モータNVH解析（電磁加振力）
理論と物理
なぜEVでモータ騒音が問題になるのか
EVは静かなはずなのに、モーターの音が問題になるんですか？
エンジン音がないぶん、モーター起因の高周波音（ヒューン、キーン）が逆に目立つ。特に電磁加振力の48次成分が2〜4kHzに来ると不快。テスラModel 3の初期型もこれで問題になった。
48次成分って何ですか？ なんでそんな特定の次数が問題になるんですか？
「次数」というのはモータの回転に対する電磁力の周波数の倍率だ。例えば48スロットのモータでは48次が最も支配的な加振力成分になる。モータが5,000 rpmで回転していたら、48次成分は $48 \times \frac{5000}{60} = 4{,}000\ \text{Hz}$ になる。人の耳が最も敏感な2〜4kHz帯にちょうど入ってくるわけだ。
なるほど、高速道路の巡航速度でちょうど不快な周波数になるわけですね。これってエンジン車だと気にならなかったんですか？
エンジン音は広帯域のブロードバンドノイズだから、特定周波数の純音を「マスク」してくれていた。EVではそのマスキング効果がなくなるから、モーターのトーナル音（純音性騒音）がドライバーの耳に突き刺さるように感じるんだ。日産リーフの初期型でも、加速時の「キーン」音が話題になったことがある。
径方向電磁力の基礎
モーターの中で、どうやってそんな振動力が生まれるんですか？
エアギャップ中の磁束密度 $B_r(\theta, t)$ がステータの歯面に径方向の圧力を与える。これがMaxwell応力テンソルの径方向成分で、以下の式で表される：
ここで $B_r$ はエアギャップの径方向磁束密度、$\mu_0 = 4\pi \times 10^{-7}$ H/m は真空の透磁率。磁束密度の2乗に比例するため、磁束密度に含まれる基本波 $f_1$ の成分同士の積から $2f_1$ の加振力が生まれ、異なる高調波 $f_m$ と $f_n$ の積から $f_m + f_n$ および $|f_m - f_n|$ の成分が発生する。
2乗するから周波数成分が増えるんですね。ギターのディストーションみたいな感じですか？
いい例えだね。まさに非線形処理で倍音が生まれるのと同じ原理だ。実務的に重要なのは、エアギャップ磁束密度を空間高調波と時間高調波の積として展開すること。磁束密度は一般に次のように書ける：
ここで $\nu$ は空間高調波の次数、$\omega_\nu$ は対応する角周波数。$\nu = p$（極対数）が基本波で、$\nu = p \pm kQ_s$（$Q_s$: スロット数、$k$: 正の整数）がスロット起因の高調波だ。
力の次数（Force Order）とモード
その高調波がどうステータの振動につながるのか、もう少し詳しく教えてください。
異なる空間次数 $\nu_1$ と $\nu_2$ の磁束密度成分が掛け合わさると、径方向力の力モード次数（circumferential mode order） $r$ が決まる：
この $r$ がステータの振動モード（周方向の変形パターン）に対応する。$r = 0$ は呼吸モード（ステータが一様に膨張・収縮）、$r = 2$ は楕円モード、$r = 3$ は三角形モード、と続く。
ステータの $r$ 次モードの固有振動数は、薄肉円筒シェルの近似式で見積もれる：
ここで $D_s$ は曲げ剛性、$\rho$ は密度、$h$ はヨーク厚さ、$R$ はステータ内径の代表半径。$r = 0$ と $r = 2$ のモードが最も音響放射効率が高いため、これらのモードと電磁加振力の周波数が一致する運転条件で大きな騒音が発生する。
低いモード次数ほど危険なんですね。$r = 0$ の呼吸モードってそんなに大きいんですか？
$r = 0$ はステータ全体が一様に膨張収縮するから音響放射効率がほぼ100%。実際には $r = 0$ 成分は比較的小さいことが多いけど、$r = 2$ の楕円モードは力のレベルも放射効率もともに大きくて、実務上最も問題になるパターンだ。例えば8極48スロットのIPMSMでは、$\nu_1 = 4$, $\nu_2 = 4$ の組み合わせで $r = 0$ と $r = 8$ が生じ、$\nu_1 = 4$, $\nu_2 = 44$ の組み合わせで $r = 48$ が生じる。
キャンベル線図（Campbell Diagram）
キャンベル線図ってよく聞くんですけど、何を見るための図なんですか？
横軸に回転数、縦軸に周波数をとって、電磁加振力の各次数成分を直線で描く図だ。$n$ 次成分の周波数は $f = n \cdot N / 60$ だから、傾き $n/60$ の直線になる。ステータの固有振動数を水平線で入れれば、交点が共振点になる。
交点が多いと危険ということですか？
そうだね。ただし全ての交点が問題になるわけじゃない。重要なのは、（1）その次数の加振力レベルが大きいか、（2）共振するモードの音響放射効率が高いか、（3）常用回転数域にあるか、の3条件。自動車用では2,000〜10,000 rpmの範囲で2kHz以下に交点がないのが理想だけど、実際にはほぼ不可能だから、交点での力レベルをできるだけ下げる設計を目指す。
構造-音響連成
ステータが振動したとして、それがどうやって車内の騒音になるんですか？
ステータ表面の振動速度 $v_n$ が空気を押し退けて音波を放射する。放射音のパワーは以下で表される：
ここで $\sigma_{\text{rad}}$ は放射効率（モード次数と周波数に依存）、$\rho_0 c_0$ は空気の特性インピーダンス（約415 Pa・s/m）、$S$ はステータ外表面積、$\langle v_n^2 \rangle$ は面平均の法線方向振動速度の2乗平均。
ポイントは放射効率 $\sigma_{\text{rad}}$ だ。低次モード（$r = 0, 2$）ではほぼ1に近いけど、高次モードでは周波数が低いと放射効率が劇的に下がる。つまり同じ振動レベルでも、モード次数によって実際の音の大きさがまったく違う。これが「力のレベルが大きくても音にならない」ケースと「小さい力でもうるさい」ケースを生む理由だ。
ステータの振動がハウジングに伝わって、さらにそこから音が出るルートもありますよね？
その通り。実際の伝達経路はステータ → ハウジング → マウント → 車体 → 車室内空気、という構造伝播経路（Structure-borne path）が主経路だ。ステータから直接空気に放射される経路（Air-borne path）もあるが、通常は構造伝播の方が支配的。だからNVH設計では、マウントのゴム硬度やハウジングの剛性設計も含めた「システムレベル」の評価が不可欠になる。
テスラModel 3の「ヒューン音」とOTA修正の裏側
テスラModel 3の初期型（2017年〜）で有名になった高周波の「ヒューン音」は、IPMSMの48次加振力成分がステータの楕円モード（$r=2$）と中速域で共振したことが原因だとされている。興味深いのは、テスラがOTA（Over-the-Air）ソフトウェアアップデートでインバータのPWMキャリア周波数を変更し、体感的な不快感を軽減したこと。ハードウェアを変えずにソフトで「聞こえ方」を変えるアプローチは、従来のNVH設計の常識を覆すものだった。ただし根本原因である電磁加振力自体は変わっていないため、加速時の特定回転域では依然としてトーナル音が残る。
Maxwell応力テンソルの完全形と接線方向力
- 径方向成分 $\sigma_r = \frac{B_r^2 - B_\theta^2}{2\mu_0}$：NVHの主要加振力。径方向磁束密度 $B_r$ の2乗に支配されるが、接線方向磁束密度 $B_\theta$ も差し引きで寄与する。通常 $B_r \gg B_\theta$ なため $B_r^2 / (2\mu_0)$ で近似される。
- 接線方向成分 $\sigma_\theta = \frac{B_r B_\theta}{\mu_0}$：トルク生成に直結する成分。NVHへの寄与はトルクリプルを通じた回転方向の振動（torsional vibration）として現れる。ギヤを介してドライブシャフトに伝わり、ギヤノイズとして車室内に入ることもある。
仮定条件と適用限界
- 2D断面解析の仮定：端部効果（コイルエンドの軸方向磁場）を無視。軸方向力が問題になる場合（スキュー付きモータ、アキシャルフラックス型）は3D解析が必要。
- 線形材料仮定の破綻：鉄心の磁気飽和により $B$-$H$ カーブが非線形になると、高調波成分が大幅に増加する。特に高負荷時（高トルク時）のNVH解析では非線形解析が必須。
- 剛体ロータ仮定：通常のNVH解析ではロータ変形を無視するが、高速回転ロータでは遠心力による変形がエアギャップ長を変化させ、磁束密度分布に影響する。
- 定常状態仮定：急加減速時の過渡状態では、電流制御の遅れやPWMスイッチングの影響が加わり、定常解析より複雑な加振パターンが生じる。
次元解析：加振力レベルの目安
| 物理量 | 典型値 | 単位 | 備考 |
|---|---|---|---|
| エアギャップ磁束密度 $B_r$ | 0.7〜1.2 | T | IPMSMの通常動作範囲 |
| 径方向力密度 $F_r = B_r^2/(2\mu_0)$ | $2 \times 10^5$〜$6 \times 10^5$ | Pa | エアギャップ面に作用する圧力 |
| ステータ外表面の振動加速度 | 0.1〜10 | m/s² | モード・周波数による |
| 1m距離の音圧レベル | 60〜90 | dB(A) | 共振時はさらに増加 |
数値解法と実装
電磁-構造-音響の連成フロー
モータのNVH解析って、1つのソルバーで全部できるんですか？
基本的には3ステップの弱連成（one-way coupling）だ。それぞれ異なる物理を解くから、通常は別々のソルバーを使い分ける：
- 電磁界FEM（JMAG, Maxwell等）：ロータ回転に伴うエアギャップ磁束密度の時間・空間分布を計算し、ステータ歯面の径方向力を抽出
- 構造FEM（Nastran, Abaqus等）：ステータ+ハウジングのモーダル解析（固有振動数・モード形）、および電磁力を外力とした周波数応答解析
- 音響解析（Virtual.Lab, Actran等）：ステータ/ハウジング表面の振動速度から放射音圧を計算（BEM or FEM）
「弱連成」ってことは、振動が電磁場に影響を及ぼすフィードバックは無視するんですか？
いい質問だ。ステータの振動でエアギャップ長が数μm変化するけど、エアギャップ自体が0.5〜1.0mmあるから、変位の影響は0.1%程度。通常は無視して問題ない。ただし超薄型モータやフレキシブルロータの場合は双方向連成（two-way coupling）が必要になることもある。
電磁界FEMと力の抽出
電磁界FEMで力を計算するとき、Maxwell応力法とVirtual Work法のどっちを使えばいいですか？
NVH解析ではMaxwell応力法（MST）が標準だ。エアギャップ中の閉曲面上で磁束密度からノードごとの力を計算する。Virtual Work法はトルク計算には良いけど、空間分布力を出すのには不向き。
Maxwell応力テンソルからの径方向力密度は：
この力を周方向位置 $\theta$ と時間 $t$ の2次元で離散化し、2D FFTによって空間次数 $r$ と時間高調波次数 $n$ に分解する。結果は $r$-$n$ マップとして可視化され、どの力モード・どの次数が支配的かが一目で分かる。
実務上のポイントは、力のFFTを取るときのサンプリング点数だ。空間方向は最低でもスロット数の4倍（48スロットなら192点以上）、時間方向は電気角1周期あたり最低360ステップ欲しい。これをケチると高次成分がエイリアシングで潰れる。
構造FEMによるモーダル解析
ステータのモーダル解析って、何か特殊なことはありますか？
ステータ単体ではなく、ステータ+ハウジング+ベアリング支持のアセンブリとして解くのが正解だ。ステータ単体のFree-Free固有振動数と、ハウジングに焼嵌めた状態の固有振動数は大きく異なる。例えばステータ単体の楕円モード（$r=2$）が2,500 Hzだったのが、ハウジング焼嵌め後に3,200 Hzに上がることもある。
構造の運動方程式は：
ここで $[M]$, $[C]$, $[K]$ はそれぞれ質量、減衰、剛性マトリクス。$\{F_{\text{em}}(t)\}$ は電磁力のノード力ベクトル。モーダル座標に変換すると、各モード $i$ の応答は：
$\zeta_i$ はモード減衰比（構造体では通常0.01〜0.05）。加振力の空間分布がモード形 $\{\phi_i\}$ と直交していれば、力レベルが大きくても当該モードは励起されない。これをモード参加因子（Modal Participation Factor）で評価する。
音響解析と放射音予測
振動から音圧を計算するのって、FEMでやるんですか？ BEMでやるんですか？
モータ単体の放射音はBEM（境界要素法）が主流だ。外部の無限空間を自動的に扱えるから、音響用のボリュームメッシュを切る必要がない。Helmholtz方程式を境界積分で解く：
ここで $G$ はHelmholtzのGreen関数、$c(\mathbf{r})$ は立体角係数。車室内の音場（閉空間）を解く場合はFEMの方が適している。
車室内の音まで計算しようとすると、相当大規模になりそうですね。
その通り。実務ではモータ単体の放射音パワーを計算して、伝達関数（Transfer Path Analysis: TPA）で車室内音圧を予測するアプローチが一般的だ。モータ→マウント→車体→車室の各伝達経路をFRFで繋ぐ。車体FEMの結果からTPA用の伝達関数を抽出して使う。
次数分析（Order Analysis）
次数分析って、FFTとは違うんですか？
FFTは「時間→周波数」の変換だけど、次数分析は「回転角→次数」の変換だ。回転速度が変化するスイープ運転で使う。角度ベースのリサンプリング（回転に同期した等角度刻みで再サンプリング）をしてからFFTすると、各次数成分がきれいに分離できる。
結果はカラーマップ（横軸：回転数、縦軸：周波数、色：振動/音圧レベル）で表示する。これをスペクトログラムとかウォーターフォール図と呼ぶ。各次数成分が直線で見えるから、「48次がこの回転数で固有振動数と交差している」というのが視覚的に分かる。キャンベル線図の実測版と思ってもらえばいい。
電磁力マッピングのイメージ
電磁力を構造モデルに渡す工程は、「大量の圧力センサーのデータを1つずつ貼り付ける」作業に似ている。電磁FEMのメッシュと構造FEMのメッシュは一般にノード配置が異なるため、エアギャップ面の力分布を構造メッシュのノードに「補間（マッピング）」する必要がある。このとき、力の合力（N単位）が保存されていることを必ずチェックすること。マッピング前後で合力が5%以上ずれていたら、補間設定を見直す必要がある。
実践ガイド
解析ワークフロー（ステップバイステップ）
最初から通して、モータNVH解析の手順を教えてください。何から手をつければいいですか？
実務のフローはこうだ。大きく5ステップに分かれる：
- 2D電磁過渡解析：モータ断面の電磁FEMモデルを作り、ロータを回転させながら数電気角周期分の過渡解析を実行。エアギャップの磁束密度時刻歴を取得。
- 径方向力の抽出と周波数分解：Maxwell応力法でステータ歯面のノード力を計算。2D FFT（空間×時間）で力モード次数と周波数の分解を行い、$r$-$n$ マップを作成。
- ステータ構造のモーダル解析：ステータ+ハウジングの3D構造FEMモデルでモーダル解析。$r = 0, 2, 3, 4$ モードの固有振動数を特定。
- キャンベル線図の作成と共振チェック：電磁力の主要次数線とステータ固有振動数の交点を確認。常用回転域での共振リスクを評価。
- 周波数応答解析と音響予測：問題のある次数成分を外力として構造FEMに入力し、振動応答を計算。表面振動速度からBEMで放射音を予測。
2Dでいいんですか？ 3Dの方が正確じゃないですか？
NVH解析では2D電磁FEMが主流だ。理由は計算時間。3Dだと1運転点あたり数時間〜数十時間かかるけど、2Dなら数分〜数十分で済む。スロットスキューの影響は2D解析を複数スライスで行って重ね合わせる「マルチスライス法」で近似できる。3Dが必要になるのはアキシャルフラックスモータや、コイルエンドの軸方向振動が問題になるケースだ。
メッシュと時間刻みの設定
電磁FEMのメッシュで特に気をつけることはありますか？
NVH解析で最も重要なのはエアギャップ領域のメッシュだ。ここの磁束密度分布から力を計算するから、メッシュが粗いと高調波成分が数値的に消えてしまう。
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| エアギャップ径方向分割数 | 3〜5層 | 最低3層。5層で高次成分が安定 |
| エアギャップ周方向要素数 | スロットピッチあたり6〜10 | 48スロットなら全周で288〜480要素 |
| 時間ステップ | 電気角1°以下 | 6極: 機械角0.33°/step以下 |
| 解析周期 | 最低2電気角周期 | 初期過渡を除去するため |
| 構造メッシュ（ステータ） | 最大要素サイズ 2〜3mm | 歯先とヨーク接合部を細分化 |
電磁力マッピングの実務
電磁FEMから構造FEMへの力の受け渡しって、自動でできるんですか？
ツール間の連携機能が充実してきているから、JMAG→Nastranとか、Maxwell→ANSYS Mechanicalの連成はほぼ自動化されている。ただし気をつけるポイントが3つある：
- メッシュ一致性：電磁メッシュと構造メッシュでステータ歯面のノード配置が異なる場合、補間が必要。最近接ノード法よりバリセントリック補間の方が精度が高い。
- 力の保存チェック：マッピング前後で各方向の合力と合モーメントが一致しているか必ず確認。1%以上のずれがあれば設定を見直す。
- 2D→3D押し出し：2D電磁解析の歯面力を3D構造モデルに適用する際、軸方向に一様分布として押し出す。スキュー付きの場合は軸方向位置ごとに位相をずらす。
データ形式はどうしているんですか？ CSVですか？
ツール間の直接連携が使えればそれがベストだけど、汎用的にはUNVフォーマットやCSVでノード番号・座標・力ベクトルを受け渡すことが多い。JMAGは構造ソルバー向けのforce export機能が充実していて、Nastran用のBDF形式やAbaqus用のINP形式で直接出力できる。Maxwellの場合はACT（Application Customization Toolkit）拡張でMechanicalへのマッピングを自動化する。
NVH低減のための設計指針
解析で問題が見つかった後、どうやってNVHを改善するんですか？ ハードウェアの変更が必要ですか？
設計変更には電磁側と構造側の2つのアプローチがある。まず電磁側：
- スロットスキュー：ステータまたはロータのスロットを軸方向にねじると、特定次数の力成分を低減できる。1スロットピッチ分のスキューで基本スロット次数成分が理論上ゼロになる。ただしトルクも若干低下する（通常2〜5%）。
- 極スロット数の組み合わせ最適化：$r = |p \pm kQ_s|$ で低次モード（$r = 0, 2$）が生じにくい組み合わせを選ぶ。例えば8極48スロットより8極36スロットの方が $r = 2$ 成分が小さい場合がある。
- 磁石形状の最適化：磁石エッジの面取りや偏心配置で磁束密度の高調波を低減。
- PWMキャリア周波数の調整：インバータのスイッチング周波数を可聴帯域外（20kHz以上）に上げると電流高調波起因のNVHが大幅に改善する。SiCインバータの普及で実現しやすくなった。
構造側のアプローチもある：
- ステータヨーク厚の増加：固有振動数を上げて共振回避。ただし質量・コスト増。
- ハウジングのリブ追加：特定モードの固有振動数を狙って上げる。
- マウントの最適化：ゴムマウントの硬度・配置を調整して構造伝播経路を遮断。
- ダンピング材の適用：ステータ外周に制振材を貼付して振動エネルギーを熱に変換。
電磁側と構造側のどっちを先にやるべきですか？
「源流対策」の原則で、まず電磁側で加振力自体を下げるのが第一。加振力が大きいまま構造でなんとかしようとすると、ハウジングが重くなったりコストが嵩む。実務的には、概念設計段階で極スロット数の組み合わせを選び、詳細設計で磁石形状とスキューを最適化し、最後に構造側の微調整という順番が効率的だ。
聞こえ方の評価——dB(A)だけでは足りない
モータの騒音は音圧レベル（dB(A)）だけでは不快度を正しく評価できない。モータ音は純音成分が強いため、同じdB(A)でもエンジンの広帯域音より不快に感じられる。欧州の自動車メーカーでは、ラウドネス（sone）、シャープネス（acum）、トーナリティ（tu）、ラフネス（asper）といった心理音響指標を使って「音質」を定量化するのが標準的なプラクティスだ。特にトーナリティ（純音の突出度）はモータNVHの評価に不可欠で、ISO 1996-2に基づく評価手法が使われる。
ソフトウェア比較
NVH連成ワークフロー比較
結局、モータNVH解析をやるにはどんなツールの組み合わせがあるんですか？
主要な連成ワークフローを整理すると、こうなる：
| ワークフロー | 電磁 | 構造 | 音響 | 特徴 |
|---|---|---|---|---|
| JMAG連成 | JMAG-Designer | Nastran / Abaqus | Actran / Virtual.Lab | 日本の自動車メーカーで主流。JMAG-RTで高速力計算も可能 |
| Ansys統合 | Maxwell 2D/3D | ANSYS Mechanical | ANSYS Sound (旧Actran) | Ansysエコシステム内で完結。ACT拡張で連成自動化 |
| COMSOL一体型 | COMSOL Multiphysics（AC/DC + 構造 + 音響モジュール） | 1ツールで完結。研究用途に強いが大規模は重い | ||
| Motor-CAD簡易 | Motor-CAD（Ansys傘下） | 概念設計向け。解析的手法で高速NVH評価。精度は有限要素に劣る | ||
| Siemens統合 | Simcenter SPEED / Motorsolve | Simcenter Nastran | Simcenter Acoustics | Teamcenter PLM連携に強み |
ツール別の強みと弱み
JMAGとMaxwell、どっちがNVH解析に向いてますか？
それぞれ特徴がある。ざっくりまとめると：
- JMAG-Designer：NVH解析に特化した機能が充実。径方向力の空間高調波分解、キャンベル線図の自動生成、Nastran向けのforce export機能がある。日本語サポートの手厚さも強み。ただしAnsysエコシステムとの連携は弱い。
- Ansys Maxwell：ANSYS Mechanicalへのforce mappingがACT拡張で自動化されている。2025年以降はActran（音響解析）も統合されて、Ansysプラットフォーム内でEnd-to-EndのNVHワークフローが組める。電磁解析の精度自体はJMAGとほぼ同等。
- COMSOL Multiphysics：1ツールで電磁→構造→音響の全てを解けるのが最大の利点。ツール間のデータ変換エラーがない。ただし大規模3Dモデル（100万ノード超）では計算効率がJMAG/Maxwellに劣る。研究用途や小型モータの設計最適化に適する。
- Motor-CAD：解析的手法ベースで数秒〜数分でNVH予測ができる。概念設計での大量パラメータスタディ向き。ただしFEMほどの精度はない。
選定の指針
うちの会社でこれから導入するなら、何を基準に選べばいいですか？
3つの観点で考えるといい：
- 既存の構造解析環境：すでにNastranを使っていればJMAG連成が自然。ANSYS Mechanicalを使っていればMaxwell連成が楽。ツール間のインターフェースコストが最も大きい。
- 解析の目的と頻度：概念設計で大量のケースを回すならMotor-CAD。詳細設計での精密NVH予測ならJMAGかMaxwell。研究開発でカスタム物理を入れたいならCOMSOL。
- チームのスキルと教育コスト：JMAGは日本語ドキュメント・サポートが充実していて立ち上げが速い。COMSOLは「何でもできる」反面、正しい設定に到達するまでの学習曲線が急。
まずは2D電磁＋構造モーダルのキャンベル線図チェックから始めて、問題が見つかった運転点だけ詳細な音響解析に進む、というのが効率的ですね。
その通り。全回転域で詳細なBEM音響解析を回すのは計算コスト的に現実的じゃない。キャンベル線図で「危険な交差点」をまず特定して、その回転数±500rpmの範囲だけ詳細解析をするのが実務のスタンダードだ。
Motor-CADの台頭と「シフトレフト」の流れ
近年のEV開発では「シフトレフト」（設計初期段階での問題発見）の重要性が増している。詳細な有限要素NVH解析には1運転点あたり数時間かかるが、Motor-CADの解析的NVH手法なら数秒で結果が出る。Toyota、BMW、Hyundaiなどが概念設計段階でMotor-CADによるNVHスクリーニングを実施し、有望な設計候補だけを詳細FEM解析に回すワークフローを採用している。2023年のAnsysによるMotor Design Ltd.買収により、Motor-CADとMaxwellの統合がさらに進んでいる。
トラブルシューティング
電磁力が構造モデルに正しく伝わらない
先生、電磁力をNastranにマッピングしたんですけど、振動レベルが実測の1/10しかなくて...
それは典型的な「力のマッピングミス」だ。チェックポイントを順に確認しよう：
- 合力の保存：電磁側の全径方向力の合計（N）と構造側に入力された力の合計が一致しているか。座標変換で符号が反転したり、2Dの力を3Dに押し出す際に軸長を忘れていたりする。
- 力の適用面：力をステータ歯面（エアギャップ側）のノードだけに与えているか。ステータ外周面に与えてしまうと伝達が全く変わる。
- 単位系：電磁ソルバーがN/m（2D線荷重）で出力しているのに、構造ソルバーがN（3D点荷重）として読み込んでいないか。2Dから3Dへの変換で積層長 $L_{\text{stk}}$ を掛ける必要がある。
- 時間/周波数のサンプリング：力の時刻歴をFFTした後に構造に入力する場合、Nyquist周波数以下のデータだけが有効。サンプリングが粗すぎると高調波が消える。
キャンベル線図で共振が見えない
キャンベル線図を作ったんですが、交差点があるのに実測では問題になっていません。逆に、交差点がないのに実測で騒音が出ている回転域があります。
キャンベル線図はあくまでスクリーニングツールで、交差点の「有無」だけでは判断できない。よくあるケースを整理する：
- 交差点あるが問題なし：その次数の力モード $r$ が高次（$r > 6$ 等）で、ステータモードへの励振力が弱いか、音響放射効率が低い場合。モード参加因子が小さければ交差しても振動は小さい。
- 交差点ないのに問題あり：PWMキャリア高調波による加振が原因の可能性。キャリア周波数 $f_c$ とその側帯波 $f_c \pm nf_1$ はモータ回転数に依存しない（または緩やかにしか変化しない）ため、キャンベル線図では水平線に近くなり見落としやすい。インバータの電流波形を含んだ電磁解析をしないと出てこない。
- 温度の影響：運転中にステータ温度が上昇すると、珪素鋼板のヤング率が低下して固有振動数が下がる。冷間でのモーダル解析と実動条件で乖離が生じる。
音圧レベルが実測と合わない
BEMで計算した放射音が実測より15dBも小さいんです。何が原因でしょう？
15dBは振動振幅にして約5.6倍の差だから、単一の原因ではなく複合要因の可能性が高い：
- 構造減衰の過大評価：モーダル減衰比 $\zeta$ を大きく設定しすぎると共振ピークが抑えられる。実測のモーダル減衰比を使っているか確認。ステータの減衰比は通常0.005〜0.02程度。
- ハウジング接合部のモデル化：焼嵌め部の接触剛性やボルト締結部の摩擦減衰をどうモデル化しているか。線形バネで代用していると実際の伝達特性と乖離する。
- 磁歪（Magnetostriction）の無視：珪素鋼板の磁歪変形（磁場による直接的な変形）は通常のMaxwell応力法では計算されない。低周波域では磁歪の寄与が全体の20〜30%に達することがある。
- ベアリングの伝達特性：ベアリングを剛結合として扱っていると、実際にはバネ特性を持つベアリングを通過する際の減衰が反映されない。
高次高調波成分の取りこぼし
96次以上の高調波が解析に出てこないんですが、実測では明らかにあります。
高次成分が消える原因は主に3つ：
- 時間刻みが粗すぎる：電気角1°刻みだと360ステップ/電気周期で、理論的にはNyquist限界の180次まで取れるはずだけど、実際にはそれより遥かに低い次数でエイリアシングの影響が出る。安全のために最高対象次数の5倍の時間分解能を確保する。96次を見たいなら電気角 $360/(96 \times 5) \approx 0.75$°以下の時間刻みが必要。
- 空間メッシュが粗い：96次の空間力成分を解像するには、周方向に最低192要素（Nyquist）、実用的には500要素以上必要。
- B-Hカーブの飽和未考慮：線形材料で解くと磁束密度の高調波が小さくなり、力の高次成分も過小評価される。非線形B-Hカーブを使って鉄心の飽和を正しくモデル化すること。
なるほど、NVH解析って「細かいところをサボると全然合わない」解析なんですね。
そうだ。トルク計算なら基本波だけでもそこそこ合うけど、NVHは高調波の積み重ねだから、細部の精度が結果に直結する。「NVH解析が得意な人」と「電磁解析が得意な人」は必ずしも同じではなくて、NVHではFFT処理やサンプリング理論の知識、構造力学の感覚が欠かせない。電磁屋と構造屋と音響屋のトリプルスキルが求められる、ある意味でCAEの中でも最も総合力が問われる分野だよ。
「解析が合わない」と思ったら——NVH解析のデバッグ手順
- まず力を疑え：電磁力の合力レベルと周波数成分が妥当か、解析的な概算値と比較する。$F_r \approx B^2 S / (2\mu_0)$（$S$: 歯面面積）でオーダー合っているか。
- 次にモーダルを疑え：FEMのステータ固有振動数がハンマリング試験の実測値と合っているか。10%以上ずれていたら接触条件や材料定数を見直す。
- 次に減衰を疑え：共振ピーク高さは減衰比に反比例。$\zeta$ を半分にするとピークが2倍（+6dB）になる。
- 最後に伝達経路を疑え：マウント、ブラケット、ボルト締結のモデル化が現実を反映しているか。
なった
詳しく
報告