プログレッシブ損傷解析 — トラブルシューティングガイド

PDAはFEMの中でも最も複雑な解析の一つ。トラブルも多岐にわたる。

確認項目（優先度順）：

1. 破壊エネルギー $G_c$ は正しいか — 最も影響が大きい。文献値ではなく試験値を使うべき

2. メッシュサイズは適切か — 0.5〜2 mmの要素。粗すぎると損傷が広がりすぎる

3. 材料強度は適切か — B-basis vs. 平均値で20%の差

4. 積層の定義は正しいか — 繊維角、積層順序、材料座標系

5. 境界条件は試験と一致するか — クランプの拘束、荷重の作用位置

応力集中による局所化。対策：

• 破壊エネルギーの正則化が正しく機能しているか確認
• メッシュを細かくして損傷が複数要素に分散するか確認
• 粘性正則化のパラメータを調整

要素削除の基準が厳しすぎるか、損傷モデルが不安定。

対策：

• 要素削除の損傷閾値を $d = 0.99$ ではなく $d = 0.999$ に（完全損傷に近いときだけ削除）
• 繊維損傷でのみ要素削除する（マトリクス損傷では削除しない）
• 最大ひずみベースの削除基準を追加（異常変形の防止）

PDAは非線形解析＋損傷更新＋接触（CZMの場合）が組み合わさるため、計算コストが大きい。

効率化：

• 陽解法＋質量スケーリング — 収束問題を回避しつつ計算を高速化
• 対称条件の活用 — 対称な問題なら1/2, 1/4モデル
• 局所的なPDA — 全体モデルは弾性、損傷予想領域のみPDAを適用
• MPI並列計算 — 計算ノード数を増やす

• 最終破壊荷重のずれ → $G_c$、メッシュサイズ、材料強度を確認
• 損傷の局所化 → 正則化の確認、メッシュの細分化
• 大量の要素削除 → 削除基準の緩和、繊維損傷でのみ削除
• 計算速度 → 陽解法＋質量スケーリング、局所PDA
• 試験との比較が全ての検証の基本 — 荷重-変位曲線の一致

PDAの精度は材料パラメータの品質に直結する。試験データなしにPDAを行うのは「地図なしで登山する」ようなものだ。