プログレッシブ損傷解析
プログレッシブ損傷の理論基礎
プログレッシブ損傷とは
先生、「プログレッシブ損傷解析（PDA）」って何ですか？
複合材は金属と違って局所的な破壊が起きても構造全体は崩壊しないことがある。マトリクスクラックが入っても繊維が荷重を担い続ける。この段階的な損傷の進展と荷重の再配分をシミュレーションするのがPDAだ。
破壊判定（Tsai-Wu, Hashin）だけでは不十分ということですか？
Tsai-WuやHashinは「初期破壊（First Ply Failure）」を予測する。しかし複合材構造の最終破壊荷重（Last Ply Failure）は初期破壊荷重の2〜3倍になることがある。PDAはこの「初期破壊から最終崩壊まで」の全過程を追跡する。
PDAの3要素
PDAは3つの要素で構成される：
1. 損傷開始基準（Damage Initiation）
「いつ損傷が始まるか」。Hashin基準、Puck基準、LaRC基準など。各破壊モードの開始を判定。
2. 損傷進展則（Damage Evolution）
「損傷がどう進むか」。損傷開始後の剛性低減をどう表現するか。
- 瞬間的低減 — 損傷開始で剛性を一気にゼロに（突然劣化モデル）。シンプルだがメッシュ依存性大
- 漸進的低減 — 破壊エネルギーに基づいて徐々に剛性を低減。メッシュ依存性が小さい
3. 要素削除（Element Deletion）
損傷変数が1.0（完全損傷）に達した要素を解析から除去。材料が完全に破壊した領域を表現。
漸進的低減のほうが物理的に正確ですか？
そう。破壊エネルギー $G_c$ を使った漸進的低減は、メッシュサイズに依存しない（正則化された）結果を出す。AbaqusのHashin Damage Evolutionがこの手法だ。
損傷変数
損傷の程度を損傷変数 $d$（0〜1）で表現：
- $d = 0$: 健全（損傷なし）
- $0 < d < 1$: 部分的損傷
- $d = 1$: 完全破壊
Hashinの4モードそれぞれに損傷変数があるんですか？
そう。繊維引張 $d_{ft}$、繊維圧縮 $d_{fc}$、マトリクス引張 $d_{mt}$、マトリクス圧縮 $d_{mc}$ の4つの独立した損傷変数が定義される。
損傷後の剛性マトリクス：
各モードの損傷が独立に進行し、剛性マトリクスの対応する成分が低減する。
これがCDM（Continuum Damage Mechanics）ベースのPDAだ。連続体力学の枠組みで損傷を扱うため、FEMの標準的な枠組みに自然に組み込める。
まとめ
PDAの理論を整理します。
要点：
- 初期破壊から最終崩壊までの全過程をシミュレーション — FPFからLPFまで
- 3要素 — 損傷開始基準 + 損傷進展則 + 要素削除
- 損傷変数 $d$（0〜1） — 各破壊モードの損傷度を追跡
- CDM（連続体損傷力学）ベース — FEMに自然に組み込み可能
- 破壊エネルギーによる正則化 — メッシュ依存性の排除
PDAは「複合材の強度を最後まで搾り出す」ための解析手法ですね。
複合材の後荷重強度（post-failure strength）を評価できるPDAは、軽量化の限界を押し上げるための不可欠な技術だ。
Fracture Mechanicsとプログレッシブ破損の統合
複合材のプログレッシブ破損（Progressive Damage）とは荷重増加に伴い局所的な材料破損が順次拡大していく過程だ。1990年代にChaboche・Ladevezらが連続体損傷力学（CDM）をCFRPに適用し、損傷変数d（0=未損傷・1=完全破損）の発展方程式を定式化した。CDMベースのプログレッシブ解析はFEMにおける最終破壊荷重予測を可能にし、設計の安全余裕の定量化に使われる。
プログレッシブ損傷の数値計算手法
PDAの実装
PDAの具体的な実装方法を教えてください。
2つの主要な実装アプローチ：
1. Abaqusの組み込みHashin Damage
AbaqusのDAMAGE INITIATION (HASHIN) + DAMAGE EVOLUTION。前のページで説明した設定で、損傷開始→漸進的低減→要素削除の全プロセスが自動で動く。
2. ユーザーサブルーチン（UMAT/VUMAT）
より高度な損傷モデル（Puck、LaRC05、CDMカスタム）を使いたい場合は、ユーザーサブルーチンを自作する。Abaqusの場合：
VUMATのほうが実装しやすいんですか？
VUMATは接線剛性マトリクスの導出が不要で、与えられたひずみ増分から応力を更新するだけ。PDAの初回実装にはVUMAT + Explicit が推奨だ。陰解法のUMATは接線剛性の導出が難しく、収束性にも影響する。
収束性の問題
PDAは収束が難しいんですよね。
損傷による剛性低減は局所的な軟化を引き起こし、荷重-変位経路がスナップバックすることがある。対策：
| 手法 | 特徴 |
|---|---|
| 粘性正則化 | 微小な粘性で軟化を「なまらせる」。$\eta \approx 10^{-5}$ |
| 陽解法 | 収束問題なし（反復なし）。計算コスト大 |
| 弧長法（Riks） | スナップバックを追跡可能。設定が複雑 |
| 安定化法 | *STATIC, STABILIZE。エネルギー比で検証 |
陽解法が最も安定ですか？
収束の心配がないという意味では最も安定。ただし準静的問題を陽解法で解くには質量スケーリングが必要で、慣性効果に注意。実務では陽解法＋質量スケーリングがPDAの主流になりつつある。
LS-DYNAでのPDA
LS-DYNAではMAT54（Chang-Chang基準）とMAT58（連続損傷力学）がPDAの標準：
| モデル | 特徴 |
|---|---|
| MAT54 | 突然劣化モデル。シンプルだがメッシュ依存性大 |
| MAT58 | CDMベース。漸進的低減。MAT54より正確 |
自動車の衝突解析ではMAT54/58が広く使われているんですよね。
CFRPのクラッシュボックスやバンパーの衝突解析はLS-DYNA MAT54/58が事実上の標準。衝突時のエネルギー吸収量を予測する。
まとめ
PDAの数値手法、整理します。
要点：
- Abaqusの組み込みHashin Damageが最も手軽 — 設定のみで動作
- VUMAT（陽解法）が自作PDAの推奨 — 接線剛性不要
- 陽解法＋質量スケーリングがPDAの主流 — 収束問題を回避
- LS-DYNA MAT54/58が自動車業界の標準 — 衝突時のCFRP破壊
- 粘性正則化でメッシュ依存性を軽減 — $\eta \approx 10^{-5}$
剛性低減と損傷変数の更新スキーム
プログレッシブ損傷解析では①荷重増分②Hashin等の破損基準で初期破損判定③破損モードに応じた剛性低減（e.g. Eidax法でE11を20%・E22を90%低減）④変形更新⑤残余応力再分配⑥次ステップの順で反復する。剛性低減の速度（突然低減 vs 指数的低減）が最終破壊荷重に20〜40%影響するため、実験キャリブレーションが精度の鍵だ。
プログレッシブ損傷の実務適用
PDAの実務適用
PDAは実務でどう使われていますか？
航空宇宙と自動車で主要な適用がある。
航空宇宙: OHT/OHC解析
Open Hole Tension/Compression（穴あき板の引張/圧縮）はPDAの最も基本的な適用。ASTM D5766（OHT）, D6484（OHC）の試験をFEMで再現する。
手順：
1. 穴あき板のFEMモデル（シェル＋Hashin損傷）
2. 端面に変位制御で引張/圧縮
3. 荷重-変位曲線のピークが最終破壊荷重
4. 損傷の進展パターン（繊維破壊、マトリクスクラック）を可視化
5. 試験結果と比較
OHT/OHCの試験とFEMの一致精度はどの程度ですか？
適切にキャリブレーションされたPDAモデルで10〜15%以内の精度が達成可能。ただし破壊エネルギーのキャリブレーションが必要。「ブラインド予測」（キャリブレーションなし）では20〜30%の誤差が出ることがある。
自動車: CFRPクラッシュ解析
自動車のCFRPクラッシュボックスの衝突解析。エネルギー吸収量（SEA: Specific Energy Absorption）を予測する。
CFRPのクラッシュは金属と全く異なる：
- 金属: 塑性変形で座屈 → 安定したエネルギー吸収
- CFRP: 破砕（fragmentation）→ 不安定だが高い比エネルギー吸収
CFRPのクラッシュは「壊れ方」が複雑なんですね。
実務チェックリスト
PDAのチェックリストをお願いします。
「荷重-変位曲線の試験との比較」が最終的な検証ですね。
PDAの目的は「最終破壊荷重の予測」。荷重-変位曲線が試験と一致するかどうかが、PDAモデルの信頼性を示す唯一の指標だ。
炭素繊維プレッシャーベッセルの最終破壊
タイプIV圧力容器（CFRP強化・ポリマーライナー）はCNG（圧縮天然ガス）車両に使われ、プログレッシブ損傷解析で最終破裂圧力を予測する。70MPa（Hワンパス水素タンク）の最終破壊解析では繊維破断（σ11/XT≥1）の発生点から容器全体の破裂までの経路をFEMで追跡し、設計破裂圧力（使用圧力の2.25倍以上）を確保する。トヨタFCEVのMirai用タンクはこの解析で設計最適化されている。
プログレッシブ損傷のソフトウェア比較
PDAのツール
PDAに使えるツールを比較してください。
航空宇宙はAbaqus、自動車はLS-DYNA。
この棲み分けは複合材PDAで特に明確。Abaqusの論文実績とLS-DYNAの衝突実績は各分野で圧倒的だ。
専用ツール
| ツール | 特徴 |
|---|---|
| CompDam（NASA） | Abaqus VUMAT。LaRC05基準＋CDM。NASAが開発・公開 |
| Helius:MCT（Autodesk/旧Firehole） | 多くのソルバーに対応。MCT（Multi-Continuum Theory）ベース |
| Genoa（AlphaSTAR） | マルチスケールPDA。ミクロ→マクロの連成 |
CompDamはNASAが公開しているんですか？
そう。CompDamはNASAランリー研究所が開発したAbaqus VUMATで、GitHubで公開されている。LaRC05基準ベースのPDAモデルで、研究レベルの最先端がオープンソースで利用可能。
選定ガイド
CompDamが無料で使えるのは素晴らしいですね。
NASAの複合材研究成果がオープンソースで公開されていることは、学術・産業の両方にとって大きな財産だ。
Helius Composite（Autodesk）の自動化破損解析
Autodesk Helius Compositeは複合材の多スケール解析とプログレッシブ損傷解析を統合したソフトで、Abaqus・Nastran・ANSYSのプラグインとして機能する。繊維・マトリクスレベルの破損をマクロスケールの損傷変数にスケールアップする「Multi-continuum Theory（MCT）」を使い、初期破損から最終破壊まで自動計算する。Lockheed Martinの宇宙機CFRP構造の設計解析に採用されている。
プログレッシブ損傷の先端研究
PDAの先端研究
PDAの最前線を教えてください。
3つの方向がある。
マルチスケールPDA
マイクロスケール（繊維-マトリクスのRVE）の破壊力学とマクロスケール（積層板）のPDAを連成する。マイクロスケールの破壊メカニズム（繊維/マトリクス界面の剥離、マトリクスの亀裂）がマクロスケールの損傷進展にどう影響するかを直接シミュレーションする。
計算コストが膨大そうですね。
FE²法やROMで効率化する。各積分点でRVEの応答を計算するFE²は最も正確だが、RVEの応答をニューラルネットワークで代替するData-Driven PDAが注目されている。
Phase-Field法によるPDA
Phase-Field法を複合材の面内損傷に適用する研究が急増。Hashinのようなモード分離を事前に仮定せず、亀裂の経路が自然に決定される。繊維方向と直交方向で異方的なPhase-Field定式化が開発されている。
デジタルツインとPDA
損傷が「いつ危険なレベルに達するか」をリアルタイムで予測できるんですね。
CBM（Condition-Based Maintenance）の基盤技術だ。従来の定期検査（TBM）から、状態ベースの保全に移行するための鍵がPDAベースのデジタルツインだ。
まとめ
PDAの先端研究、まとめます。
PDAは複合材構造の「寿命予測」の核心技術であり、航空宇宙から自動車までの構造安全性を支えている。
NLPFA：非線形プログレッシブ疲労破損
疲労荷重下での複合材プログレッシブ破損（NLPFA）は非常に複雑で、各サイクルの損傷蓄積と剛性低下を追跡する必要がある。Shokrieh・Lessard（2000年）が提案した疲労破損モデルは繰返し数ΔNごとに損傷変数を更新し、剛性低下曲線と最終破壊を予測できる。風力タービンブレードの設計寿命20年（10⁸サイクル以上）の予測にこの手法が2010年代から採用されている。
プログレッシブ損傷のトラブル対応
PDAのトラブル
PDAでよくあるトラブルを教えてください。
PDAはFEMの中でも最も複雑な解析の一つ。トラブルも多岐にわたる。
最終破壊荷重が試験と合わない
荷重-変位曲線のピークが試験値と大きくずれます。
確認項目（優先度順）：
1. 破壊エネルギー $G_c$ は正しいか — 最も影響が大きい。文献値ではなく試験値を使うべき
2. メッシュサイズは適切か — 0.5〜2 mmの要素。粗すぎると損傷が広がりすぎる
3. 材料強度は適切か — B-basis vs. 平均値で20%の差
4. 積層の定義は正しいか — 繊維角、積層順序、材料座標系
5. 境界条件は試験と一致するか — クランプの拘束、荷重の作用位置
損傷が1要素に集中する
損傷が特定の1要素にだけ集中します。
応力集中による局所化。対策：
- 破壊エネルギーの正則化が正しく機能しているか確認
- メッシュを細かくして損傷が複数要素に分散するか確認
- 粘性正則化のパラメータを調整
要素が大量に削除される
解析の途中で要素が大量に消えて、構造がバラバラになります。
要素削除の基準が厳しすぎるか、損傷モデルが不安定。
対策：
- 要素削除の損傷閾値を $d = 0.99$ ではなく $d = 0.999$ に（完全損傷に近いときだけ削除）
- 繊維損傷でのみ要素削除する（マトリクス損傷では削除しない）
- 最大ひずみベースの削除基準を追加（異常変形の防止）
計算が遅い
PDAの計算時間が非常に長いです。
PDAは非線形解析＋損傷更新＋接触（CZMの場合）が組み合わさるため、計算コストが大きい。
効率化：
- 陽解法＋質量スケーリング — 収束問題を回避しつつ計算を高速化
- 対称条件の活用 — 対称な問題なら1/2, 1/4モデル
- 局所的なPDA — 全体モデルは弾性、損傷予想領域のみPDAを適用
- MPI並列計算 — 計算ノード数を増やす
まとめ
PDAのトラブル対処、整理します。
PDAは「パラメータのキャリブレーション」が避けて通れないんですね。
PDAの精度は材料パラメータの品質に直結する。試験データなしにPDAを行うのは「地図なしで登山する」ようなものだ。
プログレッシブ解析で最終破壊前に解が発散する場合
プログレッシブ損傷FEM解析が途中で発散する場合、主原因は損傷要素での急激な剛性低減（突然低減法）だ。指数的低減法（E'=E×exp(-d/dc)）に変えるか、人工粘性（Viscous regularization）を追加することで収束性が大幅に改善する。Abaqusでのviscous stabilization係数η=1e-4〜1e-3が多くの複合材解析で有効な範囲だ。収束後には安定化エネルギーが全ひずみエネルギーの5%未満であることを確認する。
関連トピック
なった
詳しく
報告