損傷変数Dってどうやって進展させるんですか？

損傷則（Damage Evolution Law）で記述する。代表的なのがLemaitre（1985年）の延性損傷モデルで、dD/dt = (Y/S)^s * epsilon_dot_p という形をとる——Yは損傷エネルギー解放率、Sは材料定数、epsilon_dot_pは塑性ひずみ速度だ。要するに「塑性変形が起きると損傷が蓄積して剛性が下がる」という連成だ。他にMazarsモデル（コンクリートの引張損傷）、Gursonモデル（ボイド成長に基づく延性破断）など材料・破壊モードに合わせた多様なモデルがある。

FEMで損傷力学を使うときの難しさは何ですか？

メッシュ依存性が最大の難関だ。損傷が進展すると局所的にひずみが集中して、メッシュを細かくするほど損傷域が狭くなって破断変位が変わってしまう——これを「病的なメッシュ依存性」と呼ぶ。対策として非局所損傷モデル（Nonlocal Damage）やひずみ勾配正則化（Gradient-Enhanced Damage）が使われる。これらは損傷の広がり方に「材料の特性長（characteristic length）」を導入して、メッシュを細かくしても収束する解を保証する。実装は複雑だが、き裂進展を物理的に正確に予測するには必須の処理だ。

Gursonモデルって延性破断の解析で有名ですよね？

そう。Arthur Gurson（1977年）が多孔質金属（ボイドを含む材料）の降伏面を理論的に導いた。孔食率f（ボイドの体積分率）を内部変数として、von Misesの古典的な降伏面を拡張した形になっている。Tvergaard-Needlemanが実験フィットのパラメータを加えたGTNモデルが今もよく使われる。アルミや高強度鋼の延性破断——衝突解析での変形・破断シミュレーション——で採用されることが多い。AbaqusではPorous Metal Plasticityとして実装されている。

損傷力学 — CAE用語解説

Q: 先生、損傷力学ってプラスチックモデルとどう違うんですか？ 「損傷」っていう言葉がポイントですか？

そう、キーワードは「材料の劣化」だ。塑性力学はどれだけ変形しても材料の剛性（ヤング率）は変わらない——降伏後も剛性は同じで、ひずみが残るだけだ。損傷力学（Damage Mechanics）は材料内部のマイクロクラック・ボイドの蓄積を「損傷変数D」という内部変数で表現して、Dが増えるにつれて弾性係数が低下する——E_effective = E * (1 - D)——という形で剛性の劣化を記述する。D=0が無損傷、D=1が完全破断（材料としての機能喪失）だ。

カテゴリ: 用語集 | 2026-01-15

CAE visualization for damage mechanics - technical simulation diagram

損傷力学

🧑‍🎓

先生、損傷力学ってプラスチックモデルとどう違うんですか？「損傷」っていう言葉がポイントですか？

理論と物理

損傷力学の基本概念

🧑‍🎓

損傷力学でいう「損傷変数」って、具体的に何を表しているんですか？ひび割れの本数みたいなものですか？

🎓

いいえ、本数ではなく、材料の有効断面積の減少率を表す無次元量です。例えば、初期断面積を

$$ A_0 $$

、微小空洞やき裂で負荷を支えられない面積を

$$ A_D $$

とすると、損傷変数

$$ D $$

は

$$ D = A_D / A_0 $$

と定義されます。D=0は無傷、D=1は完全破断を意味します。実際の金属では、D=0.2〜0.3で巨視的な破断が始まることが多いです。

🧑‍🎓

「有効応力」という言葉も出てきますが、通常の応力とどう違うんですか？

🎓

名目応力

$$ \sigma $$

は初期断面積で割った値ですが、実際に材料が負担しているのは、損傷で失われた部分を除いた有効断面積あたりの応力です。これが有効応力

$$ \tilde{\sigma} $$

で、

$$ \tilde{\sigma} = \sigma / (1 - D) $$

の関係があります。つまり、Dが増大すると、同じ名目応力でも材料が実際に感じる応力は急激に増加し、加速的に破壊が進みます。

🧑‍🎓

等価応力（フォンミーゼス応力）に損傷の効果を組み込むには、どうすればいいですか？

🎓

有効応力の概念を使って置き換えます。損傷を考慮した等価応力

$$ \tilde{\sigma}_{eq} $$

は、有効応力テンソル

$$ \tilde{\boldsymbol{\sigma}} $$

から計算します。具体的には、

$$ \tilde{\sigma}_{eq} = \sqrt{ \frac{3}{2} (\tilde{\boldsymbol{\sigma}}^{dev} : \tilde{\boldsymbol{\sigma}}^{dev}) } $$

となります。ここで

$$ \tilde{\boldsymbol{\sigma}}^{dev} $$

は有効応力の偏差成分です。これにより、損傷が進んだ領域では降伏条件が早く満たされるようになります。

数値解法と実装

損傷進化則とその離散化

🧑‍🎓

損傷変数Dは時間とともにどう増えていくんですか？増え方を決める式は？

🎓

多くのモデルでは、塑性ひずみの蓄積に伴って損傷が進化すると仮定します。Lemaitreのモデルでは、損傷駆動力Yを用いて、

$$ \dot{D} = \left( \frac{Y}{S} \right)^s \dot{p} $$

という進化則が使われます。ここで、

$$ \dot{p} $$

は等価塑性ひずみ速度、Sとsは材料定数です。駆動力Yはひずみエネルギー解放率に関連し、

$$ Y = \frac{\tilde{\sigma}_{eq}^2 R_v}{2E(1-D)^2} $$

で表されます。Eはヤング率、

$$ R_v $$

は三軸応力度の関数です。

🧑‍🎓

その微分方程式をFEMで解く時、どうやって離散化するんですか？陽解法と陰解法で違いは？

🎓

時間積分には通常、陽的オイラー法か陰的後退オイラー法を使います。例えば、時間ステップ

$$ n $$

から

$$ n+1 $$

への更新を陽解法で行うと、

$$ D_{n+1} = D_n + \dot{D}_n \Delta t $$

となります。計算は簡単ですが、条件付き安定です。陰解法では、

$$ D_{n+1} = D_n + \dot{D}_{n+1} \Delta t $$

となり、

$$ \dot{D}_{n+1} $$

が

$$ D_{n+1} $$

に依存するため、ニュートン・ラフソン法などの反復解法が必要ですが、数値的に安定です。Abaqusのダメージモデルの多くは陰的積分を採用しています。

🧑‍🎓

Dが1に近づくと、剛性マトリックスが特異になって計算が発散しませんか？

🎓

その通り、重大な問題です。これを回避するため、「要素消去法」や「剛性減衰法」が使われます。要素消去法では、Dがある閾値（例: 0.99）を超えた要素の剛性を極めて小さな値に設定し、実質的に荷重を支えなくします。LS-DYNAのMAT_ADD_EROSIONなどが該当します。剛性減衰法では、Dの増加に伴って連続的に剛性を減少させ、数値的不安定性を緩和します。いずれにせよ、結果のメッシュ依存性に注意が必要です。

実践ガイド

解析ワークフローと材料定数決定

🧑‍🎓

実際に金属の破壊解析を損傷力学で行う場合、まず何から始めればいいですか？

🎓

まずは材料定数の同定です。最低限、弾性率、降伏応力、硬化則に加え、損傷パラメータ（先ほどのS, sなど）が必要です。これらは、一軸引張試験で破断まで行い、特に破断直前の応力-ひずみ関係の軟化挙動から逆算して決定します。例えば、JIS Z 2241に準拠した試験片を使い、デジタル画像相関法（DIC）で局部のひずみ集中を計測すると、損傷の進展をより正確に捉えられます。

🧑‍🎓

材料定数が決まったら、その次は？いきなり製品モデルで解析？

🎓

それは危険です。必ず「検証と校正」のステップを挟みます。同定した材料定数を使って、試験片レベルの単純なモデル（切欠き付き引張試験など）を再解析し、実験で得られた破断位置、破断荷重、荷重-変位曲線を再現できるか確認します。このステップでパラメータを微調整します。これをスキップすると、製品解析の結果は全く信用できません。

🧑‍🎓

メッシュの大きさによって、破断の予測結果が大きく変わりそうですが、どう対策しますか？

🎓

局所化問題と呼ばれる本質的な課題です。実務では2つのアプローチがあります。1つは「メッシュ依存性低減モデル」の採用。例えば、損傷進化則に材料の特性長さ（例: 0.1mm）を組み込み、ひずみ勾配を考慮する方法です。もう1つは、破壊エネルギーを基準とする方法。要素が消える（または剛性が失われる）際に消散するエネルギーが、材料の破壊エネルギー（鉄鋼で約100 kJ/m²程度）に一致するようにメッシュサイズを調整・設定します。後者はAbaqusのダメージモデルでよく使われます。

ソフトウェア比較

主要ソルバーの損傷モデル実装

🧑‍🎓

Abaqusで使える代表的な損傷・破壊モデルは何ですか？

🎓

主に3つあります。1. **延性損傷モデル**: 塑性ひずみと三軸応力度に基づき損傷を定義。破断ひずみを指定するシンプルなものから、Lemaitreモデルに近い進化則を持つものまで。2. **Johnson-Cook損傷モデル**: 高ひずみ速度・高温環境での破壊予測に強く、衝撃解析で多用。3. **連続体損傷力学（CDM）に基づくモデル**: Abaqusでは「Damage for Traction Separation Laws」と組み合わせ、コヒーシブゾンモデルと併用することが多い。いずれもStandard/Explicitの両方で利用可能です。

🧑‍🎓

Ansys MechanicalとLS-DYNAでは、どういう特徴がありますか？

🎓

**Ansys Mechanical** は、Ansys GRANTAの材料データベースと連携した「Ansys Material Designer」内で損傷パラメータを管理できるのが特徴です。複合材料の層間破壊などに強い。また、「Ansys Explicit Dynamics」では、GISSMO (Generalized Incremental Stress State dependent damage Model) モデルが利用でき、複雑な応力状態下での破壊を精度良く予測できます。
**LS-DYNA** は、破壊解析のバリエーションが非常に豊富。MAT_ADD_EROSIONをはじめ、*MAT_XXX_MODEL_FAILURE_YYYといった多数の破壊基準が用意されており、ユーザーが複数の破壊条件（塑性ひずみ、応力、時間など）をAND/ORで組み合わせて定義できる柔軟性が高いです。自動車衝突解析のデファクトスタンダードです。

🧑‍🎓

COMSOL Multiphysicsのような汎用ソフトでは、損傷モデルをどう扱うんですか？

🎓

COMSOLの強みは、ユーザーが支配方程式を比較的自由にカスタマイズできる「弱形式」や「PDEインターフェース」です。これを用いて、論文などで提案された独自の損傷進化則

$$ \dot{D} = f(\sigma, \epsilon, D, T, ...) $$

を直接実装できます。組み込みの物理インターフェースとしては、「固体力学」モジュールに「損傷」機能があり、粘弾性体の損傷や疲労損傷のモデリングが可能です。ただし、AbaqusやLS-DYNAのように産業界で徹底検証された多数の破壊モデルが最初から箱入り、というわけではないので、ユーザー自身の実装と検証責任が大きくなります。

トラブルシューティング

よくある収束問題と対策

🧑‍🎓

損傷を考慮した非線形解析で、ニュートン・ラフソン法が収束せずに中止されます。最初に確認すべき点は？

🎓

まず、**損傷の進化が急激すぎないか**を疑います。損傷進化則のパラメータ（特に指数s）が大きすぎると、1ステップでDが0から0.8などにジャンプし、剛性マトリックスの急激な変化で収束不能になります。対策として、初期のタイムステップを非常に小さく（例: 総時間の1e-5）設定し、自動増分法に任せます。また、Abaqusでは「STABILIZE」オプションで人工的な減衰を導入し、不安定な軟化挙動を数値的に安定化させる手法もあります。

🧑‍🎓

解析は最後まで回るのですが、破断が予想より極端に早く、または遅く発生します。原因は？

🎓

主に3つ考えられます。1. **材料定数の誤り**: 特に破断ひずみや損傷エネルギーが実験値と合っているか再確認。2. **三軸応力度の計算精度**: 要素形状が歪んでいると、三軸応力度

$$ \eta = \sigma_m / \sigma_{eq} $$

の計算が不正確になり、損傷進化が狂います。メッシュ品質を改善。3. **メッシュサイズ効果**: 先述の通り、エネルギー基準で調整していない場合、メッシュが粗いと破断が遅く、細かいと早くなりがち。メッシュコンバージェンススタディが必須です。

🧑‍🎓

要素が消去された後、その空洞部分を囲む要素が異常な変形を起こし、結果が物理的に不自然です。

🎓

これは「要素消去に伴う接触問題」です。消えた要素の両側にあった面同士が、互いに貫通してしまう現象です。対策として、1. **自己接触条件の設定**: 破壊が発生する可能性のある領域全体に、あらかじめ通用面自己接触を定義する（Abaqusの「General Contact」など）。2. **剛性減衰法の採用**: 要素を完全に消去するのではなく、剛性を極小値まで減衰させることで、幾何学的な「すき間」が生じないようにする。3. **コヒーシブゾンモデルの併用**: 界面に初期からコヒーシブ要素を挿入し、そこでの分離として破壊を表現する。これらはLS-DYNAやAbaqus/Explicitで標準的に対処法が用意されています。