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パラメータで算出できる画期的な式が生まれた。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
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秒で完了する。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
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ジェット推進研究所）でも標準手法として採用。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
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は出力フォーマットに互換性があり大規模モデルでの高速化に優れる。
選定で最も重要な3つの問い
- 「何を解くか」：Miles方程式による簡易評価に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
Miles方程式の先端トピック
Miles方程式の精度限界と補正
Miles方程式は入力PSDが対象固有振動数周辺で平坦（白色雑音近似）である場合にのみ厳密に成立する。実際のロケット振動環境では低周波側で急峻な立ち上がりがあるため、非平坦PSDに対して精度が最大30%劣化することが1970年代のNASA TN D-8303で報告されている。補正法としてNarrowband修正係数や有限帯域幅積分法が提案されており、Ansysのランダム振動モジュールは後者を採用している。
トラブルシューティング
Miles方程式のトラブル
Q値設定ミスによる過大評価
Miles方程式の適用で最も多いミスはQ値の誤設定。溶接構造物の実測ではζ＝2〜3%（Q≈17〜25）になることがあるが、保守的にζ＝5%（Q＝10）を一律適用するとRMS値が√(10/5)≈1.41倍、すなわち41%過大評価となる。過剰設計を避けるため、実機振動試験でのモード同定結果や、類似構造のデータベースを活用してQを適切に設定することが重要。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——Miles方程式による簡易評価の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告