C語言進(jìn)行潮流計算

C語言進(jìn)行潮流計算

ID:45097100

大?。?46.00 KB

頁數(shù):16頁

時間:2019-11-09

C語言進(jìn)行潮流計算_第1頁
C語言進(jìn)行潮流計算_第2頁
C語言進(jìn)行潮流計算_第3頁
C語言進(jìn)行潮流計算_第4頁
C語言進(jìn)行潮流計算_第5頁
資源描述:

《C語言進(jìn)行潮流計算》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、電力系統(tǒng)課程設(shè)計C語言潮流計算學(xué)院:電氣工程班級:電092班學(xué)號:0912002020學(xué)生姓名:閔凱2013.3.7-16-電力系統(tǒng)的潮流計算是對電力系統(tǒng)分析的最基本步驟也是最重要的步驟,是指在一定的系統(tǒng)結(jié)構(gòu)和運(yùn)行條件下,確定系統(tǒng)運(yùn)行狀態(tài)的計算,也即是對各母線(節(jié)點(diǎn))電壓,各元件(支路)傳輸電線或功率的計算。通過計算出的節(jié)點(diǎn)電壓和功率分布用以檢查系統(tǒng)各元件是否過負(fù)荷,各點(diǎn)電壓是否合理,以及功率損耗等。即使對于一個簡單的電力系統(tǒng),潮流計算也不是一件簡單就可以完成的事,其運(yùn)算量很大,因此如果對于一個大的、復(fù)雜的電網(wǎng)來說的話,由于其節(jié)點(diǎn)多,分支雜,其計算量可想而知,

2、人工對其計算也更是難上加難了。特別是在現(xiàn)實(shí)生活中,遇到一個電力系統(tǒng)不會像我們期望的那樣可以知道它的首端電壓和首端功率或者是末端電壓和末端功率,而是只知道它的首端電壓和末端功率,更是使計算變的頭疼萬分。為了使計算變的簡單,我們就可以利用計算機(jī),用C語言編程來實(shí)現(xiàn)牛頓-拉夫遜(Newton-Raphson)迭代法,最終實(shí)現(xiàn)對電力系統(tǒng)潮流的計算。一.用牛頓-拉夫遜迭代法進(jìn)行電力系統(tǒng)潮流計算的相關(guān)概念1.節(jié)點(diǎn)導(dǎo)納矩陣如圖所示的電力網(wǎng)絡(luò),將節(jié)點(diǎn)i和j的電壓用i和j表示,它們之間的支路導(dǎo)納表示為yij,那么有基爾霍夫電流定律可知注入接點(diǎn)I的電流i(設(shè)流入節(jié)點(diǎn)的電流為正)等

3、于離開節(jié)點(diǎn)I的電流之和,因此有jI(1-1)(1-2)如令則可將(1-2)改寫為:-16-I=1,2,…,n.(1-3)上式也可以寫為:I=YU(1-4)其中Y為節(jié)點(diǎn)導(dǎo)納矩陣,也稱為稀疏的對稱矩陣,它是n×n階方陣。對角元Yii稱為自導(dǎo)納,它等于與該節(jié)點(diǎn)I直接相連的所有支路導(dǎo)納總和;非對角元Yij(i≠j)稱為互導(dǎo)納或轉(zhuǎn)移導(dǎo)納,它等于連結(jié)節(jié)點(diǎn)I,j支路導(dǎo)納的負(fù)數(shù),且有Yij=Yji,當(dāng)節(jié)點(diǎn)I,j之間沒有支路直接相連時,Yij=Yji=0。電力系統(tǒng)的分析計算中,往往要作不同的運(yùn)行方式下的潮流計算,如果系統(tǒng)發(fā)生變化,如投切一條線路或一臺變壓器,由于改變了一條支路的

4、狀態(tài)或參數(shù)只影響該支路兩端節(jié)點(diǎn)的自導(dǎo)納和他們之間的互導(dǎo)納,因而對每一種運(yùn)行方式不必重新形成導(dǎo)納矩陣,只需對原有導(dǎo)納矩陣作相應(yīng)的修改即可。2.潮流計算的功率方程在實(shí)際的電力系統(tǒng)中,已知的條件往往不是節(jié)點(diǎn)的注入電流而是負(fù)荷和發(fā)電機(jī)功率,且這些功率一般不隨節(jié)點(diǎn)的電壓變化而變化,而節(jié)點(diǎn)的電流則是隨電壓的變化而變化的,因此在已知節(jié)點(diǎn)導(dǎo)納矩陣的情況下,必須用已知的節(jié)點(diǎn)功率來替代未知的節(jié)點(diǎn)注入電流,才能求出節(jié)點(diǎn)電壓,每一個節(jié)點(diǎn)的注入功率方程式為:GI(1-5)(1-6)(1-7)節(jié)點(diǎn)注入電流用功率和電壓表示為:(1-8)功率方程可以表示為:(1-9)3.節(jié)點(diǎn)分類對于有n個節(jié)

5、點(diǎn)的電力網(wǎng)絡(luò),可以列出n個功率方程,由圖可知一個節(jié)點(diǎn)有四個變量:注入有功功率Pi,注入無功功率Qi,節(jié)點(diǎn)電壓幅值Ui和相角。n-16-個節(jié)點(diǎn)有4n個變量,但只有2n個關(guān)系式,所以為了使潮流有確定解,必須給定其中2n個變量。根據(jù)給定節(jié)點(diǎn)變量的不同,可以有以下三種類型的節(jié)點(diǎn):(1)PQ節(jié)點(diǎn):給定注入功率Pi,Qi,即已知PGi,PLi,QGi,QLi,待求Ui,δi。例如:降壓變電所母線(負(fù)荷節(jié)點(diǎn)),固定出力的發(fā)電廠母線。(2)PV節(jié)點(diǎn):給定了注入有功功率Pi(PGi,PLi),Ui和QLi,待求QGi(Qi),δi。例如:有一定無功電源的降壓變電所母線,有一定儲

6、備的發(fā)電廠母線。(3)平衡節(jié)點(diǎn):給定了Ui,δi和PLi,QLi,待求PGi,QGi,即Pi,Qi,用來平衡全電網(wǎng)的功率,通常在一個獨(dú)立的電力系統(tǒng)中只設(shè)一個平衡節(jié)點(diǎn)。4.牛頓-拉夫遜迭代法牛頓-拉夫遜迭代法將解非線性方程組的過程轉(zhuǎn)化為反復(fù)求與之相對應(yīng)的線性方程的求解過程。對于一個n維非線性方程組:n=1,2,3,…,n假定其初值為x1(0),x2(0),…,xn(0),也即其近似解,它與真值之間的誤差為也即各變量與真解之間的修正量。將這n個方程式都在初值的附近展開成Taylor級數(shù)且忽略二次項(xiàng)及高次項(xiàng),則可得修正方程,I=1,2,…n.(1-10)將修正方程寫

7、成矩陣形式:(1-11)其中令J=,稱之為雅可比(Jacobi)方陣。它的第I行,第j列交點(diǎn)的元素為第I個函數(shù)對第j個變量xj的偏導(dǎo)數(shù)在點(diǎn)(x1(0),x2(0),…,xn(0))的值,所以方程組是線性方程,可用于求出,從而得到新的近似解,(1-12)于是得到一般迭代式:-16-(1-13)于是得到近似解:(1-14)迭代一直進(jìn)行到Max{

8、yi-fi(x1(0),x2(0),…,xn(0))

9、}<ε或Max{

10、Δxi(k)

11、}<ε為止。一.用牛頓-拉夫遜迭代法進(jìn)行潮流計算設(shè)網(wǎng)絡(luò)中除參考節(jié)點(diǎn)外有n個節(jié)點(diǎn),其中1個平衡節(jié)點(diǎn)(并令第n個節(jié)點(diǎn)為平衡節(jié)點(diǎn)),m個PQ節(jié)

12、點(diǎn)(為第1~m個節(jié)點(diǎn)),有n-m-1個

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動畫的文件,查看預(yù)覽時可能會顯示錯亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負(fù)責(zé)整理代發(fā)布。如果您對本文檔版權(quán)有爭議請及時聯(lián)系客服。
3. 下載前請仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時可能由于網(wǎng)絡(luò)波動等原因無法下載或下載錯誤,付費(fèi)完成后未能成功下載的用戶請聯(lián)系客服處理。