短路电流计算程序(共12页).doc

《短路电流计算程序(共12页).doc》由会员分享，可在线阅读，更多相关《短路电流计算程序(共12页).doc（12页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、精选优质文档-倾情为你奉上短路电流计算程序2 OPEN BH2.txt FOR OUTPUT AS #13 REM B44 INPUT I96 REM GREETING-1 PROGREM8 READ N, L, B, L1, G10 o = N + 112 DIM W(o), Z(o), a(o), r(B + 1), x(B + 1), i(N), J(N), t(N)14 DIM L(B + 1), k(B + 1), N(B + 1), M(B + 1), s(B + 1)16 FOR i = 1 TO N18 READ W(i)20 NEXT i22 M = W(N): m9 = M

2、 + N + 124 DIM C(m9), D(m9), H(m9), y(m9), E(o), F(o), P(o)26 FOR i = 1 TO m928 H(i) = 030 NEXT i32 FOR i = 1 TO L34 READ J, k, r, x, b136 s = W(J): t = W(k)38 m1 = s - J + k40 IF J k THEN 4442 m1 = t - k + J44 IF r = 0 THEN GOSUB 6846 IF r 0 THEN GOSUB 7848 NEXT i50 IF I9 = 1 THEN 6652 FOR J = 1 TO

3、 N54 J0 = W(J) - J: j1 = J - W(J) + 156 IF J 1 THEN j1 = j1 + W(J - 1)58 FOR k = j1 TO J60 PRINT Y; J; k; =; H(J0 + k)62 NEXT k64 NEXT J66 GOTO 9068 H(s) = H(s) + b1 - 1 / x70 H(t) = H(t) + b1 - 1 / x72 IF i * k = 0 THEN 7574 H(m1) = H(m1) + 1 / x75 IF J = (N - G) THEN 7776 t(J) = x77 RETURN78 IF AB

4、S(r) = 2 THEN 8480 H(s) = H(s) - 1 / x: H(t) = H(t) - 1 / (b1 * b1 * x)82 GOTO 8684 H(s) = H(s) - 1 / (b1 * b1 * x): H(t) = H(t) - 1 / x86 H(m1) = H(m1) + 1 / (b1 * x)88 RETURN90 FOR i = 1 TO m992 C(i) = 0: D(i) = 094 NEXT i96 FOR i = 1 TO N98 READ E1, F1100 E(i) = E1: F(i) = F1102 NEXT i104 FOR o =

5、 1 TO L106 READ J, k, I1, I2107 IF J * k = 0 THEN 113108 IF J k THEN 144142 m1 = t - k + J144 IF r = 0 THEN GOSUB 168146 IF r 0 THEN GOSUB 178148 NEXT i150 IF I9 = 1 THEN 166152 FOR J = 1 TO B154 J0 = r(J) - J: j1 = J - r(J) + 1156 IF J 1 THEN j1 = j1 + r(J - 1)158 FOR k = j1 TO J160 PRINT Y0; J; k;

6、 U(J0 + k)162 NEXT k164 NEXT J166 GOTO 196168 U(s) = U(s) + b1 - 1 / x170 U(t) = U(t) + b1 - 1 / x172 IF J * k = 0 THEN 176174 U(m1) = U(m1) + 1 / x176 RETURN178 IF J * k = 0 THEN 188180 IF ABS(r) = 2 THEN 184182 U(s) = U(s) - 1 / x: U(t) = U(t) - 1 / (b1 * b1 * x): GOTO 185184 U(s) = U(s) - 1 / (b1

7、 * b1 * x): U(t) = U(t) - 1 / x185 U(m1) = U(m1) + 1 / (b1 * x)186 RETURN188 IF k = 0 THEN 192190 U(t) = U(t) - 1 / x: t(k) = x: GOTO 194192 U(s) = U(s) - 1 / x: t(J) = x194 GOTO 186196 PRINT +-7-+198 INPUT X2200 IF X2 = 0 THEN 252202 PRINT J,K,L: INPUT J, k, L204 E1 = (1 - L) * E(J) + L * E(k): F1

8、= (1 - L) * F(J) + L * F(k)206 E(N + 1) = E1: F(N + 1) = F1208 IF J k THEN 220210 J0 = W(J) - J: Y0 = H(J0 + k)212 H(J0 + k) = 0: C1 = C(J0 + k): D1 = D(J0 + k)214 C(J0 + k) = 0: D(J0 + k) = 0216 C(W(N) + J) = -C1: C(W(N) + k) = C1: D(W(N) + J) = -D1: D(W(N) + k) = D1218 GOTO 226220 K0 = W(k) - k: Y

9、0 = H(K0 + J): H(K0 + J) = 0222 C1 = C(K0 + J): D1 = D(K0 + J): C(K0 + J) = 0: D(K0 + J) = 0224 C(W(N) + J) = C1: C(W(N) + k) = -C1: D(W(N) + J) = D1: D(W(N) + k) = -D1226 FOR i = 1 TO N228 IF i J THEN 232230 H(M + i) = Y0 / L: GOTO 238232 IF i k THEN 236234 H(M + i) = Y0 / (1 - L): GOTO 238236 H(M

10、+ i) = 0238 NEXT i240 H(M + N + 1) = -Y0 * (1 / L + 1 / (1 - L)242 J2 = W(J): K2 = W(k)244 H(J2) = H(J2) + Y0 - Y0 / L246 H(K2) = H(K2) + Y0 - Y0 / (1 - L)248 W(N + 1) = W(N) + N + 1: N = N + 1250 GOTO 256252 PRINT AT WHICH BUS: INPUT J254 E1 = E(J): F1 = F(J)256 FOR i = 1 TO N258 IF X2 = 0 THEN 262

11、260 s = N: GOTO 264262 s = J264 IF i s THEN 268266 Z(i) = -1: GOTO 270268 Z(i) = 0270 NEXT i272 GOSUB 404274 Z0 = Z(s)276 PRINT +-2-+: INPUT X3278 IF X3 = 0 THEN 332280 IF X2 = 0 THEN 312282 J2 = r(J): K2 = r(k)284 IF J k THEN 288286 Y2 = U(J2 - J + k): U(J2 - J + k) = 0: GOTO 290288 Y2 = U(K2 - k +

12、 J): U(K2 - k + J) = 0290 I2 = r(B)292 FOR i = 1 TO B294 IF i J THEN 298296 U(I2 + i) = Y2 / L: GOTO 304298 IF i k THEN 302300 U(I2 + i) = Y2 / (1 - L): GOTO 304302 U(I2 + i) = 0304 NEXT i306 U(I2 + B + 1) = -Y2 * (1 / L + 1 / (1 - L)308 U(J2) = U(J2) + Y2 - Y2 / L: U(K2) = U(K2) + Y2 - Y2 / (1 - L)

13、310 r(B + 1) = r(B) + B + 1: B = B + 1312 FOR i = 1 TO B314 IF X2 = 0 THEN 318316 s = B: GOTO 320318 s = J320 IF i s THEN 324322 x(i) = -1: GOTO 326324 x(i) = 0326 NEXT i328 GOSUB 404330 X0 = x(s): GOTO 346332 PRINT TTT334 INPUT X4336 IF X4 = 0 THEN 372338 Z4 = 0: m1 = 1340 W1 = 0: W2 = 0342 PRINT R

14、ESULT OF 3 SHORT-CURCUIT344 GOTO 378346 PRINT TTI348 INPUT X5350 IF X5 = 1 THEN 360352 Z4 = Z0 + X0: m1 = 3354 W1 = 1: W2 = 1356 PRINT RESULT OF SINGLE CONNECTING THE EARTH358 GOTO 378360 Z4 = Z0 * X0 / (Z0 + X0)362 m1 = SQR(3 * (1 - Z0 * X0 / (Z0 + X0) * (Z0 + X0)364 W1 = -X0 / (Z0 + X0)366 W2 = -Z

15、0 / (Z0 + X0)368 PRINT RESULT OF DOUBLE CONNECTING THE EARTH370 GOTO 378372 Z4 = Z0: m1 = SQR(3)374 W1 = -1: W2 = 0376 PRINT RESULT OF DOUBLE SHORT-CURCUIT378 Z5 = Z4 + Z0: I1 = F1 / Z5: I2 = -E1 / Z5380 I3 = W1 * I1: I4 = W1 * I2382 I5 = W2 * I1: I6 = W2 * I2384 I7 = m1 * SQR(I1 * I1 + I2 * I2)386

16、E1 = E1 + I2 * Z0: F1 = F1 - Z0 * I1388 E2 = Z0 * I4: F2 = -Z0 * I3390 E3 = X0 * I6: F3 = -X0 * I5391 PRINT #1, IK+=; I1; +J; I2392 PRINT #1, IK-=; I3; +J; I4393 PRINT #1, IK0=; I5; +J; I6394 PRINT #1, IMAX=; I7396 PRINT #1, UK+=; E1; +J; F1398 PRINT #1, UK-=; E2; +J; F2399 PRINT #1, UK0=; E3; +J; F

17、3400 GOTO 516404 N0 = N406 FOR i = 1 TO N408 IF X3 = 1 THEN 412410 a(i) = W(i): P(i) = Z(i): GOTO 416412 IF i B THEN 416414 a(i) = r(i): P(i) = x(i)416 NEXT i418 IF X3 = 1 THEN 430420 I7 = a(N)422 FOR i = 1 TO I7424 y(i) = H(i)426 NEXT i428 GOTO 440430 I7 = a(B)432 FOR i = 1 TO I7434 y(i) = U(i)436

18、NEXT i438 IF X3 = 1 THEN N = B440 FOR o = 1 TO N442 o0 = a(o) - o: O1 = o - a(o) + 1444 IF o 1 THEN O1 = O1 + a(o - 1)446 FOR q = O1 TO o448 Q0 = a(q) - q: Q1 = q - a(q) + 1450 IF q 1 THEN Q1 = Q1 + a(q - 1)452 IF Q1 = O1 THEN 456454 Q1 = O1456 FOR i = Q1 TO q - 1458 IF i q - 1 THEN 462460 y(o0 + q)

19、 = y(o0 + q) - y(o0 + i) * y(Q0 + i) / y(a(i)462 NEXT i464 NEXT q466 NEXT o468 FOR o = 1 TO N470 o0 = a(o) - o: O1 = o - a(o) + 1472 IF o 1 THEN O1 = O1 + a(o - 1)474 FOR q = O1 TO o - 1476 IF q o - 1 THEN 480478 P(o) = P(o) - y(o0 + q) * P(q)480 NEXT q482 P(o) = P(o) / y(a(o)484 NEXT o486 FOR o = N

20、 TO 1 STEP -1488 o0 = a(o) - o: O1 = o - a(o) + 1490 IF o 1 THEN O1 = O1 + a(o - 1)492 FOR q = O1 TO o - 1494 IF q o - 1 THEN 498496 P(q) = P(q) - y(o0 + q) * P(o) / y(a(q)498 NEXT q500 NEXT o502 IF X3 = 1 THEN N = B504 FOR i = 1 TO N506 IF X3 = 1 THEN 510508 Z(i) = P(i): GOTO 512510 x(i) = P(i)512

21、NEXT i514 RETURN516 FOR i = 1 TO N0518 E1 = E(i) + Z(i) * I2: F1 = F(i) - Z(i) * I1520 E2 = Z(i) * I4: F2 = -Z(i) * I3522 IF X2 = 1 THEN 530524 IF i B THEN 528526 E3 = x(i) * I6: F3 = -x(i) * I5: GOTO 542528 E3 = 0: F3 = 0: GOTO 542530 IF i B THEN 526532 IF i B THEN 536534 E4 = x(B) * I6: F4 = -x(B)

22、 * I5536 E3 = 0: F3 = 0538 IF i N0 THEN 542540 E3 = E4: F3 = F4542 A1 = E1 + E2 + E3: A2 = F1 + F2 + F3544 A4 = (F1 + F2) * COS(3.14159 * 2 / 3) + (E2 - E1) * SIN(3.14159 * 2 / 3) + F3546 A3 = (E1 + E2) * COS(3.14159 * 2 / 3) + (F1 - F2) * SIN(3.14159 * 2 / 3) + E3548 A5 = (E1 + E2) * COS(3.14159 *

23、2 / 3) + (F2 - F1) * SIN(3.14159 * 2 / 3) + E3550 A6 = (F1 + F2) * COS(3.14159 * 2 / 3) + (E1 - E2) * SIN(3.14159 * 2 / 3) + F3551 PRINT #1, U; i; +=; E1; +J; F1552 PRINT #1, U; i; -=; E2; +J; F2553 PRINT #1, U; i; 0=; E3; +J; F3554 PRINT #1, UA=; A1; +J; A2555 PRINT #1, UB=; A3; +J; A4556 PRINT #1,

24、 UC=; A5; +J; A6557 NEXT i558 FOR J = 2 TO N0560 j1 = W(J) - J562 IF J = B THEN J2 = r(J) - J564 J3 = J - W(J) + W(J - 1) + 1566 IF J = B THEN J4 = J - r(J) + r(J - 1) + 1568 FOR k = J3 TO J - 1569 IF H(j1 + k) = 0 THEN 624570 C1 = -I1 * (Z(J) - Z(k) * H(j1 + k): C2 = -I2 * (Z(J) - Z(k) * H(j1 + k)5

25、72 C3 = -I3 * (Z(J) - Z(k) * H(j1 + k): C4 = -I4 * (Z(J) - Z(k) * H(j1 + k)574 IF X3 = 0 THEN 604576 IF X2 = 1 THEN 586578 IF J = B THEN 582580 C5 = 0: C6 = 0: GOTO 604582 C5 = -I5 * (x(J) - x(k) * U(J2 + k): C6 = -I6 * (x(J) - x(k) * U(J2 + k): GOTO 604583 IF t(k) 0 THEN 585584 C5 = L(k): C6 = s(k)

26、: GOTO 604585 C5 = I5 * x(k) / t(k): C6 = I6 * x(k) / t(k): GOTO 604586 IF J B THEN 582588 IF J B THEN 596590 FOR q = J4 TO J - 1592 L(q) = -I5 * (x(J) - x(q) * U(J2 + q): s(q) = -I6 * (x(J) - x(q) * U(J2 + q)594 NEXT q596 C5 = 0: C6 = 0598 IF J N0 THEN 604600 IF k B - 1 THEN 604602 GOTO 584604 C1 =

27、 C1 + C(j1 + k): C2 = C2 + D(j1 + k): GOTO 606606 A1 = C1 + C3 + C5: A2 = C2 + C4 + C6608 b1 = (C1 + C3) * COS(3.14159 * 2 / 3) + (C2 - C4) * SIN(3.14159 * 2 / 3) + C5610 B2 = (C2 + C4) * COS(3.14159 * 2 / 3) + (C3 - C1) * SIN(3.14159 * 2 / 3) + C6612 D1 = (C1 + C3) * COS(3.14159 * 2 / 3) + (C4 - C2

28、) * SIN(3.14159 * 2 / 3) + C5614 D2 = (C2 + C4) * COS(3.14159 * 2 / 3) + (C1 - C3) * SIN(3.14159 * 2 / 3) + C6616 PRINT #1, IF+(; J; -; k; )=; C1; +J; C2618 PRINT #1, IF-(; J; -; k; )=; C3; +J; C4620 PRINT #1, IF0(; J; -; k; )=; C5; +J; C6621 PRINT #1, IA(; J; -; k; )=; A1; +J; A2622 PRINT #1, IB(;

29、J; -; k; )=; b1; +J; B2623 PRINT #1, IC(; J; -; k; )=; D1; +J; D2624 NEXT k626 IF X3 = 1 THEN 644628 IF X4 = 0 THEN 644630 IF X2 = 0 THEN 634632 IF J (N0 - G - 1) THEN 640634 IF J = (N0 - G) THEN 644636 i(J) = i(J) + Z(J) * I1 / t(J): J(J) = J(J) + Z(J) * I2 / t(J)638 PRINT #1, IG; J; =; i(J); +J; J

30、(J): GOTO 644640 IF J = N0 THEN 644642 IF J N0 THEN 636644 NEXT J646 END注：以下为数据输入部分，此为一个电网一种运行方式下的例子。700 DATA 72,104,4,7,10702 DATA 1,3,6,10,13,19,21,23,32,35,40,46,53,61,70,80,91,92,93,94,108,110,113,131,150,170,191,213,236,260,285,311,338,366,395,425,456,488,521,555,592,594,597,601,606,612,619,627

31、,672,718,765,813,862,912,960,1009,1059,1110,1162,1215,1269,1324,1381,1439,1498,1558,1619,1681,1744,1805,1870,1941706 DATA 1,2,0,.0728,0708 DATA 1,3,.0292,.0368,0710 DATA 1,4,.0145,.0183,0712 DATA 1,6,0,.3270,0714 DATA 1,9,0,.3270,0716 DATA 2,72,0,.0381,0718 DATA 3,4,.0145,.0183,0720 DATA 3,5,0,.985,

32、0722 DATA 5,41,.8314,2.255,0724 DATA 5,49,.4157,1.1275,0726 DATA 5,50,.03495,.01178,0728 DATA 5,51,.20972,.0707,0730 DATA 5,52,.0373,87.004,0732 DATA 5,53,.3263,32.0375,0734 DATA 5,54,.3263,32.0375,0736 DATA 6,7,0,.5819,0738 DATA 7,8,0,.1082,0740 DATA 7,11,.0799,3.8492,0742 DATA 7,12,.0307,.9681,074

33、4 DATA 7,13,.2305,3.9307,0746 DATA 7,14,.2305,3.9307,0748 DATA 7,15,.0091,3.8108,0750 DATA 7,16,.00152,.6351,0752 DATA 7,17,.1383,23.8808,0754 DATA 7,63,.03535,.0191,0756 DATA 7,68,.068,.0366,0758 DATA 7,67,.229,.1239,0760 DATA 7,66,.2244,.1214,0762 DATA 7,65,.0499,.0333,0764 DATA 7,64,.03535,.0191,

34、0766 DATA 7,24,.0707,.0382,0768 DATA 7,25,.0814,.0441,0770 DATA 7,26,.0615,.0333,0772 DATA 7,27,.027,.0145,0774 DATA 7,28,.106,.0574,0776 DATA 7,29,.1844,.0998,0778 DATA 7,30,.0692,.0374,0780 DATA 7,31,.2228,.1205,0782 DATA 7,32,.04225,.02286,0784 DATA 7,33,.1022,.0553,0786 DATA 7,34,.0745,.0403,078

35、8 DATA 7,35,.02305,.01247,0790 DATA 7,36,.15369,.0831,0792 DATA 7,37,.0307,.0166,0794 DATA 7,38,.0499,.02702,0796 DATA 7,39,.04666,.02411,0798 DATA 7,40,.05533,.0299,0799 DATA 7,71,.1783,.1135,0800 DATA 8,9,0,.5819,0801 DATA 8,0,0,.2114,0802 DATA 8,10,0,.4687,0803 DATA 8,21,.2619,.1746,0804 DATA 8,5

36、5,0,.5174,0806 DATA 8,56,0,1.6236,0808 DATA 8,57,1.222,10.3761,0810 DATA 8,58,0,.31874,0812 DATA 8,59,0,.8135,0814 DATA 8,60,0,87.1915,0816 DATA 8,61,0,.4779,0820 DATA 8,62,.9769,.373,0822 DATA 8,69,.0135,.00608,0826 DATA 10,70,.012,.0076,0828 DATA 18,62,0,2.64007,0830 DATA 19,62,0,40.126,0832 DATA

37、20,62,0,4.1958,0834 DATA 21,22,0,1.125,0836 DATA 21,23,0,2.94033,0838 DATA 24,0,0,69.1155,0840 DATA 25,0,0,90.099,0842 DATA 26,0,0,70.96,0844 DATA 27,0,0,38.573,0846 DATA 28,0,0,95.2934,0848 DATA 29,0,0,110.542,0850 DATA 30,0,0,111.37,0852 DATA 31,0,0,111.37,0854 DATA 32,0,0,169.9445,0856 DATA 33,0,

38、0,156.966,0858 DATA 34,0,0,146.87,0860 DATA 35,0,0,73.0382,0862 DATA 36,0,0,190.2789,0864 DATA 37,0,0,38.0558,0866 DATA 38,0,0,83.7852,0868 DATA 39,0,0,97.8652,0870 DATA 40,0,0,135.5427,0872 DATA 42,0,0,17.778,0874 DATA 41,49,.4157,1.1275,0876 DATA 41,42,.03495,.01178,0878 DATA 41,43,.20972,.0707,08

39、80 DATA 41,44,.0932,29.851,0882 DATA 41,45,.0932,29.851,0884 DATA 41,46,1.8141,14.1602,0886 DATA 41,47,1.5873,7.048,0888 DATA 41,48,1.7,7.1493,0889 DATA 43,0,0,106.668,0890 DATA 63,0,0,4.68955,0892 DATA 64,0,0,9.22,0894 DATA 65,0,0,7.9,0896 DATA 66,0,0,7.9,0898 DATA 67,0,0,7.9,0900 DATA 68,0,0,7.9,0907 DATA 69,0,0,4.4643,0908 DATA 70,0,0,.363,0910 DATA 71,0,0,.88,0912 DATA 72,0,0,.04605,0914 DATA 1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0916 DATA 1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0918 DATA 1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1