感应电动机的电磁场FEM分析

问得好。感应电机是通过滑差在次级侧感应出涡流来产生转矩的。FEM分析需要准确捕捉转子导条的涡流分布——这是与PMSM分析根本不同的关键点。

简单来说，流程是这样的：

1. 定子绕组通入三相交流电产生旋转磁场（同步转速 $n_s$）
2. 转子以低于旋转磁场的速度旋转（实际转速 $n < n_s$），导致转子导体产生相对的磁通变化
3. 根据法拉第定律感应出电动势，在转子导条（鼠笼式）或转子绕组中产生涡流
4. 该涡流与旋转磁场相互作用产生转矩

也就是说，旋转磁场与转子速度之差 = 滑差为零时，涡流也为零，转矩也为零。这是感应电机最大的特点，也是分析的难点所在。

在工业领域它可是份额第一哦。泵、风机、压缩机、传送带等旋转机械的七成以上都是感应电机。它不使用稀土，所以成本低，坚固耐用且维护方便。特斯拉的Model S早期也用的是鼠笼式感应电机。电磁场FEM分析的需求依然非常高。

滑差 $s$ 是用比率表示同步转速 $n_s$ 与实际转速 $n$ 之差：

同步转速由电源频率 $f$ 和极对数 $p$ 决定：

例如，给4极电机提供50Hz电源，则 $n_s = 1500$ rpm。额定运行时的滑差通常在 $s = 0.02 \sim 0.05$ 左右，转速大约在1425〜1470 rpm。

由等效电路推导出的转矩-滑差特性如下：

这里 $V_1$ 是相电压，$R_1$ 是定子电阻，$R_2'$ 是折算到一次侧的转子电阻，$X_1, X_2'$ 是各自的漏电抗，$\omega_s = 2\pi n_s / 60$ 是同步角速度。

从这个公式可以计算出起动转矩（$s=1$）、最大转矩（失步转矩，$s = R_2'/\sqrt{R_1^2 + (X_1+X_2')^2}$）和额定转矩。这是验证FEM分析结果的基准值，所以首先必须掌握这个等效电路模型。

实际工作中，传统方法是通过空载试验和堵转试验（转子锁定试验）这两个实验来求取。但如果是FEM分析，则可以通过仿真重现这些试验来提取参数。

通过FEM的2D瞬态分析模拟这些试验，然后根据端电压、电流、功率反算等效电路参数，这是目前标准的做法。

感应电机的电磁场FEM分析，出发点是以磁矢量势 $\mathbf{A}$ 为未知数的涡流方程：

这里 $\nu = 1/\mu$ 是磁阻率（磁导率的倒数），$\sigma$ 是电导率，$\mathbf{J}_s$ 是定子绕组的外部电流密度。

很敏锐。在2D分析（假设轴向无限长）中，设 $\mathbf{A} = A_z(x,y,t) \hat{z}$，则归结为一个分量的标量问题：

转子导条的电导率 $\sigma$ 约为 $3.5 \times 10^7$ S/m（铝的情况），这一项 $\sigma \partial A_z / \partial t$ 就表示转子导条的涡流。定子铁心和转子铁心通常假设 $\sigma \approx 0$（叠片钢抑制涡流），但如果要详细评估铁损，则需要在铁心中也设置有限的 $\sigma$，或者在后处理中进行损耗分离计算。

从2D涡流方程的Galerkin弱形式出发。使用权函数 $w$：

这里 $\Omega_c$ 仅对导体区域（转子导条）积分。用形函数 $N_i$ 近似 $A_z$：

代入后，得到半离散化的常微分方程组：

$[K]$ 是磁刚度矩阵（依赖于 $\nu$），$[M]$ 是质量矩阵（依赖于 $\sigma$，仅在导体区域非零），$\{f\}$ 是由外部电流产生的右端向量。

没错。而且 $[K]$ 由于 $\nu(B)$ 是非线性的（磁饱和），需要根据 $B$ 的值在每一步更新。这就是为什么需要牛顿-拉夫森法进行非线性迭代的原因。

在2D分析中，$A_z$ 是标量，所以使用普通的节点单元（拉格朗日单元）没问题。另一方面，在3D分析中，$\mathbf{A}$ 是矢量场，因此边单元（Nedelec单元）是必需的。

边单元能在单元边界上自动保证磁通密度法向分量的连续性，因此在3D电磁场分析中已成为事实上的标准。JMAG、Maxwell、COMSOL在3D分析中都采用了边单元。

时间步长 $\Delta t$ 的设置直接关系到解的精度。参考标准如下：

• 电气角每度一个步长是基本。50Hz时，1周期 $= 20$ ms，360步/周期，则 $\Delta t \approx 55.6 \, \mu$s
• 若要捕捉槽谐波，则需要更细：$\Delta t = T_e / (360 \times 2 \sim 3)$
• 若是PWM逆变器驱动，则需要开关频率10倍以上的采样：例如10kHz开关频率，则 $\Delta t \leq 10 \, \mu$s

时间积分多采用Crank-Nicolson法（$\theta = 0.5$）。$\theta = 0$ 是前向欧拉法（容易不稳定），$\theta = 1$ 是后向欧拉法（稳定但精度为一阶）。

实际工作中常用的是滑动网格法。在气隙中央设置滑动面（接触面），独立保持定子侧网格和转子侧网格。根据旋转角度重新连接滑动面上的节点。

JMAG、Maxwell、COMSOL都标配了滑动网格法，所以用户几乎不需要显式地编写重划网格程序。不过，为了让滑动面上的单元分割角度与转子旋转角度匹配，将网格分割角度设为旋转步长的整数倍是个技巧。

铁心材料（电工钢片）的B-H曲线具有很强的非线性。例如35A300，在 $B = 1.0$ T附近 $\mu_r \approx 5000$，但超过 $B = 1.8$ T后 $\mu_r$ 会急剧下降到100以下。这就是磁饱和。

在FEM中，每个单元的 $\nu = 1/\mu(B)$ 成为磁通密度的函数，因此每个时间步内都需要进行牛顿-拉夫森法的非线性迭代：

收敛判据通常为 $\| \Delta a \| / \| a \| < 10^{-4}$ 左右。实际工作中，通常3〜8次迭代即可收敛，但当齿尖或轭背进入深度饱和时，迭代次数会增加，最坏情况会发散。此时，增加B-H曲线的数据点（特别是膝点附近更密集）或引入欠松弛（$\omega = 0.3 \sim 0.7$）是常用手段。

主要有三种方法：

从实用性来说，虚功法最为稳定。麦克斯韦应力张量法在网格较粗时会产生较大的振荡成分，所以在对结果进行FFT分析时需要特别注意。

这是最佳实践。如果两种方法的结果相差5%以上，就应该重新审视网格。特别是气隙网格较粗时，麦克斯韦应力张量法的结果会剧烈波动，所以应该先确认网格收敛性，然后再做判断。

典型的分析流程是这样的：

1. 相关主题

   耦合分析誘導加熱解析电磁分析PMSM电机设计分析耦合分析誘電加熱解析电磁分析BLDC电机分析热分析板级热分析

形状関数 $N_i$ を用いて未知量を近似:

これを数式で表すとこうなるよ。

これを数式で表すとこうなるよ。

連続体の支配方程式を離散化すると、以下の代数方程式系が得られる:

ここで $[K]$ は全体剛性マトリクス（または同等のシステムマトリクス）、$\{u\}$ は未知節点変数ベクトル、$\{F\}$ は外力ベクトルなんだ。

• 完全積分: 全ての項を正確に積分。剛性過大評価の傾向（ロッキング）
• 低減積分: 積分点数を削減。計算効率向上だが、アワーグラスモード発生のリスク
• 選択的低減積分 (B-bar法): 体積項と偏差項を分離して積分。ロッキング回避

• h-refinement: メッシュを細分化（要素サイズ h を小さく）して精度向上
• p-refinement: 要素の多項式次数を上げて精度向上
• hp-refinement: h と p を同時に最適化

収束速度: 二次要素で $O(h^2)$ のオーダーで誤差が減少（滑らかな解の場合）

誘導モータ解析の実務的な解析フローと注意点を解説する。

1. 前処理 (Pre-processing)

• CADデータのインポートと形状簡略化
• 材料特性の定義
• メッシュ生成（要素タイプ・サイズの決定）
• 境界条件と荷重条件の設定

2. 求解 (Solving)

• ソルバー設定（解法、収束基準、出力制御）
• ジョブ投入と計算実行
• 収束モニタリング

3. 後処理 (Post-processing)

• 結果の可視化（変位、応力、その他の物理量）
• 結果の検証と妥当性確認
• レポート作成

• 応力集中部: 最低3層以上の要素を配置
• 応力勾配の大きい領域: 要素サイズを周囲の1/3〜1/5に
• 荷重印加点近傍: 局所細分化
• 遠方領域: 粗いメッシュで計算効率を確保

• 過拘束に注意: 剛体移動の拘束は6自由度のみ
• 対称条件の活用: 計算規模の削減
• 荷重の等価分配: 集中荷重 vs. 分布荷重の選択

日本のJSOL Corporationが開発。電気機器設計に特化した電磁場解析ツール。

現在の所属: JSOL Corporation

Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。

現在の所属: Ansys Inc.

• メッシュ収束性を3水準以上で確認したか
• 力の釣り合い（反力合計）を検証したか
• 結果が物理的に妥当な範囲か確認したか
• 既知の理論解またはベンチマーク問題と比較したか

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

誘導モータ解析に対応する主要な商用CAEツールの機能比較と、各製品の歴史的背景を詳述する。

日本のJSOL Corporationが開発。電気機器設計に特化した電磁場解析ツール。

現在の所属: JSOL Corporation

Ansoft Maxwell。低周波電磁場解析。2008年Ansysに統合。

現在の所属: Ansys Inc.

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

現在の所属: COMSOL AB

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

現在の所属: Ansys Inc.

• 要素タイプの非互換: ソルバー固有要素は中立フォーマットで表現不可
• 材料モデルの差異: 同名でも内部実装が異なる場合がある
• 境界条件の再定義: 多くの場合、手動での再設定が必要
• 結果データの比較: 出力変数の定義（節点値 vs. 要素値、積分点値）に差異

誘導モータ解析のツール選定においては以下を考慮:

• 解析規模: 数万〜数億DOFへのスケーラビリティ
• 物理モデル: 必要な構成則・要素タイプの対応状況
• ワークフロー: CADとの連携、自動化の容易さ
• コスト: 初期投資 + 年間保守 + 教育コスト
• サポート: 技術サポートの質とレスポンス

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

誘導モータ解析における最新の研究動向と先進的手法を見ていこう。

• 等幾何解析 (IGA): NURBS基底関数を直接使用し、CAD-CAE間のシームレスな連携を実現
• 粒子法 (SPH, MPM): メッシュフリー手法による大変形・破壊の追跡
• 位相場法 (Phase-Field): 界面の暗示的表現による複雑な界面追跡
• 機械学習支援: サロゲートモデル、物理インフォームドニューラルネットワーク (PINN)

症状: ソルバーが指定反復回数内に収束せず異常終了

考えられる原因:

• メッシュ品質の不足（過度に歪んだ要素）
• 材料パラメータの不適切な設定
• 不適切な初期条件
• 非線形性が強すぎる（荷重ステップの不足）

対策:

• メッシュ品質チェックを実施（アスペクト比、ヤコビアン）
• 材料パラメータの単位系を確認
• 荷重を複数ステップに分割（サブステップ数の増加）
• 収束判定基準の緩和（ただし精度に注意）

症状: 応力/変位/温度等が物理的に非現実的な値

考えられる原因:

• 境界条件の誤設定
• 単位系の混在（SI単位と工学単位の混同）
• 不適切な要素タイプの選択
• 応力特異点の存在

対策:

• 反力の合計を確認（力の釣り合い）
• 単位系の一貫性を確認
• 要素タイプの適切性を再検討
• 特異点除去またはサブモデリング

症状: 計算が想定時間の何倍もかかる

対策:

• メッシュの粗密分布の最適化
• 対称性の活用（1/2, 1/4モデル）
• ソルバー設定の最適化（反復法、前処理の選択）
• 並列計算の活用

症状: Out of Memory エラー

対策:

• アウトオブコア解法の使用
• メッシュ規模の削減
• 64bit版ソルバーの使用確認
• メモリ割り当ての増加

• FATAL 2012: 特異剛性マトリクス → 拘束条件の見直し
• USER WARNING 5291: 要素品質不良 → メッシュ修正
• SYSTEM FATAL 3008: メモリ不足 → MEM設定の調整

• Excessive distortion: 要素の過大変形 → NLGEOM確認、メッシュ改善
• Zero pivot: 拘束不足 → 境界条件追加
• Time increment too small: 収束失敗 → ステップ設定見直し