电力电子与电力系统报告-兰州交通大学

（电力电子技术与电力系统分析）课程实训报告

附录I同步发电机突然三相短路暂态过程的仿真

MATLAB/Simulink下系统模型图及运行结果

程序

%%****************************************************************** N=48;

t1=(0:0.02/N:1.00);

fai=0*pi/180;%空载短路全电流表达式

Ia=(-cos(2*pi*50*t1+fai)-1.43*exp(-2.97*t1).*cos(2*pi*50*t1+fai)-2.34*exp(-0.608*t1).*cos(2*pi*50*t1+fai)+4*exp(-6.3*t1).*cos(-fai*pi/180)+0.77*exp(-6.3*t1).*cos(2*2*pi*50*t1+fai));%基频分量

Ia1=-cos(2*pi*50*t1+fai)-1.43*exp(-2.97*t1).*cos(2*pi*50*t1+fai)-2.34*exp(-0.608*t1).*cos(2*pi*50*t1+fai);%倍频分量

Ia2=0.77*exp(-6.3*t1).*cos(2*2*pi*50*t1+fai);%非周期分量 Iap=4*exp(-6.3*t1).*cos(-fai*pi/180);

subplot(4,1,1); %绘制空载短路全电流波形图 plot(t1,Ia); grid on;

axis([0 1 -10 10]); ylabel('Ia(p.u.)');

subplot(4,1,2); %绘制基频分量波形图 plot(t1,Ia1); grid on;

axis([0 1 -10 10]); ylabel('Ia1(p.u.)');

subplot(4,1,3); %绘制倍频分量波形图 plot(t1,Ia2); grid on; axis([0 1 -1 1]); ylabel('Ia2(p.u.)');

subplot(4,1,4); %绘制非周期分量波形图

- 12 -

（电力电子技术与电力系统分析）课程实训报告

plot(t1,Iap); grid on;

axis([0 1 -10 10]); ylabel('Iap(p.u.)'); xlabel('t/s');

程序运行结果：

(1)发电机端突然三相短路时的a相定子电流仿真波形：

图1发电机端突然发生三相短路时的a相定子电流波形图

(2)发电机端突然三相短路时的定子电流仿真波形：

图2 发电机端突然三相短路时的定子电流仿真波形图

- 13 -

（电力电子技术与电力系统分析）课程实训报告

(3)发电机端突然三相短路时id、iq以及if 的电流仿真波形：

图3 发电机端突然三相短路时id、iq以及if 的电流仿真波形图

- 14 -