Math ClubMath Club
v1 · padrão canônico

Lição 96 — Vibrações mecânicas: massa-mola-amortecedor

m x'' + c x' + k x = F(t). Frequência natural, amortecimento, ressonância. Subamortecido, crítico, superamortecido.

Used in: Spécialité Maths (France, Terminale) · Leistungskurs alemão Klasse 12 · University Physics (global)

mx¨+cx˙+kx=F(t)x¨+2ζω0x˙+ω02x=F(t)mm\ddot x + c\dot x + kx = F(t) \quad\Longleftrightarrow\quad \ddot x + 2\zeta\omega_0\dot x + \omega_0^2 x = \frac{F(t)}{m}

Oscilador massa-mola-amortecedor. Frequência natural ω0=k/m\omega_0 = \sqrt{k/m} [rad/s]; fator de amortecimento ζ=c/(2mk)\zeta = c/(2\sqrt{mk}). Quatro regimes: não-amortecido (ζ=0\zeta=0), subamortecido (ζ<1\zeta<1), criticamente amortecido (ζ=1\zeta=1), superamortecido (ζ>1\zeta>1).

Choose your door

Rigorous notation, full derivation, hypotheses

Oscilador completo — quatro regimes

Equação de movimento

Equacao caracteristica e regimes

λ2+2ζω0λ+ω02=0\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0. Discriminante Δ=4ω02(ζ21)\Delta = 4\omega_0^2(\zeta^2 - 1).

"The most important case is b24km<0b^2 - 4km < 0, which occurs when the damping is small... In this case the solution oscillates with exponentially decaying amplitude." — Lebl, Notes on Diffy Qs, §2.4

Resposta forçada harmonica

Para F(t)=F0cosωtF(t) = F_0\cos\omega t: solução particular (regime permanente)

xp(t)=F0/m(ω02ω2)2+4ζ2ω02ω2cos(ωtϕ)x_p(t) = \frac{F_0/m}{\sqrt{(\omega_0^2-\omega^2)^2 + 4\zeta^2\omega_0^2\omega^2}}\cos(\omega t - \phi)

onde tanϕ=2ζω0ω/(ω02ω2)\tan\phi = 2\zeta\omega_0\omega/(\omega_0^2 - \omega^2).

Ressonancia

Diagrama qualitativo dos regimes

txsubamortecido (zeta menor que 1)superamortecidocritico

Resposta livre (F=0F = 0, x(0)=x0>0x(0) = x_0 > 0, x˙(0)=0\dot x(0) = 0): subamortecido oscila enquanto decai; crítico e super convergem monotonicamente.

Exemplos resolvidos

Exercise list

24 exercises · 6 with worked solution (25%)

Application 12Understanding 2Modeling 5Challenge 3Proof 2
  1. Ex. 96.1Application

    Mola sem amortecimento: m=1kgm = 1\,\text{kg}, k=4N/mk = 4\,\text{N/m}. Calcule ω0\omega_0, período e escreva a solução geral.

    Show solution
    ω0=k/m=4=2\omega_0 = \sqrt{k/m} = \sqrt{4} = 2 rad/s. Período: T=2π/2=πT = 2\pi/2 = \pi s. Solução: x(t)=C1cos2t+C2sin2tx(t) = C_1\cos 2t + C_2\sin 2t.
  2. Ex. 96.2Application

    m=1kgm = 1\,\text{kg}, k=16N/mk = 16\,\text{N/m}. Classifique o regime para: (a) c=2c = 2, (b) c=8c = 8, (c) c=10c = 10.

    Show solution
    ζ=c/(2mk)\zeta = c/(2\sqrt{mk}). Com m=1m = 1, k=16k = 16: cc=216=8c_c = 2\sqrt{16} = 8. (a) c=2c = 2: ζ=0,25\zeta = 0{,}25 — subamortecido. (b) c=8c = 8: ζ=1\zeta = 1 — crítico. (c) c=10c = 10: ζ=1,25\zeta = 1{,}25 — superamortecido.
    Show step-by-step (with the why)
    1. ω0=16/1=4\omega_0 = \sqrt{16/1} = 4 rad/s.
    2. Amortecimento crítico: cc=2mω0=8c_c = 2m\omega_0 = 8 N.s/m.
    3. Calcule ζ=c/cc\zeta = c/c_c para cada caso.
    4. Compare com 1: subamortecido se <1< 1, crítico se =1= 1, superamortecido se >1> 1.
    Macete: o amortecimento crítico cc=2mkc_c = 2\sqrt{mk} é o valor que zera o discriminante da equação característica.
  3. Ex. 96.3ApplicationAnswer key

    x¨+6x˙+9x=0\ddot x + 6\dot x + 9x = 0, x(0)=1x(0) = 1, x˙(0)=1\dot x(0) = -1. Classifique e resolva.

    Show solution
    ω0=3\omega_0 = 3, ζ=6/(23)=1\zeta = 6/(2 \cdot 3) = 1. Criticamente amortecido: x=(A+Bt)e3tx = (A+Bt)e^{-3t}. CI: A=1A = 1, B3=1B=2B - 3 = -1 \Rightarrow B = 2. x=(1+2t)e3tx = (1+2t)e^{-3t}.
  4. Ex. 96.4Application

    x¨+4x˙+8x=0\ddot x + 4\dot x + 8x = 0, x(0)=3x(0) = 3, x˙(0)=0\dot x(0) = 0. Resolva e calcule ωd\omega_d.

    Show solution
    λ2+4λ+8=0λ=2±2i\lambda^2 + 4\lambda + 8 = 0 \Rightarrow \lambda = -2 \pm 2i. Subamortecido: x=e2t(Acos2t+Bsin2t)x = e^{-2t}(A\cos 2t + B\sin 2t). CI: A=3A = 3, 2A+2B=0B=3-2A + 2B = 0 \Rightarrow B = 3. x=3e2t(cos2t+sin2t)x = 3e^{-2t}(\cos 2t + \sin 2t).
  5. Ex. 96.5Application

    x¨+6x˙+5x=0\ddot x + 6\dot x + 5x = 0, x(0)=2x(0) = 2, x˙(0)=4\dot x(0) = -4. Superamortecido — resolva.

    Show solution
    λ2+6λ+5=(λ+1)(λ+5)=0\lambda^2 + 6\lambda + 5 = (\lambda+1)(\lambda+5) = 0. Superamortecido: x=Aet+Be5tx = Ae^{-t} + Be^{-5t}. CI: A+B=2A+B = 2, A5B=4B=0,5-A-5B = -4 \Rightarrow B = 0{,}5, A=1,5A = 1{,}5. x=1,5et+0,5e5tx = 1{,}5e^{-t} + 0{,}5e^{-5t}.
  6. Ex. 96.6Application

    x¨+4x=cost\ddot x + 4x = \cos t (sem amortecimento). Calcule a amplitude do regime permanente.

    Show solution
    Amplitude do regime permanente: X=F0/(kmω2)2+c2ω2|X| = F_0/\sqrt{(k - m\omega^2)^2 + c^2\omega^2}. Com m=1m = 1, k=4k = 4, c=0c = 0, ω=1\omega = 1, F0=1F_0 = 1: X=1/(41)2=1/3|X| = 1/\sqrt{(4-1)^2} = 1/3.
  7. Ex. 96.7Application

    Resolva x¨+4x=2cos3t\ddot x + 4x = 2\cos 3t.

    Show solution
    Solução particular: xp=Acos3t+Bsin3tx_p = A\cos 3t + B\sin 3t. Substituindo em x¨+4x=2cos3t\ddot x + 4x = 2\cos 3t: 9A+4A=2A=2/5-9A + 4A = 2 \Rightarrow A = -2/5, B=0B = 0. x=C1cos2t+C2sin2t25cos3tx = C_1\cos 2t + C_2\sin 2t - \tfrac{2}{5}\cos 3t.
  8. Ex. 96.8ApplicationAnswer key

    Ressonância pura: resolva x¨+4x=cos2t\ddot x + 4x = \cos 2t. O que acontece com a amplitude?

    Show solution
    Ressonância: ω=ω0=2\omega = \omega_0 = 2. Solução particular: xp=t(Acos2t+Bsin2t)x_p = t(A\cos 2t + B\sin 2t). Determinando: B=1/4B = 1/4, A=0A = 0. xp=t4sin2tx_p = \tfrac{t}{4}\sin 2t. Amplitude cresce sem limite.
    Show step-by-step (with the why)
    1. Forma: x¨+4x=cos2t\ddot x + 4x = \cos 2t. A frequência ω=2=ω0\omega = 2 = \omega_0.
    2. Ressonância: tentativa xp=t(Acos2t+Bsin2t)x_p = t(A\cos 2t + B\sin 2t).
    3. xp=4Atcos2t4Btsin2t4Asin2t+4Bcos2tx_p'' = -4A t\cos 2t - 4B t\sin 2t - 4A\sin 2t + 4B\cos 2t.
    4. xp+4xp=4Asin2t+4Bcos2t=cos2tx_p'' + 4x_p = -4A\sin 2t + 4B\cos 2t = \cos 2t. Logo A=0A = 0, B=1/4B = 1/4.
    5. xp=t4sin2tx_p = \frac{t}{4}\sin 2t. O fator tt faz a amplitude crescer linearmente.
    Curiosidade: a Ponte Tacoma Narrows colapsou em 1940 por fenômeno análogo: o vento entrou em ressonância com a frequência natural da ponte.
  9. Ex. 96.9Application

    Resolva x¨+2x˙+5x=5cost\ddot x + 2\dot x + 5x = 5\cos t.

    Show solution
    xh=et(C1cos2t+C2sin2t)x_h = e^{-t}(C_1\cos 2t + C_2\sin 2t). Particular para q=5costq = 5\cos t: xp=Acost+Bsintx_p = A\cos t + B\sin t. Substituindo: (4A+2B)cost+(4B2A)sint=5cost(4A+2B)\cos t + (4B-2A)\sin t = 5\cos t. Resolvendo: A=1A = 1, B=1/2B = 1/2. x=et(C1cos2t+C2sin2t)+cost+12sintx = e^{-t}(C_1\cos 2t+C_2\sin 2t) + \cos t + \tfrac{1}{2}\sin t.
  10. Ex. 96.10ApplicationAnswer key

    Em ensaio de vibração, dois picos consecutivos medem x1=1,20x_1 = 1{,}20 m e x2=0,89x_2 = 0{,}89 m. Calcule o decremento logarítmico e o fator de amortecimento ζ\zeta.

    Show solution
    Decremento logarítmico: mede-se δ=ln(xn/xn+1)=2πζ/1ζ2\delta = \ln(x_n/x_{n+1}) = 2\pi\zeta/\sqrt{1-\zeta^2}. Com δ=0,3\delta = 0{,}3: ζ=δ/4π2+δ20,3/6,280,0477\zeta = \delta/\sqrt{4\pi^2 + \delta^2} \approx 0{,}3/6{,}28 \approx 0{,}0477. Subamortecido leve.
  11. Ex. 96.11Modeling

    Suspensão automotiva: m=100kgm = 100\,\text{kg}, k=50000N/mk = 50000\,\text{N/m}, c=2000N.s/mc = 2000\,\text{N.s/m}. Calcule ω0\omega_0, ccc_c e ζ\zeta. Está sub ou superamortecida?

    Show solution
    ω0=50000/100=50022,4\omega_0 = \sqrt{50000/100} = \sqrt{500} \approx 22{,}4 rad/s. cc=2100×50000=25×1064472c_c = 2\sqrt{100 \times 50000} = 2\sqrt{5 \times 10^6} \approx 4472 N.s/m. ζ=2000/44720,447\zeta = 2000/4472 \approx 0{,}447 — subamortecido.
  12. Ex. 96.12Modeling

    m=1kgm = 1\,\text{kg}, k=100N/mk = 100\,\text{N/m}, c=2c = 2 N.s/m. Calcule ω0\omega_0, ζ\zeta, frequência de pico e fator de amplificação.

    Show solution
    Frequência de ressonância: ωr=ω012ζ2\omega_r = \omega_0\sqrt{1-2\zeta^2}. Fator de amplificação máximo: Hmax=1/(2ζ1ζ2)|H|_{\max} = 1/(2\zeta\sqrt{1-\zeta^2}). Com ω0=10\omega_0 = 10 rad/s, ζ=0,1\zeta = 0{,}1: ωr100,989,9\omega_r \approx 10\sqrt{0{,}98} \approx 9{,}9 rad/s. Hmax1/(0,2×0,995)5,03|H|_{\max} \approx 1/(0{,}2 \times 0{,}995) \approx 5{,}03.
    Show step-by-step (with the why)
    1. ω0=k/m\omega_0 = \sqrt{k/m}. Para m=1m = 1, k=100k = 100: ω0=10\omega_0 = 10.
    2. ζ=c/(2mk)=2/(2×10)=0,1\zeta = c/(2\sqrt{mk}) = 2/(2\times 10) = 0{,}1.
    3. Frequência de pico: ωr=ω012ζ2=100,989,9\omega_r = \omega_0\sqrt{1-2\zeta^2} = 10\sqrt{0{,}98} \approx 9{,}9 rad/s.
    4. Fator Q = 1/(2ζ)=51/(2\zeta) = 5. Amplificação ≈ 5×.
    Curiosidade: um fator de qualidade Q=5Q = 5 significa que em ressonância a amplitude é 5 vezes a estática. Estruturas com Q>50Q > 50 são perigosas em sismos.
  13. Ex. 96.13Modeling

    Pêndulo de comprimento L=1mL = 1\,\text{m}. Calcule ω0\omega_0 e o período TT. (Use g=9,8m/s2g = 9{,}8\,\text{m/s}^2.)

    Show solution
    Pêndulo simples linearizado: θ¨+(g/L)θ=0\ddot\theta + (g/L)\theta = 0. ω0=g/L=9,8/13,13\omega_0 = \sqrt{g/L} = \sqrt{9{,}8/1} \approx 3{,}13 rad/s. Período: T=2π/ω02,01T = 2\pi/\omega_0 \approx 2{,}01 s.
  14. Ex. 96.14Modeling

    Isolamento de vibração: para isolar uma máquina de vibração a 4 Hz (do piso), qual deve ser a frequência natural máxima do suporte?

    Show solution
    Para isolamento efetivo: operar em r=ω/ω0>2r = \omega/\omega_0 > \sqrt{2}. Logo ω0<ωestrada/2\omega_0 < \omega_{estrada}/\sqrt{2}. Frequência de excitação típica de estrada a 80 km/h com irregularidades espaçadas 5 m: f=80000/(3600×5)4,4f = 80000/(3600 \times 5) \approx 4{,}4 Hz = 27,7 rad/s. Necessário ω0<27,7/219,6\omega_0 < 27{,}7/\sqrt{2} \approx 19{,}6 rad/s, ou seja, k/m<384k/m < 384 N/(kg m).
  15. Ex. 96.15Modeling

    Circuito RLC em série: L=0,01HL = 0{,}01\,\text{H}, C=100μFC = 100\,\mu\text{F}, R=10ΩR = 10\,\Omega. Calcule ω0\omega_0 e QQ.

    Show solution
    Circuito RLC série: Lq¨+Rq˙+q/C=V(t)L\ddot q + R\dot q + q/C = V(t). Analogia: LmL \leftrightarrow m, RcR \leftrightarrow c, 1/Ck1/C \leftrightarrow k. Frequência de ressonância: ω0=1/LC\omega_0 = 1/\sqrt{LC}. Fator de qualidade: Q=ω0L/R=1/(RC/L)Q = \omega_0 L/R = 1/(R\sqrt{C/L}).
  16. Ex. 96.16Understanding

    Como se compara a frequência amortecida ωd\omega_d com a natural ω0\omega_0 no regime subamortecido?

    Select the correct option
    Select an option first
    Show solution
    ωd=ω01ζ2<ω0\omega_d = \omega_0\sqrt{1-\zeta^2} < \omega_0 para qualquer ζ>0\zeta > 0. O amortecimento sempre reduz a frequência de oscilação. Em ζ=1\zeta = 1 (crítico): ωd=0\omega_d = 0 — sem oscilação.
  17. Ex. 96.17Understanding

    Em projeto de controle, quando se prefere amortecimento crítico versus subamortecido?

    Select the correct option
    Select an option first
    Show solution
    Subamortecido: atinge o valor alvo mais rápido (tempo de subida menor) mas ultrapassa (overshoot). Criticamente amortecido: sem overshoot e tempo de assentamento mínimo. Superamortecido: sem overshoot mas mais lento que o crítico.
  18. Ex. 96.18Application

    Duas molas com k1=200N/mk_1 = 200\,\text{N/m} e k2=300N/mk_2 = 300\,\text{N/m} conectadas em série com massa m=5kgm = 5\,\text{kg}. Calcule keqk_{\text{eq}} e ω0\omega_0.

    Show solution
    Mola em série: rigidez equivalente keq=k1k2/(k1+k2)=200×300/500=120k_{eq} = k_1 k_2/(k_1+k_2) = 200 \times 300/500 = 120 N/m. Molas em paralelo: keq=k1+k2=500k_{eq} = k_1 + k_2 = 500 N/m. ω0\omega_0 muda conforme a configuração.
  19. Ex. 96.19Application

    Para o oscilador amortecido com F=F0cosωtF = F_0\cos\omega t, escreva a fórmula da amplitude e fase do regime permanente.

    Show solution
    Regime permanente (transitório decai): xp=Acosωt+Bsinωtx_p = A\cos\omega t + B\sin\omega t. Amplitude: R=A2+B2=F0/(kmω2)2+c2ω2R = \sqrt{A^2+B^2} = F_0/\sqrt{(k-m\omega^2)^2 + c^2\omega^2}. Fase: ϕ=arctan(cω/(kmω2))\phi = \arctan(c\omega/(k-m\omega^2)).
  20. Ex. 96.20ChallengeAnswer key

    Compare a resposta em ω=ω0\omega = \omega_0 para (a) ζ=0\zeta = 0 e (b) ζ=0,05\zeta = 0{,}05. Qual a amplitude máxima em cada caso?

    Show solution
    Com F=F0cosω0tF = F_0\cos\omega_0 t e ζ=0\zeta = 0: solução particular via regra de modificação é xp=F02mω0tsinω0tx_p = \tfrac{F_0}{2m\omega_0}t\sin\omega_0 t. Para ζ>0\zeta > 0 pequeno: amplitude máxima finite F0/(2ζmω02)=F0/(2ζk)F_0/(2\zeta m\omega_0^2) = F_0/(2\zeta k).
  21. Ex. 96.21ChallengeAnswer key

    Batimento: x¨+4x=cos2,1t\ddot x + 4x = \cos 2{,}1t, x(0)=0x(0) = 0, x˙(0)=0\dot x(0) = 0. Calcule a frequência de batimento e esboce qualitativamente a solução.

    Show solution
    Batimento: soma de dois cosenos com frequências próximas. cosω1t+cosω2t=2cos ⁣(ω1+ω22t)cos ⁣(ω1ω22t)\cos\omega_1 t + \cos\omega_2 t = 2\cos\!\left(\frac{\omega_1+\omega_2}{2}t\right)\cos\!\left(\frac{\omega_1-\omega_2}{2}t\right). Frequência de batimento: ωbat=ω1ω2\omega_{bat} = |\omega_1 - \omega_2|. Para ω1=2\omega_1 = 2, ω2=2,1\omega_2 = 2{,}1: ωbat=0,1\omega_{bat} = 0{,}1 rad/s, período de batimento =2π/0,163= 2\pi/0{,}1 \approx 63 s.
  22. Ex. 96.22Challenge

    Aplique variação de parâmetros ao oscilador subamortecido x¨+2x˙+5x=etsint\ddot x + 2\dot x + 5x = e^{-t}\sin t.

    Show solution
    Variação de parâmetros: y1=etcos2ty_1 = e^{-t}\cos 2t, y2=etsin2ty_2 = e^{-t}\sin 2t, W=2e2tW = 2e^{-2t}. xp=y1y2q/Wdt+y2y1q/Wdtx_p = -y_1\int y_2 q/W\,dt + y_2\int y_1 q/W\,dt. Para q=et/sin2tq = e^{-t}/\sin 2t — integral não elementar. Se q=etsintq = e^{-t}\sin t: utilizável com integração por partes.
  23. Ex. 96.23ProofAnswer key

    Demonstre que a energia total E=12mx˙2+12kx2E = \frac{1}{2}m\dot x^2 + \frac{1}{2}kx^2 do oscilador amortecido (c>0c > 0) é estritamente decrescente.

    Show solution
    Dada a equação mx¨+cx˙+kx=0m\ddot x + c\dot x + kx = 0, compute E˙\dot E onde E=12mx˙2+12kx2E = \tfrac{1}{2}m\dot x^2 + \tfrac{1}{2}kx^2: E˙=mx˙x¨+kxx˙=x˙(mx¨+kx)=x˙(cx˙)=cx˙20\dot E = m\dot x\ddot x + kx\dot x = \dot x(m\ddot x + kx) = \dot x(-c\dot x) = -c\dot x^2 \leq 0. Logo a energia total é não-crescente; igualdade apenas em x˙=0\dot x = 0.
  24. Ex. 96.24Proof

    Use o teorema de Abel para mostrar que o Wronskiano de y1=eζω0tcosωdty_1 = e^{-\zeta\omega_0 t}\cos\omega_d t e y2=eζω0tsinωdty_2 = e^{-\zeta\omega_0 t}\sin\omega_d t é sempre não-nulo (0<ζ<10 < \zeta < 1).

    Show solution
    A equação de modo normal satisfaz λ2+2ζω0λ+ω02=0\lambda^2 + 2\zeta\omega_0\lambda + \omega_0^2 = 0. Abel: W+2ζω0W=0W=W0e2ζω0tW' + 2\zeta\omega_0 W = 0 \Rightarrow W = W_0 e^{-2\zeta\omega_0 t}. Se W00W_0 \neq 0, então W(t)0W(t) \neq 0 para todo tt — par é linearmente independente e forma base.

Fontes

Updated on 2026-05-06 · Author(s): Clube da Matemática

Found an error? Open an issue on GitHub or submit a PR — open source forever.