電力系統(tǒng)分析課程設(shè)計---電力系統(tǒng)短路故障的計算機算法程序設(shè)計

1、<p>　　電力系統(tǒng)分析課程設(shè)計報告書</p><p>　　題目: 電力系統(tǒng)短路故障的計算機算法程序設(shè)計 </p><p>　　專 業(yè)：電氣工程及其自動化 </p><p>　　班 級： </p><p>　　學 號： </

2、p><p>　　學生姓名： </p><p>　　指導(dǎo)教師： </p><p>　　2013年 3月 14日</p><p>　　信息工程學院課程設(shè)計任務(wù)書</p><p><b>　　目錄</b></p><p&

3、gt;　　1.任務(wù)提出與方案論證1</p><p>　　1.1計算機計算的由來和要求1</p><p>　　1.2論證方案的可行性1</p><p><b>　　2.總體設(shè)計3</b></p><p>　　2.1程序語言及論述3</p><p>　　2.2輔助知識及參數(shù)4</p&

4、gt;<p><b>　　3.詳細設(shè)計5</b></p><p>　　3.1程序流程圖及習題5</p><p>　　3.2程序及解題7</p><p><b>　　4.總結(jié)12</b></p><p><b>　　參考資料13</b></p>

5、<p>　　1.任務(wù)提出與方案論證</p><p>　　1.1計算機計算的由來和要求</p><p>　　隨著電子計算機的發(fā)展與完善，運算能力已日新月異，人工運算的時代慢慢隱形化，</p><p>　　程序化計算設(shè)計必將成為主體。為此，運用計算機計算短路電流已成為可能，也會帶來很大的方便。</p><p>　　依照所給的電力系統(tǒng)

6、，編制短路電流計算程序，通過計算機進行調(diào)試，最后完成一個切</p><p>　　實可行的電力系統(tǒng)計算應(yīng)用程序。通過自己設(shè)計電力系統(tǒng)計算程序使同學對電力系統(tǒng)分析有</p><p>　　進一步理解，同時加強計算機實際應(yīng)用能力的訓(xùn)練。</p><p>　　所謂短路，是指電力系統(tǒng)正常運行情況以外的相與相之間或相與地之間的“短接”。在電力系統(tǒng)正常運行時，除中性點外，相與相或相

7、與地之間是絕緣的。如果由于某種原因使其絕緣破壞而構(gòu)成了通路，我們就稱電力系統(tǒng)是發(fā)生了短路故障。在三相系統(tǒng)中，短路故障可分為兩大類：即對稱短路（三相短路）和不對稱短路（兩相短路、兩相接地短路、單相接地短路）。其中三相短路雖然發(fā)生的機會較少，但情況嚴重，又是研究其它短路的基礎(chǔ)。所以我們先研究最簡單的三相短路電流的暫態(tài)變化規(guī)律。</p><p>　　1.2論證方案的可行性</p><p>　　在

8、大量實際計算中，短路計算有著比較固定的方式，這就為計算機計算墊定了基礎(chǔ)，根據(jù)電力系統(tǒng)三相對稱短路的特點，建立了合理的三相短路的數(shù)學模型，在此基礎(chǔ)上， 形成電力系統(tǒng)短路電流實用汁算方法；節(jié)點阻抗矩陣的支路追加法．編制了對任意一個電力系統(tǒng)在任意點發(fā)生短路故障時三相短路電流及其分布的通用計算程序。 該辦法適用予各種復(fù)雜結(jié)構(gòu)的電力系統(tǒng)．從一個側(cè)面展示了計算機應(yīng)用于電力系統(tǒng)的廣闊前景．</p><p><b>　

9、　其原理如下：</b></p><p>　　1.節(jié)點導(dǎo)納矩陣計算電壓初始值的原理</p><p>　　在電力系統(tǒng)分析的第四章中學到了矩陣方程:YV=I,故由矩陣性質(zhì)得：V=ZI</p><p>　　2.節(jié)點阻抗矩陣計算短路電流的原理</p><p>　　利用節(jié)點阻抗矩陣計算短路電流如圖所示，假定系統(tǒng)中的節(jié)點f經(jīng)過過渡阻抗發(fā)生短路。

10、對于正常狀態(tài)的網(wǎng)絡(luò)而言，發(fā)生短路相當于在故障節(jié)點f增加了一個注入電流。因此，網(wǎng)絡(luò)中任一節(jié)點i的電壓可表示為：</p><p><b>　　公式一 ：</b></p><p>　　由式可見，任一節(jié)點i的電壓都由兩項疊加而成。第一項是當時由網(wǎng)絡(luò)內(nèi)所有電源在節(jié)點i產(chǎn)生的電壓，也就是短路前瞬間正常運行狀態(tài)下的節(jié)點電壓，記為 。第二項是當網(wǎng)絡(luò)中所有電流源都斷開，電勢源都

11、短接時，僅僅由短路電流 在節(jié)點i產(chǎn)生的電壓。這兩個分量的疊加，就等于發(fā)生短路后節(jié)點i的實際電壓，即</p><p>　　公式二也適用于故障節(jié)點f，于是有：</p><p>　　是故障節(jié)點f的自阻抗，也稱輸入阻抗。 根據(jù)故障的邊界條件：</p><p>　　由以上兩個方程式解出：</p><p><b>　　2.總體

12、設(shè)計</b></p><p>　　2.1程序語言及論述</p><p>　　MATLAB是一門用途廣泛的語言，它提供了專業(yè)水平的符號計算，文字處理，可視化建模仿真和實時控制等功能。</p><p>　　MATLAB的基本數(shù)據(jù)單位是矩陣，它的指令表達式與數(shù)學、工程中常用的形式十分相似,故用MATLAB來解決問題要比用C,FORTRAN等語言方便的多。在新的

13、版本中也加入了對C,FORTRAN,c++ ,JAVA的支持，可以直接調(diào)用,用戶也可以將自己編寫的實用程序?qū)氲組ATLAB函數(shù)庫中方便自己以后調(diào)用。</p><p>　　MATLAB的程序極其簡短且可移植性很好，基本上不做修改就可以在各種型號的計算機和操作系統(tǒng)上運行。MATLAB既具有結(jié)構(gòu)化的控制語句（如for循環(huán)，while循環(huán)，break語句和if語句），又有面向?qū)ο缶幊痰奶匦?。但較之高級語言，由于其程序不

14、用編譯等預(yù)處理，也不生成可執(zhí)行文件，而是邊解釋邊執(zhí)行，所以速度較慢。</p><p>　　綜合起來看，MATLAB是一個高度集成的系統(tǒng)，集科學計算、圖象處理、聲音處理于一體，具有極高的編程效率。近年來，MATLAB已經(jīng)從最初的“矩陣實驗室”，滲透到科學與工程計算的多個領(lǐng)域，在自動控制、信號處理、神經(jīng)網(wǎng)絡(luò)、模糊邏輯、小波分析等多個方向，都有著廣泛的應(yīng)用。因此，我選擇的計算機語言是MATLAB。</p>

15、<p>　　2.2輔助知識及參數(shù)</p><p>　　計算時選定一個基準容量(Sjz)和基準電壓(Ujz).將短路計算中各個參數(shù)都轉(zhuǎn)化為和該參數(shù)的基準量的比值(相對于基準量的比值),稱為標么值(這是短路電流計算最特別的地方,目的是要簡化計算). </p><p><b>　　(1)基準 </b></p><p>　　基準容量 Sjz

16、 =100 MVA </p><p>　　基準電壓 UJZ規(guī)定為8級. 230, 115, 37, 10.5, 6.3, 3.15 ,0.4, 0.23 KV </p><p>　　有了以上兩項,各級電壓的基準電流即可計算出,例: UJZ (KV)3710.56.30.4 </p><p>　　因為 S=1.73*U*I 所以 IJZ (KA)1.565.59.16

17、144 </p><p><b>　　(2)標么值計算 </b></p><p>　　容量標么值 S* =S/SJZ.例如:當10KV母線上短路容量為200 MVA時,其標么值容量 </p><p>　　S* = 200/100=2. </p><p>　　電壓標么值 U*= U/UJZ ; 電流標么值 I* =I/IJ

18、Z </p><p>　　有了標幺值的概念后，短路電流計算的過程就大大減少了，簡化了運算。</p><p><b>　　2.3 等值電路</b></p><p>　　計算前先要做出等值電路，并計算等值參數(shù)。在詳細設(shè)計中會將習題6-10及習題6-11的等值電路及參數(shù)列寫出來，這里只做簡單說明：Sd三相短路容量 (MVA)簡稱短路容量校核開關(guān)分斷容

19、量 Id三相短路電流周期分量有效值(KA)簡稱短路電流校核開關(guān)分斷電流 和熱穩(wěn)定 IC三相短路第一周期全電流有效值(KA) 簡稱沖擊電流有效值校核動穩(wěn)定 ic三相短路第一周期全電流峰值(KA) 簡稱沖擊電流峰值校核動穩(wěn)定 x電抗(Ω) 其中系統(tǒng)短路容量Sd和計算點電抗x 是關(guān)鍵。</p><p><b>　　3.詳細設(shè)計</b></p><p>　　3.1程序流程圖及

20、習題</p><p><b>　　一 流程圖</b></p><p><b>　　二 習題</b></p><p>　　1習題6-10 （題目見書）</p><p><b>　　等值網(wǎng)絡(luò)及參數(shù)：</b></p><p>　　選SB=120MV.A,VB=

21、Vav,各元件標幺值參數(shù)如下：</p><p>　　X10=0.105 X20=0.28 X24=0.21 X43=0.43554 X35=0.29036 X45=0.25406</p><p>　　Y340=Y430=j0.01852 Y350=Y530=j0.01234 Y450=Y540=j0.01080</p><p><b>　

22、　故節(jié)點導(dǎo)納為：</b></p><p>　　Y11=-j13.8716 Y13=Y31=j9.5238 Y22=-j8.33343 Y24=Y42=j4.7619 Y33=-j15.2329 </p><p>　　Y34=Y43=j2.2960 Y35=Y53=j3.4440 Y44=-j10.9646 Y45=Y54=j3.9360 Y55=-j7

23、.3569</p><p>　　2習題6-11（題目見書）</p><p><b>　　等值網(wǎng)絡(luò)及參數(shù)：</b></p><p><b>　　節(jié)點導(dǎo)納為：</b></p><p>　　Y11=-j10 Y12=Y21=j5 Y22=-j14.5 Y23=Y32=j10 Y33=-j17

24、.6667</p><p><b>　　3.2程序及解題</b></p><p><b>　　程序如下：</b></p><p><b>　　1.通用程序為：</b></p><p>　　function dl63</p><p><b>　

25、　text63;</b></p><p><b>　　Z=inv(Y)</b></p><p>　　m=size(Z,1)</p><p>　　I=I' </p><p><b>　　V=Z*I </b></p><

26、p>　　i=input(‘請輸入節(jié)點i=')</p><p>　　Ii=1/Z(i,i)</p><p><b>　　for k=1:m</b></p><p><b>　　k</b></p><p>　　Vk=V(k,1)-Z(k,i)*Ii</p><p>

27、<b>　　end</b></p><p><b>　　for p=1:m</b></p><p>　　for q=(1+p):m</p><p>　　if(Y(p,q)~=0)</p><p><b>　　p,q</b></p><p>　　Vp=V(

28、p,1)-Z(p,i)*Ii;</p><p>　　Vq=V(q,1)-Z(q,i)*Ii;</p><p>　　Ipq=-(Vp-Vq)*Y(p,q)</p><p><b>　　end</b></p><p><b>　　end </b></p><p><b>

29、;　　end</b></p><p>　　1習題6-10 解：</p><p><b>　　輸入數(shù)據(jù)為：</b></p><p>　　Y=[-13.8716j 0 9.5238j 0 0 </p><p>　　0 -8.3333j

30、 0 4.7619j 0</p><p>　　9.5238j 0 -15.2329j 2.2960j 3.4440j</p><p>　　0 4.7619j 2.2960j -10.9646j 3.9360j</p><p>　　0 0

31、3.4440j 3.9360j -7.3569j]</p><p>　　I=[4.5652j 3.75j 0 0 0]</p><p><b>　　m=5</b></p><p><b>　　程序為：</b></p><p>　　function dl610

32、 %定義一個函數(shù)名為dl610的函數(shù)%</p><p>　　text610; %調(diào)用函數(shù)名為text610的函數(shù)%</p><p>　　Z=inv(Y) %求節(jié)點阻抗矩陣%</p><p>　　I=I'

33、 %求I的轉(zhuǎn)置矩陣%</p><p>　　V=Z*I %求各節(jié)點的初始電壓%</p><p>　　SB=120; %選定基準功率為120MV.A%</p><p>　　VB=115; %選定基準電壓為115KV%&l

34、t;/p><p>　　IB=SB/(VB*sqrt(3)) %求基準電流%</p><p>　　i=input('ÇëÊäÈë¶Ì·½Úµãi=') %輸入短路節(jié)點i%</p><p>　　I

35、i=V(i,1)/Z(i,i); %求短路節(jié)點電流%</p><p>　　Ii1=Ii*IB </p><p><b>　　for p=1:m</b></p><p>　　for q=(1+p):m</p><p>　　if(Y(p,q)~=0)

36、 %如果p,q節(jié)點間的導(dǎo)納不為%</p><p>　　Vp=V(p,1)-Z(p,i)*Ii; %分別求p,q的節(jié)點電壓%</p><p>　　Vq=V(q,1)-Z(q,i)*Ii;</p><p>　　Ipq=-(Vp-Vq)*Y(p,q); %求p,q節(jié)點間的支路電流%</p><p>　　Ipq1=

37、Ipq*IB; %求短路電流的實際值 %</p><p>　　if(p==3&&q==5) %求支路L-2的短路電流%</p><p><b>　　p</b></p><p><b>　　q</b></p><p><b>　

38、　Ipq</b></p><p><b>　　Ipq1</b></p><p><b>　　end</b></p><p>　　if(p==4&&q==5) %求支路L-3的短路電流%</p><p><b>　　p</b></p&g

39、t;<p><b>　　q</b></p><p><b>　　Ipq</b></p><p><b>　　Ipq1</b></p><p><b>　　end</b></p><p><b>　　end</b><

40、/p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　輸出結(jié)果為：</b></p><p><b>　　Y =</b></p><p>　　0 -13.8716i

41、0 0 + 9.5238i 0 0 </p><p>　　0 0 - 8.3333i 0 0 + 4.7619i 0 </p><p>　　0 + 9.5238i 0

42、 0 -15.2329i 0 + 2.2960i 0 + 3.4440i</p><p>　　0 0 + 4.7619i 0 + 2.2960i 0 -10.9646i 0 + 3.9360i</p><p>　　0 0

43、 0 + 3.4440i 0 + 3.9360i 0 - 7.3569i</p><p><b>　　I =</b></p><p>　　0 + 4.5652i 0 + 3.7500i 0 0 0 </p>

44、;<p><b>　　m =</b></p><p><b>　　5</b></p><p><b>　　Z =</b></p><p>　　0 + 0.1819i 0 + 0.0616i 0 + 0.1600i 0 + 0.1079i

45、 0 + 0.1326i</p><p>　　0 + 0.0616i 0 + 0.2078i 0 + 0.0898i 0 + 0.1536i 0 + 0.1242i</p><p>　　0 + 0.1600i 0 + 0.0898i 0 + 0.2330i 0 + 0.1571i

46、 0 + 0.1932i</p><p>　　0 + 0.1079i 0 + 0.1536i 0 + 0.1571i 0 + 0.2689i 0 + 0.2174i</p><p>　　0 + 0.1326i 0 + 0.1242i 0 + 0.1932i 0 + 0.2174i

47、 0 + 0.3427i</p><p><b>　　I =</b></p><p>　　0 - 4.5652i</p><p>　　0 - 3.7500i</p><p>　　0 </p><p>　　0 </p><p>　　0

48、 </p><p><b>　　V =</b></p><p><b>　　1.0617</b></p><p><b>　　1.0607</b></p><p><b>　　1.0671</b></p><p>&

49、lt;b>　　1.0686</b></p><p><b>　　1.0713</b></p><p><b>　　IB =</b></p><p><b>　　0.6025</b></p><p>　　請輸入短路節(jié)點i=5</p><p&g

50、t;<b>　　i =</b></p><p><b>　　5</b></p><p><b>　　Ii1 =</b></p><p>　　0 - 1.8835i</p><p><b>　　p =</b></p><p><

51、b>　　3</b></p><p><b>　　q =</b></p><p><b>　　5</b></p><p><b>　　Ipq =</b></p><p>　　0 - 1.5954i</p><p><b>　　I

52、pq1 =</b></p><p>　　0 - 0.9611i</p><p><b>　　p =</b></p><p><b>　　4</b></p><p><b>　　q =</b></p><p><b>　　5</

53、b></p><p><b>　　Ipq =</b></p><p>　　0 - 1.5310i</p><p><b>　　Ipq1 =</b></p><p>　　0 - 0.9223i</p><p><b>　　2習題6-11解：</b>&

54、lt;/p><p><b>　　輸入數(shù)據(jù)：</b></p><p>　　Y=[-10.0j 5.0j 0.0</p><p>　　5.0j -14.5j 10.0j</p><p>　　0.0 10.0j -17.667j]</p><p>　　I=[5.25j 0

55、 7j]</p><p><b>　　程序為：</b></p><p>　　function dl611 %定義一個函數(shù)名為dl610的函數(shù)%</p><p>　　text611; %調(diào)用函數(shù)名為text610的函數(shù)%</p><p>　　Z=

56、inv(Y) %求節(jié)點導(dǎo)納矩陣%</p><p>　　I=I' %求I的轉(zhuǎn)置矩陣%</p><p>　　V=Z*I %求節(jié)點電壓的初始值%</p><p>　　i=input(請輸入節(jié)點i=')

57、 %輸入短路節(jié)點%</p><p>　　Ii=V(i,1)/Z(i,i) %計算短路電流%</p><p><b>　　輸出結(jié)果為：</b></p><p><b>　　Y =</b></p><p>　　0 -10.0000i 0 + 5.0000i

58、 0 </p><p>　　0 + 5.0000i 0 -14.5000i 0 +10.0000i</p><p>　　0 0 +10.0000i 0 -17.6670i</p><p><b>　　I =</b></p>&

59、lt;p>　　0 + 5.2500i 0 0 + 7.0000i</p><p><b>　　Z =</b></p><p>　　0 + 0.1394i 0 + 0.0789i 0 + 0.0446i</p><p>　　0 + 0.0789i 0

60、+ 0.1577i 0 + 0.0893i</p><p>　　0 + 0.0446i 0 + 0.0893i 0 + 0.1071i</p><p><b>　　I =</b></p><p>　　0 - 5.2500i</p><p>　　0 </p&

61、gt;<p>　　0 - 7.0000i</p><p><b>　　V =</b></p><p><b>　　1.0445</b></p><p><b>　　1.0390</b></p><p><b>　　0.9843</b><

62、/p><p>　　請輸入短路節(jié)點i=3</p><p><b>　　i =</b></p><p><b>　　3</b></p><p><b>　　Ii =</b></p><p>　　0 - 9.1875i</p><p>&l

63、t;b>　　4.總結(jié)</b></p><p>　　在這次課程設(shè)計當中，在學到一定東西知識的同時，也暴露了更多以前不曾注意的問題。例如這次做短路故障分析，看書覺得都懂了，等到真正開始做的時候才發(fā)現(xiàn)很多學的都不是很透徹，這也算是一次查漏補缺的機會，把以前學過的知識再鞏固一遍。</p><p>　　這次編程選用MATLAB語言，因為之前只是了解，掌握和運用的并不是很好，所以再做

64、的時候出現(xiàn)了很多問題，現(xiàn)在還有很多疑問不清楚，深深的理解到要學會一種軟件是很必要的，我也在繼續(xù)努力，雖然本課程設(shè)計已告于此，但是對MATLAB的學習還只是個開始，我還要進一步去學習，做到真正的學會學活。總之，再這次課程設(shè)計中，我清楚的意識到理論與實踐有很大的差別，理論只有與實踐相結(jié)合才能更好的發(fā)揮理論的作用，只有在實踐中才能找到實用的東西，僅有理論是遠遠不夠的。</p><p><b>　　參考資料

65、 </b></p><p>　　1.電力系統(tǒng)分析（上、下） 何仰贊 華中科技大學出版社</p><p>　　2.MATLAB電子仿真與應(yīng)用教程 王華 國防工業(yè)出版社</p><p>　　3.孫麗華 電力工程基礎(chǔ) 機械工業(yè)出版社</p><p>　　4.傅知蘭 電力系統(tǒng)電氣設(shè)備選擇與實用計算 中國電力出版