Projet d’Analyse Numérique
Nom1:TAGBA
Nom2:AGBELA
Année Universitaire 2023-2024
0.1 Introduction
Ce rapport présente la résolution de deux problèmes d’analyse numérique liés
à un couloir coudé et à un système mécanique articulé. Chaque problème est
résolu en détail et les solutions sont implémentées en Python.
0.2 Problème 1: Couloir coudé et optimisa-
tion
0.2.1 Énoncé
On veut faire passer un grand tableau dans un couloir coudé selon un angle
γ, avec des largeurs ℓ1 et ℓ2 avant et après le coude. Quelle est la longueur
maximale L∗ d’un tableau pouvant passer le coude en glissant le long des
murs ?
0.2.2 Justification de la formule de L(α)
Pour chaque angle α ∈ [0, π − γ], la longueur maximale L(α) d’un tableau
tenant dans le couloir selon un angle α est donnée par:
ℓ1 ℓ2
L(α) = +
sin(α) sin(π − γ − α)
soit D1 et D2 telque:
L(α) = D1 + D2
Dans la première partie de la figure on a:
ℓ1
sin(α) = Dℓ11 ⇒ D1 = sin(α)
ℓ2
Dans la deuxième partie de la figure on a : sin(β) = D2
avec β + γ + α = π ⇒ β = π − γ − α
1
ℓ2
Donc D2 = sin(π−γ−α)
D’ou
ℓ1 ℓ2
L(α) = +
sin(α) sin(π − γ − α)
0.2.3 le minimum de L(α)
Nous allons chercher la dérivée séconde de L(α) et résoudre l’équation L′ (α)=0
∀α ∈ [0, π − γ], L(α) est continue et on a :
L′ (α) = −ℓsin
1 cos(α)
2 (α) + ℓsin
2 cos(π−γ−α)
2 (π−γ−α)
0.2.4 Résolution de l’équation L′ (α) = 0
Pour trouver les valeurs de α qui minimisent L(α), nous devons résoudre
l’équation L′ (α) = 0 :
−ℓ1 cos(α) ℓ2 cos(π − γ − α)
+ =0
sin2 (α) sin2 (π − γ − α)
Cela équivaut à :
−ℓ1 cos(α) sin2 (π − γ − α) = ℓ2 cos(π − γ − α) sin2 (α)
Pour résoudre l’équation
−ℓ1 cos(α) sin2 (π − γ − α) = ℓ2 cos(π − γ − α) sin2 (α)
avec ℓ1 = 1 m, ℓ2 = 1.5 m et γ = 120◦ , nous allons d’abord simplifier
l’équation en utilisant les identités trigonométriques suivantes:
1. cos(π − x) = − cos(x) 2. sin(π − x) = sin(x)
L’équation devient donc :
−ℓ1 cos(α) sin2 (γ + α) = −ℓ2 cos(γ + α) sin2 (α)
En simplifiant les signes négatifs, nous obtenons :
ℓ1 cos(α) sin2 (γ + α) = ℓ2 cos(γ + α) sin2 (α)
Substituons les valeurs de ℓ1 , ℓ2 et γ :
1 · cos(α) sin2 (120◦ + α) = 1.5 · cos(120◦ + α) sin2 (α)
Nous savons que cos(120◦ + α) = cos(120◦ ) cos(α) − sin(120◦ ) sin(α) et
sin(120◦ + α) = sin(120◦ ) cos(α) + cos(120◦ ) sin(α).
2
√
Utilisons les valeurs suivantes : - cos(120◦ ) = − 21 - sin(120◦ ) = 23
Nous obtenons donc :
√
1 3
cos(120◦ + α) = − cos(α) − sin(α)
2 2
√
◦ 3 1
sin(120 + α) = cos(α) − sin(α)
2 2
En substituant ces valeurs dans l’équation initiale, nous avons :
√ !2 √ !
3 1 1 3
cos(α) cos(α) − sin(α) = 1.5 − cos(α) − sin(α) sin2 (α)
2 2 2 2
Cette équation est complexe, nous allons simplifier en supposant α = 0
et résoudre pour vérifier :
Pour α = 0 :
√ !2 √ !
3 1 1 3
cos(0) · 1 − · 0 = 1.5 − · 1 − ·0 ·0
2 2 2 2
√ !2
3
1· =0
2
3
̸= 0
4
α = 0 n’est donc pas une solution. Nous devons donc résoudre l’équation
numériquement. Pour cela, nous pouvons utiliser des méthodes numériques
comme l’algorithme de Newton-Raphson pour trouver une solution approchée
pour α.
0.2.5 Résolution et implémentation en Python
Pour ℓ1 = 1 m, ℓ2 = 1.5 m et γ = 120◦ , nous utilisons la méthode de Newton
pour trouver le minimum de L(α).
3
4
0.2.6 Résultats
Le minimum de la fonction L(α) est obtenu pour α ≈ 0.47 radians, et la
longueur maximale admissible est L ≈ 4.96 m. Un tableau de 5 m ne peut
donc pas passer le coude.
0.3 Problème 2: Position d’un système mécanique
0.3.1 Énoncé
Un système formé de quatre barres rigides articulées entre elles est représenté
ci-dessous.
5
Si on impose l’angle β, alors le système se déforme et définit l’angle α (et
γ). En autorisant les angles α négatifs, on devine qu’il peut y avoir plusieurs
solutions. Selon les longueurs a1 , a2 , a3 , a4 , tous les angles β ne sont pas
réalisables. Par exemple, β = 0 n’est possible que si a1 + a4 < a2 + a3 .
−→ −−→ −−→ −−→
1. En exprimant l’égalité AB + BC + CD + DA = 0 en coordonnées,
écrire deux équations reliant α, β et γ. En exprimant sin2 γ + cos2 γ =
1, obtenir une relation entre α et β. En utilisant une identité pour
cos(α − β), vérifier que :
a1 a1 a22 + a23 + a24 − a21
cos β − cos α − cos(α − β) + =0
a2 a4 2a2 a4
2. On prend a1 = 10 cm, a2 = 13 cm, a3 = 8 cm, a4 = 10 cm. Pour
chaque valeur de β entre 0 et π (on prendra par exemple 20 valeurs),
utiliser la méthode de Newton pour estimer α. On partira de x0 = 2π 3
et de x0 = −1 pour tenter de chercher deux solutions différentes, que
l’on stockera dans deux vecteurs va1 et va2; lorsque la méthode de
Newton échoue, on posera va1(i)=NaN (not a number). Représenter
graphiquement va1 et va2 en fonction de β sur le même graphe (les
points NaN ne seront pas représentés).
0.3.2 Relations entre α et β
−→ −−→ −−→ −−→
Cherchons les vecteurs AB , BC , CD , DA
D’après les projections sur la figure on a:
6
−→ −a1 −−→ a2 cos(α) −−→ a3 cos(γ)
AB = , BC = , CD =
0 a2 sin(α) a3sin(γ)
−−→ a4 cos(π − β) −
−→ −a4 cos(β)
, DA = ou DA =
a4 sin(π − β) −a4 sin(β)
on a :
Cx = −a1 + a2 cos(α) + a3 cos(γ) − a4 cos(β) = 0
Cy = a2 sin(α) + a3 sin(γ) − a4 sin(β) = 0
on sait que : sin2 (γ) + cos2 (γ) = 1 (1)
et
cos(α − β) = cos(α) cos(β) + sin(α) sin(γ)(2)
Commençons par tirer a3 cos(γ) et a3 sin(γ)
a3 cos(γ) = a1 − a2 cos(α) + a4 sin(β)
a3 sin(γ) = −a2 sin(α) + a4 sin(β)
a23 cos2 (γ) = (a1 − a2 cos(α) + a4 sin(β))2
a23 sin(γ) = (−a2 sin(α) + a4 sin(β))2
a23 cos2 (γ)+a23 sin(γ) = (a1 −a2 cos(α)+a4 sin(β))2 +(−a2 sin(α)+a4 sin(β))2
D’après (1) on a:
⇒ a23 = a21 − a22 − a24 − 2a2 a1 cos(α) + 2a1 a4 cos(β) − 2a2 a4 (cos(α) cos(β) +
sin(α) sin(β))
D’après (2) on a: a23 = a21 −a22 −a24 −2a2 a1 cos(α)+2a1 a4 cos(β)−2a2 a4 cos(α−
β)
Divison le tout par 2a2 a4
on aura :
a1 a1 a2 + a23 + a24 − a21
cos β − cos α − cos(α − β) + 2 =0
a2 a4 2a2 a4
D’ou la réponse
7
0.3.3 Résolution et implémentation en Python
Pour a1 = 10 cm, a2 = 13 cm, a3 = 8 cm, a4 = 10 cm, nous utilisons la
méthode de Newton pour estimer α pour chaque valeur de β entre 0 et π.
8
0.3.4 Représenter graphiquement
0.3.5 Conclusion
Les deux problèmes du projet 2MTH1422 d’analyse numérique explorent
des concepts mathématiques et des techniques de résolution numérique
dans le contexte d’applications concrètes.Les deux problèmes illustrent la
puissance des outils mathématiques et des techniques numériques pour
résoudre des problèmes concrets et analyser des situations complexes par la
méthode de NEWTON RAPHSON.
