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
|
BEGIN NEW DATA CASE
C BENCHMARK DC-19
C Test of the control of TACS sources by node voltages and switch currents
C of the electric network. The electric network has no dynamics, with both
C the node voltage and switch current equal to one half, always. TACS
C supplemental variable GENT adds these to give unity. The output signal
C is then passed through a function block 1/(1+S) to give ELEMT, and
C through a S/(1+S) block to give GT/ELT. There also is a disconnected
C validation of the Type-17 electric network source, which provides for
C modulation of the following source function by a TACS variable. For use
C here, the TACS variable is the constant TWO = 2.0 --- easy to verify.
.05 2.0
1 1 1 1 1
TACS HYBRID
1ELEMT +GENT
1.0
1.0 1.0
1GT/ELT +GENT
1.0
1.0 1.0
90BUS2
91BUS1
11TWO 2.0 { Constant for modulation of electric network COSINE
99GENT = BUS1 + BUS2
33ELEMT GENT GT/ELTBUS1 BUS2
BLANK card ending all TACS data
BUS3 BUS2 1.0
BUS2 BUS1 1.0
COSINE 1.0 { Dummy branch connects source of interest
SURGE 1.0 { Dummy branch connects Type-15 surge function
TAKUSG 1.0 { Dummy branch for Taku Noda's Type-15 surge
BLANK card ending electric network branches
BUS1 MEASURING 1
BLANK card ending switches
11BUS3 1.0
17TWO { The constant 2.0 of TACS multiplies following 1/2 to give unity:
14COSINE 0.5 0.5 { Cosine starts at peak 0.5, and also ends there
15SURGE 2.0 -1.0 -3.0 { Surge function is a smooth pulse > 0
15TAKUSG 2.0 -1.0 -3.0 -6666. { Same except 2.0 = peak valu
BLANK card ending electric network source cards.
C Step Time BUS1 BUS2 BUS3 COSINE SURGE TAKUSG BUS1 TACS TACS TACS
C TERRA ELEMT GENT GT/ELT
C *** Switch "BUS1 " to " " closed before 0.00000000E+00 sec.
C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C 1 .05 0.0 0.5 1.0 0.0 .181042896 .470363242 0.5 .024390244 1.0 .975609756
1
C 39 1.95 0.0 0.5 1.0 .987688341 .278788345 .724313367 0.5 .854137166 1.0 .145862834
C 40 2.0 0.0 0.5 1.0 1.0 .265713062 .690342786 0.5 .861252427 1.0 .138747573
C maxima : 0.0 0.5 1.0 1.0 .769799804 1.99999856 0.5 .861252427 1.0 .975609756
C maxima : 0.0 .05 .05 2.0 .55 .55 .05 2.0 .05 .05
PRINTER PLOT
143 .4 0.0 2.0 SURGE COSINE { Axis limits: (-1.000, 1.000)
193 .4 0.0 2.0 TACS GT/ELTTACS ELEMT { Axis limits: (0.000, 9.756)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 2nd of 5 subcases of DC-19 tests important supplemental variable functions
C Unlike 1st, there is no electric network (in theory). In practice, ATP
C will automatically supply a dummy electric network. This is since around
C 1993 when TACS STAND ALONE was internally converted to TACS HYBRID.
.100 6.0
1 1 0 1 0 0 0 0
TACS STAND ALONE
X2 -X3 +TQ 6.6667 -100. 0.2
1X3 +X2 0.0 1.05
1.0
0.0 1.
11TQ 1.0 0.0
98SUPVAR = 0.4 * ( ABS(X3) + 1.E-6 ) ** 3.0 * 2.0
26T { Type-26 is given value equal to ATP variable in columns 3-8
33X2 X3 T SUPVAR
C Step Time TACS TACS TACS TACS
C X2 X3 T SUPVAR
C 0 0.0 0.0 0.0 0.0 0.0
C 1 0.1 0.2 .01 0.1 .80024E-6
C 2 0.2 0.2 .03 0.2 .216022E-4
BLANK card ending all TACS data cards
C 58 5.8 .203206E-4 .999997968 5.8 .799997523
C 59 5.9 .677346E-5 .999999323 5.9 .800000774
C 60 6.0 .22578E-5 .999999774 6.0 .800001858
PRINTER PLOT
143 1. 0.0 6.0 X2 SUPVAR { Axis limits: (0.000, 8.000)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 3rd of 5 subcases consists only of Type-15 lightning functions, including
C both old (the 2-exponential surge function) and new (Bernd Stein of FGH)
C models. Bernd's FGH improvement is documented in a letter dated 30.7.86
C Smoothing of plot data points after 50 ups and downs is also illustrated.
1.0E-8 3.E-6
1 1 1 1 1 -1 1
5 5 20 20
LIGHTN 1.0
ALGHTN 1.0
LIGHT1 1.0
ALGHT1 1.0
ALGHT2 1.0
ALGHT3 1.0
STAN1 1.0 { Resistor for Standler surge added 29 Sept 00
STAN2 1.0 { Resistor for 2nd Standler surge added 27 Oct 00
TWOE1 1.0 { Resistor for TWO EXP surge added 19 Jan 01
TWOE2 1.0 { Resistor for 2nd TWO EXP surge added 19 Jan 01
HEID1 1.0 { Resistor for Heidler surge added 19 Jan 01
HEID2 1.0 { Resistor for 2nd Heidler surge added 19 Jan 01
C Orlando Hevia mentions legality and danger of negative L of series R-L-C.
C Add a branch to illustrate as mentioned in newsletter dated January, 2001
HEVIA 10.E3 -3.0 { Series R-L with negative L } 1
BLANK card ending electric network branches
BLANK card ending switches
15ALGHTN-1 1.0 2.E-6 50.E-6 5.0
15LIGHTN-1 1.0 -20000. -500000.
15ALGHT2-1 1.0 2.E-6 50.E-6 10.0
15ALGHT1-1 1.0 2.E-6 50.E-6 5.0
15LIGHT1-1 1.0 -20000. -500000.
15ALGHT3-1 1.0 2.E-6 50.E-6 10.0
C 345678901234567890123456789012345678901234567890
C <--- Off-line Standler 3F10.0 ->
15STAN1 -1 0.6 20.E-6 0.04 Standler { "Standler" is same as -7777
C C A B { Identify of 3 Standler parameters
C <--- In-line Standler 4F8.0 ---><-Request word-> <T-start >< T-stop >
C T1 T2 Perc peak
15STAN2 -1 1.0E-6 20.0E-6 0.0 1.000 Standler in-line 0.0 1.0
C
C <--- Off-line Two exp 3F10.0 ->
15TWOE1 -1 1.20 -36000.0 -5.0E+6 { it is an old two exp source
C C A B { Identify of 3 Two exp parameters
C <--- In-line Two exp 4F8.0 ---><-Request word-> <T-start >< T-stop >
C T1 T2 Perc peak
15TWOE2 -1 1.0E-6 20.0E-6 0.0 1.0 Two exp in-line 0.0 1.0
C
C <--- Off-line Heidler 3F10.0 -->
15HEID1 -1 0.6 1.0E-06 3.0E-5 5.0 { it is an old Heidler source
C C A B N { Identify of 4 Heidler parameters
C <-Heidler in line 4F8.0 ---><-key----------><><T-start >< T-stop >
C T1 T2 Perc peak N
15HEID2 -1 1.0E-6 20.0E-6 0.0 1.000 Heidler in-line 5 0.0 1.0
C
14HEVIA 1.0 60. { Voltage source is shorted at 1 us } 1.E-6
BLANK card ending electric network source cards.
STAN1 STAN2 ALGHT2ALGHT1LIGHT1ALGHTNLIGHTNALGHT3HEID1 HEID2 TWOE1 TWOE2
C First 12 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 STAN1 STAN2 ALGHT2 ALGHT1 LIGHT1 ALGHTN LIGHTN ALGHT3 HEID1 HEID2
C
C
C TWOE1 TWOE2 HEVIA
C TERRA
C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C 0.0 0.0 0.0
C 1 .1E-7 .442479244 .855069065 .25405E-20 .812971E-4 957.508161 .406485E-9 .004787541 .50811E-15 .930195E-8 .140651E-7
C .058092768 .04993186 -.16949E-5
C 2 .2E-7 .454691626 .88074615 .2601E-17 .002438392 -4.9670744 .130049E-7 .009550246 .51918E-12 .297563E-6 .449919E-6
C .113331409 .097441938 -.51422E-5
C 3 .3E-7 .461895211 .895961187 .14996E-15 .01470791 952.565899 .987364E-7 .01428824 .28952E-10 .225886E-5 .34153E-5
C .165855128 .142646809 -.87063E-5
C 4 .4E-7 .467007506 .906801471 .26623E-14 .048743046 -9.8846491 .415991E-6 .019001647 .473526E-9 .951554E-5 .143866E-4
C .21579634 .185657436 -.12391E-4
C 5 .5E-7 .470959015 .91521116 .2479E-13 .121908619 947.672889 .126925E-5 .023690588 .395202E-8 .290285E-4 .43887E-4
C .263281003 .226579442 -.16201E-4
BLANK card ending output requests
C 260 .26E-5 .48556872 .972472279 .991750958 -15.735499 -362.71992 .993881314 .676797074 -18.724315 .572623725 .949657569
C 1.09277363 .94980214 -.55047746
C 280 .28E-5 .482164403 .967184298 .987942471 -17.163302 -375.61152 .990578487 .698942172 -19.269102 .568828878 .942698718
C 1.08493561 .942824605 -1.0722499
C 300 .3E-5 .478686003 .961711993 .984066679 -17.995991 -387.2694 .987055255 .718634373 -19.452928 .565055727 .935784527
C 1.07715275 .935897248 -2.0885866
C Variable maxima : .506829483 1.00000002 1.0 1757.33053 957.508161 1.0 .718634373 2231.07839 0.6 1.0
C 1.14948506 .999999945 0.0
C Times of maxima : .8E-6 .1E-5 .2E-5 .7E-6 .1E-7 .2E-5 .3E-5 .113E-5 .1E-5 .1E-5
C .99E-6 .1E-5 0.0
C Variable minima : 0.0 0.0 0.0 -17.995991 -387.2694 0.0 0.0 -19.452928 0.0 0.0
C 0.0 0.0 -2.0885866
C Times of minima : 0.0 0.0 0.0 .3E-5 .3E-5 0.0 0.0 .3E-5 0.0 0.0
C 0.0 0.0 .3E-5
PRINTER PLOT
145 .4 0.0 3.0 LIGHTNALGHTNALGHT2 { Plot limits: (0.000, 1.000)
145 .4 0.0 3.0 LIGHT1ALGHT1ALGHT3 { Plot limits: (-0.110, 2.231)
C 145 .4 0.0 3.0 STAN1 STAN2
C 145 .4 0.0 3.0 HEID1 HEID2
C 145 .4 0.0 3.0 TWOE1 TWOE2
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 4th of 5 subcases was added 17 September 1999 to illustrate delayed
C opening of a Type-13 TACS-controlled switch. There are two parallel,
C identical circuits on the electrical side. A sinusoidal source drives
C current through a series R-L branch. The old logic is illustrated by
C the QUICK alternative whereas the new is the DELAY alternative.
C The new optional DO NOT OPEN UNTIL CURRENT 0 is illustrated by the
C latter of these two: opening will be delayed until current passes
C through zero. Both circuits oscillate following opening because both
C have a hanging inductor. But QUICK opens first, and it oscillates
C more. The graph is distinctive.
PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns)
.0004 .010
1 1
TACS HYBRID
14CLAMP 100. 50. 30. 0. -1.
33CLAMP
BLANK card ending TACS
GEN QUICK 10. 4. { Series R-L will oscillate when opened
GEN DELAY 10. 4.
BLANK card ending branches
13DELAY Do not open until current 0 CLOSED CLAMP 11
$DISABLE
C The preceding involves a high-level, English-language declaration. It
C is equivalent to the following low-level, numeric declation.
13RESIS -77333. CLOSED CLAMP 11
C I.e., the data field of otherwise-unused columns 25-34 can be used to
C hold the special numeric value. About the English, this must begin
C to the right of column 24. If 55-60 are used as here for CLOSED, the
C beginning of "do not open until current 0" can be delayed at most 3
C columns (note 3 blanks between ending 0 and start of CLOSED). Case is
C arbitrary (note 1st letter has been made upper case on switch card).
$ENABLE
13QUICK CLOSED CLAMP 11
BLANK card ending switches
14GEN 100. 50. 0.0 0. -1.
BLANK card ending sources
GEN QUICK DELAY
BLANK card ending voltage printout
CALCOMP PLOT { Needed to cancel PRINTER PLOT of preceding subcase
SUPERIMPOSE 2 { Superimpose plots of the following 2 cards
144 1. 0.0 10. GEN DELAY QUICK { First, three node voltages
194 1. 0.0 10.-200.200.TACS CLAMP { 2nd of 2 in overlay is I-branch
BLANK card ending plot
BEGIN NEW DATA CASE
C 5th of 5 subcases was added 28 July 2002 to illustrate 3 new source
C types: CIGRE, CESI, and USRFUN. These are all special sub-types
C within the Type-15 source. The supporting code comes from Orlando P.
C Hevia of UTN in Santa Fe, Argentina. For an alternate and different
C illustration of the USRFUN alternative, see DC-5. For comments
C about required precision, see note following this subcase.
.05E-6 10.E-6 { Hevia used dT = 1.0E-9, but such detail is not needed
1 1 1 1 -1
5 5 20 20
CESI1 1.0
CESI2 1.0
CIGRE1 1.0
CIGRE4 1.0
USER2 1.0
USER8 1.0 { Add another USRFUN test on 10 December 2002
BLANK card ending branch cards
BLANK card ending switch cards (none for this data)
C CORRENTE DI CRESTA -0.001 [kA]
C TEMPO DI CRESTA 1 [us]
C TEMPO EMIVALORE 10 [us]
C GENERATORI EQUIVALENTI FRONTE
15CESI1 -1 1.5 1.0E-6 10.0-6 CESI
15CESI2 -1 -1.2 1.0E-6 10.0-6 CESI
15CIGRE1 1.0 1.0E-6 10.0E-6 CIGRE 20.0E+5 0.5E-6 9.5E-6
15CIGRE4 1.0 1.0E-6 10.0E-6 CIGRE 50.0E+5 0.5E-6 9.5E-6
15USER2 2.0 USRFUN
15USER8 8.0 USRFUN
BLANK card ending electric network sources
1 { Request for all possible node voltage outputs (here, just 4)
CALCOMP PLOT
145 .5 0.0 5.0 CESI1 CESI2 USER2 CIGRE1
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
EOF
Warning about the preceding. The CIGRE1 output depends greatly on arithmetic
precision. Even 64 bits seems inadequate since turning the Salford debugger
on and off (minor update using SUBR5) results in the following differences.
Note that the 3rd row is what will be published in a future newsletter. The
CIGRE1 value on step 40 rises from .661 to .668 when the debugger is
turned on.
FC version 2.0 - Copyright (c) 1990 Mike Albert Wed Aug 21 14:08:13 2002
Options are: -c1 -ds5 -t8
Compare files:
dc19.lis 71885 8-21-102 2:07p
dc19.sal 71885 8-21-102 1:48p
Changed lines 859-868
> 15 .75E-6 1.000594 -.800475 .0648095 .0662431 .0783034
> 20 .1E-5 1.477403 -1.18192 .1299899 .1324862 .1044063
> 40 .2E-5 1.338601 -1.07088 .668485 .4280678 .2093254
> 60 .3E-5 1.23813 -.990504 .9529333 .9505637 .3267037
> 80 .4E-5 1.167656 -.934125 .8811803 .8800206 .5541328
> 100 .5E-5 1.104946 -.883957 .8148301 .8147126 .9983771
> 120 .6E-5 1.044245 -.835396 .7534758 .7542513 .9931261
> 140 .7E-5 .9840634 -.787251 .6967413 .6982769 .9835613
> 160 .8E-5 .9240164 -.739213 .6442788 .6464564 .9740871
> 180 .9E-5 .8640042 -.691203 .5957665 .5984817 .9647042
To
> 15 .75E-6 1.000594 -.800475 .0644145 .0662431 .0783034
> 20 .1E-5 1.477403 -1.18192 .1291901 .1324862 .1044063
> 40 .2E-5 1.338601 -1.07088 .661331 .4280678 .2093254
> 60 .3E-5 1.23813 -.990504 .9533223 .9505637 .3267037
> 80 .4E-5 1.167656 -.934125 .88154 .8800206 .5541328
> 100 .5E-5 1.104946 -.883957 .8151627 .8147126 .9983771
> 120 .6E-5 1.044245 -.835396 .7537834 .7542513 .9931261
> 140 .7E-5 .9840634 -.787251 .6970257 .6982769 .9835613
> 160 .8E-5 .9240164 -.739213 .6445418 .6464564 .9740871
> 180 .9E-5 .8640042 -.691203 .5960097 .5984817 .9647042
Changed line 872
>Variable maxima : 1.497392 0.0 .999673 .9989941 .9999945
To
>Variable maxima : 1.497392 0.0 .9997968 .9989941 .9999945
Comparison complete
21 August 2002, remove from before the 5th, which was not being executed:
BLANK ----------------- Protect following 4th data subcase from execution
C Note: Following is copy of 3rd subcase between October, 1990, and
C July, 1991. It was wrong during 9 months. The solution
C prior to October, 1990, was kept below here, and discovered
C to be correct once again on July 22, 1991. So, it would
C seem that Guido's "correction" of Oct, 1990, was wrong, and
C during the spring or summer of 1991 he merely corrected the
C correction. Well, as history, we keep the wrong data:
C 3rd of 3 subcases consists only of Type-15 lightning functions, including
C both old (the 2-exponential surge function) and new (Bernd Stein of FGH)
C models. Bernd's FGH improvement is documented in a letter dated 30.7.86
C Smoothing of plot data points after 50 ups and downs is also illustrated.
C Comment cards altered 23 Oct 1990. Original contents preserved below
C as 4th data subcase until Bernd Stein approves of changes.
1.E-8 300.E-8
1 1 1 1 1 -1
5 5 20 20
LIGHTN 1.0
ALGHTN 1.0
LIGHT1 1.0
ALGHT1 1.0
ALGHT2 1.0
ALGHT3 1.0
BLANK card ending electric network branches
BLANK card ending switches
15ALGHTN-1 1.0 2.E-6 50.E-6 5.0
15LIGHTN-1 1.0 -20000. -500000.
15ALGHT2-1 1.0 2.E-6 50.E-6 10.0
15ALGHT1-1 1.0 2.E-6 50.E-6 5.0
15LIGHT1-1 1.0 -20000. -500000.
15ALGHT3-1 1.0 2.E-6 50.E-6 10.0
BLANK card ending electric network source cards.
C Step Time ALGHT2 ALGHT1 LIGHT1 ALGHTN LIGHTN ALGHT3
C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C 1 .1E-7 .26503E-20 .853443E-4 957.508161 .426722E-9 .004787541 .53005E-15
C 2 .2E-7 .27133E-17 .002559784 -4.9670744 .136524E-7 .009550246 .54161E-12
C 3 .3E-7 .15643E-15 .015440124 952.565899 .103652E-6 .01428824 .30202E-10
C 4 .4E-7 .27773E-14 .051169655 -9.8846491 .436701E-6 .019001647 .493979E-9
1
C 300 .3E-5 1.02657121 -18.891898 -387.2694 1.03619451 .718634373 -20.293153
C maxima : 1.04319273 1844.81694 957.508161 1.0497837 .718634373 2327.44477
C Times of max : .2E-5 .7E-6 .1E-7 .2E-5 .3E-5 .113E-5
C minima : 0.0 -18.891898 -387.2694 0.0 0.0 -20.293153
C Times of min : 0.0 .3E-5 .3E-5 0.0 0.0 .3E-5
PRINTER PLOT
145 .4 0.0 3.0 LIGHTNALGHTNALGHT2 { Plot limits: (0.000, 1.050)
145 .4 0.0 3.0 LIGHT1ALGHT1ALGHT3 { Plot limits: (-0.110, 2.327)
BLANK card ending plot cards
BEGIN NEW DATA CASE
|