【正文】 暫態(tài)穩(wěn)定。lim=。從總體上看，三種方法求出的結(jié)果較為接近。=%5%（6）通過求解故障后功角隨時間的變化曲線。圖3tc=同理。x2]。b=。%給定微分方程組的初值a=。代入式13后得：（18）當(dāng)節(jié)點3的故障切除后，通過網(wǎng)絡(luò)變換法同理可得，轉(zhuǎn)移阻抗Z17= + =，輸入阻抗Z77= + =。由于采用Matlab仿真時需要計算各種狀態(tài)下的轉(zhuǎn)移阻抗與輸入阻抗，而且為復(fù)數(shù)運算，所以存在一點的誤差。點擊圖中交點后，輸出x =，y =，即==*180/pi=。c=*pi/180。其中的Matlab求解程序如下：（子程序與穩(wěn)定計算程序一樣）functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。plot(t,x(:,1),39。c=*pi/180。b=。假設(shè)在5s時刻節(jié)點3發(fā)生了三相對稱短路，則采用Matlab編制的程序如下所示：子函數(shù):functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。5 故障方案仿真 選定故障1：節(jié)點3發(fā)生三相對稱短路 參數(shù)的求取慣性時間常數(shù)TJN=7，歸算到全網(wǎng)后，三臺機組等效為一臺機組后。等效后的機組出力P=，Xd=，Xq=，Xd’=發(fā)電機采用二階經(jīng)典模型，即E’恒定模型。所以，兩種計算潮流得出的節(jié)點電壓具有很好的一致性。圖接在平衡節(jié)點處的發(fā)電機參數(shù)設(shè)置節(jié)點6為PV節(jié)點，P=，V=，為了方便后面的穩(wěn)定計算，按題目中要求，節(jié)點6接3臺參數(shù)一樣的機組，其中正常運行時各機組發(fā)出的有功功率P1=P2=P3=。該仿真器是一個集成的產(chǎn)品，其核心是一個全面、強大的潮流計算程序。基于MATLAB平臺編寫程序，我們對待求解電網(wǎng)進行潮流計算，從而得出各PQ節(jié)點的節(jié)點電壓以及各線路的線路損耗。 潮流計算結(jié)果（潮流圖）178。二、 設(shè)計的主要內(nèi)容：進行參數(shù)計算和標(biāo)幺化，形成潮流計算參數(shù)；用Matlab編制潮流計算程序，要求采用PQ分解法潮流計算方法。發(fā)電機組主要參數(shù)如下表：機組臺數(shù)單臺容量（MW）額定電壓（EV）功率因數(shù)升壓變?nèi)萘縈VAXdXd’XqTd0’TJ=2H330035087當(dāng)發(fā)電機采用三階模型時，勵磁環(huán)節(jié)（含勵磁機和勵磁調(diào)節(jié)器）模型如下:圖中參數(shù)如下：TR=0，KA=20，TA=，Te=，KE=1，SE=0，KF=，TF=不考慮PSS，即圖中的附加信號Vs=0。電力系統(tǒng)穩(wěn)定分析和計算設(shè)計報告電力系統(tǒng)穩(wěn)定分析和計算設(shè)計報告專業(yè)：電氣工程及其自動化班級： 09電氣（4）班學(xué)號：32號、34號姓名：禤培正、陳子程華南理工大學(xué)電力學(xué)院2012229目錄課程設(shè)計題目：電力系統(tǒng)穩(wěn)定分析和計算 11 潮流的參數(shù)計算和標(biāo)幺化 4 線路參數(shù)計算與標(biāo)幺化 4 節(jié)點信息統(tǒng)計 52 采用MATLAB計算電網(wǎng)潮流 6 PQ分解法潮流計算 6 MATLAB平臺潮流計算結(jié)果分析 73 采用powerworld仿真軟件計算潮流 7 PowerWorld仿真軟件簡介 7 PowerWorld軟件建模與潮流計算 8 PowerWorld仿真結(jié)果與Matlab計算結(jié)果對比分析 114 網(wǎng)絡(luò)變換法求解輸入阻抗與轉(zhuǎn)移阻抗 12 負(fù)荷、發(fā)電機、變壓器、線路等效電容的處理 12 采用網(wǎng)絡(luò)變換法的具體步驟 15 求解輸入阻抗和轉(zhuǎn)移阻抗 19 求解有功功率傳輸特性 195 故障方案仿真 20 選定故障1：節(jié)點3發(fā)生三相對稱短路 20 選定故障2：線路L23中點處發(fā)生三相對稱短路 266 采用等面積法求解極限切除角 31 故障方案一 31 故障方案二 327 基于PowerWorld，分析三階模型的發(fā)電機的暫態(tài)過程 33 故障方案一：節(jié)點3發(fā)生三相對稱短路 33 故障方案二：線路L23中點處發(fā)生三相對稱短路 358 發(fā)電機模型及勵磁調(diào)節(jié)參數(shù)對穩(wěn)定計算結(jié)果的影響 36 發(fā)電機模型對穩(wěn)定計算結(jié)果的影響 36 TJ對穩(wěn)定計算結(jié)果的影響 38 勵磁參數(shù)對穩(wěn)定計算結(jié)果的影響 39參考文獻 40附錄： 40[鍵入文字]課程設(shè)計題目：電力系統(tǒng)穩(wěn)定分析和計算姓名：指導(dǎo)教師：一、 一個220kV分網(wǎng)結(jié)構(gòu)和參數(shù)如下：500kV站（1）的220kV母線視為無窮大母線，電壓恒定在230kV。電廠220kV側(cè)采用雙母接線。負(fù)荷采用恒阻抗模型。 所編制的潮流和穩(wěn)定計算源程序（主要語句應(yīng)加注釋）；178。綜上，各節(jié)點的參數(shù)如表2所示：2 采用MATLAB計算電網(wǎng)潮流通過對待求解電網(wǎng)的等效線路參數(shù)的計算與標(biāo)幺化，得出待求解電網(wǎng)的各PQ節(jié)點的有功功率和無功功率，各PV節(jié)點的有功功率和節(jié)點電壓。而且由于其可交互性和可繪圖性，它也可以用于向非專業(yè)用戶解釋電力系統(tǒng)的運行操作。具體參數(shù)設(shè)置如圖3所示，機組的有功出力是不確定的，可在框內(nèi)暫時填為100MW。接下來，我們將兩種方法得出的各節(jié)點電壓和各支路的潮流狀態(tài)進行對比：（1）各節(jié)點電壓的對比如表9所示：表各節(jié)點電壓的對比節(jié)點編號標(biāo)幺電壓實際電壓(kV)相角 (度)PowerWorldMatlabPowerWorldMatlabPowerWorldMatlab10023456比較兩種計算結(jié)果的各節(jié)點電壓可以發(fā)現(xiàn)，%。則歸算到全網(wǎng)基準(zhǔn)后的發(fā)電機參數(shù)可按式10來計算： （10）代入?yún)?shù)，求得xd=，xq=，xd’=.為了簡化網(wǎng)絡(luò)，將三臺發(fā)電機并聯(lián)等值成一臺機組。由于穩(wěn)定時，P=，代入式13后求得功角初值。此外，為了加快功角的收斂速度，添加了阻尼部分，阻尼系數(shù)D設(shè)為10。%給定微分方程組的初值a=。b=。x2]。圖tc=（穩(wěn)定）圖2tc=（不穩(wěn)定）（2）采用Matlab求極限切除角通過求解故障時發(fā)電機轉(zhuǎn)子運動方程來確定功角隨時間變化的特性，進而可以求出極限切除角。b=。在主窗口中輸入：[x,y]=ginput(1)。（5）節(jié)點3發(fā)生三相對稱短路時，用Matlab軟件求出的極限切除時間tc=，用PowerWorld軟件求出的極限切除時間tc’=。 參數(shù)的求取線路L23中點處發(fā)生三相對稱短路，通過網(wǎng)絡(luò)變換法同理可得，轉(zhuǎn)移阻抗Z17==，輸入阻抗Z77= + =。w0]。%給第二個計算區(qū)間賦初值a=。x=[x1。=。綜上，對于故障方案2，=，=。臨界角為，滿足式，即（22）求得設(shè)極限切除角為，由等面積定則有（23）將參數(shù)代入后解方程得，極限切除角幾種方法求出來的結(jié)果統(tǒng)計如下：表1幾種方法求出來的結(jié)果統(tǒng)計計算方法MatlabPowerWorld等面積法（s）極限切除角（度）124從表14可知，采用等面積法求出的極限切除角是最小的，其次是Matlab仿真法，最大的是PowerWorld仿真法。圖43綜上，取tc。圖45（3）5s時線路L23中點處發(fā)生三相對稱短路。 不同發(fā)電機模型對暫定過渡時間的影響暫定過渡時間是指從發(fā)生故障時刻算起到功角穩(wěn)定于一個定值后的這一段時間，暫態(tài)過渡時間越短意味功角穩(wěn)定地越快，即系統(tǒng)的暫態(tài)穩(wěn)定性越好。計算結(jié)果如表16所示：表1TJ取不同值時系統(tǒng)的極限切除角TJ20406080100(s)(s)隨TJ變化特性如圖30所示，可以發(fā)現(xiàn)，當(dāng)TJ增大時，極限切除角也變大。此外，當(dāng)Tr較小時，極限切除角下降較為明顯；當(dāng)Tr，極限切除角基本不變。0 +i* 1 0 0 2。4 0。)。請輸入由支路參數(shù)形成的矩陣：B1=39。%輸入Xna=input(39。O=zeros(1,n)。end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。end %求導(dǎo)納矩陣G=real(Y)。f(i)=imag(B2(i,3))。 for k=i+1:n for j1=i+1:n B(k,j1)=B(k,j1)B(k,i)*B(i,j1)。 else k=i+1。K=1。DP(i)=DP1(i)./V(i)。 IC4=L1。for k=1:n C(i)=C(i)+V(k)*(G(i,k)*sin(O(i)O(k))BI(i,k)*cos(O(i)O(k)))。fori=1:n if B2(i,6)==2C(i)=0。fori=1:naDQ(i)=A(i,i)*DQ(i)。fori=1:nif B2(i,6)==2 L=L+1。迭代次數(shù)39。每次沒有達(dá)到精度要求的無功功率個數(shù)為39。disp(E)。disp(O)。disp(39。disp(Si(p,q))。 end Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。q=B1(i,1)。 plot(Cs,Dy),xlabel(