summaryrefslogtreecommitdiffstats
path: root/benchmarks/dc4.dat
blob: ccead6d59de57d9300021541a39ea050c4305088 (plain)
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
BEGIN NEW DATA CASE
C     BENCHMARK DC-4
C        Solution agrees (key variables, anyway) with DC-5 saturable TRANSFORMER
C        Here,  there is no transformer,  since the secondary impedance has been
C        reflected to the primary.   Also, the pseudo-nonlinear reactance of the
C        saturable  TRANSFORMER  is  here  represented by a "true" nonlinearity.
C        2nd subcase will illustrate residual flux for Type-93 nonlinear induct.
C        Etc.  5 more subcases were added 21 March 2006 bringing total to 7.
PRINTED NUMBER WIDTH, 13, 2,  { Request maximum precision (for 8 output columns)
    .010     6.0
       1       1       1       1       1      -1
       5       5      20      20     100     100
  BRANCH NAME:First   { Even though name could go on next card, use this instead
  GEN   TRAN                 5.0  5.E4                                         3
  TRAN        NAME  R-mag   1.E4
93TRAN        NAME  Magnet  .005   30.                                         3
C            -5.0           -100.         { No longer needed 3-rd quadrant point
C             -.1            -50.         { No longer needed 3-rd quadrant point
C            -.02            -45.         { No longer needed 3-rd quadrant point
C            -.01            -40.         { No longer needed 3-rd quadrant point
C           -.005            -30.         { No longer needed 3-rd quadrant point
             0.0             0.0  { 1st point being origin is request to reflect
            .005             30.
             .01             40.
             .02             45.
             .10             50.
             5.0            100.
            9999
  TRAN  LOADG               255.  5.E4                                         3
C   Note: original fixed-format card (next comment) converted to free-format:
C LOADG                    1.E-6
0,LOADG, , , ,             1.E-6,  0.0, 0.0, , , , , ,
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14GEN            70.  .1591549                                     -1.
BLANK card terminating program source cards.
C       Total network loss  P-loss  by summing injections =   8.286714400785E+00
C Inject:  GEN     70.    70.  .23676326859385  .25769284993889  8.2867144007848
C Inject:          0.0    0.0  -.1017288531066      -23.2514964  3.5605098587304
C          ---- Initial flux of coil  "TRAN  "  to  "      "  =  -1.13295190E+01
C Step    Time        GEN          TRAN         TRAN         GEN          TRAN
C                     TRAN         TERRA        LOADG
C    0     0.0   6.270257621  63.72974238  63.72974215          70.  63.72974238
C    1     .01   6.156651781  63.83984825  63.83984802  69.99650003  63.83984825
  GEN   TRAN
BLANK card ending program output-variable requests.
C Last step begins: 600  6.0  12.26276725  54.94912105  54.94912087  67.21188829
C Last step continued .... :  54.94912105  -.004557817  .1742561973  .1733191018
C Variable max : 35.83743967  64.72841365   64.7284134          70.  64.72841365
C Times of max :        5.37          .18          .18          0.0          .18
C Variable min : -34.5302967  -63.9279756  -63.9279753  -69.9999111  -63.9279756
C Times of min :         2.2         3.35         3.35         3.14         3.35
  PRINTER PLOT
 193 .5 0.0 6.0         TRAN  LOADG GEN   TRAN   { Axis limits: (-3.494,  3.289)
 193 1. 0.0 6.0         TRAN                     { Axis limits: (-3.253,  3.031)
BLANK card ending all plot cards
BEGIN NEW DATA CASE
C          2nd of 8 subcases illustrates residual flux for the Type-93 nonlinear
C          inductor. The solution is a little different because the segment that
C          passes through the origin has been modified to draw no current  (this
C          is the poor man's approximation to hysteresis).   Usage of a residual
C          flux that  is a little different than the phasor solution also alters
C          the solution somewhat.  Appearance  of  the output column of inductor
C          current is altered since whenever operation  is on the first segment,
C          the current now is identically zero.
PRINTED NUMBER WIDTH, 13, 2,  { Request maximum precision (for 8 output columns)
C         4 May 2007,  add following  $PREFIX  after modification of  OVER1  and
C         CIMAGE  that moves  CHAR*80 PREFIX  into  BLKCOM  for use by  $INSERT.
C         This allows remote execution, from \UNFORM, without data modification.
$PREFIX, []    { $INCLUDE files are located in same place as this main data file
    .010     6.0
       1       1       1       1       1      -1
       5       5      20      20     100     100
  GEN   TRAN                 5.0  5.E4                                         3
  TRAN        NAME  R-mag   1.E4
C      Next, move entire nonlinear element into file to demonstrate new feature:
$INSERT, dc4bincl.dat       { Like $INCLUDE only dynamic (no burden on LIMCRD)
  TRAN  LOADG               255.  5.E4                                         3
C   Note: original fixed-format card (next comment) converted to free-format:
C LOADG                    1.E-6
0,LOADG, , , ,             1.E-6,  0.0, 0.0, , , , , ,
BLANK card ending program branch cards.
BLANK card terminating program switch cards (none, for this case)
14GEN            70.  .1591549                                     -1.
BLANK card terminating program source cards.
C   20 July 2009, convert the following N.L. element I.C. card from the original
C   use of terminal node names  (TRAN, TERRA)  to the use of branch name Magnet.
C   Yes, within $INSERT file dc4bincl.dat this NL elem. branch name is declared.
C   Just today, such alternative use of N.L. names is being added to ATP.  WSM
C  The following is an initial-condition card for the Type-93 nonlinear inductor
C 4TRAN                    0.0    -11.5508593 { --- Used to reaffirm phasor flux
C 4TRAN                    0.0          -13.0    { Residual flux differs somewhat
 4Name: Magnet            0.0          -13.0    { Residual flux differs somewhat
  GEN   TRAN                                { Request for 2 node voltage outputs
C  Step    Time         GEN          TRAN         TRAN         GEN          TRAN
C                       TRAN         TERRA        LOADG
C     0    0.0   5.816303518  64.18369648  64.18369625          70.  64.18369648
C     1    .01   5.700507029    64.295993  64.29599277  69.99650003    64.295993
C     2    .02    5.58414044  64.40186003   64.4018598  69.98600047  64.40186003
BLANK card ending program output-variable requests.
C   600    6.0   12.62786557  54.58402272  54.58402255  67.21188829  54.58402272
C Variable max :  37.1903484  65.21465971  65.21465947          70.  65.21465971
C Times of max :         5.4          .18          .18          0.0          .18
C Variable min : -32.6413814  -64.4282966  -64.4282964  -69.9999111  -64.4282966
C Times of min :        2.16         3.35         3.35         3.14         3.35
  PRINTER PLOT
 193 .5 0.0 6.0         TRAN  LOADG GEN   TRAN   { Axis limits: (-3.701,  3.055)
 193 1. 0.0 6.0         TRAN                     { Axis limits: (-3.474,  2.777)
BLANK card ending all plot cards
BEGIN NEW DATA CASE
C          3rd of 8 subcases begins illustration of enhanced  $INSERT  which
C          is added to the UTPF on  21 March 2006.  Arguments and "/"-card
C          sorting and nesting (2 or more levels) all are allowed.  The data
C          now being added was named  DC17NEW.DAT  during early testing because
C          initial testing involved the configuration of  DC-17.  WSM.
$INSERT,  dc4ins.dat,   9,   .005,   TRAN##,
    .005     4.0
       1       1       1       1       1      -1
       5       5      20      20
  TRAN  LOADG               255.  5.E4                                         3
  LOADG                    1.E-6
BLANK card ending  BRANCH  cards
BLANK card ending  SWITCH  cards (none exists, for this case)
BLANK card ending  SOURCE  cards
  GEN   TRAN
BLANK card ending  OUTPUT  variable requests
 193 1. 0.0 7.0         GEN   TRAN               { Axis limits: (-2.544,  3.289)
BLANK card ending  PLOT  cards
BEGIN NEW DATA CASE
C          4th of 8 subcases continues illustration of enhanced  $INSERT.  This
C          illustrates  /LOAD FLOW  within  $INSERT  using data from the 5th
C          of 5 data subcases of DC-26.  This is the former  DC17A.DAT
$INSERT, DC4AINS.DAT, { Move all data except misc. data into this file
 .000200    .000     60.   { T-max = 0  means that no transient solution follows
       1       1       1       0       1     
BLANK card ending branch cards.
BLANK card ending switch cards.
BLANK card terminating EMTP source cards.
BLANK card ending requests for output variables
BLANK card ending plot cards
BEGIN NEW DATA CASE
C          5th of 8 subcases continues illustration of enhanced  $INSERT.  This
C          illustrates  /TACS  but no  /MODELS.  The configuration is taken
C          from the 1st subcase of  DC-33.  This is former  DC17B.DAT
$INSERT, DC4BINS.DAT, { All data except misc. data has been moved to here
 .000050    .020
       1       1       1       1       1      -1
       5       5      20      20
BLANK card terminates all TACS data
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
C       Total network loss  P-loss  by summing injections =   7.106450000000E-06
C Inject:  GEN       377.             377.             .377E-7   2.6786345332877
C Inject:            0.0              0.0      -2.678634533288       -89.9999992
C           ---- Initial flux of coil  "GROUND"  to  "GEN   "  =  1.00000000E-09
C
C  Step     Time        GEN          GROUND       GROUND       TACS         TACS
C                                     TERRA        GEN          FLUX         GEN
C ***      Phasor I(0) =  0.3770000E-07   Switch  "GROUND"  to  "      "  closed
C   0      0.0          377.      .377E-7          0.0          0.0          0.0
C   1    .5E-4   376.9330268  .5235646391  -.523564601  .0094233257  376.9330268
C   2    .1E-3   376.7321312   1.04694322  -1.04694318  .0282649546  376.7321312
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE
C          6th of 8 subcases continues illustration of enhanced  $INSERT.  This
C          illustrates  /MODELS  but no  /TACS.  The configuration is taken
C          from the 2nd of 3 subcases of DC-33.  This is former  DC17C.DAT
$INSERT, DC4CINS.DAT, { All data except misc. data has been moved to here
 .000050    .020
       1       1       1       1       1      -1
       5       5      20      20
BLANK card terminates all MODELS data  { 9 July 1995,  MODELS  replaces old
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
C       Total network loss  P-loss  by summing injections =   7.106450000000E-06
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
BEGIN NEW DATA CASE                                                             
C          7th of 8 subcases continues illustration of enhanced  $INSERT.  This
C          illustrates both  /MODELS  and  /TACS.  The configuration is taken
C          from the 6th of 6 subcases of DC-30.  This is former  DC17D.DAT
$INSERT, DC4DINS.DAT,
   1.E-6   5.E-6         { Only take 5 steps;  these are plenty for illustration
       1      -1
BLANK card terminates all TACS data
BLANK card ends all MODELS data
ENDMODELS
BLANK card ending all BRANCH cards
BLANK card ending all SWITCH cards
BLANK terminates the last SOURCE card
BLANK card ends OUTPUT variable requests
BLANK card ending all batch-mode PLOT cards
C   Comment about 6 lines above.  Note  ENDMODELS  has like the initiation
C   word MODELS:  there is only one of them,  and has nothing to do with
C   the actual modeling.  It is part of the structure in which actual data
C   is carried.  Normally,  ENDMODELS  would precede the blank card ending
C   MODELS,  but with sorting that is not possible.  /-card sorting will
C   put the actual data there.  If we raise  ENDMODELS  line by one row,
C   it would precede all real  MODELS  data,  and that would be completely
C   wrong.  So,  we put it after the blank.  Recall  MODELS  ignores blank
C   lines,  so the blank card ending MODELS is ignored.  It is essential
C   to the sorting,  but then is ignored by MODELS itself.  Since it was
C   optional, anyway, this works well.
BEGIN NEW DATA CASE
C          8th of 8 subcases continues illustration of enhanced  $INSERT.  This
C          illustrates nesting: a $INSERT line within a $INSERT file.  The
C          use is quite artificial.  WSM simply took a single-level example
C          from Orlando Hevia and modified it to illustrate the 2nd level.
C          Use is realistic in that  $INSERT  is required.  $INCLUDE  can
C          not be used in its place because a PCVP loop defines one of the
C          arguments.  In fact,  this was the data that inspired WSM to add
C          the service of arguments to $INSERT capability during January, 2006.
C          Arguments came first.  "/"-card sorting by class came later.
C          Misc. data parameter  ICAT = 1  will create one .PL4 file for each
C          pass of the loop.  But because this is not the first subcase within
C          the disk file,  they will not be named  DC4.xxx  as might be hoped.
C          Instead,  they will be renamed using the date and time.  Remember
C          that  KOMPAR = 4  in startup results in the WW I Armistice name of
C          8B111100.  So,  there will be a  .001  file for the first pass,  a
C          .002  file for the second,  etc. through  MAXKNT  passes.
C   Controls of following request card:   MAXKNT  IOPCVP  NOSTAT { Loop 3 times with
POCKET CALCULATOR VARIES PARAMETERS            3       0       1 { full printout
$PARAMETER
RHORHO = 10.000*KNT
BLANK card ends parameter block.
$DUMMY, XYZ000
   1.E-6    .001                
       1       1       1       1       1      -1       0       1
       5       5      20      20     100     100
  XX0001XX0003              1.E3                                               0
  XX0001                  RHORHO
$INSERT, dc4drtt.dat, XX0003, 10.0, RHORHO
BLANK card ending BRANCH cards
BLANK card ending SWITCH cards
15XX0001 0      1.E4    1.2E-6     5.E-5        2.                            1.
BLANK card ending SOURCE cards
  XX0003
BLANK card ending node voltage OUTPUTs
BLANK card ending PLOT cards
BLANK card to be read and ignored by MNTINC.  This avoids 2-line EOF warning.
BEGIN NEW DATA CASE
BLANK