ARTICLE DETAIL

资讯详情

深耕编程入门与网站建设的一线实战洞察。

基于SLSPC的高阶PT-WPT系统Simulink仿真复现指南

基于SLSPC的高阶PT-WPT系统Simulink仿真复现指南 你是不是也有这种经历导师甩给你一篇SCI一区论文说“把它仿真复现一下”你满怀信心打开PDF两页公式看完就懵了。今天说的这个项目——基于SLSPC系列的高阶PT-WPT无线电能传输系统就是我曾经被“甩”过的一篇。前后折腾了将近三周把Simulink模型从零搭到能跑把论文里的补偿网络一步一步翻译成电路模块最后总算把波形和效率曲线给对上了。这篇文章专门写给正在做无线电能传输仿真复现的人尤其是刚接触PT对称、高阶补偿这些概念的同学。我会把SLSPC拓扑怎么理解、参数怎么算、Simulink里具体怎么搭、中途踩了哪些坑原原本本讲清楚。1. 项目背景与论文思路拆解1.1 这个课题到底在做什么很多人第一次看到“基于SLSPC系列的高阶PT-WPT无线电能传输系统”这个标题第一反应是查缩写。PT-WPT在近几年的无线电能传输文献里通常指Parity-Time对称无线电能传输也就是宇称-时间对称WPT。这个思想最早来自量子力学里的PT对称理论后来被引入到电磁学和电路领域核心是让系统的增益和损耗刚好达到平衡使系统在小扰动下表现出非常稳定的输出特性。对WPT来说这意味着什么意味着负载变化、线圈偏移、距离改变时输出功率不用靠复杂的控制算法去补偿也能基本保持不变。这个特性对无线电能传输来说太有吸引力了。常规的磁耦合谐振WPT系统只要负载一变化或者发射线圈和接收线圈稍微错位输出功率和效率就会剧烈波动必须额外加阻抗匹配网络、频率跟踪算法或者DC-DC后端来稳压。而PT对称结构相当于把“自适应”做进了电路本身在特定区间内天然具备恒功率输出能力。这篇论文的贡献就在于把PT对称机制和高阶补偿网络结合起来用SLSPC系列拓扑去构造这种有源-无源平衡关系从而在较宽的耦合系数和负载范围内保持传输特性平稳。再说SLSPC。这个缩写原论文里给出的是一个拓扑族逐字母拆开看S是串联SeriesL是电感支路InductorSPC是串并联补偿Series-Parallel Compensation合在一起可以理解为“原边串联加电感支路、副边串并联混合补偿”的一类高阶网络。不同论文里这个缩写可能指代不同结构复现之前一定要回到原文看拓扑图。我这里讲的是我当时复现时采用的解读整个Simulink模型也是按这个思路建的。1.2 为什么值得花时间复现这篇论文一区论文不一定都复杂但通常把一个关键问题讲透了。这篇论文的价值不在某个花哨的新元件而是把两个相对成熟的概念——高阶补偿网络和PT对称原理——有机地结合到了一起。复现这种论文收获是双重的既要吃透高阶谐振网络的设计逻辑又要把PT对称的实现条件在仿真里落地这比单纯搭一个SS型WPT模型要深入得多。从实操角度看这篇论文的可复现性其实不错。它给出的系统结构明确补偿网络和线圈参数都有据可循Simulink仿真模型作为论文验证手段关键波形和效率曲线也比较清晰。我复现的时候目标是做到两点一是把论文里提出的SLSPC高阶网络搭成可运行的仿真模型二是通过参数扫描验证PT对称区间内的恒功率输出特性。如果你正在做类似课题或者需要复现别的WPT论文这套方法完全可以平移过去。整个系统的信号链是这样的直流电源 → 高频全桥逆变器 → 原边SLSPC补偿网络 → 发射线圈 → 接收线圈 → 副边SLSPC补偿网络 → 整流桥 → 滤波电容 → 负载电阻。和常规WPT系统相比差别全在补偿网络那一块。二阶系统只有一个谐振电容高阶系统则有多级LC谐振单元每个单元对频率、增益、谐波抑制都有各自的贡献。2. SLSPC拓扑与高阶补偿网络设计要点2.1 从二阶补偿到高阶补偿到底多了什么理解SLSPC的前提是先搞清楚为什么二阶补偿不够用。最基础的SS型WPT系统原边和副边各一个补偿电容谐振条件要求工作频率等于线圈自感和补偿电容的谐振频率。它的优点是结构简单、恒流特性好但缺点也很明显当耦合系数变化或者负载变化时系统的传输功率和效率会明显波动。这是因为二阶谐振网络只有一个谐振点对外部参数变化几乎没有“缓冲”能力。高阶补偿网络加进去的本质上是额外的谐振零点和极点。以LCC补偿为例它在原边线圈前面又插入了一个电感支路和电容支路这些支路和线圈自身电感、补偿电容共同作用会在目标频率附近形成更陡峭的阻抗特性同时让系统的等效输入阻抗在一定范围内对负载变化不敏感。换句话说高阶网络是用更多的电抗元件换取系统对参数漂移的容忍度。把二阶SS、常规LCC-LCC和SLSPC放在一起对比区别会更直观特性二阶SS补偿LCC-LCC补偿SLSPC系列谐振网络阶数2阶4阶以上高阶多谐振点对负载变化的敏感度高中低接近PT平衡时几乎恒定对耦合系数变化的容忍度低中高恒流/恒压输出能力单一特性可通过参数调节结合串并联结构灵活设计仿真复杂度低中较高谐波抑制能力弱中强2.2 SLSPC网络拓扑是怎么推算出来的复现的时候我先把论文给的拓扑图翻译成电路级描述。常见的二阶SS网络原边就是电感Lp串联电容Cp副边是Ls串联Cs。SLSPC系列在此基础上做扩展原边保留串联谐振支路并在输入端加入一个谐振电感形成一个L型输入匹配结构副边不再只是单一串联电容而是拆成一条串联谐振支路和一条并联补偿支路两条支路组合后再通过一个电感到整流桥。用语言描述比较抽象我拆成三个特征来记原边输入端有一个附加电感Lf和原边补偿电容一起构成输入谐振网络它的作用是让逆变器的输出电流接近正弦并抑制高频谐波。原边线圈和副边线圈各自保留基础串联谐振电容保证能量能在线圈之间有效耦合。副边除了串联谐振电容还并联了一个补偿电容支路用来调节整流桥入口的等效阻抗使系统整体更容易落在PT对称的工作区间。这种结构的好处是多了两个可调参数维度原边和副边的阻抗特性不再锁死在单一谐振频率上。我在仿真里验证过只要Lf和副边并联电容取值合适系统的输出功率随负载变化会明显变平这正是PT对称特性的体现。2.3 关键参数计算与匹配思路参数设计是复现中最容易翻车的一步但逻辑其实很固定。先把工作频率定下来WPT仿真里常用的是50kHz到150kHz这一档我用的100kHz。然后根据线圈自感算基础谐振电容。以原边线圈Lp50μH为例谐振电容计算公式是C 1 / (ω²L)其中ω 2πf。在100kHz下ω约等于628318.5 rad/s代入后得到C约等于50.7nF。副边线圈如果自感相同补偿电容也取50.7nF。附加谐振电感Lf我取20μH对应补偿电容Cf也按同一公式算大约126.7nF。耦合系数k取0.3互感M k√(Lp×Ls)也就是0.3×50μH15μH。负载电阻先按20Ω设置。我把整套初始参数整理成了一张表仿真搭建时直接往参数脚本里填就行参数符号数值输入直流电压Uin100 V开关频率fs100 kHz原边线圈自感Lp50 μH副边线圈自感Ls50 μH耦合系数k0.3互感M15 μH原边补偿电容Cp50.7 nF副边补偿电容Cs50.7 nF附加谐振电感Lf20 μH附加谐振电容Cf126.7 nF负载电阻RL20 Ω这里要特别说一句这套参数的核心验证逻辑是让原边和副边的基础谐振频率都落在100kHz而附加谐振网络的谐振点稍微偏移让系统在工作频率附近形成两个相邻的谐振模式。PT对称系统恰好需要这种接近简并的模式结构。如果所有支路都谐振在同一个频率PT特性反而出不来。3. Simulink仿真建模全流程实操3.1 环境准备与模块选型我用的是MATLAB R2021b以上版本Simscape Electrical库里的Specialized Power Systems子库老版本叫SimPowerSystems界面略有不同但元件用法基本一致。搭建WPT功率回路核心用到的模块就四类电源和测量、开关器件、无源元件、耦合线圈。直流电压源用DC Voltage Source。逆变器我建议直接用Universal Bridge模块设成两桥臂2 arms开关器件选MOSFET。这样省去自己搭四个开关管的麻烦而且内部自带缓冲电路仿真不容易发散。无源元件用Series RLC Branch把电阻值设成0或者一个很小的寄生电阻设成电感或电容模式。耦合线圈用Mutual Inductance双绕组模式直接填Lp、Ls、M三个参数。所有参数不要手填在模块对话框里改成用MATLAB脚本统一初始化模型变量名引用脚本变量。这样后面做参数扫描不用反复改模块参数非常省事。初始化脚本内容如下% 系统参数初始化 fs 100e3; % 开关频率 100 kHz w 2*pi*fs; % 角频率 Lp 50e-6; % 原边线圈自感 Ls 50e-6; % 副边线圈自感 k 0.3; % 耦合系数 M k*sqrt(Lp*Ls); % 互感 Cp 1/(w^2*Lp); % 原边补偿电容 Cs 1/(w^2*Ls); % 副边补偿电容 Lf 20e-6; % 附加谐振电感 Cf 1/(w^2*Lf); % 附加谐振电容 RL 20; % 负载电阻 Uin 100; % 输入直流电压3.2 功率主回路搭建步骤搭建顺序我建议从后往前先搭负载和整流再搭副边补偿然后是线圈接着是原边补偿最后才是逆变器和电源。这种从后往前的顺序有个好处每一级模块插上后都能有一个明确的测试点方便定位问题。整流部分比较简单用四个二极管搭全桥整流后面加一个滤波电容和负载电阻。滤波电容一般取10μF到100μF之间太小了输出电压纹波大太大了启动电流冲击明显。我测试下来20μF在这个功率等级下比较合适纹波和启动速度都能接受。副边补偿网络按之前说的结构连接收线圈的一端接串联谐振电容Cs然后接整流桥接收线圈的另一端接一个并联补偿电容支路再汇合到整流桥的另一端。之后在整流桥入口处串入附加谐振电感Lf。这个电感的位置很关键它同时承担谐振和滤波两个功能不能省略。线圈部分用Mutual Inductance模块原边绕组和副边绕组没有电气连接各自引出两个端子。设置的时候自感分别填Lp和Ls互感填M注意同名端符号。同名端接反会导致耦合振荡相位完全不对输出功率会异常低落这是非常容易踩的坑。原边补偿网络稍微简单一些一端接发射线圈另一端接逆变器输出的一个桥臂中点。原边补偿电容Cp与电感Lp串联附加输入端再加一个Lf就行。需要说明的是这里Lf和Cp的节点位置如果接错输入阻抗特性会完全改变判断标准是逆变器输出电流应该接近正弦如果看到明显的尖峰电流多半就是网络形态不对。逆变器用Universal Bridge直流侧接100V电压源交流侧两个输出端分别接原边补偿网络的输入端和电源地。3.3 PWM驱动与死区配置WPT系统里的逆变器工作在固定频率方波模式不需要SPWM或SVPWM用50%占空比方波驱动全桥即可。关键是两个桥臂的驱动信号要相差180度对角桥臂同时导通。最简单的方法是用两个Pulse Generator一个初始相位0度另一个初始相位180度频率都设成fs。脉冲宽度设为50%。如果用的是Universal Bridge四个开关管的门极信号在模型里要通过一个Demux拆开或者使用内部的PWM Generator模块。死区时间不能省。虽然纯仿真里开关管是理想器件但如果你后面要接硬件没有死区习惯的话很容易出问题。在Universal Bridge里可以直接设死区时间我通常设0.5μs到1μs对应100kHz开关频率死区占开关周期的5%到10%是合理的。有一点要特别注意Pulse Generator是连续时间模块而Simscape Electrical里很多元件用离散求解器两者混用会导致数据采样不匹配。解决方法是把所有PWM信号通过Rate Transition模块统一到同一个采样时间推荐设置成1μs也就是开关周期的十分之一既能分辨占空比和死区又不会让仿真速度慢到难以接受。3.4 求解器配置与波形提取Simscape Electrical的仿真步长配置直接影响模型能不能稳定跑起来。我试过默认的ode45效果很差开关瞬间的电流尖峰会让求解器频繁报错。后来改成ode23tb问题少了很多。这个求解器对刚性系统有更好的适应性电力电子电路正属于典型的刚性系统。具体设置上相对误差1e-4最大步长限制为开关周期的百分之一也就是1e-7秒。如果不限制最大步长仿真会在开关切换时刻跳过关键细节波形看起来平滑实际上全是错的。仿真时间先跑2ms大概200个开关周期足够系统进入稳态。波形测量用Current Measurement和Voltage Measurement模块分别接在需要观测的支路上。输出功率不要用瞬时功率直接看那样波形毛刺很大建议把电压电流信号送到To Workspace然后在MATLAB里用rms函数计算有效值再用有效值乘积计算输出功率。效率计算则是在输入端加一个功率测量模块把直流侧输入功率和交流侧输出功率相除。最后处理的时候我习惯把节点电压、原边电流、副边电流、输出电压这四路信号存下来然后画在一张图里和论文里的波形图对比。论文里的波形通常是在稳态下截取的所以仿真数据也要取后半段的稳态区间不要包含启动暂态。4. 踩坑记录与问题排查实录4.1 波形毛刺和系统发散仿真最容易遇到的现象是波形前面一段剧烈震荡然后求解器报错停止或者波形布满密密麻麻的毛刺。这种情况九成是步长问题。高频开关电路的最大步长必须足够小我之前用1μs最大步长跑100kHz系统波形还能看但一提高到200kHz就崩了。后来设置最大步长不超过0.1μs稳定多了。另一种常见原因是电感和电容的初始状态。Simscape里的电感默认初始电流是0这没问题但电容要特别注意如果整流后的滤波电容初始电压是0启动瞬间会有一个很大的充电电流冲击可能导致求解器发散。可以在电容模块里设置初始电压也可以接受这个启动暂态只取稳态数据。我的做法是模型里加一个很小的预充电路电阻启动后1ms断开仿真更稳定。4.2 输出功率远低于预期值的排查顺序如果系统跑通了但输出功率比论文低很多不要急着调参数按这个顺序排查先看逆变器输出端波形确认输出电压是±100V方波没有异常压降。再看原边电流正常应该是接近正弦的交流波形如果电流很小说明原边补偿网络没有谐振起来。接着看副边电流和整流后电压如果所有这些都正常但输出功率低问题多半在耦合线圈或负载参数。互感同名端接反是很多人的盲区。Mutual Inductance模块如果同名端接反两个线圈的感应电压方向相反耦合效果几乎为0输出功率会掉到正常值的十分之一以下。判断方法是看原边电流和副边电流的相位关系正常应该接近同相或反相如果相差90度以上多半同名端不对。4.3 谐振频率和驱动频率不匹配的问题很多人觉得“谐振频率100kHz开关频率100kHz两个一样不就行了”。实际不是这么简单。SLSPC高阶网络里有多个谐振支路它们的谐振频率并不完全一样而是分布在100kHz左右形成一个工作带。如果所有支路都严格谐振在100kHz系统可能会表现出非常陡峭的频率选择性稍微偏离就掉功率。反之如果各支路谐振点稍微错开系统在100kHz附近的响应会更平缓PT对称特性更容易出现。参数设置里Lf和Cf是调节这个“频带宽度”的关键。Lf取20μH、Cf取126.7nF时谐振频率其实还是100kHz但由于和线圈支路参数不同整个网络形成两个相邻模式实际工作点处于两个模式之间的“最佳平衡点”。这个平衡点需要通过参数扫描来找到不是一眼能看出来的。4.4 与论文数据对比时的注意事项复现论文最忌讳的是拿着仿真结果追求和论文图一模一样这不现实也没有必要。论文里的参数可能经过归一化处理波形可能有测量误差。正确做法是看趋势一致性论文说负载从5Ω变到50Ω时输出功率恒定你的仿真曲线也应该是基本水平的论文说某个频率点效率最高你的效率曲线峰值位置应该能对上具体数值差两个百分点问题不大。如果趋势都不对那问题就大了。这时候检查论文的补偿网络结构是否和你搭建的一致特别是原边和副边的串并联关系。有的论文把副边分成两个支路其中一个支路经整流桥后接负载另一个支路是空闲谐振器看漏一个支路整个模型就废了。常见现象可能原因处理方法波形毛刺多求解器报错最大步长过大设置最大步长为开关周期的1/100启动冲击电流过大滤波电容初始电压为0允许启动暂态只取稳态数据输出功率异常低互感同名端接反检查Mutual Inductance同名端输出功率随负载变化剧烈高阶网络参数失配查看Lf、Cf参数重新扫描匹配输出电压纹波很大滤波电容偏小增大滤波电容PWM驱动信号不完整采样时间不匹配用Rate Transition统一采样时间5. 扩展思路与实操心得5.1 用参数扫描验证PT对称特性复现论文的高阶PT-WPT系统仿真模型跑通只是第一步真正验证PT特性需要做参数扫描。最直观的实验是扫描负载电阻从5Ω逐步增加到50Ω记录每个负载下的输出功率。普通二阶SS系统在这个范围内输出功率曲线会明显下坠而SLSPC系统如果PT对称条件满足输出功率应该呈现出一个“平台”波动控制在10%以内。扫描耦合系数也有意义。固定负载20Ω把k从0.1扫到0.5观察输出功率。PT对称区间内耦合系数在一定范围内变化时输出功率基本不变超出区间后功率迅速下降形成一个明显的工作窗口。这个窗口的宽度直接反映了系统抗偏移能力也是论文里最常见的对比曲线。5.2 从Simulink仿真走向实验验证仿真模型的一个重要用途是给硬件实验指明方向。比如Lf和Cf的取值可以直接对应到实际电感电容选型逆变器的开关频率和死区时间可以直接照搬到驱动板。我自己做硬件时仿真里的稳态电流有效值用来选功率器件和导线线径电容电压峰值用来确定器件耐压等级。但仿真和硬件之间有明显的鸿沟。仿真里可以忽略线圈的交流电阻和电容的ESR硬件里这些寄生参数对效率有显著影响。我做实验时线圈用利兹线绕制工作频率100kHz下的Q值明显高于普通漆包线。电容也尽量选C0G或薄膜电容避免陶瓷电容的直流偏压特性导致容值偏移。5.3 复现论文的经验与通用方法这套复现过程走下来我最大的体会是高阶WPT系统的仿真搭建核心不在于会连电路图而在于知道每一级网络在干什么。先搭二阶SS系统跑通看波形然后加入附加电感Lf看波形有什么变化再加入副边并联补偿电容继续观察。每加一级系统行为都会变只有每一级的变化都在预期之内最后整个SLSPC网络才能稳定工作。我在实际复现中还养成了一个习惯每次修改参数前先把当前模型的波形截图保存写上参数和现象做成一份简单的实验记录。SLSPC这种多参数系统调参空间很大靠记忆很容易乱。有了记录回头排查问题时一眼就能看出是哪次改动引入了问题。后面再复现类似的LCC-LCC、LCL-LCL或者其他高阶WPT论文时这套方法可以直接平移。万变不离其宗先把基础谐振单元理解透再谈高阶拓扑的组合你也能在一个星期内把一篇一区论文的模型搭得明明白白。
返回列表