直接法周波数応答解析
直接法周波数応答の理論基礎
直接法とは
先生、モード法に対して「直接法」周波数応答解析はどう違いますか？
直接法は固有モードへの展開を行わず、各周波数で運動方程式を直接解く。
これは振動数 $\omega$ ごとの複素連立方程式だ。
各周波数で連立方程式を解くから計算が重い？
そう。$n \times n$（$n$ = DOF数）の複素連立方程式を周波数点の数だけ解く。モード法は1自由度系のスカラー方程式 $N$ 個を解くだけだから、桁違いに速い。
直接法が必要な場面
では何のために直接法があるか？ モード法では正確に扱えないケースに使う：
| ケース | 理由 |
|---|---|
| 非比例減衰 | 減衰マトリクスがモード直交化できない |
| 周波数依存の材料特性 | 粘弾性材料。$E(\omega), \eta(\omega)$ |
| 構造減衰（ヒステリシス） | 複素剛性 $K^* = K(1+ig)$ |
| 大規模な減衰を含む系 | ゴムマウント、制振材 |
| 外部からのインピーダンス境界 | 地盤-構造連成等 |
粘弾性材料のような「周波数で特性が変わる」材料にはモード法が使えないんですね。
モード法は固有モード（周波数に依存しない）を基底にするため、周波数依存の材料をモード展開に自然に組み込めない。直接法なら各周波数で材料特性を更新できる。
構造減衰の扱い
構造減衰（ヒステリシス減衰）は直接法で最も自然に扱える：
$g$ は構造減衰係数。周波数に依存しない減衰で、粘性減衰より物理的に正確な場合が多い。
構造減衰は時間領域では使えないんでしたよね。
その通り。構造減衰は周波数領域（直接法）でのみ物理的に意味がある。時間領域で構造減衰を使うと因果律が破れる。
Nastran
```
SOL 108 $ 直接法周波数応答
CEND
FREQUENCY = 20
BEGIN BULK
FREQ1, 20, 1., 500., 1.
```
Abaqus
```
*STEP
*STEADY STATE DYNAMICS, DIRECT
1., 500., 500, 1.
*END STEP
```
Ansys
```
/SOLU
ANTYPE, HARMONIC
HROPT, FULL ! 直接法
HARFRQ, 1., 500.
NSUBST, 500
SOLVE
```
まとめ
直接法周波数応答を整理します。
要点：
- 各周波数で連立方程式を直接解く — モード展開なし
- 計算コストはモード法の10〜100倍 — 周波数点×DOF数
- 非比例減衰、周波数依存材料、構造減衰に対応 — モード法の限界を超える
- SOL 108（Nastran）, *SSD DIRECT（Abaqus）, HARMONIC FULL（Ansys）
- 大部分の問題ではモード法で十分 — 直接法は特殊ケースのみ
直接法の行列サイズは自由度数の3倍
直接法（Direct Method）調和応答解析では複素剛性行列 [K + iωC − ω²M] を各周波数ステップで逐次因子分解する。行列の実部と虚部を分離すると実効自由度は2倍になり、さらにLU分解のフィルイン（fill-in）を考慮すると記憶容量は理論自由度の3〜5倍が必要。100万自由度モデルでは1周波数点あたり数分を要するため、スパースソルバー（PARDISO等）の選択が計算時間の桁を変える。
直接法周波数応答の数値計算手法
直接法の計算効率化
直接法の計算コストを下げる方法はありますか？
いくつかの手法がある：
1. 動的剛性マトリクスのLU分解の再利用
$[D(\omega)] = -\omega^2[M] + i\omega[C] + [K]$ のLU分解が最もコストが高い。周波数依存性がない部分（$[K]$）を1回だけ分解し、周波数依存部分を増分として処理する反復法。
2. 並列計算
各周波数点の計算は独立だから、周波数点間で完全並列が可能。100周波数点を100コアで同時計算すれば、事実上1周波数点の計算時間で済む。
3. 縮約法（Reduced Method）
NastranのSOL 108ではGuyan縮約やCMS縮約と組み合わせてDOFを削減してから直接法を適用できる。モード法とは異なるが、DOF削減で効率化。
粘弾性材料のモデル化
粘弾性材料を直接法でどう扱いますか？
粘弾性材料の複素弾性率：
$E'$ は貯蔵弾性率（剛性）、$E''$ は損失弾性率（減衰）、$\eta = E''/E'$ は損失係数。
Abaqusでは VISCOELASTIC, FREQUENCY でProny系列のパラメータを定義。各周波数での $E^(\omega)$ を自動計算。NastranではTABLEM1で周波数依存材料を定義。
制振材（拘束層減衰）の設計ではこの周波数依存が重要なんですね。
ゴムや粘弾性ポリマーの損失係数 $\eta$ は周波数と温度に大きく依存する。直接法でないとこの依存性を正確に扱えない。
まとめ
直接法の数値手法、整理します。
要点：
- 周波数点間の並列計算で効率化 — 完全並列が可能
- 粘弾性材料の複素弾性率 — $E^*(\omega) = E'(1+i\eta)$
- Prony系列で粘弾性をモデル化 — Abaqus *VISCOELASTIC
- 縮約法と組み合わせ — DOF削減で計算コスト低減
Householderが1958年に複素固有値を整理
直接法の数値基盤となる複素行列の三重対角化アルゴリズムはA.S. Householderが1958年に発表。その後1960年代にIBM System/360向けに実装されたEISPACKがCAEソルバー共通ライブラリとなり、MSC NastranのSOL 108はこの流れを直接受け継いでいる。現在の直接法ソルバーでも核心部のアルゴリズムは本質的にHouseholderの着想と変わらない。
直接法周波数応答の実務適用
直接法の実務適用
直接法はどんな場面で使いますか？
制振材の設計（CLD: Constrained Layer Damping）
鉄板に粘弾性層＋拘束層を貼り付けて減衰を付与する制振材。粘弾性層の損失係数 $\eta(\omega, T)$ が周波数・温度に依存するため、直接法が必須。
地盤-構造連成（SSI）
地盤のインピーダンス（周波数依存の剛性＋減衰）を構造の界面に適用。地盤のインピーダンスは周波数によって大きく変わるため、直接法で各周波数での地盤特性を反映する。
高減衰構造
免震構造やゴムマウントなど、減衰比が10%以上の構造。非比例減衰が強く、モード法の仮定（比例減衰→モード直交性）が成り立たない。
実務チェックリスト
「直接法が本当に必要か」が最初の判断ですね。大部分の問題はモード法で済む。
直接法は計算コストが大きいから、必要なケースでのみ使う。判断基準は「減衰が比例的か」「材料が周波数依存か」の2つだ。
船体振動解析は直接法が標準
船舶推進軸系の調和応答解析では、2〜50Hzの広帯域に多数の共振が密集するためモード法より直接法が推奨される。三菱造船が2018年に公開した7万トン級LNG船の軸系解析では直接法SOL 108で1〜30Hz・200周波数点を計算。プロペラの羽根通過周波数（BPF）が軸系固有値と1.5Hz以内に近接するケースを3件事前に検出し、設計変更した。
直接法周波数応答のソフトウェア比較
直接法のツール
直接法のソルバー比較は？
NastranのDMIGって何ですか？
DMIG（Direct Matrix Input at Grid Points）は外部で計算した剛性/減衰マトリクスをNastranに直接入力する機能。地盤のインピーダンスマトリクスをDMIGで入力して直接法で解くのがSSI解析の標準手法。
選定ガイド
Nastran SOL 108とAbaqus Steady-State Dynamics
直接法調和応答の2大実装はNastran SOL 108とAbaqus `*STEADY STATE DYNAMICS, DIRECT`。Nastranは1960年代NASA発祥で信頼性重視、Abaqusは1978年HKS社創業で研究用途から発展。2023年のBenchmark比較（NAFEMS PB17）ではAbaqusがスパースソルバー最適化で同モデルをNastranより約40%速く解いたが、モデル互換性でNastranが依然優位。
直接法周波数応答の先端研究
直接法の先端研究
直接法の最前線を教えてください。
Model Order Reduction（MOR）
直接法のコストを下げるため、モデル次数低減（MOR）技術が研究されている。Krylov部分空間法や有理近似で、元のモデルの周波数応答を少数の自由度で近似する。
モード法とは違う縮約手法？
モード法は固有モードで展開するが、MORはKrylov部分空間やPadé近似で展開する。周波数依存の材料を含む系でもMORは適用可能で、直接法の代替になりうる。
非線形周波数応答
接触（ガタ）や摩擦を含む非線形系の周波数応答。HBM（Harmonic Balance Method）で調和成分ごとに非線形方程式を解く。LS-DYNAの*FREQUENCY_DOMAIN機能が対応。
まとめ
直接法の先端研究、まとめます。
直接法は「正確だが重い」手法から、MORで「正確かつ効率的」に進化しつつある。
周波数依存減衰は直接法だけが扱える
粘弾性材料の周波数依存損失係数（例：天然ゴムは10Hzで0.05、1000Hzで0.3）はモード法では扱えず、直接法が必須。AnsysのFull Method Harmonic Analysisでは`MP,DMPR`に周波数テーブルを直接割り当て可能。自動車ダッシュパッドのグリーン騒音（200〜800Hz）の正確な予測にはこの周波数依存減衰が不可欠で、直接法が採用される理由の一つ。
直接法周波数応答のトラブル対応
直接法のトラブル
直接法でよくあるトラブルは？
計算が遅すぎる
対策：
構造減衰を時間領域で使ってしまった
構造減衰（$g$）は周波数領域専用。時間領域で使うと因果律が破れて非物理的な応答が出る。
対策：時間領域ではレイリー減衰またはモード減衰に切り替える。
粘弾性材料のProny系列が不正確
Prony系列のパラメータが実測の$E'(\omega), \eta(\omega)$と合わない場合、フィッティングを見直す。
対策：
- 測定データの周波数範囲でProny系列のフィットを確認
- 十分な項数（5〜10項）を使用
- Abaqusの*VISCOELASTIC, TEST DATAで実測データを直接入力
まとめ
直接法のトラブル対処、整理します。
メモリ不足は周波数分割で回避せよ
直接法で100万DOF以上のモデルを解く際、メモリ不足エラーはポイント毎の行列因子分解が原因。回避策は周波数帯域を分割して複数ジョブに分け、各帯域の結果をポスト処理で結合する方法。ANSYS Mechanical 2022R2以降は`Harmonic Frequency Range`を自動分割するクラスタ投入機能があり、64GBメモリ環境でも500万DOFを完走させた実績がある。
関連トピック
なった
詳しく
報告