$$\frac{da}{dN}= C(\Delta K)^m$$
$\Delta K = \Delta\sigma \cdot Y\sqrt{\pi a}$
$Y = \sec\!\left(\frac{\pi a}{2W}\right)^{1/2}$(有限幅補正)
臨界き裂:$a_c = \frac{1}{\pi}\!\left(\frac{K_{Ic}}{\sigma_{max} Y}\right)^2$
初期き裂寸法・応力振幅・材料定数を操作して、き裂が臨界寸法に達するまでの伝播挙動と残余寿命をリアルタイムで計算・可視化しよう。
$$\frac{da}{dN}= C(\Delta K)^m$$
$\Delta K = \Delta\sigma \cdot Y\sqrt{\pi a}$
$Y = \sec\!\left(\frac{\pi a}{2W}\right)^{1/2}$(有限幅補正)
臨界き裂:$a_c = \frac{1}{\pi}\!\left(\frac{K_{Ic}}{\sigma_{max} Y}\right)^2$
航空機の定期点検間隔の設定:機体の主要構造部材(ウィングスキンなど)に想定される最大の初期欠陥サイズから出発し、実際の飛行荷重履歴を考慮してパリス則でき裂進展を予測します。き裂が臨界寸法に達する前に点検で発見できるよう、その「残余寿命」よりも短い間隔で検査が設定されます。
発電プラント・橋梁などの長寿命構造物の健全性評価:数十〜百年単位で使用されるインフラ構造物では、経年劣化としての疲労き裂進展が問題となります。保守計画策定のため、現在検知されているき裂サイズから将来の進展を予測し、補修や取替えの時期を判断する材料データとして活用されます。
自動車のシャシー・エンジン部品の耐久設計:設計段階で、想定される荷重条件下での部品寿命をパリス則を用いて予測します。特に、軽量化が求められる部品では、き裂進展寿命を正確に見積もることで、必要以上に肉厚にすることなく、安全で信頼性の高い設計が可能になります。
金属アディティブマニュファクチャリング(3Dプリント)製品の信頼性評価:積層造形により内部に生じうる微小な気孔や不完全融合部を「初期き裂」とみなして疲労寿命を評価します。造形パラメータや熱処理が材料定数C, mに与える影響を調査し、製品の信頼性保証に役立てられます。
このシミュレーターを使い始める際、特にCAE初心者が陥りがちな落とし穴がいくつかあります。まず「パリス則は万能ではない」という点。この法則は、き裂進展の中間領域(いわゆる「領域II」)でよく成り立ちますが、進展が非常に遅い領域(下限界)や、破壊直前の高速領域では別のモデルが必要です。シミュレーターの結果をそのまま絶対的な寿命として信じるのは危険で、あくまで目安と考えてください。
次に、初期き裂サイズa₀の設定の重要性。例えば、a₀を0.1mmと1mmで計算すると、残余寿命が数倍も変わることがあります。実務では、非破壊検査で検出可能な最小サイズをa₀に設定することが多いです。「傷はないはず」とa₀を極端に小さく設定すると、現実とはかけ離れた楽観的な寿命が出てしまうので注意。
最後に、「応力範囲Δσは一定ではない」という現実。シミュレーターでは固定値ですが、実際の機械や構造物は様々な大きさの荷重がランダムに加わります。このような変動振幅荷重に対しては、ミナー則などと組み合わせた累積損傷計算が必要になります。ツールで基本を学んだら、この「荷重履歴の扱い」が次の実践的な課題です。
鋼部材の疲労解析:初期き裂a₀=2.0mm、応力範囲Δσ=100MPa、鋼プリセット(C=6.9×10⁻¹²、m=3.0)、板幅W=100mm、応力比R=0.1、破壊靭性K_Ic=50MPa√mで計算すると、臨界き裂a_c≈47.4mm、残余寿命N≈0.90×10⁶サイクル、ΔK_max≈7.9MPa√m、安全率a_c/a₀≈23.7が得られます。アルミ合金プリセット(C=5.6×10⁻¹¹、m=2.9)に切り替えると伝播が速まり、N≈0.14×10⁶サイクルまで短くなります。