電力系統(tǒng)暫態(tài)上機計算課程設計報告

上傳人：奔*** IP屬地：河北更新時間：2024-03-01 格式：doc 頁數(shù)：24 大?。?57.50KB 人氣指數(shù)：12 舉報 版權申訴

已閱讀1頁，還剩23頁未讀，繼續(xù)免費閱讀

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、　　課程設計報告　　( 2010—2011年度第二學期)　　名稱：電力系統(tǒng)暫態(tài)上機 　　院系：電氣與電子工程學院 　　班級：電氣0803班 </p&g

2、t;　　日期： 2011年 7月6　　一、課程設計的目的與要求　　鞏固電力系統(tǒng)暫態(tài)分析理論知識，使學生掌握采用計算機對電力系統(tǒng)電磁暫態(tài)過程和機電暫態(tài)過程進行計算的基本方法，并進一步鞏固計算機編程能力，為將來從事相關的技術工作打下必要的基礎。　　二、主要內容</b&g

3、t;　　1．電磁暫態(tài)過程計算　　1.1 三相短路計算　　1.1.1 計算機計算三相短路流程圖（程序見附錄）　　附：式3-17 If = Uf|0| / ( Zff + zf )≈1/( Zff + zf ) ≈1/Zff　　式3

4、-19 U1= U1|0| +△U1 = U1|0| – Z1fIf ≈1 – Z1fIf　　Uf = Uf|0| +△Uf = ZfIf　　Un= Un|0|+△Un = Un|0| – ZnfIf ≈1 – ZnfIf　　式3-20 Iij = (Ui – Uj) / zij ≈(△Ui –△Uj )

5、/ zij = (△Ui –△Uj) yij　　1.1.2 計算教材《電力系統(tǒng)暫態(tài)分析(第三版)》P68例（3-2）　　例3-2：下圖所示一環(huán)形網(wǎng)絡，已知各元件參數(shù)為：　　發(fā)電機：G1～G3 100MW,10.5KV,cosφN=0.86，xd”=0.183。　　變壓器

6、：T1～T3 120MVA,115/10.5KV,Us(%)=10.5。　　線路：三條線路完全相同，長50km,電抗0.44歐姆/㎞　　試計算母線③三相短路后時刻的各節(jié)點電壓、各支路電流以及發(fā)電機端電壓。　　電路圖和等值電路圖如下：　　故障分量網(wǎng)絡：&l

7、t;/b>　　1.1.3節(jié)點發(fā)生三相短路接地計算結果得：（下面各數(shù)值均為標幺值）　　各節(jié)點的三相電壓　　各支路電流　　發(fā)電機的端電壓

8、　　1.2.不對稱短路計算過程流程圖和計算結果　　1.2.1不對稱短路計算過程流程圖：　　1.2.2計算教材《電力系統(tǒng)暫態(tài)分析(第三版)》P127例（5-1）　　例5-1：在例3-2的系統(tǒng)中又已知三臺發(fā)電機中性點均不接地；三臺變壓器均為YNd11接線（發(fā)電機側為三角形）；經試驗得三條輸電線路

9、的零序電抗均為0.20（以60MVA為基準值）。　　要求計算節(jié)點③分別發(fā)生a相接地短路、a相經10歐姆的電阻接地短路、bc兩相相間短路、bc兩相接地短路時各節(jié)點的電壓、各支路電流以及發(fā)電機端電壓。　　三個序網(wǎng)圖：　　正序：</p&g

10、t;　　負序： 　　零序：　　3.2.1 計算結果（標幺值） 　?。?）A相接地短路　　各節(jié)點電壓：</b&g

11、t;　　各支路電流：　　發(fā)電機端電壓　?。?）A相經10歐姆的電阻接地短路　　各節(jié)點電壓：　　各支

12、路電流　　發(fā)電機端電壓　?。?）BC兩相相間短路　　各節(jié)點電壓　　各支路電流

13、　發(fā)電機端電壓　?。?）BC兩相接地短路　　各節(jié)點電壓　　各支路電流　　發(fā)電機端電壓<b&

14、gt;　　2．靜態(tài)穩(wěn)定計算　　2.1靜態(tài)穩(wěn)定性計算過程流程圖：　　2.2 計算教材《電力系統(tǒng)暫態(tài)分析(第三版)》P191例（7-2）中，采用試探法選擇放大倍數(shù)　　例7-2：一簡單電力系統(tǒng)參數(shù)為：xd=xq=0.982,xd’=0.344,xe=0.504,xd∑=1.486, xd∑’=0

15、.48。已知正常運行時PE0=1，Eq0=1.972, σ0=49°,U=1。若發(fā)電機裝有按端電壓偏差調節(jié)的比例式勵磁調節(jié)器，采用試探法求不同勵磁放大倍數(shù)情況下，系統(tǒng)的靜態(tài)穩(wěn)定性。　　2.3 計算結果：　　Ke delta P　　將三者關系兩兩繪成圖線：</

16、p>　　用試探法可知，隨著勵磁放大倍數(shù)Ke的增大，功角delta先增大再減小，功率極限P也是先增大再減小。分析數(shù)據(jù)可知：當Ke=1.8時，系統(tǒng)可達靜態(tài)穩(wěn)定極限，即有功角deltasl≈105°、P=1.613924。　　3.暫態(tài)穩(wěn)定計算　　3.1暫態(tài)穩(wěn)定性計算過程流程圖：<

17、/p>　　3.2 已知：系統(tǒng)的額定頻率f=50，歸算后的發(fā)電機慣性時間常數(shù)Tj=8.47，正常運行時發(fā)電機向無窮大系統(tǒng)傳輸?shù)挠泄β蔖T=1，發(fā)生故障時發(fā)電機的最大功率P2M=0.48，故障切除后發(fā)電機的最大功率 P3M=1.38。用改進歐拉法編程，分析切除時間為0.15s和0.25s兩種情況下?lián)u擺曲線以及分析系統(tǒng)的暫態(tài)穩(wěn)定性。　　3.3 計算

18、結果　　對于不同的切除時間其結果是不一樣的。　?。?）切除時間為0.15秒時搖擺曲線：　　可以看出由于及時切除故障，系統(tǒng)會保持暫態(tài)穩(wěn)定。　　（2）切除時間為0.25秒時搖擺曲線：　　可以看出由于切除故障時間

19、過晚，系統(tǒng)無法保持暫態(tài)穩(wěn)定。　　4．思考題　　4.1計算短路電流，書中給出的手算方法與計算機編程方法有何區(qū)別？　　答：手算方法不用形成導納矩陣，通過星角變換化簡，求得最后的等值阻抗，但是如果要求非故障點的電壓或電流量還得倒回去求。計算機編程方法要形成導納矩陣，主要是利用公式來

20、求故障點和非故障點的量。　　4.2用計算機方法，進行電力系統(tǒng)潮流計算和短路電流計算，有哪些區(qū)別？　　答：電力系統(tǒng)常規(guī)的潮流計算的任務是根據(jù)給定的電網(wǎng)結構，發(fā)電規(guī)劃及負荷分布情況，求出整個電網(wǎng)的運行狀態(tài)，如各母線上的電壓（幅值及相角），網(wǎng)絡中的功率分布以及功率損耗等，潮流計算的結果是電力系統(tǒng)穩(wěn)定計算和故障分析的基礎，是正常時候，系統(tǒng)里功率和電流的分布。而短

21、路電流計算則是在系統(tǒng)發(fā)生故障時計算短路后系統(tǒng)里各點電壓，各支路流過的電流。它要先算出正常時電壓和電流，即潮流分布。再算故障分量，兩者相加才是最后系統(tǒng)里各點電壓，流過的電流。可以說潮流計算是短路電流計算的一部分。　　4.3如果交給你一個任務，請你用已有的短路電流計算軟件計算某地區(qū)電網(wǎng)某點發(fā)生短路后的短路電流，你應該做哪些工作？（收集哪些數(shù)據(jù)，如何整理，計算結果如何分析）&l

22、t;p>　　答：首先要收集該地區(qū)電網(wǎng)的電氣連接圖，收集各元件的參數(shù)，有發(fā)電機、變壓器和線路的正序、負序和零序的阻抗參數(shù)。將電氣連接圖畫出電網(wǎng)的等值電路的圖，將有名值參數(shù)裝換成在同一基準下的標幺值。利用保護裝置可以檢測出是哪里出現(xiàn)何種故障，根據(jù)故障點位置以及故障類型可以計算出短路電流?？炊搪冯娏髦档拇笮?，如果較小，電網(wǎng)可以持續(xù)運行1～2個小時，但應發(fā)出報警信號，讓工作人員去把故障切除；如果短路電流較小，應迅速切除故障，可以采用自動重

23、合閘來減低故障帶來的危害。　　4.4靜態(tài)穩(wěn)定的概念，暫態(tài)穩(wěn)定的概念？　　答：靜態(tài)穩(wěn)定：是指電力系統(tǒng)受到小擾動后，不發(fā)生非周期性失步，自動恢復到初始運行狀態(tài)的能力。　　暫態(tài)穩(wěn)定：是指電力系統(tǒng)受到大擾動后，各同步電機保持同步運行并過渡到新的或是恢復到原來穩(wěn)定運行方式的能力。

24、　4.5提高靜態(tài)穩(wěn)定、暫態(tài)穩(wěn)定的措施？　　答：提高靜態(tài)穩(wěn)定措施：采用自動調節(jié)勵磁裝置；減小元件的電抗：采用分裂導線、提高線路的額定電壓等級、采用串聯(lián)電容補償；改善系統(tǒng)的結構和采用中間補償設備。　　提高暫態(tài)穩(wěn)定措施：故障的快速切除和自動重合閘裝置的應用；提高發(fā)電機輸出的電磁功率：對發(fā)電機施行強行勵磁、電氣制動、變壓器中性點經小電阻接地；減少原動機輸出的機械

25、功率。　　三.課程設計總結或結論　　為期兩周的課程設計使我獲益良多：首先初步掌握了Matlab的基本使用方法；之后對自己學過的穩(wěn)態(tài)暫態(tài)知識有了新的理解和感悟，并在整個過程中發(fā)現(xiàn)了很多自己對知識掌握的不足之處；另外，基本掌握了采用計算機對電力系統(tǒng)電磁暫態(tài)過程和機電暫態(tài)過程進行計算的基本方法，并進一步鞏固了計算機編程的能力。<p

26、>　　對系統(tǒng)導納矩陣和阻抗矩陣：首先一定要理解好自導納、互導納以及自阻抗、互阻抗的物理意義，以便其更好地服務于對短路電流計算；其次要掌握形成節(jié)點導納矩陣的正確方法，對節(jié)點導納矩陣求逆就可以形成阻抗矩陣；此外注意嵌套循環(huán)的邏輯順序，盡量使循環(huán)清楚明朗便于執(zhí)行。　　對于短路電流計算部分：首先要理解好用對稱分量法計算不對稱故障短路電流的基本原理；其次，要會畫三序網(wǎng)圖，以及不同故障類型系統(tǒng)復合序網(wǎng)

27、圖；最后，除了會計算故障點的各參數(shù)外，還要會算非故障處電壓、電流值。對稱故障是不對稱故障的一種特殊情況，可以一起討論。在編程中藥注意的是，計算發(fā)電機端電壓時，對稱分量經變壓器后的相位變化。　　靜態(tài)穩(wěn)定部分：首先要熟練掌握并使用公式，了解公式的物理意義；其次要學會利用所學的知識，解決實際的靜態(tài)穩(wěn)定問題，了解影響靜態(tài)問的相關因素，并對改善靜態(tài)穩(wěn)定性的措施有所了解；另外，還需理解試探法的精神，學會使用

28、這一方法解決其他類似問題。　　暫態(tài)穩(wěn)定部分：學會使用改進的歐拉法，通過迭代算出下時段的delta和omega，繪制搖擺曲線，并通過搖擺曲線能判斷出在相應的切除時間下，系統(tǒng)能否保持暫態(tài)穩(wěn)定，可求取極限切除時間；另外，需掌握影響暫態(tài)穩(wěn)定的基本因素，并了解改善暫態(tài)穩(wěn)定的基本措施，學會解決實際問題。　　感謝老師的指導和同學的熱心幫助，才得以解決一個又一個的編程問

29、題、補足知識的缺陷，順利完成這次課程設計。　　四.參考文獻　　《電力系統(tǒng)暫態(tài)分析》，李光琦，中國電力出版社，2006年，第三版　　《電力系統(tǒng)分析》（上、下），何仰贊，華中科技大學出版社，1998年，第二版　　《電力系統(tǒng)故障的計算機輔助

30、分析》重慶大學出版社米麟書等 　　《電力系統(tǒng)潮流計算》天津大學出版社宋文南等 　　《電力系統(tǒng)故障分析》清華大學出版社周榮光　　《短路電流實用計算方法》電力工業(yè)出版社西安交通大學等　　《精通Matlab6.5》北京航空航天大學出版社，張志涌

31、;　　附錄：　　1.Matrix.m，網(wǎng)絡原始數(shù)據(jù)以及節(jié)點導納矩陣、節(jié)點阻抗矩陣的形成程序　　clc;%清空屏幕顯示　　clear;%清空所有變量　　N1=3; %輸入網(wǎng)絡的節(jié)點數(shù)<

32、;p>　　N2=5; %輸入包括所有發(fā)電機節(jié)點的網(wǎng)絡的節(jié)點數(shù)　　k=3; %輸入短路點的節(jié)點號　　%輸入支路各序阻抗，ZZ1表示支路正序阻抗，ZZ2表示支路負序阻抗，ZZ0表示支路零序阻抗　　ZZ1(1,1)=j*0.15;　　ZZ1(1,2)=j*0.1;&

33、lt;/p>　　ZZ1(1,3)=j*0.1;　　ZZ1(2,1)=j*0.1;　　ZZ1(2,2)=j*0.075;　　ZZ1(2,3)=j*0.1;　　ZZ1(3,1)=j*0.1;　　ZZ1(3,

34、2)=j*0.1;　　ZZ1(3,3)=j*0.0;　　ZZ2=ZZ1;　　YYY1=zeros(N1);　　YYY2=zeros(N1);　　YYY0=zeros(N1);

35、;　　for p=1:N1　　for q=1:N1　　if ZZ1(p,q)==0　　YYY1(p,q)=0;　　else 　　YYY1(p,q)=1/ZZ1(p,q);

36、　　end　　end　　end　　YYY2=YYY1;　　ZZ0(1,1)=j*0.0525;

37、　　ZZ0(1,2)=j*0.2;　　ZZ0(1,3)=j*0.2;　　ZZ0(2,1)=j*0.2;　　ZZ0(2,2)=j*0.00263;　　ZZ0(2,3)=j*0.2;　　ZZ0(3,1)=j*0.2

38、;　　ZZ0(3,2)=j*0.2;　　ZZ0(3,3)=j*0.0;　　for p=1:N1　　for q=1:N1　　if ZZ0(p,q)==0　　YYY0(p,q)=0;</

39、p>　　else 　　YYY0(p,q)=1/ZZ0(p,q);　　end　　end　　end

40、;　　%求節(jié)點導納矩陣，Y1表示不計發(fā)電機節(jié)點的正序網(wǎng)絡節(jié)點導納陣，Y2表示不計發(fā)電機節(jié)點的負序網(wǎng)絡節(jié)點導納陣，Y0表示不計發(fā)電機節(jié)點的零序網(wǎng)絡節(jié)點導納陣，　　Y1=zeros(N1);　　for p=1:N1　　for q=1:N1<p&g

41、t;　　if p==q　　for a=1:N1　　Y1(p,q)=Y1(p,q)+YYY1(p,a);%???　　end　　else

42、　　Y1(p,q)=(-1)*YYY1(p,q);　　end　　end　　end　　Y2=zeros(N1);<

43、;p>　　Y2=Y1;　　Y0=zeros(N1);　　for p=1:N1　　for q=1:N1　　if p==q　　for a=1:N1<

44、;/p>　　Y0(p,q)=Y0(p,q)+YYY0(p,a);%???　　end　　else 　　Y0(p,q)=(-1)*YYY0(p,q);

45、end　　end　　end　　Y1　　Y2

46、;　　Y0　　%節(jié)點導納矩陣參考結果，YY1表示計及發(fā)電機節(jié)點的正序網(wǎng)絡節(jié)點導納陣，YY2表示計及發(fā)電機節(jié)點的負序網(wǎng)絡節(jié)點導納陣　　%輸入支路各序阻抗，ZZ1表示支路正序阻抗，ZZ2表示支路負序阻抗，ZZ0表示支路零序阻抗　　ZZ1(1,1)=0;

47、ZZ1(1,2)=j*0.1;　　ZZ1(1,3)=j*0.1;　　ZZ1(1,4)=j*0.0525;　　ZZ1(1,5)=0;　　ZZ1(2,1)=j*0.1;　　ZZ1(2,2)=0;

48、;　　ZZ1(2,3)=j*0.1;　　ZZ1(2,4)=0;　　ZZ1(2,5)=j*0.0263;　　ZZ1(3,1)=j*0.1;　　ZZ1(3,2)=j*0.1;　　ZZ1(3,3)=j*0.0;&

49、lt;p>　　ZZ1(3,4)=0;　　ZZ1(3,5)=0;　　ZZ1(4,1)=ZZ1(1,4);　　ZZ1(4,2)=ZZ1(2,4);　　ZZ1(4,3)=ZZ1(3,4);　　ZZ1(4,4)=j*0.095;&l

50、t;/p>　　ZZ1(4,5)=0;　　ZZ1(5,1)=ZZ1(1,5);　　ZZ1(5,2)=ZZ1(2,5);　　ZZ1(5,3)=ZZ1(3,5);　　ZZ1(5,4)=0;　　ZZ1(5,5)

51、=j*0.0475;　　ZZ2=ZZ1;　　YYY1=zeros(N2);　　YYY2=zeros(N2);　　for p=1:N2　　for q=1:N2<p&

52、gt;　　if ZZ1(p,q)==0　　YYY1(p,q)=0;　　else 　　YYY1(p,q)=1/ZZ1(p,q);　　end

53、　　end　　end　　YYY2=YYY1;　　YY1=zeros(N2);　　for p=1:N2　　for q=1:N2<

54、;b>　　if p==q　　for a=1:N2　　YY1(p,q)=YY1(p,q)+YYY1(p,a);　　end　　else <p&g

55、t;　　YY1(p,q)=(-1)*YYY1(p,q);　　end　　end　　end　　YY2=zeros(N2);

56、;　　YY2=YY1;　　YY1　　YY2　　%正序節(jié)點阻抗　　Z1=eye(N1)/Y1;% eye(N1)是3階的單

57、位矩陣,求節(jié)點導納矩陣Y1的逆陣,即節(jié)點阻抗矩陣Z　　%負序節(jié)點阻抗　　Z2=eye(N1)/Y2;　　%零序節(jié)點阻抗　　Z0=eye(N1)/Y0;

58、2.shortcircuit.m，短路電流計算的程序　　clc;%清空屏幕顯示　　clear;%清空所有變量　　Matrix %計算節(jié)點導納矩陣、節(jié)點阻抗矩陣，形成全局變量　　fault=input('輸入短路類型 ;(0)fault=0為三相短路;\n(1)f

59、ault=1為a相接地短路;\n(2)fault=2為a相經10歐姆電阻接地短路;\n(3)fault=3為bc兩相相間短路;\n(4)faul=4為bc兩相短路接地;\nfault=')　　%第一部分：計算所有節(jié)點的a,b,c三相電壓　　%計根據(jù)故障類型選擇不同的計算公式，計算故障點各序電流　　if f

60、ault==0　　Ik1=1/Z1(k,k);%fault=0為三相短路 　　Ik2=0;　　Ik0=0;　　elseif fault==1 　　Ik1

61、=1/(Z1(k,k)+Z2(k,k)+Z0(k,k));%fault=1為a相接地短路　　Ik2=Ik1;　　Ik0=Ik1;　　elseif fault==2　　Zf=10*60/115^2;

62、%換成標幺值　　Ik1=1/(Z1(k,k)+Z2(k,k)+Z0(k,k)+3*Zf);%fault=2為a相經10歐姆電阻接地短路　　Ik2=Ik1;　　Ik0=Ik1;　　elseif fa

63、ult==3　　Ik1=1/(Z1(k,k)+Z2(k,k));%fault=3為bc兩相相間短路 　　Ik2=(-1)*Ik1;　　Ik0=0; 　　elseif fault==4　　Ik1=1/(Z1(k,k)+Z2

64、(k,k)*Z0(k,k)/(Z2(k,k)+Z0(k,k))); %fault=4為bc兩相短路接地　　Ik2=(-1)*Ik1*Z0(k,k)/(Z0(k,k)+Z2(k,k));　　Ik0=(-1)*Ik1*Z2(k,k)/(Z0(k,k)+Z2(k,k));　　end;

65、;　　%計算各序電流的故障分量　　for p=1:N1 　　if p==k　　I1(p)=-1*Ik1;　　I2(p)=-1*Ik2;　　I0(p)

66、=-1*Ik0;　　else　　I1(p)=0;　　I2(p)=0;　　I0(p)=0;<

67、p>　　end;　　end;　　%計算各點各序電壓的故障分量　　uu1(:,k)=Y1\I1.'; %“\”表示左除 A\B=inv(A)*B inv是求逆陣。即AX=B。 I.'表示轉置

68、　　uu2(:,k)=Y2\I2.'; %uu(:,k)取uu2的第k列 P89　　uu0(:,k)=Y0\I0.';　　%假設故障前空載，各節(jié)點電壓正序分量為1 　　for p=1:N1　　U1(p)=1;<

69、/b>　　end;　　%計算發(fā)生故障后的電壓值:正常運行的分量+故障分量　　u1=U1'+uu1(:,k); %計算所有節(jié)點正序電壓　　u2=uu2(:,k); %計算所有節(jié)點負序電壓

70、　　u0=uu0(:,k); %計算所有節(jié)點零序電壓　　a=-0.5+j*sqrt(3)/2;　　T=[1 1 1 % T為對稱分量法的變換矩陣　　a^2 a 1　　a a^2 1];&

71、lt;/p>　　for p=1:N1　　U=[u1(p) u2(p) u0(p)];　　p　　Uabc=T*U.' %Uabc表示p節(jié)點的a,b,c三相電壓向量　　UUabc

72、=abs(Uabc) %UUabc表示p節(jié)點的a,b,c三相電壓有效值　　end;　　%第二部分：計算支路電流　　for m=1:N1　　for n=1:N1　　if

73、m<n　　mn(1)=m; %mn是什么？？　　mn(2)=n;　　mn　　I1(m,n)=(u1(m)-u1(n))/ZZ1(m,n);%支路正序

74、電流　　I2(m,n)=(u2(m)-u2(n))/ZZ2(m,n);%支路負序電流　　I0(m,n)=(u0(m)-u0(n))/ZZ0(m,n);%支路零序電流　　Iabc=T*[I1(m,n) I2(m,n) I0(m,n)].';%由序分量合成相分量<p&g

75、t;　　Iabc　　abs(Iabc) 　　end; 　　end;　　end;<

76、;p>　　%第三部分：計算發(fā)電機的端電壓，只寫了部分程序，余下的工作由同學們完成。要考慮變壓器繞組是11點接線的　　for p=1:N2　　if p==k　　II(p)=-Ik1;　　else

77、;　　II(p)=0;　　end;　　end;　　vv1(:,k)=YY1\II.';　　vv2(:,k)=

78、YY2\II.'; 　　%假設故障前空載，各節(jié)點電壓正序分量為1 　　for p=1:N2　　V1(p)=1;　　end;　　vv0=zeros(

79、N2);　　%計算故障點發(fā)生故障后的電壓值:正常運行的分量+故障分量　　v1=V1'+vv1(:,k); %計算所有節(jié)點正序電壓　　v2=vv2(:,k); %計算所有節(jié)點負序電壓　　v0=vv0(:,k); %計算所有節(jié)點零序電壓

80、 　　for p=1:N2　　if p==4|p==5　　V=[v1(p)*(exp(j*pi/6)) v2(p)*(exp(-j*pi/6)) v0(p) ];%對稱分量經變壓器后的相位變化　　p<p&

81、gt;　　Vabc=T*V.' %Vabc表示p節(jié)點的a,b,c三相電壓向量 　　VVabc=abs(Vabc) %VVabc表示p節(jié)點的a,b,c三相電壓有效值　　end;　　end;&

82、lt;p>　　3.stability_smallsignal.m，靜態(tài)穩(wěn)定計算的程序　　clear;　　clc;　　U=1; %系統(tǒng)電壓　　Eq0=1.972;

83、%空載電動勢　　UG0=1.21; %機端電壓　　Xe=0.504; %線路電抗值　　Xd=0.982; %同步電抗　　Xdd=0.344; %暫態(tài)電抗　　Xd_all=1.486; %系統(tǒng)

84、電抗(包含同步電抗)　　Xdd_all=0.848; %系統(tǒng)電抗(包含暫態(tài)電抗)　　Ke=0.1;　　P_result=1:56;　　K=0.1:0.1:5.6;　　delta_result=

85、1:56;　　i=1;　　q=1;　　while(Ke<5.7)　　for delta=49:110　　%同學補充程序，計算極限功角及校驗放大倍數(shù)&

86、lt;/p>　　a=(Ke*Ke*Xe*Xe)/(Xd_all*Xd_all)-1;%192頁　　b=2*(Ke*Ke*Xd*Xe*U*cos(delta*pi/180)/(Xd_all*Xd_all)+Eq0+Ke*UG0);　　c=Ke*Ke*Xd*Xd*U*U/(Xd_all*Xd_all)-(Eq0+Ke*UG0)*(

87、Eq0+Ke*UG0);　　Eq=(-b+sqrt(b*b-4*a*c))/(2*a);　　Eqq=Eq*Xdd_all/Xd_all+(Xd_all-Xdd_all)*U*cos(delta*pi/180)/Xd_all;%191頁　　UGd=U*Xd*sin(delta*pi/180)/Xd_all;<

88、/p>　　UGq=Eq*Xe/Xd_all+Xd*U*cos(delta*pi/180)/Xd_all;　　UG=sqrt(UGd*UGd+UGq*UGq); K1=Eqq*U*cos(delta*pi/180)/Xdd_all+U*U*(Xdd_all-Xd_all)*cos(2*delta*pi/180)/(Xd_all*Xdd_all);%191頁</p

89、>　　K2=U*sin(delta*pi/180)/Xdd_all;　　K3=Xdd_all/Xd_all;　　K4=(Xd_all-Xdd_all)*U*sin(delta*pi/180)/Xdd_all;　　K5=UGd*Xd*cos(delta*pi/180)/(UG*Xd_al

90、l)-UGq*Xdd*sin(delta*pi/180)/(UG*Xdd_all);　　K6=UGq*(Xd_all-Xd)/(UG*Xdd_all); 　　Kemax=-K4/K5;　　Kemin=-(Eq*cos(delta*pi/180))/(Xd_all*K3*(K1*K6-K2*K5));</

91、p>　　if ((K1>0)&(Ke>=Kemin)&(Ke<=Kemax))　　Ke_result=Ke;　　delta_result1=delta;　　P_result1=Eq*U*sin(delta*pi/180)/Xd_all; </p&

92、gt;　　else　　s=sprintf('Ke=%f, delta=%d, P=%f\n',Ke_result,delta_result1,P_result1);　　P_result(i)=P_result1;　　delta_res

93、ult(i)=delta_result1;　　i=i+1;　　disp(s);　　break;　　end;</p&

94、gt;　　end;　　K(q)=Ke;　　q=q+1;　　Ke=Ke+0.1; 　　end;&l

95、t;/p>　　subplot(2,1,1);　　plot(K,P_result), grid on;　　xlabel('Ke');　　ylabel('P');　　title('Ke和P之間的關系'

96、)　　subplot(2,1,2);　　plot(K,delta_result), grid on;　　xlabel('Ke');　　ylabel('delta');　　title('Ke和

97、delta之間的關系')　　4．sability_transient.m，暫態(tài)穩(wěn)定計算的程序　　clear;　　clc;　　f=50; %系統(tǒng)額定頻率&l

98、t;p>　　Tj=8.47; %歸算后的發(fā)電機慣性時間常數(shù)　　PT=1; %正常運行時發(fā)電機向無窮大系統(tǒng)傳輸?shù)挠泄β?lt;/p>　　P2M=0.48; %故障存在時發(fā)電機的最大功率　　P3M=1.38; %故障切除后發(fā)電機的最大功率

99、　　%下面是利用改進歐拉法進行逐段計算　　%需要注意故障切除前后電磁功率有躍變　　h=0.05; %設置步長0.05s　　Duration=2; %設置計算時段長度2s 　　CutTime=input('輸入故

100、障切除時間\nCutTime=');　　%故障發(fā)生時的功角變化過程　　Delta(1)=33.92; %初始功角　　Omega(1)=1; %初始轉速　　t(1)=0;<p

101、>　　for i=1:round(CutTime/h) 　　Delta0(i)=(Omega(i)-1)*360*f;　　Omega0(i)=(PT-P2M*sin(Delta(i)*pi/180))/Tj;　　Delta(i+1)=Delta(i)+Delta0(i)*h;<

102、;p>　　Omega(i+1)=Omega(i)+Omega0(i)*h;　　Delta0(i+1)=(Omega(i+1)-1)*360*f;　　Omega0(i+1)=(PT-P2M*sin(Delta(i+1)*pi/180))/Tj;　　Delta0_aver(i+1)=(Delta0(i)+Del

103、ta0(i+1))/2;　　Omega0_aver(i+1)=(Omega0(i)+Omega0(i+1))/2;　　Delta(i+1)=Delta(i)+Delta0_aver(i+1)*h;　　Omega(i+1)=Omega(i)+Omega0_aver(i+1)*h;<p&

104、gt;　　t(i+1)=i*h;　　end;　　%故障切除后的功角變化過程　　for i=round(CutTime/h)+1:round(Duration/h)　　Delta0(i)=(Omega(i)-1)*360*f;</

105、p>　　Omega0(i)=(PT-P3M*sin(Delta(i)*pi/180))/Tj;　　Delta(i+1)=Delta(i)+Delta0(i)*h;　　Omega(i+1)=Omega(i)+Omega0(i)*h;　　Delta0(i+1)=(Omega(i+1)-1)

106、*360*f;　　Omega0(i+1)=(PT-P3M*sin(Delta(i+1)*pi/180))/Tj;　　Delta0_aver(i+1)=(Delta0(i)+Delta0(i+1))/2;　　Omega0_aver(i+1)=(Omega0(i)+Omega0(i+1))/2;

107、　　Delta(i+1)=Delta(i)+Delta0_aver(i+1)*h;　　Omega(i+1)=Omega(i)+Omega0_aver(i+1)*h;　　t(i+1)=i*h;　　end;　　plot

眾賞文庫> 全部分類> 畢業(yè)設計

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
6. 下載文件中如有侵權或不適當內容，請與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

電力系統(tǒng)暫態(tài)上機計算課程設計報告

文檔簡介

溫馨提示

最新文檔

評論

電力系統(tǒng)暫態(tài)上機計算課程設計報告

文檔簡介

溫馨提示

最新文檔

評論

免費下載