4. Vienna整流器数学模型

Vienna整流器是一种三电平、三开关桥臂的Boost型PWM整流器,每相桥臂由两个串联二极管和一个双向开关构成,输出侧分裂为上下两个电容,形成正母线P、中性点O、负母线N三个电平。要设计电流环、电压环和锁相环,必须先建立能反映其开关周期平均行为的数学模型。本节从三相静止坐标系abc出发,依次给出坐标变换理论、αβ坐标系模型和dq旋转坐标系模型,并说明各模型的适用场合与推导要点。

一、建模前的约定与假设

为避免符号混乱,先统一参照约定:三相电网电压记为 ea、eb、ec,三相入端电流记为 ia、ib、ic,均以从电网流入整流器为正方向;交流侧电感为 L,等效串联电阻为 R;上下电容各为 C,负载电流为 iL;开关函数 Sx(x=a,b,c)取值为 1、0、-1,分别对应该相接入正母线、中性点、负母线。建模采用开关周期平均法,假设开关频率远高于电网频率,忽略开关纹波;三相电网为理想对称正弦,且中性点电位平衡。

二、三相静止坐标系abc模型

根据基尔霍夫电压定律,每相回路可写为:

L·di_a/dt = e_a - R·i_a - v_aO
L·di_b/dt = e_b - R·i_b - v_bO
L·di_c/dt = e_c - R·i_c - v_cO

其中 vaO、vbO、vcO 为各相桥臂中点对中性点O的电压。由于 Vienna 整流器每相输出三电平,桥臂中点电压由开关函数和电容电压决定。设上电容电压为 uC1、下电容电压为 uC2,在电容电压均衡时 uC1=uC2=Udc/2,则:

v_aO = S_a·U_dc/2
v_bO = S_b·U_dc/2
v_cO = S_c·U_dc/2

直流侧电流由三相开关函数与相电流共同决定,对中性点列写节点电流方程可得:

C·du_C1/dt = (S_a·i_a + S_b·i_b + S_c·i_c)/2 - i_L
C·du_C2/dt = -(S_a·i_a + S_b·i_b + S_c·i_c)/2 - i_L

abc 模型物理意义直观,但三相变量均为时变正弦量,且相互耦合,难以直接用于控制器设计。因此需要引入坐标变换,把交流量变为直流量或两相正交量。

三、坐标变换理论

坐标变换的核心是保持功率或幅值不变的正交变换。常用 Clarke 变换将 abc 三相量映射到两相静止 αβ 坐标系,Park 变换再将 αβ 量旋转到与电网同步的 dq 坐标系。采用等幅值变换时,Clarke 变换矩阵为:

[x_α]   [ 1   -1/2   -1/2 ] [x_a]
[x_β] = [ 0   √3/2  -√3/2 ] [x_b]
                              [x_c]

Park 变换矩阵为:

[x_d]   [ cosθ   sinθ ] [x_α]
[x_q] = [-sinθ   cosθ ] [x_β]

其中 θ 为电网电压矢量与 α 轴的夹角,由锁相环给出。变换前后需保持同一参照方向,否则 dq 轴电流的物理含义会反转。若采用等功率变换,矩阵前需乘 √(2/3) 系数,本节统一采用等幅值变换,便于与电流环给定值直接对应。

Vienna整流器坐标系变换与模型关系 三相静止 abc e_a, e_b, e_c i_a, i_b, i_c S_a, S_b, S_c ∈ {1,0,-1} 时变正弦、相互耦合 Clarke 等幅值 两相静止 αβ e_α, e_β i_α, i_β v_α, v_β 正交、仍为交流量 Park θ 由PLL给出 两相旋转 dq e_d, e_q i_d, i_q v_d, v_q 稳态为直流量 dq 旋转坐标系电压方程(ω 为电网角频率) L·di_d/dt = e_d − R·i_d + ωL·i_q − v_d L·di_q/dt = e_q − R·i_q − ωL·i_d − v_q 耦合项 +ωL·i_q 与 −ωL·i_d 需前馈解耦 直流侧与中性点方程 C·du_C1/dt = (S_a·i_a + S_b·i_b + S_c·i_c)/2 − i_L C·du_C2/dt = −(S_a·i_a + S_b·i_b + S_c·i_c)/2 − i_L 反变换用于生成调制波

四、αβ坐标系模型推导

将 abc 三相电压方程两边同时左乘 Clarke 变换矩阵,并利用三相无中线系统中 ia+ib+ic=0 的条件,可得 αβ 坐标系下的电压方程:

L·di_α/dt = e_α - R·i_α - v_α
L·di_β/dt = e_β - R·i_β - v_β

其中 vα、vβ 为桥臂中点电压在 αβ 轴上的分量。αβ 模型把三相耦合降为两相独立方程,结构简洁,适合用于滞环电流控制、模型预测控制以及空间矢量调制。但 αβ 轴上的电流给定仍是正弦量,若要在同步旋转坐标系下做 PI 无静差调节,还需进一步变换到 dq 坐标系。

五、dq旋转坐标系模型推导

将 αβ 方程代入 Park 变换,并注意旋转坐标系下微分算子会引入旋转电动势项。设 d 轴与电网电压矢量 e 对齐,则 ed=Em、eq=0。经过推导可得 dq 坐标系电压方程:

L·di_d/dt = e_d - R·i_d + ωL·i_q - v_d
L·di_q/dt = e_q - R·i_q - ωL·i_d - v_q

式中 ωL·iq 和 −ωL·id 为 d、q 轴之间的耦合项,来源于坐标旋转。稳态时 did/dt=diq/dt=0,方程退化为代数关系,d、q 轴电流为直流量,PI 调节器可实现无静差跟踪。这也是 Vienna 整流器电流内环普遍采用 dq 模型的原因。

六、三种模型的对比与选用

模型变量形式耦合程度典型用途
abc 模型三相时变正弦量三相强耦合物理机理分析、仿真验证
αβ 模型两相正交交流量两相解耦滞环控制、模型预测、SVPWM
dq 模型两相直流量存在旋转耦合项PI 电流环、前馈解耦、锁相环
推导要点:Clarke 变换采用等幅值形式时,dq 轴电流幅值与相电流幅值相等,便于直接设定电流给定;若采用等功率形式,需在控制器中同步调整给定值系数,否则会出现 1.5 倍或 √(2/3) 倍的稳态偏差。
注意:Vienna 整流器开关函数取值为 1、0、-1,与两电平变换器的 0、1 不同。在推导 vα、vβ 时,必须把三电平特性代入,否则 dq 模型无法反映中性点电位波动对电流的影响。

以一组教学示例参数说明模型的使用:设电网相电压幅值 Em=311 V,频率 50 Hz,L=2 mH,R=0.1 Ω,Udc=700 V,开关频率 20 kHz。在 dq 坐标系下,若要求单位功率因数运行,则 iq*=0,id* 由电压外环给出。将上述参数代入 dq 方程,可估算出稳态时 vd≈ed−R·id+ωL·iq,vq≈−R·iq−ωL·id,其中 ωL=2π×50×0.002≈0.628 Ω,耦合项量级较小,但在高功率或低开关频率场合不可忽略。该示例仅用于说明参数代入方式,实际数值需结合具体样机重新核算。

综上,abc 模型是理解 Vienna 整流器物理过程的起点,αβ 模型是简化分析与调制实现的桥梁,dq 模型则是控制器设计的核心。三者通过 Clarke 和 Park 变换相互联系,推导过程中需始终保持参照方向、变换系数和开关函数定义的一致性,才能保证后续电流环与电压环设计的正确性。