REPUBLIQUE DU CAMEROUN REPUBLIC OF CAMEROON
******** ********
Paix-travail-patrie Peace-Work-Fatherland
******** ********
UNIVERSITE DE YAOUNDE I THE UNIVERSITY OF YAOUNDE I
******** ********
FACULTE DES SCIENCES FACULTY OF SCIENCES
******** ********
DEPARTEMENT DE PHYSIQUE DEPARTEMENT OF PHYSICS
******** ********
COMPTE RENDU DES TRAVAUX PRATIQUES DE PHYSIQUE 372
TRAVAUX PRATIQUES VII : ETUDE D’UN PENDULE
SIMPLE
GROUPE DE MARDI 10H – 12H55
PRESENTE PAR : NKONGO OLANGUENA JULIENNE BRINDA
MATRICULE : 23S2187
OPTION : FONDA
ENCARDREUR : DR. PETER NYIFEH
M. BALENG SONE
RESPONSSABLE DE UE : PR. MBINACK CLEMENT
ANNEE ACADEMIQUE : 2025 –2026
SOMMAIRE
INTRODUCTION ................................................................................................................................... 3
• OSCILLATEUR LIBRE AMORTIE .............................................................................................. 4
1) Adimensionnons l’équation différentielle régissant de teta (𝜃𝜃) .................................................. 4
2) Algorithme de résolution d’une équation différentielle par la méthode d’Euler ......................... 4
3) Evolution temporelle de la solution analytique ........................................................................... 5
4) Testons les conditions initiales et du paramètre Q. ..................................................................... 6
5) Représentation du portrait de phase ............................................................................................ 6
6) Modification de la valeur de delta t et commentaire. .................................................................. 7
7) Résolution de l’équation différentielle en construisant le diagramme bloc correspondant sur
Simulink (Sim 2) ................................................................................................................................. 8
• OSCILLATEUR ENTRETENU ..................................................................................................... 9
8) Adimensionnons l’équation différentielle régissant de teta (𝜃𝜃) .................................................. 9
9) Représentation de la solution numérique 𝜃𝜃 (t) ........................................................................... 9
10) Représentation de la solution numérique du plan de phase (𝜃𝜃,V) ......................................... 10
11) Graphe du nombre d’attracteurs périodique des pendules..................................................... 10
12) Dédoublement de période ...................................................................................................... 10
13) Etude de l’etat chaotique du portrait de phase....................................................................... 11
CONCLUSION ..................................................................................................................................... 12
INTRODUCTION
L’objectif de ce travail pratique est de mettre en évidence la nature des oscillations d’un
pendule simple et d’étudier leur comportement en fonction de différents paramètres physiques.
Le pendule simple est un système oscillant constitué d’une masse ponctuelle suspendue à une
tige rigide, qui oscille sous l’effet de la gravité. Son étude permet de mieux comprendre les
phénomènes d’oscillation, d’amortissement et de résonance. Le pendule est utilisé comme
modèle de base dans l’analyse des systèmes dynamiques linéaires et non linéaire. Il est
également très utile en physique pour illustrer des concepts comme la périodicité et l’énergie
mécanique. Ses domaines d’application son vastes :il intervient en horlogerie (pendule
d’horloge), en robotique, en instrumentation et même dans certains capteur d’orientation.
Le principe repose sur la modélisation, l’analyse numérique et la visualisation de
l’évolution angulaire du pendule.
Application du travaux pratique
OSCILLATEUR LIBRE AMORTIE
Soit l’équation différentielle régissant de l’angle θ
𝑑𝑑 2 𝜃𝜃 𝑑𝑑𝑑𝑑
𝑚𝑚𝑚𝑚 = −𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝜃𝜃 – KL (1)
𝑑𝑑𝑡𝑡 2 𝑑𝑑𝑑𝑑
1) Adimensionnons l’équation différentielle régissant de teta (𝜃𝜃 )
En divisant l’équation (1) par mL on obtient :
𝑑𝑑2 𝜃𝜃 𝑔𝑔 𝐾𝐾 𝑑𝑑𝑑𝑑
= − 𝑠𝑠𝑠𝑠𝑠𝑠𝜃𝜃 –
𝑑𝑑𝑡𝑡 2 𝐿𝐿 𝑚𝑚 𝑑𝑑𝑑𝑑
On introduit les grandeurs adimensionnées suivantes :
• 𝜏𝜏=t/T0 avec T0 = √𝐿𝐿/𝑔𝑔 qui équivaut à dire que t= 𝜏𝜏 T0 => dt= T0 𝑑𝑑𝜏𝜏
• Q=m/(KT0) facteur de qualité (inverse de l’amortissement)
En remplacent dans (1) on obtient :
𝑑𝑑2 𝜃𝜃 1 𝑑𝑑𝑑𝑑
𝑑𝑑 𝜏𝜏2
+ 𝑄𝑄 𝑑𝑑 𝜏𝜏
+ 𝜃𝜃 = 0
2) Algorithme de résolution d’une équation différentielle par la méthode d’Euler
Debut
Q← ? ; h← ? ; g←? ……………..
Pour i allant de 0 à N-1 faire
teta(i+1) = teta(i)+h*V(i);
V(i+1) = V(i)+h*(-teta(i)-(1/Q) *V(i)) ;
Fin pour
Fin
Résolution numérique sur MATLAB
3) Evolution temporelle de la solution analytique
Comparaison de la solution analytique et numérique
COMPARAISON ANALYTIQUE ET THEORIQUE
1
0.8
0.6
0.4
0.2
teta
-0.2
-0.4
-0.6
-0.8
0 10 20 30 40 50 60 70 80 90 100
Temps(s)
On constate que la forme de la trajectoire est une sinusoïde et le mouvement est mouvement
oscillatoire libre amortie dont les amplitudes décroisent au cour du temps et finissent par
s’annulés.
4) Testons les conditions initiales et du paramètre Q.
Q 5 10 100 1000
REPRESENTATION de teta nalytique pour différentes valeurs de Q
a
1
0.8
0.6
0.4
0.2
0
teta
-0.2
-0.4
-0.6
-0.8
-1
0 10 20 30 40 50 60 70 80 90 100
Temps(s)
5) Représentation du portrait de phase
Pour 𝜃𝜃(0) = 1 𝑒𝑒𝑒𝑒 𝑣𝑣(0) = 4
Portrait de phase pour teta0=1 et V0=4
4
0
V
-1
-2
-3
-4
-3 -2 -1 0 1 2 3 4
theta
Pour 𝜃𝜃(0) = 1 𝑒𝑒𝑒𝑒 𝑣𝑣(0) = 4.5
Portrait de phase pour teta0=1 et V0=4.5
5
1
V
-1
-2
-3
-4
-4 -3 -2 -1 0 1 2 3 4 5
theta
On constate que la trajectoire décrite par le portrait de phase pour différentes valeurs de teta et
v sont des spirales presque identiques.
6) Modification de la valeur de delta t et commentaire.
Commentaire : d’après nos trois figures on remarque que plus le pas est grand plus on observe
une courbe linéaire telle est le cas de la dernière figure qui a un pas h=1 et moins le pas est petit
on observe des oscillations libre amorti (des amplitudes qui décroissent au cour du temps) et
qui sont continue tels sont le cas de la première figure avec un pas h=0.001. De plus lorsque le
pas est moyen les oscillations se forme en petit intervalle de temps puis disparaissent.
7) Résolution de l’équation différentielle en construisant le diagramme bloc correspondant
sur Simulink (Sim 2)
Graphe obtenue su Simulink (Sim2)
Conclusion : on constate que la courbe obtenue par Simulink est une courbe exacte des
oscillations libre amorti de l’angle teta en fonction du temps. D’où Simulink nous permet de
mieux visualiser le phénomène d’oscillation effectuer par le pendule.
OSCILLATEUR ENTRETENU
Soit l’équation différentielle régissant de l’angle θ
𝑑𝑑 2 𝜃𝜃 𝑑𝑑𝑑𝑑
𝑚𝑚𝑚𝑚 = −𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝜃𝜃 – KL + Csin𝜔𝜔𝜔𝜔 (2)
𝑑𝑑𝑡𝑡 2 𝑑𝑑𝑑𝑑
8) Adimensionnons l’équation différentielle régissant de teta (𝜃𝜃 )
En divisant l’équation (1) par mL on obtient :
𝑑𝑑2 𝜃𝜃 𝑔𝑔 𝐾𝐾 𝑑𝑑𝑑𝑑 𝐶𝐶
= − 𝑠𝑠𝑠𝑠𝑠𝑠𝜃𝜃 – + sin𝜔𝜔𝜔𝜔
𝑑𝑑𝑡𝑡 2 𝐿𝐿 𝑚𝑚 𝑑𝑑𝑑𝑑 𝑚𝑚𝑚𝑚
On introduit les grandeurs adimensionnées suivantes :
• 𝜏𝜏=t/T0 avec T0 = √𝐿𝐿/𝑔𝑔 qui équivaut à dire que t= 𝜏𝜏 T0 => dt= T0 𝑑𝑑𝜏𝜏
• Q=m/(KT0) facteur de qualité (inverse de l’amortissement)
En remplacent dans (2) on obtient :
𝑑𝑑2 𝜃𝜃 1 𝑑𝑑𝑑𝑑
𝑑𝑑 𝜏𝜏2
+ 𝑄𝑄 𝑑𝑑 𝜏𝜏
+ 𝜃𝜃 + Asin𝜔𝜔𝜔𝜔 = 0
9) Représentation de la solution numérique 𝜃𝜃 (t)
10) Représentation de la solution numérique du plan de phase (𝜃𝜃 ,V)
11) Graphe du nombre d’attracteurs périodique des pendules
On constate que pour un premier temps lorsque lorsque 𝜃𝜃(0) = −0.03 𝑒𝑒𝑒𝑒 𝑣𝑣(0) = −2.01 on
observe une droite linéaire et lorsque 𝜃𝜃(0) = −0.78 𝑒𝑒𝑒𝑒 𝑣𝑣(0) = 0.15 on observe des
oscillations un peu plus retardé( appelé oscillation entretenu) qui son continue au cour du temps
12) Dédoublement de période
On constate que dans les deux cas la courbe du dédoublement de la période décrite à la
forme des oscillations entretenu qui évolue instantanément au cour du temps.
13) Etude de l’etat chaotique du portrait de phase
On constate que à l’instant 𝜃𝜃(0) = 𝑣𝑣(0) = 0 obtenu au 10) le plan de phase n’était pas
vraiment visible mais au fur et mesure qu’on augmente l’angle le plan de phase essaye de
prendre la forme tel est le cas observé ci haut.
CONCLUSION
Ce travail pratique nous a permis de modéliser et de simuler le mouvement d’un pendule
simple, d’analyser les oscillation libres et amorties, et d’évaluer l’influence des paramètres sur
sa dynamique. La méthode numérique utilisée (Euler sous MATLAB) s’est montrée efficace
pour observer l’évolution temporelle de l’angle, les portraits de phase et la stabilité du système.
L’approche expérimentale et théorique conjointe renforce notre compréhension des systèmes
oscillants et prépare à l’’étude des systèmes plus complexes en physique et en ingénierie.
ANNEXE
close all;
clear;
clc;
%% OXILLATEUR LIBRE AMORTI %%
%%Declaration des variables
Q=5; g=9.8; L=1; h=0.01; t=0:h:100; N=length(t); omega0=sqrt(g/L);
%% Condition initial
teta0=1; V0=0;
%% allocation mémoire
teta=zeros(1,N); V=zeros(1,N); teta(1)=teta0; V(1)=V0; t(1)=0;
%% Résolution par la methode d'euler
for i=1:N-1
teta(i+1)= teta(i)+h*V(i);
V(i+1)= V(i)+h*(-teta(i)-(1/Q)*V(i));
end
figure
plot(t,teta ,'r','linewidth',2)
title('oxcilation libre amorti')
grid on
%% solution analytique
omega=omega0*sqrt(1-(1/(4*Q^2)));
teta_analy=teta0*exp(-omega*t/(2*Q)).*cos(omega*t);
figure
plot(t,teta_analy,'m','linewidth',2)
grid on
%% Superposition
figure
plot(t,teta ,'r','linewidth',2)
hold on
plot(t,teta_analy,'m','linewidth',1.5)
xlabel('Temps(s)')
ylabel('teta')
title('COMPARAISON ANALYTIQUE ET THEORIQUE')
grid on
%% Test des conditions initiales et du paarametre Q
Qv=[5 10 100 1000]; couleur=['c','r','m','g'];
figure
for j=1:4
Q=Qv(j);
teta_analy=teta0*exp(-omega*t/(2*Q)).*cos(omega*t); hold on;
plot(t,teta_analy,[couleur(j) ''],'linewidth',1.5)
xlabel('Temps(s)')
ylabel('teta')
title('REPRESENTATION de teta_analytique pour différentes valeurs de Q')
grid on
end
%% Représentation du portrait de phase pour teta0=1 et V0=4
%% Conditions initiales
teta0=1; V0=4; Q=5; h=0.01;
%% allocation mémoire
teta=zeros(1,N); V=zeros(1,N); teta(1)=teta0; V(1)=V0; t(1)=0;
%% Résolution par la methode d'euler
for k=1:N-1
teta(k+1)= teta(k)+h*V(k);
V(k+1)= V(k)+h*(-teta(k)-(1/Q)*V(k));
end
figure
plot(teta,V,'r','linewidth',2)
title('Portrait de phase pour teta0=1 et V0=4 ')
xlabel('theta')
ylabel('V')
grid on
%% Représentation du portrait de phase pour teta0=1 et V0=4.5
%% Conditions initiales
teta0=1; V0=4.5; Q=5; h=0.01;
%% allocation mémoire
teta=zeros(1,N); V=zeros(1,N); teta(1)=teta0; V(1)=V0; t(1)=0;
%% Résolution par la methode d'euler
for k=1:N-1
teta(k+1)= teta(k)+h*V(k);
V(k+1)= V(k)+h*(-teta(k)-(1/Q)*V(k));
end
figure
plot(teta,V,'c','linewidth',2)
title('Portrait de phase pour teta0=1 et V0=4.5 ')
xlabel('theta')
ylabel('V')
grid on
%% Variation du pas pour h=0.001
%% Condition initial
teta0=1; V0=0; h=0.001;
%% allocation mémoire
teta=zeros(1,N); V=zeros(1,N); teta(1)=teta0; V(1)=V0; t(1)=0;
%% Résolution par la methode d'euler
for i=1:N-1
teta(i+1)= teta(i)+h*V(i);
V(i+1)= V(i)+h*(-teta(i)-(1/Q)*V(i));
end
figure
subplot(3,1,1)
plot(t,teta ,'r','linewidth',2)
xlabel('Temps(s)')
ylabel('teta')
title('solution numérique')
legend('h=0.001')
grid on
%% Variation du pas pour h=0.1
%% Condition initial
teta0=1; V0=0; h=0.1;
%% allocation mémoire
teta=zeros(1,N); V=zeros(1,N); teta(1)=teta0; V(1)=V0; t(1)=0;
%% Résolution par la methode d'euler
for i=1:N-1
teta(i+1)= teta(i)+h*V(i);
V(i+1)= V(i)+h*(-teta(i)-(1/Q)*V(i));
end
subplot(3,1,2)
plot(t,teta ,'m','linewidth',2)
xlabel('Temps(s)')
ylabel('teta')
title('solution numérique')
legend('h=0.1')
grid on
%% Variation du pas pour h=1
%% Condition initial
teta0=1; V0=0; h=1;
%% allocation mémoire
teta=zeros(1,N); V=zeros(1,N); teta(1)=teta0; V(1)=V0; t(1)=0;
%% Résolution par la methode d'euler
for i=1:N-1
teta(i+1)= teta(i)+h*V(i);
V(i+1)= V(i)+h*(-teta(i)-(1/Q)*V(i));
end
subplot(3,1,3)
plot(t,teta ,'r','linewidth',2)
xlabel('Temps(s)')
ylabel('teta')
title('solution numérique')
legend('h=1')
grid on
%% OXILLATEUR ENTRETENU %%
%%Declaration des variables
Q=10; A=3; We=3; g=9.8; L=1; h=0.01; t=0:h:100; N=length(t);
%% Condition initial
teta0=0; V0=0;
%% allocation mémoire
teta=zeros(1,N); V=zeros(1,N); teta(1)=teta0; V(1)=V0;
%% Résolution par la methode d'euler
for j=1:N-1
teta(j+1)= teta(j)+h*V(j);
V(j+1)= V(j)+h*(-teta(j)-(1/Q)*V(j)+A*sin(We*t(j)));
end
figure
plot(t,teta ,'g','linewidth',2)
title('oxcilateur entretenu')
grid on
%% Représentation du plan de phase
for i=1:N-1
teta(i+1)= teta(i)+h*V(i);
V(i+1)= V(i)+h*(-teta(i)-(1/Q)*V(i)+A*sin(We*t(i)));
end
figure
plot(teta,V,'g','linewidth',2)
title('Plan de phase oxillateur entretenu')
xlabel('theta')
ylabel('V')
grid on
%% les attracteurs periodique dependent des amplitudes
%% Condition initial
tetab0=-0.03; Vb0=-2.01;
%% allocation mémoire
tetab=zeros(1,N); Vb=zeros(1,N); tetab(1)=tetab0; Vb(1)=Vb0;
%% Résolution par la methode d'euler
for i=1:N-1
tetab(i+1)= tetab(i)+h*Vb(i);
Vb(j+1)= Vb(j)+h*(-tetab(j)-(1/Q)*Vb(j)+A*sin(We*t(j)));
end
figure
subplot(2,1,1)
plot(t,tetab,'r','linewidth',2)
xlabel('Temps(s)')
ylabel('tetab')
title('Oxillateur entretenu pour tetab0=-0.03 et Vb0=-2.01')
grid on
%% Condition initial
tetac0=-0.78; Vc0=0.15;
%% allocation mémoire
tetac=zeros(1,N); Vc=zeros(1,N); tetac(1)=tetac0; Vc(1)=Vc0; t(1)=0;
%% Résolution par la methode d'euler
for i=1:N-1
tetac(i+1)= tetac(i)+h*V(i);
Vc(j+1)= Vc(j)+h*(-tetac(j)-(1/Q)*Vc(j)+A*sin(We*t(j)));
end
subplot(2,1,2)
plot(t,tetac,'c','linewidth',2)
xlabel('Temps(s)')
ylabel('tetac')
title('Oxillateur entretenu pour tetab0=-0.78 et Vb0=0.15')
grid on
%% phénomène de dedoublement de période
%%Declaration des variables
Q=10; A=4.9; omega=1.4; We=2*omega; g=9.8; L=1; h=0.01; t=0:h:180; N=length(t);
%% Condition initial
tetab0=-2.03; Vb0=-2.03;
%% allocation mémoire
tetab=zeros(1,N); Vb=zeros(1,N); tetab(1)=tetab0; Vb(1)=Vb0;
for i=1:N-1
tetab(i+1)= tetab(i)+h*Vb(i);
Vb(i+1)= Vb(i)+h*(-tetab(i)-(1/Q)*Vb(i)+A*sin(We*t(i)));
end
figure
subplot(3,1,1)
plot(t,tetab,'g','linewidth',2)
xlabel('Temps(s)')
ylabel('tetab')
title('dedoublement pour tetab0=-2.03 et Vb0=-2.03')
grid on
%% Condition initial
tetad0=0.42; Vd0=-3.01; omega=1.4; We=2*omega;
%% allocation mémoire
tetad=zeros(1,N); Vd=zeros(1,N); tetad(1)=tetab0; Vd(1)=Vd0;
for i=1:N-1
tetad(i+1)= tetad(i)+h*Vd(i);
Vd(i+1)= Vd(i)+h*(-tetad(i)-(1/Q)*Vd(i)+A*sin(We*t(i)));
end
subplot(3,1,2)
plot(t,tetad,'m','linewidth',2)
xlabel('Temps(s)')
ylabel('tetad')
title('dedoublement pour tetad0=0.42 et Vd0=-3.01')
grid on
%% Etude de l'etat chaotiquee
%% phénomène chaotique
%% Condition initial
tetad0=0.001; Vd0=0; omega=0.8; We=2*omega; Q=10; A=1.6;
%% allocation mémoire
tetad=zeros(1,N); Vd=zeros(1,N); tetad(1)=tetab0; Vd(1)=Vd0;
for i=1:N-1
tetad(i+1)= tetad(i)+h*Vd(i);
Vd(i+1)= Vd(i)+h*(-tetad(i)-(1/Q)*Vd(i)+A*sin(We*t(i)));
end
subplot(2,1,1)
plot(tetad,Vd,'b','linewidth',2)
xlabel('Temps(s)')
ylabel('tetad')
title('etude chaotique pour tetad0=0.001; Vd0=0')
grid on
%% Condition initial
tetaa0=1e-6; Va0=0; omega=0.8; We=2*omega; Q=10; A=1.6;
%% allocation mémoire
tetaa=zeros(1,N); Va=zeros(1,N); tetaa(1)=tetaa0; Va(1)=Va0;
for i=1:N-1
tetaa(i+1)= tetaa(i)+h*Va(i);
Va(i+1)= Va(i)+h*(-tetaa(i)-(1/Q)*Va(i)+A*sin(We*t(i)));
end
subplot(2,1,2)
plot(tetaa,Va,'c','linewidth',2)
xlabel('Temps(s)')
ylabel('tetaa')
title('etude chaotique pour tetaa0=1e-6; Vd0=0')
grid on