減衰振動における自由振動の計算-過減衰の形と双曲線関数について解説

この記事は下記の減衰振動で出てきた2階の線形微分方程式を解いた結果を、もう少し詳しく解説します。

[box class="blue_box" title="参考記事"]

振動・波動の基礎②-振動の減衰

[/box]

まずは復習で2階の線形微分方程式について、説明し、双曲線関数の説明をします。

1 2階の線形微分方程式
- 1.1 解の求め方（簡易版）
2 双曲線関数
- 2.1 関数の形
- 2.2 グラフ
3 一般解の求め方
- 3.1 双曲線関数での表記の導出
4 まとめ
5 参考文献

2階の線形微分方程式

下記のような式を2階の線形微分方程式を考えます。

$$\ddot{x}+2ζω\dot{x}+ω^2x = 0$$

$$ζ = \frac{c}{2\sqrt{mk}}$$

$$ω = \sqrt{\frac{k}{m}}$$

これは減衰を考慮した振動の運動方程式ですね。

ζは減衰比、ωは固有振動数を示します。

解の求め方（簡易版）

詳細は参考記事を見てくださいね。

[box class="blue_box" title="参考記事"]

振動・波動の基礎②-振動の減衰

[/box]

これを解くためには、微分しても関数系が変わらない解を仮定してから、一般解を求めます。

$$x = Ae^{λt}　Aは定数$$

このように指数関数の解を仮定して、運動方程式に代入すると、

$$A(λ^2+2ζωλ+ω^2)e^{λt} = 0$$

となるので、λの部分の2次方程式がゼロになると、すべての$t$で成立しますね。

$$λ^2+2ζωλ+ω^2 = 0$$

この2次方程式を特性方程式と言い、この解は

$$λ = -ζω±ω\sqrt{ζ^2-1}$$

となります。

このルートの中の$ζ^2-1$が1以上のとき、どのようになるか、について解説します。

双曲線関数

双曲線関数というものを使用しますので、この関数について説明します。

関数の形

[box class="glay_box" title="双曲線関数"]

$sinhx = \frac{e^x-e^{-x}}{2}$
$coshx = \frac{e^x+e^{-x}}{2}$
$tanhx = \frac{e^x-e^{-x}}{e^x+e{-x}}$

[/box]

グラフ

$sinhx、coshx、tanhx$のグラフを下記に示します。

[voice icon="https://deltapower.site/wp-content/uploads/2020/04/男性悩み.png" name="物理苦手君" type="l"]これって振動と関係あるの？？[/voice]

[voice icon="https://deltapower.site/wp-content/uploads/2020/09/にっこり右側.jpg" name="デルタ先生" type="r"]振動している領域では使わないんだけど、減衰を大きくした過減衰の状態で少し使うんだよ[/voice]

一般解の求め方

線形微分方程式の一般解は、任意定数を含む2つの解の足し合わせ、となります。

$$x(t) = C_1 x_1+C_2 x_2　C_1、C_2は定数$$

$x_1 = e^{λ_1 t}$、$x_2 = e^{λ_2 t}$としますと、下記のようになりますね。

$$x(t) = C_1 e^{λ_1 t}+C_2 e^{λ_2 t}$$

特性方程式の解は下記の通りで、$ζ>1$のときを考えます。

$$λ = -ζω±ω\sqrt{ζ^2-1}$$

二つの解を$λ_1$、$λ_2$とすると

$$λ_1 = -ζω＋ω\sqrt{ζ^2-1} = ω(-ζ+\sqrt{ζ^2-1}) $$

$$λ_2 = -ζω-ω\sqrt{ζ^2-1} =ω(-ζ-\sqrt{ζ^2-1}) $$

となり、これを一般解の式に代入して、整理します。

$$x(t) = e^{-ωζt}(C_1 e^{-\sqrt{ζ^2-1}・ωt}+C_2 e^{\sqrt{ζ^2-1}・ωt})$$

これが過減衰の式ということなのですが、本記事ではこれを双曲線関数で示します。

双曲線関数での表記の導出

まず、$ω\sqrt{ζ^2-1} = γ$と置いて変形すると

$$x(t) = e^{-ωζt}(C_1 e^{-γt}+C_2 e^{γt})$$

この式を双曲線関数で示そうと思います。

適当な定数$A$、$B$を用意し、$Acoshx$と$Bsinhx$の和の式を変形していきます。

$$Acoshx + Bsinhx = A\frac{e^x+e^{-x}}{2} + B\frac{e^x-e^{-x}}{2}$$

これらを$e^x$と$e^{-x}$で整理しますと、

$$Acoshx + Bsinhx =\frac{A+B}{2}e^x +\frac{A-B}{2}e^{-x}$$

となります。

運動方程式の一般解の式に似てきましたね。

$C_1= \frac{A+B}{2}、C_2=\frac{A-B}{2}、x = γt$とすると

$$Acoshx + Bsinhx = C_1e^{γt} + C_2e^{-γt}$$

この式から一般解は双曲線関数を使って表すと下記のようになります。

$$x(t) = e^{-ωζt}(Acosh(γt) + Bsinh(γt))$$

さて、この式がどのような関数になるか、についてが重要ですね。

eの指数部分ωζと双曲線関数のγのどちらが大きいかで、関数の形は変わります。

（双曲線関数もeの指数関数なので、どちらが強いかは指数関数の指数部分の大きさ次第ということです）

しかし、$ω\sqrt{ζ^2-1} = γ$と置いていますので、

$ωζ > γ =ω\sqrt{ζ^2-1}$　ですね。

この場合のグラフを下記に示します。

なお、$fx = e^{2x}(coshx+sinhx)$として黒線で示しています。

黒線は単調減少のグラフになっていることがわかりますね。

ちなみに緑の線はζ=1で双曲線関数がゼロになったとき、すなわち臨界減衰の形と一致します。

最後にζで場合分けした減衰の振動のグラフを示して終わります。

まとめ

過減衰の部分の計算について詳細を説明しました。

双曲線関数を$exp$関数で表記すれば、一般解からグラフの形は想像しやすくなりますね。

どのような関数か、イメージをもっていただければ理解しやすくなりますので、関数の形は覚えておきましょう。

参考文献

振動工学の基礎：岩壺卓三、松久寛、森北出版株式会社
機械力学-振動の基礎から制御まで：日高照晃、小田哲、川辺尚志、曽我部雄次、吉田和信、朝倉書店
構造と連続体の力学基礎：熊でもわかる変形できる物体の力学：岩熊哲夫、小山茂 web版