2. 功率电子热源建模

液冷系统仿真的精度，首先取决于热源是否被正确描述。功率电子器件的损耗并不是一个固定数值，它随电流、母线电压、开关频率、结温、驱动电阻和调制方式持续变化。若把 IGBT 模块简单当作一个恒定发热块，液冷侧的流量、温升和结温预测都会出现系统性偏差。本节围绕三类主流器件——硅基 IGBT、SiC MOSFET 与功率二极管——说明损耗的物理来源、可计算的分量，以及如何把计算结果转成仿真可用的热源边界条件。

损耗的物理来源与分量划分

功率器件的总损耗通常拆成三部分：导通损耗、开关损耗和截止/驱动损耗。导通损耗来自器件在通态下的等效压降或通态电阻；开关损耗来自每次开通、关断过程中电压与电流交叠产生的能量；截止损耗在阻断状态下一般很小，但高温漏电流不可忽略。驱动损耗主要消耗在栅极回路，通常不进入主散热路径，建模时可单独处理或忽略。

IGBT 是双极型器件，通态压降由 PN 结压降和漂移区电阻压降共同决定，常用线性化模型 Vce = Vce0 + Rc·Ic 描述，其中 Vce0 是门槛压降，Rc 是通态等效电阻。SiC MOSFET 是单极型器件，通态特性更接近纯电阻，Vds = Rds(on)·Id，且 Rds(on) 随结温上升明显增大。二极管则用 Vf = Vf0 + Rd·If 描述正向压降，反向恢复过程会额外产生开关损耗。

导通损耗的计算

导通损耗的基本形式是电流与压降的乘积在一个基波周期内的平均值。对 IGBT，若已知某工况下的集电极电流有效值 Ic,rms 和平均值 Ic,avg，则

P_cond_IGBT = Vce0 · Ic,avg + Rc · Ic,rms²
P_cond_MOS  = Rds(on) · Id,rms²
P_cond_Diode= Vf0 · If,avg + Rd · If,rms²

在正弦电流、调制比为 m、功率因数为 cosφ 的典型逆变工况下，可用解析式快速估算。以 IGBT 为例，设峰值电流 Î，则

Ic,avg = Î/π · [1 + m·cosφ·π/4]
Ic,rms² = β/4 · [1 + 8·m·cosφ/(3π)]

这些公式的价值在于：不必先跑详细开关波形，就能得到与工况直接挂钩的导通损耗估计。实际建模时，更稳妥的做法是从器件手册的 Vce-Ic 或 Rds(on)-Tj 曲线中提取参数，并注意这些参数对应的结温条件。

开关损耗的计算

开关损耗由每次开关的能量 Eon、Eoff 和 Err（二极管反向恢复）乘以开关频率得到。手册通常给出在参考母线电压 Vref、参考电流 Iref、参考结温 Tref 和参考栅极电阻 Rg,ref 下的单脉冲能量。实际工况下需按比例缩放：

Eon(Ic, Vdc, Tj) ≈ Eon,ref · (Ic/Iref) · (Vdc/Vref) · (1 + kT·(Tj - Tref))
P_sw = fsw · (Eon + Eoff + Err) · (1/π)   // 正弦调制下的平均修正

其中 kT 是结温修正系数，SiC 器件的开关能量对温度相对不敏感，而硅 IGBT 的 Eoff 随温度上升较明显。栅极电阻 Rg 对开关损耗影响很大：增大 Rg 会降低 dv/dt 和 EMI，但开关能量近似线性增加。因此同一器件在不同驱动配置下，热源强度可能相差 30% 以上。

SiC MOSFET 的开关损耗通常显著低于同规格硅 IGBT，尤其在高温、高频下优势更明显；但其体二极管反向恢复虽小，仍会在桥式电路中产生额外损耗。二极管的反向恢复损耗 Err 与正向电流、di/dt 和结温相关，在硬开关桥臂中不可忽略。

从损耗到热源边界条件

算出总损耗后，下一步是把它分配到仿真模型中的热源位置。IGBT 模块内部通常包含多个芯片和二极管芯片，它们焊接在 DBC 基板上，再通过基板、导热硅脂和冷板散热。因此热源不能简单加在模块外壳，而应加在芯片层或等效热阻网络的功率注入节点。

常用做法是建立 Cauer 或 Foster 热阻网络：从芯片结到壳、壳到散热器、散热器到冷却液，逐级串联热阻和热容。仿真中把 P_total 注入结节点，通过热阻网络计算结温 Tj，再把 Tj 反馈回损耗模型，形成电-热耦合迭代。若只做单向热仿真，则需给定一个保守的结温假设，例如 Tj = 125°C 或 150°C，并说明该假设对结果的影响。

冷却液侧的边界条件同样关键：入口温度、流量、冷却液物性（乙二醇水溶液浓度）、流道几何和表面粗糙度共同决定对流换热系数。仿真中通常把冷板壁面设为第三类边界条件，即给定对流换热系数 h 和流体温度 Tf；若做共轭传热，则直接求解流固耦合，此时只需给定入口流量和入口温度。

典型示例：三相逆变桥单桥臂损耗估算

假设某 800V 母线、SiC MOSFET 三相逆变桥，单管有效值电流 Id,rms = 40A，峰值 Î = 56.6A，开关频率 fsw = 20kHz，调制比 m = 0.9，cosφ = 0.85，结温假设 Tj = 125°C。器件参数示例：Rds(on) = 8mΩ，Eon = 0.35mJ，Eoff = 0.25mJ，Err = 0.08mJ（均为参考条件下的合理假设值）。

P_cond = 0.008 × 40² = 12.8 W
P_sw   = 20000 × (0.35 + 0.25 + 0.08)e-3 / π ≈ 4.3 W
P_total ≈ 17.1 W / 管

一个桥臂上下两管合计约 34W，三相桥共约 103W。这个量级与液冷系统热设计直接相关：若冷板热阻为 0.15K/W，则仅器件损耗就会带来约 15K 的温升。实际建模时还需叠加二极管、驱动和母排损耗，并按占空比和电流波形修正。

建模中的常见偏差与处理

第一类偏差来自参数条件不匹配。手册中的 Eon、Eoff 通常标注在特定 Rg、Vdc 和 Tj 下，若实际驱动电阻或母线电压不同，必须缩放。第二类偏差来自电流波形假设：解析式假设正弦电流，而实际相电流可能含谐波或处于过调制区，此时应改用数值积分或从仿真波形中提取 Ic,avg 和 Ic,rms²。第三类偏差来自热耦合：多芯片之间通过 DBC 和基板存在热耦合，单管独立热阻网络会低估相邻芯片的温升，高功率密度模块建议使用有限元提取的耦合热阻矩阵。

边界条件设定还需注意单位与参考系一致性：热阻用 K/W，热容用 J/K，功率用 W，流量用 L/min 或 m³/s，对流换热系数用 W/(m²·K)。若冷板流道内出现局部沸腾或流量分配不均，单值 h 假设会失效，此时应改用共轭传热或分区边界条件。完成建模后，建议用一组已知工况做校核：例如给定总损耗和入口温度，检查仿真结温是否落在热阻网络手算值的 ±10% 以内，以此判断热源映射和边界条件是否合理。