Abaqus 中充气轮胎稳态滚动的应力应变建模

基于 CST Studio Suite 与 SIMULIA Abaqus 的两个耦合场案例解读电磁场涡流 · 趋肤 温度场焦耳热 · 对流 应力场热膨胀 · 断裂 电气设备在小型化，功率密度在攀升，先进复合材料不断进入工程现场。工程师已经很难再把电磁、传热、流体和结构力学分开看待——多物理场建模（耦合现象仿真）正在成为设计的基本范式。本文用两个反差很大的案例，说明 CST Studio Suite 与 SIMULIA Abaqus 如何把这条“电 → 热 → 力”的链条完整算下来。本文两个案例案例一：额定电流幅值 1 kA 的空心电抗器，在 CST 中完成热–电磁全耦合分析。案例二：嵌入陶瓷基体的磁铁矿棒，经高频磁场循环加热后的热–电磁–力学分析，电磁部分由 CST 完成，热与断裂部分在 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 专家团队的工程经验，科研人员和工业工程师可以在设计阶段就预测产品行为，远远超出传统计算的范围。把断裂力学与电磁现象结合起来，也为智能材料、混合结构以及能源和航空航天领域的新型部件打开了新的设计路径——在这些领域，减重必须与更高的工作载荷同步实现。耦合现象仿真不再只是可选的验证工具，而是前沿技术在探索、优化和验证阶段的基础环节，能够保证可靠性与安全性，并降低传统样机试制的成本。世界级的数值建模，是当今技术进步所使用的语言。来源：ABAQUS仿真世界