延性破壊 — CAE用語解説

大きな違いは、破壊前に材料が「塑性変形」するかどうかです。延性破壊では、例えば軟鋼（S45C）を引っ張ると、局所的にくびれ（ネッキング）が生じ、断面積が大幅に減少します。破断面は繊維状で、肉眼でも変形がわかります。一方、脆性破壊（鋳鉄など）はほとんど変形なしに突然割れ、破断面は粒状で平坦です。

良い着眼点だ。一軸引張りでは、公称応力は下がっても、真応力は上がり続けます。しかし、材料の「加工硬化」能力が、局所的な断面積減少による応力集中に追いつかなくなった点でネッキングが開始します。Considèreの条件として知られています。真応力-真ひずみ曲線で、

多くの金属では、「ボイドの発生・成長・合体」という微視的プロセスです。材料中の介在物（MnSなど）や第二相粒子の界面で微細な空洞（ボイド）が発生し、大きな塑性ひずみによってこれらが成長、隣接するボイド同士がつながって（合体）、最終的な破断に至ります。破断面を電子顕微鏡で見ると、ディンプル（小さなくぼみ）と呼ばれる特徴的な模様が観察されます。

単純な塑性ひずみの閾値だけでは不十分です。代表的なアプローチは「損傷力学」に基づくモデルです。例えば、LemaitreやGursonのモデルがあります。Gursonモデルは、ボイド体積分率

その通りです。これが「応力三軸度」の影響です。応力三軸度

「要素消去法」が古典的ですが、応力が不連続にゼロになるため動解析で問題が出ます。現在は「粘性正則化」や「弾性剛性低減法」が一般的です。例えば、Abaqusの延性損傷モデルでは、破断ひずみに達した後も、指定したひずみ（例：0.05）の範囲で徐々に応力を減衰させ、要素の剛性を滑らかにゼロに近づけます。これにより、過度な収束問題や不安定なエネルギー解放を防ぎます。

その通りで、これが最大のハードルです。最低限、異なる応力三軸度下での実験データが必要です。実務では、平滑材の引張試験（η≈0.33）、ノッチ付き丸棒試験（η≈0.6-1.0）、純せん断試験（η≈0）などを組み合わせ、破断時のひずみデータを取得します。例えば自動車鋼板の成形限界図（FLD）作成では、JIS Z 2254に準拠した Nakazima試験が行われ、各ひずみ経路での破断データが収集されます。

局所的なひずみ集中を追うため、メッシュ依存性は本質的な問題です。対策は主に二つ。一つは「メッシュサイズに依存したパラメータ調整」です。例えば、破断エネルギー密度が一定となるように、メッシュの特性長さ

破壊位置の一致は必要条件ですが、十分条件ではありません。荷重-変位曲線のピーク後の減衰挙動、破壊に至るエネルギー吸収量、そして複数の異なる試験形状（例えば、平面ひずみ状態と平面応力状態）での破壊挙動を再現できて初めて、モデルの信頼性が高まります。自動車のクラッシュ解析なら、コンポーネント単体の三点曲げ試験と、実車衝突試験の両方で整合性を取る必要があります。

大きな違いがあります。Abaqus/Explicitは「延性損傷」モデルが標準で組み込まれており、応力三軸度とひずみ速度を考慮した破壊ひずみを定義できます。また、Gurson-Tvergaard-Needleman（GTN）モデルも利用可能で、ボイド成長を直接モデル化できます。一方、Ansys Mechanical（APDL/LS-DYNA以外）の標準的な塑性モデルには破壊判定がなく、ユーザーサブルーチン（UPF）で独自モデルを実装するか、LS-DYNA連携が一般的です。

LS-DYNAの強みは、破壊モデルの豊富さと実績です。特に「GISSMO（Generalized Incremental Stress State dependent damage Model）」は、応力三軸度とロード角（Lode角）の両方を考慮した高度なモデルで、鋼板の複雑な破壊モード（せん断破壊など）を精度良く再現できます。また、MAT_ADD_EROSION（要素消去条件）で、塑性ひずみ、応力、時間など、様々な単純な閾値を設定する簡易的な方法も用意されており、ユーザーの目的とデータ量に応じて選択できます。

COMSOL Multiphysicsでは、「構造力学モジュール」の「損傷」インターフェースで延性損傷を扱えます。デフォルトでLemaitreの損傷モデルが実装されており、損傷変数Dの進化則を定義できます。また、ユーザーが偏微分方程式（PDE）インターフェースを使って独自の損傷進化則を組み込む柔軟性が高いのが特徴です。ただし、AbaqusやLS-DYNAのような自動車衝突解析向けに最適化された多数の実験フィッティング済みモデルは少なく、ユーザー自身が理論と実験データに基づいてモデルを構築する必要があります。

これには主に二つの原因があります。第一に、要素消去や剛性低下により、モデルの質量が急激に失われたり、剛性マトリックスが悪条件になったりするためです。第二に、破壊が進展する先端で非常に大きなひずみ速度が発生し、陽解法の安定条件であるクーラン条件を厳しくするためです。対策としては、先ほど話した「粘性正則化」パラメータを調整して破壊をなだらかにすること、および「質量スケーリング」を適用して、破壊領域以外の要素のタイムステップを維持することが有効です。

これは「損傷の局所化帯」が一つの要素幅に集中し、隣接要素に伝播しない典型的な問題です。非局所化モデルが解決策です。もし使っているソルバーにその機能がない場合は、応力三軸度などの破壊パラメータを計算する際に、要素中心の値ではなく、周囲数要素の平均値を使うようなフィルタリングをユーザーサブルーチンで実装する必要があります。あるいは、メッシュをある程度粗く保ち、破壊エネルギーを正しく散逸させる「メッシュサイズ調整」手法に徹する選択肢もあります。

まず、入力パラメータが「本当に同じ」か確認してください。例えば、破断ひずみが真ひずみか工学ひずみか、応力三軸度の定義（η = p/q など、分母分子が逆の定義もある）、デフォルトの要素積分法（完全積分/低減積分）による応力・ひずみの計算精度の違いが原因です。信頼するのは、単純な検証モデル（例えば、ノッチ付き引張試験）で実験結果と最も整合性の高いソルバーです。ソルバーのマニュアルで、使用しているモデルの理論的詳細を必ず確認し、比較する際は「ブラックボックス」として扱わないことが重要です。