Electromagnetic Field FEM Analysis of Induction Motors

Good question. In an induction motor, torque is generated by eddy currents induced in the secondary side due to slip. In FEM analysis, it's crucial to accurately capture the eddy current distribution in the rotor bars—this is a fundamentally different point from PMSM analysis.

Roughly speaking, the flow is like this:

1. Applying three-phase AC to the stator winding creates a rotating magnetic field (synchronous speed $n_s$)
2. The rotor rotates slower than the rotating field (actual speed $n < n_s$), causing a relative change in magnetic flux in the rotor conductors
3. An electromotive force is induced by Faraday's law, causing eddy currents to flow in the rotor bars (squirrel cage) or rotor windings
4. Torque is generated by the interaction between these eddy currents and the rotating magnetic field

In other words, if the speed difference between the rotating field and the rotor—the slip—is zero, the eddy current is zero, and the torque is zero. This is the defining characteristic of induction motors and also the source of analytical difficulty.

They are number one in market share for industrial use. Over 70% of rotating machinery like pumps, fans, compressors, and conveyors use induction motors. They are cost-effective as they don't use rare earth elements, robust, and easy to maintain. Tesla's Model S also initially used a squirrel-cage induction motor. The demand for electromagnetic field FEM analysis remains very high.

Slip $s$ is the ratio of the difference between synchronous speed $n_s$ and actual speed $n$:

Synchronous speed is determined by the power supply frequency $f$ and the number of pole pairs $p$:

For example, a 4-pole motor supplied with 50Hz power gives $n_s = 1500$ rpm. The slip during rated operation is typically around $s = 0.02 \sim 0.05$, so the speed is roughly 1425 to 1470 rpm.

The torque-slip characteristic derived from the equivalent circuit is:

Here, $V_1$ is the phase voltage, $R_1$ is the stator resistance, $R_2'$ is the rotor resistance referred to the primary side, $X_1, X_2'$ are the respective leakage reactances, and $\omega_s = 2\pi n_s / 60$ is the synchronous angular velocity.

From this equation, you can calculate the starting torque ($s=1$), maximum torque (breakdown torque, $s = R_2'/\sqrt{R_1^2 + (X_1+X_2')^2}$), and rated torque. These serve as reference values for verifying FEM analysis results, so you should always have this equivalent circuit model at hand.

In practice, the traditional method is to obtain them from two experiments: the no-load test and the locked-rotor test. However, with FEM analysis, you can simulate these tests and extract the parameters.

The current standard approach is to simulate these using 2D transient FEM analysis and back-calculate the equivalent circuit parameters from the terminal voltage, current, and power.

For electromagnetic field FEM analysis of induction motors, the starting point is the eddy current equation with the magnetic vector potential $\mathbf{A}$ as the unknown:

Here, $\nu = 1/\mu$ is the magnetic reluctivity (reciprocal of permeability), $\sigma$ is the electrical conductivity, and $\mathbf{J}_s$ is the external current density from the stator winding.

Sharp observation. In 2D analysis (assuming infinite length in the axial direction), with $\mathbf{A} = A_z(x,y,t) \hat{z}$, it reduces to a scalar problem with one component:

The conductivity $\sigma$ of the rotor bars is about $3.5 \times 10^7$ S/m (for aluminum), and the term $\sigma \partial A_z / \partial t$ represents the eddy current in the rotor bars. It's common to assume $\sigma \approx 0$ for the stator and rotor cores (laminated steel suppresses eddy currents), but if you want to evaluate core losses in detail, you either set a finite $\sigma$ for the core or perform post-processing loss separation calculations.

We start from the Galerkin weak form of the 2D eddy current equation. Using a test function $w$:

Here, $\Omega_c$ is the integral only over the conductor region (rotor bars). Approximating $A_z$ with shape functions $N_i$:

Substituting this yields a semi-discretized system of ordinary differential equations:

$[K]$ is the magnetic stiffness matrix (depends on $\nu$), $[M]$ is the mass matrix (depends on $\sigma$, non-zero only in conductor regions), and $\{f\}$ is the right-hand side vector from the external current.

Exactly. Moreover, $[K]$ is nonlinear due to $\nu(B)$ (magnetic saturation), so it needs to be updated every step depending on the value of $B$. This is why nonlinear iteration via the Newton-Raphson method is necessary.

形状関数 $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: 収束失敗 → ステップ設定見直し