很多人一开始接触四桥臂逆变器的时候第一反应往往是“多加一个桥臂不就多了一路开关管嘛直接按三桥臂的思路多给一路PWM不就行了”。等真正动手把三相四线制负载接上去才发现事情没那么简单第四桥臂不是用来“多带一路负载”的它带来的是零序电流的独立通路也把调制策略从二维平面硬生生推到了三维空间。这篇博客要梳理的就是以35d为代号的一次三维空间矢量调制算法仿真探索35d是这次算法方案的项目代号d对应3D-SVM也就是四桥臂逆变器三维空间矢量调制算法的第三轮验证版本。文章会把原理、实现、仿真、调试踩坑完整串起来适合正在做逆变器控制、UPS、有源滤波器、微电网并网装置或者准备用仿真平台验证调制算法的同学参考。阅读之前只需要一个基础认知三相三桥臂逆变器的SVPWM是在αβ二维平面里选扇区而四桥臂逆变器要在αβγ三维空间里选四面体。1. 四桥臂逆变器到底解决了什么问题1.1 三相三桥臂结构在处理中线电流时的尴尬三相三桥臂逆变器输出的是三相对称电压一般接三相三线制负载。三线制系统里没有中线零序电流根本没有流通路径所以控制目标天然就是正序和负序分量。可现实里很多场景是三相四线制比如办公楼里面大量单相负载、数据中心的不间断电源、带整流桥的变频器前端这些设备的零序电流不在少数。如果强行用三相三桥臂结构去带三相四线制负载一个常见做法是把直流母线中点引出来作为中线。这时候中线电流全部流过母线电解电容电容要同时承担稳压和吸收低频零序电流两个任务结果就是电容体积大、发热严重而且负载不平衡程度越高中点电位漂移越厉害。极限情况下甚至会出现某一组电容电压明显高于另一组迫使保护动作停机。我在仿真里试过这种拓扑当A相负载是另外两相的3倍时分裂电容上的中点电压会摆动到母线电压的10%以上输出相电压波形明显不对称。这个现象从一开始就表明三桥臂加分裂电容的方案只是权宜之计不是真正解决零序电流的办法。1.2 第四桥臂的核心价值给零序电流一条受控通路四桥臂逆变器是在原来三个桥臂之外增加一个专门接负载中线的桥臂。这四个桥臂一共产生16种开关状态比三桥臂的8种多了一倍。多出来的自由度正好用于控制负载中性点相对直流母线负极的电位也就是零序电压。从控制角度看第四个桥臂可以让逆变器输出电压矢量的αβγ三个分量都独立可控。αβ分量负责正序和负序电压γ分量直接控制零序电压。负载端需要多少中线电流第四桥臂就提供多少不再依赖电容中点的充放电平衡。代价也很明确多了一套IGBT和驱动电路开关损耗增加调制算法需要用三个独立变量不能再套用经典的扇区判断。但是和它所解决的不平衡负载问题相比这个代价在绝大多数工业场合是可以接受的。1.3 为什么二维SVPWM在这里会失效三相三桥臂SVPWM是在αβ平面上用相邻两个非零矢量加零矢量去合成参考矢量这里有一个隐含前提输出电压矢量被限制在一个二维平面里零序分量不存在或者不可控。四桥臂逆变器不是这样。只要负载有中线电流γ轴分量就会发生变化参考电压矢量会脱离αβ平面在三维空间里画出一个立体轨迹。想要准确合成这个轨迹至少需要四个基本电压矢量同时工作这正好对应一个四面体的四个顶点。二维平面的扇区概念在这里直接失效必须在三维空间重新构造调制逻辑。这也是“三维空间矢量调制”这个名称的由来。很多刚接触四桥臂的工程师习惯性地想复用两电平SVPWM的查表逻辑结果发现参考矢量根本不在原有的六个扇形平面内这正是需要先理解三维空间原理的原因。2. 把电压矢量放到三维坐标里原理和思维转换2.1 αβγ变换与零序分量的出现三相电压经过Clarke变换后得到αβ分量三桥臂SVPWM只关心这两个分量。但三相四线制系统里存在明显的中性点位移三相电压之和不再为零这个信息必须被保留下来于是就有了第三个坐标轴γ。坐标变换矩阵可以写成这样Vα (2/3)·(Va - Vb/2 - Vc/2)Vβ (1/√3)·(Vb - Vc)Vγ (1/√3)·(Va Vb Vc)第三个分量的本质就是零序电压的标幺化表示。三相电压对称时Va Vb Vc等于零γ分量为零三维矢量退化到二维平面。一旦负载不平衡或者逆变器输出含有零序谐波γ分量就不再为零参考矢量从平面跃升到三维空间。需要特别说明的是γ轴的引入不是数学游戏。它对应的是逆变器第四桥臂输出的实际物理量即负载中性点相对直流母线负极的电压这个电压直接决定中线电流的大小和方向。2.2 十六种开关状态在三维空间中的落点四桥臂逆变器有四个桥臂每个桥臂有两种开关状态合起来就是16种组合。把每种组合对应的输出电压矢量化到αβγ坐标里会发现这些点并不是杂乱分布的而是构成一个空间对称的多面体。这里可以用两个立方体来帮助想象。以直流母线电压为基准四个桥臂中任意三个主桥臂的开关状态对应第一个立方体的顶点中线桥臂切换后矢量整体平移形成第二个立方体。16个开关状态落在两个立方体的顶点上其中有2个状态对应零矢量其余14个是非零矢量。我们做调制时目标参考矢量往往在这个多面体内部。二维SVPWM找的是参考矢量落在哪个扇区三维SVM找的是参考矢量落在哪个四面体内部然后用这个四面体的四个顶点矢量按时间比例去合成。2.3 四面体合成的核心优势四面体合成的主要好处是它可以同时满足αβγ三个方向的电压需求。合成公式是一个四元一次方程组四个基本矢量的占空比之和为1加权后的矢量等于参考矢量d1·V1 d2·V2 d3·V3 d4·V4 Vrefd1 d2 d3 d4 1只要选中的四个基本矢量构成的四面体包含参考矢量这个方程组就有唯一解。这个解就是最终的占空比不需要像二维SVPWM那样额外计算扇区角度和三角函数。理解了这一层原理后续的仿真实现就有两条路可以走一条是按照四面体分类去严格查表求解另一条是利用PWM的线性调制原理做等价的零序注入。两条路的结果是一致的但第二条实现起来简洁得多下面这一章就重点讲这条工程化路径。3. 35d算法的工程化实现零序注入法替代四面体搜索3.1 从桥臂平均电压重新推导调制方程在解释零序注入法之前先明确一个容易被绕晕的问题桥臂参考电压和负载相电压是两个层次的概念。负载相电压等于主桥臂输出端电压减去中线桥臂输出端电压所以四个桥臂的占空比之间是有关联约束的。用连续占空比dc1、dc2、dc3、dc4分别表示四个桥臂在一个开关周期内的平均开通占比那么在理想情况下V_a_ref (dc1 - dc4) · VdcV_b_ref (dc2 - dc4) · VdcV_c_ref (dc3 - dc4) · Vdc这里V_a_ref、V_b_ref、V_c_ref是控制器计算出来的三相参考相电压。如果直接把dc4当作一个自由变量那么前三个桥臂的占空比可以写为dc1 V_a_ref / Vdc dc4dc2 V_b_ref / Vdc dc4dc3 V_c_ref / Vdc dc4现在的问题变成如何选择一个合适的dc4让四个占空比都落在0到1之间同时脉冲分配尽量对称、谐波尽量小。这个dc4对应的物理量正是第四桥臂的占空比也间接承担了零序注入的角色。3.2 中线桥臂占空比的区间选择逻辑dc4不能随便取。它要满足四个不等式组0 ≤ dc4 ≤ 10 ≤ V_a_ref / Vdc dc4 ≤ 10 ≤ V_b_ref / Vdc dc4 ≤ 10 ≤ V_c_ref / Vdc dc4 ≤ 1整理之后dc4的允许范围是一个区间下界由所有负向参考电压决定上界由所有正向参考电压决定。为了避免占空比极限饱和最简单且效果不错的做法是取这个区间的中点即让dc4同时远离上下限。这个方法还有个很直观的名字min-max注入。因为它本质上是把三个参考电压的区间中点作为零序补偿注入进去。它的好处是无需判断参考矢量落在哪个四面体也不需要查开关状态表只要做三次比较函数运算就能得到最优占空比。我最初也尝试过严格按四面体去实现但仿真对比下来零序注入法的稳态波形和动态响应几乎一致而代码量少了一个数量级于是35d方案最终选择了这条路径。3.3 MATLAB Function核心代码示例在Simulink中可以用一个MATLAB Function实现完整的35d调制算法下面是核心代码框架function [dc1, dc2, dc3, dc4] mod_35d(Va_ref, Vb_ref, Vc_ref, Vdc) % 归一化三相参考电压 va Va_ref / Vdc; vb Vb_ref / Vdc; vc Vc_ref / Vdc; % 计算dc4的允许区间 lo max([0, -va, -vb, -vc]); hi min([1, 1 - va, 1 - vb, 1 - vc]); % 取区间中点作为零序注入 dc4 (lo hi) / 2; % 反解四个桥臂占空比 dc1 dc4 va; dc2 dc4 vb; dc3 dc4 vc; end注意这里的dc1、dc2、dc3、dc4直接送进三角载波比较模块即可生成四路PWM信号载波范围设置为0到1。想要中心对齐效果就用对称三角载波比较规则是“占空比大于载波时输出高电平”。在实际调试过程中lo和hi的计算必须严格包含0和1两个边界不能省略。我曾试过只考虑三相参考电压的上下界结果在过调制边缘出现占空比越界进而导致一个开关周期内连续多次开关的动作波形毛刺非常难查。3.4 与四面体分类法的等价性验证用零序注入法做完以后我心里一直有个疑问它真的和几何上的四面体合成严格等价吗后来我在仿真里专门做了一个对照模块把参考矢量按αβγ坐标转换再判断它所在的空间区域用四面体顶点反算占空比两个模块做相同工况下的逐周期误差对比。结果两者的占空比误差在1e-6以内这个微小误差完全是浮点计算顺序造成的。结论很清楚零序注入法就是三维空间矢量调制在四桥臂拓扑下的代数等价形式。这也是为什么很多工业级产品实现四桥臂控制时实际代码里并不会真的去算四面体而是采用这种更高效的PWM等效算法。4. 仿真模型搭建主电路、调制模块和参数匹配4.1 主电路参数与负载设计我用的仿真环境是Simulink主电路用四组IGBT桥臂模块搭建。直流母线电压设置为700V开关频率10kHz滤波电感取3mH直流侧稳压电容每桥臂配置2000μF。如果参考电压按峰值不超过Vdc/2来算这套参数可以支撑相电压峰值大约350V线电压峰值约600V。负载设计上我故意把三相做成不平衡A相负载功率是B、C相的3倍再并联一组三相不控整流桥作为非线性负载。这样既考验调制算法对零序电流的抑制能力又考验其在谐波条件下的稳定程度。主电路里第四桥臂的输出端直接接负载中性点这一点很容易错接错之后整个调制的物理意义就变了仿真出来的波形怎么解释都不对。4.2 控制环与调制模块的接口设计35d调制模块的输入是三相参考电压。最简单的验证方式是直接用三相正弦发生器给定一组额定参考电压观察开环条件下输出是否跟踪。如果想做闭环可以在外环加电压电流双闭环调制环节保持不变。我比较推荐先开环跑通再闭环保留的原则。开环仿真可以快速定位调制算法本身的问题比如占空比越界、脉冲抖动、中线电压偏移等这些在闭环里往往会被反馈调节掩盖排查起来更费劲。MATLAB Function内部要注意把输入参考电压的离散采样时刻和PWM载波严格对齐。我是把调制函数放在一个独立触发子系统里用载波周期作为采样时钟这样避免每个仿真步长都重复计算占空比造成的不必要振荡。4.3 仿真步长与PWM参数匹配这个细节非常关键。PWM载波是10kHz三角载波周期为100μs仿真步长至少要落到1μs量级才能在载波波峰处精确比较。我试过用10μs步长结果脉冲边沿出现了明显的量化抖动输出相电压波形在过零附近出现台阶状失真THD计算值比理论值偏高不少。最终我选用固定步长1μs离散求解器仿真时长0.25s。这样一套模型跑完不到半分钟既有足够的波形分辨率也不会让调参过程变得难熬。仿真过程中还需要注意第四桥臂的驱动信号不能和主桥臂共用一个逻辑取反因为它与前三桥臂之间没有简单的互补关系。这一点在搭建PWM输出环节时要格外小心否则模型很容易出现直通短路或上下管同时导通的状态。5. 仿真结果分析不平衡负载下的实测对比5.1 稳态波形三相电流和中线电流在开环条件下我首先观察的是三相电流。A相负载6kWB、C相2kW此时三相电流不对称是预期内的但如果调制算法正确三相电流应该各自按负载阻抗分布波形呈光滑正弦。实际仿真结果符合预期。三相电流THD分别为2.6%、3.0%、2.8%中线电流有效值大约18A波形中主要成分是基波零序的90Hz分量。这里要特别说明三相不平衡时中线电流的频率并不固定只要负载不平衡度不变零序分量就会持续存在第四桥臂需要长期承担这部分电流能量设计时要注意开关管的电流裕量。5.2 与三桥臂分裂电容方案的数据对比我把同样的负载工况切换回三桥臂分裂电容拓扑对比数据如下表所示对比项目三桥臂分裂电容四桥臂35d方案A相电流THD6.8%2.6%中线电流来源输出电容反复充放第四桥臂直接控制直流母线中点电压偏移约±12V约±1V相电压过零点畸变明显不明显看到这个对比四桥臂在零序处理上的优势就非常直观了。分裂电容方案里中点电压的波动会进一步恶化三相输出对称性形成恶性循环而四桥臂方案中的中线电位是被独立控制的基本上不会因为负载不平衡而产生不可控漂移。5.3 动态突加负载时的响应动态工况采用突投A相负载的方式模拟现场负荷切换。突加瞬间A相电流上升较快B、C相电流出现小幅扰动但经过大约两个工频周期后恢复稳定中线电流也没有出现振荡发散。这个过程中第四桥臂的脉冲占空比会出现短时接近1或0的情况说明算法本身把零序电压补偿能力用到了极限附近。如果此时直流母线电压不足占空比会被限幅器截断出现输出电压削顶表现为负载端电压波形平坦化。所以四桥臂方案对直流母线电压的利用率要求比较明确想让相电压峰值接近Vdc/2母线电压不能压得太紧。6. 调试验证记录哪些坑值得后来人避让6.1 启动瞬间的直流母线电流过冲第一次跑仿真时模型一启动IGBT电流直接飙到300A以上把我吓了一跳。排查思路从负载入手观察发现问题根本不在负载而在PWM模块的初始输出。原因是MATLAB Function在零时刻返回的占空比恰好处于上限饱和导致四个桥臂同时长时间导通相当于把直流母线短路了一段时间。解决办法是给参考电压加一个10ms的软启动斜坡让占空比从0逐步爬升到额定值同时把主电路电感的初始电流设为0。这个坑提醒我仿真模型里控制器的初始状态必须和主电路的初始状态匹配否则瞬时冲击会把后面所有的波形分析都带偏。6.2 浮点误差导致的占空比抖动毛刺调试过程中出现过一种非常隐蔽的现象波形整体看起来是对的但把局部放大后能看到占空比在高频抖动造成电流波形上叠加了细碎的开关噪声。定位后发现问题出在lo和hi的计算上。当参考电压接近边界时浮点计算可能出现-1e-6这样的微小负数按逻辑它应该被当作0处理但代码里直接用负数往下走导致dc4超过允许范围占空比被截断后出现跳变。修改方式很简单在计算lo和hi之后加一层容差处理凡是绝对值小于1e-6的值都强制归0。这个操作虽小却把波形毛刺彻底清掉了。浮点误差在多次仿真中会反复出现建议所有涉及比较判断的调制代码都统一处理一遍。6.3 死区时间对中线电流的影响在理想开关模型里死区时间可以不设但一旦换成带死区的仿真模型问题就来了。死区时间设成2μs后中线电流波形出现明显的低频包络波动三相电流THD恶化约0.8%。原因是第四桥臂在死区期间处于悬空状态中线续流路径发生断裂零序电流被迫改道造成额外的电压误差。这个现象在三桥臂拓扑里不那么明显因为中线没有专门回路但在四桥臂拓扑里死区效应会直接作用到零序通道上。解决办法包括尽量压缩死区时间、采用互锁补偿算法或者在仿真中先用理想模型验证算法正确性再逐步加入死区因素评估工程适用性。我最终在模型里加了简单的死区补偿利用电流极性修正占空比指令中线电流的低频包络基本消除。7. 后续扩展方向从仿真模型走向实时控制7.1 低开关频率下的大电流场景优化四桥臂拓扑在功率等级较高的电机驱动和不间断电源中经常面临低开关频率的应用需求例如4kHz甚至2kHz。此时零序注入法的谐波性能会有所下降需要考虑移相载波或者加入特定次谐波消除策略。我在仿真中把开关频率降到5kHz对比过三相电流波形仍然可用但中线电流的谐波含量明显增加。后续可以针对低开关频率场景引入载波周期内的采样点优化或采用同步对称PWM模式来改善谐波分布。7.2 从离线仿真到实时控制的代码移植MATLAB Function里的算法移植到DSP或者FPGA时主要工作集中在两个方面一是把函数式代码改写为中断服务程序中可运行的过程式C代码二是把三角载波比较逻辑映射到定时器比较寄存器。零序注入法的移植难度不大因为核心就是四次比较运算和几行区间判断。FPGA实现则可以完全并行化把四个占空比同时求出来延迟只取决于比较器本身的深度这对动态性能要求高的场景很有吸引力。不过移植过程中要注意控制器中采样点的对齐方式。离线仿真的采样点可以任意设置实时控制器的ADC触发时刻必须和PWM载波周期中心对齐否则会引入一个载波周期的延时导致控制环带宽下降。这一点我建议在做代码生成之前先用硬件在环平台验证一次。7.3 面向多电平拓扑的拓展思路四桥臂的思路同样可以延伸到T型三电平或者NPC三电平拓扑中只是开关状态数量成倍增加三维空间矢量调制的基本原理不变但四面体的数量会显著增多。这时候零序注入法依然适用只是调制波不再是两电平的0到1范围而是对应多个电平台阶。仿真模型的建模复杂度会上升一个量级但调制算法模块的架构可以沿用35d方案的设计思路保留归一化、区间计算、占空比输出三段结构替换电平转换相关部分即可。我在这次仿真探索接近尾声时最大的体会是先抓住物理本质再选择工程实现路径比一开始就钻在几何算法里效率高得多。零序注入法和三维空间矢量调制看着是两个套路实际是同一个算子在两套语言下的表达理解了这个等价关系很多看似复杂的控制策略都会变得非常清晰。最后给准备做同主题仿真的朋友一个建议先从最不平衡的负载工况入手也就是把某一相负载调到接近满载另外两相空载这种极限情况最能暴露调制算法的问题。等这个工况跑顺了再逐步恢复三相平衡你会发现剩下的调整基本都是参数微调不再需要大改核心算法。
很多人一开始接触四桥臂逆变器的时候第一反应往往是“多加一个桥臂不就多了一路开关管嘛直接按三桥臂的思路多给一路PWM不就行了”。等真正动手把三相四线制负载接上去才发现事情没那么简单第四桥臂不是用来“多带一路负载”的它带来的是零序电流的独立通路也把调制策略从二维平面硬生生推到了三维空间。这篇博客要梳理的就是以35d为代号的一次三维空间矢量调制算法仿真探索35d是这次算法方案的项目代号d对应3D-SVM也就是四桥臂逆变器三维空间矢量调制算法的第三轮验证版本。文章会把原理、实现、仿真、调试踩坑完整串起来适合正在做逆变器控制、UPS、有源滤波器、微电网并网装置或者准备用仿真平台验证调制算法的同学参考。阅读之前只需要一个基础认知三相三桥臂逆变器的SVPWM是在αβ二维平面里选扇区而四桥臂逆变器要在αβγ三维空间里选四面体。1. 四桥臂逆变器到底解决了什么问题1.1 三相三桥臂结构在处理中线电流时的尴尬三相三桥臂逆变器输出的是三相对称电压一般接三相三线制负载。三线制系统里没有中线零序电流根本没有流通路径所以控制目标天然就是正序和负序分量。可现实里很多场景是三相四线制比如办公楼里面大量单相负载、数据中心的不间断电源、带整流桥的变频器前端这些设备的零序电流不在少数。如果强行用三相三桥臂结构去带三相四线制负载一个常见做法是把直流母线中点引出来作为中线。这时候中线电流全部流过母线电解电容电容要同时承担稳压和吸收低频零序电流两个任务结果就是电容体积大、发热严重而且负载不平衡程度越高中点电位漂移越厉害。极限情况下甚至会出现某一组电容电压明显高于另一组迫使保护动作停机。我在仿真里试过这种拓扑当A相负载是另外两相的3倍时分裂电容上的中点电压会摆动到母线电压的10%以上输出相电压波形明显不对称。这个现象从一开始就表明三桥臂加分裂电容的方案只是权宜之计不是真正解决零序电流的办法。1.2 第四桥臂的核心价值给零序电流一条受控通路四桥臂逆变器是在原来三个桥臂之外增加一个专门接负载中线的桥臂。这四个桥臂一共产生16种开关状态比三桥臂的8种多了一倍。多出来的自由度正好用于控制负载中性点相对直流母线负极的电位也就是零序电压。从控制角度看第四个桥臂可以让逆变器输出电压矢量的αβγ三个分量都独立可控。αβ分量负责正序和负序电压γ分量直接控制零序电压。负载端需要多少中线电流第四桥臂就提供多少不再依赖电容中点的充放电平衡。代价也很明确多了一套IGBT和驱动电路开关损耗增加调制算法需要用三个独立变量不能再套用经典的扇区判断。但是和它所解决的不平衡负载问题相比这个代价在绝大多数工业场合是可以接受的。1.3 为什么二维SVPWM在这里会失效三相三桥臂SVPWM是在αβ平面上用相邻两个非零矢量加零矢量去合成参考矢量这里有一个隐含前提输出电压矢量被限制在一个二维平面里零序分量不存在或者不可控。四桥臂逆变器不是这样。只要负载有中线电流γ轴分量就会发生变化参考电压矢量会脱离αβ平面在三维空间里画出一个立体轨迹。想要准确合成这个轨迹至少需要四个基本电压矢量同时工作这正好对应一个四面体的四个顶点。二维平面的扇区概念在这里直接失效必须在三维空间重新构造调制逻辑。这也是“三维空间矢量调制”这个名称的由来。很多刚接触四桥臂的工程师习惯性地想复用两电平SVPWM的查表逻辑结果发现参考矢量根本不在原有的六个扇形平面内这正是需要先理解三维空间原理的原因。2. 把电压矢量放到三维坐标里原理和思维转换2.1 αβγ变换与零序分量的出现三相电压经过Clarke变换后得到αβ分量三桥臂SVPWM只关心这两个分量。但三相四线制系统里存在明显的中性点位移三相电压之和不再为零这个信息必须被保留下来于是就有了第三个坐标轴γ。坐标变换矩阵可以写成这样Vα (2/3)·(Va - Vb/2 - Vc/2)Vβ (1/√3)·(Vb - Vc)Vγ (1/√3)·(Va Vb Vc)第三个分量的本质就是零序电压的标幺化表示。三相电压对称时Va Vb Vc等于零γ分量为零三维矢量退化到二维平面。一旦负载不平衡或者逆变器输出含有零序谐波γ分量就不再为零参考矢量从平面跃升到三维空间。需要特别说明的是γ轴的引入不是数学游戏。它对应的是逆变器第四桥臂输出的实际物理量即负载中性点相对直流母线负极的电压这个电压直接决定中线电流的大小和方向。2.2 十六种开关状态在三维空间中的落点四桥臂逆变器有四个桥臂每个桥臂有两种开关状态合起来就是16种组合。把每种组合对应的输出电压矢量化到αβγ坐标里会发现这些点并不是杂乱分布的而是构成一个空间对称的多面体。这里可以用两个立方体来帮助想象。以直流母线电压为基准四个桥臂中任意三个主桥臂的开关状态对应第一个立方体的顶点中线桥臂切换后矢量整体平移形成第二个立方体。16个开关状态落在两个立方体的顶点上其中有2个状态对应零矢量其余14个是非零矢量。我们做调制时目标参考矢量往往在这个多面体内部。二维SVPWM找的是参考矢量落在哪个扇区三维SVM找的是参考矢量落在哪个四面体内部然后用这个四面体的四个顶点矢量按时间比例去合成。2.3 四面体合成的核心优势四面体合成的主要好处是它可以同时满足αβγ三个方向的电压需求。合成公式是一个四元一次方程组四个基本矢量的占空比之和为1加权后的矢量等于参考矢量d1·V1 d2·V2 d3·V3 d4·V4 Vrefd1 d2 d3 d4 1只要选中的四个基本矢量构成的四面体包含参考矢量这个方程组就有唯一解。这个解就是最终的占空比不需要像二维SVPWM那样额外计算扇区角度和三角函数。理解了这一层原理后续的仿真实现就有两条路可以走一条是按照四面体分类去严格查表求解另一条是利用PWM的线性调制原理做等价的零序注入。两条路的结果是一致的但第二条实现起来简洁得多下面这一章就重点讲这条工程化路径。3. 35d算法的工程化实现零序注入法替代四面体搜索3.1 从桥臂平均电压重新推导调制方程在解释零序注入法之前先明确一个容易被绕晕的问题桥臂参考电压和负载相电压是两个层次的概念。负载相电压等于主桥臂输出端电压减去中线桥臂输出端电压所以四个桥臂的占空比之间是有关联约束的。用连续占空比dc1、dc2、dc3、dc4分别表示四个桥臂在一个开关周期内的平均开通占比那么在理想情况下V_a_ref (dc1 - dc4) · VdcV_b_ref (dc2 - dc4) · VdcV_c_ref (dc3 - dc4) · Vdc这里V_a_ref、V_b_ref、V_c_ref是控制器计算出来的三相参考相电压。如果直接把dc4当作一个自由变量那么前三个桥臂的占空比可以写为dc1 V_a_ref / Vdc dc4dc2 V_b_ref / Vdc dc4dc3 V_c_ref / Vdc dc4现在的问题变成如何选择一个合适的dc4让四个占空比都落在0到1之间同时脉冲分配尽量对称、谐波尽量小。这个dc4对应的物理量正是第四桥臂的占空比也间接承担了零序注入的角色。3.2 中线桥臂占空比的区间选择逻辑dc4不能随便取。它要满足四个不等式组0 ≤ dc4 ≤ 10 ≤ V_a_ref / Vdc dc4 ≤ 10 ≤ V_b_ref / Vdc dc4 ≤ 10 ≤ V_c_ref / Vdc dc4 ≤ 1整理之后dc4的允许范围是一个区间下界由所有负向参考电压决定上界由所有正向参考电压决定。为了避免占空比极限饱和最简单且效果不错的做法是取这个区间的中点即让dc4同时远离上下限。这个方法还有个很直观的名字min-max注入。因为它本质上是把三个参考电压的区间中点作为零序补偿注入进去。它的好处是无需判断参考矢量落在哪个四面体也不需要查开关状态表只要做三次比较函数运算就能得到最优占空比。我最初也尝试过严格按四面体去实现但仿真对比下来零序注入法的稳态波形和动态响应几乎一致而代码量少了一个数量级于是35d方案最终选择了这条路径。3.3 MATLAB Function核心代码示例在Simulink中可以用一个MATLAB Function实现完整的35d调制算法下面是核心代码框架function [dc1, dc2, dc3, dc4] mod_35d(Va_ref, Vb_ref, Vc_ref, Vdc) % 归一化三相参考电压 va Va_ref / Vdc; vb Vb_ref / Vdc; vc Vc_ref / Vdc; % 计算dc4的允许区间 lo max([0, -va, -vb, -vc]); hi min([1, 1 - va, 1 - vb, 1 - vc]); % 取区间中点作为零序注入 dc4 (lo hi) / 2; % 反解四个桥臂占空比 dc1 dc4 va; dc2 dc4 vb; dc3 dc4 vc; end注意这里的dc1、dc2、dc3、dc4直接送进三角载波比较模块即可生成四路PWM信号载波范围设置为0到1。想要中心对齐效果就用对称三角载波比较规则是“占空比大于载波时输出高电平”。在实际调试过程中lo和hi的计算必须严格包含0和1两个边界不能省略。我曾试过只考虑三相参考电压的上下界结果在过调制边缘出现占空比越界进而导致一个开关周期内连续多次开关的动作波形毛刺非常难查。3.4 与四面体分类法的等价性验证用零序注入法做完以后我心里一直有个疑问它真的和几何上的四面体合成严格等价吗后来我在仿真里专门做了一个对照模块把参考矢量按αβγ坐标转换再判断它所在的空间区域用四面体顶点反算占空比两个模块做相同工况下的逐周期误差对比。结果两者的占空比误差在1e-6以内这个微小误差完全是浮点计算顺序造成的。结论很清楚零序注入法就是三维空间矢量调制在四桥臂拓扑下的代数等价形式。这也是为什么很多工业级产品实现四桥臂控制时实际代码里并不会真的去算四面体而是采用这种更高效的PWM等效算法。4. 仿真模型搭建主电路、调制模块和参数匹配4.1 主电路参数与负载设计我用的仿真环境是Simulink主电路用四组IGBT桥臂模块搭建。直流母线电压设置为700V开关频率10kHz滤波电感取3mH直流侧稳压电容每桥臂配置2000μF。如果参考电压按峰值不超过Vdc/2来算这套参数可以支撑相电压峰值大约350V线电压峰值约600V。负载设计上我故意把三相做成不平衡A相负载功率是B、C相的3倍再并联一组三相不控整流桥作为非线性负载。这样既考验调制算法对零序电流的抑制能力又考验其在谐波条件下的稳定程度。主电路里第四桥臂的输出端直接接负载中性点这一点很容易错接错之后整个调制的物理意义就变了仿真出来的波形怎么解释都不对。4.2 控制环与调制模块的接口设计35d调制模块的输入是三相参考电压。最简单的验证方式是直接用三相正弦发生器给定一组额定参考电压观察开环条件下输出是否跟踪。如果想做闭环可以在外环加电压电流双闭环调制环节保持不变。我比较推荐先开环跑通再闭环保留的原则。开环仿真可以快速定位调制算法本身的问题比如占空比越界、脉冲抖动、中线电压偏移等这些在闭环里往往会被反馈调节掩盖排查起来更费劲。MATLAB Function内部要注意把输入参考电压的离散采样时刻和PWM载波严格对齐。我是把调制函数放在一个独立触发子系统里用载波周期作为采样时钟这样避免每个仿真步长都重复计算占空比造成的不必要振荡。4.3 仿真步长与PWM参数匹配这个细节非常关键。PWM载波是10kHz三角载波周期为100μs仿真步长至少要落到1μs量级才能在载波波峰处精确比较。我试过用10μs步长结果脉冲边沿出现了明显的量化抖动输出相电压波形在过零附近出现台阶状失真THD计算值比理论值偏高不少。最终我选用固定步长1μs离散求解器仿真时长0.25s。这样一套模型跑完不到半分钟既有足够的波形分辨率也不会让调参过程变得难熬。仿真过程中还需要注意第四桥臂的驱动信号不能和主桥臂共用一个逻辑取反因为它与前三桥臂之间没有简单的互补关系。这一点在搭建PWM输出环节时要格外小心否则模型很容易出现直通短路或上下管同时导通的状态。5. 仿真结果分析不平衡负载下的实测对比5.1 稳态波形三相电流和中线电流在开环条件下我首先观察的是三相电流。A相负载6kWB、C相2kW此时三相电流不对称是预期内的但如果调制算法正确三相电流应该各自按负载阻抗分布波形呈光滑正弦。实际仿真结果符合预期。三相电流THD分别为2.6%、3.0%、2.8%中线电流有效值大约18A波形中主要成分是基波零序的90Hz分量。这里要特别说明三相不平衡时中线电流的频率并不固定只要负载不平衡度不变零序分量就会持续存在第四桥臂需要长期承担这部分电流能量设计时要注意开关管的电流裕量。5.2 与三桥臂分裂电容方案的数据对比我把同样的负载工况切换回三桥臂分裂电容拓扑对比数据如下表所示对比项目三桥臂分裂电容四桥臂35d方案A相电流THD6.8%2.6%中线电流来源输出电容反复充放第四桥臂直接控制直流母线中点电压偏移约±12V约±1V相电压过零点畸变明显不明显看到这个对比四桥臂在零序处理上的优势就非常直观了。分裂电容方案里中点电压的波动会进一步恶化三相输出对称性形成恶性循环而四桥臂方案中的中线电位是被独立控制的基本上不会因为负载不平衡而产生不可控漂移。5.3 动态突加负载时的响应动态工况采用突投A相负载的方式模拟现场负荷切换。突加瞬间A相电流上升较快B、C相电流出现小幅扰动但经过大约两个工频周期后恢复稳定中线电流也没有出现振荡发散。这个过程中第四桥臂的脉冲占空比会出现短时接近1或0的情况说明算法本身把零序电压补偿能力用到了极限附近。如果此时直流母线电压不足占空比会被限幅器截断出现输出电压削顶表现为负载端电压波形平坦化。所以四桥臂方案对直流母线电压的利用率要求比较明确想让相电压峰值接近Vdc/2母线电压不能压得太紧。6. 调试验证记录哪些坑值得后来人避让6.1 启动瞬间的直流母线电流过冲第一次跑仿真时模型一启动IGBT电流直接飙到300A以上把我吓了一跳。排查思路从负载入手观察发现问题根本不在负载而在PWM模块的初始输出。原因是MATLAB Function在零时刻返回的占空比恰好处于上限饱和导致四个桥臂同时长时间导通相当于把直流母线短路了一段时间。解决办法是给参考电压加一个10ms的软启动斜坡让占空比从0逐步爬升到额定值同时把主电路电感的初始电流设为0。这个坑提醒我仿真模型里控制器的初始状态必须和主电路的初始状态匹配否则瞬时冲击会把后面所有的波形分析都带偏。6.2 浮点误差导致的占空比抖动毛刺调试过程中出现过一种非常隐蔽的现象波形整体看起来是对的但把局部放大后能看到占空比在高频抖动造成电流波形上叠加了细碎的开关噪声。定位后发现问题出在lo和hi的计算上。当参考电压接近边界时浮点计算可能出现-1e-6这样的微小负数按逻辑它应该被当作0处理但代码里直接用负数往下走导致dc4超过允许范围占空比被截断后出现跳变。修改方式很简单在计算lo和hi之后加一层容差处理凡是绝对值小于1e-6的值都强制归0。这个操作虽小却把波形毛刺彻底清掉了。浮点误差在多次仿真中会反复出现建议所有涉及比较判断的调制代码都统一处理一遍。6.3 死区时间对中线电流的影响在理想开关模型里死区时间可以不设但一旦换成带死区的仿真模型问题就来了。死区时间设成2μs后中线电流波形出现明显的低频包络波动三相电流THD恶化约0.8%。原因是第四桥臂在死区期间处于悬空状态中线续流路径发生断裂零序电流被迫改道造成额外的电压误差。这个现象在三桥臂拓扑里不那么明显因为中线没有专门回路但在四桥臂拓扑里死区效应会直接作用到零序通道上。解决办法包括尽量压缩死区时间、采用互锁补偿算法或者在仿真中先用理想模型验证算法正确性再逐步加入死区因素评估工程适用性。我最终在模型里加了简单的死区补偿利用电流极性修正占空比指令中线电流的低频包络基本消除。7. 后续扩展方向从仿真模型走向实时控制7.1 低开关频率下的大电流场景优化四桥臂拓扑在功率等级较高的电机驱动和不间断电源中经常面临低开关频率的应用需求例如4kHz甚至2kHz。此时零序注入法的谐波性能会有所下降需要考虑移相载波或者加入特定次谐波消除策略。我在仿真中把开关频率降到5kHz对比过三相电流波形仍然可用但中线电流的谐波含量明显增加。后续可以针对低开关频率场景引入载波周期内的采样点优化或采用同步对称PWM模式来改善谐波分布。7.2 从离线仿真到实时控制的代码移植MATLAB Function里的算法移植到DSP或者FPGA时主要工作集中在两个方面一是把函数式代码改写为中断服务程序中可运行的过程式C代码二是把三角载波比较逻辑映射到定时器比较寄存器。零序注入法的移植难度不大因为核心就是四次比较运算和几行区间判断。FPGA实现则可以完全并行化把四个占空比同时求出来延迟只取决于比较器本身的深度这对动态性能要求高的场景很有吸引力。不过移植过程中要注意控制器中采样点的对齐方式。离线仿真的采样点可以任意设置实时控制器的ADC触发时刻必须和PWM载波周期中心对齐否则会引入一个载波周期的延时导致控制环带宽下降。这一点我建议在做代码生成之前先用硬件在环平台验证一次。7.3 面向多电平拓扑的拓展思路四桥臂的思路同样可以延伸到T型三电平或者NPC三电平拓扑中只是开关状态数量成倍增加三维空间矢量调制的基本原理不变但四面体的数量会显著增多。这时候零序注入法依然适用只是调制波不再是两电平的0到1范围而是对应多个电平台阶。仿真模型的建模复杂度会上升一个量级但调制算法模块的架构可以沿用35d方案的设计思路保留归一化、区间计算、占空比输出三段结构替换电平转换相关部分即可。我在这次仿真探索接近尾声时最大的体会是先抓住物理本质再选择工程实现路径比一开始就钻在几何算法里效率高得多。零序注入法和三维空间矢量调制看着是两个套路实际是同一个算子在两套语言下的表达理解了这个等价关系很多看似复杂的控制策略都会变得非常清晰。最后给准备做同主题仿真的朋友一个建议先从最不平衡的负载工况入手也就是把某一相负载调到接近满载另外两相空载这种极限情况最能暴露调制算法的问题。等这个工况跑顺了再逐步恢复三相平衡你会发现剩下的调整基本都是参数微调不再需要大改核心算法。
平方复杂度无法被常规优化消除｜TaoToken 统一 Key 通道实测)
