磁気トルク計算
理論と物理
トルクの基本原理
先生、回転機のトルクはどう定式化されるんですか？
磁気トルクはマクスウェル応力テンソルの角度成分から求まる。エアギャップにおけるトルク：
$r$: エアギャップ半径、$L_{stk}$: 積厚、$B_r$: 半径方向磁束密度、$B_\theta$: 接線方向磁束密度。
トルクは$B_r$と$B_\theta$の積なんですね。
そう。永久磁石同期モータ(PMSM)ではトルクをマグネットトルクとリラクタンストルクに分離できる：
第1項がマグネットトルク、第2項がリラクタンストルク。IPMモータでは$L_d \neq L_q$でリラクタンストルクも利用する。
まとめ
マクスウェル応力テンソル——「磁場が面に加える力」の物理
電磁トルクをFEMで計算するとき「マクスウェル応力テンソル法」が最もよく使われる。磁場が境界面に加える機械的応力（法線成分：吸引力、接線成分：せん断力）を面積分することでトルク・力を求める。理論的にはマクスウェルの電磁場理論（1865年）から直接導かれ、任意形状の物体に適用できる汎用性が高い。ただし積分面の位置が空気中のメッシュ密度に依存するため、メッシュが粗いとトルク計算精度が低下する。JMAG・ANSYSはこの問題に対して「平均マクスウェル法」や「仮想仕事法」と組み合わせる改良実装を採用している。
各項の物理的意味
- 電場項 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：ファラデーの電磁誘導法則。時間変動する磁束密度が起電力を生じさせる。【日常の例】自転車のダイナモ（発電機）は、磁石を回転させることで近くのコイルに電圧が発生する——磁場が時間的に変化すると電場が誘起されるというこの法則の直接的応用。IHクッキングヒーターも同じ原理で、高周波磁場の変化が鍋底に渦電流を誘起し、ジュール熱で加熱する。
- 磁場項 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：アンペア-マクスウェルの法則。電流と変位電流が磁場を生成する。【日常の例】電線に電流を流すと周囲に磁場が生じる——これがアンペアの法則。電磁石はこの原理で動作し、コイルに電流を流して強力な磁場を作る。スマートフォンのスピーカーも、電流→磁場→振動板の力というこの法則の応用。高周波（GHz帯のアンテナ等）では変位電流 $\partial D/\partial t$ が無視できなくなり、電磁波の放射を記述する。
- ガウスの法則 $\nabla \cdot \mathbf{D} = \rho_v$：電荷が電束の発散源であることを示す。【日常の例】下敷きで髪の毛をこすると静電気で髪が逆立つ——帯電した下敷き（電荷）から電気力線が放射状に広がり、軽い髪の毛に力を及ぼす。コンデンサ（キャパシタ）の設計では、電極間の電場分布をこの法則で計算する。ESD（静電気放電）対策もガウスの法則に基づく電場解析が基盤。
- 磁束保存 $\nabla \cdot \mathbf{B} = 0$：磁気単極子が存在しないことを表す。【日常の例】棒磁石を半分に割っても、N極だけ・S極だけの磁石は作れない——必ずN極とS極がペアで存在する。これは磁力線が「始点も終点もない閉じたループ」を描くことを意味する。数値解析では、この条件を満たすためにベクトルポテンシャル $\mathbf{B} = \nabla \times \mathbf{A}$ という定式化を用い、磁束保存を自動的に保証する。
仮定条件と適用限界
- 線形材料仮定：透磁率・誘電率が磁場・電場強度に依存しない（飽和領域では非線形B-Hカーブが必要）
- 準静的近似（低周波）：変位電流項を無視可能（$\omega \varepsilon \ll \sigma$）。渦電流解析で一般的
- 2D仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
FEMでのトルク算出
FEMでトルクを計算する手法の比較を教えてください。
| 手法 | 精度 | メッシュ依存性 | 備考 |
|---|---|---|---|
| マクスウェル応力法 | 中 | 高 | 積分面の位置に依存 |
| Arkkio法 | 高 | 中 | エアギャップ体積積分で平均化 |
| 仮想仕事法 | 最高 | 低 | 2回の計算が必要 |
| ノーダルフォース法 | 高 | 中 | 構造連成に便利 |
実務ではArkkio法が多いですか？
JMAGではArkkio法が標準。Ansys MaxwellではVirtual Workがデフォルト。トルクの空間次数分析（FFT）でコギングトルクやトルクリプルの次数成分を分離することも重要。
まとめ
仮想仕事法——「エネルギー微分」によるトルク計算の代替手法
仮想仕事法はトルクを「磁場エネルギーの変位微分」として計算する手法で、マクスウェル応力法の代替として使われる。小さな変位δを与えたときのエネルギー変化δWからトルク=δW/δθを計算する。長所はメッシュ品質への感度がマクスウェル法より低いこと、欠点は2回の計算（変位前後）が必要なことだ。実装上の注意点は変位量δが大きすぎると非線形誤差が生じ、小さすぎると数値誤差が増大すること。最適なδは概ね1電気角度の0.01〜0.1倍程度で、モータ解析の指針として定着している。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務でのトルク評価
モータの性能評価ではN-T特性（速度-トルク特性）の作成が基本。
実務チェックリスト
「計算トルクと実測が10%違う」——電磁トルク計算精度の落とし穴
モータの電磁トルク計算結果が実測値より10〜15%高い場合、モデル上で考慮されていない損失成分がある可能性が高い。主な候補は①機械損（ベアリング損失・風損）、②追加鉄損（端板や締め付けボルトの渦電流）、③巻線端部の3D効果だ。2Dシミュレーションでは巻線端部の磁場が計算されず、端部インダクタンスが過小評価される。精度向上のために2D+端部等価回路、または3D全体モデル解析が推奨される。JMAGの「端部補正機能」は2Dモデルに端部の等価インピーダンスを自動追加して精度を高める。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| JMAG | トルク計算の自動化。N-T特性マップ生成。コギングトルク解析 |
| Ansys Maxwell | 仮想仕事法によるトルク。パラメトリック最適化連携 |
| MotorCAD | 電磁-熱-機械の統合モータ設計。効率マップ作成 |
| SPEED (Siemens) | 解析式ベースの高速トルク計算。初期設計に最適 |
電磁トルク計算ツール——JMAG vs Motor-CAD vs ANSYS Maxwell
電磁トルク計算の主要ツール比較：JMAGはコギングトルク・トルクリップルの高精度FEM計算と鉄損計算の充実で、日本の自動車・電機業界でデファクト。Motor-CAD（Ansys製）は電磁-熱-NVHの統合評価をルミテッドパラメータで高速実行でき、概念設計フェーズの多案検討が得意。ANSYS MaxwellはFEMに加えて電磁-回路-機械の三連成（Motion Solver）でリニアアクチュエータや可変リラクタンスモータの動特性計算に強みを持つ。OpenFOAMとElmer FEMのオープンソース組み合わせは研究機関での無料解析環境として普及中だ。
選定で最も重要な3つの問い
- 「何を解くか」：磁気トルク計算に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術
コギングトルクの精密計算——スロットハーモニクスと計算精度の限界
コギングトルクはスロットと磁石の相互作用で生じる周期的な無励磁トルク変動で、1電気角周期を「最小公倍数」の1周期とする。コギングトルクのFEM計算精度は「回転角ステップ」に強く依存し、角度刻みが粗いとピーク値が過小評価される。実用的な精度（5%以内）には1°以下のステップが必要で、36スロット6極モータでは最低1周期60ステップが目安だ。平均化フラックス法やFourier変換法で高調波成分を抽出し、スキューや磁石配置最適化の効果をシミュレートするのが設計の定石だ。
トラブルシューティング
トラブル
「トルクリップルが予測より大きい」——解析設定と製造誤差の相互作用
モータのトルクリップルが設計シミュレーションの2〜3倍になる場合、磁石着磁の非均一性・スロット形状の製造誤差・電流波形の高調波が主な原因だ。設計シミュレーションは理想的な正弦波電流と完全な幾何対称を前提とするが、実機では必ずしもそうではない。トラブルシューティングの手順：①実機電流波形をフーリエ解析し高調波成分を特定、②その高調波電流をFEMに入力してトルクを再計算、③差異を「製造誤差」「電流波形」「磁石ばらつき」に分離する。JMAG-RTモデルをMATLAB/Simulinkに取り込んだシステム連成解析が診断効率を大幅に上げる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——磁気トルク計算の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告