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
|
BEGIN NEW DATA CASE
C BENCHMARK DC-29
C Test of "SEMLYEN SETUP" for 138-mile, single-circuit, 500-kV overhead
C line that connects "John Day" with "Lower Monumental" (BPA substations).
C There are 2 ground wires (both eliminated), and each phase consists of a
C bundle of 2 conductors. Punched LUNIT7 branch cards contribute to DC-31.
C 1st of 3 subcases. 2nd is for Pi-equiv. output (fast). 3rd is bypassed.
C Solution changed on 23 Oct 1990 due to Mustafa Kizilcay's modifications
C (preceding commented data has been copied to bottom for preservation).
SEMLYEN SETUP
$ERASE
BRANCH JDA LMA JDB LMB JDC LMC
TOLERANCES 10 5000. { Illustration only; value of FMED actually unchanged
C 345678901234567890123456789012345678901234567890
C 60 80 230 1 7 7777 0 { Semlyen miscellaneous data
C 60 80 1 230 1 7 7777 0 { Semlyen miscellaneous data
C Guido: whenever n4=0(col 26), n1 (col 21) should be 1 !! avoid kill=152
200 150 10 230 1 7 7777 0 { Semlyen miscellaneous data
C Guido: overhead lines; rule of thumb: NOO=100-300, NPOINT=150
C Note: Guido modifies Semlyen misc. data card above on 23 Oct 1990 following
C the addition of Mustafa Kizilcay's latest changes. The old one gave
C an early KILL code. Later (November 7th), WSM and THL convert the
C "1 " of columns 21-22 to "10" to make it universal. Although the
C VAX had no trouble reading 10 from the original as I2 information,
C Sun (which agreed with Apollo, so Apollo most likely also was wrong)
C erroneously read unity! Sun seems to ignore blanks in "I" FORMATs.
LINE CONSTANTS
1.3636 .05215 4 1.602 -20.75 50. 50.
1.3636 .05215 4 1.602 -19.25 50. 50.
2.3636 .05215 4 1.602 - 0.75 77.5 77.5
2.3636 .05215 4 1.602 0.75 77.5 77.5
3.3636 .05215 4 1.602 19.25 50. 50.
3.3636 .05215 4 1.602 20.75 50. 50.
0.5 2.61 4 0.386 -12.9 98.5 98.5
0.5 2.61 4 0.386 12.9 98.5 98.5
BLANK card ending conductor cards within "LINE CONSTANTS" data
C 27. 5000. 1 138. { Transient frequency
C 27. 60.00 1 138. { Phasor solution frequency
C 27. 6.00 1 138. 6 20 { log looping
27. 5000. 138. { Transient frequency
27. 60.00 138. { Phasor solution frequency
27. 6.00 138. 6 20 { log looping
BLANK card ending frequency cards of "LINE CONSTANTS" data
C Frequency = 5.00000000E+03
C --------------------------
C Impedance matrix Admittance matrix
C (ohm/mile) (mho/mile)
C
C 3.2463317E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
C 6.3079678E+01 5.2869139E-04 -8.2803449E-05 -6.2048083E-05
C
C 2.3065452E+00 2.0246396E+00 0.0000000E+00 0.0000000E+00
C 1.4273670E+01 5.9110509E+01 5.3917160E-04 -8.2803449E-05
C
C 2.9595110E+00 2.3065452E+00 3.2463317E+00 0.0000000E+00
C 1.4114577E+01 1.4273670E+01 6.3079678E+01 5.2869139E-04
BLANK card ending "LINE CONSTANTS" data cases
C The following values are computed on the basis of exact diagonalization of the
C Mode to phase voltage Mode to phase current
C transformation matrix transformation matrix
C
C 1.000000 -1.000000 -0.266985 0.420272 -0.500000 -0.299082
C 0.000000 0.000000 0.019744 0.007566 0.000000 0.006219
C
C 0.711142 0.000000 1.000000 0.224711 0.000000 0.840545
C -0.027600 0.000000 0.000000 -0.012556 0.000000 0.015131
C
C 1.000000 1.000000 -0.266985 0.420272 0.500000 -0.299082
C 0.000000 0.000000 0.019744 0.007566 0.000000 0.006219
C
C Alpha Beta Attenuation
C Mode (neper/mi) (radians/mi) (db/mi)
C 1 0.8245225702687424E-02 0.1861008889722703E+00 0.7161712049448019E-01
C 2 0.4981192766022477E-03 0.1700760589248700E+00 0.4326609063159914E-02
C 3 0.4573686287515477E-03 0.1691150685564419E+00 0.3972653433249082E-02
C Number of frequencies = 121 Starting and ending Hz = 6.0000E+00 6.000E+06
C Fitting of the propagation step response for mode number 1 :
C Iteration Amplitude X(1) X(2) X(3) Error
C 1 7.26461134E-1 1.32035374E+1 1.17295449E+0 2.21313604E-2 7.35218983E-04
C 1 7.13648369E-1 1.45238912E+1 1.17295449E+0 2.21313604E-2 1.23450128E-03
C 1 7.02961062E-1 1.58442449E01 1.17295449E+0 2.21313604E-2 1.78841949E-03
C < < Etc. until final (23-rd) such entries, which is:
C 3 7.67204181E-1 1.05695813E+1 1.17465822E+0 2.53215686E-2 3.72675596E-05
C Successful convergence .... SUBROUTINE TDFIT. Tolerance EPS = 5.00000000E-05
C ++++ Warning! TDFIT can not adjust time-domain step response to fit exactly
C Instead, the initial fitting will be used.
C Amplitude 1 Time const. 1 Amplitude 2 Time const. 2 Delay time
C 7.6720418E-01 7.0400012E-05 2.3279582E-01 6.3345970E-04 7.7782236E-04
C Mode 1, time step = 3.72049322E-06 Symbols : "." = Semlyen approximation
C Error Time 0.0 8.93193367E-01
C .1........1........1........1........1........1..
C 0.00000E+00 7.62701E-04 .*
C 0.00000E+00 7.66422E-04 ..0
C 0.00000E+00 7.70142E-04 ..0
C 0.00000E+00 7.73863E-04 .. 0
C 0.00000E+00 7.77583E-04 .. 0
C -2.07761E-02 7.81304E-04 . . 0
C -4.84345E-03 7.85024E-04 . *
C 6.50870E-03 7.88745E-04 . 0.
C 1.41364E-02 7.92465E-04 . 0.
C 1.88276E-02 7.96186E-04 . 0.
C 2.12721E-02 7.99906E-04 . 0 .
C 2.20463E-02 8.03627E-04 . 0 .
C 2.16177E-02 8.07347E-04 . 0 .
C 2.03607E-02 8.11068E-04 . 0 .
C 1.85647E-02 8.14788E-04 . 0 .
C 1.64473E-02 8.18509E-04 . 0.
C 1.41747E-02 8.22229E-04 . 0 .
C 1.18700E-02 8.25949E-04 . 0. C 9.61526E-03 8.29670E-04 . 0.
C 7.46654E-03 8.33390E-04 . 0.
$PUNCH
C -1JDA LMA 5.72862E-03 7.77822E-04 1 1 2 2 3
C 1.76406006E+01 9.04231135E+01-1.97906661E-04 1.44978522E-03 6.00000000E+01
C 0.00000E+00 1.42045E+04 7.67204E-01 0.00000E+00 1.57863E+03 2.32796E-01
C 0.00000E+00 1.70053E+05-7.75536E-04 0.00000E+00 5.92740E+02-1.43454E-03
C -1JDB LMB 7.00175E-03 7.40982E-04 2 2 2 2 3
C 5.19844042E-01 4.16341654E+01-1.59314154E-05 1.96071950E-03 6.00000000E+01
C 0.00000E+00 1.82512E+05 9.69335E-01 0.00000E+00 2.58428E+03 3.06652E-02
C 0.00000E+00 9.60136E+03-1.29672E-04 0.00000E+00 5.91248E+01-4.28783E-04
C -1JDC LMC 4.11495E-03 7.40975E-04 3 3 2 2 3
C 2.42272029E+00 7.27673245E+01-1.71561794E-05 1.15316857E-03 6.00000000E+01
C 0.00000E+00 4.30709E+05 9.63717E-01 0.00000E+00 1.32518E+03 3.62828E-02
C 0.00000E+00 2.12922E+03-1.21998E-04 0.00000E+00 8.97571E+01-3.32339E-04
C Prior to improvement on 27 November 2000, zero imaginary parts were
C included, and this forced low precision of 6E12.5:
C 1.00000E+00 0.00000E+00 1.00000E+00 0.00000E+00-2.66985E-01 0.00000E+00
C 7.11142E-01 0.00000E+00-1.02636E-14 0.00000E+00 1.00000E+00 0.00000E+00
C 1.00000E+00 0.00000E+00-1.00000E+00 0.00000E+00-2.66985E-01 0.00000E+00
C 4.20272E-01 0.00000E+00 5.00000E-01 0.00000E+00-2.99082E-01 0.00000E+00
C 2.24711E-01 0.00000E+00-5.36721E-15 0.00000E+00 8.40545E-01 0.00000E+00
C 4.20272E-01 0.00000E+00-5.00000E-01 0.00000E+00-2.99082E-01 0.00000E+00
C Beginning 27 November 2000, the zero imaginary parts are implied,
C giving the following 3E26.18 output following a new request word. A
C big paragraph about this is planned for the April, 2001, newsletter:
C NO IMAGINARY PART
C 1.0000000000000000E+00 1.0000000000000000E+00 -2.6698449851193501E-01
C 7.1114175640039812E-01 -1.1764114701044076E-14 1.0000000000000000E+00
C 9.9999999999999945E-01 -9.9999999999999467E-01 -2.6698449851194972E-01
C 4.2027250134843613E-01 5.0000000000000333E-01 -2.9908213367540470E-01
C 2.2471123211965521E-01 -6.1556662545037000E-15 8.4054500269687149E-01
C 4.2027250134843563E-01 -4.9999999999999922E-01 -2.9908213367541620E-01
BLANK card ending "SEMLYEN SETUP" data cases
BEGIN NEW DATA CASE
C BENCHMARK DC-29
C 2nd of 3 subcases illustrates punching of Pi-equivalents on LUNIT8.
C Feature added by Mustafa Kizilcay of University of Hannover, Oct 88.
C There is no rational function fitting, so execution is fast. Also,
C the time spent in LINE CONSTANTS is minimal: 2 decades, 1 pt/dec.
C { Running interactively ("SPY DC29.") gave debug on second subcase:
DIAGNOSTIC { Somehow needed to erase overlay-25 diagnostic. ???? 9 July 1989
SEMLYEN SETUP
$OPEN, UNIT=8 FILE=dc29out.lis STATUS=UNKNOWN FORM=FORMATTED ! { Pi-equivalents will go here
C 345678901234567890123456789012345678901234567890
60 80 230 3 1 7 4100 4100 { Semlyen miscellaneous data
LINE CONSTANTS
1.3636 .05215 4 1.602 -20.75 50. 50.
1.3636 .05215 4 1.602 -19.25 50. 50.
2.3636 .05215 4 1.602 - 0.75 77.5 77.5
2.3636 .05215 4 1.602 0.75 77.5 77.5
3.3636 .05215 4 1.602 19.25 50. 50.
3.3636 .05215 4 1.602 20.75 50. 50.
0.5 2.61 4 0.386 -12.9 98.5 98.5
0.5 2.61 4 0.386 12.9 98.5 98.5
BLANK card ending conductor cards within "LINE CONSTANTS" data
27. 5000. 138. { Transient frequency
27. 50.00 138. { Phasor solution frequency
27. 6.00 138. 2 1 { log looping
BLANK card ending frequency cards of "LINE CONSTANTS" data
BLANK card ending "LINE CONSTANTS" data cases
$CLOSE, UNIT=8 STATUS=KEEP { Disconnect file now containing Pi-equivalents
BLANK card ending "SEMLYEN SETUP" data cases
C Transfer [Y] of equivalent Pi for F = 5.0000E+01 Hz. 23-Oct-90 19.45.28
C 1.151519555275E-3 -1.137441619542E-02 -1.060634758717E-04 3.491463556430E-03
C 3.537076770591E-4 3.210624087768E-03 -1.060634758717E-04 3.491463556430E-03
C 2.050900981651E-3 -1.162860843144E-02 -1.060634758717E-04 3.491463556430E-03
C 3.537076770591E-4 3.210624087768E-03 -1.060634758717E-04 3.491463556430E-03
C 1.151519555275E-3 -1.137441619542E-02
C Shunt [Y] of equivalent Pi for F = 5.0000E+01 Hz.
C -8.541869714622E-6 3.665747845891E-04 6.607667284137E-06 -5.660142619185E-05
C -8.648581530269E-6 -4.295696271522E-05 6.607667284137E-06 -5.660142619185E-05
C -2.342599469890E-5 3.722080025785E-04 6.607667284137E-06 -5.660142619185E-05
C -8.648581530269E-6 -4.295696271522E-05 6.607667284137E-06 -5.660142619185E-05
C -8.541869714623E-6 3.665747845891E-04
C Transfer [Y] of equivalent Pi for F = 6.0000E+00 Hz. 23-Oct-90 19.45.28
C 3.262709023195E-2 -7.809118174262E-02 -1.341276675346E-02 2.522775392732E-02
C -1.150007192882E-2 2.343312809664E-02 -1.341276675346E-02 2.522775392732E-02
C 3.436831351543E-2 -7.888248987903E-02 -1.341276675347E-02 2.522775392732E-02
C -1.150007192882E-2 2.343312809664E-02 -1.341276675347E-02 2.522775392732E-02
C 3.262709023195E-2 -7.809118174262E-02
C Shunt [Y] of equivalent Pi for F = 6.0000E+00 Hz.
C -3.427096070385E-8 4.377964661972E-05 2.594787689802E-08 -6.855484667815E-06
C -3.570745329639E-8 -5.136885722574E-06 2.594787689802E-08 -6.855484667815E-06
C -9.461439943596E-8 4.464722732424E-05 2.594787689809E-08 -6.855484667815E-06
C -3.570745329639E-8 -5.136885722574E-06 2.594787689809E-08 -6.855484667815E-06
C -3.427096070393E-8 4.377964661972E-05
C Transfer [Y] of equivalent Pi for F = 6.0000E+01 Hz. 23-Oct-90 19.45.28
C 9.441668531285E-4 -9.566451344743E-03 -7.083376044220E-05 2.864853469711E-03
C 3.766916044259E-4 2.651337447173E-03 -7.083376044220E-05 2.864853469711E-03
C 1.809395347023E-3 -9.800403930488E-03 -7.083376044219E-05 2.864853469711E-03
C 3.766916044259E-4 2.651337447173E-03 -7.083376044219E-05 2.864853469711E-03
C 9.441668531284E-4 -9.566451344743E-03
C Shunt [Y] of equivalent Pi for F = 6.0000E+01 Hz.
C -1.202671448681E-5 4.408181056666E-04 9.416194169739E-06 -6.772347154278E-05
C -1.218469464166E-5 -5.164200772177E-05 9.416194169739E-06 -6.772347154278E-05
C -3.311677630903E-5 4.467905957272E-04 9.416194169739E-06 -6.772347154278E-05
C -1.218469464166E-5 -5.164200772177E-05 9.416194169739E-06 -6.772347154278E-05
C -1.202671448681E-5 4.408181056666E-04
C Transfer [Y] of equivalent Pi for F = 6.0000E+02 Hz. 23-Oct-90 19.45.28
C -3.712845487162E-3 -5.603787237254E-03 -1.561079196979E-03 4.482090715102E-03
C -3.184914860379E-3 5.584231439014E-03 -1.561079196979E-03 4.482090715102E-03
C -4.822417739649E-4 -7.822757815584E-03 -1.561079196979E-03 4.482090715102E-03
C -3.184914860379E-3 5.584231439014E-03 -1.561079196979E-03 4.482090715102E-03
C -3.712845487162E-3 -5.603787237254E-03
C Shunt [Y] of equivalent Pi for F = 6.0000E+02 Hz.
C 7.516929477229E-3 1.093259393216E-02 3.157367799464E-03 -8.756571161718E-03
C 6.427094062687E-3 -1.089709886987E-02 3.157367799464E-03 -8.756571161718E-03
C 1.100477005468E-3 1.520008342172E-02 3.157367799464E-03 -8.756571161718E-03
C 6.427094062687E-3 -1.089709886987E-02 3.157367799464E-03 -8.756571161718E-03
C 7.516929477229E-3 1.093259393216E-02
BEGIN NEW DATA CASE
C 3rd of 3 subcases illustrates SEMLYEN SETUP for two of the actual
C 3 phases of the 500-kV line data of the 1st subcase. Results will
C be used in the 4th subcase of DC-41 to illustrate selective branch
C current output. This addition is made 11 December 2000.
SEMLYEN SETUP
$ERASE
BRANCH JDA LMA JDB LMB { Note only 2 of the actual 3 phases are used
TOLERANCES 10 5000. { Illustration only; value of FMED actually unchanged
200 150 10 230 1 7 7777 0 { Semlyen miscellaneous data
LINE CONSTANTS
1.3636 .05215 4 1.602 -20.75 50. 50.
1.3636 .05215 4 1.602 -19.25 50. 50.
2.3636 .05215 4 1.602 - 0.75 77.5 77.5
2.3636 .05215 4 1.602 0.75 77.5 77.5
C Following two conductor cards for phase "c" are omitted to create 2-phase line
C 3.3636 .05215 4 1.602 19.25 50. 50. ---- omit phase "c"
C 3.3636 .05215 4 1.602 20.75 50. 50. ---- omit phase "c"
0.5 2.61 4 0.386 -12.9 98.5 98.5
0.5 2.61 4 0.386 12.9 98.5 98.5
BLANK card ending conductor cards within "LINE CONSTANTS" data
27. 5000. 138. { Transient frequency
27. 60.00 138. { Phasor solution frequency
27. 6.00 138. 6 20 { log looping
BLANK card ending frequency cards of "LINE CONSTANTS" data
BLANK card ending "LINE CONSTANTS" data cases
$PUNCH
BLANK card ending "SEMLYEN SETUP" data cases
BLANK { Note this extra blank will terminate execution, to avoid 3rd subcase
BLANK { 14 Dec 94, we add more, to make 5 even when comments are destroyed
BLANK { 14 Dec 94, we add more, to make 5 even when comments are destroyed
BLANK { 14 Dec 94, we add more, to make 5 even when comments are destroyed
BLANK { 14 Dec 94, we add more, to make 5 even when comments are destroyed
BLANK { 14 Dec 94, we add more, to make 5 even when comments are destroyed
BEGIN NEW DATA CASE
C 4th of 4 subcases, not normally executed, is Rule Book example of XXII-D.
C This was 3rd of 3 until 11 December 2000 when preceding subcase was added
SEMLYEN SETUP
$ERASE
BRANCH SORCA LINEA SORCB LINEB SORCC LINEC
60 80 2 1 7 7777 1637
LINE CONSTANTS
1.3159 .192 4 .858 -10.67 31.5 31.5
2.3159 .192 4 .858 0.00 31.5 31.5
3.3159 .192 4 .858 10.67 31.5 31.5
0.5 9.48 4 0.313 -5.33 40.5 40.5
0.5 9.48 4 0.313 5.33 40.5 40.5
BLANK card ending conductor cards within "LINE CONSTANTS" data
1000. 5000. 18.65 { Transient frequency
1000. 60.00 18.65 { Phasor solution frequency
1000. 60.00 18.65 4 21 { log looping
BLANK card ending frequency cards of "LINE CONSTANTS" data
BLANK card ending "LINE CONSTANTS" data cases
$PUNCH
BLANK card ending "SEMLYEN SETUP" data cases
BEGIN NEW DATA CASE
BLANK
|
