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
|
BEGIN NEW DATA CASE
C 1st of 3 subcases is added 18 October 2001. It is unrelated to following
C 2nd. Since a DO KNT loop is involved, data must be in 1st subcase.
C 1st of 2 subcases illustrates "simultaneous usage of $PARAMETER
C blocks and DO loops, which are used to serialize power network
C branches." This is the way Prof. Juan Martinez of UPC in Barcelona,
C Spain, described his latest observation of inadequacy of the new ATP
C tools, in an E-mail message dated 16 October 2001. The attached small
C example NEWDIS2.DAT served as the model for the present addition.
C Among other changes, /-card data sorting was removed for simplicity.
C It had nothing to do with the program error, which was corrected the
C following day. Before correction, execution died early in branch input.
$PARAMETER
AAAAAAAA='$ENABLE '
BBBBBBBB='$DISABLE'
BLANK card ends block of parameter definitions
.000100 .010 { dT and T-max will result in 100 steps (Martinez had 400)
20 -1 { Minimal printout and no plotting demonstrate simulation
DO KNT=1, 5 /OUTPUT
$PARAMETER
BU1___ = KNT - 1. SERIALIZE 'BUS000'
BU2___ = KNT SERIALIZE 'BUS000'
BLANK card ending symbols definitions
BU1___BU2___ 0.075 0.5570
BU2___ 1000. 2387.3
ENDDO KNT
C The following 3 data cards will be ignored since BBBBBBBB will be replaced
C by $DISABLE (see the first $PARAMETER block). Yet, these 3 lines
C produce substantial output in the .LIS file because the $DISABLE is
C dynamic. Each line of the disk file must be read since ATP has no way of
C knowing that one of them does not end the discarding.
BBBBBBBB
$INCLUDE dofile.dat
$ENABLE
BLANK card ending BRANCH cards
SOURCEBUS000 MEASURING 1
C The following data card will be ignored since AAAAAAAA will be
C replaced by $ENABLE by the first $PARAMETER block:
AAAAAAAA
BLANK card ending SWITCH cards
14SOURCE 20412.4 50.0 0.0 -1.0
BLANK card ending SOURCE cards
BUS000BUS001BUS004BUS005
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 BUS000 BUS001 BUS004 BUS005 SOURCE
C BUS000
C *** Phasor I(0) = 6.5184598E+01 Switch "SOURCE" to "BUS000" closed in the steady-state.
C 0 0.0 20412.4 20398.9415 20374.7213 20372.0305 65.1845978
C 20 .002 16513.9785 16507.6359 16496.2181 16494.9494 81.5202569
C 40 .004 6307.7785 6310.97466 6316.72113 6317.3591 66.7198316
C 60 .006 -6307.7785 -6296.2642 -6275.5483 -6273.2473 26.4355142
C 80 .008 -16513.978 -16498.544 -16470.771 -16467.686 -23.945918
C 100 .01 -20412.4 -20398.941 -20374.72 -20372.029 -65.18067
BLANK card ending OUTPUT cards
BLANK card ending PLOT cards
BEGIN NEW DATA CASE
C BENCHMARK DCNEW-2
C 2nd of 2 subcases was the 1st prior to 18 October 2001.
C Test of automatic U.M. initialization using 3-phase synchronous mode of
C operation. The steady-state lasts 50 msec until a TACS "fault" occurs.
C The machine data comes from the text by Fouad and Anderson (real units)
C Note: WSM and Yin change comment cards in Leuven, October 22, 1990,
C following modifications of code that change answers in 3rd or 4th
C decimal place (peak values of PRINTER PLOTs are barely affected).
C In response to a request for specifics, G. Papagiannis of Aristotle
C University of Thessaloniki in Greece posted the following on February
C 18, 1998 (public E-mail from Fargo): A. A. Fouad and P. M. Anderson;
C "Power System Control and Stability," Iowa State University Press in
C Ames, Iowa, USA, 1982. "An excellent treatment of synhronous machine
C modelling" was the recommendation. Neil Fitzgerald of Electrowatt
C Engineering (UK) Ltd concurred later that same day, adding that the
C book "was reprinted in 1993 by the IEEE Press and I think it costs
C around $39.00 from the IEEE." He agreed on the recommendation. As
C for the "P. M.", this is Dr. Paul M. of PMA (Power Math Associates)
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
ABSOLUTE U.M. DIMENSIONS, 20, 2, 50, 100
POWER FREQUENCY, 60, { Corrects possible 50-Hz declaration of European STARTUP
.000200 .070
1 1 1 1 1 -1
5 5 20 20 100 100
TACS HYBRID
STEP1 +BAT1 -TFAULT
STEP2 +BAT2 -TCLEAR
FIRE -STEP0 +STEP1 -STEP2
C --------- TACS UNITY SIGNAL SOURCE
11STEP0 1.0 -1.0
C --------- EMTP OUTPUT TO TACS
90BAT1
90BAT2
90TFAULT
90TCLEAR
BLANK card ends all TACS data cards
C --------- EMTP CONNECTIVITY
BUSA0 1.0E+8
BUSB0 1.0E+8
BUSC0 1.0E+8
BUSA5 BUSA1 0.001 1
BUSB5 BUSB1 0.001 1
BUSC5 BUSC1 0.001 1
C --------- TRANSMISSION LINES (0.001 OHMS AND 1.0 MILLI HENRY)
BUSA1 BUSA0 1.0E-2 1.0 1
BUSB1 BUSB0 BUSA1 BUSA0 1
BUSC1 BUSC0 BUSA1 BUSA0 1
C --------- RESISTANCES FROM BUS1 TO SWITCH (BUS2 TO GROUND)
BUSA1 BUSA2 0.0001 1
BUSB1 BUSB2 0.0001
BUSC1 BUSC2 0.0001
C --------- EMTP CONNECTIVITY REQUIREMENT
BUSFS 1.0E+6
BUSFS BUSF0 1.0E-5
C --------- RESISTANCES IN PARALLEL TO SWITCHES (BUS2 TO GROUND)
BUSA2 1.0E+9
BUSB2 1.0E+9
BUSC2 1.0E+9
C --------- CIRCUIT FOR SWITCH CONTROL
BAT1 1.0
BAT2 1.0
BAT1 TFAULT 1.0
BAT2 TCLEAR 1.0
C --------- MECH NETWORK
BUSMSM 1.E+09 1
BUSMS BUSMSM 1.0E-3 1
BLANK card bounds the last branch card
C --------- TIME OF FAULT AND TIME OF CLEARING
TFAULT .05001 { Fault is applied 1/2 DELTAT after 50 msec
TCLEAR 3600.0
C --------- TACS CONTROLLED SWITCHES
11BUSA2 FIRE
11BUSB2 FIRE
11BUSC2 FIRE
BUSA9 BUSA5 -1.0 1.0 1
BUSB9 BUSB5 -1.0 1.0 1
BUSC9 BUSC5 -1.0 1.0 1
BLANK card ends the last electric network switch card
C --------- EXCITATION SOURCE VOLTAGE (ACTUAL VALUE SET BY SS)
14BUSFS 0.001 0.00001 -1.0
C --------- INFINITE BUS SOURCE VOLTAGES
14BUSA0 3000.0 60.0 0.0 -1.0
14BUSB0 3000.0 60.0 -120.0 -1.0
14BUSC0 3000.0 60.0 +120.0 -1.0
C --------- BATTERIES FOR SWITCH CONTROL
14BAT1 2.0 0.00001 -1.0
14BAT2 3.0 0.00001 -1.0
C --------- MECH INPUT TORQUE (ACTUAL VALUE SET BY SS) :
14BUSMS -1 0.0000001 0.00001 -1.0
C --------- SYNCHRONOUS MACHINE SPECIFICATIONS
C --------- UM DATA
19 UM
1
BLANK CARD ENDING CLASS 1 UM DATA
C --------- UM-1 MACH-TABLE
1 2 1111BUSMSM 1 0.0373
0.0057815
0.0055850
3030.0 15.0 BUSFS BUSMS
C --------- UM-1 COIL-TABLE
BUSA9 1
0.001542 0.0005595 BUSB9 1
0.001542 0.0005595 BUSC9 1
0.0010 0.0003767 BUSF0 1
0.01834 0.0002075 1
0.07560 0.0001335 1
BLANK card ends all U.M. data cards
C --------------+------------------------------
C From bus name | Names of all adjacent busses.
C --------------+------------------------------
C BUSA0 |TERRA *BUSA1 *
C BUSB0 |TERRA *BUSB1 *
C BUSC0 |TERRA *BUSC1 *
C BUSA5 |BUSA1 *BUSA9 *
C BUSA1 |BUSA0 *BUSA5 *BUSA2 *
C BUSB5 |BUSB1 *BUSB9 *
C BUSB1 |BUSB0 *BUSB5 *BUSB2 *
C BUSC5 |BUSC1 *BUSC9 *
C BUSC1 |BUSC0 *BUSC5 *BUSC2 *
C BUSA2 |TERRA *TERRA *BUSA1 *
C BUSB2 |TERRA *TERRA *BUSB1 *
C BUSC2 |TERRA *TERRA *BUSC1 *
C BUSFS |TERRA *BUSF0 *
C BUSF0 |BUSFS *
C < < Etc. (many more not shown) > >
BLANK card ending all electric network sources
C Total network loss P-loss by summing injections = 7.248075223468E+04
C Total network loss P-loss by summing injections = 1.461915999075E+15
C Total network loss P-loss by summing injections = 1.461915980400E+15
C Total network loss P-loss by summing injections = 1.461915980711E+15
C Gen 1: BUSMSM 376.99111843078 376.99111843078 -24959.40530999 24959.416549762
C Gen 1: 0.0 0.0 23.687050562614 179.9456250
BUSA1 BUSA0 BUSMSM
C Step Time BUSA1 BUSA0 BUSMSM BUSA9 BUSB9
C BUSA5 BUSB5
C
C BUSC5 BUSA1 BUSB1 BUSC1 BUSA1
C BUSC1 BUSA0 BUSB0 BUSC0 BUSA2
C
C UM-1 UM-1 UM-1 UM-1 UM-1
C OMEGM THETAM IPA IPB IPC
C *** Phasor I(0) = 2.0727726E+03 Switch "BUSA9 " to "BUSA5 " closed
C *** Phasor I(0) = -8.1575149E+02 Switch "BUSB9 " to "BUSB5 " closed
C *** Phasor I(0) = -1.2570211E+03 Switch "BUSC9 " to "BUSC5 " closed
C 0 0.0 2924.682776 3000. 376.9911184 2072.772625 -815.751494
C -1257.02113 2072.772622 -815.751493 -1257.02113 .2928573E-5
C 376.9911184 1.174129748 2072.761251 -815.747018 -1257.01423
C 1 .2E-3 2857.218771 2991.476701 376.9911184 2047.694642 -668.660479
C -1379.03416 2047.69464 -668.660478 -1379.03416 .2833076E-5
C 376.9911184 1.174129748 2047.694642 -668.660479 -1379.03416
C 2 .4E-3 2773.75849 2965.955234 376.9911185 2010.99049 -517.772984
C -1493.21751 2010.990487 -517.772984 -1493.2175 .2755769E-5
C 376.9911184 1.174129748 2010.99049 -517.772984 -1493.21751
BLANK card ending selective variable output request (node voltages, here)
C 350 .07 .7006894039 927.0509831 376.7455301 2408.782669 10873.00302
C -13281.7857 -4598.11137 -807.600679 5405.712048 7006.894039
C 376.7757426 1.169603437 2408.782669 10873.00302 -13281.7857
C Variable max: 3027.62917 3000. 376.9914104 7182.444042 21231.90574
C 2087.531085 11218.73045 2087.574976 12006.72539 11206.93469
C 377.0003725 1.174129748 7182.444042 21231.90574 2087.531085
C Times of max: .016 0.0 .0504 .0532 .0578
C .0108 .0624 .0386 .0568 .0534
C .0652 .0352 .0532 .0578 .0108
PRINTER PLOT
194 8. 0.0 96. UM-1 TQGEN { Plot limits: (-0.913, 1.466)
19415. 0.0 75. UM-1 OMEGM { Plot limits: (0.000, 3.770)
194 8. 0.0 96. UM-1 IPA { Plot limits: (-1.335, 0.718)
BLANK card ending batch-mode plot cards
BEGIN NEW DATA CASE
C 3rd of 3 subcases is added 12 November 2001. It is unrelated to preceding.
C Data to be used to test counting of actual (minimal) table sizes. This is
C ESTIMATE ACTUAL TABLE SIZES (EATS) which entered the UTPF 12 November 2001.
C Data is simple, and was obtained from DC-7 by canibalization. Prior to
C addition here, data was stored separately as COUNT.DAT This is just the
C first illustration of EATS so keep it simple by omission of batch-mode
C plot requests and other unrelated features.
PRINTED NUMBER WIDTH, 9, 2, { Full precision on each of 8 columns of printout
.000025 .0007 60. 60.
1 1 1 1 -1
5 5 10 10 20 1
ESTIMATE ACTUAL TABLE SIZES { New request must follow immediately the misc. data
2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0
2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0
2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0
10 10 10 10
0 0 10
10 10
C The preceding data is similar to NEW LIST SIZES in that there is one entry
C for each of 30 lists, with 10 per card. There are two such sets: first
C the multiplication factors, and second the additive constants.
GENA SWA 89352. 3
BRANCH NAME:2nd LIN
GENB SWB GENA SWA 1
GENC SWC GENA SWA
93RECA NAME 1st NL .559 300. 3
0.0 0.0
.5590 300.
.9344 400.
9999
BRANCH NAME:2nd NL ! { Hold lower case of this 2nd, but not following 3rd
93RECB RECA .559 300. 4
-1SENDA RECA .137681.07755.6806124.27
-2SENDB RECB .03009.450289.5000124.27
-3SENDC RECC
-1RECA ENDA SENDA RECA
-2RECB ENDB
-3RECC ENDC
BLANK card ends all branch cards
SWA SENDA -1. 1.0
SWB SENDB 0.0 1.0
SWC SENDC 0.0 1.0
BLANK card ends all switch cards
14GENA 188000. 60. 0.0 -1.
14GENB 188000. 60. 120. -1.
14GENC 188000. 60. -120. -1.
BLANK card ending source cards
RECA RECB RECC
BLANK card ending output variable requests (just node voltages, here)
PRINTER PLOT
BLANK card terminating plot cards
BEGIN NEW DATA CASE
BLANK
EOF
10 June 2002, WSM adds output to the screen in case of DISK use.
Without any EATS, this is simple as should be illustrated in the October
(or later) newsletter. But with EATS, there are variations depending upon:
1) the subcase number; and 2) whether the order for EATS comes from
STARTUP, or is buried in data, or both. Consider first the case of no
STARTUP request. I..e., assume FLZERO > 0. Then the screen will show:
---- Begin next subcase number KNTSUB = 2
---- Begin next subcase number KNTSUB = 3
Explicit EATS request now is being processed.
I.e., the exceptional 3rd subcase will have two lines rather than
the usual one. On the other hand, if FLZERO < 0, the screen
should show:
---- Begin EATS for subcase number KNTSUB = 1
---- Begin EATS for subcase number KNTSUB = 2
---- Both implied and explicit EATS. KNTSUB = 3
There is yet another variation. Suppose FLZERO > 0 with
EATS data embedded in the 1st subcase. Output will differ from the
preceding as follows:
---- Explicit EATS request in subcase KNTSUB = 1
|
