感应加热 — CAE术语解说
感应加热
老师，感应加热和IH电磁炉的原理一样吗？CAE也能处理这个问题吗？
感应加热的理论基础
基本概念和支配方程
感应加热是线圈通电时工件被加热的现象，对吧？但为什么没有直接接触却会变热呢？
很好的问题。核心在于"电磁感应"和"涡流"。交流电流流过线圈周围会产生随时间变化的磁场。根据法拉第电磁感应定律，这个变化的磁场在导体（工件）内部产生感应电动势。这个电动势使导体内部形成涡形电流，即"涡流"。这个涡流在导体的电阻中转化为焦耳热，工件本身就会发热。用数学表示，感应电动势是磁通随时间的变化率。
那么，工件的哪个部分会最热呢？表面似乎比内部热的印象。
你的感觉是对的。这是由"表皮效应（skin effect）"引起的。交流电流及其产生的磁场的影响在导体表面从强到弱，向内部急剧衰减。这个衰减的参考深度称为"表皮深度（skin depth）"δ，由频率f、导体的磁导率μ和导电率σ决定。
例如，钢（μ_r≈100，σ=7×10^6 S/m）中通入10kHz电流时，表皮深度约为0.7mm。也就是说，几乎所有的发热都发生在表面以下1mm的区域内，内部通过导热升温。频率越高，加热越集中在表面。
表皮深度的式子中有磁导率μ，听说钢的磁导率随温度变化很大。这对计算有什么影响呢？
这是非常重要的一点。铁钢材料的磁导率在居里温度（约770°C）处急剧下降，变为常温值的百分之一。代入前面的表皮深度公式，加热进行到透磁率下降时，表皮深度δ就会增加。也就是说，加热进行过程中"加热模式本身会改变"这种非线性现象会发生。要在CAE中准确模拟，必须有温度相关的磁导率、比热和电导率数据。
感应加热的数值计算方法
电磁-热耦合分析的离散化
用CAE解感应加热时，需要电磁场分析和热传导分析耦合，对吧？具体要解什么方程，按什么顺序解呢？
完全正确。通常采用弱耦合（loosely coupled）方法。首先使用从麦克斯韦方程组推导出的低频近似（忽略位移电流）方程，比如A-φ法（磁矢量势A和电位φ）来求解涡流分布。
从这里得到的焦耳发热密度
每个时间步都用更新后的温度对应的材料物性（σ、μ）重新计算电磁场，进行迭代计算。
电磁场分析的网格需要考虑前面说的表皮深度来细分吗？
必须要做。在表皮深度δ范围内，至少要放置3～5层以上的单元。比如δ=1mm时，工件表面的网格尺寸应该在0.2～0.3mm左右。更粗的网格无法准确捕捉涡流分布，会低估发热量。相反，工件内部和空气区域不需要那么细。很多求解器都有生成表面边界层网格的功能。
听说耦合计算有时会发散。有什么办法保持数值稳定吗？
主要原因有两个。一是温度上升导致材料物性急剧变化。特别是磁导率在居里点处不连续下降。缓解办法是对物性数据进行平滑处理，或者设置足够小的时间步长（比如0.01秒以下）。二是电磁场和温度场的更新时机。不仅仅用简单的逐次代入法，在每个时间步内进行数次电磁→热的计算"子步进"或选择更加强韧的耦合算法。在Ansys中，"Magnetic Transient"和"Thermal Transient"耦合时的"Coupling Control"设置很重要。
感应加热的实际应用
工作流和验证
实际开始感应加热仿真时，最初应该定义什么？
首先是"明确加热条件"。如果有实机规格书，从中提取频率（比如中频30kHz）、逆变器额定功率（比如100kW）、线圈圈数（5圈）和形状、目标温度（淬火则850°C以上）、加热时间（3秒）。这是所有输入条件的基础。从材料物性数据库（比如日本钢铁联盟的"热物性计算系统"或软件内置库）中准备对应材料（比如S45C）的温度相关数据（比热、热传导率、电导率、磁导率）。透磁率数据缺失是最大的障碍。
线圈建模要注意什么？实际线圈是铜管，里面还有冷却水流动。
实际上，线圈自身的发热和温度升高也不能忽视。建模有两种方法。一种是"实心线圈模型"：定义铜的物性，在线圈截面内设置均匀励磁电流密度。要考虑冷却影响，就在线圈表面应用对流边界条件（水冷约5000 W/m²K）。另一种是"绕线模型"：将线圈简化为通电的"导体"，省去发热计算，改用基于实测的等效冷却条件。后者计算成本低，但在线圈设计阶段应该用前者。
结果的验证怎么做才能保证可信度？和实验比较时有什么要点？
最直接的是用热电偶或红外温度计测温度履历来比较。要点是把温度测量位置和仿真模型上的特定点（节点）严格对应。还必须验证输入电功率。把仿真计算出的线圈阻抗（电阻和电抗）和用LCR表等测实际值比较。比如在100kHz时计算阻抗5mΩ + j20mΩ，实测5.2mΩ + j19mΩ，模型的可信度就很高。还要检查发热量累积值（焦耳积分）和从工件温度升高反推的吸热量之间的平衡。
感应加热软件对比
各软件的特点和应用案例
感应加热分析中经常听到Ansys、COMSOL、JSOL的"JMAG"等，它们在方法上有什么区别吗？
电磁场求解的核心都是有限元法，但在耦合的"便利性"和"专业性"上有差异。COMSOL Multiphysics的最大优势是能在一个界面内直接耦合"电磁热"。用户也可以修改支配方程，适合研发。Ansys用"Maxwell"解电磁场，用"Mechanical"或"Twin Builder"解热传导，在Workbench环境中自动传递数据。适合大规模问题和已有Ansys工作流的集成。JMAG是电磁专用，特别是旋转机类强项，但感应加热也有丰富的模板和材料库（特别是钢材），工业界应用实绩丰富。
OpenFOAM或Elmer等免费或低成本软件能做感应加热分析吗？
可以，但门槛很高。Elmer原生支持电磁热耦合，可以耦合"MagnetoDynamics"求解器和"HeatSolve"求解器。但GUI比商用软件简洁，材料物性温度相关性设置和复杂网格生成需要花时间。对于OpenFOAM，需要用户自己耦合电磁场分析功能扩展（比如"foam-extend"的electroMagnetic模块）和热传导求解器，要求相当高的知识和开发能力。学习基本原理时不错，但设计工作用它就不太现实了。
选商用软件时，从材料库的角度应该注意什么？
这是购买决策的关键之一。感应加热主要对象是钢材，各钢种（碳钢、合金钢、不锈钢）的温度相关数据，特别是"相对磁导率"曲线的丰富度很重要。Ansys可以和Granta MI等数据库联动，但默认库有限。JMAG和日本厂商联系紧密，JIS规格钢材数据较多。COMSOL自有库外，用户易用.csv文件自定义物性表。如果没有数据，还要从厂商目录（比如"大同特殊钢"的技术资料）手动输入。
感应加热问题排查
常见错误和解决方案
计算能跑，但工件的温度明显低于实测。可能的原因是什么？
首先怀疑的是"表皮深度区域网格不足"。如前所述，表皮深度内单元层数不够会低估涡流密度。重新细化网格再计算。其次是"材料物性，特别是常温磁导率设置"。目录上的"相对磁导率"不一定是初始磁导率。比如"相对磁导率500"，那是某个磁场强度下的值。实际感应加热时强磁场作用，饱和状态的磁导率（数十）可能更合适。需要用实测阻抗反推"参数辨识"，调整磁导率。
反过来，温度异常高，上升过快，这是什么问题？
是热损失考虑不周。感应加热中，工件的热传导、对流、放射冷却都不能忽视。特别是高温区（500°C以上）放射传热占优。确认发射率ε是否设置为约0.8（氧化钢板）。还有，如果工件和治具接触，那部分的热传导放热也要建模。另一个原因是输入功率过大估计。是否忽视了线圈电阻发热（铜损）和电源-线圈匹配损耗，而直接把电源功率当成工件发热？要重新审视考虑功率效率的输入条件。
非线性计算经常出现"无法收敛"错误。初始设置中应该审查哪些参数？
改善收敛性，按以下顺序尝试。1. **时间步长**：最初极小（0.001秒），然后自动增大。2. **求解器松弛因子**：非线性强时，放松电磁场求解的收敛判定（比如残差范数从1e-4改为1e-3）。3. **材料数据平滑**：磁导率在居里点附近的急变，用窄温度范围线性插值平滑化。4. **初始条件**：工件初始温度不用室温（20°C），改高点（比如100°C），物性变化平缓，易收敛。Ansys中调整"Step Control"，COMSOL中调整"Time-Dependent Solver"的详细设置。
3D模型计算成本太大。有什么效率化方法？
首先是"利用对称性"。线圈和工件轴对称就降到2D模型。这最有效。不对称但有周期对称就只建其中一部分。再次是"网格自适应"。只把发热区（工件表面）网格极细，空气区和工件内部粗化。还要考虑"频域分析"。不做过渡分析，假设为定常正弦波用复数求解"谐波"分析，求涡流分布，把发热量传给定常热传导分析。虽然不能考虑温度相关物性，但对加热初期或温度上升小的情况有效，计算时间会大幅缩短。
细节
错误