CAEソリューション
振動解析の数値計算
第三回:1自由度減衰系の自由振動
今回は前回無視した減衰を考慮した自由振動について考えていきたいと思います。
減衰は振動を抑制する働きを持ちますが、これは系に存在するエネルギーを散逸させていることになります。物理モデルでは下図のようにダンパーで表現され、減衰係数と速度を乗じた力として表現されます。イメージとしては物体の周囲が粘度の高い流体で満たされていると振動が止まることや、物体の速度が速いと抵抗力が大きくなることが想像できるかと思います。また、減衰係数がある値より大きくなると物体は振動せずに止まりますが、その閾値を臨界減衰係数といいます。
![ダンパーで表現された減衰係数と速度を乗じた力の図](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_001.jpg)
それでは、減衰を考慮した自由振動の解を求めていきましょう。
減衰を考慮した1自由度自由振動の運動方程式は下記式①になります。
![減衰を考慮した1自由度自由振動の運動方程式](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_002.jpg)
上記は、ある関数とその関数を1回微分および2回微分したものに係数(m,c,k)を乗じた和が0になる必要があるため、これを満足する関数は下記のような複素指数関数に限定されます。
![複素指数関数](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_003.jpg)
式②を微分すると
![式②を微分した式](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_004.jpg)
式②と式③を式①に代入すると
![式②と式③を式①に代入の式](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_005.jpg)
式④は jω に関する2次方程式ですので下記解が得られます。
![式④は jω に関する2次方程式](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_006.jpg)
ここで、CC=2mω0 , ζ=C/CC と新たに定義すると
![C_C=2mω_0 , ζ=C/C_C と新たに定義](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_007.jpg)
となります。
ω0は不減衰固有角振動数、CCは臨界減衰係数、ζは減衰比と呼ばれます。
また、jωは2通りあるため、jω1,jω2として式②は下記のようになります。
![jωは2通りあるため、jω_1,jω_2として式②は下記の図になる](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_008.jpg)
式⑥には平方根があり、ζ2の値によって平方根の中身の正負が変わりますので、場合分けを見ていきましょう。
![式⑥には平方根があり、ζ^2の値によって平方根の中身の正負が変わりますので、場合分けの図](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_009.jpg)
ここで、
![式⑨と定義](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_012.jpg)
と定義し、式⑧に代入すると、
![式⑧に代入](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_010.jpg)
となります。
式⑩は、実部(-σ)が減衰の速さ、虚部( jωd)が振動の速さを表現しています。
![減衰を考慮した自由振動のグラフ](./landingpage/common/imgs/cae_column_ichijiyuudogensuikei_011.jpg)
今回は減衰を考慮した1自由度バネマスダンパーモデルの自由振動について説明いたしました。次回は加振力が作用する際の強制振動について説明する予定です。
疑問やお困りごとがございましたら、以下お問い合わせより、ご相談ください。