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/dc10.dat | 378 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 378 insertions(+) create mode 100644 benchmarks/dc10.dat (limited to 'benchmarks/dc10.dat') diff --git a/benchmarks/dc10.dat b/benchmarks/dc10.dat new file mode 100644 index 0000000..7e37b09 --- /dev/null +++ b/benchmarks/dc10.dat @@ -0,0 +1,378 @@ +BEGIN NEW DATA CASE +C BENCHMARK DC-10 +C Companion to "CASCADED PI" test of preceding data case. Same answer +C Use batch-mode solution only, since SPY sets TMAX to nearly infinity, +C and then the 0 DELTAT will result in an error stop. 6 subcases total. +C 20 June 2007, add the following comment merely to satisfy a new 6th +C subcase which will illustrate CONCATENATE MONTE CARLO PL4 ( CMCP ) . +C Since subcases 1 through 5 do not use .PL4 files significantly, the +C switch from C-like to widexx format should go unnoticed : +C FMTPL4 L4BYTE NEWPL4 Next card FORMAT ( 16X, A6, 2X, 2I8 ) +C CONCATENATE. wide08 0 0 key text anywhere +C The preceding overrides .PL4 choices of STARTUP. Service begins 17 June 2007 +C A second illustration of just this declaration can be found in DCNEW-22. + 0.0 0.0 60. + 1 1 + USE OLD { Extraneous request to undo preceding USE RL or AB will be rejected +C Note about preceding, which was added 5 November 2001. This first +C occurence will illustrate the warning message. For a legal use of +C USE OLD, see the 2nd subcase. + 1RA1 1A1 .877 8.40 .1628 + 2RB1 1B1 .747 4.14-.0252 .852 8.43 .1559 + 3RC1 1C1 .735 3.47-.0067 .723 4.17-.0277 .829 8.46 .1571 + 11A1 2A1 RA1 1A1 + 21B1 2B1 + 31C1 2C1 + 02B1 3B1 13.1449.071 3 + 12C1 4C1 RA1 1A1 + 23A1 4A1 + 33B1 4B1 + 04C1 5C1 13.14 1 + 04A1 5A1 13.1449.071 2 + 04B1 5B1 4C1 5C1 + 05C1 5G1 2B1 3B1 + 05A1 5G1 2B1 3B1 + 05B1 5G1 2B1 3B1 + 5G1 5.0 13.14 + 15B1 6B1 RA1 1A1 + 25C1 6C1 + 35A1 6A1 + 6A1 7A1 4B1 5B1 + 6B1 7B1 4C1 5C1 + 6C1 7C1 4A1 5A1 + 7A1 7G1 2B1 3B1 + 7B1 7G1 2B1 3B1 + 7C1 7G1 2B1 3B1 + 7G1 5G1 + 17A1 GA1 1.658 16.92 .3142 + 27B1 GB1 1.446 8.34-.0554 1.704 16.86 .3118 + 37C1 GC1 1.470 6.94-.0134 1.494 8.28-.0504 1.754 16.80 .3256 +BLANK card ending branch cards +BLANK card ending switch cards (none, for this problem) +14GA1 424.35 60. 0.0 -.1 +14RA1 424.35 60. 10.0 -.1 +14GB1 424.35 60. -120.0 -.1 +14RB1 424.35 60. -110.0 -.1 +14GC1 424.35 60. 120.0 -.1 +14RC1 424.35 60. 130.0 -.1 +C --------------+------------------------------ +C From bus name | Names of all adjacent busses. +C --------------+------------------------------ +C RA1 |1A1 * +C 1A1 |RA1 *2A1 * +C RB1 |1B1 * +C 1B1 |RB1 *2B1 * +C RC1 |1C1 * +C 1C1 |RC1 *2C1 * +C 2A1 |1A1 * +C 2B1 |1B1 *3B1 * +C 2C1 |1C1 *4C1 * +C 3B1 |2B1 *4B1 * +C 4C1 |2C1 *5C1 * +C 3A1 |4A1 * +C Etc. (many more connections!) ..... +BLANK card ending all source cards +C Total network loss P-loss by summing injections = 9.311041032866E+03 +C End injection: -12.95674346101 44.419110587004 -6432.468410934 9424.6247887975 +C End injection: -42.48741206788 -106.9593405 -6888.171205186 -0.6825172 +C Selective branch outputs follow (for column-80 keyed branches only). Any req +C augmented to include branch voltage. But the converse is not true (a request +C From To (======== Branch voltage Vkm = Vk - Vm =========) (== +C bus K bus M Magnitude Degrees Real part Imag part Magni +C 2B1 3B1 1.6180264E+03 -139.359766 -1.2277813E+03 -1.0538323E+03 3.9545 +C 4C1 5C1 3.9925372E+02 -138.975087 -3.0120669E+02 -2.6206500E+02 3.0384 +C 4A1 5A1 2.9415954E+00 92.386571 -1.2249239E-01 2.9390439E+00 +BLANK card ending output variable requests + PRINTER PLOT +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 2 of 6 subcases is unrelated to the first. It shows "BCTRAN" usage. +C Illustration of 9x9 matrix representation of 3-phase, 3-winding +C transformer as punched by the 2nd subcase of DCNEW-8. Note the use +C of "USE AR" for better representation via [A] and [R] matrices. +C To validate, compare with the same solution in 3rd subcase which has +C no "USE AR" (it uses the more common "USE RL" for [R], [X]). +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .000200 .010 + 1 1 1 1 1 -1 + 5 5 + USE AR +$VINTAGE, 1, + 1H-1 26.512692374898 .2054666 + 2L-1 -59.57848438329 0.0 + 180.85474339841 .0742333 + 3T-1 T-2 5.1245421614045 0.0 + -71.06950226953 0.0 + 76.560711306734 .0822 + 4H-2 1.317410104146 0.0 + -1.044760156768 0.0 + -2.174181664407 0.0 + 26.512692374898 .2054666 + 5L-2 -1.044760156768 0.0 + .10024670965949 0.0 + 2.6475868140585 0.0 + -59.57848438329 0.0 + 180.85474339841 .0742333 + 6T-2 -2.174181664407 0.0 + 2.6475868140585 0.0 + 2.4174362478358 0.0 + 5.1245421614045 0.0 + -71.06950226953 0.0 + 76.560711306734 .0822 + 7H-3 1.317410104146 0.0 + -1.044760156768 0.0 + -2.174181664407 0.0 + 1.317410104146 0.0 + -1.044760156768 0.0 + -2.174181664407 0.0 + 26.512692374898 .2054666 + 8L-3 -1.044760156768 0.0 + .10024670965949 0.0 + 2.6475868140585 0.0 + -1.044760156768 0.0 + .10024670965949 0.0 + 2.6475868140585 0.0 + -59.57848438329 0.0 + 180.85474339841 .0742333 + 9 T-1 -2.174181664407 0.0 + 2.6475868140585 0.0 + 2.4174362478358 0.0 + -2.174181664407 0.0 + 2.6475868140585 0.0 + 2.4174362478358 0.0 + 5.1245421614045 0.0 + -71.06950226953 0.0 + 76.560711306734 .0822 +C USE RL This old, absolute form is replaced by the following universal form: + USE OLD { Beginning 5 November 2001, this will undo preceding USE AB or RL + GENA H-1 5.0 + GENB H-2 5.0 + GENC H-3 5.0 + L-1 1.0 +BLANK card ending all branch cards +BLANK card terminating the last switch card +C Balanced 3-phase source of 400 Hz is connected to transformer delta side: +14GENA 100. 60. -30. -1. +14GENB 100. 60. -150. -1. +14GENC 100. 60. 90. -1. +BLANK card ending all source cards +C Total network loss P-loss by summing injections = 1.731473944895E+02 +C Last gen: GENC .6123031769E-14 100. .11363284501067 .62844772817969 +C Last gen: 100. 90.0000000 .61808909033326 79.5827588 + H-1 L-1 H-2 L-2 H-3 L-3 +C Step Time H-1 L-1 H-2 L-2 H-3 +C 0 0.0 83.84032026 .9255817033 -87.1526014 -35.2685815 -.568164225 +C 1 .2E-3 85.28151988 1.894879197 -82.9051812 -33.2151021 -7.86654662 +C 2 .4E-3 86.23807533 2.853434433 -78.1866756 -30.9728926 -15.1202281 +C 3 .6E-3 86.7045542 3.795799431 -73.0238961 -28.5546937 -22.2879919 +BLANK card ending program output-variable requests. +C 50 .01 -80.9365963 -8.33136565 39.29673143 13.29079634 57.4223614 +C Variable max : 86.7045542 12.93833086 102.0061131 43.76610174 57.4223614 +C Times of max : .6E-3 .004 .0068 .0066 .01 +C Variable min : -86.757459 -8.33136565 -87.1526014 -35.2685815 -96.8960673 +C Times of min : .009 .01 0.0 0.0 .0042 + PRINTER PLOT + 144 2. 0.0 10. L-1 L-2 L-3 { Axis limits: (-4.603, 4.377) +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C 3rd of 6 subcases has the same solution as the 2nd, but it is to be +C obtained using [R], [L] branch cards (punched by the 1st subcase +C of DCNEW-8) rather than [A], [R] branch cards of preceding case. +C Note XOPT = 60, so we use not L but rather wL in ohms at 60 Hz. +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .000200 .010 60. { Note XOPT = 60 to agree with DCNEW-8 punched cards + 1 1 1 1 1 -1 + 5 5 + USE RL { Should not be needed, since RL units are the default starting choice +$VINTAGE, 1, + 1H-1 .2054666 41432.097487177 + 2L-1 0.0 19771.027634356 + .0742333 9437.8794714269 + 3T-1 T-2 0.0 15579.567888173 + 0.0 7437.5502771391 + .0822 5866.2181589463 + 4H-2 0.0 -.0533106335013 + 0.0 .95599943703728 + 0.0 1.7517639180077 + .2054666 41432.097487177 + 5L-2 0.0 .95599943703728 + 0.0 .73764739804435 + 0.0 .83627686643027 + 0.0 19771.027634356 + .0742333 9437.8794714269 + 6T-2 0.0 1.7517639180077 + 0.0 .83627686643027 + 0.0 .65959655490825 + 0.0 15579.567888173 + 0.0 7437.5502771391 + .0822 5866.2181589463 + 7H-3 0.0 -.0533106335013 + 0.0 .95599943703728 + 0.0 1.7517639180077 + 0.0 -.0533106335013 + 0.0 .95599943703728 + 0.0 1.7517639180077 + .2054666 41432.097487177 + 8L-3 0.0 .95599943703728 + 0.0 .73764739804435 + 0.0 .83627686643027 + 0.0 .95599943703728 + 0.0 .73764739804435 + 0.0 .83627686643027 + 0.0 19771.027634356 + .0742333 9437.8794714269 + 9 T-1 0.0 1.7517639180077 + 0.0 .83627686643027 + 0.0 .65959655490825 + 0.0 1.7517639180077 + 0.0 .83627686643027 + 0.0 .65959655490825 + 0.0 15579.567888173 + 0.0 7437.5502771391 + .0822 5866.2181589463 + GENA H-1 5.0 + GENB H-2 5.0 + GENC H-3 5.0 + L-1 1.0 +BLANK card ending all branch cards +BLANK card terminating the last switch card +C Balanced 3-phase source of 400 Hz is connected to transformer delta side: +14GENA 100. 60. -30. -1. +14GENB 100. 60. -150. -1. +14GENC 100. 60. 90. -1. +BLANK card ending all source cards +C Total network loss P-loss by summing injections = 1.731473944910E+02C 100. .11363284500541 .6284477281075 30.904454513041 31.422386405375 +C 90.0000000 .61808909026083 79.5827588 5.6816422502705 0.9835171 +C Step Time H-1 L-1 H-2 L-2 H-3 +C 0 0.0 83.84032026 .9255817034 -87.1526014 -35.2685815 -.568164225 +C 1 .2E-3 85.28151988 1.894879197 -82.9051812 -33.2151021 -7.86654662 +C 2 .4E-3 86.23807533 2.853434433 -78.1866756 -30.9728926 -15.1202281 +C 3 .6E-3 86.7045542 3.795799431 -73.0238961 -28.5546937 -22.2879919 + H-1 L-1 H-2 L-2 H-3 L-3 +BLANK card ending program output-variable requests. +C 50 .01 -80.9365963 -8.33136565 39.29673143 13.29079634 57.4223614 +C Variable max : 86.7045542 12.93833086 102.0061131 43.76610174 57.4223614 +C Times of max : .6E-3 .004 .0068 .0066 .01 +C Variable min : -86.757459 -8.33136565 -87.1526014 -35.2685815 -96.8960673 +C Times of min : .009 .01 0.0 0.0 .0042 + PRINTER PLOT + 144 2. 0.0 10. L-1 L-2 L-3 { Axis limits: (-4.603, 4.377) +$WIDTH, 79, +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C 4th of 6 subcases is unrelated to the preceding. It shows "SSOMIT" +C usage as introduced 7 November 1998. More precisely, this 4th case +C does _not_ involve SSOMIT to omit a linear branch during the +C phasor solution. The switch (LEFT, RIGHT) is open, isolating 50 Hz +C on the left from 60 Hz on the right. But, because this switch is +C paralleled by a leakage branch (snubber circuits in case of hvdc), +C phasor solution logic will reject the mixing of frequencies. Use of +C an all-resistive network allows a large time step for illustration. +$WIDTH, 132, { Back to full width to make new 4th and 5th subcases look better +PRINTED NUMBER WIDTH, 9, 2, + .001 .020 + 1 1 1 1 1 + GEN1 LEFT 0.5 1 + RIGHT GEN2 0.5 1 + LEFT RIGHT 1.E4 +BLANK card ending branch cards + LEFT RIGHT .010 1.0 1 +BLANK card ending switch cards +14GEN1 1.0 50. 0. -1. +14GEN2 1.0 60. 0. -1. +BLANK card ending source cards + 1 +C PRINTER PLOT + 194 2. 0.0 20. LEFT RIGHT +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 5th of 6 subcases is the same as the preceding except for the name +C of the snubber branch that parallels the switch. It shows "SSOMIT" +C has the effect of omitting this branch during the phasor solution, +C thereby avoiding the illegal mixing of 50- and 60-Hz excitation (a +C problem of the preceding data subset). See January, 1999, newsletter. +PRINTED NUMBER WIDTH, 9, 2, + .001 .020 + 1 1 1 1 1 + GEN1 LEFT 0.5 1 + RIGHT GEN2 0.5 1 + LEFT RIGHT NAME SSOMIT 1.E4 { Snubber to be ignored during phasor solution +BLANK card ending branch cards + LEFT RIGHT .010 1.0 { Switch is open during phasor solution } 1 +BLANK card ending switch cards +14GEN1 1.0 50. 0. -1. +14GEN2 1.0 60. 0. -1. +BLANK card ending source cards +C First 4 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 RIGHT LEFT GEN1 GEN2 LEFT GEN1 RIGHT +C RIGHT LEFT GEN2 +C 0 0.0 1.0 1.0 1.0 1.0 0.0 0.0 0.0 +C 1 .1E-2 .929778 .951055 .951057 .929776 0.0 .213E-5 .213E-5 +C 2 .002 .728973 .809013 .809017 .728969 0.0 .8E-5 .8E-5 +C 3 .003 .425787 .587777 .587785 .425779 0.0 .162E-4 .162E-4 + 1 { Output all node voltages +C 19 .019 .79424 .79424 .951057 .637424 .313633 .313633 .313633 +C 20 .02 .654508 .654508 1.0 .309017 .690983 .690983 .690983 +C Variable max : 1.0 1.0 1.0 1.0 .690983 .690983 .690983 +C Times of max : 0.0 0.0 0.0 0.0 .02 .02 .02 +C Variable min : -.99211 -.99999 -1. -.99211 -.84484 -.84484 -.84484 +C Times of min : .008 .01 .01 .008 .014 .014 .014 + PRINTER PLOT + 194 3. 0.0 20. LEFT RIGHT { Axis limits: (-8.448, 6.910) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 6th of 6 subcases is unrelated to preceding. It illustrates use of +C CONCATENATE MONTE CARLO PL4 ( CMCP ) which becme available with the +C UTPF update of 20 June 2007. Data is copied from DC-24, so the +C solution is well known. What is new is the .PL4 file, which will +C contain all 3 signal sets, for the 3 energizations, in 1 disk file. +C For illustration, use formatted .PL4 as per the +C declaration of the 1st subcase (mandatory location if present). +CONCATENATE MONTE CARLO PL4 { New form of IPLOT= 54321 frees IPLOT for frequency +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) +REPEATABLE RANDOM NUMBERS { Same as NSEED=1 (this overrides zero value below) +OMIT BASE CASE { Comment out this card, if a base case solution is to be added + 100.E-6 20.E-3 60. + 1 1 1 2 1 -1 1 1 3 +C ISW ITEST IDIST IMAX IDICE KSTOUT NSEED + 1 1 0 0 1 { KSTOUT is blank, not 0! } 0 + 2 2 10 10 { Printout frequency change only if base case + 0GENA A1 7. 2 + 0GENB B1 7. 1 + 0GENC C1 7. 2 + 0ENDA A10 7. 3 + 0ENDB B10 7. + 0ENDC C10 7. +-1ASW1 A5 .3 2.1146 0.645 50. 0 +-2BSW1 B5 .0268 .5397 0.021 50. 0 +-3CSW1 C5 + 0A5 A5F 1. + 0B5 B5F 1. + 0C5 C5F 1. +-1A5F ASW10 ASW1 A5 +-2B5F BSW10 +-3C5F CSW10 +BLANK card ending branch cards +76A1 ASW1 2.E-3 .1E-3 STATISTICS + B1 BSW1 4.E-3 0.5E-3 STATISTICS 4 + C1 CSW1 6.E-3 1.E-3 STATISTICS + A10 ASW10 7.95E-3 1.0 + B10 BSW10 9.95E-3 1.0 + C10 CSW10 11.95E-3 1.0 +BLANK card ending switch cards +14GENA 303. 60. 0.0 -1. +14GENB 303. 60. -120.0 -1. +14GENC 303. 60. 120.0 -1. +14ENDA 303. 60. - 10.0 -1. +14ENDB 303. 60. -130.0 -1. +14ENDC 303. 60. 110.0 -1. +BLANK card ending source cards +$OPEN, UNIT=LUNIT9 FILE=dc24at40.ext STATUS=UNKNOWN FORM=UNFORMATTED RECL=1000 ! + ASW10 BSW10 CSW10 { Request for these node voltage outputs +BLANK card ending the specification of program outputs (node voltages, here) +$OPEN, UNIT=LUNIT2 STATUS=SCRATCH FORM=UNFORMATTED { Disconnect *.BIN on LUNIT2 +-1 100. GENA A1 +BLANK card ending statistical tabulation requests +BEGIN NEW DATA CASE +BLANK -- cgit v1.2.3