École Polytechnique de Ouagadougou Année Académique 2023 - 2024
MENTION : CPEI
PARCOURS : PC/L2
TP de Simulation Numérique
Durée : 02h00 - Date : Avril 2024
TP : Simulation d’un pendule simple non amorti
Objectif : Écrire un code Python pour simuler le mouvement d’un pendule simple non amorti en
utilisant les équations mathématiques du pendule et en résolvant numériquement l’équation
différentielle.
Équations du pendule simple non amorti : les équations du pendule simple non amorti peuvent être
décrites par les équations différentielles suivantes :
𝒅𝟐 𝜽 𝒈
+ 𝒔𝒊𝒏(𝜽) = 𝟎
𝒅𝒕𝟐 𝑳
Où :
- 𝜽 est l'angle entre le pendule et la verticale,
- 𝒕 est le temps,
- 𝒈 est l'accélération due à la gravité,
- 𝑳 est la longueur du pendule.
Questions I:
1. Écrivez une fonction Python pour définir l’équation différentielle du pendule simple.
2. Spécifiez les conditions initiales et les paramètres du pendule.
3. Utilisez la fonction solve_ivp de SciPy pour résoudre numériquement l’équation différentielle.
4. Tracez l’évolution de l’angle du pendule en fonction du temps à l’aide de Matplotlib.
Questions II:
1. Modifier votre code python pour le cas d’une variante avec une équation différentielle légèrement
différente, utilisant la loi de Hooke pour modéliser un pendule simple. L’équation différentielle
pour un pendule simple avec une force de rappel linéaire de la loi de Hooke est :
𝒅𝟐 𝜽 𝒌
+ 𝒔𝒊𝒏(𝜽) = 𝟎
𝒅𝒕𝟐 𝒎
Où 𝜽 est l'angle entre le pendule et la verticale, 𝒕 est le temps, 𝒌 est la constante de raideur du
ressort (force de rappel) en N/m et 𝒎 est la masse du pendule en kg.
Remarques :
1. Utilisez des unités cohérentes pour les paramètres tels que la longueur du pendule et l’accélération
due à la gravité.
2. Assurez-vous de convertir les angles en radians lorsque vous utilisez les fonctions
trigonométriques de NumPy.
3. Choisissez judicieusement l’intervalle de temps et le nombre de points de données pour la
simulation.
EPO – CPEI – PC – Simulation Numérique - Année 2023 - 2024 Page 1 sur 1
