線形座屈（固有値座屈）解析

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

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

2段階の手順になる。

そう。ここで重要なのは、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]$ だ。

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

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

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

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

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

構造解析の固有値問題は典型的には一般化固有値問題だ。座屈の場合、標準形に変換すると：

Lanczos法は大規模疎行列の固有値問題に特化したアルゴリズムで、全固有値を求めずに指定した範囲の固有値だけを効率的に抽出できる。座屈解析では下位数モード（$\lambda$ が小さい方）だけが必要だから、最適なんだ。

アルゴリズムの骨格はこうだ：

いい質問だ。$[K_0]^{-1}$ は明示的には作らない。代わりに $[K_0]\{x\} = [K_\sigma]\{q_j\}$ という連立方程式を解く。これなら $[K_0]$ の疎行列構造を維持したまま計算できる。事前に $[K_0]$ をLU分解（またはLDLT分解）しておけば、各反復での連立方程式の求解は前進後退代入だけで済む。

まさにその通り。だから固有値問題のコストの大部分は行列分解が占める。これは通常の静解析とほぼ同じコストだ。つまり固有値座屈解析は、静解析の1.5〜2倍程度の計算時間で済むことが多い。

そういう場合はシフト-反転Lanczos法が有効だ。注目する固有値 $\sigma$ の近くを重点的に探索する：

実務的には、まずシフトなしで粗い固有値を求め、特定のモードを詳細に調べたいときにシフトを設定する。Nastranでは EIGRL カードの V1, V2 パラメータで固有値の探索範囲を指定できるし、Abaqusでも同様の機能がある。

サブスペース反復法（Bathe法）もまだ使われている。特に固有値が密集している場合にLanczos法よりロバストなことがある。

違いを整理すると：

99%のケースではLanczosで問題ない。サブスペース法に切り替えるのは、Lanczosが収束しない場合（まれにある）の代替手段と考えておけばいい。

シェル要素の場合で説明しよう。膜応力状態 $\{N_{xx}, N_{yy}, N_{xy}\}$ を持つシェル要素の幾何剛性は：

そう。逆に圧縮やせん断が支配的だと $[K_\sigma]$ が負に効いて剛性が下がる。面白いのは面内せん断 $N_{xy}$ だけで座屈を引き起こすケースもあること。せん断座屈は航空機の翼のウェブやI桁のウェブで典型的だ。

いくつかのテクニックがある：

例えば対称面でz方向変位をゼロにすると、対称モードしか出ない。反対称モードは別途、反対称境界条件（z方向自由、x方向回転拘束など）で解く必要がある。対称条件を使うなら、対称・反対称の両方を解くのが鉄則だ。

要点を整理すると：

設計フェーズによって使い方が変わる。整理しよう。

そう。固有値座屈解析は単独で設計判断に使う場合と、非線形解析の前処理として使う場合がある。後者では座屈モード形状を初期不整として利用するのが典型的だ。

ここは初心者が混乱しやすいポイントだ。座屈荷重係数 $\lambda$ は参照荷重に対する倍率だから、参照荷重が変われば $\lambda$ も変わる。でも臨界荷重 $P_{cr} = \lambda \times P_{ref}$ は同じだ。

そう。実務的には設計荷重（使用荷重）を参照荷重にすると便利だ。$\lambda > 1.0$ なら設計荷重の範囲で座屈しない、$\lambda < 1.0$ なら座屈する。設計コードの安全率と直接比較できる。

そういうことだ。ただし、これは固有値座屈が上界値であるという前提のもとでの話。不整敏感性の高い構造では、さらにノックダウンファクターを掛ける必要がある。

これは固有値座屈解析の重要な制約だ。全ての荷重が同じ比率でスケーリングされる前提だから、独立にスケーリングされる複数荷重系には直接適用できない。

対策は3つある：

その場合はプリロード＋座屈荷重の手法を使う。まず自重だけで静解析（プリロード）し、その応力状態を初期条件として、風荷重の座屈解析を行う。Nastranでは SOL 105 の STATSUB + BUCKLE サブケースで実現できる。Abaqusでは Staticステップの後に Buckleステップを置く。

モード形状からは座屈の物理的メカニズムが見える。チェックポイントは：

1. 座屈の種類の判別

2. 変形が集中する部位

3. モード間の差異

極めて危険だ。例えば $\lambda = 5.0$ で安全に見えても、モード形状が局所的な板の波打ちだったら、メッシュが粗くて本当の局所座屈を捕捉していない可能性がある。メッシュを細かくしたら $\lambda = 1.2$ に急落した、という事例は珍しくない。

設計コードによってアプローチが異なる：

そう。ユーロコード3の6.3.4項「General method」では、FEMによる弾性座屈解析の固有値を直接使って座屈耐力を評価できる。最もモダンなアプローチだ。

そう。固有値座屈解析は答えの数字よりも、モード形状から読み取れる物理に価値がある。$\lambda$ が何倍かを報告するだけなら、解析をやる意味は半分もない。

基本的な考え方（静解析→固有値解析の2段階）は同じだけど、設定の書き方やオプションは結構違う。ソルバーごとに見ていこう。

Nastranの線形座屈はSOL 105で、Executive Control と Case Control の組み合わせで設定する。

これがプリロード付き座屈解析のキモだ。自重や恒荷重をサブケース1で与え、変動荷重（風荷重など）の座屈をサブケース2で評価する。METHOD = 10 は EIGRL カードを参照していて、そこで固有値ソルバーの設定（モード数、手法）を指定する。

Nastranの座屈解析で注意すべきテクニカルポイント：

Abaqusでは2つのステップを順に設定する：

eigensolver, number_of_vectors, max_iterations, blocksize の順だ。デフォルトの Lanczos で通常は問題ない。サブスペース法を使いたければ明示的に指定する。

Abaqusのポイント：

Workbenchの手順は直感的だ：

APDLの場合：

注意点: APDLでは PSTRES, ON を忘れると応力剛性マトリクスが計算されず、座屈解析が無意味な結果を出す。Workbenchでは自動的にONになるが、APDLでは手動で指定が必要。これは非常によくあるミスだ。

理想的には同じだが、実際には微妙に異なる。要因は：

そう。ベンチマーク問題（NAFEMSの標準問題など）では5%以内の一致が得られれば十分。10%以上ずれたら、モデルの解釈に差がある可能性を疑うべきだ。クロスチェックの際は、まずモード形状が同じかどうかを視覚的に確認し、その上で $\lambda$ の値を比較すること。

OptiStructの最大の特徴は座屈制約付きトポロジー最適化に対応していること。通常のトポロジー最適化は剛性最大化が目的だが、OptiStructでは「座屈荷重係数 $\lambda$ が指定値以上」という制約を追加できる。

その通り。薄肉構造の軽量化では、材料を削りすぎると座屈してしまうジレンマがある。座屈制約付き最適化はこのトレードオフを自動で解いてくれる。航空宇宙や自動車のパネル設計で特に有用だ。ただし、座屈固有値の感度計算が数値的に不安定になりやすいので、計算パラメータの調整が必要なこともある。

判断のポイント：

その通りだ。固有値座屈解析自体はどのソルバーでも成熟した技術。差別化要因はワークフロー全体の効率にある。

手法自体は1960年代に確立しているけど、使い方と解釈の方法にはまだ大きな研究の余地がある。特に以下の3つの方向が活発だ。

従来のFEMの固有値座屈は、離散化された行列の固有値問題として解く。一方、連続体ベース座屈は微分方程式レベルで座屈条件を記述し、連続的な座屈場を求める。

具体的には、シェルの座屈では Donnell 方程式や Sanders-Koiter 方程式を数値的に解く。FEMの離散化誤差なしに「連続体としての座屈荷重」が得られる。メッシュ依存性がない点が大きな利点だ。

そう。現状では回転シェル（円筒、円錐、球）のような単純形状が対象。だけど、これらはまさに不整敏感性が高い構造であり、FEMのメッシュ依存性が問題になりやすい構造でもある。補完的に使うと効果的だ。

固有値 $\lambda$ の設計変数 $s$ に対する感度は、随伴法（adjoint method）で効率的に計算できる：

その通り。これが座屈最適化の基盤になる。板厚を0.1 mm変えたとき $\lambda$ がどれだけ変わるか…この感度情報を使って勾配法で最適化を回す。

ただし落とし穴がある。固有値が近接（クラスタリング）しているとき、感度が不連続になる。モードの入れ替わり（mode switching）が起きて、最適化が振動的になる。

対策としては：

歪み座屈（distortional buckling）だね。薄板リップ溝形鋼やZ形鋼で、断面形状自体が変わる座屈モードだ。

一般化梁理論（GBT）は、薄肉断面の座屈を「変形モード」に分解して解析する手法だ：

そう。GBTの大きな利点は座屈モードの参加度（modal participation）がわかること。例えば「この座屈モードは60%が歪み座屈、30%が全体曲げ座屈、10%がねじり座屈」と定量化できる。これは設計の判断に非常に有用だ。

オープンソースのGBTUL（リスボン大学開発）やCUFSM（コーネル大学開発、有限ストリップ法ベース）で手軽に試せる。設計コード（ユーロコード3、AISI S100）のDSM（Direct Strength Method）にも直結する。

従来の決定論的な座屈評価（ノックダウンファクター）に代わる手法として、確率論的座屈評価が注目されている。

手順はこうだ：

確かにコストは大きい。だからサロゲートモデル（Kriging、ニューラルネットワーク等）で非線形解析を代替する研究が盛んだ。数十〜百回の非線形解析結果からサロゲートを構築し、数万回のモンテカルロサンプリングをサロゲート上で実行する。

サロゲートモデルの入力変数を設計する段階で使う。固有値座屈のモード形状を初期不整の基底として利用し、各モードの振幅をランダム変数とする。つまり固有値座屈は「不整のパラメトリゼーション」の道具として活躍するんだ。

3つの方向が見える：

固有値座屈は「古い手法」と侮れない。座屈問題のスクリーニング、モード分類、最適化の勾配情報、確率論的解析の基底…どの先端手法にも固有値座屈が組み込まれている。基礎がしっかりしているからこそ、応用が広がるんだ。

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

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

原因チェック：

そう。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ケースが最低限。最も低い崩壊荷重を設計値とする。

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

対策：

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

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

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

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

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

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

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