オイラー座屈 — 数値解法と実装

座屈問題のFEM定式化は、通常の静解析とは本質的に異なる。普通の静解析は $[K]\{u\} = \{F\}$ を解くだけだけど、座屈解析は固有値問題になるんだ。

そう、振動解析と数学的には同じ構造だ。振動では $([K] - \omega^2 [M])\{\phi\} = \{0\}$ だけど、座屈では：

ここで $[K]$ は通常の剛性マトリクス、$[K_\sigma]$ は幾何剛性マトリクス（応力剛性マトリクス）、$\lambda$ が座屈荷重係数、$\{\phi\}$ が座屈モード形状だ。

通常の剛性マトリクス $[K]$ は材料の弾性と形状から決まる。一方 $[K_\sigma]$ は現在の応力状態に依存する。圧縮応力があると $[K_\sigma]$ は負になり、全体の剛性を下げる方向に働く。

完璧な理解だ。$[K] + \lambda [K_\sigma]$ が特異になる（行列式がゼロになる）ときの $\lambda$ が座屈荷重係数。参照荷重に $\lambda$ をかけたものが臨界荷重になる。

梁要素の場合で考えよう。軸力 $N$ を受ける梁要素の幾何剛性マトリクスは：

そう。これが意味するのは、座屈解析ではまず静解析で応力分布を求めてから、その応力で $[K_\sigma]$ を組み立てて固有値問題を解く、という2段階の手順になるということだ。

その通り。Nastranなら SOL 105、AbaqusならBuckle ステップの前にStatic ステップを置く。この前段の静解析を忘れると、応力状態が不明なので $[K_\sigma]$ が作れない。

もちろん。座屈解析で使われる主な手法を整理しよう。

NastranのEIGRL、AbaqusのLanczosオプション、AnsysのBlock Lanczos法…主要ソルバーはすべてLanczos系を実装している。ポイントは求めるモード数の指定だ。1次モードだけでなく、少なくとも5〜10モードは求めておくべき。

2つ理由がある。1つは、1次と2次の固有値が近い場合（クラスタリング）、モード間の相互作用で実際の座屈挙動が変わる。もう1つは、局所座屈モード（フランジの座屈など）が全体座屈モードより低い固有値を持つことがあり、それを見逃すと危険だからだ。

整理するとこうだ：

計算コストが桁違いだ。線形座屈解析は固有値問題1回で済むけど、非線形解析は荷重増分ごとに平衡反復を回す。100倍以上の計算時間がかかることもある。

実務的なアプローチは：

1. まず線形座屈解析で座屈荷重の目安とモード形状を把握

2. 必要に応じて、線形座屈のモード形状を初期不整として与えて非線形解析を実施

3. 非線形解析では弧長法（Riks法）を使って、荷重のピーク（限界点）を通過した後座屈パスも追跡

通常のNewton-Raphson法は荷重を制御量にするから、荷重-変位曲線のピーク（限界点）を通過できない。Riks法は弧長をパラメータにして、荷重と変位を同時に変化させる。これにより荷重の折り返し（スナップスルー、スナップバック）も追跡できるんだ。

非常に重要だ。座屈モード形状を正しく表現するには、座屈波長に対して十分なメッシュ密度が必要になる。目安として：

• 座屈の半波長あたり最低4〜6要素（2次要素の場合）
• 薄板・シェル構造では板厚方向にも注意（ソリッド要素なら最低2層）
• 局所座屈を捕捉するには、フランジやウェブに十分な要素数

まさにそこが実務の難しさだ。まず粗めのメッシュで線形座屈解析を回し、モード形状を確認してから、座屈変形が集中する部位のメッシュを細分化して再解析する。メッシュ収束性の確認は必須だよ。

その通り。もう一つ注意点がある。座屈解析では要素タイプの選択が結果に大きく影響する。特に薄肉構造では、ソリッド要素よりシェル要素のほうが効率的かつ正確なことが多い。ソリッド要素で薄板をモデル化すると、シアロッキングで剛性が過大評価され、座屈荷重も過大になる。

要点をまとめると：

• 座屈解析は固有値問題 — 幾何剛性マトリクス $[K_\sigma]$ がカギ
• 2段階の手順 — まず静解析で応力を求め、次に固有値解析
• Lanczos法が実務標準 — 複数モードを必ず確認すること
• 線形→非線形の段階的アプローチ — 線形で概要把握、非線形で精密評価
• Riks法（弧長法）で後座屈まで追跡可能
• メッシュは座屈波長を解像できる密度が必要

まさにそう。振動の質量マトリクス $[M]$ が座屈の幾何剛性マトリクス $[K_\sigma]$ に入れ替わっただけ。この対応関係を掴んでおくと、ソルバーのマニュアルを読むときにも理解が早くなるよ。