振動疲労PSD応力計算 戻る
振動解析

振動疲労PSD応力計算 — Steinberg 3σ法

PSD(パワースペクトル密度)入力から1自由度応答を計算し、RMS応力・ピーク応力・Steinberg法による疲労損傷・期待寿命をリアルタイムに求める。

構造パラメータ
固有振動数 fₙ (Hz)
Hz
減衰比 ζ (%)
%
モーダル応力係数 cσ (MPa/g)
MPa/g
PSD入力
PSD種別
PSDレベル W₀ (g²/Hz)
g²/Hz
周波数下限 f₁ (Hz)
Hz
周波数上限 f₂ (Hz)
Hz
S-N パラメータ
S-N指数 b
基準応力 σref (MPa)
MPa
基準寿命 Nref (×10⁶サイクル)
×10⁶
試験時間 T (hours)
h
計算結果
計算結果
RMS加速度 (g)
RMS応力 (MPa)
ピーク応力 3σ (MPa)
等価静荷重 (N)
疲労損傷 D
期待寿命 (h)
PSDスペクトル
PSD スペクトル
理論・主要公式
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」を変えてみると、寿命が大きく変わることがわかるよ。材料の疲労特性の影響を実感できるはずだ。

よくある質問

cσ は、対象部品の固有モードにおいて、単位加速度(1g)入力時に発生する応力値です。FEM解析のモーダル過渡応答解析や、実験モーダル解析で同定した値を入力してください。経験値や類似設計のデータで代用することも可能ですが、精度は低下します。
評価周波数帯域 f1〜f2 は、対象部品の固有振動数 fn を中心に、減衰比 ζ に応じて拡がる共振ピークを十分に含む範囲に設定してください。目安として fn ± 3fnζ 程度の範囲を推奨します。範囲が狭すぎるとRMS応力を過小評価する恐れがあります。
この手法は応力振幅が正規分布する前提に基づく簡易評価であり、実際のランダム振動疲労試験の平均的な傾向を推定します。ただし、非線形挙動や平均応力効果、材料のばらつきは考慮されないため、設計の一次スクリーニングとしてご利用ください。
一般的な電子部品実装基板では ζ=0.01〜0.05、金属構造物では ζ=0.001〜0.01程度が目安です。不明な場合は安全側として ζ=0.01(1%)を仮設定し、感度解析で複数の値を試すことを推奨します。減衰比が小さいほど応答が大きくなるため、過小評価を避けられます。

実世界での応用

車載電子機器:エンジンや走行路からのランダム振動に曝される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法などの確率密度関数ベースの手法と結果を比較することをお勧めする。

使い方ガイド

  1. PSD入力データ(CSV形式)を準備:周波数[Hz]、パワースペクトル密度[g²/Hz]の2列で記述。例:20Hz~2kHzの帯域で加速度PSDを取得
  2. 構造パラメータを入力:質量m[kg]、固有振動数fn[Hz]、減衰比ζ[%]。アルミ部品の場合fn=150Hz、ζ=3~5%が標準
  3. シミュレータが1自由度応答スペクトラムを計算し、RMS加速度・RMS応力・Steinberg 3σピーク応力を自動算出。疲労寿命曲線(Basquin則)から期待寿命を直接判定

具体的な計算例

PCB搭載アルミフレーム(質量0.5kg、E=70GPa、断面係数Z=5cm³)が振動環境(0~500Hz帯域、最大PSD密度0.08g²/Hz@100Hz)に曝露される場合:RMS加速度2.8g→RMS応力12.5MPa、Steinberg 3σ法でピーク応力37.5MPa。S-N曲線(k=7、応力基準値σ₀=120MPa@10⁷サイクル)適用により、推定疲労損傷D=0.18、期待寿命1200時間。

実務での注意点

  1. 減衰比ζの過小評価は危険:実測値がない場合、金属材料は3~8%、複合材料は5~12%を設定。ζ低下でピーク応力は3倍以上増加
  2. PSD統合範囲の確認が必須:共振周波数±3倍帯域外のデータは応答に寄与しないため、計算前に帯域制限を実施。例えば150Hzの部品なら50~450Hzを確認
  3. Steinberg係数(通常3σ)は非ガウス加速度に対して保守的。実装時に試験結果と照合し、係数2.5~3.5の範囲で補正検証を行う