一个空心电抗器，一块磁铁矿陶瓷：当电、热、力同时登场，单一物理场仿真为何不够用

电气设备在小型化，功率密度在攀升，先进复合材料不断进入工程现场。工程师已经很难再把电磁、传热、流体和结构力学分开看待——多物理场建模（耦合现象仿真）正在成为设计的基本范式。

本文用两个反差很大的案例，说明 CST Studio Suite 与 SIMULIA Abaqus 如何把这条“电 → 热 → 力”的链条完整算下来。

耦合背后的物理链条

要正确设置并解读仿真结果，先要理解这些现象遵循的基本定律：麦克斯韦方程组、傅里叶–毕奥热传导方程，以及广义胡克定律。

低频电磁现象建立在法拉第电磁感应定律、安培定律和高斯定律之上：时变磁场感生电场，流动的电流产生磁通。在有损介质中，电磁能会不可逆地转化为焦耳热，这是通往热力学的桥梁。热量按热传导规律扩散，温度梯度造成材料体积热膨胀，并在内部引起局部和整体应力场。

更关键的是，这些过程互相依赖：温度改变材料属性，又反过来影响电磁场分布，最终形成非线性的强耦合动力系统。

案例一：1 kA 空心电抗器的热–电磁耦合

为什么空心电抗器值得细算

空心电抗器没有铁磁芯，是高压电力系统、高次谐波滤波电路和直流换流站中的关键部件。没有饱和磁芯，电感随电流的线性度很好，对抑制短路电流很有利。

代价是磁通无法被约束在闭合磁路中，会在线圈周围空间大量弥散：对外电磁辐射增强，绕组内部的电流密度分布也变得极其复杂。

第一步：电磁场——抓住趋肤与邻近效应

仿真使用 CST 的低频频域求解器（LF Frequency Domain），基于有限元法对工作空间离散。它能准确捕捉两个决定大电流交流导体焦耳损耗的关键效应：

在这台大电流电抗器中，两种效应以高度非线性的方式叠加。计算得到的磁通密度矢量分布如图 1 所示，最大值如预期出现在线圈几何轴线的中央区域，向外扩展出典型的偶极场形态。

图 1　电抗器内磁通密度矢量分布（CST）

但从设计角度看，更重要的是电流密度的空间分布（图 2）：电流集中在最靠近对称轴的导体边缘。这一分布直接决定了功率损耗——焦耳损耗等于电流密度矢量与电场强度的标量积。

图 2　电抗器导线中的电流密度矢量分布（CST）

由此得到的体积功率损耗（图 3）就是热分析的输入数据——一张明确的热源分布图，需要导入热求解器。

图 3　电抗器线圈中的焦耳功率损耗分布（CST）

第二步：热与流体——不用经验换热系数

在 CST Multiphysics Studio 中，电磁求解器与热求解器之间的数据传递是自动、无损的，避免了人工插值带来的误差。

稳态温度的求解没有采用基于经验换热系数的简化方法，而是使用共轭传热（CHT）求解器：同时求解铜导体内部的热扩散方程，以及导体内外空气流动的完整纳维–斯托克斯方程。由于没有强制冷却，分析的是纯浮力驱动的自然对流，并用布辛尼斯克（Boussinesq）近似处理线圈表面受热引起的空气密度变化，把温度场与动量守恒方程耦合起来。

结果：最热的地方并不在你以为的位置

稳态温度分布给出了一个反直觉的结论：温度最高的区域不在外层线匝，也不在线圈最底部，而是集中在内部的中间层（图 4）。

图 4　所研究电抗器的温度分布（CST）

原因在于气流：较冷的环境空气从下方被吸入线匝形成的通道，一路上升一路被强烈加热，于是线圈上部被明显升温的空气包围，散热能力大幅下降，出现热节流。这些带有精确温度梯度的结果，对绝缘设计非常有价值：可以准确定位关键热点，并优化绕组饼间的间距来提升自然冷却效率。

案例二：磁铁矿棒加热陶瓷，直到它开裂

第一个案例是稳态问题，几何中没有铁磁元件。第二个案例则进入了热冲击下脆性材料的非线性断裂力学：热–电磁–力学耦合分析。

模型中，磁铁矿（magnetite）制成的细铁磁棒嵌入大体积陶瓷基体，再由外部高频电磁场对其循环加热。这类复合结构可用于催化反应器、体积加热系统和航空航天领域。

热只产生在铁磁棒里

磁铁矿磁导率高，与其他铁氧化物相比导电性也较好，在时变磁场中是出色的能量吸收体。陶瓷基体具有介电特性，对电磁波“透明”，因此发热只发生在铁磁夹杂物内部。

这恰恰是把变压器、电抗器中的“寄生”现象——涡流损耗和磁畴翻转引起的磁滞损耗——有意加以利用。

电磁部分：极端的趋肤效应

电磁仿真同样使用频域求解器，对包围复合体的感应线圈建模，并考虑磁铁矿的等效磁化特性。计算表明，时变磁场强烈集中在铁磁棒内部，绕过绝缘陶瓷（图 5）。

图 5　含磁性棒系统中的磁通密度矢量分布（CST）

高激励频率叠加磁铁矿的材料特性，造成极端的趋肤效应：涡流几乎只在棒的近表面薄层内感应产生（图 6）。电流密度和体积损耗在这里达到极值，棒芯部分的直接发热则可以忽略（图 7）。这部分分析的结果，是一张局部功率密度极高的热源分布图。

图 6　棒内电流密度矢量分布（CST）

图 7　棒内焦耳功率损耗分布（CST）

跨软件的场映射：能量必须守恒

与电抗器案例不同，要评估系统的时域响应、温度演化和最终的力学失效，需要把数据转入专门的固体力学环境——SIMULIA Abaqus。

不同软件网格之间的联合仿真与场映射是复杂的数值问题，必须严格遵守能量守恒。棒区域内精细致密的电磁网格上的标量和矢量场，需要插值到 Abaqus 的六面体网格上，同时还要考虑两个程序可能采用不同的形函数插值方法。

热：悬殊的升温速度

热载荷施加完成后，在 Abaqus 中启动瞬态分析，在时域内耦合求解热传导方程和弹性力学方程。

这种复合材料的加热过程极不均衡：磁铁矿迅速体积升温，而陶瓷基体只能靠缓慢的热传导吸热。陶瓷热导率极低，两种材料的交界面上因此形成极陡的温度梯度，热量来不及向基体深处扩散，出现局部过热区。加热过程中的温度分布见图 8。

图 8　棒加热过程中的温度分布（Abaqus）

力：热膨胀失配撑裂陶瓷

快速的温度变化直接成为力学方程的载荷。关键的驱动因素是两种材料热膨胀系数的巨大差异：热膨胀系数较高的磁铁矿迅速膨胀，却受到体积上更冷、更硬的陶瓷基体的约束，形成典型而强烈的热载荷接触问题。

膨胀的棒对周围陶瓷施加巨大的径向压力，按连续介质力学规律，在陶瓷中产生很大的环向拉应力。结构陶瓷抗等静压缩的能力很强，但缺少通过塑性变形松弛应力的机制，单轴或双轴载荷下的抗拉强度很低。基体结构中的微小缺陷会立即成为应力集中点；一旦超过线弹性断裂力学定义的临界应力强度因子，裂纹便开始萌生（图 9）。

图 9　陶瓷中的机械应力分布与开裂（Abaqus）

用 XFEM 追踪热疲劳损伤

仿真的目的不只是确认许用应力被超过——那种做法已经过时，信息量也不足——而是完整分析反复加热、冷却循环造成的材料退化过程。这种现象称为脆性材料的热疲劳，需要先进的损伤扩展算法。

Abaqus 中使用了扩展有限元法（XFEM）：通过在裂纹尖端附近用跳跃函数和渐近函数增强标准形函数，裂纹的萌生与扩展可以不受网格单元边界的限制。

外磁场的周期性开关，在基体中产生脉动应力。每个循环中，应变能在磁铁矿–陶瓷界面应力梯度最大的区域积累，并以微裂纹的形式释放。后处理结果呈现出清晰的损伤图案：裂纹网络从铁磁棒周边萌生，并迅速沿径向向陶瓷试样外缘扩展。冷却过程中材料的多次微观楔入效应进一步加剧了这一过程，最终复合材料完全丧失结构完整性，灾难性地碎裂成多块。

小结：孤立地看物理场，会得出危险的结论

两个案例清楚表明，把物理现象孤立处理，会得到失真、甚至危险的工程结论：

结合 CST Studio Suite 与 SIMULIA Abaqus 的分析能力，再加上 TECHNIA 专家团队的工程经验，科研人员和工业工程师可以在设计阶段就预测产品行为，远远超出传统计算的范围。把断裂力学与电磁现象结合起来，也为智能材料、混合结构以及能源和航空航天领域的新型部件打开了新的设计路径——在这些领域，减重必须与更高的工作载荷同步实现。

耦合现象仿真不再只是可选的验证工具，而是前沿技术在探索、优化和验证阶段的基础环节，能够保证可靠性与安全性，并降低传统样机试制的成本。