2. 电磁场与热场基础参数回顾
电磁热耦合分析的本质,是把“电磁损耗从哪来”和“热量往哪去”这两条链路用同一套材料参数串起来。电磁场求解器关心介电常数 ε、磁导率 μ、电导率 σ,它们决定场分布与损耗密度;热场求解器关心热导率 k、比热容 cp、密度 ρ 和对流换热系数 h,它们决定温升幅度与温度梯度。若两侧参数各自独立取值、忽略温度依赖性,耦合结果往往在高温工况下明显偏离实际。本节先把六个关键参数的定义、物理意义和温度行为逐一厘清,为后续损耗映射与热网络建模打基础。
一、电磁侧三个参数:场如何分布、损耗如何产生
介电常数 ε描述介质在电场作用下的极化能力,宏观上体现为电位移与电场强度的比值。工程中更常用相对介电常数 εr = ε/ε0,其中 ε0 ≈ 8.854×10-12 F/m。对时谐场,复介电常数写作 ε = ε′ − jε″,虚部与实部之比即损耗角正切 tanδ = ε″/ε′。介质损耗密度可写为 pd = ωε0εr′tanδ·E2(单位 W/m³),可见 εr 与 tanδ 直接决定单位体积发热量。温度升高时,极性介质(如环氧、聚酰亚胺)的 εr 通常略降,而 tanδ 在玻璃化转变温度附近可能先升后降,是高频绝缘件热失控的重要诱因。
磁导率 μ描述材料对磁场的响应,相对磁导率 μr = μ/μ0,μ0 = 4π×10-7 H/m。铁磁材料(硅钢、铁氧体、纳米晶)的 μr 可达数百至数万,但存在饱和与非线性:B = μ0μrH 仅在弱场线性区成立。磁滞损耗与涡流损耗都与 μ 相关,且铁磁材料的 μ 随温度上升而下降,居里点附近急剧跌落。做耦合分析时,若磁芯工作温度从 25 ℃ 升到 100 ℃,μr 可能下降 10%~30%,电感量随之变化,励磁电流与损耗分布都会偏移。
电导率 σ是电场驱动自由电荷迁移的能力,单位 S/m,其倒数为电阻率 ρe。导体涡流损耗密度为 pe = J2/σ = σE2。金属的 σ 随温度升高而下降,常用线性近似 σ(T) = σ20/[1 + α(T − 20)],铜的 α ≈ 0.00393 /℃。这意味着高温下导体电阻增大、涡流损耗反而略降,但直流铜损 I2R 上升,两者方向相反,需在耦合迭代中分别处理。
二、热侧三个参数:热量如何传导、储存与散出
热导率 k(W/(m·K))表征材料传导热量的能力,进入傅里叶定律 q = −k∇T。金属 k 高(铜约 400、铝约 237),导热绝缘材料低(环氧约 0.2~0.5、空气约 0.026),这一数量级差异决定了电机槽绝缘、功率模块衬板往往是热瓶颈。金属的 k 随温度略降,而部分陶瓷与半导体在高温下 k 下降更明显。
比热容 cp(J/(kg·K))表示单位质量材料升高 1 K 所需热量,与密度 ρ 共同决定体积热容 ρcp,控制瞬态温升速度。铜的 cp ≈ 385 J/(kg·K),环氧约 1000~1500 J/(kg·K)。稳态分析中 cp 不出现,但短时过载、脉冲工况下它决定温升斜率,忽略它会高估瞬时温升。
对流换热系数 h(W/(m²·K))不是材料属性,而是表面与流体换热条件的集总参数,进入牛顿冷却定律 q″ = h(Ts − Tf)。自然对流空气约 5~10,强制风冷约 20~100,水冷可达 1000~5000。h 强烈依赖流速、流态与表面状态,工程上常由经验关联式或实测反推,其不确定性往往是热模型误差的主要来源。
三、温度依赖性:耦合迭代中必须显式处理的部分
把六个参数放在一起看,会发现一个关键事实:电磁参数的温度漂移会改变损耗,热参数的温度漂移会改变散热能力,两者互为因果。工程上常见的处理方式有两种。其一是单向顺序耦合:先按常温参数算电磁损耗,再把损耗作为体热源导入热模型求温度,适用于温升较小(如 ΔT < 30 K)的工况。其二是双向迭代耦合:每轮用当前温度场更新 ε(T)、μ(T)、σ(T)、k(T),重算损耗与温度,直到相邻两轮最高温度差小于设定容差(例如 0.5 K)。对功率密度高、温升大的对象,双向迭代几乎是必需的。
以一台示例性高频变压器为例(以下数值为教学假设,非实测):磁芯初始 μr = 2000,磁芯损耗 8 W;绕组铜损 12 W;散热面 h = 25 W/(m²·K),环境 40 ℃。第一轮按常温参数求得热点约 95 ℃。若磁芯 μr 在 95 ℃ 时降至 1700,励磁电流上升约 15%,磁芯损耗可能增至 9.5 W,第二轮热点升至约 99 ℃。此时温度变化已小于 5 K,可认为收敛。这个例子说明:温度依赖不是修正项,而是可能改变结论的一阶效应。
四、参数取值与建模的实操建议
实际建模时,建议按以下顺序整理参数表:先确认每个参数是常数、分段线性表还是解析函数;再标注数据来源(datasheet、实测、经验值)与适用温度区间;最后检查单位与量纲一致性,尤其是 σ 用 S/m 还是 S/cm、h 用 W/(m²·K) 还是 mW/(mm²·K)。对于 tanδ、μr 这类强非线性参数,宁可给出 3~5 个温度点的表格并线性插值,也不要用单一常温值贯穿全温区。对于 h,应记录其对应的流速与冷却介质,因为同一个“风冷”在不同风速下可能相差数倍。
掌握这六个参数的定义与温度行为,后续的损耗计算、热网络搭建和耦合迭代才有可靠的输入基础。参数取值阶段的谨慎程度,往往比求解器设置更直接地决定最终分析结果的可信度。