定常熱伝導 — CAE用語解説
定常熱伝導
熱解析で「定常」と「非定常」を選ぶ場面がありますが、定常熱伝導ってどういう状態を指すんですか？
理論と物理
定常熱伝導の支配方程式
「定常熱伝導」の「定常」って、具体的にどういう状態を指すんですか？時間が止まっているってことですか？
いい質問だ。時間が止まるわけではなく、温度場が時間に対して変化しなくなる状態だ。エンジンが十分に暖機された後や、電子基板が通電してから十分に時間が経過した状態をイメージしてほしい。非定常の熱伝導方程式
時間微分項がゼロになる条件は、具体的にはどんな時ですか？例えば、熱源Qが突然消えたらどうなりますか？
熱源Qが突然消えれば、それは境界条件が変わったことになり、システムは新たな定常状態を目指して非定常変化を始める。定常状態が成立するのは、「境界条件（温度や熱流束）」と「内部発熱Q」が全て時間的に一定で、かつ初期の過渡現象が完全に収束した後だ。例えば、CPUのTDP（熱設計電力）が65Wで一定に保たれ、冷却ファンの風量も一定なら、数分から十数分後にはチップの温度分布はほぼ定常とみなせる。
熱伝導率kが温度によって変わる材料、例えば炭素鋼だと定常解析はどう扱うんですか？kが一定じゃないと、先ほどの方程式は単純なラプラス方程式にならないですよね。
その通りで、その場合は非線形問題になる。支配方程式は
数値解法と実装
離散化とソルバー設定
定常熱伝導の有限要素法（FEM）では、質量マトリックスは作らないんですよね？時間項がないから。
正解だ。非定常解析では
では、ソルバーは直接法と反復法、どちらを選ぶのが一般的ですか？また、収束判定の「残差」は何を見ればいいんですか？
問題規模による。10万自由度以下の比較的小規模で、特に複数の荷重ケースを解く場合は、AnsysのSPARSE直接法（Intel MKL PARDISOベース）が安定して速い。100万自由度を超える大規模問題では、反復法（例えばCG法）に前処理（ICCGやAMG）を組み合わせるのが主流だ。収束判定は、無次元化された熱流束の残差ノルムを見る。Ansys Mechanicalではデフォルトで1.0e-4に設定されているが、熱流束が極端に小さい領域に注目する場合は1.0e-6まで厳しくする場合もある。
熱伝導率が異なる材料が接する界面（例えばアルミニウムと樹脂）は、FEMではどう扱われるんですか？特別な要素が必要ですか？
特別な「界面要素」を挿入する方法と、節点を共有させる「完全接合」とする方法がある。現実的には接触熱抵抗を無視できるほど密着していると仮定し、異なる材料要素が節点を共有するモデルが最も一般的だ。熱流束の連続性は、要素の積分式の中で自然に満たされる。ただし、実際の剥離やエアギャップを考慮するなら、Abaqusの「Gap Conductance」や、接触面に熱伝達率（例えば 1000 W/(m²·K)）を定義する必要がある。
実践ガイド
解析ワークフロー
定常熱伝導解析を始める時、最初に確認すべき材料物性値は何ですか？全部揃えようとすると膨大になりそうです。
最低限必要なのは「熱伝導率 k」だけだ。密度ρや比熱c_pは定常解析では不要（非定常や熱応力連成では必要）。ただし、熱伝導率が温度依存する場合は、想定温度範囲（例えば20°Cから最大到達温度まで）で数点のデータを用意する。アルミニウムA6061なら、室温で約167 W/(m·K)だが、150°Cでは約180 W/(m·K)と少し上がる。このデータは「JSMEデータブック」や「JAHM熱物性データベース」で確認できる。
境界条件で、「対流」と「放射」はどちらも熱を逃がす効果ですが、定常解析では両方入れるべきですか？計算が複雑になりませんか？
影響度を見極めることが重要だ。まず、強制対流（ファン冷却）があれば、それが支配的で放射は無視できることが多い。自然対流のみの場合、放射の寄与は無視できない。目安として、物体表面温度が周囲温度より50K以上高く、表面がアルミ裸面（放射率ε≈0.1）なら対流が主役だ。同じ条件で表面が塗装黒（ε≈0.9）なら、放射による放熱は対流に匹敵するかそれを上回る。COMSOLやAnsys Steady-State Thermalでは、対流と放射を同じ境界面に簡単に追加できる。
メッシュの細かさは、どう決めればいいですか？温度勾配が大きいところを細かくするのはわかりますが、具体的な指標は？
熱流束（温度勾配×熱伝導率）を追うといい。まず粗いメッシュで一度解き、要素ごとの熱流束をコンター表示する。隣接する要素間の熱流束の差が大きい領域、特に発熱体の角や異種材料界面、対流境界付近を重点的に細分化する。実用的には、メッシュを段階的に細かくして、注目点の温度変化が1%以内（例えば目標温度85°Cに対して0.85°C以内）に収束することを確認する「メッシュ収束性検証」を行う。これがCAEの基本だ。
ソフトウェア比較
各ソフトウェアの特徴
Ansys MechanicalとAbaqus/Standardの定常熱伝導解析で、根本的に違うところはありますか？どちらもFEMですよね。
解法のコアは同じだが、ユーザーインターフェースと連成解析への拡張性が大きく異なる。Ansys Mechanicalは「Steady-State Thermal」モジュールが独立しており、設定が直感的で、後段の静構造解析への温度荷重の受け渡しがワンクリックでできる。一方、Abaqus/Standardでは「*HEAT TRANSFER, STEADY STATE」というキーワードで定義し、同じINPファイル内で「*COUPLED TEMPERATURE-DISPLACEMENT」解析に拡張できる柔軟性が高い。また、Abaqusのユーザーサブルーチン（DFLUX, FILM）は高度な境界条件を自作できる点が強みだ。
無償で使えるCAEソフト、例えばCalculixやOpenFOAMの熱伝導解析は、有償ソフトと比べて何が足りないんですか？
Calculix（Abaqus互換のFEMソルバー）は定常熱伝導の基本機能は備えるが、大きく2点不足する。第一に、放射の角係数（View Factor）の自動計算機能が貧弱で、複雑形状では実用が難しい。第二に、熱物性の温度依存データをグラフィカルに設定・管理するGUIがなく、テキスト入力を要求される。OpenFOAMはCFDが主なので、固体のみの熱伝導解析には過剰で、メッシュ生成から学習コストが高い。有償ソフトはこれらの「実務効率化」の部分で価値を生んでいる。
COMSOL Multiphysicsは「マルチフィジックス」が売りですが、単純な定常熱伝導だけを解く場合の利点は何ですか？
「弱形式」や「係数型PDE」を直接編集できる柔軟性が最大の利点だ。例えば、熱伝導率が場所の関数（
トラブルシューティング
よくあるエラーと対策
解析を実行したら、「ソルバーが収束しません」というエラーが出ました。まず何を疑えばいいですか？
まず、熱伝導率kにゼロまたは極端に小さい値（1e-10など）が設定されていないか確認せよ。次に、境界条件の矛盾だ。例えば、全ての表面を「断熱」にした上で内部発熱Qを与えると、熱が逃げ場を失い、温度が無限大に発散する（数値的にはオーバーフロー）。少なくとも1か所は、温度固定条件か、熱を逃がす対流/放射条件が必要だ。材料の単位系（W/(m·K)）と発熱Qの単位（W/m³）の不一致もよくある原因だ。
温度依存の熱伝導率を入れたら、解が振動したり、現実離れした高温（数千℃）になりました。これはなぜですか？
それは非線形ソルバーの発散だ。原因は2つ考えられる。第一に、熱伝導率が温度上昇で急激に減少する材料（一部の断熱材）の場合、負のフィードバックがかかり不安定になる。第二に、初期推定値が悪い。対策として、まず熱伝導率を一定とした線形解析で得た温度場を、非線形解析の初期値として与える。それでもダメなら、Ansysでは「Nonlinear Controls」で荷重ステップを細かく分割し、一度に全負荷をかけずに段階的に載荷する「Ramped loading」を適用する。
異方性熱伝導（例えばグラファイトシート）を設定したら、熱流束の向きがおかしくなりました。異方性の軸定義で気をつけることは？
これも頻出問題だ。異方性は「材料座標系」に対して定義される。多くの学生はグローバル座標系（X,Y,Z）で定義したつもりが、ソフトウェアのデフォルトでは「要素座標系」や「局所座標系」が割り当てられていて、要素の向きによって軸がバラバラになる。対策は、Abaqusなら「*ORIENTATION」で、Ansysなら「Coordinate System」を明示的に作成し、それを材料プロパティに紐付けることだ。グラファイトシートでは面内方向（x,y）の熱伝導率が約1500 W/(m·K)、厚さ方向（z）は10 W/(m·K)程度と、3桁も違うので、設定を間違えると結果は全く違うものになる。
関連トピック
なった
詳しく
報告