TD de méthode des éléments finis

Ecole hassania des travaux public Annee academique : 2015/2016

Departement ponts chaussees et tranports Le 29/10/2015
Calcul des structures en element finis

Devoir libre en Calcul des structures en element finis

Prepare par : AL ABDALI Abdelhamid - EL FATHI Izdin
Classe : 2GC1

EXERCICE 1 :

Soit le treillis de la figure suivante:

Figure 1: Shéma du treillis

L’objectif de cet exercice est de :

• Calculer les déplacements aux noeuds.

• Calculer les réactions aux noeuds 1 et 3.

• Calculer les forces intérieurs dans les barres du treillis .

• Et enfin,calculer tous éléments avec la Résistance des matériaux .

1
1. Méthode des déplacements :

Elément 1-2:
La matrice de rigidité correspond à cet élément est :
   
ES 1 −1 u1
ke1 = , sur
l1 −1 1 u2
Elément 1-3 :
La matrice de rigidité correspond à cet élément est :
   
ES 1 −1 v1
ke2 = , sur
l2 −1 1 v3
Elément 2-3 :
La matrice de rigidité correspond à cet élément est :
 
  u2
ES [A] −[A]  v2 
ke3 = , sur  
l3 −[A] [A] u3 
v3
Avec :  
cos2 φ cos φ sin φ 3π
[A] = et φ =
cos φ sin φ sin2 φ 4
Puisque les barres ont une section circulaire :

πd2 30.109 .π.0.0362

ES = E = = 2.1375.108 N
4 4
Tout calculs fait ,on trouve :
   
0.250 −0.250 u1
ke1 = ES , sur
−0.250 0.250 u2
   
0.200 −0.200 v1
ke2 = ES , sur
−0.200 0.200 v3
   
0.061 −0.076 −0.061 0.076 u2
 −0.076 0.095 0.076 −0.095   v2 
ke3 = ES 
 −0.061 0.076
 , sur  
0.061 −0.076  u 3 
0.076 −0.095 −0.076 0.095 v3
En tenant compte les conditions aux limites, les vecteurs déplacement nodal et force nodal sont
respectivement :

       
u1 0 Fx1 X1
 v1   0  Fy1   Y1 
       
→ u2  u2  → Fx2   Fx2 
δ= 
 v2  =  v2  et F = Fy2  = −Fy2 
      
       
u 3   0  Fx3   X3 
v3 v3 Fy3 0

2
 → →
En assemblant les trois matrices ke1 , ke2 et ke3 ,puis écrivons le système total [K] δ =F :

    
0.250 0 −0.250 0 0 0 u1 Fx1

 0 0.200 0 0 0 −0.200   v1  Fy1 
   
 −0.250 0 0.250 + 0.061 0.076 0.061 0.076  u2  Fx2 
ES    =  

 0 0 0.076 0.095 0.076 0.095   v2  Fy2 
   
 0 0 0.061 0.076 0.061 0.076  u3  Fx3 
0 −0.200 0.076 0.095 0.076 0.20 + 0.095 v3 Fy3
ou:
      
0.250 0 −0.250 0 0 0 0 X1 X1

 0 0.200 0 0 0 −0.200 
  0   Y1   Y1 
     
 −0.250 0 0.311 0.076 0.061 0.076  u2   Fx2   180 
     
ES   v2  = −Fy2  = −240


 0 0 0.076 0.095 0.076 0.095       
 0 0 0.061 0.076 0.061 0.076   0   X3   X3 
0 −0.200 0.076 0.095 0.076 0.295 v3 0 0

Le système réduit s’écrit :

    
0.311 0.076 0.076 u2 180
ES  0.076 0.095 0.095   v2  = −240
0.076 0.095 0.295 v3 0
Les réactions inconnues sont donnés par le système :
    
X1 −0.250 0 0 u2
 Y1  = ES  0 0 −0.200   v2 
X3 0.061 0.076 0.076 v3

Donc :
−50×103
 
 u2 = ES
= −0.0234 cm  X1 = 12.5 kN
−3807.681×103
v = = −1.78 cm et Y1 = 240.5 kN
 2 ES
−1247.507×103
= −0.58 cm X3 = −387.24 kN

v3 = ES

L’éffort inetrne de chaque barre est donné par la formule :

ES
Nij = [(uj − ui ) cos φ + (vj − vi ) sin φ]
le
Soit :

ES ES
• Pour la barre 1-2: N12 = l1
(u2 − u1 ) = u
l1 2
= −12.5 kN

ES ES
• Pour la barre 1-3: N13 = l2
(v3 − v1 ) = v
l2 3
= −240.5 kN

ES
• Pour la barre 2-3: N23 = l3
[(u3 − u2 ) cos φ + (v3 − v2 ) sin φ] = 307.350kN
,avec φ = 128◦ .65

3
2. Calcul avec la rdm :

Figure 2: Equilibre des noeuds

Equilibre du noeud 2 :
−→ −→ →
−
N21 + N23 + 180→
−
ex − 240→
−
ey = 0
Soit après la projection on trouve :

−N21 − N23 cos θ + 180 = 0
N23 sin θ − 240 = 0
Ce qui donne ,après le calcul:
240

N23 = sin θ
= 307.350kN
N21 = 12kN
Equilibre du noeud 2 :
−→ −→ →
−
N12 + N13 + X1 →
−
ex + Y1 →
−
ey = 0
Soit après la projection on trouve :

X1 + N12 = 0
Y1 + N13 = 0
Ce qui donne ,après le calcul:

X1 = 12kN
Y1 = 240kN
Equilibre du noeud 3 :
−→ −→ →
−
N32 + N31 + X3 →
−
ex = 0
Soit après la projection on trouve :

X3 + N32 cos θ = 0
N31 + N32 sin θ = 0

4
Ce qui donne ,après le calcul:

X3 = 192kN
Y1 = 240kN
En résumé :  
 X1 = 12 kN  N12 = −12 kN
Y1 = 240 kN et N13 = −240 kN
X3 = 192 kN N23 = 307.350 kN
 

Calculs des déplacements u2 , v2 et v3 :

On d’après la formule de Bresse en cas de la traction/compression:
Z x1
−−−→ −−−→ −−−→ −−−→ N (x) →
−
u(x1 ) = u(x0 ) + ω(x0 ) ∧ G0 G1 + ex dx
x0 ES

Après la projection ,on trouve :

NL
(v(x1 ) − v(x0 )) sin φ + (u(x1 ) − u(x0 )) cos φ =
ES
Tout calculs fait , on trouve :

 u2 = −0.0234 cm
v2 = −1.78 cm
v3 = −0.58 cm


3. Vérification avec le logiciel Robot :

On lance le logiciel Robot ,on choisit comme envirenement de travail : Les treillis plans,puis
on introduit notre treillis qui est composé de trois barre d’acier de diamètre d = 36 mm et
de module de Young E = 210 GP a,on trouve des résultats très proche de ce qu’on calculé
précedement soit avec la Mef ou la Rdm ,les figures suivantes illustrent les résultats obtenus
avec le même logiciel:

Figure 3: Diagramme de l’effort normal

5
 Ainsi le déplacement aux nœuds,les efforts normaux et les réactions inconnus:

Réactions inconnus :

N oeud Xi (kN ) Yi (kN )

1 12 240
3 −192 0

Déplacement inconnus :
N oeud ui vi
1 0 0
2 −0.02 −1.8
3 0 −0.6

Efforts normaux :

Barre Nij
1−→ 2 12
2−→ 3 −307.35
1−→ 3 240

Conclusion:
D’après ce qui précède,on trouve que les résultas sont identiques ,cela due à la modélisation
simple par la méthode des éléments finis .En effet ,le champs déplacement approché par la Mef
est identique au champs déplacement réel(Résultats confirmé aussi par le logiciel Robot ).
NB : vous pouvez vérifier les résultats par un autre logiciel comme le SAP2000 .

EXERCICE 2 :

Soit la poutre continue de la figure suivante:

Figure 4: Schéma de la poutre continue

6
1. Modélisation par la Mef :Modèle à deux élément

Elément 1-2 :
     
12 6l1 −12 6l1 0.188 0.375 −0.188 0.375 v1
ES  6l1 4l12 −6l1 2l12   0.375 1 −0.375 0.5  
θ1 
 
ke1 = 3   = ES  , sur
l1  −12 −6l1 12 −6l1   −0.188 −0.375 0.188 −0.375  v2 
6l1 2l12 −6l1 4l12 0.375 0.5 −0.375 1 θ2
Elément 2-3 :
     
12 6l2 −12 6l2 0.444 0.667 −0.444 0.667 v2
ES  6l2 4l22 −6l2 2l22   0.667 1.333 −0.667 0.667 
θ2 
 
ke2 = 3   = ES  , sur
l2  −12 −6l2 12 −6l2   −0.444 −0.667 0.444 −0.667  v3 
6l2 2l22 −6l2 4l22 0.667 0.667 −0.667 1.333 θ3
Le vecteur déplacement est donné par(En utilisant les conditions aux limites) :
→
T T



δ= v1 θ1 v2 θ2 v3 θ3 = 0 0 0 θ2 0 θ3

Le vecteur force nodal connu et le vecteur force inconnu sont donnés respactivement par :
−
→  2 2

FC = T − ql21 − ql121 − ql21 − F2 ql121 − F8l2 − F2 F8l2 = T −120 80 −200 20 −80 60

et
−
→ T


FI = Y1 Z1 Y2 0 Y3 0
→ −
→ − →
En assemblant les deux matrices ke1 et ke2 ,puis écrivons le système total [K] δ = FC + FI :
      
0.188 0.375 −0.188 0.375 0 0 0 −120 Y1
 0.375
 1 −0.375 0.5 0 0   0   −80  Z1 
      
 −0.188 −0.375 0.632 0.292 −0.444 0.667    0  = −200 +  Y2 
     
ES  0.375
 0.5 0.292 2.333 −0.667 0.667  θ2   20 
   
 0
 
 0 0 −0.444 −0.667 0.444 −0.667   0   −80   Y3 
0 0 0.667 0.667 −0.667 1.333 θ3 60 0
Le système réduit s’écrit :
    
2.333 0.667 θ2 20
ES =
0.667 1.333 θ3 60
Les réactions inconnues et les moments sont donnés par le système :
   
Y1 − 120 0.375 0  
 Z1 − 80 
 = ES  0.5 0
  θ2
 
Y2 − 200  0.292 0.667  θ3
Y3 − 80 −0.667 −0.667
Donc :

 Y1 = 118.120 kN
−5.013×103
 
θ2 = Z1 = 77.500 kN.m

ES et
47.520×103 Y2 = 230.231 kN
θ3 = ES


Y3 = 51.648 kN


Diagramme du moment fléshissant :

7
 Un calcul de Rdm rapide donne l’expression du moment fléshissant :

 −77.5 − 118.31x + 30.02x2 , si 0 ≤ x ≤ 4 P outre1
M (x) = −85.10 − 108.76x , si 0 ≤ x ≤ 1.5 P outre2
77.54 + 51.74x , si 1.5 ≤ x ≤ 3 P outre2


L’effort tranchant est donné par la formule V (x) = − ∂M∂x(x) ,soit :


 118.31 − 60.08x , si 0 ≤ x ≤ 4 P outre1
V (x) = 108.76 , si 0 ≤ x ≤ 1.5 P outre2
−51.74 , si 1.5 ≤ x ≤ 4 P outre2


2. Modélisation par le logiciel Robot :

En choisit comme envirennement de travail dans le logiciel Robot ,étude d’une portique
plane ,on trouve les résultats suivants:

Figure 5: Diagramme du moment fléshissant obtenu par Robot

Figure 6: Diagramme de l’effort tranchant obtenu par Robot

Réaction : valeur
N oeud Yi ( kN ) Mi ( kN.m)
1 118.29 −77.61
2 230.52 0
3 51.75 0

8
 Figure 7: Diagramme de la déformé obtenu par Robot

Remarque : ce document est préparé par AL ABDALI Abdelhamid - EL FATHI Izdin , Ecole Hassania des traveux public .