鉄心飽和解析
理论与物理
铁心饱和
老师，铁心饱和是一种什么现象？
即使增大磁场强度$H$，磁通密度$B$也几乎不再增加的现象。当铁的磁矩（磁畴）全部排列整齐后，就无法进一步磁化。
B-H曲线
非线性的$B$-$H$关系：
| 区域 | $\mu_r$ | 特征 |
|---|---|---|
| 初始区域 | 100〜500 | 低磁场。可逆 |
| 急剧上升区 | 1000〜10000 | 磁畴壁移动占主导 |
| 饱和附近（膝点） | 10〜100 | 磁畴旋转 |
| 饱和区 | ≈1 | $B \approx \mu_0 H + B_{sat}$ |
饱和后$\mu_r \approx 1$…和空气一样吗？
是的。饱和的铁心不再导磁。磁通会泄漏出去，导致效率下降、损耗增加、发热。
典型的饱和磁通密度
总结
- 饱和 → $\mu_r ≈ 1$ — 铁心不再导磁
- B-H曲线是非线性的 — FEM中需要Newton-Raphson迭代
- 电工钢: $B_{sat} ≈ 2.0$ T — 电机设计的极限
铁心饱和的物理——磁畴结构造就BH曲线的“膝点”
电工钢磁化曲线上出现的“膝点”，标志着磁畴壁移动完成、只剩下磁化旋转阶段的过渡。软钢的Bsat（饱和磁通密度）约为2.0〜2.1 T，铁氧体则为0.3〜0.5 T，差异巨大。变压器或电机铁心设计的基本准则是将工作点保持在膝点以下，但为了轻量化、高输出功率的要求，将B值设定在膝点附近的案例正在增加。将饱和区的非线性BH曲线精确地纳入FEM，决定了铁损和磁通分布的计算精度，而选择Langrange单元还是Nedelec单元也变得至关重要。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度产生电动势。【日常例子】自行车发电机（磨电灯）通过旋转磁铁在附近的线圈中产生电压——这是磁场随时间变化会感应出电场这一原理的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流产生磁场。【日常例子】电线通电后周围产生磁场——这就是安培定律。电磁铁根据此原理工作，通过线圈通电产生强磁场。智能手机的扬声器也应用了此原理：电流→磁场→振膜的力。在高频（如GHz频段天线等）下，位移电流 $\partial D/\partial t$ 不可忽略，它描述了电磁波的辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电通量的发散源。【日常例子】用垫板摩擦头发后，静电会使头发竖立——带电的垫板（电荷）向外辐射出电力线，对轻质的头发施加力。电容器设计中，电极间的电场分布就依据此定律计算。ESD（静电放电）防护也基于高斯定律的电场分析。
- 磁通连续性 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半，也无法得到只有N极或只有S极的磁铁——N极和S极总是成对出现。这意味着磁力线描绘的是“没有起点和终点的闭合回路”。在数值分析中，为了满足此条件，采用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 的公式化方法，自动保证磁通连续性。
假设条件与适用范围
- 线性材料假设：磁导率、介电常数不依赖于磁场、电场强度（饱和区需要非线性B-H曲线）
- 准静态近似（低频）：位移电流项可忽略（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀、可忽略端部效应时有效
- 各向同性假设：各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的情形：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
非线性磁场FEM
$\nu = \nu(|\mathbf{B}|)$（$\nu = 1/\mu$）依赖于磁通密度。用Newton-Raphson迭代求解：
1. 初始估计（用$\mu_r = 1000$等求线性解）
2. 根据$|\mathbf{B}|$从B-H曲线更新$\nu$
3. 用切线磁导率$d\nu/d(B^2)$构成雅可比矩阵（切线刚度）
4. 迭代直至残差为零
和结构的弹塑性分析用的是同一种非线性求解器呢。
正是如此。B-H曲线对应着应力-应变曲线。
B-H曲线的输入
将钢板制造商的数据表（磁化曲线）以表格形式输入。注意：
- 饱和区的延伸 — 如果数据只到$B = 1.8$ T，则按$\mu_0$的斜率延伸
- 插值方法 — 样条或线性插值。急剧的转折点对收敛有不良影响
总结
- Newton-Raphson迭代 — 与结构弹塑性分析相同
- B-H曲线的表格输入 — 注意饱和区的延伸
非线性磁场分析的迭代法——Newton-Raphson法与BH曲线的处理
考虑铁心饱和的磁场分析是非线性问题，Newton-Raphson法（NR法）是标准的迭代解法。NR法收敛速度快（二次收敛，误差平方递减），但缺点是初始值不好会导致发散。在BH曲线的“陡峭部分（膝点附近）”，磁导率变化剧烈，雅可比矩阵的计算精度直接关系到收敛性。实用中常采用先给线性初始解再转入NR的两阶段法，或与Picard迭代（定点法）结合以提高稳定性，这些方法在ANSYS Maxwell等商用软件中也有采用。
边单元（Nedelec单元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM：适用于非线性材料、非均匀介质。BEM：自然处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉夫逊法处理B-H曲线的非线性。残差标准：通常为 $||R||/||R_0|| < 10^{-4}$。
频域分析
通过时间谐波假设归结为稳态问题。需要进行复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中可使用更大的步长，但需注意精度。
频域与时域的使用区分
频域分析类似于“将收音机调谐到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实务
电机、变压器、电抗器的铁心设计。设计时要确保工作点不进入饱和区。
检查清单
- [ ] B-H曲线是否基于实际材料的数据表
- [ ] 是否包含饱和区（$B > 1.8$ T）的数据
- [ ] Newton-Raphson是否收敛（残差在$10^{-4}$以下）
- [ ] 铁心的最大$B$是否超过饱和值（若超过则需修改设计）
- [ ] 是否评估了饱和导致的磁通泄漏影响
“变压器比设计时更热”——局部饱和引发的意外损耗
在变压器或电抗器设计中，当出现“实际样机的温升比计算值高10〜15 K”的问题时，往往是因为忽略了铁心的“局部饱和”。在拐角处、接缝处，磁通集中导致局部饱和，涡流损耗急剧增加。即使截面平均磁通密度低于Bsat，也可能存在局部超出的情况。通过FEM分析查看磁通密度的彩色云图，可以一目了然地发现拐角处的高密度区域。有效的对策包括改变拐角处的叠片方向（斜切）或优化接缝位置，采用FEM预先验证后再制造的流程已成为标准。
分析流程的比喻
电机的电磁场分析，感觉上类似于“给吉他调音”。调整琴弦的粗细（线圈匝数）和琴桥的位置（磁铁布置），以引出最美妙的音色（高效的扭矩特性）。改变一个参数，整体的平衡就会改变——所以参数化研究很重要。
初学者容易掉入的陷阱
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会扩散到铁心之外”。如果将分析区域设置得紧贴铁心，无处可去的磁通就会“撞上”边界并反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙上不断弹跳的状态。
边界条件的思考方式
远场边界条件看似不起眼，但极其重要。需要在数值上表达“从这里开始是无限延伸的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被反射回来。
软件比较
工具
全磁场FEM求解器都支持非线性B-H曲线。JMAG、Maxwell、COMSOL、FEMM都支持饱和分析。JMAG和Maxwell的电工钢材料数据库非常丰富。
铁心饱和分析工具——ANSYS Maxwell vs JMAG vs FEMM
包含铁心饱和的磁场分析主流工具是ANSYS Maxwell、JMAG和FEMM（开源）。Maxwell擅长大规模3D非线性磁场分析和电路耦合，用于变压器、电机的集成设计。JMAG凭借其丰富的日本电工钢BH数据库和铁损计算精度，在日本国内汽车、电机行业占据压倒性份额。FEMM虽然仅限于2D轴对称和平面问题，但免费，广泛应用于教育、研究和简易设计。仿真精度的关键与其说是计算方法，不如说是“BH曲线数据的质量”，这一事实暗示着投资于材料测量比工具选型更重要。
选型时最重要的三个问题
なった
詳しく
報告