简化时域模型与前馈控制：LLC谐振变换器动态响应优化

简介这是一份面向电力电子工程师与研究者的PDF技术文档围绕LLC谐振变换器传统线性控制动态响应不佳的问题提出基于简化时域分析并融合频率前馈的优化控制方法。作者按谐振点附近、高于谐振点且远离、高于谐振点但接近三种工作状态分别建立简化时域迭代算法可在数字控制器中实时计算理想开关频率并前馈至线性控制器同时输出电压增益、谐振电流与谐振电压波形实验验证显著改善了动态响应。资源包仅1个PDF文件大小777KB内含论文全文、数学推导、算法实现代码示例与解释方便读者直接复现和调试。已有147人学习下载适合从事LLC变换器控制算法开发或电力电子数字控制研究的工程师作为技术参考。 做电力电子的人大概都绕不过LLC。不管是数据中心电源、车载充电机还是储能里的双向变换器LLC几乎成了中大功率隔离DC-DC的默认选择。可真正上手做控制的时候很多人会卡在同一个地方小信号模型太麻烦时域分析又怕算力扛不住最后只能用基波近似法FHA凑合调环动态响应一塌糊涂。我这次想聊聊一套相对“接地气”的做法——简化时域算法Simplified Time-Domain Model配合频率前馈控制用一张查表表把谐振变换器的非线性工作点“录”下来再用前馈把动态响应从“PI慢慢爬”变成“一步到位”。整套思路我已经在DSP上跑过也在PSIM里做过联合验证下面会给出可复现的代码框架和完整的推导思路。这篇文章适合谁如果你已经在用LLC做电源但觉得FHA建模精度不够或者正在为动态响应头疼再或者你想知道“时域算法到底怎么落到实时控制里”那这篇应该能帮你省下不少试错时间。就算你是刚入门LLC的工程师读懂里面的归一化思路和查表法也会对整个拓扑的非线性特性有很大帮助。1. 为什么选时域算法FHA的短板和时域模型的真正价值1.1 基波近似法在什么场景下会“翻车”先说说大家最熟悉的FHA。它的思路很简单把方波激励分解成基波用基波等效电路来分析谐振腔然后算出增益曲线。在谐振频率附近、满载条件下FHA算出来的增益和工作点基本靠谱这也是它成为“入门第一课”的原因。但问题是LLC的实际工况很少一直待在谐振点附近。轻载的时候励磁电感的分压效应非常明显FHA的等效负载电阻已经不准宽范围输出的时候开关频率可能跑到1.5倍谐振频率以上高频段的谐波分量占比上升基波近似误差会迅速拉大。我见过不少项目FHA算出来明明能软开关实测却丢失了ZVS根源就在建模精度不够。更关键的一点FHA只能告诉你“某个频率下增益大概是多少”它给不出谐振电流的瞬时波形、死区时间内的励磁电流方向、以及状态变量在每个开关周期的具体轨迹。而这些东西恰恰是设计死区、设计同步整流时序、判断ZVS边界所必须的。1.2 时域算法在做一件什么事所谓简化时域算法本质上是把LLC的开关周期按照工作模态切分成若干段。每个模态下谐振腔等效成一个线性二阶电路状态变量谐振电感电流、谐振电容电压、励磁电流可以用一组解析表达式描述。然后通过模态间的边界条件递推迭代到稳态就能得到精确到开关周期级别的波形信息。相比全数值仿真比如Simulink里搭非线性磁件模型跑瞬态时域算法的优势在于计算量可控。不需要跑到微秒级的小步长而是一次开关周期内解几十次代数方程即可收敛。这给了它两条路一条是离线扫频建表把“负载-输入电压-频率-增益”的关系提前算好另一条是在线轻量迭代在控制周期里直接算出当前工况下的目标频率和软开关边界。对我个人来说时域算法更大的价值在控制侧它能提供一个“频率前馈”所需的全工况映射。前馈控制的核心前提是你得知道“输入电压变了多少频率应该跟着变多少”的定量关系这个问题FHA给不准但时域表可以给得很细。2. 简化时域算法的数学建模与归一化推导2.1 工作阶段划分从半个周期来看LLC先限定一个常见场景全桥LLCPFM调频工作在断续模式DCM和连续模式CCM的交界附近这也是中大功率电源最常用的区间。简化时域分析最常用的是“半周期对称法”——因为LLC在正负半周的工作过程完全对称我们只用分析从开关管Q1关断到Q2关断这半个周期另一半周期做镜像处理即可。在这半个周期里典型的工作阶段可以划分为三个励磁电感参与谐振阶段Lr与Lm串联谐振变压器原边电流不再向副边传递负载由输出电容供电谐振电感与谐振电容谐振阶段Lm被输出电压钳位只有Lr和Cr参与谐振能量从原边传递到副边死区阶段所有开关管关断励磁电流给结电容充放电实现ZVS。每个阶段都是一个线性时不变电路可以写出解析解。以第2阶段为例谐振槽路的微分方程为[ L_r \frac{di_{Lr}}{dt} v_{Cr} v_p v_{in} ]其中 (v_p) 是变压器原边电压在电流传递阶段等于 (nV_o)n为变压器原副边匝比(V_o) 为输出电压。把状态方程写成矩阵形式后解出来就是标准的二阶振荡解[ i_{Lr}(t) \frac{V_{in} - nV_o}{Z_r}\sin(\omega_r t) I_{L0}\cos(\omega_r t) ][ v_{Cr}(t) V_{in} - nV_o - (V_{in} - nV_o - V_{Cr0})\cos(\omega_r t) I_{L0}Z_r\sin(\omega_r t) ]其中[ \omega_r \frac{1}{\sqrt{L_r C_r}}, \quad Z_r \sqrt{\frac{L_r}{C_r}} ](I_{L0}) 和 (V_{Cr0}) 是当前阶段的起始电流和起始电容电压。看起来简单但真正麻烦的是如何把三个阶段的边界条件串联起来迭代求稳态。2.2 归一化处理把麻烦的参数压缩到最小集直接对着实际参数亨、法拉、伏特、安培去迭代当然也能算但效率太低而且不利于通用化。我建议先把所有量纲做归一化让算法只和几个无量纲数打交道。定义归一化基准[ V_{base} V_{in}, \quad Z_{base} Z_r, \quad I_{base} \frac{V_{base}}{Z_{base}}, \quad \omega_{base} \omega_r ]于是得到归一化状态变量[ i_n \frac{i_{Lr} \cdot Z_r}{V_{in}}, \quad v_{cn} \frac{v_{Cr}}{V_{in}} ]还有两个极其关键的归一化系数增益 (M \frac{n V_o}{V_{in}})励磁电感比 (k \frac{L_m}{L_r})这两个参数一出来整个LLC的稳态特性基本就被锁定了。M决定工作点是高于还是低于谐振频率k决定励磁电流大小也就是决定ZVS的能力和环流损耗。用归一化之后的方程迭代你会发现所有曲线都能塌缩到少数几张无量纲图表上这比每次改参数都重跑一遍仿真高效得多。归一化后谐振电流的表达式可以重新写成[ i_n(\theta) (1 - M)\sin\theta I_{n0}\cos\theta ][ v_{cn}(\theta) 1 - M - (1 - M - V_{cn0})\cos\theta I_{n0}\sin\theta ]这里的 (\theta \omega_r t) 成为自变量也就是我们常说的“角度域”分析。用角度域有一个好处在PFM调频下开关周期的变化体现在总角度 (\theta_{sw} 2\pi f_r / f_s) 上不同负载条件只改变各阶段占用的角度长度而不再需要反复处理时间尺度上的变化。2.3 从边界条件到稳态解一遍迭代的完整过程稳态求解的思路是用半个周期的末态作为下半个周期的初态不断迭代直到状态变量不再变化。具体步骤如下给定输入参数 (V_{in}, V_o, P_o, f_s)算出 M、等效负载电阻、励磁电感电流幅值假定一个初态通常取 (i_{Lr}0, v_{Cr}0)或者上一次的末态按时间顺序依次求解三个阶段的解析表达式每个阶段算完后把末态更新为下一阶段的初态半个周期结束后利用对称性构造另半个周期的初始条件(i_{Lr}(T_s/2) -i_{Lr}(0), v_{Cr}(T_s/2) -v_{Cr}(0))检查当前末态和由对称性推导出的状态是否一致如果误差大于容忍度把末态作为新初态回到步骤3继续迭代迭代收敛后记录这个工况下的稳态增益、ZVS边界信息、各阶段占用的角度。实际编写这段代码时迭代通常15~30次就能收敛比Simulink跑一段上百微秒的瞬态仿真快几个量级。我在一台普通PC上用MATLAB做全工况扫描输入电压20个点、负载30个点、频率100个点大约几分钟就能生成完整的时域稳态表。这就是“简化”二字的来源不去追求每个瞬态毛刺的精确还原而是聚焦于稳态工作点的快速求解。对控制工程师来说稳态工作点的精度远比瞬态波形细节重要。3. 动态响应为什么要靠前馈环路增益的固有矛盾3.1 LLC小信号模型的物理难点很多人在LLC上直接用电压环PI动态响应调到某一带宽后就再也推不上去了。表面上看是补偿器参数没调好根子上其实在小信号模型。LLC是一个高阶非线性系统谐振腔里有两个储能元件再加上励磁电感和输出电容在频率域里至少是四阶系统。它的控制到输出的传递函数 (G_{vd}(s) \hat{v_o}/\hat{f_s}) 在不同负载条件下变化极大轻载时增益曲线平滑重载时在谐振频率附近存在一个“山谷”相位变化剧烈。用一个固定参数的PI去覆盖整个工作区间只能折中不可能让每个工况都达到最优。另一个麻烦是LLC的开关频率本身就很高通常100kHz到1MHz数字控制的采样频率虽然能到几十kHz但相位延迟依然明显。如果完全依赖反馈环去纠正输入电压突变带来的输出跌落至少要经过一个“采样—PI计算—PWM更新”的延迟周期动态响应天然受限。3.2 频率前馈控制把扰动“计算”掉而不是“等待”掉前馈控制的思想很朴素既然我们知道扰动源输入电压、负载电流也知道稳态下频率和输出的映射关系那就可以在扰动出现的瞬间直接把频率推向新工作点附近剩下的细微误差再由PI赔偿。关键问题是这个“映射关系”从哪来FHA给的关系在宽范围内不够准而时域算法恰恰能提供一组精确的静态映射表或拟合函数。用查表的方式做前馈等于把一个“经过精确校准的非线性补偿器”并联在反馈环上PI只需要处理模型误差和未建模的动态带宽压力大幅下降。我实测下来的效果是在突加50%负载的工况下输出电压的跌落幅度能从300mV左右压到100mV以内恢复时间也能缩短到原来的三分之一左右。这是单纯的“把PI调激进”做不到的因为带宽一旦推高噪声和相位裕度问题马上会反噬。3.3 前馈量的数学表达时域模型给出的映射可以写成在某一输入电压和负载条件下为了维持输出电压恒定需要的开关频率为[ f_s F_{map}(V_{in}, P_o) ]前馈量可以定义为[ \Delta f_{ff} \left.\frac{\partial f_s}{\partial V_{in}}\right|{P_o, V_o} \cdot \Delta V{in} \left.\frac{\partial f_s}{\partial P_o}\right|{V{in}, V_o} \cdot \Delta P_o ]其中两个偏导项就是从时域表中数值差分得到的增益系数。实际写代码时我会把这两个系数本身也做成二维表或者用多项式拟合这样在线计算只需要几次乘加运算占用资源极小。数字实现上有两个地方要注意前馈和反馈必须做解耦。前馈输出的是“参考频率的粗调值”反馈输出的是“参考频率的精修值”两者相加后送给VCO或PWM模块。前馈增益不能过大否则扰动方向上“过冲”会和反馈环形成振荡。我一般会给前馈通道加一个一阶低通滤波截止频率放在开关频率的1/10以下既能快速响应负载突变又能滤除开关纹波对采样值的污染。4. 实时控制实现代码架构与核心模块4.1 控制系统的整体框架整个实时控制程序可以分成三层底层ADC采样、PWM比较器更新、故障保护过流、过压、过温中间层频率前馈计算、时域状态查表、软启动/掉电检测等模式管理上层电压环PI调节器、环路状态监控、通信接口。核心思路是每个开关周期或者每两个开关周期做一次控制计算ADC采样和PWM更新严格同步在载波的峰谷点避免采样点落在开关噪声上。我实现的架构示意如下文字版便于你在自己的工程里复刻/* LLC控制结构体 */ typedef struct { float fs; /* 当前开关频率Hz */ float fs_ref; /* 频率参考值前馈 反馈合成 */ float vo_ref; /* 输出电压参考值V */ float vin_meas; /* 输入电压采样V */ float po_est; /* 输出功率估算W */ float kff_vin; /* 输入电压前馈增益Hz/V */ float kff_load; /* 负载前馈增益Hz/W */ float f_min; /* 最低开关频率限制ZVS边界 */ float f_max; /* 最高开关频率限制死区占比 */ } LLC_Control_t;4.2 频率前馈的代码实现前馈计算的代码很简单但有两个细节值得注意一是前馈增益必须根据当前工作点查表而不是用固定常数二是前馈生效的前提是输入电压和负载功率的估算足够及时所以ADC采样的顺序和滤波窗口要设计好。我通常在中断服务函数里这样写float LLC_FeedForward(LLC_Control_t *ctrl) { float df_vin 0.0f; float df_load 0.0f; /* 查表获取当前工作点的前馈增益 */ float k_vin Table_GetKvin(ctrl-vin_meas, ctrl-po_est); float k_load Table_GetKload(ctrl-vin_meas, ctrl-po_est); /* 输入电压前馈电压升高频率应增大 */ df_vin k_vin * (ctrl-vin_meas - VIN_NOMINAL); /* 负载前馈负载增大增益需求增大频率应降低 */ df_load k_load * (ctrl-po_est - PO_NOMINAL); return ctrl-fs df_vin df_load; }这里查表函数返回的是 ( \partial f_s/\partial V_{in} ) 和 ( \partial f_s/\partial P_o )这两个表是在上电初始化时从存储区加载到RAM的。如果MCU资源充足也可以把整个 ( f_s F_{map}(V_{in},P_o) ) 做成二维表前馈时直接查频率值连增益计算都省了。我推荐后一种做法因为它在负载突变时更直接只是初始化时表会稍微大一点比如21×21个float折合约1.7KB绝大多数MCU都放得下。4.3 主控制循环从采样到PWM更新的完整逻辑主控制循环采用经典的“状态机 周期中断”结构。ADC在PWM载波峰值触发的同步信号下启动转换转换完成后进入中断服务函数中断里执行所有控制计算最后把更新好的占空比/频率值写入PWM比较寄存器。核心伪代码如下interrupt void ISR_Control(void) { float vout ADC_GetValue(CH_VOUT); float vin ADC_GetValue(CH_VIN); float iout ADC_GetValue(CH_IOUT); /* 1. 更新功率估算 */ float po vout * iout; /* 2. 频率前馈计算 */ float fs_ff LLC_FeedForward(ctrl); /* 3. 电压环PI调节器 */ float verr ctrl.vo_ref - vout; float pi_out PI_Update(verr); /* PI输出作为频率修正量 */ /* 4. 合成参考频率并做限幅 */ ctrl.fs_ref fs_ff pi_out; if (ctrl.fs_ref ctrl.f_max) ctrl.fs_ref ctrl.f_max; if (ctrl.fs_ref ctrl.f_min) ctrl.fs_ref ctrl.f_min; /* 5. 查时域表获取死区时间和ZVS边界信息 */ Lookup_DeadTime(ctrl.fs_ref, ctrl.po_est, dead_time); /* 6. 更新PWM模块频率、死区、互补相位 */ PWM_UpdateFrequency(ctrl.fs_ref); PWM_UpdateDeadBand(dead_time); }这段代码的精华在第2步和第4步的配合前馈先把频率“预置”到接近目标值的位置PI只在残差上做功。这样一来即使PI带宽调得不太高动态响应也能保持优秀。用MATLAB Simulink做控制参数验证时我会把前馈通道先断开单独调好PI的稳定裕度再接上前馈检查动态指标。4.4 死区时间的动态查表时域算法的“赠品”很多LLC控制实现里死区时间是固定值。但时域模型可以直接输出每个工作点下的励磁电流峰值进而估算出最小所需死区时间。死区不够ZVS可能丢失死区太长又增加有效占空比损失和环流损耗。这个问题在宽范围应用中很容易成为效率杀手。简化时域表里我额外存了一列“安全死区时间”定义基于励磁电流在死区阶段对开关管结电容的充放电时间估算。代码里用一个线性插值查表来读取数据和频率前馈表共用同一个索引结构几乎不增加额外开销。实测下来动态死区比固定死区在中载区间能提升约0.3%~0.5%的效率在轻载区间提升更明显因为轻载时励磁电流较小固定死区往往偏大。5. 离线建表与在线查表的工程落地5.1 离线扫描脚本的设计思路时域表不是手动算出来的而是用一个离线脚本预先扫描生成。我习惯用MATLAB写扫描脚本因为它语法自由矩阵运算方便。核心逻辑就是遍历输入电压、负载功率、开关频率三个维度在每个组合下调用时域稳态迭代函数得到该点的工作状态。关键代码框架%% 参数定义 Lr 18e-6; % 谐振电感 Cr 22e-9; % 谐振电容 Lm 110e-6; % 励磁电感 n 4; % 匝比 fr 1/(2*pi*sqrt(Lr*Cr)); %% 扫描网格 Vin_list 250:25:450; Po_list 50:50:1000; fs_list 0.3*fr:2e3:1.6*fr; %% 逐一扫描并记录 for i 1:length(Vin_list) for j 1:length(Po_list) for k 1:length(fs_list) [M, zvs_ok] LLC_TimeDomainSteady(...); if zvs_ok table(i,j,k) M; else table(i,j,k) NaN; end end end end扫描完以后为了满足实时控制对存储空间的限制我会做两步压缩第一步是只保留满足ZVS条件的区域把不可行点的表项直接置空第二步是降采样并保存为C语言的const数组或者二进制文件在MCU初始化时加载到内存。对21×21×96的表用float存储大约16KB如果空间紧张可以改成uint16存储索引再配合线性插值恢复精度。5.2 在线查表与插值技巧在线查表的核心是二维插值。对输入的 (V_{in}) 和 (P_o)先找到它们在表中最接近的四个邻域点然后做双线性插值。实际测试中双线性插值的误差远小于表分辨率本身足够工程使用。插值函数的核心代码如下float Table_Interp(const float *tab, int xIdx, int yIdx, float fx, float fy) { float v00 tab[ yIdx * XSIZE xIdx]; float v10 tab[ yIdx * XSIZE xIdx 1]; float v01 tab[(yIdx1) * XSIZE xIdx]; float v11 tab[(yIdx1) * XSIZE xIdx 1]; float v0 v00 (v10 - v00) * fx; float v1 v01 (v11 - v01) * fx; return v0 (v1 - v0) * fy; }查表的效率高到几乎可以忽略了一次二维插值只需要四次浮点乘加也就是几十个CPU周期。这也是我敢在控制中断里每次都查表的原因。5.3 表数据校验仿真和实验的两道关时域表建出来之后不能直接拿去用第一次必须做校验。第一道关是仿真校验。我用PSIM搭了LLC的详细开关模型在若干个采样点把仿真的稳态频率和时域表预测的频率做对比。一般误差在1%以内就说明表的结构和迭代代码没有原则性错误。如果误差偏大优先检查Lm被钳位的边界条件是否处理正确这是最常见的建模错误来源。第二道关是实验校验。样机上跑一个输入电压斜坡同时记录输出电压和开关频率反推实际工作点的增益曲线跟表数据对照。这样能发现器件寄生参数如变压器漏感、开关管结电容造成的偏差。通常我会在时域表的基础上加一个“标定修正系数”把实测偏差吸收入系数里而不用推翻整个模型。6. 调试过程中踩过的坑常见问题与排查实录6.1 轻载启动时查表无解频率指令乱跳第一版程序上电后轻载空载启动时PWM频率经常跳变输出电压有时间歇性跌落。查了半天发现是表里轻载区域出现了“无解区”——当负载很轻时为了维持目标增益频率必须拉得非常高甚至超出我设定的频率上限导致查表返回无效值。解决方法是两个一是把表的频率上限放宽到2.5倍谐振频率让算法在轻载时能找到合法的过渡工作点二是在查表函数里加了一个“回退逻辑”如果目标点不可用就沿着负载功率方向逐步回退到最近的可用点。6.2 前馈增益过强突加负载后电压过冲反而不降前馈确实大幅压缩了电压跌落但有一轮测试发现恢复阶段之后出现了明显的过冲大概超过目标电压2%左右。分析下来是前馈增益过大导致频率被推到了低于最终稳态值的位置随后PI再把它拉回来时已经产生了超调。对策是给前馈通道加了一阶低通时间常数设为50微秒左右。同时在代码里限制了频率变化率dFs/dt不让频率在单次控制周期里跳变超过2kHz。这样牺牲了极少量的响应速度换来了更平稳的动态过程整体上是划算的。6.3 高频区ZVS丢失实测效率突然掉了一个点还有一次在高压轻载区间测试效率曲线在某个频率点附近突然掉了将近一个百分点。用示波器看开关管Vds波形发现是ZVS丢失了。查时域表里的励磁电流估算值发现它确实在那一带变得很小已经不足以在死区时间内完成结电容充放电。这个问题只靠动态死区只能缓解不能根治。后来我在控制里加了“ZVS保护频率下限”——一旦查表发现励磁电流校核不过就直接把频率钳制到一个更低的值不进入ZVS失效区。虽然会牺牲一些轻载电压调整率但保证了硬开关不出现对可靠性的收益远大于效率微损。6.4 ADC采样点和开关噪声打架导致PI输出抖动硬件调试时还发现一个特隐蔽的问题控制中断里的输出电压采样值有周期性抖动PI输出跟着抖导致开关频率出现微小但明显的音频噪声。查了半天发现是采样触发点没和PWM载波对齐导致采样窗口落入了死区附近的振铃尖峰上。修法是把ADC同步触发信号改为PWM载波的峰值点并在采样后做3次滑动平均。这样可以在不降低环路带宽的前提下把采样噪声压到几乎可以忽略的程度。6.5 常见问题速查表问题现象可能原因排查手段与对策启动时频率乱跳查表区域不连续或无解放宽频率表范围增加回退逻辑突加负载过冲明显前馈增益过大前馈通道加低通限制频率变化率高频区效率骤降ZVS丢失校核励磁电流估算设置ZVS保护频率下限PI输出抖动采样点落在开关噪声上同步ADC触发滑动平均采样表数据与实测偏差大器件寄生参数影响实测标定修正系数更新表数据7. 一点实操体会把这套简化时域算法和频率前馈从“论文里的公式”搬到“样机上的代码”整个过程比我预想的要曲折但最后的收益非常直接。时域算法让我第一次在LLC上不再“对着空气调参”每一个频率指令背后都有明确的物理映射调试时用示波器抓一次波形就能定位问题。如果你准备在自己的项目里复现这套思路我建议先别急着写完整代码。先花一天时间把时域迭代的MATLAB脚本跑通生成一张增益表对照FHA画出的曲线看差异再用这张表在仿真里验证前馈效果最后才移植到MCU。很多人图快直接写固件结果表数据的边界条件没搞清越调越乱。后续如果你想在这套基础上做三相LLC或者CLLC双向变换器时域建模的思路大体一致只是每个阶段的状态方程会多一些耦合项查表维度也可能从二维变三维。但只要把“归一化—分段解析—迭代稳态—离线建表—在线查表”这个方法论吃透扩展过去不会太折腾。反正在我这里用过时域表之后我是真不愿意再退回纯FHA去调LLC了。本文还有配套的精品资源点击获取