1. 三相整流器双闭环控制仿真实战
做电力电子仿真这么多年，三相PWM整流器始终是个经典课题。这次接到个车载充电机的预研项目，要求用Simulink搭建一套15kW整流器模型，正好把双闭环PI控制的全套设计过程记录下来。这个仿真最迷人的地方在于，它完美展现了如何用数学工具驯服电力电子系统的"野性"——既要让750V直流电压稳如泰山，又要保证网侧电流正弦度满足THD<5%的严苛标准。
先看系统关键参数：
- 交流输入：三相220V/50Hz
- 直流输出：750V/15kW
- 滤波电感：1.8mH（每相）
- 开关频率：20kHz
- 控制方式：电压外环+电流内环双PI控制
- 调制策略：SPWM异步调制
提示：车载系统对体积敏感，1.8mH电感是权衡开关损耗与谐波抑制后的折中选择，后文会详细分析其设计依据。
2. 主电路建模与参数设计
2.1 拓扑结构选择
采用典型的三相两电平电压型PWM整流器结构，主电路包含：
- 三相交流电压源（220Vrms/50Hz）
- L型滤波器（简化设计，实际产品多用LCL）
- IGBT全桥模块（6个开关管）
- 直流侧电容（取2000μF）
matlab复制% Simulink中主电路关键参数设置
AC_Voltage = 220*sqrt(2); % 峰值电压311V
L_filter = 1.8e-3; % 每相滤波电感
C_dc = 2000e-6; % 直流支撑电容
R_load = 750^2/15e3; % 等效负载电阻37.5Ω
2.2 关键器件参数计算
电感选型公式：
[ L_{min} = \frac{V_{dc}}{4 \cdot f_{sw} \cdot \Delta I_{pp}} ]
其中ΔIpp允许纹波电流取额定电流20%：
[ I_{rated} = \frac{P}{\sqrt{3} \cdot V_{ac}} = \frac{15kW}{\sqrt{3} \times 220V} ≈ 39.4A ]
[ \Delta I_{pp} = 0.2 \times 39.4 ≈ 7.88A ]
[ L_{min} = \frac{750V}{4 \times 20kHz \times 7.88A} ≈ 1.19mH ]
实际取1.8mH留有设计余量，仿真验证发现：
- 1.0mH时THD=8.2%（超标）
- 1.8mH时THD=4.7%（达标）
- 2.5mH时THD=4.1%（改善有限但体积增大）
电容选择经验公式：
[ C \geq \frac{2P_o}{\omega V_{dc} \cdot \Delta V_{pp}} ]
取电压纹波ΔVpp<1%即7.5V：
[ C \geq \frac{2 \times 15kW}{314 \times 750V \times 7.5V} ≈ 1695μF ]
选用2000μF电容后，实测纹波仅4.3V。
3. 控制策略实现细节
3.1 双闭环控制架构
采用经典的dq旋转坐标系控制：
- 外环：直流电压PI控制，输出d轴电流参考
- 内环：交流电流PI控制，生成PWM占空比
mermaid复制graph TD
A[电压外环] -->|Vdc_ref - Vdc| B(PI_V)
B -->|Id_ref| C[电流内环]
D[电网电压] --> E(PLL)
E -->|θ| F[abc/dq变换]
C -->|Vd_Vq| G[PWM生成]
3.2 PI参数整定实战
电流内环设计（响应速度要求高）：
- 计算临界增益法初值：
[ K_{p,i} = L \cdot 2\pi f_c \cdot 0.5 ]
取穿越频率fc=1kHz（开关频率1/20）：
[ K_{p,i} = 1.8mH \times 6283 \times 0.5 ≈ 5.65 ] - 考虑线路电阻R≈0.1Ω：
[ K_{i,i} = K_{p,i} \cdot R/L ≈ 314 ] - 实际调试最优值：
- Kp=0.35（较理论值减小16倍）
- Ki=2000（较理论值增大6倍）
注意：理论计算值通常偏大，实际需考虑数字控制延迟、PWM非线性等因素。
电压外环设计（带宽取内环1/10）：
- 经验公式：
[ K_{p,v} = \frac{C}{2T_s} ]
取Ts=100μs（控制周期）：
[ K_{p,v} = \frac{2000μF}{2 \times 100μs} = 10 ] - 实测稳定参数：
- Kp=0.01（减小1000倍）
- Ki=5
参数调整秘诀：
- 先调电流环：固定Ki=0，增大Kp至临界振荡
- 再调Ki值：取Kp的1/10开始递增
- 电压环同理，但调整幅度要小一个数量级
4. SPWM调制优化技巧
4.1 异步调制实现
采用S函数避免过零点畸变，核心算法：
matlab复制function sys=mdlOutputs(t,x,u)
carrier = sawtooth(2*pi*20e3*t, 0.5); % 20kHz三角载波
mod_wave = u(1)*sin(2*pi*50*t + phase_shift);
pwm = (mod_wave > carrier) - (mod_wave < -carrier);
end
4.2 死区补偿策略
添加0.5μs死区时间后，需在调制波中注入补偿电压：
[ V_{comp} = \frac{T_{dead}}{T_{sw}} \cdot V_{dc} ]
对于750V系统：
[ V_{comp} = \frac{0.5μs}{50μs} \times 750V = 7.5V ]
实测效果：
- 无补偿时THD恶化至6.3%
- 补偿后THD恢复至4.9%
5. 仿真结果分析
5.1 稳态性能
| 指标 | 实测值 | 标准要求 |
|---|---|---|
| 直流电压 | 749.8V | 750V±1% |
| 网侧THD | 4.7% | <5% |
| 功率因数 | 0.998 | >0.99 |
| 效率（估算） | 97.2% | >96% |
5.2 动态响应测试
- 突加负载（5kW→15kW）：
- 电压跌落：750V→700V（6.7%）
- 恢复时间：0.18秒
- 电网电压骤降（220V→180V）：
- 直流电压波动：±15V
- 恢复时间：0.25秒
5.3 关键波形展示
- 电网电压与电流（完美同相位）
- 直流侧电压纹波（4.3Vpp）
- 相电流频谱分析（50Hz基波占比95.3%）
6. 工程经验总结
-
电感饱和问题：实际产品要选用降额30%的磁芯，仿真中1.8mH电感在40A电流时磁通密度应留有余量
-
散热设计要点：20kHz开关频率下，IGBT损耗估算：
[ P_{sw} = \frac{1}{2} V_{dc} I_o (t_r + t_f) f_{sw} ]
以FF300R12KE3为例：
[ P_{sw} ≈ 0.5 \times 750V \times 40A \times 100ns \times 20kHz = 30W ]
需配150W散热器 -
参数敏感度测试：
- 电感±20% → THD变化±1.5%
- 开关频率±5kHz → 效率变化∓0.8%
- PI参数±30% → 恢复时间变化±50%
这个仿真模型最让我惊喜的是，用最传统的PI控制也能达到高性能指标。那些盲目追求先进控制算法的同行可能需要反思——有时候把经典方法吃透，比追逐新技术更有效。下次可以尝试在电流环加入前馈补偿，应该能把动态响应再提升一个档次。