モード減衰の同定と設定
モード減衰の同定と設定の理論基礎
減衰とは
先生、構造力学で「減衰」はどう扱いますか?
減衰は振動のエネルギーが散逸する効果だ。減衰がなければ構造は永遠に振動し続ける。実構造には必ず減衰があり、これが共振時の振幅を有限に保つ。
運動方程式
減衰付きの運動方程式:
$[C]$ が減衰マトリクス。$[M]$ と $[K]$ は物理的に明確だが、$[C]$ は一般に不確かさが大きい。
減衰マトリクスをどう決めるかが問題なんですね。
その通り。減衰のモデル化はFEMの動的解析で最も不確かなパラメータだ。
減衰のモデル
主要な減衰モデル:
1. モード減衰(Modal Damping)
各モードに減衰比 $\zeta_i$ を直接指定。最も一般的な手法。
$q_i$ はモード座標、$\omega_i$ は固有角振動数。
2. レイリー減衰(Rayleigh Damping)
$[C] = \alpha [M] + \beta [K]$。$\alpha$ と $\beta$ は2つの振動数で $\zeta$ を合わせて決める。
3. 構造減衰(Structural Damping)
ヒステリシス減衰とも呼ばれる。周波数に依存しない減衰で、複素剛性 $[K^*] = K$ で表現。$g$ が構造減衰係数。
3つのモデル、どう使い分けますか?
減衰比の典型値
| 構造 | 減衰比 $\zeta$ |
|---|---|
| 鋼構造(溶接) | 0.5〜1% |
| 鋼構造(ボルト接合) | 1〜2% |
| RC構造 | 3〜5% |
| 免震構造 | 10〜30% |
| 機械構造 | 1〜3% |
| 複合材構造 | 0.5〜2% |
鋼構造で0.5〜1%…非常に小さいですね。
鋼は内部減衰が小さい。だから鋼構造は共振で大きな振幅が出やすく、減衰の設定が結果に大きく影響する。
まとめ
モード減衰の理論を整理します。
要点:
- 減衰は動的解析で最も不確かなパラメータ — 感度分析が必須
- 3つのモデル — モード減衰、レイリー減衰、構造減衰
- モード減衰が最も一般的 — 各モードに $\zeta_i$ を指定
- レイリー減衰は時間領域用 — $\alpha, \beta$ で2周波数を合わせる
- 減衰比の典型値 — 鋼: 1%、RC: 5%。用途で異なる
減衰の設定で結果が何倍も変わることがあるんですね。共振時の振幅は $1/(2\zeta)$ に比例するから、$\zeta = 1\%$ と $\zeta = 2\%$ で振幅が2倍違う。
だから減衰は「最も影響が大きく、最も不確かなパラメータ」。減衰に対する感度分析なしに動的解析の結果を信じてはいけない。
減衰比2%という「魔法の数字」
構造減衰比ζ=2%は設計上の慣例として広く使われるが、実際の鋼構造物では0.5〜5%と大幅に異なる。この値はLankford(1954年)が建築物の実測データから統計的に提案した中央値だ。1970年代のUBC(米国建築基準法)に採用されてから「標準値」として定着したが、溶接構造は1%未満、ボルト締結は3〜5%とまったく異なることを忘れてはならない。
モード減衰の同定と設定の数値計算手法
減衰の設定方法
FEMで減衰をどう設定しますか?
Nastranでのモード減衰
```
TABDMP1, 1, CRIT
, 0., 0.02, 100., 0.02, ENDT
```
全モードに $\zeta = 2\%$ を設定。振動数範囲ごとに異なる $\zeta$ も可能。
Abaqusでのモード減衰
```
*MODAL DAMPING
1, 50, 0.02
```
モード1〜50に $\zeta = 2\%$ を一括設定。
レイリー減衰の設定
$\alpha$ と $\beta$ の決定方法:
2つの振動数 $f_1, f_2$ で $\zeta_1 = \zeta_2 = \zeta$ とする場合:
$f_1$ と $f_2$ はどう選びますか?
着目する振動数範囲の下限と上限。例えば地震応答で1〜10 Hzを対象とするなら $f_1 = 1$ Hz, $f_2 = 10$ Hz。この範囲外では減衰比がずれるので注意。
Abaqusでのレイリー減衰
```
*DAMPING, ALPHA=0.5, BETA=0.001
```
減衰の同定
実構造の減衰比はどう測定しますか?
実験モード解析で測定:
1. ハンマー加振法 — インパルスハンマーで加振し、加速度を計測
2. 加振器法 — 正弦波/ランダム加振で周波数応答関数(FRF)を取得
3. 減衰の同定 — FRFの半値幅法、または曲線フィッティング
半値幅法:共振ピークの振幅が $1/\sqrt{2}$ になる2つの振動数 $f_1, f_2$ から:
簡単に測定できるんですね。
原理はシンプルだが、実験の精度(加振点、計測点、ノイズ処理)が結果に影響する。複数の方法で交差検証するのが安全だ。
まとめ
減衰設定の数値手法、整理します。
要点:
- TABDMP1(Nastran), *MODAL DAMPING(Abaqus) — モード減衰の設定
- $\alpha, \beta$ の決定 — 2振動数での合わせ込み
- 実験モード解析で減衰比を測定 — 半値幅法が基本
- 減衰比の範囲に注意 — レイリー減衰は指定範囲外でずれる
- 実験データがなければ文献値を使い、感度分析する
半値幅法と対数減衰率の実測手順
実験から減衰比を求める半値幅法では、FRF(周波数応答関数)の共振ピーク周波数から±3dB点の幅Δfとピークfrからζ=Δf/(2fr)で計算する。精度は周波数分解能に依存し、Δfの5分の1以下の分解能が必要だ。インパルスハンマ試験では自由減衰波形の包絡線の対数から対数減衰率を直接求める方が精度が高い場合もある。
モード減衰の同定と設定の実務適用
減衰の実務適用
減衰の設定は実務でどう行いますか?
設計コードの減衰値
設計コードで減衰比が規定されている場合がある:
| コード | 構造 | 減衰比 $\zeta$ |
|---|---|---|
| 建築基準法(日本) | RC造 | 5% |
| 建築基準法(日本) | S造 | 2% |
| ユーロコード8 | RC造 | 5% |
| ASCE 7 | 全構造 | 5%(デフォルト) |
| API 617 | 回転機械 | 実測または1% |
地震応答ではRC造5%がデフォルトなんですね。
RC造の5%は「ひび割れを含むRC構造の等価減衰」を意味する。無損傷のRC構造の減衰は1〜2%だが、地震で微小なひび割れが入ることでエネルギー散逸が増加し、等価的に5%になる。
実験データがない場合
実務チェックリスト
減衰設定のチェックリストをお願いします。
減衰比を $\zeta/2$ にすると共振振幅が2倍…感度分析は必須ですね。
減衰はFEMの動的解析で最大の不確定要素。感度分析を怠ると、結果に意味がない。
免震建物のダンパー最適化
超高層ビルの制振設計では粘性ダンパーや摩擦ダンパーをモーダル減衰比として換算し、目標減衰比(通常3〜5%)を達成するダンパー配置を最適化する。六本木ヒルズ森タワー(238m)では粘性制振ダンパー設置後に低次モードの減衰比が1%から4%に上昇し、地震応答加速度を40%低減させた。
モード減衰の同定と設定のソフトウェア比較
減衰設定のツール
減衰設定に関するソルバー間の差は?
材料ごとに異なる減衰を設定できますか?
NastranのMAT1カードのGEフィールドで材料ごとの構造減衰係数を設定できる。鋼とゴムで異なる減衰を持つ構造で有用。AbaqusやAnsysでも同様の機能がある。
実験モード解析ツール
| ツール | 特徴 |
|---|---|
| LMS Test.Lab | FRFからの減衰同定。FEMとの相関 |
| ME'scope | モード形状と減衰の可視化 |
| ModalVIEW | Excelベース。手軽なモード解析 |
選定ガイド
解析手法によって使える減衰モデルが異なるんですね。
モード減衰はモード重畳法、レイリー減衰は直接積分法、構造減衰は周波数応答。これを間違えると物理的に意味のない結果が出る。
MSC Nastranモーダル減衰入力形式
Nastran SOL 108(周波数応答)ではTABDMP1カードでモーダル減衰比の周波数依存性を定義できる。各モードに個別の減衰比を設定するCRITDAMP=MODEオプションも使え、実測値をそのまま入力可能だ。LMS(現Siemens Simcenter)社はこの機能を活用したテスト-解析相関フローを構築し、試験から1時間以内にFEMの減衰モデルを更新するシステムを自動車向けに提供している。
モード減衰の同定と設定の先端研究
減衰の先端研究
減衰モデルの最前線を教えてください。
非粘性減衰(Non-viscous Damping)
レイリー減衰は粘性減衰($[C]\{\dot{u}\}$、速度に比例)を仮定するが、実構造の減衰は非粘性であることが多い。摩擦減衰(接合部の微小すべり)、ヒステリシス減衰(材料の非線形履歴)は粘性モデルでは正確に表現できない。
一般化減衰モデル:
$[G(t)]$ は減衰核関数。粘性減衰は $[G(t)] = [C]\delta(t)$ の特殊ケース。
接合部の減衰
ボルト接合部の微小すべり(フレッティング)が構造全体の減衰の大部分を占めることがある。接合部減衰モデル(Iwan要素、Jenkins要素)が研究されている。
接合部の微小すべりが減衰の主因?
溶接の鋼構造では材料自体の内部減衰が非常に小さい($\zeta < 0.1\%$)。実測の1〜2%の減衰は接合部のフレッティングと非構造材(塗装、ケーブル等)の摩擦が主因。
能動減衰(Active Damping)
圧電素子やMRダンパーで能動的に減衰を制御する。FEMではセンサー-アクチュエータのフィードバック制御と構造解析を連成する。
まとめ
減衰の先端研究、まとめます。
減衰は構造力学の「最もソフトな話題」だが、振動制御の実効性を決める最も重要なパラメータだ。
材料損失係数とモーダル減衰の関係
構造物全体のモーダル減衰比ζnは各部材の材料損失係数ηと部材がモードで蓄えるひずみエネルギーの比率で決まる(モーダルひずみエネルギー法)。高減衰ゴム(η=0.3〜1.0)を少量挟むだけで構造全体の減衰比を大幅に高められるのはこの原理による。制振鋼板(η=0.1〜0.3)を自動車フロアに使うと室内騒音を3〜6dB低減できる。
モード減衰の同定と設定のトラブル対応
減衰設定のトラブル
減衰設定でよくあるトラブルは?
共振振幅が非現実的
共振での応答が非常に大きくなります。
減衰が小さすぎる。減衰比を確認。$\zeta = 0$(減衰なし)だと共振で振幅が無限大。
確認:
- 減衰が設定されているか(設定忘れが最多)
- 減衰比の値が適切か(文献値・コード値と比較)
- $\zeta = 0.01$ で振幅が $1/(2 \times 0.01) = 50$ 倍。これが妥当か
レイリー減衰で低周波/高周波の応答がおかしい
レイリー減衰は $\alpha$ 項が低周波で支配的(減衰大)、$\beta$ 項が高周波で支配的(減衰大)。指定した2周波数の間では $\zeta$ が適切だが、範囲外では減衰が過大になる。
どう対処しますか?
構造減衰を時間領域で使ってしまう
時間領域ではモード減衰かレイリー減衰を使うべきですね。
そう。周波数領域 = 構造減衰OK。時間領域 = レイリー or モード減衰。これを間違えると因果律が破れて非物理的な応答が出る。
まとめ
減衰設定のトラブル対処、整理します。
モーダル減衰比が実測より低く解析される場合
FEM解析でのモーダル減衰比は実測より低く見積もられることが多い。これはFEMが接合部(ボルト・リベット・溶接)での摩擦・微小すべりによる散逸エネルギーを表現できないためだ。接合部の減衰を集中減衰要素(C要素)として付加するか、実測減衰比をモーダル解析にそのまま入力するのが実務的な対処法だ。
なった
詳しく
報告