電力系統(tǒng)暫態(tài)上機(jī)計(jì)算課程設(shè)計(jì)報(bào)告-資料下載頁(yè)
2025-08-03 06:10本頁(yè)面
【正文】 nch*exp(j*pi/6)I2_G2=I2_branch*exp(j*pi/6)I2_G0=0Iabc=T*[I1_branch*exp(j*pi/6) I2_branch*exp(j*pi/6) 0].39。%由序分量合成相分量 考慮Yd11聯(lián)結(jié)組別Iabc_effective=abs(Iabc)l stability_smallsignalclear。clc。 U=1。 %系統(tǒng)電壓Eq0=。 %空載電動(dòng)勢(shì)UG0=。 %機(jī)端電壓Xe=。 %線路電抗值Xd=。 %同步電抗Xdd=。 %暫態(tài)電抗Xd_all=。 %系統(tǒng)電抗(包含同步電抗)Xdd_all=。 %系統(tǒng)電抗(包含暫態(tài)電抗)Tj=10。 %慣性時(shí)間常數(shù)Td0=10。 %勵(lì)磁繞組時(shí)間常數(shù)Ke=。while(Ke) for delta=49:110 %各個(gè)參數(shù) a=(Ke*Ke*Xe*Xe)/(Xd_all*Xd_all)1。 b=2*(Ke^2*Xd*Xe*U*cos(delta*pi/180)/(Xd_all^2)+Eq0+Ke*UG0)。 c=Ke^2*Xd^2*U*U/Xd_all^2(Eq0+Ke*UG0)^2。 Eq=(b+sqrt(b^24*a*c))/(2*a)。 Eqq=Eq*Xdd_all/Xd_all+(Xd_allXdd_all)*U*cos(delta*pi/180)/Xd_all。 UGd=U*Xd*sin(delta*pi/180)/Xd_all。 UGq=Eq*Xe/Xd_all+Xd*U*cos(delta*pi/180)/Xd_all。 UG=sqrt(UGd^2+UGq^2)。 A=zeros(3)。 %P186公式741狀態(tài)方程各個(gè)參數(shù) K1=Eqq*U*cos(delta*pi/180)/Xdd_all+U*U*(Xdd_allXd_all)*cos(2*delta*pi/180)/(Xd_all*Xdd_all)。 K2=U*sin(delta*pi/180)/Xdd_all。 K3=Xdd_all/Xd_all。 K4=(Xd_allXdd_all)*U*sin(delta*pi/180)/Xdd_all。 K5=UGd*Xd*cos(delta*pi/180)/(UG*Xd_all)UGq*Xdd*sin(delta*pi/180)/(UG*Xdd_all)。 K6=UGq*(Xd_allXd)/(UG*Xdd_all)。 Kemax=K4/K5。 Kemin=Eq0*cos(delta/180*pi)/(K3*Xd_all*(K1*K6K2*K5))。 A(1,1)=0。 A(1,2)=1。 A(1,3)=0。 A(2,1)=K1/Tj。 A(2,2)=0。 A(2,3)=K2/Tj。 A(3,1)=(K4+Ke*K5)/Td0。 A(3,2)=0。 A(3,3)=(1/K3+Ke*K6)/Td0。 V=eig(A)。 V1=real(V(1))。 V2=real(V(2))。 V3=real(V(3))。 if ((V10)amp。(V20)amp。(V30)) Ke_result=Ke。 delta_result=delta。 PE_result=Eq*U*sin(delta*pi/180)/Xd_all。 else s=sprintf(39。Ke=%f, delta=%d, P=%f\n39。,Ke_result,delta_result,PE_result)。 disp(s)。 break。 end。 end。 Ke=Ke+。 end。l stability_transientclear。clc。 f=50。 %系統(tǒng)額定頻率Tj=。 %歸算后的發(fā)電機(jī)慣性時(shí)間常數(shù)PT=1。 %正常運(yùn)行時(shí)發(fā)電機(jī)向無(wú)窮大系統(tǒng)傳輸?shù)挠泄β蔖2M=。 %故障存在時(shí)發(fā)電機(jī)的最大功率P3M=。 %故障切除后發(fā)電機(jī)的最大功率 %下面是利用改進(jìn)歐拉法進(jìn)行逐段計(jì)算%需要注意故障切除前后電磁功率有躍變h=。 %Duration=2。 %設(shè)置計(jì)算時(shí)段長(zhǎng)度2s CutTime=input(39。輸入故障切除時(shí)間\n39。)。 %故障發(fā)生時(shí)的功角變化過(guò)程Delta(1)=。 %初始功角Omega(1)=1。 %初始轉(zhuǎn)速t(1)=0。for i=1:round(CutTime/h) Delta_1=(Omega(i)1)*360*f。 %第一時(shí)段開(kāi)始時(shí)功角和轉(zhuǎn)速變化率% Omega_1=(PTP2M*sin(Delta(i)*pi/180))/Tj。 Delta(i+1)=Delta(i)+Delta_1*h。 %時(shí)段末功角和轉(zhuǎn)速估計(jì)值% Omega(i+1)=Omega(i)+Omega_1*h。 Delta_2=(Omega(i+1)1)*360*f。 %時(shí)段末估計(jì)值的變化率% Omega_2=(PTP2M*sin(Delta(i+1)*pi/180))/Tj。 Delta_z=(Delta_1+Delta_2)/2。 %時(shí)段中功角、轉(zhuǎn)速平均變化率% Omega_z=(Omega_1+Omega_2)/2。 Delta(i+1)=Delta(i)+Delta_z*h。 %第一時(shí)段末功角、轉(zhuǎn)速的值% Omega(i+1)=Omega(i)+Omega_z*h。 t(i+1)=i*h。end。%故障切除后的功角變化過(guò)程for i=round(CutTime/h)+1:round(Duration/h) Delta_1=(Omega(i)1)*360*f。 %第一時(shí)段開(kāi)始時(shí)功角和轉(zhuǎn)速變化率% Omega_1=(PTP3M*sin(Delta(i)*pi/180))/Tj。 Delta(i+1)=Delta(i)+Delta_1*h。 %時(shí)段末功角和轉(zhuǎn)速估計(jì)值% Omega(i+1)=Omega(i)+Omega_1*h。 Delta_2=(Omega(i+1)1)*360*f。 %時(shí)段末估計(jì)值的變化率% Omega_2=(PTP3M*sin(Delta(i+1)*pi/180))/Tj。 Delta_z=(Delta_1+Delta_2)/2。 %時(shí)段中功角、轉(zhuǎn)速平均變化率% Omega_z=(Omega_1+Omega_2)/2。 Delta(i+1)=Delta(i)+Delta_z*h。 %第一時(shí)段末功角、轉(zhuǎn)速的值% Omega(i+1)=Omega(i)+Omega_z*h。 t(i+1)=i*h。end。 plot(t,Delta)。
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1