3. 热仿真与实验对标的核心指标体系

热仿真与实验对标之所以困难,根本原因在于双方比较的往往不是同一个物理量。仿真输出的是网格单元上的温度场、节点上的结温、由温差与热流反推的热阻;实验测到的是热电偶电压、红外辐射亮度、瞬态电压曲线。要让两者可比,必须先把“温度场、结温、热阻、热流密度、对流换热系数、辐射角系数、瞬态热响应时间常数”这七类关键量各自的物理定义、测量方式与仿真提取方式对齐,再谈误差。本节按“定义—测量—仿真提取”三条线逐一梳理,并给出对标时的口径约定。

温度场:空间分布量,对标要固定切片与坐标

温度场是温度在空间与时间上的分布 T(x, y, z, t),单位 ℃ 或 K。仿真中它天然存在,是求解结果本身;实验中它只能被离散采样。对标温度场时最容易犯的错误是拿仿真最高温去比某几个测点温度,二者不是同一对象。正确做法是约定切片:例如在 PCB 表面上方 0.5 mm 处取一个平行平面,或沿某条过芯片中心的截面线,把仿真结果插值到与实验测点相同的坐标上再比较。

实验获取温度场的手段包括热电偶阵列、红外热像仪和热敏涂料。热电偶给出点温,精度高但空间分辨率受布点限制;红外热像仪给出面分布,但读数是辐射温度,需用发射率修正:T_true 由 ε·σ·T_ir⁴ = ε·σ·T_true⁴ 关系反推,实际按 T_true = T_ir / ε^(1/4)(近似,单位 K)修正,金属表面 ε 低至 0.05~0.2,不修正会带来几十度误差。仿真提取温度场时,应输出与实验同坐标、同高度的切片数据,并记录网格尺寸,因为网格越粗,峰值温度越容易被“抹平”。

结温:芯片内部最高温度,测量与仿真口径差异最大

结温 T_j 是半导体芯片有源区的温度,是可靠性判据的核心。仿真中 T_j 取芯片体(die)内的最高节点温度,需注意是否包含金属层、是否取体积平均还是面平均。实验中 T_j 无法直接接触测量,常用方法有三类:一是电学法(如 VBE 或体二极管正向压降作为温度敏感参数 TSP),先标定 TSP-温度曲线,再在通电状态下测 TSP 反推 T_j;二是红外显微测温,需去除封装或使用红外透明封装;三是热敏电阻贴装,测到的是壳温或板温而非真正结温。

对标时务必统一“结温”的定义位置:是 die 最高点、die 平均温度,还是结-壳路径上的等效节点。工程上常约定 T_j 为 die 有源区面积平均温度,仿真按同一区域做面积加权平均,实验用电学法标定,二者才具备可比性。若仿真取最高节点而实验取平均,误差可能被高估 5~15 ℃。

热阻:从温差与热流定义的派生量

热阻 R_th = ΔT / P,单位 K/W,其中 ΔT 是两点温差,P 是流经该路径的热功率。它不是一个独立测量量,而是由温度和功率算出来的。结-壳热阻 R_thJC 用 T_j 与壳温 T_c 之差除以芯片功耗;结-板热阻 R_thJB 用 T_j 与板温之差除以功耗。仿真提取时,需要在芯片下方设置一个等温参考面作为“壳”,读取该面平均温度,再除以施加的体热源功率。实验则按 JEDEC 标准(如 JESD51 系列)在规定冷板上测量。

对标热阻的关键是功率口径一致:仿真中施加的体热源功率是否等于实验中电功率减去光功率或其它输出功率。对 LED 而言,约 30%~50% 的电功率转化为光,若仿真把全部电功率当热源,热阻会被系统性低估。建议在仿真中显式设置热功率 P_heat = P_elec − P_optical,并在报告中注明。

热流密度与对流换热系数:边界条件的对标

热流密度 q″ = P / A,单位 W/m²,描述单位面积通过的热功率。仿真中它由温度梯度与导热系数决定:q″ = −k·∂T/∂n,可在后处理中沿指定面提取。实验中热流密度用热流计或由已知热阻的辅助结构反推。对标热流密度能直接暴露边界条件设置是否合理。

对流换热系数 h = q″ / (T_s − T_∞),单位 W/(m²·K),是仿真中最不确定的输入之一。自然对流下 h 通常为 5~15 W/(m²·K),强制风冷可达 20~100 W/(m²·K)。仿真中 h 作为边界条件直接给定,或用 CFD 求解得到;实验中 h 需通过加热平板测温差与功率反算。对标时若仿真用固定 h、实验是真实流场,二者差异往往来自风速分布与遮挡,应优先核对风速测点位置与仿真入口条件。

辐射角系数与瞬态时间常数

辐射角系数 F_ij 表示表面 i 发出的辐射能中直接落到表面 j 的比例,取值 0~1,满足 ΣF_ij = 1。仿真中由几何与遮挡关系自动计算,可在辐射模块中输出角系数矩阵;实验中难以直接测量,通常通过辐射热流与表面温度反推等效辐射换热。对标时重点检查发射率设置与遮挡面是否被忽略。

瞬态热响应时间常数 τ 描述温度随时间的指数变化,T(t) = T_∞ + (T_0 − T_∞)·e^(−t/τ)。仿真中可由阶跃功率下的温升曲线拟合得到,也可由热容与热阻乘积估算 τ ≈ R_th·C_th。实验中用阶跃功率激励,记录 TSP 或热电偶响应,拟合出 τ。对标 τ 能有效检验热容、材料密度与比热设置是否正确,因为稳态温度对热容不敏感,而瞬态对热容非常敏感。

热仿真与实验对标:七类核心量的定义—测量—仿真提取 仿真侧提取 温度场 T(x,y,z,t) 切片插值 / 网格节点 结温 T_j die 体最高/面积平均 热阻 R_th = ΔT/P 等温参考面 + 体热源 热流密度 q″ = −k∂T/∂n 后处理面通量 对流换热系数 h 边界条件 / CFD 求解 辐射角系数 F_ij 几何遮挡自动计算 瞬态时间常数 τ 实验侧测量 热电偶 / 红外热像 点温 / 面温(需 ε 修正) 电学法 TSP 标定 V_BE / 体二极管 JESD51 冷板法 T_j、T_c、P 同步测 热流计 / 辅助热阻 反推 q″ 加热平板 + 风速仪 h = q″/(T_s−T_∞) 辐射热流 + 表面温度 反推等效辐射 阶跃功率 + 响应拟合 对标口径约定 同坐标、同高度切片 避免峰值 vs 点温错配 T_j 定义位置统一 最高点 / 面积平均 热功率口径一致 P_heat = P_elec − P_optical 边界条件可追溯 h、ε、风速入口记录 发射率修正 金属 ε 低,必须标定 瞬态激励形式一致 阶跃幅值、起始温度 误差带与判据预先约定

对标实施:一份可填写的口径核对表

把上述七类量落到一次具体对标中,建议先填一张口径核对表,再开始跑仿真或做实验。下表以某 15 W 功率器件、自然对流散热为例,参数为教学示例,非实测结果。

关键量 仿真提取方式 实验测量方式 对标口径 示例值
温度场 PCB 上方 0.5 mm 切片插值 红外热像,ε=0.92 同坐标、同高度 峰值 78.5 ℃
结温 T_j die 面积加权平均 V_BE 电学法标定 面积平均,非最高点 92.3 ℃
热阻 R_thJB (T_j−T_board)/P_heat JESD51 冷板法 P_heat = 15 W×0.7 4.2 K/W
热流密度 q″ 底面 −k∂T/∂n 面积分 热流计贴装 同一投影面积 1.8×10⁴ W/m²
对流换热系数 h 边界条件 8 W/(m²·K) 加热平板反算 同一环境温度 7.5 W/(m²·K)
辐射角系数 F_ij 辐射模块输出矩阵 等效辐射热流反推 ε 与遮挡一致 F=0.31
时间常数 τ 阶跃功率拟合 阶跃功率 + TSP 响应 同幅值、同起始温度 42 s
对标前先问三个问题:比的是同一个物理量吗?定义位置一致吗?功率与边界条件口径一致吗?三问过关,误差才有讨论意义。
常见陷阱:用仿真最高节点温度对比实验平均结温,会把误差放大数倍;用全部电功率当热源计算热阻,会系统性低估 R_th;红外测温不做发射率修正,金属表面读数可能偏低 20 ℃ 以上。

完成一次对标后,建议按量分别记录偏差:温度场看切片最大偏差与均方根偏差,结温看绝对偏差,热阻看相对偏差,瞬态看 τ 的相对偏差。若结温偏差小于 5 ℃、热阻偏差小于 10%、τ 偏差小于 15%,通常认为模型可用于该工况下的工程判断;若某一项显著超差,优先回到该量的口径定义与边界条件,而不是盲目调整材料参数。