永磁电机退磁分析

问得好。退磁分为“可逆退磁”和“不可逆退磁”两种。可逆退磁在温度下降后可以恢复。问题在于不可逆退磁。当B-H曲线的工作点降到knee point（拐点）以下时，永磁体的磁力会永久性降低。即使移除外部原因也无法恢复。

简单来说，永磁体内部存在称为“磁畴”的小区域。正常状态下这些磁畴的取向一致，因此能产生强磁力。但是，当暴露在强反向磁场或高温下时，部分磁畴会反转成杂乱的方向。一旦被推压到拐点以下，反转的磁畴就无法恢复原来的方向——这就是不可逆退磁。

最危险的是高温与大电流的组合。具体来说，比如电动汽车爬陡坡时（大扭矩=大电流）、逆变器短路故障时，或者低速高扭矩长时间运行等情况。当线圈铜损导致磁体温度上升到180°C或200°C左右时，矫顽力会急剧下降，工作点会一举跌破拐点。

永磁体的退磁曲线（第二象限的B-H曲线）通常是直线下降，但当超过某个磁场强度时，曲线会急剧弯曲，磁通密度骤降。这个“弯曲角”就是拐点。用公式表示如下：

没错。而且温度升高时，拐点会向右上方移动——也就是说，即使反向磁场较弱也更容易进入不可逆退磁状态。可以想象成高温下磁体变得“脆弱”。因此，退磁分析中必须使用温度依赖的B-H曲线。

NdFeB（钕铁硼）磁体的温度依赖性可以用以下两个公式表示：

典型的温度系数如下：

是的。所以在EV电机设计中，如果磁体温度会上升到150°C这样的条件，就必须选择矫顽力高的SH系或UH系等级。如果为了省钱使用便宜的N系磁体，现实中确实会发生夏天爬坡时退磁导致输出功率骤降的情况。

你注意到了关键点。铁氧体在高温下会变强，但反过来低温下矫顽力会下降。冬季冷启动时有退磁风险。所以铁氧体电机的退磁分析需要考虑“最低温度×最大电流”的情况。这与NdFeB的分析思路正好相反。

量化这个的就是退磁裕度（Demagnetization Margin: DM）。定义如下：

实际工作中的标准如下：

这正是退磁分析的难点所在。假设温度是多少度、电流加到多少安培、如何考虑磁体的离散性（公差），这些都会导致结果大不相同。因此，“最坏情况的叠加”非常重要。

按危险程度从高到低排列如下：

• 逆变器短路故障：流过最大电流5〜10倍的过电流，产生巨大的d轴方向退磁场。是最严苛的条件。
• 低速·高扭矩连续运行：EV爬陡坡、牵引。铜损大，磁体温度急剧上升。
• 弱磁控制（Field Weakening）：高速旋转时施加负向d轴电流以削弱磁通的制御。直接增大退磁场。
• 铁氧体在低温环境：如前所述，-40°C冷启动时。

是的。高速巡航时电机温度也升高了，所以是高温+弱磁的双重打击。实际设计中，需要在最高速度×连续运行×最高冷却液温度这样的条件下确认退磁裕度。

基础是使用磁矢势 $\mathbf{A}$ 的麦克斯韦方程组公式化。对于2D电机截面分析，控制方程如下：

没错。退磁后 $\mathbf{M}_r$ 会变小。但磁体并非整体均匀退磁，而是会发生磁体内不同位置退磁程度不同的“局部退磁”。磁体边缘或角落附近退磁场集中，退磁从那里开始。

FEM会计算每个单元的磁通密度 $B$ 和磁场强度 $H$ 的工作点，并在每个单元判断其是否低于拐点。

使用牛顿-拉弗森法（NR法）进行迭代计算。每次迭代更新磁导率 $\mu$。磁体区域的本构关系可以写成：

很敏锐。实际上，由于拐点附近B-H曲线的斜率急剧变化，NR法的收敛性容易变差。对策包括：

• 使用样条插值平滑连接B-H曲线的拐点附近
• 对NR法应用阻尼（欠松弛：0.3〜0.7）
• 将收敛判据设为残差范数 $\|R\|/\|R_0\| < 10^{-4}$ 左右
• 将磁场增量分割成小步长

要准确评估退磁，需要进行热-电磁耦合分析。流程如下：

1. 通过电磁场分析计算扭矩·损耗分布（铁损+铜损+磁体涡流损）
2. 将损耗分布作为热分析的热源输入，计算温度分布
3. 根据温度分布更新每个单元的磁体B-H曲线
4. 使用更新后的B-H曲线再次进行电磁场分析 → 回到步骤1

重复此过程直到温度收敛。实际工作中，通常2〜5次迭代就足以收敛。

确实很耗时。所以实际工作中常用“单向耦合”。首先通过稳态热分析预测磁体的最高温度，然后仅使用该温度的B-H曲线进行一次电磁场分析。严格的双向耦合通常只在最终确认时进行。

局部退磁建模主要有两种方法：

• 单元剩磁更新法：判断每个有限单元的工作点，仅减少跌破拐点的单元的 $B_r$。JMAG的“退磁分析”功能采用此方式。
• 回复线追踪法：退磁后的磁体沿回复线（斜率 $\mu_0 \mu_{\text{rec}}$）从拐点返回。新回复线的截距即为退磁后的 $B_r'$。

要正确评估局部退磁后的电机特性，需要使用退磁单元的 $B_r'$ 再运行一次电磁场分析。因为退磁前后转矩脉动或齿槽转矩的波形会改变，这也会影响NVH（噪声振动）。

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