Salta al contenuto
Note per Studenti Esercizio 174 · oscillatore forzato con smorzamento, transitorio e regime permanente

Esercizio 174oscillatore forzato con smorzamento, transitorio e regime permanente

In questa pagina 5
In questa pagina 4

Testo. Una massa m=1m = 1 kg, collegata a una molla di costante k=5k = 5 N/m e a uno smorzatore di coefficiente c=2c = 2 N⋅\cdots/m, è spinta da una forza esterna F(t)=10cos⁡(2t)F(t) = 10\cos(2t) N. La posizione x(t)x(t) soddisfa x′′+2x′+5x=10cos⁡(2t)x'' + 2x' + 5x = 10\cos(2t)

  1. Determinare l'integrale generale.
  2. Determinare ampiezza e sfasamento del moto a regime.
  3. Determinare la soluzione con la massa inizialmente ferma nell'equilibrio, x(0)=0x(0) = 0 e x′(0)=0x'(0) = 0.
  4. Dopo quanto tempo la parte transitoria è minore dell'1%1\% dell'ampiezza a regime?

Teoria: Oscillatore massa-molla-smorzatore e transitoriMassa m, molla k e smorzatore viscoso c con forza esterna F(t): m x'' + c x' + k x = F(t), cioè x'' + (c/m) x' + (k/m) x = F/m. Grandezze: ω0 = √(k/m), smorzamento critico c_cr = 2√(km), ζ = c/c_cr. Libero: c > c_cr due esponenziali negativi (sovrasmorzato, nessuna oscillazione), c = c_cr (c1 + c2 t) e^(−ω0 t), c < c_cr e^(−ct/2m)(c1 cos ωd t + c2 sin ωd t) con ωd = ω0√(1−ζ²); con c > 0 ogni moto libero tende a 0 (transitorio). Forzato F0 cos Ωt con c > 0: x = transitorio + regime permanente X cos(Ωt − φ), X = F0/√((k−mΩ²)² + c²Ω²), tan φ = cΩ/(k−mΩ²); massimo di X per Ω² = k/m − c²/(2m²). Primo ordine m v' + c v = F: v = F/c + (v0 − F/c) e^(−ct/m), costante di tempo τ = m/c.Oscillatore massa-molla-smorzatore e transitori →, Equazioni lineari a coefficienti costanti - equazione completaPer y'' + a1 y' + a0 y = b(t) a coefficienti costanti l'integrale generale è y0 + yp: y0 dal polinomio caratteristico, yp con il metodo di somiglianza quando b(t) = Qm(t) e^(αt) cos(ωt) (o sin). Si pone λ̄ = α + iω: se p(λ̄) ≠ 0, yp = e^(αt)[Pm(t) cos ωt + Rm(t) sin ωt] con polinomi di grado m da determinare; se λ̄ è radice semplice si moltiplica per t, se è doppia per t² (risonanza). Se b = b1 + b2 si sommano le soluzioni particolari (principio di sovrapposizione).Equazioni lineari a coefficienti costanti - equazione completa →.


1. Integrale generale

Omogenea. p(λ)=λ2+2λ+5=0p(\lambda) = \lambda^2 + 2\lambda + 5 = 0, Δ=4−20=−16<0\Delta = 4 - 20 = -16 < 0, λ=−1±2i\lambda = -1 \pm 2i, quindi α=−1\alpha = -1, ω=2\omega = 2: x0(t)=e−t(c1cos⁡2t+c2sin⁡2t)∀c1,c2∈Rx_0(t) = e^{-t}\big(c_1\cos 2t + c_2\sin 2t\big) \qquad \forall c_1, c_2 \in \mathbb{R} (è il sottosmorzato: c=2<ccr=25≈4,47c = 2 < c_{cr} = 2\sqrt{5} \approx 4{,}47).

Soluzione particolare. b(t)=10cos⁡2tb(t) = 10\cos 2t, quindi λˉ=0+2i=2i\bar\lambda = 0 + 2i = 2i. Non è radice di pp (le radici hanno parte reale −1-1, quella di 2i2i è 00): niente risonanza, e xp(t)=Acos⁡2t+Bsin⁡2tx_p(t) = A\cos 2t + B\sin 2t (anche il seno, perché derivando il coseno compare il seno). Le derivate sono xp′=−2Asin⁡2t+2Bcos⁡2t,xp′′=−4Acos⁡2t−4Bsin⁡2tx_p' = -2A\sin 2t + 2B\cos 2t, \qquad x_p'' = -4A\cos 2t - 4B\sin 2t Sostituendo in xp′′+2xp′+5xp=10cos⁡2tx_p'' + 2x_p' + 5x_p = 10\cos 2t e raccogliendo seno e coseno: cos⁡2t:  −4A+4B+5A=A+4B,sin⁡2t:  −4B−4A+5B=B−4A\cos 2t:\ \ -4A + 4B + 5A = A + 4B, \qquad \sin 2t:\ \ -4B - 4A + 5B = B - 4A Il termine 2xp′2x_p' contribuisce con 2⋅2Bcos⁡2t=4Bcos⁡2t2\cdot 2B\cos 2t = 4B\cos 2t e con 2⋅(−2A)sin⁡2t=−4Asin⁡2t2\cdot(-2A)\sin 2t = -4A\sin 2t. L'uguaglianza vale per ogni tt se {A+4B=10−4A+B=0\begin{cases} A + 4B = 10 \\ -4A + B = 0 \end{cases} Dalla seconda B=4AB = 4A; nella prima A+16A=17A=10A + 16A = 17A = 10, quindi A=1017,B=4017,xp(t)=1017cos⁡2t+4017sin⁡2tA = \frac{10}{17}, \qquad B = \frac{40}{17}, \qquad x_p(t) = \frac{10}{17}\cos 2t + \frac{40}{17}\sin 2t Integrale generale: x(t)=e−t(c1cos⁡2t+c2sin⁡2t)+1017cos⁡2t+4017sin⁡2t∀t∈Rx(t) = e^{-t}\big(c_1\cos 2t + c_2\sin 2t\big) + \frac{10}{17}\cos 2t + \frac{40}{17}\sin 2t \qquad \forall t \in \mathbb{R}

2. Regime permanente

Si scrive xp=Xcos⁡(2t−φ)x_p = X\cos(2t - \varphi) con X=A2+B2=100+1600289=170017=1017≈2,425 m,tan⁡φ=BA=4 ⇒ φ=arctan⁡4≈1,326 radX = \sqrt{A^2 + B^2} = \sqrt{\frac{100 + 1600}{289}} = \frac{\sqrt{1700}}{17} = \frac{10}{\sqrt{17}} \approx 2{,}425 \ \text{m}, \qquad \tan\varphi = \frac BA = 4 \ \Rightarrow\ \varphi = \arctan 4 \approx 1{,}326 \ \text{rad} (φ\varphi è circa 7676 gradi, in ritardo di φΩ=1,3262≈0,66\frac{\varphi}{\Omega} = \frac{1{,}326}{2} \approx 0{,}66 s rispetto alla forza.) La stessa ampiezza si ottiene dalla formula della nota: D=(k−mΩ2)2+(cΩ)2=(5−4)2+42=17D = (k - m\Omega^2)^2 + (c\Omega)^2 = (5 - 4)^2 + 4^2 = 17 e X=F0D=1017X = \frac{F_0}{\sqrt D} = \frac{10}{\sqrt{17}} ✓.

3. Condizioni iniziali

Serve x′x': x′(t)=e−t[(−c1+2c2)cos⁡2t+(−2c1−c2)sin⁡2t]−2017sin⁡2t+8017cos⁡2tx'(t) = e^{-t}\big[(-c_1 + 2c_2)\cos 2t + (-2c_1 - c_2)\sin 2t\big] - \frac{20}{17}\sin 2t + \frac{80}{17}\cos 2t In t=0t = 0: {x(0)=c1+1017=0x′(0)=−c1+2c2+8017=0\begin{cases} x(0) = c_1 + \dfrac{10}{17} = 0 \\[2mm] x'(0) = -c_1 + 2c_2 + \dfrac{80}{17} = 0 \end{cases} Dalla prima c1=−1017c_1 = -\frac{10}{17}. Nella seconda: 1017+2c2+8017=0\frac{10}{17} + 2c_2 + \frac{80}{17} = 0, cioè 2c2=−90172c_2 = -\frac{90}{17} e c2=−4517c_2 = -\frac{45}{17}. Quindi x(t)=117[e−t(−10cos⁡2t−45sin⁡2t)+10cos⁡2t+40sin⁡2t]x(t) = \frac{1}{17}\Big[e^{-t}\big(-10\cos 2t - 45\sin 2t\big) + 10\cos 2t + 40\sin 2t\Big] Verifica (Python, dsolve): coincide con la soluzione trovata.

4. Quando il transitorio diventa trascurabile

Il transitorio è x0(t)=e−t(c1cos⁡2t+c2sin⁡2t)x_0(t) = e^{-t}\big(c_1\cos 2t + c_2\sin 2t\big) con c1=−1017c_1 = -\frac{10}{17}, c2=−4517c_2 = -\frac{45}{17}. Per la disuguaglianza ∣c1cos⁡2t+c2sin⁡2t∣≤c12+c22|c_1\cos 2t + c_2\sin 2t| \le \sqrt{c_1^2 + c_2^2} (Cauchy-Schwarz, oppure scrivendo la combinazione come Ccos⁡(2t−ψ)C\cos(2t - \psi)): ∣x0(t)∣≤C e−t,C=100+202517=212517≈2,712|x_0(t)| \le C\,e^{-t}, \qquad C = \frac{\sqrt{100 + 2025}}{17} = \frac{\sqrt{2125}}{17} \approx 2{,}712 Si vuole C e−t≤0,01 X≈0,02425C\,e^{-t} \le 0{,}01\,X \approx 0{,}02425: e−t≤0,024252,712≈0,00894  ⟺  t≥ln⁡2,7120,02425≈ln⁡111,8≈4,72 se^{-t} \le \frac{0{,}02425}{2{,}712} \approx 0{,}00894 \iff t \ge \ln\frac{2{,}712}{0{,}02425} \approx \ln 111{,}8 \approx 4{,}72 \ \text{s} Dopo circa 4,74{,}7 s il moto coincide, a meno dell'1%1\% dell'ampiezza, con l'oscillazione a regime. Il tempo è legato al decadimento e−te^{-t}: la costante di tempo è 1∣α∣=2mc=1\frac{1}{|\alpha|} = \frac{2m}{c} = 1 s, e servono circa 4,74{,}7 costanti di tempo.

Grafico interattivo: Esercizio 174: la soluzione x(t) (partenza da fermo) si avvicina al regime permanente x_p(t)


Errori comuni

  • Cercare xp=Acos⁡2tx_p = A\cos 2t senza il seno: con c≠0c \ne 0 i due coefficienti sono entrambi diversi da 00 (B=4AB = 4A).
  • Imporre le condizioni iniziali alla sola x0x_0 (dimenticando xp(0)=1017x_p(0) = \frac{10}{17} e xp′(0)=8017x_p'(0) = \frac{80}{17}): le costanti c1,c2c_1, c_2 vanno trovate sull'integrale generale completo.
  • Pensare che c1c_1 e c2c_2 cambino il regime permanente: xpx_p non dipende dalle condizioni iniziali.
  • Calcolare l'ampiezza come A+BA + B anziché A2+B2\sqrt{A^2 + B^2}.

Versione ripasso

Testo. x′′+2x′+5x=10cos⁡2tx'' + 2x' + 5x = 10\cos 2t. Integrale generale; ampiezza e sfasamento a regime; soluzione con x(0)=x′(0)=0x(0) = x'(0) = 0; tempo perché il transitorio sia sotto l'1%1\% dell'ampiezza di regime.

TeoriaMassa m, molla k e smorzatore viscoso c con forza esterna F(t): m x'' + c x' + k x = F(t), cioè x'' + (c/m) x' + (k/m) x = F/m. Grandezze: ω0 = √(k/m), smorzamento critico c_cr = 2√(km), ζ = c/c_cr. Libero: c > c_cr due esponenziali negativi (sovrasmorzato, nessuna oscillazione), c = c_cr (c1 + c2 t) e^(−ω0 t), c < c_cr e^(−ct/2m)(c1 cos ωd t + c2 sin ωd t) con ωd = ω0√(1−ζ²); con c > 0 ogni moto libero tende a 0 (transitorio). Forzato F0 cos Ωt con c > 0: x = transitorio + regime permanente X cos(Ωt − φ), X = F0/√((k−mΩ²)² + c²Ω²), tan φ = cΩ/(k−mΩ²); massimo di X per Ω² = k/m − c²/(2m²). Primo ordine m v' + c v = F: v = F/c + (v0 − F/c) e^(−ct/m), costante di tempo τ = m/c.Oscillatore massa-molla-smorzatore e transitori →

1. Integrale generale

  • Omogenea: λ=−1±2i\lambda = -1 \pm 2i, x0=e−t(c1cos⁡2t+c2sin⁡2t)x_0 = e^{-t}(c_1\cos 2t + c_2\sin 2t).
  • λˉ=2i\bar\lambda = 2i non è radice: xp=Acos⁡2t+Bsin⁡2tx_p = A\cos 2t + B\sin 2t.
  • Sostituzione: A+4B=10A + 4B = 10, −4A+B=0-4A + B = 0 ⇒\Rightarrow A=1017A = \frac{10}{17}, B=4017B = \frac{40}{17}.

x(t)=e−t(c1cos⁡2t+c2sin⁡2t)+1017cos⁡2t+4017sin⁡2tx(t) = e^{-t}(c_1\cos 2t + c_2\sin 2t) + \tfrac{10}{17}\cos 2t + \tfrac{40}{17}\sin 2t

2. Regime permanente

X=A2+B2=1017≈2,425X = \sqrt{A^2 + B^2} = \frac{10}{\sqrt{17}} \approx 2{,}425 m, tan⁡φ=4\tan\varphi = 4, φ≈1,326\varphi \approx 1{,}326 rad.

3. Condizioni iniziali

c1+1017=0c_1 + \frac{10}{17} = 0, −c1+2c2+8017=0-c_1 + 2c_2 + \frac{80}{17} = 0 ⇒\Rightarrow c1=−1017c_1 = -\frac{10}{17}, c2=−4517c_2 = -\frac{45}{17}.

4. Transitorio

∣x0∣≤C e−t|x_0| \le C\,e^{-t} con C=212517≈2,712C = \frac{\sqrt{2125}}{17} \approx 2{,}712. C e−t≤0,01X  ⟺  t≥ln⁡111,8≈4,72C\,e^{-t} \le 0{,}01X \iff t \ge \ln 111{,}8 \approx 4{,}72 s (costante di tempo 2mc=1\frac{2m}{c} = 1 s).

Errori tipici: omettere il seno in xpx_p; imporre le condizioni iniziali solo a x0x_0; ampiezza come A+BA + B.

Teoria collegata