• 2
    回复
  • 收藏
  • 点赞
  • 分享
  • 发新帖

LLC matlab仿真

一全桥LLC闭环仿真,仅供参考。LLC matlab

     针对全桥LLC拓扑,利用Matlab软件搭建模型,分别对轻载,满载进行仿真实验,对比使用增益曲线,以及matlab siulink搭建LLC模型的相关工作频率和输出电压关系,仅供参考。

案例1,输入Vin=275V,输出电压Vo=20V,谐振电感Lr=20uH,谐振电容Cr=88nF,励磁电感Lm=66uH,变压器匝比n=13,额定功率P=2kW。

一、LLC增益曲线

全桥LLC增益表达式(s域):

相关谐振电路参数定义如下:,特征阻抗,品质因数。其中等效负载。为实际负载。

画出LLC不同负载下增益曲线图:

图1 LLC不同负载下增益曲线图

图2 LLC不同负载下ZVS以及ZCS划分示意图

图2中可见,在满载情况下,为了保持在ZVS区LLC最低工作频率大约在70.4kHz左右。

二、LLC仿真模型

于MATLAB中搭建全桥LLC仿真模型如下:

图3 全桥LLC仿真模型

       控制部分设计,PI输出的物理意义是LLC的开关频率,最高设置为200kHZ,最低为80kHZ。软起动的时候LLC从低增益开始,实现电压软起动,所以PI的积分初始值设置为200kHZ,本仿真中是带载启动,为了加快软起,所以设置为150kHz。控制以及PWM产生模块如下:

图4 全桥LLC控制部分

2.1 R=5轻载相关仿真波形

1)Vin=275V,R=5轻载谐振腔部分的仿真波形

图5 R=5全桥LLC桥臂电压,谐振电流,以及输出电压

观察上图5可以发现,谐振电流ir滞后于原边桥臂电压VAB,可见谐振腔阻抗在当前开关频率下呈感性。

图6 R=5全桥LLC第一桥臂上管S1的应力,驱动,电流

图7 R=5全桥LLC第二桥臂上管S3的应力,驱动,电流

由图6,图7观察S1,S3的VDS下降为0时,流过S1,S3的电流为负,可见开关管是在VDS下降为0而且对应管子的反并联二极管导通时开通的,因此它们都是零电压开通,而当它们关断时,其结电容的存在,使得它们是零电压关断。

2)R=5轻载输出整流侧相关波形

图8 R=5变压器原副边电压波形

观察图8可以发现,变压器原副边电压符合变压器基本规律,即:

Vsec1,2=Vp/n

图9 R=5变压器原副边相关电流波形

观察图9可以发现,变压器原副边电流符合变压器基本规律,即:

isec1+isec2=n*ip

       同时,在图9也可以发现,只有在谐振电流绝对值ir大于励磁电流im,才会有能量从原边传输到副边(图中画虚线部分)。虚线范围外,对应原边电流ip=0区域的能量将留在谐振腔体之内,这应该是影响效率的关键点之一。

3)根据增益曲线计算输出电压

对比增益曲线计算以及matlab仿真电压,观察matlab模型中在R= 5,最终输出Vo=20V时对应PI工作频率约为128.1Khz。

根据前面增益表达式,我们可以得到对应R=5的增益曲线:

图11 matlab仿真模型R=5对应的增益曲线

已知Vin=275V,n=13,由上图中可知,在频率为128kHZ时,其增益M=0.968

根据:M=(n*Vo)/Vin

则:Vo=M*Vin/n= 20.47V

Matlab中输出电压平均值约为19.994V,PI工作闭环计算的工作频率约为128.1Khz与增益曲线计算结果基本保持一致。

如果输入电压变化范围为:225V至325V,输出20V,M=(n*Vo)/Vin则理论计算需要的增益分别为:1.1556,以及0.8。根据满载的增益曲线图12可知在输入电压为225V时,需要的工作频率为96.8k,在输入电压为325V时,工作频率为178kHz。

图12 matlab仿真模型R=0.2即满载对应的增益曲线

 

2.2 Vin=225V且R=0.2即满载相关仿真波形

1)Vin=225V ,R=0.2满载谐振腔部分的仿真波形

图13 Vin=225V,R=0.2全桥LLC桥臂电压,谐振电流,以及输出电压

由图13可知,此时VAB超前谐振电流,谐振腔在当前工作频率呈感性。

图14 Vin=225V,R=0.2,S1管子应力以及其驱动,电流波形

与R=5类似,如图14所示,S1的VDS下降为0时,流过S1的电流为负,可见开关管是在VDS下降为0而且对应管子的反并联二极管导通时开通的,因此S1是零电压开通,而当S1关断时,其结电容的存在,使得S1是零电压关断,其他管子驱动应力分析同理。

2)Vin=225V ,R=0.2输出侧波形

图15 Vin=225V,R=0.2,变压器原副边电压波形

观察图15同样可以发现,变压器原副边电压符合变压器基本规律,即:

Vsec1,2=Vp/n

图16 Vin=225V,R=0.2,整流侧电流以及原边对应

由图16变压器原副边电流符合变压器基本规律,即:

isec1+isec2=n*ip

在图16也可以发现,只有在谐振电流绝对值ir大于励磁电流im,才会有能量从原边传输到副边,对应原边电流ip=0区域的能量将留在谐振腔体之内,同时由图16可知副边二极管是零电流开关。

在Vin=225V,R=0.2时,控制环路PI输出的工作频率约为99.1kHz,与图12增益曲线对应的96.8k大体一致。

2.2 Vin=325V且R=0.2即满载相关仿真计算

Vin=325V且R=0.2即满载相关仿真(波形略),观察控制环路PI输出的工作频率约为151.6kHz左右。与图12增益曲线得到的178k相差较大。

综合Vin=225V,325V,满载R=0.2,结合图12满载增益曲线图得到下表一:

表一 满载时增益曲线以及matlab仿真数据对比

注:

       可见,matlab仿真闭环的频率与增益曲线所计算的工作频率有所出入,而且在频率高的时候,出入偏得更多,原因是增益曲线求取是基于基波近似的,必然存在误差,而该误差与所在工作频率也存在一定关系。同时matlab仿真模型的元件均为理想元件,所以跟实际项目依旧是必然存在误差的。

 

三、结语

通过matlab,计算LLC增益曲线,以及Matlab搭建模型,方便分析LLC的相关波形,后续可根据实际项目继续补充。

附:llc仿真模型(直接复制,matlabR2016b之后版本均可打开):

 

matlab附增益曲线计算m函数源代码:

 

Lr=20e-6;

Lm=66e-6;

n=13;

Cr=88e-9;

w0=1/sqrt(Lr*Cr);

f0=w0/2/pi

wp=1/sqrt((Lr+Lm)*Cr);

fp=wp/2/pi

 

% R=pi^2*sqrt(Lr/Cr)/8/2

R=20^2/2000;

Re=8*R*n^2/pi^2;

Q=sqrt(Lr/Cr)/Re

Hs=tf([Re*Lm 0 0],[Lr*Lm Re*Lm+Re*Lr Lm/Cr Re/Cr]);

 

% R=pi^2*sqrt(Lr/Cr)/8/1

R=20^2/1000;

Re=8*R*n^2/pi^2;

Q=sqrt(Lr/Cr)/Re

Hs1=tf([Re*Lm 0 0],[Lr*Lm Re*Lm+Re*Lr Lm/Cr Re/Cr]);

 

R=5;

Re=8*R*n^2/pi^2;

Q=sqrt(Lr/Cr)/Re

Hs2=tf([Re*Lm 0 0],[Lr*Lm Re*Lm+Re*Lr Lm/Cr Re/Cr]);

 

P1 = bodeoptions;

P1.FreqUnits = 'Hz';

P1.Grid = 'on';

P1.Xlim={[1e4,1e6]};

 

bodemag(Hs,P1);

% hold on;

% bodemag(Hs1,P1);

% bodemag(Hs2,P1);

legend('满载');

 

 

 

matlab计算llc增益方法二

Lm=66e-6;

Lr=20e-6;

Cr=88e-9;

Np=13;

Ns=1;

Uo=20;

Io=100;

 

K=Lm/Lr;

fr=1/(2*pi*(Lr*Cr)^0.5);

fs=1/(2*pi*((Lr+Lm)*Cr)^0.5);

Rac=Uo/Io*8/pi^2*(Np/Ns)^2;

Q=2*pi*fr*Lr/Rac;

 

Uo=20;

Io=50;

Rac=Uo/Io*8/pi^2*(Np/Ns)^2;

Q1=2*pi*fr*Lr/Rac;

 

Uo=20;

Io=4;

Rac=Uo/Io*8/pi^2*(Np/Ns)^2;

Q2=2*pi*fr*Lr/Rac;

 

X=0.01:0.01:2;

for num=1:1:length(X)

   G(num)=K*X(num)^2/(((1+K)*X(num)^2-1)^2+(Q*K*X(num)*(X(num)^2-1))^2)^0.5;

   

   G1(num)=K*X(num)^2/(((1+K)*X(num)^2-1)^2+(Q1*K*X(num)*(X(num)^2-1))^2)^0.5;

   

   G2(num)=K*X(num)^2/(((1+K)*X(num)^2-1)^2+(Q2*K*X(num)*(X(num)^2-1))^2)^0.5;

   X(num)=X(num)*fr/1000;

end

% [Gmax Gindex]=max(G);

 

axes1=axes();

plot(X,G,'b');

hold on;

plot(X,G1,'r');

plot(X,G2,'y');

grid on;

全部回复(23)
正序查看
倒序查看
zanlu_2
LV.1
2
2020-09-29 15:53
0
回复
兰波
LV.8
3
2020-09-29 22:15
你做的吗?好厉害
0
回复
flymaco
LV.1
4
2020-09-30 11:57
@兰波
你做的吗?好厉害
只是仿真,仅供参考
0
回复
flyergu
LV.1
5
2020-12-28 11:41
看不到?
0
回复
2020-12-29 10:42
看不到仿真楼主大神
0
回复
2020-12-30 23:56
支持!
0
回复
flymaco
LV.1
8
2021-01-07 16:13
@flyergu
看不到?
挨着 在表情后面
0
回复
flymaco
LV.1
9
2021-01-07 16:14
@dy-1atn4AQL
看不到仿真楼主大神
挨着表情后面
0
回复
ygtjcng1
LV.1
10
2021-01-09 14:44
0
回复
dy-cLy8SYla
LV.1
11
2022-10-25 14:49

杠杠滴

0
回复
dy-cLy8SYla
LV.1
12
2022-10-25 15:06
@dy-cLy8SYla
杠杠滴

sin函数作用是啥

0
回复
听听1234
LV.3
13
2022-10-25 18:34

学习学习

0
回复
yingjieso
LV.5
14
2022-10-26 08:13

这水平没得说

0
回复
2023-04-24 23:04

感谢分享,受益匪浅

0
回复
南雁
LV.2
16
2023-10-08 16:09

附件呢

0
回复
南雁
LV.2
17
2023-10-08 16:11

找到了谢谢大神

0
回复
dy-Qdm76yft
LV.4
18
2023-10-08 23:16

48V  35A用半桥还是全桥做好呢

0
回复
dy-Qdm76yft
LV.4
19
2023-10-08 23:16

48V  35A用半桥还是全桥做好呢

0
回复
greenbug
LV.4
20
2023-10-18 15:11

mark下,学习到了,

0
回复
orangi
LV.1
21
2023-12-14 17:06

大佬,表情后面的连接打不开呀

0
回复
orangi
LV.1
22
2023-12-14 17:09
@orangi
大佬,表情后面的连接打不开呀

打开了,谢谢大佬

0
回复
05-15 09:48

值得学习,谢谢作者分享

0
回复
菱奇
LV.4
24
06-12 17:55

0
回复