熱座屈解析
理論と物理
熱座屈とは何か
熱座屈って温度が上がるだけで座屈するんですか？ 荷重がなくても？
そうだ。ざっくり言うと、拘束された構造体が加熱されると、膨張したくても膨張できないから圧縮応力が溜まる。その圧縮応力が臨界値を超えた瞬間、外部荷重ゼロでも座屈が起きるんだ。
え、荷重ゼロで座屈って直感に反しませんか？ どんな場面で実際に起きるんですか？
身近な例を3つ挙げよう：
- 鉄道レールの「太陽くねり」：夏場にレールが加熱されるが、枕木にボルト固定されて長手方向に伸びられない。圧縮力が蓄積して、ある温度を超えると横方向にウネウネと座屈する。運転見合わせの原因になるやつだ。
- 航空機の外板パネル：超音速飛行時の空力加熱で外板温度が上がるが、フレーム（骨格）で拘束されているため圧縮応力が発生し、パネル座屈が起きる。
- 配管の熱伸び：化学プラントの高温配管がサポート間で拘束されて圧縮され、座屈で配管が曲がる事故がある。
どれも「拘束 + 温度上昇 → 圧縮応力 → 座屈」という同じメカニズムだよ。
なるほど！ 「拘束」がキーワードなんですね。自由に伸びられれば座屈しないと。
その通り。自由端なら $\Delta L = \alpha L \Delta T$ だけ伸びて終わり。でも両端固定だと伸びがゼロに拘束されるから、その分がそっくり圧縮応力 $\sigma = E \alpha \Delta T$ に変わる。これが熱座屈の出発点だ。
熱応力の発生メカニズム
熱応力の式をもう少し詳しく教えてもらえますか？
一軸拘束の棒を考えよう。温度が $\Delta T$ 上がると自由熱ひずみは：
でも両端固定だと全ひずみがゼロ（$\varepsilon_{\text{total}} = 0$）だから、弾性ひずみが熱ひずみを打ち消す方向に生じる：
フックの法則から、拘束された構造に発生する圧縮応力は：
マイナスだから圧縮ですね。例えば鋼（$E = 200\,\text{GPa}$、$\alpha = 12 \times 10^{-6}\,/\text{K}$）で $\Delta T = 50\,\text{K}$ だと…
$\sigma = 200 \times 10^3 \times 12 \times 10^{-6} \times 50 = 120\,\text{MPa}$。たった50度の温度上昇で120 MPaもの圧縮応力が発生する。薄板ならこの応力で十分に座屈するよ。
臨界温度差の理論式
じゃあ、具体的に「何度上がったら座屈するか」はどうやって求めるんですか？
最も基本的なケースとして、四辺単純支持の矩形板が一様加熱される場合を考えよう。オイラー座屈理論の熱版だ。
臨界座屈応力はプレートの座屈公式から：
ここで $k$ は座屈係数（四辺単純支持・一軸圧縮で $k = 4$）、$t$ は板厚、$b$ は短辺幅だ。これと熱応力 $\sigma = E \alpha \Delta T$ を等置すると：
おお、ヤング率 $E$ が消えるんですね！
いい着眼点だ。熱応力が $E \alpha \Delta T$ で、臨界応力も $E$ に比例するから、$E$ が約分されて消える。つまり臨界温度差は板厚・幅・熱膨張係数・ポアソン比だけで決まり、ヤング率に依存しない。これは熱座屈特有の性質だよ。
具体的な数値で感覚をつかみたいです。例えばアルミ板（$\alpha = 23 \times 10^{-6}$、$\nu = 0.33$）で $t = 2\,\text{mm}$、$b = 200\,\text{mm}$ だとどうなりますか？
$k = 4$（四辺単純支持）として計算すると：
$\Delta T_{cr} = \dfrac{4 \times \pi^2 \times (0.002)^2}{12 \times 23 \times 10^{-6} \times (1 - 0.33^2) \times (0.2)^2} \approx 1.6\,\text{K}$
たった1.6度で座屈する。薄いアルミパネルがいかに熱座屈に弱いかが分かるだろう。航空機の外板設計で熱座屈が重要視される理由がこれだ。
各項の物理的意味
- $\alpha$（線膨張係数）：温度上昇1Kあたりの伸び率。アルミは鋼の約2倍の値を持つため、アルミ構造は熱座屈に弱い。
- $t/b$（板厚/短辺幅比）：板の細長比に相当。この比が小さい（薄くて幅広い板）ほど $\Delta T_{cr}$ は低くなり、座屈しやすい。$t^2$ に比例するので、板厚を2倍にすると $\Delta T_{cr}$ は4倍になる。
- $\nu$（ポアソン比）：$(1-\nu^2)$ の項は二軸応力状態の効果。ポアソン比が大きいと $\Delta T_{cr}$ は若干小さくなる。
- $k$（座屈係数）：境界条件とアスペクト比で変わる。四辺固定だと $k \approx 10.67$ で単純支持の約2.7倍になる。
仮定条件と適用限界
- 線形弾性材料（降伏前の範囲で有効）
- 小たわみ理論（たわみが板厚に対して十分小さい）
- 一様な温度分布（板面内・板厚方向とも均一）
- 初期不整なし（完全平板を仮定）
- 温度依存の材料物性は無視
- 実構造ではノックダウンファクター（0.5〜0.8程度）の適用が必要
次元解析と単位系
| 変数 | SI単位 | 代表値・メモ |
|---|---|---|
| $\alpha$（線膨張係数） | 1/K | 鋼: 12×10⁻⁶、アルミ: 23×10⁻⁶、チタン: 8.6×10⁻⁶ |
| $E$（ヤング率） | Pa | 鋼: 200 GPa、アルミ: 70 GPa（$\Delta T_{cr}$ には直接影響しない） |
| $t$（板厚） | m | 航空機外板: 1〜3 mm、圧力容器: 10〜50 mm |
| $\sigma_{cr}$（臨界座屈応力） | Pa | $E \alpha \Delta T_{cr}$ で逆算可能 |
固有値問題としての定式化
FEMではどうやって定式化するんですか？
線形座屈解析は固有値問題として定式化される。熱応力による応力剛性マトリクス $[K_\sigma]$ を使って：
各記号の意味はこうだ：
- $[K_0]$：通常の弾性剛性マトリクス
- $[K_\sigma(T)]$：温度 $T$ による熱応力から計算された応力剛性マトリクス（幾何剛性マトリクスとも呼ぶ）
- $\lambda$：固有値＝臨界温度荷重倍率
- $\{\phi\}$：座屈モード形状
$\lambda = 1$ なら与えた温度荷重ちょうどで座屈。$\lambda = 0.5$ なら半分の温度で座屈するから危険。$\lambda > 1$ なら安全側ということだ。
なるほど、$\lambda$ が安全率みたいな役割を果たすんですね。
境界条件と臨界温度の関係
境界条件を変えると臨界温度はどのくらい変わるんですか？
座屈係数 $k$ が境界条件で大きく変わるから、臨界温度も比例して変わる。正方形板（$a/b = 1$）の場合で比較しよう：
| 境界条件 | 座屈係数 $k$ | $\Delta T_{cr}$ の倍率（SS基準） |
|---|---|---|
| 四辺単純支持 (SSSS) | 4.0 | 1.0 |
| 四辺固定 (CCCC) | 10.67 | 2.67 |
| 2辺固定+2辺自由 (CCFF) | ≈2.2 | 0.55 |
| 全辺自由端（面内拘束あり） | 1.0 | 0.25 |
四辺固定にするだけで臨界温度が2.7倍になる。設計で固定端を増やすのが熱座屈対策の基本だよ。逆に自由端が多い構造は弱い。
固定条件を変えるだけで臨界温度が4倍近く違うのは設計上大きいですね。
数値解法と実装
FEM 2ステップ解析手順
FEMで熱座屈を解析するとき、具体的にはどういう手順になるんですか？
基本は2ステップで進める：
Step 1：熱応力解析（Static Analysis with Thermal Load）
- 温度分布 $\Delta T(\mathbf{x})$ を荷重として入力
- 熱ひずみ $\varepsilon_{\text{th}} = \alpha \Delta T$ を外力等価荷重に変換
- $[K]\{u\} = \{F_{\text{th}}\}$ を解いて応力場 $\sigma_{ij}$ を求める
Step 2：線形座屈固有値解析（Linear Buckling Eigenvalue Analysis）
- Step 1 の応力場から応力剛性マトリクス $[K_\sigma]$ を構築
- 固有値問題 $([K_0] + \lambda [K_\sigma]) \{\phi\} = \{0\}$ を解く
- 最小正の固有値 $\lambda_1$ が臨界温度倍率
Step 1 で温度を100度で入れて、$\lambda_1 = 0.6$ が出たら、60度で座屈するってことですか？
その通り。$\Delta T_{cr} = \lambda_1 \times \Delta T_{\text{ref}}$ だから、$0.6 \times 100 = 60$ 度が臨界温度差になる。実務では安全率1.5〜2.0を掛けるから、$\Delta T_{\text{allow}} = 30 \sim 40$ 度という設計になるね。
温度分布が非一様な場合
一様加熱でなく、片面加熱や局所加熱の場合は、Step 1 で定常/非定常の熱伝導解析を先に実行して温度分布を求め、その結果をStep 1の構造解析に温度荷重としてマッピングする。つまり実質3ステップ（熱伝導 → 熱応力 → 座屈固有値）になる。
要素選択と定式化
熱座屈解析ではどんな要素を使えばいいですか？
構造の形状によって使い分ける：
| 構造タイプ | 推奨要素 | 注意点 |
|---|---|---|
| 薄板パネル | シェル要素（S4R, SHELL181等） | 板厚方向5〜7積分点で曲げ勾配を捕捉 |
| 厚板・ブロック体 | ソリッド要素（C3D20R, SOLID186等） | 板厚方向に最低3〜4層のメッシュ |
| レール・梁 | ビーム要素（B31, BEAM188等） | 薄肉部の局所座屈は捕捉不可、シェルと併用 |
| 配管 | パイプ/シェル要素 | 周方向の座屈モードを捕捉するには十分な周方向分割が必要 |
シェル要素が多いんですね。板厚方向の積分点数って何に効くんですか？
板厚方向に温度勾配がある場合、積分点が少ないと曲げ応力の分布を正確に捕捉できない。例えば片面加熱では板厚方向に線形の温度分布が生じるが、積分点が3点だけだと曲げ熱応力を過小評価することがある。実務では最低5点、できれば7点を推奨している。
非線形ポスト座屈解析
線形座屈解析だけで十分ですか？ 座屈後の挙動も知りたい場合はどうするんですか？
いい質問だ。線形座屈は「いつ座屈するか」を教えてくれるが、「座屈後にどうなるか」は分からない。ポスト座屈の挙動を追うには非線形静解析（幾何学的非線形）が必要だ。
手順はこうなる：
- 線形座屈解析で座屈モード形状 $\{\phi\}$ を取得
- そのモード形状を初期不整として幾何形状にスケーリングして加える（典型的には板厚の10〜50%）
- 温度荷重を段階的に増加させる非線形静解析（NLGEOM=ON）を実行
- 荷重-変位曲線の追跡（Riks法 or Arc-length法）
初期不整を入れないとどうなるんですか？
完全な幾何形状だと分岐点で解が一意に決まらず、ソルバーが座屈を捕捉できない。非線形解析では「完璧な構造は現実に存在しない」という前提で初期不整を入れるのがセオリーだ。実構造にも製造誤差による微小な凹凸があるから、むしろ初期不整入りのモデルの方がリアルだよ。
ソルバー設定の要点
ソルバーの設定で気をつけるべきポイントはありますか？
| 設定項目 | 推奨値 | 理由 |
|---|---|---|
| 求める固有値数 | 10〜20個 | 最低モードだけでなく高次モードの影響も確認 |
| 固有値ソルバー | Lanczos法 | 大規模モデルで安定・高速 |
| 温度参照温度 | 応力フリー温度 | 設定ミスで全応力がオフセットする事故が多い |
| NLGEOM（非線形） | ON | ポスト座屈追跡時は必須 |
| Riks法の初期弧長 | 荷重の5〜10% | 大きすぎると分岐点を飛び越す |
実践ガイド
実務での熱座屈事例
実務ではどんな場面で熱座屈解析が必要になるんですか？
代表的な分野と事例をまとめよう：
| 分野 | 事例 | 温度範囲 | 特徴 |
|---|---|---|---|
| 鉄道 | ロングレールの太陽くねり | ΔT ≈ 30〜50 K | 管理基準温度で徐行規制を発動 |
| 航空宇宙 | 超音速機の外板パネル | ΔT ≈ 100〜300 K | 空力加熱、スティフナー配置で対策 |
| プラント | 高温配管の座屈 | ΔT ≈ 200〜500 K | エクスパンションループで熱伸び吸収 |
| 電子機器 | PCB基板の反り | ΔT ≈ 150〜250 K | はんだリフロー時、異材接合のCTE差 |
| 原子力 | 原子炉容器の熱過渡座屈 | ΔT ≈ 50〜200 K | 急速昇温時の熱衝撃と座屈の複合 |
| 土木 | 鋼橋の路面鋼板座屈 | ΔT ≈ 20〜40 K | 日射による不均一加熱 |
鉄道レールの場合、具体的にはどう管理しているんですか？
日本のJRでは「レール温度管理基準」があって、敷設温度（通常25〜30℃で応力ゼロに設定）からの温度上昇量を監視している。レール温度が管理値（例えば60℃）に近づくと徐行規制を発令し、超えると運転見合わせになる。最近はレール内部にセンサーを埋め込んで軸力をリアルタイム監視するシステムも導入されているよ。
モデリングワークフロー
実際の解析フローを一通り教えてもらえますか？
一般的な熱座屈解析のフローはこうだ：
- 形状モデル作成：CADから形状をインポート、解析に不要なフィレット・穴の簡略化
- 材料物性の定義：$E(T)$、$\alpha(T)$、$\nu(T)$ の温度依存性を入力。高温域で降伏応力が低下する場合は塑性データも必要
- メッシュ生成：座屈が予想される領域（薄板部、拘束近傍）を細かく分割。板厚方向の積分点数を確認
- 境界条件の設定：拘束条件（固定、対称、周期）+ 温度荷重（一様 or 分布）+ 参照温度の設定
- Step 1 実行：熱応力解析（静解析）→ 応力分布を確認し、物理的に妥当か検証
- Step 2 実行：線形座屈固有値解析 → $\lambda_1$ と座屈モード形状を確認
- 結果の検証：理論解との比較、メッシュ収束性の確認、モード形状の妥当性チェック
- （必要に応じて）Step 3：初期不整入りの非線形ポスト座屈解析
商用ツール別の実装
各ソルバーでの具体的な設定方法を教えてもらえますか？
Abaqusの場合：
- Step 1：*STEP, *STATIC + *TEMPERATURE で温度荷重を適用
- Step 2：*STEP, *BUCKLE, EIGENSOLVER=LANCZOS
- 初期不整：*IMPERFECTION で座屈モードをスケーリング
- ポスト座屈：*STEP, NLGEOM=YES, *STATIC, RIKS
Ansys Mechanicalの場合：
- Step 1：Static Structural + Thermal Condition
- Step 2：Eigenvalue Buckling（pre-stress環境をON）
- APDL：PSTRES,ON → SOLVE → /SOLU → ANTYPE,BUCKLE → BUCOPT,LANB,10 → SOLVE
Nastranの場合：
- SOL 105：線形座屈（TEMPERATURE(LOAD)で温度荷重指定）
- SOL 106：非線形静解析でポスト座屈追跡
- EIGRL カードで求める固有値数とシフト値を指定
検証と妥当性確認
FEM結果が正しいかどうか、どうやって確かめるんですか？
検証には3段階のアプローチを取る：
- 理論解との比較：単純な形状（矩形板、円筒殻）では閉形式解が存在するから、まずそれで誤差5%以内であることを確認
- メッシュ収束性：メッシュを2倍に細かくして $\lambda$ の変化が2%以内であることを確認。3段階以上のメッシュ密度で Richardson外挿を適用できるとベスト
- 実験データとの照合：可能であれば熱座屈試験の実験値と比較。FEMは理想構造を仮定するため、通常は実験値よりも10〜30%高い $\Delta T_{cr}$ を出す（ノックダウンファクターの根拠）
FEMが実験より高く出るのは、完全な形状を仮定しているからですね。
その通り。実構造には製造時のわずかな凹凸・残留応力・板厚のばらつきがある。だからNASAのガイドラインではノックダウンファクター 0.65 を推奨している。つまり FEM の $\Delta T_{cr}$ に 0.65 を掛けた値を許容値とする。保守的だけど安全側に倒す設計思想だ。
ソフトウェア比較
主要ツールの機能比較
熱座屈解析ができるソフトって、どれがいいんですか？
主要4ツールを比較しよう：
| 機能 | Abaqus | Ansys Mechanical | Nastran | COMSOL |
|---|---|---|---|---|
| 線形座屈（固有値） | ◎ BUCKLE | ◎ BUCKLE | ◎ SOL 105 | ◎ |
| 熱-構造連成 | ◎ 順次連成 | ◎ Workbench連成 | ○ TEMP(LOAD) | ◎ 完全連成 |
| 非線形ポスト座屈 | ◎ Riks法 | ◎ Arc-Length | ○ SOL 106 | ○ |
| 温度依存物性 | ◎ テーブル入力 | ◎ テーブル入力 | ◎ MATTi | ◎ 関数入力 |
| 初期不整の導入 | ◎ *IMPERFECTION | ○ UPGEOM | ○ 手動 | △ |
| 大規模並列 | ◎ MPI | ◎ MPI | ◎ MPI | ○ |
Abaqusは初期不整の導入とRiks法が使いやすく、熱座屈のポスト座屈解析では最も実績がある。Nastranは航空宇宙分野で圧倒的シェアで認証実績が豊富。COMSOLはマルチフィジクス連成が得意だが、大規模モデルではやや不利だ。
Nastran SOL 105 vs SOL 106
NastranのSOL 105とSOL 106はどう使い分けるんですか？
ざっくり言うと：
- SOL 105（線形座屈）：固有値問題を解いて臨界荷重倍率と座屈モードを高速に求める。設計初期段階のスクリーニングに最適
- SOL 106（非線形静解析）：幾何学的非線形を考慮して、初期不整入りのモデルのポスト座屈挙動を追跡する。最終的な耐荷力評価に使う
熱座屈の場合、Case Control Section に TEMPERATURE(LOAD) = SET_ID を指定して温度荷重を取り込むのがポイント。TEMPERATURE(INITIAL) と混同すると参照温度の設定が狂って結果がメチャクチャになるから注意だ。
夏の線路がウネウネ曲がる理由——太陽くねりの物理
真夏の炎天下、鉄道の線路が大きくS字にうねって運転見合わせになるニュースを見たことがある人も多いでしょう。あれが「太陽くねり」、英語では "sun kink" と呼ばれる熱座屈現象です。日本のロングレール（200m以上の連続レール）では、レール温度が敷設基準温度（通常30℃前後）から33℃以上上昇すると座屈リスクが高まるとされます。レールの断面積が約77cm²、鋼の $E\alpha$ が約2.4 MPa/K なので、ΔT=40Kで約96 MPa、軸力にして約740 kNもの圧縮力が発生します。これはレールの横方向座屈荷重を超えうる値です。対策としては、バラスト（砂利）の横抵抗力確保、レールの定期的な応力開放（レール切断・再溶接）、温度センサーによるリアルタイム監視が行われています。
先端技術
傾斜機能材料（FGM）の熱座屈
最近の研究ではどんなテーマがホットなんですか？
ここ10年で特に注目されているのが傾斜機能材料（FGM: Functionally Graded Materials）の熱座屈だ。FGMは板厚方向に材料組成が連続的に変化する材料で、例えばセラミックス面（耐熱側）から金属面（構造強度側）へグラデーションする。
FGMでは板厚方向に $E(z)$、$\alpha(z)$、$\rho(z)$ が連続的に変化するため、従来の均質板の $\Delta T_{cr}$ 式がそのまま使えない。高次せん断変形理論（HSDT）やレイヤーワイズ理論で定式化する必要がある。研究によれば、セラミックスリッチ面を高温側にすると $\Delta T_{cr}$ が20〜40%向上するケースが報告されている。
確率論的初期不整感度解析
初期不整の影響をもっと厳密に評価する方法はありますか？
最先端のアプローチは確率論的初期不整感度解析だ。従来は「1次モードを板厚の何%で入れる」という決定論的な方法だったが、最近は以下の手法が使われている：
- モンテカルロ法：初期不整の形状・振幅をランダムに生成して数千ケースの座屈解析を実行し、$\Delta T_{cr}$ の確率分布を求める
- カルーネン-レーヴェ展開：実測データから初期不整のランダム場を統計的にモデル化
- PINN（物理インフォームドニューラルネットワーク）：FEMの代わりにニューラルネットで座屈荷重を高速予測するサロゲートモデルを構築
特にロケットや航空機の薄肉シェルでは、初期不整に対する感度が極めて高いため、このアプローチの重要性が増している。
トラブルシューティング
よくある失敗パターン
熱座屈解析で「やりがちな間違い」って何がありますか？
現場で最も多い失敗トップ5を挙げよう：
1. 参照温度（ストレスフリー温度）の設定ミス
これが断トツで多い。参照温度を室温（20℃）にすべきところを0℃にしていたり、設定し忘れたりする。結果として全節点に余計な熱応力がオフセットされ、$\lambda$ がまったく変わってくる。
2. 拘束条件の過不足
剛体モードが残っていると固有値解析が発散する。逆に拘束しすぎると「拘束点の反力で座屈が抑制される」という非現実的な結果になる。特に対称モデルで対称面の拘束を間違えると、反対称モードが出なくなる。
3. 温度依存材料物性の無視
高温域では $E(T)$ が低下するため、定温物性で計算すると $\sigma_{cr}$ を過大評価する。特にアルミは200℃で $E$ が約15%低下する。
4. メッシュが粗すぎる
座屈モードの半波長に対して要素が4個未満だと、モード形状を捕捉できずに $\lambda$ が過大評価される。
5. 初期不整なしのポスト座屈解析
完全形状で非線形解析をすると、座屈分岐点でソルバーが収束失敗する。初期不整を入れ忘れているケースが非常に多い。
ソルバー別エラーと対策
具体的にどんなエラーメッセージが出るんですか？
Nastran：
- FATAL 2012（特異剛性マトリクス）：剛体モード未拘束 → 座屈解析の前段（静解析）の拘束条件を確認。AUTOSPC=YES で微小拘束を自動追加する手もある
- USER WARNING 5291（要素品質不良）：歪んだ要素が座屈モードを歪める → メッシュ修正
- INFO 4158（負の固有値）：圧縮ではなく引張状態で固有値解析している → 温度荷重の符号・参照温度を確認
Abaqus：
- Zero pivot in factorization：拘束不足で剛性マトリクスが特異 → 境界条件を追加
- Too many attempts（Riks法）：ポスト座屈で収束失敗 → 初期弧長を小さくする、初期不整の振幅を調整
- Excessive element distortion：ポスト座屈で変形が大き過ぎる → NLGEOM=ONを確認、要素の歪み許容値を緩和（ただし結果の妥当性に注意）
デバッグチェックリスト
「解析結果がおかしい」と思ったとき、何からチェックすればいいですか？
以下の順序でチェックすると効率的だ：
- Step 1（熱応力解析）の結果を先に確認：応力分布が物理的に妥当か？ 拘束点の反力の合計がゼロに近いか？ 一様加熱なら面内応力が一様になっているか？
- 参照温度の確認：温度荷重と参照温度の差が $\Delta T$ になっているか？ 「参照温度＝応力フリー温度」の設定を確認
- 座屈モード形状の可視化：物理的に意味のあるモードか？ 局所的な要素の歪みがモードとして出ていないか？
- 固有値の符号：負の固有値は「圧縮ではなく引張で座屈する温度」を意味する（通常は不要）。正の最小固有値のみを評価
- 単位系の整合性：$\alpha$ の単位が 1/K か 1/℃ か、$E$ の単位系とメッシュの寸法単位が一致しているか
熱座屈の全体像がかなりクリアになりました。まずは単純な矩形板モデルで理論解と照合するところから始めてみます！
いいね。四辺単純支持の正方板は理論解が明確だから最高の練習材料だよ。まずStep 1で応力が $E \alpha \Delta T$ と一致することを確認して、Step 2 で $\lambda$ が理論値の $\pm 5\%$ 以内に収まれば、解析環境のセットアップは合っている。そこから実モデルに展開していけばいい。頑張れ！
なった
詳しく
報告