【正文】 電力系統(tǒng)暫態(tài)上機課程設(shè)計報告 12 12 ? 計算結(jié)果 切除故障時的搖擺曲線 切除故障時的搖擺曲線 切除故障時，系統(tǒng)為暫態(tài)穩(wěn)定，最大切 此時系統(tǒng)出現(xiàn)失穩(wěn)現(xiàn)象 除時間大約為 。在分析不對稱故障時，各相之間電磁量的耦合使問題的分析更為復(fù)雜，此時常用的分析方法是采用對稱分量法將不對稱問題轉(zhuǎn)化為對稱問題來求解 。 1 課程設(shè)計報告 ( 2020— 2020 年度第二學(xué)期 ) 名 稱： 電力系統(tǒng)暫態(tài)上機 院 系： 電氣與電子工程學(xué)院 班 級： 學(xué) 號： 學(xué)生姓名： 指導(dǎo)教師： 設(shè)計周數(shù)： 兩周 成 績： 日期： 2020 年 7 月 1 日 課程 電 電力系統(tǒng)暫態(tài)上機課程設(shè)計報告 2 2 一、課程 設(shè) 計 的目的與要求 鞏固電力系統(tǒng)暫態(tài)分析理論知識，使學(xué)生掌握采用計算機對電力系統(tǒng)電磁暫態(tài)過程和機電暫態(tài)過程進行計算的基本方法，并進一步鞏固計算機編程能力，為將來從事相關(guān)的技術(shù)工作打下必要的基礎(chǔ)。在分析過程中通常盡量避免對微分方程直接求解，而是采用一定的工具和假設(shè)使問題得以簡化，即把 “ 微分方程代數(shù)化，暫態(tài)分析穩(wěn)態(tài)化 ” 。 Ia= 0 Ib= Ic= + 、三相電壓： 節(jié)點 U1 U2 U0 Ua Ub Uc 1 0 + 2 0 + 3 0 4 0 + + 5 0 + + 、三相電流： 支路 I1 I2 I0 Ia Ib Ic 12 0 + 0 0 0 13 0 0 + 0 0 課程 電 電力系統(tǒng)暫態(tài)上機課程設(shè)計報告 8 8 23 0 0 + 0 0 14 0 + 0 + 0 0 25 0 + 0 0 0 發(fā)電機 1 + 0 發(fā)電機 3 + 0 (4) 節(jié)點 3 發(fā)生 b、 c 兩相短路接地故障 a. 故障點三序電流、三相電流： I1= 0 I2= 0 + I0= 0 + 。 潮流計算需要選定 PQ 和 PV 節(jié)點，短路計算不需要區(qū)分節(jié)點類型。 %節(jié)點數(shù) bus_Num2=5。 z1_branch(2)=j*。 Y1(2,2)=1/(j*)。 end Y2=Y1。 Y0(b,c)=a。 YY1(5,5)=1/(j*)。 end YY2=YY1。 c=branch2_No2(i)。 ? Shortcircuit clc。\n(1)Fault_Type=0為三相短路 。 elseif Fault_Type==1 %a相接地短路各序電流計算 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node))。 I_Fault0=0。短路相電流 39。 else %非故障點各序電流各序分量 I1fault_node1(m,1)=0。 U1fault_node0=Y0\I1fault_node0。 U1_node2=U1fault_node2+U1normal_node2。 for i=1:bus_Num1 m=branch1_No1(i)。各支路相電流 39。 else %非故障點短路電流各序分量 I2fault_node1(m,1)=0。 %正常分量 U2normal_node1=ones(bus_Num2,1)。 U2_Node0=U2fault_node0+U2normal_node0。%由序分量合成相分量 Iabc_effective=abs(Iabc) G=39。 G I2_G1=I1_branch*exp(j*pi/6) I2_G2=I2_branch*exp(j*pi/6) I2_G0=0 Iabc=T*[I1_branch*exp(j*pi/6) I2_branch*exp(j*pi/6) 0].39。 %同步電抗 Xdd=。 c=Ke^2*Xd^2*U*U/Xd_all^2(Eq0+Ke*UG0)^2。 K2=U*sin(delta*pi/180)/Xdd_all。 A(1,2)=1。 V=eig(A)。 PE_result=Eq*U*sin(delta*pi/180)/Xd_all。 Ke=Ke+。 %故障存在時發(fā)電機的最大功率 P3M=。 %初始轉(zhuǎn)速 t(1)=0。 %時段中功角、轉(zhuǎn)速平均變化率 % Omega_z=(Omega_1+Omega_2)/2。 %時段末功角和轉(zhuǎn)速估計值 % Omega(i+1)=Omega(i)+Omega_1*h。 end。 %時段中功角、轉(zhuǎn)速平均變化率 % Omega_z=(Omega_1+Omega_2)/2。 end。 %時段末功角和轉(zhuǎn)速估計值 % Omega(i+1)=Omega(i)+Omega_1*h。輸入故障切除時間 \n39。 f=50。 disp(s)。 if ((V10)amp。 A(2,3)=K2/Tj。 K6=UGq*(Xd_allXd)/(UG*Xdd_all)。 UGq=Eq*Xe/Xd_all+Xd*U*cos(delta*pi/180)/Xd_all。 %慣性時間常數(shù) Td0=10。 %系統(tǒng)電壓 Eq0=。2539。 %考慮 Yd11聯(lián)結(jié)組別 U4abc_effective=abs(U4abc) m=5 U5abc=T*[U2_Node1(5,1)*exp(j*pi/6) U2_Node2(5,1)*exp(j*pi/6) 0].39。 U2_Node2=zeros(bus_Num2,1)。 U2fault_node2=zeros(bus_Num2,1)。 I2fault_node0=zeros(bus_Num2,1)。 I_branch0(i,1)=(U1fault_node0(m,1)+U1normal_node0(m,1)U1fault_node0(n,1)U1normal_node0(n,1))/z0_branch(i)。 str U1=U1_node1 U2=U1_node2 U0=U1_node0 課程 電 電力系統(tǒng)暫態(tài)上機課程設(shè)計報告 18 18 Uabc123=T*[U1_node1 U1_node2 U1_node0].39。 U1_node1=zeros(bus_Num1,1)。 U1fault_node2=zeros(bus_Num1,1)。 I1fault_node0=zeros(bus_Num1,1)。 end。 I_Fault2=I_Fault1。 T=[1 1 1 % T為對稱分量法的變換矩陣，見 P87公式（ 44） a^2 a 1 a a^2 1]。\n Fault_Node=39。 YY0(c,c)=YY0(c,c)+a。 YY0(4,4)=1/(j*)。 YY1(b,b)=YY1(b,b)+a。 z1_branch(5)=j*。 %零序節(jié)點導(dǎo)納矩陣 for i = 1:branch_Num1 a=1/z0_branch(i)。 Y1(b,b)=Y1(b,b)+a。 z0_branch(2)=j*。 branch1_No2=[2,3,3]。 暫態(tài)穩(wěn)定是指電力系統(tǒng)受到大的干擾后，各發(fā)電機組保持同步運行并且過渡到新的或恢復(fù)到原來穩(wěn)定運行狀態(tài)的能力。 b. 用計算機方法，進行電力系統(tǒng)潮流計算和短路電流計算，有哪些區(qū)別？ 答： 計算機潮流計算是根據(jù)發(fā)電機和線路參數(shù)算出系統(tǒng)中各線路的功率和電壓分布，所計算的數(shù)值實際各節(jié)點電壓和各個線路流過的功率。 3. 《電力系統(tǒng)故障的計算機輔助分析》 重慶大學(xué)出版社 米麟書等 4. 《電力系統(tǒng)潮流計算》 天津大