変動荷重スペクトル疲労累積損傷 戻る
破壊力学

変動荷重スペクトル疲労累積損傷計算

S-N曲線とMiner則を用いた変動振幅疲労解析。鋼・アルミ・CFRPの材料パラメータ、Goodman平均応力補正、残存寿命と安全率をリアルタイム計算。

材料 & 荷重スペクトル
材料
荷重スペクトル種別
応力振幅 S₁ (MPa)200
繰返し数 n₁ (×10⁶)1.0
応力振幅 S₂ (MPa)120
繰返し数 n₂ (×10⁶)5.0
平均応力 σₘ (MPa)50
引張強さ Sᵤ (MPa)600
稼働周波数 f (Hz)10

理論メモ

S-N曲線: $N = \left(\frac{\sigma_f'}{S_a}\right)^{1/b}$
Miner則: $D = \sum_i \frac{n_i}{N_i}$
Goodman補正: $S_{a,eq}= \dfrac{S_a}{1 - \sigma_m / S_u}$
計算結果
許容繰返し数 N₁
許容繰返し数 N₂
部分損傷 D₁
部分損傷 D₂
累積損傷度 D
安全率 1/D
推定残存寿命(年)

累積損傷度メーター(D = 1 で破損)

変動荷重スペクトル疲労累積損傷計算とは

🧑‍🎓
「疲労累積損傷」って何ですか? 材料が疲れて壊れるまでの「疲労度」みたいなものですか?
🎓
ざっくり言うと、その通りだね。例えば、飛行機の翼は離着陸のたびに荷重が変動するよね。大きな荷重が1回かかるだけでは壊れないけど、小さな荷重を何百万回も繰り返すと、目に見えない傷が蓄積して突然壊れる。その「傷の蓄積度合い」を数値で表したのが累積損傷Dだ。このシミュレーターでは、上の「材料」や「荷重スペクトル」を選ぶと、リアルタイムでDの値が計算されて、メーターで危険度がわかるよ。
🧑‍🎓
え、そうなんですか! で、メーターのD=1で破損って書いてありますけど、実際に1になったら絶対壊れるんですか?
🎓
実はそこがポイントで、絶対じゃないんだ。Miner則という単純な理論上はD=1で破損だけど、実際の実験では0.3から3までバラつく。だから実務では「安全率」が超重要で、シミュレーターの右側に出てくる「安全率 = 1/D」を見るんだ。例えばDが0.5なら安全率は2.0。試しに「応力振幅S₁」のスライダーを大きくしてみて、Dが0.8を超えると安全率が1.25を切って、メーターが黄色や赤に変わるのがわかるよ。
🧑‍🎓
なるほど!「平均応力」と「Goodman補正」ってのもありますけど、これは何のためにあるんですか?
🎓
いいところに気が付いたね。例えば、回転するシャフトは引っ張り方向に常に力がかかっている(平均応力が正)。この状態で振動が加わると、平均応力がゼロの時より疲労寿命が短くなるんだ。Goodman補正はその影響を考慮して「等価な応力振幅」を計算するための式だ。ツールで「平均応力σₘ」を0から100MPaくらいに増やしてみて。同じ応力振幅でも、Dの値がどんどん大きくなって、残存寿命が短くなるのが確認できるはずだよ。

物理モデルと主要な数式

材料の疲労寿命(破断までの繰返し数N)と応力振幅S_aの関係を表すS-N曲線です。両対数グラフでほぼ直線になることが多く、ここではBasquinの式を用いています。

$$N = \left(\frac{\sigma_f'}{S_a}\right)^{1/b}$$

ここで、$N$: 破断までの繰返し数、$S_a$: 応力振幅 [MPa]、$\sigma_f'$: 疲労強度係数 [MPa](材料定数)、$b$: 疲労強度指数(材料定数、通常は負の値)。例えば、鋼は$\sigma_f'=1000, b=-0.085$という値がツールに設定されています。

異なる応力レベルでの疲労損傷を足し合わせるためのMinerの線形累積損傷則です。また、平均応力の影響を考慮する修正Goodman則も示します。

$$D = \sum_i \frac{n_i}{N_i}, \quad S_{a,eq}= \dfrac{S_a}{1 - \sigma_m / S_u}$$

ここで、$D$: 累積損傷度($D \geq 1$で破損)、$n_i$: 応力レベル$i$での実際の繰返し数、$N_i$: 応力レベル$i$での破断繰返し数(S-N曲線から算出)。$S_{a,eq}$: 平均応力$\sigma_m$を考慮した等価応力振幅、$S_u$: 材料の引張強さ。Goodman補正により、$S_a$を$S_{a,eq}$に置き換えて寿命$N$を計算します。

実世界での応用

自動車・鉄道車両:路面やレールの凹凸による繰返し荷重が車体やサスペンションに加わります。走行距離(繰返し数)と想定される荷重スペクトルから、定期点検や部品交換の目安となる残存寿命を計算するために用いられます。

航空宇宙機体:離着陸時の圧力変化や乱気流による荷重変動が主翼や胴体に加わります。非常に高い信頼性が要求されるため、Miner則による計算結果に大きな安全率を乗じて設計し、実際にはフルスケールの疲労試験で検証します。

発電プラント・回転機械:タービンブレードや発電機のシャフトは、遠心力(平均応力)に加え、振動(応力振幅)が常時作用しています。Goodman補正を適用した疲労評価は、これらの機械の保全計画において不可欠です。

橋梁・建築構造物:風や交通量の変動による荷重が長期にわたって作用します。特に溶接部などのき裂進展起点となりやすい箇所に対して、保守管理の優先度を判断する材料として累積損傷度が参照されます。

よくある誤解と注意点

このツールを使い始める際、特に初心者の方が陥りがちな落とし穴がいくつかあります。まず第一に、「計算結果が絶対的な寿命を保証する」という誤解です。先輩エンジニアからよく言われるのは、「シミュレーションはあくまで『机上の計算』。現実はもっと厳しい」ということ。例えば、ツールで鋼材を選び、安全率が1.5と出たからといって、そのまま設計してはいけません。実際の製品には、計算に入れていない腐食製造上の欠陥(溶接ビードの微細なき裂など)、想定外の過大荷重が必ず存在します。実務では、計算で得られた安全率にさらに「経験係数」を乗じるのが普通です。

第二に、荷重スペクトルの入力の質です。ツールではあらかじめ用意されたスペクトルを選べますが、実際の設計では、実測データや規格に基づいて自分で作成する必要があります。ここで、低応力レベルの繰り返しを「無視してしまう」のは大きな間違い。例えば、10MPaの小さな振動を100万回与えると、Dの値は微増に見えても、材料の内部では微小き裂が確実に発生・進展しています。この「下地」ができた状態で大きな荷重が一度かかると、予想より早く破断に至る「下地効果」が起こり得ます。

最後に、材料定数の信頼性。ツールに登録されているσ_f'やbの値は代表値です。実際の材料では、熱処理のバラつきやロット差によって疲労強度は変動します。例えば、同じ「S45C」という鋼材でも、焼き入れ・焼き戻しの条件次第で疲労寿命が数倍変わることは珍しくありません。信頼性の高い設計のためには、自社で調達する材料の疲労試験データを取得し、ツールのパラメータをカスタマイズすることが理想です。

関連する工学分野

この疲労累積損傷の計算は、単体で完結するものではなく、CAEの広大な世界への「入り口」です。まず直接的に繋がるのが「振動解析」と「耐久性設計」です。例えば、自動車のエンジンマウントの耐久性を評価する場合、まず路面からの振動を多体動力学シミュレーションで計算し、その応力時系列データを「レインフロー計数法」などの手法で荷重スペクトルに変換します。このツールで扱っているのは、まさにその最終ステップの計算部分なのです。

さらに発展させると、「き裂進展解析」の分野と深く関連します。Miner則は「き裂が無い状態からの」累積損傷を扱いますが、既存のき裂(初期欠陥)がある場合の寿命予測には「破壊力学」の考え方、特にパリの法則が使われます。ここでは応力拡大係数ΔKが主役となり、き裂長さが時間とともにどのように進展するかを計算します。疲労設計の実務では、Miner則による初期寿命とき裂進展による残存寿命を組み合わせて、総合的な安全性を評価します。

また、複合材料(CFRP)を扱う場合、その異方性(方向によって強度が異なる性質)から「多軸応力疲労」の考え方が必要になります。ツールでは一軸(引張圧縮)の応力を仮定していますが、現実の部品は複雑な応力状態にあります。これを評価するには、等価応力(フォン・ミーゼス応力など)に変換したり、より高度な多軸疲労基準(臨界面アプローチなど)を用いたりする必要があり、それが次の学習ステップとなります。

発展的な学習のために

もしこのツールの計算に興味を持ち、もっと深く知りたいと思ったら、次のステップを踏むことをお勧めします。まずは数学的な背景を固めましょう。S-N曲線のBasquinの式 $$N = \left(\frac{\sigma_f'}{S_a}\right)^{1/b}$$ は、両対数グラフ上で直線になることを前提としています。この「両対数プロットで線形」という関係を理解するために、対数計算($$ \log N = \frac{1}{b} \log \sigma_f' - \frac{1}{b} \log S_a $$)に慣れることが第一歩です。

次に、「確率的」な視点を取り入れてください。疲労現象は本質的にバラつきが大きいため、信頼性工学の考え方が不可欠です。例えば、S-Nデータを確率分布(ワイブル分布など)で表現し、「99%の生存確率を持つS-N曲線」を導出する方法を学びましょう。これが「P-S-N曲線」と呼ばれるもので、航空宇宙分野など超高信頼性が要求される設計では標準的に用いられます。

実践的な次のトピックとしては、「レインフロー計数法」を強く推奨します。これは、実際の変動応力の時系列データ(例えば、テストコースを走行する車両の車軸に計測されたひずみデータ)を、このツールが入力として想定する「応力振幅とその繰り返し数」のペア(荷重スペクトル)に変換するアルゴリズムです。この手法を学べば、シミュレーションと実測の橋渡しができるようになり、設計の精度が格段に向上します。まずは、単純なサイン波が重なった複雑な波形から、どのようにサイクルを抽出するのか、その基本原理を理解することから始めてみてください。