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
|
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 <Force .PL4 type> 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 <Force .PL4 type>
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
|
