diff options
Diffstat (limited to 'benchmarks/dc38.dat')
| -rw-r--r-- | benchmarks/dc38.dat | 885 |
1 files changed, 885 insertions, 0 deletions
diff --git a/benchmarks/dc38.dat b/benchmarks/dc38.dat new file mode 100644 index 0000000..7bed8fb --- /dev/null +++ b/benchmarks/dc38.dat @@ -0,0 +1,885 @@ +BEGIN NEW DATA CASE
+C BENCHMARK DC-38
+C ZnO simulation similar to DC-37, only here a 3-phase network is used.
+C The same arrester having characteristic i = 2500 * ( v / V-ref ) ** 26
+C is used, only here the coefficient has been cut in four (to COEF =625)
+C so that the column multiplier COL = 4.0 can be used: 4 * 625 = 2500.
+C Also, the usual, recommended (and more accurate) exponential modeling
+C (Type-92 nonlinear R(i) requested by "5555.") is only used for two of
+C the three phases. In order to illustrate the piecewise-linear alter-
+C native (requested by "4444."), such less-accurate modeling (for the
+C highly-nonlinear ZnO, anyway) has been placed in the 3rd phase ("c").
+C There are a total of 11 subcases.
+ZO, 20, , , , 0.9, ,{ To improve ZnO convergence,control Newton ZnO iteration
+ .000050 .020000
+ 1 1 1 0 1 -1
+ 5 5 20 1 30 5 50 50
+-1SENDA RECA .305515.8187.01210 200. 0 { 200-mile, constant-
+-2SENDB RECB .031991.5559.01937 200. 0 { parameter, 3-phase
+-3SENDC RECC { transmission line.
+92RECA 5555. { 1st card of 1st of 3 ZnO arresters
+C VREF VFLASH VZERO COL
+ 778000. -1.0 0.0 4.0
+C COEF EXPON VMIN
+ 625. 26. 0.5
+ 9999.
+92RECB RECA 5555. { Phase "b" ZnO is copy of "a"
+92RECC 4444. { Phase "c" ZnO is piecewise-linear
+C VREF VFLASH VZERO
+ 0.0 -1.0 0.0
+ 1.0 582400. { First point of i-v curve.
+ 2.0 590800. { Data is copied from DC-39
+ 5.0 599200. { which was used to create
+ 10. 604800. { the ZnO branch cards that
+ 20. 616000. { are used in phases "a" &
+ 50. 630000. { "b". But there is some
+ 100. 644000. { distortion due to the use
+ 200. 661920. { of linear rather than the
+ 500. 694400. { more accurate exponential
+ 1000. 721280. { modeling, of course.
+ 2000. 756000.
+ 3000. 778400. { Last point of i-v curve.
+ 9999. { Terminator for piecewise-linear characteristic
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 408000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 408000. 60. -120. { three-phase, sinusoidal excitation
+14SENDC 408000. 60. 120. { with no phasor solution.
+C --------------+------------------------------
+C From bus name | Names of all adjacent busses.
+C --------------+------------------------------
+C SENDA |RECA *
+C RECA |TERRA *SENDA *
+C SENDB |RECB *
+C RECB |TERRA *SENDB *
+C SENDC |RECC *
+C RECC |TERRA *SENDC *
+C TERRA |RECA *RECB *RECC *
+C --------------+------------------------------
+BLANK card follows the last source card
+C Step Time RECC RECB RECA SENDA SENDB SENDC
+C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+C 1 .5E-4 .47358E-13 .15692E-13 .15692E-13 407927.52 -197303.88 -210623.64
+C 2 .1E-3 -.4736E-13 -.1569E-13 -.1569E-13 407710.105 -190537.66 -217172.44
+C 3 .15E-3 .47358E-13 .15692E-13 .15692E-13 407347.832 -183703.75 -223644.08
+ 1
+C Last step: 400 .02 -601371.07 152342.824 295692.924 126078.934 273005.287
+C Variable maxima : 651691.033 676288.521 709562.656 407991.946 407999.105
+C Times of maxima : .00985 .00455 .00115 .01665 .00555
+C Variable minima : -669507.52 -663771. -717417.08 -407991.95 -407996.42
+C Times of minima : .00325 .01435 .0085 .00835 .0139
+C To appreciate the distortion that is involved in the use of piecewise-linear
+C representation for phase "c", I also show the result for exponential "c".
+C The following are derived from a simulation where RECC is a copy of RECA:
+C Last step: 400 .02 -600972.73 179505.6 299541.296 126078.934 273005.287
+C Variable maxima : 680201.783 671644.425 709538.839 407991.946 407999.105
+C Times of maxima : .0098 .00455 .00115 .01665 .00555
+C Variable minima : -704350.77 -664092.88 -718634.71 -407991.95 -407996.42
+C Times of minima : .00325 .01435 .00855 .00835 .0139
+ PRINTER PLOT
+ 144 3. 0.0 20. RECA { Axis limits: (-7.174, 7.096)
+ CALCOMP PLOT
+ 144 2. 0.0 20. RECB
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 2nd of 11 subcases. This one uses the same ZnO arrester as the second
+C of DC-37, only here the gap has been omitted by V-flash < 0. The line
+C is the same as the 1st subcase, too, although here we illustrate the
+C specialized request for modal output. The first six branches are very
+C large resistors that have been added to reserve outputs for this usage.
+STEP ZERO COUPLE { No reason for this, other than illustration of the feature
+MODE VOLTAGE OUTPUT
+ZO { Needed to restore default values that were upset by first subcase?
+ .000050 .020000
+ 1 1 1 0 1 -1
+ 5 5 20 1 30 5 50 50
+ SENDA 1.E18 { 1st of 6 high-R branches that serve } 1
+ SENDB 1.E18 { only to reserve output variables in } 1
+ SENDC 1.E18 { the output vector for modal voltages } 1
+ RECA 1.E18 1
+ RECB 1.E18 1
+ RECC 1.E18 { 6th of 6 high-R branches } 1
+92RECA 5555. 1
+C VREF VFLASH VZERO COL
+ 0.778000000000000E+06 -0.100000000000000E+03
+C COEF EXPON VMIN
+ 0.294795442961157E+05 0.265302624185338E+02 0.545050636122854E+00
+ 9999
+92RECB RECA 5555. { Phase "b" ZnO is copy of "a"
+92RECC RECA 5555. { Phase "c" ZnO is copy of "a"
+-1SENDA RECA .305515.8187.01210 200. 0
+-2SENDB RECB .031991.5559.01937 200. 0
+-3SENDC RECC
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 408000. 60. 0.0
+14SENDB 408000. 60. -120.
+14SENDC 408000. 60. 120.
+BLANK card follows the last source card
+C Step Time RECC RECB RECA SENDA SENDB SENDC
+C
+C RECA RECB RECC
+C TERRA TERRA TERRA
+C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+C 0.0 0.0 0.0
+C 1 .5E-4 0.0 0.0 0.0 407927.52 -197303.88 -210623.64
+C 0.0 0.0 0.0
+ 1
+C Last step: 400 .02 -592984.99 209476.019 234551.89 126078.934 273005.287
+C Last step: -86000.409 17731.3181 665443.85
+C Maxima: 639771.795 647447.415 658163.048 407991.946 407999.105 407996.421
+C Maxima: 181518.845 784118.005 883946.319
+C Associated times: .00985 .00455 .00115 .01665 .00555 .0111
+C Associated times: .00995 .01435 .00325
+ PRINTER PLOT
+ 194 1. 0.0 10. RECA { Axis limits: (-3.030, 3.485)
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 3rd of 11 subcases. This one uses the same ZnO arrester as subcase three
+C of DC-37 -- a single characteristic (no gap, since V-flash < 0) that
+C consists of two exponentials. 3-phase line, sources remain unchanged.
+C An important addition is a 4th nonlinear element, a TACS-controlled R(t)
+C that is practically disconnected from the 3 ZnO surge arresters and the
+C line that they terminate. But to test the logic, we couple the TACS-
+C controlled R(t) with the ZnO by means of the high resistance R = 1.E8.
+C The TACS control and electrical use is copied from the second subcase of
+C DC-22. Note the small EPSILN to ensure all 4 NL elements are coupled.
+PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 9 output columns)
+ .000050 .020000 1.E-10 { Small EPSILN to couple 2 subnetworks
+ 1 1 1 0 1 -1
+ 5 5 21 -1 30 -5 50 50
+TACS HYBRID { In a real case, arcs are on electric side, and equations in TACS
+99RESIS = 1.0 + SIN ( 300 * TIMEX ) { 1st R(t) signal -- constant + sine wave
+33RESIS { Output the only TACS variable: the resistance function R(t)
+77RESIS 1.0 { Initial condition on 1st R(t) insures smooth start
+BLANK card ending all TACS data
+ BUS1 BUS2 1.0 { Master copy of two 1-ohm resistors } 1
+ BUS2 BUS1 BUS2 { 2nd of 2 linear branches in second subnetwork
+91BUS2 TACS RESIS { R(t) controlled by TACS variable "RESIS" } 1
+ RECA BUS1 1.E+8 { Near-infinite R couples ZnO and TACS R(t)
+-1SENDA RECA .305515.8187.01210 200. 0
+-2SENDB RECB .031991.5559.01937 200. 0
+-3SENDC RECC
+92RECA 5555. 1
+C VREF VFLASH VZERO COL
+ 0.778000000000000E+06 -1.0
+C COEF EXPON VMIN
+ 0.505584788677197E+07 0.464199973324622E+02 0.632754084797274E+00
+ 0.122767153039007E+05 0.166775903445228E+02 0.816748018907843E+00
+ 9999
+92RECB RECA 5555. { Phase "b" ZnO is copy of "a"
+92RECC RECA 5555. { Phase "c" ZnO is copy of "a"
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 408000. 60. 0.0 { 1st of 3 sources for transmission
+14SENDB 408000. 60. -120. { line that is terminated by the ZnO
+14SENDC 408000. 60. 120.
+11BUS1 1.0 { 1-volt battery excites ladder network of TACS R(t)
+BLANK card follows the last source card
+ RECC RECB RECA SENDA BUS2
+C Note immediate voltage at RECC, RECB, RECA, due to near-infinite R coupling:
+C Step Time RECC RECB RECA SENDA BUS2
+C 0 0.0 0.0 0.0 0.0 0.0 0.0
+C 1 .5E-4 .141238E-5 .141238E-5 .426252E-5 407927.52 .333333333
+C 2 .1E-3 .141238E-5 .141238E-5 .426252E-5 407710.105 .334983437
+BLANK card terminating selective output variables
+C Last step: 400 .02 -600366.6 218884.325 220028.774 126078.934 .292739915
+C Variable maxima : 625666.843 631501.662 634878.442 407991.946 .399999647
+C Times of maxima : .00985 .00455 .00115 .01665 .0053
+C Variable minima : -633451.49 -613164.81 -630265.01 -407991.95 0.0
+C Times of minima : .0032 .01435 .0088 .00835 0.0
+ PRINTER PLOT
+ 144 3. 0.0 20. RECA { Axis limits: (-6.303, 6.349)
+ 194 3. 0.0 20. BUS2 TACS RESIS { Axis limits: (0.000, 2.000)
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 4th of 11 subcases is unrelated to the preceding, although it does use
+C a Type-91 TACS-defined R(t) as the preceding subcase does. But the
+C subject is quite different as should be explained in the April, 2003,
+C newsletter: corona modeling by TACS control of series R-L-C branches.
+C 27 December 2002, combine 3 separate, disconnected demonstrations of
+C TACS CONTROL of series R-L-C branches. The 3 disconnected subnetworks are:
+C 1) Series R-L with L fixed; R is ramped to a limiting value;
+C 2) Series R-L with R fixed; L is ramped to a limiting value;
+C 3) Series R-C with R fixed; C is is stepped (cut in half);
+C In each case, unit current at radian frequency 1.0 will be forced through
+C the branch, and voltage will be measured. There are 3 disconnected
+C networks, and there will be 3 screen plots to demonstrate reasonableness
+C of the answers.
+PRINTED NUMBER WIDTH, 10, 2, { Limited precision (not needed) & good separation
+TACS POCKET CALCULATOR { Required for use of IF-THEN-ELSE-ENDIF below
+ .10 20.0 { 200 steps over 3 cycles at radian frequency equal to unity
+ 1 1 1 1 1 -1
+ 5 5
+TACS HYBRID { TACS is required to define R of the series R-L branch
+C The first 2 of 3 problems each can use a discontinuity at T = 15 seconds:
+IF( TIMEX .LE. 15.0 ) THEN { If simulation time T does not exceed 15 sec:
+88OHMS = 0.5 + TIMEX / 2.5 { R is ramped linearly from 0.5 to 6.5 at end
+88HENRY = 0.5 + TIMEX / 10.0 { L increases linearly from 0.5 through 2.0
+ELSE { Alternatively (if simulation time T does exceed 15 sec):
+88OHMS = 6.5 { Limiting R in ohms for 15 or more seconds.
+88HENRY = 2.0 { Limiting L in Henries for 15 or more seconds.
+ENDIF { Terminate 5-line block that chooses among 2 formula for inductance HENRY
+C The 3rd of 3 problems requires discontinuity at T = 11 seconds:
+IF( TIMEX .LE. 11.05 ) THEN { If simulation time T is 11 sec or less:
+88FARAD = 2.0 { C is fixed for first 11 of 20 seconds of simulation
+ELSE { Alternatively (if simulation time T is 11 or more):
+88FARAD = 1.0 { Half the capacitance corresponds to switch being open
+ENDIF { Terminate 5-line block that chooses among 2 formulas for supplemental X1
+33HENRY { Bring out just 1 of 3 TACS signals to show it is not necessary
+77HENRY 0.5 { Initial condition on L(t) avoids jump from 0 on step 1
+77OHMS 0.5 { Initial condition on R(t) avoids jump from 0 on step 1
+77FARAD 2.0 { Initial condition on C(t) avoids jump from 0 on step 1
+BLANK card ending all TACS data
+C First comes the R-L test where R is varied and L is held fixed. We have
+C 3 signals of interest: a) old Type-91 model; b) new TACS CONTROL;
+C and c) limiting value (for large times, this agrees with a and b):
+ TYP91 COMP 1000. { Inductance of 1 Henry is fixed half
+91COMP TACS OHMS { TACS-defined R(t) is the variable half
+ COMP 1.E+7 { Leakage path avoids floating subnetwork
+ RAMPR 0.5 1000. { New modeling begins with R-L branch
+ TACS CONTROL OHMS { TACS signal "OHMS" will define R of series R-L
+C A TOLERANCE= tag could be added to any TACS CONTROL card such as
+C the preceding if the tolerance EPSRLC for the application of parameter
+C changes should be different from EPSILN of the miscellaneous data card.
+C Location is arbitrary, so typically will be to the right of column 44
+C (end of the 3rd of 3 TACS names). For example, TOLERANCE=1.E-5 will
+C serve to define EPSRLC = 1.E-5 In that case, any relative parameter
+C change in excess of this value will order re-triangularization whereas
+C any smaller change will not. For this data, there would be no change,
+C however, since all changes are large. dT is artificially large.
+C The phasor solution of the Type-91 branch is wrong because Type-91 content
+C is ignored prior the the dT loop. Using SSONLY of STARTUP, we can add
+C a branch that will correct this problem. The following branch will be
+C present only during the phasor solution; it will draw the current that
+C really should be going through the Type-91 branch. This will avoid a very
+C high voltage spike (e.g., 1.E7) at time zero. It also demonstrates that
+C use of SSONLY is compatible with TACS CONTROL of a series R-L-C:
+ COMP NAME PHASOR 0.5 { Branch that will be erased as dT loop begins
+ LIMIT 6.5 1000.
+C Next comes R-L test where L is varied and R is held fixed. We have 2
+C signals of interest: a) assumptote (for large T): b) new TACS CONTROL:
+ ASSYM 0.5 2000. { Assymptote (where variation will end)
+ RAMPL 0.5 500. { Branch to be varied begins at 1/2 Henry
+ TACS CONTROL HENRY { TACS signal will define L of series R-L
+C Finally (3rd of 3), we have an R-C test where R fixed and C is is stepped
+C to correspond exactly to electric network switching (breaker opening). The
+C answer seems believable; it agrees by eyeball with switching.
+ GEN CAP 1.0 { Inductance of 1 Henry is fixed half
+ CAP 1.0E6 { This capacitance always is used
+ CAP2 1.0E6 { This capacitance is switched
+ NEWRC 1.0 2.0E6 { For comparison, begin with R-C 1
+ TACS CONTROL FARAD { TACS defines C of series R-C
+BLANK card ending electric network branches
+ CAP CAP2 -1. 8.0 { Switch will open on current 0 at T = 11.0
+BLANK card ending switches
+C Each of the branches is to be driven by the same current source having
+C radian frequency equal to unity. I.e., 1 / frequency = 2 * Pi. Excite
+C the three networks in order:
+C 1) Series R-L with L fixed; R is ramped to a limiting value;
+14TYP91 -1 1.0 .1591549 -1.
+14RAMPR -1 1.0 .1591549 -1.
+14LIMIT -1 1.0 .1591549 -1.
+C 2) Series R-L with R fixed; L is ramped to a limiting value;
+14ASSYM -1 1.0 .1591549 -1.
+14RAMPL -1 1.0 .1591549 -1.
+C 3) Series R-C with R fixed; C is is stepped (cut in half);
+14GEN -1 1.0 .1591549 -1.
+14NEWRC -1 1.0 .1591549 -1.
+BLANK card ending electric network source cards.
+C Total network loss P-loss by summing injections = 5.249999987500E+00
+C Output for steady-state phasor switch currents.
+C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive
+C CAP CAP2 5.00000000E-01 0.00000000E+00 5.00000000E-01 0.0000 0.00000000E+00 -1.25000034E-01
+C Node voltage outputs will be grouped by network for easy visual comparison:
+C <---- Test a ---->< Test b ><---- Test c ---->
+ TYP91 RAMPR LIMIT COMP RAMPL ASSYM GEN NEWRC CAP CAP2
+C First 10 output variables are electric-network voltage differences (upper voltage minus lower voltage);
+C Next 1 output variables belong to TACS (with "TACS" an internally-added upper name of pair).
+C Step Time TYP91 RAMPR LIMIT COMP RAMPL ASSYM GEN NEWRC CAP CAP2 TACS
+C HENRY
+C *** Phasor I(0) = 5.0000000E-01 Switch "CAP " to "CAP2 " closed in the steady-state.
+C 0 0.0 0.5 0.5 6.5 0.5 0.5 0.5 1.0 1.0 0.0 0.0 0.5
+C 1 0.1 .3975854 .3975854 6.36761 .4975021 .4475438 .2976688 1.044879 1.044879 .0498751 .0498751 .51
+C 2 0.2 .330401 .3709362 6.171598 .5292359 .3953563 .0923634 1.079318 1.079318 .0992519 .0992519 .52
+C 3 0.3 .2583286 .2576778 5.913921 .5540951 .3218644 -.113865 1.102973 1.102973 .147637 .147637 .53
+C *** Open switch "CAP " to "CAP2 " after 1.10000000E+01 sec.
+BLANK card ending names of nodes for node voltage output
+C 200 20. 1.738859 1.773888 1.738861 2.652564 -1.6153 -1.62337 1.819737 1.833053 1.41165 -.499468 2.0
+C Variable maxima: 6.576562 6.580848 6.576566 6.492014 2.065211 2.062599 1.912249 1.925565 1.497945 .4993704 2.0
+C Times of maxima: 18.7 18.6 18.7 18.8 17.6 17.5 19.6 19.6 14.1 1.6 15.1
+C Variable minima: -6.57 -6.60161 -6.57619 -6.49979 -2.00396 -2.06317 -1.11783 -1.11783 -.499579 -.499579 0.5
+C Times of minima: 15.6 15.5 3.0 15.7 14.4 8.1 3.6 3.6 11. 11. 0.0
+ CALCOMP PLOT { Switch to screen plot from printer plot of preceding subcase
+C 1) Series R-L with L fixed; R is ramped to a limiting value;
+C Plot the 3 branch voltages that result from 1 amp of current being driven
+C through each branch. Note TYP91 should lie on top of with RAMPR, and
+C this common signal should be close to the limiting value LIMIT for times
+C in excess of 12 seconds :
+ 143 2. 0.0 20. TYP91 RAMPR LIMIT Ramp R of R-L
+C 2) Series R-L with R fixed; L is ramped to a limiting value;
+C Plot the 2 branch voltages that result from 1 amp of current being driven
+C through each branch. Note RAMPL should be close to the limiting value
+C ASSYM for times in excess of 12 seconds :
+ 143 2. 0.0 20. RAMPL ASSYM Ramp L of R-L
+C 3) Series R-C with R fixed; C is is stepped (cut in half);
+C Plot the 2 branch voltages that result from 1 amp of current being driven
+C through each branch. Note NEWRC should agree with GEN for all time.
+C Following removal of capacitance, the curves are offset significantly:
+ 143 2. 0.0 20. NEWRC GEN Step C of R-C
+BLANK card ending plot cards
+BEGIN NEW DATA CASE
+C 5th of 11 subcases illustrates a practical (although oversimplied)
+C application of the preceding. Data comes from Orlando Hevia of UTN
+C in Santa Fe, Argentina, as originally named TACSCAPA.DAT Data is
+C being added to this test case on 30 December 2002. Whereas the first
+C such example from Mr. Hevia involved 200 cascaded line sections, this
+C more manageable illustration involves just 2. TACS is used to vary the
+C shunt capacitance of the line as an approximation to corona. Note that
+C comment cards below are machine-produced (Mr. Hevia seems to have a
+C separate program to create such cascaded data automatically). Numerical
+C burden of the simulation has been reduced by a factor of 20 without much
+C loss to the plot or extrema. A factor of 2 was gained by shortening
+C the simulation from 20 to 10 usec, and a factor of 10 was gained by
+C increasing the time step from the original 5 nanoseconds (5.E-9 sec).
+C The surge (lightning) is fast, so very high frequencies are involved.
+C Note Mr. Hevia's use of JMARTI frequency-dependent line modeling.
+C Warning. 7 September 2003, the answer changes substantially following
+C the correction of an error in OVER12 (introduce new variable N7).
+PRINTED NUMBER WIDTH, 10, 2, { Limited precision (not needed) & good separation
+TACS POCKET CALCULATOR OFF { End use of pocket calculator (preceding subcase)
+C The preceding probably is necessary because of complex definition of VAR002
+ 5.0E-08 20.E-06 { Hevia's dT increased by a factor of 10; cut Tmax in half
+ 1 1 1 0 1 -1
+ 5 5 10 10 134 1 170 10
+C $INCLUDE, CORONA1.PCH
+C FIRST STEP CAPACITY 3.000000E-06 uF/KM
+C SLOPE 3.000000E-12 uF/KVKM
+C CORONA INCEPTION VOLTAGE 3.600000E+05 V
+C LENGTH OF LINE SEGMENT 1.000000E+00 KM
+C NUMBER OF SEGMENTS 2.000000E+00
+TACS HYBRID
+90BUS002
+88DER00259+BUS002
+C DV/DT MUST BE POSITIVE, BUT THIS TEST PRODUCES
+C OSCILLATIONS ON CAPACITANCE
+88VAR002 = BUS002 .GT. 360000.00 { .AND. DER002 .GT. 0.0
+88CAP002 = 1.0E-08+VAR002*((BUS002- 360000.00)*0.3000E-11+0.3000E-05)
+33CAP002DER002BUS002VAR002 { Output all TACS signals including control C(t)
+BLANK card ending TACS data
+-1BUS000BUS001 2. 0.00 -2 1
+ 14 3.9461680140762559000E+02
+ 7.68954468040036890E+02 1.09493340867763940E+03 2.77331232270879630E+03
+ 1.24494695098279860E+04 4.87585677225587210E+04 1.94958822722845510E+05
+ 7.82012894548635460E+05 3.09109899381158690E+06 1.48401963798197680E+07
+ 3.34339104652340860E+07 1.56456366517231150E+07 4.10038300055303800E+07
+ 2.60359793110293930E+07 4.14639643816612140E+07
+ 6.61711924983759210E+00 1.43260235003813180E+01 1.39885566693366850E+02
+ 6.43953575180861780E+02 2.62156097340268890E+03 1.08866830412747530E+04
+ 4.53734562567173710E+04 1.87083684125800150E+05 9.33229189322630060E+05
+ 4.32016631012824080E+06 8.22729460640732390E+06 2.11696813048871940E+07
+ 1.36911740150641220E+07 2.35107210671712680E+07
+ 15 3.3528019962850977000E-06
+ 1.48107642189314750E+01 8.18386897856797330E+01 1.07718234528722760E+02
+ 1.39846901178167800E+02 1.72162896702735680E+02 2.28340646958654700E+02
+ 3.44414362842715720E+02 1.63690212466734790E+04 8.08859081632825200E+03
+ 5.81880629665730960E+04 8.57432646874608240E+05 5.25522742751047830E+05
+ 3.86658063350409460E+06 -1.10174538112164120E+07 1.31615353212200510E+07
+ 7.09885572628100910E+03 3.82889067640842040E+04 5.07666709936286210E+04
+ 6.42506894830861860E+04 7.75796172424984980E+04 1.09484865717845850E+05
+ 8.28324070221879670E+04 4.32415798449636150E+05 4.46749849274677110E+05
+ 8.54778751513758670E+05 3.29015425966867800E+06 4.21774707623910620E+06
+ 1.25365993856300990E+07 2.57506236853497770E+07 2.21885826483384670E+07
+ 1.00000000
+ 0.00000000
+-1BUS001BUS002BUS000BUS001
+C
+C THE OLD FILE HAD THE CAPACITANCE IN AN ISOLATED BUS!
+C
+ BUS002CAP002 0.1 { Capacitance is to be made voltage-dependent
+ TACS CONTROL CAP002 TOLERANCE=1.E-2
+C Note preceding card includes optional definition of the tolerance for use
+C of the TACS signal CAP002. Without this declaration, EPSRLC = EPSILN =
+C 1.E-8, and this results in 159 triangularizations to [Y] as seen in case-
+C summary statistics when KOMPAR = 0 (see STARTUP):
+C Size List 5. Storage for [Y] and triangularized [Y]. No. times = 159 ...
+C Using 1.E-3, this is reduced slightly to 148. This is the effect of not
+C making a change if the change to C is less than 1/10 of 1%. This ignores 11
+C of the 159 changes. Using 1.E-2, the "No. times" drops to 29; and using
+C 1/10, it drops to 6. So 1.E-2 is practical. Using 29 steps to approximate
+C C(t) should be plenty good (see plot of C). Yet 29 of 159 is only 18%, so
+C simulation is a lot faster (82% of triangularization is avoided).
+ BUSXXXBUS000 0.0001
+ BUSXXX 394.61
+ BUS002 394.61
+BLANK
+C TACS CONTROLLED SWITCH TO CONNECT/DISCONNECT THE SOURCE
+13CAP002CP1002 VAR002
+BLANK
+C DC SOURCE
+11CP1002 360000.0
+15BUSXXX 9.0 USRFUN { Hevia's own user-supplied so
+C Recall USRFUN sources are a family of user-supplied sources as first
+C described in the October, 2002, newsletter. Alternative sources that might
+C interest the reader include the following:
+C 15BUSXXX 8.0 usrfun
+C 15BUSXXX 1 0.3E-6 7.00E-6 10.01.000E06Heidler in-line 5
+C 15BUSXXX-1 0.3E-6 7.0E-6 30.05.000E03Heidler in-line 5
+C 15BUSXXX 1 1.2E-6 10.0E-6 10.01.000E06TWO EXP in-line
+BLANK card ending electric-network source cards
+ BUSXXXBUS000BUS001BUS002 { List of nodes for node-voltage output
+C First 4 output variables are electric-network voltage differences (upper voltage minus lower voltage);
+C Next 4 output variables belong to TACS (with "TACS" an internally-added upper name of pair).
+C Step Time BUSXXX BUS000 BUS001 BUS002 TACS TACS TACS TACS
+C CAP002 DER002 BUS002 VAR002
+C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+C 1 .5E-7 144899.2 72286.18 0.0 0.0 .1E-7 0.0 0.0 0.0
+C 2 .1E-6 523890.1 334255.7 0.0 0.0 .1E-7 0.0 0.0 0.0
+C 3 .15E-6 922280. 724128.3 0.0 0.0 .1E-7 0.0 0.0 0.0
+C 4 .2E-6 .1199E7 .10635E7 0.0 0.0 .1E-7 0.0 0.0 0.0
+C 5 .25E-6 .13008E7 .1254E7 0.0 0.0 .1E-7 0.0 0.0 0.0
+BLANK card ending output variable requests
+C 400 .2E-4 83987.68 85044.93 222821.1 191231.3 .1E-7 -.237E11 191231.3 0.0
+C Extrema of output variables follow. Order and column positioning are the same as for the preceding time-step loop output.
+C Variable maxima : .13008E7 .128E7 839549.4 588280.4 .3695E-5 .2214E13 588280.4 1.0
+C Times of maxima : .25E-6 .3E-6 .375E-5 .815E-5 .815E-5 .7E-5 .815E-5 .705E-5
+C Variable minima : 0.0 0.0 0.0 0.0 0.0 -.567E11 0.0 0.0
+C Times of minima : 0.0 0.0 0.0 0.0 0.0 .149E-4 0.0 0.0
+C 145 2. 0.0 20. BUS000BUS001BUS002 { Not enough space for Y-max
+C Replace the preceding normal plot card by following alternative wide format:
+ 145 BUS000BUS001BUS002 Voltage on line Volts
+C Zero units/inch in columns 5-7 means that another card carries the info:
+C Units/inch Minimum time Maximum time Bottom Y-axis Top of Y-axis
+ 2.0 0.0 20.0 0.0 1.4E6
+ 195 2. 0.0 20. TACS CAP002 Capacitance C(t)Farads
+BLANK card ending plot cards
+BEGIN NEW DATA CASE
+C 6th of 11 subcases illustrates a practical (although oversimplied)
+C application of the preceding. Data comes from Orlando Hevia of UTN
+C in Santa Fe, Argentina, as originally named TIDDHHC.DAT Data is
+C being added to this test case on 10 September 2003.
+C A SAMPLE OF CORONA WITH TACS CONTROLLED CAPACITORS
+C THE OUTPUT LOOKS BELIEVABLE
+C AN AVERAGE CAPACITANCE IS CALCULATED BETWEEN TIME STEPS
+PRINTED NUMBER WIDTH, 11, 1, { Restore default settings as if no declaration
+ 2.0E-08 40.E-06 { Orlando used dT = 1.E-8 for more realistic looking plots
+ 1 1 0 0 1 -1
+ 5 5 20 20 100 100 500 500
+TACS HYBRID
+90BUS001
+90BUS002
+88DER00159+BUS001
+88DER00259+BUS002
+88VOLTA1 = BUS001.GT.270000.0
+88VOLTA2 = BUS002.GT.270000.0
+88DELTA1 = (BUS001-270000.0)*1.0E-5
+88DELTA2 = (BUS002-270000.0)*1.0E-5
+88CAP011 = 1.0E-08+(DER001.GT.0.0)*1.0E-8*VOLTA1*DELTA1
+88CAP02153+CAP011 1.0E-8
+88CAP001 =(CAP021+CAP011)/2.0
+88CAP012 = 1.0E-08+(DER002.GT.0.0)*1.0E-8*VOLTA2*DELTA2
+88CAP02253+CAP012 1.0E-8
+88CAP002 =(CAP022+CAP012)/2.0
+33CAP001CAP002
+BLANK
+$VINTAGE, 1
+-1BUS000BUS001 7.88076E+01 4.80104E+02 2.93720E+05 1.00000E+00 1
+$VINTAGE, -1,
+-1BUS001BUS002BUS000BUS001
+-1BUS002BUS003BUS000BUS001
+C TACS CONTROLLED CAPACITANCES
+ BUS001 10.0 1
+ TACS CONTROL CAP001
+ BUS002 10.0 1
+ TACS CONTROL CAP002
+ BUSXXXBUS000 1.0 1
+ BUS003 468.82
+C
+C CONSTANT CAPACITANCES
+C
+-1VUS000VUS001BUS000BUS001
+-1VUS001VUS002BUS000BUS001
+-1VUS002VUS003BUS000BUS001
+ VUS003 468.82
+C
+ BUSXXXVUS000 1.0 1
+C CONSTANT CAPACITANCES
+ VUS001 10.0 1.0E-2 1
+ VUS002 10.0 1.0E-2 1
+C
+C NO CAPACITANCES
+C
+C CONSTANT CAPACITANCES
+C
+-1XUS000XUS001BUS000BUS001
+-1XUS001XUS002BUS000BUS001
+-1XUS002XUS003BUS000BUS001
+ XUS003 468.82
+C
+ BUSXXXXUS000 1.0 1
+C
+BLANK
+BLANK
+C ------==--------========--------========
+15BUSXXX 1 1.0E-6 7.0E-6 0.01.770E06TWO EXP in-line
+C ------==----------==========----------
+BLANK
+ BUS000
+ BUS001
+ BUS002
+ BUS003
+ VUS000
+ VUS001
+ VUS002
+ VUS003
+ XUS000
+ XUS001
+ XUS002
+ XUS003
+C Step Time BUS000 BUS001 BUS002 BUS003 VUS000 VUS001 VUS002 VUS003 XUS000 XUS001
+C
+C
+C XUS002 XUS003 BUS001 BUS002 BUSXXX BUSXXX VUS001 VUS002 BUSXXX TACS
+C TERRA TERRA BUS000 VUS000 TERRA TERRA XUS000 CAP001
+C
+C TACS
+C CAP002
+C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+C 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+C 0.0
+C 1 .2E-7 133766.113 0.0 0.0 0.0 133766.113 0.0 0.0 0.0 133766.113 0.0
+C 0.0 0.0 0.0 0.0 267.636123 267.636123 0.0 0.0 267.636123 .75E-8
+C .75E-8
+C 2 .4E-7 258225.657 0.0 0.0 0.0 258225.657 0.0 0.0 0.0 258225.657 0.0
+C 0.0 0.0 0.0 0.0 516.651878 516.651878 0.0 0.0 516.651878 .1E-7
+C .1E-7
+BLANK card ending nodes for node-voltage outputs
+C 2000 .4E-4 17081.5133 -23694.238 -37514.031 75051.7774 16473.2675 -93546.663 -101202.23 51712.6484 16072.4173 18461.1635
+C 21456.0213 25168.3439 -430.45364 -154.33355 -1000.651 -392.40513 507.97555 -115.34624 8.44498556 .1E-7
+C .1E-7
+C Variable maxima : .1766463E7 814156.947 655233.802 577911.165 .1766463E7 .1131763E7 .1118349E7 967992.864 .1766463E7 .1629911E7
+C .15039E7 .1343678E7 6074.55465 4854.86191 5930.26754 5930.26754 5072.07156 1967.01281 3534.29804 .640067E-7
+C .479668E-7
+C Times of maxima : .1E-5 .1322E-4 .2316E-4 .2692E-4 .1E-5 .1142E-4 .1676E-4 .2012E-4 .1E-5 .44E-5
+C .78E-5 .1118E-4 .132E-4 .1668E-4 .742E-5 .742E-5 .404E-5 .908E-5 .1E-5 .1322E-4
+C .2304E-4
+C Variable minima : 0.0 -23694.238 -38329.779 0.0 0.0 -166322.75 -101202.23 0.0 0.0 0.0
+C 0.0 0.0 -3733.6936 -2735.0324 -1337.4272 -1647.1189 -2170.2496 -1613.5916 -331.24485 0.0
+C 0.0
+C Times of minima : 0.0 .4E-4 .3998E-4 0.0 0.0 .3744E-4 .4E-4 0.0 0.0 0.0
+C 0.0 0.0 .14E-4 .2642E-4 .3788E-4 .3112E-4 .1392E-4 .2462E-4 .205E-4 0.0
+C 0.0
+ 145 2. 0.0 40. 0.02.E6VUS000VUS001VUS002VUS003 Constant C
+ 145 2. 0.0 40. 0.02.E6BUS000BUS001BUS002BUS003 TACS CONTROL
+ 145 2. 0.0 40. 0.02.E6BUS001BUS002VUS001VUS002 Both
+ 145 2. 0.0 40. 0.02.E6XUS000XUS001XUS002XUS003 No capacitors
+ 195 2. 0.0 40.-2.E36.E3BUSXXXBUS000BUSXXXVUS000 Currents
+BLANK card ending plot cards
+BEGIN NEW DATA CASE
+C 7th of 11 subcases illustrates a true delta connection of nonlinear
+C elements that use compensation. Prior to November of 2006, ATP would
+C have halted with a complaint that the Thevenin impedance matrix [Z-thev]
+C was singular as follows:
+C KILL code number Overlay number Nearby statement number
+C 209 18 3471
+C KILL = 209. ZnO solution by Newton`s method of 3 coupled ...
+C Order is critical. For the delta to be recognized, the 3 N.L. elements
+C must be contiguous and must have triplets of (BUS1, BUS2) names ordered
+C as NAMEA to NAMEB first, then NAMEB to NAMEC 2nd, and finally
+C NAMEC to NAMEA. Data appended 15 December 2006. WSM.
+PRINTED NUMBER WIDTH, 11, 2, { Deliberately reduce 9 output columns by 1 digit
+ZO, 20, , , , 0.9, ,{ To improve ZnO convergence,control Newton ZnO iteration
+ .000050 .020000
+ 1 1 1 0 1 -1
+ 5 5 20 1 30 5 50 50
+-1SENDA RECA .305515.8187.01210 200. 0 { 200-mile, constant-
+-2SENDB RECB .031991.5559.01937 200. 0 { parameter, 3-phase
+-3SENDC RECC { transmission line.
+92RECA RECB 5555. { 1st card of 1st of 3 ZnO arrest} 3
+C VREF VFLASH VZERO COL
+ 778000. -1.0 0.0 4.0
+C COEF EXPON VMIN
+ 625. 26. 0.5
+ 9999.
+92RECB RECC RECA RECB 5555. { Phase "bc" ZnO is copy of "ab" } 3
+92RECC RECA RECA RECB 5555. { Phase "ca" ZnO is copy of "ab" } 3
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 236000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 236000. 60. -120. { three-phase, sinusoidal excitation
+14SENDC 236000. 60. 120. { with no phasor solution.
+BLANK card follows the last source card
+BLANK card ending node voltage outputs
+ PRINTER PLOT
+ 194 2. 0.0 20. BRANCH { Axis limits (-1.829, 0.525)
+ RECA RECB RECB RECC RECC RECA
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 8th of 11 subcases unites the 1st with the 7th. Both the Y & the delta
+C connections are present with the Y of the 1st subcase having node names
+C as follows: SEND ---> LINE REC ---> END The two subnetworks
+C are physically disconnected but mathematically coupled by one very high
+C resistance branch (see comment cards) that makes the difference between
+C two 3x3 matrices [Z-thev] and one 6x6 matrix. See (RECA, ENDA). Data
+C is added 15 December 2006. WSM.
+PRINTED NUMBER WIDTH, 11, 2, { Deliberately reduce 9 output columns by 1 digit
+ZO, 20, , , , 0.9, ,{ To improve ZnO convergence,control Newton ZnO iteration
+ .000050 .020000
+ 1 1 1 0 1 -1
+ 5 5 20 1 30 5 50 50
+C Begin with branches of the 1st subcase:
+-1LINEA ENDA .305515.8187.01210 200. 0 { 200-mile, constant-
+-2LINEB ENDB .031991.5559.01937 200. 0 { parameter, 3-phase
+-3LINEC ENDC { transmission line.
+92ENDA 5555. { 1st card of 1st of 3 ZnO arresters
+C VREF VFLASH VZERO COL
+ 778000. -1.0 0.0 4.0
+C COEF EXPON VMIN
+ 625. 26. 0.5
+ 9999.
+92ENDB ENDA 5555. { Phase "b" ZnO is copy of "a"
+92ENDC 4444. { Phase "c" ZnO is piecewise-linear
+C VREF VFLASH VZERO
+ 0.0 -1.0 0.0
+ 1.0 582400. { First point of i-v curve.
+ 2.0 590800. { Data is copied from DC-39
+ 5.0 599200. { which was used to create
+ 10. 604800. { the ZnO branch cards that
+ 20. 616000. { are used in phases "a" &
+ 50. 630000. { "b". But there is some
+ 100. 644000. { distortion due to the use
+ 200. 661920. { of linear rather than the
+ 500. 694400. { more accurate exponential
+ 1000. 721280. { modeling, of course.
+ 2000. 756000.
+ 3000. 778400. { Last point of i-v curve.
+ 9999. { Terminator for piecewise-linear characteristic
+C Done with branches of the 1st subcase; follow by branches of 7th subcase:
+-1SENDA RECA .305515.8187.01210 200. 0 { 200-mile, constant-
+-2SENDB RECB .031991.5559.01937 200. 0 { parameter, 3-phase
+-3SENDC RECC { transmission line.
+92RECA RECB 5555. { 1st card of 1st of 3 ZnO arrest} 3
+C VREF VFLASH VZERO COL
+ 778000. -1.0 0.0 4.0
+C COEF EXPON VMIN
+ 625. 26. 0.5
+ 9999.
+92RECB RECC RECA RECB 5555. { Phase "bc" ZnO is copy of "ab" } 3
+92RECC RECA RECA RECB 5555. { Phase "ca" ZnO is copy of "ab" } 3
+C Remove the following large resistance to solve each 3-phase bank of surge
+C arresters separately. With this branch present, the 6 N.L. elements all
+C are in the same subnetwork, so 6 N.L. equations in 6 unknowns are solved
+C by Newton's method at each time step. Without the branch, there will be
+C two sequential solutions of 3 N.L. equations each. The difference can be
+C seen in Lists 24 and 26 of the case-summary statistics:
+C With R : Size 21-30: 9 0 13 6 -9999 36 -9999 ...
+C Without: Size 21-30: 9 0 12 3 -9999 9 -9999 ...
+C Of course, the latter should simulate faster than the former. Resistance
+C is high enough so the solution changes little. For example, the two
+C printer plots are identical.
+ RECA ENDA 1.E+8 { Leakage resistanc ties 2 subnetworks together
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+C Begin with sources of the 1st subcase:
+14LINEA 408000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14LINEB 408000. 60. -120. { three-phase, sinusoidal excitation
+14LINEC 408000. 60. 120. { with no phasor solution.
+C Done with sources of the 1st subcase; follow by sources of 7 subcase:
+14SENDA 236000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 236000. 60. -120. { three-phase, sinusoidal excitation
+14SENDC 236000. 60. 120. { with no phasor solution.
+C --------------+------------------------------
+C From bus name | Names of all adjacent busses.
+C --------------+------------------------------
+C LINEA |ENDA *
+C ENDA |TERRA *LINEA *RECA *
+C LINEB |ENDB *
+C ENDB |TERRA *LINEB *
+C LINEC |ENDC *
+C ENDC |TERRA *LINEC *
+C SENDA |RECA *
+C RECA |ENDA *SENDA *RECB *RECC *
+C SENDB |RECB *
+C RECB |RECA *SENDB *RECC *
+C SENDC |RECC *
+C RECC |RECA *RECB *SENDC *
+C TERRA |ENDA *ENDB *ENDC *
+C --------------+------------------------------
+BLANK card terminates the last source card
+ ENDA ENDB ENDC { Arrester voltages of Y-connected 1st subcase
+C Column headings for the 9 EMTP output variables follow. These are divided among the 5 possible classes as follows ....
+C First 6 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 Step Time RECA RECB RECC ENDA ENDB ENDC RECA RECB RECC
+C RECB RECC RECA RECB RECC RECA
+C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+C 1 .5E-4 .3494E-21 -.47E-37 -.349E-21 .1226E-15 .1226E-15 .37E-15 -.334E-31 .9006E-47 .3341E-31
+C 2 .1E-3 -.349E-21 .4702E-37 .3494E-21 -.123E-15 -.123E-15 -.37E-15 .3341E-31 -.135E-46 -.334E-31
+C 22 .0011 28277.756 622.32762 -28900.08 32953.961 -15932.9 -17000.47 .2708E-5 .59598E-7 -.2768E-5
+C 23 .00115 658750.92 10864.876 -669615.8 709562.41 -421760.1 -448957.3 33.053897 .10405E-5 -50.5753
+C 400 .02 94123.287 379306.97 -473430.3 295693.53 152343.78 -601370.3 .90138E-5 .36325E-4 -.0061536
+C Variable maxima : 667041.09 609376.16 670599.37 709562.41 676288.65 651690.5 45.754917 4.3600592 52.54267
+C Times of maxima : .0152 .0036 .00985 .00115 .00455 .00985 .0152 .0036 .00985
+C Variable minima : -694529.7 -467073.7 -703549.2 -717416.9 -663770.9 -669507.6 -130.7469 -.00433 -182.8636
+C Times of minima : .00775 .0143 .0025 .0085 .01435 .00325 .00775 .0143 .0025
+BLANK card ending node voltage outputs
+ PRINTER PLOT
+ 144 2. 0.0 20. ENDA { Axis limits: (-7.174, 7.096)
+ 194 2. 0.0 20. BRANCH { Axis limits: (-1.829, 0.525)
+ RECA RECB RECB RECC RECC RECA
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 9th of 11 subcases is like the 1st except that exponential ZnO modeling
+C is used for all 3 surge arresters. The piecewise-linear 3rd arrester of
+C the 1st subcase has been replaced by a copy of the 1st arrester. Also,
+C the alternative [Z]-based Newton iteration replaces the default choice
+C of [Y]-based iteration. The request for [Z] is made by the line that
+C immediately follows the first line of the first arrester. Unlike the
+C 11th subcase of DC-37, here the request for [Z] is active. Because
+C of a lack of NO Y-BASED NEWTON declaration, that request that has
+C been added to the first arrester is active and necessary (to obtain Z).
+PRINTED NUMBER WIDTH, 11, 2, { Deliberately reduce 9 output columns by 1 digit
+ .000050 .020000
+ 1 1 1 0 1 -1
+ 5 5 20 1 30 5 50 50
+-1SENDA RECA .305515.8187.01210 200. 0 { 200-mile, constant-
+-2SENDB RECB .031991.5559.01937 200. 0 { parameter, 3-phase
+-3SENDC RECC { transmission line.
+92RECA 5555. { 1st card of 1st of 3 ZnO arrest} 3
+ [Z]-based Newton iteration { Column and case matter. Declare not use of [Y]
+C VREF VFLASH VZERO COL
+ 778000. -1.0 0.0 4.0
+C COEF EXPON VMIN
+ 625. 26. 0.5
+ 9999.
+92RECB RECA 5555. { Phase "b" ZnO is copy of "a" } 3
+92RECC RECA 5555. { Phase "c" ZnO is copy of "a" } 3
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+C SENDA 208000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14SENDA 408000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 408000. 60. -120. { three-phase, sinusoidal excitation
+14SENDC 408000. 60. 120. { with no phasor solution.
+BLANK card follows the last source card
+ SENDA SENDB SENDC
+BLANK card ending node voltage outputs
+ PRINTER PLOT
+ 194 2. 0.0 20. BRANCH
+ RECA RECB RECC
+ CALCOMP PLOT
+ 184 2. 0.0 20. BRANCH
+ RECA RECB RECC
+ 194 2. 0.0 20. BRANCH
+ RECA RECB RECC
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 10th of 11 subcases is like the 1st except that exponential ZnO modeling
+C is replaced by piecewise-linear modeling for all 3 surge arresters. Such
+C modeling became available 2 February 2007 for [Z]-based Newton iteration
+C which continues to be used in place of the default [Y]-based iteration.
+PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
+C Demonstrate that the following request for [Z]-based Newton iteration is a
+C binary toggle. Note that 3 uses has the same effect as a single use:
+NO Y-BASED NEWTON { Every subnetwork is to be solved using [Z] rather than [Y]
+NO Y-BASED NEWTON { 2nd use cancels the 1st. At this point, use [Y] not [Z]
+NO Y-BASED NEWTON { Every subnetwork is to be solved using [Z] rather than [Y]
+C ZINC OXIDE STARTUP 20 1.D-8 1.D-3 0.1 1.0 1.5
+ .000050 .020
+ 1 1 1 0 1 -1
+ 5 5 20 1 30 5 50 50
+-1SENDA RECA .305515.8187.01210 200. 0 { 200-mile, constant-
+-2SENDB RECB .031991.5559.01937 200. 0 { parameter, 3-phase
+-3SENDC RECC { transmission line.
+92RECA 4444. { 1st card of 1st of 3 ZnO arres } 1
+C VREF VFLASH VZERO
+ 0.0 -1.0 0.0
+ 0.0 0.0 { Origin. 3rd quadrant copy
+ 1.0 582400. { First point of i-v curve.
+ 2.0 590800. { Data is copied from DC-39
+ 5.0 599200. { which was used to create
+ 10. 604800. { the ZnO branch cards that
+ 20. 616000. { are used in phases "a" &
+ 50. 630000. { "b". But there is some
+ 100. 644000. { distortion due to the use
+ 200. 661920. { of linear rather than the
+ 500. 694400. { more accurate exponential
+ 1000. 721280. { modeling, of course.
+ 2000. 756000.
+ 3000. 778400. { Last point of i-v curve.
+ 9999. { Terminator for piecewise-linear characteristic
+92RECB RECA 4444. { Phase "b" ZnO is copy of "a" } 1
+92RECC RECA 4444. { Phase "c" ZnO is copy of "a" } 1
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 408000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 408000. 60. -120. { three-phase, sinusoidal excitation
+14SENDC 408000. 60. 120. { with no phasor solution.
+BLANK card follows the last source card
+ RECA RECB RECC { Names of nodes for voltage output
+C First 3 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 Step Time RECA RECB RECC RECA RECB RECC
+C TERRA TERRA TERRA
+C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
+C 1 .5E-4 .615267E-15 .615327E-15 .615387E-15 -.868E-18 -.86821E-18 -.86842E-18
+C 2 .1E-3 -.61527E-15 -.61533E-15 -.61539E-15 .867995E-18 .868206E-18 .868416E-18
+C 21 .00105 .615267E-15 .615327E-15 .615387E-15 -.868E-18 -.86821E-18 -.86842E-18
+C 22 .0011 32933.78136 -15929.2096 -17004.5718 .0565483883 -.027350978 -.02919741
+C 23 .00115 674022.6244 -433252.07 -460655.512 311.7853238 -.743908087 -.790960701
+BLANK card ending requests for node voltage output
+C 400 .02 248862.4504 197049.3444 -599097.73 .4273050316 .3383402205 -4.96347489
+C Variable maxima : 674022.6244 656282.6288 651788.8183 311.7853238 168.5414551 143.4643879
+C Times of maxima : .00115 .00455 .0098 .00115 .00455 .0098
+C Variable minima : -675779.414 -635023.773 -669662.422 -328.011827 -67.9420472 -271.51252
+C Times of minima : .00865 .01435 .00325 .00865 .01435 .00325
+ PRINTER PLOT
+ 194 3. 0.0 20. BRANCH { Axis limits: ( -3.280, 3.118 )
+ RECA RECB RECC
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+C 11th of 11 subcases is like the 1st. But the 1st was solved by [Y]-based
+C Newton iteration. Here, use [Z]-based iteration. Answer is the same.
+C Note that there is no NO Y-BASED NEWTON request because the one used
+C by the preceding subcase remains in effect. The choice was set to [Z].
+C NO Y-BASED NEWTON { If data is removed as separate subcase, activate this card
+ .000050 .020000
+ 1 1 1 0 1 -1
+ 5 5 20 1 30 5 50 50
+-1SENDA RECA .305515.8187.01210 200. 0 { 200-mile, constant-
+-2SENDB RECB .031991.5559.01937 200. 0 { parameter, 3-phase
+-3SENDC RECC { transmission line.
+92RECA 5555. { 1st card of 1st of 3 ZnO arresters
+C VREF VFLASH VZERO COL
+ 778000. -1.0 0.0 4.0
+C COEF EXPON VMIN
+ 625. 26. 0.5
+ 9999.
+92RECB RECA 5555. { Phase "b" ZnO is copy of "a"
+92RECC 4444. { Phase "c" ZnO is piecewise-linear
+C VREF VFLASH VZERO
+ 0.0 -1.0 0.0
+ 1.0 582400. { First point of i-v curve.
+ 2.0 590800. { Data is copied from DC-39
+ 5.0 599200. { which was used to create
+ 10. 604800. { the ZnO branch cards that
+ 20. 616000. { are used in phases "a" &
+ 50. 630000. { "b". But there is some
+ 100. 644000. { distortion due to the use
+ 200. 661920. { of linear rather than the
+ 500. 694400. { more accurate exponential
+ 1000. 721280. { modeling, of course.
+ 2000. 756000.
+ 3000. 778400. { Last point of i-v curve.
+ 9999. { Terminator for piecewise-linear characteristic
+BLANK card follows the last branch card
+BLANK line terminates the last (here, nonexistent) switch
+14SENDA 408000. 60. 0.0 { 1st of 3 sources. Note balanced,
+14SENDB 408000. 60. -120. { three-phase, sinusoidal excitation
+14SENDC 408000. 60. 120. { with no phasor solution.
+BLANK card follows the last source card
+ 1
+ PRINTER PLOT
+ 144 3. 0.0 20. RECA { Axis limits: (-7.174, 7.096)
+BLANK termination to plot cards
+BEGIN NEW DATA CASE
+BLANK
|