熱座屈解析 — トラブルシューティングガイド

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

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

1. 材料の線膨張係数 $\alpha$ が設定されているか — $\alpha = 0$ だと温度が変化しても熱応力がゼロ

2. 参照温度（応力フリー温度）が正しいか — 温度と参照温度の差が0だと熱応力がゼロ

3. Ansysで PSTRES, ON を忘れていないか — 応力剛性が計算されない

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

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

確認：

• 辺の境界条件で面内変位を拘束しているか — 面外変位（座屈変位方向）だけ自由にする
• 拘束が対称条件になっていないか — 一辺だけ拘束だと全体が片側に膨張して応力が低下
• 実構造の拘束度を正しくモデル化しているか

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

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

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

• Riks法に切り替える — 温度制御→弧長制御
• 安定化法を使う — *STATIC, STABILIZE
• 初期不整を導入する — 分岐点を滑らかに通過

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

典型的な問題は：

• 温度-材料テーブルの外挿 — 定義範囲外の温度でソルバーが外挿し、非現実的な材料特性になる
• 温度ステップの粗さ — 温度が急変する区間で増分が大きすぎると追従できない
• 応力フリー温度の扱い — 多段階の温度履歴で参照温度の管理が複雑になる

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

• $\alpha$ と参照温度の設定を確認 — ゼロ応力の原因はほぼこれ
• 面内拘束＋面外自由が基本 — 拘束不足だと座屈荷重が無限大に
• 一様温度でも応力は不均一になり得る — 境界条件次第
• 座屈点通過にはRiks法または安定化法 — 通常のNewton-Raphsonでは止まる
• 温度依存材料は外挿に注意 — テーブルの定義範囲を十分広く

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