【正文】 ],x02,h_opt,a,b,c)。b39。 h_opt=odeset。,[0,tc],x0,h_opt,a,b,c)。 x03=x2(end,:)39。 a=。t3。,t,x(:,2),39。 圖中曲線為發(fā)電機(jī)的功角變化曲線，直線為平衡機(jī)的功角變化曲線，由此可以看出發(fā)電機(jī)的初始功角為 度。因而，計算結(jié)果是正確的。由此可得出發(fā)電機(jī)采用三階模型時系統(tǒng)的極限切除時 間 =。 % 1 3 + 1 0 0。 % 0 + 0 1 0 2。 % 0 + 0 1 0 2。 % 4 5 + 1 0 0。 圖 26. 線路 24 中點 時單相故障， 時切除， 時重合閘， 時再次切除故障線路 由此可以看出，單相接地故障時發(fā)電機(jī)的功角擺動幅度較小，極限切除時間（ ）也較三相短路時（ ）長。 33 圖 18. 線路 24 發(fā)生 三相接地故障時各節(jié)點的電壓 圖 19. 線路 24 發(fā)生單相接地故障時各節(jié)點的電壓 采用 Power World 仿真軟件分析故障（三階） 三相短路時的故障分析 在 Power World 軟件里給發(fā)電機(jī)建立暫態(tài)穩(wěn)定模型，其中平衡節(jié)點用一臺大發(fā)電機(jī)代替，選定 GENPWTwoAxis 模型，把其中的 H 值設(shè)得非常大（如 ），其他都用默認(rèn)參數(shù)。 圖 15. 線路 發(fā)生三相短路故障， 切除時發(fā)電機(jī)的暫態(tài)圖像 3）用試探法找出發(fā)電機(jī)的極限切除時間， 線路 24 中點發(fā)生三相對稱短路故障， 切除故障（圖 16 左）時系統(tǒng)穩(wěn)定， 切除故障（圖 16 右）時系統(tǒng)不穩(wěn)定，如圖 。) 由以上程序仿真發(fā)電機(jī)搖擺曲線的結(jié)果如圖 10 所示。 x=[x1。c=()*pi/180。b=。 a=。w0]。 0 5 10 15 20 25 30 35 40 0 . 4 0 . 200 . 20 . 40 . 60 . 811 . 21 . 41 . 6 圖 8. 切除 24 三相短路故障時發(fā)電機(jī)的功角 （紅） 、角速度 （藍(lán)） 變化曲線 當(dāng)故障切除時間改為 tc=（左）和 tc=（右），修改相應(yīng)的故障切除時間 tc，得到的發(fā)電機(jī)暫態(tài)搖擺曲線如圖 。t2]。%求解第一個區(qū)間 [0,tc] x02= x1 (end,:)39。 x0=[*pi/180。慣性時間 TJN=11，歸算到全網(wǎng)后 300= 1 1 = 3 8 . 8 2 40 . 8 5 1 0 0GNJ J N BSTT S?? ?，三臺機(jī)組等效為一臺后TJΣ= 3=。時電磁功率最大， PE39。 穩(wěn)態(tài)無故障 +(176。 11PR Q XV V??? （ 11） 11PX Q RV V? ?? （ 12） 2239。φ， φ為輸入阻抗或轉(zhuǎn)移阻抗的阻抗角。處理后的網(wǎng)架圖如圖 7 所示。39。因而可以使用一個純電抗等效三臺變壓器，其中電抗值 2 2()* 22()% 1 0 . 5 1 0 . 5 1 0 0 0 . 0 3 6 31 0 0 1 0 0 3 5 0 9 . 5 4 5 5TNS BT T N BVV SX SV? ? ? ? ? ? ? （ 10） 三臺變壓器等效為一臺變壓器時的標(biāo)幺阻抗為 XT*=。等效后的出力為837MW， Xd*=， Xq*=， Xd39。兩種計算潮流得出的支路潮流具有很好的一致性。 表 8. 各支路的狀態(tài) 線路名稱 首端有功功率（ MW） 首端無功功率（ Mvar） 末端有功功率（ MVA） 末端無功功率（ Mvar） 有功損耗（ MW） 無功損耗（ Mvar） L12 L13 L24 L45 0 L46 L37 L57 L67 注： Lij 表示 i 為線路首端， j 為線路末端 4）電網(wǎng)的網(wǎng)損率計算 電網(wǎng)的網(wǎng)損為發(fā)電機(jī)發(fā)出的功率與平衡節(jié)點 1與外界的交換功率的之和減去負(fù)荷功率之和。其余各節(jié)點只設(shè)置基準(zhǔn)電壓 VB=220kV 即可。 與其它商業(yè)潮流計算軟件包不同，該軟件可以讓用戶通過生動詳細(xì)的全景圖來觀察電力系統(tǒng)。π，表 的 功角是角度的形式表示的。計算不平衡功率 ΔQ 判斷 ΔQ 是否小于可接受誤差。本題形成的 B1 矩陣見附錄。節(jié)點 6 均為 PQ 節(jié)點，并且?guī)в幸欢ǖ呢?fù)荷。 例如，求解如下微分方程： 12212sin( )xx?? 初值為 x1=1,x2=0 從 0s積分到 3s，步長： 則，首先定義函數(shù) myfunc，計算微分方程右端項的值： function dx=myfunc(x) dx=[x(2) 2*sin(x(1) ] 。 編制課程設(shè)計報告。 3. 考慮該電廠開機(jī)三臺，均為額定容量輸出的方式，用編制的程序計算潮流。電廠 220kV 側(cè)采用單母分段接線，正常運行時分段開關(guān)閉合。圖表整潔，布局合理，文字注釋必須使用工程字書寫，不準(zhǔn)用徒手畫 3）畢業(yè)論文須用 A4 單面打印，論文 50 頁以上的雙面打印 4）圖表應(yīng)繪制于無格子的頁面上 5）軟件工程類課題應(yīng)有程序清單，并提供電子文檔 1）設(shè)計（論文） 2）附件：按照任務(wù)書、開題報告、外文譯文、譯文原文（復(fù)印件）次序裝訂 3）其它 5 目 錄 課 程 設(shè)計題 目 C： 電 力系 統(tǒng)穩(wěn) 定分析和 計 算 ............................................................................... 6 算 .......................................................................................................................... 11 線 路 參數(shù) 的 計 算和 標(biāo) 幺化 ................................................................................................ 11 節(jié) 點信息 統(tǒng)計 .................................................................................................................. 14 潮流 計 算 .......................................................................................................................... 14 采用 Matlab 計 算 電網(wǎng) 潮流 .............................................................................................. 14 網(wǎng)絡(luò) 信息 處 理 ................................................................................................................................................... 14 基于 Matlab 的 PQ 分解法程序 ................................................................................................................. 15 PQ 分解法潮流 計 算的 結(jié) 果分析 ............................................................................................................. 15 采用 Power World 計 算潮流 ............................................................................................ 16 Power world 仿 真軟 件 簡 介 ........................................................................................................................... 16 基于 Power World 軟 件 構(gòu) 建潮流模型 ..................................................................................................... 17 采用 power world 仿 真軟 件的潮流 計 算 結(jié) 果 ........................................................................................ 18 Matlab 仿 真計 算 結(jié) 果 與 Power World 仿 真軟 件 計 算 結(jié) 果 對 比 ........................................ 20 4. 轉(zhuǎn) 移阻抗的求解 ..................................................................................................................... 21 負(fù) 荷、 發(fā)電 機(jī)、 變壓 器、 線 路的 處 理 .............................................................................. 21 負(fù) 荷的 處 理 ....................................................................................................................................................... 21 發(fā)電 機(jī)的 處 理 ................................................................................................................................................... 21 變壓 器的 處 理 ................................................................................................................................................... 21 輸電線 路的 處 理 .............................................................................................................................................. 22 網(wǎng)絡(luò)變換 （ 穩(wěn)態(tài)時 ） ....................................................................................................... 23 消去 節(jié) 點 6 ................................................................................................................................................. 23 消去 節(jié) 點 7 ................................................................................................................................... 24 各 種運 行方式下 轉(zhuǎn) 移阻抗的求解 ..................................................................................... 24 功角特性 計 算 ....................................................................................