REPUBLIQUE DU CAMEROUN REPUBLIC OF CAMEROUN

Paix-Travail-Patrie Peace-Work-Fatherland
MINISTERE DE L’ENSEIGNEMENT SUPERIEUR
MINISTRY OF HYER EDUCATION
**********************************************
UNIVERSITE DE DOUALA
UNIVERSITY OF DOUALA
**************************
ECOLE NATIONALE SUPERIEUR POLYTECHNIQUE DE DOUALA
NATIONANAL POLYTECHNIC OF DOUALA

DEPARTEMENT DE TECHNOLOGIE ET CONSTRUCTION

INDUSTRIEL
TCI 4

EC : CALCUL DYNAMIQUE

THEME :
Discrétisation de l’équation des déplacements
transversales d’une poutre ainsi que l’équation des
déplacements d’une plaque en vue d’une résolution
numérique

Rédigé par : CHEGUEGUIM Fleury Noel (19G00378)

Enseignant : Dr ANAGUE Année académique : 2020/2021

TAF
I. Discrétisation de l’équation de vibration transversale d’une poutre.
Sans perturbation et dissipation la dynamique du système est donné par :

Page 1
 2 4
∂ y 2∂ y
2
+ω o 4 =0(1)
∂t ∂x
EI 2
Avec ω o=
ρS
La discrétisation se fait ainsi :

( )
k+1 k k−1
2
∂ y y i −2 y i + y i
2
= 2
(2)
∂t i ∆t

( )
k k k k k k
4
∂ y y i+2−4 y i+1 +6 y i −4 y i−1 + y i−2
4
= 4
(3)
∂x ∆x
(2) et (3) dans (1) nous permet d’avoir connaissant les conditions aux bords et le déplacement
initial, le déplacement aux autres nœuds de la poutre :

2
∆t
4 ( i+2
y −4 y i+1 +6 y i −4 y i−1 + y i−2 ) ( 4)
k +1 k k−1 2 k k k k k
y i =2 y − y
i i −ω o
∆x

En considérant la présence de la dissipation et des excitations extérieures, l’équation générale

des vibrations transversales de la poutre se met sous la forme :
2 2 4
∂ y ∂y ∂ y
2
+2 ξ + ω 4
=F ( x ,t ) (5)
∂t ∂t o ∂x
λ
Avec ξ= où λ est le coefficient de dissipation.
2 ρs
k+ 1 k−1
∂y
∂t i
= ( )
yi − yi
2 Δt
(6)

(2), (3) et (6) dans (5) donne :

2 2
2 ξΔ t 2∆t 2∆t
y
k +1
i =
2+ ξ Δ t
( 2 y i − y i )+
k k−1
2+ξ Δ t
k−1 2
y i −ω o 4 ( i +2
y −4 y i+1 +6 y i −4 y i−1 + y i−2 ) +
k k k k k
F
( 2+ ξ Δ t ) i ,k
( 2+ξ Δ t ) ∆ x
(7)
Avec F i ,k =F (x i ,t k )

II. Discrétisation de l’équation d’une plaque

L’équation générale est donnée par :

( )
2 2 4 4 4
∂ w ∂w ∂ y ∂ w ∂ w
2
+2 ξ +ω 4
+ 2 2 + 4 =F ( x , y , t ) (8)
∂t ∂t o ∂x ∂x ∂ y ∂ y

Page 2
 w ki+1 , j+1−2 wki , j+1 +w ki−1 , j +1−2 ( w ki+1 , j−2 w ki−1, j +w ki−1 , j ) +wik+1 , j−1−2 wki , j−1+ wi−1
( )
4 k k
∂ w , j−1
2 2
= 2 2
( 9)
∂x ∂ y Δx Δ y
En substituant (2), (3), (6) et (9) dans (8) et après arrangement nous obtenons :

{
k k k k k k k
2 ξ Δt 2 2∆t
2
wi+ 2, j−4 wi+1 , j +6 w i , j−4 wi−1 , j +w i−2 , j wi+1 , j+ 1−2 wi , j+1
w
k+1
i, j =
2+ξ Δ t
( 2 w i , j−w i , j ) +
k k−1
2+ ξ Δ t
k−1
wi , j −ωo
2+ ξ Δ t 4
+
∆x

III. Discrétisation de l’équation d’une membrane

( )
2 2 2
∂u 2 ∂ u ∂ u
2
=c 2
+ 2 (11)
∂t ∂x ∂ y

En substituant (2) dans (11) nous obtenons après arrangement que :

( )
k k k k k k
c u −2ui , j +ui−1 , j ui , j+1−2ui , j +ui , j−1
2
k +1 k k−1
u i, j =2 u −u
i, j i, j + 2 i+ 1, j +
Δt ∆ x2 ∆ y2

Page 3