スイッチトリラクタンスモータ（SRM）の電磁界シミュレーション
理論と物理
SRMの動作原理
スイッチトリラクタンスモーターって永久磁石を使わないんですか？ レアアース不要？
その通り。ロータは単純な鉄の突極構造で磁石もコイルもない。レアアースフリーでコスト安・高温耐性が利点。ただしトルクリプルと騒音が大きいのが課題で、最近のNidec製洗濯機用SRMではAI制御で騒音を60%低減した例がある。
磁石がないのにどうやってトルクが出るんですか？ 全然イメージ湧かないです…
ざっくり言うと「鉄が磁石に引き寄せられる力」を利用している。ステータのコイルに電流を流して電磁石を作ると、ロータの鉄の突極がそこに吸い寄せられる。これがリラクタンストルクだ。磁気回路のインダクタンスが最大になる位置にロータが動こうとする性質を利用しているんだ。
なるほど！ じゃあコイルの電流をタイミングよく切り替えていけば連続回転するってことですか？
まさにそう。典型的な8/6構造（ステータ8極・ロータ6極）では、A相→B相→C相→D相と順番に通電すると、ロータの突極が次々にステータの突極に引き寄せられて回転する。ただ、相の切り替えタイミングと電流波形が命で、ここをミスるとトルクが出ないどころか逆回転することもある。
SRMの基本的な構造的特徴をまとめると以下の通りになる。
| 項目 | SRM | PMSM（参考） |
|---|---|---|
| ロータ構造 | 突極鉄心のみ（磁石・コイルなし） | 永久磁石を埋込/表面貼付 |
| レアアース使用 | なし | ネオジム磁石使用 |
| 耐熱性 | 250°C以上でも動作可能 | 磁石減磁温度で制約（約150°C） |
| トルクリプル | 大きい（20〜40%） | 小さい（3〜8%） |
| 電磁騒音 | 大きい（課題） | 小さい |
| コスト | 低い | 磁石分高い |
| 主用途 | 家電、産業用、EV（研究段階） | EV、サーボ |
トルク発生メカニズムと支配方程式
SRMのトルクの式ってどうなるんですか？ PMSMとは全然違う？
根本的に違う。PMSMのトルクは「電流×磁石の磁束」だけど、SRMはインダクタンスの変化率からトルクが出る。線形領域（磁気飽和なし）では、瞬時トルクはこうなる：
ここで $I$ は相電流 [A]、$L$ はインダクタンス [H]、$\theta$ はロータ角度 [rad] である。この式からわかる重要なポイントが3つある。
- トルクは電流の2乗に比例 — 電流の向きに無関係なので、単方向電流（ユニポーラ）で駆動できる。PMSMのような正弦波交流は不要。
- $dL/d\theta > 0$ のときのみ正トルク — インダクタンスが増加する区間で通電し、減少区間に入る前に電流を切る必要がある。
- 磁気飽和を無視した近似式 — 実機では鉄心が飽和するため、この線形式だけでは不十分。
飽和するとどうなるんですか？ 式が変わる？
そう、飽和領域ではインダクタンスが電流にも依存するから $L(\theta, I)$ になる。この場合、正確なトルク計算にはコエナジー（共エネルギー）からの導出が必要だ：
$\psi(\theta, i)$ は鎖交磁束 [Wb]。この関数は解析的に求まらないため、FEMで磁束-電流-角度の3次元マップ $\psi(\theta, I)$ を作成し、数値微分でトルクを求めるのが実務の標準アプローチとなる。
さらに、SRMの各相の電気回路方程式は以下の通り：
右辺第2項が自己インダクタンス項（$L \cdot dI/dt$に相当）、第3項が逆起電力項（$e = \omega \cdot d\psi/d\theta$、回転速度に比例）である。
インダクタンスプロファイル L(θ,I)
$L(\theta, I)$ のマップってどんな形になるんですか？ イメージが掴めなくて…
低電流ではきれいな台形波に近い形をしている。ロータ突極がステータ突極と向かい合う位置（整列位置、aligned）で最大 $L_{\max}$、突極がずれた位置（非整列、unaligned）で最小 $L_{\min}$ になる。
ところが電流が大きくなると、整列位置付近で鉄心が磁気飽和して $L_{\max}$ が大幅に低下する。一方 $L_{\min}$ は空気のギャップで磁気回路が支配されるからほとんど変わらない。つまり飽和すると $dL/d\theta$ が小さくなってトルクが出にくくなる。これがSRM設計の最大の悩みだ。
| パラメータ | 典型値（8/6 SRM, 1kW級） | 備考 |
|---|---|---|
| $L_{\max}$（低電流） | 40〜80 mH | 整列位置、線形領域 |
| $L_{\max}$（定格電流） | 15〜30 mH | 飽和により大幅低下 |
| $L_{\min}$ | 5〜10 mH | 非整列位置、電流依存性小 |
| $L_{\max}/L_{\min}$比 | 3〜8（線形）/ 2〜4（飽和） | この比が大きいほど高トルク |
磁気飽和と非線形B-Hカーブ
磁気飽和って具体的にどの程度の磁束密度で起きるんですか？ 材料によって違う？
使用する電磁鋼板によって全然違う。例えばSRM頻出の材料でいうと：
- 35A300（日本製薄板）: 飽和磁束密度 $B_s \approx 1.7$ T。鉄損少ないが飽和早め
- 50A470: $B_s \approx 1.8$ T。コスト安だが鉄損やや大きい
- 10JNEX900（JFEの6.5%Si鋼）: $B_s \approx 1.5$ T だが高周波鉄損が極端に小さい。高速SRM向き
実務では1.5〜1.8 T付近からB-Hカーブが寝始め、$\mu_r$（比透磁率）が数千から数十に急落する。FEMではこの非線形B-Hカーブを忠実に入力することが精度の要だ。
SRMに特有の磁気飽和の影響
- 局所飽和：ステータ突極の根元（ヨーク接合部）とロータ突極先端に磁束が集中し、局所的に2.0 T以上に達することがある。この部分のメッシュが粗いと、トルク計算に大きな誤差が生じる。
- 相互インダクタンスの顕在化：飽和状態では磁気回路が非線形になるため、隣接相の電流が通電相のインダクタンスに影響を与える。線形モデルでは無視できた相互インダクタンスがトルク計算に影響する。
- 鉄損増大：飽和領域では磁束密度の時間変化 $dB/dt$ が急峻になり、渦電流損失・ヒステリシス損失ともに急増する。高速回転時には鉄損が銅損を上回ることもある。
トルクリプルとNVH（騒音振動）
SRMの騒音って具体的にどういうメカニズムで発生するんですか？ PMSMより本当にそんなにうるさいんですか？
うるさい。体感で「ギャーン」とか「ビーッ」という金属的な音が出る。メカニズムは明快で、ラジアル方向の電磁力がステータを楕円に変形させるんだ。
SRMでは相の通電ON/OFFでラジアル力が急激に変化する。この力がステータ鉄心の固有振動モード（特に楕円モード = 呼吸モード）を加振する。加振力の周波数成分がステータの固有振動数と一致すると共振し、大きな振動・騒音になる。FEMでやるべきことは：
- 電磁界解析でラジアル力の時間波形を求める
- FFTで加振力のスペクトルを分析する
- 構造モード解析でステータの固有振動数を求める
- 加振周波数と固有振動数の一致（共振）を回避する設計にする
なるほど、電磁界と構造の連成解析が必要ってことですね。かなり大変そう…
その通り。だからSRMのNVH解析は「電磁界FEM → ラジアル力抽出 → 構造FEM（モード解析）→ 音響解析」という多段階のワークフローになる。JMAGやMaxwellには構造解析への力データエクスポート機能があるから、それを活用するのが一般的だ。
SRMは「磁石ゼロのEVモータ」として注目されている
SRMは永久磁石を一切使わないため、中国のレアアース輸出規制リスクと完全に無縁だ。ロータは積層鋼板を打ち抜いて重ねるだけなので、製造コストはPMSMの1/2〜1/3と言われる。Nidecは2023年に家庭用洗濯機向けSRMを量産開始し、AIベースの電流波形最適化で従来SRM比60%の騒音低減を実現した。一方、Land RoverのJaguar I-PACEの初期プロトタイプではSRMが検討されたが、NVH要件を満たせずPMSMに戻った経緯がある。「静かなSRM」は今なおモータ業界のホットトピックで、電磁界FEMと制御シミュレーションの精度が製品化の鍵を握っている。
数値解法と実装
電磁界FEMの定式化
SRMの電磁界解析って、基本はマクスウェル方程式をFEMで解くってことですよね？ でも具体的にどう定式化するんですか？
SRMの解析周波数帯（〜数十kHz）では変位電流を無視できるから、準静的近似（渦電流方程式）を使う。2D解析ではベクトルポテンシャル $A_z$ の1成分だけの問題になる：
ここで $\nu = 1/\mu$ はリラクティビティ（透磁率の逆数）、$\sigma$ は導電率、$\mathbf{J}_s$ は印加電流密度である。$\nu$ が $B$（すなわち $A$）に依存するため、方程式は非線形になる。
2Dで十分なんですか？ 3Dは必要ない？
SRMの基本的なトルク・インダクタンス計算は2D-FEMで十分な精度が出る。軸方向に構造が一様だからね。ただし以下の場合は3Dが必要になる：
- スキュー（ねじり）構造のロータ — トルクリプル低減のためにロータをねじった設計
- 端部効果の評価 — コイルエンドの漏れインダクタンスや端部渦電流
- 軸方向の磁束成分が無視できないケース — アキシャルフラックス型SRM
FEMの弱形式への変換では、ガレルキン法を適用し、三角形要素（2D）または四面体/プリズム要素（3D）で空間を離散化する。SRMの電磁界FEMで特に重要な要素技術を以下にまとめる。
| 要素技術 | 用途 | SRMでの重要度 |
|---|---|---|
| 節点要素（スカラー $A_z$） | 2D断面解析 | 最も基本的。計算コスト低 |
| 辺要素（Nedelec要素） | 3D解析（ベクトル $\mathbf{A}$） | スプリアスモード回避に必須 |
| 2次要素（6節点三角形） | エアギャップの磁束精度向上 | トルク精度に直結。推奨 |
| スライディングメッシュ | ロータ回転の取り扱い | 必須（後述） |
非線形ニュートン・ラフソン法
非線形のB-Hカーブを扱うには、各時間ステップでニュートン・ラフソン法を回すんですか？ 収束は大丈夫ですか？
その通り。各時間ステップの中でニュートン・ラフソン反復を行う。残差方程式はこうだ：
ここで $[S(\mathbf{A})]$ は非線形剛性マトリクス（$\nu(B)$ に依存）、$[M]$ は質量マトリクス（渦電流項）、$\{\mathbf{F}\}$ は電流源ベクトルである。ヤコビアン行列は：
収束のコツは3つある：
- B-Hカーブのスプライン補間を滑らかにする — 折れ線近似だと $d\nu/dB^2$ が不連続になりヤコビアンが振動する
- 初期値に前ステップの解を使う — ゼロ初期値から始めると収束が遅い
- ダンピング（減衰係数）を導入する — $\Delta \mathbf{A}^{(k+1)} = \alpha \cdot [J]^{-1}\mathbf{R}$ で $\alpha=0.5\sim0.8$ 程度
通常3〜8回の反復で $\|R\|/\|R_0\| < 10^{-4}$ に収束する。それ以上かかるならB-Hデータかメッシュに問題がある。
時間ステッピングと回転運動の取り扱い
ロータが回転するのに、メッシュはどうするんですか？ 毎ステップ作り直し？
いい質問だね。実務で使われる手法は主に2つある：
1. スライディングメッシュ法（Band法）
エアギャップ内に円周方向の「スライディングバンド」を設け、ロータ側メッシュとステータ側メッシュを独立に作成する。回転ごとにバンド上の節点対応関係だけを更新する。JMAGとAnsys Maxwellの標準手法。メッシュの再生成は不要だが、バンド上の要素サイズと回転角度ステップが整合していないと、補間誤差がトルクに乗る。
2. エアギャップ・リメッシュ法
各ステップでエアギャップ領域だけを再メッシュする。スライディングバンドより一般的だが計算コストが高い。COMSOLで使われることがある。
時間ステップ幅 $\Delta t$ の決め方は、SRMでは回転速度と通電角度から以下で目安を決める：
$\theta_{\text{step}}$ は1相の通電角度 [rad]、$\omega$ は角速度 [rad/s]。例えば6000 rpmの8/6 SRMでは $\theta_{\text{step}} = 15°$、$\omega = 200\pi$ rad/s なので $\Delta t \leq 4.2 \times 10^{-6}$ s（約4 μs）程度が必要になる。
回路連成解析
実際のSRMはインバータで駆動するんですよね？ FEMだけだと電流波形がわからないのでは？
鋭い。SRMの電磁界FEMを本気でやるなら、駆動回路（インバータ）との連成は避けて通れない。典型的なSRMの非対称ハーフブリッジ回路では、上下のスイッチ素子とフリーホイーリングダイオードの状態によって端子電圧が $+V_{dc}$、$0$、$-V_{dc}$ と切り替わる。
連成のアプローチは2つ：
- 直接連成：FEMソルバー内に回路方程式を組み込む。JMAG・Maxwellともに標準で対応。電圧源駆動でFEMが電流を計算。最も精度が高い
- テーブルベース連成：FEMで $\psi(\theta, I)$ マップを事前計算し、Simulinkの制御モデルに埋め込む。計算は圧倒的に速いが、相互インダクタンスや過渡的な飽和効果を近似している
初期設計にはテーブルベース、詳細設計・NVH解析には直接連成、という使い分けが現場では一般的だ。
直接FEM連成 vs テーブル連成の使い分け
直接連成は「フルオーケストラのリハーサル」のようなもの——全パートが同時に演奏するから最も忠実だが、1回のリハに膨大な時間がかかる。テーブル連成は「パート譜を個別録音してDAWで合成」——各パートは正確だが、即興的な相互作用（飽和時の相間干渉など）は再現できない。設計段階に応じて使い分けるのが実務の知恵だ。
実践ガイド
SRM解析の全体フロー
SRMの電磁界解析を最初から最後まで通すと、どんな流れになりますか？ 全体像を知りたいです。
SRMの設計解析は、だいたいこういうステップで進む：
- 形状パラメータの決定 — ステータ/ロータの極数（例: 8/6, 12/8）、突極幅、エアギャップ長、スタック長、コイル巻数
- 2D断面モデル作成 — CADまたはFEMツール内のパラメトリックモデラーで形状定義
- 材料設定 — B-Hカーブ（非線形）、鉄損係数（Steinmetz定数）、銅の抵抗率
- メッシュ生成 — エアギャップと突極先端を最優先で細分化
- 静磁場解析（まず簡単なもの） — 各ロータ角度で $L(\theta, I)$ マップを作成
- 過渡解析 — 回路連成で実際の電流波形・トルク波形・鉄損を計算
- 後処理 — トルクリプル率、効率、ラジアル力スペクトルを評価
- パラメトリックスタディ / 最適化 — 突極幅・巻数・制御角度の最適化
いきなり過渡解析から始めちゃダメなんですか？
ダメ。静磁場解析で $L(\theta, I)$ マップを取って、まず手計算やMATLABで「だいたいの性能」を確認するのが鉄則。いきなり過渡解析を回すと1ケースに何時間もかかって、間違いに気づくのが遅れる。現場で多いのは「静磁場で2日検証 → 過渡解析で3日 → パラスタで1週間」というスケジュール感だ。
メッシュ戦略（エアギャップ・突極）
SRMのメッシュで特に注意すべき場所はどこですか？
SRMのメッシュでミスると一番痛いのはエアギャップだ。ここの磁束密度分布がトルク計算の精度を直接左右する。目安はこうだ：
| 領域 | 推奨要素サイズ | 理由 |
|---|---|---|
| エアギャップ（半径方向） | ギャップ長の1/3〜1/5（例: 0.5mmギャップなら0.1〜0.17mm） | マクスウェル応力テンソルの精度確保 |
| エアギャップ（周方向） | 突極先端の円弧を40分割以上 | ロータ回転に伴うトルク波形の滑らかさ |
| 突極先端（ステータ/ロータ） | 0.2〜0.5 mm | 局所飽和の正確な捕捉 |
| ヨーク部 | 1〜3 mm | 磁束密度勾配が緩やかなので粗くてOK |
| コイル領域 | 1〜2 mm | 電流密度は均一なので粗くてよい |
| 外部空気領域 | 5〜20 mm（漸増） | 遠方では精度要求が低い |
外部空気領域ってどこまで広げればいいんですか？ 無限遠って言っても…
実務ではステータ外径の2〜3倍の半径まで空気領域を取り、最外周にゼロ・ディリクレ条件 $A_z = 0$ を課すのが一般的だ。それより小さいと「見えない壁」で磁束が反射してトルクがおかしくなる。もっとエレガントにやりたければ、Kelvin変換（無限要素）を使う手もある。JMAGもMaxwellもこの機能を持っている。
境界条件と対称性の活用
8/6 SRMだと、全周モデルを作るのは無駄が多い気がするんですけど、対称性って使えますか？
使える。8/6 SRMではステータ極数とロータ極数の最大公約数 = GCD(8,6) = 2 だから、1/2モデル（180°分）が最小対称単位になる。ただし注意点がある：
- 周期境界条件（反対称/同型対称）を正しく設定すること。SRMでは一般に反対称（アンチペリオディック）条件 $A_z(\theta) = -A_z(\theta + 180°)$ を使う
- 相互インダクタンスの影響を見たい場合は1/2モデルで十分だが、ロータの偏心や製造偏差を考慮する場合はフルモデルが必要
- 12/8 SRMなら GCD(12,8) = 4 で、1/4モデル（90°分）まで減らせる
後処理とトルク計算の検証
FEMでトルクを計算する方法って何種類かあるって聞いたんですけど、どれを使えばいいですか？
主に3つの方法がある。信頼性を確保するために少なくとも2つの方法で一致を確認するのが鉄則だ：
| 方法 | 原理 | 精度 | SRMでの注意点 |
|---|---|---|---|
| マクスウェル応力テンソル法 | エアギャップ中の面積分 | メッシュ依存性が高い | 積分面をギャップ中心に設定。メッシュが粗いと大きな誤差 |
| 仮想仕事の原理 | コエナジーの角度微分 | 高い | 2つのロータ角度での解が必要。摂動量 $\delta\theta$ の選び方に注意 |
| エイリング・トンクス法 | 要素レベルでの局所力計算 | 非常に高い | JMAG標準。メッシュ依存性が小さく推奨 |
マクスウェル応力テンソル法がメッシュに敏感なのは困りますね…
そう。だからメッシュ収束性の確認は必須だ。トルクの値がメッシュを細かくしても変わらなくなる（収束する）まで、少なくとも3段階のメッシュ密度で計算するのがプロの仕事。特にエアギャップの要素層数を2層→4層→6層と変えてトルクの変化が1%以内に収まれば合格だ。
SRMの騒音問題 — 「ガラガラ音」をどう消すか
SRMを実際に回すと、切り替えのたびに「ギャーン」「バコバコ」と耳に刺さる電磁騒音が発生する。これがSRMのEV搭載を阻んできた最大の障壁だ。騒音の元は相切り替え時のラジアル電磁力によるステータの楕円変形（呼吸モード）。対策として (1) 極数を増やして加振次数を上げる（12/8→16/12等）、(2) 電流波形を台形→正弦波風にソフトスイッチングする、(3) ステータ外周にリブを追加して固有振動数をシフトする、などがある。実践では複数の対策を組み合わせ、FEMで加振力スペクトルを予測しながら最適解を探るプロセスが不可欠になる。
初心者が陥りやすい落とし穴 — 空気領域の省略
「空気領域？ なんで空気をメッシュで切るの？」 — 初めて電磁界解析に触れた人がほぼ全員抱く疑問だ。答えは「磁力線は鉄心の外にも広がるから」。特にSRMはエアギャップの磁束が全てなので、ここのメッシュを省略したり粗くしたりするとトルクが数十%ずれる。もう一つ多いのが「B-Hカーブを線形（一定μ）で入れてしまう」ミス。SRMは定格付近で確実に飽和するので、非線形B-Hカーブの入力は絶対に省略してはいけない。
ソフトウェア比較
SRM設計対応ツール比較
SRMの解析に使えるソフトってどのくらいあるんですか？ やっぱり高いんですか？
主要なところだと4つ。それぞれSRM解析での得意・不得意がかなり違う。特にSRMは「回路連成」と「非線形収束の安定性」がツール選びのポイントになる：
| ツール | 開発元 | SRM向けの強み | 弱点 | ライセンス |
|---|---|---|---|---|
| JMAG-Designer | JSOL（日本） | SRMテンプレート搭載、非線形収束が安定、鉄損モデルが豊富、NVH連携が充実 | 海外文献はMaxwellが多い | ノードロック/フローティング（年額300万〜） |
| Ansys Maxwell | Ansys（米国） | Simplorer/Twin Builderとの制御連成が強力、適応メッシュ自動生成 | 初期メッシュが粗い場合トルクリプルに誤差 | Ansysスイート内（年額500万〜） |
| COMSOL Multiphysics | COMSOL AB（スウェーデン） | 電磁-構造-音響の完全連成が1ツール内で可能。カスタマイズ性が高い | モータ専用テンプレートがない。回路連成がやや煩雑 | モジュール制（AC/DCモジュール必須、年額200万〜） |
| Motor-CAD | Ansys/Motor Design Ltd（英国） | SRM含む各種モータの設計特化。高速パラメトリックスタディ | FEMは内蔵だが詳細メッシュ制御不可。NVH非対応 | Ansysスイートまたは単体（年額150万〜） |
オープンソースの選択肢はないんですか？ 大学の研究室だと予算が厳しくて…
あるよ。FEMM（Finite Element Method Magnetics）は2D磁界解析に特化した無料ツールで、SRMの静磁場解析・インダクタンス計算なら十分使える。Lua/Pythonスクリプトで自動化もできる。ただし過渡解析や回路連成はないから、自分でMATLABやPythonで外部ループを組む必要がある。もう一つ、Elmer FEM（フィンランドCSC開発）は3D電磁界の過渡解析も可能なオープンソースだけど、SRM向けのスライディングメッシュの実装がまだ成熟していない。
用途別ツール選定ガイド
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| 大学研究（論文レベル） | FEMM + Python | 無料、スクリプト駆動で再現性確保 |
| 初期設計・パラスタ（産業） | Motor-CAD | 数分で100ケース回せる高速性 |
| 詳細設計（トルク精度重視） | JMAG / Maxwell | 非線形FEM+回路連成の信頼性 |
| NVH解析（騒音対策） | JMAG → Nastran/Actran連携 | ラジアル力→構造モード→音響の連成 |
| マルチフィジクス（発熱も） | COMSOL | 電磁-熱-構造を1ツールで完結 |
| 制御最適化 | Maxwell + Simplorer / JMAG-RT + Simulink | 制御ループとFEMの連成 |
SRM制御シミュレーションに「Simulink + FEM連携」が欠かせない理由
SRMの設計では電磁界FEMだけでは不十分で、電力変換器（インバータ）と制御アルゴリズムを含めたシステムレベルのシミュレーションが必要になる。なぜなら、相電流波形の最適化が騒音・効率の両方を同時に左右するからだ。JMAGの「JMAG-RT」やAnsys Maxwellの「Simplorer」はMatlab/Simulinkとの連携インターフェースを持っており、FEMで計算した$\psi(\theta, I)$テーブルをSimulinkの制御モデルに埋め込むワークフローが普及している。この「電磁界FEM + 制御シミュレーション連成」がSRMの制御最適化には今や標準的なアプローチだ。
トラブルシューティング
非線形収束の失敗
先生、SRMの過渡解析を回したらニュートン・ラフソンが収束しなくて止まりました… 何が悪いんでしょう？
SRMで非線形収束が失敗する原因の95%は以下の4つだ：
| 原因 | 症状 | 対策 |
|---|---|---|
| B-Hカーブの不連続・非単調 | 残差が振動して収束しない | B-Hデータを3次スプラインで滑らかに補間。$dB/dH > 0$ を全域で確認 |
| 時間ステップが大きすぎる | 電流が急変するタイミングで発散 | $\Delta t$ を半分にして再実行。特に通電ON/OFF時に自動ステップ制御を有効化 |
| エアギャップメッシュが粗すぎる | トルクが暴れる、負値が出る | ギャップ半径方向に最低3層の要素を配置 |
| 初期磁化状態の不整合 | 最初の数ステップで発散 | 無通電状態（$I=0$）から徐々に電流を立ち上げる「ランプアップ」を挿入 |
B-Hカーブって、メーカーのデータシートをそのまま入れちゃダメなんですか？
データシートの点数が少ない（10点程度）場合はそのまま入れると折れ線になって微分が不連続になる。最低30〜50点に細分化して滑らかなカーブにするのが安全だ。あと、高磁場側（$B > 2.0$ T）のデータがないことも多い。その場合は $\mu_0$ で外挿するのが標準的なやり方。JMAGにはB-Hカーブの自動外挿・平滑化機能があるから、それを活用するといい。
トルク波形の異常振動
計算は収束したんですけど、トルク波形がギザギザで全然滑らかじゃないんです。これって正しい結果ですか？
SRMはもともとトルクリプルが大きいけど、「物理的なリプル」と「数値的なノイズ」を見分ける必要がある。見分け方は簡単で、メッシュを2倍に細かくして再計算する。ギザギザが減ればメッシュ起因、変わらなければ物理的リプルだ。
数値ノイズの主因は：
- マクスウェル応力テンソル法のメッシュ依存性 — 仮想仕事法やエイリング・トンクス法に切り替える
- スライディングバンドの要素不整合 — バンド上の要素サイズと回転角度ステップを揃える。回転角度 $\Delta\theta$ が要素の1つ分に一致するのが理想
- 時間ステップが粗い — 電気角で0.5°以下のステップにする
初心者が陥る典型的ミス
先生が今まで見てきた中で、SRM解析の初心者が一番やりがちなミスって何ですか？
ベスト5を挙げるとこんな感じだ：
| # | ミスの内容 | 結果 | 正しいやり方 |
|---|---|---|---|
| 1 | B-Hカーブを線形（一定$\mu_r$）で入力 | トルクが実測の2〜3倍になる | 必ず非線形B-Hカーブを使用 |
| 2 | 空気領域が狭すぎる（ステータ外径ぎりぎり） | 磁束が外周で反射しトルクが異常値に | 外径の2〜3倍まで空気領域を確保 |
| 3 | 通電角度（turn-on/turn-off角）の設定ミス | 負トルクが発生し平均トルクが激減 | $dL/d\theta > 0$ の区間でのみ通電 |
| 4 | 対称境界条件を間違える（同型⇔反対称） | 磁束分布が物理的にありえないパターンに | 通電パターンに応じて対称/反対称を使い分け |
| 5 | 鉄損を無視して効率を過大評価 | 実機テストで効率が10〜20%低い | Steinmetzモデルまたは鉄損FEM（JMAGのiron lossソルバー等）を使用 |
SRMは構造が単純な分、シミュレーションの「設定の正しさ」がダイレクトに精度に効く。まずはFEMMで8/6 SRMの静磁場解析を動かしてみるところから始めてみるといい。手を動かすと一気に理解が深まるよ。
「解析が合わない」と思ったら — SRM専用デバッグチェックリスト
- B-Hカーブの確認 — 非線形データが正しく入っているか。単位（A/m vs kA/m、T vs mT）は合っているか
- メッシュ収束性 — エアギャップのメッシュを倍にしてトルクが変わらないことを確認
- トルク計算法の比較 — マクスウェル応力テンソル法と仮想仕事法の結果が5%以内で一致するか
- エネルギーバランス — 入力電力 = 銅損 + 鉄損 + 機械出力 + 蓄積エネルギー変化 が成り立つか
- 既知問題との比較 — 文献の12/8 SRMベンチマーク問題（例: TU Darmstadtのベンチマーク）で自分の手法を検証
なった
詳しく
報告