非線形後座屈解析 — 実践ガイドとベストプラクティス

典型的なフローはこうだ：

1. 線形座屈解析 — 座屈荷重の概算とモード形状の取得

2. 初期不整の導入 — 1次モード形状を初期不整として与える

3. 非線形解析（Riks法または安定化法） — 荷重-変位経路を追跡

4. 崩壊荷重の特定 — 荷重のピーク、または許容変形量を超えた点

5. 感度確認 — 不整振幅を変えて崩壊荷重の感度を確認

6. 結果の検証 — メッシュ収束性、エネルギーバランス

一般的には板厚 $t$ の 10%、50%、100%、200% の4〜5ケース。グラフにすると不整振幅 vs. 崩壊荷重の曲線（ノックダウン曲線）が得られる。この曲線が急峻なら不整敏感性が高く、なだらかなら低い。

Abaqusが後座屈で最もよく使われるので、詳しく見よう。

順番に：初期増分, 全弧長, 最小増分, 最大増分。

• 初期増分 0.01 — 参照荷重の1%から始める。小さめが安全
• 全弧長 1.0 — 参照荷重の100%まで（実際はピーク後も自動で続く）
• 最小増分 1e-10 — 座屈点近傍で極小の増分を許容
• 最大増分 0.05 — 1回の増分で5%以上進まないように制限

座屈点近傍では荷重-変位曲線の曲率が非常に大きくなるため、増分を極端に小さくしないと追跡できないことがある。ただし実際に $10^{-10}$ まで小さくなったら、それは収束困難な兆候だから、不整振幅やメッシュを見直すべきだ。

重要な違いが2つある。

1. 荷重が出力されない — 通常の解析では荷重は入力値だが、Riks法では荷重も未知数。各増分での実際の荷重は $\lambda \times F_{ref}$ で、$\lambda$ はLPF（Load Proportionality Factor）として出力される。

2. 増分番号が時間に対応しない — Riks法の「時間」は弧長に対応する仮想的な値。物理的な時間とは関係ない。

History OutputでLPFと注目点の変位を出力し、LPF × 参照荷重を横軸（または縦軸）、変位を縦軸（または横軸）にプロットする。Abaqusの結果ファイル（.odb）からPythonスクリプトで抽出するのが一般的だ。

これは明確な定義がない場合もある。一般的な基準：

ユーロコードでは終局限界状態（ULS）として「LPFのピーク」を用いることが多い。航空宇宙ではMIL-HDBK-5等で許容ひずみが規定されていて、「ひずみが規定値を超えた時点」を崩壊とする場合もある。

後座屈はベンチマーク問題が限られるので、検証が難しい。以下の多角的なアプローチが必要だ：

1. エネルギーバランスの確認 — 外力仕事 ≈ 内部ひずみエネルギー + 散逸エネルギー。安定化を使っている場合、ALLSD/ALLIE < 5%

2. メッシュ収束性 — 2水準以上のメッシュで崩壊荷重の変化が5%以内

3. 不整敏感性の妥当性 — ノックダウン曲線が物理的に妥当か（不整が増えたら荷重は下がるはず）

4. 変形パターンの物理的妥当性 — 座屈モード形状が予想と一致するか

5. 既知の実験データとの比較 — 可能であれば最も信頼性が高い

そう。特にエネルギーバランスとメッシュ収束性は最低限チェックすべき項目だ。この2つが合っていなければ、結果は信用できない。