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/dc45.dat | 288 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 288 insertions(+) create mode 100644 benchmarks/dc45.dat (limited to 'benchmarks/dc45.dat') diff --git a/benchmarks/dc45.dat b/benchmarks/dc45.dat new file mode 100644 index 0000000..d06f9be --- /dev/null +++ b/benchmarks/dc45.dat @@ -0,0 +1,288 @@ +BEGIN NEW DATA CASE +C BENCHMARK DC-45 +C 1st half of split version of DC-33 . Rather than use TACS to integrate +C flux during the simulation, that operation is now posponed until the post- +C processing of compansion case DC-46. Connection is via LUNIT4 plot file. +C Total of 4 subcases. +C DISK PLOT DATA { Toggle the Apollo default of LUNIT4 = -4 to +4 (use disk) +C $CLOSE, UNIT=4 STATUS=DELETE { Destroy empty date/time plot file of "SYSDEP" +C $OPEN, UNIT=4 FILE=DC45TO46.pl4 ! { Expected to be C-like (for L4BYTE = 1 ) +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .000050 .050 + 1 1 1 1 1 -1 0 2 + 5 5 20 20 100 100 +96 GEN 8888. 1.E-9 1 + 1.0 -0.7 + 2.0 0.9 + 3.5 1.0 + 9999. +BLANK card ending branch cards +BLANK card ending (here non-existent) switch cards +14GEN 377. 60. -1. +C --------------+------------------------------ +C From bus name | Names of all adjacent busses +C --------------+------------------------------ +C GEN |TERRA *TERRA * +C TERRA |GEN *GEN * +C --------------+------------------------------ +BLANK card ending source cards +C Total network loss P-loss by summing injections = 7.106450000000E-06 +C ---- Initial flux of coil " " to "GEN " = 1.00000000E-09 +C Inject: GEN 377. 377. .377E-16 2.6786345332877 .710645E-14 +C Inject: 0.0 0.0 -2.678634533288 -90.0000000 504.92260952473 +C Step Time GEN TERRA +C GEN +C 0 0.0 377. 0.0 +C 1 .5E-4 376.9330268 -.523564601 +C 2 .1E-3 376.7321312 -1.04694318 +C 3 .15E-3 376.3973844 -1.02374791 +C 4 .2E-3 375.9289054 -1.04551661 +C 5 .25E-3 375.3268607 -1.06725433 +C 10 .5E-3 370.3222935 -1.17520854 + 1 +C 800 .04 -304.999407 1.070136441 +C 900 .045 -116.499407 2.765761443 +C 1000 .05 377. -1.4375 +C Variable maxima : 377. 3.499909256 +C Times of maxima : 0.0 .0125 +C Variable minima : -377. -3.49990923 +C Times of minima : .025 .0375 + PRINTER PLOT + 194 4. 0.0 20. GEN { Axis limits: (-3.500, 3.500) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 2nd of 4 subcases is unrelated to 1st (no hysteresis). But there is a +C nonlinear inductor (Type-93, true-nonlinear element). This example +C illustrates TACS control of the optional shunt current source that is +C declared by the second data field (columns 21-26) of the TACS CONTROL +C declaration. Since the controlling TACS signal is identically unity, +C this injects an extra amp of current into node "TRAN" beginning on +C the first time step. This would produce a discontinuity in the branch +C current (GEN, TRAN) at time zero, were it not for exact cancellation +C by a one-amp Type-11 current source that is connected to node "TRAN". +C It will be noted that the current (GEN, TRAN) continues smoothly out +C of the phasor solution. The current of the NL reactor (TRAN, TERRA) +C shows a one-amp jump beginning on the first step, due to the internal, +C TACS-controlled source (remember, the inductor current is continuous). +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .000200 .020 + 1 1 1 1 1 -1 0 2 + 5 5 10 10 +TACS HYBRID { We use TACS only to produce signal that controls shunt current +99SHUNT = 1.0 { Unit current (dc source) is inserted in parallel with ZnO +33SHUNT { Output only this one TACS variable that controls NL shunt current +77SHUNT 1.0 { Initial condition required for smooth electrical step 1 +BLANK card ending all TACS data + GEN TRAN 2.0 1 +93TRAN 1.0 1.0 1 + 0.0 0.0 + 0.9 0.9 + 2.0 1.1 + 10. 1.2 + 9999 + TACS CONTROL SHUNT { Use 2nd of 3 names: for TACS-controled shunt current +BLANK card ending branch cards +BLANK card ending (here non-existent) switch cards +14GEN 377. 60. -1. +11TRAN -1 -1.0 { 1-amp external dc current cancels internal NL source +BLANK card ending source cards +C Total network loss P-loss by summing injections = 1.000018973436E+00 +C +C GEN 377. 377. .00530514044263 1.000009486673 1.0000189734357 +C 0.0 0.0 -.9999954144499 -89.6960393 188.49913562381 +C +C ---- Initial flux of coil "TRAN " to " " = 5.30514044E-03 +C +C Step Time TRAN GEN TRAN GEN TACS +C TERRA TRAN SHUNT +C 0 0.0 376.9893897 377. .0053051404 .0053051404 1.0 +C 1 .2E-3 375.7677437 375.9289054 -.919419146 .0805808538 1.0 + 1 +C 100 .02 114.1251632 116.4994069 .1871218422 1.187121842 1.0 +C Variable max : 376.9893897 377. .4424771753 1.442477175 1.0 +C Times of max : 0.0 0.0 .0042 .0042 0.0 +C Variable min : -376.817325 -376.880939 -2.44858983 -1.44858983 1.0 +C Times of min : .0084 .0084 .0124 .0124 0.0 + PRINTER PLOT +C -0.653 -0.353 -0.054 0.245 0.545 0.844 1.143 1.442 +C ------+---------+---------+---------+---------+---------+---------+---------+ +C BBBBBBBBBBBBBBBBBBBBBBBBBBBB*AA +C BBBBB | AAAAAA +C BBBBBBB | AAAAAAAA +C BBBBB | AAAAA +C BBBB | AAAA +C BBBBBBBBBBBBBB AAAAAAAAAAAAAA +C | BBB AAA +C | BBBBBBBBBB AAAAAAAAA +C BBBBBB AAAAAAA +C BBBBB | AAAA +C BBBBB | AAAAA +C BBBBBBBB | AAAAAAAA +C BB |AAAAAA + 194 4. 0.0 20. GEN TRAN TRAN { Axis limits: (-2.449, 1.442) + 144 4. 0.0 20. TRAN { Axis limits: (-3.768, 3.770) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 3rd of 4 subcases is related to 2nd. There is the same nonlinear L(i) +C inductor (Type-93, true-nonlinear element). This example illustrates +C optional TACS control of the (i, PSI) characteristic by multiplicative +C scaling. The TACS multiplier of "i" is declared by the third data +C field (columns 27-32) of the "TACS CONTROL" card. For this example, +C the TACS variable is FACTOR, which varies smoothly (via one half a +C cosine cycle) from 2.0 to 1.0, thereby reducing the current that is +C drawn by approximately a factor of two. The series resistor is small, +C to it does not drop much voltage. As a result, the voltage, and hence +C the peak flux of the inductor, are nearly constant. Only the current +C increases by about a factor of two over the 10 msec "ramp" time. The +C presence of resistance merely reduces this ideal result: the factor is +C less than two --- maybe 1.5 or so. +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .000200 .030 + 1 1 1 1 1 -1 0 2 + 5 5 49 1 55 5 +TACS HYBRID { We use TACS only to produce signal that controls shunt current +99FACTOR = 1.5 + COSINE + STEP { Scaling factor changes from 3 to 1 over 0-10ms +14COSINE 0.5 50. { Smooth transition over 1st 10 msec } .010 +11STEP -0.5 50. { Hold COSINE value after 10 ms} .010 +33FACTORSTEP COSINE { Output only the 3 user-defined TACS variables +77FACTOR 2.0 { Initial condition required for smooth electrical step 1 +BLANK card ending all TACS data + GEN TRAN 20. 1 +93TRAN 2.0 1.0 1 + 0.0 0.0 + 0.9 0.9 + 2.0 1.1 { Execution does reach this 2nd seg of (i, PSI) + 9999 + TACS CONTROL FACTOR { Use 3rd of 3 names: TACS-scaled curve (x,y) +BLANK card ending branch cards +BLANK card ending (here non-existent) switch cards +14GEN 350. 50. -1.0 +BLANK card ending source cards +C Total network loss P-loss by summing injections = 4.885536661314E+01 +C ---- Initial flux of coil "TRAN " to " " = 1.39586762E-01 +C GEN 350. 350. .27917352350367 2.2103250125975 48.855366613143 +C 0.0 0.0 -2.19262372629 -82.7439171 383.7091521007 +C +C Step Time TRAN GEN TRAN GEN TACS +C TERRA TRAN FACTOR +C +C 0 0.0 344.4165295 350. .2791735235 .2791735235 2.0 +C 1 .2E-3 340.9842812 349.3093549 .4162536857 .4162536857 1.999013364 +C 2 .4E-3 336.2117309 347.2401455 .5514207281 .5514207281 1.996057351 +C +C 49 .0098 -347.838306 -349.309355 -.073552423 -.073552423 1.000986636 +C 50 .01 -347.141951 -350. -.14290245 -.14290245 1.0 +C 51 .0102 -345.0697 -349.309355 -.211982761 -.211982761 1.0 + 1 +C 150 .03 -347.789928 -350. -.110503595 -.110503595 1.0 +C Variable maxima : 349.3458152 350. 3.110373844 3.110373844 2.0 +C Times of maxima : .0198 0.0 .0042 .0042 0.0 +C Variable minima : -348.491847 -350. -2.24512833 -2.24512833 1.0 +C Times of minima : .0298 .01 .0146 .0146 .01 + PRINTER PLOT + 194 3. 0.0 30. GEN TRAN TACS FACTOR { Axis limits: (-2.245, 3.110) +C CALCOMP PLOT { Repeat this same plot, only in high-resolution mode +C 194 3. 0.0 30. GEN TRAN TACS FACTOR +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 4th of 4 subcases is related to 3rd. That nonlinear inductor L(i) +C (Type-93, true-nonlinear element) has been made linear by reduction of +C the flux-current characteristic to a single point. The result is a +C TACS-controlled linear inductor L(t). Except that here dynamics are +C involved (because L stores energy, and has history), this is the +C inductive equivalent of commonly-used TACS-controlled resistance R(t). +C Between time 0 and 100 msec, inductance is to be varied sinusoidally +C between 125 mH and 500 mH. Disconnected branches BEGIN and END +C correspond to these two limiting values, and they define assymptotes +C of the transition between the two. The graph is smooth, and it shows +C a smooth transition between the two steady-state assymptotes. +PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) + .000400 .120 + 1 1 1 1 1 -1 + 5 5 +TACS HYBRID { We use TACS only to produce signal that controls shunt current +99FACTOR = 2.5 + COSINE + STEP { Scale factor changes from 4 to 1 over 0-100 ms +14COSINE 1.5 5.0 { Smooth transition over 1st 100 msec } .100 +11STEP -1.5 { Hold COSINE valu after 100 ms} .100 +33FACTORSTEP COSINE { Output only the 3 user-defined TACS variables +77FACTOR 4.0 { Initial condition to make output signal smooth at t = 0 +77COSINE 1.5 { Initial condition to make output signal smooth at t = 0 +BLANK card ending all TACS data + GEN TRAN 35. { R of series R-L branch to ground } 1 +93TRAN 8.0 1.0 { i-steady and PSI-steady for phasors + 0.0 0.0 { Symmetric curve always begins at origin 0,0 + 2.0 1.0 { Single (i, PSI) point defines L b4 FACTOR + 9999 { Terminate the characteristic after just the one point + TACS CONTROL FACTOR { Use 3rd of 3 names: TACS-scaled curve (x,y) +C That is all we need for the TACS-controlled L of series R-L circuit. But +C for a nice plot, we want assymptotic curves. The first in in the steady +C state, corresponding to i-steady and PSI-steady of the Type-93 branch +C card. This is L = PSI / i = 1.0 / 8.0 = 125 mH. Note that this matches +C the initial signal. Initial FACTOR = 4.0 and the characteristic gives +C L = PSI / i = 1.0 / 2.0 = 500 mH. Dividing by FACTOR gives 125 mH as the +C initial inductance within the dT loop. So, show phasor flow: + GEN BEGIN 35. 125. { Branch flow is phasor solution +C That was for the start. Next, we want the termination, which corresponds +C to FACTOR = 1.0 so no scaling. I.e., L = 500 mH (the characteristic): + GEN END 35. 500. { Branch flo is assymptote for infinite t +BLANK card ending branch cards +C Switches are used only to measure the currents through the preceding two +C branches. These are the two assymptotic signals of following plot: + BEGIN MEASURING 1 + END MEASURING 1 +BLANK card ending switch cards +C --------------+------------------------------ +C From bus name | Names of all adjacent busses. +C --------------+------------------------------ +C GEN |TRAN *BEGIN *END * +C TRAN |TERRA *GEN * +C BEGIN |TERRA *GEN * +C END |TERRA *GEN * +C TERRA |TRAN *BEGIN *END * +C --------------+------------------------------ +14GEN 420. 50. -1.0 +C TRAN 234.06699294777 313.5412844237 5.3123716300638 7.9842632447061 -.677791157E-13 -.677791157E-13 +C 208.61634607684 41.7095865 -5.960467030767 -48.2904135 1251.698076461 1251.6980765 +C +C TERRA 0.0 0.0 -5.312371630064 7.9842632447061 0.0 +C 0.0 0.0 5.9604670307669 131.7095865 0.0 +C Total network loss P-loss by summing injections = 2.350389825256E+03 +C Output for steady-state phasor switch currents. +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C BEGIN 5.31237163E+00 -5.96046703E+00 7.98426324E+00 -48.2904 0.00000000E+00 0.00000000E+00 +C END 5.67589241E-01 -2.54733456E+00 2.60980285E+00 -77.4387 0.00000000E+00 0.00000000E+00 +BLANK card ending source cards + GEN TRAN { Output node voltages for the controlled branch L(t) only +C First 2 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 3 output variables are branch currents (flowing from the upper node to the lower node); +C Next 3 output variables belong to TACS (with "TACS" an internally-added upper name of pair). +C Step Time GEN TRAN BEGIN END GEN TACS TACS TACS +C TERRA TERRA TRAN FACTOR STEP COSINE +C *** Phasor I(0) = 5.3123716E+00 Switch "BEGIN " to " " closed in the steady-state. +C *** Phasor I(0) = 5.6758924E-01 Switch "END " to " " closed in the steady-state. +C 0 0.0 420. 234.06699 5.3123716 .56758924 5.3123716 4.0 0.0 1.5 +C 1 .4E-3 416.68817 206.10551 6.0166476 .88197067 6.0166476 3.9998816 0.0 1.4998816 +C 2 .8E-3 406.80493 174.8932 6.6262352 1.1824658 6.6260494 3.9995263 0.0 1.4995263 +C 3 .0012 390.50612 140.93094 7.1314997 1.464335 7.1307196 3.9989342 0.0 1.4989342 +BLANK card ending node voltage outputs +C 290 .116 129.78714 208.41029 -4.028677 -2.244816 -2.246376 1.0 -1.5 0.0 +C 295 .118 339.78714 376.14579 .78854743 -1.037462 -1.038819 1.0 -1.5 0.0 +C 300 .12 420. 400.22543 5.3045736 .56616677 .5649877 1.0 -1.5 0.0 +C Variable maxima : 420. 409.93775 7.9767416 2.6063662 7.9698488 4.0 0.0 1.5 +C Times of maxima : 0.0 .0992 .0028 .0044 .0028 0.0 0.0 0.0 +C Variable minima : -420. -409.6697 -7.973184 -2.60532 -7.845698 1.0 -1.5 -1.499882 +C Times of minima : .01 .1092 .1128 .1144 .0128 0.1 0.1 .0996 + CALCOMP PLOT +C Orlando Hevia in Santa Fe, Argentina, supplied the following one-line +C correction in E-mail dated 10 January 2007. He learned of the bad ATP data +C during experimentation with the new GNU F95 compiler, which choked on this: +C NOTE A MISTAKE IN DATA 120. IS ONE COLUMN TOO TO THE RIGHT +C DATA FM1180 / '( I2, 2I1, E3.0, 2E4.0, E5.0, E4.0, 4A6, 2A16 )' / +C -=---====----=====---- +C 19412. 0.0 120. BRANCH Vary linear L Inductor current + 19412. 0.0120. BRANCH Vary linear L Inductor current + BEGIN END GEN TRAN TACS FACTOR +BLANK card ending plot cards +BEGIN NEW DATA CASE +BLANK + -- cgit v1.2.3