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
|
BEGIN NEW DATA CASE
C BENCHMARK DC-62 (DC-66 of "M39." vintage)
C Same four-phase South African line simulation as in DC-61, only here a
C high-order Hauer model ("HAUER SETUP", no longer in program) is used.
C 25 November 2000, WSM converts [Tv] and [Ti] from original use of
C free-format for complex elements to wide format for real parts only.
C This involves use of NO IMAGINARY PART as described in a newsletter
C (April 2001 or later).
PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 9 output columns)
PEAK VOLTAGE MONITOR
POWER FREQUENCY, 50.,
.000010 .020 50. 50.
1 5 1 1 1 -1
5 5 20 20 100 100 1000 200
SA-3 BUSA-3 .1260.75
SA-4 BUSA-4 .1260.75
BUSA-1 2. .03
BUSA-2 2. .03
BUSA-3 FIL1 .1 47.12
FIL1 .1 90.76
FIL1 3. 74.24
FIL1 6.166.94 94.25
FIL1 12.314.16 11.06
FIL1 12.314.16 5.87
FIL2 .01 1021.
FIL2 .1 14.14
BUSA-4 FIL3 .1 47.12
FIL3 FIL4 .1 90.76
FIL3 FIL4 3. 74.24
FIL3 FIL4 6.166.94 94.25
FIL3 FIL4 12.314.16 11.06
FIL3 FIL4 12.314.16 5.87
FIL4 .01 1021.
FIL4 .01 14.14
C &&&&&&&&&&&&&&&&&&&&&& BEGIN UNIT-13 CARDS PUNCHED BY DC-65 .
C #DATA FOR LINE SECTION 1: LENGTH = 0.2330096E+03
C METRIC
C 1 .2 .273 4 1.7653 0. 39. 20.
C 2 .2 .273 4 1.7653 30. 39. 20.
C 3 .392 .0585 4 3.1750 2.3 27.5 8.5 50.8 45.
C 4 .392 .0585 4 3.1750 32.3 27.5 8.5 50.8 45.
C
C 100. 5000. 1 375.
C 100. 50. 1 375.
C 100. 5. 1 375. 9 21
C
-1BUSA-1BUSB-1 0.74300E-02 0.13695E-02 1 1 14 6 4
0.24524677E+02 0.13280541E+03 0.00000000E+00 0.29592165E-02 0.50000000E+02 0
0., 0.301017131E+03, 0.240474368E-01, 0., 0.911921587E+04, 0.866135515E+01
0., 0.142378188E+03, 0.467439158E-01, 0., 0.597875358E+05, 0.157341296E+05
0., 0.129882539E+04, 0.158087056E+00, 0., 0.602061446E+05, 0.350951915E+05
0., 0.108960413E+05,-0.108480154E+02, 0., 0.600607627E+05,-0.503745837E+05
1., 0.289467030E+05, 0.211443554E+00, -1.,-0.537514515E+04,-0.421452862E+01
0., 0.178074955E+04, 0.254168820E+00, 0., 0.315149545E+04, 0.732001459E+00
0., 0.239325862E+04,-0.862687875E-01, 0., 0.618967874E+05,-0.453102444E+03
0., 0.619168960E+01,-0.420104801E-02, 0., 0.100000000E+01, 0.000000000E+00
0., 0.157926984E+03,-0.556405973E-03, 0., 0.179349724E+04,-0.616894886E-03
0., 0.154682388E+05,-0.740819400E-03, 0., 0.298676381E+06,-0.131478773E-02
-1BUSA-2BUSB-2 0.77523E-02 0.12698E-02 2 2 14 6 4
0.51203735E+01 0.54502236E+02 0.00000000E+00 0.31183848E-02 0.50000000E+02 0
0., 0.313224674E+02,-0.663044870E-02, 0., 0.210626498E+06, 0.432723088E+01
0., 0.182100421E+02, 0.260951895E-01, 0., 0.849632630E+06, 0.182628050E-02
0., 0.227829005E+05, 0.478630198E+00, 0., 0.102901020E+07,-0.212641132E-02
0., 0.124533203E+06, 0.211022159E+02, 0., 0.157231418E+07, 0.281738783E-03
0., 0.142701493E+06,-0.267953188E+02, 0., 0.181740797E+07,-0.719862890E-04
1., 0.398474670E+05, 0.418428806E+00, -1.,-0.214066736E+05, 0.125874270E+00
1., 0.151593680E+06, 0.515504958E+00, -1.,-0.147087464E+06, 0.267749108E+00
0., 0.365121200E+01,-0.350144906E-02, 0., 0.100000000E+01, 0.000000000E+00
0., 0.249321282E+04,-0.431759968E-04, 0., 0.182369337E+06,-0.242606205E-04
0., 0.239163482E+02,-0.406028697E-02, 0., 0.207918724E+06,-0.123123552E-03
-1BUSA-3BUSB-3 0.32619E-02 0.12582E-02 3 3 14 6 4
0.53536833E+02 0.12389964E+03 0.00000000E+00 0.12908552E-02 0.50000000E+02 0
0., 0.134598777E+06,-0.150159946E-01, 0., 0.140301928E+07,-0.753366163E-01
0., 0.284966044E+02, 0.819451304E-01, 0., 0.666744070E+07, 0.157059296E+04
0., 0.272105388E+05, 0.213806667E+00, 0., 0.782363486E+07, 0.971629566E+04
0., 0.272559469E+06,-0.223283773E+00, 0., 0.799352013E+07,-0.800863365E+06
0., 0.361659615E+06, 0.100462472E+01, 0., 0.799483818E+07, 0.792080675E+06
0., 0.491172730E+05,-0.165633449E+00, 0., 0.599908733E+07,-0.177051614E+03
0., 0.792280597E+05, 0.366175685E+00, 0., 0.698913494E+07,-0.232733391E+04
0., 0.115763918E+02,-0.117805658E-02, 0., 0.100000000E+01, 0.000000000E+00
0., 0.197441131E+04,-0.373045199E-04, 0., 0.181295123E+06, 0.445567784E-04
0., 0.976839246E+02,-0.202948282E-02, 0., 0.208740876E+06,-0.616217575E-04
-1BUSA-4BUSB-4 0.30402E-02 0.12587E-02 4 4 6 6 4
0.64988120E+02 0.13246385E+03 0.00000000E+00 0.12025589E-02 0.50000000E+02 0
0., 0.169333103E+03, 0.332803226E-01, 0., 0.566559087E+06, 0.730448580E+00
0., 0.680693506E+02, 0.708523592E-01, 0., 0.490942013E+07,-0.111902264E+01
0., 0.602799468E+05, 0.326699789E+00, 0., 0.528347745E+07, 0.957741589E+00
0., 0.646378180E+01,-0.792079833E-03, 0., 0.100000000E+01, 0.000000000E+00
0., 0.142573596E+03,-0.117927457E-02, 0., 0.182764058E+06, 0.449937719E-04
0., 0.578852181E+02,-0.104480206E-02, 0., 0.209333115E+06,-0.690541063E-04
C ** TX MATRICES: FTX= 0.50000E+04
C For historical reasons, let's document data that existed prior to 25 Nov
C 2000. This former use of free-format is being replaced by the following
C use of fixed format that was copied from newer JMarti [T] : 3E26.0
C 0.87142042642E+00, 0.,-0.68364313402E+00, 0., 0.99967944506E+00, 0.
C -0.42048418843E+00, 0.,
C 0.88407538153E+00, 0., 0.55548365361E+00, 0., 0.41119257618E+00, 0.
C 0.99968059968E+00, 0.,
C 0.99996548754E+00, 0.,-0.95428680565E+00, 0.,-0.36537856597E+00, 0.
C 0.81709827065E-01, 0.,
C 0.99535953259E+00, 0., 0.99993400837E+00, 0.,-0.14508479886E+00, 0.
C -0.28684188937E+00, 0.,
C 0.13749143884E+00, 0.,-0.13487437424E+00, 0., 0.66382281216E+00, 0.
C -0.32097010954E+00, 0.,
C 0.13713603525E+00, 0., 0.10502753896E+00, 0., 0.25667806258E+00, 0.
C 0.72423739519E+00, 0.,
C 0.38040412741E+00, 0.,-0.42207400039E+00, 0.,-0.56925331264E+00, 0.
C 0.13593347464E+00, 0.,
C 0.38065940953E+00, 0., 0.44863255553E+00, 0.,-0.23815031667E+00, 0.
C -0.49588955641E+00, 0.,
NO IMAGINARY PART: { This allows us to omit the imaginary part, which is zero:
0.87142042642E+00 -0.68364313402E+00 0.99967944506E+00
-0.42048418843E+00
0.88407538153E+00 0.55548365361E+00 0.41119257618E+00
0.99968059968E+00
0.99996548754E+00 -0.95428680565E+00 -0.36537856597E+00
0.81709827065E-01
0.99535953259E+00 0.99993400837E+00 -0.14508479886E+00
-0.28684188937E+00
0.13749143884E+00 -0.13487437424E+00 0.66382281216E+00
-0.32097010954E+00
0.13713603525E+00 0.10502753896E+00 0.25667806258E+00
0.72423739519E+00
0.38040412741E+00 -0.42207400039E+00 -0.56925331264E+00
0.13593347464E+00
0.38065940953E+00 0.44863255553E+00 -0.23815031667E+00
-0.49588955641E+00
C &&&&&&&&&&&&&&&&&&&&&& END UNIT-13 CARDS PUNCHED BY DC-65 .
TOWER 2.
BUSB-1 2. .03
BUSB-2 2. .03
BLANK card ends all branch cards
BUSB-3TOWER .000995 1. 3
BLANK cards ends all switch cards (here, just one, representing the fault)
14SA-3 533000. 50. 0.0 -1.
14SA-4 533000. 50. 180. -1.
C --------------+------------------------------
C From bus name | Names of all adjacent busses.
C --------------+------------------------------
C SA-3 |BUSA-3*
C BUSA-3 |SA-3 * FIL1*BUSB-3*
C SA-4 |BUSA-4*
C BUSA-4 |SA-4 * FIL3*BUSB-4*
C BUSA-1 |TERRA *BUSB-1*
C BUSA-2 |TERRA *BUSB-2*
C FIL1 |TERRA *TERRA *TERRA *TERRA *TERRA *BUSA-3*
C FIL2 |TERRA *TERRA *
C FIL3 |BUSA-4* FIL4* FIL4* FIL4* FIL4* FIL4*
C FIL4 |TERRA *TERRA * FIL3* FIL3* FIL3* FIL3* FIL3*
C BUSB-1 |TERRA *BUSA-1*
C BUSB-2 |TERRA *BUSA-2*
C BUSB-3 |BUSA-3*TOWER *
C BUSB-4 |BUSA-4*
C TOWER |TERRA *BUSB-3*
C TERRA |BUSA-1*BUSA-2* FIL1* FIL1* FIL1* FIL1* FIL1* FIL2* FIL2*
C --------------+------------------------------
BLANK card ends source cards
C Total network loss P-loss by summing injections = 1.427244164239E+07
C 1st inject: SA-3 533000. 533000. 27.206416840329 1511.633587503
C 1st inject: 0.0 0.0 1511.3887368079 88.9687335
C Last inject: SA-4 -533000. 533000. -26.3487112737 1515.8296473689
C Last inject: 6527332453E-10 180.0000 -1515.600628549 -90.9959866
C
C Step Time BUSB-3 FIL2 TOWER FIL4 FIL3
C TOWER
C
C BUSB-1 BUSA-2 BUSA-1 BUSA-4 BUSA-3
C
C 0 0.0 .1012112E7 0.0 0.0 630.150593 3939.80139
C 161.608804 147.910277 -155.53173 -928190.23 927091.892
C 1 .1E-4 .1011754E7 0.0 0.0 638.819613 4046.20822
C 186.778615 22.6361792 -114.24647 -927914.88 927042.528
1
C Last step begins: 2000 .02 0.0 0.0 -3274.1701 1954.38677 -19609.364
C Last step begins: 250.022661 -61.141145 257.948281 -951899.68 138624.063
C Overall simulation peak node voltage = 1.31694423E+06 occurred at bus "BUSB-4"
C Variable maxima : .1012112E7 0.0 6501.36843 40746.7902 124200.262
C Variable maxima : 8139.84047 3396.88809 18096.1397 .1056252E7 927119.833
C Times of maxima : 0.0 0.0 .00101 .00951 .0024
C Times of maxima : .00102 .0023 .00227 .01074 .9E-4
C Final 5 maxima : .1056252E7 927119.833 533000. 533000. 3250.68422
C Associated times : .01074 .9E-4 0.0 .01 .00101
PRINTER PLOT
2One section of 375 km.
144 .6 0.0 7.0 BUSB-3BUSB-4 { Axis limits: (-1.087, 1.012)
CALCOMP PLOT
144 .7 0.0 7.0 BUSB-3BUSB-4
BLANK card ending all plot cards
BEGIN NEW DATA CASE
C 2nd of 2 subcases halts immediately following the illustration of 5-phase
C Semlyen branch card input. To the data of the preceding real case we add
C an artificial 5th phase that is a copy of the 4th. But both terminal
C nodes are grounded, to illustrate that such use should be possible. As
C for [Tv] and [Ti], they were artificially constructed by making column 5
C a copy of column 4 and row 5 a copy of row 4. It is the interpretation of
C these matrices that are being illustrated (it looks OK). WSM, 25 Nov 2000
PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 9 output columns)
PEAK VOLTAGE MONITOR
POWER FREQUENCY, 50.,
.000010 .020 50. 50.
1 5 1 1 1 -1
5 5 20 20 100 100 1000 200
SA-3 BUSA-3 .1260.75
SA-4 BUSA-4 .1260.75
BUSA-1 2. .03
BUSA-2 2. .03
BUSA-3 FIL1 .1 47.12
FIL1 .1 90.76
FIL1 3. 74.24
FIL1 6.166.94 94.25
FIL1 12.314.16 11.06
FIL1 12.314.16 5.87
FIL2 .01 1021.
FIL2 .1 14.14
BUSA-4 FIL3 .1 47.12
FIL3 FIL4 .1 90.76
FIL3 FIL4 3. 74.24
FIL3 FIL4 6.166.94 94.25
FIL3 FIL4 12.314.16 11.06
FIL3 FIL4 12.314.16 5.87
FIL4 .01 1021.
FIL4 .01 14.14
C &&&&&&&&&&&&&&&&&&&&&& BEGIN UNIT-13 CARDS PUNCHED BY DC-65 .
C #DATA FOR LINE SECTION 1: LENGTH = 0.2330096E+03
C METRIC
C 1 .2 .273 4 1.7653 0. 39. 20.
C 2 .2 .273 4 1.7653 30. 39. 20.
C 3 .392 .0585 4 3.1750 2.3 27.5 8.5 50.8 45.
C 4 .392 .0585 4 3.1750 32.3 27.5 8.5 50.8 45.
C
C 100. 5000. 1 375.
C 100. 50. 1 375.
C 100. 5. 1 375. 9 21
C
-1BUSA-1BUSB-1 0.74300E-02 0.13695E-02 1 1 14 6 5
0.24524677E+02 0.13280541E+03 0.00000000E+00 0.29592165E-02 0.50000000E+02 0
0., 0.301017131E+03, 0.240474368E-01, 0., 0.911921587E+04, 0.866135515E+01
0., 0.142378188E+03, 0.467439158E-01, 0., 0.597875358E+05, 0.157341296E+05
0., 0.129882539E+04, 0.158087056E+00, 0., 0.602061446E+05, 0.350951915E+05
0., 0.108960413E+05,-0.108480154E+02, 0., 0.600607627E+05,-0.503745837E+05
1., 0.289467030E+05, 0.211443554E+00, -1.,-0.537514515E+04,-0.421452862E+01
0., 0.178074955E+04, 0.254168820E+00, 0., 0.315149545E+04, 0.732001459E+00
0., 0.239325862E+04,-0.862687875E-01, 0., 0.618967874E+05,-0.453102444E+03
0., 0.619168960E+01,-0.420104801E-02, 0., 0.100000000E+01, 0.000000000E+00
0., 0.157926984E+03,-0.556405973E-03, 0., 0.179349724E+04,-0.616894886E-03
0., 0.154682388E+05,-0.740819400E-03, 0., 0.298676381E+06,-0.131478773E-02
-1BUSA-2BUSB-2 0.77523E-02 0.12698E-02 2 2 14 6 5
0.51203735E+01 0.54502236E+02 0.00000000E+00 0.31183848E-02 0.50000000E+02 0
0., 0.313224674E+02,-0.663044870E-02, 0., 0.210626498E+06, 0.432723088E+01
0., 0.182100421E+02, 0.260951895E-01, 0., 0.849632630E+06, 0.182628050E-02
0., 0.227829005E+05, 0.478630198E+00, 0., 0.102901020E+07,-0.212641132E-02
0., 0.124533203E+06, 0.211022159E+02, 0., 0.157231418E+07, 0.281738783E-03
0., 0.142701493E+06,-0.267953188E+02, 0., 0.181740797E+07,-0.719862890E-04
1., 0.398474670E+05, 0.418428806E+00, -1.,-0.214066736E+05, 0.125874270E+00
1., 0.151593680E+06, 0.515504958E+00, -1.,-0.147087464E+06, 0.267749108E+00
0., 0.365121200E+01,-0.350144906E-02, 0., 0.100000000E+01, 0.000000000E+00
0., 0.249321282E+04,-0.431759968E-04, 0., 0.182369337E+06,-0.242606205E-04
0., 0.239163482E+02,-0.406028697E-02, 0., 0.207918724E+06,-0.123123552E-03
-1BUSA-3BUSB-3 0.32619E-02 0.12582E-02 3 3 14 6 5
0.53536833E+02 0.12389964E+03 0.00000000E+00 0.12908552E-02 0.50000000E+02 0
0., 0.134598777E+06,-0.150159946E-01, 0., 0.140301928E+07,-0.753366163E-01
0., 0.284966044E+02, 0.819451304E-01, 0., 0.666744070E+07, 0.157059296E+04
0., 0.272105388E+05, 0.213806667E+00, 0., 0.782363486E+07, 0.971629566E+04
0., 0.272559469E+06,-0.223283773E+00, 0., 0.799352013E+07,-0.800863365E+06
0., 0.361659615E+06, 0.100462472E+01, 0., 0.799483818E+07, 0.792080675E+06
0., 0.491172730E+05,-0.165633449E+00, 0., 0.599908733E+07,-0.177051614E+03
0., 0.792280597E+05, 0.366175685E+00, 0., 0.698913494E+07,-0.232733391E+04
0., 0.115763918E+02,-0.117805658E-02, 0., 0.100000000E+01, 0.000000000E+00
0., 0.197441131E+04,-0.373045199E-04, 0., 0.181295123E+06, 0.445567784E-04
0., 0.976839246E+02,-0.202948282E-02, 0., 0.208740876E+06,-0.616217575E-04
-1BUSA-4BUSB-4 0.30402E-02 0.12587E-02 4 4 6 6 5
0.64988120E+02 0.13246385E+03 0.00000000E+00 0.12025589E-02 0.50000000E+02 0
0., 0.169333103E+03, 0.332803226E-01, 0., 0.566559087E+06, 0.730448580E+00
0., 0.680693506E+02, 0.708523592E-01, 0., 0.490942013E+07,-0.111902264E+01
0., 0.602799468E+05, 0.326699789E+00, 0., 0.528347745E+07, 0.957741589E+00
0., 0.646378180E+01,-0.792079833E-03, 0., 0.100000000E+01, 0.000000000E+00
0., 0.142573596E+03,-0.117927457E-02, 0., 0.182764058E+06, 0.449937719E-04
0., 0.578852181E+02,-0.104480206E-02, 0., 0.209333115E+06,-0.690541063E-04
C -1FIVEA FIVEB 0.30402E-02 0.12587E-02 5 5 6 6 5
-1 0.30402E-02 0.12587E-02 5 5 6 6 5
0.64988120E+02 0.13246385E+03 0.00000000E+00 0.12025589E-02 0.50000000E+02 0
0., 0.169333103E+03, 0.332803226E-01, 0., 0.566559087E+06, 0.730448580E+00
0., 0.680693506E+02, 0.708523592E-01, 0., 0.490942013E+07,-0.111902264E+01
0., 0.602799468E+05, 0.326699789E+00, 0., 0.528347745E+07, 0.957741589E+00
0., 0.646378180E+01,-0.792079833E-03, 0., 0.100000000E+01, 0.000000000E+00
0., 0.142573596E+03,-0.117927457E-02, 0., 0.182764058E+06, 0.449937719E-04
0., 0.578852181E+02,-0.104480206E-02, 0., 0.209333115E+06,-0.690541063E-04
C NO IMAGINARY PART: { This allows us to omit the imaginary part, which is zero:
C The preceding declaration is not needed because the one in the 1st subcase
C remains in effect. This is like a PRINTER PLOT declaration that way.
C First of two transformation matrices is 5x5 [Tv]:
0.87142042642E+00 -0.68364313402E+00 0.99967944506E+00
-0.42048418843E+00 -0.42048418843E+00
0.88407538153E+00 0.55548365361E+00 0.41119257618E+00
0.99968059968E+00 0.99968059968E+00
0.99996548754E+00 -0.95428680565E+00 -0.36537856597E+00
0.81709827065E-01 0.81709827065E-01
0.99535953259E+00 0.99993400837E+00 -0.14508479886E+00
-0.28684188937E+00 -0.28684188937E+00
0.99535953259E+00 0.99993400837E+00 -0.14508479886E+00
-0.28684188937E+00 -0.28684188937E+00 { End 5th row of [Tv]
C 2nd of two transformation matrices is 5x5 [Ti]:
0.13749143884E+00 -0.13487437424E+00 0.66382281216E+00
-0.32097010954E+00 -0.32097010954E+00
0.13713603525E+00 0.10502753896E+00 0.25667806258E+00
0.72423739519E+00 0.72423739519E+00
0.38040412741E+00 -0.42207400039E+00 -0.56925331264E+00
0.13593347464E+00 0.13593347464E+00
0.38065940953E+00 0.44863255553E+00 -0.23815031667E+00
-0.49588955641E+00 -0.49588955641E+00
0.38065940953E+00 0.44863255553E+00 -0.23815031667E+00
-0.49588955641E+00 -0.49588955641E+00 { End 5th row of [Ti]
YES, IMAGINARY PART: { Illustrate optional return to default use of imaginary
C The preceding has no effect for this data. But if the user wants to
C cancel the NO IMAGINARY PART declaration of the first subcase, this
C is the way to do it. Note such use would be required if old Semlyen
C branch cards followed these newer ones having implied zeros.
TOWER 2.
$STOP { Having illustrated 5-phase Semlyen data input (particularly [T]), halt
BEGIN NEW DATA CASE
BLANK
|