1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
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
|
