非線形後座屈解析 — トラブルシューティングガイド

後座屈は非線形解析の中でも最も収束が難しい分野だ。よくあるトラブルを症状別に整理しよう。

最も多いトラブルだ。原因チェック：

1. 初期不整が入っていない（または極めて小さい）

完璧に対称な構造では、分岐点の手前で接線剛性がゼロに近づき、増分がどんどん小さくなる。初期不整を導入して対称性を崩すこと。

2. 最小増分が大きすぎる

最小増分（Abaqusの *STATIC, RIKS の3番目のパラメータ）を $10^{-10}$ 程度まで小さくしてみる。

3. 最大増分数が少なすぎる

デフォルトの増分数（100程度）では足りないことがある。1000〜10000に増やす。

その場合は安定化法（*STATIC, STABILIZE）に切り替えるのが現実的だ。Riks法は万能ではなく、複雑な座屈モード（複数の限界点が連鎖するケースなど）では追跡に失敗することがある。安定化法はRiksより頑健だ。

Riks法は最初に見つけた分岐に沿って進む。意図した座屈モードではなく、別の安定な経路に乗ってしまうことがある。

対策：

• 初期不整のモード選択を見直す — 目的のモードに対応する不整を与えているか
• 不整振幅を大きくする — 小さすぎると分岐点を「乗り越えて」しまうことがある
• 初期増分を小さくする — 分岐点付近の挙動を細かく追跡

まさにそう。不整が大きいと鋭い分岐が滑らかなカーブに変わり、Riks法が追跡しやすくなる。ただし不整が大きすぎると物理的に非現実的になるから、振幅の感度分析が重要だ。

人工粘性が大きすぎる。安定化エネルギーが実際の変形エネルギーを上回っている状態だ。

確認手順：

1. Abaqusの場合、History Outputで ALLSD（安定化散逸エネルギー）と ALLIE（内部エネルギー）を出力

2. ALLSD / ALLIE < 0.05 (5%) であることを確認

3. 超えていたら安定化係数を下げる（*STATIC, STABILIZE, FACTOR=0.0001 など）

Abaqusのデフォルトの安定化係数は自動計算で、多くの場合は適切だが、座屈問題では大きすぎることがある。特に急激なスナップスルーでは粘性が効きすぎて、本来の荷重ドロップが抑制される。必ずエネルギー比で検証すること。

その感覚は正しい。後座屈では以下の理由でメッシュ依存性が強くなる：

• 局所座屈の波長がメッシュで解像できていない — 後座屈ではモード遷移が起きるため、複数のモード波長を正しく捕捉する必要がある
• 塑性ヒンジの局所化 — 弾塑性座屈では変形がメッシュ1要素に集中し、メッシュ依存の結果になる

対策：

• 座屈波長に対して最低6要素（2次要素の場合）
• 弾塑性の局所化が起きる場合は正則化手法（粘塑性モデル、非局所モデル）を検討
• 最低3水準のメッシュで崩壊荷重の収束を確認

後座屈解析でエネルギーバランスが崩れる主な原因：

鉄則として、後座屈解析では全てのエネルギー成分をHistory Outputに出力しておくこと。何か問題があったとき、エネルギーのどの成分が異常かで原因が特定できる。

1. まず線形座屈が正しいか確認 — 非線形の前段が間違っていたら、全部狂う

2. 初期不整の有無と振幅を確認 — 入れ忘れ、または不適切な振幅が最多の原因

3. 収束しない場合 — 最小増分を下げる → 安定化法に切り替え → 陽解法を検討

4. 結果がおかしい場合 — モード形状を可視化 → エネルギーバランスを確認 → メッシュ収束性

5. 最小モデルで再現テスト — 1枚の板や1本の柱で、理論解と合うことを確認してから実モデルへ

その通り。後座屈解析はパラメータが多くて問題の切り分けが難しい。だからこそ、一つずつ変数を固定して段階的に複雑化するアプローチが最速の解決策なんだ。