構造減衰（ヒステリシス減衰） — 複素剛性・損失係数・FEM実装
理論と物理
複素剛性と損失係数の定義
先生、振動解析で「構造減衰」って出てきたんですけど、普通の減衰と何が違うんですか？ ダンパーみたいなものとは別物ですよね？
いい質問だ。ダッシュポット（粘性ダンパー）の減衰力は速度に比例する — $F_d = c\dot{x}$。これが粘性減衰だ。一方、構造減衰（ヒステリシス減衰）は、材料内部の微視的な摩擦やすべりによるエネルギー散逸を表す。数学的には複素剛性で表現する：
ここで $\eta$ が損失係数（loss factor）だ。実部 $K$ が通常の弾性剛性、虚部 $K\eta$ がエネルギー散逸の大きさを表している。
複素数で剛性を表すって、物理的にはどういう意味ですか？ 虚数のばねってイメージしにくいんですが…
こう考えるといい。調和振動 $x = X e^{i\omega t}$ を仮定すると、復元力は：
実部は変位と同位相の弾性復元力、虚部は変位から90°進んだ力 — つまり速度成分に対応する散逸力だ。粘性減衰力 $c\dot{x} = ci\omega x$ と比較すると、構造減衰の散逸力は $\omega$ に依存しない。これが決定的な違いだよ。
粘性減衰との本質的な違い
周波数に依存しないっていうのは、具体的にどういう差が出るんですか？ 例えば自動車の振動解析とかだと。
粘性減衰では等価減衰比が周波数に比例する：
だから20 Hzと200 Hzでは減衰比が10倍も違ってしまう。自動車のダッシュボード振動解析で50〜500 Hzの広帯域を扱うとき、粘性減衰だけでは低周波で減衰不足、高周波で過減衰になる。
構造減衰なら等価減衰比は：
全周波数で一定。実際の金属やFRP構造の減衰特性は、構造減衰モデルのほうが圧倒的に近い。だから周波数応答解析（NVH解析）では構造減衰が標準なんだ。
じゃあ $\eta$ と $\zeta$ の関係は $\eta \approx 2\zeta$ ってことですね。モード解析の減衰比がわかっていれば変換できる？
その通り。共振周波数での等価関係として：
例えば実験モード解析で $\zeta = 1\%$ と同定されたら、$\eta = 0.02$ と設定する。ただし注意点があって、この等価関係は共振周波数でのみ厳密に成り立つ。共振から離れた周波数では粘性減衰と構造減衰の応答は異なってくる。
エネルギー散逸とヒステリシスループ
「ヒステリシス減衰」って呼ばれるのは、応力-ひずみ曲線のヒステリシスループと関係があるんですか？
まさにそうだ。材料に繰り返し荷重をかけると、応力-ひずみ曲線が閉じたループを描く。このループの面積がエネルギー散逸量 $\Delta W$ に対応する。損失係数は1サイクルあたりの散逸エネルギーと最大弾性エネルギーの比で定義できる：
ここで $W_{\max} = \frac{1}{2}K X^2$ は最大弾性ひずみエネルギー。例えば鋼材のループ面積はごく小さいから $\eta \approx 0.001$〜$0.01$、ゴムのループは大きいから $\eta \approx 0.2$〜$1.0$ になる。文字通り「ヒステリシスの大きさ」が減衰の強さなんだ。
半値幅法による損失係数の同定
実験で損失係数を測定するにはどうすればいいですか？ インパルスハンマー試験から求められますか？
最も一般的なのが半値幅法（Half-Power Bandwidth Method）だ。FRFの共振ピークから損失係数を読み取る。手順はこうだ：
- FRFの共振ピーク周波数 $f_n$ と最大振幅 $|H|_{\max}$ を特定
- 振幅が $|H|_{\max}/\sqrt{2}$（= -3 dB）になる2つの周波数 $f_1, f_2$ を求める
- 損失係数を計算：
例えばインパルスハンマー試験で共振が200 Hzにあり、-3 dBの帯域幅が2 Hzなら $\eta = 2/200 = 0.01$。これは $\zeta = 0.5\%$ に相当する。鋼構造としては溶接構造の典型的な値だね。
半値幅法って、ピークが近接しているモードだとうまくいかなそうですね。
鋭いね。モードが密集していると-3 dBラインが隣のピークと重なって正確に読めない。その場合はCircle Fit法（Nyquistプロット上で円弧フィッティング）やカーブフィッティング法（RFP法、LSCE法など）を使う。特にCFRP積層板のように減衰が大きくモードが重なりやすい材料では、単純な半値幅法は避けたほうがいい。
複素剛性は航空機のフラッター解析から生まれた
構造減衰を複素剛性 $K^* = K(1+i\eta)$ で表すモデルは、1930年代の航空機翼フラッター研究で誕生した。翼の弾性変形と空気力学的な力の連成振動を解くために、材料の散逸特性を簡潔に表現する方法が必要だった。1960年にCaltech の T.K. Caughey が複素剛性と粘性減衰の等価関係を厳密に整理し、$\zeta = \eta/2$ という実用式が確立された。現在もNastranのGEパラメータとしてそのまま使われている。
数値解法と実装
直接法とモード法での構造減衰
構造減衰って周波数応答解析で使うんですよね。直接法とモード法のどちらで使うんですか？
どちらでも使える。直接法（Nastran SOL 108）では、複素剛性を運動方程式に直接代入する：
モード法（Nastran SOL 111）では、固有モードに展開してからモードごとの減衰を適用する：
ここで $q_r$ はモード座標、$\phi_r$ はモードベクトル。実務ではモード法のほうが圧倒的に多い。計算コストが直接法の1/10〜1/100で済むからね。
モード法だとモードごとに $\eta_r$ を変えられるんですか？ 実験で各モードの減衰比が違うことは多いですよね。
NastranのTABDMP1カードを使えば、周波数に対する減衰のテーブルを定義できる。各モードの固有振動数に応じた $\eta$ を自動補間してくれる。例えば実験モード解析で1次モード $\zeta_1=0.5\%$、2次モード $\zeta_2=1.2\%$ と同定されたなら、$\eta_1=0.01$、$\eta_2=0.024$ とテーブルに入れる。
周波数応答関数（FRF）の定式化
構造減衰を含むFRFの式を教えてください。共振点での振幅がどう変わるか知りたいです。
1自由度系の構造減衰入りFRF（レセプタンス）はこうなる：
ここで $r = \omega/\omega_n$（振動数比）。共振点 $r=1$ では：
つまり共振振幅は $1/(K\eta)$ で、損失係数に反比例する。$\eta$ を1桁間違えると共振振幅が10倍になるから、設定値の桁は絶対に確認すること。自動車のステアリングコラム振動解析で $\eta=0.02$ のところを $\eta=0.002$ と入力して、FRFピークが10倍になって大騒ぎ — 実際にあったケースだよ。
レイリー減衰への変換手法
時刻歴応答解析では構造減衰が使えないんですよね。レイリー減衰に変換するとき、どうやって $\alpha$ と $\beta$ を決めるんですか？
レイリー減衰は $C = \alpha M + \beta K$ で、減衰比は：
構造減衰の $\zeta = \eta/2$ が一定なので、2つの周波数 $\omega_1, \omega_2$ で $\zeta(\omega_i) = \eta/2$ を満たすように連立方程式を解く：
例えば$\eta = 0.02$ で、関心帯域が50〜300 Hz（$\omega_1 = 100\pi, \omega_2 = 600\pi$）なら、この2点で $\zeta = 1\%$ になるように $\alpha, \beta$ を計算する。注意すべきは、この2点の間は減衰が過小に、外側は過大になること。帯域の選び方で結果がかなり変わるから慎重に。
因果律の制約と時間領域での扱い
構造減衰が時間領域で使えない理由を、もう少し詳しく教えてください。「因果律が破れる」って具体的にどういうことですか？
複素剛性 $K(1+i\eta)$ をフーリエ逆変換して時間領域に戻すと、インパルス応答関数 $h(t)$ が $t < 0$ でゼロにならない。つまり入力がまだ来ていないのに応答が出る — 因果律（原因→結果の時間順序）が崩れるんだ。
これは損失係数 $\eta$ が全周波数で一定という仮定自体が、Kramers-Kronig関係（因果律を満たすための実部と虚部の拘束条件）を破っているため。物理的に厳密な減衰モデルは必ず周波数依存性を持つ。構造減衰は「周波数領域での近似モデル」として割り切って使うべきだ。
半値幅法は60年以上使われ続けている
半値幅法（Half-Power Bandwidth Method）は、FRFピークの $1/\sqrt{2}$ 振幅の2周波数から損失係数 $\eta = \Delta f / f_n \approx 2\zeta$ を算出する手法。D.J. Ewins が1960年代に体系化して以来、現在も振動試験の標準手法として使われている。鉄鋼材料の $\eta$ は約0.001〜0.01、CFRP複合材は0.005〜0.02、制振鋼板は0.05〜0.2と材料によって2桁以上の差がある。
実践ガイド
材料・構造別の典型的な損失係数
損失係数 $\eta$ の値ってどのくらいが普通ですか？ 解析で入力するとき、何を参考にすればいいですか？
代表的な値をまとめておくよ。まずは材料単体の損失係数から：
| 材料 | 損失係数 $\eta$ | 等価減衰比 $\zeta$ | 備考 |
|---|---|---|---|
| 純アルミニウム | 0.0001〜0.001 | 0.005〜0.05% | 純度・結晶構造に依存 |
| アルミ合金（A6061等） | 0.001〜0.005 | 0.05〜0.25% | 航空機構造で一般的 |
| 軟鋼（SS400等） | 0.001〜0.006 | 0.05〜0.3% | 材料単体の値 |
| 鋼構造（溶接） | 0.005〜0.015 | 0.25〜0.75% | 接合部の摩擦を含む |
| 鋼構造（ボルト締結） | 0.01〜0.05 | 0.5〜2.5% | ボルト面の微小すべりが支配的 |
| CFRP積層板 | 0.005〜0.02 | 0.25〜1% | 積層構成・マトリクスに依存 |
| コンクリート | 0.02〜0.1 | 1〜5% | ひび割れ状態で大きく変化 |
| ゴム（天然ゴム） | 0.1〜0.5 | 5〜25% | 温度・周波数に強く依存 |
| 制振ゴム・ブチルゴム | 0.3〜1.5 | 15〜75% | 制振材料として使用 |
| 制振鋼板（サンドイッチ） | 0.05〜0.3 | 2.5〜15% | 拘束型制振材 |
ボルト締結の構造は材料単体より減衰が大きいんですね。なぜですか？
いいところに気づいたね。Beards（1983）の研究によると、機械構造物の全減衰の60〜90%は接合部の摩擦に起因する。ボルト締結面では振動時に微小すべり（マイクロスリップ）が発生し、クーロン摩擦でエネルギーが散逸する。だから同じ鋼材でも、溶接構造（$\eta \approx 0.01$）よりボルト締結構造（$\eta \approx 0.03$）のほうが3倍近く減衰が大きい。FEMで材料減衰だけ入力して実測と合わないのは大抵このせいだよ。
自動車NVH解析での適用事例
自動車のNVH解析で構造減衰を使う場合、具体的にどう設定するんですか？
自動車の車体（BIW: Body In White）の周波数応答解析は構造減衰の代表的な適用事例だ。典型的な設定をまとめると：
- 鋼板ボディパネル：$\eta = 0.01$〜$0.02$（スポット溶接部の摩擦を含む）
- ダッシュボード制振材貼付部：$\eta = 0.02$〜$0.05$（拘束型制振材の効果）
- ゴムマウント（エンジン/サブフレーム）：$\eta = 0.1$〜$0.3$
- ガラス：$\eta = 0.002$〜$0.005$
NastranのSOL 111でモード法解析し、MAT1のGEフィールドで材料ごとに $\eta$ を設定するのが標準的なやり方だ。全体にPARAM, Gで一律の値を入れるのは粗い方法で、精度が必要なら材料別設定は必須だよ。
航空宇宙・プラント配管での適用
航空宇宙ではどうですか？ 自動車とは違いますか？
航空宇宙では減衰が小さいぶん、設定値がさらにクリティカルだ。具体例を挙げると：
- 人工衛星のパネル（CFRP/アルミハニカム）：$\eta = 0.005$〜$0.01$。ランチャーの振動環境解析で使う。$\eta$ を過大に設定するとピーク応力を過小評価してしまい、軌道上で破損する危険がある
- ジェットエンジンのブレード：$\eta = 0.001$〜$0.003$。高サイクル疲労の評価には共振振幅の正確な予測が必須
- プラント配管系：$\eta = 0.01$〜$0.03$（サポート・クランプの摩擦を含む）。ASME/JSME規格では設計用減衰比としてパイプ系に $\zeta = 1$〜$2\%$（$\eta = 0.02$〜$0.04$）が推奨されている
実務チェックリスト
構造減衰の設定でチェックすべき項目をまとめたよ：
- 解析タイプの確認：周波数応答解析（SOL 108/111, *STEADY STATE DYNAMICS, HROPT）であること。時刻歴解析では使えない
- 損失係数の桁：金属構造で $\eta > 0.1$ になっていたら入力ミスを疑う。ゴム以外で $\eta > 0.05$ は稀
- 材料別 vs 全体：異種材料（鋼＋ゴム等）では材料ごとに $\eta$ を設定。PARAM,G の一律値で済ませない
- GEとPARAM,Gの関係：NastranではMAT1のGEフィールドの値にPARAM,Gが加算される（$\eta_{\text{total}} = \text{GE} + G$）。二重定義に注意
- 実験との突合：可能なら半値幅法やモード解析の結果と比較して $\eta$ を検証する
- $\eta = 2\zeta$ の変換：モード解析の減衰比 $\zeta$ から変換する場合、$\eta = 2\zeta$ を使う（共振点での等価関係）
超高層ビルの減衰比は高さとともに低下する
高層建築物の減衰比は階数の増加とともに低下する傾向がある。ATC-72（2010年）によると、20階建てで約2%、50階建てで約1%、100階超では0.5%前後。これはジョイント（接合部）の数に対する構造全体の質量の比率が変化するためと考えられている。東京スカイツリー（634m）では内筒と外筒の間に摩擦型制振壁を設置し、等価的に損失係数 $\eta \approx 0.1$（$\zeta \approx 5\%$）を実現して強風時の居住性基準を満たしている。
ソフトウェア比較
Nastranでの設定
Nastranで構造減衰を設定する方法を具体的に教えてください。SOL 111で使います。
Nastranには3つの設定方法がある。使い分けが重要だ。
1. 全体に一律の構造減衰（PARAM, G）
$ 全体構造減衰 η=2%
PARAM, G, 0.02
$ W3：直接法SOL 108 で等価粘性減衰に変換する基準周波数
PARAM, W3, 628.3 $ = 2π×100 Hz2. 材料ごとの構造減衰（MAT1のGEフィールド）
$ 鋼材 η=1%
MAT1, 1, 210000., , 0.3, 7.85E-9, , , 0.01
$ ゴム η=30%
MAT1, 2, 10., , 0.49, 1.1E-9, , , 0.33. 周波数依存テーブル（TABDMP1）
$ SOL 111 用：モード減衰をG（構造減衰）形式で定義
TABDMP1, 100, G
+, 10.0, 0.01, 100.0, 0.015, 500.0, 0.02, ENDT注意：GEフィールドとPARAM,Gは加算される。MAT1のGE=0.01とPARAM,G=0.02を同時定義すると、その材料は $\eta=0.03$ になる。意図しない二重定義に要注意だ。
Abaqusでの設定
Abaqusだとどう設定しますか？ Nastranとの違いはありますか？
Abaqusでは*DAMPINGキーワードで設定する。
材料レベルの構造減衰
*MATERIAL, NAME=STEEL
*ELASTIC
210000., 0.3
*DAMPING, STRUCTURAL=0.01*STEADY STATE DYNAMICS解析ステップ
*STEP
*STEADY STATE DYNAMICS, DIRECT
20., 500., 1, 1.0
** 周波数範囲20-500Hz、1Hz刻みNastranとの大きな違いは、Abaqusには「全体に一律の構造減衰」を1パラメータで設定する機能がないこと。材料ごとに *DAMPING, STRUCTURAL を定義する必要がある。また Abaqus の *MODAL DAMPING では構造減衰を直接指定できず、粘性の減衰比 $\zeta$ で入力して内部で変換する。
Ansysでの設定
Ansysではどうですか？ DMPSTRコマンドで設定するんですよね？
Ansys Mechanical APDLでの設定方法を整理するよ。
全体構造減衰
! 全体に η=2% の構造減衰
DMPSTR, 0.02材料ごとの構造減衰（MP,DMPS）
! 材料1（鋼）：η=1%
MP, DMPS, 1, 0.01
! 材料2（ゴム）：η=30%
MP, DMPS, 2, 0.3Harmonic解析ステップの定義
/SOLU
ANTYPE, HARMIC
HROPT, MSUP ! モード重ね合わせ法
HARFRQ, 20, 500 ! 20-500 Hz
NSUBST, 480 ! サブステップ数Ansys Workbenchでは、Engineering Data の材料プロパティで "Structural Damping Coefficient" を設定する GUI がある。APDL と違って直感的に操作できるが、内部的には同じ MP,DMPS コマンドが生成される。
3ソルバー比較表
| 機能 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 全体構造減衰 | PARAM, G | （機能なし、材料別のみ） | DMPSTR |
| 材料ごと | MAT1 GEフィールド | *DAMPING, STRUCTURAL | MP, DMPS |
| モードごと | TABDMP1（G形式） | *MODAL DAMPING（ζで入力） | MDAMP |
| 周波数依存テーブル | TABDMP1 | *DAMPING, DEPENDENT | TB, SDAMP |
| GE + G の加算 | あり（要注意） | N/A | あり |
| 直接法での利用 | SOL 108 | *STEADY STATE DYNAMICS, DIRECT | HROPT, FULL |
| モード法での利用 | SOL 111 | *STEADY STATE DYNAMICS, MODAL | HROPT, MSUP |
結局、どのソルバーを使うのが一番楽ですか？
構造減衰の柔軟性ではNastranが最も優れている。PARAM,G で全体設定、MAT1 GE で材料別、TABDMP1 でモード別・周波数依存テーブルと3段階の制御ができる。航空宇宙業界でNastranが標準なのは、この減衰モデリングの充実度も理由の一つだよ。Abaqusは非線形解析には強いが、線形周波数応答での構造減衰の設定はやや不便。Ansysはバランスが良いが、APDL のコマンド体系に慣れが要る。
NastranのPARAM,Gは「加算」される落とし穴
NastranのPARAM,Gは全要素に一律に加算される全体構造減衰であり、MAT1カードのGEフィールド（材料構造減衰）とは独立して加算される。つまり PARAM,G=0.02 かつ MAT1のGE=0.01 と設定すると、その材料の実効損失係数は $\eta = 0.03$ になる。これを知らずに「GEで材料別に設定したのに全体が大きくなっている」と悩むエンジニアは多い。解決策は PARAM,G=0.0 にしてGEだけで制御するか、GEを0にしてPARAM,Gだけで制御するか、明確に使い分けること。
トラブルシューティング
構造減衰を時間領域で使ってしまった
先生、Nastranの過渡応答解析（SOL 112）にPARAM,Gを設定したら、応答が入力前から始まっている異常な結果になりました。何が起きたんですか？
典型的な初心者ミスだね。構造減衰は周波数領域専用のモデルだ。時間領域で使うと因果律が崩れて「未来の入力に対する応答」が出てしまう。NastranのSOL 109/112では PARAM,G を設定しても実質的に無視されるか、不正確な結果になる。
対策：時間領域ではレイリー減衰に変換する。$\eta = 0.02$ なら $\zeta = 0.01$ で、関心帯域の2点で $\zeta(\omega) = 0.01$ を満たすように $\alpha, \beta$ を算出する。
$ SOL 112 用：レイリー減衰で代替
$ 50-300 Hz 帯域で ζ≈1% になるように設定
PARAM, ALPHA1, 5.984 $ 質量比例減衰
PARAM, ALPHA2, 2.274E-6 $ 剛性比例減衰損失係数の桁を間違えた
FRFの共振ピークが実験の10倍の振幅になってしまいました。何が原因ですか？
共振振幅は $|H| = 1/(K\eta)$ で損失係数に反比例するから、$\eta$ を1桁小さく入力すると振幅が10倍になる。よくあるパターンは：
- $\zeta$ と $\eta$ の混同：実験で得た $\zeta = 1\%$ をそのまま $\eta = 0.01$ と入力してしまう。正しくは $\eta = 2\zeta = 0.02$
- 百分率と小数の混同：$\eta = 2\%$ のつもりで $\eta = 2.0$ と入力。実質200%の減衰で、ピークが消えてしまう
- Abaqusでの*DAMPING入力ミス：*MATERIAL の STRUCTURAL DAMPING ではなく MATERIAL DAMPING（レイリー減衰）を設定してしまう。パラメータの意味が全く異なる
2019年に国内プラントメーカーで配管系解析において $\eta = 0.01$ を $\eta = 0.001$ と設定し、疲労寿命を100倍過小評価してしまった事例がある。桁の確認は解析レビューの最重要項目だ。
材料GEとPARAM,Gの二重定義
Nastranで各材料のGEフィールドに $\eta$ を設定したのに、全体的に減衰が大きすぎます。共振ピークが低すぎる感じで…
PARAM, G が残っていないか確認してみて。NastranではPARAM,GとMAT1のGEが加算される。例えばPARAM,G=0.02 でMAT1のGE=0.02 だと、実効 $\eta = 0.04$ になる。意図の倍の減衰が入っている状態だ。
対策：材料別にGEを設定する場合は、必ず PARAM,G=0.0 にするか、PARAM,Gの行を削除する。逆に全体一律で十分なら、MAT1のGEフィールドは空白（ゼロ）にしてPARAM,Gだけで制御する。
FRFが実験と合わない
構造減衰の値は文献値を使ったのに、FRFのピーク振幅が実験より大きいんです。共振周波数は合っているんですが。
共振周波数が合っていてピーク振幅だけ違うなら、減衰の過小評価が原因だ。考えられるのは：
- 接合部減衰の未考慮：FEMでは材料の $\eta$ しか入れていないが、実構造のボルト結合面・溶接部の摩擦による減衰が全体の60〜90%を占める。材料の $\eta = 0.005$ に加えて接合部効果で実効 $\eta = 0.02$ 程度になることが多い
- 空気抵抗の影響：薄板パネルの面外振動では空気の粘性抵抗が無視できない減衰を提供する。特に低周波の面外1次モード
- 境界条件の減衰：実験のジグ固定部が完全拘束ではなく、微小なすべりで減衰を生んでいる
対策：まず実験のFRFから半値幅法で各モードの $\eta$ を実測し、それをFEMに入力する。文献値はあくまで出発点で、最終的には実験同定値を使うのが正道だよ。
減衰の1桁ミスは疲労寿命の100倍ミスにつながる
共振点での応力振幅は $\sigma \propto 1/\eta$ に比例する。S-N曲線の疲労寿命は一般に $N \propto \sigma^{-m}$（$m \approx 3$〜$5$）で応力のべき乗に反比例するため、$\eta$ を1/10に間違えると応力が10倍、疲労寿命は $10^3$〜$10^5$ 倍も変わってしまう。つまり「10年持つ」つもりの部品が「数日で破損する」計算になる（またはその逆）。減衰値のレビューは構造物の安全性に直結する最重要チェック項目だ。
関連トピック
なった
詳しく
報告