力の残差発散 — トラブルシューティング

どの荷重増分から発散が始まるか特定する:

• 最初の増分で発散 → 初期条件、拘束条件、材料定数の問題
• 途中の増分で発散 → 非線形現象（座屈、接触変化、塑性局所化）
• 特定の荷重レベルで発散 → 物理的な限界（座屈荷重、破壊荷重）

最大残差が発生している節点と自由度を確認する。Abaqusでは.msgファイルに詳細が出力される。残差が大きい節点の周辺を重点的にチェック:

• その節点の拘束条件は適切か
• 材料定義に問題がないか
• メッシュ品質は十分か

モデルを段階的に簡略化して問題を切り分ける:

1. 線形解析で結果の妥当性を確認

2. 非線形性を1つずつ追加（大変形→材料非線形→接触）

3. 荷重を小さくして線形範囲で確認

4. メッシュを粗くして基本挙動を確認

モデルが物理的に正しい場合の数値安定化:

• 粘性正則化（Viscous Regularization）: 材料の急激な軟化に対して*DAMAGE STABILIZATIONで粘性を追加
• 自動安定化: *STATIC, STABILIZE=DISSIPATED ENERGYで散逸エネルギーベースの安定化
• Arc-Length法: *STATIC, RIKSで荷重制御から弧長制御に切り替え

荷重-変位曲線にスナップスルーやスナップバックがある構造（薄肉シェルの座屈、ゴム部品のスナップフィット等）で使う。通常のNewton-Raphson法は荷重を制御パラメータにするが、Arc-Length法は荷重と変位を同時に制御パラメータにするため、荷重が減少する経路も追跡できる。