Sommaire …………………………………………………………………………………………………….
1
Introduction……………………………………………………………………………………..............3
I-Détermination de l’équation d’équilibre par la méthode vibrationnelle…….6
1-Equation d’équilibre………………………………………………………………………….………6
2-Détermination du système à résoudre………………………………………………….…..7
II-Programmation sur
matlab……………………………………………………………………………………………….…………
9,10,11,12
III-modèle numérique et résultat du calcul………………………………….……………….13
Conclusion……………………………………………………………………………………………………15
Bibliographie………………………………………………………………………………………………..16
1
SCHEMA
Formulation du modele variationnelle ………………………………………………………………………………4
2
INTRODUCTION
La mécanique non linéaire désigne le comportement des systèmes mécaniques dont la
réponse n’est pas proportionnelle à la force appliquée ce qui désigne que les lois de Hooke
(déformation proportionnelle à la force) ne s’appliquent plus. Selon les conditions limites aux liées
aux appuis et encastrement que peut être la contrainte σ et le module de Young μ dans une plaque
grande ponctuation. A l’aide de la méthode vibrationnelle et la méthode des éléments finis puis la
programmation sur Matlab il s’agira tout au long de notre devoir de calculer puis trouver la
contrainte σ et le module de Young μ.
3
SCHEMA
FORMULATION DU MODELE VARIATIONNELLE
4
5
I-Détermination de l’équation d’équilibre par la méthode variationnelle :
1-Equation d’équilibre :
div T0 = f dans Γ (1)
T0x= P sur Γ1 (2)
U0=0 sur Γ3
T0=A (x, F) =Fσ(x, F)
F=I+𝛁U de (I) on a :
div T0 = f <=> ∫ fvdx=∫ ¿ Tv
On pose :
On pose ∇=JU on a :
6
On pose : σ =a ( x ) I +b ( x ) +c ( x ) E avec ;p=2 KN/m2 , u=0 sur Γ0
∫ a ( x ) I + b ( x ) trEI + c ( x ) E=¿ ¿
Or a ( x ) I + c ( x )=0 d’où
2-Détermination du système à résoudre
On pose :
D’où σ =f(ε) en élasticité non linéaire .
σij=(σ1 ,σ2,σ3,σ4,σ5,σ6)=
= (σ11 ,σ22,σ33,σ44,σ55,σ66)
= (σ11 ,σ22,σ33,σ44,σ55,σ66)
ε ij=(ε 1 ,ε 2,ε 3,ε 4,ε 5,ε 6)=
=(ε 11 ,ε 22,ε 33,ε 44,ε 55,ε 66)
On a :
et avec
On a : e1=cosθe 1+ sinθe2 ; e 2=-sinθe 1+ sinθe2
σ11=λ(ε 11 +ε 22+ε 33)+ 2με 11
σ22=λ(ε 11 +ε 22+ε 33)+2με 11) + 2με 22
σ33=λ(ε 11 +ε 22+ε 33)+2με 11) + 2με 33
σ12=2με 12
σ13=2με 13
σ23=2με 23
7
Par conséquent notre équation à résoudre devient :
8
II-PROGRAMMATION SUR MATLAB :
9
10
11
12
III-MODELE NUMERIQUE ET RESULTAT DU CALCUL
13
14
CONCLUSION
Tout au long de notre devoir il était question de déterminer la contrainte σ et le module
d’Young μ. L’analyse par la mécanique non linéaire nous a permis de retrouver un système
d’équation complexe en passant par la formulation variationnelle et matricielle jusqu’à la
programmation sur Matlab et en fin de compte nous avions trouvé la contrainte σ et le module de
Young μ.
15
Bibliographie :
Optimisation non linéaire : Stéphane Mottelet, Mohamed Elbagdouri.
Analyse non linéaire en grande rotations et grands déplacements des structures mince métalliques
de AHMED RAHEM.
Géométrie différentielle de surface du Professeur Robert NZENGWA, premier Doyen de l’ex-Faculté
de Génie industriel de l’université de Douala devenue école nationale supérieur de polytechnique de
Douala.
16
