構造パラメータ
PSD入力
PSD種別
S-N パラメータ
理論メモ
RMS応力: $\sigma_{rms}= c_\sigma \cdot \sqrt{\int W_{resp}(f)df}$
ピーク応力: $\sigma_{peak}= 3\sigma_{rms}$
Steinberg: $D = f_n T \sum p_i / N(\sigma_i)$
振動疲労PSD応力計算 — Steinberg 3σ法とは
🧑🎓
「振動疲労」って、普通の疲労と何が違うんですか？ あと、PSDって何ですか？
🎓
ざっくり言うと、時間的にランダムで不規則な振動による疲労のことだね。PSD（パワースペクトル密度）は、そのランダム振動の強さを周波数ごとに分解して表した「レシピ」みたいなものだよ。例えば自動車の走行振動やロケットの打ち上げ振動は、特定の周波数だけじゃなく広い帯域で揺れるから、PSDで定義するんだ。このツールの「PSD種別」をホワイトノイズからピンクノイズに変えてみると、周波数による強さの違いがわかるよ。
🧑🎓
え、そうなんですか！じゃあ、そのPSDからどうやって「疲労」が計算できるんですか？ 上の「モーダル応力係数 cσ」ってのもよくわからないです。
🎓
いい質問だね。まずcσは、振動の加速度（g）を部品に発生する応力（MPa）に変換する「倍率」だ。実務では有限要素解析（FEA）のモーダル解析で求めることが多いよ。この値とPSDから、応力の実効値（RMS応力）が計算できる。疲労寿命は、このランダムな応力波形の中に含まれる様々な大きさの応力振幅を、全部足し合わせて評価するんだ。ツールで「減衰比 ζ」を大きくすると、共振が抑えられてRMS応力が下がるのが確認できるはずだよ。
🧑🎓
なるほど！で、「Steinberg 3σ法」って名前がすごいですけど、あの「3σ」は統計のやつですか？ どうやって計算してるんですか？
🎓
その通り、統計の正規分布（ガウス分布）の3σ（3シグマ）から来てるんだ。ランダム振動による応力はほぼ正規分布すると仮定する。すると、応力振幅の約68%は±1σ_rms以内、約27%は±2σ_rms、約4.5%だけが±3σ_rmsの大きな振幅を持つことになる。Steinberg法は、この3つのバンドに分けて、それぞれの発生回数と材料のS-N曲線から疲労損傷を足し上げる、実務で非常によく使われる簡便法なんだ。ツールで「S-N指数 b」を変えてみると、寿命が大きく変わることがわかるよ。材料の疲労特性の影響を実感できるはずだ。
物理モデルと主要な数式
入力PSD $W_{in}(f)$ が、固有振動数 $f_n$、減衰比 $\zeta$ の1自由度系（部品の主要な共振モード）に作用した時の応答PSD $W_{resp}(f)$ を計算し、その面積から応力の実効値（RMS）を求めます。
$$ \sigma_{rms}= c_{\sigma}\cdot \sqrt{\int_{f_1}^{f_2}W_{resp}(f) \, df}$$
ここで、$c_{\sigma}$: モーダル応力係数 [MPa/g]、$W_{resp}(f)$: 応力の応答PSD [MPa²/Hz]、$f_1$, $f_2$: 評価周波数帯域。積分範囲内のPSDの面積（パワー）の平方根がRMS値となります。
Steinberg 3σ法では、応力振幅が正規分布すると仮定し、RMS応力 $\sigma_{rms}$ を基準として3つの応力レベルを定義し、Miner則を用いて累積損傷率 $D$ を計算します。
$$ D = f_n T \left( \frac{p_1}{N(\sigma_1)}+ \frac{p_2}{N(\sigma_2)}+ \frac{p_3}{N(\sigma_3)}\right) $$
ここで、$f_n$: 固有振動数 [Hz]、$T$: 曝露時間 [s]、$p_1=0.683, p_2=0.271, p_3=0.045$は各応力レベル発生確率、$\sigma_1=1\sigma_{rms}, \sigma_2=2\sigma_{rms}, \sigma_3=3\sigma_{rms}$、$N(\sigma_i)$はS-N曲線から求める応力レベル$\sigma_i$での破断寿命。$D=1$で破損と予測します。S-N曲線は $N \sigma^b = N_{ref}\sigma_{ref}^b$ でモデル化されています。
実世界での応用
車載電子機器：エンジンや走行路からのランダム振動に曝されるECU（エンジンコントロールユニット）やセンサーの寿命予測に使用されます。実車の測定データをPSD化し、筐体や基板のはんだ接続部の疲労を評価します。
航空宇宙機器：航空機の搭載機器やロケットのペイロードは、打ち上げ・飛行中の極めて過酷なランダム振動環境に耐えなければなりません。仕様としてPSDが規定され、Steinberg法による予備設計検証が広く行われています。
産業用機械・プラント配管：ポンプやコンプレッサーなどの回転機械から発生する振動が、支持構造や接続配管に伝わり、振動疲労破壊を引き起こすことがあります。振動測定データをPSD解析し、危険個所の寿命を評価します。
民生電子機器：スマートフォンやデジタルカメラなど、輸送中や使用中の偶発的衝撃・振動に対する耐久性を評価する際の参考として利用されます。特に落下時の内部回路基板の応答評価に応用されます。
よくある誤解と注意点
この手法を使い始める際、いくつかハマりやすいポイントがあるよ。まず「PSD入力は加速度だけど、応力に変換するcσは本当に定数？」という点。cσ（モーダル応力係数）は、線形で比例関係が成り立つ範囲でしか一定だ。例えば、大きな変形で幾何非線形性が出たり、材料が塑性域に入ったりすると、cσは変わってしまう。実務では、想定される振動レベルでFEAを走らせ、応力と加速度の関係が本当に直線的か確認するのが鉄則だ。
次に「1自由度系（1モード）の近似で本当に大丈夫？」という根本的な疑問。確かに、部品の共振周波数が広く離れていて、支配的なモードが1つだけならこれでOK。しかし、例えば2つのモードが近接している（例えば100Hzと110Hz）と、それらが相互作用してRMS応力が単純な足し算以上に大きくなる「モーダルカップリング」が起きる。ツールで計算した寿命が極端に短い/長い時は、モードの重なりを疑ってみよう。
最後に「Steinberg法の3つのバンド分けは万能ではない」こと。材料のS-N曲線の傾き（指数b）が非常に大きいor小さい場合、あるいは非ガウス性が強い振動（衝撃を多く含むなど）では、精度が落ちる。あくまで「簡便で保守的な見積もり法」と心得て、重要な設計では、より精緻なDirlik法などの確率密度関数ベースの手法と結果を比較することをお勧めする。
関連する工学分野
この「振動疲労PSD応力計算」の考え方は、実は様々な工学分野に応用されているんだ。まず「音響疲労」だ。航空機の機体やエンジンカウルは、乱気流やエンジン騒音による広帯域の音圧変動（これもPSDで定義される）で疲労損傷を受ける。音圧をPSDで入力し、構造の音響応答を求めて疲労寿命を予測する流れは、振動疲労と数学的にほぼ同じだ。
次に「海上構造物の波浪荷重疲労」。石油プラットフォームや洋上風車の支柱は、不規則な波浪の力で揺れる。波浪のエネルギー周波数分布（波浪スペクトル）がPSDに相当し、これが構造物の固有振動数で応答を増幅させる。波浪データからPSDを作り、同じくモーダル応答解析とMiner則を組み合わせて耐用年数を評価する。
さらに「半導体の信頼性試験」にも繋がる。基板実装されたBGA（ボールグリッドアレイ）のようなパッケージは、輸送中のランダム振動ではんだ接合部が疲労破壊する。ここでは、はんだのクリープ挙動を考慮した修正S-N曲線（例えばCoffin-Manson則）と組み合わせて、Steinberg法が適用されることがある。このように、「広帯域不規則励振→線形応答→累積損傷評価」という枠組みは、分野を超えた共通言語になっているんだ。
発展的な学習のために
このツールの背後にある理論を深めたいなら、まずは「ランダム振動理論」と「疲労強度学」の2本柱を勉強しよう。最初のステップとして、フーリエ変換とパワースペクトルの関係（Wiener-Khinchinの定理）を理解すると、PSDが「時間信号のパワーの周波数分布」であることが腹落ちするよ。例えば、ホワイトノイズPSDが周波数によらず一定なのは、全ての周波数成分が同じ強さで含まれるから、とイメージできる。
次に、ツールで省略されている「PSDから応力の確率密度関数を求める」部分を学ぶと、Steinberg法の位置づけがわかる。より一般的な手法であるDirlik法や、シンプルなレインフロー計数法の基礎となる「レインフロー法」を調べてみよう。これらのアルゴリズムは、ランダム波形の中から疲労損傷を与える「サイクル」をどう抽出して数えるか、という核心に迫っている。
実務に直結する次のトピックとしては、「モーダル解析とモーダル座標系への変換」がおすすめだ。cσを求めるために行うFEAのモーダル解析では、複雑な構造を複数の1自由度系（モード）の重ね合わせとして表現する。この「モーダル座標」の考え方を理解すれば、なぜツールの計算が1自由度系を前提としているのか、そして複数モードをどう扱えばいいのかが見えてくる。まずは、質量・剛性・減衰行列から固有値とモード形状を求める、という線形代数の計算に触れてみるのが良い第一歩だ。