振動・波動の基礎-⑫2自由度系の減衰自由振動について

[voice icon="https://deltapower.site/wp-content/uploads/2020/04/男性悩み.png" name="物理苦手君" type="l"]2自由度の減衰ありで運動方程式を解けないんだけど…[/voice]

[voice icon="https://deltapower.site/wp-content/uploads/2020/09/にっこり右側.jpg" name="デルタ先生" type="r"]減衰ありの状態で計算するのは、凄く難しくなってくるんだ。考え方はこれまでと同じでもいいんだけど、計算力が必須だね。[/voice]

本記事では、2自由度の減衰自由振動について解説します。

注意事項ですが、過去記事のように、本記事では方程式を最後まで解いていません。

というか 手計算ではそうそう簡単に解けません。

ということを解説します。

[box class="green_box" title="本記事をおススメする人"]

2自由度の振動について減衰がある場合を理解したい人
物理を勉強していて興味がある人、仕事で振動の知識を使う人

[/box]

1 モデルと運動方程式
- 1.1 運動方程式を解く
2 どのような挙動をするのか
3 まとめ
4 参考文献

モデルと運動方程式

モデルは下記にようになります。

前回までは、減衰部分の$c_1$、$c_2$が無い非減衰の場合を考えていましたが、今回はこの減衰項が入った状態で考えます。

1自由度のときと同じく、速度に比例した減衰力が発生すると考えた場合、運動方程式は、下記のようになります。

$$m_1\ddot{x_1}+(c_1+c_2)\dot{x_1}+(k_1+k_2)x_1-c_2\dot{x_2}-k_2x_2=0・・・①$$

$$m_2\ddot{x_2}+c_2\dot{x_2}+k_2x_2-c_2\dot{x_1}-k_2x_1=0・・・②$$

連立の線形微分方程式ですね。

運動方程式を解く

さて、1自由度の減衰あり運動方程式の解き方は覚えていますでしょうか？

下記の記事を参考にしてくださいね。

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

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

[/box]

1回微分の項があるので、特殊解を三角関数ではなく、eの指数関数で置くのでしたね。

この解法を、2自由度の非減衰自由振動の解き方に応用してみましょう。

特殊解を$x_1=A_1e^{λt}、x_2=A_2e^{λt}$とすると、

$\dot{x_1}=A_1λe^{λt}、\dot{x_2}=A_2λe^{λt}、\ddot{x_1}=A_1λ^2e^{λt}、\ddot{x_2}=A_2λ^2e^{λt}$となりますので、

これらをすべて①、②式の運動方程式に代入しましょう。

①式は

$$m_1A_1λ^2e^{λt}+(c_1+c_2)A_1λe^{λt}+(k_1+k_2)A_1e^{λt}-c_2A_2λe^{λt}-k_2A_2e^{λt}=0$$

$e^{λt}$を消去できるので、消去すると、

$$m_1A_1λ^2+(c_1+c_2)A_1λ+(k_1+k_2)A_1-c_2A_2λ-k_2A_2=0・・・③$$

同様に②式は、

$$m_2A_2λ^2e^{λt}+c_2A_2λe^{λt}+k_2A_2-e^{λt}c_2A_1λe^{λt}-k_2A_1e^{λt}=0$$

これも$e^{λt}$を消去して、

$$m_2A_2λ^2+c_2A_2λ+k_2A_2-c_2A_1λ-k_2A_1=0・・・④$$

③と④を連立させ、行列で書くと、

$$\left(\begin{array}{cc}m_1λ^2+(c_1+c_2)λ+(k_1+k_2)&-c_2λ-k_2\\-c_2λ-k_2&m_2λ^2+c_2λ+k_2)\end{array}\right)\left(\begin{array}{c}A_1\\A_2\end{array}\right)=0$$

ここから、$A_1、A_2$がともにゼロでないという、1次従属の条件から、λを求めます。

行列式=0となることが、その条件ですので、

$$det\left|\begin{array}{cc}m_1λ^2+(c_1+c_2)λ+(k_1+k_2)&-c_2λ-k_2\\-c_2λ-k_2&m_2λ^2+c_2λ+k_2)\end{array}\right|=0$$

整理すると、

$$(m_1λ^2+(c_1+c_2)λ+(k_1+k_2))(m_2λ^2+c_2λ+k_2)-(-c_2λ-k_2)(-c_2λ-k_2)=0$$

[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"]λの4次方程式になっているね。しかも$λ^3$の項も残るから、高度な数学を使わないと難しいんだ。[/voice]

結論、解けないことはないけど、非常に計算が難しい、ということです。

実は、減衰がある2自由度の自由振動や強制振動について、詳細に解いた教科書も少ないのが現状です。

計算は複雑なので略、みたいな形で書いてある本が多いですね。