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
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
|
BEGIN NEW DATA CASE
C BENCHMARK DC-5
C Saturable TRANSFORMER problem with the same basic solution as DC-4.
C 1st of 7 subcases illustrates narrow (79-column) output.
$WIDTH, 79, { Request narrow, 80-column LUNIT6 output as an illustration
BEGIN PEAK VALUE SEARCH, -1, { Time of -1 is request to read following card
0.5 1.5 2.5 3.5 { Limit extrema to these 2 subranges
PRINTED NUMBER WIDTH, 10, 2, { Request max precision of 6 numbers in 80 bytes
.010 6.0
1 1 1 1 1 -1
5 5 20 20 100 100
C The following $VINTAGE restoration is added 3 November 2001 to illustrate
C the associated rejection message. The request is meaningless since there is
C no preceding $VINTAGE request. For correct use, see later subcase.
$VINTAGE, -1 { Cancel preceding $VINTAGE request, returning the value to zero
GEN GENT 1.E-3 3
LOAD 1.0E3 3
TRANSFORMER .005 30.TRANFF 1.E4 3
.005 30.
.01 40.
.02 45.
0.1 50.
5.0 100.
9999
1GENT OPEN 5.0 5.E4 50.
2LOADFF 20. 2.E5 100.
OPEN 1.E7
TRANSFORMER TRANFF TRAN 3
1GENT
2LOAD
BLANK card finishing all branch cards
BLANK card ending non-existent switch cards
14GEN 70. .1591549 -1.
C --------------+------------------------------
C From bus name | Names of all adjacent busses.
C --------------+------------------------------
C GEN |GENT *
C GENT |GEN *TRANFF*TRAN *
C LOAD |TERRA *TERRA *
C TRANFF |GENT *OPEN *OPEN *OPEN *
C OPEN |TERRA *TRANFF*TRANFF*TRANFF*
C LOADFF |TERRA *
C TRAN |TERRA *TERRA *TERRA *GENT *
C TERRA |LOAD *LOAD *OPEN *LOADFF*TRAN *TRAN *TRAN *
C |
C --------------+------------------------------
BLANK card terminating program source cards
C Total network loss P-loss by summing injections = 8.286936503209E+00
C Inject: GEN 70. .25769841013987 8.2869365032093 9.0194443548954
C Inject: 0.0 -23.2507996 3.5604859026308 0.9187857
C ---- Initial flux of coil "TRANFF" to "OPEN " = 3.08658228E-02
C ---- Initial flux of coil "TRAN " to " " = -1.13293881E+01
C Step Time GEN LOAD TRANFF TRAN LOAD TRAN
C GENT TERRA OPEN TERRA
C
C GENT LOADFF OPEN TRANFF GEN TRANFF
C OPEN
C
C TRAN GEN LOAD
C TERRA GENT TERRA
C 0 0.0 .2368E-3 116.1389 .0185382 63.72954 116.1389 63.72954
C 69.99976 .0370763 69.98119 69.99973 70. .5144E-5
C -.001888 .2367696 .1161389
C 1 .01 .2378E-3 116.583 .0182286 63.83965 116.583 63.83965
C 69.99626 .0364571 69.978 69.99623 69.9965 .5175E-5
C -.001782 .237775 .116583
1 { Request for all node voltage outputs
C 600 6.0 .1743E-3 86.66574 .0264244 54.92308 86.66574 54.92308
C 67.21171 .0528488 67.18516 67.21158 67.21189 .4076E-5
C -.004558 .1742727 .0866657
C Variable maxima : .3281E-3 123.4893 .0360051 61.35953 123.4893 61.35953
C 67.21171 .0720101 67.18516 67.21158 67.21189 .6001E-5
C .2904095 .3281234 .1234893
C Times of maxima : 1.49 0.5 5.25 0.5 0.5 0.5
C 6.0 5.25 6.0 6.0 6.0 .54
C 1.49 1.49 0.5
C Variable minima : -.349E-3 -122.626 -.033186 -63.9283 -122.626 -63.9283
C -69.9997 -.066372 -69.9811 -69.9996 -69.9999 -.6E-5
C -.32526 -.349375 -.122626
C Times of minima : 4.67 3.56 2.51 3.35 3.56 3.35
C 3.14 2.51 3.14 3.14 3.14 3.68
C 4.76 4.67 3.56
PRINTER PLOT
C Axis limits: (-2.544, 2.624)
193 .5 0.0 4.0 LOAD GEN GENT
BLANK card terminating plot cards
$WIDTH, 132, { Done with 80 columns, so return to wide, 132-column LUNIT6 output
BEGIN NEW DATA CASE
C 2nd of 7 subcases. This is the former (before January, 1987) DC-43.
C Test of [A], [R] usage ("USE AR"). For DIAGNOSTIC, see DCPRINT-14.
C [A],[R] data cards used here were punched by 3rd subcase of DCNEW-8.
C Solution can be checked by hand, since it is really just 1 big series
C circuit, with an ideal transformer separating the primary & secondary.
C The turns ratio is 2:1, with the primary being the high-voltage side.
C Transformer has bottoms of primary and secondary grounded. The primary
C has a 1-volt, 1-[rad/sec] sinusoidal voltage source, that is connected
C to the transformer through a reactance X = j0.1. Values X = 0.1 ohm
C and R = 0.2 ohm are next, coming from the [A], [R] matrices when
C reflected entirely to the primary of the ideal transformer. Finally,
C the secondary is closed through a 1-ohm resistor. Reflecting all to
C the primary, R = 4.2 ohm, X = j0.2, or Z-tot = 4.2047592083257 at
C an angle of 2.726226 degrees. Reciprocating gives phasor current
C I-source = .23782576065627 at an angle of -2.726226 degrees. This
C can be checked exactly. Then, correctness of the time-step solution
C follows by smoothness (there is no discontinuity on the first step).
PRINTED NUMBER WIDTH, 25, 2, { Request 25-column precision for 4 output columns
C Apollo will not deliver more than about 16 digits of output, since
C this is the precision limit of REAL*8. Hence the preceding width
C request will in fact have more than two blank separator columns.
.100 6.0.1591549 { Note XOPT = 1 / (2 * Pi ), and w = 1.0 [rad/sec]
1 1 1 1 1 -1
5 5
GEN PRIM NAME Excite 0.1 { Branch named for later current output
USE AR { Col. 3-8 request switches to [A], [R] usage
1SEC 40. .025
2PRIM -20. 10. 0.1
USE RL { Col. 3-8 request returns to [R], [L] usage
SEC 1.0
BLANK card ending branches
BLANK card ending non-existent switches
14GEN 1.0 .15915 0.0 -1.
C The following forward and reverse flows are for 1st phasor branch:
C GEN 1.0 1.0 .23755659127612 .23782576065627
C 0.0 0.0 -.0113118686787 -2.7262268
C
C PRIM .99886884795867 .99915127584094 -.2375565912761 .23782576065627
C -.0237549277475 -1.3623416 .01131186867874 177.2737732
BLANK card ending source cards
C Total network loss P-loss by summing injections = 1.187782956381E-01
C Inject: GEN 1.0 1.0 .23755659127612 .23782576065627 .11877829563806
C Inject: 0.0 0.0 -.0113118686787 -2.7262268 .00565593433937
C
C Step Time SEC PRIM GEN GEN
C PRIM
C 0 0.0 .4751131825522436 .9988688479586737 1.0 .2375565912761218
C 1 0.1 .4749983247294277 .9962504852347109 .9950044753394702 .2374991623647137
C 2 0.2 .470139617191015 .9836805108498354 .9800678118911468 .2350698085955073
C 3 0.3 .460583675647767 .961282487768105 .9553392425962317 .2302918378238831
1 XXXX { Cols. 1-2 requests all node voltage outputs; A6 nonblank ==> another
-1Excite { Cols 1-2 type code "-1" ===> request for branch current outputs
BLANK card terminating selective (and all node voltages) output requests
C The secondary current need not be monitored, since it is twice
C the primary current, and is equal to the voltage of node SEC.
C 55 5.5 .3206618611460799 .690969577591206 .708549242853354 .1603309305730407
C 60 6.0 .4498323687101901 .9523830926143848 .9601182008522722 .2249161843550953
C Max : .4751131825522436 .9988688479586737 1.0 .2375565912761218
C Times of max: 0.0 0.0 0.0 0.0
C Min -.4756234626742181 -.9985549285218509 -.999131142217714 -.2378117313371089
C T-min 3.200000000000002 3.200000000000002 3.100000000000001 3.200000000000002
PRINTER PLOT
C Following plot shows in-phase V-prim and V-sec, with ratio about 2 to 1:
143 1. 0.0 6.0 PRIM SEC { Axis limits: (-9.986, 9.989)
C Following plot shows gen current, I = .23782576065627 /__ -2.726226 degrees
193 1. 0.0 6.0 GEN PRIM { Axis limits: (-2.378, 2.376)
BLANK card ending non-existent plot cards
BEGIN NEW DATA CASE
C 3rd of 7 subcases. Test of [A], [R] usage ("USE AR"), with the
C branch cards punched by the 3rd subcase of DCNEW-8. Trivial 1-phase
C problem. For identical network using [R], [L], see the 4th subcase
PRINTED NUMBER WIDTH, 15, 2, { Request more than enough precision for 3 outputs
.0005 .050 50.
1 1 1 1 2 -1
5 5
S1 25.782 1
$VINTAGE, 1
$UNITS, 0., 0. {XOPT, COPT
USE AR
1P1 609.14725451514 .05780196523054
2S1 -355.0194488376 0.0
207.17315316786 .16973833711262
C $VINTAGE, 0
C 3 November 2001, replace the preceding by just-added cancellation of
C the preceding $VINTAGE,1. Of course, the starting value is zero:
$VINTAGE, -1 { Cancel preceding $VINTAGE request, returning the value to zero
BLANK card ending branches
BLANK card ending non-existent switches
14P1 311. 50. 000.0 0. -1.
BLANK card ending source cards
C Total network loss P-loss by summing injections = 5.418483002156E+03
C P1 311. 311. 34.845549853093 34.95612473617 5418.483002156
C 0.0 0.0 -2.77818502017 -4.5584667 432.0077706365
C Step Time S1 P1 S1
C TERRA
C 0 0.0 524.252284827 311. 20.3340425424
C 1 .5E-3 522.548672234 307.171073925 20.2679649458
C 2 .1E-2 507.987759921 295.778576568 19.7031944737
C 3 .0015 480.917065702 277.103029023 18.6532102126
1
C 95 .0475 -349.18551236 -219.91020895 -13.543771327
C 100 .05 -524.24501851 -311. -20.333760706
C Variable maxima : 524.252284827 311. 20.3340425424
C Times of maxima : 0.0 0.0 0.0
C Variable minima : -524.24501851 -311. -20.333760706
C Times of minima : .01 .01 .01
PRINTER PLOT
144 4. 0.0 20. S1 { Axis limits: (-5.242, 5.243)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 4th of 7 subcases. Test of [R], [wL] usage ("USE RL"), with the
C branch cards punched by the 3rd subcase of DCNEW-8. This is the same
C simulation as the preceding subcase, only [R],[L] replaces [A],[R]
PRINTED NUMBER WIDTH, 15, 2, { Request more than enough precision for 3 outputs
.0005 .050 50.
1 1 1 1 2 -1
5 5
S1 25.782 1
$VINTAGE, 1,
$UNITS, 50., 0., {XOPT, COPT
USE RL
1P1 .05780196523054 406.40657142442
2S1 0.0 696.43307921366
.16973833711262 1194.949458531
$VINTAGE, 0
BLANK card ending branches
BLANK card ending non-existent switches
14P1 311. 50. 000.0 0. -1.
BLANK card ending source cards
C Total network loss P-loss by summing injections = 5.418483002158E+03
C P1 311. 311. 34.845549853106 34.956124736168 5418.483002158
C 0.0 0.0 -2.778185019985 -4.5584667 432.00777060763
C Step Time S1 P1 S1
C TERRA
C 0 0.0 524.252284827 311. 20.3340425424
C 1 .5E-3 522.548672234 307.171073925 20.2679649458
C 2 .001 507.987759921 295.778576568 19.7031944737
C 3 .0015 480.917065702 277.103029023 18.6532102126
1
C 95 .0475 -349.18551236 -219.91020895 -13.543771327
C 100 .05 -524.24501851 -311. -20.333760706
C Variable maxima : 524.252284827 311. 20.3340425424
C Times of maxima : 0.0 0.0 0.0
C Variable minima : -524.24501851 -311. -20.333760706
C Times of minima : .03 .01 .03
PRINTER PLOT
144 4. 0.0 20. S1 { Axis limits: (-5.242, 5.243)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 5th of 7 subcases is used to illustrate RENUMBER WITHOUT COUPLING
C Same data as 1st subcase, but with two important changes. First, output is
C wide (not 80-column). Second and critically important, the transient node
C renumbering is different because of the following request that is introduced
C in the January, 1999, newsletter. Output variables are the same, but the
C "1" request for all node voltage outputs will change the order of occurance
C of these variables in the output vector.
RENUMBER WITHOUT COUPLING { Upset node renumbering; change output of all V-node
PRINTED NUMBER WIDTH, 10, 2, { Request max precision of 6 numbers in 80 bytes
.010 .030 { Take only 3 steps (plenty to illustrate the effect)
1 1 { Note no extrema (values would not agree, anyway)
GEN GENT 1.E-3 3
LOAD 1.0E3 3
TRANSFORMER .005 30.TRANFF 1.E4 3
.005 30.
.01 40.
.02 45.
0.1 50.
5.0 100.
9999
1GENT OPEN 5.0 5.E4 50.
2LOADFF 20. 2.E5 100.
OPEN 1.E7
TRANSFORMER TRANFF TRAN 3
1GENT
2LOAD
BLANK card finishing all branch cards
BLANK card ending non-existent switch cards
14GEN 70. .1591549 -1.
BLANK card terminating program source cards
1 { Request for all node voltage outputs
BLANK card terminating plot cards
BEGIN NEW DATA CASE
C 6th of 7 subcases is identical to the 2nd. But because KOMPAR > 4
C (see the $DEPOSIT use immediately below), KOMPAR is interpreted
C to be the output width of phasor and dT-loop numbers. It works. Not
C only are the numbers reduced in precision, the dT-loop headings also
C are adjusted to match. I.e., the 25 of the PRINTED NUMBER WIDTH
C declaration below is overridden. Although operation was described in
C the April, 1998, newsletter, an illustration (this data) did not
C exist prior to 5 December 2001.
C Add Type-15 USRFUN source on 28 July 2002 to illustrate this new
C feature from Orlando Hevia. The USRFUN source provides a higher-
C level replacement for ANALYTIC SOURCES USAGE that is illustrated
C in DC-6. The same node name NODE1 will be used. The same signal
C is generated.
$DEPOSIT, KOMPAR=11, { Illustrate reduced field width of phasor & dT-loop #-s
PRINTED NUMBER WIDTH, 25, 2, { Request 25-column precision for 4 output columns
C Apollo will not deliver more than about 16 digits of output, since
C this is the precision limit of REAL*8. Hence the preceding width
C request will in fact have more than two blank separator columns.
.100 6.0.1591549 { Note XOPT = 1 / (2 * Pi ), and w = 1.0 [rad/sec]
1 1 0 1 0 -1
10 5
GEN PRIM NAME Excite 0.1 { Branch named for later current output
USE AR { Col. 3-8 request switches to [A], [R] usage
1SEC 40. .025
2PRIM -20. 10. 0.1
USE RL { Col. 3-8 request returns to [R], [L] usage
SEC 1.0
NODE1 1.0 { Dummy resistor serves to anchor USRFUN source
BLANK card ending branches
BLANK card ending non-existent switches
14GEN 1.0 .15915 0.0 -1.
C Each Type-15 USRFUN source must be provided by the user in the user-supplied
C FUNCTION USRFUN. Here, illustrate type code 5 (CREST in columns 11-20) that
C is identical to the ANALYT function of DC-6. There is a difference, however
C Whereas ANALYT internally zeroed the function after T = .75 seconds, here
C no such logic has been built into USRFUN. Instead, use T-stop of 71-80 to
C zero the function beginning on time step number 8. Speaking of time steps,
C modify the printout frequency to show the first 10 steps. About columns,
C the request word "USRFUN " (note 2 trailing blanks) must occupy 43-50.
C < BUS> < Index > Request T-start T-stop
15NODE1 5.0 USRFUN 0.0 0.75
BLANK card ending source cards
C SEC .475113189 .475651527 .475113189 .475651527 .113122188 .113122188
C -.02262374 -2.7262268 -.02262374 -2.7262268 0.0 0.0
C
C TERRA 0.0 0.0 -.47511319 .475651527 0.0
C 0.0 0.0 .022623738 177.2737732 0.0
C Total network loss P-loss by summing injections = 1.187782971625E-01
C Note all of the white space surrounding the preceding phasor valus.
C The precision is reduced but column positioning and headings are
C unchanged. This is unlike dT-loop output, for which the headings,
C too, are modified so as to maintain the same inter-# separation.
1 XXXX { Cols. 1-2 requests all node voltage outputs; A6 nonblank ==> another
-1Excite { Cols 1-2 type code "-1" ===> request for branch current outputs
C First 4 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 SEC PRIM GEN NODE1 GEN
C PRIM
C 0 0.0 .47511319 .99886885 1.0 0.0 .23755659
C 1 0.1 .47499833 .99625049 .99500448 40. .23749917
C 2 0.2 .47013962 .98368051 .98006781 80. .23506981
C 3 0.3 .46058368 .96128249 .95533924 95. .23029184
C 4 0.4 .44642603 .92928025 .92106583 85. .22321301
C 5 0.5 .42780811 .88799352 .87759001 75. .21390405
C 6 0.6 .40491594 .8378348 .82534614 65. .20245797
C 7 0.7 .37797823 .77930524 .76485619 55. .18898912
C 8 0.8 .34726412 .7129896 .69672453 0.0 .17363206
BLANK card terminating selective (and all node voltages) output requests
C 60 6.0 .44983237 .95238309 .9601182 0.0 .22491619
PRINTER PLOT
C Following plot shows in-phase V-prim and V-sec, with ratio about 2 to 1:
143 1. 0.0 6.0 PRIM SEC { Axis limits: (-9.986, 9.989)
BLANK card ending non-existent plot cards
BEGIN NEW DATA CASE
C 7th of 7 subcases is fundamentally different from preceding subcases
C in that it involves a frequency scan rather than time simulation. But
C there is similarity to the preceding subcase in the user-supplied
C source code is involved. This tests USER10 as built into the UTPF
C beginning 26 October 2003.
PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns)
C FMIMFS--DELFFS--FMAXFS--NPD-----
C 11111111222222223333333344444444
FREQUENCY SCAN 25.0 25.0 200.0 0 { F = 25, 50, ... 200
0.00010 0.0
1 1
RRR 10. { 10 ohm resistor so V = 10 * source curr } 1
BLANK card ending branches
BLANK card ending switches (none)
BOTH POLAR AND RECTANGULAR { Request for (in order): mag, angle, real, imag
C Preceding is one of 3 alternatives. The other two are, after commented:
C POLAR OUTPUT VARIABLES { 2nd of 3 alternatives for output gives mag, angle
C RECTANGULAR OUTPUT VARIABLES { 3rd of 3 alternative outputs gives real, imag
C Following ANALYTIC source is user-defined. The function name HYPERB must
C be defined in USER10. Note T-start of cols. 61-70 is left blank as an
C illustration that ATP automatically will set this to value -1.0 Amplitude
C (columns 11-20) and frequency (columns 21-30) likewise will be defined
C internally. Only columns 1-10, plus the text in 43-57, is required:
C Node V? Request> -Name-
14RRR -1 ANALYTIC HYPERB
C Preceding is defined as follows: Real part = 100. * Pi / Omega
C Imaginary part = ( Real part + 1.0 ) / 2
BLANK card ending source cards
RRR { Node voltage will be 10 times the source current, note.
BLANK card ending requests for node voltage output
C Step F [Hz] RRR RRR RRR RRR RRR RRR RRR RRR
C TERRA TERRA TERRA TERRA
C 1 25. 25. 36.8699 20. 15. 2.5 36.8699 2.0 1.5
C 2 50. 14.14214 45. 10. 10. 1.414214 45. 1.0 1.0
C 3 75. 10.67187 51.34019 6.666667 8.333333 1.067187 51.34019 .6666667 .8333333
C 4 100. 9.013878 56.30993 5.0 7.5 .9013878 56.30993 0.5 .75
C 5 125. 8.062258 60.25512 4.0 7.0 .8062258 60.25512 0.4 0.7
C 6 150. 7.45356 63.43495 3.333333 6.666667 .745356 63.43495 .3333333 .6666667
C 7 175. 7.034898 66.03751 2.857143 6.428571 .7034898 66.03751 .2857143 .6428571
C 8 200. 6.731456 68.19859 2.5 6.25 .6731456 68.19859 .25 .625
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
|
