M. V - DYNAMIQUE - OSCILLATIONS FORCÉES


Équation du mouvement

• On considère un pendule à ressort vertical avec amortissement visqueux, mis en oscillations forcées (de pulsation ωω quelconque) par un mouvement sinusoïdal permanent du support :  zA=zA0+ZAmcos(ωt)z_A=z_{A0}+Z_{Am} \; \cos(ω \,t) .

En posant  ω02=km\displaystyle ω_0^{\:2}=\frac{k}{m}   et   α=λ2m\displaystyle α=\frac{λ}{2 \,m}  et en prenant pour origine la position d’équilibre obtenue pour  zA=zA0z_A=z_{A0}  l’équation du mouvement  peut s’écrire :
z̈+2αz˙+ω02z=ω02ZAmcos(ωt)\ddot{z}+2 \,α \:\dot{z}+ω_0^{\:2} \: z=ω_0^{\:2} \: Z_{Am} \; \cos(ω \,t) .

dynOscForc_Im/dynOscForc_Im1.jpg

◊ remarque : l'effet du fluide est décrit par   m=m+mem′=m+m_e   et   g=gmmdm+me\displaystyle \overset{→}{g}{}′=\overset{→}{g} \:\frac{m-m_d}{m+m_e}  (avec mdm_d pour le fluide “déplacé” et mem_e pour le fluide “oscillant”).

• D’une façon analogue à ce qui a été étudié en électricité, la solution générale de cette équation est la somme :

Bien que le régime transitoire soit moins vite amorti en mécanique qu’en électricité, on se limite ici à l’étude du régime forcé qu’on cherche sous la forme :   z=Zmcos(ωt+ϕ)z=Z_m \; \cos(ω \,t+ϕ)  ;  donc en notations complexes :
z_=Zmej(ωt+ϕ)=Z_mejωt\underline{z}=Z_m \: \mathrm{e}^{\mathrm{j}\,(ω\,t+ϕ)} =\underline{Z}_m \: \mathrm{e}^{\mathrm{j}ω\,t}   avec   Z_m=Zmejϕ\underline{Z}_m=Z_m \: \mathrm{e}^{\mathrm{j}ϕ} .

◊ remarque : l’amplitude ZmZ_m n'est pas une impédance… (l’analogie électromécanique permet de définir des “impédances” en mécanique ; par exemple λλ ).

◊ remarque : l’analogue électromécanique de zz est la charge  q=CuC q=C \:u_C  du condensateur ; les résultats sont donc un peu différents de ceux obtenus pour le courant  i=uRR\displaystyle i=\frac{u_R}{R}  dans un circuit “RLC” (il faudrait utiliser la variable z˙\dot{z} ).

• L’équation  z̈+2αz˙+ω02z=ω02ZAmcos(ωt)\ddot{z}+2 \,α \:\dot{z}+ω_0^{\:2} \: z=ω_0^{\:2} \: Z_{Am} \; \cos(ω \,t)   correspond à l’équation complexe  (ω2+2jαω+ω02)Z_mejωt=ω02ZAmejωt\left(-ω^2+2 \,\mathrm{j} α \:ω+ω_0^{\:2} \right) \:\underline{Z}_m \: \mathrm{e}^{\mathrm{j}ω\,t}=ω_0^{\:2} \: Z_{Am} \: \mathrm{e}^{\mathrm{j}ω\:t}   d’où on déduit par identification :  Z_m=ZAmω02ω02ω2+2jαω\displaystyle \underline{Z}_m=\frac{Z_{Am} \: ω_0^{\:2}}{ω_0^{\:2}-ω^2+2 \,\mathrm{j} α \:ω}   c'est-à-dire :
Zm=ZAmω02(ω02ω2)2+(2αω)2\displaystyle Z_m=\frac{Z_{Am} \: ω_0^{\:2}}{\sqrt{\left(ω_0^{\:2}-ω^2 \right)^2+(2 \,α \:ω)^2}}   et   sin(ϕ)=2αω(ω02ω2)2+(2αω)2\displaystyle \sin(ϕ)=-\frac{2 \,α \:ω}{\sqrt{\left(ω_0^{\:2}-ω^2 \right)^2+(2 \,α \:ω)^2}} .


Amplitude du mouvement et résonance

dynOscForc_Im/dynOscForc_Im2.jpg

• La variation de Zm(ω)Z_m (ω) donne une résonance si l’amortissement n’est pas trop grand, mais avec un maximum pour  ωrω0ω_r≠ω_0 .

En fait :  dZmdω=2ZAmω02ω.(ω02ω22α2)[(ω02ω2)2+(2αω)2]3/2\displaystyle \frac{dZ_m}{dω}=\frac{2 \,Z_{Am} \: ω_0^{\:2} \: ω .\left(ω_0^{\:2}-ω^2-2 \,α^2 \right)}{\left[\left(ω_0^{\:2}-ω^2 \right)^2+(2 \,α \:ω)^2 \right]^{3/2}}   s’annule :

• On peut définir la “bande passante” comme la largeur de la courbe de résonance pour  Zm=ZmM2\displaystyle Z_m=\frac{Z_{mM}}{\sqrt{2}} :

dynOscForc_Im/dynOscForc_Im3.jpg

On obtient ainsi :  ω42ω2ωr2+2ωr4ω04=0ω^4-2 \,ω^2 \:ω_r^{\:2}+2 \,ω_r^{\:4}-ω_0^{\:4}=0   et   ω2=ωr2±ω04ωr4ω^2=ω_r^{\:2}±\sqrt{ω_0^{\:4}-ω_r^{\:4}} .

Le calcul exact est ensuite compliqué, mais pour  X=αωr\displaystyle X=\frac{α}{ω_r}  petit :
ω2=ωr2.(1±2X1+X2)ω^2=ω_r^{\:2}.\left(1±2 \,X \:\sqrt{1+X^2}\right)  ;  ωωr.(1±X12X2±X3)ω≈ω_r .\left(1±X-\frac{1}{2} X^2±X^3\:⋯\right)  ;
ω2α.(1+X2)∆ω≈2 \,α .\left(1+X^2\:⋯\right) .

• On peut alors définir le “facteur de qualité” :  𝒬=ωrωω02α=kmλ\displaystyle 𝒬=\frac{ω_r}{∆ω}≈\frac{ω_0}{2 \,α}=\frac{\sqrt{k \:m}}{λ}  (résonance plus “aiguë” si l’amortissement est faible).

On peut en outre noter que  ZmMZAm=ω022αω02α2𝒬\displaystyle \frac{Z_{mM}}{Z_{Am}} =\frac{ω_0^{\:2}}{2 \,α \:\sqrt{ω_0^{\:2}-α^2}}≈𝒬   pour αα petit.

Déphasage

• Le déphasage varie (entre 00 et π ) de façon analogue à celle obtenue en électricité pour l’intensité du courant dans un circuit RLC (entre  π2-\frac{π}{2}  et  π2\frac{π}{2} ) :

dynOscForc_Im/dynOscForc_Im4.jpg

📖 exercices n° I et II

Puissance

• La puissance mécanique d'une force f\overset{→}{f} est :  𝒫=fv𝒫=\overset{→}{f}∙\overset{→}{v}  ;  pour un mouvement rectiligne, l'expression algébrique  𝒫=f.v𝒫=f.v  est l'analogue mécanique de la puissance instantanée électrique  p(t)=u(t)i(t)p(t)=u(t) \:i(t) .

◊ remarque : si on utilise en électricité le “produit scalaire complexe” pour noter la puissance moyenne :  P_=U_.I_*\underline{P}=\underline{U}.\underline{I}^*  on peut être tenté de considérer qu'il s'agit de l'analogue logique du produit scalaire envisagé précédemment ; cette pseudo-analogie est toutefois “bâtarde” car les puissances considérées (instantanée et moyenne) sont des grandeurs différentes.

• Pour calculer la puissance moyenne reçue par le système, on peut considérer que c'est celle fournie par le moteur qui actionne le support :
fA(t)=k.(zAz𝓁0)f_A (t)=k .(z_A-z-𝓁_0)
=k.[zA0+ZAmcos(ωt)zeqZmcos(ωt+ϕ)𝓁0]=k .\left[z_{A0}+Z_{Am} \; \cos(ω \,t)-z_{eq}-Z_m \; \cos(ω \,t+ϕ)-𝓁_0\right]  ;
vA(t)=z˙A=ωZAmsin(ωt)v_A (t)=\dot{z}_A=-ω \:Z_{Am} \; \sin(ω \,t)  ;
𝒫(t)=fA.vA=12kωZAmZm.[sin(2ωt+ϕ)sin(ϕ)]𝒫(t)=f_A .v_A=\frac{1}{2} k \:ω \:Z_{Am} \: Z_m .\left[\sin(2 \,ω \,t+ϕ)-\sin(ϕ) \right] 
12kωZAm2sin(2ωt)mgωZAmsin(ωt)-\frac{1}{2} k \:ω \:Z_{Am}^{\:2} \; \sin(2 \,ω \,t)-m \:g \:ω \:Z_{Am} \; \sin(ω \,t)  ;
P=𝒫(t)=12kωZAmZmsin(ϕ) = k α Z Am 2 ω0 2 ω2 (ω02ω2)2+(2αω)2 \displaystyle P=〈 𝒫(t) 〉=-\frac{1}{2} k \:ω \:Z_{Am} \: Z_m \; \sin(ϕ)=\frac{k \:α \:Z_{Am}^{\:2} \: ω_0^{\:2} \: ω^2)}{\left(ω_0^{\:2}-ω^2 \right)^2+(2 \,α \:ω)^2} .

◊ remarque : on considère ici que le ressort “fait partie” du système décrit par un point matériel ; son énergie cinétique est négligeable ; son énergie potentielle élastique est considérée comme propriété du “point”.

◊ remarque : avec les valeurs efficaces  P=FA𝒻𝒻VA𝒻𝒻cos(ϕA)P=F_{Aℯ𝒻𝒻} \: V_{Aℯ𝒻𝒻} \; \cos(ϕ_A )  où le déphasage entre fAf_A et vAv_A est  ϕAϕϕ_A≠ϕ .

• On peut aussi considérer que cette puissance est celle fournie par le ressort à MM ,  ou dissipée par les frottements (équivalant à l'effet Joule en électricité) :
v(t)=z˙=ωZmsin(ωt+ϕ)v(t)=\dot{z}=-ω \:Z_m \; \sin(ω \,t+ϕ)  ;
f(t)=λz˙=λωZmsin(ωt+ϕ)f(t)=-λ \:\dot{z}=λ \:ω \:Z_m \; \sin(ω \,t+ϕ)  ;
P(t)=f.v=λωZm2sin2(ωt+ϕ)P(t)=-f.v=-λ \:ω \:Z_m^{\:2} \; \sin^2(ω \,t+ϕ)  ;
P=𝒫(t)=12λω2Zm2=kαZAm2ω02ω2(ω02ω2)2+(2αω)2=P0x2(1x2)2+4α2x2\displaystyle P=〈 𝒫(t) 〉=-\frac{1}{2} λ \:ω^2 \: Z_m^{\:2}=\frac{k \:α \:Z_{Am}^{\:2} \: ω_0^{\:2} \: ω^2}{\left(ω_0^{\:2}-ω^2 \right)^2+(2 \,α \:ω)^2}=\frac{P_0 \: x^2}{\left(1-x^2 \right)^2+4 \:α^2 \: x^2} .

Ceci correspond ici encore à une courbe de résonance :

dynOscForc_Im/puissance.jpg