線形座屈（固有値座屈）解析 — トラブルシューティングガイド

固有値座屈特有のトラブルを整理しよう。オイラー座屈のトラブルシューティングとは違う、固有値解析ならではの問題がある。

これは前段の静解析で有意な応力が発生していない場合に起きる。$[K_\sigma]$ がほぼゼロ行列なので、固有値が極端に大きいか小さい値になる。

原因チェック：

1. 参照荷重がゼロまたは極小 — 荷重の値や方向を確認

2. 荷重が全て拘束点に作用している — 反力だけで平衡し、構造内部に応力が生じない

3. 剛体移動が拘束されていない — 静解析が特異になっている

4. Ansys APDLで PSTRES, ON を忘れている — 応力剛性が計算されない

そう。Nastranの SOL 105 やAbaqusの *BUCKLE は自動的に応力剛性を計算するが、AnsysのAPDLでは PSTRES, ON を明示的に指定しないと幾何剛性マトリクスが計算されない。Workbenchでは自動だが、APDLスクリプトでは忘れがちだ。

参照荷重の方向が意図と逆になっている可能性が高い。例えば圧縮したいのに引張荷重を与えている場合だ。

そう。あるいは、参照荷重はそのままで $\lambda$ を負で読むこともできる。$\lambda = -3.0$ は「反対方向の3倍の荷重で座屈」という意味だ。ただし解釈が紛らわしいので、最初から荷重方向を合わせる方が安全だ。

これが固有値のクラスタリングだ。物理的には「ほぼ同時に2つの座屈モードが発生する」状況で、モード相互作用による耐力低下が懸念される。

対処法：

• モード形状を比較 — 同じ形ならメッシュの対称性による重複（問題なし）。異なる形なら要注意
• 対称モードと反対称モード — 対称構造では対称/反対称がペアで出やすい
• 非線形解析で検証 — 両モードを初期不整として重ね合わせ、相互作用の影響を確認

$\{u_{impf}\} = \alpha_1 \{\phi_1\} + \alpha_2 \{\phi_2\}$ として、$\alpha_1, \alpha_2$ の比率を変えた複数ケースを実施する。例えば $\alpha_1 : \alpha_2 = 1:0, 1:1, 0:1$ の3ケースが最低限。最も低い崩壊荷重を設計値とする。

これは薄肉構造でよくある問題だ。フランジやウェブの板要素が局所的に座屈するモードが、全体座屈より低い固有値を持つことがある。

対策：

1. 求めるモード数を大幅に増やす（50〜100モード） — 全体モードが上位に埋もれている

2. シフト-反転法で特定の固有値域を探索 — 全体座屈の概算値付近にシフトを設定

3. サブモデリング — 全体モデル（梁要素）で全体座屈を評価、局所モデル（シェル要素）で局所座屈を別途評価

その通り。大規模な構造（ビル全体、橋梁全体など）では、全体を詳細なシェルモデルにするのは現実的でない。梁モデルで全体座屈を確認し、クリティカルな部材だけシェルで局所座屈を確認するマルチスケールアプローチが実務的だ。

プリロードによる応力が $[K_\sigma]$ に反映されるからだ。自重で柱に圧縮が入っていれば、追加荷重に対する座屈余裕は当然小さくなる。

逆に、プリロードが引張（例：プレストレス）だと座屈荷重が上がる。プレストレスコンクリート柱が座屈に強い理由の一つはこれだ。

必ずしも過大評価とは限らない。プリロードが引張方向に効く部材もあるから。ただし自重のような永続荷重を無視する理由はないので、常にプリロードを含めるのが安全だ。

NAFEMSが公開しているベンチマーク問題が最も信頼性が高い：

そう。新しいソルバーや新しい要素タイプを使い始めるときは、必ずベンチマーク検証をしてから実問題に入ること。座屈解析は設定のミスが結果に直結する分野だから、検証のステップを省いてはいけない。

その通り。数字だけ見て判断するのではなく、なぜその固有値が出たのかを常に考える習慣をつけてほしい。そうすればトラブルの原因は自然と見えてくる。