板の座屈 — 理論と支配方程式

カテゴリ: 構造解析 | 2026-01-15
plate-buckling-theory
理論と物理の世界へ

柱の座屈との違い

🧑‍🎓

オイラー座屈は1次元の柱でしたが、板(2次元)の座屈は何が違うんですか?


🎓

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


🎓

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


🎓

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


🧑‍🎓

座屈しても荷重を担えるって、すごく実用的ですね。航空機の外板がそうだと聞きました。


🎓

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


板の座屈の支配方程式

🧑‍🎓

板の座屈を支配する方程式は何ですか?


🎓

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


$$ D\nabla^4 w + N_x \frac{\partial^2 w}{\partial x^2} + 2N_{xy} \frac{\partial^2 w}{\partial x \partial y} + N_y \frac{\partial^2 w}{\partial y^2} = 0 $$

ここで $D = Et^3 / 12(1-\nu^2)$ は板の曲げ剛性、$w$ は面外たわみ、$N_x, N_y, N_{xy}$ は面内応力合力だ。


🧑‍🎓

$D$ の中に板厚 $t$ の3乗が入っている! 板厚が倍になると曲げ剛性は8倍…板の座屈は板厚に極めて敏感なんですね。


🎓

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


座屈係数 $k$

🧑‍🎓

教科書で $\sigma_{cr} = k \cdot \pi^2 D / (b^2 t)$ という式を見ましたが、$k$ は何ですか?


🎓

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


🎓

代表的な $k$ の値:


荷重境界条件$k$備考
一様圧縮四辺単純支持4.0基準ケース($a/b \geq 1$)
一様圧縮荷重辺単純支持+非荷重辺固定6.97固定辺の拘束効果
一様圧縮荷重辺単純支持+一辺自由0.425フランジの突出板座屈
純せん断四辺単純支持5.34 + 4.0$(b/a)^2$せん断座屈
純曲げ四辺単純支持23.9曲げ圧縮の座屈
🧑‍🎓

$k = 4.0$ と $k = 0.425$ で10倍近い差! 境界条件がこんなに効くんですか。


🎓

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


座屈モード形状

🧑‍🎓

板の座屈モードはどんな形ですか?


🎓

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


$$ w(x,y) = A \sin\frac{m\pi x}{a} \sin\frac{\pi y}{b} $$

ここで $m$ は荷重方向の半波数だ。$m$ が変わると座屈係数 $k$ も変わる:


$$ k = \left(\frac{m}{a/b} + \frac{a/b}{m}\right)^2 $$

🧑‍🎓

$k$ が最小になる $m$ が実際の座屈モードですよね。


🎓

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


🧑‍🎓

どんなに長くしても $k = 4.0$ に収束するんですか!


🎓

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


有効幅の概念

🧑‍🎓

板が座屈した後、どうやって荷重を担うんですか?


🎓

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


🎓

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


$$ b_e = b \sqrt{\frac{\sigma_{cr}}{\sigma}} = t \sqrt{\frac{k \pi^2 E}{12(1-\nu^2)\sigma}} $$

🧑‍🎓

有効幅は作用応力が大きくなると狭くなる…つまり荷重が増えるほど板の有効断面が減っていくんですね。


🎓

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


まとめ

🧑‍🎓

板座屈の理論、整理します。


🎓

要点:


  • 板は座屈後も荷重を担える — 柱と本質的に異なる
  • $\sigma_{cr} = k \pi^2 D / (b^2 t)$ — 座屈係数 $k$ が境界条件と荷重で決まる
  • 板厚 $t$ が支配的 — $D \propto t^3$ なので板厚感度が非常に高い
  • 長い板の $k$ は幅 $b$ だけで決まる — $b/t$ が設計の鍵
  • 有効幅 — 座屈後の荷重伝達能力を表す設計概念

🧑‍🎓

板の座屈は「壊れる」のではなく「荷重の流れが変わる」イメージですね。


🎓

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


Coffee Break よもやま話

タコマナローズ橋の崩壊(1940年)

完成からわずか4ヶ月で崩壊した吊り橋。風速わずか65km/hで起きた空力弾性フラッター(共振)が原因でした。この事故は「振動解析を怠るとどうなるか」の最も有名な教訓として、今でも構造力学の教科書に載っています。現代のCAEは、この種の問題を設計段階で発見できます。もし当時にCAEがあれば、橋は今も架かっていたかもしれません。

各項の物理的意味
  • 慣性項(質量項):$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか? あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
  • 剛性項(弾性復元力):$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね? あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか? 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い:「剛性が高い=強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
  • 外力項(荷重項):体積力 $f_b$(重力など)と表面力 $f_s$(圧力、接触力など)。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」(体積力)、タイヤが路面を押す力は「表面だけにかかる力」(表面力)。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗:荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
  • 減衰項:レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか? いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら? 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
  • 連続体仮定:材料を連続的な媒質として扱い、ミクロな不均質性を無視する
  • 微小変形仮定(線形解析の場合):変形が初期寸法に比べて十分小さく、応力-歪み関係が線形
  • 等方性材料(特に指定がない場合):材料特性が方向に依存しない(異方性材料では別途テンソル定義が必要)
  • 準静的仮定(静解析の場合):慣性力・減衰力を無視し、外力と内力の釣り合いのみを考慮
  • 適用外ケース:大変形・大回転問題では幾何学的非線形性が必要。塑性・クリープ等の非線形材料挙動では構成則の拡張が必要
次元解析と単位系
変数SI単位注意点・換算メモ
変位 $u$m(メートル)mm入力時は荷重・弾性率もMPa/N系に統一すること
応力 $\sigma$Pa(パスカル)= N/m²MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意
歪み $\varepsilon$無次元(m/m)工学歪みと対数歪みの区別に注意(大変形時)
弾性率 $E$Pa鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意
密度 $\rho$kg/m³mm系ではtonne/mm³(= 10⁻⁹ tonne/mm³ for 鋼)
力 $F$N(ニュートン)mm系ではN、m系ではNで統一

数値例:片持ち梁の先端荷重(L=1m, 断面50×100mm, 鋼材E=210GPa, P=1kN)

最大たわみ δ = PL³/(3EI) = 1000×1000³/(3×210000×4,166,667) ≈ 0.381 mm 最大応力 σ = PL×(h/2)/I ≈ 12.0 MPa(降伏応力235MPaに対して安全率19.6)

メッシュ密度を変えた収束性の確認:

粗いメッシュ(500要素)0.362 mm
-5.0%
中程度(2,000要素)0.378 mm
-0.8%
細かいメッシュ(8,000要素)0.380 mm
-0.3%
理論解0.381 mm
基準

ポイント:要素数を4倍にしても結果は0.5%しか変わらない→8,000要素で十分収束。これが「メッシュ収束性」の確認です。

簡易計算ツール:構造力学基礎

片持ち梁の先端集中荷重における最大たわみ・最大応力を計算します。

単軸応力状態における応力・歪み・伸びの相互換算。

CAE実務でよく使う単位の換算。

構造解析の収束問題や計算コストに課題を感じていませんか? — Project NovaSolverは、実務者が日々直面するこうした課題の解決を目指す研究開発プロジェクトです。

CAEの未来を、実務者と共に考える

Project NovaSolverは、板の座屈における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。

プロジェクトの最新情報を見る →