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
|
BEGIN NEW DATA CASE
C BENCHMARK DC-52
C "FREQUENCY SCAN" use, with subsequent plotting of phasors vs. frequency
C Unlike DC-51, this case uses geometric spacing: f-min, max = 1.0, 15 Hz;
C two points/decade. Also, there is full phasor solution output, & CalComp
PRINTED NUMBER WIDTH, 9, 2,
FREQUENCY SCAN, 1.0, 0.0, 15., 2,
1.0 0.0
1 1 1 1 { Note request for phasor branch flows
SWIT LOAD 10.
LOAD 1000.
-1SWIT OPEN .3055 5.82 .012 138.
BLANK card ending all branches
GEN SWIT -1.
BLANK card ending all switch cards
BOTH POLAR AND RECTANGULAR { Request for (in order): mag, angle, real, imag
14GEN 1.0 0.5 -1.
BLANK card ending source cards
C Column headings for the 2 output variables follow. These are divided among the 3 possible FS variable classes as follows ....
C First 2 output variables are electric-network voltage differences (upper voltage minus lower voltage);
C For each variable, magnitude is 1st, angle is 2nd, real part is 3rd, and imaginary is 4th. All 4 are labeled identically, note.
C Step F [Hz] GEN GEN GEN GEN LOAD LOAD LOAD LOAD
C 1 1.0 1.0 0.0 1.0 0.0 .532018 57.8581 .283043 .450477
C 2 3.1622777 1.0 0.0 1.0 0.0 .893248 26.7158 .797892 .401573
C 3 10. 1.0 0.0 1.0 0.0 .98757 9.04306 .975295 .155223
C 4 31.622777 1.0 0.0 1.0 0.0 .998736 2.88122 .997473 .050202
C Preceding was new output. Following comments are left from pre-1998 output:
C Total network loss P-loss by summing injections = 3.584784087208E-02
C Inject: SWIT 1.0 1.0 .07169568174416 .08466776778383 .03584784087208
C Inject: 0.0 0.0 -.0450373192002 -32.1359449 .02251865960008
C ----- Output vector for step number 1. Frequency = 1.00000000E+00 Hz.
C 1.0000000E+00 0.0000000E+00 5.3201804E-01 5.7858092E+01
C 1.0000000E+00 0.0000000E+00 2.8304320E-01 4.5047724E-01
GEN LOAD
BLANK card ends output requests (just node voltages, for FREQUENCY SCAN)
C Total network loss P-loss by summing injections = 1.272220336009E-04
C ----- Output vector for step number 4. Frequency = 3.16227766E+01 Hz.
C 1.0000000E+00 0.0000000E+00 9.9873589E-01 2.8812203E+00
C 1.0000000E+00 0.0000000E+00 9.9747337E-01 5.0202049E-02
F-SCAN COMPONENTS ANGLE { Access a single "angle" next
C 197 .2 0.0 1.6 -20. 60.LOAD angle
147 .2 0.0 1.6 -20. 60.LOAD
PRINTER PLOT
C 197 .4 0.0 1.6 LOAD angle { Axis limits: (0.000, 5.786)
147 .4 0.0 1.6 LOAD { Axis limits: (0.000, 5.786)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 2nd of 3 subcases illustrates LINE MODEL FREQUENCY SCAN (LMFS) for
C a cable (CORE ONLY, 3 PHASE). Branch cards generated by JMARTI SETUP.
C DISK PLOT DATA { Has no effect as explained in Can/Am EMTP News of April, 1992
C DIAGNOSTIC 0 0 0 0 0 0 0 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9
LINE MODEL FREQUENCY SCAN, 6, , 60000, 10, , , { f-begin, f-end, points/decade
C LINE MODEL FREQUENCY SCAN, 6, , 599.9, 1, , , { f-begin, f-end, points/decade
40.E-9 -.0600 60. 1.E-15
1 1 0 0 0 0 0 2
C IOUT KSSOUT ICAT
$INCLUDE, []dc52fs3.dat { 1st of 2 $INCLUDE contains line section to be tested
$INCLUDE, []dc51fs3z.dat, JDA###, LMA###, JDB###, LMB###, JDC###, LMC###,
C []dc51fs3p.dat
C Unused! ----- 2nd of 2 file names is not used for this single-circuit case
CABLE CONSTANTS
C Dimensioning of CABLE PARAMETERS (CP) is under manual control as the
C following declaration illustrates. Normally, the user will omit such
C information of columns 33-48, and the program will use available working
C space to provide for the most burdensome possible case. While consistent,
C this is wasteful. The resulting limit LV on phases will be smaller than
C necessary. For data that is in danger of overflowing, an intelligent user
C should consider the following manual allocation. In the following, LVPLUS
C is the number of extra phases for use within SUBROUTINE CYMTRX. A value in
C excess of 2 * LV, or non-positive, will result in the default 2 * LV.
C But most data requires nothing like (3*LV,3*LV) for storage. In the
C following, use of (39,39) is illustrated. No matter what is used, the
C program should be protected. If inadequate, an error message should halt
C execution within CYMTRX. Of course, the present data requires nothing
C like (33,33) and (39,39). Rather, these numbers come from CABLE32.DAT
C as used for testing by BPA's Dr. Tsu-huei Liu on 9 May 2000. The following
C declaration is adequate to handle that 33-conductor data:
C 345678901234567890123456789012<---- LV LVPLUS
CABLE PARAMETERS 33 6
2 -1 3 0 1 1 1
1 1 1
25.4E-3 45.6E-3
6.8912E-8 1. 1. 3.52
25.4E-3 45.6E-3
6.8912E-8 1. 1. 3.52
25.4E-3 45.6E-3
6.8912E-8 1. 1. 3.52
0.75 0.0 0.75 0.3 0.75 0.15
C EARTH RESISTIVITY AND FREQUENCY CARDS
20. 50. 90.1 { 1st of 3 is Xform matrix [T]
BLANK card ending frequency cards of "CABLE CONSTANTS" data
PRINTER PLOT
186 .5 0.0 5.0 GENIA MAG Marti Section Z in Ohms
196 .5 0.0 5.0 GENIA MAG Equivalent Pi Z in Ohms
196 .5 0.0 5.0 GENIB % MAG Error Percent
BLANK card ending batch-mode (CalComp) plot cards
BEGIN NEW DATA CASE
C 3rd of 3 subcases illustrates the use of MODELS with FREQUENCY SCAN.
C Uses same circuit as first subcase. From Laurent Dube, 15 December 1997
C For practical usage, there would be many more frequency points, and the
C printout of phasor branch flows would be suppressed (1 --> 0 in column 32).
FREQUENCY SCAN, 1.0, 0.0, 15., 2,
1.0 0.0
1 1 1 1 { Note request for phasor branch flows
C ==============================================================================
MODELS
INPUT v_Re_swit {v(swit)} -- real part of voltage at node SWIT
v_Im_swit {imssv(swit)} -- imag part of voltage at node SWIT
v_Re_load {v(load)} -- real part of voltage at node LOAD
v_Im_load {imssv(load)} -- imag part of voltage at node LOAD
i_Re_swit {i(swit)} -- real part of current at switch SWIT
i_Im_swit {imssi(swit)} -- imag part of current at switch SWIT
state_swit {switch(swit)} -- state of switch SWIT (0=open, 1=closed)
MODEL m1 -----------------------------------------------------------------------
VAR pass -- to verify that its value is kept between passes
INPUT v_Re_1 {dflt: 0} -- to verify that all types of inputs are
v_Im_1 {dflt: 0} -- accessed correctly
v_Re_2 {dflt: 0} -- (machine quantities are not tested here)
v_Im_2 {dflt: 0}
i_Re_1 {dflt: 0}
i_Im_1 {dflt: 0}
state_1 {dflt: 0}
INIT pass:=0 ENDINIT
EXEC
pass:=pass+1
write('************** In models, pass= ', pass)
write(' freq= ', t) -- <<=== <<====
write(' Re{v1}, Im{v1} = ', v_Re_1, ', ', v_Im_1 )
write(' Re{v2}, Im{v2} = ', v_Re_2, ', ', v_Im_2 )
write(' Re{i1}, Im{i1} = ', i_Re_1, ', ', i_Im_1 )
write(' state_1 = ', state_1 )
ENDEXEC
ENDMODEL
USE m1 as m1
INPUT v_Re_1 := v_Re_swit
v_Im_1 := v_Im_swit
v_Re_2 := v_Re_load
v_Im_2 := v_Im_load
i_Re_1 := i_Re_swit
i_Im_1 := i_Im_swit
state_1 := state_swit
ENDUSE
RECORD v_Re_swit AS vre1
v_Im_swit AS vim1
v_Re_load AS vre2
v_Im_load AS vim2
i_Re_swit AS ire1
i_Im_swit AS iim1
m1.state_1 AS state1
ENDMODELS
C Begin branches ===============================================================
SWIT LOAD 10.
LOAD 1000.
-1SWIT OPEN .3055 5.82 .012 138.
BLANK card ending all branches
GEN SWIT -1.
BLANK card ending all switch cards
14GEN 1.0 0.5 -1.
BLANK card ending source cards
GEN LOAD
BLANK card ends output requests (just node voltages, for FREQUENCY SCAN)
PRINTER PLOT
C 197 .5 0.0 1.6 0. 1.LOAD mag
147 .5 0.0 1.6 0. 1.LOAD
C C 197 .2 0.0 1.6 -20. 60.LOAD angle
C C 197 .2 0.0 1.6 0. 1.LOAD real
C C 197 .2 0.0 1.6 0. 1.VRE2 MODELS
C C 197 .2 0.0 1.6 0. 1.LOAD imag
C C 197 .2 0.0 1.6 0. 1.VIM2 MODELS
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
|
