孟祥萍电力系统分析教案20.pptVIP

第20章 电力系统稳定计算程序及其实现 20.1 简化模型的暂态稳定计算 20.1.1程序及框图 20.1.2 上机说明 20.1.3 实例 20.2 简化模型的静态稳定计算 20.2.1 程序及框图 20.2.2 上机说明 20.2.3 实例 小　　结 电力系统分析 本章提示 20.1简化模型的暂态稳定计算 20.2简化模型的静态稳定计算 小 结 本章提示 电力系统暂态稳定的计算机算法； 电力系统静态稳定的计算机算法。 20.1.1 程序及框图 20.1.2 上机说明 20.1.3 实例 开 始 输入原始数据 计算出功角及初始电势E0 判断此时系统有无扰动，根据扰动性质选择特性曲线 K=2 打印 K是否是要求的时段 K=K+1 no yes 结束 简化模型暂 态稳定计算 框图 由框图可编写程序如下： %本程序是用分段计算法计算电力系统的暂态稳定。 S0=input(请输入初始功率：S0=); V0=input(请输入无限大系统母线电压：V0=); B=input(请输入系统等值电抗矩阵：B=); Tj=input(请输入惯性时间常数：Tj=); N=input(请输入时段数：N=); Ni=input(请输入哪个时段发生故障：Ni=); dt=input(请输入每段间隔的时间：dt=); for i=1:Ni C(i,1)=1;C(i,2)=B(1); end for i=Ni:N C(i,1)=0;C(i,2)=B(2); 20.1.1程序及框图 end E0=sqrt((V0+imag(S0)*B(1)./V0)^2+(real(S0)*B(1)./V0)^2); dtj0=atan(real(S0)*B(1)./(V0*(V0+imag(S0)*B(1)./V0))); Q(1)=dtj0*180./pi; K1=360*50*dt^2./Tj; P2m=E0*V0./B(2); P3m=E0*V0./B(3); dtk=pi-asin(real(S0)./P3m); dtjx=acos((real(S0)*(dtk-dtj0)+P3m*cos(dtk)-P2m*cos(dtj0))./(P3m-P2m)); dtjx=dtjx*180./pi; for K=2:N+1 P(K-1)=E0*V0./C(K-1,2); if C(K-1,1)==1 %故障发生之前 Peq(K-1)=P(K-1)*sin(Q(K-1)*pi./180); DP(K-1)=real(S0)-Peq(K-1); if K==2 %第一时段 20.1.1程序及框图 DQ(K)=K1*DP(K-1)./2; else DQ(K)=DQ(K-1)+K1*DP(K-1)./2; end Q(K)=Q(K-1)+DQ(K); elseif C(K-1,1)==0%切除故障时 P1(K-1)=P(K-2); P2(K-1)=P(K-1); Peq(K-1)=P1(K-1)*sin(Q(K-1)); DP(K-1)=real(S0)-Peq(K-1); DP2(K-1)=real(S0)-P2(K-1)*sin(Q(K-1)); DQ(K)=DQ(K-1)+K1*(DP(K-1)+DP2(K-1))./2; Q(K)=Q(K-1)+DQ(K); end end 20.1.1程序及框图 disp(程序计算结果如下：) disp(输出不平衡功率：DP=) disp(DP) disp(输出功角增量：DQ=) disp(DQ); disp(输出功角：Q=) disp(Q) disp(输出电磁功率：Peq=) disp(Peq) disp(输出发电机转子摇摆曲线如图Figure No.1所示:) t=0:dt:dt*N; plot(t,Q,t,dtjx);xlabel(时间（秒）);ylabel(功角(度数));title(发电机摇摆曲线); 20.1.1程序及框图 1请输入初始功率S0,形如a+jb 2请输入无限大系统母线电压U0 3请输入系统等值电抗矩阵B 矩阵B是由以下元素组成的行矩阵： 正常运行时的系统直轴等值电抗Xd 故