电力系统暂态稳定性仿真

在图(1)中,发电机(电抗Xd = 0.25 pu)连接到母线(4)。 发电机的端电压调节为1.2 pu。 该发电机通过升压电源变压器(漏电抗Xtr = 0.2 pu)连接到母线(1)。 总线(1)通过两条相同的平行传输线(每单位电抗长度为0.3174每公里,长为100 Km)连接到无限总线(2)。 该系统的基本电压和MVA为345 KV和1500 MVA

分别。 最初,发电机的功率为1500 MW。 一根传输线上距母线(1)20公里的距离处发生三相对称故障。 一段时间后,该故障将清除。 注意,在这个问题中,由于饱和效应(由于电枢电流的上升),在故障期间反电动势E和Xd的大小变化可以忽略。 摆动方程使用以下参数:H = 3.5 pu,!syn = 2 60,阻尼系数D = 7.5.1

用matlab和powerworld分别进行了暂态仿真,已经分析了一个简单的问题,以便了解电源系统中的瞬态稳定性。 这是可以分析解决的最简单的问题。

了解详情请下载附件http://www.apollocode.net/a/448.html

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
clear; clc; f=50; %系统额定频率 Tj=8.47; %归算后的发电机惯性时间常数 PT=1; %正常运行时发电机向无穷大系统传输的有功功率 P2M=0.48; %故障存在时发电机的最大功率 P3M=1.38; %故障切除后发电机的最大功率 %下面是利用改进欧拉法进行逐段计算 %需要注意故障切除前后电磁功率有跃变 h=0.05; %设置步长0.05s Duration=2; %设置计算时段长度2s Delta(1)=33.92; %初始功角 Omega(1)=1; %初始转速 t(1)=0; Delta_h=pi-asin(1/1.38); Delta_cm=acos((PT*(Delta_h-Delta(1)*pi/180)+P3M*cos(Delta_h)-P2M*cos(Delta(1)*pi/180))/(P3M-P2M))*180/pi; d=Delta(1); for i=1:round(Duration/h) if d<Delta_cm d_Delta(i)=(Omega(i)-1)*360*f; d_Omega(i)=(PT-P2M*sin(Delta(i)*pi/180))/Tj; Delta0(i+1)=Delta(i)+d_Delta(i)*h; Omega0(i+1)=Omega(i)+d_Omega(i)*h; d_Delta0(i+1)=(Omega0(i+1)-1)*f*360; d_Omega0(i+1)=(PT-P2M*sin(Delta0(i+1)*pi/180))/Tj; d_Deltaa(i+1)=(d_Delta(i)+d_Delta0(i+1))/2; d_Omegaa(i+1)=(d_Omega(i)+d_Omega0(i+1))/2; Delta(i+1)=Delta(i)+d_Deltaa(i+1)*h; Omega(i+1)=Omega(i)+d_Omegaa(i+1)*h; d=Delta(i+1); t(i+1)=i*h; T=t(i); end; end; s=sprintf(' 最大摇摆角 Delta_h=%f\n 最大切除角 Delta_cm=%f\n 最大稳定切除时间 Tmax=%f \n',Delta_h*180/pi,Delta_cm,T); disp(s); CutTime=input('输入故障切除时间\n'); %故障发生时的功角变化过程 Delta(1)=33.92; %初始功角 Omega(1)=1; %初始转速 t(1)=0; for i=1:round(CutTime/h) d_Delta(i)=(Omega(i)-1)*360*f; d_Omega(i)=(PT-P2M*sin(Delta(i)*pi/180))/Tj; Delta0(i+1)=Delta(i)+d_Delta(i)*h; Omega0(i+1)=Omega(i)+d_Omega(i)*h; d_Delta0(i+1)=(Omega0(i+1)-1)*f*360; d_Omega0(i+1)=(PT-P2M*sin(Delta0(i+1)*pi/180))/Tj; d_Deltaa(i+1)=(d_Delta(i)+d_Delta0(i+1))/2; d_Omegaa(i+1)=(d_Omega(i)+d_Omega0(i+1))/2; Delta(i+1)=Delta(i)+d_Deltaa(i+1)*h; Omega(i+1)=Omega(i)+d_Omegaa(i+1)*h; t(i+1)=i*h; end; %故障切除后的功角变化过程 for i=round(CutTime/h)+1:round(Duration/h) d_Delta(i)=(Omega(i)-1)*360*f; d_Omega(i)=(PT-P3M*sin(Delta(i)*pi/180))/Tj; Delta0(i+1)=Delta(i)+d_Delta(i)*h; Omega0(i+1)=Om
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值