板の座屈

根本的な違いが2つある。

1. 板は座屈後も荷重を担える。 柱は座屈すると荷重が急落するが、板は座屈後に荷重を再配分して追加荷重を担える。これが後座屈強度（post-buckling strength）だ。

2. 板の座屈は2次元の問題。 柱は1方向のたわみだけ考えればよいが、板は面内2方向の応力状態と面外たわみの相互作用を扱う。

その通り。航空機の翼外板は設計荷重の60〜70%で局所座屈するけど、スティフナー（桁、リブ）が荷重を再配分して全体としては持つ。後座屈強度を活用しない設計は過剰に重くなるから、航空機設計では板座屈の理論が必須なんだ。

薄板の座屈はvon Kármán方程式で記述される。一様圧縮を受ける矩形板の場合：

その通り。座屈応力は板厚の2乗に比例する。だから薄板の座屈設計では板厚が最も重要なパラメータだ。

$k$ は座屈係数（buckling coefficient）で、境界条件と荷重条件とアスペクト比 $a/b$ で決まる無次元パラメータだ。板の座屈問題の本質は、この $k$ を求めることに帰着する。

代表的な $k$ の値：

一辺が自由（拘束されていない）だと、その辺が自由に変形できるから座屈しやすい。H形鋼のフランジの突出部はまさにこれで、$k = 0.425$ という低い値になる。逆に両辺が固定されるウェブは $k = 6.97$ で座屈しにくい。

四辺単純支持の矩形板の座屈変位は：

そう。正方形板（$a/b = 1$）では $m = 1$ で $k = 4.0$。$a/b = 2$ なら $m = 2$ で $k = 4.0$。長い板ほど多くの半波で座屈する。

その通り。これは板座屈の重要な特徴で、幅 $b$ だけで座屈応力が決まる（長さ $a$ には依存しない）ということ。だから板の座屈設計では「幅/厚さ比」$b/t$ が最も重要なパラメータになる。

板が座屈すると中央部がたわんで剛性が低下するが、支持辺近くはまだ平面を保っていて荷重を担える。この「まだ有効に荷重を担える部分の幅」が有効幅 $b_e$ だ。

von Kármánの有効幅の式：

まさにそう。設計基準（ユーロコード3、AISI S100など）はこの有効幅の概念に基づいている。板の全幅ではなく有効幅で断面性能を計算し、その有効断面で応力チェックを行う。

要点：

いい表現だ。特に航空宇宙や薄板鋼構造では、座屈を「許容する」設計が標準だから、後座屈強度と有効幅の理解が不可欠なんだ。

板座屈は柱座屈に比べてFEM特有の問題が多い。最大の課題は要素タイプの選択だ。

基本的にはそうだ。薄板の座屈解析ではシェル要素が標準で、ソリッド要素は特殊な場合にのみ使う。

最低でも2次要素で2層は必要。1次要素だとシアロッキングで座屈荷重が著しく過大になる。ただしソリッドで薄板をモデル化するとアスペクト比が大きくなって要素品質が悪化するから、$b/t > 20$ 程度の薄板ではシェル要素一択だ。

有限ストリップ法（Finite Strip Method, FSM）が板座屈には非常に有効だ。板をストリップ（細長い短冊状の要素）に分割し、ストリップ方向にはフーリエ級数を仮定する。

FSMの利点：

FSMは一様断面の長尺部材に限定される。穴、補強、接合部、荷重の変化がある実構造には適用できない。だから実務ではFSMで初期スクリーニングし、FEMで詳細解析するのが最適な使い分けだ。

オープンソースのCUFSM（コーネル大学開発）でFSMを手軽に試せる。冷間成形鋼の設計基準（AISI S100）のDSM（Direct Strength Method）はCUFSMの結果を直接利用する前提で書かれている。

座屈の半波長に対して十分な要素数が必要だ。具体的には：

四辺単純支持の圧縮板なら、座屈半波長は板幅 $b$ とほぼ等しい。だからまず $b$ を6〜8等分するメッシュで始める。実際の問題ではまず粗いメッシュで座屈モードを確認し、波長に合わせてメッシュを細分化する。

重要なポイントだ。板の座屈解析では荷重の与え方が結果に大きく影響する。

違う。一様変位（端面が平面を保つ拘束）のほうが座屈荷重が高くなる。端面が自由に変形する一様分布荷重より、拘束されている分だけ座屈しにくい。差は5〜15%程度。実構造の端面条件がどちらに近いかを考えて選択する必要がある。

手順は非線形後座屈解析と同じだが、板特有のポイントがある：

そう。設計基準の有効幅公式はかなりの近似を含んでいるから、複雑な形状や荷重条件ではFEMの後座屈解析のほうが正確だ。ただし計算コストが高いので、標準的な問題では設計式を使い、非標準的な問題だけFEMで確認するのが実務的だ。

要点：

板座屈が設計の中心になる分野は多い：

分野ごとに異なる設計基準がある：

EN 1993-1-5 の手順はこうだ：

そう。$\bar{\lambda}_p \leq 0.673$ なら $\rho = 1.0$ で板全体が有効。これはClass 3以上の断面に相当する。$\bar{\lambda}_p > 0.673$ ならClass 4（薄肉断面）で有効幅による低減が必要だ。

スティフナー付きパネルの座屈は3つのレベルで考える必要がある：

ユーロコード3 Part 1-5 ではカラムモデル法を使う。スティフナー付きパネル全体を等価な柱として扱い、有効断面での柱座屈チェックを行う。内部では局所座屈による有効幅低減が組み込まれている。

FEMでは固有値座屈解析で3つのモードを全て抽出できる。モード形状を見れば、どのレベルの座屈が支配的かがわかる。設計式で扱いにくい非標準的なスティフナー配置では、FEMが威力を発揮する。

そう。I桁のウェブが支配的にせん断力を受ける場合、せん断座屈が生じる。ウェブに斜め45°方向のしわが入る変形パターンだ。

せん断座屈応力：

担える。これがテンションフィールド（張力場）理論だ。せん断座屈後、ウェブは引張の対角場を形成して荷重を伝達する。フランジとスティフナーがこの引張力のアンカーとなる。橋梁のプレートガーダーで後座屈のせん断強度を設計に活用するのは一般的だ。

その通り。どちらを採用するかは分野と構造によるが、軽量化が求められる場面では後座屈強度の活用が避けられない。FEMはどちらの設計フィロソフィーでも定量的な裏付けを与えてくれる。

板座屈は汎用FEMだけでなく、専用ツールもある。用途別に整理しよう。

コーネル大学のBen Schafer教授が開発したフリーソフト。MATLAB版とスタンドアローン版がある。

特徴：

その通り。米国のAISI S100だけでなく、オーストラリアのAS/NZS 4600でもCUFSMの使用が想定されている。

フランスCTICM開発のフリーソフト。ユーロコード3 Part 1-5に準拠した板座屈評価が可能。スティフナー付きパネルの座屈係数を直接計算できる。ユーロコード圏の設計者にとっては実質的な標準ツール。

航空宇宙分野ではESDU（Engineering Sciences Data Unit）のデータシートが座屈係数の権威的な参照先。板の圧縮座屈（ESDU 72019）、せん断座屈（ESDU 71005）、複合荷重（ESDU 81047）など、多数のシートがある。

チャートと式で座屈係数を求めるもので、FEMの代替ではない。ただしFEM結果の検証に非常に有用。航空宇宙の認証ではESDUとの比較が求められることが多い。

板座屈は基本的な問題だから、どの汎用ソルバーでも正確に解ける。差が出るのはワークフローの効率だ。

NastranのPCOMP/PCOMPG（複合材積層シェル定義）とSOL 200（最適化）の組み合わせが航空宇宙で圧倒的に使われている。各層の繊維角度・厚さを設計変数にして、座屈荷重を制約条件にした最適化が標準的なワークフローだ。この分野のノウハウの蓄積はNastranが頭一つ抜けている。

専用ツールは「速くて安い」が適用範囲が限定される。汎用FEMは「何でもできる」が手間とコストがかかる。両方を使いこなせるのが理想だ。

板座屈は古典的な問題だが、新材料と新しい設計手法で研究が活発化している。

本質的に異なる点が3つある：

1. 異方性。 繊維方向とそれに直交する方向で剛性が全く違う。$[0/90]$ 積層と $[\pm 45]$ 積層では座屈荷重が2倍以上変わる。

2. 座屈と破壊の関係。 金属板は座屈しても塑性変形で壊れるが、CFRP板は座屈時の局所曲げで層間剥離が発生し、急激に強度が低下する。後座屈強度の活用はより慎重になる。

3. 積層順序の自由度。 繊維角と積層順序の組み合わせは事実上無限にあり、最適設計の余地が大きい。

2層のアプローチがある：

まさにそこが研究のフロンティアだ。IJsselmuiden（デルフト工科大学）らのバイレベル最適化フレームワークや、遺伝的アルゴリズムによるスタッキングシーケンス最適化が活発に研究されている。

AFP（Automated Fiber Placement）で製造可能な可変剛性複合材だね。従来の一様な繊維角の積層に比べて、座屈荷重を30〜50%向上できるとの研究結果がある。

例えば圧縮パネルで、支持辺付近は $0°$ 繊維（荷重方向に強い）、中央は $\pm 45°$ 繊維（面内せん断で荷重を再配分）とすることで、応力を最適に分布させる。

要素ごとに異なる積層構成を定義する。汎用ソルバーでも可能だが、前処理が非常に煩雑になる。HyperMeshやDrapingシミュレーション（PAM-FORM等）との連携が実務的には必要だ。

座屈拘束ブレース（Buckling-Restrained Brace, BRB）は、鋼板の芯材を拘束材（コンクリート充填管など）で囲んで座屈を防止する制振部材だ。芯材が座屈せずに塑性圧縮できるため、引張と圧縮で対称的なヒステリシスループが得られる。

そう。BRBの設計では「拘束材のギャップ」と「芯材の座屈モード」の関係が重要で、FEMの固有値座屈解析で拘束の効果を検証する。ギャップが大きすぎると芯材が局所座屈してしまう。

コルゲート板は等方性板とは座屈挙動が全く異なる。波の方向とそれに直交する方向で曲げ剛性が大幅に違う直交異方性板として扱う。

面白い特性として、コルゲート板のせん断座屈荷重は波の振幅に強く依存する。波高が小さいコルゲートはほぼ平板に近い座屈をするが、波高が大きくなると「全体せん断座屈」ではなく「局所座屈」（波一つ一つが座屈する）に遷移する。

そう。コルゲートウェブ橋梁では、ウェブの面内剛性がほぼゼロのため、プレストレス力がウェブに伝わらずフランジに集中する。座屈の観点では有利だが、せん断座屈の評価が設計のポイントになる。

板座屈は「解けた問題」ではなく、新材料（CFRP、可変剛性）と新工法（AFP、3Dプリント）で常に新しい課題が生まれている。基礎理論（座屈係数、有効幅、Koiter理論）を押さえた上で、これらの応用に取り組んでほしい。

板座屈は柱座屈に比べて、FEM特有のトラブルが多い。代表的なものを見ていこう。

考えられる原因を順にチェック：

1. メッシュが粗い — 板座屈は最低6×6要素（2次シェル）が必要。粗いメッシュは剛性を過大評価し、座屈荷重が高くなる方向にずれる。

2. 境界条件が厳しすぎる — 「単純支持」のつもりが回転を拘束してしまっている。シェル要素の場合、面外変位（$w = 0$）だけ拘束すべきところ、面外回転（$\theta = 0$）まで拘束すると固定支持になる。

3. 面内拘束の影響 — 面内変位の拘束条件が理論の仮定と合っていない。面内自由な辺で面内を拘束すると剛性が上がる。

圧倒的に多い。板座屈の「単純支持」は：

荷重辺の拘束条件を確認。荷重辺が面内に完全に拘束（一様変位）されていると、短辺方向の変形も拘束されて高次モードが抑制される。

もう一つの可能性はアスペクト比に対してモード数が求められている固有値の中に入っていないこと。$a/b = 3$ では $m = 3$ のモードの $k$ 値は $m = 1$ よりわずかに高いが、$m = 2$ のモードはさらに高い。固有値を多めに（20モード程度）求めて、各モードの半波数を確認しよう。

応力集中部やコーナー部で局所座屈モードが出現している可能性がある。穴の周辺、荷重集中点、形状の不連続部では、メッシュを細かくするほど局所的なモードが低い固有値で出る。

対策：

設計基準では相互作用式を使う：

そう。FEMの固有値座屈解析で圧縮とせん断を同時に参照荷重として与えれば、相互作用を自動的に考慮した座屈荷重が得られる。ただし荷重比が固定される点に注意。荷重比が変わる場合は複数ケースの解析が必要だ。

これは実務で本当に起きるミスだ。H形鋼のフランジは片側がウェブに接続（支持辺）、もう片側が自由端。$k = 0.425$ で評価すべきところを $k = 4.0$ にしたら、座屈応力が10倍になってしまう。断面のどの部分が「内部板」でどの部分が「突出板」かを正しく認識することが板座屈設計の第一歩だ。