Vienna 整流器系统双PI闭环仿真

三相三电平六开关 VIENNA 整流器的电路拓扑结构如下图所示：

在该电路拓扑结构中，A、B、C为三相对称输入的交流电压源，电压峰值设置为220*1.414V，频率设为50Hz，通过采样得到三相电路各自的瞬时电压 ua、ub、uc和各自的瞬时电流ia、ib、ic。RL1、RL2、RL3为三相电路上的电感和等效的电阻，电阻大小设为0.05Ω，电感设为425e-6H。D1~D6是6个箝位二极管，导通电阻设为0.001Ω，正向电压设为0.8V，缓冲电阻大小设为1e5Ω，缓冲电容大小设为inf；Sa、Sb、Sc是各相桥臂的双向开关管，均由2个参数相同的IGBT组成，IGBT内阻为1e-3，缓冲电阻大小为1e5Ω，缓冲电容大小设为inf，通过在每个周期内对各开关段通断的控制，可以实现交流侧单位功率因数运行，并且稳定直流侧输出电压。C1、C2是2个等值滤波电容，用于减小输出电压纹波，从而稳定直流侧输出电压，电容容值为1100e-6F，初始电压设为200V，通过采样得到上下两个均压电容的电压分别为Vc1、Vc2。R代表直流侧的负载，其大小为10Ω，并对输出电压和输出电流进行采样，得到io和vo，并计算出输出的功率Po。

在该仿真工程中，采样电压电流双闭环控制方法。电压外环采用传统的 PI 控制，电流内环先采用 PI 前馈解耦，再用 SVPWM 进行调制。三相三电平 VIENNA 整流器控制回路结构框图如下所示。

电压外环仿真结构框图如下所示。

在该仿真工程中，电压内环采用传统的PI控制器，参考电压为600V。用参考电压V_ref减去直流侧采样得到的电压Vdc得到差值，经过离散形式的PI调节器，得到电流内环给定值Id*。锁相环仿真结构框图如下所示。

该模块的输入量为A、B、C三相的电压，选定电压基准值为100V，通过乘以系数0.01，将三相输入电压转换为标幺值送到simulink中的PLL(3ph)模块中，三相锁相环的参数为：最小频率设为45Hz，表示输入信号的频率不会小于该频率；初始输入相位和频率设为0°和50Hz；调节器增益[Kp,Ki,Kd]设为[180,3200,1]；微分作用时间常数为1e-4；频率最大变化率为12Hz/s；滤波器频率测量的截止频率为25Hz；采样时间为1e-4。三相锁相环输出的值有频率和ωt，此处只需使用输出的相角值。从三相锁相环的功能可知，输出的ωt是一个和A相正序同步的三角波，其周期为0.02s，纵轴取值范围为0到2pi，将该信号分别取正弦和余弦后得到周期为0.02s，幅值为1的正弦波和余弦波，将两个波形一起输出为sin-cos。电流内环仿真结构框图如下所示。

前馈解耦策略控制结构图如下所示。

利用输入的三相电压、三相电流已经相角，经过Park变换，将三相静止坐标系下的值Vabc转换成两相同步旋转坐标系下的值Vdq0。设置离散PI控制器的比例系数为0.5，积分系数为1000。使用基于PI的电流前馈解耦，使得id和iq不再相互耦合，实现了电流在dq坐标系下的解耦。因为其是对有功分量和无功分量的单独控制，只要令无功电流参考值为 0，则系统可以实现单位功率因数。经过前馈解耦之后，得到dq坐标系下的参考电压，经过Park逆变换得到在三相静止坐标系下的电压参考值Vabc*。SVPWM模块仿真结构框图和大扇区划分模块仿真结构框图分别如下所示。

SVPWM调制仿真模块结构图中，首先将采样得到的三相电流经过大扇区划分模块得到此时所处的区间。具体原理如下：

判断ia、ib、ic电流的正负情况，若A相为正，则计数值N加4；若B相为正，则计数值N加2；若C相为正，则计数值N加1。得到的计数值N的表达式为N=4A+2B+C，计数值与大扇区的对应关系为：

将电流内环得到的三相静止坐标系下的电压参考值Vabc*通过Clarke变换转换到两相静止坐标系下，得到Vα和Vβ。两电平大扇区到小扇区坐标变换模块仿真结构框图和大扇区向小扇区的转换关系表如下所示。

得到两电平下小扇区的目标矢量后，根据下图进行小扇区的判断。

将大扇区分成6个正三角形小扇区，判断条件如下表所示：

具体实现过程是首先将得到的vα*和vβ*进行相应的运算，得到：

分别与阈值0进行比较，其中Uref1比较得到的值记作A，Uref2比较得到的值记作B，Uref3比较得到的值记作C，得到的小扇区编号Vqy=A+2B+4C。

在判断完电压矢量处于何小扇区后，就需要确定起始矢量、终止矢量和零矢量各自的作用时间Ti、Tj、T0，实现该过程的仿真框图如下图所示。

经过计算可知各扇区内基本矢量作用时间对应表如下表所示：

得出Ti和Tj表达式后，再判断Ti与Tj之和是否大于Ts=0.00001，若Ti+Tj<Ts，则t1=Ti，t2=Tj；否则t1=Ts*Ti/(Ti+Tj)，t2=Ts*Tj/(Ti+Tj)。

通过表达式可得：

各小扇区占空比的计算公式如下表所示：

在进行最后的开关管触发信号的生成前，需要对均压电容中点电压的平衡进行控制。零矢量、大矢量对中点电位不影响，中矢量对中点电位的影响无法控制，只有小矢量对中点电位的影响可控，而且存在正小矢量和负小矢量，对中点电位的影响正好相反，因此可以通过控制正负小矢量的工作时间来调整中点电位电压。当中点电位偏低时，提高正小矢量的作用时间；当中点电位偏高时，提高负小矢量的作用时间。

上图为中点电压平衡在simulink中实现的仿真结构框图，将上下电容的差值netural与参考值Vref做差，并通过连续时间形式的PID控制器，经过限幅后乘以开关周期Ts=0.00001，得到中点电压平衡的调整时间time。

上图为各大扇区中电压正负值判断模块仿真结构框图，根据大扇区号N的不同来区分所处的大扇区中三相电压的正负值：2表示对应的电压为正，1表示对应的电压为负。

上图为3个开关管占空比转化触发信号模块仿真结构框图，得到各开关管的触发脉冲。主要工作原理是，首先将得到的时间Ta、Tb、Tc以及中点平衡时间netural同时叠加到生成的周期与开关周期Ts=0.00001相同且幅值为5e-6的对称锯齿波上，因为时间常数值不同，锯齿波会产生不同程度的上下平移。锯齿波的波形位于时轴上方时，比较器输出0，当位于时轴下方时，输出1，首先PWM的调制。最后根据上述表7判断得到的三相电压的正负值来输出不同的触发信号。当电压为负，即X(X=A,B,C)=1时，直接输出比较器得到的值；否则当电压为正，即X=2时，取反后输出比较器得到的值。