【正文】 functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。此外兩個軟件采用的算法不同，對精確的要求也不一樣，所以也會造成一定的計算差距。 采用 PowerWorld軟件分析結(jié)果在PowerWorld軟件里給發(fā)電機建立暫態(tài)穩(wěn)定模型，其中平衡節(jié)點用一臺大發(fā)電機代替，選定GENPWTwoAxis模型，把其中的H值設(shè)得非常大（），其他都用默認(rèn)參數(shù)。%短路后的電磁功率[t,x]=ode45(39。 1/*(*sin(x(1)+c)10*(x(2)1))]。r39。%短路后的P3，如式15[t2,x2]=ode45(39。c=*pi/180。 1/*(*sin(x(1)+c)10*(x(2)1))]。由文中第4部分知，E’=，V=，穩(wěn)定運行時的轉(zhuǎn)移阻抗Z17= +=，輸入阻抗Z77= +j =。在網(wǎng)絡(luò)變換中，用電壓為E’、內(nèi)阻為Xd’的電壓源來等效原來的三臺機組。說明兩種計算方法都是準(zhǔn)確可行的。具體參數(shù)設(shè)置如圖4所示。它能夠有效地計算高達10,0000個節(jié)點的電力網(wǎng)絡(luò)，因此當(dāng)它作為一個獨立的潮流分析軟件包時，性非常實用。我們采用PQ分解法對待求解電網(wǎng)進行潮流計算，大致流程如下：1）編程思想是將所有支路的狀態(tài)信息整合成一個矩陣進行計算、將所有節(jié)點的狀態(tài)信息整合成另一個矩陣進行計算，具體如下：構(gòu)造支路信息矩陣B1，每一行表示一條支路，一行中的各元素表示該支路的信息：支路首端號；末端號；支路阻抗；支路對地電納 支路的變比；支路首端處于K側(cè)為1，首端處于1側(cè)則為0；構(gòu)造節(jié)點信息矩陣B2，每一行表示一個節(jié)點，一行中的各元素表示該節(jié)點的信息：該節(jié)點發(fā)電機功率；該節(jié)點負荷功率；節(jié)點電壓初始值。 穩(wěn)定計算的功角曲線等；178。用PowerWorld軟件進行潮流計算并與自己編制的軟件計算結(jié)果進行校核和分析；用Matlab編制穩(wěn)定計算程序（三臺機可并聯(lián)等值成一臺機），發(fā)電機采用二階經(jīng)典模型，要求給出網(wǎng)絡(luò)變換法求解輸入阻抗和轉(zhuǎn)移阻抗的變換過程圖；自行選擇23種故障方案，給出搖擺曲線，并計算故障的極限切除時間和極限切除角。發(fā)電廠升壓變參數(shù)均為Vs%=%。圖中，各變電站參數(shù)如下表：編號類型220kV最大負荷，MVA1500kV站平衡節(jié)點2220kV站250+j403220kV站380+j904220kV站280+j855220kV站420+j95各線路長度如圖所示。燃煤發(fā)電廠G有三臺機組，均采用單元接線。穩(wěn)定仿真中不考慮發(fā)電廠的調(diào)速器和原動機模型。 以上設(shè)計任務(wù)每一部分的計算過程和結(jié)果分析；178。此外，節(jié)點6為PV節(jié)點，其中總發(fā)出的有功功率（5）。該仿真軟件能夠進行專業(yè)的工程分析。 PowerWorld軟件建模與潮流計算 采用PowerWorld仿真軟件建立潮流模型在軟件的編輯模式下，節(jié)點1為平衡節(jié)點，為了反映平衡節(jié)點這一特性，節(jié)點1應(yīng)接一臺容量無窮大的發(fā)電機。則發(fā)電機輸出功率PD=3*+=總負荷PL =250MW+380MW+280MW+420MW=1330MW PowerWorld仿真結(jié)果與Matlab計算結(jié)果對比分析以上，我們分別用PowerWorld仿真軟件和Matlab平臺對待求解電網(wǎng)進行了潮流計算，得出各節(jié)點的節(jié)點電壓和各支路的潮流狀態(tài)。在本例中，正常運行下負荷所在節(jié)點電壓近似為VN，由式9求得各負荷的等效阻抗如表6所示：表1各負荷的等效阻抗負荷所在節(jié)點節(jié)點電壓V（kV）負荷大小S（MVA）等效阻抗ZL的標(biāo)幺值2220250+j40+3220380+j90+4220280+j85+5220420+j95+ 發(fā)電機的處理接在節(jié)點6的有三臺發(fā)電機，發(fā)電機的容量SG(N)=300/=，VG(N)=，而基準(zhǔn)容量SB=100MW，為了使變壓器的標(biāo)幺變比k*=1，在10kV電壓等級下的基準(zhǔn)電壓VB2=*220/242kV=。（13）代入數(shù)據(jù)后得（14）則有。由于Pe的表達式在不同的狀態(tài)下（故障前、故障后、切除后）是不一樣的，所以需要Matlab定義多個時間段進行仿真。w0]。%給第二個計算區(qū)間賦初值a=。x=[x1。圖1tc=(1) 采用試湊法不斷改變切除時間，即給tc賦值，找到穩(wěn)定與不穩(wěn)定的臨界時間，如圖21所示，取tc=。%給定微分方程組的初值a=。)%畫出搖擺曲線line([,],[0,3]) %畫t=圖2故障時功角隨時間變化的特性已知極限切除時間tc=，則可由圖22的曲線上找到對應(yīng)的切除角，即極限切除角。圖28綜合上面4點的仿真結(jié)果，=。所以，為了檢驗系統(tǒng)的暫態(tài)穩(wěn)定特性，文本所提到的線路故障皆為雙回路在同一點發(fā)生短路，即線路L23的兩個回路同時在中點處發(fā)生三相對稱短路。x0=[*pi/180。%求解第一個區(qū)間[0,tc]x02= x1 (end,:)39。t2]。圖tc=，修改相應(yīng)的Matlab程序，得到的搖擺曲線如下，即此時發(fā)電暫態(tài)不穩(wěn)定。圖36綜合上面4點的仿真結(jié)果，=。 故障方案二由第5部分知，L23中點處發(fā)生三相對稱短路，當(dāng)故障切除后。圖42（4）在5s時節(jié)點3發(fā)生三相對稱短路，所得搖擺曲線如下，即此時系統(tǒng)暫態(tài)不穩(wěn)定。所得搖擺曲線如下，此時系統(tǒng)暫態(tài)不穩(wěn)定。但總的來說兩種發(fā)電機模型的極限切除時間相差并不大。在PowerWorld軟件里，發(fā)電機采用三階模型，勵磁系統(tǒng)采用BPA_EA模型，TJ分別取50時，分別求出系統(tǒng)的極限切除角，其中故障選取故障方案1：節(jié)點3發(fā)生三相對稱短路。計算結(jié)果如表18：表1定量分析TR對穩(wěn)定計算結(jié)果的影響Tr012(s)圖5(s)隨勵磁參數(shù)Tr的變化特性從圖30中可知，當(dāng)勵磁參數(shù)Tr變大，極限切除角也變大。B2=[0 0 1 0 1。3 0。請輸入支路數(shù)：nl=39。 B1=input(39。)。V=zeros(1,n)。q=B1(i,1)。YI(p,p)=YI(p,p)+1./B1(i,3)。fori=1:n e(i)=real(B2(i,3))。endB(i,i)=1./B(i,i)。endendendendfori=1:na ifi==naA(i,i)=1./A(i,i)。kq=1。endDP1(i)=P(i)V(i)*C(i)。 endelsefor LZ=3:i L=i+3LZ。L=L+1。 L=0。if ICT3~=0 L=0。endendend L=0。endenddisp(39。disp(39。)。)。disp(S)。 end Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。q=B1(i,1)。else p=B1(i,2)。)。47。迭代次數(shù)39。 endDS(i)=Si(p,q)+Sj(q,p)。disp(Sj(q,p))。enddisp(39。各條支路的首端功率Sj（順序同您輸入B1時一樣）為：39。for p=1:nC(p)=0。disp(39。)。)。V(i)=V(i)DQ(L)。ifi==nafor LZ=2:i L=i+2LZ。L=L+1。endDQ1(i)=Q(i)V(i)*C(i)。for MZ=2:IC4 I=IC4+2MZ。 DET=abs(DP1(i))。DET=0。for j1=k:naA(i,j1)=A(i,j1)./A(i,i)。endendendendp=0。V(i)=B2(i,4)。B=imag(YI)。YI(p,q)=YI(p,q)1./B1(i,3)。fori=1:n if X(i,2)~=0。請輸入PQ節(jié)點數(shù)na=39。)。 isb=input(39。5 0。0 +*i 1 0 0 2。參考文獻[1] 《電力系統(tǒng)分析》（上下冊）華中科技大學(xué)出版[2] 《發(fā)電廠電氣部分》中國電力出版社[3] PowerWorld 16使用手冊[4] 基于Matlab/Simulink的系統(tǒng)仿真技術(shù)與應(yīng)用附錄：Matlab編制PQ分解法潮流計算程序：本文中的實例數(shù)據(jù)如下：節(jié)點數(shù)為6；支路數(shù)為6；平衡母線節(jié)點號為1；；PQ節(jié)點數(shù)為4；B1=[1 2 +i** 1 0。當(dāng)TJ足夠大后，再增大TJ時，極限切除角變化幅度不大。表16為兩種不同的發(fā)電機模型在不同故障方案下的暫態(tài)過渡時間統(tǒng)計；表1不同的發(fā)電機模型在不同故障方案下的暫態(tài)過渡時間發(fā)電機模型故障方案一故障方案二故障切除時間（s）二階模型暫態(tài)過渡時間（s）454748424344三階模型暫態(tài)過渡時間（s）141516131415從表中可以看出，對于同一模型，故障切除時間越大，暫態(tài)過渡時間越大。所得搖擺曲線如下，此時系統(tǒng)