減衰 (Damping) — CAE用語解説
減衰(Damping)とは
減衰って結局なんですか? 振動がだんだん小さくなることですか?
そう、振動系のエネルギーが熱や音として散逸して、振幅が徐々に小さくなる現象のこと。車のサスペンションのショックアブソーバーが典型例で、段差を乗り越えてもすぐ揺れが収まるのは減衰のおかげだ。
減衰がないとどうなるんですか?
減衰ゼロだと一度振動が始まったら永遠に止まらない。さらに共振が起きると振幅が理論上は無限大になる。実際には材料の降伏や破壊で止まるけど、橋が風で崩壊した1940年のタコマナローズ橋事故は、減衰不足と共振が重なった有名な事例だよ。
粘性減衰と構造減衰
粘性減衰と構造減衰ってどう違うんですか? 教科書で両方出てきて混乱してます。
粘性減衰(Viscous Damping)は速度に比例する減衰力で、ダッシュポットのイメージだ。運動方程式では
と書く。$c$ は減衰係数で、単位は N·s/m。油圧ダンパーのように流体の粘性抵抗で減衰させる機構がこれに近い。
じゃあ構造減衰はどういうメカニズムなんですか?
構造減衰(Structural Damping / Hysteretic Damping)は、材料内部の微視的な塑性変形や粒界すべりによるエネルギー散逸をモデル化したものだ。減衰力は変位に比例するけど位相が90度ずれていて、複素剛性で表現する:
ここで $\eta$ は損失係数(Loss Factor)と呼ばれる。構造減衰の大きな特徴は、周波数に依存しないこと。鉄鋼構造物の $\eta$ は 0.001〜0.01 程度、ゴムだと 0.1〜1.0 にもなる。実構造物の挙動は粘性減衰よりこちらに近いことが多いよ。
どっちを使えばいいんですか?
臨界減衰と減衰比
臨界減衰ってどういう状態ですか? 「減衰比 $\zeta = 1$」って書いてあるんですけど、ピンときません。
1自由度系の運動方程式
の特性方程式の判別式がちょうどゼロになる減衰を臨界減衰という。臨界減衰係数は
$$c_{cr} = 2\sqrt{km} = 2m\omega_n$$で、減衰比はこれとの比だ:
$$\zeta = \frac{c}{c_{cr}} = \frac{c}{2m\omega_n}$$$\zeta$ の値で振動の仕方がどう変わるんですか?
3つのケースがある:
- $\zeta < 1$(不足減衰):振動しながら減衰する。工学的な構造物のほとんどがこれ。鋼構造で $\zeta \approx 0.01$〜$0.02$、コンクリートで $0.03$〜$0.07$ が目安だ
- $\zeta = 1$(臨界減衰):振動せずに最速で平衡位置に戻る。ドアクローザーはこの付近に設計されている
- $\zeta > 1$(過減衰):振動しないけど戻りが遅い。精密計測器のダンパーなど、オーバーシュートを絶対に避けたい場合に使う
不足減衰系の自由振動の解は
$$x(t) = X e^{-\zeta\omega_n t} \cos(\omega_d t - \phi)$$ここで $\omega_d = \omega_n\sqrt{1-\zeta^2}$ が減衰固有振動数だ。$\zeta$ が小さければ $\omega_d \approx \omega_n$ だけど、$\zeta = 0.3$ を超えると差が無視できなくなるよ。
Rayleigh減衰
CAEの動解析で「Rayleigh減衰」ってよく出てきますけど、これはどういう減衰モデルなんですか?
多自由度系の運動方程式は
で、問題は減衰マトリクス $[C]$ をどう作るかだ。Rayleigh減衰は
$$[C] = \alpha [M] + \beta [K]$$と、質量マトリクスと剛性マトリクスの線形結合で構成する。$\alpha$(質量比例項)と $\beta$(剛性比例項)の2つのパラメータで定義されるんだ。
$\alpha$ と $\beta$ はどうやって決めるんですか?
第 $i$ モードの減衰比は
だから、2つの角周波数 $\omega_1, \omega_2$ でそれぞれ目標の減衰比 $\zeta_1, \zeta_2$ を決めると、連立方程式で $\alpha, \beta$ が求まる。実務でよくやるのは、$\zeta_1 = \zeta_2 = \zeta$(同じ減衰比)と仮定して
$$\alpha = \frac{2\omega_1 \omega_2}{\omega_1 + \omega_2}\zeta, \quad \beta = \frac{2}{\omega_1 + \omega_2}\zeta$$とする方法だ。
$\omega_1$ と $\omega_2$ の選び方にコツはありますか?
これがRayleigh減衰の一番のポイントだ。$\omega_1$ と $\omega_2$ の間の周波数帯では減衰比が目標値より小さくなり(減衰不足)、外側では大きくなる(過減衰)。だから 応答に寄与する主要な周波数帯をカバーするように選ぶ のが鉄則。例えば地震応答解析なら、1次固有振動数と、入力地震動の卓越周波数を含む範囲にするのが一般的だ。$\omega_1$ を1次固有振動数の0.5〜1倍、$\omega_2$ を3〜5次固有振動数あたりにする実務者が多いよ。
モーダル減衰比
「モード解析で各モードに減衰比を直接指定する」って聞いたんですけど、Rayleigh減衰との違いは何ですか?
モーダル減衰(Modal Damping)は、モード重ね合わせ法を使うときに、各モードに個別の減衰比 $\zeta_i$ を直接与える方法だ。Rayleigh減衰のようなU字型の周波数依存性に縛られないから、実験データに基づいて各モードに正確な減衰比を設定できる。
それなら全部モーダル減衰でいいのでは?
実験からモーダル減衰比を求めるにはどうすればいいですか?
代表的な方法は3つある:
- 半値幅法:周波数応答関数のピークから$-3$ dB下がった帯域幅 $\Delta\omega$ を使って $\zeta = \Delta\omega / (2\omega_n)$
- 対数減衰率:自由振動波形の連続するピーク比から $\delta = \ln(x_n / x_{n+1})$、$\zeta = \delta / \sqrt{4\pi^2 + \delta^2}$
- カーブフィッティング:周波数応答関数を理論式にフィッティングして同定。多モードが近接している場合に有効
自動車のボディだとモードによって $\zeta$ が 0.01〜0.05 くらいばらつくから、一律の値を使うより実験同定する価値は大きいよ。
CAE実務での減衰設定
実際にCAEで減衰を設定するとき、実験データがない場合はどうするんですか?
文献値や規格値を使うのが現実的だ。参考までに代表的な値を挙げると:
- 鋼構造(溶接):$\zeta = 0.02$〜$0.03$
- 鋼構造(ボルト接合):$\zeta = 0.03$〜$0.05$(接合部の摩擦で増える)
- 鉄筋コンクリート:$\zeta = 0.03$〜$0.07$
- アルミ合金:$\zeta = 0.002$〜$0.01$
- ゴム・高分子材料:$\zeta = 0.05$〜$0.15$
ただし、ボルト接合部のすべりや非構造部材の影響で実測値はこれより大きくなることが多い。迷ったら複数の減衰比でパラメトリックスタディを回して、結果の感度を確認するのが良い実務だよ。
減衰を大きく設定しすぎるとどうなりますか?
応答が過小評価されて危険側の設計になる。特にRayleigh減衰で $\omega_2$ を低く取りすぎると、高次モードに過大な減衰がかかって、衝撃荷重のような高周波成分を含む応答が不当に抑えられてしまう。逆に減衰を小さくしすぎると応答が過大になるから保守側だけど、不経済な設計になる。いずれにしても、減衰は解析結果に大きな影響を与えるのに不確かさが最も大きいパラメータだという認識が大事だ。
CAE用語の正確な理解は、チーム内のコミュニケーションの基盤です。 — Project NovaSolverは実務者の学習支援も視野に入れています。
減衰 (Damping)の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
お問い合わせ(準備中)関連トピック
なった
詳しく
報告