Miles方程式による簡易評価
Miles方程式による簡易評価の理論基礎
Miles方程式とは
先生、Miles方程式って何ですか？
John Miles（1954）が導出した1自由度系のランダム振動の簡易公式だ。入力PSDが一定（ホワイトノイズ近似）の場合、応答のRMS値を1行の式で求められる。
または加速度RMS：
ここで $Q = 1/(2\zeta)$ は共振の品質係数、$S_{\ddot{u}}(f_n)$ は固有振動数での入力加速度PSD。
FEMなしで1行で計算できる！
Miles方程式はスクリーニング評価に最適。FEMのPSD解析の前に概算を出し、オーダーを把握する。
仮定と限界
Miles方程式の仮定：
1. 1自由度系 — 多自由度系には直接適用不可
2. 入力PSDが固有振動数付近で一定 — ホワイトノイズ近似
3. 小減衰 — $\zeta < 0.1$ 程度
入力PSDが一定でない場合はどうなりますか？
固有振動数 $f_n$ での入力PSD値を使えば、多くの場合で10〜20%の精度が得られる。入力PSDが $f_n$ 付近で急変する場合は不正確。
実務での使い方
Miles方程式の実務的な使い方：
1. 機器の1次固有振動数 $f_n$ を推定
2. 振動環境のPSD値 $S(f_n)$ を読み取る
3. 減衰比 $\zeta$ を仮定（通常 $\zeta = 0.02 \sim 0.05$）
4. RMS応答を計算
5. 3σ（3×RMS）で最大応答を推定
6. 許容値と比較
FEMの代わりに使える場面があるんですね。
概念設計段階のスクリーニングには十分。詳細設計ではFEMのPSD解析に進む。
まとめ
要点：
- $a_{rms} = \sqrt{\pi f_n Q S(f_n) / 2}$ — 1行の公式
- スクリーニング評価に最適 — FEM前の概算
- ホワイトノイズ近似 — 入力PSDが$f_n$付近で一定
- 3σで最大応答 — 設計値として使用
- 多自由度系には直接適用不可 — 各モードに個別適用は可能
Miles方程式の誕生秘話
John W. Milesが1954年にJournal of the Aeronautical Sciencesへ発表した論文「On Structural Fatigue Under Random Loading」が原点。当時の米空軍が航空機疲労破壊の予測に苦しんでいた時代背景があり、白色雑音近似＋1自由度系という大胆な仮定のもと、応答RMS値をわずか3パラメータで算出できる画期的な式が生まれた。
Miles方程式による簡易評価の数値計算手法
Miles方程式の計算例
具体的な計算例を見せてください。
電子機器のプリント基板（PCB）の振動評価：
- $f_n = 200$ Hz（PCBの1次固有振動数）
- $S_{\ddot{u}} = 0.04$ g²/Hz（MIL-STD-810の入力PSD）
- $Q = 20$（$\zeta = 2.5\%$）
67 Gの加速度！ PCB上のBGAはんだが持つかどうか心配ですね。
BGAの典型的な耐衝撃加速度は50〜100 G。Miles方程式で67 Gと出たら、FEMで詳細評価すべき。PCBの補強やマウント変更を検討する。
多モード系への拡張
多自由度系では各モードにMilesを適用し、SRSS（Square Root of Sum of Squares）で合成：
各モードの寄与を二乗和の平方根で合成。応答スペクトル法のSRSSと同じですね。
モードが十分離れていれば（$f_{i+1}/f_i > 1.2$ 程度）SRSS合成は正確。密集モードではCQC（Complete Quadratic Combination）が必要。
まとめ
Miles式の3ステップ計算法
Miles方程式の適用手順は①固有振動数fn（Hz）の確認、②入力PSDのfnにおける値G²/Hz（W(fn)）の読み取り、③応答RMS＝√(π/2 × fn × Q × W(fn))の計算、の3ステップ。Q値（≈1/(2ζ)）は構造減衰比ζから決まり、航空宇宙構造では通常Q＝10（ζ＝5%）が仮定される。計算はExcelでも10秒で完了する。
Miles方程式による簡易評価の実務適用
Miles方程式の実務
電子機器、宇宙機器、軍用機器のランダム振動評価でスクリーニングとして広く使われている。
感度パラメータ
Miles方程式で応答に最も影響するパラメータ：
| パラメータ | 応答への影響 | 備考 |
|---|---|---|
| $Q$（品質係数） | $a_{rms} \propto \sqrt{Q}$ | $Q$が2倍→RMSが$\sqrt{2}$倍 |
| $f_n$（固有振動数） | $a_{rms} \propto \sqrt{f_n}$ | 高い$f_n$→大きいRMS |
| $S(f_n)$（入力PSD） | $a_{rms} \propto \sqrt{S}$ | 入力が2倍→RMSが$\sqrt{2}$倍 |
$Q$（減衰の逆数）が最も不確かなパラメータですね。
$Q = 10$ と $Q = 50$ で応答が$\sqrt{5} \approx 2.2$倍変わる。減衰の推定精度がMiles方程式の精度を支配する。
実務チェックリスト
Miles方程式はFEMの「サニティチェック」としても使えますね。
FEMのPSD解析の結果がMilesの概算と桁で一致するか確認する。大きくずれていたらFEMの設定ミスを疑う。
ロケット搭載機器への適用実例
SpaceX Falcon 9の搭載電子機器設計では、打上げ時の音響・振動環境としてMIL-STD-810G Method 514.6の振動プロファイルが使われる。典型的な入力PSDは20〜2000Hzで0.04 G²/Hz、設計担当者はMiles方程式でRMS加速度を算出し3σ値（＝3×RMS）を設計荷重として疲労評価に用いる。実績ある宇宙機器メーカーJPL（NASAジェット推進研究所）でも標準手法として採用。
Miles方程式による簡易評価のソフトウェア比較
Miles方程式のツール
Miles方程式は手計算の公式だから、特別なツールは不要。
| ツール | 特徴 |
|---|---|
| 手計算（電卓） | 最速。1分で結果 |
| Excel/Python | パラメトリック評価。$Q$の感度分析 |
| Vibration Research社のアプリ | Miles方程式＋PSD解析の統合 |
| FEMソルバーのPSD解析 | 詳細評価。Miles方程式はスクリーニング用 |
選定ガイド
Nastranのランダム振動解析実装
MSC NastranのSOL 111（Modal Frequency Response）ではランダム振動解析をCASECC制御文RANDOM＋PSD入力カードDLOADで実行。出力はRMSオプションのPESEFまたはPSDFを指定する。一方、Ansys Mechanical（Workbench 2023 R2以降）ではランダム振動解析をGUI上で3クリック設定でき、Miles方程式との比較検証ツールも内蔵。SIMcenter Nastranは出力フォーマットに互換性があり大規模モデルでの高速化に優れる。
Miles方程式による簡易評価の先端研究
Miles方程式の先端トピック
Miles方程式の精度限界と補正
Miles方程式は入力PSDが対象固有振動数周辺で平坦（白色雑音近似）である場合にのみ厳密に成立する。実際のロケット振動環境では低周波側で急峻な立ち上がりがあるため、非平坦PSDに対して精度が最大30%劣化することが1970年代のNASA TN D-8303で報告されている。補正法としてNarrowband修正係数や有限帯域幅積分法が提案されており、Ansysのランダム振動モジュールは後者を採用している。
Miles方程式による簡易評価のトラブル対応
Miles方程式のトラブル
Q値設定ミスによる過大評価
Miles方程式の適用で最も多いミスはQ値の誤設定。溶接構造物の実測ではζ＝2〜3%（Q≈17〜25）になることがあるが、保守的にζ＝5%（Q＝10）を一律適用するとRMS値が√(10/5)≈1.41倍、すなわち41%過大評価となる。過剰設計を避けるため、実機振動試験でのモード同定結果や、類似構造のデータベースを活用してQを適切に設定することが重要。
関連トピック
なった
詳しく
報告