プログレッシブ損傷解析
理論と物理
プログレッシブ損傷とは
先生、「プログレッシブ損傷解析（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における最終破壊荷重予測を可能にし、設計の安全余裕の定量化に使われる。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
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%影響するため、実験キャリブレーションが精度の鍵だ。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
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用タンクはこの解析で設計最適化されている。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
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構造の設計解析に採用されている。
選定で最も重要な3つの問い
- 「何を解くか」：プログレッシブ損傷解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
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%未満であることを確認する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——プログレッシブ損傷解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告