多自由度系のランダム振動
多自由度系のランダム振動の理論基礎
多自由度系のランダム振動
先生、多自由度系のランダム振動はどう扱いますか？
多自由度系ではモード間の相関（クロスモード項）が重要。応答PSDは：
モードが十分離れていればSRSS合成で十分。密集モードではCQC（Complete Quadratic Combination）が必要。
多点入力
複数の支持点に相関のあるランダム入力：
$[S_{in}]$ はクロススペクトルマトリクス。対角=オートPSD、非対角=クロスPSD。
まとめ
要点：
- クロスモード項 — 密集モードで重要。CQCで合成
- クロススペクトルマトリクス — 多点入力の相関を記述
- SRSS vs. CQC — モードが離れていればSRSS、密集ならCQC
- FEMのPSD解析は自動的にクロスモード項を含む
多自由度ランダム振動の基礎理論
多自由度（MDOF）ランダム振動は、モーダルスーパーポジション法で各固有モードに分解し、それぞれのモーダル座標でのランダム応答を求めた後に重ね合わせる。モード間の相関（クロスPSD）を考慮するSRSS（Square Root of Sum of Squares）法とCQC（Complete Quadratic Combination）法があり、固有振動数が接近している場合（比率<10%）はCQCが必須。1970年代にE.L. Wilsonらがカリフォルニア大学バークレー校で定式化した。
多自由度系のランダム振動の数値計算手法
FEMでの多自由度PSD
NastranのSOL 111 + RANDOMは全モードの自己項＋クロス項を自動計算。
多点入力：
```
RANDPS, 1, 1, 1, 100, 0.0 $ 入力1のオートPSD
RANDPS, 1, 2, 1, 101, 0.5 $ 入力1-2のクロスPSD
RANDPS, 2, 2, 1, 102, 0.0 $ 入力2のオートPSD
```
相関が不明なら完全相関＋無相関の2ケースで範囲を評価。
まとめ
クロスPSD行列を用いた解法
MDOFランダム振動の厳密解法では、入力クロスパワースペクトル密度行列S_FF(ω)と周波数応答関数行列H(ω)から出力クロスPSD行列S_XX(ω)＝H(ω)S_FF(ω)H*(ω)^Tを計算し、応答分散を周波数積分で求める。入力点が複数ある場合（例：4点サスペンション）はこの定式化が必須。Pythonのscipyライブラリではsignal.coherenceやnumpy.fftを組み合わせて数値実装できる。
多自由度系のランダム振動の実務適用
実務チェックリスト
自動車シャシーへの実際の適用
自動車業界では路面入力を複数タイヤ（4点入力）のランダム荷重として扱い、シャシーの疲労耐久評価を行う。Ford社は1990年代にベルギーのポーベルーテストコースの路面PSDデータをベースにロード・シミュレーション・テーブル（LST）試験規格を整備した。解析ではMSC Nastranのランダム振動SOL 111で最大10,000自由度モデルを解き、ストレスRMSから疲労損傷量をRainflow計数法で評価する。
多自由度系のランダム振動のソフトウェア比較
ツール
選定ガイド
ソルバー別MDOF対応比較
MDOF ランダム振動の主要ソルバー比較：Ansys Mechanicalは大規模モデルでのGPU並列化に対応（RTX 4090使用で従来比8倍速）、MSC Nastranは50年以上の実績と豊富なNASAの検証事例、Abaqusはランダム振動に非対応（過渡陰解法＋時刻歴で代替）、SIMcenter NastranはDAKOTAと連携した確率論的設計最適化が可能。HPC環境ではNastranとAnSys Mechanicalがほぼ同等のスケーリング性能を持つ。
多自由度系のランダム振動の先端研究
先端研究
非線形MDOFへの等価線形化
ゴムブッシュや摩擦ダンパーを含む非線形MDOFランダム振動では、等価線形化法（Equivalent Linearization, EL法）が有効。1954年にBogoliubovが提案した確率的EL法は、非線形復元力を等価な線形剛性と等価減衰で置き換え、反復計算で収束させる。Simuliaの研究チームが2018年に発表した論文では、自動車サスペンション非線形モデルで応答RMSの誤差を5%以内に抑えたことを報告している。
多自由度系のランダム振動のトラブル対応
トラブル
モード打ち切りが招く過小評価
1997年打ち上げのNASAカッシーニ探査機では振動解析に200モード以上を使用したが、150モードで打ち切った初期モデルでは特定周波数帯の応答が18%過小評価された。多自由度ランダム解析では、対象周波数の2倍まで固有値を含めることがMIL-STD-810Gでも推奨されている。
関連トピック
なった
詳しく
報告