弱め界磁制御の電磁界シミュレーション
理論と物理
弱め界磁制御とは何か
弱め界磁って、わざと磁石の力を弱めるんですか？ もったいなくないですか？
一見もったいないように見えるけど、基底速度以上で回すにはインバータの電圧制限に引っかかるんだ。d軸電流を負方向に流して磁石の磁束を打ち消すことで、逆起電力を抑えて高速回転を実現する。EVの高速道路走行には必須の技術だよ。
え、EVってギアで速度を変えるんじゃないんですか？
多くのEVは変速機を持たない、いわゆるシングルギアだ。エンジン車なら5速や6速に入れて高速走行するところを、EVのモータは弱め界磁制御で定格回転数の3〜4倍まで回して対応する。つまり、弱め界磁制御はモータに「仮想的なギアチェンジ」をさせる技術とも言えるね。
なるほど、ギアの代わりに電流の使い方を変えるってことですね！ でも磁石の力を弱めるとトルクは下がりませんか？
そのとおり。弱め界磁領域ではトルクは回転数に反比例して低下する。つまり出力（=トルク×回転数）はほぼ一定の「定出力領域」になるんだ。高速道路で100km/h巡航するくらいなら必要トルクは少ないから、これで十分なんだよ。
弱め界磁制御（Field Weakening Control）は、永久磁石同期モータ（IPMSM/SPMSM）において基底速度を超えた高速動作を実現するための電流制御技術である。d軸に負方向の電流（$i_d < 0$）を通流して永久磁石磁束 $\psi_m$ を電気的に打ち消し、逆起電力を低減することで、インバータの出力電圧制約内での高速回転を可能にする。
dq軸電圧方程式と制約条件
具体的に、どんな数式で制御が制約されるんですか？
IPMSMのdq軸電圧方程式から出発しよう。定常状態（$di/dt = 0$）を仮定すると、こうなる：
ここで $R_s$ は巻線抵抗、$L_d, L_q$ はdq軸インダクタンス、$\omega_e$ は電気角速度、$\psi_m$ は永久磁石磁束鎖交数である。高速域では $R_s$ による電圧降下は小さいため無視すると、電圧制限条件は次のように整理できる：
同時に、インバータの電流定格による電流制限条件が課される：
つまり、電圧と電流の両方に上限があって、その範囲内でモータを動かさないといけないんですね。
そう。そしてここが重要なんだけど、電圧制限式を $i_d$-$i_q$ 平面上で見ると楕円になる。回転数が上がるほど $\omega_e$ が増えて楕円が縮小する。だから高速域では動作可能な電流領域がどんどん狭くなっていくんだ。
電圧制限円と電流制限円
電圧制限条件を $i_d$-$i_q$ 平面で書き直すと、次のようになる：
これは中心 $(-\psi_m/L_d, \, 0)$、長軸・短軸が $V_{\max}/(\omega_e L_d)$, $V_{\max}/(\omega_e L_q)$ の楕円である。IPMSMでは $L_d < L_q$ のためd軸方向が短軸となる。回転数の上昇に伴い楕円は縮小し、電流制限円 $i_d^2 + i_q^2 = I_{\max}^2$ との交差領域が変化する。
楕円の中心が原点じゃなくて $(-\psi_m/L_d, 0)$ にズレてるのはなぜですか？
永久磁石磁束 $\psi_m$ が常にd軸方向に存在するから、その分だけ電圧が「食われて」いるんだ。d軸電流がゼロでも磁石磁束による逆起電力は発生する。だから楕円の中心は $i_d = -\psi_m/L_d$ の位置にシフトする。弱め界磁制御とは、動作点をこの楕円中心に向かって移動させることに他ならないんだよ。
MTPA制御とMTPV制御
MTPAとかMTPVとかよく聞くんですけど、何が違うんですか？
簡単に言うと、速度域で使い分ける2つの最適電流戦略だ。
MTPA（Maximum Torque Per Ampere）： 電流制限円の内側で、与えられた電流振幅に対して最大トルクを発生する電流ベクトル角 $\beta$ を選ぶ戦略。IPMSMのトルク式：
において、リラクタンストルク項 $(L_d - L_q) i_d i_q$ を活用するため、$i_d < 0$ の領域で動作する。MTPA軌跡は $i_d$-$i_q$ 平面上で双曲線状のカーブとなる。
MTPV（Maximum Torque Per Voltage）： 電圧制限楕円上で最大トルクを発生する動作点を選ぶ戦略。高速域で電圧制限に達した後、楕円の縁に沿って動作点を移動させる。MTPV軌跡の条件式は：
じゃあ速度によって、MTPA軌跡 → 電圧制限楕円上 → MTPV軌跡、って動作点が移り変わるんですか？
そのとおり。低速域ではMTPA軌跡上で最高効率運転、基底速度付近で電圧制限に到達してからは楕円上を移動、さらに高速域ではMTPV軌跡に遷移する。この3つの領域を滑らかに切り替えるのが制御の腕の見せ所で、CAEでは各速度・トルクの組合せに対して最適電流点を全部計算して「効率マップ」を作るんだ。
脱磁リスクの物理
d軸電流を大きくしすぎると磁石が壊れると聞いたんですが、本当ですか？
壊れるというより「脱磁」する。ネオジム磁石には減磁曲線（B-H曲線の第2象限）があって、動作点がニーポイント（曲線の折れ点）を超えると不可逆的に磁力が低下する。弱め界磁制御では大きな負のd軸電流を流すから、磁石内部の磁束密度がニーポイント近くまで下がる危険性があるんだ。
温度も関係しますよね？ 夏場とかヤバそう...
まさにそこが一番危ない。ネオジム磁石の保磁力は温度上昇で急激に低下する。20°Cで2000 kA/mあった保磁力が150°Cでは半分以下になることもある。だからCAEでは、最悪条件（最高温度 × 最大d軸電流 × 過渡的な電流オーバーシュート）での磁石動作点を必ず確認するんだ。
脱磁判定は、FEMで求めた磁石内部の磁束密度分布と、温度に対応した減磁曲線を照合して行う。磁石コーナー部（角部）は反磁場の集中により局所的に磁束密度が低下しやすいため、特に注意が必要である。
テスラModel 3のモータ設計が「弱め界磁の教科書」である理由
テスラModel 3のリアモータ（IPMSMタイプ）は、定格約5,000rpmに対して最高回転数が約18,000rpm。つまり基底速度の3.6倍まで弱め界磁で回している。これを可能にしているのは、磁石埋込み構造（V字配置）によるリラクタンストルクの活用と、SiCインバータによる高い電圧利用率だ。設計段階ではFEMによる弱め界磁時の脱磁余裕解析が膨大な回数行われたと推察される。日本の自動車OEMでも同様の設計プロセスが標準化されており、JMAGやAnsys MaxwellでのFEM-回路連成解析が中心的な役割を果たしている。
dq軸電圧方程式の導出過程
- d軸電圧 $v_d = R_s i_d + L_d \frac{di_d}{dt} - \omega_e L_q i_q$：d軸の磁束鎖交数 $\lambda_d = L_d i_d + \psi_m$ の時間微分に抵抗降下を加えたもの。$-\omega_e L_q i_q$ の項はq軸磁束の回転による速度起電力（干渉項）。
- q軸電圧 $v_q = R_s i_q + L_q \frac{di_q}{dt} + \omega_e(L_d i_d + \psi_m)$：q軸の磁束鎖交数 $\lambda_q = L_q i_q$ の時間微分に、d軸磁束回転による逆起電力 $\omega_e(L_d i_d + \psi_m)$ を加えたもの。この逆起電力項が弱め界磁制御の本質的制約。
- トルク式 $T = \frac{3}{2}p(\psi_m i_q + (L_d - L_q)i_d i_q)$：第1項は磁石トルク（永久磁石磁束とq軸電流の相互作用）、第2項はリラクタンストルク（突極性 $L_d \neq L_q$ による磁気的非等方性を利用）。IPMSMでは $L_q > L_d$ なので $i_d < 0$ でリラクタンストルクが正に寄与する。
弱め界磁制御の適用限界と仮定
- 定常状態仮定：$di/dt = 0$ の定常状態で電圧制限楕円を描く。過渡時は楕円が変動するため、瞬時的に電圧超過が発生する可能性がある
- 磁気飽和の影響：高負荷時には $L_d$, $L_q$ が電流に依存して変化する。FEMで得た非線形インダクタンスマップを使わないと誤差が大きくなる
- 温度依存性：$\psi_m$ は磁石温度の上昇で減少し、$R_s$ は巻線温度で増加する。広い温度範囲で効率マップを作る必要がある
- 空間高調波：集中巻モータでは巻線配置に起因する空間高調波磁束が存在し、dq軸モデルの精度が低下する場合がある
- 鉄損の影響：d軸電流を増やすと鉄心の磁束密度変動幅が変わり、鉄損分布が変化する。等価回路モデルでは精度不足な場合がある
主要パラメータの典型値（EV用IPMSM）
| パラメータ | 典型値 | 備考 |
|---|---|---|
| 永久磁石磁束鎖交数 $\psi_m$ | 0.05〜0.15 Wb | 温度上昇で約-0.1%/°C低下 |
| d軸インダクタンス $L_d$ | 0.1〜0.5 mH | 電流依存（飽和で低下） |
| q軸インダクタンス $L_q$ | 0.3〜1.5 mH | IPMSMでは $L_q > L_d$ |
| 突極比 $L_q/L_d$ | 2〜8 | 大きいほどリラクタンストルク大 |
| 基底速度 | 3,000〜6,000 rpm | 電圧利用率100%到達点 |
| 最高速度 / 基底速度比 | 3〜5 | 弱め界磁定出力領域の幅 |
| DC母線電圧 | 300〜800 V | $V_{\max} = V_{dc}/\sqrt{3}$（SVPWM） |
数値解法と実装
FEMモータモデリングの基礎
弱め界磁制御をFEMでシミュレーションするって、具体的にはどういう計算をするんですか？
大きく分けて2つのアプローチがある。1つは2D断面FEMで磁場解析を行い、各回転角度位置でのトルク・磁束鎖交数・損失を計算する方法。もう1つは、FEMと駆動回路（インバータ）を連成させて、弱め界磁制御アルゴリズムそのものをシミュレートする方法だ。
モータの電磁界解析では、ベクトルポテンシャル $\mathbf{A}$ を用いた準静的磁場方程式を解く：
ここで $\nu$ は磁気抵抗率（$\nu = 1/\mu$）、$\mathbf{J}_s$ はコイル電流密度、$\sigma \partial \mathbf{A}/\partial t$ は渦電流項、$\mathbf{M}_r$ は永久磁石の残留磁化ベクトルである。
2D解析で大丈夫なんですか？ モータって3次元の物体ですよね？
鋭いね。積層珪素鋼板を使うモータでは、軸方向の電流が支配的だから2D断面解析で十分な精度が得られることが多い。スキュー（捩り）やエンドコイルの影響は補正係数で扱う。ただし3D特有の端部効果、例えばスキューによるトルクリプル低減やエンドコイル漏れインダクタンスを正確に評価したいなら3D解析が必要だ。
回路連成解析
FEMと回路を連成させるって、どうやるんですか？
FEMの各コイル領域から磁束鎖交数 $\psi$ と逆起電力を抽出し、外部回路のKirchhoffの法則と連立させる。コイル電流が磁場に影響し、磁場がコイルの逆起電力に影響するから、毎時間ステップで反復的に解く必要があるんだ。
FEM-回路連成の基本方程式系は以下のようになる：
ここで $K_{ff}$ はFEMの磁場剛性マトリクス、$K_{cc}$ は回路のインピーダンスマトリクス、$K_{fc}$, $K_{cf}$ はFEMと回路の結合項、$A$ は節点ベクトルポテンシャル、$I$ はコイル電流、$V_s$ はインバータ出力電圧である。
インバータのPWMスイッチングも含めて解くんですか？ それだとものすごい計算量になりそう...
目的次第だね。弱め界磁領域の平均的なトルクや効率を見たいだけなら、基本波電圧源でOK。トルクリプルやEMC解析をしたいなら、PWMキャリアの時間刻み（数μs）で解く必要がある。実務では多くの場合、まず基本波解析で効率マップを作り、気になる動作点だけPWM連成で詳細検討するよ。
損失計算（鉄損・銅損・磁石渦電流損）
弱め界磁領域での効率を正しく評価するには、以下の損失成分を正確に計算する必要がある：
| 損失成分 | 計算式 | 弱め界磁での変化 |
|---|---|---|
| 銅損 $P_{cu}$ | $\frac{3}{2}R_s(i_d^2 + i_q^2)$ | $|i_d|$ 増加で増大（最大の損失因子） |
| ヒステリシス鉄損 $P_h$ | $k_h f B_m^{\alpha}$ | 周波数 $f$ 増加で増大、$B_m$ 低下で減少 |
| 渦電流鉄損 $P_e$ | $k_e f^2 B_m^2$ | $f^2$ に比例し急増する主要損失 |
| 磁石渦電流損 | FEM時間領域解析で直接計算 | 高調波磁束増加で増大 |
| 機械損 | ベアリング摩擦 + 風損 | 回転数の2〜3乗に比例 |
鉄損の計算式で $B_m$ は何を使うんですか？ 磁石の $B_r$ とは違いますよね？
$B_m$ は各要素での磁束密度の振幅だよ。FEM解析で1電気角周期の過渡解析をして、各有限要素の磁束密度波形を取得する。そこからFFTで基本波成分と高調波成分を分離して、それぞれの周波数で鉄損を計算する。これを全要素について積算するのがFEMベースの鉄損計算だ。
効率マップの構築手法
効率マップって、どうやって作るんですか？ 全部の速度×トルクの組合せでFEM解析するんですか？
それだと計算量が膨大になるから、実務ではFEMを使って「dq軸の電流をパラメトリックに振った磁場解析」を先に実行する。例えば $i_d$ を -200A〜0A、$i_q$ を 0A〜300Aまで20点ずつ振れば400ケースの静磁場解析になる。各ケースでトルク、磁束鎖交数、鉄損が得られるから、それをルックアップテーブル（LUT）にして、任意の速度・トルク条件で最適電流を探索するんだ。
効率マップ構築の手順は以下のとおりである：
- FEMパラメトリック解析： $(i_d, i_q)$ の格子点で静磁場解析 → トルク $T(i_d, i_q)$, $\psi_d(i_d, i_q)$, $\psi_q(i_d, i_q)$, 鉄損 $P_{fe}(i_d, i_q)$ のマップを作成
- 最適電流探索： 各速度・トルク指令に対して、電圧制限・電流制限を満たしつつ損失最小の $(i_d^*, i_q^*)$ を探索
- 効率計算： $\eta = P_{out} / (P_{out} + P_{cu} + P_{fe} + P_{mech})$ をプロット
非線形収束と磁気飽和の取り扱い
珪素鋼板のB-Hカーブって非線形ですよね。FEMでどう扱うんですか？
ニュートン・ラフソン法で非線形反復を行う。各反復で、現在の磁束密度からB-Hカーブを参照して透磁率を更新し、剛性マトリクスを再構築して解き直す。弱め界磁領域では、d軸電流によって磁石近傍の鉄心が局所的に脱飽和する場合があり、収束が改善することもあるし、逆にq軸の高電流で歯部が深い飽和に入ると収束が悪化することもある。
非線形収束の判定基準は残差ノルムの相対値が一般的で、$||R||/||R_0|| < 10^{-4}$ 程度が推奨される。磁気飽和が深い場合は、ライン探索や適応的緩和を併用してロバスト性を確保する。
FEM解析とLPM（集中定数モデル）の使い分け
「FEMで全部やるべき」と思いがちだが、効率マップの構築に1ケースあたり数時間かかるFEM解析を数百ケース回すのは現実的でない。実務ではLPM（等価回路モデル）で概略設計を行い、最適解の候補を絞ってからFEMで精密検証する2段階アプローチが標準だ。Motor-CADのようなLPMツールは1つの動作点を秒単位で計算できるため、パラメトリック探索に圧倒的に有利である。
実践ガイド
弱め界磁解析のワークフロー
実際に弱め界磁のシミュレーションをするとき、どういう手順で進めるんですか？
典型的なワークフローを順を追って説明するね。
Step 1: 幾何モデル作成
- ステータ・ロータ・永久磁石・空気領域（ギャップ含む）の断面形状を作成
- 対称性の活用：8極48スロットなら1/8モデル（1極分）で解析可能
- スライディングメッシュ用の回転界面を定義
Step 2: 材料設定
- 珪素鋼板：非線形B-Hカーブ + 鉄損係数（Steinmetz式またはBertotti式）
- 永久磁石：残留磁化 $B_r$、保磁力 $H_{cJ}$、温度係数、減磁曲線（第2象限）
- 導体：銅の導電率（温度依存）
Step 3: パラメトリック励磁解析
- d軸電流 $i_d$ と q軸電流 $i_q$ の格子点（例：20×20 = 400点）を定義
- 各点で1電気角周期のトランジェント解析を実行（回転角度刻み：電気角1〜2°）
- トルク平均値、磁束鎖交数（$\psi_d$, $\psi_q$）、鉄損を抽出
Step 4: 効率マップ構築
- LUTデータから最適電流探索（損失最小化）を速度×トルクの各点で実施
- MTPA/弱め界磁/MTPV領域の自動判別
- 等効率線（例：95%、90%、85%）のコンターマップを作成
Step 5: 脱磁リスク評価
- 最大d軸電流 × 最高温度条件でのFEM解析
- 磁石各部の磁束密度分布と減磁曲線の照合
- 脱磁余裕（ニーポイントまでの距離）の定量評価
モータ電磁界解析のメッシュ戦略
メッシュで特に気をつけるところはどこですか？
モータ特有の注意点がいくつかある。
- エアギャップ： 最低3〜5層の要素を配置。ギャップ長0.5mmに対して0.1mm以下の要素サイズが必要。トルク精度に直結する最も重要な領域
- 磁石エッジ： 局所的な減磁判定のため、磁石コーナー部は細かくメッシュを切る
- スロット開口部： 磁束密度の変動が大きい領域。コギングトルク精度に影響
- ロータブリッジ： 磁気飽和が深い薄肉部。最低2〜3層の要素が必要
- 回転界面： スライディングメッシュの分割数は、精度と回転角度ステップの整合性に注意
| 領域 | 推奨要素サイズ | 要素タイプ | 理由 |
|---|---|---|---|
| エアギャップ | 0.05〜0.15 mm | 三角形1次 | トルク精度・マクスウェル応力テンソル |
| 磁石内部 | 0.3〜1.0 mm | 三角形1次 | 脱磁評価用に十分な分解能 |
| 歯部・ヨーク | 0.5〜2.0 mm | 三角形1次 | 飽和・鉄損計算に必要 |
| 外部空気領域 | 5〜20 mm | 三角形1次 | 漏れ磁束の減衰を表現 |
境界条件と対称性の活用
8極48スロットのモータで1/8モデルを使うとき、境界条件はどうするんですか？
周方向の両端に周期境界条件（反周期境界条件）を適用する。隣り合う極の磁極が反転するから、$A(0°) = -A(45°)$ のような反周期条件を使うんだ。外周境界は $A = 0$（ディリクレ境界）で磁束が外に漏れないようにする。外部空気領域はステータ外径の1.5〜2倍程度取れば十分。
脱磁余裕チェックの実務手順
脱磁チェックの実務的な手順を教えてください。判定基準は何ですか？
手順はこうだ。まず最悪条件を定義する：
- 磁石温度：最高使用温度（例：150°C〜180°C）
- d軸電流：最大（短絡故障時を含む）
- q軸電流：ゼロ（d軸電流の効果を最も厳しく評価）
この条件でFEM解析を回して、磁石内部の最小磁束密度 $B_{\min}$ を確認する。そして対応温度の減磁曲線のニーポイント $B_{knee}$ と比較して：
- $B_{\min} > B_{knee} + \Delta B_{margin}$ であればOK
- 余裕 $\Delta B_{margin}$ は一般に 0.05〜0.1T を確保する
- $B_{\min}$ が $B_{knee}$ を下回る場合は、磁石グレードの変更（例：NからSHへ）、d軸電流の上限制限、磁石形状の変更（フラックスバリア追加）で対策する
「弱め磁束で脱磁した」——EV開発現場の実失敗事例
ある国産EV開発プロジェクトで弱め磁束制御の電流制限値の設定ミスにより、過大なd軸電流が流れてネオジム磁石が部分脱磁した事例がある。脱磁はトルク低下・効率悪化として現れ、製品回収寸前まで原因特定に時間がかかった。根本原因は「弱め磁束領域での磁石動作点がシミュレーションに含まれていなかった」ことだ。現在の標準プロセスでは弱め磁束時の磁石動作点をFEM上で全温度範囲にわたり確認し、減磁曲線との余裕（demagnetization margin）を設計仕様として管理する。
初心者が陥りやすい落とし穴：エアギャップメッシュ不足
「メッシュ数を増やしてもトルクが収束しない」という相談を受けると、大抵はエアギャップのメッシュが粗い。ステータ歯部やヨークをいくら細かく切っても、トルクを決めているのはギャップの磁束密度分布だ。例えると、プールの水量を知りたいのに排水口の計測器だけ精密にしているようなもの。ギャップに3〜5層の要素を入れるだけでトルクの計算精度が劇的に改善する。
ソフトウェア比較
弱め界磁解析対応ツール比較
弱め界磁のシミュレーションに使えるソフトって、どんなものがありますか？
モータ電磁界解析の主要ツールを弱め界磁解析の観点で比較してみよう。
| 機能 | JMAG | Ansys Maxwell | Motor-CAD | COMSOL |
|---|---|---|---|---|
| 2D FEM過渡解析 | ◎ | ◎ | -（LPM） | ○ |
| 3D FEM解析 | ◎ | ◎ | - | ○ |
| FEM-回路連成 | ◎ | ◎（Twin Builder） | - | ○ |
| 効率マップ自動生成 | ◎（標準搭載） | ○（スクリプト要） | ◎（主機能） | △ |
| 脱磁解析 | ◎（温度連成） | ○ | △（簡易） | ○ |
| MTPA/MTPV自動計算 | ◎ | ○ | ◎ | △ |
| 磁石温度依存 | ◎ | ○ | ◎ | ○ |
| 計算速度（効率マップ） | 中（数時間〜） | 中（数時間〜） | ◎（数分） | 低 |
| 日本語サポート | ◎（国産） | ○ | ○ | ○ |
JMAGでの弱め界磁解析
日本の自動車メーカーだとJMAGが多いと聞きますが、弱め界磁解析の機能は充実してますか？
JMAGはモータ設計に特化しているだけあって、弱め界磁関連の機能が非常に充実しているよ。JMAG-Designerの「効率マップ」機能を使えば、dq軸電流のパラメトリック解析から効率マップの生成まで半自動で実行できる。特に脱磁解析は温度連成で行えるのが強みで、磁石の温度分布を考慮した上で各部位の減磁余裕を評価できる。
JMAGの弱め界磁解析のワークフロー：
- テンプレートモデル作成： モータテンプレート機能で幾何形状を定義
- dq軸パラメトリック解析： 「特性解析」モジュールで $(i_d, i_q)$ 格子の自動掃引
- 効率マップ生成： 損失計算結果を用いた最適電流探索と効率コンター作成
- 脱磁チェック： 「減磁耐力解析」で磁石動作点と減磁曲線の余裕を自動評価
Ansys Maxwellでの弱め界磁解析
Ansys Maxwellだとどんな流れになりますか？
MaxwellはRMxprt（解析設計）→ Maxwell 2D/3D（FEM）→ Twin Builder（システムシミュレーション）→ Ansys Motor-CAD（効率マップ）という連携が強みだ。RMxprtで概略設計したモデルをMaxwellに自動エクスポートしてFEM検証できるし、Twin BuilderでインバータのPWM制御まで含めた時間領域シミュレーションが可能。ただし効率マップの自動生成はJMAGほど簡単ではなく、Pythonスクリプト（PyAEDT）での自動化が実務では必要になることが多い。
Motor-CADによるシステムレベル設計
Motor-CADはFEMじゃないんですよね？ それでも弱め界磁の解析はできるんですか？
Motor-CADは集中定数モデル（LPM）ベースだけど、効率マップの高速計算という点では圧倒的に優れている。1つの速度×トルク点の計算が秒単位で終わるから、冷却条件を変えながら効率マップを即座に更新するとか、100パターンのロータ形状をスクリーニングするとか、概念設計段階では無敵の武器だよ。ただし磁石コーナー部の局所脱磁評価はFEMでないと精度が出ないから、最終検証はJMAGやMaxwellに任せるのが一般的だね。
JMAG vs Motor-CAD——「精密さ」と「速さ」の二刀流
弱め磁束領域を含むモータ解析ツールとして国内ではJMAG（JSOL製、FEM）が定評がある。弱め磁束時の鉄損・銅損・脱磁余裕を温度依存で計算するワークフローが充実しており、自動車OEMと一次サプライヤで広く採用されている。Motor-CAD（Ansys製）は熱-電気-機械の連成解析を迅速に行えるルミテッドパラメータモデルが強みで、システムレベルの効率マップ計算速度が圧倒的に速い。詳細設計はJMAG、概念設計はMotor-CADという使い分けが多く見られる。近年はJMAGの解析結果をMotor-CADにインポートして、FEM精度のLUTとLPM速度を組合せるハイブリッドワークフローも一般化しつつある。
ツール選定のポイント
- 概念設計・パラメトリック探索 → Motor-CAD： 数千パターンの形状・制御条件を高速にスクリーニング。冷却系との連成も得意。
- 詳細設計・精密損失計算 → JMAG / Ansys Maxwell： 非線形FEMによる磁気飽和・鉄損・脱磁の高精度評価。日本のOEMではJMAGの使用率が高い。
- マルチフィジックス連成 → COMSOL： 電磁-熱-構造の一貫した連成解析が必要な研究開発向け。ただしモータ専用GUIは他ツールに劣る。
- 制御アルゴリズム検証 → MATLAB/Simulink + FEMエクスポート： FEMで作成したLUTをSimulinkの制御モデルに組み込み、HILSで実装検証する流れが主流。
トラブルシューティング
弱め界磁シミュレーションの典型的トラブル
先生、弱め界磁のシミュレーションでハマりやすいポイントってどこですか？
いい質問だね。弱め界磁特有の落とし穴をまとめよう。
| 症状 | 原因 | 対策 |
|---|---|---|
| 効率マップの高速域でトルクがゼロにならない | 電圧制限条件の実装ミス。$V_{\max}$ にSVPWMの電圧利用率 $1/\sqrt{3}$ を考慮していない | $V_{\max} = V_{dc}/\sqrt{3}$ を正しく設定 |
| FEMと等価回路のトルクが一致しない | FEM側の $L_d$, $L_q$ が電流依存で変化するのに等価回路では定数としている | FEMで得た $L_d(i_d, i_q)$, $L_q(i_d, i_q)$ のマップを等価回路に反映 |
| 脱磁余裕が厳しすぎる | 磁石温度を一律最高温度で評価している（実際は温度分布がある） | 熱解析結果の温度分布を反映。磁石中心部は周辺より低温 |
| 弱め界磁領域の鉄損が過小評価 | 基本波成分のみで鉄損を計算し、高調波鉄損を無視 | FFT分解して各次数の鉄損を個別計算し積算 |
| 非線形収束が弱め界磁条件で悪化 | ロータブリッジの深い飽和と磁石近傍の脱飽和が同時に発生し、透磁率の空間変動が大きい | ニュートン・ラフソンのライン探索を有効化、初期解に低電流ケースの解を使用 |
| 効率マップでMTPV領域が欠落 | $\psi_m/L_d > I_{\max}$ の条件（楕円中心が電流制限円の外）でMTPV軌跡が存在しないモータ設計 | 設計パラメータの見直し。$\psi_m/L_d \leq I_{\max}$ にすれば広い定出力域が実現 |
トルクリプルと高調波の問題
弱め界磁領域でトルクリプルが増大するという話を聞いたんですが、なぜですか？
弱め界磁制御ではd軸電流が大きくなるため、ステータ歯部の飽和パターンが変わる。通常運転時とは異なる高調波成分が発生して、コギングトルクやトルクリプルの波形が変化するんだ。特にIPMSMではリラクタンストルクの空間高調波成分が変わるため、基底速度以下ではリプルが小さかったのに弱め界磁域で急増することがある。
それってシミュレーションで予測できるんですか？
もちろん。FEMの過渡解析で弱め界磁条件（$i_d$ を大きくした状態）のトルク波形を計算すればリプルが直接出る。ただし精度を出すには回転角度刻みを十分小さくすること（電気角0.5〜1°以下）と、エアギャップのメッシュを十分細かくすることが条件だ。スキュー角の最適化による対策を検討する場合は3D解析が必要になることもある。
デバッグチェックリスト
シミュレーション結果がおかしいときの確認手順をまとめてもらえますか？
弱め界磁シミュレーションの系統的なデバッグ手順を教えるよ。
- 単位系の確認： $\psi_m$ [Wb], $L_d$/$L_q$ [H], $i_d$/$i_q$ [A], $\omega_e$ [rad/s] が整合しているか
- 無負荷逆起電力の検証： FEMの無負荷解析で得た逆起電力が、$E = \omega_e \psi_m$ と一致するか確認。一致しなければ磁石材料定義に問題あり
- インダクタンスの妥当性： FEMで計算した $L_d$, $L_q$ が教科書的な範囲（$L_q/L_d$ = 2〜8）に入っているか
- トルクの理論値比較： FEMトルクと $T = \frac{3}{2}p[\psi_m i_q + (L_d - L_q)i_d i_q]$ を比較。大きなズレは飽和の影響か定式化の問題
- 電圧制限の確認： 基底速度での電圧が $V_{\max}$ に到達するか確認。到達しなければ $V_{\max}$ の定義か巻線抵抗の問題
- 脱磁チェック温度： 減磁曲線が正しい温度のものか。20°Cのデータで150°C条件を評価していないか
この手順をひとつずつ潰していけば、大抵の問題は解決できそうですね。ありがとうございます！
そうだね。弱め界磁解析は電気・磁気・熱の3つの物理が絡み合う複合問題だから、1つずつ原因を切り分けるのが鉄則だ。「全部一度に直そう」とすると泥沼にハマるよ。
よくある間違い：$V_{\max}$ の定義
DC母線電圧が400Vのとき、$V_{\max}$ はいくつか？ 答えは制御方式による。SVPWMなら $V_{\max} = 400/\sqrt{3} \approx 231$V（線間電圧実効値の相電圧換算）。過変調やパルスモードを使えばもっと上げられるが、高調波が増える。この1つの数値の違いで効率マップの高速域が大きく変わるので、制御チームとの認識合わせが極めて重要だ。
なった
詳しく
報告