線形座屈（固有値座屈）解析 — 理論と支配方程式

同じだ。FEMの文脈では「線形座屈解析」（linear buckling analysis）、「固有値座屈解析」（eigenvalue buckling analysis）、「線形化前座屈解析」（linearized pre-buckling analysis）は全て同じ手法を指す。基準状態の応力から幾何剛性マトリクスを作り、固有値問題を解いて座屈荷重係数を求める手法だ。

その通り。オイラー座屈は「柱」という特定の構造に限定した話だったけど、固有値座屈解析は板、シェル、フレーム、さらにはそれらの組み合わせ構造に適用できる汎用的な方法だ。

2段階の手順になる。

Step 1: 参照荷重 $\{F_{ref}\}$ による静解析を実行し、応力分布 $\{\sigma_{ref}\}$ を求める：

Step 2: 得られた応力から幾何剛性マトリクス $[K_\sigma]$ を構成し、固有値問題を解く：

そう。ここで重要なのは、Step 1の静解析が線形静解析であることだ。つまり大変形効果を考慮しない。荷重に比例して応力が増えるという仮定のもとで、「いつ不安定になるか」を予測している。

まさに。この「線形性の仮定」が固有値座屈解析の本質的な限界であり、同時に計算の速さの源泉でもある。

本質的な問いだね。座屈とは「荷重が増加しても変形が一定だった状態から、突然別の変形パターンに遷移する」現象だ。数学的には、平衡経路の分岐（bifurcation）と捉える。

荷重を $\lambda \{F_{ref}\}$ とスケーリングすると、全体の剛性は $[K_0] + \lambda [K_\sigma]$ になる。ここで $[K_\sigma]$ は参照応力に対する幾何剛性だから、$\lambda$ に比例してスケールする。この全体剛性が特異になる（det = 0）ときが分岐点で、その $\lambda$ が座屈荷重係数だ。

そう。振動解析の $([K] - \omega^2 [M])\{\phi\} = \{0\}$ と全く同じ数学構造だ。振動では質量マトリクス $[M]$ の位置に、座屈では幾何剛性マトリクス $[K_\sigma]$ が入る。

$[K_\sigma]$ は「現在の応力状態が、微小な変位摂動に対してどれだけ仕事をするか」を表す。圧縮応力下で横方向に少したわませると、圧縮力がたわみを増大させる方向に仕事をする。これが負の剛性として効く。

具体的に言うと、軸力 $N$ を受ける梁要素で、微小な横たわみ $\delta v$ が生じたとき、軸力は $N \cdot (\delta v')^2 / 2$ の追加仕事をする。この「応力による追加仕事」を行列化したものが $[K_\sigma]$ だ。

その直感は正しい。ただし注意点がある。例えばプレストレスされたケーブルのような引張部材でも、横方向の座屈（本質的には引張材の振動モードに近い現象）が問題になることはある。基本的には「圧縮応力が支配する部分が座屈の起点」と理解しておけばいい。

以下の前提条件が満たされているときに信頼性が高い：

1. 座屈前の変形が微小 — 形状が実質的に変わらない状態

2. 材料が弾性範囲内 — 座屈点で降伏していない

3. 荷重が比例的 — 全ての荷重が同じ比率で増減する

4. 初期不整が小さい — 理想形状からの偏差が無視できる程度

固有値座屈解析の結果は上界値（unconservative estimate）になる。つまり実際の崩壊荷重は固有値座屈荷重より低い。具体的にどの程度低いかは構造タイプによる：

だから円筒シェルの座屈設計には、固有値解析だけでは不十分で、初期不整を導入した非線形解析が不可欠なんだ。一方、柱の全体座屈は固有値解析でかなり良い推定が得られる。構造タイプに応じた「固有値解析の信頼度」を把握しておくことが実務では極めて重要だよ。

• 固有値問題 $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$ を解く — 2段階手順（静解析 → 固有値解析）
• 振動解析と双対的 — $[M]$ の代わりに $[K_\sigma]$
• $[K_\sigma]$ は応力状態に依存する負の剛性 — 圧縮が座屈を駆動する
• 結果は上界値 — 実崩壊荷重は構造タイプに応じて低下する
• 適用限界を理解した上で使う — 信頼度が低い構造には非線形解析を併用

その通り。設計の初期段階で「どこが危ないか」を素早く見つけるには最適な手法だ。ただし、それだけで設計OKを出せるかどうかは、構造の不整敏感性次第ということを忘れないでほしい。