Near East University

Department of Civil Engineering

Fall 2023-2024

CIV581-Matrix Methods of Structural Analysis

Assignment No.1
 DEMONSTRATION

1. Elément finis en résistance des matériaux

Considérons cet élément de la structure :

F FL L
En résistance des matériaux : σ = et dl= =F
A AE AE

2. La structure unidimensionnelle

 Equation d’une droite passant par 2 points : u = α 1+ α 2 x (1)

α1
Ecrivons l’équation (1) sous forme matricielle : u = ( 1 x ) ( ) (2)
α2

Avec ; A : aire, E : module d’élasticité et u : une fonction linéaire passant par deux points.

On sait que : {
si , X=0 →u=u1 → u1=α 1 +0=α 1
Si, X =L→ u=u 2 → u2=α 1+ α 2 L

Ecrivons ces équations sous forme matricielle :

u1 1 0 α1 α1 u1 1 0 −1 1 0 −1 1
( )=( ) ( ) →( ) = ( ) * ( ¿ ; or ( ¿ = ¿)
u2 1 L α2 α2 u2 1 L 1 L L

1 0 −1
Démonstration de ( ¿
1 L

−1 1
A = adj A ; Avec Adj A = t com A
⌈ A⌉

+ L −1 L −1 L 0
Com A = ( )=( ) → t com A = ( ) = Adj A
−0 +1 0 1 −1 1

1 0
⌈ A ⌉=⌈ ⌉ = (1*L) – (1*0) = L
1 L
A
−1
= (
1 L 0 1
= ¿)
L −1 1 L )
α1 u1 1 α1 1 0 u
Donc : ( ) = ( ) * ¿) → ( ) = ( −1 1 ) ( 1) (3)
α2 u2 L α2 u2
L L

1 0 u
L’équation (3) dans (2) : u = ( 1 x ) ( −1
L L
u2 (x
1 ) ( 1) (4) →u = (1− )
L )
x u1
( )
L u2

1 0 u
∂u
Dérivons l’équation (4) par rapport à x → ε x = =¿ ( 0 1 ) ( −1 1 ) ( 1) (5)
∂x u2
L L
1 0
Posons : G = ( 1 x ) ; A −1
= ( −1
L L
u
u2 ( x
1 ) ; δ = ( 1) ; N = (1− )
L
x
L )
and Bx = ( 0 1 )

Nous aurons : ε x = Bx * A−1 * δ

 La loi de Hooke : σ =¿E * ε (6)

1 0 u
L’équation (5) dans (6) : σ =¿E ( 0 1 ) ( −1 1 ) ( 1) (7)
u2
L L

Nous aurons : σ =¿E * Bx * A−1 * δ

 Théorème de Castigliano de l’énergie de la déformation : δU

1 1
∫ ε T σ dv = 2 ∫ (B x∗A
−1 T
δU = ∗δ) E∗B x∗A−1∗δ dv
2

1
δU =
2 ∫ B xT ¿ A−1T ¿ δT * E∗B x∗A−1∗δ dv
On sait que : dv = A dx

1
δU =
2 ∫ B xT ¿ A−1T ¿ δT * E∗B x∗A−1∗δ∗A dx
1
¿ A ¿ δ ¿ δ∗A E ∫ Bx ¿ Bx ∗A dx
T −1
δU = −1T T
2

( )
−1
1 0 1
L u1
On sait que : A = ( −1
−1
1)→ A =
−1T
; δ = ( ) → δ T = (u1 ¿ u2 )
1 u2
L L 0
L
¿ Bx = ( 0 1 ) → Bx T = (01)
1 0 0 0
0
= ( ) ( 0 1 ) ( −1
T −1
- Bx ∗Bx ∗A 1 ) = ( −1 1)
1
L L L L

1 0 1 0 0 0
0 0 0
Démonstration : Bx ∗Bx ∗A T −1
= ( ) ( 0 1 ) ( −1 1)=( ) ( −1 1 ) = ( −1 1)
1 0 1
L L L L L L

Nous introduisons les bornes : (0 ; L)

( )
L L 0 0 0 0 0 0
∫B dx = ∫ −1
T −1
- x ∗B x∗A 1 dx = ( −1 1 ) ¿ = ( −1 1 ) (L 0)
0
0 L L L L L L
0 0
=( )
−1 1

( )
−1
1
1 0 0 1 L 0 0
δU = ¿A
−1T
¿ δ ¿ δ∗A E ( −1 1) =
T
( ) T ¿ δ∗A E
2 2 1 −1 1 ¿ δ
0
L

( ) ( )
−1 1 −1
1
L 0 0 L L 1 1 −1
- ( )= = ( )
1 −1 1 −1 1 L −1 1
0
L L L

δU =
1 1 1 −1
*
2 L −1 1
¿(δ
T
¿ δ∗A E = )
AE 1 −1
2 L −1 1
u
(u1 ¿ u2 ) ( 1)
u2 ( )
δU =
AE
2L
u −u
(u1 ¿ u2 ) 1 2 =
(
−u 1+u 2
AE
)
(u ( u −u ) +u 2 (−u1 +u2) )
2L 1 1 2

AE
δU = (u12−2 u1 u2+ u22)
2L

δU AE δU AE
= (2u1−2 u2) = Q1 → = (u1−u 2)
δ u1 2 L δ u1 L

δU AE δU AE
= (- 2u1 +2 u2) = Q2 → = (- u1 +u2)
δ u2 2 L δ u2 L

- Sous forme matricielle

()
δU
δ u1
δU
=
AE 1 −1 u1
( ¿ =
L −1 1 u2
Q1
Q2( ) ( )
δ u2

L’équation finale est :

( )
Q1
Q2
=
AE 1 −1
L −1 1 ( ) (uu )
1

3. Nœud à 2 degré de liberté

x’ et y’ : système de coordonnée locale

X et Y : système de coordonnée globale

( ) )( )
'
xi
y
'
j
=¿ (
AE 1 −1 u ' i
L −1 1 u' j
(8)

Enfin de permettre la transformation des coordonnées locales en coordonnées globales,

il faut que la matrice de l’équation (8) qui est 2X2 soit exprimée en matrice de 4X4 en
incluant les déplacements v ' i et v ' j ,on aura :
Transformation de la machine de rigidité à partir de l’équation (8)

x i = x i cosθ+ y i sinθ
'

y i = - x i sinθ+ y i cosθ
'

x j = x j cosθ + y j sinθ
'

y j = - x j sinθ+ y j cosθ
'

Posons : cosθ = C et sinθ = S

Nous aurons :

x i = xi C + yi S
'

y i = - x i S+ y i C (11)
'

x j = x j C+ y j S
'

y j = - x j S + y j C (12)
'

Ecrivons l’équation (11) et (12) sous forme matricielle, on aura :

Donc la transposée de la matrice T, T T sera :

On sait que, si une matrice est symétrique, alors M −1=M T ; or la matrice T est
symétrique, doncT −1=T T , et nous aurons :

[ X ' ]= [ T −1 ] [ X ] ¿ [ T T ] [ X ] (16)

Similairement, les vecteurs déplacements se transforment de coordonnées locales en

coordonnées globales de la même façon que les forces : [ ∆ ' ]= [ T ][ ∆ ] (17)

L’équation (14) et (17) dans (10) : [ T ][ X ] =[ K ' ][ T ] [ ∆ ] → [ X ] = [ T −1 ] [ K ' ][ T ] [ ∆ ]

[X]
→ =[ T −1 ] [ K ' ][ T ]
[ ∆]
[X]
D’après l’équation (10) : [ K ] =¿ et posons T = R, on aura :
[ ∆]

[ K ] =¿ [ R−1 ] [ K ' ][ R ] (18)

[ K ] : La matrice de rigidité en coordonnée globale

[ K ' ] : La matrice de rigidité en coordonnée locale
[ T ] = [ R ] : La matrice de transformation.

4. La démonstration de la matrice [ K ]

[ K ] =¿ [ R−1 ] [ K ' ][ R ]

Substituons l’équation (13) et (15) dans l’équation (18), on aura :