RÉPUBLIQUE DÉMOCRATIQUE DU CONGO
MINISTÈRE DE L’ENSEIGNEMENT SUPÉRIEUR ET UNIVERSITAIRE
Université Loyola du Congo
Faculté d’Ingénierie ULC-Icam
Département de Maintenance et Génie Industriel
Mission catholique KIMWENZA
Commune de Mont-Ngafula
BP 3724 KINSHASA – GOMBE
Cours de math
Travail pratique numéro 4 du 24 novembre 2025
Rédigé par :
MABILA NSOLE Bavon
MANZAMBI KYANZILA Stéphane
MUHIYA KABIYA Nicolas
MUSHAGALUSA KOKO Josué
Étudiants en A3/I3/L3 Maintenance et Génie industriel
Encadré par : Dr. Demco MUKONGO
Année académique : 2025 – 2026
Exercice 1 : Two Masses on Springs (masses-ressorts)
(a) Modélisation du système
On considère deux masses verticales m1 et m2 reliées par des ressorts. Les déplacements verti-
caux vers le bas sont notés :
y1 (t), y2 (t),
mesurés par rapport à la position d’équilibre. Cette convention élimine le poids des équations
puisque les efforts gravitaires sont compensés à l’équilibre.
Données :
m1 = m2 = 1, k1 = 3, k2 = 2.
Forces agissant sur m1
Ressort fixé au plafond :
Fk1 = −k1 y1 .
Ressort reliant m1 et m2 :
Fk2 →m1 = k2 (y2 − y1 ).
En appliquant la deuxième loi de Newton :
y1′′ = −k1 y1 + k2 (y2 − y1 ) = −5y1 + 2y2 .
Forces agissant sur m2
Seul le ressort k2 agit :
Fk2 →m2 = −k2 (y2 − y1 ).
Donc :
y2′′ = 2y1 − 2y2 .
Système final
!
y ′′ = −5y + 2y , −5 2
1 1 2
y′′ = Ay, A= .
y ′′ = 2y − 2y ,
2 1 2
2 −2
(b) Recherche des modes propres
On cherche des solutions harmoniques :
y(t) = xeiωt .
1
Substitution :
(A + ω 2 I)x = 0, det(A + ω 2 I) = 0.
Calcul du déterminant :
!
−5 + ω 2 2
det = (−5 + ω 2 )(−2 + ω 2 ) − 4.
2 −2 + ω 2
Développement :
(−5 + ω 2 )(−2 + ω 2 ) = 10 − 7ω 2 + ω 4 .
⇒ ω 4 − 7ω 2 + 6 = 0.
On pose λ = ω 2 . L’équation devient :
λ2 − 7λ + 6 = 0 ⇒ (λ − 1)(λ − 6) = 0.
On en déduit les pulsations :
√
ω1 = 1, ω2 = 6 .
Vecteurs propres
Mode 1 : ω 2 = 1 ! !
−4 2 x1
(A + I)x = 0, = 0.
2 −1 x2
La relation devient :
!
1
x2 = 2x1 ⇒ x(1) = .
2
Mode 2 : ω 2 = 6 ! !
1 2 x1
(A + 6I)x = 0, = 0.
2 4 x2
On obtient :
!
(2) −2
x1 = −2x2 ⇒ x = .
1
Solution générale
√ √
y(t) = x(1) (C1 cos t + C2 sin t) + x(2) C3 cos( 6 t) + C4 sin( 6 t) .
Explicitement :
2
√ √
y1 (t) = C1 cos t + C2 sin t − 2C3 cos( 6t) − 2C4 sin( 6t),
√ √
y2 (t) = 2C1 cos t + 2C2 sin t + C3 cos( 6t) + C4 sin( 6t).
Détermination des constantes
Conditions initiales :
y1 (0) = a, y2 (0) = b, y1′ (0) = u, y2′ (0) = v.
Àt=0:
C1 − 2C3 = a, 2C1 + C3 = b.
Résolution :
a + 2b b − 2a
C1 = , C3 = .
5 5
Vitesses initiales :
√ √
C2 − 2 6 C4 = u, 2C2 + 6C4 = v.
u + 2v v − 2u
C2 = , C4 = √ .
5 5 6
(c) Analyse des effets des conditions initiales
Mode 1 : déplacement en phase
x(1) = (1, 2) ⇒ y2 = 2y1 .
Mode 2 : déplacement en opposition de phase
x(2) = (−2, 1) ⇒ y2 = − 21 y1 .
Cas particuliers
— Si b = 2a et v = 2u, alors C3 = C4 = 0 seul le mode 1 apparaît.
— Si b = − a et v = − u, alors C = C = 0 seul le mode 2 est excité.
1
2
1
2 1 2
— Dans la plupart des cas : superposition des deux modes battements et échanges d’énergie.
3
Exercice 2 : Mixing Problem
Deux cuves de 200 gallons échangent du fluide. Conditions initiales :
y1 (0) = 100, y2 (0) = 200.
Concentrations :
y1 y2
c1 = , c2 = .
200 200
Cuve T1
Entrée provenant de T2 :
y2
4c2 = .
50
Sorties :
y1 4y1
4c1 = , 16c1 = .
50 50
Débit total sortant :
y1
.
10
Équation différentielle :
y2 y1
y1′ = − .
50 10
Cuve T2
Entrée en provenance de T1 :
4y1
16c1 = .
50
Sorties :
3y2 y2
, .
50 50
Débit total sortant :
4y2
.
50
Équation :
4y1 4y2
y2′ = − .
50 50
Système final
y1′ = − 10
1
y1 + 1
y,
50 2
y ′ =
4
y − 4
y.
2 50 1 50 2
4
Matrice associée : !
−0.1 0.02
A= .
0.08 −0.08
Valeurs propres
Équation caractéristique :
λ2 + 0.18λ + 0.0064 = 0.
Discriminant :
∆ = 0.0068.
λ1 = −0.12, λ2 = −0.06.
Vecteurs propres
! ! ! !
1 1 1 1
v1 = ∼ , v2 = ∼ .
−1 −2 2 1
Solution générale
! ! !
y1 (t) 1 1
= C1 e−0.12t + C2 e−0.06t .
y2 (t) −2 1
En imposant les conditions initiales :
C1 = −100, C2 = 200.
Solutions finales :
y1 (t) = −100e−0.12t + 200e−0.06t
y2 (t) = 200e−0.06t + 200e−0.12t