Paris則 — CAE用語解説

カテゴリ: 用語集 | 2026-01-15
CAE visualization for paris law - technical simulation diagram

Paris則

🧑‍🎓

先生、Paris則の式 da/dN = C(ΔK)^m ってよく教科書に出てくるんですけど、この式は何を意味してるんですか？


理論と物理

Paris則の基本概念

🧑‍🎓

Paris則って、教科書には「亀裂進展速度を表す経験則」と書いてありますが、具体的に何を「経験」して得られた式なんですか？

🎓

良い質問だ。1960年代にP.C. Parisらが、アルミニウム合金や鋼材の疲労試験データを大量にプロットした結果、亀裂進展速度

$$ da/dN $$
と応力拡大係数範囲
$$ \Delta K $$
の対数グラフ上で、中間領域が直線になることを発見したんだ。つまり、
$$ \log(da/dN) $$
と
$$ \log(\Delta K) $$
の関係が線形になる、という「経験」に基づいている。これがParis則の基本形
$$ da/dN = C (\Delta K)^m $$
につながる。

🧑‍🎓

応力拡大係数範囲

$$ \Delta K $$
はどうやって求めるんですか？ 最大応力と最小応力だけで決まるんですか？

🎓

基本的には、最大応力と最小応力から求める応力範囲

$$ \Delta \sigma = \sigma_{max} - \sigma_{min} $$
と、現在の亀裂長さ
$$ a $$
、そして幾何形状を表す補正係数
$$ Y $$
を使って計算する。式は
$$ \Delta K = Y \Delta \sigma \sqrt{\pi a} $$
だ。この
$$ Y $$
は、無限板中心き裂なら1.0だが、実際の構造物では有限幅や自由端の影響で1.0より大きくなる。例えば、幅50mmの板の端に10mmのき裂がある場合、
$$ Y $$
は約1.12になる。

🧑‍🎓

材料定数の

$$ C $$
と
$$ m $$
は、例えば普通の構造用鋼（SM490）だとどのくらいの値なんですか？

🎓

ASTM E647などの標準試験法で測定される。大まかな目安として、空気中の鉄鋼材料では、指数

$$ m $$
は2から4の範囲にあることが多い。例えば、あるSM490鋼のデータでは、
$$ C = 6.9 \times 10^{-12} $$
(単位は m/cycle, MPa√m)、
$$ m = 3.0 $$
という報告がある。ただし、これはあくまで一例で、実際の設計では対象材料の試験データか、信頼できるデータベース（例えば、日本材料学会の疲労データシート）から選定する必要がある。

数値解法と実装

CAEソフトウェアでの扱い方

🧑‍🎓

CAEソフトで疲労き裂解析をするとき、Paris則は具体的にどのモジュールや機能の中で使われるんですか？

🎓

主に「き裂進展解析」や「疲労き裂解析」専用のモジュールで使われる。例えば、Abaqusでは「Abaqus/Standard」の「Contour Integral」計算機能で

$$ \Delta K $$
を求め、「Direct Cyclic」分析ステップなどと組み合わせて、Paris則に基づくき裂進展シミュレーションを実行できる。Ansysでは「Ansys Mechanical」の「Fracture」ツールボックス内の「Fatigue Crack Growth」機能が該当する。

🧑‍🎓

解析では、き裂が進むたびにメッシュを更新する必要がありますか？ それとも自動でやってくれるんですか？

🎓

手法による。主流は「メッシュリメッシュ法」と「拡張有限要素法（XFEM）」だ。メッシュリメッシュ法は、き裂先端が進むと、その周辺のメッシュを削除して新しいき裂形状に合わせて再メッシュを作成する。多くのソフトで（例えばFRANC3Dなど）自動化されているが、計算コストがかかる。一方、AnsysやAbaqusに実装されているXFEMは、き裂が要素を横切るように進むため、メッシュを更新せずに済む。ただし、進展方向の予測アルゴリズム（最大主応力方向など）の設定が重要になる。

🧑‍🎓

Paris則の式

$$ da/dN = C (\Delta K)^m $$
は、ソルバー内部でどう積分されて、き裂長さの増分が計算されるんですか？

🎓

通常、増分計算にオイラー法のような単純な数値積分が使われる。例えば、ある荷重サイクルブロック

$$ \Delta N $$
でのき裂進展量
$$ \Delta a $$
は、現在のき裂長さ
$$ a_i $$
から求まる
$$ \Delta K_i $$
を使って、
$$ \Delta a = (da/dN)_i \times \Delta N = C (\Delta K_i)^m \times \Delta N $$
と計算する。そして、き裂長さを
$$ a_{i+1} = a_i + \Delta a $$
で更新する。この
$$ \Delta N $$
の取り方は解析の精度と効率に直結するので、ソフトではユーザーが「進展増分」を設定できるようになっている。

実践ガイド

解析ワークフローと注意点

🧑‍🎓

実際に航空機部品のき裂進展寿命をParis則で評価したい場合、最初に何から始めればいいですか？

🎓

まずは「初期き裂サイズの定義」だ。これは規格で決まっている。例えば、航空宇宙分野では損傷許容設計（Damage Tolerance）が求められ、初期欠陥サイズは製造・検査能力に基づいて決める。ミリ波探傷検査（NDI）が可能なら、初期き裂長さを1.27mm（0.05インチ）と仮定することがFAAのアドバイザリー・サーキュラーで示されていることもある。この初期値を間違えると、寿命予測が現実とかけ離れたものになる。

🧑‍🎓

荷重データはどう準備するんですか？ 実際の飛行データはランダム振動みたいな不規則な波形ですが。

🎓

その通り。実稼働荷重履歴（Service Loading History）は不規則だ。そこで「雨流計数法」などの手法で、応力範囲

$$ \Delta \sigma $$
と平均応力の組として統計化する。その後、平均応力の影響を考慮するために、
$$ \Delta K $$
の代わりに有効応力拡大係数範囲を使うか、Walkerの式のような平均応力補正をParis則に組み込む。ソフトウェア（例えばnCode DesignLife）では、この雨流計数と疲労き裂進展計算を連携させるワークフローが組まれている。

🧑‍🎓

解析結果として「あと何サイクルで破壊する」という寿命が出ますが、その結果をどう検証すれば信頼できますか？

🎓

絶対的な検証は実物試験しかない。ただし、解析の妥当性を確認するチェックリストはある。1) 使用した材料定数

$$ C, m $$
が、解析対象の環境（湿気、温度）と合っているか。2) き裂進展中に、
$$ \Delta K $$
が材料のしきい値
$$ \Delta K_{th} $$
を下回ったり、破壊靭性
$$ K_{IC} $$
に近づいたりしていないか（Paris則が適用できない領域）。3) き裂進展経路が物理的に妥当か（例えば、モードI（開口型）が主か）。これらのチェックをせずに寿命の数値だけを信じてはいけない。

ソフトウェア比較

主要ソフトの機能と特徴

🧑‍🎓

Ansys、Abaqus、NASGRO（旧NASAのコード）で、Paris則を使ったき裂進展解析のアプローチはどう違うんですか？

🎓

核となるParis則の式は同じだが、統合度とデータベースが大きく異なる。Ansys Mechanicalは汎用構造解析に「Fatigue Crack Growth」ツールが統合されており、XFEMによるメッシュ非依存の解析が売りだ。Abaqusも同様で、Pythonスクリプトによるカスタマイズ性が高い。一方、NASGRO（現在は商用ソフト）は、航空宇宙分野に特化して開発された経緯があり、その最大の強みは膨大な材料のき裂進展データベースと、Paris則を含むより高度な進展則（NASGRO方程式）を内蔵している点だ。規格（ESDU, MMPDS）に準拠したデータが使える。

🧑‍🎓

COMSOL Multiphysicsでも疲労き裂解析はできますか？ その場合、Paris則の実装はどうなっていますか？

🎓

COMSOLでは「Structural Mechanics Module」の「Fracture Mechanics」機能でJ積分や応力拡大係数を計算できる。しかし、き裂が自動進展する組み込みの「疲労き裂進展ソルバー」は、私の知る限りAnsysやAbaqusほど直接的ではない。代わりに、「変数

$$ da/dN $$
をParis則で定義し、時間（またはサイクル）に対してき裂長さパラメータを更新する」というプロセスを、COMSOLの「Events」インターフェースや「ODE（常微分方程式）インターフェース」を利用してユーザー自身が構築するアプローチが取られる。柔軟性は高いが、専門知識が必要だ。

🧑‍🎓

無料やオープンソースのソフト（CalculiX、Code_Asterなど）ではParis則による解析は可能ですか？

🎓

可能だが、作業の多くは手動またはスクリプトに依存する。例えば、Code_Aster（Salome-MECAプラットフォーム）は非常に強力なき裂力学機能を持ち、

$$ K $$
の計算はできる。しかし、Paris則に基づいてき裂を自動進展させ、メッシュを更新する一連のプロセスを実行するには、専用のコマンド（DEFI_MATERIAUで進展則を定義し、PROPAGATION_MAILLAGEなどを使用）を正しく組み合わせる必要がある。商用ソフトのようなGUIベースのワンクリック操作とは程遠く、理論とソフトの両方に対する深い理解が必須だ。

トラブルシューティング

よくあるエラーと対策

🧑‍🎓

解析を実行したら、き裂が予想外の方向に急に曲がって進んでしまいました。考えられる原因は何ですか？

🎓

主に二つの原因が考えられる。第一に、き裂先端の応力場の計算精度が不十分。特に、き裂先端のメッシュが粗すぎると、応力拡大係数

$$ K_{II} $$
（面内せん断モード）の計算値が不安定になり、進展方向予測が狂う。第二に、進展方向の判定基準の問題。多くのソフトは最大主応力方向（
$$ \theta = 2 \arctan\left(\frac{1}{4}\frac{K_I}{K_{II}} \pm \sqrt{(\frac{K_I}{K_{II}})^2 + 8}\right) $$
のような式）を使うが、複雑な多軸応力状態下ではこの基準自体が不適切な場合がある。まずはき裂先端のメッシュを細かくし、特異要素（シンギュラリティ要素）を使っているか確認せよ。

🧑‍🎓

寿命予測の結果が、過去の実験データと比べて極端に短く（または長く）出ます。材料定数は実験データから取ったのに、なぜこんなにずれるんですか？

🎓

Paris則の式

$$ da/dN = C (\Delta K)^m $$
は対数グラフ上の直線部分だけを表す。ずれの原因として最も多いのは、解析で計算される
$$ \Delta K $$
が、実験で使われた標準試験片（CT試験片など）の
$$ \Delta K $$
と「同じ計算方法で求められているか」の見落としだ。特に、幾何形状補正係数
$$ Y $$
が正しいか。また、解析対象の構造物では、き裂が進むに連れて
$$ Y $$
が大きく変化する。この
$$ Y $$
の変化を正しく反映した「き裂進展曲線」と、定数
$$ Y $$
を使った単純計算とでは、寿命が数倍違ってくる。

🧑‍🎓

XFEMを使って解析しているのですが、き裂進展の計算が特定のサイクル数で突然停止して「収束しません」というエラーが出ます。

🎓

XFEMでの収束失敗はよくある。直接的な原因は、き裂先端が要素内の特定の位置（例えばガウス積分点の近く）に来た時に、レベルセット関数の更新で数値的不安定性が生じることだ。対策はいくつかある。1) き裂進展増分

$$ \Delta a $$
を小さくする。デフォルトが0.001mなら、0.0005mにしてみる。2) ソルバーの設定を変える。Abaqus/Standardなら、「ステップ」設定で「より大きな増分サイズを許可」をオンにしたり、イテレーション行列の更新頻度を上げたりする。3) 根本的には、き裂進展経路付近のメッシュを均一で細かいものに変更する。XFEMはメッシュ非依存と言われるが、荒いメッシュでは進展挙動が不安定になりやすい。

関連トピック

用語集Coffin-Manson則 — CAE用語解説用語集Miner則 — CAE用語解説用語集SIF — CAE用語解説用語集損傷許容設計 — CAE用語解説用語集応力拡大係数 — CAE用語解説用語集CTOD — CAE用語解説
この記事の評価
ご回答ありがとうございます！
参考に
なった
もっと
詳しく
誤りを
報告
参考になった
0
もっと詳しく
0
誤りを報告
0
Written by NovaSolver Contributors
Anonymous Engineers & AI — サイトマップ
プロフィールを見る