【文章內(nèi)容簡介】 0)J1=0 (21)從Ｉ＝１步長１到L作 (22)IY1(K)=I N (23) YJ3=IY2(K) (24)從Ｉ＝１步長１到J1作 (25)K2=K0+K1 (26)J3=IY(K2) Y (27) NJ=J+1:Y1(J)=YZ1(K):Y2(J)=YZ2(K):IY(J)=J3:J1=J1+1 (28)Y1(K2)=Y1(K2)+YZ1(K):Y2(K2)=Y2(K2)+YZ2(K) (29)IN(I)=J1:K0=K0+J1 (30)由于導(dǎo)納矩陣的自導(dǎo)納是在追加支路的過程中累加而造成的，所以儲存自導(dǎo)納的D1D12數(shù)組要預(yù)先置零?？颍?）（2）是 完成這一工作的?？颍?）中L是非零互導(dǎo)納元素的記數(shù)單元，開始置零?？颍?）K循環(huán)表示逐次追加一條支路?？颍?）中 IG是當(dāng)前支路狀態(tài)數(shù)的臨時記存單元?？颍?）判別當(dāng)前支路是否處于停 運狀態(tài)，停運狀態(tài)不影響導(dǎo)納矩陣。框（7） 中I、J分別為當(dāng)前支路兩端節(jié)點號的臨時記存單元，R、X分別為當(dāng)前支路正序電阻、電抗的臨時記存單元，B為當(dāng)前支路正序電納或非標(biāo)準(zhǔn)變比的臨時 記存單元。有時會出現(xiàn)這樣的支路，只有零序阻抗，而無正序阻抗，即正序支路斷開，框（9）A單元為0，就是正序支路情況。框（10）計算支路導(dǎo)納。它等于支路阻抗的倒數(shù)。其實部存于 GIJ單元虛部存于BIJ單元?？颍?1）到（14） 分別計算輸電線支路或變壓器支路I節(jié)點自導(dǎo)納，J節(jié)點自導(dǎo)納和它們之間的互導(dǎo)納。I節(jié)點自導(dǎo)納的實部與虛部分別存于GI 和BI單元。J節(jié)點；；自導(dǎo)納實部虛部分別記存與GJ與BJ單元，互導(dǎo)納的實部與虛部去負(fù)號分別記存于GIJ、BIJ單元。如果當(dāng)前支路是對地支路，則執(zhí)行框（15），將當(dāng)前支路導(dǎo)納累加到I節(jié)點的自導(dǎo)納上?？颍?6）將當(dāng)前非對地支路導(dǎo)納累加到I、J節(jié)點的自導(dǎo)納上。每增加一條非對地支路，在導(dǎo)納矩陣中要增加兩個非零互導(dǎo)納元素。每增加一個非零互導(dǎo)納，非零互導(dǎo)納元素的記數(shù)單元L增一?？颍?7）記存互導(dǎo)納，框（18）記存 導(dǎo)納互。導(dǎo)納互的實部、虛部分別記存于YZYZ 2 數(shù)組，其行號和列號分別記存于IYIY2數(shù)組。對所有的支路掃描一遍，便形成正序自導(dǎo)納和不規(guī)則的互導(dǎo)納。應(yīng)注意，L單元存放的是非互 導(dǎo)納元素的總個數(shù)?？颍?9）中J是有規(guī)則非零互導(dǎo)納元素的記事單元。它是逐個累加的。K0也是有規(guī)則非零互導(dǎo)納的記數(shù)單位，但它是逐個累加的。它們開始均置零?？颍?0）——（30）I循環(huán)實現(xiàn)按行號由小到大將非零互導(dǎo)納元素排列在 YY2數(shù)組中。具體做法如下：框 （21）中JI是當(dāng)前行I非零互導(dǎo)納元素的記數(shù)單元，開始置零?？颍?2）、（23）K循環(huán)逐個檢查不規(guī)則非零互導(dǎo)納，其行號IY1（K）是否等于當(dāng)前行號I。若等于，則表明該非零互導(dǎo)納YZI（K）是第I行的元素?？颍?4）將該非零互導(dǎo)納的列好從I Y2（K）單元取出，存于J3單元?？颍?5）、（27） KI循環(huán)對當(dāng)前行I已有規(guī)則 排列在YY2數(shù)組中的互導(dǎo)納元素（它在 YY2數(shù)組中的位置 為K0+K1）逐個進(jìn)行檢查，是否有列號等于J3的元素。若有，就是前述并列支路的情況，執(zhí)行框（29）將非零互導(dǎo)納YZ1（K）、YZ2（K）累加到Y(jié)1（K2）、Y2（K2）單元中。若對當(dāng)前行I已有規(guī)則排列在 YY2數(shù)組中的所有元素掃描一遍，都找不到列號等于J3的元素，則表明YZ1（K）、YZ2（K）非并聯(lián)支路產(chǎn)生的互導(dǎo)納，執(zhí)行框（28）?？颍?8）首先使用規(guī)則排列的非零互導(dǎo)納元素總個數(shù)增加一個，即J=J+1。然后將非零互導(dǎo)納YZ1（K）、YZ2（K）存入Y1（J）、 Y2（J）單元，其列號存于IY（J）單元。同時，使當(dāng)前行I非零互導(dǎo)納元素的個數(shù)增加一個，即J1=J1+1。對不規(guī)則存放在YZYZ2數(shù)組的非零互導(dǎo)納掃描一遍，便找出當(dāng)前行I所有的非零互導(dǎo)納，且有規(guī)則的排列在Y1=Y2數(shù)組中?？颍?0）先將當(dāng)前行I非零元素的個數(shù)J1存于IN數(shù)組。IN數(shù)組是用來存放正序?qū)Ъ{矩陣每行非零互導(dǎo)納元素的個數(shù)的。然后按 行累計有規(guī)則排列的非零互導(dǎo)納元素的總個數(shù)，即K0=K0+J1。用K0=J語句代替也行。對所有的行掃描一遍，便得到導(dǎo)納矩陣有規(guī)則排列的非零互導(dǎo)納元素。 故障程序簡介短路電流計算是電力系統(tǒng)最常用的計算之一。分單相接地短路，兩相短路,兩相接地短路和三相短路四種。根據(jù)電力系統(tǒng)暫態(tài)分析課程可知，故障計算是用對稱分量法進(jìn)行的。首先，假設(shè)網(wǎng)絡(luò)是線性的，即系統(tǒng)各元件的參數(shù)是恒定的。應(yīng)用迭加原理，將三相網(wǎng)絡(luò)分解為正負(fù)零序三個網(wǎng)絡(luò)。并假定正常情況下網(wǎng)絡(luò)是對稱的，即三個網(wǎng)絡(luò)是各自獨立的。再應(yīng)用迭加原理，將各序網(wǎng)的電壓，電流分解為正常分量和故障分量。最后根據(jù)故障點故障類型的邊界條件，將三個序網(wǎng)連成一個完整的網(wǎng)絡(luò)，應(yīng)用線性交流電路理論，計算三個序網(wǎng)電壓電流的故障分量，在與正常分量相加，便求得三序電壓電流的實際值。由三序電壓和電流就可計算出三相電和電流。故障點三序電流：單相接地電流： （）兩相短路： （）三相短路： （）三序網(wǎng)各節(jié)點的實際電壓 （）實用短路計算有下列基本假設(shè)：1） 不計系統(tǒng)元件的電阻，只計元件的電抗，或用阻抗的模替代電抗。2） 不計輸電線對地電納。3） 不壓器的非標(biāo)準(zhǔn)變比。4） 不計負(fù)荷，或負(fù)荷用恒定電抗表示。5） 發(fā)電機次暫態(tài)電動勢的幅值標(biāo)么均為一，幅角均為零。所以正序網(wǎng)所有節(jié)點正常電壓均為一。短路電流計算主程序SCP包括九個子程序，SCGP、SCOP、SCKP、SCDP、SCCP、SCJP、SCWP、SCRP、SCOM。（1）SCGP子程序的功能是輸入原始數(shù)據(jù)。實用短路計算程序需要的原始數(shù)據(jù)有：支路參數(shù)；發(fā)電機參數(shù)。建立一個數(shù)據(jù)文件，名為SCDATA。DAT。該數(shù)據(jù)文件第一行有三個數(shù)據(jù)；網(wǎng)絡(luò)節(jié)點數(shù)N，網(wǎng)絡(luò)支路數(shù)M，發(fā)電機臺數(shù)IQ。第三行開始為支路數(shù)據(jù)。每條支路占據(jù)一行，每行有五個數(shù)據(jù)：1）支路狀態(tài)數(shù)IZA；2）支路一端的節(jié)點號IZ1；3）支路另一端的節(jié)點號IZ2；4）支路正序電抗Z2；5）支路零序電抗Z5。緊接著是發(fā)電機數(shù)據(jù)，每臺發(fā)電機占一行。每行有四個數(shù)據(jù)：1發(fā)電機狀態(tài)數(shù)IWGA；2發(fā)電機連接的節(jié)點號IWG；3發(fā)電機的負(fù)序電抗W3；4發(fā)電機的次暫態(tài)電抗W4。（2）SCOP子程序的功能是節(jié)點優(yōu)化，目的是尋找一種節(jié)點編號方式，使在解網(wǎng)絡(luò)方程中注入新的非零元素最少，能夠節(jié)省內(nèi)存的開銷和提高計算速度。本程序采用的是半動態(tài)優(yōu)化程序。這種方法的基本思想是將連接支路最少的節(jié)點先編號，然后消去該節(jié)點，每消一個節(jié)點，尚未編號的節(jié)點所連接支路數(shù)就隨之變化，再在未編號的節(jié)點中查找連接之路數(shù)最少的節(jié)點挨編號，如此反復(fù)，直至消去所有節(jié)點。 （3）SCKP子程序的功能是將支路兩端舊節(jié)點號轉(zhuǎn)換為新節(jié)點號。 （4）SCDP子程序的功能是形成正序?qū)Ъ{矩陣，SCCP子程序的功能是形成零序?qū)Ъ{矩陣，這兩個子程序與LDP和LCP子程序基本一樣，不同之處是，這里不計電阻，所以導(dǎo)納只有虛部，無實部。還有不計輸電線的對地電納和變壓器的非標(biāo)準(zhǔn)變比。（5）SCJP子程序的功能是形成導(dǎo)納矩陣因子表。計算三個序網(wǎng)阻抗矩陣第f列元素。因子表的下三角矩陣的元素可以由上三角矩陣的元素求得，所以在形成因子表時，只需要貯存因子表的上三角矩陣的非零元素。對角元素與非對角元素分開存放。同時要記存非對角非零元素的列號和每行非對角非零元素的個數(shù)。三個序網(wǎng)的導(dǎo)納矩陣都要形成因子表。并編成程序供各程序調(diào)用。（6）SCWP子程序的功能是對常數(shù)項向量進(jìn)行前代和回代運算。（7）SCRP子程序的功能是輸出短路電流計算結(jié)果。短路電流計算結(jié)果有四部分：1短路點電壓；2分量，即正，負(fù)，零序分量和a,b,c三相量，且輸出其幅值。如何判斷短路電流計算結(jié)果是否正確，從三方面校驗：1)是否滿足故障點的邊界條件；2)檢查三序電壓分布情況。正確的分布是，正序網(wǎng)故障電壓最低，越靠近電源，電壓越高；負(fù)序網(wǎng)與零序網(wǎng)故障點電壓的幅值最高，離開故障點越遠(yuǎn)的，電壓越低。3)不論是序電流還是相電流，流入一個節(jié)點的電流之和應(yīng)等于零。 復(fù)雜故障計算的主要程序及其實現(xiàn)復(fù)雜計算的主程序FAP有十一個子程序，現(xiàn)將各個子程序的功能及具體容敘述如下。（1）FGP子程序的功能是輸入原始數(shù)據(jù)。復(fù)雜故障計算所需要的原始數(shù)據(jù)包括支路參數(shù)，發(fā)電機參數(shù)和，正常運行狀態(tài)下發(fā)電機的有功、無功功率和節(jié)點電壓，負(fù)荷的有功、無功功率和節(jié)點電壓。（2）FBP子程序的功能是輸入故障信和分解故障線路故障信息輸入方式有兩種：一是調(diào)用磁盤中的數(shù)據(jù)文件。二是人機對話輸入故障信息。不論選用那種方式，故障信息輸入后都可以進(jìn)行修改。屏幕上提問修改哪一重故障后，若鍵入2，則表示要修改2重故障的信息。若鍵入0，則表示修改結(jié)束。 （3）FOP子程序的功能是節(jié)點編號優(yōu)化，與潮流計算程序中的子程LOP子程序基本一樣。（4）FKP子程序的功能是將舊節(jié)點號轉(zhuǎn)換為新節(jié)點號。（5）FDP子程序的功能是形成正序?qū)Ъ{矩陣，與潮流計算中的子程序LDP基本一樣。（6）FCP子程序的功能是形成零序?qū)Ъ{矩陣。（7）FIP子程序的功能是進(jìn)行出值計算。（8）FYP子程序的功能是計算修改導(dǎo)支路納矩陣IV=1：B=SQRT（3）/2 （1）從 K=1步長1到N5作 （2）IG=FA（K、1）：R=FA（K、4）：X=FA（K、5） （3）IV=1^IG≥7^IG≤0 N （4） YIV=2 （5）IG≥1^IG≤6 N （6） YIV=3 （7）IG≥10^IG≤15 N （8） YIV=3 （9）A=3*（R*F+X*X） （10）A=0 （11） N YR=0：：A=3*（R*R+X*X） （12） R1=R/A：X1= —X/A （13）R2=B**R1：X2=*X1B*R1 （14）R3=*R1B*X1：X3=B**X1