Methode des elements finis :

treillis plans a nuds articules

Yves Debard
Institut Universitaire de Technologie du Mans
Departement Genie Mecanique et Productique

http://iut.univ-lemans.fr/ydlogi/index.html

24 mars 2006 29 mars 2011

Table des matieres
Introduction 3

1 Matrices elementaires 3

2 Exemple 1 : treillis soumis a une force nodale 6

2.1 Enonce . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.2 Partition des degres de liberte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.3 Etude elementaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.4 Assemblage et calcul des deplacements inconnus . . . . . . . . . . . . . . . . . . . . . 7
2.5 Efforts normaux dans les elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.6 Actions de liaison . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

3 Exemple 2 : treillis soumis a une force nodale 8

3.1 Enonce . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3.2 Partition des degres de liberte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3.3 Etude elementaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3.4 Assemblage et calcul des deplacements inconnus . . . . . . . . . . . . . . . . . . . . . 9
3.5 Efforts normaux dans les elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.6 Actions de liaison . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

4 Exemple 3 : treillis soumis a une variation de temperature 10

4.1 Enonce . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
4.2 Partition des degres de liberte . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
4.3 Etude elementaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
4.4 Premier cas de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4.4.1 Assemblage et calcul des deplacements inconnus . . . . . . . . . . . . . . . . . 12
4.4.2 Efforts normaux dans les elements . . . . . . . . . . . . . . . . . . . . . . . . . 12
4.4.3 Actions de liaison . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4.5 Deuxieme cas de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
4.5.1 Assemblage et calcul des deplacements inconnus . . . . . . . . . . . . . . . . . 13
4.5.2 Efforts normaux dans les elements . . . . . . . . . . . . . . . . . . . . . . . . . 13
4.5.3 Actions de liaison . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

5 Programmes Maple 14
5.1 tre mat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
5.2 exemple 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
5.3 exemple 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
5.4 exemple 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

References 16
Treillis plans a nuds articules 3

Introduction

Un treillis est un ensemble de poutres droites (elements) reliees entre elles par des rotules (nuds).
Les liaisons exterieures sont des rotules et des appuis simples. Les charges sont des forces portees par
les rotules, des gradients thermiques et des deplacements dappui. La force interieure dans une
section droite se reduit a leffort normal.

Le treillis est plan si :

Le plan {O; x, y} est un plan de symetrie pour toutes les sections droites.
Les forces appliquees sont situees dans le plan {O; x, y}.

On suppose que les deplacements sont petits.

1 Matrices elementaires
Soit (i j) un element de treillis plan de section droite constante (figure 1).

Figure 1 Element i j

L est la longueur de lelement et A laire de sa section droite.

(xi , yi ) et (xj , yj ) sont les coordonnees des nuds de lelement.

Le vecteur unitaire ~n porte par laxe de la poutre est defini par :

nx 1 xj xi cos
= = , L2 = (xj xi )2 + (yj yi )2 (1.1)
ny L yj yi sin

ou est langle que fait ~n avec laxe x.

E et sont respectivement le module de Young et le coefficient de dilatation du materiau.

Lelement est soumis a un effort normal N (positif : traction, negatif : compression) et a une variation
de temperature T constante.
4 Methode des elements finis

(ui , vi) et (uj , vj ) sont les deplacements nodaux (figure 2).

Figure 2 Deplacements elementaires

Les efforts aux extremites de lelement sont (figure 3) :

N ~n en i , N ~n en j (1.2)

Figure 3 Efforts elementaires

Differentions la relation :
L2 = (xj xi )2 + (yj yi )2 (1.3)
Il vient :
2 L dL = 2 (xj xi ) (dxj dxi ) + 2 (yj yi ) (dyj dyi ) (1.4)
dou lexpression de lallongement unitaire suivant ~n :

dL 1 (xj xi ) (yj yi )
n = = (dxj dxi ) + (dyj dyi ) (1.5)
L L L L
soit :
1
n = ( nx (uj ui ) + ny (vj vi ) ) (1.6)
L
Cet allongement unitaire est du a leffort normal (loi de Hooke) et a la variation de temperature :
N
n = + T (1.7)
EA
Leffort normal secrit en fonction des deplacements nodaux :

N = EA (n T )
EA (1.8)
= ( nx (uj ui ) + ny (vj vi ) ) EA T
L
soit :
ui

EA vi
N= nx ny nx ny EA T (1.9)
L
uj

vj
Treillis plans a nuds articules 5

On en deduit :
{fnod } = [ k ] {u} {fth } (1.10a)
avec :

nx

ui

ny vi
{fnod } = N , {u} = (1.10b)

n u
x j

ny vj
2

nx
nx nx ny n2x nx ny
EA ny EA
n2y nx ny n2y
[k] = nx ny nx ny = 2 (1.10c)
L
nx
L n x nx n y
2
ny sym. ny


nx


ny
{fth } = EA T (1.10d)

n
x
ny
{fnod } est le vecteur force nodal (N).

{u} est le vecteur deplacement elementaire (m).

[ k ] est la matrice de rigidite (N/m).

{fth } est le vecteur force equivalent au gradient thermique (N).

Remarque 1 : la matrice de rigidite peut se mettre sous la forme :

[k] [k] EA n2x nx ny
[k] = avec [k] = (1.11)
[k] [k] L nx ny n2y

Remarque 2 : la matrice de rigidite est egale a :

nx 0
ny 0 EA 1 1 nx ny 0 0
[k] = (1.12)
0 nx L 1 1 0 0 nx ny
0 ny
ou :

nx ny 0 0 1 0 1 0 nx ny 0 0
ny nx 0 0 EA 0 0 0 0 ny nx 0 0
[k] =
0
(1.13)
0 nx ny L 1 0 1 0 0 0 nx ny
0 0 ny nx 0 0 0 0 0 0 ny nx

Remarque 3 : lenergie de deformation est egale a (a un coefficient pres independant des deplacements
et de leurs derivees) :

1
Edef = EA 2n L EA n T L
2 (1.14)
1
= {u}T [ k ] {u} {u}T {fth }
2

Le travail des forces exterieures se reduit au travail des forces nodales :

Wext = {u}T {fnod } (1.15)

6 Methode des elements finis

Lenergie potentielle est egale a :

Epot = Edef Wext (1.16)

La matrice de rigidite est la matrice hessienne (ou matrice de Hess) de lenergie de deformation par
rapport aux deplacements nodaux (programme tre mat) :

2 Edef
kij = (= kji ) (1.17)
ui uj

Le vecteur des efforts aux nuds est le gradient de lenergie de deformation par rapport aux depla-
cements nodaux :
Edef
fnod,i = (1.18)
ui

2 Exemple 1 : treillis soumis a une force nodale

2.1 Enonce
Le treillis plan a noeuds articules represente sur la figure 4 est com-
pose de trois poutres de meme nature et de meme section droite.

Soient E le module de Young du materiau et A laire des sections

droites.

Le noeud 1 est articule et le nud 3 repose sur un appui simple dont

la normale est horizontale.

Le noeud 2 porte une charge de composantes (0, P ).

Application numerique : on donne :

A = 100 mm2 , L = 0.2 m , E = 200000 MPa , P = 10000 N Figure 4 Exemple 1

2.2 Partition des degres de liberte

Effectuons une partition des degres de liberte en deplacements connus et inconnus ([1], [12]) :


u2

v2
u2 u1



{UL } v3
{UL } = v2 , {US } = v1 dou {U } = =
{US } u1

v3 u3




v 1

u3

On en deduit la localisation des degres de liberte dans les matrices globales :

u1 0

v 0


1


u2 1
{DDL} =

v2 2



u3 0


v3 3
Treillis plans a nuds articules 7

2.3 Etude elementaire

coordonnees nodales :
nud x y
1 0 L
2 L 0
3 0 L
element 1 2 :
1 1
caracteristiques : L 2 , A , E , nx = , ny =
2 2


u1 0 1 1 1 1
EA
v1 0 1 1 1 1
{ddl12 } = , [k12 ] =

u 1
2 2L 1 1 1 1
2
v2 2 1 1 1 1
element 3 1 :
caracteristiques : 2 L , A , E , nx = 0 , ny = 1


u3 0 0 0 0 0
EA
v3 3 0 1 0 1

{ddl31 } = , [k31 ] =

u 0 2 L 0 0 0 0
1
v1 0 0 1 0 1
element 3 2 :
1 1
caracteristiques : L 2 , A , E , nx = , ny =
2 2


u3 0 1 1 1 1

v3 3 EA 1 1 1 1
{ddl32 } = , [k32 ] =

u 1
2 2L 1 1 1 1
2
v2 2 1 1 1 1

2.4 Assemblage et calcul des deplacements inconnus

Les deplacements inconnus sont les solutions de lequation [KLL ]{UL } = {Fnod,L } :

2 0 1 u2 0
EA
0 2 1 v2 = P
2 2 L 1 1 1 + 2 v 0
3

dou (programme exemple 1) :

PL PL
u2 = = 0.050 mm , v2 = (1 + 2 2) = 0.191 mm
2 EA 2 EA
PL
v3 = = 0.100 mm
EA

2.5 Efforts normaux dans les elements

Ils sont calcules a laide de la formule (1.8) :

EA 1 1 P 2
N12 = u2 v2 = = 7071 N
L 2 2 2 2
EA P
N31 = ( v3 ) = = 5000 N
2L 2

EA 1 1 P 2
N32 = u2 + ( v2 v3 ) = = 7071 N
L 2 2 2 2
8 Methode des elements finis

2.6 Actions de liaison

Les actions de liaisons sont calculees a partir des efforts normaux (equation 1.2) :

nud 1 :
F~1 = N12 ~n12 + N31 ~n31 dou

1 P 1
F1x = N12 = = 5000 N , F1y = N12 + N31 = P = 10000 N
2 2 2
nud 3 :
F~3 = N31 ~n31 N32 ~n32 dou

1 P
F3x = N32 = = 5000 N
2 2
Remarque : lequilibre de la structure est verifie :

F1x + F2x + F3x = 0 , F1y + F2y + F3y = 0 , 2 LF1x LF2x + LF2y = 0

3 Exemple 2 : treillis soumis a une force nodale

3.1 Enonce
Le treillis plan represente sur la figure 5 est compose de trois
poutres de meme section.

Soient E le module de Young du materiau et A laire des sec-

tions droites.

Le nud 1 est articule et le nud 2 repose sur un appui simple

dont la normale est horizontale.

Le nud 3 porte une charge dintensite (P, 3 P, 0).

Figure 5 Exemple 2
Application numerique : on donne :

A = 100 cm2 , E = 200000 MPa , L = 0.7 m , P = 120 kN

3.2 Partition des degres de liberte

Effectuons une partition des degres de liberte en deplacements connus et inconnus ([1], [12]) :


v2
u

3
v2 u1
{UL } v3
{UL } = u3 , {US } = v1 dou {U } = =
{US } u1

v3 u2


v1


u2
Treillis plans a nuds articules 9

On en deduit la localisation des degres de liberte dans les matrices globales :

u1 0

v 0


1


u2 0
{DDL} =

v2 1




u3 2


v3 3

3.3 Etude elementaire

coordonnees nodales :
nud x y
1 0 0
2 0 L
3 L 0
element 1 2 :
caracteristiques : L , A , E , nx = 0 , ny = 1


u1 0 0 0 0 0
EA
v1 0 0 1 0 1
{ddl12 } = , [k12 ] =

u 0

L 0 0 0 0
2
v2 1 0 1 0 1
element 1 3 :
caracteristiques : L , A , E , nx = 1 , ny = 0


u 1 0
1 0 1 0

v1 0 EA 0 0 0 0

{ddl13 } = , [k13 ] =

u 2
L 1 0 1 0
3
v3 3 0 0 0 0
element 2 3 :
1 1
caracteristiques : L 2 , A , E , nx = , ny =
2 2


u2 0 1 1 1 1
EA
v2 1 1 1 1 1

{ddl23 } = , [k23 ] =

u 2 2 2 L 1 1 1 1
3
v3 3 1 1 1 1

3.4 Assemblage et calcul des deplacements inconnus

Les deplacements inconnus sont les solutions de lequation [KLL ]{UL } = {Fnod,L } :

1+2 2 1 1 v2 0
EA
1 1 + 2 2 1 u3 = P
2 2L
1 1 1 v3 3P

dou (programme exemple 2) :

3PL
v2 =
= 0.126 mm
EA

4PL (7 + 6 2) P L
u3 = = 0.168 mm , v3 = = 0.650 mm
EA EA
10 Methode des elements finis

3.5 Efforts normaux dans les elements

Ils sont calcules a laide de la formule (1.8) :
EA
N12 = v2 = 3 P = 360 kN
L
EA
N13 = u3 = 4 P = 480 kN
L

EA 1 1
N23 = u3 ( v3 v2 ) = 3 2 P = 509 kN
L 2 2 2

3.6 Actions de liaison

Les actions de liaisons sont calculees a partir des efforts normaux (equation 1.2) :

nud 1 :
F~1 = N12 ~n12 N13 ~n13 dou
F1x = N13 = 4 P = 480 kN , F1y = N12 = 3 P = 360 kN
nud 2 :
F~2 = N12 ~n12 N23 ~n23 dou
1
F2x = N23 = 3 P = 360 kN
2
Remarque : lequilibre de la structure est verifie :

F1x + F2x + F3x = 0 , F1y + F2y + F3y = 0 , LF3y LF2x = 0

4 Exemple 3 : treillis soumis a une variation de temperature

4.1 Enonce
Le treillis plan a noeuds articules represente sur la figure 6 est compose de trois poutres de meme
materiau et de meme section droite.

Figure 6 Exemple 3

Soient E et respectivement le module de Young et le coefficient de dilatation du materiau.

Soit A laire des sections droites.

Treillis plans a nuds articules 11

Les nuds 1, 2 et 4 sont lies a lexterieur par une rotule.

Premier cas de charge : la structure est soumise a une variation de temperature T .

Deuxieme cas de charge : la poutre (2 3) est soumise a une variation de temperature T .

Application numerique : on donne :

A = 100 mm2 , L = 0.1 m , E = 200000 MPa , = 105 K1 , T = 100 K

4.2 Partition des degres de liberte

Effectuons une partition des degres de liberte en deplacements connus et inconnus ([1], [12]) :


u3
u1
v



3


v1

u


1

u3 u2 {UL } v1
{UL } = , {US } = dou {U } = =
v3 v2
{US } u2








u4

v2


v4
u4

v4

On en deduit la localisation des degres de liberte dans les matrices globales :

u1 0
v 0



1


u 0

2



v2 0
{DDL} =

u3 1



v3 2




u4 0


v4 0

4.3 Etude elementaire

coordonnees nodales :
nud x y
1 0 0
2 L 0
3 L L
4 2L L
element 1 3 :
1 1
caracteristiques : L 2 , A , E , , nx = , ny =
2 2


u1 0 1 1 1 1
1
EA
v1 0 1 1 1 1 {fth,13 } = 1 EA T 1
{ddl13 } = [k13 ] =

u 1
2 2L 1 1 1 1 2 1

3
v3 2 1 1 1 1 1
element 2 3 :
caracteristiques : L , A , E , , nx = 0 , ny = 1
12 Methode des elements finis


u 2 0
0 0 0 0
0
EA
v2 0 0 1 0 1
1
{ddl23 } = [k23 ] = {fth,23 } = EA T

u 1 L 0 0 0 0
0
3
v3 2 0 1 0 1 1
element 3 4 :
caracteristiques : L , A , E , , nx = 1 , ny = 0


u3 1 1 0 1 0
1

v3 2 EA 0 0 0 0 {fth,34 } = EA T 0
{ddl34 } = [k34 ] =

u 0
L 1 0 1 0
1
4
v4 0 0 0 0 0 0

4.4 Premier cas de charge

4.4.1 Assemblage et calcul des deplacements inconnus
Les deplacements inconnus sont les solutions de lequation [KLL ]{UL } = {Fth,L } :

EA 1 + 2 2 1 u3 1 2 2
= EA T
2 2L 1 1 + 2 2 v3 2 2+2
dou (programme exemple 3) :

u3 = ( 2 2) L T = 0.0586 mm , v3 = 2 L T = 0.1414 mm

4.4.2 Efforts normaux dans les elements

Ils sont calcules a laide de la formule (1.8) :

EA 1 1
N13 = u3 + v3 EA T = ( 2 2) EA T = 11716 N
L 2 2 2
EA
N23 = v3 EA T = ( 2 1) EA T = 8284 N
L
EA
N34 = (u3 ) EA T = (1 2) EA T = 8284 N
L

4.4.3 Actions de liaison

Les actions de liaisons sont calculees a partir des efforts normaux (equation 1.2) :

nud 1 :
F~1 = N13 ~n13dou :
1
F1x = N13 = ( 2 1) EA T = 8284 N , F1y = F1x = 8284 N
2
nud 2 :
F~2 = N23 ~n23 dou :

F2x = 0 , F2y = N23 = (1 2) EA T = 8284 N
nud 4 :
F~4 = N34 ~n34 dou :

F4x = N34 = (1 2) EA T = 8284 N , F4y = 0
Remarque : lequilibre de la structure est verifie :
F1x + F2x + F4x = 0 , F1y + F2y + F4y = 0 , L F2y L F4x = 0
Treillis plans a nuds articules 13

4.5 Deuxieme cas de charge

4.5.1 Assemblage et calcul des deplacements inconnus
Les deplacements inconnus sont les solutions de lequation [KLL ]{UL } = {FL } :

EA 1 + 2 2 1 u3 0
= EA T
2 2L 1 1 + 2 2 v3 1

dou (programme exemple 3) :

1 2 3 2
u3 = L T = 0.0207 mm , v3 = L T = 0.0793 mm
2 2

4.5.2 Efforts normaux dans les elements

Ils sont calcules a laide de la formule (1.8) :

EA 1 1 2 2
N13 = u3 + v3
= EA T = 5858 N
L 2 2 2 2

EA 1 2
N23 = v3 EA T = EA T = 4142 N
L 2

EA 21
N34 = (u3 ) = EA T = 4142 N
L 2

4.5.3 Actions de liaison

Les actions de liaisons sont calculees a partir des efforts normaux (equation 1.2) :

nud 1 :
F~1 = N13 ~n13 dou :

1 1 2
F1x = N13 = EA T = 4142 N , F1y = F1x = 4142 N
2 2
nud 2 :
F~2 = N23 ~n23 dou :

21
F2x = 0 , F2y = N23 = EA T = 4142 N
2
nud 4 :
F~4 = N34 ~n34 dou :

21
F4x = N34 = EA T = 4142 N , F4y = 0
2
Remarque : lequilibre de la structure est verifie :

F1x + F2x + F4x = 0 , F1y + F2y + F4y = 0 , L F2y L F4x = 0

14 Methode des elements finis

5 Programmes Maple
Les programmes suivant se trouvent dans le fichier treillis.txt.

5.1 tre mat

# calculs elementaires
restart:with(linalg):

# allongement unitaire
eps:=(nx*(uj-ui)+ny*(vj-vi))/L;

# energie de deformation
Edef:=EA*eps^2*L/2-eps*EA*alpha*DT*L;

# matrice de rigidite
k:=hessian(Edef,[ui,vi,uj,vj]);

# efforts nodaux
fnod:=grad(Edef,[ui,vi,uj,vj]);

# remarque
k:=jacobian(fnod,[ui,vi,uj,vj]);

# vecteur d^
u au gradient thermique
fth:=-jacobian(fnod,[DT]);

5.2 exemple 1
restart:with(linalg):

# application numerique

#L:=200;E:=200000;A:=Pi*30^2/4;P:=-10000;

# matrice de rigidite

KL:=matrix([[2,0,-1],[0,2,-1],[-1,-1,1+sqrt(2)]]):
KL:=scalarmul(KL,E*A/2/sqrt(2)/L);

# vecteur FL

FL:=vector([0,P,0]);

# calcul des deplacements nodaux

UL:=linsolve(KL,FL);

#evalf(%);
Treillis plans a nuds articules 15

5.3 exemple 2
restart:with(linalg):

# application numerique

# L:=700;E:=200000;A:=10000;P:=-120e3;

# matrice de rigidite
KL:=matrix([[1+2*sqrt(2),1,-1],[1,1+2*sqrt(2),-1],[-1,-1,1]]):
KL:=scalarmul(KL,E*A/2/sqrt(2)/L);

# vecteur FL
FL:=vector([0,P,3*P]);

# calcul des deplacements nodaux

UL:=linsolve(KL,FL);
#evalf(%);

5.4 exemple 3
restart:with(linalg):

# application numerique

#L:=100;E:=200000;A:=100;alpha:=1e-5;DT:=100;

# matrice de rigidite

x:=1+2*sqrt(2):
KL:=matrix([[x,1],[1,x]]):
KL:=scalarmul(KL,E*A/2/sqrt(2)/L);

# vecteurs F
x:=E*A*alpha*DT/2:
FLcas1:=vector([(sqrt(2)-2)*x,(sqrt(2)+2)*x]);
FLcas2:=vector([0,E*A*alpha*DT]);

# calcul des deplacements nodaux

ULcas1:=linsolve(KL,FLcas1);
ULcas2:=linsolve(KL,FLcas2);

# map(evalf,ULcas1);
# map(evalf,ULcas2);
16 Methode des elements finis

References
[1] J. H. Argyris et H.-P. Mlejnek Die methode der finiten elemente, Band I. Verschiebung-
smethode in der statik, Vieweg, 1986.
[2] J.-L. Batoz et G. Dhatt Modelisation des structures par elements finis, Volume 1. Solides
elastiques, Hermes, 1990.
[3] , Modelisation des structures par elements finis, Volume 2. Poutres et plaques, Hermes, 1990.
[4] A. Bazergui, T. Bui-Quoc, A. Biron, G. McIntyre et C. Laberge Resistance des
materiaux, 3 ed., Editions de lEcole Polytechnique de Montreal, 2002.
[5] L. Chevalier Mecanique des systemes et des milieux deformables. Cours, exercices et pro-
blemes corriges, Ellipses, 2004.
[6] G. Dhatt, G. Touzot et E. Lefrancois Methode des elements finis, Hermes, 2005.
[7] F. Frey Traite du genie civil, Volume 1. Analyse des structures et milieux continus. Statique
appliquee, Presses Polytechniques et Universitaires Romandes, 1998.
[8] , Traite du genie civil, Volume 2. Analyse des structures et milieux continus. Mecanique des
structures, Presses Polytechniques et Universitaires Romandes, 2000.
[9] F. Frey et J. Jirousek Traite du genie civil, Volume 6. Methode des elements finis, Presses
Polytechniques et Universitaires Romandes, 2001.
[10] D. Gay et J. Gambelin Une approche simple du calcul des structures par la methode des
elements finis, Hermes, 1989.
[11] , Dimensionnement des structures. Une introduction, Hermes, 1999.
[12] J.-F. Imbert Analyse des structures par elements finis, 3 ed., Cepadues, 1995.
[13] S. Laroze Mecanique des structures, Tome 2. Theorie des poutres, 2 ed., Eyrolles/Masson,
1988.
[14] A. Portela et A. Charafi Finite elements using Maple. A Symbolic Programming Approach,
Springer, 2002.
[15] J. S. Przemieniecki Theory of matrix structural analysis, Dover, 1986.
[16] W. Weaver et J. M. Gere Matrix analysis of framed structures, 3 ed., Van Nostrand Rein-
hold, 1990.
[17] C. Wielgoz Cours et exercices de resistance des materiaux : elasticite, plasticite, elements
finis, Ellipses, 1999.
[18] W. Wunderlich et W. D. Pilkey Mechanics of structures. Variational and computational
methods, 2 ed., CRC PRESS, 2003.