BEGIN NEW DATA CASE C BENCHMARK DCNEW-16 C 1st of 12 subcases that confirm various aspects of switches that C touch compensation-based elements (bugs removed 17 December 1995). C E-mail from Laurent Dube on Date: Sat, 15 Jul 1995 13:14:55 -0700 (PDT) C "Here is a trimmed-down version of Janko Kosmac's data case showing C wrong switch current values (which in his case were picked up by C MODELS to drive some other logic)." Test here is even more simplified: C C SRCE N1 Type-91 N2 1.0 C o------_------o----/\/\/\----o----/\/\/\----------- C ^ | | | C closed | |<-- switch | C switch | | closed | C |----/\/\/\----o | C 1.E-6 N12 | C ___ C _ 1.E-5 1.E-5 1 1 1 1 TACS HYBRID 99RESIS = 1.E-6 77RESIS 1.E-6 33RESIS BLANK card ending TACS 91N1 N2 TACS RESIS 1 N1 N12 1.E-6 N2 1. BLANK card ending branches N12 N2 MEASURING SRCE N1 MEASURING 1 BLANK card ending switches 14SRCE 100. 60. 0.0 0. -1. BLANK card ending sources C Output for steady-state phasor switch currents. C Node-K Node-M I-real I-imag I-magn Degrees C N12 N2 9.99999000E+01 0.00000000E+00 9.99999000E+01 0.0000 C SRCE N1 9.99999000E+01 0.00000000E+00 9.99999000E+01 0.0000 N1 N12 N2 C Step Time N1 N12 N2 SRCE N1 TACS C N1 N2 RESIS C Phasor I(0) = 9.999990E1 Switch "N12 " to "N2 " closed in steady-st C Phasor I(0) = 9.999990E1 Switch "SRCE " to "N1 " closed in steady-st C 0 0.0 100. 99.9999 99.9999 99.9999 0.0 1.E-6 C 1 .1E-4 99.9992894 99.9992394 99.9992394 99.9992394 49.9996197 1.E-6 BLANK card ending voltage printout BLANK card ending plot BEGIN NEW DATA CASE C 2nd of 12 subcases began as the big, 7th subcase of DC-16. Eventually, C nearly everything could be thrown away. After closing, the switch C voltage had better be zero (one of the improvements of 17 Dec 95). This C was not a problem 2 or more days earlier, however. This problem resulted C after making just the single OVER16 correction. Correction of it C resulted in the two changes of SUBTS3 (see WSM95DEC idents). .000100 .0004 60. 1 -1 1 1 CR20A CR30A 93.40 1 92CR30A CR20A 5555. 3 147500. -1. 1.0 40. .80 9999. RAVBA CR30A 0.5 14.0 GRCBA CR20A 0.4 9.0 1 BLANK card ending program branch cards. CR20A 0.00015 10. { Fault switch, phase "a" to ground } 3 BLANK card terminating program switch cards 14RAVBA 440000. 60. -20.0 -1. 14GRCBA 440000. 60. 0.0 -1. BLANK card terminating program source cards. C Total network loss P-loss by summing injections = 3.505855948763E+08 C Node Source node voltage Injected source current C name Rectangular Polar Rectangular Polar C RAVBA 413464.7531458 440000. 26317.186674075 27911.988458816 C -150488.8630633 -20.0000000 -9299.719634859 -19.4619099 BLANK card ending output variables requests (none here, since all column 80) C First 2 output variables are electric-network voltage differences (upper C Next 4 output variables are branch currents (flowing from the upper nod C Step Time CR30A CR20A CR20A CR30A CR20A GRCBA C CR20A TERRA TERRA CR20A CR30A CR20A C 0 0.0 -264114.27 534224.351 0.0 0.0 -26317.187 -26317.187 C 1 .1E-3 -195321.89 509234.084 0.0 -75571.13 -102186.97 -26615.843 C *** Close switch "CR20A " to " " after 2.00000000E-04 sec. C 2 .2E-3 -126265.42 483522.274 0.0 -.00199387 -26810.516 -26810.514 C 3 .3E-3 -97153.756 0.0 1625.62684 -.10944E-3 -27570.069 -25944.443 C 4 .4E-3 -66918.302 0.0 4833.9446 -.75381E-4 -28909.76 -24075.815 BLANK card ending plot cards BEGIN NEW DATA CASE C 3rd of 12 subcases that confirm various aspects of switches that C touch compensation-based elements. C IMTESTA3.DAT --- Name used by Prof. Juan Martinez of Barcelona, Spain C Based on I. Bonfanti's case : EMTP News, Vol. 2, no. 3, September 1989 C Steady State Initialization C Wye-connected ungrounded armature - Resistor between neutral and ground C Switches permanently closed are connected to the armature terminals C Until correction on 26 Jan 96, TQGEN was completely wrong on step 2. POWER FREQUENCY, 50 1.E-4 .0002 1.E-15 1 1 1 1 C ----- Network description FEMR CBR .37024 1.196 1 FEMS CBS .37024 1.196 1 FEMT CBT .37024 1.196 1 C ----- Grounding resistor CSMT 1000. C ----- Motor parasitic capacitances M1.R .02 1 M1.S .02 1 M1.T .02 1 C ----- Mechanical network ROTORD 1.E-08 ROTORQ 1.E-08 COPPIA 9.16E6 2 COPPI2COPPIA 1.0E-6 BLANK ENDING BRANCHES CBR M1.R -1. 15.00 1 CBS M1.S -1. 15.00 1 CBT M1.T -1. 15.00 1 BLANK ENDING SWITCHES 14FEMR 4899. 50. 0. -1.0 9999. 14FEMS 4899. 50. -120. -1.0 9999. 14FEMT 4899. 50. 120. -1.0 9999. 14COPPI2-1 -1. 1.E-9 0. -1.0 9999. C ----- UM specification 19 UM 01 0 - Compensation; change 0 to 1 if prediction is wanted BLANK 3 1 1111COPPIA 1 0.0 0.3964 0.0 0.3964 4. COPPI2 .1674 .001 M1.R CSMT 1 .1674 .01 M1.S CSMT 1 .1674 .01 M1.T CSMT 1 .7819 .00453 ROTORD 1 .7819 .00453 ROTORQ 1 BLANK card ending U.M. data BLANK card ending sources C Total network loss P-loss by summing injections = 5.005262005074E-01 C Total network loss P-loss by summing injections = 2.721172542328E+04 C Total network loss P-loss by summing injections = 3.083034016021E+04 C Total network loss P-loss by summing injections = 3.209218252624E+04 C Output for steady-state phasor switch currents. C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive C CBR M1.R 2.18879785E+02 -8.91358682E+01 2.36333585E+02 -22.1580 5.25806422E+05 2.07845254E+05 C CBS M1.S -1.86633819E+02 -1.44987520E+02 2.36333585E+02 -142.1580 5.25806422E+05 2.07845254E+05 C CBT M1.T -3.22459663E+01 2.34123388E+02 2.36333585E+02 97.8420 5.25806422E+05 2.07845254E+05 C 1st gen: COPPIA 301.59289474462 301.59289474462 4868.1875 4868.1875 734105.38014229 734105.38014229 C 0.0 0.0 .17357870652E-4 0.0000002 -.0026175052283 1.0000000 C C 2nd gen: FEMR 4899. 4899. 218.87978507236 236.33358480802 536146.03353474 578899.11598725 C 0.0 0.0 -89.13586817262 -22.1579539 218338.30908884 0.9261476 BLANK card ending node voltage output requests (none here) C Step Time COPPIA CBR CBS CBT FEMR FEMS FEMT M1.R M1.S M1.T C TERRA M1.R M1.S M1.T CBR CBS CBT TERRA TERRA TERRA C C 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 C *** Phasor I(0) = 2.1887979E+02 Switch "CBR " to "M1.R " closed in the steady-state. C *** Phasor I(0) = -1.8663382E+02 Switch "CBS " to "M1.S " closed in the steady-state. C *** Phasor I(0) = -3.2245966E+01 Switch "CBT " to "M1.T " closed in the steady-state. C 0 0.0 301.592895 218.879785 -186.63382 -32.245966 218.879785 -186.63382 -32.245966 .309378E-3 .02587952 -.0261889 C -4976.4349 301.592895 1.57079633 -218.87948 186.659698 32.2197774 65.266856 -275.17162 C 1 .1E-3 301.593554 221.571447 -181.98783 -39.583621 221.571447 -181.98783 -39.583621 -.63563E-3 .026350533 -.02571491 C -4976.4333 301.592895 1.60095562 -221.57208 182.014176 39.5579061 64.9214972 -275.25321 C 2 .2E-3 301.594872 224.044463 -177.16225 -46.882216 224.044463 -177.16225 -46.882216 -.00157782 .026791347 -.02521353 C -4976.4315 301.592895 1.63111491 -224.04604 177.189039 46.8570021 64.5760272 -275.33436 BLANK card ending plot cards BEGIN NEW DATA CASE C 4th of 12 subcases that confirm various aspects of switches that C touch compensation-based elements. C Ivano Bonfanti of CESI in Milano, Italy, published this troublesome data C in LEC's EMTP News, Vol. 2, no. 3, September, 1989. Although the answer C here is slightly different than the published "M39." EMTP answer shown by C Bonfanti, either is believed to be correct for engineering purposes. The C machine obviously is in the steady state as it should we. WSM, 28 Jan 96. C Note EPSILN = 1.E-8 a few lines below is needed to restore default value C of STARTUP following distortion of the 3rd subcase. Without restoration, C the answer will be quite different because U.M. uses this to build network C during the 4 phasor solutions. PRINTED NUMBER WIDTH, 10, 1, { 10-digit col width including 1 blank separators POWER FREQUENCY, 50 2.E-4 .020 50. 0.0 1.E-8 1 1 1 1 1 -1 5 5 20 20 100 100 TACS HYBRID 11CARICO -4797. -1. 1. 91COPPI1 -1. 1. 99COPPIA =CARICO+COPPI1 33COPPIA BLANK card ending TACS FEMR CBR .37024 1.196 1 FEMS CBS .37024 1.196 1 FEMT CBT .37024 1.196 1 C ----- Motor parasitic capacitances M1.R .02 1 M1.S .02 1 M1.T .02 1 CSMT 1000. { Motor neutral grounding resistance C ----- Mechanical network ROTORD 1.E-08 ROTORQ 1.E-08 COPPIA 9.16E6 2 COPPI1COPPIA 1.0E-6 BLANK ENDING BRANCHES COPPI1COPPI2 MEASURING 1 CBR M1.R -1. 15.00 1 CBS M1.S -1. 15.00 1 CBT M1.T -1. 15.00 1 BLANK ENDING SWITCHES 14FEMR 1 4898.98 50. 0. -1.0 1. 14FEMS 1 4898.98 50. -120. -1.0 1. 14FEMT 1 4898.98 50. 120. -1.0 1. 14COPPI2-1 -1. 1.E-9 0. -1.0 1. 60COPPIA-1 19 UM 01 0 { Compensation; use "1" rather than "0" for prediction BLANK card ends Class-1 cards 3 1 1111COPPIA 1 0.0 0.3964 0.0 0.3964 4. COPPI2 .1674 .001 M1.R CSMT 1 .1674 .01 M1.S CSMT 1 .1674 .01 M1.T CSMT 1 .78187 .00453 ROTORD 1 .78187 .00453 ROTORQ 1 BLANK card ending U.M. coils BLANK ENDING SOURCES C Total network loss P-loss by summing injections = 5.005262016348E-01 C Total network loss P-loss by summing injections = 3.007635104268E+04 C Total network loss P-loss by summing injections = 3.007585104268E+04 C Total network loss P-loss by summing injections = 3.008749127839E+04 C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive C COPPI1 COPPI2 4.82497291E+03 0.00000000E+00 4.82497291E+03 0.0000 7.27577133E+05 0.00000000E+00 C CBR M1.R 2.12219079E+02 -9.54916543E+01 2.32713544E+02 -24.2262 5.09803228E+05 2.01520807E+05 C CBS M1.S -1.88807738E+02 -1.36041286E+02 2.32713544E+02 -144.2262 5.09803228E+05 2.01520807E+05 C CBT M1.T -2.34113411E+01 2.31532941E+02 2.32713544E+02 95.7738 5.09803228E+05 2.01520807E+05 C 1 of 4 gen: COPPIA 301.59289474462 301.59289474462 4824.9729075506 4824.9729075506 727588.77312628 727588.77312628 C 0.0 0.0 .17357870652E-4 0.0000002 -.0026175052283 1.0000000 C C 2 of 4 gen: FEMR 4898.98 4898.98 212.21907905668 232.71354399216 519828.51195855 570029.49887335 C 0.0 0.0 -95.49165430404 -24.2261891 233905.85230121 0.9119327 C C FEMS -2449.49 4898.98 -188.807738005 232.71354399216 519828.51195855 570029.49887335 C -4242.641132632 -120.0000000 -136.0412864788 -144.2261891 233905.85230121 0.9119327 C C FEMT -2449.49 4898.98 -23.41134105164 232.71354399216 519828.51195855 570029.49887335 C 4242.6411326319 120.0000000 231.53294078285 95.7738109 233905.85230121 0.9119327 C C Step Time COPPIA COPPI1 CBR CBS CBT FEMR FEMS FEMT M1.R M1.S M1.T C TERRA COPPI2 M1.R M1.S M1.T CBR CBS CBT TERRA TERRA TERRA C C TACS UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 C COPPIA TQGEN OMEGM THETAM IPA IPB IPC IE1 IE2 C *** Phasor I(0) = 4.8249729E+03 Switch "COPPI1" to "COPPI2" closed in the steady-state. C *** Phasor I(0) = 2.1221908E+02 Switch "CBR " to "M1.R " closed in the steady-state. C *** Phasor I(0) = -1.8880774E+02 Switch "CBS " to "M1.S " closed in the steady-state. C *** Phasor I(0) = -2.3411341E+01 Switch "CBT " to "M1.T " closed in the steady-state. C 0 0.0 301.59289 4824.9729 212.21908 -188.8077 -23.41134 212.21908 -188.8077 -23.41134 .00137262 .024922 -.0262946 C 0.0 -4824.973 301.59289 1.5707963 -212.2177 188.83266 23.385046 74.005629 -268.4606 C 1 .2E-3 301.59289 4824.9729 217.79488 -179.8952 -37.89965 217.79488 -179.8952 -37.89965 -.4897E-3 .02588688 -.0253972 C 27.972911 -4824.959 301.59289 1.6311149 -217.7954 179.92112 37.874249 73.334597 -268.6438 C 2 .4E-3 301.5932 4824.9729 222.51143 -170.273 -52.23846 222.51143 -170.273 -52.23846 -.002339 .02673945 -.0244004 C 27.972911 -4824.941 301.59289 1.6914335 -222.5138 170.29971 52.21406 72.662955 -268.8252 BLANK ENDING OUTPUT REQUEST C 100 .02 301.64888 4824.9729 211.99353 -188.6458 -23.34774 211.99353 -188.6458 -23.34774 .00136234 .02495568 -.026318 C 27.972911 -4819.913 301.64979 7.6030487 -211.9922 188.67075 23.321418 5.0277011 -278.1347 C Variable maxima : 301.64888 4824.9729 232.67504 232.63158 232.51922 232.67504 232.63158 232.51922 .02962796 .02960618 .02960763 C 27.972911 -4819.913 301.65916 7.6030487 232.6071 232.54547 232.60707 74.005629 -268.4606 C Times of maxima : .02 .0118 .0014 .008 .0146 .0014 .008 .0146 .0152 .002 .0084 C .0118 .02 .0122 .02 .0114 .018 .0046 0.0 0.0 C Variable minima : 301.59289 4824.9729 -232.5957 -232.5347 -232.5969 -232.5957 -232.5347 -232.5969 -.029647 -.0296246 -.0296779 C 0.0 -4824.973 301.59289 1.5707963 -232.6864 -232.6424 -232.5295 5.0277011 -278.1347 C Times of minima : .2E-3 0.0 .0114 .018 .0046 .0114 .018 .0046 .0052 .0118 .0184 C 0.0 0.0 0.0 0.0 .0014 .008 .0146 .02 .02 PRINTER PLOT 194 5. 0.0 20. UM-1 TQGEN COPPI2COPPI1 { Limits: (-4.825, 0.000) 194 5. 0.0 20. TACS COPPIAUM-1 THETAM { Limits: (0.000, 2.797) 194 4. 0.0 20. BRANCH { Limits: (-2.327, 3.017) M1.R CBR M1.S CBS M1.T CBT UM-1 OMEGM 194 4. 0.0 20. BRANCH { Limits: (-2.327, 2.326) UM-1 IPA UM-1 IPB UM-1 IPC BLANK card ending plot cards BEGIN NEW DATA CASE C 5th of 12 subcases that confirm various aspects of switches that C touch compensation-based elements. This is like the 1st except C that here we have 3 phases. As explained in the July, 1996, C newsletter, Marjan Popov in Macedonia first complained of trouble. C Prior to 1 June 1996, step 2 would never be reached because the C error occurred following opening on step 1. Here, we go 3 steps: .000050 .000150 50. 50. 1.E-9 1 1 TACS HYBRID C Note. Mr. Popov used MODELS, and he had complicated dynamics. But C none of that complexity was required to demonstrate the trouble C (KILL = 209). So, instead use simpler TACS model of 1st subcase: 99RESIS = 1.E-6 77RESIS 1.E-6 33RESIS BLANK card ending TACS $DEPOSIT, MATFUL=1 { Internal conversion of Type-51, 52, ... to unsymmetric data C Add preceding definition of MATFUL on 6 August 2008 as the update to allow C unsymmetric branch [R], [L] is being packaged. Of course, the following three C coupled branches are symmetric, so by pretending that they are not we merely C waste both storage (the size of List 3 will increase by 3 cells) and computer C time (the unsymmetric phasor solution requires more work). As for location, C there is no magic. If MATFUL has value unity at the start of a Type-51, 52, C ... branch group, ATP will convert the R-L data to unsymmetric (full matrix) C storage. To turn off this "service," use $DEPOSIT, MATFUL=0 which will be C found near the start of the following subcase (it could be moved to this one, C to any point below the Type-53 branch card, just as well). WSM. 51SRCR BUSR 21.07 { ZERO SEQUENCE 52SRCS BUSS 2.39 { POSITIVE SEQUENCE 53SRCT BUST $VINTAGE, 1 1LINBR LINER 3.98507949E+00 1.16971333E+01 5.00880327E+01 2LINBS LINES 1.64595110E+00 3.84463816E+00 -7.10131829E+00 3.81191239E+00 1.21114617E+01 4.86330090E+01 3LINBT LINET 1.57601253E+00 4.21324332E+00 -8.66201631E+00 1.50186508E+00 4.04356685E+00 -6.51737363E+00 3.69294935E+00 1.24065061E+01 4.91034827E+01 $VINTAGE, 0 91ARCR LINBR TACS RESIS 2 91ARCS LINBS TACS RESIS 2 91ARCT LINBT TACS RESIS 2 C NAME NAME REFERENCE R(OHM)X(OHM)C(miF) BUSR 0.314 { Bus bar capacitances BUSS 0.314 BUST 0.314 C Provide linear elements connected to ARCR, ARCS, ARCT (leakage to ground): ARCR 1.E8 ARCS 1.E8 ARCT 1.E8 BLANK Card ending branch cards BUSR ARCR -0.1 2.60 1 BUSS ARCS -0.1 2.60 1 BUST ARCT -0.1 2.60 1 ARCR LINBR -0.1 0.000 50.E3 ARCS LINBS -0.1 0.001 50.E3 ARCT LINBT -0.1 0.001 50.E3 BLANK Card ending switch cards 14SRCR 89814. 50. 0. 0. -1. 14SRCS 89814. 50. 120. 0. -1. 14SRCT 89814. 50. -120. 0. -1. BLANK Card ending source cards LINER LINES LINET BUSR BUSS BUST BLANK Card ending output cards C Step Time LINER LINES LINET BUSR BUSS BUST BUSR BUSS BUST TACS C ARCR ARCS ARCT RESIS C *** Phasor I(0) = -1.2037996E-01 Switch "BUSR " to "ARCR " closed in the steady-state. C *** Phasor I(0) = -4.2912571E+00 Switch "BUSS " to "ARCS " closed in the steady-state. C *** Phasor I(0) = 4.3267803E+00 Switch "BUST " to "ARCT " closed in the steady-state. C *** Phasor I(0) = -1.2127823E-01 Switch "ARCR " to "LINBR " closed in the steady-state. C *** Phasor I(0) = -4.2908080E+00 Switch "ARCS " to "LINBS " closed in the steady-state. C *** Phasor I(0) = 4.3272295E+00 Switch "ARCT " to "LINBT " closed in the steady-state. C 0 0.0 89846.7243 -44918.879 -44929.061 89826.4331 -44913.159 -44913.549 -.12037996 -4.2912571 4.32678032 1.E-6 C *** Open switch "ARCR " to "LINBR " after 5.00000000E-05 sec. C 1 .5E-4 89835.7272 -46135.534 -43701.414 89815.3391 -46129.508 -43686.133 -.20216547 -4.2510535 4.36857684 1.E-6 C 2 .1E-3 89802.5653 -47340.817 -42462.973 89782.085 -47334.497 -42447.917 -.28390448 -4.2098234 4.40932077 1.E-6 C 3 .15E-3 89747.2461 -48534.443 -41214.032 89726.6777 -48527.808 -41199.225 -.3655754 -4.1675019 4.44892539 1.E-6 BLANK CARD ending plot BEGIN NEW DATA CASE C 6th of 12 subcases illustrates trouble with compensation-based elements C even though no switch is involved. The data came from Gabor Furst on C 14 Sept 2000 as mentioned in the January, 2001, newsletter. Execution C should be stopped in overlay 16 because the Type-91 element R(t) is not C isolated in a separate subnetwork from the U.M. Case-summary statistics C show List 24 = 4 as the number of phases of compensation. This is C entirely wrong: Size 21-30: 0 0 9 4 337 ... C STEP ZERO COUPLE { If activated, the correct KILL = 9 of overlay 15 results C 20 May 2001, the answer changes following correction to SOLVUM to C solve a problem from Dr. Michael Steurer of CAPS at Florida State. C No, answer is not right. It remains wrong, but is a different wrong C because compensation logic has been modified. C 23 May 2001, introduce request for special verification of possible C overlap between U.M. compensation and that of List-9 elements. This is C in OVER16 (to find the code, search for KOMPUM). Add new request: VERIFY U.M. COMPENSATION { If removed, KILL = 9 error termination will disappear $DEPOSIT, MATFUL=0 { Cancel conversion of Type-51, 52, ... to unsymmetric data UM TO TACS 5.0E-4 5.E-4 50.0 { Only 1 time step is needed to demonstrate the error 1 1 TACS { The Rule Book recommends TACS HYBRID or TACS STAND ALONE, but this works 92TQGEN 98TPER = 1/FREQHZ 1TAVG1 +TQGEN 1 1.0 1.0 0.020 1TAVG +TAVG1 1.0 1.0 1.0 0.020 98ROTRES = 1000.0*(1.0 - TIMEX/8.0) 33IM TAVG TQGEN ROTRES BLANK card ending TACS data SRA MOTA 1.00 { source impedance ohm } 1 SRB MOTB 1.00 1 SRC MOTC 1.00 1 C Tops of 3 rotor coils are grounded through 1/2-ohm resistors: XOTA 0.50 { Fixed resistance for rotor phase a } 1 XOTB 0.50 { Fixed resistance for rotor phase b } 1 XOTC 0.50 { Fixed resistance for rotor phase c } 1 C Bottoms of 3 rotor coils are connected to neutral NEUT, which then is C grounded by the following TACS-defined, time-varying resistor: 91NEUT TACS ROTRES { Variable resistance in rotor circuit } 1 C But node NEUT is connected to no linear branch. To avoid a warning C message, parallel the preceding R(t) by fixed high reistance: NEUT 1.E6 C the anlogue network records follow C the separator from the 14 source INERS INER 1.E-6 INER 3.3E7 {inertia uF} C the damping term in ohms INER 2.26 {damping 1/mho} BLANK ending BRANCHes BLANK ending SWITCHes 14SRA 3400.00 50.0 0.0 -1 14SRB 3400.00 50.0 240.0 -1 14SRC 3400.00 50.0 120.0 -1 C the source for the analogue network 14INERS -1 0.000001 .0000001 -1 C the source of any additional load applied 14INERS -1 -0.000001 0.0000001 5.00 19 { Begin U.M. data for single 3-phase induction motor with wound rotor 0 0 { Note no auto-initialization (no multiple phasor solutions) BLANK 4 111INER 2 0.15700 0.465131 0.465131 C armature coils MOTA 1 0.3393571 0.0109180 MOTB 1 0.3393571 0.0109180 MOTC 1 C rotor coils 0.5785651 0.0109180 XOTB NEUT 1 0.5785651 0.0109180 XOTC NEUT 1 XOTA NEUT 1 BLANK card ending U.M. data BLANK card ending sources BLANK card ending output variables requests (none here, since all column 80) BLANK card ending plot cards BEGIN NEW DATA CASE C 7th of 12 subcases illustrates new TACS R-thev option on rotor coil C cards to define connected resistance using a TACS variable. It does C not solve the preceding problem, but it avoids it, and is superior C for cases of common interest. See January, 2000, newsletter. Data is C from preceding subcase, although here we energize the motor. To speed C the simulation, double Gabor Furst's dT and shorten his T-max (was 8): .001 2.0 50.0 1 3 1 0 1 -1 5 5 20 20 100 100 200 200 TACS HYBRID 98ROTRES = 10. * ( 1.0 - TIMEX / 2.0 ) { Linear ramp of R from 10 ohms to zero 33ROTRES { Output this one and only TACS variable of interest 77ROTRES 10. { Initial condition is needed for smooth start (not R=0) BLANK card ending TACS cards SRA MOTA 1.00 { Source impedance = 1 ohm reactive. SRB MOTB 1.00 { This separates the machine armature SRC MOTC 1.00 { MOT from the infinite bus SR. XOTA 100. { Resistance across each rotor coil is XOTB 100. { arbitrary since it will be replaced XOTC 100. { by value of TACS signal ROTRES, anyway C the anlogue network records follow C the separator from the 14 source INERS INER 1.E-6 INER 3.3E7 {inertia uF} C the damping term in ohms INER 2.26 {damping 1/mho} BLANK card ending BRANCH cards BLANK card ending SWITCH cards 14SRA 3400.00 50.0 0.0 -1. 14SRB 3400.00 50.0 240.0 -1. 14SRC 3400.00 50.0 120.0 -1. C the source for the analogue network 14INERS -1 0.000001 .0000001 -1. C the source of any additional load applied 14INERS -1 -0.000001 0.0000001 5.00 19 { Begin U.M. data, which will consist of a single 3-phase induction motor C No autoinitialization. Also, this data does use compensation on both sides. 0 0 BLANK 4 111INER 2 0.15700 0.465131 0.465131 C 3 armature coils come first: MOTA 1 0.3393571 0.0109180 MOTB 1 0.3393571 0.0109180 MOTC 1 C 3 rotor coils follow: 0.5785651 0.0109180 XOTB 1 TACS R-thev ROTRES 0.5785651 0.0109180 XOTC 1 TACS R-thev ROTRES XOTA 1 TACS R-thev ROTRES C Note about preceding. The request word is "TACS R-thev " in columns 63 C through 74. It is to be followed by the A6 name of a TACS variable in C columns 75-80. Each coil can have a different TACS variable (the three C being the same means that the 3-phase resistance is balanced. BLANK card ending U.M. data cards BLANK card ending source cards BLANK card ending output specifications (none here) 193 .4 0.0 2.0 0.0 200.UM-1 OMEGM C 193 .2 0.0 2.0 UM-1 IE1 { A vector plot really is needed for this BLANK card ending plot cards BEGIN NEW DATA CASE C 8th of 12 subcases is similar to 1st (Janko's data), but is slightly C different. Prior to correction on 14 October 2000, the solution was C wrong following closure of the switch to ground. Data was contributed C by Steve Nurse of Reyrolle in England on 11 Oct 2000. The following C data is a smaller, simpler illustration of the problem: C C GEN SEND Type-91 LOAD 1.0 C o------_-------||------/\/\/\-------||------/\/\/\-------|| E C ^ || R = 1 || R = 1 || a C MEASURING || || || r C switch || || || t C (always ||--------/----------||--------/----------|| h C closed) 1st, switch 2nd, switch C will open will close C C There are 3 switches, with one permanently closed. All 3 touch the C compensation-based Type-91 element. Simulation begins with the Type-91 C element shorted, so the source GEN feeds the 1-ohm resistor from C LOAD to ground. The current is 6 volts / 1 ohm = 6 amps. But then C the switch from SEND to LOAD opens, inserting the Type-91, adding C another ohm. This drops the current to 3 amps. Finally, the ordinary C 1-ohm resistor from LOAD to ground is shorted, restoring the current C to 6 amps. Prior to correction, the final change did not happen. See C also the 8th (final) subcase of DC-68, for which answers changed C slightly as a result of this change to OVER16. PRINTED NUMBER WIDTH, 10, 1, { Request maximum precision (for 8 output columns) .001 .009 1 -1 1 91SEND LOAD 3333. 1 C -------R(tr)---------><-----tr----------------> 0.001 { V of flashover gap 1.0 0.0 1.0 .020 9999 SEND 10.E8 { Avoid ATP warning and such internal addition LOAD 1.0 1 BLANK card ending branch cards C Switch data: T-close T-open I-epsiln LOAD SEND -1.0 .0025 1.E8 1 LOAD .0065 1.0 1 GEN SEND MEASURING 1 BLANK card ending switch cards 11GEN 6.0 BLANK card ending source cards GEN SEND LOAD C First 3 output variables are electric-network voltage differences (upper voltage minus lower voltage); C Next 5 output variables are branch currents (flowing from the upper node to the lower node); C Step Time GEN SEND LOAD LOAD LOAD GEN SEND LOAD C SEND TERRA SEND LOAD TERRA C *** Switch "LOAD " to "SEND " closed before 0.00000000E+00 sec. C *** Switch "GEN " to "SEND " closed before 0.00000000E+00 sec. C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 C 1 .1E-2 6.0 6.0 6.0 -6. 0.0 6.0 0.0 6.0 C 2 .002 6.0 6.0 6.0 -6. 0.0 6.0 0.0 6.0 C *** Open switch "LOAD " to "SEND " after 3.00000000E-03 sec. C 3 .003 6.0 6.0 6.0 0.0 0.0 6.0 0.0 6.0 C 4 .004 6.0 6.0 0.0 0.0 0.0 .6E-8 0.0 0.0 C 5 .005 6.0 6.0 3.0 0.0 0.0 3.0 3.0 3.0 C 6 .006 6.0 6.0 3.0 0.0 0.0 3.0 3.0 3.0 C *** Close switch "LOAD " to " " after 7.00000000E-03 sec. C 7 .007 6.0 6.0 3.0 0.0 0.0 3.0 3.0 3.0 C 8 .008 6.0 6.0 0.0 0.0 6.0 6.0 6.0 0.0 C 9 .009 6.0 6.0 0.0 0.0 6.0 6.0 6.0 0.0 BLANK card ending output variables requests (just node voltages, here) BLANK card ending plot cards BEGIN NEW DATA CASE C 9th of 12 subcases is unrelated to preceding subcases. C Illustrate that the 3-phase compensation of DC-34 can coexist with single- C phase compensation of a disconnected Type-93 nonlinear reactor. Prior to C correction to SOLVUM on 20 May 2001, the solution was obviously wrong. C The first report of such trouble came from Dr. Michael Steurer of CAPS at C Florida State University in Tallahassee. His data involved a Type-91 C TACS-controlled resistor, but in fact TACS had nothing to do with the C trouble, so is not being involved in the present illustration. C Comment about 6th subcase: It, too, involves both 3-phase compensation C of the U.M. and 1-phase compensation of a separate element. But that C data was illegal because the two were not isolated in disconnected C subnetworks. POWER FREQUENCY, 60.0, { Europeans need this (LEC letter dated 6 Jan 89, page 2) PRINTED NUMBER WIDTH, 11, 1, { Restore default STARTUP column width, separation .001 .500 1 1 1 1 1 -1 5 5 20 20 100 100 0 BUS-F0 1.0 0 BUS-A0 1.0 0 BUS-B0 1.0 0 BUS-C0 1.0 0BUS-A1BUS-A0 0.02 1.0610 0BUS-B1BUS-B0BUS-A1BUS-A0 0BUS-C1BUS-C0BUS-A1BUS-A0 0BUS-M1 2.00E6 1 0BUS-M0 2.00E6 1 0BUS-M0BUS-M1 1. 1 C End of DC-34 branches; begin branches from DC-45b. Note that these C involve different names, so are completely disconnected from the machine. C This single-phase subnetwork will illustrate an energization transient: GEN TRAN 2.0 1 93TRAN 1.0 1.0 1 0.0 0.0 0.9 0.9 2.0 1.1 10. 1.2 9999 BLANK card ends the last branch card BLANK card ends (in this case nonexistent) switch cards 11BUS-F0 0.002091 { 18 Oct 90, remove ineffective T-start = -1.0 from card 11BUS-M0-1 1.02 11BUS-M0-1 -.4 .019500 14BUS-A0 1.41421356 60.0 0.0 14BUS-B0 1.41421356 60.0 -120.0 14BUS-C0 1.41421356 60.0 120.0 C End of DC-34 static sources. B4 the U.M., insert the source from DC-45b: 14GEN 377. 60. 19 UM { Beginning of U.M. data (Type-19 source) 1 BLANK card ending class-1 U.M. data 1 2 1111BUS-M1 1 1786.98 { 1st card of U.M. machine table 1.0 1.550 0 0.3 1.5 1.0 0.93787 1.490 0 0.0 0.0 BUS-A1 1 1.3860 { 1st card of coil table 0.001096 0.150 BUS-B1 1 -0.95877 0.001096 0.150 BUS-C1 1 -0.42721 0.00074 0.101 BUS-F0 1 -2.826 0.0131 0.055 1 0.0540 0.036 1 BLANK card ending all U.M. data BLANK card ending EMTP source cards C Total network loss P-loss by summing injections = 2.999999989932E+00 2BUS-A1 1.56413 { First of many initial condition cards for the 2BUS-B1 -0.30745 { electric network. Since the U.M. is not a 2BUS-C1 -1.25677 { part of the phasor solution (see DCNEW-1 for 2BUS-A0 1.41421356 { such a more modern problem), synchronous 2BUS-B0 -0.70710678 { operation can begin smoothly only if the 2BUS-C0 -0.70710678 { initially conditions are manually applied. 2BUS-F0 0.002091 2BUS-M0 1.0 2BUS-M1 1.0 { Final card of node voltage initial conditions 3BUS-A1BUS-A0 1.38494 { 1st card of branch current initial conditions 3BUS-B1BUS-B0 -0.95793 3BUS-C1BUS-C0 -0.42701 3 BUS-A0 -1.41421356 3 BUS-B0 +0.70710678 3 BUS-C0 +0.70710678 3 BUS-F0 -0.002091 3BUS-M0BUS-M1 +1.01 3BUS-M0 +1.0 3BUS-M1 +1.0 { Last card of branch current init. condit. TRAN BUS-M1BUS-M0 { Request for selective node voltage output C First 3 output variables are electric-network voltage differences (upper voltage minus lower voltage); C Next 5 output variables are branch currents (flowing from the upper node to the lower node); C Final 9 output variables pertain to Type-19 U.M. components (names are generated internally); C Step Time TRAN BUS-M1 BUS-M0 TRAN BUS-M1 BUS-M0 BUS-M0 GEN UM-1 UM-1 C TERRA TERRA TERRA BUS-M1 TRAN TQGEN OMEGM C UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 UM-1 C THETAM IPA IPB IPC IE1 IE2 IE3 C 0 0.0 0.0 1.0 1.0 0.0 0.0 0.0 1.01 0.0 1.00124132 1.0 C .93787 1.386 -.95877 -.42721 -2.826 0.0 0.0 C 1 .1E-2 350.17556 1.00000302 1.0000025 .17508778 .012094107 .010000262 1.00999974 .17508778 .997905631 1.0 C .93787 1.40049668 -.50888659 -.89163009 -2.8259503 .85956E-4 -.00481942 C 2 .002 273.846974 1.00000988 1.0000075 .487099047 .015338089 .010001714 1.00999829 .487099047 .994660197 1.0 C .93787 1.22097729 .010653511 -1.2316108 -2.8258183 .309395E-3 -.0092568 C 3 .003 159.111636 1.00001841 1.0000125 .703578352 .018762645 .010005857 1.00999414 .703578352 .991231498 1.0 C .93787 .870655506 .528248848 -1.3989244 -2.8248425 .001999525 -.01265851 BLANK card ending the specification of output variables C 500 0.5 377.026496 1.00472639 1.00854631 -.01324809 -.13068128 .049259392 .570740608 -.01324809 .701421892 1.00469204 C .660491474 .974954484 -.86183035 -.11310413 -2.4109229 -.04254056 .118833141 C Variable max: 377.260976 1.00779145 1.00854631 1.25309712 .260847165 .297988719 1.01 1.25309712 1.0019948 1.00877611 C .93787 1.41537866 1.40559809 1.40738077 -2.0788757 .133544596 .152562874 C Times of max: .05 .442 0.5 .471 .378 .296 0.0 .471 .015 .439 C .064 .034 .056 .045 .358 .176 .471 C Variable min: -376.94611 .989230544 .990404468 -6.3622515 -.29505746 -.3892114 .322011281 -6.3622515 .29918812 .988252614 C .361026488 -1.4043962 -1.4122326 -1.411928 -2.826 -.04419925 -.18596365 C Times of min: .475 .158 .217 .012 .102 .02 .296 .012 .346 .159 C .358 .059 .031 .02 0.0 .485 .188 PRINTER PLOT 193 .1 0.0 1.0 UM-1 THETAM { Plot limits: (0.000, 9.379) C The preceding plot is identical to DC-34, so provides the simplest way to C validate the solution of the U.M. As for the Type-93 reactor, inrush C current is shown by the following vector plot. It seems believable, and C is unchanged by deletion of all data related to DC-34. CALCOMP PLOT 193.05 0.0 0.5 -7.0 1.0TRAN BLANK card ending all plot cards BEGIN NEW DATA CASE C 10th of 12 subcases is derived from 1st to demonstrate that the TACS dc C solution in fact represents the superposition of the network solution C within TACS (driven by Type-11 sources, which are batteries) and the C user-supplied initial conditions. 1.E-5 1.E-5 1 1 1 2 TACS HYBRID C 99RESIS = 1.E-6 C 77RESIS 1.E-6 C Replace the preceding 2 lines (see 1st subcase) in attempt to demonstrate C that the TACS dc solution printout is the sum of the TACS network solution C value and the initial condition value: -37.E-6 + 38.E-6 = 1.E-6 (same as C the 1st subcase). Unfortunately, this has not yet been accomplished. But C different treatment of the initial condition is easily demonstrated. The C following shows 4 different ways to define a TACS variable. All will be C given the same initial condition value 38.E-6. But two of the four will C have one value for dc printout, and the remaining two will have another. 11SOURCE -37.E-6 { Note T-start < 0 so TACS phasor solve } -1. 1. C Following RESA and RESB will show source value (-37.E-6) in dc printout: 0RESA +SOURCE 1RESB +SOURCE 1. 1. 1. C Following RESC & RESD will show initial condition valu (+38) in dc printout: 1RESC +SOURCE 1. 1. 88RESD = 1.E-6 77RESA 38.E-6 77RESB 38.E-6 77RESC 38.E-6 77RESD 38.E-6 33RESA RESB RESC RESD BLANK card ending TACS 91N1 N2 TACS RESA 1 N1 N12 1.E-6 N2 1. BLANK card ending branches N12 N2 MEASURING SRCE N1 MEASURING 1 BLANK card ending switches 14SRCE 100. 60. 0.0 0. -1. BLANK card ending sources C Zero-frequency (dc) steady-state solution for TACS follows. C (Name) TACS value (Name) TACS value (Name) TACS value (Name) TACS value (Name) TACS value C RESA -3.70000000E-05 RESB -3.70000000E-05 RESC 3.80000000E-05 RESD 3.80000000E-05 BLANK card ending output variables (none) BLANK card ending plot BEGIN NEW DATA CASE C 11th of 12 subcases always should be last since it demonstrates a halt C to execution using CALL STOPTP. This is added 29 August 2002 to C illustrate what previously was a tight loop within Z-thev computation C of OVER16. There is similarity to the 6th subcase. But whereas that C data would simulate (giving the wrong answer), this data would die. C The original complaint was in E-mail of the EEUG list server dated C 22 August 2002. In this, Alejandro Montenegro at the University of C Florida asked: "The problem appears when I combine TACS-controlled C Type 13 switches, Type 94 components and UM type 3." But the present C data is much simpler. Only the Type-3 U.M. is used. No TACS, no C MODELS, and no switches at all. The fundamental problem was an overlap C of compensation, and it led to an unbounded index prior to the trap C that this subcase illustrates. .000100 .15 1 -1 1 C Begin with 3 single-phase nonlinear reactors in the same subnetwork as the C U.M. It is not clear why all 3, rather than just 1, is required for C the tight loop. But this seems to be the case. 93VPA VA .559 300. 1 0.0 0.0 .5590 300. .9344 400. 9999 93VPB VB VPA VA .559 300. 1 93VPC VC VPA VA .559 300. 1 C Nodes VA, VB, and VC have no connected linear branch. To prevent warning C messages about floating subnetwork, & automatic additions to ground, add C: VA 1.0 { Shunt capacitance avoids floating VB 1.0 { Shunt capacitance avoids floating VC 1.0 { Shunt capacitance avoids floating C The following 2 lines are rotor mass and damping from Bonfanti's motor: COPPIA 9.16E6 2 COPPI2COPPIA 1.0E-6 BLANK card ending branch cards BLANK card ending switches 14VPA 10777.75 60. -1. 10. 14VPB 10777.75 60. -120. -1. 10. 14VPC 10777.75 60. 120. -1. 10. C The following is Bonfanti's motor from the 3rd subcase. 14COPPI2-1 -1. 1.E-9 0. -1.0 9999. 19 UM 01 0 - Compensation; change 0 to 1 if prediction is wanted BLANK general UM specification 3 1 1111COPPIA 1 0.0 0.3964 0.0 0.3964 4. COPPI2 .1674 .001 VA 1 .1674 .01 VB 1 .1674 .01 VC 1 .7819 .00453 1 .7819 .00453 1 BLANK card ending U.M. data BLANK card ending sources VA VB VC VPA VPB VPC { Names of nodes for node voltage output BLANK card ending output variables BLANK card ending plot C Document the end of output of this 11th data subcase. This begins C exactly as in years past. An error is recognized after the dT-loop heading: C 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 C Error. Halt execution above S.N. 2322 in OVER16. Compensation is in error. Bad subscript N1. If U.M. is involved, try adding C VERIFY U.M. COMPENSATION or STEP ZERO COUPLE as illustrated in DCN16. C Temporary error stop in ENTRY STOPTP of "WINDOW". NCHAIN, LASTOV = 16 15 C 12345678901234567890123456789012345678901234567890123456789012345678901234567890 C ABUFF(1:80) = C C So much for the old. The new adds explanation of the recovery. C Before output of the 12th subcase begins, there is this: C C ------------------------------------------------------------------------------------------------------------------------------------ C ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ C ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ C ------------------------------------------------------------------------------------------------------------------------------------ C C This is not actually a KILL error termination. But ATP now will try to continue as if it were. A recursive CALL is involved. C C ------------------------------------------------------------------------------------------------------------------------------------ C ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ C ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ERROR/ C ------------------------------------------------------------------------------------------------------------------------------------ BEGIN NEW DATA CASE C 12th of 12 subcases demonstrates that the STOPTP termination of the C preceding subcase is not fatal. By means of a recursive CALL, ATP C can be restarted to continue processing following data subcases such C as this one, which is a copy of the 1st subcase. Output is identical. C For brevity, connectivity and the phasor output have been omitted. C This modification is made 23 January 2011. WSM. 1.E-5 1.E-5 1 1 TACS HYBRID 99RESIS = 1.E-6 77RESIS 1.E-6 33RESIS BLANK card ending TACS 91N1 N2 TACS RESIS 1 N1 N12 1.E-6 N2 1. BLANK card ending branches N12 N2 MEASURING SRCE N1 MEASURING 1 BLANK card ending switches 14SRCE 100. 60. 0.0 0. -1. BLANK card ending sources N1 N12 N2 BLANK card ending voltage printout BLANK card ending plot BEGIN NEW DATA CASE BLANK