永久磁石同期モータ（PMSM）の設計シミュレーション

ざっくり言うと、ロータに永久磁石を使っているのがPMSM、ロータにも電流を流すのが誘導モータだ。PMSMはロータの銅損がゼロだから効率が高い。トヨタの初代プリウス（1997年）がIPMSMを採用して以来、EV/HEVの駆動モータの主流になった。

磁石をロータ表面に貼るのがSPM（Surface Permanent Magnet）、鉄心の中に埋め込むのがIPM（Interior PM）だ。IPMは「突極性」を持つから、磁石トルクに加えてリラクタンストルクも使える。EV用は大半がIPMだね。

初期設計は磁気等価回路で素早く回し、詳細設計でFEMを投入する。特にd-q軸インダクタンスの正確な値はFEMでないと出せない。磁気飽和の影響が大きいIPMSMでは、等価回路だけだと10〜20%のトルク誤差が平気で出るからね。

3相の交流量をロータに固定された回転座標系に変換すると、定常状態では直流量になる。制御もFEMとの連携も格段にやりやすくなるんだ。Park変換（dq変換）と呼ばれる手法で、PMSMの制御・解析の基礎中の基礎だよ。

$\omega_e = p \cdot \omega_m$ は電気角速度で、$p$ は極対数、$\omega_m$ は機械角速度だ。$\psi_m$ は永久磁石による磁束鎖交数で、FEMで求める場合は磁石のみで励磁した状態のコイル鎖交磁束を積分する：

3相PMSMのトルクはこう表せる：

第1項が永久磁石による磁石トルク（マグネットトルク）で、q軸電流 $i_q$ に比例する。第2項がリラクタンストルクで、$L_d \neq L_q$ の突極性がないと発生しない。SPMSMでは $L_d \approx L_q$ だからリラクタンストルクはほぼゼロ。IPMSMは $L_d < L_q$ だから $i_d < 0$ を流すことでリラクタンストルクを上乗せできる。

典型的なEV用IPMSMだと、最大トルク運転点でリラクタンストルクが全トルクの30〜50%を占めることもある。磁石使用量を減らしつつトルクを確保する設計の肝がこの突極比 $L_q / L_d$ の最適化だ。最近は希土類フリーのフェライトIPMSMもリラクタンストルク主体で実用化が進んでいるよ。

逆起電力 $e = \omega_e \psi_m$ は回転速度に比例する。高速域では逆起電力がインバータの出力電圧上限に達するから、弱め磁束制御（$i_d < 0$）で $\psi_m$ を見かけ上減らす必要がある。FEMではBack-EMFの波形を計算して高調波成分（THD）を評価するのが重要だ。THDが大きいとトルクリプルや振動騒音の原因になる。

d軸は磁石が並んでいる方向で、磁石の比透磁率は空気とほぼ同じ（$\mu_r \approx 1.05$）だから磁気抵抗が大きく $L_d$ は小さくなる。一方、q軸は磁石と磁石の間——鉄心を通る磁路だから磁気抵抗が小さく $L_q$ が大きくなる。だからIPMSMでは常に $L_q > L_d$ になる。

そこが重要なポイントだ。鉄心が磁気飽和すると透磁率が下がるから、$L_d$ も $L_q$ も電流依存になる。さらにクロスカップリング（d軸電流がq軸インダクタンスに影響する）もある。だからFEMで $i_d$-$i_q$ マップ全体にわたって $L_d(i_d, i_q)$, $L_q(i_d, i_q)$ を計算するんだ。これがPMSM設計でFEMが不可欠な最大の理由と言ってもいい。

モータの動作周波数（数百Hz〜数kHz）では準静的近似が使える。つまり変位電流項を無視して、磁気ベクトルポテンシャル $\mathbf{A}$ を用いた定式化が標準だ。2D断面解析なら $A_z$ のスカラー方程式に帰着する：

$\mathbf{J}_s$ はステータ巻線の強制電流密度、$-\sigma \partial \mathbf{A}/\partial t$ は渦電流項（鉄心やロータ導体の渦電流）、$\nabla \times \mathbf{M}$ は永久磁石の等価電流源だ。2D解析では辺要素（Nedelec要素）または節点要素が使われる。モータの場合は2D三角形の節点要素で十分なケースが多い。

まさにそこがPMSM解析の肝だ。珪素鋼板のB-Hカーブは強い非線形性を持つ。磁束密度が1.5Tを超えると透磁率が急激に下がる。これをニュートン・ラフソン法で処理する。各反復ステップで $\nu(B)$ を更新しながらヤコビアンマトリクスを再構築する：

収束判定は通常 $||R|| / ||R_0|| < 10^{-4}$ 程度。IPMSMでは歯先やブリッジ部が2.0T近くまで飽和するから、B-Hカーブのデータが2.0T以上まで入っていないと「飽和領域の外挿」になって結果が大きく狂う。実測データの品質がFEM結果の信頼性を左右するんだ。

エアギャップにスライディングインターフェース（スリップ面）を設ける方法が標準だ。ロータ側メッシュを剛体回転させ、エアギャップの境界でステータ側と補間接続する。JMAGでは「スライドメッシュ法」、Maxwellでは「Band」オブジェクトがこれに相当する。再メッシュ不要で計算効率が高い。

電気角1度あたり1ステップ（つまり1電気周期360ステップ）が出発点だ。コギングトルクを正確に評価したいなら0.5度以下が必要。回転速度×時間刻み＝角度刻みだから、例えば3000rpm・4極のモータなら電気角周波数は100Hzで、0.5度刻みだと $\Delta t = 0.5/(360 \times 100) \approx 13.9 \mu s$ になる。

基本的な手法は「凍結透磁率法」（Frozen Permeability法）だ。手順はこうだ：

その通り。$i_d$ を -200A〜0A、$i_q$ を 0A〜400A くらいの範囲で格子状に計算すると、$L_d$-$L_q$ マップが得られる。これが制御シミュレーション（Simulink等）に渡すテーブルデータになる。JMAGではこの計算を自動で回す「効率マップ解析」機能がある。

鉄損は一般にSteinmetzの経験式をベースにした分離モデルで計算する：

FEMでは各要素の磁束密度波形を周波数分解し、各高調波成分に対して上記の式を適用する。これをポスト処理で積分して全体の鉄損を得る。最近のソフトはPWMキャリア高調波による追加鉄損も考慮できるようになっている。

典型的なフローはこうだ：

2D解析の計算時間は3Dの1/100以下だから、設計初期のパラメータ探索には2Dが断然効率的だ。端部効果やスキューの影響は2D結果に補正係数を掛けることで対応できる。3Dが本当に必要なのは、端部漏れ磁束が大きいアキシャルギャップモータや、複雑なロータ形状を持つモータだけだよ。

PMSM解析で最も重要なメッシュ領域はエアギャップだ。ここに最低3〜5層の要素を配置する。エアギャップ長が0.5mmのモータなら、要素サイズは0.1〜0.15mm程度になる。

コギングトルクのピーク値で判断するのが最も厳しい（感度が高い）。メッシュを半分に細かくしてコギングトルクの変化が5%以内なら十分。平均トルクは割と早く収束するから、コギングが一番シビアなインジケータだ。

8極48スロットのモータなら、最小周期が1/8モデル（1極分、45度）になる。ただし周期境界条件の設定に注意がいる。PMSMは磁石の極性が交互だから、奇数倍の周期では「反周期境界条件」（$A_z$ の符号が反転）を使う必要がある。偶数倍なら通常の周期境界でOKだ。

外周の境界条件は $A_z = 0$（ディリクレ条件）が標準。ステータ外周を十分遠方と仮定する。実際のモータではケースがあるから、ケースまで含めたモデルの方が正確だけど、設計初期段階では $A_z = 0$ で十分だよ。

パラメータを変えながらFEMを自動で回すパラメトリックスタディが基本だ。例えば磁石厚、磁石幅、フラックスバリア形状をパラメータ化して、トルク最大化・コギング最小化・材料コスト最小化を目的関数にする。最適化アルゴリズムは：

そうだね。JSOL社が開発していて、トヨタ・日産・ホンダなど国内の自動車OEMでのシェアが圧倒的に高い。最大の強みはモータ設計に特化したワークフローで、テンプレートからIPMSMやSPMSMを数分で定義できる。効率マップ解析、$L_d$-$L_q$マップ自動生成、NVH連携（音圧レベルの評価）まで一気通貫で回せる。

汎用電磁界ソルバーとしては範囲が限られる。アンテナやEMC解析には向かない。あくまで「回転機・変圧器・アクチュエータ専用」だ。あと海外ではAnsys Maxwell系のシェアが大きいから、グローバルチームではツール統一が難しい場合もある。

Ansys Maxwell（旧Ansoft Maxwell）は低周波電磁界の汎用ソルバーで、モータだけでなく変圧器、リレー、ソレノイド、磁気シールドなど幅広い用途に対応する。最大の武器はAnsysエコシステムとの連携だ。Mechanical（構造）、Fluent/Icepak（熱流体）、Twin Builder（システムシミュレーション）との双方向連成が自然に組める。

もう一つの特徴はアダプティブメッシュ。解の精度が不十分な領域を自動で細分化してくれる。初心者がメッシュ設計に悩まなくても、ある程度の精度が出せるのは大きい。ただ、細分化が過剰になって計算時間が爆発するケースもあるから、上限は設定しておいた方がいい。

AC/DCモジュールでPMSM解析は可能だ。COMSOLの真価はマルチフィジックス。電磁-熱-構造を同一メッシュで直接連成できるから、例えば「高速回転時の渦電流損→温度上昇→磁石の特性変化→トルク低下」をワンストップで計算できる。研究開発やカスタム物理の実装には最適だけど、量産設計の効率マップを何百ケースも回すような用途ではJMAGやMaxwellの方がワークフローが洗練されている。

いくつかあるよ：

PMSM解析での非線形収束失敗は、ほぼ100%がB-Hカーブ関連だ。チェックすべき点：

それが原因だ。応急処置として、飽和域は $\mu_0$ に漸近する外挿をかける。具体的には1.8T以降のデータ点を $B = \mu_0 H + J_s$（$J_s$は飽和磁化）で延長する。35A300のような一般的な珪素鋼板なら $J_s \approx 2.0$ T だ。

コギングトルクは全トルクの1%以下のオーダーだから、わずかなモデル誤差で大きく変動する。よくある原因：

最悪条件でFEMを回して、磁石内部の各要素で動作点が膝点（不可逆減磁が始まるH値）を超えていないか確認する。条件は：

3つの鉄則を覚えておけ：