ジュール発熱の電磁-熱連成解析
理論と物理
ジュール発熱って、電流が流れると熱くなるだけですよね？ なぜわざわざ連成解析が必要なんですか？
いい質問だね。「電流が流れると熱くなる」はその通りなんだけど、話はそこで終わらないんだ。温度が上がると金属の電気抵抗率が増加する。例えば銅は1°Cあたり約0.4%抵抗が増える。すると同じ電流でも発熱量がさらに増え、温度がもっと上がり、抵抗がもっと増え...という正のフィードバックループが回り始める。
え、それって無限に温度が上がっちゃうんですか？
実際には放熱（対流・輻射・伝導）があるから、どこかで発熱と放熱が釣り合って定常温度に落ち着く。でも、その定常温度が何°Cになるかは「電磁場→発熱→温度変化→抵抗率変化→電磁場」というループを繰り返し計算しないと求められない。だからジュール発熱の精度を出すには電磁-熱の連成解析が必須なんだ。
具体的にどんな場面で問題になるんですか？
身近な例でいうと、EV（電気自動車）の800Vインバータの接続端子だね。400A以上の大電流が流れるバスバーで、局所温度が150°Cを超えないか確認する必要がある。ここを超えると絶縁材が劣化して短寿命になる。抵抗率の温度依存を無視して計算すると、実際より20〜30°C低い温度が出てしまい、設計不良を見逃す。あとはヒューズの溶断特性の予測、パワー半導体の基板設計、送電線のサグ（たわみ）計算なんかも典型例だよ。
なぜ連成解析が必要か
ジュール発熱解析が単純な「P = I²R」の計算で済まない理由は、電気伝導率 $\sigma$ が温度 $T$ の関数であることに集約される。定電流条件下でも、温度上昇により抵抗率が変化し、電流密度分布が再配分される。この双方向のカップリングを無視した解析は、以下の点で実測と乖離する。
- 定常温度の過小評価: 温度依存性を無視すると、銅バスバーで20〜30°Cの温度過小評価が生じ得る
- 電流集中の見落とし: 断面変化部・接続部で局所的に抵抗が増加し、ホットスポットが形成される
- 過渡応答の誤り: 短絡電流など過渡的な大電流での温度上昇曲線は、抵抗率変化により非線形になる
支配方程式
じゃあ実際にどんな数式で表されるんですか？
まず電磁場側。定常電流の場合、電荷保存則から連続の式が出る。そして発熱量は電流密度と電場の内積で決まる:
電場の支配方程式（定常電流近似）:
単位体積あたりのジュール発熱量:
$Q = J^2/\sigma$ と $Q = \sigma E^2$ は同じものなんですか？ なんか逆っぽく見えるんですけど...
するどいね。$\mathbf{J} = \sigma \mathbf{E}$ の関係を代入すると同じ式になる。使い分けとしては、電流密度が既知なら $J^2/\sigma$、電位が既知なら $\sigma E^2$ が便利だよ。実務では電位を未知数として解くことが多いから、$\sigma |\nabla V|^2$ の形で実装されることが多い。
熱場の支配方程式（エネルギー保存）:
ここで $\rho$ は密度、$c_p$ は定圧比熱、$k$ は熱伝導率、$Q$ はジュール発熱量（熱源項）。
各項の物理的意味
- $Q = \mathbf{J} \cdot \mathbf{E}$（ジュール発熱項）: 電場が電荷を加速し、格子振動（フォノン）との衝突で運動エネルギーが熱に変換される。電子ドリフト速度が高いほど、また散乱頻度が高いほど発熱量が増加する。【日常の例】スマホの充電ケーブルを触ると温かい——ケーブル内部の抵抗でジュール熱が発生している。急速充電ほど電流が大きく発熱も増える。
- $\nabla \cdot \mathbf{J} = 0$（電流連続の式）: 定常状態では電荷が蓄積しないことを表す。断面積が変わる部分では電流密度が反比例して変化し、狭い部分で局所的に発熱が集中する。【実務の例】バスバーのボルト穴周辺で電流密度が増大し、ホットスポットが形成される典型的なパターン。
- $\rho c_p \partial T / \partial t$（蓄熱項）: 材料の熱容量による温度変化の遅れ。銅のような熱容量の大きい材料は過渡的な温度上昇が緩やかになる。短絡電流の持続時間が数秒以下なら、この項の影響が大きい。
- $\nabla \cdot (k \nabla T)$（熱伝導項）: 温度勾配に比例した熱の拡散。銅は熱伝導率 $k \approx 400$ W/(m·K) と非常に高く、局所的な温度上昇を素早く均一化する。ステンレス鋼（$k \approx 16$）だとホットスポットが顕著に残る。
次元解析と単位系
| 変数 | SI単位 | 典型値（銅, 20°C） |
|---|---|---|
| 電気伝導率 $\sigma$ | S/m | $5.96 \times 10^7$ |
| 電気抵抗率 $\rho_e = 1/\sigma$ | $\Omega$·m | $1.68 \times 10^{-8}$ |
| 温度係数 $\alpha$ | 1/K | $3.93 \times 10^{-3}$（銅） |
| 電流密度 $J$ | A/m² | $10^5$〜$10^7$（バスバー） |
| 発熱密度 $Q$ | W/m³ | $10^4$〜$10^8$ |
| 熱伝導率 $k$ | W/(m·K) | 401（銅）, 16（SUS304） |
抵抗率の温度依存性
金属の電気抵抗率は温度に対してほぼ線形に増加する（金属の場合、室温付近で）:
ここで $\rho_0$ は基準温度 $T_0$ における抵抗率、$\alpha$ は抵抗温度係数。
$\alpha$ の値って材料ごとにどれくらい違うんですか？
代表的な値を表にまとめると一目瞭然だよ。
| 材料 | $\rho_0$ [$\mu\Omega$·cm] | $\alpha$ [1/K] | 備考 |
|---|---|---|---|
| 銅 (Cu) | 1.68 | $3.93 \times 10^{-3}$ | 最も使用頻度が高い導体 |
| アルミ (Al) | 2.65 | $4.29 \times 10^{-3}$ | 軽量化目的でバスバーに使用 |
| 銀 (Ag) | 1.59 | $3.80 \times 10^{-3}$ | 接点材料、高コスト |
| ニクロム (NiCr) | 108 | $1.7 \times 10^{-4}$ | ヒーター線、αが極小 |
| SUS304 | 72 | $9.4 \times 10^{-4}$ | 構造材料、高抵抗 |
ニクロムの $\alpha$ がすごく小さいですね。ヒーター線に使われる理由がわかりました——温度が上がっても抵抗があまり変わらないから、安定した発熱が得られるんですね？
その通り！ 逆に銅やアルミは $\alpha$ が大きいから、バスバーやケーブルの設計では温度依存性を無視できないんだ。100°Cの温度上昇で銅の抵抗は約40%も増加する。これを無視した解析は全く使い物にならないよ。
連成のフィードバック構造
ジュール発熱の電磁-熱連成には、明確な正のフィードバックループが存在する:
- 電流密度 $\mathbf{J}$ による発熱 $Q = J^2/\sigma(T)$
- 温度 $T$ が上昇
- 抵抗率 $\rho(T)$ が増加（$\sigma(T)$ が減少）
- 同じ電位差に対して電流分布が変化し、局所的な発熱密度が変化
- 1に戻る
定電流源と定電圧源で振る舞いが違ったりしますか？
鋭い視点だね。定電圧源の場合、温度上昇→抵抗増加→電流減少→発熱減少、と自己制限的に働く。逆に定電流源の場合、温度上昇→抵抗増加→発熱 $Q = I^2 R$ が増加→さらに温度上昇、と暴走方向に働く。パワーエレクトロニクスではスイッチング素子の駆動条件によって両者が混在するから、境界条件の設定を間違えると結果が全く変わってしまう。
身近な例で理解するフィードバック
白熱電球のフィラメント（タングステン）は室温で約5.3 $\mu\Omega$·cmの抵抗率だが、2500°Cの動作温度では約15倍に跳ね上がる。電源ON直後の突入電流が定常時の10倍以上になるのはこのためだ。CAE解析でこの過渡挙動を再現するには、抵抗率の温度依存性と電磁-熱連成が不可欠になる。
ジュールの法則から180年——P = I²R が現代ICを支配する
1841年、ジェームズ・プレスコット・ジュールは鉄線に電流を流す実験で「発熱量は電流の2乗と抵抗に比例する」ことを発見した。当時は電気の本質すら謎だった時代だ。180年後の今日、最先端の3nm半導体チップでは数十億個のトランジスタが1 cm²あたり100 W以上のジュール熱を発生させている。この熱密度は太陽表面の約1/6に相当する。ジュールが実験室で見出した法則は、EV・5G基地局・データセンターの設計を左右する現代のボトルネックとなっている。
数値解法と実装
弱連成と強連成
「弱連成」と「強連成」って聞いたことあるんですけど、どう使い分けるんですか？
ジュール発熱解析の連成方法は大きく3つに分類できる:
| 手法 | 概要 | 適用条件 | 精度 |
|---|---|---|---|
| 非連成（片方向） | 電磁場を解く→発熱量を熱解析に渡す（1回のみ） | $\Delta T < 50$°C、$\alpha \Delta T \ll 1$ | 低 |
| 弱連成（逐次反復） | 電磁場→熱→電磁場→熱...を収束まで反復 | $\Delta T = 50$〜$200$°C | 中〜高 |
| 強連成（モノリシック） | 電磁場と熱を1つの連立方程式系で同時に解く | $\Delta T > 200$°C、非線形性が強い | 最高 |
実務ではどれが一番よく使われますか？
圧倒的に弱連成（逐次反復法）だね。理由は3つある。まず既存の電磁場ソルバーと熱ソルバーをそのまま使い回せること。次にデバッグしやすいこと——各物理場を個別に検証できる。そして大規模モデルでもメモリが節約できること。モノリシック法は自由度が倍になるからメモリが厳しい。
弱連成の反復アルゴリズム（Gauss-Seidel型）:
- 初期温度 $T^{(0)}$ を設定（通常は環境温度）
- $\sigma(T^{(k)})$ を用いて電場方程式を解き、$V^{(k+1)}$ を得る
- $Q^{(k+1)} = \sigma(T^{(k)}) |\nabla V^{(k+1)}|^2$ を計算
- $Q^{(k+1)}$ を熱源項として熱伝導方程式を解き、$T^{(k+1)}$ を得る
- 収束判定: $\|T^{(k+1)} - T^{(k)}\|_\infty < \varepsilon$ ならば終了、そうでなければ2へ戻る
収束判定の $\varepsilon$ はどれくらいに設定すればいいですか？
温度の絶対値基準で0.1°C〜1°Cが実務的な目安だよ。EVバスバーのように許容温度マージンが10°C程度しかない場合は0.1°C、ヒーター設計のように数十°Cの余裕がある場合は1°Cで十分。反復回数は通常5〜15回で収束する。
有限要素定式化
電場の弱形式:
ここで $N_i$ は形状関数、$J_n$ は境界上の法線電流密度。
熱場の弱形式:
右辺第2項は対流境界条件（$h$: 熱伝達係数、$T_\infty$: 環境温度）。
電場と熱場で同じメッシュを使うんですか？
ジュール発熱解析の場合は、電場と熱場が同じ領域に存在するから、基本的に同一メッシュを使う。これが流体-構造連成（FSI）と大きく違う点で、データ転写の補間誤差を心配する必要がないんだ。COMSOL やAnsys Workbench では自動的に同一メッシュで連成計算される。
時間積分と収束制御
過渡解析における時間積分は、電場の応答時間（$\tau_e \sim \varepsilon/\sigma \sim 10^{-18}$ s）と熱場の応答時間（$\tau_T \sim \rho c_p L^2/k \sim 1$〜$100$ s）が極端に異なるため、工夫が必要。
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| 電場解析 | 定常解 | $\tau_e \ll \tau_T$ なので各タイムステップで定常解を使用 |
| 熱場の時間刻み | $\Delta t \leq \tau_T / 10$ | 後退Euler法（無条件安定）推奨 |
| 連成反復の収束基準 | $\Delta T < 0.1$〜$1$ °C | 温度の絶対値基準 |
| 連成反復の緩和係数 | 0.5〜0.8 | 発散防止、Aitken緩和も有効 |
| 最大連成反復回数 | 20〜30 | 超過時は緩和係数を下げる |
時間スケールの分離を利用した計算効率化
電磁場の応答速度（フェムト秒オーダー）は熱場（秒〜分オーダー）に比べて桁違いに速い。つまり熱場の各タイムステップにおいて、電磁場は瞬時に定常状態に達したと仮定できる。この「準定常近似」により、電磁場の過渡解析を省略し、各熱タイムステップで電場の定常解を1回だけ求めればよい。計算コストが劇的に削減される。交流電流の場合は、電磁場の1周期平均の発熱量を求めて熱場に渡す。
実践ガイド
解析フローと境界条件
実際にジュール発熱の連成解析を始めるとき、最初に何をすればいいですか？
実務的な解析フローはこうなる:
- 形状モデリング: CADデータから導体部分を抽出。ボルト穴・面取り・フィレットは省略せず残す（電流集中に影響）
- 材料定義: $\rho(T)$、$k(T)$、$c_p(T)$ の温度依存テーブルを入力。特に $\rho(T)$ は必須
- 電気境界条件: 電流入口面に電流密度 or 全電流値、出口面にGND（$V = 0$）を設定
- 熱境界条件: 外表面に自然対流（$h = 5$〜$15$ W/m²K）、必要に応じて輻射（$\varepsilon = 0.1$〜$0.9$）
- 連成設定: 弱連成の反復回数・収束基準を設定
- ポスト処理: 最高温度・電流密度分布・熱流束分布を確認
対流の熱伝達係数 $h$ ってどうやって決めるんですか？ 教科書には5〜25って書いてあるけど、幅が広すぎて...
いい悩みだね。現場では以下のように使い分ける:
- 密閉筐体内の自然対流: $h = 3$〜$8$ W/m²K
- 開放空間の自然対流（水平面上向き）: $h = 5$〜$15$ W/m²K
- 開放空間の自然対流（垂直面）: $h = 8$〜$20$ W/m²K
- 強制空冷: $h = 25$〜$250$ W/m²K（風速による）
実務的には、初回解析で $h = 10$ として感度分析し、$\pm 50$%変えても結論が変わらないか確認するのが安全だ。それで結論が変わるなら、CFD連成か実測が必要になる。
メッシュ戦略
メッシュの切り方で注意することはありますか？
ジュール発熱特有のメッシュ戦略がある。ポイントは電流密度の勾配が大きい箇所を細かくすること:
- 断面変化部: バスバーの幅や厚さが変わる部分。電流が集中するためメッシュを3層以上配置
- ボルト穴周辺: 電流が穴を迂回するため密度が急変。穴径の1/5以下の要素サイズが目安
- 接触面: 接触抵抗を評価する場合、接触面に少なくとも2層の薄い要素を配置
- 表皮効果（AC電流）: 表皮深さ $\delta = \sqrt{2/(\omega \mu \sigma)}$ 内に最低3要素
一方、導体中央部や電流が一様に流れる長い直線部は粗いメッシュでOK。計算コストの8割はメッシュの賢い粗密配分で決まるよ。
| 領域 | 要素サイズの目安 | 要素タイプ | 理由 |
|---|---|---|---|
| 電流集中部 | 代表長さの1/10以下 | 六面体2次推奨 | 電流密度の精度確保 |
| 接触面近傍 | 0.1〜0.5 mm | 六面体、薄い層状 | 接触抵抗の解像 |
| 表皮深さ領域 | $\delta/3$ 以下 | 六面体（境界層メッシュ） | AC電流分布の再現 |
| 一様電流部 | 5〜10 mm | 四面体でも可 | 計算コスト節約 |
接触抵抗のモデリング
バスバーのボルト接続部って、接触抵抗が重要だと聞いたんですけど、どうモデリングするんですか？
これは実務で最も頭を悩ませるポイントだ。接触抵抗は接触面圧・表面粗さ・酸化膜の状態に依存して、値が1〜2桁変動する。モデリング方法は3つある:
- 薄層要素: 接触面に厚さ数$\mu$mの仮想薄層を挿入し、その電気伝導率で接触抵抗を表現
- 界面抵抗条件: COMSOLの「Electric Contact」やAnsysの「Contact Resistance」で面の電気抵抗を直接指定
- 等価回路モデル: 接触部を集中定数の抵抗として扱い、発熱量を手計算で加える
いずれの方法でも、入力する接触抵抗値は実測が必須。典型的には銅-銅接触で $10^{-5}$〜$10^{-3}$ $\Omega$·cm² のオーダーだけど、表面処理やボルト締め付けトルクで大きく変わる。
実務事例: EVバスバー設計
実際のEV開発で、ジュール発熱の連成解析ってどう使われてるんですか？
典型的なケースを紹介しよう。800Vバッテリーパックのバスバー設計で、解析の目的は「連続400A通電時にバスバー温度が120°C（絶縁材の耐熱温度-30°Cマージン）を超えないこと」の確認だ。
解析条件:
- 材料: 無酸素銅 C1020（$\rho_0 = 1.72 \times 10^{-8}$ $\Omega$·m、$\alpha = 3.93 \times 10^{-3}$ /K）
- 断面: 30 mm × 3 mm（バスバー）、M6ボルト4本で接続
- 電流: 400 A DC 連続
- 環境: 60°C（筐体内温度）、$h = 8$ W/m²K（自然対流）
解析結果（抵抗率温度依存あり/なし比較）:
| 条件 | 最高温度 | ボルト穴周辺温度 | 電力損失 |
|---|---|---|---|
| 温度依存なし | 95°C | 102°C | 18.2 W |
| 温度依存あり（連成） | 118°C | 131°C | 24.7 W |
| 差分 | +23°C | +29°C | +36% |
29°Cも違うんですか！ 温度依存なしだと合格だけど、連成解析だと不合格に変わるケースですね。
まさにそう。この差が設計判断を完全にひっくり返す。対策としては断面を3 mmから4 mmに増厚するか、ボルト穴のレイアウトを変えて電流集中を緩和する。パラメトリック解析で最適な断面形状を探索するのがセオリーだ。実際にティア1サプライヤーでは、断面形状の最適化で重量15%削減と温度5°C低減を同時に達成した例がある。
EV開発の裏側——バスバー温度が1°C下がるたびに設計者が笑顔になる理由
EVのバッテリーパック設計で、バスバー温度が1°C下がるということは、絶縁材料のグレードを下げられる可能性を意味する。高耐熱絶縁材（PI: ポリイミドフィルム）は安価な絶縁材（PET: ポリエチレンテレフタレート）の3〜5倍のコストがかかる。数百万台規模の量産車では、絶縁材のコスト差が数億円のインパクトになる。だからバスバー温度の「最後の1°C」を詰めるためにジュール発熱の連成解析が使い倒されているのだ。
初心者が陥りやすい落とし穴: 輻射の無視
バスバー設計で「対流だけモデリングすれば十分でしょ」と輻射を省略するケースが多い。しかし100°C以上の高温では輻射放熱が全放熱の30〜50%を占めることがある。特に酸化した銅表面（$\varepsilon \approx 0.7$）や塗装面（$\varepsilon \approx 0.9$）では輻射が支配的になる場合もある。逆に研磨された銅面（$\varepsilon \approx 0.05$）では無視しても問題ない。放射率の見積もりを間違えると温度予測が10〜20°C狂う。
ソフトウェア比較
対応ツールと機能比較
ジュール発熱の連成解析ができるソフトって、具体的にどれですか？
主要なツールを比較するとこうなる。ジュール発熱は比較的シンプルな連成なので、対応ソフトは多い:
| ツール | 連成方式 | AC対応 | 接触抵抗 | 特徴 |
|---|---|---|---|---|
| COMSOL Multiphysics | 強連成（自動） | 周波数領域 | 薄層/界面 | GUIで連成設定が最も容易 |
| Ansys Workbench | 弱連成（Electric-Thermal） | Maxwell連携 | 界面抵抗 | 大規模モデルに強い |
| Abaqus | 強連成（*COUPLED THERMAL-ELECTRICAL） | 限定的 | GAPCON/GAPELEC | 非線形性に強い |
| JMAG-Designer | 弱連成（熱-電磁場） | 得意 | 限定的 | モータ・変圧器に特化 |
| MSC Marc | 強連成 | 限定的 | 接触テーブル | 非線形材料に強い |
| Simcenter STAR-CCM+ | 弱連成（CHT） | 電磁連携 | 界面 | 熱流体連成が得意 |
初心者がまず試すならどれがいいですか？
COMSOLが圧倒的に入門しやすい。「Electric Currents」と「Heat Transfer in Solids」の2つのモジュールを選ぶだけで、自動的に連成が設定される。チュートリアルに「Busbar Heating」という名前のまさにそのものの例題がある。ただしCOMSOLはライセンスがモジュール単位で高額（基本パッケージ+追加モジュール×2）なので、会社として大規模に展開するならAnsys Workbenchの方がコスパがいい場合もある。
ツール別設定手順
COMSOL Multiphysics:
- 「モデルウィザード」→ 3D → 「Electric Currents (ec)」+「Heat Transfer in Solids (ht)」+「Electromagnetic Heating (emh)」マルチフィジクス
- 材料に「Copper」を選択（温度依存性はデフォルトで含まれる）
- ec: 電流入口=「Normal Current Density」、出口=「Ground」
- ht: 外表面に「Heat Flux」→「Convective heat flux」
- 「Study」→「Stationary」で実行
Ansys Workbench:
- 「Steady-State Thermal-Electric」解析システムを選択
- Engineering Data で銅の $\rho(T)$ テーブルを入力
- 電流入口: Voltage/Current → 片端にVoltage=0、他端にCurrent
- 熱境界: Convection を外表面に適用
- Solve → 電流密度と温度の結果を確認
Abaqus:
- Step: *COUPLED THERMAL-ELECTRICAL, STEADY STATE
- 要素タイプ: DC3D8E（8節点六面体、熱-電気連成要素）
- 材料: *ELECTRICAL CONDUCTIVITY に $\sigma(T)$ テーブルを入力
- 電気境界: *CECURRENT（電流）+ *BOUNDARY, TYPE=ELECTRICAL（電位）
- 熱境界: *SFILM（対流）
COMSOLの「Busbar」チュートリアル——世界一有名な電磁-熱連成の例題
COMSOLを学んだことのある人なら、ほぼ全員が「Busbar」チュートリアルを解いた経験があるだろう。L字型バスバーにDC電流を流し、温度分布を求めるこの例題は、COMSOLのマルチフィジクス連成の「Hello World」とも言える存在だ。モデルの自由度はわずか数万だが、ジュール発熱の本質——電流集中・温度依存抵抗・放熱バランス——がすべて詰まっている。多くのエンジニアが「この例題を通じてマルチフィジクスの考え方を初めて理解した」と語る。
先端技術
マルチスケール連成
ジュール発熱の解析で、最新の研究トレンドってどんなものがありますか？
注目されているのはマルチスケール連成だね。例えばパワーモジュールの設計では:
- チップスケール（$\mu$m）: 半導体デバイス内部の電流分布とジュール発熱
- パッケージスケール（mm）: ワイヤボンド・はんだ層の接触抵抗と発熱
- システムスケール（cm〜m）: バスバー・筐体を含む全体の温度分布
これらを階層的にモデリングし、各スケールの情報を受け渡す手法が開発されている。全スケールを一括で解くのは計算コスト的に不可能だから、サブモデリングや均質化理論が使われる。
機械学習によるサロゲートモデル
最近はAI・機械学習も解析に使われてるって聞いたんですけど？
バスバー設計のパラメトリック最適化では、サロゲートモデルが威力を発揮する。数百パターンのジュール発熱連成解析を実行して訓練データを生成し、ニューラルネットワークで「形状パラメータ→最高温度」の関係を近似する。1回の連成解析に30分かかるとして、1000パターンの探索にはFEMで500時間かかるが、サロゲートモデルなら数秒で完了する。
またPINN（Physics-Informed Neural Network）で、支配方程式を損失関数に組み込んだ学習を行うアプローチも研究されている。ただし実務で使えるレベルにはまだ課題が多い。
その他の先端技術:
- デジタルツイン: リアルタイムセンサーデータ（電流値・温度）を連成解析モデルに入力し、運用中のホットスポットを予測
- トポロジー最適化: 導体の形状を最適化し、最小重量で温度制約を満たすバスバー設計を自動探索
- 確率論的解析: 接触抵抗のばらつき（$\pm 50$%以上）をモンテカルロ法で考慮し、信頼性設計に反映
トラブルシューティング
先生、ジュール発熱の連成解析で計算が上手くいかないケースって、どんなパターンがありますか？
典型的なトラブルパターンを3つ挙げよう。
1. 連成反復が収束しない
| 原因 | 診断方法 | 対策 |
|---|---|---|
| 緩和係数が大きすぎる | 反復ごとの温度変化が振動 | 緩和係数を0.3〜0.5に下げる |
| $\alpha \Delta T > 1$（抵抗変化が100%超） | $\rho(T)$ が基準値の2倍以上に | 温度範囲を確認、融点超えなら解析設定を見直す |
| メッシュ粗すぎ（電流集中部） | メッシュ細分化で結果が大きく変化 | 電流集中部のメッシュ細分化 |
| 境界条件の不整合 | 全電流のイン/アウトが不一致 | 電流保存を確認 |
「温度が振動する」って言われても、どう見ればいいかわからないんですけど...
COMSOLなら「Convergence Plot」で反復ごとの残差をグラフ化できる。残差が単調減少せずにジグザグしていたら振動の兆候。Ansysでは「Solution Information」タブのForce Convergenceグラフで確認できる。振動が見えたら、まず緩和係数を0.5に設定してみよう。それでもダメなら0.3まで下げる。
2. 温度が非物理的に上昇する
| 原因 | 診断方法 | 対策 |
|---|---|---|
| 放熱境界条件の未設定 | 全外表面が断熱になっている | 対流・輻射条件を追加 |
| 電流値の単位ミス（A↔mA） | 発熱量が1000倍以上になっている | SI単位で再入力 |
| 接触抵抗が過大 | 接触面で温度が突出 | 接触抵抗値を実測値で再設定 |
| 過渡解析の時間刻みが大きすぎる | 温度が階段状に急上昇 | $\Delta t$ を1/5に減らして再計算 |
一番多いミスは放熱境界条件の付け忘れだね。ジュール発熱では、熱の行き場がないと温度が際限なく上がり続ける。定常解析なら発散するし、過渡解析なら非現実的な高温（数千°C）が出る。まず外表面すべてに対流条件が設定されているか確認しよう。
3. 電流分布が偏る
| 原因 | 診断方法 | 対策 |
|---|---|---|
| 電気境界条件の誤り | 一部の面だけに電流が集中 | 入口/出口面の選択を再確認 |
| 材料の電気伝導率誤り | 異なる材料間で抵抗率の差が10桁以上 | 材料定義を確認（空気の伝導率を0にしていないか） |
| メッシュの不整合 | 部品間の接続面でメッシュが接合されていない | 共有トポロジー or 接触条件を設定 |
いろんなトラブルがあるんですね... 最初に確認すべきことを一言でまとめると？
ジュール発熱解析のデバッグで最初にやるべきことはエネルギー収支の確認だ。「投入電力 = $I^2R$ or $VI$ 」と「放熱量 = 対流+輻射+伝導の合計」が一致しているか。これが合っていれば、あとはメッシュ精度と局所的な問題に絞り込める。合っていなければ、境界条件の根本的な間違いがある。
「解析が合わない」と思ったら
- エネルギー収支を確認: 投入電力と放熱量が一致しているか
- 単位系を統一: SI単位で全入力を再確認。特に電気伝導率 [S/m] と抵抗率 [$\Omega$·m] の取り違いに注意
- 最小モデルで検証: 単純な直方体バスバーで解析解 $\Delta T = Q L / (2 k A)$（1次元近似）と比較
- メッシュ収束性: 要素数を2倍にして最高温度が5%以上変化するなら、メッシュ不足
- 温度依存性を一旦OFF: $\rho(T)$ を定数にして解析し、結果の妥当性を確認してから温度依存を戻す
なった
詳しく
報告