パリス則 き裂伝播 戻る
破壊力学・疲労

パリス則 疲労き裂伝播シミュレーター

初期き裂寸法・応力振幅・材料定数を操作して、き裂が臨界寸法に達するまでの伝播挙動と残余寿命をリアルタイムで計算・可視化しよう。

材料・荷重設定
材料
Paris定数 C (m/cycle)
Paris指数 m
初期き裂 a₀ (mm)
mm
板幅 W (mm)
mm
応力範囲 Δσ (MPa)
MPa
応力比 R
破壊靭性 K_Ic (MPa√m)
MPa√m
計算結果
臨界き裂 a_c (mm)
残余寿命 N (×10⁶)
ΔK_max (MPa√m)
安全率 a_c/a₀
き裂寸法 a-N 曲線
Paris則 log(da/dN)-log(ΔK)
き裂進展アニメーション
理論・主要公式

$$\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$

パリス則による疲労き裂伝播とは

🙋
「疲労き裂伝播」って何ですか? き裂が進むスピードがわかるということですか?
🎓
その通り!大まかに言うと、金属に小さな傷(き裂)があって、繰り返し力がかかると、その傷が少しずつ大きくなっていく現象だ。その「1回の荷重でどれだけ傷が進むか」を予測するのがパリス則だよ。このシミュレーターでは、上のスライダーで「応力範囲Δσ」や「初期き裂a₀」を変えると、き裂が進んでいく様子がグラフでリアルタイムに見えるんだ。
🙋
え、そうなんですか!材料定数の「C」と「m」って何を変えてるんですか?数字が変わるとグラフの形が大きく変わりますね。
🎓
いいところに気づいたね。「C」と「m」は材料ごとに決まる定数で、き裂の進みやすさを決める。例えば、鋼はmが3前後、アルミは2.9前後だ。mの値を大きくすると、き裂が進むスピードが応力に非常に敏感になる。実務では、材料のデータシートからこの値を取ってくることが多いよ。右の「パリス則プロット」を見ると、この関係が直線で表されるのがわかる。
🙋
「臨界き裂寸法」って、これ以上大きくなると一気に壊れてしまうということですか?「破壊靭性K_Ic」を変えるとどうなりますか?
🎓
その解釈でバッチリだ!「破壊靭性K_Ic」は材料がき裂を受け止める強さだ。このシミュレーターでK_Icの値を大きくしてみてごらん。臨界き裂寸法a_cが大きくなるだろう?例えば、非常に靭性の高い鋼はK_Icが大きいので、大きな傷があっても突然破壊しにくい。逆に、もろい材料は小さな傷で壊れる。残余寿命の計算は、き裂がこのa_cに達するまでのサイクル数を積分で求めてるんだ。

よくある質問

Cとmは材料固有の値で、実験データから決定されます。一般的な金属材料の参考値は文献や材料データベースで確認できます。本シミュレーターでは初期値としてアルミニウム合金の代表値(C=1e-11, m=3)を設定していますが、解析対象に合わせて変更してください。
いいえ、本シミュレーターはパリス則に基づき、き裂が臨界寸法(入力された限界き裂長さ)に達するまでの伝播挙動と残余寿命のみを計算します。臨界以降の急速破壊や破面解析は対象外です。あくまで疲労き裂伝播寿命の予測ツールとしてご利用ください。
パリス則は線形弾性破壊力学に基づくため、応力振幅が大きすぎるとき裂先端の塑性域が無視できなくなり、適用範囲を超えます。目安として、最大応力が材料の降伏応力の30%を超える場合は結果の精度が低下します。その場合は弾塑性破壊力学の適用を検討してください。
有限幅の板における中央き裂の厳密解として、Feddersenの補正式(sec関数)を用いています。この式は板幅に対するき裂長さの比が大きくなるほど、き裂進展が加速される効果(有限幅補正)を正確に表現します。他の形状を解析したい場合は、対応するYの式に置き換えてご利用ください。

実世界での応用

航空機の定期点検間隔の設定:機体の主要構造部材(ウィングスキンなど)に想定される最大の初期欠陥サイズから出発し、実際の飛行荷重履歴を考慮してパリス則でき裂進展を予測します。き裂が臨界寸法に達する前に点検で発見できるよう、その「残余寿命」よりも短い間隔で検査が設定されます。

発電プラント・橋梁などの長寿命構造物の健全性評価:数十〜百年単位で使用されるインフラ構造物では、経年劣化としての疲労き裂進展が問題となります。保守計画策定のため、現在検知されているき裂サイズから将来の進展を予測し、補修や取替えの時期を判断する材料データとして活用されます。

自動車のシャシー・エンジン部品の耐久設計:設計段階で、想定される荷重条件下での部品寿命をパリス則を用いて予測します。特に、軽量化が求められる部品では、き裂進展寿命を正確に見積もることで、必要以上に肉厚にすることなく、安全で信頼性の高い設計が可能になります。

金属アディティブマニュファクチャリング(3Dプリント)製品の信頼性評価:積層造形により内部に生じうる微小な気孔や不完全融合部を「初期き裂」とみなして疲労寿命を評価します。造形パラメータや熱処理が材料定数C, mに与える影響を調査し、製品の信頼性保証に役立てられます。

よくある誤解と注意点

このシミュレーターを使い始める際、特にCAE初心者が陥りがちな落とし穴がいくつかあります。まず「パリス則は万能ではない」という点。この法則は、き裂進展の中間領域(いわゆる「領域II」)でよく成り立ちますが、進展が非常に遅い領域(下限界)や、破壊直前の高速領域では別のモデルが必要です。シミュレーターの結果をそのまま絶対的な寿命として信じるのは危険で、あくまで目安と考えてください。

次に、初期き裂サイズa₀の設定の重要性。例えば、a₀を0.1mmと1mmで計算すると、残余寿命が数倍も変わることがあります。実務では、非破壊検査で検出可能な最小サイズをa₀に設定することが多いです。「傷はないはず」とa₀を極端に小さく設定すると、現実とはかけ離れた楽観的な寿命が出てしまうので注意。

最後に、「応力範囲Δσは一定ではない」という現実。シミュレーターでは固定値ですが、実際の機械や構造物は様々な大きさの荷重がランダムに加わります。このような変動振幅荷重に対しては、ミナー則などと組み合わせた累積損傷計算が必要になります。ツールで基本を学んだら、この「荷重履歴の扱い」が次の実践的な課題です。

使い方ガイド

  1. 初期き裂寸法a₀(mm)を設定します。鋳鉄部品の場合0.5~2.0mm、溶接部は2.0~5.0mmが一般的です
  2. 応力振幅ΔS(MPa)を入力します。回転曲げ試験機での標準値は50~200MPa範囲です
  3. 材料定数C値とm値を設定します。アルミ合金2024-T3ではC=1.62×10⁻¹¹、m=3.16が標準値です
  4. 板厚w(mm)を指定後、シミュレーション実行ボタンで伝播挙動を計算開始します
  5. 臨界き裂寸法a_c、残余寿命N、最大応力拡大係数ΔK_maxを確認します

具体的な計算例

航空機翼パネル(Al合金)の疲労解析:初期き裂a₀=1.5mm、応力振幅ΔS=120MPa、C=1.62×10⁻¹¹、m=3.16、板厚w=3.0mmで計算すると、臨界き裂a_c≈15.2mm、残余寿命N≈8.7×10⁶サイクル、ΔK_max≈18.5MPa√mが得られます。溶接鋼(E=200GPa)の場合はC=3.0×10⁻¹²、m=2.8でより低速伝播を示します。

実務での注意点