座屈解析（Buckling Analysis） — CAE用語解説

座屈解析は、構造物が圧縮荷重を受けたときに突然横方向にバコッと変形する——いわゆる座屈現象を予測する解析だ。普通の線形静解析は荷重に比例して変形が増えることを前提にしているけど、座屈はある荷重に達すると急に別の変形モードに飛び移る不安定現象だから、静解析だけでは捉えられないんだよ。

一番わかりやすいのは定規だね。長い定規の両端を手で押すと、最初は縮むだけだけど、ある力を超えた瞬間にグニャッと横に曲がるだろう？ あれが座屈だ。実務だと、薄肉の航空機パネル、建築の柱、圧力容器の胴体、自動車のフレームなど、圧縮を受ける薄い構造では常に座屈が設計上の主要リスクになる。

そうだ。理想的な真っ直ぐで一様な断面の柱に対する、古典的な座屈荷重の公式がこれだ：

その通り。$K$ が小さいほど有効長さ $KL$ が短くなって、座屈荷重 $P_{cr}$ は大きくなる。つまり両端固定は両端ピンの4倍の座屈荷重に耐えられる計算になる。ただし実際のボルト接合やフランジ接合が「完全固定」なのかは別問題で、現実には固定とピンの中間的な拘束になることが多い。そこが設計の難しさだね。

できるよ。座屈荷重を断面積 $A$ で割れば座屈応力になる：

大きく分けるとこうだ。線形座屈解析（固有値座屈解析とも呼ぶ）は、剛性行列 $[K]$ と幾何剛性行列（応力剛性行列）$[K_\sigma]$ を使って、次の固有値問題を解く：

線形座屈にはいくつか致命的な限界があるんだ。まず、初期不整を考慮できない。完璧にまっすぐな柱を仮定しているから、実際より高い座屈荷重を出してしまう。次に、材料非線形（塑性）を含められない。降伏が始まったあとの剛性低下を無視している。そして大変形の影響を追跡できない。だから線形座屈の結果はあくまで「上界推定」と考えるべきだ。

非線形座屈解析は、荷重を少しずつ段階的に増やしながら、毎ステップで幾何学的非線形（大変形）や材料非線形を含む平衡方程式を反復法（Newton-Raphson法など）で解いていく。荷重-変位曲線の勾配がゼロになる点、あるいは荷重が極大値を示す点が座屈荷重に相当する。座屈後の挙動（後座屈）を追跡するには、弧長法（Riks法）を使って荷重の減少も追えるようにする必要がある。

そうだ。通常の荷重制御だと、座屈点で荷重を増やせなくなって計算が発散する。弧長法は荷重と変位の両方をパラメータにして、荷重-変位空間上の「弧」に沿って追跡するから、スナップスルーやスナップバックのような不安定な平衡経路も追えるんだ。AbaqusではStatic, Riks、NastranではSOL 106のARC-LENGTH法として実装されている。

これは本当に重要なポイントだ。初期不整敏感性（imperfection sensitivity）といって、構造の種類によって影響度が全く違う。例えば柱の座屈は比較的鈍感で、少々曲がっていてもEuler値に近い荷重で座屈する。しかし薄肉円筒シェルの軸圧縮座屈は極端に敏感で、理論値の20〜40%の荷重で座屈してしまうことがある。有名な話で、NASAの実験データが理論値から大幅にばらついたのがきっかけで、ノックダウンファクタという経験的な補正係数が導入されたんだ。

実務的なアプローチはこうだ。まず線形座屈解析で座屈モード形状を求める。次に、その第1座屈モード形状を微小な振幅（板厚の10〜50%程度）でFEモデルの節点座標に重ね合わせて初期不整を導入する。その初期不整付きモデルで非線形座屈解析を実行すれば、より現実的な座屈荷重が得られる。複数のモード形状を組み合わせたり、不整の振幅をパラメトリックに変えて敏感性を調査することも多い。

業界によって異なるけど、一般的には製造公差に対応する値を使う。例えば航空宇宙ではNASA SP-8007が薄肉シェルのノックダウンファクタを規定しているし、鋼構造ではEurocodeがEN 1993-1-6で初期不整の振幅を板厚の割合として規定している。CAE解析では「板厚の10%」「板厚の50%」「板厚の100%」など複数ケースで感度解析するのが安全だよ。

いくつかあるよ。まずメッシュの細かさ。座屈モードは局所的な波形を含むことが多いから、波長の1/4以下のメッシュサイズが必要だ。特にシェル構造では、シェル要素の選択（S4R vs S4など）も結果に影響する。次に境界条件の忠実な再現。先ほど話した $K$ 値のように、実際の拘束条件が結果を大きく左右する。

負の固有値は、基準荷重の方向を逆にした（つまり引張にした）ときに座屈するモードを意味する。圧縮で座屈するモードは正の固有値に対応する。ソルバーは正負両方の固有値を出すから、圧縮の座屈を調べたいなら正の値に注目すればいい。ただし、荷重の方向設定を間違えていないかの確認にもなるから、負の値が出たこと自体は無視しなくていいよ。

座屈点の近傍は解が分岐するから、Newton-Raphson法が収束しなくなるのは当然ありえる。対策としては、(1) 弧長法（Riks法）を使う、(2) 初期不整を導入して完全な分岐を避ける（対称性が崩れることで荷重-変位曲線が滑らかになる）、(3) 増分ステップを細かくする、(4) 人工的な減衰（stabilize）を加える——ただし減衰が大きすぎると結果が歪むから注意が必要だ。実務では(1)と(2)の組み合わせが最も信頼性が高い。

完璧にそのとおり。線形座屈で「どのくらいの荷重で、どんな形で座屈しそうか」を把握して、その情報を使って非線形解析モデルを構築する。これが座屈設計の王道パターンだ。特に安全性が重要な航空宇宙・原子力・圧力容器の分野では、このプロセスが設計コードで明確に要求されている。