Vibrations Mécaniques des Systèmes Forcés à 1 Degré de Liberté

Exercice n°2 (5 points) : Pendule forcé.

Image centrée

Un pendule constitué d'une tige de masse négligeable, portant à ses deux extrémités deux masses `m_1` et `m_2,` peut tourner dans le plan vertical autour d'un axe horizontal `(Δ)` fixe passant par le point `O` de la tige. Ces deux masses, `m_1` et `m_2,` sont distantes de `O` de ##{ℓ}_1## et ##{ℓ}_2## respectivement. La masse `m_1` est reliée au bâti `B_1` par un ressort de constante de raideur `K` et la masse `m_2` à un bâti `B_2` par un amortisseur de coefficient de frottement visqueux `\alpha ` (voir la figure ci-dessus).

A l'équilibre, la tige est verticale et le ressort non déformé. La position du pendule est repérée par l'angle θ que fait la tige avec la verticale. Le bâti `B_2` est animé d'un mouvement horizontal harmonique `s(t)=S_0 sin(ωt).`

Données :
`ℓ_2=2ℓ_1=2ℓ  ,  m_2=2m_1=2m  , m=1 kg,K=3500 N//m, ℓ= 50 cm ` et ` g=10 m//s^2.`

1/Etablir l'équation du mouvement du système dans le cas des mouvements de faible amplitude.

Réponse

Image centrée

Les coordonnées cartésiennes et la vitesse de la masse `m_1` :
## \begin{cases} x_1=-{ℓ}_1sin(θ) \\ \\ y_1={ℓ}_1cos(θ) \end{cases} \implies## ## \begin{cases} \dot{x}_1=-{ℓ}_1 \dotθ cos(θ) \\ \\ \dot{y}_2=-{ℓ}_1\dotθ sin(θ) \end{cases} \implies## ## v_1^2=\dot{x}_1+\dot{y}_1={ℓ}_1^2 \dot{θ}^2 ##

et de la masse `m_2`:
## \begin{cases} x_2={ℓ}_2sin(θ) \\ \\ y_2=-{ℓ}_2cos(θ) \end{cases} \implies## ## \begin{cases} \dot{x}_1={ℓ}_2 \dotθ cos(θ) \\ \\ \dot{y}_2={ℓ}_2\dotθ sin(θ) \end{cases} \implies## ## v_2^2=\dot{x}_2+\dot{y}_2={ℓ}_2^2 \dot{θ}^2 ##

Energie cinétique :

## E_c=E_{c1}+E_{c2}\implies## ##=\dfrac 12 m_1 v_1^2+\dfrac 12 m_2 v_2^2\implies## ##=\dfrac 12 m_1{ℓ}_1^2 \dot{θ}^2+\dfrac 12 m_2{ℓ}_2^2 \dot{θ}^2\implies## ## E_c=\dfrac 12 (m_1{ℓ}_1^2 +m_2{ℓ}_2^2 )\dot{θ}^2 ##


Energie potentielle :
Comme θ est petit, le déplacement du ressort est supposé pratiquement horizontal.
A l'équilibre, le ressort n'est pas déformé, le ##Δ {ℓ}_{equ}## est nul.
## E_P=\dfrac 12 K(-x_2+s )^2+m_1 gy_1-m_2 gy_2 \implies## ## E_P= \dfrac 12 K(- {ℓ}_2 sin(θ)+s)^2+m_1 g{ℓ}_1 cos(θ)-m_2 g{ℓ}_2 cos(θ) ##

On peut encore faire les approximations suivantes :
`sin(θ)=θ` et ##cos(θ)=1-\dfrac {θ^2}{2}##
## E_P= \dfrac 12 K (-{ℓ}_2 θ+s)^2- \dfrac 12 m_1 g{ℓ}_1 θ^2+\dfrac 12 m_2 g{ℓ}_2 θ^2+ constante\implies ## ## E_P= \dfrac 12 ( K (-{ℓ}_2 θ+s)^2- m_1 g{ℓ}_1 θ^+m_2 g{ℓ}_2 θ^2)+ constante ##


Lagrangien :
## ℒ=E_c-E_P \implies## ## ℒ=\dfrac 12 (m_1 {ℓ}_1^2+m_2 {ℓ}_2^2 ) θ ̇^2- \dfrac 12 (K (-{ℓ}_2 θ+s)^2+m_2 g{ℓ}_2-m_1 g {ℓ}_1 ) θ^2 ##


Fonction de dissipation :
## 𝒟=\dfrac 12 \alpha v_1^2 \implies## ##𝒟=\dfrac 12 \alpha ( {ℓ}_1 \dotθ )^2\implies## ## 𝒟=\dfrac 12 \alpha {ℓ}_1^2 \dotθ^2 ##


Equation du mouvement :
## \dfrac {d}{dt} \left(\dfrac {∂ ℒ}{∂\dotθ}\right) -\dfrac {∂ℒ} {∂θ}=-\dfrac {∂{𝒟}} {∂\dotθ}\implies## ## (m_1{ℓ}_1^2 +m_2{ℓ}_2^2 )\ddot{θ}-(K {ℓ}_2({ℓ}_2 θ+s)+m_2 g{ℓ}_2-m_1 g {ℓ}_1 ) θ=-\alpha {ℓ}_1^2 \dotθ\implies## ## (m_1{ℓ}_1^2 +m_2{ℓ}_2^2 )\ddot{θ}+\alpha {ℓ}_1^2 \dotθ+(K {ℓ}_2^2+m_2 g{ℓ}_2-m_1 g {ℓ}_1 ) θ=K {ℓ}_2s ##


On divise l'équation par le coefficient de ##\ddot{θ}##:
## \ddot{θ}+\dfrac {\alpha {ℓ}_1^2}{m_1{ℓ}_1^2 +m_2{ℓ}_2^2 } \dotθ+\dfrac {K {ℓ}_2^2+m_2 g{ℓ}_2-m_1 g {ℓ}_1 }{m_1{ℓ}_1^2 +m_2{ℓ}_2^2 } θ= \dfrac {K {ℓ}_2}{m_1{ℓ}_1^2 +m_2{ℓ}_2^2 }s\implies## ## \ddot{θ}+2\dfrac {\alpha {ℓ}_1^2}{2(m_1{ℓ}_1^2 +m_2{ℓ}_2^2) } \dotθ+\dfrac {K {ℓ}_2^2+m_2 g{ℓ}_2-m_1 g {ℓ}_1 }{m_1{ℓ}_1^2 +m_2{ℓ}_2^2 } θ= \dfrac {K {ℓ}_2}{m_1{ℓ}_1^2 +m_2{ℓ}_2^2 }s ##


Comme `ℓ_2=2ℓ_1=2ℓ  ,  m_2=2m_1=2m  `
## \ddot{θ}+2\dfrac {\alpha {ℓ}^2}{2(m{ℓ}^2 +8m{ℓ}^2) } \dotθ+\dfrac {4K {ℓ}^2+4m g{ℓ}-m g {ℓ} }{m{ℓ}^2 +8m{ℓ}^2 } θ= \dfrac {2K {ℓ}}{m{ℓ}^2 +8m{ℓ}^2 }s\implies## ## \ddot{θ}+2\dfrac {\alpha {ℓ}^2}{2(m{ℓ}_1^2 +8m{ℓ}^2) } \dotθ+\dfrac {4K {ℓ}^2+4m g{ℓ}-m g {ℓ} }{m{ℓ}^2 +8m{ℓ}^2 } θ= \dfrac {2K {ℓ}}{m{ℓ}^2 +8m{ℓ}^2 }s\implies## ## \ddot{θ}+2\dfrac {\alpha }{18m } \dotθ+ \left(\dfrac {4}{9}\dfrac {K}{m}+\dfrac {1}{3}\dfrac {g}{{ℓ} }\right) θ= \dfrac {2}{9}\dfrac {K}{m{ℓ} }s\implies## ## \ddot{θ}+2\times\dfrac {13,4}{18\times 1} \dotθ+ \left(\dfrac {4}{9}\times\dfrac {3500}{1}+\dfrac {1}{3}\times\dfrac {10}{{0,5} }\right) θ= \dfrac {2}{9}\dfrac {3500}{1\times 0,5}s \implies##
## \ddot{θ}+2\times0,744\dotθ+39,5^2θ=1556S_0 cos(ωt) ##


Pulsation propre :
##ω_0=39,5 rd/s##

et
Coefficient d'amortissement :
##δ=0,744kg/s ##

2/ Donner l'équation horaire `θ_p=g(t)` en régime permanent.

Réponse

## θ_p (t)=\dfrac{\dfrac{2}{9}\dfrac{K}{m{ℓ} }S_0 }{\sqrt{(ω_0^2-ω^2 )^2+4δ^2 ω^2 }}cos(ωt-ϕ)## où ## ϕ=arccos\left(\dfrac{ω_0^2-ω^2}{\sqrt{(ω_0^2-ω^2 )^2+4δ^2 ω^2 }}\right) (et sin(ϕ)>0)\implies##

## θ_p (t)=\dfrac{1556 S_0 }{\sqrt{(1560-ω^2 )^2+2,21ω^2 }}cos(ωt-ϕ) ## où ## ϕ=arccos\left(\dfrac{1560-ω^2}{\sqrt{(1560-ω^2 )^2+2,21ω^2 }}\right) ## ##(et sin(ϕ)>0) ##

3/Calculer l'amplitude de la déviation maximale du pendule en fonction de `S_0.`

Réponse

L'amplitude de la déviation maximale est obtenue à la pulsation de résonance :

## ω_R=\sqrt{ω_0^2-2δ^2} \implies## ##ω_R=\sqrt{39,5^2-2\times0,744^2 }\implies## ##ω_R≅39,5 rd/s ##

## θ_{max}=\dfrac{1556 S_0 }{\sqrt{(39,5^2-ω_R^2 )^2+2,21ω_R^2 }}\implies## ## θ_{max}=\dfrac{1556 S_0 }{\sqrt{(39,5^2-39,5^2 )^2+2,21\times 39,5^2 }}\implies##
##θ_{max}=26,5 S_0##