熱座屈解析 — 理論と支配方程式

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

熱座屈とは

🧑‍🎓

先生、熱で座屈が起きるんですか? 力がかかっていないのに?


🎓

いい質問だ。熱座屈(thermal buckling)は温度変化によって生じる圧縮応力が原因で起きる座屈だ。構造が自由に膨張できれば温度応力は生じないが、膨張が拘束されていると圧縮応力が発生する。


🧑‍🎓

線路のレールが暑い日に曲がるのも熱座屈ですか?


🎓

まさにそう。レールは枕木で軸方向の膨張が拘束されている。温度上昇 $\Delta T$ で生じる軸応力は:


$$ \sigma_{th} = E \alpha \Delta T $$

ここで $\alpha$ は線膨張係数。鋼($\alpha \approx 12 \times 10^{-6}$ /°C, $E = 200$ GPa)で $\Delta T = 40$ °Cなら $\sigma_{th} = 96$ MPa。これが座屈応力を超えるとレールが横に曲がる。


🧑‍🎓

40°Cの温度差で96 MPa…思ったより大きいですね。


🎓

そう。熱応力は拘束度に比例する。完全拘束なら $E\alpha\Delta T$ だが、部分的に膨張が許されれば応力は下がる。実構造の拘束度を正しく評価することが熱座屈解析の第一歩だ。


熱座屈の支配方程式

🧑‍🎓

熱座屈の定式化は機械的な座屈と違いますか?


🎓

固有値座屈の枠組みは同じだ。ただし参照荷重が温度荷重になる:


1. 温度分布 $\Delta T(x,y,z)$ を与えて熱応力解析を実施

2. 得られた熱応力で $[K_\sigma]$ を構成

3. 固有値問題 $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$ を解く


$\lambda$ が座屈温度荷重係数。$\lambda \cdot \Delta T$ が臨界温度上昇。


🧑‍🎓

温度分布が一様でない場合もありますよね。


🎓

実構造では温度分布は通常不均一だ。例えば火災時の梁は下面が高温、上面が比較的低温。板の片面だけ加熱されると温度勾配による曲げモーメントも発生する。


🎓

温度の効果は2つある:

  • 膜応力(面内の圧縮/引張) — 座屈の駆動力
  • 曲げモーメント(板厚方向の温度勾配) — 追加の変形を誘発

板の熱座屈

🧑‍🎓

板の熱座屈の理論解はありますか?


🎓

四辺拘束の矩形板が一様に加熱される場合の臨界温度上昇:


$$ \Delta T_{cr} = \frac{k \pi^2}{12(1+\nu)} \cdot \frac{1}{\alpha} \left(\frac{t}{b}\right)^2 $$

🧑‍🎓

板厚 $t$ と幅 $b$ の比の2乗が入っている。薄い板ほど低い温度上昇で座屈する。


🎓

しかもこの式でヤング率 $E$ が入っていない! 理想的な拘束条件では臨界温度は材料の剛性に依存せず、幾何学的パラメータ($t/b$)と線膨張係数 $\alpha$ だけで決まる。


🧑‍🎓

材料に依存しないって直感に反しますが…。


🎓

完全拘束の場合、温度応力が $E\alpha\Delta T$、座屈応力が $E \cdot f(t/b)$ だから、比をとると $E$ が消える。ただし実構造は完全拘束ではないので、$E$ の影響は出てくる。


実構造での熱座屈問題

🧑‍🎓

熱座屈が問題になる実例を教えてください。


🎓

代表的な例:


構造状況特徴
鉄道レール夏季の温度上昇連続溶接レールの横座屈
パイプライン高温流体の輸送海底パイプラインの蛇行座屈(lateral buckling)
航空機外板超音速飛行時の空力加熱薄板パネルの座屈
宇宙構造日向/日影の温度差太陽電池パネル、アンテナ
溶接構造溶接時の局所加熱溶接変形(角変形、座屈変形)
火災時の建物鉄骨梁の加熱端部が拘束された梁の軸力座屈
🧑‍🎓

海底パイプラインの座屈って面白いですね。


🎓

海底パイプラインの熱座屈は石油・ガス業界の重要問題だ。高温の原油を通すと管が膨張しようとするが、海底の摩擦で拘束される。ある程度以上温度が上がると管が横方向に蛇行する(lateral buckling)。これを意図的に制御する「designed lateral buckling」という設計手法がある。


まとめ

🧑‍🎓

熱座屈の理論を整理します。


🎓

要点:


  • 熱座屈は「拘束された熱膨張」による圧縮応力が原因 — 自由膨張なら座屈しない
  • $\sigma_{th} = E\alpha\Delta T$ — 完全拘束時の熱応力
  • 拘束度の評価が鍵 — 実構造は完全拘束でも完全自由でもない
  • FEMでは熱応力解析→固有値座屈解析の2段階 — 温度分布が不均一でも対応可
  • 温度勾配は膜応力と曲げの両方を発生 — 片面加熱が特に危険

🧑‍🎓

「力をかけなくても座屈する」というのが熱座屈の怖いところですね。温度管理がそのまま構造安全性に直結する。


🎓

その通り。熱座屈は「荷重」が目に見えないから設計で見落とされやすい。温度変化が大きい環境で使われる構造は、常に熱座屈の可能性を検討すべきだ。


Coffee Break よもやま話

タイタニック号と安全率の教訓

「不沈」と謳われたタイタニック号は、低温でのリベット材の脆性破壊が沈没の一因とされています。現代の破壊力学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は、実務者が日々直面するこうした課題の解決を目指す研究開発プロジェクトです。

熱座屈解析の実務で感じる課題を教えてください

Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。

実務課題アンケートに回答する →