スナップバック解析 — トラブルシューティングガイド

スナップバックは非線形解析で最も収束が困難な問題だ。

これはスナップバック点でRiks法の弧長方向が定まらない状態だ。

対策（段階的に試す）：

1. 最大増分を小さくする — 0.01 → 0.005 → 0.001

2. 最小増分をさらに小さくする — $10^{-12}$ → $10^{-15}$

3. cylindrical arc-length に切り替える — 変位空間のみで弧長制約

4. 安定化法に切り替える — *STATIC, STABILIZE

5. 動的解析に切り替える — *DYNAMIC, QUASI-STATIC

通常のRiks法の設定調整（1, 2）で解決することも多い。それでダメなら手法の切り替え（3〜5）を検討する。一つの手法に固執せず、柔軟に切り替えるのがスナップバック解析のコツだ。

安定化係数が大きすぎる。人工粘性がスナップバックのエネルギー解放を吸収してしまっている。

対策：

• 安定化係数を下げる（FAQTORを0.001 → 0.0001 → 0.00001と段階的に）
• ALLSD/ALLIE比を確認（1%以下を目標）
• Riks法の結果と比較検証

そう。安定化係数の「スイートスポット」を見つける必要がある。小さすぎると不安定で収束しない、大きすぎるとスナップバックが抑制される。パラメトリックスタディで適切な範囲を特定するのが実務的だ。

スナップバック時のエネルギー解放が運動エネルギーに変換されて振動が起きる。対策：

• 減衰を追加 — Rayleigh減衰（α, β）またはAbaqusの*DYNAMIC, APPLICATION=MODERATE DISSIPATIONで数値減衰を増やす
• 段階的荷重除荷 — スナップバック点の直前で荷重を一度除荷し、再載荷
• 質量スケーリング — 質量を増やして固有振動数を下げ、時間増分を大きくする（ただし慣性効果に注意）

局所化に対する正則化が不足している。

確認と対策：

• 亀裂帯モデルを使っているか — 軟化曲線をメッシュサイズに応じて調整
• 特性長さ（characteristic length）は正しいか — Abaqusでは自動計算されるが、特殊な要素形状では不正確になることがある
• 破壊エネルギー $G_f$ が一定か — メッシュを変えても $G_f$ が保存されているか確認

理想的にはそう。ただし実際には完全にメッシュ独立にはならないことが多い。2〜3水準のメッシュで結果を比較し、差が許容範囲内（10〜20%）であることを確認するのが実務的だ。

• Riks法の振動 — 増分を小さく → cylindrical arc-length → 安定化 → 動的
• 安定化でスナップバック消失 — 安定化係数を下げる。ALLSD/ALLIE < 1%
• 動的解析の振動 — 減衰追加。MODERATE DISSIPATION
• 破壊のメッシュ依存性 — 亀裂帯モデルで正則化。$G_f$ 保存を確認
• 一つの手法に固執しない — 複数の手法を試して結果を比較

その通り。スナップバックを正しく追跡できたら、非線形構造力学の上級者だ。焦らず、一つずつ手法を試して経験を積んでいこう。