Projet Mef

UNIVERSITE DES SCIENCES ET TECHNIQUES DE MASUKU
ECOLE POLYTECHNIQUE DE TECHNIQUES DE MASUKU
ÉC O LE
P O LYTEC HN IQ UE
D E M ASUKU
Département de Génie Civil
Projets de
METHODE
DES
ELEMENTS FINIS
Présenté par :
- BOHILE Sophie
- EYEYE NGANGA Vianney
- NZENGUE Claille Darlin
- LIPOBO MINGONDZA Sarah Anne
- MONIKA NGOUNDI Cruz Lunel
- PAMBOU MALOUBA Lenny Van
- PENDI MOUELE Ulrich
- ZINGA PAMBOU Greg Pulgervie
Dirigé par :
Pr MOUTOU PITTI Rostand, enseignant à l’École Polytechnique de Masuku
PROJETS DE METHODE DES ELEMENTS FINIS
TABLE DES MATIÈ RES
PROJET 1 : RÉ SOLUTION D’UN PROBLÈ ME D’É LASTICITÉ __________________________________2

I- Résolution par la méthode de la Théorie des poutres (RDM)_____________________2
II- Application du principe des travaux virtuels (PTV)____________________________4
III- Solution trigonométrique_________________________________________________11
PROJET 2 : FORMULATION DE LA MEF EN DÉ PLACEMENT________________________________14
Exemple 1 : Traitement en parallèle du problème d’une barre élastique verticale soumise à
son poids propre________________________________________________________________14
Exemple 2 : Maillage de la barre____________________________________________________15
Exemple 3______________________________________________________________________24
PROJET 3 : MÉ THODE DES É LÉ MENTS FINIS : É TUDE D’UN É LÉ MENT DE POUTRE DE
BERNOULLI_____________________________________________________________________________________25
I. Donnons l’expression du champ axial u(x) en tout point de l’élément en fonction des
valeurs nodales u1 et u2.___________________________________________________________25
II. Déterminons le degré du champ v pour lequel nous pouvons représenter un champ du
moment fléchissant linéaire_______________________________________________________26
PROJET 4 : MÉ THODE DES É LÉ MENTS FINIS : É TUDE DES É LÉ MENTS FINIS
TRIANGULAIRES________________________________________________________________________________29
I- É tude des éléments finis à trois nœuds__________________________________________29
II- É tude des éléments finis triangulaires à 6 nœuds_________________________________38
1
École Polytechnique de Masuku _ Ingénieur 3 GC _ 2023-2024
PROJET 1 : RÉSOLUTION D’UN PROBLÈME D’ÉLASTICITÉ
Exemple d’utilisation de la méthode de Ritz : Étude d’une poutre
x u
→ →
R1 R2 y v
L
z
h
y
Soi
t une poutre uniformément chargée sur toute sa longueur et reposant sur deux appuis.
I- Résolution par la méthode de la Théorie des poutres (RDM)
Déterminons l’expression de la flèche v ( x ) et du moment fléchissant M(x).
- Calcul des réactions aux appuis
D’après le principe fondamental de la statique :
∑ F H =0 → H =0
1
∑ M A =0 → R 2 L+ 2 q L2=0
−1
R 2= qL
2
∑ F V =0=¿ R 1+ R 2=−qL
1
R1=−qL+ qL
2
−1
R 1= qL
2
2
- Calcul des efforts internes
Considérons la partie gauche de la coupure : 0 ≤ x ≤ L
-N(x) +
M(x)
T(x)
R1
1 1
∑ M ∑=0=¿ M ( x )+ 2 q x2 − 2 qx=0
1
M ( x ) = qx (L−x)
2
Équation de la déformée d’après TIMOSCHENKO

''
EI V ( x )=−M (x )
'' −1
V ( x )= M (x )
EI
'' −q
V ( x )= x (L−x )
2 EI
' −q 2 q 3
V (x )= x+ x +α
4 EI 6 EI
−qL 3 q 4
V ( x )= x+ x + αx+ β
12 EI 24 EI
Le champ étant cinématiquement admissible, on peut alors appliquer les conditions aux
limites.
 V ( 0 )=0=¿ β=0
4 4
−q L qL
 V ( L )= + +αL=0
12 EI 24 EI
3
qL
→α=
24 EI
3
q 4 qL 3 q L
→ V ( x )= x − x + x
24 EI 12 EI 24 EI
q 4 3 3
V ( x )= (x −2 L x + L x)
24 EI
()
4
L L qL
Pour : x= →V ≈ 0.0130
2 2 EI
3
II- Application du principe des travaux virtuels (PTV)
On cherche maintenant à obtenir une solution approchée de ce problème par la méthode

de Ritz (cela serait le cas si l’on ne savait pas déterminer la solution analytique. Le fait de
connaitre ici cette solution, nous permettra d’apprécier la véracité de notre
approximation). Pour cela, on se propose d’utiliser le PTV.
δ W i+ δ W e=0 ∀ δu ϵ C . A
 δ W e =∫ qδv ( x ) dx
0
Le matériau étant élastique linéaire et isotrope, nous pouvons alors écrire la loi de
HOOKE réduite : σ x =E ε x
−N M
Par ailleurs, σ = + y comme la flexion est la seule composante prise en compte, on
A I
M
a : σ= y
I
M
Connaissant l’expression de σ on a y=E ε xy
I
Or M ( x ) =−EIV ' '( x )
''
→−E V ( x ) y=E ε xy
→ ε xy =− yv ' '( x)
2
d v (x )
→ E ε xy =−E y
d x2
1
d 2 v (x ) d 2 δv ( x )
 δ wi=∫ (−Ey 2
¿)(− y 2
)dV ¿
V dx dx
On sait que : dV =dSdx
2
d v ( x) 2
L L ∗d δv ( x )
d x2
→ δ W i =EI ∫ y dS ∫ (¿
2
) dx ¿
0 0 d x2
4
2
d v (x) 2
L ∗d δv ( x )
d x2
δ W i=EI ∫ (¿ )dx ¿
0 d x2
on a : −δ W i=δ W e
2
d v(x) 2
L 2
∗d δv ( x ) L
dx
→ δJ=EI ∫ (¿ )dx =∫ qδv ( x ) dx ¿
0 d x2 0
Solution approchée :
La courbe de la déformée étant une parabole, on peut donc écrire :
V ( x )=α x 2 + βx +γ
Conditions aux limites :
 v ( 0 )=0 → γ =0
 v ( L )=0→ β=−αL
2
→ V ( x )=α x −αLx
2
δV ( x ) =δα (x −Lx)
L L
∫ qδα (x 2¿−Lx )dx =qδα ∫ (x 2 ¿−Lx)dx ¿ ¿

0 0
1 3 1 3
¿ qδα ¿ ¿ qδα ( L − L )
3 2
∫ qδv ( x ) dx= −1
6
3
qδα L 1¿ ¿
0
( )
L 2 2 L
d v ( x ) d δv ( x )
∫ d x 2 d x 2 dx=∫ 2 α∗2 δαdx
0 0
( )
L 2 2
d v ( x ) d δv ( x )
∫ d x 2 d x 2 dx=4 αδαL 2¿
0
3
qL
δ W e =−δ W i=¿− δα=EI 4 Lαδα
6
5
2
−q L
Par identification : α =
24 EI
2 2 3
−q L x q L x
D’où v ( x )= +
24 EI 24 EI
2
( ) qL
v x= x (L−x)
24 EI
2
d v ( x ) dv ( x ) q L3 q L2
M ( x ) =−EI ∨ = − x
d x2 dx 24 EI 12 EI
M ( L2 )=0.083
M RDM ( L2 )=0.125
Err =
M RDM ( L2 )−M ( L2 ) = 0.125−0.083
M RDM( L2 ) 0.125
Err =33.6 %
Divisons la poutre en 4 éléments de même longueur
1 2
2
É lément 1 : v 1 ( x )=α 1 x + β 1 x +γ 1
2
É lément 2 : v 2 ( x )=α 2 x + β 2 x +γ 2
On peut appliquer les conditions aux limites :
- Condition en déplacement :
v 1 ( 0 ) =0 → γ 1=0
v1 ( L4 )=v ( L4 ) → 161 α L + 14 β L= 161 α L + 14 β L+ γ 1 ¿

2 1
2
1 2
2
2 2
- Continuité en rotation :
6
v'1 ( L4 )=v ' ( L4 )=¿ 12 α L+ β = 12 α L+ β

2 1 1 2 2
1
→ L(α ¿ ¿ 1−α 2 )+ β 1−β 2=0 2 ¿
2
v'2 ( L2 ) → α L+ β =0 3 ¿
2 2
On obtient le système suivant :

1 2 1
L ( α 1 −α 2 ) + L ( β 1−β 2 )−γ 2=0 1 ¿
16 4
1
L(α ¿ ¿ 1−α 2)+ β 1−β 2=0 2 ¿
2
α 2 L+ β 2=0 3 ¿
D’où 3 ¿ → β 2=−α 2 L
1 −1
2 ¿ → L ( α 1−α 2) + β1 +α 2 L=0=¿ β 1= L(α 1 +α 2)
2 2
1 2 1 2
1¿→ L ( α 1−α 2) + L ( α 1−α 2−2 α 2 )−γ 2=0
16 8
−1 2
→ γ 2= L ( α 1−α 2 )
16
On a alors :
2 1
v 1 ( x )=α 1 x − L( α 1 +α 2) x
2
2 1 2
v 2 ( x )=α 2 x −α 2 Lx− L (α 1 −α 2 )
16
Application du PTV : δ W e =−δ W i
- É lément 1) :
7
L
4
 δ W =∫ qδ v ( x ) dx
e1 1
0
L
4
−1
¿ q ∫ [δ α 1 x 2 ¿ L ( δ α 1 +δ α 2 ) x ]dx ¿ ¿ q ¿
0 2
¿q¿
3 3
−L L
¿q( δ α 1− δ α 2)
96 64
3 δ α1 δ α2
δ W e1 =−q L ( + )
96 64
L
4 2 2
 δ W =−EI ∫ d v 1 ( x ) ( d δv1 ( x ) )dx
i1
0 d x2 d x2
L
4
¿−EI ∫ (2 α 1)(2 δ α 1 ¿ )dx ¿

0
L
4
¿−4 α 1 δ α 2 EI ∫ dx
0
δ W i 1=−α 1 δ α 1 EIL
- É lément 2) :
L
2
 δ W = ∫ qδ v ( x ) dx
e2 2
L/ 4
L
2
1 2
¿ q ∫ [δ α 2 x 2¿−δ α 2 Lx− L ( δ α 1 +δ α 2 ) ]dx ¿
L/ 4 16
¿q¿
¿q¿
3 3
−L L
¿q( δ α 2− δ α 1)
24 64
3 δ α 2 δ α1
δ W e2 =−q L ( − )
24 64
L
2
 δ W i 2=−EI ∫ (2 α 2)(2 δ α 2 ¿ )dx ¿

L
4
8
L L
−4 EI α 2 δ α 2 ( − )
2 4
δ W i 2=−EI α 2 Lδ α 2
PTV :
−q L3 ( −δ α 1 δ α 2
96
+
64
−q L3 ) (
−δ α 2 δ α 1
24
+
64
=EI ¿ )
−5 3 11 3
L qδ α 1− L qδ α 2=EI α 2 Lδ α 2 + EI α 1 δ α 1
192 192
Par identification :
−5 3
α 1= L q
192 EI
−11 2
α 2= L q
192 EI
 β 2=−α 2 L
2
11 L
→ β2 = q
192 EI
3 2
−L −L −5 L 11 L
 β 1= (α 2 +α 2)= ( q− q)
2 2 192 EI 192 EI
3
L
β 1= q
24 EI
2 2 2 2
−L L −5 L 11 L
 γ 2= (α 1−α ¿ ¿ 2)= ( q+ )¿
16 16 192 EI 192 EI
4
−L
¿> γ 2 = q
512 EI
2 3
−5 L 2 L
D’où ( )
v1 x = qx + qx
192 EI 24 EI
2 3 4
−11 L 2 11 L L
v 2 ( x )= qx + qx− q
192 EI 192 EI 512 EI
()
2
L 5L
→ M1 = q
2 96
9
M1 ( L2 )≈ 0.052 q L 2
()
2
L 11 L
→ M2 = q
2 96
M2 ( L2 )≈ 0.114 q L 2
Calcul de l’erreur :
- Calcul de la flèche :
4 4 4 4
−11 L 11 L L L
 v 2 ( L/2 )= q+ q− q=0.0124 q
768 EI 384 EI 512 EI EI
L
V −V 2 ( )
RDM
L
2 ( ) 2
Err =
V
RDM ( L2 )
0.0130−0.0124
Err =
0.0130
Err ≈ 4.6 %
- Calcul du moment :
L
M −M 2 ( )
RDM
L
2 ( ) 2
Err =
M
RDM ( L2 )
0.125−0.114
Err =
0.125
Err ≈ 8.8 %
III- Solution trigonométrique
n
iπx
V (x )=∑ α i sin( ) 𝑃𝑜𝑢𝑟 𝑡𝑜𝑢𝑡 𝑉(𝑥) 𝐶𝐴
i=0 L
10
πx
Pour n=1, 𝑉(𝑥) = 𝛼1 sin( )
L
D’après le PTV
𝑤𝑖 +𝑤𝑒 = 0 ⇒ 𝑤𝑖 = − 𝛿𝑤𝑒
Mais −𝑤𝑖 = 𝐽
L
On a : 𝑤𝑒 = ∫ q δ V ( x )dx
0
πx
Or 𝑉(𝑥) = 𝛿𝛼1 sin( )
L
L L
πx πx
Alors : 𝑤𝑒 = ∫ q δα 1 sin ( ) d x = 𝑞𝛿𝛼1 ∫ sin ( ) d x
0 L 0 L
πx πx πx πx πx πx πx
Or sin( )= sin( + ) = sin( ¿ cos( ¿ + sin( ¿ cos( ¿
L 2L 2 L 2L 2L 2L 2L
πx πx πx
sin( )=2 sin( ¿ cos( ¿
L 2L 2L
L
πx πx
⇒ 𝑤𝑒 = ∫ q δα 1 sin ( )cos( )d x
0 2L 2L
πx πx π L πx πx
Soit sin( ¿ cos( ¿ = ∗¿ sin( ¿ cos( ¿
2L 2L L π 2L 2L
L πx πx
= ¿sin( ¿ cos( ¿
π 2L 2L
πx π πx πx
On sait que :( sin( ¿ ²)′ = 2∗ sin( ¿ cos( )
2L 2L 2L 2L
π
∗L L
D’où : 𝑤𝑒 = 𝑞1 2 L πx πx
∫ π sin ( 2 L )cos ( 2 L )d x
0
L
2L π πx πx
= 𝑞1 ∫ sin( )cos( )d x
π 0 L 2L 2L
2L π
Donc 𝑤𝑒 = 𝑞 (sin ² ( −0)
π 1 2
2L
𝑤𝑒 = 𝑞
π 1
11
πx
On a: 𝑉(𝑥) = 𝛿𝛼1 sin( )
L
2
d V (x )
L ∗d δ V ( x )
−𝑤𝑖 =EIZ d x2
∫( d x2
)d x
0
dV ( x ) π πx
= 𝛼1 cos( )
dx
2
L L
2
d V (x) π² πx
=¿ - 𝛼1 sin( )
dx 2
L² L
Par analogie
2
d V (x) π² πx
=¿ - δ 𝛼1 sin( )
dx 2 L² L
L
( ) ( )
2 2
π πx π πx
Donc −𝑤𝑖 = EIZ ∫ [−α 1 2
sin ][−δα 1 2 sin ]d x
0 L L L L
[ ( )]
4 L
π πx 1
= EIZ 𝛼1 ²δ 4∫ Or 𝑠𝑖𝑛²(𝑥) = (1−cos(2x))
sin ² dx
L 0 L 2
πx 1 πx
D’où 𝑠𝑖𝑛²( ) = (1−cos(2 ))
L 2 L
[ ]
4 L
π 2 πx
⇒ −𝑤𝑖 = EIZ 𝛼1 ²δ 4∫
1−cos ( ) dx
2L 0 L
Calcul de l’intégral
[ ]
L L L
2 πx 2 πx
∫ 1−cos (
L
) d x=∫ d x−∫ cos(
L
)d x
0 0 0
L
=L - [ (0−0)¿ = L
2π
[ ]
L
2 πx
Donc ∫ 1−cos (
L
) d x=L
0
4
π
Alors 𝐽 = −𝑤𝑖 = EIZ 𝛼1 ²δ 3
2L
Système d’équilibre
12
𝐽 = 𝑤𝑒
4
π 2L
⇒ EIZ 𝛼1 ²δ 3= 𝑞
2L π 1
Par identification :
4
4 qL
𝛼1 = 5
EI z π
Finalement
4
4 qL πx
𝑉(𝑥) = 5 𝑠𝑖𝑛( )
EI z π L
Pour 𝑥 = 𝐿 ; On a :
4
L 4 qL πL
𝑉( )= 5 𝑠𝑖𝑛( )
2 EI z π 2L
4
L 4 qL πL
Donc 𝑉( )= Car 𝑠𝑖𝑛( )=1
2 EI z π 5 2L
4
L 4 qL L 0.0130 qL4
Pour n=1, 𝑉( )= = 𝑉( )=
2 EI z π 5 2 EIz
Calcul de l’erreur
L
V −V trigo ( )
RDM(L
2 ) 2
Err =
V
RDM ( L2 )
0.0130−0.0134
Err =
0.0130
Err =0 %
13
PROJET 2 : FORMULATION DE LA MEF EN DÉPLACEMENT

Exemple 1 : Traitement en parallèle du problème d’une barre élastique verticale
soumise à son poids propre
Le système étant linéaire et la barre isotrope homogène on a :
- Efforts externes :
∑ F ¿ x =0=¿ V + LAρg=0 𝐻
→ V =− AρgL V
∑ F ¿ y =0=¿ H =0
ρg
∑ M =0=¿ C=0
𝑥
- Effort internes :
𝑀1
∑ F ¿ x =0 → N ( x )−V + Aρgx=0
x 𝑉1
→ N ( x ) + Aρgx− AρgL=0
→ N ( x )= Aρg(L−x )
T
- Calcul de la contrainte : −𝑁(𝑥)
N ( x)
→ σ ( x )= −𝑀(𝑥)
A
→ σ ( x )=ρg( L−x)
- Calcul de la déformation :
14
σ (x)
→ ε x=
E
ρg
→ ε x= (L−x)
E
- Calcul du déplacement :
du
→ ε x=
dx
→ du=ε x dx
x
→ u=∫ ε x dx
0
x
ρg
→ u( x)=
E 0
∫ (L−x )dx
2
ρg x
→ u( x)= (Lx− )
E 2
Exemple 2 : Maillage de la barre
On se propose de considérer 2 maillages différents :
- Un maillage avec 2 éléments unidimensionnels à 2 nœuds de longueur

L
(maillage 1)
2
- Un maillage à un élément unidimensionnel à 3 nœuds de longueur L
(maillage 2)
Maillage 1 Maillage 2
1 x=0 1
I
2 L I 2
x=
I 2
I 3 x=L 3
 Les champs sont cinématiquement admissibles.
15
 On a deux éléments de la barre constituées de deux nœuds ceci entraine qu’on a une
fonction de forme linéaire.
N i ( x )=α i+ β i x
 É lément 1 :
N 1 ( x ) =α 1+ β 1 x et N 2 ( x )=α 2 + β 2 x
 Conditions limites :
Nœuds 1 : N 1 ( 0 )=1 → α 1=1
N 2 ( 0 )=0→ α 2=0
Nœuds 2 : N 1 ( L2 )=0 → β =−2L

1
N2 ( L2 )=0 → β = 2L
2
2 2
D’où N 1 ( x ) =1− x et N 2 ( x )= x
L L
 É lément 2 :
N 2 ( x )=α 2+ β 2 x et N 3 ( x )=α 3 + β 3 x
Nœuds 2 : N2 ( L2 )=1 → α + L2 β =1
2 2
N3 ( L2 )=0→ α + L2 β =0
3 2
Nœuds 3 : N 2 ( L )=0 → α 2 + L β2 =0
N 3 ( L )=1 α 3 + L β 3=1
−2
D’où : β 2= , α 2=2
L
2 2
→ N 2 ( x )=2− x et N 3 ( x )=−1+ x
L L
- Maillage 2 :
16
 Les champs sont cinématiquement admissibles.
 On a un élément de barre constitué de trois nœuds ceci entraine qu’on a une fonction
de forme quadratique.
2
N 1 ( x ) =α 1+ β 1 x +C 1 x
2
N 2 ( x )=α 2+ β 2 x+C 2 x
2
N 3 ( x )=α 3+ β3 x+ C3 x
N 1 ( 0 )=1 → α 1=1
()
2
L L L
N1 =0 → 1+ β1 + C1=01 ¿
2 2 4
2
N 1 ( L )=0 →1+ L β 1+ L C1 =0
N 2 ( 0 )=0→ α 2=0
()
2
L L L
N2 =1 → β 2 + C 2=1 2 ¿
2 2 4
N 2 ( L )=0 → L β2 + L2 C1=0
N 3 ( 0 )=0→ α 3=0
()
2
L L L
N3 =0→ β3 + C 3=0 3 ¿
2 2 4
17
N 3 ( L )=1 → L β3 + L2 C3=0
−3
D’où : β 1=
L
2
C 1= 2
1¿
L
α 1=1
4
β 2=
L
−4
C 2= 2
2¿
L
α 2=0
−1
β 3=
L
2
C 3= 2
3¿
L
α 3=0
3 2 4 2 2 −1 2 2
→ N 1 ( x )=1− x + 2 x ; N 2 ( x )= x − 2 x et N 3 ( x )= x+ 2 x
L L L L L L
- Matrice de rigidité globale :

1
On a [ K e ]=∫ ¿ ¿ar le matériau est homogène linéaire et isotrope

V
L 1
¿ E ∫∫ ¿ ¿ ¿
0 S
[ K e ]=ES ∫ ¿ ¿
0
18
 Déterminons [ Be ]
On a { ε }=[D]{u }
Or { u } =[ N e ]{q }e
→ { ε }=[ D][ N e ]{q }e
{ ε }=[ Be ] {q }e
Avec : [ Be ]=[D][N e ]
- Maillage 1 :
 É lément 1 :
¿
¿
[ ]
∂
∂x
[N 1 N 2 ]
→¿
→¿
↦¿
 Élément 2
¿¿ ¿
Assemblage :
Matrice de rigidité globale :
[ ]
1 −1 0
2 ES
[ K ]= −1 2 −1
L
0 −1 1
Calcul des déplacements aux nœuds
19
 É lément 1 :
[ K e ] {qe }={F e }
1
{F e }1= {F }+ { F }=∫ ¿ ¿
V S
e e
V
L /2
¿ ∫∫ ρgds ∫ ¿ ¿
0
{1−22 x /xL/ L}dx

L/2
¿ ρgs ∫
0
2
x −x / L L/ 2
¿ ρgs { 2 }0
x /L
¿ ρgs {L/4
L/4 }
{F e }1=
ρgLs 1
4 1 {}
 É lément 2 :
1
{F e }2=∫ ¿ ¿
V
{ }
2
2 x−x / L L ¿ ρgs L/4
¿ ρgs { }
− x+ x 2 / L L/2 L/4
{F e }2={ Fe }1=
ρgLs 1
4 1 {}
- Force génératrice globale :
{ } {}
1 1
ρgLs ρgLs
{F e }= 1+1 = 2
4 4
1 1
D’où [ K ] { q }= { F }
- Application de la méthode directe :
20
Les déplacements aux nœuds sont uni-axiales.
Soit, u1 ,u 2 , u3 avec u1=0 car il y ' a encastrement
[ ]{ } { }
1 −1 0 0 1
2 ES ρgLs
u
−1 2 −1 2 = 2 , donc :
L 4
0 −1 1 u3 1
2 2
ρg L ρg L
2 u2−u3 = =¿ 2u 2−u3=
4E 4E
2 2
ρg L ρg L
−u2 +u3 = →−2u 2+2 u3=
8E 4E
2 2
ρg L ρg L
u3 = →u 3=
2E 2E
2 2 2
ρg L ρg L 3 ρg L
2 u2= + u2 =
4E 2E 8E
- Calcul des contraintes pour chaque élément :
 É lément 1 :
U 1 ( x )=[ N 1 N 2 ]
[] u1
u2
¿ N 1 ( x ) u1 + N 2 ( x ) u2
¿ N 2 ( x ) u2
2x 2
∗3 ρg L
L
¿
8E
3 ρgL
u1 ( x ) = x
4E
Par ailleurs, la déformation ε x =ε est donnée par la première dérivée car la fonction
de forme est linéaire.
∂ u1 ( x ) 3 ρgL
ε x= =
∂x 4E
21
3 ρgL
σ 1=Eε x =
4
 É lément 2 :
U 2 ( x )=[ N 2 N 3 ]
[]u2
u3
¿ N 2 ( x ) u2 + N 3 ( x ) u3
(¿ 2− 2Lx )∗3 ρg L + (−1+ 2Lx )∗ρg L

2 2
8E 2E
1 1
∗ρg L
2
ρg L 4
¿ + x
4E E
1
∗ρgL
∂ u2 ( x ) 4
ε x= =
∂x E
ρgL
D’où σ 2=E ε x =
4
- Matrice de rigidité globale du maillage 2 :
[ K e ]=∫ ¿ ¿
V
¿ [ H ] =E car ≤matériau est homogene et isotrope
[ K e ]=E ∫ ¿ ¿
V
Avec [ Be ] =[ D ] [ N e ] = [ d N1 d N 2 d N3
dx
;
dx
;
dx ]
¿
[ −3 4 x 4 8 x 1 4x
+ 2 ; − 2 ;− + 2
L L L L L L ]
¿
22
¿ L
[ K e ]=ES ¿ 0
[ ]
7 −8 1
3L
[ K e ]=ES 16 3L
¿
7
=
[
3 L ES 7 −8 1
3 L 16 ¿ 7 ]
3L 3L
{ F e }=∫ ¿ ¿
V
L
¿ ρgs ∫ ¿ ¿
0
¿ ρgs ¿
{ } { }
L /6 1/6
{ F e }=ρgs 2 L /3 =ρgsL 2/3
L /6 1/6
- Application de la méthode directe :
[ K ] = {q }= { F }
[ ]{ } { }
7 −8 1 u1 1/6
ES
−8 16 −8 u2 =ρgsL 2/3
3L
1 −8 7 u3 1/6
[
ES 16 −8 u2
3 L −8 7 u3
=ρgsL
2 /3
1/6 ]{ } { }
2 2
2 ρg L 3 ρg L
16 u2−8u 3= =¿ u2=
E 8E
2 2
ρg L ρg L
−8 u2 +7 u3= =¿ u3=
2E 2E
23
{}
0
{ U ( x ) }= [ N 1 N 2 N 3 ] u2
u3
{ U ( x ) }=N 2 ( x ) u 2+ N 3 (x)u3
(¿ 4L x− L4 x )∗3 ρg L + ( −1L x + L2 x )∗ρg L

2
2 2
2
2 2
8E 2E
ρg
{ U ( x ) }= 2
(2 Lx−x )
2E
dU ρg ρg
ε x= = ( 2 L−2 x )= ( L−x )
dx 2 E E
- Calcul des contraintes :
σ =E ε x =ρg(L−x)
24
PROJET 3 : MÉTHODE DES ÉLÉMENTS FINIS : É TUDE D’UN É LÉ MENT DE

POUTRE DE BERNOULLI
Y
U1 U2 X
ϴ1 ϴ2
V1 V2
Données :
Caractéristiques géométriques L, S, I sont des constantes
Caractéristiques mécaniques
E, f sont des constantes
I- Expression du champ axial u(x) en tout point de l’élément en fonction des

valeurs nodales u1 et u2.
U(x) = ax + b
CA →U(0) = U1
U(L) = U2
U(0)= a * 0 + b = U1
U(L)= aL + b = U2
U 2−U 1
→ b=U 1 → aL+ U1 = U2→ a=
L
U(x)= ( U 2−U 1
L )
x +U 1
U(x)= 1−( Lx ) U + Lx U
1 2
25
II- Degré du champ v
III- Expressions de V(x) et U
EIV’’(x)= - M(x)
Afin de pouvoir représenter exactement avec cet élément un champ de moment linéaire,
nous devons avoir un champ de déplacement V(x) de degré 3.
Soit V ( x )=α 1+ α 2 x+ α 3 x 2+ α 4 x 3
[]
α1
α2
V ( x )=[ 1 x x ² x ]
3
α3
α4
Conditions aux limites :
{
1
V y ( 0 )=V 1
1
'
V y ( 0 )=θ1
1
1
V ( L ) =V 2
y2
'
V ( L )=θ 2
y2
V 1=α 1
θ1=α 2
3
V 2=α 1 +α 2 L+α 3 L ²+ α 1 L
'
V 2=α 2 +2 α 3 L+3 α 4 L ²
[ ][ ][ ]
V1 1 0 0 0 α1
θ1 0 1 0 0 α2
= 2 3
V2 1 L L L α3
2
θ2 0 1 2L 3 L α4
[ ][ ][ ]
α1 1 0 0 0 V1
α2 = 0 1 0 0 θ1
2 3
α3 1 L L L V2
2
α4 0 1 2L 3L θ2
26
[ ]
1 0 0 0
0 1 0 0
A= 2
1 L L L3
0 1 2 L 3 L2
−1 1
A = ∗com [ A ] avec detA =L4
detA
[ ]
[ ] [ ][ ][ ]
1 0 0 0 0 0 0 1 0 0 1 0
L L ² L3 − 1 L 2 L3 1 L L 3
− 1 L L²
1 2 L 3 L² 0 2 L 3 L2 0 1 3 L² 0 1 2L
[ ][ ] [ ][ ]
0 0 0 1 0 0 1 0 0 1 0 0
− L L ² L3 1 L ² L3 −1 3
L L 1 L L²
2
1 L 2 L 3 L² 0 2L 3 L² 0 1 3L 0 1 2L
com A=
[ ] [ ] [ ][ ]
0 0 0 1 0 0 1 0 0 1 0 0
1 0 0 −0 0 0 0 1 0 −0 1 0
2
1 2 L 3 L² 0 2L 3L² 0 1 3L 0 1 2L
[ ] [ ] [ ][ ]
0 0 0 1 0 0 1 0 0 1 0 0
−1 0 0 0 0 0 0 0 0 0 0 0
3 3 3 3
L L² L 1 L² L 1 L² L 1 L² L
[ ]
4 2
L 0 3L −2 L
4 4
0 L −2 L L²
→ comA=
0 0 3 L ² −2 l
3
0 0 −L L²
[ ]
4
L 0 0 0
4
T 0 L 0 0
→❑ comA= 2 4 3
3 L −2 L 3 L ² −L
−2 L L² −2l L²
[ ]
1 0 0 0
0 1 0 0
−1 3 3 −1
→ A = 2 −2 2
L L L
−2 1 −2 1
3 2 3 2
L L L L
27
[] []
α1 V1
α2 θ1
=[ A ]
−1
α3 V2
α4 θ2
[ ][ ]
1 0 0 0
0 1 0 0 V1
3 3 −1 θ1
V ( x )=[ 1 x x ² x ] −2
3
2 2
L L L V2
−2 1 −2 1 θ2
3 2 3 2
L L L L
][ ]
U1
[
x x V1
1− 0 0 0 0
D’où | |
u
{ u } = [x ] =
v (x) 0
L
1+
2
3x 2x
− 3
3
x−2 x 2 +
x
3
L
0
2
3x 2x
− 3
3 2
−x x U 2
+
3
θ1
L2 L L2 L2 L L L V2
θ2
IV- Expression de [Be] et [Ke]
2
d d
{ ε }=[ B1e ] { q } on a : { ε }= U ( x )− y 2 V ( x)
dx dx
[]
U1
V1
−1
(
6 12 x
L
12 x 1
L L
6 12 x
) (
−2 6 x
Or { ε x } =[ L − 2 − 3 y− −4+ 2 y L − 2 − 3 y− L + L ² y ]
L L
θ1
U2 ) ( ) ( )
V2
θ2
Par ailleurs, { ε }=[ Be ] { q } . Il vient par identification que :

1
[]
U1
V1
−1 −6
( 12 x 12 x 1 −6 12 x
)(2 6 x
) (
[ B❑e ]=[ L L2 + L3 y 4− L2 y L L2 + L3 y L − L² y ] Uθ1
2
)( )
V2
θ2
28
IL manque une partie
V- Expression du vecteur des forces surfaciques { F se }
On a { F }=∫ ¿ ¿ ]T{ q s❑}dx

s
e
s
Notre étude étant ramenée à la fibre moyenne on a ds=1.dx

L
→ { F }=∫ ¿ ¿]T{ q 0❑ }dx
s
e
0
29
[ ]
x
1− 0
L
3 x 2 2 x3
0 1− 2 + 3
L L
2 x x3
2
L 0 x− ❑ + 2
{ F }=∫
s
e
0 x
L L
0
{0q}dx
L
3 x2 2 x3
0 − 3
L2 L
2 3
−x x
0 ❑ + 2
L L ❑
F1=0
L 2 3
3x 2x
F2=∫ (1− 2
+ 3 )dx
0 L L
3 4
x x L
¿ q [x − 2
+ 3]
L 2L 0
qL
F2=
2
L 2 3
2x x
F3=∫ (x− ❑ + 2 )dx
0 L L
2 3 4
x 2x x L
¿q[ − + 2]
2 3L 4L 0
qL ²
F3=
12
F4=0
L 2 3
3 x 2x
F5=∫ ( 2 − 3 )dx
0 L L
3 4
x x L
¿q[ 3
− 3]
L 2L 0
qL
F5=
2
L 2 3
−x x
F6=∫ ( + ) dx
0 L ❑ L2
30
3 4
−x x L
¿q[ ❑+ ]
3 L 4 L2 0
−qL ²
F6=
12
PROJET 4 : MÉTHODE DES ÉLÉMENTS FINIS : É TUDE DES É LÉ MENTS FINIS

TRIANGULAIRES
Y
O
2L x
2e
2L y
X
Figure : géométrie étudiée
I- É tude des éléments finis à trois nœuds
31
1) L’hypothèse pour le calcul bidimensionnel par élément finis est que nous
sommes en contraintes planes.
2) Proposition de maillage
1 4
II
2 3
3) Degré de liberté
Ly
1 4
Lx
II
2 3
 Chargement suivant ⃑
X → v 2=v 1=0
 Symétrie de chargement par rapport OY → u1=u 4=0
'
D où nos inconnus sont :u1 ,u3 , v 3 et v 4
Ainsi, la matrice de rigidité globale [K] est de dimension 4*4
- 4) Expression de [B e ]
[ Be ] =[ D ] [ N e ]
| |
∂
0
∂x
[ Be ] = 0
∂
∂ N1 0
∂ y 0 N2
∂
[ N2 0 N 3 0
0 N2 0 N3 ]
∂y ∂x
32
[ ]
∂ N1 ∂N2 ∂ N3
0 0 0
∂x ∂x ∂x
∂ N1 ∂ N2 ∂ N3
[ Be ] = 0
∂y
0
∂y
0
∂y
∂ N1 ∂ N1 ∂ N 2 ∂ N2 ∂ N3 ∂ N3
∂y ∂x ∂ y ∂x ∂y ∂x
On donne :
1 ∂Ni 1 ∂ Ni 1
N i ( x , y )=
2∆
( ai +b i x +Ci y ) ; ( x , y )= bi ; ( x , y )= c
∂x 2∆ ∂y 2∆ i
a i=x j y k −x k y J ; bi= y J − y K ; c i=x k −x j
[ ]
b 0 b2 0 b3 0
1 1
[ Be ] = 2 ∆ 0 C 1 0 C 2 0 C 3
C 1 b1 C 2 b2 C 3 b3
Comme [ Be ] est constant → l état de déformation est constant { ε }=[ B e ] { F e }

'
- Expression de [ K e ]
[ K e ]=∫ ¿ ¿
V
Nous sommes en contraintes planes donc :
| |
1 ν 0
E ν 1 0
[ H ]= 2
1−v 0 0 1−ν
2
D’où : [ K e ]=¿
-Calcul de K1 et K2
 É lément 1 :
On pose Ly= L2 et Lx=L1

L2
1(0 ; 0) 4 (0 ; L2)
I 33
L1
- Détermination des coefficients aréolaires
 b i= y J − y K
j=1+1=2
k = j+1=2+ 1=3 → k=4
b 1= y 2− y 4 =− y 4 =−l 2
 b 2 : j=2+1=3=¿ j=4
k =4 +1=5=¿ k =1
b 2= y 4 − y 1=¿ L2
 b 4=0
 C 1=x 4−x 2=−l 1
 C 2=x 1−x 4=0
 C 4 : j=4+1=5 ≡4 → j=1
k =1+1=2
→ C4 =x 2−x 1
C 4=l 1
Donc :
34
- Appliquons la méthode directe :
On obtient :
- Calcul de [ K 2 ¿
4 (0 ; L2)
L1
II
2 (L1 ; 0) 3 (L1 ; L2)

L2
On a :
[ ]
b2 0 b 3 0 b4 0
[B¿ ¿ e]= 0 C 2 0 C 3 0 C 4 ¿
C 2 b2 C 3 b3 C 4 b 4
[ K 2 ]=¿
 b 2 : j=2+1=3 avec k =3+1=4
b 2= y 4 − y 3=0
 b 3 : j=3+1=4 avec k =4 +1=5 ≡ [ 4 ] =2
b 3= y 4 − y 3=L2
b 4 : j=4+1=5≡ [ 4 ] (car 1 n existe pas)avec k =3+1=4

'

35
b 4= y 2 − y 3=−L2
 C 2 : j=3 avec k=4
C 2=x 4−x 3=− L1
 C 4 : j=2 et k=3
C 4=x 3−x 2=0
Donc :
[ ]
0 0 L2 0 −L2 0
[B¿ ¿ e]= 0 −L1 0 L1 0 0 ¿
−L1 0 L1 L2 0 −L2
Ainsi :
[ ]
0 0 −L1
| |[ ]
0 −L1 0 1 ν 0 0 0 L2 0 −L2 0
1 L2 0 L 1 E ν 1 0 1
[ K ¿¿ 2]= 0 −L1 0 L1 0 0 2e∆¿
2 ∆ 0 L1 L2 1−ν 2 1−ν 2 ∆
0 0 −L1 0 L1 L2 0 −L2
−L2 0 0 2
0 0 −L2
[ ][
0 0 −L1
]
0 −L1 0 0 −L1 ν L2 ν L1 −L2 0
eE L2 0 L1 0 −L1 L2 ν L1 −ν L2 0
¿
( )
2
L1 L2 (1−ν ) 0 L 1 L2 1−ν 1−ν L ( 1−ν ) 1−ν
0 −L2 ( )
− L2 0 0 −L1 ( 2 ) 0 L1
2
2
2 2
0 0 −L2
eE
[ K ¿¿ 2]= 2
¿¿
l 1 l 2 (1−ν )
36
- Appliquons la méthode directe
eE
[ K ¿¿ 2]= 2
¿¿
l 1 l 2 (1−ν )
On obtient la matrice symétrique :
eE
[ K ¿¿ 2]= 2
¿¿
l 1 l 2 (1−ν )
- Assemblage :
 Tableau de connectivité
Eléments Nœuds
I 1 2 4
II 2 3 4
[ ]
00 0 0 0 0
00 0 0 0 0
1 1
0 0 k 33 0 0 k 36
[ 1] 0 0 0 0
k =
0 0
00 0 0 0 0
0 0 k 163 0 0 k 166
[ ]
2 2 2 2
k 11 0 k 13 k 14 0 k 16
0 0 0 0 0 0
k2 k2 k2 k 236
[ k 2 ]= 231 00 332 234 0
0
k 41 k 43 k 44 k 246
0 0 0 0 0 0
k 61 0 k 63 k 64 0
2 2 2 2
k 66
37
[ ][
2 2 2 2
]
k 11 0 k 13 k 14 0 k 16
0 0 0 0 0 0 k 211 k 213 k 214 k 216
k 31 0 k 33+ k 33 k 234
2 2 1 2
0 k 36+ k 36
1
k 231 k 233 +k 133 k 234 k 236 +k 136
[ G] 2 0 2 2
k =
0
=
k 41 k 43 k 44 k 246 k 241 k 243 k 244 k 246
0 0 0 0 0 0 k 261 k 263 +k 163 k 264 k 266 +k 166
2 0 2 1 2 0 k 2 + k1
k 61 k 63+ k 63 k 64 66 66
Donc :
[ k G ]=¿
5) Force surfacique :
1
{ F e }=∫ ¿ ¿
S
[ ]
S
F xi
{ F Se }= S
F yi
{ F xi }=∫ ¿ ¿
S
{ F }=∫ ¿ ¿
S
yi
S
- Elément 1 :
L
1(0 ; 0) 4 (0 ; L2)
2
I
L
1
2 (L1 ; 0)
S ' '
{F xi }1={0 0 0 0 0 0 }Car aucune charge n estappliquée à l élément .
- Elément 2 :
4 (0 ; L2)
38
École Polytechnique de Masuku _ Ingénieur
L1 3 GC _ 2023-2024
II
2 (L1 ; 0) 3 (L1 ; L2)

1
{F xi }2=∫ ¿ ¿
S
L2
¿ 2 e q 0∫ N ( x , y ) dl
0
1∗L2
¿ 2 e q0( )
2
S
{F xi }2=e q 0 L2
S T
{F e }2 ={e q0 L2 0 e q0 L2 0 0 0 }
D’où :
S T
{F } ={e q0 L2 0 e q0 L2 0 0 0
6) Expression du déplacement
On sait que :
{}
u1
v1
{u}=
[ N1 0 N2 0 N3 0
0 N 1 0 N 2 0 N3 ] u2
v2
u3
v3
- Elément 1 :
{}
0
0
{u }1= { }[
u(x , y) N 1 0 N 2 0 N 4 0
v (x , y) 0 N 1 0 N 2 0 N 4 ] u2
0
0
v4
1
u ( x , y )=N 2 u2= (a + b x+C 2 y)u2
2∆ 2 2
39
1
v ( x , y )=N 4 v 4= (a +b x +C 4 y )v 4
2∆ 4 4
a 2=0 et a4 =0
b 2=L2 b 4=0
C 2=0 C 4=L1
1
u ( x , y )= (L x)u2
2∆ 2
1
v ( x , y )= (L y )v 4
2∆ 1
L1 L2
Avec ∆=
2
x
u ( x , y )= u
L1 2
y
⇨ v ( x , y )= v
L2 4
Or :
q0 −2 νq0
u2 = et v 4 =
E E
q0 q0
u ( x , y )= x= x
EL1 E
−2 νq0 −0.3 q0
v ( x , y )= y= y
EL2 E
- Elément 2 :
40
{}
u2
0
{u }4 =
[ N2 0 N3 0 N4 0
0 N2 0 N3 0 N4 ] u3
v3
0
v4
u ( x , y )=N 2 u2 + N 3 u3
v ( x , y )=N 3 v 3 + N 4 v 4
u ( x , y )=(N ¿ ¿ 2+ N 3 )u 2 ¿
v ( x , y )=(N 3 + N 4 ) v 3
1
 N 3= (a + b x+C 3 y)u2
2∆ 3 3
a 3=−L1 L2
b 3=L2
C 2=L1
1
¿> N 3= (−L1 L2 + L2 x + L1 y )
2∆
1 1
N 3=−1+ x+ y
L1 L2
1
 N 2= (a + b x +C 3 y )u2
2∆ 2 2
a 2=L1 L2
b 3=0
C 2=−L1
1
¿> N 2= (L L −L y )
2∆ 1 2 1
1
¿> N 2=1− y
L2
41
1
 N4= (a + b x +c 4 y )
2∆ 4 4
a 4=L1 L2
b 4=−L2
C 4=0
1
¿> N 4=1− x
L1
x
u ( x , y )= u
L1 2
1
v ( x , y )= y v3
L2
q0
u ( x , y )= x
E
−y
∗2 ν q 0
2 −0.3 ν
( )
v x, y = = y
E E
q0
u ( x , y )= x
E
−y
∗2 ν q 0
2
v ( x , y )=
E
D’ou:
{ }
q0
x
E
[ ]
{ U }= u( x ) =¿
v ( x) −y
1
∗2 ν q 0
2
E
- Calcul des déformations et contrainte

 Elément 1 :
42
{ ε }=[ D ] {u }
[ ]
∂
0
∂x
¿ 0
{
∂ u( x , y )
∂ y v ( x , y) }
∂ ∂
∂y ∂x
[ ε ]=¿
{ ε }=
q0 1
E { }
−0.3
0
{ σ }=[ H ] { ϵ }
| |{
1 ν 0
}
1
E ν 1 0 q0
¿ −0.3
2
1−ν 0 0 1−ν E
0
2
|
q 0 1 0.3
|{ }
0 1
¿ 0.3 1 0 −0.3
0.91
0 0 0.35 0
¿
q 0 0.91
0.91
0
0 { }
{}
1
{ σ }=q 0 0
0
- Elément 2 :
{ ε }=
q0 1
E { }
−0.3
0
|
q 0 1 0.3
|{ }
0 1
{ σ }= 0.3 1 0 −0.3
0.91
0 0 0.35 0
43
q 0 0.91
{ } {}
1
¿> { σ }= 0 = q0 0
0.91
0 0
Partie 2
II- É tude des éléments finis triangulaires à 6 nœuds
v3
u 3
3
v5
Y v6 u 5
u 6 v2
6 5
X
v4
v1 u 2
2
u 4
4
1 u 1
Figure : élément fini triangulaire à 6 nœuds
Y
L1
L2
X
M
L3
2
1
Figure : coordonnées aréolaires
44

Projet Mef

Transféré par

Droits d'auteur :

Formats disponibles

Projet Mef

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Projet Mef

Transféré par

Droits d'auteur :

Formats disponibles

UNIVERSITE DES SCIENCES ET TECHNIQUES DE MASUKU

ECOLE POLYTECHNIQUE DE TECHNIQUES DE MASUKU

Département de Génie Civil

TABLE DES MATIÈ RES

PROJET 1 : RÉ SOLUTION D’UN PROBLÈ ME D’É LASTICITÉ __________________________________2

PROJET 1 : RÉSOLUTION D’UN PROBLÈME D’ÉLASTICITÉ

Exemple d’utilisation de la méthode de Ritz : Étude d’une poutre

I- Résolution par la méthode de la Théorie des poutres (RDM)

Déterminons l’expression de la flèche v ( x ) et du moment fléchissant M(x).

- Calcul des réactions aux appuis

D’après le principe fondamental de la statique :

- Calcul des efforts internes

Considérons la partie gauche de la coupure : 0 ≤ x ≤ L

Équation de la déformée d’après TIMOSCHENKO

II- Application du principe des travaux virtuels (PTV)

On cherche maintenant à obtenir une solution approchée de ce problème par la méthode

Or M ( x ) =−EIV ' '( x )

Conditions aux limites :

∫ qδα (x 2¿−Lx )dx =qδα ∫ (x 2 ¿−Lx)dx ¿ ¿

Divisons la poutre en 4 éléments de même longueur

On peut appliquer les conditions aux limites :

v1 ( L4 )=v ( L4 ) → 161 α L + 14 β L= 161 α L + 14 β L+ γ 1 ¿

v'1 ( L4 )=v ' ( L4 )=¿ 12 α L+ β = 12 α L+ β

On obtient le système suivant :

Application du PTV : δ W e =−δ W i

¿−EI ∫ (2 α 1)(2 δ α 1 ¿ )dx ¿

 δ W i 2=−EI ∫ (2 α 2)(2 δ α 2 ¿ )dx ¿

III- Solution trigonométrique

PROJET 2 : FORMULATION DE LA MEF EN DÉPLACEMENT

Le système étant linéaire et la barre isotrope homogène on a :

Exemple 2 : Maillage de la barre

On se propose de considérer 2 maillages différents :

- Un maillage avec 2 éléments unidimensionnels à 2 nœuds de longueur

 Les champs sont cinématiquement admissibles.

Nœuds 2 : N 1 ( L2 )=0 → β =−2L

 Les champs sont cinématiquement admissibles.

- Matrice de rigidité globale :

On a [ K e ]=∫ ¿ ¿ar le matériau est homogène linéaire et isotrope

→ { ε }=[ D][ N e ]{q }e

Matrice de rigidité globale :

Calcul des déplacements aux nœuds

{1−22 x /xL/ L}dx

- Application de la méthode directe :

Les déplacements aux nœuds sont uni-axiales.

Soit, u1 ,u 2 , u3 avec u1=0 car il y ' a encastrement

- Calcul des contraintes pour chaque élément :

(¿ 2− 2Lx )∗3 ρg L + (−1+ 2Lx )∗ρg L

- Matrice de rigidité globale du maillage 2 :

¿ [ H ] =E car ≤matériau est homogene et isotrope

- Application de la méthode directe :

(¿ 4L x− L4 x )∗3 ρg L + ( −1L x + L2 x )∗ρg L

- Calcul des contraintes :

PROJET 3 : MÉTHODE DES ÉLÉMENTS FINIS : É TUDE D’UN É LÉ MENT DE

I- Expression du champ axial u(x) en tout point de l’élément en fonction des

II- Degré du champ v

III- Expressions de V(x) et U

Conditions aux limites :

IV- Expression de [Be] et [Ke]