コギングトルク解析と低減手法

そう、永久磁石とスロットの間の磁気的な引力で生じる脈動トルクだ。EVの静粛性やロボットの位置決め精度に直結する。

その通り。コイルに電流を流さなくても、永久磁石の磁力線がステータのスロット開口部を通過するとき、磁気抵抗（リラクタンス）が角度によって変わるんだ。磁力線は「通りやすい道」を好むから、磁石がスロット開口に差し掛かると引っかかるような力が生まれる。これがコギングトルクだ。

いい例えだね。実際、EVのステアリングモータでコギングトルクが大きいと、ハンドルを切るたびにドライバーが振動を感じてしまう。産業用ロボットでは位置決め精度が落ちるし、直動モータなら停止精度に直結する。だから「いかにコギングを減らすか」がモータ設計者の腕の見せどころなんだ。

コギングトルクの本質はエアギャップ中の磁気エネルギー $W_{mag}$ が回転角度 $\theta$ に依存して変動することにある。永久磁石が作る起磁力分布と、スロット開口によるパーミアンス（磁気コンダクタンス）分布の「掛け算」で磁束密度が決まるんだ。

ざっくり言うと、パーミアンスはスロット開口部で急激に下がる。磁石がこの「穴」の上を通過するたびにエネルギーが変動し、その角度微分がコギングトルクとして現れるわけだ。

正確にはもう少し複雑で、1つの磁石が複数のスロットに同時にまたがっている場合もある。だからコギングトルクの周期は単純に「スロット数分」ではなく、極数とスロット数の最小公倍数（LCM）で決まるんだ。

コギングトルクは周期関数だから、フーリエ級数で展開するのが自然だ。

大正解。一般的に高調波次数が上がると振幅は急速に減衰する。だからLCMが大きい極スロット組合せほど、コギングトルクは小さくなる。

もう一つ重要な式がマクスウェル応力テンソルによるトルク計算だ。エアギャップ中の閉曲面 $S$ で積分すると：

典型的な組み合わせで比較してみよう。

その通り。LCM=60は1回転あたり60回の微小な脈動に分散されるから、体感的にはほぼゼロに近い。ただし分数スロットには巻線係数が下がるデメリットもあるから、トルク密度とのトレードオフは必要だ。「スロット数と極数を決める段階でコギング問題の7割は決まる」と言われるくらい重要な設計パラメータだよ。

主要な低減手法を体系的にまとめると以下の通りだ。

理論的には1スロットピッチ分（$360°/N_s$）のスキューでコギングの基本波がゼロになる。ただし平均トルクも sinc 関数的に低下するから、実務では半スロットピッチ程度に抑えることが多い。いずれにせよ「どの手法をどこまで適用するか」の定量的な判断にFEM解析が不可欠なんだ。

主に3つの方法がある。それぞれ精度と計算コストにトレードオフがあるんだ。

JMAGは内部で改良型の法線力法を使っていて精度が高い。Ansys Maxwellは仮想仕事法がデフォルトだ。COMSOLではユーザーが選択できる。MST法を使う場合は積分面をエアギャップの中央に設定し、複数の積分面で結果を平均するのがベストプラクティスだ。

モータの電磁場解析では、磁束保存 $\nabla \cdot \mathbf{B} = 0$ を自動的に満たすためにベクトルポテンシャル $\mathbf{A}$ を導入する。$\mathbf{B} = \nabla \times \mathbf{A}$ と定義すれば、2D解析での支配方程式は：

その通り。だからコギングトルク解析は比較的シンプルな静磁場問題になる。ただし鉄心の非線形B-H特性を考慮するとニュートン・ラフソン反復が必要になるし、ロータの回転位置を変えながら何十ステップも計算するから計算時間は意外にかかるんだ。

大きく影響する。ティース先端部は磁束が集中して飽和しやすく、飽和すると透磁率が急激に下がる。飽和した領域は「空気に近い」状態になるから、実質的にスロット開口が広がったような効果が生じてコギングが増大する場合がある。

だからコギングトルクの解析では線形材料では不十分。必ず非線形B-Hカーブを使った解析が必要だ。収束判定は残差ノルム $||R||/||R_0|| < 10^{-4}$ が一般的で、反復回数は通常10〜30回程度で収束する。

いい質問だ。メッシュの再生成は計算コストが大きいので、実務ではエアギャップ中にスライディング境界（回転バンド）を設定する方法が標準的だ。

その通り。特にコギングトルクのような微小量は、この接合処理の精度に敏感だ。JMAGのスライディングメッシュは内部でAZフレーム法を使っていて高精度とされるが、Maxwellではバンド上のメッシュ密度を十分確保しないと数値ノイズが混入することがある。

コギングトルクの周期性を利用すれば計算量を大幅に削減できる。

2Dならそうだね。3Dでスキューを考慮する場合は計算量が一気に増えるから、2Dのスライス積分近似（多スライス法）を使うことも多い。これは軸方向を数枚のスライスに分けて、各スライスの2D解析結果を足し合わせる手法だ。

コギングトルク解析の標準的なワークフローはこうだ：

そう。ただし低減手法を検討する段階では1電気角周期（$360°/p$）分を計算して、コギングトルクとトルクリップルの両方を確認するのが望ましい。コギングだけ減らしてもトルクリップルが増えたら意味がないからね。

コギングトルク解析で最も重要なメッシュ設計のポイントを教えるよ。

そうなんだ。メッシュ収束確認は必須で、少なくとも3水準（例えば3層→5層→7層）でコギングトルクのピーク値が5%以内に収束することを確認する。足りないと数値ノイズがコギングに混入して、本来の値が見えなくなる。

3Dでやるのが最も正確だけど、実務では多スライス2D法がよく使われる。軸方向を $M$ 枚のスライスに分け、各スライスを位相をずらして2D解析し、結果を平均する手法だ。

基本的には同じだね。3段ステップスキューなら $M=3$ でそのまま計算できる。ただしステップスキューでは各段のエッジ効果（軸方向端部の漏れ磁束）が無視できない場合があり、その場合は3D解析が必要になる。

磁石形状のパラメトリック最適化は、コギング低減の最も実践的な手法の一つだ。主なパラメータは：

一般的にコギングトルクは embrace に対してV字型や波打つ曲線を描く。特定の値（例えば 0.833 = 5/6 など分数値）でコギングが極小になる。この最適値はスロット数と極数の組み合わせで変わるから、FEMのパラメトリックスイープで見つけるのが確実だ。実務ではまず0.05刻みの粗い探索で傾向を掴み、有望な範囲を0.01刻みで詰める2段階アプローチが効率的だよ。

コギングトルク解析の妥当性検証には以下のチェックが有効だ。

十分にメッシュを切った2D解析で、実測比10〜30%程度の誤差が一般的だ。3D解析＋端部効果考慮で10%以内に迫れる。ただしコギングトルクは「量産品のバラツキ」も大きくて、磁石の着磁ムラや鉄心の組み付け精度で個体差が20〜50%出ることもある。だからFEM結果は「設計の方向性」の確認に使い、最終的な品質保証は実測で行うのが現実的だよ。

モータ設計で使われる主要な電磁場解析ツールを紹介するよ。

できるよ。ただしFEMMは2D限定でスキュー解析機能がないから、Pythonスクリプトで多スライス法を自前実装する必要がある。学生の研究や初期検討には十分だけど、製品設計にはJMAGやMaxwellの方が効率的だ。

用途別に整理するとこうだ：

これはコギングトルク解析で最も多い相談だね。チェックすべき項目を優先度順に挙げるよ。

変わる。エアギャップ長がパーミアンスに反比例するから、ギャップが0.5mmから0.45mmに変わるだけでコギングが10〜20%増大することもある。実測との照合では「公差範囲でのバラツキ」を必ず感度解析で評価することが重要だ。

数値ノイズが原因の「偽の高調波」だ。主な原因は3つ：

メッシュを2倍に細かくして再計算したときに振幅が大きく変わる成分は数値ノイズ。物理的な成分はメッシュに依存しない。また、理論的にコギングトルクは $n \cdot N_{LCM}$ 次（$n=1,2,3,...$）の成分しか持たないから、それ以外の次数に成分があれば偽物と判断できるよ。

よくあるケースだ。チェックすべきは：

磁石がスロット開口の真上に来た角度で、ティース先端の磁束密度が急変するから、非線形のB-Hカーブ上で「飽和点の乗り越え」が起きてニュートン・ラフソン法が振動しやすくなるんだ。

その通り。コギングトルクは「小さいのに厄介な量」だから、メッシュの詰めと結果の検証を手抜きすると全く意味のない数字が出てくる。でも正しくやれば、設計段階でコギングを最小化して、NVH（振動・騒音）の問題を未然に防げる。実際に手を動かして、まずは2Dの単純モデルから試してみるといいよ。