Exii 1
Exii 1
Exii 1
h h
{U } vecteur des dplacements nodaux (ui , vi ) {FD } vecteur force gnralis associ aux efforts donns {FI } vecteur force gnralis associ aux efforts inconnus
Calcul de la matrice raideur Exprimer la matrice raideur de chaque lment sur ses variables nodales. En dduire l'expression de la matrice raideur assemble complte. Extraire la matrice raideur rduite. Rsolution statique - Efforts aux appuis Dterminer la dforme statique, et reprsenter l'allure de la dforme. Calculer les efforts aux appuis, et vrifier l'quilibre global de la structure. Post traitement Calculer les efforts (contraintes) sur chaque lment, puis vrifier l'quilibre du nud qui est charg. Isoler une des barres 45 de la structure, et calculer les efforts extrieurs sur cet lment. Retrouver les rsultats prcdents. Utiliser la symtrie Prciser le nouveau maillage en tenant compte de la symtrie. Calculer la matrice raideur rduite et retrouver la solution en dplacement.
h
(1)
U T =< u1 , v1 , u 2 , v 2 , u 3 , v 3 , u 4 , v 4 >
Bilan quations inconnues Conditions aux limites (u i , v i ) = (0,0) pour i=2, 3, 4 Il reste deux dplacements nodaux inconnus : Le bilan des efforts extrieurs nous donne 1 effort donn au nud 1 6 efforts de liaison aux nuds 2, 3, 4
yo xo
Effort donn
U I T =< u1 v1 >
T
2 dplacements inconnus
6 efforts inconnus Vous devez savoir crire le travail virtuel des efforts, pour obtenir ces vecteurs de type force gnralise. Bilan des quations : quations de la statique : Nous pouvons rsoudre.
8 quations / 8 inconnues
En pratique pour limiter les calculs la main on peut se limiter la construction du systme rduit aux dplacements inconnus.
[K red ]{U I } = { D }
Calcul de la matrice raideur 2 C ES C S Rappel : [K e ] = 2 C e C S
Application la structure tudie lment 1 1 2 2
avec { I } = U
u1 0 et { D } = F v1
La dcomposition du systme par blocs est prsente dans le cours, ainsi que le calcul des efforts inconnus utilisant les 6 quations restantes.
2 C C S 2 C C S
C S 2 S C S 2 S
C S 2 S sur (ui , vi , u j , v j ) C S 2 S
ES h 2
1 2 1 2 1 2 1 2
1 2 1 2 1 2 1 2
h 2
(1)
Matrice raideur : k1 =
3 = 4 xo
1
1 2 1 [K1 ] = k1 21 2 1 2
1 2 1 2 1 2 1 2
sur (u1 , v1 , u2 , v2 )
lment 2 1 3
3 (2) 1
Matrice raideur : k 2 =
=
xo
ES h 1 [K 2 ] = k 2 1 1 sur (v1 , v3 ) 1
lment 3 1 4
Matrice raideur : k 3
ES h 2
h 2
(3)
1 2 1 2 1 2 1 2
1 2 1 2 1 2 1 2
12 0 1 2 12 [K ] = ES 202 h Vous comprenez maintenant 0 pourquoi, la main pour viter 1 des critures fastidieuses, il est 2 2 pratique de ne calculer que la matrice rduite 212
0
1 2
+1
1
2 2 1 2 2
1 2 2 1 2 2 1 2 2 1 2 2
1 2 2 1 2 2 1 2 2 1 2 2
0 1
1 2 2 1 2 2
0 0 0 0
0 0 0 0
0 0 0 0 0 0 0 0
0 1 0 0 0 1 0 0
1 2 2 1 2 2
1 2 2 1 2 2
0 0 0 0
1 2 2 1 2 2
0 0 0 0
1 2 2 1 2 2
[K red ]{U I } = { D }
0 avec { D } = F
0 2 Fh {X I } = 1 + 2 ES
La dforme respecte la symtrie du problme. Nous obtenons la solution analytique du problme statique partir dune formulation matricielle base sur la notion dapproximation. En effet pour un treillis charg aux nuds, la solution en dplacement est linaire par morceaux.
212 X 2 1 Y 2 2 2 X 3 ES 0 = Y3 h 0 1 X 4 212 Y4 2 2
Il est alors possible de vrifier les quations dquilibre de la structure Pour les rsultantes :
X 2 + X 3 + X 4 = 0 Y2 + Y3 + Y4 F = 0
La solution respecte la symtrie, lquation dquilibre des moments est donc vrifie.
Post traitement
Le calcul des efforts dans les barres utilise la loi de comportement intgre : N = ES u,x
ES ES Ne = (u j ui ) = cos e e
Application la structure tudie lment 1 1 2 2
u j ui sin v j vi
h 2
(1)
3 = 4 xo
1
N1 =
ES h 2
2 2
F
2 0 2 v1
N1 =
2+ 2
lment 2 1 3 lment 3 1 4
ES v1 h
N2 =
F 2 1+ 2
F 2+ 2
N1
N2
N3
Suivant x : Suivant y :
N1 2 N1 2
N1 2
=0
+ N2 F = 0
Du point de vue numrique, un calcul identique nous donne une information sur lerreur numrique du modle. Ici lapproximation reprsente le champ exact il ny a donc pas derreur de discrtisation.
{Fie } = [K e ]{U e }
Y4
Y13 1
1 1 1 2 2 2 X4 1 1 1 4 2 2 k3 2 (3) 1 1 1 2 2 2 1 1 1 X 13 2 2 2 X 13 1 1 Y 1 1 ES F 13 Do v1 = = X 4 2h 2 1 2(1 + 2 ) 1 Y4 1 1
1 2 0 X 1 13 v Y 2 1 = 13 1 0 X 4 2 0 Y 1 4 2
Nous retrouvons les efforts aux appuis (nud 4). Un calcul identique pour les lments 1 et 2, nous permettrait de dterminer tous les efforts agissant sur le nud 1, il serait alors possible de vrifier lquilibre du nud interne.
X 1e = 0 e Y1e F = 0 e
Le fichier MEFLAB propos sur le site vous permettra de retrouver tous ces rsultats partir dun calcul numrique. Pour utiliser ces fichiers de donnes il faut pralablement avoir tlcharg le dossier compress de lapplication pour linstaller dans votre rpertoire de travail.
Utilisation de la symtrie
Attention aux conditions de symtrie Llment 2 et la charge appartiennent au plan de symtrie Ce modle ne possde quun degr de libert v1 La matrice raideur rduite : K red = Do la solution cherche v1 =
S/2
yo h xo
1
ES 1 1 + 2h 2
F/2
Fh 1 + 2 ES
Lorsquune structure possde des symtries, il est toujours avantageux de les utiliser dans les modles lments finis.