2. 功率二极管直流特性建模

功率二极管是功率变换器中最基础的器件之一，其直流特性直接决定了导通损耗、反向阻断能力和热设计余量。在SPICE模型中，直流特性由PN结方程及其派生参数描述，理解这些参数的物理含义和提取方法，是建立可用于工程仿真的二极管模型的第一步。

PN结方程与SPICE中的直流模型

理想PN结的电流-电压关系由Shockley方程给出：

I_D = I_S · (exp(V_D / (N · V_T)) - 1)

其中 I_S 为饱和电流，N 为发射系数（理想因子），V_T = kT/q 为热电压，室温（300 K）下约为25.85 mV。SPICE中的二极管模型在此基础上增加了串联电阻 RS、反向击穿电压 BV、击穿电流 IBV 等参数，以覆盖功率二极管在大电流和反向偏置下的实际行为。

需要特别注意的是，SPICE二极管模型中的 V_D 是结电压而非端电压。当电流较大时，RS 上的压降不可忽略，端电压 V_AK = V_D + I_D · RS。这一区分在功率二极管建模中至关重要，因为功率器件的 RS 通常在毫欧到数十毫欧量级，几安培电流下就会产生几十到几百毫伏的额外压降。

正向压降与串联电阻的分离

功率二极管的正向压降 V_F 由两部分组成：结压降 V_j 和串联电阻压降 I_F · RS。在小电流下，结压降占主导，V_F 随电流对数增长；在大电流下，RS 压降线性增长并逐渐成为主要部分。建模时若只拟合小电流区，得到的 N 和 IS 在大电流下会严重低估 V_F。

实际操作中，可以从数据手册的正向特性曲线取两组数据点来初步估算。例如某肖特基二极管在25°C下，I_F = 1 A 时 V_F ≈ 0.45 V，I_F = 10 A 时 V_F ≈ 0.62 V。假设 RS = 15 mΩ，则扣除 RS 压降后：

V_j1 = 0.45 - 1 × 0.015 = 0.435 V
V_j2 = 0.62 - 10 × 0.015 = 0.470 V
ΔV_j = 0.470 - 0.435 = 0.035 V
N = ΔV_j / (V_T · ln(10/1)) = 0.035 / (0.02585 × 2.303) ≈ 0.588

该 N 值小于1，说明肖特基二极管以多数载流子导电为主，与PN结的扩散电流机制不同。对于普通PN结功率二极管，N 通常在1.0～2.0之间；若 N 接近2，则表明复合电流占比较大。

反向漏电与温度效应

反向漏电 I_R 由 IS 直接控制。SPICE中 IS 随温度的变化关系为：

IS(T) = IS(T_nom) · (T/T_nom)^(XTI/N) · exp(EG·q/(N·k) · (1/T_nom - 1/T))

其中 EG 为带隙能量（硅约1.11 eV，肖特基约0.69 eV），XTI 为饱和电流温度指数（PN结约3.0，肖特基约2.0）。温度每升高10°C，硅二极管反向漏电大约翻一倍；而正向压降则具有约 -2 mV/°C 的温度系数。这两个相反的趋势意味着：高温下正向损耗降低，但反向漏电和热失控风险增加。

理想因子提取与直流曲线拟合流程

理想因子 N 的提取应基于半对数坐标下的正向I-V曲线。在半对数坐标中，指数区的 ln(I_D) 与 V_D 呈线性关系，斜率为 1/(N·V_T)。具体步骤如下：

1. 测量或获取正向I-V数据，覆盖从微安级到额定电流的至少三个数量级。
2. 在半对数坐标中绘制 ln(I_D) vs V_D，识别低电流区的线性段。
3. 对线性段做最小二乘拟合，得到斜率 k，则 N = 1/(k · V_T)。
4. 由线性段截距得到 IS。
5. 在大电流区，用 V_F - V_j 与 I_F 的线性关系提取 RS。
6. 将 IS、N、RS 填入SPICE模型，运行直流扫描，与原始数据对比，调整参数直至误差可接受。

以下是一个可用于对比的SPICE二极管模型卡示例（以某600 V/10 A Si二极管为例，参数为教学假设值）：

.model D_PWR_600V D (
+ IS=1e-9      ; 饱和电流
+ N=1.35       ; 理想因子
+ RS=0.025     ; 串联电阻 25 mΩ
+ BV=600       ; 击穿电压
+ IBV=1e-4     ; 击穿电流
+ EG=1.11      ; 硅带隙
+ XTI=3.0      ; 温度指数
+ TT=5e-6      ; 渡越时间
+)

拟合时建议先固定 RS 和 BV，仅调整 IS 和 N 使正向指数区匹配；然后调整 RS 使大电流区匹配；最后用高温数据校验 EG 和 XTI。若反向漏电在室温下匹配但高温下偏差大，优先检查 XTI 和 EG 的取值是否与器件类型一致。