diff options
author | Angelo Rossi <angelo.rossi.homelab@gmail.com> | 2023-06-21 12:04:16 +0000 |
---|---|---|
committer | Angelo Rossi <angelo.rossi.homelab@gmail.com> | 2023-06-21 12:04:16 +0000 |
commit | b18347ffc9db9641e215995edea1c04c363b2bdf (patch) | |
tree | f3908dc911399f1a21e17d950355ee56dc0919ee /benchmarks/dc42.dat |
Initial commit.
Diffstat (limited to 'benchmarks/dc42.dat')
-rw-r--r-- | benchmarks/dc42.dat | 240 |
1 files changed, 240 insertions, 0 deletions
diff --git a/benchmarks/dc42.dat b/benchmarks/dc42.dat new file mode 100644 index 0000000..191a6e0 --- /dev/null +++ b/benchmarks/dc42.dat @@ -0,0 +1,240 @@ +BEGIN NEW DATA CASE
+C BENCHMARK DC-42
+C Test of batch-mode Fourier series and "X-Y PLOT" capability. Also, a
+C Type-97 staircase time-varying resistance element, and Type 1-10 source
+C is present. A Type-91 compensation-based nonlinear resistor was added.
+C For answer to Fourier series, see EMTP Rule Book, page 43-m, HWD example
+C 1st of 3 data subcases, of which only two are real (3rd is a near-dummy)
+CHANGE PRINTOUT FREQUENCY
+ 10 5
+ 1.0 60.
+ 1 1
+ VALUE 1.0
+ COS 1.0 1
+ SIN 1.0 1
+ GEN LOAD 1.0
+97LOAD -1. 1
+ 0.0 1.0
+ 2.01 0.5
+ 6.01 2.0
+ 7.01 1.E35
+ 9999
+ A1 A2 5.E+3 1
+ A2 5.E+3
+91A2 3333. 1
+ 300000.
+ 300. 0.0
+ 200. 3.0
+ 150. 6.0
+ 150. 1000.
+ 9999.
+BLANK card ends branch cards
+BLANK card ends switch cards (none, for this case)
+ 1VALUE
+14COS 4.0 .02 0.0
+14SIN 5.0 .02 -90.
+14GEN 100. 60. -1.
+14A1 620000. .06 -50. -1.
+BLANK card ends source cards
+C Total network loss P-loss by summing injections = 1.922250000000E+07
+C Step Time LOAD VALUE COS SIN GEN A2
+C 0 0.0 50. 0.0 0.0 0.0 100. 199264.159
+C Comment card. NUMDCD = 5527. |C Next come the data ca
+C Another input card for Type 1-10 sources. | 3.4
+C 1 1.0 50. 3.4 3.96845881 .626666168 100. 272691.058
+C Another input card for Type 1-10 sources. |.8485282
+C 2 2.0 50. .8485282 3.87433264 1.24344944 100. 307819.308
+C Another input card for Type 1-10 sources. | 0.0
+C 3 3.0 50. 0.0 3.71910594 1.84062276 100. 32112.3483
+ LOAD VALUE COS SIN GEN A2
+BLANK card ends requests for program outputs (here, just node voltages)
+C Next come the data cards for user-defined Type-1 source on node "VALUE":
+ 3.4
+.8485282
+ 0.0
+-.434314
+ -.2
+-.848528
+ -.4
+-1.56569
+ 3.41
+.8500000
+ 0.0
+-.430000
+ -.2
+ 9999
+C Last step: 60 60. 100. 0.0 1.23606798 4.75528258 100.
+C Last step continued ..... : -17025.944 1.E-33 -113.50629 -116.91148
+ PRINTER PLOT
+ FOURIER ON 10
+C Note HPI=7.6 is a special request for reduction of end-time HMAX until
+C a match with the HMIN point (for periodicity) is found.
+ 1437.6 1.5 12. VALUE
+C --- Back up the ending time from point number 11 to point number 8 in
+C 8.4852820634841920E-1 0.0E+00 -1.5656900405883790E0 3.4100000858306880E+00
+C Harmonic Cosine Sine Complex Fraction of
+C number coefficient coefficient amplitude fundamental
+C 0 1.01249525324E-01 0.00000000000E+00 1.01249525324E-01 0.086661418
+C 1 1.13313862336E+00 -2.84609199533E-01 1.16833451383E+00 1.000000000
+C 2 5.00001057982E-01 -9.02500022203E-01 1.03174965377E+00 0.883094389
+C 3 -2.84610506422E-01 -8.50297213563E-01 8.96665206061E-01 0.767473010
+C 4 -6.01250493899E-01 0.00000000000E+00 6.01250493899E-01 0.514621871
+ 143 5. 0.5 8.0 VALUE
+C Begin the Fourier series computation using 8 equidistant points. The first
+C 3.40000009536743 8.48528206348419E-1 -4.00000005960464E-1 -1.565690040588379
+C Harmonic Cosine Sine Complex Fraction of
+C number coefficient coefficient amplitude fundamental
+C 0 9.99995265156E-02 0.00000000000E+00 9.99995265156E-02 0.085748909
+C 1 9.99999101967E-01 6.00000909627E-01 1.16619007691E+00 1.000000000
+C 2 9.00000024587E-01 5.00001057982E-01 1.02956354939E+00 0.882843689
+C 3 8.00000947207E-01 4.00000906647E-01 8.94428443672E-01 0.766966262
+C 4 6.00000495091E-01 0.00000000000E+00 6.00000495091E-01 0.514496313
+ FOURIER OFF
+ 143 8. 0.0 50. COS { Axis limits: (-4.000, 4.000)
+ X-Y PLOT
+ 10. -5. 5.0
+ 13. -6.5 6.5
+ 2Arbitrary 78-character case title text of which this is an example, I hope.
+ First and only line of multi-line, 78-byte, graph subheading text.
+ 193 5. 0.0 55. COS SIN { Axis limits: (-6.500, 6.500)
+ CALCOMP PLOT { Preceding was a line printer plot; switch to vector graphics
+ X-Y PLOT Horizontal Axis label123
+ 10. -5. 5.0
+ 8. -6.0 6.0
+ 2Arbitrary 78-character case title text of which this is an example, I hope.
+ 143 5. 0.0 55. COS SIN 16-byte Heading Y-axis label
+BLANK card ending plot cards
+BEGIN NEW DATA CASE
+C 2nd of 4 subcases modifies preceding data for nice bar charts. As built
+C in this standard test case, both 5 (the original) and 13 vertical bars
+C result. Screen plotting of Salford is slowed slightly by HP-GL output
+C due to shading of hollow rectangular shells of bars by interior vertical
+C lines. This progress is made 25 December 1997 as this case is added,
+C forcing later cases down. See mention in January, 1998, newsletter.
+CHANGE PRINTOUT FREQUENCY
+ 5 5 20 20
+ 1.0 60.
+ 1 1
+ VALUE 1.0
+BLANK card ends branch cards
+BLANK card ends switch cards (none, for this case)
+ 1VALUE
+BLANK card ends source cards (here, just user-defined Type-1 --- inline below)
+ 1 { Request for all node voltages outputs (here, just one)
+C Next come the data cards for user-defined Type-1 source on node "VALUE":
+ 3.4
+.8485282
+ 0.0
+-.434314
+ -.2
+-.848528
+ -.4
+-1.56569
+ 3.41
+.8500000
+ 0.0
+-.430000
+ -.2
+ 9999
+ FOURIER ON 10 { For 5 vertical bars
+ 1437.6 1.5 12. VALUE { For 5 vertical bars
+C The preceding produces bars so wide that internal filling by 8 vertical
+C lines (drawn after the outline of a bar is drawn) still leaves lots of
+C empty space. Nonetheless, the effect in MS Word 7 for Win95 is good.
+C The following 15 bars look great in Word 7:
+ FOURIER ON 15 { For 13 vertical bars
+ 1433.0 1.5 25. VALUE { For 13 vertical bars
+ PRINTER PLOT { Toggle to line printer plotting for next subcase as b4 Jan, 97
+BLANK card ending plot cards
+BEGIN NEW DATA CASE
+C 3rd of 4 subcases tests the use of the branch name "SSNAME" (variable
+C SSONLY of the STARTUP file). The Type-91 true-nonlinear R(t) of
+C the first subcase is repeated twice. The first copy is unchanged, so it
+C provides a reference: no nonlinear element current in the phasor soution
+C and flashover so it is present only beginning on step 5. This is node
+C A2. The second copy, on node B2, has low phasor voltage because of
+C an extra 300-ohm resistor added just for the phasor solution. Voltage
+C on step 0 thus is small (21 kV). But the 300-ohm resistor is removed
+C after the phasor solution, so voltage immediately rises to 240 kV on
+C step 1. Then, since the flashover voltage was reduced to 100 kV, we
+C have a flashover, and for step 2 onward of the 2nd copy, the Type-91 R
+C is connected. With both NL elements flashed over, the 2 solutions are
+C close, but not identical. This is because of the time shift of 3 steps
+C (remember, R vs. t characteristics begin at the instant of flashover).
+C Note that the following data puts the "PHASOR" resistor as the last
+C branch, so that the EMTP will erase it completely before the time-step
+C loop, which is most efficient. Identical answers have been obtained by
+C replacing the single 300-ohm branch by two 600-ohm branches which can be
+C mixed anywhere among the branch cards.
+PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
+ 0.5 3.0 { Floating-point misc. data card
+ 1 1 1 1
+ GEN A2 5.E+3 1
+ A2 5.E+3
+91A2 3333. 1
+ 300000. { Original element will be felt on step 5 (step 4 exceeds 300 kV)
+ 300. 0.0
+ 200. 3.0
+ 150. 6.0
+ 150. 1000.
+ 9999.
+ GEN B2 5.E+3 1
+ B2 5.E+3
+ B2 NAME PHASOR 300.
+91B2 3333. 1
+ 100000. { Reduced flashover voltage means R(t) will be present for step 2
+ 300. 0.0
+ 200. 3.0
+ 150. 6.0
+ 150. 1000.
+ 9999.
+BLANK card ends branch cards
+BLANK card ends switch cards (none, for this case)
+14GEN 620000. .06 -50. -1.
+C --------------+------------------------------
+C From bus name | Names of all adjacent busses.
+C --------------+------------------------------
+C GEN |A2 *B2 *
+C A2 |TERRA *TERRA *GEN *
+C B2 |TERRA *TERRA *TERRA *GEN *
+C TERRA |A2 *A2 *B2 *B2 *B2 *
+C --------------+------------------------------
+BLANK card ends source cards
+C Total network loss P-loss by summing injections = 5.560071428571E+07
+C GEN 398528.31800565 620000. 115.28854913735 179.35714285714
+C -474947.5547338 -50.0000000 -137.3955426194 -50.0000000
+ 1 { Request the output of all node voltages
+C Step Time B2 A2 GEN A2 B2
+C TERRA TERRA
+C 0 0.0 21349.73132 199264.159 398528.318 0.0 0.0
+C 1 0.5 240232.7915 240232.7915 480465.583 0.0 0.0
+C 2 1.0 29216.89903 272691.0576 545382.1152 0.0 97.38966342
+C 3 1.5 30079.72946 295489.107 590978.214 0.0 106.163751
+C 4 2.0 29669.33084 307819.3075 615638.615 0.0 111.2599907
+C 5 2.5 28113.16869 33133.37738 618489.7112 110.4445913 112.4526748
+C 6 3.0 25585.44821 30509.93568 599430.5009 107.6821259 109.6519209
+ CALCOMP PLOT
+ 2Plot of both Type-91 voltages
+C 78901234567890123456789012345678901234567890123456789012345678901234567890
+ 143 .3 0.0 3.0 A2 B2 Type-91 nodes Voltages in [v]
+ PRINTER PLOT
+ 143 .6 0.0 3.0 A2 B2 { Axis limits: (0.000, 3.078)
+$ABORT { Illustration of request to abort this subcase but continue with next
+BLANK card ending plot cards --- needed to exit plotting where KILL is checked
+C This and any later cards of data subcase should not be executed due to abort
+$WIDTH, 79 { Switch to narrow output so case summary statistics will be compact
+XXXXX XXXXX XXXXX { 1st of several dummy data cards that should be skipped over
+YYYYY YYYYY YYYYY { 2nd of several dummy data cards that should be skipped over
+ZZZZZ ZZZZZ ZZZZZ { 3rd of several dummy data cards that should be skipped over
+AAAAA aaaaa AAAAA { 4th of several dummy data cards that should be skipped over
+bbbbb BBBBB bbbbb { 5th of several dummy data cards that should be skipped over
+CCCCC ccccc CCCCC { 6th of several dummy data cards that should be skipped over
+BEGIN NEW DATA CASE
+C 4th of 4 subcases is here just to service preceding $ABORT. Next, halt.
+$STOP { Illustrate user-requested termination of execution via input card
+CHANGE PRINTOUT FREQUENCY { 1st of otherwise-legal card that will not be used
+ 10 5 { 2nd of otherwise-legal card that will not be used
+PRINTED NUMBER WIDTH, 13, 2, { 3rd of otherwise-legal card that will not be used
+$MONITOR { 4th of otherwise-legal card that will not be used
+DISK PLOT DATA { 5th of otherwise-legal card that will not be used
+
|