ワイヤレス充電の電磁-熱連成シミュレーション
理論と物理
ワイヤレス充電って熱くなりますよね？ シミュレーションで何がわかるんですか？
いい質問だ。磁界結合で電力伝送するとき、コイルの銅損とフェライトコアの鉄損で必ず発熱する。Qi規格15Wでも伝送効率80%なら 3Wの損失 が丸ごと熱になる。スマホの薄い筐体にこの熱が閉じ込められるとバッテリー温度が45°Cを超えて寿命が縮む。
たった3Wで問題になるんですか？ EVの充電だともっと大変そう…
その通り。EV用の11kW無線充電（SAE J2954 WPT3クラス）では効率90%でも 1.1kWの損失 だ。送電パッドの表面温度を60°C以下に保つ冷却設計が必須で、電磁-熱連成FEMで温度分布を予測しないと、試作を何回も繰り返すことになる。
ワイヤレス電力伝送（WPT）の基礎
そもそもワイヤレス充電ってどういう原理なんですか？ 電磁誘導ですよね？
ざっくり言うと、送電コイル（Tx）に高周波電流を流して交番磁界を発生させ、その磁束を受電コイル（Rx）で拾って起電力を得る仕組みだ。変圧器の1次巻線・2次巻線を空間的に分離したイメージだね。
WPTの方式は主に3つある：
- 電磁誘導方式 — Qi規格のスマホ充電。コイル間距離は数mm〜1cm。周波数は100〜200kHz。
- 磁界共鳴方式 — EV充電やAirFuel。コイル間距離5〜30cm。共振回路で効率向上。85kHz（SAE J2954規格）。
- マイクロ波方式 — 遠距離給電。GHz帯。効率は低いが距離を取れる。IoTセンサー向け。
熱設計で最も需要が大きいのは電磁誘導方式と磁界共鳴方式の2つだ。
結合係数と伝送効率
「結合係数」ってよく聞くんですけど、具体的にどう効率に関わるんですか？
結合係数 $k$ は送受電コイル間の磁束結合の度合いを表す。相互インダクタンス $M$ と各コイルの自己インダクタンス $L_1, L_2$ を使って：
$k=1$ は全磁束が相手コイルを貫く理想状態（変圧器の鉄心内）。実際のワイヤレス充電では $k=0.1\sim0.6$ 程度。コイル間距離が広がったり位置がずれると $k$ が下がる。
伝送効率 $\eta$ は、各コイルの品質係数 $Q_1 = \omega L_1 / R_1$, $Q_2 = \omega L_2 / R_2$ を使って近似すると：
なるほど、$k$ が小さいと効率が落ちて、その分の損失が全部熱になるわけですね。
そういうこと。全損失は入力電力に対して：
例えばEV用11kWで $\eta=0.90$ なら $P_{\text{loss}} = 1.1\,\text{kW}$。$k$ が0.3から0.15に半減すると $\eta$ は0.9から0.7程度まで落ちて $P_{\text{loss}} = 3.3\,\text{kW}$ に跳ね上がる。だから位置ずれの許容範囲と熱設計はセットで検討しなければならない。
損失の内訳と発熱メカニズム
損失3Wとか1.1kWって言いますけど、具体的にどこが発熱してるんですか？
損失は大きく3つに分解できる。これが発熱源の内訳になる：
1. 銅損（コイル導体のジュール加熱）
コイル導線に流れる電流の $I^2 R$ 損失。ただしWPT周波数（85〜200kHz）では表皮効果と近接効果で実効抵抗がDC抵抗の数倍〜数十倍になる：
ここで $F_r$ は交流抵抗係数で、表皮深さ $\delta = \sqrt{2\rho / (\omega \mu)}$ に依存する。100kHzの銅では $\delta \approx 0.21\,\text{mm}$ だから、単線の丸線だと中心部に電流が流れず効率が悪い。リッツ線（細い素線を撚った構造）を使うのはこのためだ。
2. 鉄損（フェライトコアの磁気損失）
フェライトコアのヒステリシス損と渦電流損の合計。Steinmetz方程式で近似する：
$C_m$, $\alpha$, $\beta$ はフェライト材のデータシートから得る定数（MnZn系フェライトで典型的に $\alpha \approx 1.5$, $\beta \approx 2.5$）。$\hat{B}$ はコア内の最大磁束密度、$V_{\text{core}}$ はコア体積。
見落とされがちだが、フェライトの透磁率と飽和磁束密度は温度依存性が非常に強い。温度が上がるとキュリー温度に近づいて透磁率が急落し、磁束がコアの外に漏れて周囲の金属に渦電流が誘導される——非線形なフィードバックループだ。
3. シールド・筐体の渦電流損
コイル背面のアルミシールドやスマホ筐体の金属部品に漏れ磁束が渦電流を誘導する：
なるほど、発熱源が3つもあって、しかもフェライトは温度で特性が変わるから非線形になるんですね。だからシミュレーションが必要なんだ。
実務での損失比率の目安を表にまとめておこう：
| 損失源 | Qi 15W（スマホ） | WPT3 11kW（EV） | 主な影響因子 |
|---|---|---|---|
| 銅損（Tx+Rx） | 40〜50% | 25〜35% | リッツ線構造、周波数、電流 |
| フェライト鉄損 | 20〜30% | 30〜40% | コア材質、磁束密度、温度 |
| シールド渦電流損 | 15〜25% | 15〜20% | シールド材質・厚さ、距離 |
| 回路損失（インバータ等） | 10〜15% | 10〜20% | MOSFET/GaN素子、スイッチング周波数 |
支配方程式 — 電磁場と熱の連成
電磁場と温度場をどうやって連成させるんですか？ 支配方程式を教えてください。
電磁場はMaxwell方程式、温度場はフーリエの熱伝導方程式で支配される。WPTの場合、準静的な磁場の周波数領域解析（磁気ベクトルポテンシャル $\mathbf{A}$ を未知量とする）と、定常 or 過渡の熱伝導方程式を連成させる。
電磁場の支配方程式（A-V定式化、周波数領域）：
ここで $\mu$ は透磁率（フェライトは温度依存）、$\sigma$ は導電率（銅の $\sigma$ も温度で変化：$\sigma(T) = \sigma_0 / [1+\alpha_R(T-T_0)]$）、$\mathbf{J}_s$ はソース電流密度。
熱伝導の支配方程式：
$Q_{\text{em}}(T)$ が電磁解析から得られる体積発熱率で、銅損・鉄損・渦電流損の合計を空間分布として与える。この $Q_{\text{em}}$ は温度依存の電気抵抗率や透磁率を通じて $T$ の関数になるから、双方向連成（強連成）が本質的に必要なんだ。
電磁場が温度を変え、温度が電気特性を変え、それがまた電磁場を変える…ぐるぐる回るわけですね。
その通り。ただし実務では「温度変化に比べて電磁場の応答が十分速い」ことを利用して、弱連成（片方向反復）で解くことが多い。後で詳しく説明する。
FOD（金属異物）の渦電流発熱
ワイヤレス充電器の上にコインとか置いたら危ないって聞いたんですけど、そんなに問題なんですか？
これはFOD（Foreign Object Detection/Debris）問題と呼ばれていて、安全上非常に重要な課題だ。コイン・クリップ・アルミ箔などの金属片がTxとRxの間に挟まると、交番磁界で大きな渦電流が誘導される。
たった1枚の1円玉（アルミ、直径20mm、厚さ1.5mm）でも、磁束集中部に置かれると渦電流損で 数Wの局所発熱 が生じ、金属片自体が100°C以上に達することがある。15Wクラスでも表面に薄い布1枚しかない充電パッドでは火災リスクになる。
え、そんなに熱くなるんですか！ CAEでFODの発熱もシミュレーションできるんですか？
もちろんできる。FODの渦電流損を解析するには、異物の形状・材質・位置をパラメータとして電磁FEMで損失分布を計算し、その結果を熱FEMに入力する。異物のサイズが表皮深さ $\delta$ と同程度のときに最大損失になるから、メッシュは $\delta/3$ 以下の要素サイズが必要だ。Qi規格ではFOD検出アルゴリズムの検証にもこの解析結果を使う。
電磁-熱連成の物理的意味
- ジュール発熱フィードバック：導体の電気抵抗率は温度上昇に伴い増加する（銅の温度係数 $\alpha_R \approx 0.0039\,\text{K}^{-1}$）。抵抗増加→損失増加→さらに温度上昇、という正のフィードバックが存在する。通常のWPT条件では発散せず安定点に収束するが、極端な過負荷では熱暴走のリスクがある。
- フェライトの温度依存性：MnZnフェライトのキュリー温度は200〜300°C。100°C付近で既に透磁率が大幅に低下し、結合係数 $k$ の変化→伝送効率低下→損失増大のカスケードが生じる。Steinmetzパラメータ $C_m$, $\alpha$, $\beta$ 自体も温度依存であり、修正Steinmetzモデル（iGSE法）を用いることが望ましい。
- 表皮効果の温度依存性：表皮深さ $\delta = \sqrt{2\rho/(\omega\mu)}$ は抵抗率 $\rho$ の平方根に比例する。温度上昇で $\rho$ が増えると $\delta$ も増え、電流がやや深くまで浸透する。交流抵抗係数 $F_r$ が変化するため、銅損の温度依存性は単純な $I^2R$ の線形スケーリングでは捉えきれない。
主要な単位系と次元
| 物理量 | 記号 | SI単位 | 典型値（Qi 15W） |
|---|---|---|---|
| 結合係数 | $k$ | 無次元 | 0.4〜0.6 |
| 品質係数 | $Q$ | 無次元 | 50〜200 |
| 表皮深さ | $\delta$ | m | 0.21 mm (Cu, 100kHz) |
| 体積発熱率 | $Q_{\text{em}}$ | W/m³ | $10^5$〜$10^7$ |
| コア損失密度 | $P_v$ | kW/m³ | 100〜500 (MnZnフェライト, 100kHz, 100mT) |
| 熱伝導率 | $\lambda$ | W/(m·K) | 銅:400, フェライト:3.5〜5, 樹脂:0.2〜0.5 |
スマホが充電中に熱くなるのは「バッテリー」ではなく「コイル」が原因
ユーザーが「ワイヤレス充電中にスマホが熱い」と感じるとき、原因はバッテリー自体の発熱ではなく、受電コイルとフェライトシートの損失が背面ガラスを通じて手に伝わっている場合がほとんどだ。iPhoneのMagSafe（7.5〜15W）では、受電コイルの直下にグラファイトシートとシールドプレートが積層されており、この部分の温度がバッテリーより10〜15°C高くなることがCAE解析で確認されている。Apple・Samsung等のスマホメーカーはこの温度差を利用して「バッテリー保護のためのサーマルスロットリング」を実装している——コイル近傍のサーミスタで温度をモニタリングし、閾値を超えると充電電力を自動的に絞る仕組みだ。
数値解法と実装
連成解析の戦略 — 弱連成 vs 強連成
電磁場と温度場の連成って、どうやって解くのが効率的なんですか？ 全部同時に解く方法と、交互に解く方法があると聞いたんですけど。
WPTの電磁-熱連成では、2つの物理場の時間スケールが大きく異なることがポイントだ：
- 電磁場：周波数85〜200kHz → 周期5〜12μs
- 温度場：熱拡散の特性時間 $\tau = L^2/\alpha_{\text{th}}$ → スマホコイルで数秒〜数十秒
つまり電磁場は温度場より桁違いに速く定常に達する。このため実務では 弱連成（Sequential Coupling） が標準だ：
- 現在温度 $T^n$ での材料特性で電磁場を周波数領域で解く
- 得られた損失分布 $Q_{\text{em}}(T^n)$ を熱発生源として熱伝導方程式を解き、$T^{n+1}$ を得る
- 温度変化が収束判定基準（例: $\|T^{n+1} - T^n\| / \|T^n\| < 10^{-3}$）を満たすまで1-2を繰り返す
強連成が必要になるケースはあるんですか？
フェライトコアがキュリー温度に近い場合（透磁率が急激に変化する温度域）や、大電力WPT（50kW以上）で温度上昇が数百°Cに達する場合は、温度による特性変化が解の安定性に大きく影響するから、強連成（Monolithic）またはサブイタレーション付き弱連成が必要になる。
ただ、Qi規格の民生品やSAE J2954のEV充電では、弱連成で3〜5回の反復でほぼ収束するから、実務上は弱連成で十分だ。
電磁場FEMの定式化
電磁場のFEMって、構造解析のFEMとどこが違うんですか？
構造解析では節点変位がスカラーまたはベクトルの自由度だが、電磁場FEMでは磁気ベクトルポテンシャル $\mathbf{A}$ を未知量とする。3D問題では辺要素（Edge Element, Nedelec要素）を使うのが標準で、節点要素と違ってゲージ条件を自然に満たせる利点がある。
弱形式は：
ここで $\mathbf{w}$ は辺要素の基底関数。離散化すると複素数の連立方程式 $[K_{\text{em}}]\{\mathbf{A}\} = \{\mathbf{F}\}$ が得られ、直接法（MUMPS等）またはKrylov系反復法（BiCGSTAB, GMRES）で解く。WPTのコイルモデルは典型的に数万〜数十万DOFなので、デスクトップPCで十分解ける規模だ。
熱伝導FEMと損失マッピング
電磁解析で求めた損失を、どうやって熱解析モデルに渡すんですか？
電磁FEMの各要素で計算された時間平均損失密度 $Q_{\text{em}}$ [W/m³] を、熱FEMの対応する要素に体積発熱源として入力する。同じメッシュを使えば直接マッピングできるが、電磁場と温度場で最適なメッシュ密度が異なることが多いから、メッシュ非一致の場合は補間が必要だ。
熱伝導のFEM離散化は：
$[C]$ は熱容量マトリクス、$[K_{\text{th}}]$ は熱伝導マトリクス、$\{Q_{\text{bc}}\}$ は対流・輻射などの境界条件からの寄与。定常解析なら左辺第1項がゼロになり、$[K_{\text{th}}]\{T\} = \{Q_{\text{em}}\} + \{Q_{\text{bc}}\}$ を1回解けばよい。
メッシュ戦略 — 表皮深さとコイル形状
WPTのメッシュで特に気をつけるべきポイントはありますか？
WPT特有のメッシュ要件が3つある：
- 表皮深さの解像：導体表面から $\delta$ の深さに電流が集中するから、導体表面に $\delta/3$ 以下の厚さの境界層メッシュを最低3層入れる。100kHzの銅なら $\delta \approx 0.21\,\text{mm}$ だから、約0.07mmの要素が必要。
- コイル断面の形状：リッツ線を素線レベルでモデル化するのは非現実的（数百本の素線）なので、等価導電率モデルで均質化する。この場合、コイル断面に最低4〜6要素。
- エアギャップの解像：送受電コイル間の空気層にも磁束が分布するから、この領域も適切にメッシュを切る。TxとRxの間に最低5層の要素が欲しい。
| 領域 | 要素サイズの目安 | 要素タイプ |
|---|---|---|
| コイル導体表面 | $\delta/3 \approx 0.07\,\text{mm}$ | プリズム/六面体（境界層） |
| フェライトコア | 0.5〜2 mm | 四面体2次 |
| エアギャップ | 1〜3 mm（最低5層） | 四面体/プリズム |
| 遠方空気領域 | 5〜20 mm | 四面体1次 |
| FOD金属片 | $\delta_{\text{FOD}}/3$ | 六面体（境界層必須） |
遠方の空気領域はどこまでモデル化すればいいんですか？ 無限に広がりますよね。
実用的にはコイル外径の3〜5倍の球状領域を取り、外表面に磁気絶縁境界条件 $\mathbf{n} \times \mathbf{A} = 0$ を設定すれば十分だ。もっと精度が欲しい場合はInfinite Element（無限要素）やPML（完全吸収層）を使うが、WPTの準静的問題ではそこまで必要ないことが多い。
実践ガイド
Qi充電（5〜15W）の熱設計
スマホのQi充電で、具体的にどんな熱設計が行われているんですか？
Qi充電パッドの典型的な層構成はこうなっている（下から）：
- ベースプレート（プラスチック筐体）
- 送電コイル Tx（銅リッツ線、フラットスパイラル）
- フェライトシート（0.3〜1mm厚、MnZn系）— 磁束を上方に集中させる
- アルミシールドプレート（0.1〜0.3mm）— 下方への漏れ磁束を遮蔽
この構成で15W充電すると、定常状態でコイル表面温度が50〜65°C、フェライトシートが55〜70°Cに達する。パッド上面（ユーザーが触れる面）を45°C以下に抑えるのが設計目標だ（IEC 62368-1 安全規格）。
温度を下げるにはどんな対策があるんですか？
CAEで効果を評価できる主な対策を3つ挙げる：
- グラファイトシート追加（面方向熱伝導率 1000〜1500 W/m·K）：局所的なホットスポットを面方向に拡散。コイル直上の温度を5〜10°C下げられる。
- フェライト形状最適化：長方形からE型やI型に変更して磁束密度を均一化し、鉄損ホットスポットを低減。
- コイル巻数・線径の最適化：同じインダクタンスでも巻数と線径の組合せで銅損が大きく変わる。電磁-熱連成のパラメトリックスタディで最適点を探索。
EV用WPT（3.3〜11kW）の冷却設計
EVの無線充電はスマホとは桁が違いますよね。冷却はどうするんですか？
SAE J2954規格ではWPT1（3.7kW）からWPT3（11kW）まで定義されていて、送電パッド（GA: Ground Assembly）は地面に設置、受電パッド（VA: Vehicle Assembly）は車両底面に取り付ける。
冷却設計の課題は、GAが 地面設置で自然対流しかない 環境に置かれること。11kWで効率90%なら1.1kWの損失が密閉筐体内にこもる。
| 冷却方式 | 放熱能力 | 適用クラス | コスト |
|---|---|---|---|
| 自然対流のみ | 〜200 W | WPT1（3.7kW）の一部 | 低 |
| アルミ放熱板＋自然対流 | 〜500 W | WPT1〜WPT2 | 中 |
| 強制空冷（ファン内蔵） | 〜1.5 kW | WPT2〜WPT3 | 中 |
| 液冷（冷却水循環） | 〜5 kW | WPT3以上 | 高 |
シミュレーションでは何を予測するんですか？
EV用WPTの電磁-熱連成FEMで予測する主な項目は：
- フェライトコアの最高温度（キュリー温度に対するマージン確認）
- GA表面温度（60°C以下：SAE J2954要求）
- コイル絶縁被覆の最高温度（絶縁階級に対する裕度確認）
- 位置ずれ（X/Y方向±75mm、Z方向100〜250mm）での損失・温度変化
- FODがある場合の局所温度上昇
位置ずれと温度のマップを作成するには、パラメトリックスタディで数十〜数百ケースを回す必要がある。これが実務で計算コストが課題になる理由だ。
解析ワークフロー
実際の解析手順を最初から最後まで教えてください。
WPT電磁-熱連成解析の標準ワークフローはこうだ：
- CADモデル作成：コイル（スパイラル形状）、フェライトコア、シールド、筐体、FOD（必要に応じて）。2D軸対称で足りる場合は計算コストが激減する。
- 材料特性入力：銅の $\sigma(T)$, フェライトの $\mu_r(B,T)$ とSteinmetzパラメータ, 樹脂の $\lambda$, $c_p$ 等。温度依存テーブルが重要。
- 電磁場解析：周波数領域でAC解析。ソース電流（電流駆動）またはポート入力（電圧駆動）。$\mathbf{A}$ を求め、各領域の損失密度を計算。
- 損失マッピング：電磁メッシュから熱メッシュへ体積発熱源を転写。
- 熱解析：定常 or 過渡。外面に対流 $h$・輻射 $\varepsilon$ の境界条件。温度分布を求める。
- 収束判定：温度変化が閾値以下なら終了、そうでなければ温度依存の材料特性を更新してStep 3に戻る。
- 後処理：最高温度、温度勾配、表面温度分布の確認。安全規格への適合判定。
境界条件の設定指針
熱解析の境界条件で間違えやすいポイントはありますか？
WPT熱解析で特に注意すべき境界条件を3つ挙げよう：
- スマホ背面の接触熱抵抗：Qi充電パッドとスマホの間にはケース・空気層があり、実際の接触状態で等価熱伝達係数が大きく変わる（50〜500 W/m²K）。この値の不確実性が結果に最も影響する。実測との整合が必須。
- GA底面（地面側）：アスファルトやコンクリートの熱物性と、夏場の路面温度（60°C以上になることがある）を考慮。地面を半無限体としてモデル化するか、等価熱抵抗で近似する。
- 自然対流の熱伝達係数：密閉筐体の自然対流は $h = 5\sim10\,\text{W/m²K}$ 程度。向き（上面/下面/側面）で異なる。強制空冷の場合は $h = 30\sim100\,\text{W/m²K}$。
テスラとBMWのワイヤレス充電パッド開発裏話
EV用ワイヤレス充電の実用化では、BMW 530eが2018年に世界初のOEM純正オプションとして3.2kWのWPTシステムを市場投入した。地面設置のGAパッドはフェライトコアの発熱がボトルネックだったが、E型フェライトの形状最適化とアルミ放熱フィンの追加でパッド表面温度を55°C以下に収めた。この設計にはCOMSOLの電磁-熱連成解析が全面的に活用されたという。一方テスラは自社でWPTの開発を進め、11kWクラスの強制空冷パッドのプロトタイプでは冷却ファンの騒音が課題になっている。CAEによるファン騒音と冷却性能のトレードオフ最適化は、今後のWPT製品開発の重要テーマだ。
WPT熱解析で最もやりがちなミス
「フェライトの温度依存性を無視して解析した」——これが実務で最も多いミスだ。室温（25°C）で測定したSteinmetzパラメータをそのまま使って熱解析すると、実際の動作温度（80〜120°C）でのコア損失を20〜50%過小評価する。特にMnZnフェライトは80°C付近で損失が最小になり、それを超えると急激に増加する「バスタブカーブ」を示す。温度依存のSteinmetzパラメータ（またはiGSE: improved Generalized Steinmetz Equation）を使わないと、設計が甘くなり試作で「予想より10°C高い」という結果になる。
ソフトウェア比較
主要ツールの機能比較
WPTの電磁-熱連成をやるにはどんなソフトがあるんですか？ 選び方を教えてください。
WPT熱解析で使われる主要ツールを比較しよう：
| ツール | 電磁-熱連成 | WPTへの適性 | 特徴 |
|---|---|---|---|
| COMSOL (AC/DC + Heat Transfer) | 一体型連成 | 非常に高い | 1つのモデルで電磁＋熱。パラメトリック最適化との統合が容易 |
| Ansys Maxwell + Icepak | ツール間連携 | 高い | Maxwell→Icepakへの損失転送。Workbenchで自動化可能 |
| JMAG-Designer | 内蔵熱ソルバー | 高い | 電気機器設計に特化。コイル・コアの損失計算が得意 |
| Ansys Maxwell + Mechanical | ツール間連携 | 中程度 | 構造連成まで拡張可能（熱応力解析） |
| CST Studio Suite | 内蔵連成 | 中程度 | 高周波寄りだがLF領域も対応。3D EM+Thermal |
| FEMM（OSS）+ Python | スクリプト連携 | 2Dのみ | 無料。2D問題の予備検討に最適。自作スクリプトで熱解析 |
COMSOLとAnsysの違いって何ですか？ どっちを使えばいいか迷います。
ざっくり言うと：
- COMSOL：1つのGUI内で電磁場と温度場を自由に連成できる「マルチフィジクスの専門店」。研究開発やプロトタイピングに向く。パラメトリックスイープが直感的。ライセンスはモジュール単位（AC/DC Module + Heat Transfer Module で年間200万円程度）。
- Ansys Maxwell + Icepak：電磁場専用ソルバー（Maxwell）と電子機器冷却専用ソルバー（Icepak）の連携。それぞれの分野で成熟した機能を持つ。大企業の量産設計で実績が多い。Workbenchでの自動化・最適化が強い。ライセンスは高額だが並列計算に強い。
学生・個人研究者なら COMSOL か FEMM+Python、企業の量産設計部門なら Ansys が実務的な選択肢だ。
選定の指針
最終的にはこの3つの問いで選べばいい：
- 「2Dで足りるか3Dが必要か」：円形コイルの軸対称問題なら2D（FEMM無料）で90%の情報が得られる。位置ずれ・矩形コイル・FOD解析は3D必須。
- 「パラメトリック最適化をするか」：コイル巻数・フェライト形状・エアギャップを振って最適設計するなら、最適化ツール統合が強い COMSOL か Ansys Workbench。
- 「社内の既存ツールは何か」：電磁場チームがMaxwell、冷却チームがIcepakを使っている会社なら、新規ツール導入コストを考えてAnsys一択になることもある。
先端技術
走行中給電（Dynamic WPT）
走りながら充電できる「走行中給電」って聞いたことがあるんですけど、熱設計はどうなるんですか？
Dynamic WPT（DWPT）は道路に送電コイルを埋設し、走行中のEVに連続給電するシステムだ。韓国KAISTのOLEVプロジェクトやスウェーデンの実証路線で進んでいる。
熱設計上の新しい課題は：
- 間欠動作：車両が通過する瞬間だけ各コイルが通電される（数秒間）。定常状態に達しないため、過渡熱解析が必須。
- 路面温度の制約：アスファルト舗装の軟化点は約60°C。埋設コイルの周囲温度がこれを超えると路面が劣化する。
- 多数コイルの相互干渉：数メートル間隔で数十〜数百個のコイルが並ぶため、隣接コイルの漏れ磁束が累積し、フェライト温度が上がりやすい。
DWPTの熱解析は、単一コイルの電磁-熱連成に加えて、土壌の熱伝導や日射・降雨の影響まで考慮する必要があり、解析規模が格段に大きくなる。
機械学習サロゲートモデル
パラメトリックスタディで数百ケースも回すのは大変ですよね。AIで効率化できないんですか？
まさにそこが最近のホットトピックだ。電磁-熱連成FEMの計算を数十〜数百ケース実行してデータを蓄積し、そのデータで機械学習モデル（サロゲートモデル）を構築する。
- 入力：コイル外径、巻数、線径、エアギャップ、周波数、電流振幅、フェライト形状パラメータ
- 出力：伝送効率 $\eta$, 最高温度 $T_{\max}$, 損失内訳
Gaussian Process（Kriging）やニューラルネットワークで構築したサロゲートモデルを使えば、FEMでは1ケース数分〜数十分かかる計算が ミリ秒 で済む。これを遺伝的アルゴリズムやベイズ最適化と組み合わせて、効率と温度の多目的最適化を行う研究が活発だ。
最近ではPINN（Physics-Informed Neural Network）で物理法則をNNの損失関数に組み込み、少数のFEMデータでも高精度なサロゲートを構築する手法が注目されている。
トラブルシューティング
よくあるエラーと対策
WPTの電磁-熱連成解析で、ハマりやすいトラブルってありますか？
実務でよく遭遇するトラブルと対策をまとめておこう：
| 症状 | 原因 | 対策 |
|---|---|---|
| 損失が実測の2〜3倍 | AC抵抗を考慮せずDC抵抗で銅損を計算 | 表皮効果・近接効果を含む交流抵抗モデルを使用。または導体表面にメッシュ境界層を追加 |
| フェライト温度が異常に高い | Steinmetzパラメータの温度依存性を無視 | iGSE法または実測の温度依存損失テーブルを使用 |
| 連成反復が収束しない | 温度依存性が強すぎて解が振動 | 緩和係数 $\omega = 0.3\sim0.7$ を導入。Aitken緩和も有効 |
| $k$ が実測より低い | 2Dモデルで3D的な磁束漏れを捉えていない | 3Dモデルへ移行。または等価的な漏れインダクタンス補正 |
| 温度分布が対称でない | メッシュの非対称性、または収束不足 | 対称メッシュを強制。連成反復回数を増やす |
| COMSOLで "Failed to find a solution" | 材料特性テーブルの温度範囲外参照 | 材料特性テーブルの温度範囲を解析で想定される範囲まで拡張 |
解析結果が実測と合わないとき、まず何を確認すればいいですか？
WPT解析のデバッグ手順はこうだ：
- エネルギー収支の確認：入力電力 = 出力電力 + 全損失 になっているか。5%以上のずれがあれば、メッシュ不足か境界条件のミス。
- $k$ と $L$ の検証：電磁解析で得られた結合係数と自己インダクタンスが、LCRメーターの実測値と±10%以内で一致するか。ここがずれていると損失計算も全部ずれる。
- 単体損失の切り分け：銅損だけ、鉄損だけ、渦電流損だけを個別に確認して、どこが実測と乖離しているか特定する。
- メッシュ収束性：メッシュを2倍に細かくして結果が5%以上変わるなら、まだ収束していない。特に表皮深さ付近のメッシュ。
そうだね。ポイントをまとめると、(1) 結合係数 $k$ と品質係数 $Q$ で伝送効率が決まり、(2) 効率の裏返しが損失＝発熱、(3) 損失は銅損・鉄損・渦電流損の3つに分解でき、(4) それぞれが温度依存だから電磁-熱連成が必須——ということだ。まずは2D軸対称モデルで手を動かしてみることを勧めるよ。FEMM+Pythonなら無料で始められるからね。
なった
詳しく
報告