【正文】 end。 %時段中功角、轉速平均變化率 % Omega_z=(Omega_1+Omega_2)/2。 %時段末功角和轉速估計值 % Omega(i+1)=Omega(i)+Omega_1*h。 end。 %時段中功角、轉速平均變化率 % Omega_z=(Omega_1+Omega_2)/2。 %時段末功角和轉速估計值 % Omega(i+1)=Omega(i)+Omega_1*h。 %初始轉速 t(1)=0。輸入故障切除時間 \n39。 %故障存在時發(fā)電機的最大功率 P3M=。 f=50。 Ke=Ke+。 disp(s)。 PE_result=Eq*U*sin(delta*pi/180)/Xd_all。 if ((V10)amp。 V=eig(A)。 A(2,3)=K2/Tj。 A(1,2)=1。 K6=UGq*(Xd_allXd)/(UG*Xdd_all)。 K2=U*sin(delta*pi/180)/Xdd_all。 UGq=Eq*Xe/Xd_all+Xd*U*cos(delta*pi/180)/Xd_all。 c=Ke^2*Xd^2*U*U/Xd_all^2(Eq0+Ke*UG0)^2。 %慣性時間常數(shù) Td0=10。 %同步電抗 Xdd=。 %系統(tǒng)電壓 Eq0=。 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。2539。%由序分量合成相分量 Iabc_effective=abs(Iabc) G=39。 %考慮 Yd11聯(lián)結組別 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_Node0=U2fault_node0+U2normal_node0。 U2_Node2=zeros(bus_Num2,1)。 %正常分量 U2normal_node1=ones(bus_Num2,1)。 U2fault_node2=zeros(bus_Num2,1)。 else %非故障點短路電流各序分量 I2fault_node1(m,1)=0。 I2fault_node0=zeros(bus_Num2,1)。各支路相電流 39。 I_branch0(i,1)=(U1fault_node0(m,1)+U1normal_node0(m,1)U1fault_node0(n,1)U1normal_node0(n,1))/z0_branch(i)。 for i=1:bus_Num1 m=branch1_No1(i)。 str U1=U1_node1 U2=U1_node2 U0=U1_node0 課程 電 電力系統(tǒng)暫態(tài)上機課程設計報告 18 18 Uabc123=T*[U1_node1 U1_node2 U1_node0].39。 U1_node2=U1fault_node2+U1normal_node2。 U1_node1=zeros(bus_Num1,1)。 U1fault_node0=Y0\I1fault_node0。 U1fault_node2=zeros(bus_Num1,1)。 else %非故障點各序電流各序分量 I1fault_node1(m,1)=0。 I1fault_node0=zeros(bus_Num1,1)。短路相電流 39。 end。 I_Fault0=0。 I_Fault2=I_Fault1。 elseif Fault_Type==1 %a相接地短路各序電流計算 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node))。 T=[1 1 1 % T為對稱分量法的變換矩陣，見 P87公式（ 44） a^2 a 1 a a^2 1]。\n(1)Fault_Type=0為三相短 路 。\n Fault_Node=39。 ? Shortcircuit clc。 YY0(c,c)=YY0(c,c)+a。 c=branch2_No2(i)。 YY0(4,4)=1/(j*)。 end YY2=YY1。 YY1(b,b)=YY1(b,b)+a。 YY1(5,5)=1/(j*)。 z1_branch(5)=j*。 Y0(b,c)=a。 %零序節(jié)點導納矩陣 for i = 1:branch_Num1 a=1/z0_branch(i)。 end Y2=Y1。 Y1(b,b)=Y1(b,b)+a。 Y1(2,2)=1/(j*)。 z0_branch(2)=j*。 z1_branch(2)=j*。 branch1_No2=[2,3,3]。 %節(jié)點數(shù) bus_Num2=5。 暫態(tài)穩(wěn)定是指電力系統(tǒng)受到大的干擾后，各發(fā)電機組保持同步運行并且過渡到新的或恢復到原來穩(wěn)定運行狀態(tài)的能力。 潮流計算需要選定 PQ 和 PV 節(jié)點，短路計算不需要區(qū)分節(jié)點類型。 b. 用計算機方法，進行電力系統(tǒng)潮流計算和短路電流計算，有哪些區(qū)別？ 答： 計算機潮流計算是根據發(fā)電機和線路參數(shù)算出系統(tǒng)中各線路的功率和電壓分布，所計算的數(shù)值實際各節(jié)點電壓和各個線路流過的功率。 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)上機課程設計報告 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 + 。 3. 《電力系統(tǒng)故障的計算機輔助分析》 重慶大學出版社 米麟書等 4. 《電力系統(tǒng)潮流計算》 天津大學出版社 宋文南等 5. 《電力系統(tǒng)故障分析》 清華大學出版社 周榮光 6. 《短路電流實用計算方法》 電力工業(yè)出版社 西安交通大學等 7. 《精通 》北京航空航天大學出版社，張志涌 課程 電 電力系統(tǒng)暫態(tài)上機課程設計報告 3 3 附錄（ 設計流程圖、計算結果、思考題答案 ） 對稱短路計算流程圖及計算結果 ? 過程流程圖 輸入數(shù)據 根據 fault 值選擇故障類型 用節(jié)點導納陣求逆得到節(jié)點阻抗矩陣 選擇短路點f 計算各序電流的故障分量 計算故障點各序電流 計算故障點各序電壓的故障分量 計算故障后各節(jié)點電壓值 形成節(jié)點導納矩陣 確定 fault 值 開始 計算支路電流 計算計及發(fā) 電機端點的各節(jié)點電壓 結束 課程 電 電力系統(tǒng)暫態(tài)上機課程設計報告 4 4 ? 計算結果 節(jié)點導納矩陣 圖 1 Y Y2 的計算結果 圖 2 Y0的計算結果 圖 3 YY YY2的結算結果 (2) 節(jié)點 3 發(fā)生 三 相短路故障 a、故障點三序電流、三相電流： I1= 0 I2=0 I0=0 。在分析過程中通常盡量避免對