滞后建模
理论与物理
磁滞现象
老师，磁滞现象是指磁化强度跟不上磁场变化的滞后现象吧？
对强磁性体施加交变磁场时，B-H曲线会形成一个回线。这个回线的面积就对应着磁滞损耗。
主要模型：
- Jiles-Atherton (J-A) 模型 — 5参数。基于磁畴壁钉扎的微分方程模型
- Preisach 模型 — 磁滞元的叠加。可再现任意形状的次回线
- Play/Stop 磁滞元模型 — 与机械滞后的相似性
J-A模型是最常用的吗？
由于易于集成到FEM中，J-A模型使用较多。但高精度方面Preisach模型更优。JMAG采用了Play模型。
总结
- 磁滞损耗 — B-H回线的面积
- J-A模型 — 5参数的微分方程模型
- Preisach模型 — 高精度再现次回线
磁滞的物理——磁壁钉扎与“磁铁的记忆”
磁性材料的磁滞回线是“磁壁移动和磁化旋转”对外部磁场响应滞后的现象。磁壁被晶界、夹杂物、缺陷“钉扎”导致磁化变化不可逆，从而产生矫顽力（Hc）。软磁材料（电工钢片）磁壁钉扎小，磁滞回线面积（铁损）小；硬磁材料（钕铁硼磁铁）钉扎大，具有强矫顽力。Preisach（1935）将磁滞现象数学建模为“无限多个开关元件的叠加”，奠定了现代磁滞CAE模型的基础。
各项的物理意义
- 电场项 $\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中
如何将磁滞模型实现到FEM中呢？
计算成本似乎很高啊。
相比普通的B-H曲线（单值）会增加3〜5倍的计算时间。Preisach模型需要在每个积分点保存磁滞元的历史数据，因此内存也会增加。JMAG通过Play模型实现了高效实现。
总结
- 在每个积分点追踪历史 — 每个时间步更新B-H关系
- Newton-Raphson迭代 — 非线性收敛
- 计算成本 — 普通B-H曲线的3〜5倍
磁滞模型的数值实现——Preisach模型与Jiles-Atherton
CAE中的磁滞建模主要有“Preisach模型”和“Jiles-Atherton模型”两种方法。Preisach模型通过逆分析从测量的一次磁化曲线求取密度函数，能精确再现任意磁化过程，但计算成本高。Jiles-Atherton模型用5个参数的微分方程描述磁化，易于集成到商用FEM求解器中，计算速度快。对于磁滞损耗的高精度预测，Preisach模型更优，但在电机设计的实用分析中，Jiles-Atherton模型被广泛采用。
边单元（Nedelec单元）
专用于电磁场分析的单元。自动保证切向分量的连续性，消除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料、非均匀介质。BEM: 自然处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉弗森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 为一般标准。
频域分析
通过时间谐波假设归结为稳态问题。需要进行复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的使用区分
频域分析类似于“调收音机到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——能再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实际应用
主要用于铁损的精密评估、变压器的剩磁分析、磁屏蔽的交变磁场响应。
实践检查清单
- [ ] 是否根据实测B-H回线拟合了磁滞模型的参数
- [ ] 时间步长是否足够细（影响回线的平滑度）
- [ ] 是否需要再现次回线（需要则用Preisach模型）
- [ ] 是否存在直流偏磁情况下的非对称回线
- [ ] 是否考虑了温度依赖性（注意居里温度）
“计算铁损与实测不符”——磁滞模型的校准
使用Simation（Steinmetz公式）计算电工钢板的铁损时，与实测的差异有时可达10〜30%。原因是Steinmetz系数会因环氧处理、冲压加工的影响而变化，且模型中未包含磁通密度的高次谐波成分。提高精度的方法有：①使用从实际设备切取的E型叠片试样的测量数据拟合Jiles-Atherton参数；②引入冲压端面影响的附加损耗系数。JMAG的“磁特性输入辅助工具”具有根据测量数据自动拟合参数的功能。
分析流程的比喻
电机的电磁场分析感觉类似于“给吉他调音”。调整琴弦粗细（线圈匝数）和琴桥位置（磁铁布置），以引出最美妙的音色（高效的扭矩特性）。改变一个参数，整体的平衡就会改变——所以参数化研究很重要。
初学者容易陷入的误区
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会有的疑问。答案是“因为磁力线也会扩散到铁心之外”。如果将分析区域仅限定在铁心附近，无处可去的磁通会“撞上”边界壁并反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球在墙上不断弹跳的状态。
边界条件的思考方式
远场边界条件看似不起眼但极其重要。需要在数值上表达“从这里开始是无限广阔的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被反射回来。
软件比较
工具
| 工具 | 特点 |
|---|---|
| JMAG | 基于Play模型。在铁损分离中活用磁滞模型 |
| Ansys Maxwell | 同时支持J-A模型和Preisach模型 |
| COMSOL AC/DC | 内置J-A模型。也可实现自定义模型 |
| MagNet (Simcenter) | 支持Preisach模型。在变压器涌流分析方面有实绩 |
磁滞分析工具——JMAG vs ANSYS Maxwell vs Opera FEM
作为磁滞FEM分析工具，JMAG（JMAG钢材库丰富）、ANSYS Maxwell（电磁+电路耦合）、Cobham Technical Services Opera FEM（电机/发电机设计老牌厂商）是代表性产品。JMAG以库的形式提供日本各电工钢板制造商的测量数据，在日本汽车、家电行业采用率尤其高。Opera FEM在欧美大型电机、变压器设计方面有实绩，复杂磁滞模型的可定制性高。COMSOL的AC/DC模块最近也强化了磁滞损耗功能，持续吸引新进入的制造商。
选型时最重要的三个问题
- “要解决什么问题”：所需的物理模型、单元类型是否支持磁滞建模。例如，流体方面是否有LES支持，结构方面接触、大变形的处理能力是差异点。
- “谁来使用”：新手团队适合GUI完善的工具，有经验者适合脚本驱动的灵活工具。类似于汽车的自动挡（GUI）和手动挡（脚本）的区别。
- “未来扩展到什么程度”：考虑到未来的分析规模扩大（HPC支持）、向其他部门扩展、与其他工具的联动，这样的选择有助于长期的成本削减。
尖端技术
尖端技术
- 矢量磁滞 — 旋转磁场下的B-H矢量关系。电机铁心铁损评估不可或缺
- 应力相关磁滞 — 机械应力引起的B-H曲线变化（磁致伸缩效应）。冲压
なった
詳しく
報告