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/dcn25old.dat | 836 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 836 insertions(+) create mode 100644 benchmarks/dcn25old.dat (limited to 'benchmarks/dcn25old.dat') diff --git a/benchmarks/dcn25old.dat b/benchmarks/dcn25old.dat new file mode 100644 index 0000000..1e6123d --- /dev/null +++ b/benchmarks/dcn25old.dat @@ -0,0 +1,836 @@ +BEGIN NEW DATA CASE +C 1st of 4 subcases illustrates Robert Meredith's Type-68 TACS device as +C first described in the April, 1998, newsletter. This is AUT5 data case. +C test120a.dat (1 HERTZ AT 120% VOLTAGE) +C without resistances +C ENERGIZED AT ZERO TIME VIA 0.100 OHM RESISTORS -steady state +C FREE FORMAT TIME STEP- 1000/CYCLE FOR 1.5 SECOND+ +C .000016666666, 1.500 , 60.000 , 0.,,,,,,, { Meredith's original params + .000066666667, .100 , 60.000 , 0.,,,,,,, { Much reduced burden on computer +C PRN> PLT> NET> SS> MAX> SAVE> NENERG> > + 1 5 0 0 1 -1 0 + 5 5 20 20 100 100 500 500 +TACS HYBRID +C $INCLUDE aut5-tac.inc +C $INCLUDE aut5-tck.inc { Uses new Type-68 TACS device, so much smaller file +C TACS STEADY-STATE CORE NFLUX INITIALIZATION FOR TRANSFORMER AUT5 +C NFLUX IS PROPORTIONAL TO THE INDUCTIVE CURRENT IN THESE SWITCHES: +C TYPE 91 TACS SOURCES DERIVED FROM NETWORK SWITCH CURRENTS: +C <----A---><----A---><----A---> +91AUT5D1 -1.0 999. +91AUT5D2 -1.0 999. +91AUT5D3 -1.0 999. +91AUT5D4 -1.0 999. +91AUT5D5 -1.0 999. +C ZERO-ORDER BLOCKS TO SCALE ABOVE TO NFLUX LINKED (NREF*PHI) WITHIN CORE VOLUME +C GAIN IS BASED ON 10000-HENRY FLUX-MONITORING INDUCTANCES. +C + + + + + + 0AUT5E1 +AUT5D1 1.E4 + 0AUT5E2 +AUT5D2 1.E4 + 0AUT5E3 +AUT5D3 1.E4 + 0AUT5E4 +AUT5D4 1.E4 + 0AUT5E5 +AUT5D5 1.E4 +C NEW TACS DEVICE 68 MODELS HYSTERETIC LOSS/RESIDUAL MMF COMPONENTS. +C 68 +99AUT5T168AUT5E1 13.821 27.642 2.91727E-3 8.44316E-4 +99AUT5T268AUT5E1 27.642 62.194 2.43012E-3 2.53199E-4 +99AUT5T368AUT5E1 62.194 82.925 1.23361E-3 4.76044E-4 +99AUT5T468AUT5E1 82.925 103.66 1.22440E-3 5.90543E-4 +99AUT5T568AUT5E1 103.66 120.25 1.79633E-3 1.56980E-3 +99AUT5T668AUT5E1 120.25 131.34 1.44373E-3 3.08189E-3 +99AUT5T768AUT5E1 131.34 139.51 4.95355E-3 2.07092E-2 +99AUT5U168AUT5E2 13.821 27.642 2.91727E-3 8.44316E-4 +99AUT5U268AUT5E2 27.642 62.194 2.43012E-3 2.53199E-4 +99AUT5U368AUT5E2 62.194 82.925 1.23361E-3 4.76044E-4 +99AUT5U468AUT5E2 82.925 103.66 1.22440E-3 5.90543E-4 +99AUT5U568AUT5E2 103.66 120.25 1.79633E-3 1.56980E-3 +99AUT5U668AUT5E2 120.25 131.34 1.44373E-3 3.08189E-3 +99AUT5U768AUT5E2 131.34 139.51 4.95355E-3 2.07092E-2 +99AUT5V168AUT5E3 13.821 27.642 2.91727E-3 8.44316E-4 +99AUT5V268AUT5E3 27.642 62.194 2.43012E-3 2.53199E-4 +99AUT5V368AUT5E3 62.194 82.925 1.23361E-3 4.76044E-4 +99AUT5V468AUT5E3 82.925 103.66 1.22440E-3 5.90543E-4 +99AUT5V568AUT5E3 103.66 120.25 1.79633E-3 1.56980E-3 +99AUT5V668AUT5E3 120.25 131.34 1.44373E-3 3.08189E-3 +99AUT5V768AUT5E3 131.34 139.51 4.95355E-3 2.07092E-2 +99AUT5W168AUT5E4 14.120 28.239 2.69692E-3 7.64028E-4 +99AUT5W268AUT5E4 28.239 63.538 2.24656E-3 2.29122E-4 +99AUT5W368AUT5E4 63.538 84.717 1.14043E-3 4.30775E-4 +99AUT5W468AUT5E4 84.717 105.90 1.13191E-3 5.34391E-4 +99AUT5W568AUT5E4 105.90 122.84 1.66066E-3 1.42058E-3 +99AUT5W668AUT5E4 122.84 134.17 1.33473E-3 2.79024E-3 +99AUT5W768AUT5E4 134.17 142.42 4.58432E-3 1.91911E-2 +99AUT5X168AUT5E5 14.120 28.239 2.69692E-3 7.64028E-4 +99AUT5X268AUT5E5 28.239 63.538 2.24656E-3 2.29122E-4 +99AUT5X368AUT5E5 63.538 84.717 1.14043E-3 4.30775E-4 +99AUT5X468AUT5E5 84.717 105.90 1.13191E-3 5.34391E-4 +99AUT5X568AUT5E5 105.90 122.84 1.66066E-3 1.42058E-3 +99AUT5X668AUT5E5 122.84 134.17 1.33473E-3 2.79024E-3 +99AUT5X768AUT5E5 134.17 142.42 4.58432E-3 1.91911E-2 +C CURRENTS INJECTED INTO NETWORK FOR HYSTERESIS MODELING AND OTHER COMPENSATION. +C PHASE LEGS HAVE CURRENTS INJECTED TO COMPENSATE COUPLING XFMR. MAGNET. AMPS., +C FLUX-MONITORING INDUCTANCE CURRENT AND SMALL CORE GAPS: +99AUT5TP =+AUT5D1/( 1.30893E-2)+ 1*(+AUT5T1+AUT5T2+AUT5T3+AUT5T4+AUT5T5) +99AUT5A5 =+AUT5TP+ 1*(+AUT5T6+AUT5T7) +99AUT5TQ =+AUT5D2/( 1.30893E-2)+ 1*(+AUT5U1+AUT5U2+AUT5U3+AUT5U4+AUT5U5) +99AUT5B5 =+AUT5TQ+ 1*(+AUT5U6+AUT5U7) +99AUT5TR =+AUT5D3/( 1.30893E-2)+ 1*(+AUT5V1+AUT5V2+AUT5V3+AUT5V4+AUT5V5) +99AUT5C5 =+AUT5TR+ 1*(+AUT5V6+AUT5V7) +C YOKES ARE COMPENSATED FOR MONITORING INDUCTANCE CURRENTS: +99AUT5TS =+AUT5D4+ 1*(+AUT5W1+AUT5W2+AUT5W3+AUT5W4+AUT5W5) +99AUT501 =+AUT5TS+ 1*(+AUT5W6+AUT5W7) +99AUT5TT =+AUT5D5+ 1*(+AUT5X1+AUT5X2+AUT5X3+AUT5X4+AUT5X5) +99AUT502 =+AUT5TT+ 1*(+AUT5X6+AUT5X7) +C YOKE & PHASE HYSTERETIC CURRENTS ARE EXTRACTED, EXCEPT AT GROUND: +99AUT5B0 =-AUT501-AUT5B5 +99AUT5C0 =-AUT502-AUT5C5 +C End of $INCLUDE. File name = aut5-tck.inc +C PARTS OF HYSTERESIS COMPONENTS FOR PLOTTING +C FORTRAN STATEMENTS; 99= INPUT; 98= OUTPUT; 88= INSIDE +C =< FREE FORMAT FORTRAN TO COL 80 -----> +C TOTAL OF ALL HYSTER INJECTIONS FOR LEGS A,B,C +99ALLINA =+AUT5T1+AUT5T2+AUT5T3+AUT5T4+AUT5T5+AUT5T6+AUT5T7 +99ALLINB =+AUT5U1+AUT5U2+AUT5U3+AUT5U4+AUT5U5+AUT5U6+AUT5U7 +99ALLINC =+AUT5V1+AUT5V2+AUT5V3+AUT5V4+AUT5V5+AUT5V6+AUT5V7 +C +C TYPE 90 TACS SOURCES DERIVED FROM NETWORK VOLTAGES: +C <----A---><----A---><----A---> +C TYPE 91 TACS SOURCES DERIVED FROM NETWORK SWITCH CURRENTS (1ST NODE) +C <----A---><----A---><----A---> +91AUT5AM 0. 999. +91AUT5BM 0. 999. +91AUT5CM 0. 999. +C FORTRAN STATEMENTS; 99= INPUT; 98= OUTPUT; 88= INSIDE +C =< FREE FORMAT FORTRAN TO COL 80 -----> +C NEXT IS SUM OF COUPLING OUTPUT AND ALL INJECTED AMPS - USED BY CORE MODELS +98LEGAMP = AUT5AM + AUT5A5 +98LEGBMP = AUT5BM + AUT5B5 +98LEGCMP = AUT5CM + AUT5C5 +C NEXT IS APPARENT LEG AMPS FOR HYSTER PLOTTING +98NETAMP = LEGAMP - ALLINA +98NETBMP = LEGBMP - ALLINB +98NETCMP = LEGCMP - ALLINC +C S-BLOCKS OF ORDER 1 IN COL 2 INTEGRATE VOLTAGES TO GET FLUXES (TIMES 250.0). +C + + + + + +C N0 & D0> +C CD+< IN1> +< IN2> +< IN3> +< IN4> +< IN5> < A >< B >< C >< D >< E > +C TACS OUTPUT REQUESTS - TYPE 33 +C +33AUT5E1 ALLINA LEGAMP NETAMP +33AUT5E2 ALLINB LEGBMP NETBMP +33AUT5E3 ALLINC LEGCMP NETCMP +C 33AUT5T1 AUT5T2 AUT5T3 AUT5T4 AUT5T5 AUT5T6 AUT5 +33AUT5U1 AUT5U2 AUT5U3 AUT5U4 AUT5U5 AUT5U6 AUT5U7 +BLANK CARD ENDING ALL ATP-SORTED TACS CARDS (from blank.inc) +C THIS MUST BE EDITED TO ADD AUT5A5 >> AUT5AM, ETC METERING +C $INCLUDE aut5-brn.inc +C WSM makes non-Unix by adding a card using DOS editor. +C 3-PHASE, 3-LEG SATURABLE XFMR. MODEL FROM MODELAD.FTN - 60. HZ IMPEDANCES +C CHANGED AUT5A5, B5, C5 TO AUT5AM, BM, CM ; NOW REQUIRES METERING SWITCHES +C AUT5 = IDENTIFYING BUS CODE PREFIX OF INTERNAL NODES. +C AUT5AM = NAME OF A METERING NODE ADDED AT AUT5A5 +C EACH PHASE LEG HAS 3 WINDINGS WHOSE RATED VOLTAGES AND TURNS ARE: +C WINDING 1 (INNERMOST) = 13.130 KV. TURNS = 44.998 +C WINDING 1 INCLUDES THE EFFECT OF 0.5670 OHMS/PHASE INTERNAL (DELTA) REACTORS. +C REACTOR RESISTANCE/PHASE MUST BE INCLUDED WITH ENTERED WINDING RESISTANCES. +C WINDING 2 = 79.674 KV. TURNS = 273.052 +C WINDING 3 = 119.655 KV. TURNS = 410.072 +C INTERNAL REFERENCE = 29.179 KV. TURNS = 100.000 +C *********BEGIN USER COMMENTS ************ +C EAST GARDEN CITY AUTOTRANSFORMER*** one *** AT 345/138 KV TAP +C ABBREVIATED HYSTERESIS MODEL; SEVEN OF POSSIBLE 11 COMPONENTS +C .95 * .7274 SQUARE METER LEG; .97 * .7274 SQUARE METER YOKE +C ACTUAL DELTA IS GROUNDED BETWEEN OUTER PHASES AT Y1. +C EACH DELTA REACTOR X =.82;OR .5467 WHEN EQUALIZED; REACTOR R = .00427 EACH +C ACTUAL DELTA WINDING R = .00510; + .002847 TO EQUALIZE = .0079467 ENTERED +C AIR-CORE COUPLINGS FROM COUPCOIL, BASED ON SMIT DIMENSIONS. +C B-H CURVE BASED ON ARMCO TRAN-COR H-0 from smit +C B-H SHIFTED TO HIGHER CURRENTS IN 1.5 - 1.85 T RANGE +C MILD OPEN HEARTH STEEL FOR TANK +C TANK HAS 21.65 M CIRCUM; FLUX SHIELD = 3.15 M HIGH +C 64-36 EDDY-HYST LOSS RATIO AT RATED; +C ********** END USER COMMENTS ************ +C COUPLING, LEAKAGE IMPEDANCE & WINDING LOSS MODEL FOR TRANSFORMER AUT5: +C < R >< X >< R >< X >< R >< X > +51Y1 Y2 .01016 10124.87377 +52EGCT1A 0.000 61435.16312.04458372798.88582 +53EGCHEAEGCT1A 0.000 92263.79188 0.000559873.90936.12388840864.72923 +54AUT5A1 0.000 22499.42979 0.000136528.47263 0.000205039.49115 + 0.000 50000.88921 +55AUT5A2 0.000 22499.42071 0.000136530.26291 0.000205042.81249 + 0.000 50000.86902 0.000 50001.52468 +56AUT5A3 0.000 22499.42350 0.000136530.70112 0.000205053.04695 + 0.000 50000.87522 0.000 50001.68516 0.000 50004.18094 +57AUT5A4 0.000 22499.25174 0.000136528.91954 0.000205046.26382 + 0.000 50000.49352 0.000 50001.03270 0.000 50002.52680 + 0.000 50002.52680 +58AUT5AM 0.000 22499.02967 0.000136526.09974 0.000205035.90210 + 0.000 50000.00000 0.000 50000.00000 0.000 50000.00000 + 0.000 50000.00000 0.000 50000.00000 +51Y2 Y3 .01016 10124.87377 +52EGCT1B 0.000 61435.16312.04458372798.88582 +53EGCHEBEGCT1B 0.000 92263.79188 0.000559873.90936.12388840864.72923 +54AUT5B1AUT5B0 0.000 22499.42979 0.000136528.47263 0.000205039.49115 + 0.000 50000.88921 +55AUT5B2AUT5B0 0.000 22499.42071 0.000136530.26291 0.000205042.81249 + 0.000 50000.86902 0.000 50001.52468 +56AUT5B3AUT5B0 0.000 22499.42350 0.000136530.70112 0.000205053.04695 + 0.000 50000.87522 0.000 50001.68516 0.000 50004.18094 +57AUT5B4AUT5B0 0.000 22499.25174 0.000136528.91954 0.000205046.26382 + 0.000 50000.49352 0.000 50001.03270 0.000 50002.52680 + 0.000 50002.52680 +58AUT5BMAUT5B0 0.000 22499.02967 0.000136526.09974 0.000205035.90210 + 0.000 50000.00000 0.000 50000.00000 0.000 50000.00000 + 0.000 50000.00000 0.000 50000.00000 +51Y3 Y1 .01016 10124.87377 +52EGCT1C 0.000 61435.16312.04458372798.88582 +53EGCHECEGCT1C 0.000 92263.79188 0.000559873.90936.12388840864.72923 +54AUT5C1AUT5C0 0.000 22499.42979 0.000136528.47263 0.000205039.49115 + 0.000 50000.88921 +55AUT5C2AUT5C0 0.000 22499.42071 0.000136530.26291 0.000205042.81249 + 0.000 50000.86902 0.000 50001.52468 +56AUT5C3AUT5C0 0.000 22499.42350 0.000136530.70112 0.000205053.04695 + 0.000 50000.87522 0.000 50001.68516 0.000 50004.18094 +57AUT5C4AUT5C0 0.000 22499.25174 0.000136528.91954 0.000205046.26382 + 0.000 50000.49352 0.000 50001.03270 0.000 50002.52680 + 0.000 50002.52680 +58AUT5CMAUT5C0 0.000 22499.02967 0.000136526.09974 0.000205035.90210 + 0.000 50000.00000 0.000 50000.00000 0.000 50000.00000 + 0.000 50000.00000 0.000 50000.00000 +C PHASE LEG N-L MODELS FOR TRANSFORMER: AUT5 +C PHASE LEG CORES ARE 3.60 METERS LONG WITH AREAS OF 0.6910 SQUARE METERS +C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS: +C LINEAR:(0.20, 4.00) AND FULLY SATURATED:(2.00, 15000.00) +C FOR SS CALCS EACH LEG HAS REACT. OF 15740. OHMS (AT 100.0 TURNS) +C BASED ON SLOPE OF POINT 8 BELOW, HAVING A FLUX DENSITY OF 1.7400 TESLA. +C HYSTERESIS CURRENT OFFSET OF 3.115 AMPS IS MODELED AT SATURATION. +C LAMINATIONS ARE MODELED BY 8 LAYERS; EDDY LOSSES USE ENTERED CORE RESISTIVITY. +C

(--FOR TYPE 98 N/L REACTOR--) P +C CURRENT -----> +98AUT5A5AUT5A6 4.912234.877 + 0.14400 3.45524 + 0.39528 6.91044 + 0.81635 15.54839 + 1.12175 20.73117 + 1.41336 24.18640 + 2.00198 25.91428 + 3.13598 28.16067 + 4.23399 30.06154 + 8.81421 31.96448 + 19.73322 32.83478 + 56.33246 33.72028 + 543.11535 34.87714 + 1083.11535 35.20278 + 9999 +98AUT5A6AUT5A7AUT5A5AUT5A64.912234.877 +98AUT5A7AUT5A8AUT5A5AUT5A64.912234.877 +98AUT5A8 AUT5A5AUT5A64.912234.877 +98AUT5B5AUT5B6AUT5A5AUT5A64.912234.877 +98AUT5B6AUT5B7AUT5A5AUT5A64.912234.877 +98AUT5B7AUT5B8AUT5A5AUT5A64.912234.877 +98AUT5B8AUT5B0AUT5A5AUT5A64.912234.877 +98AUT5C5AUT5C6AUT5A5AUT5A64.912234.877 +98AUT5C6AUT5C7AUT5A5AUT5A64.912234.877 +98AUT5C7AUT5C8AUT5A5AUT5A64.912234.877 +98AUT5C8AUT5C0AUT5A5AUT5A64.912234.877 +C HIGH RESISTANCE ACROSS EACH UNPARALLELED TYPE 98 ELEMENT ABOVE: +C <----><---->RRRRRRXXXXXXCCCCCC + AUT5A6AUT5A5 .1E+11 + AUT5A7AUT5A6AUT5A6AUT5A5 + AUT5A8AUT5A7AUT5A6AUT5A5 + AUT5B6AUT5B5AUT5A6AUT5A5 + AUT5B7AUT5B6AUT5A6AUT5A5 + AUT5B8AUT5B7AUT5A6AUT5A5 + AUT5C6AUT5C5AUT5A6AUT5A5 + AUT5C7AUT5C6AUT5A6AUT5A5 + AUT5C8AUT5C7AUT5A6AUT5A5 +C YOKE MODELS FOR TRANSFORMER: AUT5 +C EACH TOP+BOTTOM YOKE LENGTH IS 3.40 METERS WITH AREA OF 0.7060 SQUARE METERS +C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS: +C LINEAR:(0.20, 4.00) AND FULLY SATURATED:(2.00, 15000.00) +C FOR SS CALCS EACH YOKE HAS REACT. OF 17026. OHMS (AT 100.0 TURNS) +C BASED ON SLOPE OF POINT 8 BELOW, HAVING A FLUX DENSITY OF 1.7400 TESLA. +C HYSTERESIS CURRENT OFFSET OF 2.943 AMPS IS MODELED AT SATURATION. +C LAMINATIONS ARE MODELED BY 8 LAYERS; EDDY LOSSES USE ENTERED CORE RESISTIVITY. +C

(--FOR TYPE 98 N/L REACTOR--) P +C CURRENT -----> +98AUT501AUT50A 4.636135.606 + 0.13600 3.52991 + 0.37332 7.05979 + 0.77100 15.88444 + 1.05943 21.17924 + 1.33484 24.70914 + 1.89076 26.47434 + 2.96176 28.76924 + 3.99877 30.71116 + 8.32457 32.65503 + 18.63694 33.54364 + 53.20871 34.44660 + 512.94277 35.60585 + 1022.94277 35.91340 + 9999 +98AUT50AAUT50BAUT501AUT50A4.636135.606 +98AUT50BAUT50CAUT501AUT50A4.636135.606 +98AUT50CAUT5B0AUT501AUT50A4.636135.606 +98AUT502AUT50FAUT501AUT50A4.636135.606 +98AUT50FAUT50GAUT501AUT50A4.636135.606 +98AUT50GAUT50HAUT501AUT50A4.636135.606 +98AUT50HAUT5C0AUT501AUT50A4.636135.606 +C HIGH RESISTANCE ACROSS EACH UNPARALLELED TYPE 98 ELEMENT ABOVE +C <----><---->RRRRRRXXXXXXCCCCCC + AUT50AAUT501AUT5A6AUT5A5 + AUT50BAUT50AAUT5A6AUT5A5 + AUT50CAUT50BAUT5A6AUT5A5 + AUT50FAUT502AUT5A6AUT5A5 + AUT50GAUT50FAUT5A6AUT5A5 + AUT50HAUT50GAUT5A6AUT5A5 +$VINTAGE, 1 +C <----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P +C SHUNT RESISTANCES FOR EDDY LOSS MODELING OF 1/8 LAMINATION THICKNESS LAYERS. +C CROSS-FLUX CORE RESISTIVITY IS 12.00 E-8 OHM-METERS. +C LAMINATIONS ARE 0.2286 MM THICK; MODELED LAYERS ARE 0.0286 MM THICK. + AUT5A5 141050.02696 + AUT5A6 70525.01348 + AUT5A7 70525.01348 + AUT5A8 70525.01348 + AUT5B5AUT5B0 141050.02696 + AUT5B6AUT5B0 70525.01348 + AUT5B7AUT5B0 70525.01348 + AUT5B8AUT5B0 70525.01348 + AUT5C5AUT5C0 141050.02696 + AUT5C6AUT5C0 70525.01348 + AUT5C7AUT5C0 70525.01348 + AUT5C8AUT5C0 70525.01348 + AUT501AUT5B0 152575.09208 + AUT50AAUT5B0 76287.54604 + AUT50BAUT5B0 76287.54604 + AUT50CAUT5B0 76287.54604 + AUT502AUT5C0 152575.09208 + AUT50FAUT5C0 76287.54604 + AUT50GAUT5C0 76287.54604 + AUT50HAUT5C0 76287.54604 +C COMBINED WINDING 1-2 AND WINDING 2-3 STRAY (EDDY CUR.) LOSS MODEL: +C 3-PHASE W1-W2 STRAY LOSSES AT RATED MVA = 636.953 KW. +C WITH W1-W2 REF.VOLTS = 2.085 KV.; AND W2-W3 REF.VOLTS = 0.359 KV. + AUT5A1AUT5A2 20.47385 + AUT5B1AUT5B2 20.47385 + AUT5C1AUT5C2 20.47385 +C 3-PHASE W2-W3 STRAY LOSSES AT RATED MVA = 127.094 KW. +C WITH W1-W2 REF.VOLTS = -0.495 KV.; AND W2-W3 REF.VOLTS = 7.203 KV. + AUT5A2AUT5A3 1706.84712 + AUT5B2AUT5B3 1706.84712 + AUT5C2AUT5C3 1706.84712 +C TANK WALL MULTIPLE CHAINED PI MODEL RESISTANCES REFERRED TO 100.000 TURNS: +C MODELED SURFACE LAYERS ARE 0.1279 MM THICK; RHOTANK = 15.00 E-8 O-M. +C <----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P + AUT581 1.67957937 + AUT582 3.35915874 + AUT583 6.71831748 + AUT584 13.43663496 + AUT585 26.87326992 + AUT586 53.74653984 + AUT587 80.61980976 + AUT504 161.23961952 +$VINTAGE, 0 +C TRANSFORMER TANK WALL MODEL FOR TRANSFORMER AUT5 +C TANK WALL AND VERTICAL BRACING HAS A CROSS-SECTIONAL AREA OF 0.354 SQ. METERS. +C ITS INSIDE PERIMETER IS 21.65 METERS. ITS AVERAGE THICKNESS IS 0.016 METERS. +C EQUIVALENT HEIGHT AFFECTED BY ZERO SEQUENCE FLUX IS ASSUMED AT 3.15 METERS. +C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS: +C LINEAR:(0.20, 70.00) AND FULLY SATURATED:(2.05, 75000.) +C STEADY-STATE CALCS USE THE SLOPES AT POINTS 3 BELOW. (1.0000 TESLA) +C

(--FOR TYPE 98 N/L REACTOR--) P +C CURRENT ----->< FLUX --------> +98AUT581 21.31036.333 + 2.20500 3.54468 + 5.04000 10.63404 + 10.39500 17.72340 + 22.05000 23.04042 + 50.40000 26.58510 + 110.25000 28.35744 + 226.80000 30.12978 + 441.00000 31.90212 + 1134.00000 35.44680 + 2362.50000 36.33297 + 4725.00000 38.03714 + 9999 +98AUT582AUT581 21.31018.166 + 2.20500 1.77234 + 5.04000 5.31702 + 10.39500 8.86170 + 22.05000 11.52021 + 50.40000 13.29255 + 110.25000 14.17872 + 226.80000 15.06489 + 441.00000 15.95106 + 1134.00000 17.72340 + 2362.50000 18.16648 + 4725.00000 19.01857 + 9999 +98AUT583AUT582 21.3109.0832 + 2.20500 0.88617 + 5.04000 2.65851 + 10.39500 4.43085 + 22.05000 5.76011 + 50.40000 6.64628 + 110.25000 7.08936 + 226.80000 7.53245 + 441.00000 7.97553 + 1134.00000 8.86170 + 2362.50000 9.08324 + 4725.00000 9.50929 + 9999 +98AUT584AUT583 21.3104.5416 + 2.20500 0.44309 + 5.04000 1.32926 + 10.39500 2.21543 + 22.05000 2.88005 + 50.40000 3.32314 + 110.25000 3.54468 + 226.80000 3.76622 + 441.00000 3.98777 + 1134.00000 4.43085 + 2362.50000 4.54162 + 4725.00000 4.75464 + 9999 +98AUT585AUT584 21.3102.2708 + 2.20500 0.22154 + 5.04000 0.66463 + 10.39500 1.10771 + 22.05000 1.44003 + 50.40000 1.66157 + 110.25000 1.77234 + 226.80000 1.88311 + 441.00000 1.99388 + 1134.00000 2.21543 + 2362.50000 2.27081 + 4725.00000 2.37732 + 9999 +98AUT586AUT585 21.3101.1354 + 2.20500 0.11077 + 5.04000 0.33231 + 10.39500 0.55386 + 22.05000 0.72001 + 50.40000 0.83078 + 110.25000 0.88617 + 226.80000 0.94156 + 441.00000 0.99694 + 1134.00000 1.10771 + 2362.50000 1.13541 + 4725.00000 1.18866 + 9999 +98AUT587AUT586 21.310.56770 + 2.20500 0.05539 + 5.04000 0.16616 + 10.39500 0.27693 + 22.05000 0.36001 + 50.40000 0.41539 + 110.25000 0.44309 + 226.80000 0.47078 + 441.00000 0.49847 + 1134.00000 0.55386 + 2362.50000 0.56770 + 4725.00000 0.59433 + 9999 +98AUT504AUT587 21.310.56770 + 2.20500 0.05539 + 5.04000 0.16616 + 10.39500 0.27693 + 22.05000 0.36001 + 50.40000 0.41539 + 110.25000 0.44309 + 226.80000 0.47078 + 441.00000 0.49847 + 1134.00000 0.55386 + 2362.50000 0.56770 + 4725.00000 0.59433 + 9999 +C FLUX SHIELD N-L MODEL FOR TRANSFORMER: AUT5 +C FLUX SHIELDS ARE 3.15 METERS LONG WITH AREA OF 0.3652 SQUARE METERS. +C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS: +C LINEAR:(0.20, 4.00) AND FULLY SATURATED:(2.00, 15000.00) +C FOR SS CALCS EACH F.S. HAS REACT. OF 9505. OHMS (AT 100.0 TURNS) +C BASED ON SLOPE OF POINT 8 BELOW, HAVING A FLUX DENSITY OF 1.7400 TESLA. +C

(--FOR TYPE 98 N/L REACTOR--) P +C CURRENT -----> +98AUT503AUT504 2.896673.034 + 0.12600 7.30340 + 0.20475 14.60680 + 0.30870 32.86530 + 0.39690 43.82040 + 0.50400 51.12380 + 0.94500 54.77550 + 1.73250 59.52271 + 2.52000 63.53958 + 6.30000 67.55645 + 15.75000 69.38230 + 47.25000 71.20815 + 472.50000 73.03400 + 945.00000 75.06272 + 9999 +C RESISTANCE UNDERESTIMATING FLUX SHIELD LOSSES, EQUAL TO LINEAR REACTANCE. +C <----><---->RRRRRRXXXXXXCCCCCC P + AUT503AUT504 9505.5 +C TANK TOP MULTIPLE CHAINED PI MODEL RESISTANCES REFERRED TO 100.000 TURNS: +C MODELED SURFACE LAYERS ARE 0.1172 MM THICK; RHOTANK = 15.00 E-8 O-M. +$VINTAGE, 1 +C <----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P + AUT591 2.00000000 + AUT592 4.00000000 + AUT593 8.00000000 + AUT594 16.00000000 + AUT595 32.00000000 + AUT596 64.00000000 + AUT597 96.00000000 + AUT509 192.00000000 +$VINTAGE, 0 +C TRANSFORMER TANK TOP MODEL FOR TRANSFORMER AUT5 +C TANK TOP EDDY CURRENT PATH IS MODELED WITH A CIRCUMFERENCE OF 15.00 METERS. +C AND A WIDTH OF 2.00 METERS; THICKNESS IS 0.015 METERS. +C THE B-H CURVE USED PASSES THROUGH THE FOLLOWING (TESLA, AMP/METER) POINTS: +C LINEAR:(0.20, 70.00) AND FULLY SATURATED:(2.05, 75000.) +C STEADY-STATE CALCS USE THE SLOPES AT POINTS 3 BELOW. +C

(--FOR TYPE 98 N/L REACTOR--) P +C CURRENT ----->< FLUX --------> +98AUT591 13.53023.062 + 1.40000 2.25000 + 3.20000 6.75000 + 6.60000 11.25000 + 14.00000 14.62500 + 32.00000 16.87500 + 70.00000 18.00000 + 144.00000 19.12500 + 280.00000 20.25000 + 720.00000 22.50000 + 1500.00000 23.06250 + 3000.00000 24.14423 + 9999 +98AUT592AUT591 13.53011.531 + 1.40000 1.12500 + 3.20000 3.37500 + 6.60000 5.62500 + 14.00000 7.31250 + 32.00000 8.43750 + 70.00000 9.00000 + 144.00000 9.56250 + 280.00000 10.12500 + 720.00000 11.25000 + 1500.00000 11.53125 + 3000.00000 12.07212 + 9999 +98AUT593AUT592 13.5305.7656 + 1.40000 0.56250 + 3.20000 1.68750 + 6.60000 2.81250 + 14.00000 3.65625 + 32.00000 4.21875 + 70.00000 4.50000 + 144.00000 4.78125 + 280.00000 5.06250 + 720.00000 5.62500 + 1500.00000 5.76562 + 3000.00000 6.03606 + 9999 +98AUT594AUT593 13.5302.8828 + 1.40000 0.28125 + 3.20000 0.84375 + 6.60000 1.40625 + 14.00000 1.82812 + 32.00000 2.10938 + 70.00000 2.25000 + 144.00000 2.39062 + 280.00000 2.53125 + 720.00000 2.81250 + 1500.00000 2.88281 + 3000.00000 3.01803 + 9999 +98AUT595AUT594 13.5301.4414 + 1.40000 0.14063 + 3.20000 0.42187 + 6.60000 0.70312 + 14.00000 0.91406 + 32.00000 1.05469 + 70.00000 1.12500 + 144.00000 1.19531 + 280.00000 1.26563 + 720.00000 1.40625 + 1500.00000 1.44141 + 3000.00000 1.50901 + 9999 +98AUT596AUT595 13.530.72070 + 1.40000 0.07031 + 3.20000 0.21094 + 6.60000 0.35156 + 14.00000 0.45703 + 32.00000 0.52734 + 70.00000 0.56250 + 144.00000 0.59766 + 280.00000 0.63281 + 720.00000 0.70312 + 1500.00000 0.72070 + 3000.00000 0.75451 + 9999 +98AUT597AUT596 13.530.36035 + 1.40000 0.03516 + 3.20000 0.10547 + 6.60000 0.17578 + 14.00000 0.22852 + 32.00000 0.26367 + 70.00000 0.28125 + 144.00000 0.29883 + 280.00000 0.31641 + 720.00000 0.35156 + 1500.00000 0.36035 + 3000.00000 0.37725 + 9999 +98AUT509AUT597 13.530.36035 + 1.40000 0.03516 + 3.20000 0.10547 + 6.60000 0.17578 + 14.00000 0.22852 + 32.00000 0.26367 + 70.00000 0.28125 + 144.00000 0.29883 + 280.00000 0.31641 + 720.00000 0.35156 + 1500.00000 0.36035 + 3000.00000 0.37725 + 9999 +C ZERO SEQUENCE MODELS ARE BASED ON AIR-CORE MODELS, PLUS DERIVED VALUE FOR +C X100T OF 3.29690 OHMS/PHASE TANK EFFECT, REFERRED TO 100 TURNS. +C THE (UNSHIELDED) TANK TOP IS MODELED AND HAS STRAY LOSSES. +C THE TANK WALL IS SATURABLE; BUT ANY STRAY LOSSES ARE NOT MEASURABLE. +C THE FLUX SHIELD IS SATURABLE AND LACKS SIGNIFICANT DISTRIBUTED GAPS. +C COUPLING REPRESENTING TOP AND SIDE GAP REACTANCES: +C PER-PHASE SIDE GAPS = 1.846 OHMS; TOP GAPS = 1.814 OHMS AT 100 TURNS. +51AUT503 0.000 10.978117 +52AUT509 0.000 5.4405450.000 5.440545 +C SHUNT FLUX PATHS BYPASSING LEGS AND COILS OF EACH PHASE FOLLOW: +C <----><---->RRRRRRXXXXXXXXXXXX +51AUT5A4AUT5B0 0.000 7.463166 +51AUT5B4AUT5C0 0.000 7.463166 +51AUT5C4AUT503 0.000 7.463166 +C ZERO-SEQUENCE FLUXES FROM EACH PHASE TO TANK ARE EQUALIZED BY THE FOLLOWING: +C <----><---->RRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXX +51AUT503AUT502 0.000 3296.896640 +52AUT502AUT501 0.000 3296.8933430.000 3296.896640 +53AUT501 0.000 3296.8933430.000 3296.8933430.000 3296.896640 +C NEXT COUPLING SUMS ZERO-SEQUENCE VOLTAGES OF OUTER WINDINGS AT NODE (0,8): +C <----><---->RRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXX +51AUT5A3 0.000 354922.42080 +52AUT506 0.000 354922.065870.000 354922.42080 +51AUT5B3AUT5B0 0.000 354922.42080 +52AUT507AUT506 0.000 354922.065870.000 354922.42080 +51AUT5C3AUT5C0 0.000 354922.42080 +52AUT508AUT507 0.000 354922.065870.000 354922.42080 +C NON-DEDICATED STRAY ZERO-SEQUENCE LOSS RESISTANCE: +C PRODUCES 158.213 KW LOSS FOR ENTERED TEST CONDITIONS. +$VINTAGE, 1 +C <----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P + AUT508 354.92242080 +$VINTAGE, 0 +C TIME-STEP-DAMPING BRANCHES ACROSS MAJOR AIR GAPS FOR HARMONIC ORDER: 500.0 +C <----><---->RRRRRRXXXXXXCCCCCC + AUT5A4AUT5B0 373.16 .00142 + AUT5B4AUT5C0 373.16 .00142 + AUT5C4AUT503 373.16 .00142 + AUT501 164.84 .00322 + AUT501AUT502 164.84 .00322 + AUT502AUT503 164.84 .00322 +C 10000-HENRY FLUX MONITORING INDUCTANCES TO INITIALIZE TACS. +C <----><---->RRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXXRRRRRRXXXXXXXXXXXX +51AUT5A5AUT5D1 3769911.185 +51AUT5B5AUT5D2 3769911.185 +51AUT5C5AUT5D3 3769911.185 +51AUT501AUT5D4 3769911.185 +51AUT502AUT5D5 3769911.185 +C End of $INCLUDE. File name = aut5-brn.inc +$VINTAGE, 1 +C <----><---->RRRRRRRRRRRRRRRRXXXXXXXXXXXXXXXXCCCCCCCCCCCCCCCC P +C CORNER GROUND DELTA + Y1 1000. +C ENERGIZE VIA 100 OHM RESISTANCES FOR DAMPING +C 0.1 ohm is now essential; 100 is too large relative to xl + SRC2A EGCT1A .1 + SRC2B EGCT1B .1 + SRC2C EGCT1C .1 +$VINTAGE, 0 +BLANK CARD ENDING ALL ATP-SORTED BRANCH CARDS (from blank.inc) +C $INCLUDE aut5-swx.inc +C WSM touches using DOS editor to make readable. +C MEAS. SWITCHES TO MONITOR CURRENT IN 10000-HENRY CORE FLUX-MONITORING INDUCTS. +C 67890123 MEASURING 56789012345678 P + AUT5D1 MEASURING + AUT5D2AUT5B0 MEASURING + AUT5D3AUT5C0 MEASURING + AUT5D4AUT5B0 MEASURING + AUT5D5AUT5C0 MEASURING +C End of $INCLUDE. File name = aut5-swx.inc +C ********************************************* STUDY-DEPENDENT CABLE SWITCHING: +C SWITCH CARD: COL 1-2 IS 0 FOR ORDINARY & GAUSS. STATISTICS SWITCHES, NOT TACS. +C CLOSING, OPENING TIMES AND STATISTICS PARAMETERS BELOW ARE IN SECONDS. +C 67890123 MEASURING 56789012345678 P +C STATISTICS P +C METERS INTERNAL CURRENT OF PHASE LEG A,B,C + AUT5AMAUT5A5 MEASURING 1 + AUT5BMAUT5B5 MEASURING 1 + AUT5CMAUT5C5 MEASURING 1 +C ENERGIZED AT ZERO + SRCA SRC2A MEASURING 1 + SRCB SRC2B MEASURING 1 + SRCC SRC2C MEASURING 1 +C ENERGIZED AT ZERO PLUS +C SRCA SRC2A .0001 999. 1. +C SRCB SRC2B .0001 999. 1. +C SRCC SRC2C .0001 999. 1. +BLANK CARD ENDING ALL ATP-SORTED SWITCH CARDS (from blank.inc) +$LISTON +C $INCLUDE aut5-src.inc +C WSM touches file with DOS editor to make it into a DOS file +C TACS-CONTROLLED SOURCES FOR HYSTERESIS/RESIDUAL MMF IN TRANSFORMER AUT5 +C *= VOLTAGE IF POSITIVE; CURRENT IF NEGATIVE +60AUT5A5-1 +60AUT5B5-1 +60AUT5C5-1 +60AUT5B0-1 +60AUT5C0-1 +60AUT501-1 +60AUT502-1 +C End of $INCLUDE. File name = aut5-src.inc +C *= VOLTAGE IF POSITIVE; CURRENT IF NEGATIVE +C +C 138 KV L-L = 112676.5 V CREST or 1877.94 at 1 hertz +C 40 percent = 751.17 +C 60 percent = 1126.76 +C 80 percent = 1502.35 +C 100 percent = 1877.94 +C 120 percent = 2253.53 +14SRCA 1 2253.53 1. 0. -1.0 999. +14SRCB 1 2253.53 1. -120. -1.0 999. +14SRCC 1 2253.53 1. -240. -1.0 999. +BLANK CARD ENDING ALL ATP-SORTED SOURCE CARDS (from blank.inc) +C < >< >< >< >< >< >< >< > + EGCT1A EGCT1B EGCT1C +BLANK CARD ENDING ALL ATP-SORTED OUTPUT CARDS (from blank.inc) + 193.01 0.0 0.1 BRANCH + TACS ALLINATACS ALLINBTACS ALLINC +BEGIN NEW DATA CASE +C 2nd of 4 subcases. +C "FREQUENCY SCAN" use, with subsequent plotting of phasors vs. frequency. +C This test case uses uniform frequency spacing. It basically like DC-51, +C but involves some variations of output and batch-mode plotting. +PRINTED NUMBER WIDTH, 8, 2, { Each column of width 8 includes 2 blank bytes +FREQUENCY SCAN, 60., 10., 280., 0, { 60 < f < 280 Hz in 10-Hz increments +CHANGE PLOT FREQUENCY { Request to vary frequency of plot points and .PL4 file + 5 5 10 3 20 1 { On step 5, change to 5, etc. + 0.1 0.0 { dT of col 1-8 is not used; TMAX of 9-16 ==> no simulation + 1 1 1 0 1 -1 + 5 5 10 3 20 1 +C Note: 1-punch in col. 80 of the following card was replaced by LIN002 below. +C This illustrates a second way to order current output: by branch name. +C When the user does not name branches himself, NMAUTO = 1 means the +C program will provide default names. The first branch is LIN001, so +C this should be used to replace the column-80 punch. But this branch +C is in series with the 2nd branch, and has no shunt leakage, so the +C same current can be found from the second branch (name LIN002). This +C alternative is more interesting since it shows current output seems to +C be possible for the 1-phase, constant-parameter distributed model. + GEN SEND 3.0 40. +-1SEND MID .306 5.82 .012 20. { 20 miles of line from DC-37 + 1MID REC 34.372457.68.15781 { 10 miles of Pi from DC-3 + EARTH 200. 1 +C The preceding branch has a 1-punch in column 80 for current output. If a 3 +C had instead been used, branch voltage would have been produced, too. Note +C it would be equal to the node voltage of node EARTH. No choice of polarity +C is possible, this way. More interesting is the "-5" request below, which +C allows polarity reversal (a minus sign will be added). +BLANK card ending branches +C Note: 1-punch in col. 80 of the following card was replaced by SWT001 below. +C This illustrates a second way to order current output: by switch name. + REC EARTH MEASURING +BLANK card ending (here, non-existent) switches +14GEN 100. 60. -1. +BLANK card ending sources +C Requests for output variables follow. About currents, there is no control +C over location in the output vector. Switch currents always precede branch +C currents, and both are in natural order (of component number). But order of +C voltages is controlled by the user. Note requests for node voltages have +C been split, and in between there is a request for a voltage difference. +C This is the order in the output vector, note. Yet, all voltage requests +C made here (after the blank card ending source cards) follow all that might +C have been made with column-80 punches on branch or switch cards. + GEN SEND MID REC +-5 EARTH { -5 ==> 2A6 name pairs for voltage differences (branch V) +-1SWT001LIN002 { -1 ==> Branch/switch current out; use A6 component names + EARTH { Node voltage. Note branch voltage 2 lines above has opposite polarity +BLANK card ending node voltage outputs + PRINTER PLOT +C Note the peculiar 60 units/inch on the following plot cards. This is not +C a mistake. With 6 lines/inch of the printer plot, & the F-scan frequency +C increment of 10 Hz, there will be exactly one line for each freq step: + 14660. 60.280. GEN SEND MID REC + 19660. 60.280. SEND MID EARTH REC +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 3rd of 4 subcases illustrates IF-THEN-ENDIF for conditional data +C assembly as described in the July, 1998, newsletter. Environment +C variable RESISTOR must be set prior to use, in order. to avoid +C an error stop. This can not be set internally, unfortunately. +PRINTED NUMBER WIDTH, 9, 2, +POWER FREQUENCY, 50., + .0001 .020 50. + 1 1 1 1 1 -1 + 5 5 20 20 +C 1st of 2 identical, disconnected networks uses manually-defined branches: +C IF ( RESISTOR .EQ. 'ONEHALF' ) SET VARIABLE { Execute OS command SET RESISTOR=ONEHALF +IF ( RESISTOR .EQ. 'ONEHALF' ) THEN { If user wants to use resistive alternative: + GEN TRAN 0.5 { 1st of 3 alternatives is resistor +ELSEIF ( RESISTOR .EQ. 'ZERO' ) THEN { If user wants to use resistive alternative: + GEN TRAN 0.0 { 2nd of 3 alternatives is illegal resistor +ELSE { Alternatively (if environment variable RESISTOR is not YES): + GEN TRAN 0.5 { 3rd of 3 alternatives is inductor +ENDIF { Terminate 5-line IF statement that uses RESISTOR to choose the model + TRAN 0.5 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 1.414 50. 0.0 -1. +BLANK card terminating program source cards. + GEN TRAN +BLANK card ending program output-variable requests. +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C BENCHMARK DC-4 +C 4th of 4 subcases is modification of DC-4 as used by SPY @5. But that +C required SPY. Beginning 23 Aug 1998, parametric studies can be run +C using higher-level batch-mode commands. See October, 1998, newsletter. +C Following request carries params MAXKNT IOPCVP +POCKET CALCULATOR VARIES PARAMETERS 5 1 { Loop five times +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .020 12.0 + 1 1 0 0 1 -1 1 + 5 5 20 20 100 100 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would be +MILLIH = 50000. - ( KNT - 1.0 ) * 5000. { Gives L =50, 45, 40, 35, & 30K mHenry +BLANK card ends $PARAMETER definitions that are processed just b4 branch cards + BRANCH NAME:First { Even though name could go on next card, use this instead + GEN TRAN 5.0MILLIH { $PARAMETER will define this inductance + TRAN NAME R-mag 1.E4 +93TRAN NAME Magnet .005 30. 3 + 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 + LOADG 1.E-6 +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. + GEN TRAN +BLANK card ending program output-variable requests. +BLANK card ending all plot cards { In fact, ignored until the final shot +BEGIN NEW DATA CASE +BLANK -- cgit v1.2.3