減衰振動(ゲーム)

上に壁・ばね・ダンパー・台車、下に台車の位置の記録(横軸は時間)。0 の線を下に越えたら行き過ぎ。

ドパリカ

CRITICAL DAMPING

減衰振動 — 臨界減衰 c = 2√(mk) で最も速く戻る

短い答え: 運動方程式 mx¨+cx˙+kx=0m\ddot x + c\dot x + kx = 0 の解は、減衰比 ζ=c/(2mk)\zeta = c/(2\sqrt{mk}) で 3 通りに分かれる。 ζ<1\zeta < 1 は行き過ぎて振動しながら減衰、ζ>1\zeta > 1 はゆっくり戻る。 行き過ぎずに最も速く戻るのは臨界減衰 c=2mkc = 2\sqrt{mk}。

前提と約束

式で納得: 特性方程式で 3 通りに分かれる

x=eλtx = e^{\lambda t} を代入すると、特性方程式と根は

mλ2+cλ+k=0  ⟹  λ=−c±c2−4mk2mm\lambda^2 + c\lambda + k = 0 \;\Longrightarrow\; \lambda = \frac{-c \pm \sqrt{c^2 - 4mk}}{2m}

根号の中の符号で運動が変わる。境目 c2=4mkc^2 = 4mk、つまり cc=2mkc_c = 2\sqrt{mk} が臨界減衰。ω0=k/m\omega_0 = \sqrt{k/m}、ζ=c/cc\zeta = c/c_c とおく。

不足減衰(ζ < 1)

根は複素数 −ζω0±iωd-\zeta\omega_0 \pm i\omega_d(ωd=ω01−ζ2\omega_d = \omega_0\sqrt{1-\zeta^2})。 x=x0e−ζω0t(cos⁡ωdt+ζω0ωdsin⁡ωdt)x = x_0 e^{-\zeta\omega_0 t}\left(\cos\omega_d t + \tfrac{\zeta\omega_0}{\omega_d}\sin\omega_d t\right)。0 を行き過ぎて、振幅が指数関数的に減る。 最初の行き過ぎの量は x0 e−πζ/1−ζ2x_0\, e^{-\pi\zeta/\sqrt{1-\zeta^2}}。

臨界減衰(ζ = 1)

重根 λ=−ω0\lambda = -\omega_0。x=x0(1+ω0t)e−ω0tx = x_0(1+\omega_0 t)e^{-\omega_0 t}。行き過ぎず、最も速く 0 に近づく。

過減衰(ζ > 1)

2 つの負の実根。遅い方の根 −ζω0+ω0ζ2−1-\zeta\omega_0 + \omega_0\sqrt{\zeta^2-1} は ζ を大きくするほど 0 に近づくので、ダンパーを強くするほど戻るのが遅くなる。

√ が効く

cc=2mkc_c = 2\sqrt{mk}。台車を 4 倍重くすると臨界のダンパーは 2 倍、ばねを 9 倍かたくすると 3 倍。 ゲームの目盛りは「2 × 重さの一辺 × ばねの一辺」。

例題

質量 1.0 kg、ばね定数 4π2≈39.5 N/m4\pi^2 \approx 39.5\ \mathrm{N/m} の振動系に、減衰係数 c=2π≈6.3 N s/mc = 2\pi \approx 6.3\ \mathrm{N\,s/m} のダンパーを付けた。減衰比・減衰振動の周期・最初の行き過ぎの割合を求めよ。

解答

cc=2mk=4πc_c = 2\sqrt{mk} = 4\pi なので ζ=0.5\zeta = 0.5。 ωd=2π1−0.25≈5.44 rad/s\omega_d = 2\pi\sqrt{1-0.25} \approx 5.44\ \mathrm{rad/s}、周期 2π/ωd≈1.15 s2\pi/\omega_d \approx 1.15\ \mathrm{s}(減衰なしの 1 秒より長い)。 行き過ぎは e−π×0.5/0.75≈0.16e^{-\pi\times0.5/\sqrt{0.75}} \approx 0.16、つまり 16%。

類題

  1. 車輪 1 本ぶんの質量 300 kg、サスペンションのばね定数 3.0×104 N/m3.0\times10^{4}\ \mathrm{N/m}。臨界減衰になるショックアブソーバーの減衰係数は?

    答え

    cc=2300×3.0×104=2×3000=6.0×103 N s/mc_c = 2\sqrt{300\times 3.0\times10^4} = 2\times 3000 = 6.0\times10^{3}\ \mathrm{N\,s/m}。

  2. 不足減衰で、隣り合う山の振幅の比の自然対数(対数減衰率)δ=ln⁡(xn/xn+1)\delta = \ln(x_n/x_{n+1}) を ζ で表せ。

    答え

    1 周期 2π/ωd2\pi/\omega_d で e−ζω0te^{-\zeta\omega_0 t} が減るので δ=2πζ1−ζ2\delta = \dfrac{2\pi\zeta}{\sqrt{1-\zeta^2}}。

  3. コイル L・抵抗 R・コンデンサ C の直列回路で、コンデンサの電荷が振動せずに最も速く放電する R は?

    答え

    Lq¨+Rq˙+q/C=0L\ddot q + R\dot q + q/C = 0 は台車の式と同じ形(m → L、c → R、k → 1/C)。R=2L/CR = 2\sqrt{L/C}。

モデルの仮定と扱わないこと

ゲームの表現について

シリーズ内の位置

力学(基礎): 斜方投射 → 単振動。大学: 減衰振動(この教材)。

更新日: ・ 誤りのご報告はお問い合わせへ ・ プライバシー