非線形後座屈解析 — 理論と支配方程式

3つの理由がある。

1. 実構造は座屈しても即座に崩壊しない場合がある。 例えば航空機の翼の外板は運用荷重下で局所座屈しているが、全体構造としては十分な強度を持つ。これを「後座屈強度の活用」という。

2. 固有値座屈は上界値であり、実崩壊荷重を過大評価する。 初期不整や材料非線形を含む実際の崩壊荷重は、非線形解析でないと正確に評価できない。

3. 崩壊モードの特定。 構造が「どのように壊れるか」を知ることは、安全設計の基本だ。荷重-変位の全経路を追跡することで、エネルギー吸収能力や延性も評価できる。

そう。薄いアルミパネルのスティフナー間は設計荷重以下で座屈するけど、スティフナーが荷重を再配分して全体としては持つ。この後座屈耐力を正しく評価できれば、構造をさらに軽量化できる。逆に後座屈を無視して全部座屈しないように設計すると、過剰に重くなる。

非線形構造力学では、荷重パラメータ $\lambda$ と変位ベクトル $\{u\}$ の関係を平衡経路（equilibrium path）と呼ぶ。この経路上に2種類の特異点がある：

分岐点（bifurcation point） — 平衡経路が分かれる点。完璧な構造でのオイラー座屈はこれ。

限界点（limit point） — 荷重が極大値に達する点。スナップスルー座屈はこれ。

そう。ただし物理的な意味は全く違う。分岐点では構造が「別の変形モードに移る」。限界点では構造が「それ以上荷重を担えなくなる」。後座屈解析では両方を正しく追跡する必要がある。

Warner T. Koiterだ。1945年のデルフト工科大学の博士論文（オランダ語で書かれていたため長く知られなかった）で、後座屈経路の安定性を系統的に分類した。

分岐点近傍の後座屈経路をべき級数展開すると：

ここで $\xi$ はモード振幅パラメータ。

• $a \neq 0$（非対称分岐） — 後座屈経路が一方向に傾く。不整に極めて敏感
• $a = 0, b > 0$（安定対称分岐） — 後座屈で荷重が上がる。平板の座屈がこれ
• $a = 0, b < 0$（不安定対称分岐） — 後座屈で荷重が落ちる。円筒シェルの座屈がこれ

まさにそう。Koiterは「後座屈経路の形状」から「不整に対する敏感性」を予測できることを示した。これは構造力学における最も重要な理論的貢献の一つだ。

大変形理論では接線剛性マトリクスは3つの項で構成される：

• $[K_0]$ — 小変位の弾性剛性（材料剛性）
• $[K_\sigma]$ — 幾何剛性（初期応力の影響）
• $[K_L]$ — 大変位剛性（変位による形状変化の影響）

その通り。線形座屈は「座屈前の変形は微小」という仮定だから $[K_L] \approx 0$。しかし後座屈では変位が大きくなるので $[K_L]$ が無視できなくなる。$[K_L]$ を含めることが「幾何学的非線形解析」の核心だ。

そう。弾塑性の後座屈解析では3つの項全てが変位と応力に依存して変化する。だからこそ増分反復法（Newton-Raphson法の荷重増分版）で少しずつ経路を追跡する必要があるんだ。

要点：

• 後座屈は「崩壊荷重の正確な評価」と「崩壊モードの特定」のために必要
• 平衡経路上の分岐点と限界点 — $[K_T]$ が特異になる2種類の点
• Koiterの理論 — 後座屈経路の形状（$a, b$ の符号）が不整敏感性を決定
• $[K_T] = [K_0] + [K_\sigma] + [K_L]$ — 大変位項 $[K_L]$ が非線形解析の核心
• 座屈後も荷重を担える構造がある — 後座屈強度の活用は軽量化の鍵

その理解は完璧だ。次回は、この経路追跡を実際にどうやるか — Riks法（弧長法）の数値アルゴリズムを詳しく見ていこう。