干并网变流器的朋友应该都有这种直觉同一个柜子、同一套控制代码在工厂强电网测试时一切正常拉到现场接到长距离架空线路后面直流母线电压开始无缘无故地波动电流波形隐约出现十几赫兹的包络振荡严重时直接报过流跳闸。这种“换个电网环境就失控”的现象往往不是硬件老化而是跟网型逆变器在小扰动下的稳定性出了问题。这篇文章围绕“跟网型逆变器小干扰稳定性分析与控制策略优化研究Simulink仿真实现”这个课题把我从失稳机理理解、分析路线选型、Simulink模型搭建到参数整定和策略优化的完整思路写下来。适合正在做并网逆变器稳定性课题的研究生、刚入行新能源变流器研发的工程师以及想弄明白“为什么弱电网下PLL会引发振荡”的人参考。整个分析链路我都会用可落地的Simulink实现方式来讲。1. 跟网型逆变器失稳从“锁不住相位”到“振荡发散”的物理画像1.1 跟网型的定位和典型控制架构先把这个概念说清楚。所谓“跟网型”Grid-FollowingGFL核心含义是逆变器本身不建立电压和频率它通过锁相环PLL去检测电网电压的相位然后把自己的输出电流相位“贴”到电网相位上。典型控制结构从上到下分三层最外层是功率控制或直流电压控制生成 dq 轴电流指令中间是电流内环通常用 dq 轴 PI 调节器生成调制波最底层是调制和驱动环节把调制波变成 PWM 脉冲。这里有个关键点跟网型逆变器所有控制都建立在“PLL 提供的角度”之上。PLL 给出一个旋转坐标系电流内环在这个坐标系里把有功电流、无功电流分别控制成给定值。换句话说逆变器没有自己的“惯性”电网是什么相位它就跟着转什么相位。这种架构在强电网下非常高效因为 PCC 电压基本由电网决定PLL 输出角度稳定、抗扰性好。但电网一弱问题就来了。1.2 弱电网下小扰动失稳的机理解读弱电网最直接的表现是线路阻抗变大我用短路比 SCR 来量化。SCR 电网短路容量 / 逆变器额定容量一般低于 5 就算比较弱了低于 2 属于极弱电网。换成电路语言就是 PCC 点到理想电压源之间串联了一个比较大的感抗 Lg。在这个条件下扰动会发生这样的链条反应假设 PCC 电压相位出现一个小偏移 → PLL 为了追踪这个偏移会调整输出角频率 → 这个角度变化直接改变 dq 旋转坐标系的朝向 → 电流内环在这个被“带偏”的坐标系下看有功电流和无功电流会出现交叉耦合 → 输出有功功率和 PCC 电压相位进一步受影响 → 又反馈给 PLL。在强电网下这条反馈链路的增益很低扰动会衰减弱电网下电网感抗很大有功功率变化对 PCC 电压相位的影响被放大反馈增益超过某个临界值原来的阻尼就变成负阻尼振荡开始发散。我给一个生活化类比跟网型逆变器就像夜路骑行的人眼睛盯着前车的尾灯保持方向。强电网时前车尾灯很亮很稳你跟着走没问题弱电网时前车尾灯的光会被你的车灯反射干扰你看到的方向里掺了自己制造的影响越修正越偏最后摔进沟里。这就是弱电网下跟网型逆变器小干扰失稳的本质——不是硬件变差了是控制环和电网阻抗形成了正反馈。2. 分析路线选型为什么把特征值法作为主线、阻抗法作为验证2.1 状态空间特征值法状态变量选取与小信号模型推导研究小干扰稳定性最正统的方法是建立系统在平衡点附近的线性化状态空间模型然后计算状态矩阵的特征值。只要所有特征值实部都为负系统在小扰动下就是渐近稳定的一旦某一对共轭复根实部穿越零点就对应某个频率的振荡从衰减变成了等幅甚至发散。跟网型逆变器的小信号模型状态变量一般这样选PLL 的积分状态和角度偏差、电流内环 PI 控制器两个积分状态、L 或 LCL 滤波器的电感电流和电容电压、外环 PI 控制器的积分状态。在这个基础上对每个状态变量求偏导得到线性化矩阵 A。手动推导这个 A 矩阵工作量很大也容易出错我在实际课题里通常先在 Simulink 中搭平均值模型再用线性化工具提取系统矩阵同时用解析推导做交叉验证两边对得上才继续往下走。在 MATLAB 里求特征值非常直接% 由解析推导或 Simulink 线性化得到状态矩阵 A sys_linear ss(A, B, C, D); eigenvalues eig(sys_linear.a); % 特征根 [wn, zeta] damp(sys_linear.a); % 自然频率与阻尼比特征值法的好处是物理意义非常清楚每一对特征根对应一个振荡模态能直接看出这个模态的频率、阻尼比以及它主要受哪些参数支配。后面做 PLL 带宽扫描、电网阻抗扫描时盯着主导特征根轨迹走就行了。2.2 阻抗法从输出阻抗相位裕度看失稳边界阻抗法是另一条成熟路线思路是把逆变器看成一个从 PCC 向逆变器侧“看进去”的阻抗 Zout(s)把电网看成理想电压源串联电网阻抗 Zg(s)稳定性就看两个阻抗的交互。跟网型逆变器在 PLL 带宽附近往往表现出负电阻特性负电阻和电网感抗串联构成的 RLC 回路等效阻尼可能是负的这个回路的特征根实部一旦变正就是失稳。阻抗法做起来比全状态空间模型轻量特别适合快速扫描参数趋势。在 Simulink 里可以用小信号扫频方法直接测阻抗在 PCC 处注入不同频率的小扰动电压测量同频率的电流响应用“电压扰动 / 电流响应”就是频域阻抗。这个结果可以直接和解析推导的 Zout(s) 做对比也方便画 Bode 图看相位裕度。特征值法和阻抗法不是二选一的关系。我自己的习惯是用解析状态空间模型做主导模态定位和参数灵敏度分析用阻抗法画稳定裕度曲线验证趋势最后全部回到 Simulink 时域仿真中用非线性模型做最终确认。三种手段的结果要能互相解释。分析方法优势局限适合场景状态空间特征值法模态归属清晰、能做参与因子和灵敏度分析建模工作量较大、模型阶数高时参数化繁琐定位失稳根源、参数协同整定阻抗法模型直观、易与测量结果对比、便于扫频获取MIMO 情形需要广义 Nyquist 判断物理意义不如特征值直观快速评估失稳边界、验证优化效果Simulink 非线性时域仿真能反映真实限幅、延时和非线性因素仅靠波形难判断稳定裕度、不能直接给出根轨迹最终验证和实验匹配2.3 Simulink 线性化能不能直接用很多刚接触这个课题的同学会直接搭一个 IGBT 开关级模型然后想说“直接在 Simulink 里线性化”。这里要泼一盆冷水开关级模型含 PWM 调制和开关器件高频通断不能直接做小信号线性化即使强行操作线性化工具也只会在一个开关周期内做平均结果往往和实际不符。正确的做法是先搭平均值模型三相桥用一个受控电压源代替控制器的输出 dq 电压指令经过反 Park 变换后直接控制受控电压源的三相电压。平均值模型保留了 PLL、电流内环、功率外环所有低频动态完全满足小干扰稳定性分析需求。对于 LCL 滤波器平均模型也可以保留全部滤波状态所以谐振峰附近的动态也跑不掉。我自己在课题里通常建两套模型一套是包含 PWM 的详细开关模型用来跑大扰动、保护动作、限幅生效等工况验证控制策略在真实工况下的表现另一套是平均值模型专门做小扰动稳定性分析、线性化提取和参数扫描。两套模型共享同一套控制结构和参数用 MATLAB 脚本批量修改参数、批量运行。3. Simulink中搭出可用于稳定性分析的跟网型逆变器模型3.1 主电路与弱电网阻抗的参数设计参数设计必须一开始就定清楚否则后面所有结果都没法比较。我以一套常用参数为例三相 380V / 50Hz逆变器额定功率 20kW直流侧电压 700VL 型滤波电感 2mH、等效电阻 0.1Ω。如果要研究 LCL 型可以设 L11mHCf10uFL20.5mHLCL 会引入谐振峰这部分动态也要进入状态空间模型。电网侧用理想三相电压源串联可调阻抗来建模。电网阻抗以感抗为主但建议保留一个较小的电阻分量否则阻尼偏小仿真结果会偏保守。电网阻抗和 SCR 的关系可以这样换算对于线电压 380V、功率 20kW 的系统如果电网阻抗模值是 Zg那么 SCR ≈ (380² / 20000) / Zg 7.22 / Zg。想要 SCR5Zg 大概取 1.44Ω想要 SCR2Zg 大概取 3.6Ω对应感抗 3.6Ω换算到 Lg ≈ 11.5mH。这是一个非常大的电感比逆变器侧滤波电感 2mH 大好几倍此时弱电网对稳定性的影响相当明显。在 Simulink 中电网阻抗用“Three-Phase Series RLC Branch”串联在理想电压源和 PCC 之间即可Rg、Lg 都做成可调变量方便用脚本循环扫描。3.2 PLL、电流内环、功率外环的模块化搭建要点控制部分的建议是尽量模块化因为后面要反复修改参数。PLL 模块在 Simulink 里可以直接用专门的“PLL (3ph)”模块也可以用 Continuous 方式自建一个同步旋转坐标系锁相环核心是 q 轴电压通过 PI 调节器锁定到零。自建的好处是你能清楚看到 PLL 的状态变量方便后续线性化和状态矩阵提取。PLL 参数按典型二阶系统设计锁相环闭环自然频率 wn_pll 取 2π×15 ≈ 94.2 rad/s阻尼比取 0.707那么 Ki_pll wn_pll²Kp_pll 2×zeta×wn_pll。注意不同文献的 PLL 模型归一化方式不同参数数值可能差异很大重点看闭环带宽而不是单个 kp、ki 的绝对值。电流内环也按带宽设计取电流环带宽 400Hz即 2π×400 ≈ 2513 rad/s。对于 L 型滤波器可以近似取 Kp_i L×ωc 2mH×2513 ≈ 5.03Ki_i R×ωc 0.1×2513 ≈ 251.3。Simulink 中电流内环要加 dq 解耦项和 PCC 电压前馈项否则 dq 轴之间会互相干扰特征值和实际系统对不上。3.3 小扰动注入实验时域波形与频谱判据模型搭好后验证稳定性的最基本手段就是小扰动注入。在电流内环的有功电流指令 i_d_ref 上叠加一个小幅值阶跃或持续几个周期的正弦扰动幅值取额定电流的 2%~5%比如额定电流约 30A 时给 1A。观察输出电流和 PCC 电压在扰动撤销后是衰减收敛、等幅振荡还是发散。光看时域波形还不够最好配合 FFT 频谱分析。弱电网失稳通常表现为一个窄带振荡频率落在 PLL 带宽附近可能是 10Hz 到 30Hz 之间。对 PCC 电压或输出电流波形做 FFT如果扰动撤销后频谱里仍然持续存在一个清晰的窄带峰值且峰值幅值不衰减那就是失稳或临界稳定的强信号。这个流程一句话总结时域看包络、频域看峰值、特征值看实部三者对上了结论才可靠。4. 参数灵敏度量化PLL带宽、电流环与电网阻抗的拉锯战4.1 PLL带宽对主导特征根的影响PLL 带宽是跟网型逆变器稳定性最敏感的参数。把 PLL 闭环带宽从 15Hz 逐步提高到 30Hz、45Hz在强电网下可能感觉不到什么差异但在 SCR2 的弱电网下主导特征根会迅速向虚轴移动甚至直接穿越到右半平面。原因是 PLL 带宽越高它对 PCC 电压相位变化的响应越快但这意味着弱电网电压降落引起的相位变化能更快地通过 PLL 进入控制环路。这条通路的增益和电网阻抗成正比电网越弱高带宽 PLL 的副作用越大。在我的仿真里SCR2 时 PLL 带宽超过 25Hz 左右主导模态阻尼比就开始接近零再往上就是失稳。这给工程整定一个关键启发跟网型逆变器不能盲目追求“PLL 响应快”。强电网下快意味着动态好弱电网下快意味着容易振荡。如果现场电网阻抗可能变化PLL 带宽需要留足裕度。4.2 电流环参数与输出阻抗曲线的耦合电流内环带宽对稳定性也有显著影响。把电流环带宽从 400Hz 提高到 1000Hz输出阻抗在中频段的幅相特性会改变特别是在几赫兹到几百赫兹这个跟 PLL 模态重合的频段。如果电流环 PI 参数整得太激进电流环自身的闭环响应出现明显超调相当于在输出阻抗曲线上多出一个凸起这个凸起如果和电网感抗曲线发生不利的交点就会制造新的负阻尼区间。反过来电流环带宽偏低会削弱逆变器的动态跟随能力外环扰动时电流指令跟踪变慢也会间接改变整个闭环系统的阻尼特性。我建议以“电流环带宽PLL 带宽 ≈ 20~40 倍”作为初始配置然后在这个基础上做特征值扫描不要直接甩一个随机参数就开始仿真。4.3 SCR下降时主导模态的迁移规律把电网阻抗从 SCR10 开始逐步调到 SCR2观察主导特征根轨迹会看到一个很典型的迁移过程在 SCR10 时主导模态通常是电流环相关的模态阻尼高、振荡频率在几百赫兹量级随着 Lg 增大PLL 相关模态的频率逐渐下降向 10~30Hz 方向移动阻尼比逐步下降并成为新的主导模态SCR 降到 2 左右时这个模态实部逼近甚至超过零。我习惯把这种迁移做成一张表格每次扫描记录主导模态频率和阻尼比。下面是基于我仿真参数的示意结果具体数值会因滤波器参数和控制参数不同而偏移但趋势是稳定的。SCR主导模态振荡频率 (Hz)主导模态阻尼比稳定性状态10850.42稳定裕度充足5460.31稳定裕度下降3240.12临界接近弱阻尼216-0.04失稳包络发散看到这个趋势后优化目标就很明确了想办法在 SCR2、3 的条件下把主导模态阻尼比抬高同时尽量不牺牲 SCR10 时的动态性能。5. 控制策略优化重塑阻抗、补偿相位与参数协同整定5.1 虚拟电阻和虚拟电感最直接的重塑手段跟网型逆变器失稳的根源是输出阻抗在特定频段呈现负电阻特性那最直接的办法就是人为改变输出阻抗形状。虚拟电阻法在电流环指令或调制波中叠加与输出电流成比例的分量等效于在输出回路串入一个电阻。这个电阻能把负阻尼拉正实现非常简单在 dq 轴调制波上分别减去 R_v × i_d 和 R_v × i_q 即可。虚拟电阻取值从电网阻抗的 0.3 倍开始试观察主导特征根变化。R_v 增大阻尼会明显改善但代价是等效输出电压被“吃掉”一块负载能力下降重载时可能达不到额定电压。所以 R_v 不是越大越好以刚好恢复正阻尼并留 5~10 度相位裕度为目标。虚拟电感类似在调制波中加入与电流微分成比例的项能在频域削掉某个特定频段的阻抗凸起。但虚拟电感在数字控制中涉及微分项对噪声敏感离散化后容易引入高频抖动实际工程里要加低通滤波或观察器处理。5.2 电网电压前馈与相位补偿在Simulink里的实现电网电压前馈是改善抗扰性的常用手段把 PCC 电压检测值经过低通滤波后直接叠加到调制波上用来抵消电网电压扰动对输出电流的影响。前馈通道越完整逆变器对外表现为“电压源跟随”但弱电网下 PCC 电压本身受电流影响大前馈增益太大会把噪声和畸变也引进来反而激发高次谐波。所以前馈通道一定要串一个低通滤波器截止频率一般取电流环带宽的两倍以上。相位补偿则更精细一些。PLL、电流环、滤波器都会在特定频段引入相位滞后如果在调制波中对某个关键频率做超前校正可以抬高相位裕度。Simulink 里实现相位补偿可以用 Continuous 传递函数模块也可以在电流内环误差信号后面串一个超前-滞后环节。做完补偿后我会回到阻抗法画一次 Bode 图看 PLL 频段附近的相位曲线是否整体上移。5.3 基于灵敏度指标的参数协同整定流程单打独斗的参数优化不够系统我在课题里跑了一套参数协同整定流程先定义目标函数比如在 SCR2、3、5 条件下所有主导特征根阻尼比不小于 0.1同时最大特征根实部留 20% 以上裕度设计变量取 PLL 带宽、PLL 阻尼比、电流环带宽、虚拟电阻、虚拟电感。然后先用解析小信号模型做批量网格扫描或粒子群优化得到候选参数再放回 Simulink 平均值模型做线性化校验最后挑一组鲁棒性最好的参数做时域验证。这个流程的核心优势是把“控制参数整定”问题转化成了“特征值约束优化”问题不会出现调了这个参数丢了那个模态的情况。粒子群优化在 MATLAB 中可以直接调 Global Optimization Toolbox目标函数里每次调用小信号状态矩阵计算特征值几分钟就能跑完一轮。优化措施实现位置主要效果副作用 / 注意点降低 PLL 带宽PLL 调节器参数减弱弱电网正反馈增益动态响应变慢需和启动需求平衡虚拟电阻调制波叠加电流比例项消除负阻尼区间等效压降增加重载能力下降电网电压前馈 低通调制波前馈通道提高抗扰性前馈过强会引入高频噪声相位超前补偿电流内环或调制波环节抬高相位裕度高频段易放大噪声需配合滤波6. 实践中容易翻车的地方假失稳、线性化误差与模型陷阱6.1 三个容易误判为失稳的实际案例第一个案例是启动瞬间的“假振荡”。PLL 初始相位和三相电压实际相位不一致时启动瞬间冲击电流很大dq 轴电流波形会出现一个非常大的衰减振荡初学者容易把它当成失稳。解决方法是设置 PLL 初始相位为 0 时同步匹配三相电网电压初始相角或者在仿真启动阶段先开环预同步等 PLL 锁定后再切入闭环。第二个案例是数值步长导致的“假发散”。开关级模型在仿真步长过大时PWM 高频分量混叠到低频波形上看就像低频振荡发散。把步长从 1e-5 改成 1e-6 后波形立刻干净了和平均值模型结果也能对上。遇到可疑振荡先检查步长和求解器别急着改控制参数。第三个案例是 dq 变换参考轴约定不一致导致的“特征值和波形对不上”。Simulink 中 Park 变换有很多种约定有的用 sin 型、有的用 cos 型有的 q 轴超前 d 轴 90 度有的滞后。控制环里的 Park、反 Park、PLL 输出角度必须遵循同一套约定。我在排错时最常用的检查手段是给模型加一个已知幅值和相位的三相正序电压看 dq 轴的输出是不是符合预期这一步 10 分钟就能定位约定错误。6.2 线性化模型与开关级仿真结果为什么会差线性化模型忽略了三类因素PWM 开关纹波、控制延时、限幅非线性。其中最容易被忽略的是数字控制延时。实际控制器从采样、计算到 PWM 更新至少有一拍到两拍延时这个延时在几十赫兹到几百赫兹频段会带来明显的相位滞后。如果线性化模型里不加入等效延时环节得到的稳定边界会比实际系统乐观也就是算出来稳定时域开关模型已经振荡了。我在状态空间模型中一般用 e^{-1.5Ts·s} 做控制延时的近似Ts 是控制周期典型值 100us。在 Simulink 平均值模型里可以用 Transport Delay 模块模拟这个延时。这样线性化结果和开关模型的一致性会好很多。限幅非线性在小扰动分析里可以忽略但大扰动验证时一定要保留因为电流限幅会让系统工作点大幅偏移线性化模型完全覆盖不了这个现象。6.3 仿真实验设计的几条建议参数扫描建议至少覆盖 SCR10、5、3、2 四档每档都记录主导特征根、阻尼比和时域扰动测试结果。保存波形时既要保存 d 轴电流、q 轴电流、PCC 电压幅值也要保存 PLL 输出角频率因为失稳早期往往先在 PLL 频率上出现振荡迹象。还要注意工作点漂移问题。电网阻抗变化不仅改变动态特性也改变了静态工作点功率角、PCC 电压幅值都变了。线性化必须在新工作点附近做不能拿 SCR10 的工作点线性化结果去判断 SCR2 的稳定性。Simulink 线性化工具支持先对模型做 trim 找到平衡点再做线性化这个流程一定要走完整。我个人的一个小经验是把参数扫描、波形提取、特征值计算全部封装成 MATLAB 脚本用结构体统一管理参数跑 10 组参数用不了一分钟。这样你可以大胆地试各种 PLL 带宽、虚拟阻抗组合不会因为手动改参数改到怀疑人生。课题做到后期这套脚本就是最值钱的资产。
跟网型逆变器小干扰稳定性分析与Simulink仿真控制优化
干并网变流器的朋友应该都有这种直觉同一个柜子、同一套控制代码在工厂强电网测试时一切正常拉到现场接到长距离架空线路后面直流母线电压开始无缘无故地波动电流波形隐约出现十几赫兹的包络振荡严重时直接报过流跳闸。这种“换个电网环境就失控”的现象往往不是硬件老化而是跟网型逆变器在小扰动下的稳定性出了问题。这篇文章围绕“跟网型逆变器小干扰稳定性分析与控制策略优化研究Simulink仿真实现”这个课题把我从失稳机理理解、分析路线选型、Simulink模型搭建到参数整定和策略优化的完整思路写下来。适合正在做并网逆变器稳定性课题的研究生、刚入行新能源变流器研发的工程师以及想弄明白“为什么弱电网下PLL会引发振荡”的人参考。整个分析链路我都会用可落地的Simulink实现方式来讲。1. 跟网型逆变器失稳从“锁不住相位”到“振荡发散”的物理画像1.1 跟网型的定位和典型控制架构先把这个概念说清楚。所谓“跟网型”Grid-FollowingGFL核心含义是逆变器本身不建立电压和频率它通过锁相环PLL去检测电网电压的相位然后把自己的输出电流相位“贴”到电网相位上。典型控制结构从上到下分三层最外层是功率控制或直流电压控制生成 dq 轴电流指令中间是电流内环通常用 dq 轴 PI 调节器生成调制波最底层是调制和驱动环节把调制波变成 PWM 脉冲。这里有个关键点跟网型逆变器所有控制都建立在“PLL 提供的角度”之上。PLL 给出一个旋转坐标系电流内环在这个坐标系里把有功电流、无功电流分别控制成给定值。换句话说逆变器没有自己的“惯性”电网是什么相位它就跟着转什么相位。这种架构在强电网下非常高效因为 PCC 电压基本由电网决定PLL 输出角度稳定、抗扰性好。但电网一弱问题就来了。1.2 弱电网下小扰动失稳的机理解读弱电网最直接的表现是线路阻抗变大我用短路比 SCR 来量化。SCR 电网短路容量 / 逆变器额定容量一般低于 5 就算比较弱了低于 2 属于极弱电网。换成电路语言就是 PCC 点到理想电压源之间串联了一个比较大的感抗 Lg。在这个条件下扰动会发生这样的链条反应假设 PCC 电压相位出现一个小偏移 → PLL 为了追踪这个偏移会调整输出角频率 → 这个角度变化直接改变 dq 旋转坐标系的朝向 → 电流内环在这个被“带偏”的坐标系下看有功电流和无功电流会出现交叉耦合 → 输出有功功率和 PCC 电压相位进一步受影响 → 又反馈给 PLL。在强电网下这条反馈链路的增益很低扰动会衰减弱电网下电网感抗很大有功功率变化对 PCC 电压相位的影响被放大反馈增益超过某个临界值原来的阻尼就变成负阻尼振荡开始发散。我给一个生活化类比跟网型逆变器就像夜路骑行的人眼睛盯着前车的尾灯保持方向。强电网时前车尾灯很亮很稳你跟着走没问题弱电网时前车尾灯的光会被你的车灯反射干扰你看到的方向里掺了自己制造的影响越修正越偏最后摔进沟里。这就是弱电网下跟网型逆变器小干扰失稳的本质——不是硬件变差了是控制环和电网阻抗形成了正反馈。2. 分析路线选型为什么把特征值法作为主线、阻抗法作为验证2.1 状态空间特征值法状态变量选取与小信号模型推导研究小干扰稳定性最正统的方法是建立系统在平衡点附近的线性化状态空间模型然后计算状态矩阵的特征值。只要所有特征值实部都为负系统在小扰动下就是渐近稳定的一旦某一对共轭复根实部穿越零点就对应某个频率的振荡从衰减变成了等幅甚至发散。跟网型逆变器的小信号模型状态变量一般这样选PLL 的积分状态和角度偏差、电流内环 PI 控制器两个积分状态、L 或 LCL 滤波器的电感电流和电容电压、外环 PI 控制器的积分状态。在这个基础上对每个状态变量求偏导得到线性化矩阵 A。手动推导这个 A 矩阵工作量很大也容易出错我在实际课题里通常先在 Simulink 中搭平均值模型再用线性化工具提取系统矩阵同时用解析推导做交叉验证两边对得上才继续往下走。在 MATLAB 里求特征值非常直接% 由解析推导或 Simulink 线性化得到状态矩阵 A sys_linear ss(A, B, C, D); eigenvalues eig(sys_linear.a); % 特征根 [wn, zeta] damp(sys_linear.a); % 自然频率与阻尼比特征值法的好处是物理意义非常清楚每一对特征根对应一个振荡模态能直接看出这个模态的频率、阻尼比以及它主要受哪些参数支配。后面做 PLL 带宽扫描、电网阻抗扫描时盯着主导特征根轨迹走就行了。2.2 阻抗法从输出阻抗相位裕度看失稳边界阻抗法是另一条成熟路线思路是把逆变器看成一个从 PCC 向逆变器侧“看进去”的阻抗 Zout(s)把电网看成理想电压源串联电网阻抗 Zg(s)稳定性就看两个阻抗的交互。跟网型逆变器在 PLL 带宽附近往往表现出负电阻特性负电阻和电网感抗串联构成的 RLC 回路等效阻尼可能是负的这个回路的特征根实部一旦变正就是失稳。阻抗法做起来比全状态空间模型轻量特别适合快速扫描参数趋势。在 Simulink 里可以用小信号扫频方法直接测阻抗在 PCC 处注入不同频率的小扰动电压测量同频率的电流响应用“电压扰动 / 电流响应”就是频域阻抗。这个结果可以直接和解析推导的 Zout(s) 做对比也方便画 Bode 图看相位裕度。特征值法和阻抗法不是二选一的关系。我自己的习惯是用解析状态空间模型做主导模态定位和参数灵敏度分析用阻抗法画稳定裕度曲线验证趋势最后全部回到 Simulink 时域仿真中用非线性模型做最终确认。三种手段的结果要能互相解释。分析方法优势局限适合场景状态空间特征值法模态归属清晰、能做参与因子和灵敏度分析建模工作量较大、模型阶数高时参数化繁琐定位失稳根源、参数协同整定阻抗法模型直观、易与测量结果对比、便于扫频获取MIMO 情形需要广义 Nyquist 判断物理意义不如特征值直观快速评估失稳边界、验证优化效果Simulink 非线性时域仿真能反映真实限幅、延时和非线性因素仅靠波形难判断稳定裕度、不能直接给出根轨迹最终验证和实验匹配2.3 Simulink 线性化能不能直接用很多刚接触这个课题的同学会直接搭一个 IGBT 开关级模型然后想说“直接在 Simulink 里线性化”。这里要泼一盆冷水开关级模型含 PWM 调制和开关器件高频通断不能直接做小信号线性化即使强行操作线性化工具也只会在一个开关周期内做平均结果往往和实际不符。正确的做法是先搭平均值模型三相桥用一个受控电压源代替控制器的输出 dq 电压指令经过反 Park 变换后直接控制受控电压源的三相电压。平均值模型保留了 PLL、电流内环、功率外环所有低频动态完全满足小干扰稳定性分析需求。对于 LCL 滤波器平均模型也可以保留全部滤波状态所以谐振峰附近的动态也跑不掉。我自己在课题里通常建两套模型一套是包含 PWM 的详细开关模型用来跑大扰动、保护动作、限幅生效等工况验证控制策略在真实工况下的表现另一套是平均值模型专门做小扰动稳定性分析、线性化提取和参数扫描。两套模型共享同一套控制结构和参数用 MATLAB 脚本批量修改参数、批量运行。3. Simulink中搭出可用于稳定性分析的跟网型逆变器模型3.1 主电路与弱电网阻抗的参数设计参数设计必须一开始就定清楚否则后面所有结果都没法比较。我以一套常用参数为例三相 380V / 50Hz逆变器额定功率 20kW直流侧电压 700VL 型滤波电感 2mH、等效电阻 0.1Ω。如果要研究 LCL 型可以设 L11mHCf10uFL20.5mHLCL 会引入谐振峰这部分动态也要进入状态空间模型。电网侧用理想三相电压源串联可调阻抗来建模。电网阻抗以感抗为主但建议保留一个较小的电阻分量否则阻尼偏小仿真结果会偏保守。电网阻抗和 SCR 的关系可以这样换算对于线电压 380V、功率 20kW 的系统如果电网阻抗模值是 Zg那么 SCR ≈ (380² / 20000) / Zg 7.22 / Zg。想要 SCR5Zg 大概取 1.44Ω想要 SCR2Zg 大概取 3.6Ω对应感抗 3.6Ω换算到 Lg ≈ 11.5mH。这是一个非常大的电感比逆变器侧滤波电感 2mH 大好几倍此时弱电网对稳定性的影响相当明显。在 Simulink 中电网阻抗用“Three-Phase Series RLC Branch”串联在理想电压源和 PCC 之间即可Rg、Lg 都做成可调变量方便用脚本循环扫描。3.2 PLL、电流内环、功率外环的模块化搭建要点控制部分的建议是尽量模块化因为后面要反复修改参数。PLL 模块在 Simulink 里可以直接用专门的“PLL (3ph)”模块也可以用 Continuous 方式自建一个同步旋转坐标系锁相环核心是 q 轴电压通过 PI 调节器锁定到零。自建的好处是你能清楚看到 PLL 的状态变量方便后续线性化和状态矩阵提取。PLL 参数按典型二阶系统设计锁相环闭环自然频率 wn_pll 取 2π×15 ≈ 94.2 rad/s阻尼比取 0.707那么 Ki_pll wn_pll²Kp_pll 2×zeta×wn_pll。注意不同文献的 PLL 模型归一化方式不同参数数值可能差异很大重点看闭环带宽而不是单个 kp、ki 的绝对值。电流内环也按带宽设计取电流环带宽 400Hz即 2π×400 ≈ 2513 rad/s。对于 L 型滤波器可以近似取 Kp_i L×ωc 2mH×2513 ≈ 5.03Ki_i R×ωc 0.1×2513 ≈ 251.3。Simulink 中电流内环要加 dq 解耦项和 PCC 电压前馈项否则 dq 轴之间会互相干扰特征值和实际系统对不上。3.3 小扰动注入实验时域波形与频谱判据模型搭好后验证稳定性的最基本手段就是小扰动注入。在电流内环的有功电流指令 i_d_ref 上叠加一个小幅值阶跃或持续几个周期的正弦扰动幅值取额定电流的 2%~5%比如额定电流约 30A 时给 1A。观察输出电流和 PCC 电压在扰动撤销后是衰减收敛、等幅振荡还是发散。光看时域波形还不够最好配合 FFT 频谱分析。弱电网失稳通常表现为一个窄带振荡频率落在 PLL 带宽附近可能是 10Hz 到 30Hz 之间。对 PCC 电压或输出电流波形做 FFT如果扰动撤销后频谱里仍然持续存在一个清晰的窄带峰值且峰值幅值不衰减那就是失稳或临界稳定的强信号。这个流程一句话总结时域看包络、频域看峰值、特征值看实部三者对上了结论才可靠。4. 参数灵敏度量化PLL带宽、电流环与电网阻抗的拉锯战4.1 PLL带宽对主导特征根的影响PLL 带宽是跟网型逆变器稳定性最敏感的参数。把 PLL 闭环带宽从 15Hz 逐步提高到 30Hz、45Hz在强电网下可能感觉不到什么差异但在 SCR2 的弱电网下主导特征根会迅速向虚轴移动甚至直接穿越到右半平面。原因是 PLL 带宽越高它对 PCC 电压相位变化的响应越快但这意味着弱电网电压降落引起的相位变化能更快地通过 PLL 进入控制环路。这条通路的增益和电网阻抗成正比电网越弱高带宽 PLL 的副作用越大。在我的仿真里SCR2 时 PLL 带宽超过 25Hz 左右主导模态阻尼比就开始接近零再往上就是失稳。这给工程整定一个关键启发跟网型逆变器不能盲目追求“PLL 响应快”。强电网下快意味着动态好弱电网下快意味着容易振荡。如果现场电网阻抗可能变化PLL 带宽需要留足裕度。4.2 电流环参数与输出阻抗曲线的耦合电流内环带宽对稳定性也有显著影响。把电流环带宽从 400Hz 提高到 1000Hz输出阻抗在中频段的幅相特性会改变特别是在几赫兹到几百赫兹这个跟 PLL 模态重合的频段。如果电流环 PI 参数整得太激进电流环自身的闭环响应出现明显超调相当于在输出阻抗曲线上多出一个凸起这个凸起如果和电网感抗曲线发生不利的交点就会制造新的负阻尼区间。反过来电流环带宽偏低会削弱逆变器的动态跟随能力外环扰动时电流指令跟踪变慢也会间接改变整个闭环系统的阻尼特性。我建议以“电流环带宽PLL 带宽 ≈ 20~40 倍”作为初始配置然后在这个基础上做特征值扫描不要直接甩一个随机参数就开始仿真。4.3 SCR下降时主导模态的迁移规律把电网阻抗从 SCR10 开始逐步调到 SCR2观察主导特征根轨迹会看到一个很典型的迁移过程在 SCR10 时主导模态通常是电流环相关的模态阻尼高、振荡频率在几百赫兹量级随着 Lg 增大PLL 相关模态的频率逐渐下降向 10~30Hz 方向移动阻尼比逐步下降并成为新的主导模态SCR 降到 2 左右时这个模态实部逼近甚至超过零。我习惯把这种迁移做成一张表格每次扫描记录主导模态频率和阻尼比。下面是基于我仿真参数的示意结果具体数值会因滤波器参数和控制参数不同而偏移但趋势是稳定的。SCR主导模态振荡频率 (Hz)主导模态阻尼比稳定性状态10850.42稳定裕度充足5460.31稳定裕度下降3240.12临界接近弱阻尼216-0.04失稳包络发散看到这个趋势后优化目标就很明确了想办法在 SCR2、3 的条件下把主导模态阻尼比抬高同时尽量不牺牲 SCR10 时的动态性能。5. 控制策略优化重塑阻抗、补偿相位与参数协同整定5.1 虚拟电阻和虚拟电感最直接的重塑手段跟网型逆变器失稳的根源是输出阻抗在特定频段呈现负电阻特性那最直接的办法就是人为改变输出阻抗形状。虚拟电阻法在电流环指令或调制波中叠加与输出电流成比例的分量等效于在输出回路串入一个电阻。这个电阻能把负阻尼拉正实现非常简单在 dq 轴调制波上分别减去 R_v × i_d 和 R_v × i_q 即可。虚拟电阻取值从电网阻抗的 0.3 倍开始试观察主导特征根变化。R_v 增大阻尼会明显改善但代价是等效输出电压被“吃掉”一块负载能力下降重载时可能达不到额定电压。所以 R_v 不是越大越好以刚好恢复正阻尼并留 5~10 度相位裕度为目标。虚拟电感类似在调制波中加入与电流微分成比例的项能在频域削掉某个特定频段的阻抗凸起。但虚拟电感在数字控制中涉及微分项对噪声敏感离散化后容易引入高频抖动实际工程里要加低通滤波或观察器处理。5.2 电网电压前馈与相位补偿在Simulink里的实现电网电压前馈是改善抗扰性的常用手段把 PCC 电压检测值经过低通滤波后直接叠加到调制波上用来抵消电网电压扰动对输出电流的影响。前馈通道越完整逆变器对外表现为“电压源跟随”但弱电网下 PCC 电压本身受电流影响大前馈增益太大会把噪声和畸变也引进来反而激发高次谐波。所以前馈通道一定要串一个低通滤波器截止频率一般取电流环带宽的两倍以上。相位补偿则更精细一些。PLL、电流环、滤波器都会在特定频段引入相位滞后如果在调制波中对某个关键频率做超前校正可以抬高相位裕度。Simulink 里实现相位补偿可以用 Continuous 传递函数模块也可以在电流内环误差信号后面串一个超前-滞后环节。做完补偿后我会回到阻抗法画一次 Bode 图看 PLL 频段附近的相位曲线是否整体上移。5.3 基于灵敏度指标的参数协同整定流程单打独斗的参数优化不够系统我在课题里跑了一套参数协同整定流程先定义目标函数比如在 SCR2、3、5 条件下所有主导特征根阻尼比不小于 0.1同时最大特征根实部留 20% 以上裕度设计变量取 PLL 带宽、PLL 阻尼比、电流环带宽、虚拟电阻、虚拟电感。然后先用解析小信号模型做批量网格扫描或粒子群优化得到候选参数再放回 Simulink 平均值模型做线性化校验最后挑一组鲁棒性最好的参数做时域验证。这个流程的核心优势是把“控制参数整定”问题转化成了“特征值约束优化”问题不会出现调了这个参数丢了那个模态的情况。粒子群优化在 MATLAB 中可以直接调 Global Optimization Toolbox目标函数里每次调用小信号状态矩阵计算特征值几分钟就能跑完一轮。优化措施实现位置主要效果副作用 / 注意点降低 PLL 带宽PLL 调节器参数减弱弱电网正反馈增益动态响应变慢需和启动需求平衡虚拟电阻调制波叠加电流比例项消除负阻尼区间等效压降增加重载能力下降电网电压前馈 低通调制波前馈通道提高抗扰性前馈过强会引入高频噪声相位超前补偿电流内环或调制波环节抬高相位裕度高频段易放大噪声需配合滤波6. 实践中容易翻车的地方假失稳、线性化误差与模型陷阱6.1 三个容易误判为失稳的实际案例第一个案例是启动瞬间的“假振荡”。PLL 初始相位和三相电压实际相位不一致时启动瞬间冲击电流很大dq 轴电流波形会出现一个非常大的衰减振荡初学者容易把它当成失稳。解决方法是设置 PLL 初始相位为 0 时同步匹配三相电网电压初始相角或者在仿真启动阶段先开环预同步等 PLL 锁定后再切入闭环。第二个案例是数值步长导致的“假发散”。开关级模型在仿真步长过大时PWM 高频分量混叠到低频波形上看就像低频振荡发散。把步长从 1e-5 改成 1e-6 后波形立刻干净了和平均值模型结果也能对上。遇到可疑振荡先检查步长和求解器别急着改控制参数。第三个案例是 dq 变换参考轴约定不一致导致的“特征值和波形对不上”。Simulink 中 Park 变换有很多种约定有的用 sin 型、有的用 cos 型有的 q 轴超前 d 轴 90 度有的滞后。控制环里的 Park、反 Park、PLL 输出角度必须遵循同一套约定。我在排错时最常用的检查手段是给模型加一个已知幅值和相位的三相正序电压看 dq 轴的输出是不是符合预期这一步 10 分钟就能定位约定错误。6.2 线性化模型与开关级仿真结果为什么会差线性化模型忽略了三类因素PWM 开关纹波、控制延时、限幅非线性。其中最容易被忽略的是数字控制延时。实际控制器从采样、计算到 PWM 更新至少有一拍到两拍延时这个延时在几十赫兹到几百赫兹频段会带来明显的相位滞后。如果线性化模型里不加入等效延时环节得到的稳定边界会比实际系统乐观也就是算出来稳定时域开关模型已经振荡了。我在状态空间模型中一般用 e^{-1.5Ts·s} 做控制延时的近似Ts 是控制周期典型值 100us。在 Simulink 平均值模型里可以用 Transport Delay 模块模拟这个延时。这样线性化结果和开关模型的一致性会好很多。限幅非线性在小扰动分析里可以忽略但大扰动验证时一定要保留因为电流限幅会让系统工作点大幅偏移线性化模型完全覆盖不了这个现象。6.3 仿真实验设计的几条建议参数扫描建议至少覆盖 SCR10、5、3、2 四档每档都记录主导特征根、阻尼比和时域扰动测试结果。保存波形时既要保存 d 轴电流、q 轴电流、PCC 电压幅值也要保存 PLL 输出角频率因为失稳早期往往先在 PLL 频率上出现振荡迹象。还要注意工作点漂移问题。电网阻抗变化不仅改变动态特性也改变了静态工作点功率角、PCC 电压幅值都变了。线性化必须在新工作点附近做不能拿 SCR10 的工作点线性化结果去判断 SCR2 的稳定性。Simulink 线性化工具支持先对模型做 trim 找到平衡点再做线性化这个流程一定要走完整。我个人的一个小经验是把参数扫描、波形提取、特征值计算全部封装成 MATLAB 脚本用结构体统一管理参数跑 10 组参数用不了一分钟。这样你可以大胆地试各种 PLL 带宽、虚拟阻抗组合不会因为手动改参数改到怀疑人生。课题做到后期这套脚本就是最值钱的资产。
