From b18347ffc9db9641e215995edea1c04c363b2bdf Mon Sep 17 00:00:00 2001 From: Angelo Rossi Date: Wed, 21 Jun 2023 12:04:16 +0000 Subject: Initial commit. --- benchmarks/dc33.dat | 711 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 711 insertions(+) create mode 100644 benchmarks/dc33.dat (limited to 'benchmarks/dc33.dat') diff --git a/benchmarks/dc33.dat b/benchmarks/dc33.dat new file mode 100644 index 0000000..4a2cb80 --- /dev/null +++ b/benchmarks/dc33.dat @@ -0,0 +1,711 @@ +BEGIN NEW DATA CASE +C BENCHMARK DC-33 +C Type-96 hysteretic inductor with TACS logic to give flux by integrating v +C Batch-mode EMTP CalComp plotting, including "X-Y PLOT", is illustrated. +C For split version with TACS postponed until postprocessing, see DC-45 and +C DC-46. As for data ordering, this is intentionally shuffled from what +C would be natural in order to illustrate the use of "/" cards for sorting. +C Initial condition & load flow classes were added during November of 1991 +C after Prof. Juan Martinez complained from Barcelona, Spain, that such +C features were not working. He was right. One year later, he complained +C again (it still did not work right. Somehow, work was lost; DC-33 had +C not been changed; the optimistic report of the January, 1992, newsletter +C was incorrect). So, sorting by class was repaired again on 28 Novemb 92. +C 1st of 9 stacked subcases total. + .000050 .020 + 1 1 1 1 1 -1 + 5 5 20 20 +TACS HYBRID +/SOURCE +14GEN 377. 60. -1. +/TACS + 1FLUX +GEN + 1.0 + 0.0 1.0 +90GEN +91GROUND +99CURR -1.0* GROUND +33FLUX GEN GROUNDCURR +/OUTPUT + GEN +/PLOT +C Last step: 400 .02 116.4994069 2.765761424 -2.76576141 .9416257618 +C Last step continued .....: 116.4994069 2.765761424 -2.76576142 +C Final 5 max. : 3.499909256 .9905492104 376.9925583 3.499613141 3.499909256 +C Associated times : .0125 .00415 .01665 .00415 .0125 + 2Arbitrary 78-character case title text of which this is an example, I hope. + First of two lines of 78-byte graph subheading text. + Second and final such line of graph subheading text. + 194 2. 0.0 20. TACS CURR ABCDEFGHIJKLMNOP1234567890123456 + X-Y PLOT Horizontal Axis label123 + 10. -1.1 1.1 + 8. -4.0 4.0 + 194 4. 0.0 20. TACS FLUX TACS CURR Graph heading---Vertical axis la + X-Y PLOT Horizontal Axis label123 + 9999. -1.1 1.1 "9999." returns to conventional (vs. time) plotting + 8. -4.0 4.0 + PRINTER PLOT + 194 4. 0.0 20. TACS CURR { Axis limits: (-3.500, 3.500) +C Now that all plotting is done with, let's illustrate the three declarations +C that will choose the destination of any subsequent "CALCOMP PLOT" use: + SCREEN PLOT { If a vector plot were to follow, it would go only to the screen + PEN PLOT { If a vector plot were to follow, ... go only to CalComp plotter + SCREEN PEN { If a vector plot were to follow, .. go to both screen & plotter +C The preceding 3 declarations really could only be tested by a user who +C had both a vector-graphic screen and a CalComp plotter. Hence no use. +/LOAD FLOW +C Comment of /LOAD FLOW. Since DC-33 does not really involve the load flow, +C Comment of /LOAD FLOW. we can not use real load flow data cards. In their +C Comment of /LOAD FLOW. place, we use these three comment cards. End test +/BRANCH +96GROUNDGEN 8888. 1.E-9 1 + 1.0 -0.7 + 2.0 0.9 + 3.5 1.0 + 9999. +/REQUEST +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) +/INITIAL + 2GROUND 0. { Node voltage initial condition in fact changes nothing +C Initial conditions really are not a part of this problem, although we do want +C to illustrate that they, too, can be sorted, and inserted after /LOADFLOW +C The preceding redefinition of node voltage at GROUND changes zero to zero. +/SWITCH + GROUND MEASURING 1 +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 BENCHMARK DC-33 +C 2nd of 9 subcases is the same as the first except that control system +C modeling is to be done using MODELS rather than the older TACS. Jerry +C Almos of BPA converted the data on 18 Nov 1991 after WSM changed the code +C to allow such sorting using the same /TACS declaration. + .000050 .020 + 1 1 1 1 1 -1 + 5 5 20 20 +MODELS { Note the change; the 1st subcase had "TACS HYBRID" here +/SOURCE +14GEN 377. 60. -1. +/MODELS { 9 July 1995, the former /TACS was changed to this new class +INPUT GEN {V(GEN)}, GROUND {I(GROUND)} +MODEL DC33 +INPUT gen, ground +VAR flux, curr +HISTORY gen {DFLT: 0}, flux {DFLT: 0}, ground {DFLT: 0} +INIT + curr:=0 +ENDINIT +EXEC + COMBINE AS INTEGRATOR + LAPLACE(flux/gen):=(1.0|S0)/(1.0|S1) + ENDCOMBINE + curr:=-ground +ENDEXEC +ENDMODEL +USE DC33 AS DC33 + INPUT GEN:=GEN, GROUND:=GROUND +ENDUSE +RECORD + DC33.FLUX AS FLUX + DC33.GEN AS GEN + DC33.GROUND AS GROUND + DC33.CURR AS CURR +ENDRECORD +ENDMODELS +/OUTPUT + GEN +/PLOT +C Last step: 400 .02 116.4994069 2.765761424 -2.76576141 .9604757618 +C Variable maxima : 377. 3.499613141 3.499909256 1.00939921 +C Times of maxima : 0.0 .00415 .0125 .00415 +C Variable minima : -376.992558 -3.49990926 -3.49961314 -.990568949 +C Times of minima : .00835 .0125 .00415 .0125 + 2Arbitrary 78-character case title text of which this is an example, I hope. + First of two lines of 78-byte graph subheading text. + Second and final such line of graph subheading text. + CALCOMP PLOT + SCREEN PLOT + 194 2. 0.0 20. MODELSCURR ABCDEFGHIJKLMNOP1234567890123456 + X-Y PLOT Horizontal Axis label123 + 10. -1.1 1.1 + 8. -4.0 4.0 + 194 4. 0.0 20. MODELSFLUX MODELSCURR Graph heading---Vertical axis la + X-Y PLOT Horizontal Axis label123 + 9999. -1.1 1.1 "9999." returns to conventional (vs. time) plotting + 8. -4.0 4.0 +/BRANCH +96GROUNDGEN 8888. 1.E-9 1 + 1.0 -0.7 + 2.0 0.9 + 3.5 1.0 + 9999. +/REQUEST +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) +/SWITCH + GROUND MEASURING 1 +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 +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 -.25696E-16 .009425 377. .377E-7 -.377E-7 +C 1 .5E-4 376.9330268 -1.42571976 1.425719797 .0282733257 376.9330268 -1.42571976 1.425719759 +C 2 .1E-3 376.7321312 -1.41394374 1.413943779 .0471149546 376.7321312 -1.41394374 1.413943741 +C 3 .15E-3 376.3973844 -1.40217609 1.40217613 .0659431925 376.3973844 -1.40217609 1.402176093 +BLANK card ends OUTPUT variable requests +BLANK card ending all batch-mode PLOT cards +BEGIN NEW DATA CASE +C 3rd of 9 subcases ends the original series. 4-9 added 5 March 2007 +C Illustration of user-supplied residual flux of Type-96 hysteretic +C inductor was supplied by Prof. Mustafa Kizilcay of FH Osnabrueck in +C Germany (see E-mail dated 23 March 2001). This data is added to show +C that non-zero (but very small) flux of the phasor solution now is +C tolerated. Prior to the change several days ago, it was not. Prior +C to the change, residual flux needed to be exactly zero for the user- +C supplied flux to be used. No longer. Phasor flux < EPSILN is all that +C now is required. +POWER FREQUENCY 50. +PRINTED NUMBER WIDTH, 13, 2, { dT loop columns have width 13 including 2 blanks +ZERO FLUX TOLERANCE 1.E-8 { Columns 33-40 carry Type-96 tolerance +C Note about preceding declaration, which was added 29 March 2001. After +C reading later discussion about trouble with the IDEAL TRANSFORMER data, +C Prof. Kizilcay agreed that definition of the flux tolerance independent of +C EPSILN would be a good idea. The user can define whatever value he wants +C although here nothing will change (1.E-8 = value of EPSILN in STARTUP). + .000200 .200 50. + 40 1 1 1 1 +C Following TACS data is not required. It is used only to integrate the +C inductor voltage, thereby indicating the change in inductor flux --- +C change with respect to the user-supplied initial value of 0.2 volt-sec. +TACS HYBRID + 1FLUSSA +MA 1. + 1. { Numerator is : 1 + 1. { Denominator is : 0 + 1 * s +90MA { Type-96 voltage of electric network defines this TACS source 1. +33FLUSSA { Output of integrator, the change of flux, is only TACS output +BLANK card ending TACS data + QA P1A .0172 .1722 0 + P3A .0172 .1722 0 + MA 1.8 0 +96MA 8888 0.20 1 + -0.1199 -0.3754 { Preceding 0.2 is residual flux in [volt-sec] + -0.06264 -0.3643 { 2nd of many current-flux pairs + -0.0353 -0.3536 { 3rd of many. These define hysteresis loop. + -0.02025 -0.3428 + -0.00732 -0.3214 + -0.002869 -0.3 + -0.000538 -0.2786 + 0.01345 0.0 + 0.02744 0.2786 + 0.02977 0.3 + 0.03422 0.3214 + 0.04715 0.3428 + 0.0622 0.3536 + 0.08954 0.3643 + 0.1468 0.3754 + 0.2254 0.3826 + 0.5494 0.3968 + 9999 { Terminate last point of Type-96 hysteresis element + P2A 1.E10 +BLANK card ending all branch cards + P2A P3A .020 1.0 { Switch closes at 20 msec } 1 +BLANK card ending switches +14QA 8165. 50. -90. -1. +11P1A 1E-10 { 1st half of ungrounded voltage source and ideal Xformer +18P2A .0025MA { 2nd half of ungrounded voltage source and ideal Xformer +C Note about preceding two source cards, which represent an ideal transformer. +C Yes, there is the ungrounded voltage source, but for engineering purposes, +C this has value zero (the 1.E-10 in columns 11-20 of the 1st source card). +C The second source card, a Type-18 source, defines turns ratio and 4 names. +C In theory, these two source cards could be replaced by the following three +C branch cards, which are more readable: +C IDEAL TRANSFORMER +C 1P1A P2A 400.0 { Rated voltage of high side +C 2MA 1.0 { Rated voltage of low side +C The resulting phasor solution would be nearly the same, too. Nearly, +C but not quite. Prof. Kizilcay's Type-11 source value is small, note +C --- small enough to make the phasor flux at time zero < EPSILN. If +C instead the IDEAL TRANSFORMER is used, it would appear that the EPSILN +C threshold is exceeded. In place of the user-desired residual flux of 0.2 +C volt-seconds, the .LIS file will show the noise from the phasor solution: +C ---- Initial flux of coil "MA " to " " = -1.87122597E-06 +C This then is carried into the dT loop. Without meaningful residual flux, +C there is no significant in-rush current. Instead of Prof. Kizilcay's +C 30 amps, the PRINTER PLOT shows only about 1.243E-2 amps. For engineering +C purposes, the simulation will be completely different (and wrong). +BLANK card ending sources +C MA .1797146122E-11 .36742498645E-8 -.998414547E-12 .9984146668E-12 .8600609455E-40 .8600609455E-40 +C -.3674249425E-8 -89.9719755 -.488343774E-15 -179.9719755 .1834212477E-20 0.0000000 +C Total network loss P-loss by summing injections = 3.333364056743E-03 +C ---- Initial flux of coil "MA " to " " = 2.00000000E-01 +C The preceding result of phasor solution confirms correct operation. The +C voltage across the Type-96 branch is not zero, but it is small enough +C compared with floating-point miscellaneous data parameter EPSILN so no +C special precautions must be taken. Understand what is happening. The +C initial flux is equal to the imaginary part of the phasor voltage divided +C by the frequency, and this is small: .367E-8 / OMEGA < EPSILN + QA MA { Names of nodes for node-voltage output +C First 2 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 2 output variables are branch currents (flowing from the upper node to the lower node); +C Next 1 output variables belong to TACS (with "TACS" an internally-added upper name of pair). +C Step Time QA MA P2A MA TACS +C P3A TERRA FLUSSA +C 0 0.0 .499946E-12 .179715E-11 0.0 -.99841E-12 0.0 +C 40 .008 4799.266585 .2154661E-8 0.0 .278303E-11 .211415E-10 +C 80 .016 -7765.37646 -.349521E-8 0.0 .438944E-12 .803576E-11 +C *** Close switch "P2A " to "P3A " after 2.02000000E-02 sec. +C 120 .024 7765.376456 70.13133651 15591.82659 .0177128714 .0990291974 +C 160 .032 -4799.26658 110.5331176 33643.29393 22.70094726 1.167639417 +BLANK card ending node voltage outputs +C 1000 0.2 .9430083E-8 -105.223163 -23383.3158 -.97653E-3 .0484350031 +C Variable maxima : 8165. 159.787186 40889.62376 30.82522183 1.523703303 +C Times of maxima : .005 .0288 .0294 .054 .054 +C Variable minima : -8165. -118.702745 -23490.9996 -.16538292 -.198576963 +C Times of minima : .015 .0792 .1996 .1794 .1846 + CALCOMP PLOT + 193.02 0.0 0.2 -10. 30.MA { Plot Type-96 current over entire time span + PRINTER PLOT + 194 4. 20. 60. MA { Axis limits: (0.000 3.083) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 4th of 9 subcases begins a series using /UM on DC-35 U.M. data +C 1st of 6 subcases that have same solution but slightly different data. +C Begin with /SOURCE (not /UM for U.M.) and full Class-1 U.M. data. +C Output has been reduced from 2 lines per dT to one which consists of +C just the U.M. variables. For the 10 dT that are taken, signals agree +C exactly with DC37.LIS and are documented on comment cards of /OUTPUT + .000400 .004 { Take only 10 time steps -- enough to begin motor rotation + 5 -1 { No plotting (which has nothing special to do with /UM use) +/OUTPUT +C The following shows dT-loop output for each of the subcases of the /UM set: +C Step Time UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 +C TQGEN OMEGM THETAM IPA IPB IPC IE1 IE2 IE3 +C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 +C 5 .002 -.6898285 0.0 0.0 -150.98163 16.0278591 134.953771 15.5970274 131.9385 -147.53553 +C 10 .004 -11.779008 -.21882419 -.42361E-3 -203.70113 -91.734138 295.43527 -90.126277 287.718333 -197.59206 +/SOURCE +14BUS-A0 180.0 60.0 0.0 +14BUS-B0 180.0 60.0 -120.0 +14BUS-C0 180.0 60.0 +120.0 +11BUS-OM-1 -6.00 +/BRANCH + 0 BUS-A0 1.0 + 0 BUS-B0 1.0 + 0 BUS-C0 1.0 + 0BUS-A1BUS-A0 0.005 1.0 + 0BUS-B1BUS-B0BUS-A1BUS-A0 + 0BUS-C1BUS-C0BUS-A1BUS-A0 + 0BUSRA1 1.E-5 + 0BUSRB1 1.E-5 + 0BUSRC1 1.E-5 + 0BUS-OM 6.00E4 + 0 BUS-OM 33.333 +/REQUEST +PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 8 output columns) +/SOURCE +19 { Data for U.M. begins with source type 19 keyed in columns 1-2 + 0 { The zero does nothing. In fact, a blank card would suffice for DC-35 +BLANK card ending U.M. header (3rd of 3 non-comment cards of Class-1 U.M. data) + 4 1 1111BUS-OM 2 0.06 0.18850 { 1st machine table +0.0 0.020291 0 +0.0 0.020291 0 +0.063 0.0001911 BUS-A1 1 0.0 { 1st winding card +0.063 0.0003925 BUS-B1 1 0.0 +0.063 0.0003925 BUS-C1 1 0.0 +0.083 0.0003925 BUSRB1 1 0.0 +0.083 0.0003925 BUSRC1 1 0.0 +0.083 0.0003925 BUSRA1 1 0.0 { 6th and final winding card +/PLOT + 2Arbitrary 78-character case title text of which this is an example, I hope. + First of two lines of 78-byte graph subheading text. + Second and final such line of graph subheading text. +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 +BEGIN NEW DATA CASE +C 5th of 9 subcases continues series using /UM on DC-35 U.M. data +C 2nd of 6 examples that have the same solution but slightly different data +C Replace /SOURCE by /UM and omit all 3 non-comment header cards. But +C then U.M. HEADER CARD (UMHC) must be added to define Class-1 data. + .000400 .004 { Take only 10 time steps -- enough to begin motor rotation + 5 -1 { No plotting (which has nothing special to do with /UM use) +/OUTPUT +C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM { Shorten output by eliminating node voltages +/SOURCE +14BUS-A0 180.0 60.0 0.0 +14BUS-B0 180.0 60.0 -120.0 +14BUS-C0 180.0 60.0 +120.0 +11BUS-OM-1 -6.00 +/BRANCH + 0 BUS-A0 1.0 + 0 BUS-B0 1.0 + 0 BUS-C0 1.0 + 0BUS-A1BUS-A0 0.005 1.0 + 0BUS-B1BUS-B0BUS-A1BUS-A0 + 0BUS-C1BUS-C0BUS-A1BUS-A0 + 0BUSRA1 1.E-5 + 0BUSRB1 1.E-5 + 0BUSRC1 1.E-5 + 0BUS-OM 6.00E4 + 0 BUS-OM 33.333 +/REQUEST +PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 8 output columns) +/UM +C Note that there is no Class-1 U.M. data; ATP will supply the Type-19 card. + 4 1 1111BUS-OM 2 0.06 0.18850 { 1st machine table +0.0 0.020291 0 +0.0 0.020291 0 +0.063 0.0001911 BUS-A1 1 0.0 { 1st winding card +0.063 0.0003925 BUS-B1 1 0.0 +0.063 0.0003925 BUS-C1 1 0.0 +0.083 0.0003925 BUSRB1 1 0.0 +0.083 0.0003925 BUSRC1 1 0.0 +0.083 0.0003925 BUSRA1 1 0.0 { 6th and final winding card +C 14BUSRA1 180.0 60.0 0.0 +/PLOT + PRINTER PLOT { Batch-mode plotting is for the line printer not vector-graphic +/REQUEST +U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1 + 0 { In fact, a blank card would suffice for DC-35 +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 +BEGIN NEW DATA CASE --- NEWSORT --- +C 6th of 9 subcases continues series using /UM on DC-35 U.M. data +C 3rd of 6 examples that have the same solution but slightly different data +C Return to use of /SOURCE and full Class-1 U.M. data as in 1st example. +C But add U.M. HEADER CARD (UMHC) which overrides the normal Class-1 data. +C See explanation of the preceding NEWSORT within DCNEW-25. There +C is old (normal) data sorting logic and there is new data sorting +C logic as summarized in the January, 1999, Can/Am newsletter. Putting +C NEWSORT on the case-separator card as shown has the effect of placing +C a minus sign on SZBED of STARTUP, and this selects the newer sorting +C alternative. Either should handle /UM sorting. The NEWSORT could +C be added to each subcase without significant change to the .LIS file. +C Text of BLANK cards will vary a little, but that should be all. + .000400 .004 { Take only 10 time steps -- enough to begin motor rotation + 5 -1 { No plotting (which has nothing special to do with /UM use) +/OUTPUT +C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM { Shorten output by eliminating node voltages +/SOURCE +14BUS-A0 180.0 60.0 0.0 +14BUS-B0 180.0 60.0 -120.0 +14BUS-C0 180.0 60.0 +120.0 +11BUS-OM-1 -6.00 +/BRANCH + 0 BUS-A0 1.0 + 0 BUS-B0 1.0 + 0 BUS-C0 1.0 + 0BUS-A1BUS-A0 0.005 1.0 + 0BUS-B1BUS-B0BUS-A1BUS-A0 + 0BUS-C1BUS-C0BUS-A1BUS-A0 + 0BUSRA1 1.E-5 + 0BUSRB1 1.E-5 + 0BUSRC1 1.E-5 + 0BUS-OM 6.00E4 + 0 BUS-OM 33.333 +/REQUEST +PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 8 output columns) +U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1 + 0 { In fact, a blank card would suffice for DC-35 +/SOURCE +C Any number of comments should be legal anywhere within Class-1 U.M. data +C as now begins. Because of UMHC (U.M. HEADER CARD) use, only the first of +C the 3 non-comment lines of Class-1 U.M. data ever will be seen. The crucial +C 2nd will be replaced by the card following the UMHC declaration, and the +C blank 3rd card is read and discarded, so will not been seen in the .LIS +C file. Yet, it is required. Even though not-used, the 3rd is checked +C for being blank, and execution will be halted if it is not so. Add "XXX" +C at the start of the 2nd to emphasize that content is arbitrary: +19 { Data for U.M. begins with source type 19 keyed in columns 1-2 +XXX 0 { 2nd of 3 non-comment header lines for U.M. will be replaced by UMHC crd +BLANK card ending U.M. header (3rd of 3 never will been seen in .LIS file) + 4 1 1111BUS-OM 2 0.06 0.18850 { 1st machine table +0.0 0.020291 0 +0.0 0.020291 0 +0.063 0.0001911 BUS-A1 1 0.0 { 1st winding card +0.063 0.0003925 BUS-B1 1 0.0 +0.063 0.0003925 BUS-C1 1 0.0 +0.083 0.0003925 BUSRB1 1 0.0 +0.083 0.0003925 BUSRC1 1 0.0 +0.083 0.0003925 BUSRA1 1 0.0 { 6th and final winding card +/PLOT + 2Arbitrary 78-character case title text of which this is an example, I hope. + First of two lines of 78-byte graph subheading text. + Second and final such line of graph subheading text. +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 +BEGIN NEW DATA CASE +C 7th of 9 subcases continues series using /UM on DC-35 U.M. data +C 4th of 6 examples that have the same solution but slightly different data +C Data of the preceding 3rd example has been massaged. First, "/" cards +C for sorting have been removed to demonstrate that UMHC use does not +C require sorting. The full 3-card Class-1 U.M. data has been replaced +C by the special Type-19 summary line which will refer to UMHC for the +C missing 2nd nonblank card of Class-1 data. Finally, since no sorting +C to insert the blank card ending U.M. data, this must be added manually. +PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 8 output columns) +U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1 + 0 { In fact, a blank card would suffice for DC-35 + .000400 .004 { Take only 10 time steps -- enough to begin motor rotation + 5 -1 { No plotting (which has nothing special to do with /UM use) + 0 BUS-A0 1.0 + 0 BUS-B0 1.0 + 0 BUS-C0 1.0 + 0BUS-A1BUS-A0 0.005 1.0 + 0BUS-B1BUS-B0BUS-A1BUS-A0 + 0BUS-C1BUS-C0BUS-A1BUS-A0 + 0BUSRA1 1.E-5 + 0BUSRB1 1.E-5 + 0BUSRC1 1.E-5 + 0BUS-OM 6.00E4 + 0 BUS-OM 33.333 +BLANK card ending all BRANCH cards +BLANK card ending all SWITCH cards +14BUS-A0 180.0 60.0 0.0 +14BUS-B0 180.0 60.0 -120.0 +14BUS-C0 180.0 60.0 +120.0 +11BUS-OM-1 -6.00 +19 { Class-1 U.M. data has been omitted WSM + 4 1 1111BUS-OM 2 0.06 0.18850 { 1st machine table +0.0 0.020291 0 +0.0 0.020291 0 +0.063 0.0001911 BUS-A1 1 0.0 { 1st winding card +0.063 0.0003925 BUS-B1 1 0.0 +0.063 0.0003925 BUS-C1 1 0.0 +0.083 0.0003925 BUSRB1 1 0.0 +0.083 0.0003925 BUSRC1 1 0.0 +0.083 0.0003925 BUSRA1 1 0.0 { 6th and final winding card +BLANK card ends U.M. data +BLANK terminates the last SOURCE card +C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM { Shorten output by eliminating node voltages +BLANK card ends OUTPUT variable requests + 2Arbitrary 78-character case title text of which this is an example, I hope. + First of two lines of 78-byte graph subheading text. + Second and final such line of graph subheading text. +BLANK card ending all batch-mode PLOT cards +BEGIN NEW DATA CASE +C 8th of 9 subcases continues series using /UM on DC-35 U.M. data +C 5th of 6 subcases that have the same solution but slightly different data +C Data is like the 2nd example except that a second and third disconnected +C machine has been added to illustrate multi-machine use for 3 machines. + .000400 .004 { Take only 10 time steps -- enough to begin motor rotation + 5 -1 { No plotting (which has nothing special to do with /UM use) +/OUTPUT +C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM { Shorten output by eliminating node voltages +/SOURCE +14BUS-A0 180.0 60.0 0.0 +14BUS-B0 180.0 60.0 -120.0 +14BUS-C0 180.0 60.0 +120.0 +11BUS-OM-1 -6.00 +/BRANCH + 0 BUS-A0 1.0 + 0 BUS-B0 1.0 + 0 BUS-C0 1.0 + 0BUS-A1BUS-A0 0.005 1.0 + 0BUS-B1BUS-B0BUS-A1BUS-A0 + 0BUS-C1BUS-C0BUS-A1BUS-A0 + 0BUSRA1 1.E-5 + 0BUSRB1 1.E-5 + 0BUSRC1 1.E-5 + 0BUS-OM 6.00E4 + 0 BUS-OM 33.333 +/REQUEST +U.M. HEADER CARD { To be followed by significant middle card of U.M. Class 1 + 0 { In fact, a blank card would suffice for DC-35 +PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 8 output columns) +/SOURCE +19 { Class-1 U.M. data has been omitted (this card is needed because of /SOURCE) + 4 1 1111BUS-OM 2 0.06 0.18850 { 1st machine table +0.0 0.020291 0 +0.0 0.020291 0 +0.063 0.0001911 BUS-A1 1 0.0 { 1st winding card +0.063 0.0003925 BUS-B1 1 0.0 +0.063 0.0003925 BUS-C1 1 0.0 +0.083 0.0003925 BUSRB1 1 0.0 +0.083 0.0003925 BUSRC1 1 0.0 +0.083 0.0003925 BUSRA1 1 0.0 { 6th and final winding card +C 14BUSRA1 180.0 60.0 0.0 +/PLOT + PRINTER PLOT { Batch-mode plotting is for the line printer not vector-graphic +C ---------------- Done with data of 1st (real) U.M. Begin data of 2nd: +C 2nd, unused U.M. has the same data as the 1st, but it is disconnected, +C so should not affect the solution. Outputs have been suppressed. For +C variety, this U.M. still uses /UM whereas the first has been converted +C to /SOURCE followed by the special Type-19 header card. Note that the +C Type-19 card is mandatory because any source might follow /SOURCE whereas +C only U.M. data can follow /UM (i.e., the "19" is needed for recognition +C as U.M. data). For the 2nd, disconnected machine, note that the original +C A3 text BUS has become GUS everywhere. To restore the 9 outputs of +C the 2nd machine, add "111" to cols. 7-9 of the card immediately following +C the /UM declaaration, and change "0" to "1" in column 47 of each of +C the 6 winding cards that end machine data. If this is done, output consists +C of 2 lines each time step, and the two will be identical. +/UM + 4 1 1 GUS-OM 2 0.06 0.18850 { 1st machine table +0.0 0.020291 0 +0.0 0.020291 0 +0.063 0.0001911 GUS-A1 0 0.0 { 1st winding card +0.063 0.0003925 GUS-B1 0 0.0 +0.063 0.0003925 GUS-C1 0 0.0 +0.083 0.0003925 GUSRB1 0 0.0 +0.083 0.0003925 GUSRC1 0 0.0 +0.083 0.0003925 GUSRA1 0 0.0 { 6th and final winding card +/BRANCH + 0 GUS-A0 1.0 + 0 GUS-B0 1.0 + 0 GUS-C0 1.0 + 0GUS-A1GUS-A0 0.005 1.0 + 0GUS-B1GUS-B0GUS-A1GUS-A0 + 0GUS-C1GUS-C0GUS-A1GUS-A0 + 0GUSRA1 1.E-5 + 0GUSRB1 1.E-5 + 0GUSRC1 1.E-5 + 0GUS-OM 6.00E4 + 0 GUS-OM 33.333 +/SOURCE +14GUS-A0 180.0 60.0 0.0 +14GUS-B0 180.0 60.0 -120.0 +14GUS-C0 180.0 60.0 +120.0 +11GUS-OM-1 -6.00 +C ---------------- Done with data of 2nd U.M. Begin data of 3rd: +C 3rd, unused U.M. has the same data as the 1st, but it is disconnected, +C so should not affect the solution. Outputs have been suppressed. Etc. +C (like the 2nd except that the original BUS ---> MUS here). +/UM + 4 1 1 MUS-OM 2 0.06 0.18850 { 1st machine table +0.0 0.020291 0 +0.0 0.020291 0 +0.063 0.0001911 MUS-A1 0 0.0 { 1st winding card +0.063 0.0003925 MUS-B1 0 0.0 +0.063 0.0003925 MUS-C1 0 0.0 +0.083 0.0003925 MUSRB1 0 0.0 +0.083 0.0003925 MUSRC1 0 0.0 +0.083 0.0003925 MUSRA1 0 0.0 { 6th and final winding card +/BRANCH + 0 MUS-A0 1.0 + 0 MUS-B0 1.0 + 0 MUS-C0 1.0 + 0MUS-A1MUS-A0 0.005 1.0 + 0MUS-B1MUS-B0MUS-A1MUS-A0 + 0MUS-C1MUS-C0MUS-A1MUS-A0 + 0MUSRA1 1.E-5 + 0MUSRB1 1.E-5 + 0MUSRC1 1.E-5 + 0MUS-OM 6.00E4 + 0 MUS-OM 33.333 +/SOURCE +14MUS-A0 180.0 60.0 0.0 +14MUS-B0 180.0 60.0 -120.0 +14MUS-C0 180.0 60.0 +120.0 +11MUS-OM-1 -6.00 +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 +BEGIN NEW DATA CASE +C 9th of 9 subcases ends the series using /UM on DC-35 U.M. data +C 6th of 6 examples that have the same solution but slightly different data +C Data is the same as the 2nd example except that /UM data has been moved +C to a $INCLUDE file. Use or non-use of $INCLUDE should have no effect, +C and this demonstrates none for this simplest, first level of $INCLUDE. +C Enhancement by Gabor Furst. WSM's single $INCLUDE was replaced by +C two for more generality. A second induction motor was added to show +C that this is possible. Mr. Furst observed the following on 9 April +C 2007: "I thought we could have this sub case with two rather than +C one $INCLUDE file. I attach my version of this, which demonstrates +C how adding another machine is simply adding another $INCLUDE file." +C Mr. Furst simulated to Tmax = 1.044 seconds but program developers +C truncate this drastically (to 10 steps) to speed execution. This is +C just enough to start the rotors and show approximate balance. To see +C the entire settling into the steady state, restore the longer T-max, +C remove the minus sign from IPLOT immediately below T-max, and +C activate the batch-mode plot card (now omitted by $DISABLE). Finally, +C control printout (every 5th step for IOUT is a lot, for realistic use). +$PREFIX, [] { $INCLUDE files are located in same place as this main data file + .000400 .004 { Take only 10 time steps -- enough to begin motor rotation + 5 -1 { No plotting (which has nothing special to do with /UM use) +/OUTPUT +C BUS-A1BUS-B1BUS-C1BUS-A0BUS-OM { Shorten output by eliminating node voltages +/SOURCE +14BUS-A0 180.0 60.0 0.0 +14BUS-B0 180.0 60.0 -120.0 +14BUS-C0 180.0 60.0 +120.0 +11BUS-OM-1 -6.00 +14BUS-A3 180.0 60.0 0.0 +14BUS-B3 180.0 60.0 -120.0 +14BUS-C3 180.0 60.0 +120.0 +11BUS-3M-1 -6.00 +/BRANCH +C the following branch records are transferred to $INCLUDE dc33inc1.dat +C 0 BUS-A0 1.0 +C 0 BUS-B0 1.0 +C 0 BUS-C0 1.0 +C 0BUS-A1BUS-A0 0.005 1.0C +C 0BUS-B1BUS-B0BUS-A1BUS-A0 +C 0BUS-C1BUS-C0BUS-A1BUS-A0 +C 0BUSRA1 1.E-5 +C 0BUSRB1 1.E-5 +C 0BUSRC1 1.E-5 +C 0BUS-OM 6.00E4 +C 0 BUS-OM 33.33 +C + 0 BUS-A3 1.0 + 0 BUS-B3 1.0 + 0 BUS-C3 1.0 + 0BUS-A4BUS-A3 0.005 1.0 + 0BUS-B4BUS-B3BUS-A4BUS-A3 + 0BUS-C4BUS-C0BUS-A4BUS-A3 + 0BUSRA4 1.E-5 + 0BUSRB4 1.E-5 + 0BUSRC4 1.E-5 + 0BUS-3M 6.00E4 + 0 BUS-3M 33.333 +$INCLUDE, dc33inc1.dat, { Separate file for /UM followed by all data for U.M. 1 +$INCLUDE, dc33inc2.dat, { Separate file for /UM followed by all data for U.M. 2 +C About the preceding data modularization, there is lack of balance. This +C is intentional. Whereas DC33INC1 has /BRANCH data, DC33INC2 has none. +C Each also involves /REQUEST which precedes a UMHC request. Since only +C one has meaning, expect rejection accompanied by a warning message for the +C 2nd of the two. +BLANK card ending all BRANCH cards +BLANK card ending all SWITCH cards +BLANK terminates the last SOURCE card +BLANK card ends OUTPUT variable requests + CALCOMP PLOT +$DISABLE { The following plot card assumes use of Gabor Furst's longer T-max + 193 .1 0.0 1.0 UM-1 OMEGM UM-2 OMEGM Induction motor w in [rad/sec] +$ENABLE +BLANK card ending all batch-mode PLOT cards +BEGIN NEW DATA CASE +BLANK -- cgit v1.2.3