熱座屈解析

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

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

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

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

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

温度の効果は2つある：

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

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

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

代表的な例：

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

要点：

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

2段階（または3段階）の手順になる。

定常熱伝導解析または非定常熱伝導解析で温度分布 $T(x,y,z)$ を求める。単純な一様加熱なら手動で温度を与えてもよい。

得られた温度場を構造解析に入力し、熱膨張による応力分布を計算。この段階で境界条件（拘束条件）が結果を大きく左右する。

Step 2の熱応力から幾何剛性マトリクスを構成し、固有値問題を解く。$\lambda$ が臨界温度荷重係数。

機械荷重と熱荷重が同時に作用する場合は、プリロードとして機械荷重を与え、追加の温度荷重に対する座屈を評価する。温度と機械荷重が独立に変動する場合は、固有値座屈の「比例荷重」の仮定が崩れるので注意が必要だ。

```

```

```

火災時のような大きな温度上昇では必須だ。鋼のヤング率は：

固有値座屈では温度依存材料を直接扱えない（線形の仮定だから）。この場合は非線形座屈解析が必要で、温度増分ごとに材料特性を更新しながら荷重-変位経路を追跡する。

シーケンシャル連成が標準的だ：

構造変形が温度場に影響を及ぼすケース（例：変形で耐火被覆が剥落する）でなければ、シーケンシャルで十分だ。火災時の鉄骨梁の座屈はほとんどの場合シーケンシャル連成で正確に評価できる。

要点：

分野ごとに典型的なケースを見ていこう。

鉄骨梁の端部がコンクリートの柱や壁に接合されている場合、梁の熱膨張が端部で拘束される。梁は軸方向に伸びられないため圧縮力が発生し、最終的にこの圧縮で座屈する。

実際のメカニズムは3段階：

Cardington火災実験（1995-96年、英国の実大鉄骨造建物の火災実験）で観察された現象だ。火災中に梁が大きくたわみ、冷却で収縮すると端部接合部に大きな引張力が発生する。接合部の回転能力が不足すると破断に至る。

できるが、計算コストが大きい。溶接の熱座屈（角変形、パネル座屈変形）は：

実際に重い。代替として固有ひずみ法がある。溶接部に事前に求めた固有ひずみ（塑性ひずみ）を入力し、1ステップの弾性解析で変形を求める手法だ。計算時間が1/100〜1/1000になる。

マッハ2以上の飛行で外板温度が200〜300°Cに達する。温度上昇による熱応力に加え、飛行荷重（気圧差、慣性力）も同時に作用するため、熱-機械複合荷重の座屈が問題になる。

SR-71ブラックバードやコンコルドの設計では熱座屈が重要な設計制約だった。現代の極超音速飛行体でも同様で、NASA X-59のような次世代機では熱座屈の評価が設計の中心にある。

最も重要であり最も難しいポイントだ。接合部の回転剛性や軸剛性は理想的な「固定」「ピン」「自由」のどれでもない。拘束度に対する感度分析（拘束剛性を変えて座屈荷重の変化を見る）を実施することを強く推奨する。

熱座屈は「熱伝導」と「構造座屈」の連成だから、両方の機能を持つソルバーが必要だ。

溶接の詳細シミュレーション（移動熱源＋弾塑性＋相変態）はAbaqusやSysweldが強い。ただし計算コストが膨大だから、専用ソルバーのSimufact Welding（旧MSC.Marc Mentat Welding）やVirfacのほうが実務的なことが多い。

溶接変形の簡易予測なら固有ひずみ法対応のツール（Virfac WELD、ASU固有ひずみ法）が圧倒的に速い。

いくつかある：

学術利用は無料だ。火災時の構造挙動に特化しているため、汎用FEMより設定が簡単。ユーロコードの耐火設計を行う場合は第一候補になる。

熱座屈は「熱」と「構造」の両方の知識が必要だから、専用ツールがユーザーをガイドしてくれるのは大きな利点だ。ただし結果の物理的解釈は最終的にエンジニアの責任だ。

3つの方向が活発だ。

Functionally Graded Materials（傾斜機能材料）は、板厚方向に材料特性が連続的に変化する材料だ。例えば片面がセラミック（耐熱）、反対面が金属（高強度）で、その間を滑らかに遷移する。

FGMの熱座屈は通常の板座屈と異なる：

もともとは宇宙往還機の耐熱タイルとして研究が始まった（日本のNAL/JAXAも初期の研究に貢献）。現在はタービンブレードの耐熱コーティング、原子力の構造材などにも応用されている。FGMの熱座屈は理論的に興味深い問題で、多数の論文が出ている。

熱弾性フラッターは、空力加熱→熱応力→剛性変化→空力弾性特性変化、というメカニズムだ。超音速パネルフラッターと熱座屈が同時に作用すると、臨界速度が大幅に低下する。

この問題は極超音速飛行体（マッハ5以上）の設計で特に重要だ。空力加熱で板が軟化し、座屈に近い状態でさらに空力弾性が不安定化する。3つの物理（熱伝導、空力、構造）の連成解析が必要になる。

だからROM（縮約モデル）の研究が盛んだ。構造の自由度を数百に縮約して空力弾性連成を解く。NASAのAFRL（空軍研究所）が中心的に研究している。

パイプラインの蛇行座屈モニタリングが実用化されている。海底パイプラインに光ファイバーセンサーを布設し、ひずみと温度をリアルタイムで計測。座屈の予兆（ひずみの急変）を検知して警報を出す。

鉄道レールでも同様のモニタリング研究がある。レールの温度と軸力をモニタリングし、座屈リスクを定量化する。気象データと連携して「座屈危険日」を予報するシステムも開発されている。

気温予報から翌日のレール温度を予測し、座屈リスクを評価する。リスクが高い日は列車の速度制限をかける。英国のNetwork Railが実運用している。

熱座屈特有のトラブルは、温度と構造の「つなぎ」に関するものが多い。

温度荷重が構造に正しく伝わっていない。確認項目：

「この温度では熱応力がゼロ」という基準温度だ。例えば施工時20°Cで組み立てた構造が100°Cに加熱されるなら、参照温度を20°Cにして荷重温度を100°Cにする。参照温度を100°Cにすると温度差がゼロになり、熱応力が出ない。

拘束が不十分で構造が自由に膨張できている。完全に自由な板を一様に加熱しても応力はゼロだから、座屈は起きない（$\lambda \to \infty$）。

確認：

その通り。熱座屈解析の境界条件の基本原則は面内を拘束して面外を自由にすることだ。面内の拘束が座屈の駆動力（圧縮応力）を生み、面外の自由度が座屈変形を許す。

境界条件が不均一だからだ。例えば板の一辺だけ面内を固定し、反対辺を自由にすると、一様温度でも応力は位置によって変わる。温度が一様でも応力が一様とは限らないのが熱座屈解析の特徴だ。

それが座屈点だ。非線形解析で座屈を通過するには：

Abaqusでは温度をRiks法の荷重パラメータにできる。$\lambda$ が温度のスケーリング係数になり、座屈点を通過した後の挙動も追跡できる。Nastranでは直接はできないので、等価な力に変換する工夫が必要。

典型的な問題は：

対策として、温度-材料テーブルは実際の使用範囲よりも広い温度で定義し、外挿を避ける。温度が急変する区間では増分を小さくする。

その通り。両方とも実構造の状態を正しく反映しているかが勝負だ。FEMの設定は簡単に見えるが、物理を理解していないと無意味な結果が出る分野だ。