From b18347ffc9db9641e215995edea1c04c363b2bdf Mon Sep 17 00:00:00 2001 From: Angelo Rossi Date: Wed, 21 Jun 2023 12:04:16 +0000 Subject: Initial commit. --- benchmarks/dc26.dat | 551 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 551 insertions(+) create mode 100644 benchmarks/dc26.dat (limited to 'benchmarks/dc26.dat') diff --git a/benchmarks/dc26.dat b/benchmarks/dc26.dat new file mode 100644 index 0000000..9793f5d --- /dev/null +++ b/benchmarks/dc26.dat @@ -0,0 +1,551 @@ +BEGIN NEW DATA CASE +C BENCHMARK DC-26 +C Test of Type-59 S.M. dynamics, with load flow ("FIX SOURCE" request) +C determination of initial conditions that observe given power constraints. +C 1st of 5 subcases involving load flow solutions; 1st of 2 machine cases. +C Answers change slightly on 10 February 1999 following the massive +C changes from TEPCO (Tokyo Electric Power Company) in Japan. See April +C newsletter. Most extrema agree to 5 or 6 decimal digits. +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns +FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections. +POWER FREQUENCY, 50.0, { Avoid warning message about suspicious XOPT, COPT. + .000200 .500 50. 50. + 1 1 0 { test } 1 -1 + 5 5 20 20 100 100 + A1 A2 0.1 1 + B1 B2 0.1 1 + C1 C2 0.1 1 +51A2 B2 .00528 571.5428 +52A3 -.029 3125.811 .1582 17115.12 +51B2 C2 A2 B2 +52B3 +51C2 A2 A2 B2 +52C3 + A2 31.416 + B2 31.416 + C2 31.416 + 1A3 A4 43.342 + 2B3 B4 12.546 43.342 + 3C3 C4 12.546 12.546 43.342 +BLANK card ending branch cards. +C The next card has T-open = 80 msec minus DELTAT/2 in order to avoid +C delayed opening (1/2 cycle) for PRIME and Burroughs. WSM. 27 FEB 1982. + A3 .01990 .0799 1 +BLANK card ending switch cards. +14A4 112.059 50. -20. -1. +14B4 112.059 50. -140. -1. +14C4 112.059 50. 100. -1. +C 59A1 11.3901 50.0 - 30.0 (before FIX SOURCE use) +$DEPOSIT, NOSM59=-1 { Special value will change Type-58 S.M. to Type-59 S.M. +C Note about preceding and following. Beginning 27 March 2000, this new +C variable NOSM59 of STARTUP can be used either to convert or to trap +C Type-58 or Type-59 S.M. type codes. Five values are recognized: +C NOSM59 = 0 ===> no change. Preserve the original meaning. +C 1 convert any Type-59 S.M. to Type-58 S.M. +C 2 trap any Type-59 S.M. (halt after all are counted). +C -1 convert any Type-58 S.M. to Type-59 S.M. +C -2 trap any Type-58 S.M. (halt after all are counted). +C So, the following machine really is Type-59. But 58 is used because ATP +C automatically will convert this to 59 as a result of preceding NOSM59 = -1 +C This preserves the old answers, which will remain unchanged. +58A1 11.0 50.0 -40.0 +58B1 +58C1 +C TOLERANCES 200. { Columns 11-20 EPSUPA = 200 } 20 { 51-60: NIOMAX=20 +C 14 October 2001, define NIOMIN of columns 41-50 so Watcom agrees with Mingw32 +C and Salford. This _does_ help. Value 8 provides better FC agreement than 5. +C Previously, Watcom was iterating a different number of times, & was producing +C a slightly different answer. So, we slow execution artificially: +TOLERANCES 200. 8 20 { 51-60: NIOMAX=20 +PARAMETER FITTING 1.0 + 1 1 2 1.0 1.0 150.0 13.8 600.0 600.0 720.0 +BLANK card for quadrature axis of machine. + 0.0014 0.175 1.85 1.76 0.2575 1.76 0.18 0.18 + 5.74757 0.051142 0.382609 0.197985 + 1 1.0 50.0 1.0 +BLANK card ending all (here, just one) mass cards. + 12 + 21 + 3 1 +BLANK card terminating all (here, three) S.M. output requests. + FINISH +BLANK card terminating all EMTP source cards. + A1 39.98755 -4.507399 10.0 12.0 -60. + B1 39.98755 -4.507399 10.0 12.0 -180. -120. + C1 39.98755 -4.507399 10.0 12.0 60.0 120. +C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER + 1 1000 20 1 0.00001 0.1 2.0 2 +C Max del-V: .7E-5 .7E-5 .6E-5 .6E-5 .5E-5 .5E-5 .5E-5 .4E-5 .4E-5 .4E-5 .3E-5 +C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 +C Exit the load flow iteration loop with counter NEKITE = 133. If no warning +C Row Node Name Voltage magnit Degrees Real power P Reactive power +C 4 8 A1 1.13901111E+01 -30.00022 3.99871502E+01 -4.50736252E+00 +C 5 9 B1 1.13901111E+01 -150.00022 3.99871502E+01 -4.50736252E+00 +C 6 10 C1 1.13901111E+01 89.99978 3.99871502E+01 -4.50736252E+00 +C Total network loss P-loss by summing injections = 8.018000511887E+00 +C Step Time C2 B2 A2 A3 C3 +C +C C1 A4 B4 C4 A3 +C TERRA +C +C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1 +C ID IQ I0 IF IKD +C +C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1 +C IB IC EFD MFORCE MANG +C 0 0.0 .0791881332 -9.29565074 9.216462609 102.5472108 -37.8516127 +C .4577064E-4 105.3010154 -85.8421743 -19.4588411 0.0 +C -6.82107489 5.325616419 .391002E-14 1.014337449 -.76501E-15 +C -5.68498411 -.791423626 -.176005017 .60792369 .887825386 +C 1 .2E-3 -.592823298 -8.9455915 9.538414794 103.3173764 -43.8432677 +C -.715886965 107.4997647 -81.1499746 -26.3497901 0.0 +C -6.82114583 5.325512789 -.4718E-14 1.014337757 .6231167E-5 +C -5.41035524 -1.23063667 -.176005017 .6079236529 .8878253793 + 1 { Request for the output of all node voltages +C Last step begin: 2500 0.5 .1013093411 -9.26844034 9.167130999 16.94173877 +C Last step begin: .0149766932 105.3010154 -85.8421743 -19.4588411 0.0 +C Last step begin: -6.73361875 5.375255555 -.34086E-14 1.004763786 -.001205227 +C Last step begin: -5.61469909 -.863326479 -.176005017 .6045428409 .889016932 +C Variable max : 11.6152661 11.51967968 12.3018062 204.4814649 119.2607541 +C 12.40458274 112.0153099 112.0562692 112.0480769 6.3761474 +C -2.00131682 9.101894661 .13086E-13 1.315246828 1.316347492 +C 12.08207638 7.893166343 -.176005017 .60792369 .8935781536 +C Times of max : .0938 .1282 .0806 .0802 .0934 +C .0938 .4812 .4878 .2344 .0252 +C .0782 .0708 .1114 .074 .0332 +C .033 .0938 0.0 0.0 .0226 + PRINTER PLOT + 19410. 0. 120. MACH 1ID MACH 1IQ { Plot limits: (-2.301, 0.910) +C Note: Former IG (next plot) now (after Prof. Juan Martinez's changes of +C October, 1989) becomes IKQ. The eddy-current coil of the Q axis +C formerly was coil "KQ", now is "G". WSM and JM, 13 October 90 + 19410. 0. 120. MACH 1IF MACH 1IKQ { Plot limits: (-0.425, 1.636) + 19410. 0. 120. BRANCH { Plot limits: (-1.874, 1.755) + A1 A2 B1 B2 C1 C2 +$DEPOSIT, NOSM59=0 { Cancel special value that was set earlier in this subcase +BLANK card terminating all plot cards. +BEGIN NEW DATA CASE +C 2nd of 6 data subcases. +C Same as preceding case, only with significant imbalance in the steady +C state due to change of generator at "A4" from -20 to -10 degrees. +C The load flow does not converge, and case is only run 5 steps. The +C idea is to force use of negative-sequence correction logic of SMINIT +C of overlay 11. No other test cases do this. Eventually, replace by +C a meaningful case. 2nd of 4 subcases (2nd of two machine cases). +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns +FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections. +POWER FREQUENCY, 50.0, { Avoid warning message about suspicious XOPT, COPT. + .000200 .001 50. 50. + 1 1 0 { test } 0 -1 + 5 5 20 20 100 100 + A1 A2 0.1 1 + B1 B2 0.1 1 + C1 C2 0.1 1 +51A2 B2 .00528 571.5428 +52A3 -.029 3125.811 .1582 17115.12 +51B2 C2 A2 B2 +52B3 +51C2 A2 A2 B2 +52C3 + A2 31.416 + B2 31.416 + C2 31.416 + 1A3 A4 43.342 + 2B3 B4 12.546 43.342 + 3C3 C4 12.546 12.546 43.342 +BLANK card ending branch cards. +C The next card has T-open = 80 msec minus DELTAT/2 in order to avoid +C delayed opening (1/2 cycle) for PRIME and Burroughs. WSM. 27 FEB 1982. + A3 .01990 .0799 1 +BLANK card ending switch cards. +14A4 112.059 50. -10. { Change -20 to -10 } -1. +14B4 112.059 50. -140. -1. +14C4 112.059 50. 100. -1. +C 59A1 11.3901 50.0 - 30.0 (before FIX SOURCE use) +59A1 11.0 50.0 -40.0 +59B1 +59C1 +TOLERANCES 200. { Columns 11-20 EPSUPA = 200 } 20 { 51-60: NIOMAX=20 +PARAMETER FITTING 1.0 + 1 1 2 1.0 1.0 150.0 13.8 600.0 600.0 720.0 +BLANK card for quadrature axis of machine. + 0.0014 0.175 1.85 1.76 0.2575 1.76 0.18 0.18 + 5.74757 0.051142 0.382609 0.197985 + 1 1.0 50.0 1.0 +BLANK card ending all (here, just one) mass cards. + 12 + 21 + 3 1 +BLANK card terminating all (here, three) S.M. output requests. + FINISH +BLANK card terminating all EMTP source cards. + A1 39.98755 -4.507399 10.0 12.0 -60. + B1 39.98755 -4.507399 10.0 12.0 -180. -120. + C1 39.98755 -4.507399 10.0 12.0 60.0 120. +C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER + 1 100 20 1 0.00001 0.1 2.0 2 +C Max del-V: .0016 .0015 .0013 .0011 .001 .9E-3 .8E-3 .6E-3 .5E-3 .4E-3 .4E-3 +C //. Warning! Divergence of the EMTP load flow has occurred. In the hope +C initialization continues, however. +C Exit the load flow iteration loop with counter NEKITE = 100. If no warning +C Row Node Name Voltage magnit Degrees Real power P Reactive power +C 4 8 A1 1.00000000E+01 -20.91363 3.45426843E+01 -2.24391750E+00 +C 5 9 B1 1.20000000E+01 -150.14225 4.22174603E+01 -6.71602291E+00 +C 6 10 C1 1.20000000E+01 90.73965 3.92028723E+01 -5.90634496E+00 +C Total network loss P-loss by summing injections = 7.625614104022E+00 + 1 { Request for the output of all node voltages +C Step Time C2 B2 A2 A3 C3 +C +C C1 A4 B4 C4 A3 +C TERRA +C +C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1 +C ID IQ I0 IF IKD +C +C MACH 1 MACH 1 MACH 1 MACH 1 MACH 1 +C IB IC EFD MFORCE MANG +C 0 0.0 -.048041917 -9.85266713 8.679836062 104.3486527 -36.7261188 +C -.154907255 110.356572 -85.8421743 -19.4588411 0.0 +C -6.92581041 5.259477223 -.445329E-4 .9859361614 .0320876279 +C -5.54502117 -1.06865338 -.170793823 .5995965889 .8894202603 +C 1 .2E-3 -1.16677225 -7.85747584 10.21232132 103.2738496 -48.1646884 +C -1.32143162 111.3606393 -81.1499746 -26.3497901 0.0 +C -6.87231021 5.289979997 .0014353419 .9858026654 .0272832017 +C -5.20835729 -1.54659375 -.170793823 .5996092398 .8894186749 +C Last step follows: +C 5 .001 -3.69592471 -6.2063594 10.61819433 99.27645368 -67.8903465 +C -4.03363622 110.9684495 -59.3822228 -52.6085139 0.0 +C -6.76271705 5.471445353 .0060095523 .9856033893 .017362355 +C -3.71587997 -3.37711507 -.170793823 .5996729546 .8894772231 + PRINTER PLOT +BLANK card terminating all plot cards. +BEGIN NEW DATA CASE +C 5-bus, 1-phase loadflow example from "Elements of Power System Analysis" +C by William D. Stevenson. McGraw-Hill Book Company, 2nd Edition (1962). +C The Gauss-Seidel solution is contained in Example 10.1, page 219 onward. +C Data is given in per unit, and we can use this directly (pretend that +C the base voltage, current, and power are all unity). We need XOPT = 60 +C since impedances are in ohms. Voltages could be multiplied by 1.414 to +C give the peak values required by the EMTP, but instead we use the VSCALE +C option that takes RMS (the square root of 2 is supplied internally) as +C input and produces RMS output (at least for the load flow solution table +C if not for EMTP branch flows). Bus number one, WHITE, is taken as the +C slack bus (the only bus without a power constraint). This is the 3rd of +C 5 subcases. It has no machine dynamics at all (unlike the preceding 2). +C This data case was 1st considered by Profs. Saul Goldberg & Bill Horton +C of Cal Poly in San Luis Obispo, California, for June, 1988, EMTP course. +C Answers change 10 February 1999 following the massive changes from +C TEPCO (Tokyo Electric Power Company) in Japan. See April newsletter. +C Convergence has been slowed substantially (when time, consider this). +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns +FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections. +$ERASE { Erase punched cards, created by preceding subcases, from punch buffer + .000200 .000 60. { T-max = 0 means that no transient solution follows + 1 1 1 0 1 + WHITE RED .10 .40 +C Stevenson's line between WHITE and BLUE appears on the final comment +C card of this comment block. For more generality, I have introduced an +C intermediate node DUMMY that is unconstrained. The solution of the +C network should be unchanged, of course. It is (original is ok, too). +C WHITE BLUE .15 .60 + WHITE DUMMY .15 .30 { 1st of 2 halves: all of R and 1/2 of X + DUMMY BLUE .30 { 2nd of 2 halves: zero R and 1/2 of X + WHITE YELLOW .05 .20 + RED GREEN .05 .20 + RED BLUE .10 .40 + GREEN YELLOW .05 .20 +BLANK card ending branch cards. +BLANK card ending switch cards. +C The following 2 cards would be used if peak rather than RMS input data: +C 14WHITE 1.4425 60. 0.0 -1. +C 14GREEN 1.4708 60. 2.0 -1. +14WHITE 1.02 60. 0.0 -1. +14GREEN 1.04 60. 2.0 -1. +C --------------+------------------------------ +C From bus name | Names of all adjacent busses. +C --------------+------------------------------ +C WHITE |RED *DUMMY *YELLOW* +C RED |WHITE *BLUE *GREEN * +C DUMMY |WHITE *BLUE * +C BLUE |RED *DUMMY * +C YELLOW |WHITE *GREEN * +C GREEN |RED *YELLOW* +C --------------+------------------------------ +BLANK card terminating all EMTP source cards. +C Next come power constraints of the load flow. There will be one +C for each non-slack generator. So, 3 of them will apply to nodes +C that have no Type-14 source as required by the algorithm. The +C program will define these internally. Yet, this is only possible +C for TMAX non-positive (no transient continuation). If data is +C modified to make TMAX > 0, a KILL = 40 error termination will +C will result. Yet, the transient simulation is possible as a +C 2nd simulation that replaces the power constraints at load nodes +C by constant-impedance loads. Branch cards for these will be +C punched by the present subcase, and the transient continuation +C will be illustrated by the following (4th of 4) subcase. +C The following 4 cards would be used if peak rather than RMS input data. +C Note that average of Vmin & Vmax = 0.5 ( 1.0 + 1.828428 ) = 1.414212 +C (power constraints RED, BLUE, and YELLOW), and 1.4708 = 1.414 * 1.02 +C RED -0.6 -0.3 1.01.828428 +C 1GREEN 1.0 1.4708 -20. 20. +C BLUE -0.4 -0.1 1.01.828428 +C YELLOW -0.6 -0.2 1.01.828428 +C 345678901234567890123456789012345678901234567890123456789012345678901234567890 + RED -0.6 -0.3 .85 1.15 + 1GREEN 1.0 1.02 -20. 20. + BLUE -0.4 -0.1 .85 1.15 + YELLOW -0.6 -0.2 .85 1.15 +C Max del-V: .0293 .0153 .0141 .0138 .0118 .0096 .0085 .0076 .0068 .0061 .0055 +C Source No. -2 -3 -4 -1 -3 -3 -3 -3 -3 -3 -3 +C Max del-V: .0023 .0021 .002 .0018 .0017 .0016 .0015 .0014 .0013 .0012 .0011 +C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 +C Max del-V: .6E-3 .5E-3 .5E-3 .5E-3 .4E-3 .4E-3 .4E-3 .4E-3 .4E-3 .3E-3 .3E-3 +C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 +C The following load-flow miscellaneous data card has two peculiarities. The +C use of VSCALE = 1.414 is the special flag requesting RMS rather than peak +C voltages. The use of KTAPER = 0 ensures constant acceleration factors +C (this works well for this problem whether RMS or peak values are used). +C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER + 1 500 20 1 .00001 0.0 0.0 1.414 0 +C Max del-V: .6E-6 .6E-6 .6E-6 .5E-6 +C Source No. -3 -3 -3 -3 +C Exit the load flow iteration loop with counter NEKITE = 164. If no warning +C Row Node Name Voltage magnit Degrees Real power P Reactive power +C 3 3 RED 9.41780959E-01 -3.81708 -5.99988408E-01 -2.99999500E-01 +C 2 7 GREEN 1.02000000E+00 2.43465 1.00001182E+00 3.82676335E-01 +C 4 5 BLUE 9.15105791E-01 -8.02053 -3.99985812E-01 -9.99988825E-02 +C 5 6 YELLOW 9.82613423E-01 -1.92907 -5.99992780E-01 -1.99999933E-01 +C Total network loss P-loss by summing injections = 5.282912652566E-02 +C Next, show the punched cards that $PUNCH (next) will display: +C $UNITS card should give X, C units XUNITS, CUNITS = 2.6526E-03 1.0000E-06 +C Following are branch cards for either loads or generator impedances: +C $VINTAGE, 1, +C Polar V, P, Q = 9.41780959E-01 -3.81708 -5.99988408E-01 -2.99999500E-01 +C RED .591308166 .295659302 +C Polar V, P, Q = 1.02000000E+00 2.43465 1.00001182E+00 3.82676335E-01 +C GREEN -.45374803 15276.6334 +C Polar V, P, Q = 9.15105791E-01 -8.02053 -3.99985812E-01 -9.99988825E-02 +C BLUE .985230493 .246313607 +C Polar V, P, Q = 9.82613423E-01 -1.92907 -5.99992780E-01 -1.99999933E-01 +C YELLOW .724153874 .241387449 +C $VINTAGE, -1, { Last of punched cardss that the following $PUNCH will display +$PUNCH { Flush punched cards: R-L branches for equivalent impedance loads +C Note about preceding card. Since no images anyway, remove. WSM. 15 Mar 98 +BLANK card ending requests for output variables +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 4th of 6 data subcases is a continuation of the 3rd into the time-step +C loop. The network differs in that constant-impedance loads (using the +C branch cards that were punched by execution of the 3rd subcase) replace +C the original power constraints at nodes without Type-14 sources. The +C time-step size DELTAT and end-time TMAX have been arbitrarily to +C selected to show that the transient continuation is smooth for just over +C a cycle. Note that this continuation involves a load flow, too, since +C otherwise angles of the Type-14 generators would not correspond to the +C network solution. The load flow could be omitted if the user were +C willing to re-key generator angles, but it generally is easier to let +C the computer solve the load flow a second time (this time with constant- +C impedance loads). +C Branch voltage output of the time-step loop, if declared by a 2-punch +C in column 80, was incompatible with load flow prior to 21 July 2002. +C So, to prove it now works, add 2-punch to 1st branch on this date. An +C older version of the program will simply ignore the column-80 punch. +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns +FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections. + .000200 .020 60. { T-max = 0 means that no transient solution follows + 1 1 1 0 1 -1 + 5 5 20 20 + WHITE RED .10 .40 2 + WHITE DUMMY .15 .30 { 1st of 2 halves: all of R and 1/2 of X + DUMMY BLUE .30 { 2nd of 2 halves: zero R and 1/2 of X + WHITE YELLOW .05 .20 + RED GREEN .05 .20 + RED BLUE .10 .40 + GREEN YELLOW .05 .20 +C <++++++> Cards punched by support routine on 20-Oct-90 05.01.52 <++++++> +$UNITS, 0.0, 0.0 { 1st of cards punched by preceding 3rd subcase: mHenry, microF +$VINTAGE, 1, { 2nd of cards punched by preceding 3rd subcase: use wide formt +C Polar V, P, Q = 9.41780959E-01 -3.81708 -5.99988408E-01 -2.99999500E-01 + RED 1.18261633E+00 1.56852131E+00 +C Polar V, P, Q = 9.15105791E-01 -8.02053 -3.99985812E-01 -9.99988825E-02 + BLUE 1.97046099E+00 1.30673427E+00 +C Polar V, P, Q = 9.82613423E-01 -1.92907 -5.99992780E-01 -1.99999933E-01 + YELLOW 1.44830775E+00 1.28060019E+00 +$UNITS, -1, -1 { Restore whatever XOPT and COPT existed before preceding $UNITS +$VINTAGE, 0, { Last of cards punched by preceding 3rd subcase +BLANK card ending branch cards. +BLANK card ending switch cards. +14WHITE 1.02 60. 0.0 -1. +14GREEN 1.04 60. 2.0 -1. +C --------------+------------------------------ +C From bus name | Names of all adjacent busses. +C --------------+------------------------------ +C WHITE |RED *DUMMY *YELLOW* +C RED |TERRA *WHITE *BLUE *GREEN * +C DUMMY |WHITE *BLUE * +C BLUE |TERRA *RED *DUMMY * +C YELLOW |TERRA *WHITE *GREEN * +C GREEN |RED *YELLOW* +C TERRA |RED *BLUE *YELLOW* +C --------------+------------------------------ +BLANK card terminating all EMTP source cards. +C Max del-V: .0282 .0015 .0013 .0011 1.E-3 .9E-3 .8E-3 .7E-3 .6E-3 .5E-3 .4E-3 +C Source No. -1 1 1 1 1 1 1 1 1 1 1 +C Max del-V: .1E-3 .9E-4 .8E-4 .7E-4 .6E-4 .5E-4 .5E-4 .4E-4 .4E-4 .3E-4 .3E-4 +C Source No. 1 1 1 1 1 1 1 1 1 1 1 +C Max del-V: .7E-5 .6E-5 .5E-5 .4E-5 .4E-5 .3E-5 .3E-5 .3E-5 .2E-5 .2E-5 .2E-5 +C Source No. 1 1 1 1 1 1 1 1 1 1 1 +C Exit the load flow iteration loop with counter NEKITE = 58. If no warning o +C Row Node Name Voltage magnit Degrees Real power P Reactive power +C 2 7 GREEN 1.02000000E+00 2.43435 9.99985664E-01 3.82679671E-01 +C Total network loss P-loss by summing injections = 1.652795946964E+00 +C 345678901234567890123456789012345678901234567890123456789012345678901234567890 + 1GREEN 1.0 1.02 -20. 20. +C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER + 1 500 20 1 .00001 0.0 0.0 1.414 0 +C Step Time YELLOW BLUE DUMMY RED WHITE +C +C 0 0.0 1.38883765 1.28149564 1.336017508 1.328924768 1.442497834 +C 1 .2E-3 1.388413609 1.291450977 1.337993527 1.331824678 1.438399553 +C 2 .4E-3 1.380101179 1.294070033 1.332367678 1.32715794 1.426128 +C 3 .6E-3 1.363947438 1.289337599 1.319171785 1.314950922 1.405752903 + 1 { Request for the output of all node voltages +C 100 .02 .4736844387 .567774976 .485772514 .4950043848 .4457563449 +C Variable max: 1.389411211 1.294070033 1.33815194 1.331824678 1.442497834 +C Times of max: .0168 .4E-3 .0168 .2E-3 0.0 +C Variable min: -1.38955118 -1.29326058 -1.33749084 -1.33106894 -1.44204228 +C Times of min: .0084 .0088 .0084 .0086 .0084 + PRINTER PLOT + 144 2. 0.0 20. RED { Axis limits: (-1.331, 1.332) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 5th of 6 data subcases is a variation of the 3rd and 4th. One line +C was omitted from the 5-bus network in order to make a radial system. +C Specifically, the connection (WHITE, BLUE) was removed, leaving just +C 4 lines connecting the 5 busses. I.e., there is no loop. The solution +C should not be much changed since loads of the two connecting nodes were +C increased or decreased by the known amount of the branch flow. Yet, a +C 2nd change must be mentioned: conversion of bus GREEN from a generator +C to a load. This was to make the present case more like NINEBUS from +C Jian Jiang of Clemson University in South Carolina. See EEUG list server +C mail dated October 7th. WSM, 10 October 2001. +FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections. + .000200 .000 60. { T-max = 0 means that no transient solution follows + 1 1 1 0 1 + WHITE YELLOW .05 .20 + YELLOWGREEN .05 .20 + GREEN RED .05 .20 + RED BLUE .10 .40 +BLANK card ending branch cards. +BLANK card ending switch cards. +14WHITE 1.02 60. 0.0 -1. +BLANK card terminating EMTP source cards. +C Next come power constraints of the load flow. There will be one +C for each non-slack generator: + RED -0.4 -.14 .85 1.15 + GREEN 1.0 0.3 .85 1.15 + BLUE -.15 .025 .85 1.15 + YELLOW -0.6 -0.2 .85 1.15 +C The following load-flow miscellaneous data card has two peculiarities. The +C use of VSCALE = 1.414 is the special flag requesting RMS rather than peak +C voltages. The use of KTAPER = 0 ensures constant acceleration factors +C (this works for this data). +C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER + 1 500 20 1 .00001 0.1 7.0 1.414 0 +C Max del-V: .0249 .0254 .0259 .0216 .0158 .0092 .0059 .006 .0039 .004 .0059 .0038 .0037 .0031 .002 .0025 .0022 .0028 .0018 .0015 +C Source No. 2 2 2 -3 2 -3 4 2 -2 -1 -2 -1 -2 -2 -1 -2 -1 -2 -1 -1 +C Max del-V: .0016 .0012 .0013 .0011 .0012 .001 1.E-3 .9E-3 .9E-3 .9E-3 .8E-3 .8E-3 .8E-3 .7E-3 .7E-3 .7E-3 .7E-3 .6E-3 .6E-3 .6E-3 +C Source No. -2 -1 -2 -1 -2 -3 -1 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 -3 ... (etc.) +C Max del-V: .2E-5 .2E-5 .2E-5 .2E-5 .2E-5 .1E-5 .1E-5 .1E-5 .1E-5 +C Source No. -3 -3 -3 -3 -3 -3 -3 -3 -3 +C Exit the load flow iteration loop with counter NEKITE = 169. If no warning on the preceding line, convergence was attained. +C Row Node Name Voltage magnit Degrees Real power P Reactive power +C 2 5 RED 9.61553517E-01 -3.15303 -4.00001772E-01 -1.39992436E-01 +C 3 4 GREEN 1.02208195E+00 2.93622 9.99997725E-01 3.00004959E-01 +C 4 6 BLUE 9.54082493E-01 -7.05955 -1.50001447E-01 2.50097905E-02 +C 5 3 YELLOW 9.83589266E-01 -1.67083 -6.00000908E-01 -1.99997137E-01 +BLANK card ending requests for output variables +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 6th of 6 data subcases is a variation of the 3rd. The solution is +C identical. But punched cards differ because T-stop of Type-14 source +C GREEN is nonzero. Previously, this data field was left blank. Nonzero +C T-stop is a flag that affects only the punched output. A nonzero value +C in columns 71-80 is taken to be internal source reactance. Instead of an +C equivalent impedance load, ATP will create an equivalent Type-14 voltage +C source behind this reactance. The A6 bus name will have "G" appended +C to it on the right to produce the name of the internal node. Thus GREEN +C yields GREENG as the internal node. Should the original name involve +C all 6 bytes, the last will be overlaid (e.g., PURPLE would be changed +C to PURPLG). This is for every Type-14 source except the one that will +C be used as the slack bus. The slack bus (here, WHITE) is the exception, +C and no output is created for it. This is why T-stop of WHITE has +C been left blank. Adding reactance to WHITE will not change the output. +C Let's summarize the change. Using X = 0.5 will change the following +C two 2 punched cards that follow the ruler: +C 12345678901234567890123456789012345678901234567890123456789012345678901234 +C Polar V, P, Q = 1.02000000E+00 2.43465 1.00001182E+00 3.82676335E-01 +C GREEN -.45374803 15276.6334 +C (one comment card followed by one wide-format R-C branch card) to: +C 14GREENG 1.84312934 60.24.5286405 -1.0 +C GREENGGREEN 0.0 0.5 +C (one Type-14 source card followed by one R-L branch card). +C Note that reactance X normally will be positive, and the R-L branch +C card carries this positive value in columns 43-58. This differs from +C the preceding equivalent impedance which might be either R-L or R-C +C depending on the power factor. Finally, note the raised voltage. Most +C of the increase from 1.02 to 1.84 is due to the need for peak rather +C than RMS on the Type-14 source card. As the reactance approaches zero, +C the voltage should approach 1.02 * 1.414 (square root of 2). This 6th +C subcase is being added 8 August 2007 following an inquiry by Prof. +C Mustafa Kizilcay of the University of Siegen in Germany. The use of +C T-stop > 0 was never documented in the ATP Rule Book. Neither was it +C illustrated in standard test case. That is why this subcase now is +C being added. WSM. +C Previous comment about the slack bus being missing is correct for normal +C data. But another undocumented option allows the user to name his slack +C bus manually at the end of power constraints. If this is done, there will +C be output for the slack generator, too. This option depends on the request +C "last is slack bus" which can appear in any case and any location to the +C right of the name field in columns 3-8. WSM. 9 August 2007 +C The preceding continues to be honored, but is believed to be clumsy. +C The more modern way to name the slack bus omits one of the 2 associated +C data cards by use of a special, in-line comment. The interpretation is +C improved at the same time. WSM. 15 August 2007 +C Final improvement was suggested by Prof. Kizilcay, who wondered why the +C user should be forced to sort the punched output manually. Prof. Kizilcay +C suggested that /SOURCE might precede each Type-14 source card & /BRANCH +C might follow it. Done. WSM. 15 August 2007. +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision for 8 output columns +FIX SOURCE { An EMTP load flow will satisfy requested phasor power injections. +$ERASE { Erase punched cards, created by preceding subcases, from punch buffer + .000200 .000 60. { T-max = 0 means that no transient solution follows + 1 1 1 0 1 + WHITE RED .10 .40 + WHITE DUMMY .15 .30 { 1st of 2 halves: all of R and 1/2 of X + DUMMY BLUE .30 { 2nd of 2 halves: zero R and 1/2 of X + WHITE YELLOW .05 .20 + RED GREEN .05 .20 + RED BLUE .10 .40 + GREEN YELLOW .05 .20 +BLANK card ending branch cards. +BLANK card ending switch cards. +14WHITE 1.02 60. 0.0 -1. 0.4 +14GREEN 1.04 60. 2.0 -1. 0.5 +BLANK card terminating all EMTP source cards. + RED -0.6 -0.3 .85 1.15 + 1GREEN 1.0 1.02 -20. 20. + BLUE -0.4 -0.1 .85 1.15 + YELLOW -0.6 -0.2 .85 1.15 + WHITE { Slack bus for punched output +C The preceding special in-line comment identifies the name of the slack bus. +C The older way, which continues to be honored, is to follow the slack bus name +C by a 2nd card which carries the special text (any case) "last is slack bus" +C NNNOUT NITERA NFLOUT NPRINT RALCHK CFITEV CFITEA VSCALE KTAPER + 1 500 20 1 .00001 0.0 0.0 1.414 0 +$PUNCH { Flush punched cards: R-L branches for equivalent impedance loads +BLANK card ending requests for output variables +BLANK card ending plot cards +BEGIN NEW DATA CASE +BLANK +EOF -- cgit v1.2.3