モータ効率マップ生成 — 電磁界FEMによる損失分離と等効率線の作成手法
理論と物理
効率マップとは何か
効率マップって何に使うんですか？ 名前だけは聞いたことがあるんですけど、いまいちピンと来なくて…
横軸に回転数、縦軸にトルクを取って、各動作点でのモータ効率を等高線で描いたものだよ。地図の等高線と同じ要領で、「効率95%の線」「90%の線」…と描いていく。
なるほど、等高線なんですね！ でもなぜそれがそんなに重要なんですか？
EVの走行パターンでの平均効率を評価するためだ。WLTPやJC08といった走行モードでどの動作点を多く使うかが分かるから、その頻出領域に高効率の「島（アイランド）」が来るよう設計を最適化できる。例えば高速道路を多く走るEVなら、8000〜12000 rpm付近に高効率域を持って行きたいし、市街地メインなら3000〜5000 rpm付近に寄せたい。
走行パターンに合わせて「島の位置」を調整するってことですか？ カーナビの経路最適化みたいですね。
まさにそう。理論上の最大効率点が1つあるだけでは不十分で、実際の使い方に合わせた「実用効率」を最大化する設計思想が、モータ設計の現場では最も重要なんだ。
損失分離の物理
効率を計算するには損失を知る必要がありますよね。モータの損失ってどういう種類があるんですか？
モータの全損失は大きく3つに分離できる。
- 鉄損 $P_{\mathrm{iron}}$：コア（鉄心）内の磁束密度変動に起因する損失。ヒステリシス損と渦電流損に分かれる
- 銅損 $P_{\mathrm{cu}}$：巻線の電気抵抗による発熱損失（ジュール損）
- 機械損 $P_{\mathrm{mech}}$：軸受摩擦と回転子の風損（ウィンデージ）
この3つを正確に分離・計算できるかが、高品質な効率マップの鍵だ。
低回転と高回転で、どの損失が支配的かって変わりますか？
いいところに気づいたね。低回転・高トルク域では銅損が圧倒的に大きい。大きな電流を流すから $I^2R$ が効く。一方、高回転域になると鉄損が急増する。鉄損は周波数（≒回転数）の1〜2乗に比例するからね。つまり回転数が上がると「銅損から鉄損へ主役が交代」するんだ。効率マップの等高線の形は、この2つのせめぎ合いを反映しているんだよ。
鉄損とSteinmetz式
鉄損の計算式って、具体的にはどうなってるんですか？
鉄損のモデルとして最も広く使われているのが修正Steinmetz式（iGSE: improved Generalized Steinmetz Equation）だ。元のSteinmetz式は正弦波磁束用だけど、モータ内の磁束波形は高調波を含むからiGSEが必要になる。
古典的なSteinmetz式は、単位体積あたりの鉄損を次のように表す：
各項の物理的意味：
- 第1項 $k_h f B_m^{\alpha}$：ヒステリシス損。磁区壁の移動に伴うエネルギー散逸。周波数 $f$ に1次比例。$\alpha$ は材料定数（典型的に1.6〜2.2）
- 第2項 $k_e (f B_m)^2$：古典渦電流損。鋼板内に誘起される渦電流によるジュール発熱。周波数の2乗に比例
- 第3項 $k_a (f B_m)^{1.5}$：異常渦電流損（excess loss）。磁区構造の動的変化に起因。Bertotti理論で説明される
高回転になると第2項と第3項がガンガン効いてくるわけですね。だから薄い鋼板を使うんですか？
その通り。渦電流損は鋼板厚 $d$ の2乗に比例するんだ。$k_e \propto d^2 / (12 \rho)$（$\rho$ は鋼板の電気抵抗率）。だから高速モータでは0.2〜0.25mm厚の電磁鋼板を使う。さらに低鉄損を追求するならアモルファス鉄心も選択肢に入る。ただしコストが跳ね上がるけどね。
銅損（I²R損失）
銅損は巻線の抵抗損失であり、dq軸電流表現では：
ここで $R_s$ は1相あたりの巻線抵抗、$i_d$, $i_q$ はdq軸電流成分である。
$R_s$ って温度で変わりますよね？ 巻線が熱くなると抵抗が上がって損失も増える…
鋭いね。銅の抵抗温度係数は約0.393%/℃だ。20℃基準で $R_s(T) = R_{s,20}[1 + 0.00393(T - 20)]$ となる。巻線温度が120℃まで上がると抵抗は約39%も増加する。だから正確な効率マップには電磁-熱連成解析が欠かせないんだよ。さらに高周波では表皮効果と近接効果で実効抵抗が増大する。これをACロスと呼ぶ。
機械損
機械損は電磁界解析では求められないですよね？ どうやって効率マップに入れるんですか？
そう、機械損は電磁FEMの直接の計算対象ではない。一般に次の経験式で見積もる：
$$ P_{\mathrm{mech}} = P_{\mathrm{bearing}} + P_{\mathrm{windage}} = k_b \omega + k_w \omega^3 $$
軸受損失 $k_b \omega$ は回転数に1次比例、風損 $k_w \omega^3$ は3乗で増加する。実測データから $k_b$, $k_w$ を同定して効率マップに組み込むのが一般的だ。
効率の定式化
各動作点（トルク $T$、角速度 $\omega$）における効率は：
式自体はシンプルですね。分母の各損失項をどれだけ正確に求められるかが勝負ってことですか。
まさにそう。式は中学生でも理解できる「出力/(出力+損失)」だけど、分母の $P_{\mathrm{cu}}$, $P_{\mathrm{iron}}$, $P_{\mathrm{mech}}$ を各動作点で精密に求めるのに電磁界FEMの全力が必要になる。特に $P_{\mathrm{iron}}$ の精度が効率マップの信頼性を左右する。
各損失項の物理的意味
- ヒステリシス損 $k_h f B_m^{\alpha}$：磁区壁が往復運動するたびに結晶格子との摩擦でエネルギーが散逸する。B-Hループの面積に相当。【日常の例】消しゴムで何度も紙をこすると熱くなる——磁区壁が「こすれて」発熱するイメージ。
- 渦電流損 $k_e (f B_m)^2$：磁束の時間変化がレンツの法則で鋼板内に渦状の電流を誘起し、その電流がジュール発熱する。IHクッキングヒーターの加熱原理と同じ。鋼板を薄くすると渦電流のループが小さくなり損失が減る。
- 銅損 $\frac{3}{2}R_s(i_d^2+i_q^2)$：電流が巻線の抵抗を通過する際のジュール熱。3/2の係数は3相交流のdq変換に由来する。ドライヤーのニクロム線が熱くなるのと同じ原理。
- 風損 $k_w \omega^3$：回転子と周囲の空気との摩擦。自転車の高速走行時に感じる空気抵抗が3乗で増えるのと同じ物理。超高速モータ（30,000 rpm以上）では無視できない。
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 磁束密度 $B$ | T（テスラ） | 1T = 1 Wb/m²。電磁鋼板の飽和磁束密度: 1.5〜2.0 T |
| 周波数 $f$ | Hz | 電気角周波数 = (極対数) × (回転数[rpm]) / 60 |
| 電流 $i_d, i_q$ | A | dq軸電流。ピーク値と実効値の使い分けに注意 |
| 巻線抵抗 $R_s$ | $\Omega$ | 温度依存性あり。20℃基準値に温度補正を適用 |
| 角速度 $\omega$ | rad/s | $\omega = 2\pi n / 60$（$n$: 回転数 [rpm]） |
| トルク $T$ | N$\cdot$m | 出力 $P = T\omega$ |
| Steinmetz係数 $k_h, k_e, k_a$ | 各々異なる | メーカーの鉄損カーブから最小二乗法で同定 |
テスラ・モデル3のモータはなぜ「IPMSynRM」なのか
テスラ・モデル3のリアモータはIPM（埋込磁石）とSynRM（シンクロナスリラクタンス）のハイブリッド構造を採用している。磁石トルクとリラクタンストルクの両方を使うことで、効率マップの高効率島を広い回転数域に拡張できるからだ。高速域では弱め磁束制御（$i_d < 0$）によって鉄損を抑え、低速域では磁石トルクでピーク効率を稼ぐ。この「いいとこ取り」戦略は、効率マップの等高線がどういう物理で決まるかを深く理解しているエンジニアだけが設計できる——だからこそ効率マップの読み書きはモータ設計者の必須スキルなんだ。
数値解法と実装
電磁界FEMの定式化
効率マップを作るために、電磁界FEMではどんな方程式を解くんですか？
モータ解析の出発点は、ベクトルポテンシャル $\mathbf{A}$ を用いた拡散方程式だ。2Dの場合、断面内の磁束密度は $z$ 成分のベクトルポテンシャル $A_z$ のみで表現できる：
$$ \nabla \times \left( \frac{1}{\mu} \nabla \times \mathbf{A} \right) + \sigma \frac{\partial \mathbf{A}}{\partial t} = \mathbf{J}_s + \nabla \times \mathbf{M} $$
左辺第1項が磁場の拡散、第2項が渦電流項。右辺の $\mathbf{J}_s$ がコイル電流源、$\mathbf{M}$ が永久磁石の磁化だ。
これを時間方向に解いて、1回転分の磁束密度の時間波形を得るわけですか？
そう。時間ステップ法（Crank-Nicolson法やBackward Euler法）で離散化して、各時刻での磁束密度分布を求める。FEM離散化後の行列方程式は：
$$ \left( [S] + \frac{1}{\Delta t} [M] \right) \{A\}^{n+1} = \{F\}^{n+1} + \frac{1}{\Delta t} [M] \{A\}^n $$
$[S]$ が剛性行列（透磁率依存）、$[M]$ が質量行列（導電率依存）、$\{F\}$ が右辺ベクトル（電流源+永久磁石）だ。磁気飽和がある場合、$[S]$ は $B$ に依存するから非線形になる。Newton-Raphson法で各時間ステップを反復的に解く必要がある。
動作点スイープの実装
効率マップ全体を作るには、トルク-回転数平面上の多数の動作点で解析を回す必要がありますよね？ それって膨大な計算量になりませんか？
まさにそこが効率マップ作成の最大の計算コスト課題だ。例えば回転数10段階 × トルク10段階 = 100点の格子を切ると、各点で電気角1周期分の過渡解析が必要になる。1点あたり2D解析で30秒かかるとして、100点なら50分。3D解析だと1点30分で100点なら50時間。
50時間…！ 何とか減らす方法はないんですか？
いくつかの高速化テクニックがある：
- DOE（実験計画法）+ サロゲートモデル：ラテン超方格法やOptimal LHSで50〜100点をサンプリングし、応答曲面（クリギング、RBF、ガウス過程回帰）で補間する方法。計算点数を1/3〜1/10に削減できる
- 周波数応答法：正弦波定常仮定が成り立つ動作点では、過渡解析の代わりに複素数演算の周波数領域解析で済ませる
- 並列計算：各動作点は互いに独立なので、embarrassingly parallelに分散処理できる。100コアあれば100点同時に走らせられる
- 2D+端部補正：2D断面解析に3Dエンド効果の補正係数を後付けする方法。3Dフル解析の1/100以下のコストで80〜90%の精度を確保
dq軸電流制御と動作点設定
各動作点の電流条件はどうやって決めるんですか？ トルクと回転数が決まれば電流は自動的に決まるわけではないですよね？
その通り。同じトルクを出すにも $i_d$（磁束弱め方向）と $i_q$（トルク生成方向）の組み合わせは無数にある。IPMSMのトルク式を見てみよう：
$$ T = \frac{3}{2} p \left[ \psi_m i_q + (L_d - L_q) i_d i_q \right] $$
第1項が磁石トルク、第2項がリラクタンストルクだ。$p$ は極対数、$\psi_m$ は磁石鎖交磁束、$L_d, L_q$ はdq軸インダクタンス。効率マップ作成では、各トルク-回転数の組み合わせに対してMTPA（Maximum Torque Per Ampere）制御または最大効率制御で最適な $(i_d, i_q)$ を決定する。
高回転で電圧制限に引っかかったらどうするんですか？
電圧制限楕円と電流制限円の交点で動作点を決める弱め磁束（Field Weakening）制御に切り替える。$i_d$ を負方向に振って磁束を弱め、逆起電力を下げることで高回転を実現する。このとき効率は下がるが、動作範囲が広がる。効率マップにはこの弱め磁束領域が必ず含まれるから、MTPA領域とFW領域の境界（基底回転数）をきちんとモデル化することが重要だ。
鉄損の後処理計算
FEM解析で磁束密度の時間波形が得られたとして、そこからどうやって鉄損を計算するんですか？
典型的な手順はこうだ：
- FEM過渡解析で各要素の $B_x(t), B_y(t)$ の時間波形を取得
- FFTで周波数分解し、基本波と各高調波の振幅 $B_k$ と周波数 $f_k$ を求める
- 各高調波成分にSteinmetz式を適用し、高調波鉄損を足し合わせる： $$ P_{\mathrm{iron}} = \sum_{k=1}^{N} \left[ k_h f_k B_k^{\alpha} + k_e (f_k B_k)^2 + k_a (f_k B_k)^{1.5} \right] $$
- 全要素について体積積分し、モータ全体の鉄損を算出
JMAGやMaxwellではこの後処理が自動化されていて、鉄損コンター図を直接可視化できるよ。
非線形収束と磁気飽和の取り扱い
モータの鉄心って磁気飽和しますよね？ FEMでそれをどう扱うんですか？
電磁鋼板のB-Hカーブは強い非線形性を持つ。$B$ が1.5T付近で透磁率が急激に低下する飽和域に入る。FEMでは各要素の透磁率 $\mu$ を前ステップの $B$ から逐次更新するNewton-Raphson反復を適用する。収束判定基準は通常 $||R||/||R_0|| < 10^{-4}$ 程度。
高トルク動作点では歯部やブリッジ部が深い飽和に入るため、反復回数が増えて計算時間が長くなる。収束を助けるテクニックとして：
- B-Hカーブの滑らかな補間（スプラインではなくAkima補間が安定）
- 緩和係数（under-relaxation factor 0.5〜0.8）の導入
- 前ステップの解を初期値として使用（continuation法）
FEM解析のたとえ
モータの電磁界FEM解析は「水彩画を描く工程」に似ている。メッシュは画布の格子、磁束密度は各マス目に塗る色の濃さだ。非線形のB-Hカーブは「色が濃くなると絵の具が固まって塗りにくくなる」ようなもの——飽和域では計算（筆の動き）が重くなる。そして効率マップは「100枚の水彩画を微妙にパラメータを変えて描き、全体の傾向を俯瞰する」作業に相当する。
実践ガイド
解析ワークフロー
実際に効率マップを作るには、最初から最後まで何をどういう順番でやるんですか？
全体の流れを5ステップで説明しよう：
- Step 1: モータ形状の作成 — 2D断面モデルをCADまたはツール内蔵エディタで作成。対称性を利用して1極対分（例：8極なら45°分）に削減
- Step 2: 材料設定 — 電磁鋼板のB-Hカーブ（非線形）、Steinmetz係数、巻線の導電率・温度係数、永久磁石の残留磁束密度 $B_r$ と保磁力 $H_{cJ}$ を設定
- Step 3: メッシュ生成＋回路連成 — エアギャップ・歯先に細かいメッシュ。外部回路（インバータ）をFEMと連成させ、dq軸電流制御を実装
- Step 4: 動作点スイープ実行 — 回転数×トルクの格子点を定義し、各点で過渡解析（電気角1〜2周期分）を実行。JMAGならEfficiency Map Study機能で自動化可能
- Step 5: 後処理 — 各動作点の鉄損・銅損・機械損を集計し、$\eta$ を計算。等効率線のコンター図を描画
メッシュ戦略
モータ解析のメッシュで特に注意すべき場所ってどこですか？
モータ特有のメッシュ要件がいくつかある：
| 領域 | 推奨要素サイズ | 理由 |
|---|---|---|
| エアギャップ | 3〜5層以上 | トルク精度に直結。磁束密度の円周方向分布を正確に捉える |
| 歯先・スロット開口部 | 鋼板厚の1/2以下 | 渦電流分布と鉄損精度に影響 |
| 磁石-鉄心境界 | 磁石厚の1/3以下 | 減磁解析の精度確保 |
| コアバック | やや粗くてOK | 磁束密度が比較的一様 |
| 外部空気領域 | 粗いメッシュ | 遠方では磁束密度が急減するため精度への影響小 |
回転するロータとステータのメッシュの接続はどうするんですか？
エアギャップ中間面にスライディングメッシュ（スリップインターフェース）を設定する。ロータメッシュが回転しても、ステータメッシュとの結合を補間で維持する仕組みだ。JMAGでは「MotionMesh」、MaxwellではBand（バンド領域）がこれに該当する。メッシュの節点がエアギャップ中間面で等間隔になるよう制約をかけると補間精度が上がるよ。
WLTP走行モードと設計最適化
WLTPモードで効率マップを使う方法をもう少し具体的に教えてください。
WLTPの速度プロファイルから、各時刻のモータ動作点（トルク, 回転数）を車両走行方程式で計算する。そして効率マップ上に動作点の軌跡をプロットすると、「この走行モードではどの領域を頻繁に使うか」が一目で分かる。
加重平均効率は次のように定義する：
$$ \bar{\eta}_{\mathrm{WLTP}} = \frac{\sum_{i} P_{\mathrm{out},i} \cdot \Delta t_i}{\sum_{i} P_{\mathrm{in},i} \cdot \Delta t_i} = \frac{\sum_{i} T_i \omega_i \Delta t_i}{\sum_{i} \frac{T_i \omega_i}{\eta_i} \Delta t_i} $$
この $\bar{\eta}_{\mathrm{WLTP}}$ を最大化するように、磁石量・スロット形状・巻数・鋼板材質を最適化するのが実務の設計フローだ。
なるほど…効率マップはモータ単体の性能図じゃなくて、「車全体のシステム効率を最適化するための道具」なんですね！
そう。だからOEMのEV開発部門では、モータ設計者とバッテリー・制御・車両運動のエンジニアが効率マップを共通言語にして議論するんだ。「この走行モードで加重平均効率を95%以上にしたい」という目標値が、設計パラメータの最適化を駆動する。
実測との突き合わせ
FEMで作った効率マップと実際のモータの効率マップって、どれくらい合うものなんですか？
適切にモデル化すれば、効率の絶対値で±1〜2%以内の精度が期待できる。ただし以下の要因で乖離が生じやすい：
- 鉄損モデルの不備：加工ストレスによるB-Hカーブの劣化（パンチング応力効果）が未反映。実測鉄損が解析値の1.3〜2倍になることも
- ACロスの無視：高速域で巻線の表皮効果・近接効果による銅損増加が見落とされている
- 温度依存性：磁石の $B_r$ は温度係数 −0.1〜−0.12%/℃で低下し、高温で大きくトルクが減少する
- 機械損の見積もり不足：シール摩擦やオイル撹拌損（油冷モータの場合）が経験式で十分カバーされていない
V&V（検証・妥当性確認）のコツ
効率マップの実測検証では、まず無負荷試験で鉄損+機械損の合計を確認し、次に拘束試験で銅損を単独で検証するのが王道だ。全損失を一度に比較しようとすると、どの損失モデルが間違っているのか切り分けられない。「一度に1つずつ検証する」という科学実験の原則を守ろう。
効率マップの「島」を読む——プロが設計で意識するポイント
効率マップには必ず「高効率の島（アイランド）」がある。島の形と位置は設計の本質を反映しており、島が細長い場合は使用回転域が限定的、丸く広い場合はロバストな設計と読める。実務では車両のWLTC走行プロファイルを効率マップに重ね、走行時間の多い動作点が高効率域に入るよう島の位置を調整する「マップ最適化」が重要だ。理論上の最大効率点より、使用状況に合わせた「実用効率」を最大化する設計思想が、現場での差を生む。
ソフトウェア比較
主要ツールの効率マップ機能比較
効率マップを作れる電磁界解析ソフトって、どれがおすすめですか？
主要なツールを効率マップ作成の観点で比較してみよう：
| 機能 | JMAG-Designer | Ansys Maxwell | Motor-CAD | COMSOL |
|---|---|---|---|---|
| 効率マップ自動生成 | ◎ 専用機能あり | ○ スクリプトで対応 | ◎ 統合機能 | △ 手動構築 |
| 鉄損モデル（Steinmetz） | ◎ iGSE + Bertotti | ○ Bertotti 3項 | ○ Steinmetz | ○ カスタム可 |
| ACロス計算 | ○ | ○ | ◎（解析的） | ○ |
| 電磁-熱連成 | ○ | ◎（Icepak連携） | ◎（内蔵熱解析） | ◎（マルチフィジクス） |
| MTPA/FW制御 | ◎ 回路連成 | ○ 外部回路 | ◎ 内蔵制御 | △ 手動実装 |
| DOE/最適化 | ○（外部連携） | ◎（optiSLang） | ○ | ○ |
| 計算速度（2D） | ◎ | ○ | ◎（解析的併用） | △ |
JMAG-Designerでのワークフロー
JMAGで効率マップを作る手順を具体的に教えてください。
JMAGには「効率マップスタディ」という専用機能がある。手順はこうだ：
- ベースモデルの作成：2D過渡解析モデルを1動作点で作成・検証する
- 効率マップスタディの設定：回転数範囲（例：500〜15000 rpm）とトルク範囲（例：0〜300 Nm）、格子点数を指定
- 制御モードの選択：MTPA + 弱め磁束制御。電圧制限・電流制限を設定
- 鉄損モデルの選択：材料DBからSteinmetz係数を自動取得、または手動入力
- 実行：JMAGが自動的に全格子点をスイープし、各動作点の最適 $(i_d, i_q)$ を探索しながら解析を実行
- 結果表示：効率マップ、損失マップ、電流ベクトル軌跡が自動生成される
JMAGの強みは、日本の自動車メーカー・サプライヤーと長年共同開発してきた実績に裏打ちされた材料DB（電磁鋼板データ）の充実度だ。
Ansys Maxwellでのワークフロー
Maxwellではどうやるんですか？ JMAGほど自動化されてない印象があるんですけど…
MaxwellにはJMAGのような専用効率マップ機能はないけど、Ansys Motor-CADと連携するのが現在の正攻法だ。Motor-CADで解析的に高速な効率マップを生成し、気になる動作点だけMaxwell 2D/3Dで精密検証する「ハイブリッドアプローチ」が実務では主流になっている。
MaxwellだけでやるならPythonスクリプト（IronPython）でパラメトリックスイープを自動化する。Ansys Electronics Desktopの「Optimetrics」でDesign of Experimentsを設定し、各動作点を並列実行する方法もある。
オープンソース代替手段
商用ソフトのライセンス料が高くて手が出ません…無料で効率マップを作る方法はありますか？
いくつかのオープンソースの選択肢がある：
- FEMM (Finite Element Method Magnetics)：David Meeker開発のフリーウェア。2D静磁場・渦電流解析が可能。Lua/Pythonスクリプトでパラメトリックスイープを組める。効率マップの自動生成にはスクリプティングが必要だが、学生・研究者には十分な精度
- Elmer FEM：CSC（フィンランド科学計算センター）開発のOSS。3D電磁場解析に対応するが、モータ設計向けのGUIや効率マップ機能はない
- PyMotor / ADEPT：Python製のモータ設計ツールキット。FEMMのPythonラッパーを使って効率マップを自動生成する研究用ツール
OSSの最大の弱点は鉄損モデルの精度と材料DBの欠如だ。商用ツールのSteinmetz係数DBに匹敵するデータを自前で整備する必要がある。
効率マップを「自動生成」する時代——DOEとサロゲートモデルの組み合わせ
かつては効率マップ1枚を作るのに数十時間の計算が必要だった。トルク×回転数の格子点を一つひとつ解析するため、100点×100点のマップなら1万回の解析が必要になる。現代ではDOE（実験計画法）でサンプル点を絞り、その結果からガウス過程回帰などのサロゲートモデルを構築することで、100〜300点の解析から全マップを補間する手法が主流になっている。Ansys OptiSLangやAltair HyperStudyがその自動化を担うが、サンプリング戦略を誤ると効率の「谷」を見落とすリスクがある。
トラブルシューティング
鉄損が実測と合わない
効率マップを作ったんですけど、高回転域で実測と全然合いません…鉄損が過小評価されているみたいです。
これは非常によくある問題だ。原因と対策を整理しよう：
| 原因 | 詳細 | 対策 |
|---|---|---|
| 加工劣化の未考慮 | パンチングやレーザー加工で鋼板端部のB-Hカーブが劣化。実効鉄損が1.3〜2倍に増加 | 加工影響係数（Building Factor: 1.3〜1.8）を適用。JMAG v23+では加工劣化モデルをサポート |
| Steinmetz係数の温度依存性 | 120℃ではヒステリシス損が20℃比で20〜30%変化する場合がある | 温度依存のSteinmetz係数を使用。複数温度点でのデータフィッティング |
| 回転磁界成分の無視 | Steinmetz式は交番磁界用。ステータ歯部では円形・楕円回転磁界が発生 | 回転鉄損モデル（回転・交番分離型）の使用 |
| DC磁化バイアス | 永久磁石のDCバイアスでヒステリシスループが非対称になる | マイナーループモデルの適用 |
高効率島の位置ずれ
効率マップの高効率島が、実測に比べて高回転側にずれてしまいます。何が原因でしょうか？
高効率島の位置ずれは、インダクタンスの推定精度に起因することが多い。具体的には：
- $L_d, L_q$ の飽和依存性：高負荷時に歯部・ブリッジ部が飽和すると $L_q$ が低下し、リラクタンストルクが減少する。結果として同じトルクを出すのにより多くの電流が必要になり、銅損が増えて効率が下がる。FEMでこの飽和依存性を正しく捉えているか確認しよう
- エンドコイルの漏れインダクタンス：2D解析では3Dのエンドコイル漏れ磁束が無視される。この分 $L_d$ が過小評価され、弱め磁束制御の効きが解析と実機でずれる
- 磁石の温度劣化：NdFeB磁石の $B_r$ は-0.11%/℃で低下。80℃→120℃の温度上昇で $B_r$ が約4.4%低下し、磁石トルクが減少する
非線形収束失敗
高トルクの動作点で解析が収束しなくなりました。Newton-Raphson反復が発散してしまいます。
高トルク動作点では鉄心が深い飽和に入るため、透磁率の急激な変化がNewton-Raphson法の収束を妨げる。以下の対策を試してみよう：
- B-Hカーブの高磁場域を延長：実測データが $B = 2.0$ T までしかない場合、その先を適切に外挿する（真空透磁率に漸近させる）。外挿が不適切だと微分（dB/dH）が不連続になり発散の原因に
- 緩和係数を0.3〜0.5に下げる：更新幅を制限して振動を抑制
- 時間刻みを半分にする：過渡解析の時間刻みが大きすぎると、前ステップの解が良い初期値にならない
- メッシュの質を確認：ブリッジ部（磁気飽和が最も厳しい箇所）に過度に歪んだ要素がないか
計算時間の爆発
効率マップの計算が3日経っても終わりません…助けてください！
落ち着いて。まず何が時間を食っているか特定しよう：
- 格子点数が多すぎないか？ — 初回は回転数5段階×トルク5段階 = 25点で十分。後からサロゲートモデルで補間すればよい
- 3D解析をしていないか？ — まず2D解析で全体傾向を掴み、気になる点だけ3Dに切り替える
- 過渡解析の電気角周期数は適切か？ — 定常に達した最後の1周期だけを使えばよい。過渡の立ち上がり期間が長すぎないか確認（初期値を前動作点の解で初期化する手もある）
- 出力データ量が過大でないか？ — 全要素・全時刻の磁束密度を保存するとディスクI/Oがボトルネックになる。鉄損後処理に必要なデータだけを保存する設定にする
- 並列計算を活用しているか？ — 各動作点は独立なので、コア数分だけ同時実行できる
ありがとうございます！ まず格子点を25に減らして2D解析からやり直してみます。
うん、それが賢明だ。効率マップ作成は「粗く速く全体を掴む → 必要な領域だけ精密化する」というアプローチが鉄則だよ。最初から完璧を目指すと計算コストで破綻する。
初心者が陥りやすい落とし穴
- 「空気領域を忘れる」：モータ外部にも磁束は漏れている。解析領域をモータ外径ぎりぎりにすると、漏れ磁束が「見えない壁」に反射して非物理的な結果になる。外径の1.5〜2倍の空気領域を確保しよう。
- 「B-Hカーブの入力ミス」：単位系の混在（CGSとSI）が最多。1 Oe = 79.577 A/m の換算を間違えるとトルクが桁違いにずれる。
- 「回転方向の符号ミス」：ロータの回転方向と電流位相の関係が逆転すると、トルクが負になったりゼロになったりする。最初の1動作点でトルク波形を確認するのを忘れずに。
- 「機械損ゼロで評価してしまう」：電磁FEMだけでは機械損は出てこない。特に高回転域で実測との乖離が大きくなる原因の筆頭。
なった
詳しく
報告