オイラー座屈 — 先端技術と研究動向

オイラー座屈は「理想化された柱の弾性座屈」という最もシンプルなケースだ。実構造の座屈問題はそこからどんどん複雑になる。大きく3つの方向に広がっている。

これは構造の種類によって劇的に異なる。Koiterの理論がこの問題の基礎を築いた。

Koiterは1945年の博士論文で、初期不整に対する座屈荷重の感度を系統的に分類した。後座屈経路の性質によって：

• 安定な後座屈（例：平板の座屈） — 座屈後も荷重を担持でき、不整敏感性は低い
• 不安定な後座屈（例：円筒シェルの軸圧縮） — 座屈後に荷重が急落し、不整敏感性が極めて高い

古典理論の座屈荷重 $\sigma_{cr} = 0.605 Et/R$（$t$: 板厚、$R$: 半径）に対して、実験値は理論値の20〜40%しかない。この乖離が1930年代から大問題になっていた。NASAが1968年に発表したSP-8007という設計ガイドラインでは、ノックダウンファクター（理論値に対する低減係数）を定めていて、大径の薄肉シェルでは0.2程度まで下がる。

現在はこのSP-8007が「過度に保守的」とされ、NASA・ESA・DLRなどが新しいノックダウンファクター手法を研究している。SBPA（Single Boundary Perturbation Approach）や確率論的不整敏感性解析が主流で、実測の不整データを使って個体ごとの座屈荷重を予測する方向に進んでいるんだ。

中間柱（intermediate column）の領域は、材料の非線形性が座屈に絡んでくる。古典的なアプローチとしては：

• 接線係数理論（Engesser-Shanley） — 座屈時の弾性係数を接線係数 $E_t$ に置き換え：

その通り。ただしShanleyが1947年に示した重要な洞察がある。座屈の瞬間に除荷が起きる繊維と載荷が続く繊維が混在するため、実際の座屈荷重は接線係数荷重と低減係数荷重の間になる。実用的にはShanleyの結果から、接線係数荷重が下限値として使えることがわかっている。

非常に重要だ。弾塑性座屈をFEMで解くには：

1. 適切な応力-ひずみ曲線を入力（真応力-真ひずみで）

2. 幾何学的非線形（NLGEOM）を有効に

3. 弧長法で平衡経路を追跡

4. 初期不整を導入（弾塑性座屈でも不整敏感性はある）

単なる弾性解析に降伏応力でカットオフする設計式（Johnsonの放物線など）は便利だけど、それはFEMの代替にはならない。

これがモード相互作用（mode interaction / mode coupling）で、座屈の中でも最も厄介な問題の一つだ。

典型例はH形鋼柱の全体座屈と局所座屈の連成だ。全体座屈（柱全体が横に倒れる）と局所座屈（フランジが波打つ）の臨界荷重が近いとき、両方のモードが相互作用して、どちらか単独の座屈荷重よりさらに低い荷重で崩壊する。

線形座屈では各モードの固有値が近いことは確認できるけど、相互作用の影響は非線形解析でないとわからない。対処法としては：

• 1次モードと2次モード（あるいはそれ以上）を重ね合わせた初期不整を与えて非線形解析
• モード振幅の比率を変えて複数ケース実施
• Koiterの理論に基づくマルチモード不整敏感性解析

だからこそ、先に線形座屈で「危ないモードの組み合わせ」を見つけてから、ピンポイントで非線形解析を回すのが効率的だ。全パラメータを非線形で回す余裕は普通ない。

動的座屈（dynamic buckling）は、荷重が急激に作用する場合の座屈だ。爆風荷重を受ける建物の柱、衝突時の車体構造、ロケット発射時の推力変動など。

静的座屈との違いは2つある：

1. 慣性効果 — 質量が座屈の発達を遅らせる。極めて短時間の荷重なら、静的座屈荷重を超えても座屈が発達しきらない

2. 座屈判定基準の曖昧さ — 静的座屈は固有値で明確に定義できるが、動的座屈は「変位がどこまで増大したら座屈と見なすか」の判定が一義的でない

陽解法（explicit法）を使うのが一般的だ。LS-DYNAやAbaqus/Explicitが代表的。陰解法の固有値座屈とは全く異なるアプローチで、時間積分で実際の変形過程を追いかける。Budiansky-Roth基準（荷重パルスの持続時間と応答振幅の関係）で臨界条件を評価する手法がよく使われる。

座屈荷重の最大化を目的関数にした構造最適化は活発な研究分野だ。

• 寸法最適化 — 板厚やフランジ幅を変数にして座屈荷重を最大化。比較的シンプル
• トポロジー最適化 — 材料配置自体を最適化。座屈制約を入れると計算が不安定になりやすく、固有値の感度計算が技術的に難しい
• 複合材のスタッキングシーケンス最適化 — 各層の繊維角を変えて座屈荷重を最大化。航空機パネルで特に重要

同じ厚さでもスタッキングシーケンス（積層順序）で座屈荷重が2倍以上変わることがある。特に $[\pm 45]$ 系は面内せん断に強く、$[0/90]$ 系は面内圧縮に強い。この最適化にはNastranのPCOMPとSOL 200の組み合わせが定番だ。

いくつかの方向が見えている：

• 確率論的座屈評価 — 決定論的な安全率ではなく、破壊確率で設計。不整の統計データが鍵
• 機械学習によるサロゲートモデル — 数千回の非線形解析をニューラルネットで代替し、不整敏感性を高速評価
• 等幾何解析（IGA） — CADのNURBS形状をそのままFEMの基底に使い、形状近似誤差をゼロに。薄肉シェルの座屈解析で特に有効
• 3Dプリンティング構造の座屈 — ラティス構造の座屈は従来の板・シェルの理論では扱えず、新しい設計指針が必要

まさに。規格やコードが追いついていない分野だから、CAEエンジニアの判断力がより一層問われる。基礎理論（オイラー座屈）をしっかり理解しておくことが、最先端の問題に取り組むときにも必ず役に立つよ。