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
|
BEGIN NEW DATA CASE
C BENCHMARK DC-45
C 1st half of split version of DC-33 . Rather than use TACS to integrate
C flux during the simulation, that operation is now posponed until the post-
C processing of compansion case DC-46. Connection is via LUNIT4 plot file.
C Total of 4 subcases.
C DISK PLOT DATA { Toggle the Apollo default of LUNIT4 = -4 to +4 (use disk)
C $CLOSE, UNIT=4 STATUS=DELETE { Destroy empty date/time plot file of "SYSDEP"
C $OPEN, UNIT=4 FILE=DC45TO46.pl4 ! { Expected to be C-like (for L4BYTE = 1 )
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
.000050 .050
1 1 1 1 1 -1 0 2
5 5 20 20 100 100
96 GEN 8888. 1.E-9 1
1.0 -0.7
2.0 0.9
3.5 1.0
9999.
BLANK card ending branch cards
BLANK card ending (here non-existent) switch cards
14GEN 377. 60. -1.
C --------------+------------------------------
C From bus name | Names of all adjacent busses
C --------------+------------------------------
C GEN |TERRA *TERRA *
C TERRA |GEN *GEN *
C --------------+------------------------------
BLANK card ending source cards
C Total network loss P-loss by summing injections = 7.106450000000E-06
C ---- Initial flux of coil " " to "GEN " = 1.00000000E-09
C Inject: GEN 377. 377. .377E-16 2.6786345332877 .710645E-14
C Inject: 0.0 0.0 -2.678634533288 -90.0000000 504.92260952473
C Step Time GEN TERRA
C GEN
C 0 0.0 377. 0.0
C 1 .5E-4 376.9330268 -.523564601
C 2 .1E-3 376.7321312 -1.04694318
C 3 .15E-3 376.3973844 -1.02374791
C 4 .2E-3 375.9289054 -1.04551661
C 5 .25E-3 375.3268607 -1.06725433
C 10 .5E-3 370.3222935 -1.17520854
1
C 800 .04 -304.999407 1.070136441
C 900 .045 -116.499407 2.765761443
C 1000 .05 377. -1.4375
C Variable maxima : 377. 3.499909256
C Times of maxima : 0.0 .0125
C Variable minima : -377. -3.49990923
C Times of minima : .025 .0375
PRINTER PLOT
194 4. 0.0 20. GEN { Axis limits: (-3.500, 3.500)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 2nd of 4 subcases is unrelated to 1st (no hysteresis). But there is a
C nonlinear inductor (Type-93, true-nonlinear element). This example
C illustrates TACS control of the optional shunt current source that is
C declared by the second data field (columns 21-26) of the TACS CONTROL
C declaration. Since the controlling TACS signal is identically unity,
C this injects an extra amp of current into node "TRAN" beginning on
C the first time step. This would produce a discontinuity in the branch
C current (GEN, TRAN) at time zero, were it not for exact cancellation
C by a one-amp Type-11 current source that is connected to node "TRAN".
C It will be noted that the current (GEN, TRAN) continues smoothly out
C of the phasor solution. The current of the NL reactor (TRAN, TERRA)
C shows a one-amp jump beginning on the first step, due to the internal,
C TACS-controlled source (remember, the inductor current is continuous).
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
.000200 .020
1 1 1 1 1 -1 0 2
5 5 10 10
TACS HYBRID { We use TACS only to produce signal that controls shunt current
99SHUNT = 1.0 { Unit current (dc source) is inserted in parallel with ZnO
33SHUNT { Output only this one TACS variable that controls NL shunt current
77SHUNT 1.0 { Initial condition required for smooth electrical step 1
BLANK card ending all TACS data
GEN TRAN 2.0 1
93TRAN 1.0 1.0 1
0.0 0.0
0.9 0.9
2.0 1.1
10. 1.2
9999
TACS CONTROL SHUNT { Use 2nd of 3 names: for TACS-controled shunt current
BLANK card ending branch cards
BLANK card ending (here non-existent) switch cards
14GEN 377. 60. -1.
11TRAN -1 -1.0 { 1-amp external dc current cancels internal NL source
BLANK card ending source cards
C Total network loss P-loss by summing injections = 1.000018973436E+00
C
C GEN 377. 377. .00530514044263 1.000009486673 1.0000189734357
C 0.0 0.0 -.9999954144499 -89.6960393 188.49913562381
C
C ---- Initial flux of coil "TRAN " to " " = 5.30514044E-03
C
C Step Time TRAN GEN TRAN GEN TACS
C TERRA TRAN SHUNT
C 0 0.0 376.9893897 377. .0053051404 .0053051404 1.0
C 1 .2E-3 375.7677437 375.9289054 -.919419146 .0805808538 1.0
1
C 100 .02 114.1251632 116.4994069 .1871218422 1.187121842 1.0
C Variable max : 376.9893897 377. .4424771753 1.442477175 1.0
C Times of max : 0.0 0.0 .0042 .0042 0.0
C Variable min : -376.817325 -376.880939 -2.44858983 -1.44858983 1.0
C Times of min : .0084 .0084 .0124 .0124 0.0
PRINTER PLOT
C -0.653 -0.353 -0.054 0.245 0.545 0.844 1.143 1.442
C ------+---------+---------+---------+---------+---------+---------+---------+
C BBBBBBBBBBBBBBBBBBBBBBBBBBBB*AA
C BBBBB | AAAAAA
C BBBBBBB | AAAAAAAA
C BBBBB | AAAAA
C BBBB | AAAA
C BBBBBBBBBBBBBB AAAAAAAAAAAAAA
C | BBB AAA
C | BBBBBBBBBB AAAAAAAAA
C BBBBBB AAAAAAA
C BBBBB | AAAA
C BBBBB | AAAAA
C BBBBBBBB | AAAAAAAA
C BB |AAAAAA
194 4. 0.0 20. GEN TRAN TRAN { Axis limits: (-2.449, 1.442)
144 4. 0.0 20. TRAN { Axis limits: (-3.768, 3.770)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 3rd of 4 subcases is related to 2nd. There is the same nonlinear L(i)
C inductor (Type-93, true-nonlinear element). This example illustrates
C optional TACS control of the (i, PSI) characteristic by multiplicative
C scaling. The TACS multiplier of "i" is declared by the third data
C field (columns 27-32) of the "TACS CONTROL" card. For this example,
C the TACS variable is FACTOR, which varies smoothly (via one half a
C cosine cycle) from 2.0 to 1.0, thereby reducing the current that is
C drawn by approximately a factor of two. The series resistor is small,
C to it does not drop much voltage. As a result, the voltage, and hence
C the peak flux of the inductor, are nearly constant. Only the current
C increases by about a factor of two over the 10 msec "ramp" time. The
C presence of resistance merely reduces this ideal result: the factor is
C less than two --- maybe 1.5 or so.
PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns)
.000200 .030
1 1 1 1 1 -1 0 2
5 5 49 1 55 5
TACS HYBRID { We use TACS only to produce signal that controls shunt current
99FACTOR = 1.5 + COSINE + STEP { Scaling factor changes from 3 to 1 over 0-10ms
14COSINE 0.5 50. { Smooth transition over 1st 10 msec } .010
11STEP -0.5 50. { Hold COSINE value after 10 ms} .010
33FACTORSTEP COSINE { Output only the 3 user-defined TACS variables
77FACTOR 2.0 { Initial condition required for smooth electrical step 1
BLANK card ending all TACS data
GEN TRAN 20. 1
93TRAN 2.0 1.0 1
0.0 0.0
0.9 0.9
2.0 1.1 { Execution does reach this 2nd seg of (i, PSI)
9999
TACS CONTROL FACTOR { Use 3rd of 3 names: TACS-scaled curve (x,y)
BLANK card ending branch cards
BLANK card ending (here non-existent) switch cards
14GEN 350. 50. -1.0
BLANK card ending source cards
C Total network loss P-loss by summing injections = 4.885536661314E+01
C ---- Initial flux of coil "TRAN " to " " = 1.39586762E-01
C GEN 350. 350. .27917352350367 2.2103250125975 48.855366613143
C 0.0 0.0 -2.19262372629 -82.7439171 383.7091521007
C
C Step Time TRAN GEN TRAN GEN TACS
C TERRA TRAN FACTOR
C
C 0 0.0 344.4165295 350. .2791735235 .2791735235 2.0
C 1 .2E-3 340.9842812 349.3093549 .4162536857 .4162536857 1.999013364
C 2 .4E-3 336.2117309 347.2401455 .5514207281 .5514207281 1.996057351
C
C 49 .0098 -347.838306 -349.309355 -.073552423 -.073552423 1.000986636
C 50 .01 -347.141951 -350. -.14290245 -.14290245 1.0
C 51 .0102 -345.0697 -349.309355 -.211982761 -.211982761 1.0
1
C 150 .03 -347.789928 -350. -.110503595 -.110503595 1.0
C Variable maxima : 349.3458152 350. 3.110373844 3.110373844 2.0
C Times of maxima : .0198 0.0 .0042 .0042 0.0
C Variable minima : -348.491847 -350. -2.24512833 -2.24512833 1.0
C Times of minima : .0298 .01 .0146 .0146 .01
PRINTER PLOT
194 3. 0.0 30. GEN TRAN TACS FACTOR { Axis limits: (-2.245, 3.110)
C CALCOMP PLOT { Repeat this same plot, only in high-resolution mode
C 194 3. 0.0 30. GEN TRAN TACS FACTOR
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 4th of 4 subcases is related to 3rd. That nonlinear inductor L(i)
C (Type-93, true-nonlinear element) has been made linear by reduction of
C the flux-current characteristic to a single point. The result is a
C TACS-controlled linear inductor L(t). Except that here dynamics are
C involved (because L stores energy, and has history), this is the
C inductive equivalent of commonly-used TACS-controlled resistance R(t).
C Between time 0 and 100 msec, inductance is to be varied sinusoidally
C between 125 mH and 500 mH. Disconnected branches BEGIN and END
C correspond to these two limiting values, and they define assymptotes
C of the transition between the two. The graph is smooth, and it shows
C a smooth transition between the two steady-state assymptotes.
PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns)
.000400 .120
1 1 1 1 1 -1
5 5
TACS HYBRID { We use TACS only to produce signal that controls shunt current
99FACTOR = 2.5 + COSINE + STEP { Scale factor changes from 4 to 1 over 0-100 ms
14COSINE 1.5 5.0 { Smooth transition over 1st 100 msec } .100
11STEP -1.5 { Hold COSINE valu after 100 ms} .100
33FACTORSTEP COSINE { Output only the 3 user-defined TACS variables
77FACTOR 4.0 { Initial condition to make output signal smooth at t = 0
77COSINE 1.5 { Initial condition to make output signal smooth at t = 0
BLANK card ending all TACS data
GEN TRAN 35. { R of series R-L branch to ground } 1
93TRAN 8.0 1.0 { i-steady and PSI-steady for phasors
0.0 0.0 { Symmetric curve always begins at origin 0,0
2.0 1.0 { Single (i, PSI) point defines L b4 FACTOR
9999 { Terminate the characteristic after just the one point
TACS CONTROL FACTOR { Use 3rd of 3 names: TACS-scaled curve (x,y)
C That is all we need for the TACS-controlled L of series R-L circuit. But
C for a nice plot, we want assymptotic curves. The first in in the steady
C state, corresponding to i-steady and PSI-steady of the Type-93 branch
C card. This is L = PSI / i = 1.0 / 8.0 = 125 mH. Note that this matches
C the initial signal. Initial FACTOR = 4.0 and the characteristic gives
C L = PSI / i = 1.0 / 2.0 = 500 mH. Dividing by FACTOR gives 125 mH as the
C initial inductance within the dT loop. So, show phasor flow:
GEN BEGIN 35. 125. { Branch flow is phasor solution
C That was for the start. Next, we want the termination, which corresponds
C to FACTOR = 1.0 so no scaling. I.e., L = 500 mH (the characteristic):
GEN END 35. 500. { Branch flo is assymptote for infinite t
BLANK card ending branch cards
C Switches are used only to measure the currents through the preceding two
C branches. These are the two assymptotic signals of following plot:
BEGIN MEASURING 1
END MEASURING 1
BLANK card ending switch cards
C --------------+------------------------------
C From bus name | Names of all adjacent busses.
C --------------+------------------------------
C GEN |TRAN *BEGIN *END *
C TRAN |TERRA *GEN *
C BEGIN |TERRA *GEN *
C END |TERRA *GEN *
C TERRA |TRAN *BEGIN *END *
C --------------+------------------------------
14GEN 420. 50. -1.0
C TRAN 234.06699294777 313.5412844237 5.3123716300638 7.9842632447061 -.677791157E-13 -.677791157E-13
C 208.61634607684 41.7095865 -5.960467030767 -48.2904135 1251.698076461 1251.6980765
C
C TERRA 0.0 0.0 -5.312371630064 7.9842632447061 0.0
C 0.0 0.0 5.9604670307669 131.7095865 0.0
C Total network loss P-loss by summing injections = 2.350389825256E+03
C Output for steady-state phasor switch currents.
C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive
C BEGIN 5.31237163E+00 -5.96046703E+00 7.98426324E+00 -48.2904 0.00000000E+00 0.00000000E+00
C END 5.67589241E-01 -2.54733456E+00 2.60980285E+00 -77.4387 0.00000000E+00 0.00000000E+00
BLANK card ending source cards
GEN TRAN { Output node voltages for the controlled branch L(t) only
C First 2 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 Next 3 output variables belong to TACS (with "TACS" an internally-added upper name of pair).
C Step Time GEN TRAN BEGIN END GEN TACS TACS TACS
C TERRA TERRA TRAN FACTOR STEP COSINE
C *** Phasor I(0) = 5.3123716E+00 Switch "BEGIN " to " " closed in the steady-state.
C *** Phasor I(0) = 5.6758924E-01 Switch "END " to " " closed in the steady-state.
C 0 0.0 420. 234.06699 5.3123716 .56758924 5.3123716 4.0 0.0 1.5
C 1 .4E-3 416.68817 206.10551 6.0166476 .88197067 6.0166476 3.9998816 0.0 1.4998816
C 2 .8E-3 406.80493 174.8932 6.6262352 1.1824658 6.6260494 3.9995263 0.0 1.4995263
C 3 .0012 390.50612 140.93094 7.1314997 1.464335 7.1307196 3.9989342 0.0 1.4989342
BLANK card ending node voltage outputs
C 290 .116 129.78714 208.41029 -4.028677 -2.244816 -2.246376 1.0 -1.5 0.0
C 295 .118 339.78714 376.14579 .78854743 -1.037462 -1.038819 1.0 -1.5 0.0
C 300 .12 420. 400.22543 5.3045736 .56616677 .5649877 1.0 -1.5 0.0
C Variable maxima : 420. 409.93775 7.9767416 2.6063662 7.9698488 4.0 0.0 1.5
C Times of maxima : 0.0 .0992 .0028 .0044 .0028 0.0 0.0 0.0
C Variable minima : -420. -409.6697 -7.973184 -2.60532 -7.845698 1.0 -1.5 -1.499882
C Times of minima : .01 .1092 .1128 .1144 .0128 0.1 0.1 .0996
CALCOMP PLOT
C Orlando Hevia in Santa Fe, Argentina, supplied the following one-line
C correction in E-mail dated 10 January 2007. He learned of the bad ATP data
C during experimentation with the new GNU F95 compiler, which choked on this:
C NOTE A MISTAKE IN DATA 120. IS ONE COLUMN TOO TO THE RIGHT
C DATA FM1180 / '( I2, 2I1, E3.0, 2E4.0, E5.0, E4.0, 4A6, 2A16 )' /
C -=---====----=====----
C 19412. 0.0 120. BRANCH Vary linear L Inductor current
19412. 0.0120. BRANCH Vary linear L Inductor current
BEGIN END GEN TRAN TACS FACTOR
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
|
