トルクリプル解析 — 電動モータのトルク脈動メカニズムと低減手法

ざっくり言うと、モータが出すトルクが回転角度によって「波打つ」現象だね。理想的なモータは一定トルクを出すはずだけど、実際には磁石の配置やコイルへの通電パターンの影響で、トルクが周期的に増減する。

EVで言えば、トルクリプルは車内の不快な振動・騒音の主原因になるんだ。ステアリング用の電動パワステだと、ハンドルが「ゴリゴリ」する感触として運転者に伝わる。産業用ロボットだと位置精度の悪化に直結する。だから、トルクリプルの大きさを数値で予測して、設計段階で低減することがモータ開発の最重要課題の一つなんだよ。

全然違うよ！ コギングトルクは無通電でも発生する磁気的な脈動で、永久磁石とスロット開口部の相互作用だけが原因だ。手でモータの軸を回したときに「カクカク」する感触、あれがコギングトルクだね。

一方、トルクリプルは通電時に発生する総合的なトルク変動だ。コギングトルクに加えて、電流高調波による起磁力の歪み、インバータのPWMスイッチングに起因するリプル電流、さらにリラクタンストルクの脈動成分まで全部含む。つまり、コギングトルクはトルクリプルの一成分に過ぎないんだ。

うまい表現だね。整理すると：

もちろん。トルク波形をフーリエ級数で展開すれば、何次の高調波がどれだけ含まれているかが分かるんだ。3相モータの場合、対称性から電気角6次の整数倍（6次、12次、18次…）が支配的になる。

もっとも広く使われるのがMaxwell応力テンソル法だ。エアギャップ内の閉曲面 $S$ 上で、磁束密度の接線成分 $B_t$ と法線成分 $B_n$ から電磁力を積分する。2D解析の場合、トルクは：

あるよ。もう一つ重要なのが仮想仕事の原理に基づく方法だ。ロータを微小角度 $\Delta\theta$ だけ回転させたときのエネルギー変化からトルクを求める：

大きく4つの発生源がある。それぞれ独立に発生して重畳するから、解析ではどれが支配的かを見極めることが重要だ。

その通り。実際のEVモータ開発では、コギングトルクは全トルクリプルの20〜30%程度にすぎないことが多い。残りはMMF高調波と飽和効果が占める。だからFEM解析では、必ず通電状態の過渡解析でトルクリプルを評価する必要があるんだ。

基本的な流れは3ステップだ。まずステップ1：モータの2D断面モデルを作り、エアギャップ周辺を十分に細かくメッシュ分割する。ステップ2：ロータを微小角度ずつ回転させながら、各ステップで磁場方程式を非線形反復で解く。ステップ3：エアギャップ上でMaxwell応力テンソルを積分してトルクを算出し、電気角1周期分のトルク波形を取得する。

そうだ。FFT分析で6次、12次…のトルク高調波振幅を特定し、どの成分が支配的かを判断する。これが低減策の方向性を決める根拠になる。

最低でも電気角1周期（360°）を120ステップ以上で計算する必要がある。6次高調波をきちんと捕らえるにはナイキストの定理から最低12ステップ/周期だが、精度を考えると20倍の余裕を持たせて240ステップ/周期が実務的な推奨値だ。

また、定常状態に達するまでに数周期の過渡期間が必要だから、解析は3〜5周期分を計算して、最後の1〜2周期のデータをFFT分析に使う。8極モータで3000 rpmの場合、電気角1周期 = 機械角の $360°/4 = 90°$ に相当するから、1周期の所要時間は5 ms程度になるよ。

トルクリプル解析で最も重要なメッシュ設計ポイントはエアギャップだ。Maxwell応力テンソル法では積分経路がエアギャップ内を通るから、ここのメッシュが粗いと結果が全く信用できない。具体的には：

そうだね。だからトルクリプル解析では、対称性を活用して計算領域を削減することが重要なんだ。8極48スロットのモータなら $1/8$ 周期対称で、断面の $45°$ 分だけ計算すれば済む。これでメッシュ数を大幅に減らせる。

鉄心の透磁率 $\mu$ は磁束密度 $B$ に依存するから、FEMではNewton-Raphson法による非線形反復が必要になる。各時間ステップで $B$-$H$ カーブを参照して $\mu(B)$ を更新し、収束するまで繰り返す。珪素鋼板（35H300など）のB-Hカーブは材料メーカーから測定データとして提供されるよ。

注意点として、B-Hカーブのデータ点が粗すぎると、補間誤差がトルクリプルの精度に直接影響する。特に膝の部分（$B = 1.2〜1.8$ T付近）は細かくデータ点を入れること。最低でも30点以上のデータで定義するのが実務の鉄則だ。

実務的なワークフローはこんな感じだ：

設計側の低減手法は大きく分けて5つある：

そうなんだ。スキューの効果はスキュー係数 $k_{\text{sk}}$ で定量化できる：

あるよ。高調波電流注入法だ。FEM解析でトルクリプルの6次成分が支配的だと分かったら、その成分を打ち消す逆位相の6次電流成分をインバータから意図的に流す。電磁界FEMで予測した高調波トルクの振幅と位相をコントローラに組み込むんだ。

まさにその通り！ ただし、電流を注入すると損失が増えるから、設計側の対策と制御側の対策を組み合わせて使うのが実務では最も効果的だ。設計で大まかに低減して、残りを制御で微調整する——「ハードとソフトの合わせ技」が最強なんだよ。

メッシュと材料データが適切なら、平均トルクは実測の±3%以内、トルクリプル率は±20〜30%程度の精度で予測できる。リプル率の精度が低いのは、製造公差（磁石寸法のばらつき、鉄心の打ち抜き歪み）がリプルに大きく影響するからだ。

実測に使う機器はトルクメーター（ロータリートルクセンサー）で、サンプリングレートは回転数の100倍以上が必要だ。3000 rpmなら50 Hz×100 = 5 kHz以上。さらにトルクメーター自体の固有振動数が十分に高いことを確認すること。安い汎用品だと高次のトルク高調波が共振して実際より大きく測れてしまう。

モータ専用の電磁場解析ツールは主に4つあるけど、トルクリプル解析に関してはJMAGとAnsys Maxwellの2強だ。自動車OEM・モータメーカーの大半がこの2つを使っている。COMSOLはマルチフィジクス連成が必要なときに選ばれるケースが多い。

JMAGはJSOL社が開発した電磁場解析ツールで、モータ・発電機の設計に特化している。トルクリプル解析で特に優れているのは：

Ansysスイートとの統合が最大の強みだ。Maxwell（電磁場）→ Mechanical（構造振動）→ Fluent（冷却流体）というマルチフィジクスのワークフローが一貫して組める。NVH（騒音・振動・ハーシュネス）解析まで含めたトータルな評価が必要なときに選ばれる。

トルクリプルから電磁加振力を算出し、それを構造解析に渡してモータケースの振動モードを予測する——こういった「電磁×構造×音響」の連成解析は、EVの車室内騒音の主要な予測手法になっている。

使えるツールはあるよ。FEMM（Finite Element Method Magnetics）は2Dの静磁場・過渡解析ができるフリーソフトで、Luaスクリプトでパラメトリックスタディもできる。Elmer FEMはフィンランドのCSCが開発した汎用FEMで、電磁場モジュールがある。ただし、どちらも商用ツールに比べるとモータ専用機能（テンプレート、スライドメッシュの自動化、FFT後処理）が弱い。学術研究には十分だけど、製品開発には厳しいね。

よくある話だね。解析＜実測の場合、以下を疑うべきだ：

あるよ。その場合はメッシュが粗すぎてスプリアスなトルク高調波が発生している可能性が高い。エアギャップメッシュを2倍に細分化して再計算してみて。あと、トルクメーターのサンプリングレートが不足して高調波成分を捕捉できていない可能性もある。

メッシュ収束性チェックは必須だ。以下の手順で確認する：

高負荷では鉄心の深い飽和領域に入るから、B-Hカーブの非線形性が急激になって収束が難しくなる。対策としては：

必ずしもエラーではないよ。低負荷・高トルクリプル条件では、瞬時トルクが負になることは物理的にあり得る。ただし以下のケースは設定ミスの可能性が高い：