オイラー座屈 — トラブルシューティングガイド

座屈解析は他の解析以上にトラブルが多い分野だ。よくある問題を症状別に整理しよう。

負の固有値は参照荷重の反対方向で座屈が起きることを意味する。つまり圧縮荷重として与えたのに引張方向の座屈が出ている状態だ。

原因と対策：

そう。ただし実務では正の最小固有値が重要だ。10モード求めて全部負だったら、荷重方向を完全に間違えている可能性が高い。

FEMが手計算より高い場合、まず以下を確認しよう：

1. 座屈方向は合っているか？

十分あり得る。境界条件によっては弱軸方向の変位を拘束してしまっていて、強軸座屈しか出ないことがある。例えば壁に面した柱を壁方向（弱軸方向）で拘束すると、弱軸座屈は抑制される。意図してそうしたなら正しいが、モデル化ミスなら深刻だ。

2. 境界条件の固定度は正しいか？

3. 要素タイプは適切か？

その場合は局所座屈が先行している可能性が高い。オイラー式は全体座屈しか評価しないが、シェルやソリッドでモデル化するとフランジの局所座屈が出ることがある。モード形状が全体の横たわみではなく、局部的な波打ちになっていないか確認してくれ。

これは非線形座屈で最もよくあるトラブルだ。原因は大きく3つ。

完璧に対称な構造に対称荷重をかけると、数値的に「分岐点」を通過できない。小さな摂動がないと、ソルバーは真っ直ぐな平衡経路に留まろうとして、分岐点で経路が分からなくなる。

対策: 線形座屈の1次モード形状を初期不整として導入する。振幅は部材長の1/1000〜1/500程度が一般的。Abaqusなら *IMPERFECTION で簡単にできる。

座屈点付近は荷重-変位曲線の勾配が急激に変化する。増分幅が大きいと、この急変に追従できない。

対策: 初期増分を小さくする。Riks法なら初期弧長を荷重の1%程度に設定。Abaqusの *STATIC, RIKS で 0.01, 1.0, 1e-8, 0.1 のように初期増分を小さく、最小増分をさらに小さく。

座屈で大変形が起きると、接触条件やバネの状態が急変して収束を乱すことがある。座屈解析のデバッグではまず境界条件を極力単純化（理想ピン・固定で）して収束を確認してから、段階的に複雑にするのが鉄則だ。

メッシュ収束性の確認が必須だ。手順はこう：

その通り。全体座屈は比較的メッシュに鈍感だけど、局所座屈（板の波状座屈など）は座屈半波長に対して要素が4つ以上ないと正しく表現できない。粗いメッシュでは局所座屈モードが「見えない」状態になる。

特に注意すべきは板の自由辺付近だ。自由辺の座屈はメッシュ密度に非常に敏感で、辺に沿って十分な要素数が必要。

経験則から言うと、この順番で調べるのが最速だ：

Step 1: モード形状を目で見る

Step 2: 最も単純なモデルで検証する

Step 3: 段階的に複雑化する

座屈解析は「ソルバーが何かおかしい」のか「自分の理解が足りないのか」の切り分けが難しい分野だ。最小モデルでオイラーの理論解を再現できれば、少なくともソルバーの使い方は合っていると確信できる。そこから一歩ずつ積み上げるのが、遠回りに見えて最速の道だよ。