Initiation à l'élasticité/Déformations

De testwiki
Version datée du 22 août 2023 à 12:52 par imported>Crochet.david.bot (Robot : remplacement de texte automatisé (-\n(==={0,3})(?: *)([^\n=\s]+)(?: *)\1(?: *)\n +\n\1 \2 \1\n))
(diff) ← Version précédente | Version actuelle (diff) | Version suivante → (diff)
Aller à la navigation Aller à la recherche

Modèle:Chapitre

Cinématique: étude du mouvement d'un corps (x,x˙,x¨,...). La cinématique est différente de la cinétique qui est l'étude des causes du mouvement.

Déplacement: Seule quantité que l’on peut mesurer.

Remarque

On ne mesure pas les forces, les contraintes ni les déformations.

[[Déformation: ϵ=ΔLL à généraliser.]]

Après déformation la nouvelle longueur est L=L0+δL

ϵ=ΔLL On pose en général L=L0

ϵ=ΔLL0+ΔL

Or on suppose ΔL<<<<LO

Remarque : L et Lo sont mesurés dans la direction d'application des forces.

Sinon :

On ne sait pas calculer ϵ

Généralisation en 3D

I.Cinématique

1.Position

Soit un point de référence, par exemple l'origine du repère, et soit un point M de l'espace occupé par un volume élémentaire de référence P à l'instant t.

La position de P=OM=(xyz).(xθz)

2. Mouvement et déplacement

Au cours d'un déplacement un VER occupe une infinité de positions.

On considère seulement:

  • l'état initial
  • l'état final

Le mouvement est causé par l’application d'une force.

Un déplacement est la différence entre les positions finales et initiales des VER.

On peut définir le déplacement du VER au cours du mouvement par:

u(P)=OMOM0=M0m

Remarque : On définit toujours la géométrie et on référence toujours les points du système dans sa position initiale.

u(d2,θ,z)=(Δd200)

3. Approximation de l'élasticité classique

Il y a beaucoup d'approximations mécaniques.

  • On considère un système statique (accélération nulle)

a=γ=0MV,(V=cste)

  • On considère des petits déplacements HPP: Hypothèse des Petites Perturbations.
Échec de l’analyse (fonction inconnue « \begin{matrix} »): {\displaystyle  \begin{matrix} d = {{unité|10|cm}} \\ si L = {{unité|10|m}} \end{matrix} \, }

Cas de petites perturbations:

-Béton

-Bois

-Les métaux (sauf domaine plastique)

On tolère Modèle:Unité maxi par mètre de poutre. La variation de la distance entre 2 points de la poutre est de quelques micromètres.

Cas où l’on est pas dans l'hypothèse des petites perturbations:

matériaux concernés :

- plastiques (mais pas tous) - liquides et fluides - domaine plastique des métaux

On a une linéarité géométrique. Par exemple une section très fine qui fait de petits efforts entraine une grande déformation (cas de géométrie particulière) Dans le cas des HPP: le déplacement u(M) est le champ de v déplacement infinitésimal à partir de l'état initial.

L'état final est infiniment proche de l'état initial.

L'intérêt est que l’on peut appliquer des règles de différenciation.

4.Tenseur de déformation

CINEMATIQUE POUR L’ELASTICITE ; LES DEFORMATIONS

Intro : « cinématique »

Décrire le mouvement d’un corps

→ cinématique du point → cinématique du solide → cinématique pour l’élasticité

≠ cinétique = (étude des causes du mouvement (les forces + PFD)) → prendre en compte le fait que le corps se déforme

mouvement solide (peut se déformer) mouvement rigide (pas de déformation)

On définit un état/configuration de référence.

Si la configuration change alors il y a déformation, distorsion ( on verra plus tard…) suppose qu'est dans la même direction que Final ? ou Initiale ?

RDM = L final → L initial → petits déplacements

Cinématique des transformations finies ne suppose pas que l’état final est proche de l’état initial (fluide - fabrication)

Pb de définition = dans une seule direction.

On veut décrire les déformations d’objet à géométrie complexe donc cette définition n’est pas suffisante.

Position et déplacement

Position

3D

M(x,y,z)
(r,θ,z)
Ajouter une légende ici

repère orthonormé

cartésien (x,y,z)
curviligne (r,θ,z)

La position de M est le vecteur OM où c’est l’origine du repère.

Cartésien

OM |xyz|

Cylindrique

OM |r00| OM=r e r M|rθz|


Notation

rp=OM

Vecteur déplacement

Soit une configuration de référence où le volume de matière considéré occupe le point M. Ce même volume occupe le point M’ pour la configuration final.

On compare seulement la position finale et initiale, on n’a pas besoin du temps.

[u(M)=OMOM]

Approximation de l’élasticité classique

On considère que le système est statique tout le temps, pour la configuration initiale, pour la configuration finale mais aussi pour toutes les étapes intermédiaires. (Il est quasi statique pour les configurations intermédiaires)

- On applique ce que dit l’approximation des « petites déformations » qui devrait être appelée « petits déplacement » :

Mathématiquement, l’approximation des petits déplacements veut dire que u(M) est le champ de vecteur déplacement infinitésimale à partir de l’état initial (OMOM )

OM |xyz| u(M)= |dxudyudzu|

Variation des déplacements autour d’un point

Variation du déplacement autour d’un point

Tenseur des déformations.

Ajouter une légende ici

B est très proche de A

AB |dxdydz| A,BV

A,BV

AA

BB On suppose que u(A) est connu. Peut-on déterminer u(B) connaissant u(A) ? [u(B)=u(A)+du(A)] u(x,y,z)
du=(dux=uxxdx+uxydy+uxzdzduy=uyxdx+uyydy+uyzdzduz=uzxdx+uzydy+uzzdz) du=(uxxuxyuxzuyxuyyuyzuzxuzyuzz) (dxdydz) [dy=grad u dx] ou OM=(dxdydz) =dx

u(B)=u(A)+grad(u(A))dx On définit : ϵ~=Sym(gradu) : tenseur des déformations

ω~=Antisym(gradu) : tenseur rotation grad(u)=12(gradu+gradu T)+12(gradugradu T) [ϵ~=12(gradu+gradu T)]

ω~=12(gradugradu T)

gradω~=ϵ~+ω~

u(B)=u(A)+(ϵ~+ω~)dx ϵ~(uxx12(uxy+uyx)12(uxz+uzx)uyy12(uyz+uxy)uzz) On a : u(B)=u(A)+ϵ~ dx+ω~ dx Pourquoi ω~ dx correspond à une rotation de corps rigide ?

ω~(012(uxyuyx)BA0CBC0) 3 valeurs scalaires suffisent pour écrire ω~.

→ un vecteur suffirait pour décrire une matrice antisym.

Soit ω~ le vecteur dual de ω~ (antisymétrique)

ω~=12 rot u On peut montrer que ω~dx=ω~1dx (cf Devoir)

u(B)=u(A)+(ϵ~+ω~)dx La décomposition du déplacement devient :

[u(B)=u(A)+ω~1dx]+ϵ~dx Soit le torseur des petits déplacements :

(ωu)A

u(B)=u(A)+ω1AB Un déplacement c’est un mouvement rigide

(translation + rotation) est une déformation.

Devoir

Q=(1a0a10001)

P : orthogonal ||P||=1=1a2

du0a0

Analyse du tenseur déformation

Exemple en traction pure

ϵ=ΔLL

Ajouter une légende ici

Champ de déplacement

U(x)(x, y, z)=ΔLLx

U(y)(x, y, z)=0

U(z)(x, y, z)=0 La façon la plus simple de décrire le déplacement de la barre ϵ~=(ΔLL00000000)

ϵxx=ΔLL

Comment analyser le tenseur

Quelles informations pratiques peut-on obtenir avec ϵ~ ?

→ Variation de longueur ?
→ Variation de volume ?
→ Variation d’angle ?
(→ Variation de surface ?)

Filière de matière :

de direction n (vecteur unitaire)

Ajouter une légende ici

A et B points V infiniment proche l’un de l’autre. AB=ABn en MV on peut définir une infinité de fibre de direction n Remarque :

→ surface dS contrainte
→ fibre
Ajouter une légende ici

Résultats mathématiques :

Ajouter une légende ici

BBAA BA+ABABBA

=ABAB du=UBUA

=BBAA

=ABAB

Variation de longueur à partir de vecteur déformation et sa décomposition en allongement et glissement

dAB=ABAB

dAB=d(ABn)=d(ABn)+(ABdn) D’autre part donne dAB=du=(ϵ~+ω~)dx=ABn

avec toujours dAB=dx=ABn donc [dABn+ABdn=(ϵ~+ω~)ABn] contient à la fois une déformation et un mouvement rigide On donne par AB dABABn+dn=(ϵ~+ω~)n dn=dnϵ~+dnω~ On regarde le mouvement rigide et la déformation : {dABABn+dnϵ~=ϵ~ndnω~=ω~n On pose : D=ϵ~n : vecteur déformation (ϵ~)

g=dnω~ : le glissement

ϵ=dABAB : l'allongement relatif D=ϵn+g ϵ~n=ϵn+g Remarque :

n.(ϵn+dnϵ~)=ϵ~n.n

(ϵ~n.n+dnϵ~n)=(ϵ~n)n [ϵ=(ϵ~n).n] On peut ainsi obtenir ϵ=dABAB en fonction de ϵ~ Interprétation de ϵn+g=D=ϵ~n

g=dnϵ~ rotation due à la déformation

Ajouter une légende ici

D=ϵn+g Remarque : parallèle avec σ~

σ~ :

C=σ~n
σ=C.n
τ=C.σn ϵ~ :

D=ϵn
ϵ=D.n
g=Dϵn

Signification des termes de ɛ̃

Si on cherche l’allongement (relatif) et le glissement dans une direction n qui correspond au vecteur de base du repère. n|ijk| n|100| n|010| n|001| ϵ=(ϵ~n).n=D.n Soit ϵ(ϵxxϵxyϵxzϵyyϵyzϵzz) ϵ~n=ϵ~|100|=|ϵxxϵxyϵxz|=D(Mi) ϵ=D.n=|ϵxxϵxyϵxz|.|100|=ϵxx n=jϵj=ϵyy n=kϵk=ϵzz Les termes diagonaux de ϵ~ correspondent aux allongements relatifs dans la direction des vecteurs de la base dans laquelle est exprimé ϵ~ . gj=Dϵn=|ϵxxϵxyϵxz|ϵxx|100|=|0ϵxyϵxz|

Variation de volume

Comment peut-on obtenir Δ V à partir de ϵ~ ?

Soit dV centré sur M

dV dV' (d(dV)=dv'-dv) Soit Δ V la variation de volume entre V et V’

V: initial

V': final avec V et V’ des éléments de vol de référence (cad infiniments petits)

Soit ϵ~ en M et ϵ~ dans son repère principal. ϵ~=(ϵI000ϵII000ϵIII) ϵI,ϵII,ϵIII déformations principales. Soit le volume V

Ajouter une légende ici

V=dX dY dZ
V=dX dY dZ dX en fonction de dX dY dZ

dX=dX+ϵIdX
dY=dY+ϵIIdY
dZ=dZ+ϵIIIdZ

V=(1+ϵI)(1+ϵII)(1+ϵIII)(dx dy dz) V=V+V(ϵI+ϵIIϵIII)+V(ϵIϵII+ϵIϵIII+ϵIIϵIII)+V(ϵIϵIIϵIII) ϵ0 ϵIϵII0 0(ϵIϵII) << 0(ϵI) VVV=ϵI+ϵII+ϵIII [dVV=tr(ϵ~)] Variation de volume local

Pour avoir la variation d’un volume V : ΔV=v tr(ϵ~)dV ϵ~ΔL ϵ~ΔV

Variation d’angle à partir de ɛ̃

Ajouter une légende ici

ϵMN = l’allongement relatif de la fibre MN


ϵMP = l’allongement relatif de la fibre MP

ϵ~ le tenseur des déformation en M

|nMNnMP| Vecteur unitaire qui représente la fibre


d φ=(ϵMN+ϵMP)cosφ2ϵ~nMNnMPsinφ ϵMN=(ϵ~.nMN).nMN Distorsion pour les directions des vecteurs de base à ϵ~ connu. Soit ϵ~(M) et on veut d φ Pour n1|100| et n2|010|

Ajouter une légende ici

φ=π2

sin φ=1/cos φ=0 Donc d φ=2ϵ~ n1 n2
=2(ϵ~ n1).n2
=2(ϵxxϵxyϵxz)n2
=2ϵxy Si d φex et ey=2ϵxy Les termes non diagonaux du tenseur des déformations donnent la variation d’angle des vecteurs de base du repère.
Conclusion analyse de ϵ~ (cf poly TD)

ϵ~=12(grad u+grad uT)

D=ϵ~n=ϵn+g

ϵ=Dn=(ϵ~n).n

g=Dϵn ϵ~ est symétrique donc il admet 3 déformations principales et des directions principales et invariants.

Avec ϵ~ on peut calculer la variation de longueur Δll=ϵ Calculer une variation de volume ΔVV Trace (ϵ~) Calculer une variation d’angle d φ=... Calculer une variation d’aire

Condition aux limites en déplacement

Soit S la surface de volume, certains points de S ont des déplacements connus (ou imposé) .

Ajouter une légende ici

Surface
x=0

n(100) C
?
u
0 Exemple : un objet posé sur une table

Ajouter une légende ici

Quel frottement ? = pas de frottement

Si on néglige strictement la gravité, pas de CL sur la surface | z = 0 Souvent on ne prend pas en compte la gravité dans l’équilibre, mais on la considère pour les CL. Surface

z = 0

n(001) C
u

Ux=?

Uz>0

Ajouter une légende ici

CL au niveau du contact ? Pour tout point tel que usphere=umassif Difficultés pour exprimer les CL

Modèle:Bas de page