ELASTICITY EN 2D WITH FREEFEM++

17

Transcript of ELASTICITY EN 2D WITH FREEFEM++

Page 1: ELASTICITY EN 2D WITH FREEFEM++

Éléments nis : Élasticité en dimension d=2.

Abdelouahab OUKHOUYA

12 juin 2012

Page 2: ELASTICITY EN 2D WITH FREEFEM++

Position du problème :

Soit Ω un ouvert régulier de de frontière Γ de classe 1 par morceaux ; soit Γ0 une partiede Γ de mesure supercielle strictement positive et soit Γ1 le complémentaire de Γ0 dans Γ.On considère un solide occupant le domaine Ω.

On note ~u le vecteur de déplacement. Si le solide est soumis à la densité de force ~f (gravitépar exemple), à la contrainte normale ~g sur une partie de son bord Γ1 et encastré sur l'autrepartie du bord Γ0, alors le vecteur déplacement satisfait aux équations de l'élasticité linéaire :

∑2j=1

∂∂xjσij(~u) + fi = 0 dans Ω

∑2j=1 σij(~u)nj = gi sur Γ1

~u = 0 sur Γ0

(1)

oú nj désigne la jième composante de la normale extérieure à Ω ; fi ∈ L2(Ω), gi ∈ L2(Γ1) aveci=1,2.(σij)1≤i,j≤N est le tenseur donné par

σij(~u) = λtrace(ε(~u))δij + 2µεij(~u) (λ et µ ∈ IR : coecients de Lamé)

εij(~u) = 12( ∂ui

∂xj+

∂uj

∂xi) (appelé tenseur des déformations)

(2)

2

Page 3: ELASTICITY EN 2D WITH FREEFEM++

0.1. ÉTUDE THÉORIQUE : 3

0.1 Étude théorique :

Cette section aborde l'approximation par éléments nis de problèmes (1).Le domaine Ω ∈ IR2 représente un milieu déformable, initialement au repos et auquel on ap-plique un chargement extérieur ~f : Ω −→ IR2.L'objectif est de déterminer ~u : Ω −→ IR2 le champ de déplacement fois que le milieu a atteintl'équilibre.

Remarque 0.1.11. Le problème (1) à des conditions limites de Van-Neumann, on peut par exemple le traité

avec des conditions de Dirichlet. Dans ce cas l'étude est plus simple.

2. Les coecients de Lamé sont des coecients λ et µ qui, pour des raisons physiques, sontcontraints par les relations µ > 0 et λ+ 2

3µ > 0.

3. Dans certaines applications, il est plus pratique d'introduire le module de Young E et lecoecient de Poisson ν tels que :

E = µ3λ+ 2µ

λ+ µ. (3)

et

ν =1

2

λ

λ+ µ. (4)

Soient les espaces L2(Ω)n et H1(Ω)n munis des normes suivantes :

‖~v‖0,Ω = (n∑

i=1

‖vi‖20,Ω)

12

et

‖~v‖1,Ω = (n∑

i=1

‖vi‖21,Ω)

12 .

Lemme 0.1.1 (Premier inégalité de Korn)Soit Ω un ouvert borné et connexe de IRn de frontière assez régulière. Alors :∫

Ω

n∑i,j=1

‖εij(~u)‖20,Ω dx ≥

1

2‖~u‖2

1,Ω ∀~u ∈ H10 (Ω)N . (5)

PreuveOn donne l'idée de la démonstration :

On suppose dans les données du problème que Γ1 = ∅1. On écrit le problème variationnel équivalent au problème (1).

2. on montre que pour tout champ ~u ∈ D(Ω)N , on a∫Ω

∂ui∂xj

∂uj∂xi

dx =

∫Ω

∂ui∂xi

∂uj∂xj

dx

3. Puis on conclut.

Page 4: ELASTICITY EN 2D WITH FREEFEM++

4

Théorème 0.1.1 (seconde inégalité de Korn)Soit Ω un ouvert borné et connexe de IRn de frontière assez régulière. Alors il existe uneconstante C(Ω) > 0 qui dépend de Ω telle que :

∀~u ∈ H1(Ω)n, C(Ω)‖~u‖21,Ω ≤

n∑i,j=1

‖εij(~u)‖20,Ω + ‖~u‖2

0,Ω. (6)

PreuveVoir P.A.Raviate et J.M.Thomas.

Le théorème précèdent nous dit tout simplement que :~u 7−→ (

∑ni,j=1 ‖εij(~u)‖2

0,Ω + ‖~u‖20,Ω)

12 est une norme équivalente à la norme : ~u 7−→ ‖~u‖1,Ω sur

H1(Ω)n.Soit V = ~u ∈ H1(Ω)N / ~u/Γ0 = 0.

On pose < = ~u ∈ H1(Ω)n; εij(~u) = 0 pour toute i, j ∈ 1, · · · , n.On a

~u ∈ < ⇔ ~u(x) = Ax+b, où A est une matrice carrée d'ordre n antisymétrique et b un vecteur de IRn.

Dans notre cas n = 2 :

~u(x) = (u1(x), u2(x)) = (a1 + bx, a2 − bx).

oú a1,a2 et b sont des réels.Si l'ouvet Ω est connexe ; alors

< = ∃a1, a2, b ∈ IR;u1(x) = a1 + bx, u2(x) = a2 − bx

Par suite ~u = 0 ; or on a supposé que Γ0 est mesure supercielle strictement positive donc :

V⋂< = 0.

Donc ~u 7→ (‖εij(~u)‖20,Ω)

12 est une norme sur V .

De plus par la première inégalité de Korne (5) on a :∫Ω

n∑i,j=1

‖εij(~u)‖20,Ω dx ≥

1

2‖~u‖2

1,Ω ∀~u ∈ V.

Donc ~u 7→ (‖εij(~u)‖20,Ω)

12 est équivalente à la norme ‖.‖1,Ω. On na le théorème :

Page 5: ELASTICITY EN 2D WITH FREEFEM++

0.1. ÉTUDE THÉORIQUE : 5

Théorème 0.1.2Soit Ω un ouvert borné et connexe de IRn, de frontière Γ assez régulière. Alors il existe uneconstante C0 strictement positive telle que :

n∑i,j=1

‖εij(~u)‖20,Ω ≥ C0‖~u‖2

1,Ω ∀~u ∈ V. (7)

De ce théorème on va montrer la V-ellipticité de la forme linéaire associée à la formulationvariationnelle de problème (1).

Formulation variationnelle et solution :

On pose

a(~u,~v) =2∑

i,j=1

∫Ω

σij(~u)εij(~v) dx. (8)

avec, σij(~u) et εij(~v) sont donnés dans l'équation (2).Par (2), on a :

a(~u,~v) = λ

∫Ω

(∂u1

∂x1

+∂u2

∂x2

)(∂v1

∂x1

+∂v2

∂x2

) dx+ 2µ2∑

i,j=1

∫Ω

εij(~u)εij(~v) dx (9)

par cette relation on a : a(~u,~v) = a(~v, ~u), donc a(., .) est symétrique.

De plus d'après le théorème (0.1.2) : il existe C0 > 0, telle que

a(~v,~v) ≥ 2µC0‖~u‖21,Ω ∀~u ∈ V. (10)

Il en découle que a(., .) est V-elliptique.Il est claire que a(., .) est continue sur V 2.

D'autre part on a la forme :

~u 7−→∫

Ω

(f1u1 + f2u2) dx+

∫Γ1

(g1u1 + g2u2) dσ (11).

est une forme linéaire continue sur V .

Conclusion :Vu les hypothèses données donc par le théorème Lax-Milgram du cours : il existe ~v ∈ V uniquesolution de

a(~u,~v) =

∫Ω

(f1u1 + f2u2) dx+

∫Γ1

(g1u1 + g2u2) dσ ∀~u ∈ V. (12)

Or la forme bilinéaire est symétrique, alors on a la fonction de l'énergie dénie par :

J(~u) =1

2

n∑i,j=1

∫Ω

σij(~u)εij(~v) dx−∫

Ω

(f1u1 + f2u2) dx−∫

Γ1

(g1u1 + g2u2) dσ. (13)

D'après un autre théorème du cours, la solution ~v de (12) est donner par :

J(~v) = minJ(~u); ~u ∈ V . (14)

Page 6: ELASTICITY EN 2D WITH FREEFEM++

6

Régularité de la solution :

Soit ~v ∈ H2(Ω)2.On a σij = σji, donc : a(~u,~v) =

∑2i,j=1

∫Ωσij(~v) ∂ui

∂xjdx (15).

Si ~u une fonction test (i.e ~u ∈ D(Ω)), alors par (12) et (15) :

2∑j=1

∂xjσij(~u) + fi = 0dans Ω, i = 1, 2 (16).

donc la formulation variationnelle (12) devient

a(~u,~v) = −2∑

i,j=1

∫Ω

(∂

∂xjσij(~v))ui dx+

2∑i=1

∫Γ1

giui dσ ∀~u ∈ V. (17)

Par la formule de Green on a,

2∑i,j=1

∫Ω

σij(~v)∂ui∂xj

dx =2∑

i,j=1

∫Ω

(∂

∂xjσij(~v))ui dx+

2∑i,j=1

∫Γ

(σij(~v)uiuj) dσ. (18)

Par suite (17) s'écrit,

2∑i=1

∫Γ1

(2∑

j=1

σij(~v)uj − gi)ui dσ = 0, ∀~u ∈ V. (19)

Finalement puisque on supposé depuis le début que la solution ~v ∈ H2(Ω)2, alors

2∑j=1

σij(~v)uj = gi sur Γ1, i = 1, 2 (20).

Par suite on a la régularité de la solution faible.

Page 7: ELASTICITY EN 2D WITH FREEFEM++

0.2. ÉTUDE PRATIQUE 7

0.2 Étude pratique

On considère un domaine Ω = ]0, 4[x ]0, 1[ un ouvert de IR2. Soit Γi, i = 1, 2, 3, 4 les cotêsdu rectangle Ω.On introduit un maillage régulier τh et Wh l'espace des fonctions P1 à valeurs vectorielles (dansIR2) associée.

Premier cas :

On va chercher à résoudre le problème de l'élasticité linéaire dans une poutre de longueur 4et de largeur 1, xée à son extrémiste gauche, soumise à des forces ~f .

D'aprés la partie théorique l'approximation du problème de l'élasticité linéaire consiste àdéterminer :

~uh = [u1h, u

2h ]∈ Vh = [w1

h, w2h

]∈ Wh : w1

h = 0 et w2h = 0 sur Γ4.

tel que ∫Ω

2µe(uh) : e(vh) + λ(Ouh)(Ovh) dx =

∫Ω

f.vh dx. (21)

pour tout vh ∈ Vh ;oú e(vh) est le tenseur linéarisé,

e(vh) =1

2(Ovh + (Ovh)>) (22).

Page 8: ELASTICITY EN 2D WITH FREEFEM++

8

Remarque 0.2.1L'achage de données vectorielles à l'aide de la fonction plot n'est pas idéal.Par contre, FreeFem++ ore la possibilité de visualiser la déformation d'un maillage τh.Dans un premier temps, on dénit le maillage déformé par :real ex=0.05 ;//coecient d'exagerationmesh Sh=movemesh(Th,[x+ex*uh1,y+ex*uh2]) ;//on déforme ou déplace le mesh

Il sut alors d'acher le maillage déformer Sh.plot(Sh) ;

Code :

Page 9: ELASTICITY EN 2D WITH FREEFEM++

0.2. ÉTUDE PRATIQUE 9

Figures :

Page 10: ELASTICITY EN 2D WITH FREEFEM++

10

Commentaires :On peut améliorerai le calcul précédent en travaillant sur un maillage de plus en plus n.Cependant, utiliser un maillage uniforme n'est pas optimal.Il plus intéressant de raner dans les régions oú il se passe réellement quelque chose,dans notrecas c'est les bords de Γ4, c'est à dire oú le gradient de la solution uh varie rapidement.Pour cela en utilise la fonction adaptmesh, qui permet de raner le maillage en l'adaptant àune fonction spéciée :real erreur=0.001 ;Th=adaptmesh(Th,uh1,uh2,err=erreur) ;permet d'adapter le maillage avec une précision inversement proportionnelle à l'erreur en fonc-tion de uh. Alors par ce possédé on améliore notre solution comme on le va montrer.

Code :

Page 11: ELASTICITY EN 2D WITH FREEFEM++

0.2. ÉTUDE PRATIQUE 11

Figures :

Page 12: ELASTICITY EN 2D WITH FREEFEM++

12

Page 13: ELASTICITY EN 2D WITH FREEFEM++

0.2. ÉTUDE PRATIQUE 13

Deuxième cas :

Cette fois, la poutre xée à son extrémiste gauche et droite.

Ce cas se traite comme le cas précédente ; sauf qu'on ajoute dans la résolution du problèmevariationnelle, une autre condition limite.Dans ce cas on trouve les singularités au bord de Γ2 et Γ4.Codes :

Page 14: ELASTICITY EN 2D WITH FREEFEM++

14

Page 15: ELASTICITY EN 2D WITH FREEFEM++

0.2. ÉTUDE PRATIQUE 15

Figures associées

Page 16: ELASTICITY EN 2D WITH FREEFEM++

16

Page 17: ELASTICITY EN 2D WITH FREEFEM++

0.2. ÉTUDE PRATIQUE 17