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