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/dc19.dat | 355 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 355 insertions(+) create mode 100644 benchmarks/dc19.dat (limited to 'benchmarks/dc19.dat') diff --git a/benchmarks/dc19.dat b/benchmarks/dc19.dat new file mode 100644 index 0000000..d803ae7 --- /dev/null +++ b/benchmarks/dc19.dat @@ -0,0 +1,355 @@ +BEGIN NEW DATA CASE +C BENCHMARK DC-19 +C Test of the control of TACS sources by node voltages and switch currents +C of the electric network. The electric network has no dynamics, with both +C the node voltage and switch current equal to one half, always. TACS +C supplemental variable GENT adds these to give unity. The output signal +C is then passed through a function block 1/(1+S) to give ELEMT, and +C through a S/(1+S) block to give GT/ELT. There also is a disconnected +C validation of the Type-17 electric network source, which provides for +C modulation of the following source function by a TACS variable. For use +C here, the TACS variable is the constant TWO = 2.0 --- easy to verify. + .05 2.0 + 1 1 1 1 1 +TACS HYBRID + 1ELEMT +GENT + 1.0 + 1.0 1.0 + 1GT/ELT +GENT + 1.0 + 1.0 1.0 +90BUS2 +91BUS1 +11TWO 2.0 { Constant for modulation of electric network COSINE +99GENT = BUS1 + BUS2 +33ELEMT GENT GT/ELTBUS1 BUS2 +BLANK card ending all TACS data + BUS3 BUS2 1.0 + BUS2 BUS1 1.0 + COSINE 1.0 { Dummy branch connects source of interest + SURGE 1.0 { Dummy branch connects Type-15 surge function + TAKUSG 1.0 { Dummy branch for Taku Noda's Type-15 surge +BLANK card ending electric network branches + BUS1 MEASURING 1 +BLANK card ending switches +11BUS3 1.0 +17TWO { The constant 2.0 of TACS multiplies following 1/2 to give unity: +14COSINE 0.5 0.5 { Cosine starts at peak 0.5, and also ends there +15SURGE 2.0 -1.0 -3.0 { Surge function is a smooth pulse > 0 +15TAKUSG 2.0 -1.0 -3.0 -6666. { Same except 2.0 = peak valu +BLANK card ending electric network source cards. +C Step Time BUS1 BUS2 BUS3 COSINE SURGE TAKUSG BUS1 TACS TACS TACS +C TERRA ELEMT GENT GT/ELT +C *** Switch "BUS1 " to " " closed before 0.00000000E+00 sec. +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 1 .05 0.0 0.5 1.0 0.0 .181042896 .470363242 0.5 .024390244 1.0 .975609756 + 1 +C 39 1.95 0.0 0.5 1.0 .987688341 .278788345 .724313367 0.5 .854137166 1.0 .145862834 +C 40 2.0 0.0 0.5 1.0 1.0 .265713062 .690342786 0.5 .861252427 1.0 .138747573 +C maxima : 0.0 0.5 1.0 1.0 .769799804 1.99999856 0.5 .861252427 1.0 .975609756 +C maxima : 0.0 .05 .05 2.0 .55 .55 .05 2.0 .05 .05 + PRINTER PLOT + 143 .4 0.0 2.0 SURGE COSINE { Axis limits: (-1.000, 1.000) + 193 .4 0.0 2.0 TACS GT/ELTTACS ELEMT { Axis limits: (0.000, 9.756) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 2nd of 5 subcases of DC-19 tests important supplemental variable functions +C Unlike 1st, there is no electric network (in theory). In practice, ATP +C will automatically supply a dummy electric network. This is since around +C 1993 when TACS STAND ALONE was internally converted to TACS HYBRID. + .100 6.0 + 1 1 0 1 0 0 0 0 +TACS STAND ALONE + X2 -X3 +TQ 6.6667 -100. 0.2 + 1X3 +X2 0.0 1.05 + 1.0 + 0.0 1. +11TQ 1.0 0.0 +98SUPVAR = 0.4 * ( ABS(X3) + 1.E-6 ) ** 3.0 * 2.0 +26T { Type-26 is given value equal to ATP variable in columns 3-8 +33X2 X3 T SUPVAR +C Step Time TACS TACS TACS TACS +C X2 X3 T SUPVAR +C 0 0.0 0.0 0.0 0.0 0.0 +C 1 0.1 0.2 .01 0.1 .80024E-6 +C 2 0.2 0.2 .03 0.2 .216022E-4 +BLANK card ending all TACS data cards +C 58 5.8 .203206E-4 .999997968 5.8 .799997523 +C 59 5.9 .677346E-5 .999999323 5.9 .800000774 +C 60 6.0 .22578E-5 .999999774 6.0 .800001858 + PRINTER PLOT + 143 1. 0.0 6.0 X2 SUPVAR { Axis limits: (0.000, 8.000) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 3rd of 5 subcases consists only of Type-15 lightning functions, including +C both old (the 2-exponential surge function) and new (Bernd Stein of FGH) +C models. Bernd's FGH improvement is documented in a letter dated 30.7.86 +C Smoothing of plot data points after 50 ups and downs is also illustrated. + 1.0E-8 3.E-6 + 1 1 1 1 1 -1 1 + 5 5 20 20 + LIGHTN 1.0 + ALGHTN 1.0 + LIGHT1 1.0 + ALGHT1 1.0 + ALGHT2 1.0 + ALGHT3 1.0 + STAN1 1.0 { Resistor for Standler surge added 29 Sept 00 + STAN2 1.0 { Resistor for 2nd Standler surge added 27 Oct 00 + TWOE1 1.0 { Resistor for TWO EXP surge added 19 Jan 01 + TWOE2 1.0 { Resistor for 2nd TWO EXP surge added 19 Jan 01 + HEID1 1.0 { Resistor for Heidler surge added 19 Jan 01 + HEID2 1.0 { Resistor for 2nd Heidler surge added 19 Jan 01 +C Orlando Hevia mentions legality and danger of negative L of series R-L-C. +C Add a branch to illustrate as mentioned in newsletter dated January, 2001 + HEVIA 10.E3 -3.0 { Series R-L with negative L } 1 +BLANK card ending electric network branches +BLANK card ending switches +15ALGHTN-1 1.0 2.E-6 50.E-6 5.0 +15LIGHTN-1 1.0 -20000. -500000. +15ALGHT2-1 1.0 2.E-6 50.E-6 10.0 +15ALGHT1-1 1.0 2.E-6 50.E-6 5.0 +15LIGHT1-1 1.0 -20000. -500000. +15ALGHT3-1 1.0 2.E-6 50.E-6 10.0 +C 345678901234567890123456789012345678901234567890 +C <--- Off-line Standler 3F10.0 -> +15STAN1 -1 0.6 20.E-6 0.04 Standler { "Standler" is same as -7777 +C C A B { Identify of 3 Standler parameters +C <--- In-line Standler 4F8.0 ---><-Request word-> < T-stop > +C T1 T2 Perc peak +15STAN2 -1 1.0E-6 20.0E-6 0.0 1.000 Standler in-line 0.0 1.0 +C +C <--- Off-line Two exp 3F10.0 -> +15TWOE1 -1 1.20 -36000.0 -5.0E+6 { it is an old two exp source +C C A B { Identify of 3 Two exp parameters +C <--- In-line Two exp 4F8.0 ---><-Request word-> < T-stop > +C T1 T2 Perc peak +15TWOE2 -1 1.0E-6 20.0E-6 0.0 1.0 Two exp in-line 0.0 1.0 +C +C <--- Off-line Heidler 3F10.0 --> +15HEID1 -1 0.6 1.0E-06 3.0E-5 5.0 { it is an old Heidler source +C C A B N { Identify of 4 Heidler parameters +C <-Heidler in line 4F8.0 ---><-key----------><>< T-stop > +C T1 T2 Perc peak N +15HEID2 -1 1.0E-6 20.0E-6 0.0 1.000 Heidler in-line 5 0.0 1.0 +C +14HEVIA 1.0 60. { Voltage source is shorted at 1 us } 1.E-6 +BLANK card ending electric network source cards. + STAN1 STAN2 ALGHT2ALGHT1LIGHT1ALGHTNLIGHTNALGHT3HEID1 HEID2 TWOE1 TWOE2 +C First 12 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 1 output variables are branch currents (flowing from the upper node to the lower node); +C Step Time STAN1 STAN2 ALGHT2 ALGHT1 LIGHT1 ALGHTN LIGHTN ALGHT3 HEID1 HEID2 +C +C +C TWOE1 TWOE2 HEVIA +C TERRA +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 +C 1 .1E-7 .442479244 .855069065 .25405E-20 .812971E-4 957.508161 .406485E-9 .004787541 .50811E-15 .930195E-8 .140651E-7 +C .058092768 .04993186 -.16949E-5 +C 2 .2E-7 .454691626 .88074615 .2601E-17 .002438392 -4.9670744 .130049E-7 .009550246 .51918E-12 .297563E-6 .449919E-6 +C .113331409 .097441938 -.51422E-5 +C 3 .3E-7 .461895211 .895961187 .14996E-15 .01470791 952.565899 .987364E-7 .01428824 .28952E-10 .225886E-5 .34153E-5 +C .165855128 .142646809 -.87063E-5 +C 4 .4E-7 .467007506 .906801471 .26623E-14 .048743046 -9.8846491 .415991E-6 .019001647 .473526E-9 .951554E-5 .143866E-4 +C .21579634 .185657436 -.12391E-4 +C 5 .5E-7 .470959015 .91521116 .2479E-13 .121908619 947.672889 .126925E-5 .023690588 .395202E-8 .290285E-4 .43887E-4 +C .263281003 .226579442 -.16201E-4 +BLANK card ending output requests +C 260 .26E-5 .48556872 .972472279 .991750958 -15.735499 -362.71992 .993881314 .676797074 -18.724315 .572623725 .949657569 +C 1.09277363 .94980214 -.55047746 +C 280 .28E-5 .482164403 .967184298 .987942471 -17.163302 -375.61152 .990578487 .698942172 -19.269102 .568828878 .942698718 +C 1.08493561 .942824605 -1.0722499 +C 300 .3E-5 .478686003 .961711993 .984066679 -17.995991 -387.2694 .987055255 .718634373 -19.452928 .565055727 .935784527 +C 1.07715275 .935897248 -2.0885866 +C Variable maxima : .506829483 1.00000002 1.0 1757.33053 957.508161 1.0 .718634373 2231.07839 0.6 1.0 +C 1.14948506 .999999945 0.0 +C Times of maxima : .8E-6 .1E-5 .2E-5 .7E-6 .1E-7 .2E-5 .3E-5 .113E-5 .1E-5 .1E-5 +C .99E-6 .1E-5 0.0 +C Variable minima : 0.0 0.0 0.0 -17.995991 -387.2694 0.0 0.0 -19.452928 0.0 0.0 +C 0.0 0.0 -2.0885866 +C Times of minima : 0.0 0.0 0.0 .3E-5 .3E-5 0.0 0.0 .3E-5 0.0 0.0 +C 0.0 0.0 .3E-5 + PRINTER PLOT + 145 .4 0.0 3.0 LIGHTNALGHTNALGHT2 { Plot limits: (0.000, 1.000) + 145 .4 0.0 3.0 LIGHT1ALGHT1ALGHT3 { Plot limits: (-0.110, 2.231) +C 145 .4 0.0 3.0 STAN1 STAN2 +C 145 .4 0.0 3.0 HEID1 HEID2 +C 145 .4 0.0 3.0 TWOE1 TWOE2 +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 4th of 5 subcases was added 17 September 1999 to illustrate delayed +C opening of a Type-13 TACS-controlled switch. There are two parallel, +C identical circuits on the electrical side. A sinusoidal source drives +C current through a series R-L branch. The old logic is illustrated by +C the QUICK alternative whereas the new is the DELAY alternative. +C The new optional DO NOT OPEN UNTIL CURRENT 0 is illustrated by the +C latter of these two: opening will be delayed until current passes +C through zero. Both circuits oscillate following opening because both +C have a hanging inductor. But QUICK opens first, and it oscillates +C more. The graph is distinctive. +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) + .0004 .010 + 1 1 +TACS HYBRID +14CLAMP 100. 50. 30. 0. -1. +33CLAMP +BLANK card ending TACS + GEN QUICK 10. 4. { Series R-L will oscillate when opened + GEN DELAY 10. 4. +BLANK card ending branches +13DELAY Do not open until current 0 CLOSED CLAMP 11 +$DISABLE +C The preceding involves a high-level, English-language declaration. It +C is equivalent to the following low-level, numeric declation. +13RESIS -77333. CLOSED CLAMP 11 +C I.e., the data field of otherwise-unused columns 25-34 can be used to +C hold the special numeric value. About the English, this must begin +C to the right of column 24. If 55-60 are used as here for CLOSED, the +C beginning of "do not open until current 0" can be delayed at most 3 +C columns (note 3 blanks between ending 0 and start of CLOSED). Case is +C arbitrary (note 1st letter has been made upper case on switch card). +$ENABLE +13QUICK CLOSED CLAMP 11 +BLANK card ending switches +14GEN 100. 50. 0.0 0. -1. +BLANK card ending sources + GEN QUICK DELAY +BLANK card ending voltage printout + CALCOMP PLOT { Needed to cancel PRINTER PLOT of preceding subcase + SUPERIMPOSE 2 { Superimpose plots of the following 2 cards + 144 1. 0.0 10. GEN DELAY QUICK { First, three node voltages + 194 1. 0.0 10.-200.200.TACS CLAMP { 2nd of 2 in overlay is I-branch +BLANK card ending plot +BEGIN NEW DATA CASE +C 5th of 5 subcases was added 28 July 2002 to illustrate 3 new source +C types: CIGRE, CESI, and USRFUN. These are all special sub-types +C within the Type-15 source. The supporting code comes from Orlando P. +C Hevia of UTN in Santa Fe, Argentina. For an alternate and different +C illustration of the USRFUN alternative, see DC-5. For comments +C about required precision, see note following this subcase. + .05E-6 10.E-6 { Hevia used dT = 1.0E-9, but such detail is not needed + 1 1 1 1 -1 + 5 5 20 20 + CESI1 1.0 + CESI2 1.0 + CIGRE1 1.0 + CIGRE4 1.0 + USER2 1.0 + USER8 1.0 { Add another USRFUN test on 10 December 2002 +BLANK card ending branch cards +BLANK card ending switch cards (none for this data) +C CORRENTE DI CRESTA -0.001 [kA] +C TEMPO DI CRESTA 1 [us] +C TEMPO EMIVALORE 10 [us] +C GENERATORI EQUIVALENTI FRONTE +15CESI1 -1 1.5 1.0E-6 10.0-6 CESI +15CESI2 -1 -1.2 1.0E-6 10.0-6 CESI +15CIGRE1 1.0 1.0E-6 10.0E-6 CIGRE 20.0E+5 0.5E-6 9.5E-6 +15CIGRE4 1.0 1.0E-6 10.0E-6 CIGRE 50.0E+5 0.5E-6 9.5E-6 +15USER2 2.0 USRFUN +15USER8 8.0 USRFUN +BLANK card ending electric network sources + 1 { Request for all possible node voltage outputs (here, just 4) + CALCOMP PLOT + 145 .5 0.0 5.0 CESI1 CESI2 USER2 CIGRE1 +BLANK card ending plot cards +BEGIN NEW DATA CASE +BLANK +EOF + +Warning about the preceding. The CIGRE1 output depends greatly on arithmetic +precision. Even 64 bits seems inadequate since turning the Salford debugger +on and off (minor update using SUBR5) results in the following differences. +Note that the 3rd row is what will be published in a future newsletter. The +CIGRE1 value on step 40 rises from .661 to .668 when the debugger is +turned on. + +FC version 2.0 - Copyright (c) 1990 Mike Albert Wed Aug 21 14:08:13 2002 +Options are: -c1 -ds5 -t8 + +Compare files: + dc19.lis 71885 8-21-102 2:07p + dc19.sal 71885 8-21-102 1:48p + +Changed lines 859-868 + > 15 .75E-6 1.000594 -.800475 .0648095 .0662431 .0783034 + > 20 .1E-5 1.477403 -1.18192 .1299899 .1324862 .1044063 + > 40 .2E-5 1.338601 -1.07088 .668485 .4280678 .2093254 + > 60 .3E-5 1.23813 -.990504 .9529333 .9505637 .3267037 + > 80 .4E-5 1.167656 -.934125 .8811803 .8800206 .5541328 + > 100 .5E-5 1.104946 -.883957 .8148301 .8147126 .9983771 + > 120 .6E-5 1.044245 -.835396 .7534758 .7542513 .9931261 + > 140 .7E-5 .9840634 -.787251 .6967413 .6982769 .9835613 + > 160 .8E-5 .9240164 -.739213 .6442788 .6464564 .9740871 + > 180 .9E-5 .8640042 -.691203 .5957665 .5984817 .9647042 +To + > 15 .75E-6 1.000594 -.800475 .0644145 .0662431 .0783034 + > 20 .1E-5 1.477403 -1.18192 .1291901 .1324862 .1044063 + > 40 .2E-5 1.338601 -1.07088 .661331 .4280678 .2093254 + > 60 .3E-5 1.23813 -.990504 .9533223 .9505637 .3267037 + > 80 .4E-5 1.167656 -.934125 .88154 .8800206 .5541328 + > 100 .5E-5 1.104946 -.883957 .8151627 .8147126 .9983771 + > 120 .6E-5 1.044245 -.835396 .7537834 .7542513 .9931261 + > 140 .7E-5 .9840634 -.787251 .6970257 .6982769 .9835613 + > 160 .8E-5 .9240164 -.739213 .6445418 .6464564 .9740871 + > 180 .9E-5 .8640042 -.691203 .5960097 .5984817 .9647042 + +Changed line 872 + >Variable maxima : 1.497392 0.0 .999673 .9989941 .9999945 +To + >Variable maxima : 1.497392 0.0 .9997968 .9989941 .9999945 + +Comparison complete + +21 August 2002, remove from before the 5th, which was not being executed: +BLANK ----------------- Protect following 4th data subcase from execution +C Note: Following is copy of 3rd subcase between October, 1990, and +C July, 1991. It was wrong during 9 months. The solution +C prior to October, 1990, was kept below here, and discovered +C to be correct once again on July 22, 1991. So, it would +C seem that Guido's "correction" of Oct, 1990, was wrong, and +C during the spring or summer of 1991 he merely corrected the +C correction. Well, as history, we keep the wrong data: +C 3rd of 3 subcases consists only of Type-15 lightning functions, including +C both old (the 2-exponential surge function) and new (Bernd Stein of FGH) +C models. Bernd's FGH improvement is documented in a letter dated 30.7.86 +C Smoothing of plot data points after 50 ups and downs is also illustrated. +C Comment cards altered 23 Oct 1990. Original contents preserved below +C as 4th data subcase until Bernd Stein approves of changes. + 1.E-8 300.E-8 + 1 1 1 1 1 -1 + 5 5 20 20 + LIGHTN 1.0 + ALGHTN 1.0 + LIGHT1 1.0 + ALGHT1 1.0 + ALGHT2 1.0 + ALGHT3 1.0 +BLANK card ending electric network branches +BLANK card ending switches +15ALGHTN-1 1.0 2.E-6 50.E-6 5.0 +15LIGHTN-1 1.0 -20000. -500000. +15ALGHT2-1 1.0 2.E-6 50.E-6 10.0 +15ALGHT1-1 1.0 2.E-6 50.E-6 5.0 +15LIGHT1-1 1.0 -20000. -500000. +15ALGHT3-1 1.0 2.E-6 50.E-6 10.0 +BLANK card ending electric network source cards. +C Step Time ALGHT2 ALGHT1 LIGHT1 ALGHTN LIGHTN ALGHT3 +C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 +C 1 .1E-7 .26503E-20 .853443E-4 957.508161 .426722E-9 .004787541 .53005E-15 +C 2 .2E-7 .27133E-17 .002559784 -4.9670744 .136524E-7 .009550246 .54161E-12 +C 3 .3E-7 .15643E-15 .015440124 952.565899 .103652E-6 .01428824 .30202E-10 +C 4 .4E-7 .27773E-14 .051169655 -9.8846491 .436701E-6 .019001647 .493979E-9 + 1 +C 300 .3E-5 1.02657121 -18.891898 -387.2694 1.03619451 .718634373 -20.293153 +C maxima : 1.04319273 1844.81694 957.508161 1.0497837 .718634373 2327.44477 +C Times of max : .2E-5 .7E-6 .1E-7 .2E-5 .3E-5 .113E-5 +C minima : 0.0 -18.891898 -387.2694 0.0 0.0 -20.293153 +C Times of min : 0.0 .3E-5 .3E-5 0.0 0.0 .3E-5 + PRINTER PLOT + 145 .4 0.0 3.0 LIGHTNALGHTNALGHT2 { Plot limits: (0.000, 1.050) + 145 .4 0.0 3.0 LIGHT1ALGHT1ALGHT3 { Plot limits: (-0.110, 2.327) +BLANK card ending plot cards +BEGIN NEW DATA CASE + -- cgit v1.2.3