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/dcn26.dat | 1334 ++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 1334 insertions(+) create mode 100644 benchmarks/dcn26.dat (limited to 'benchmarks/dcn26.dat') diff --git a/benchmarks/dcn26.dat b/benchmarks/dcn26.dat new file mode 100644 index 0000000..fdf12e5 --- /dev/null +++ b/benchmarks/dcn26.dat @@ -0,0 +1,1334 @@ +BEGIN NEW DATA CASE +C 1st of 19 subcases that illustrate POCKET CALCULATOR VARIES PARAMETERS (the +C same as introduced at the end of DCNEW-25) without simulation. The previous +C illustration involved simulation. Here, TMAX < 0 ===> none. Although +C one could save the .PL4 file, there is nothing new to illustrate in this +C usage, so it will not be done. Variable IOPCVP controls output for the +C 2nd or later pass of loop. The zero illustrated here minimizes output. If +C it is toggled to value 1, all output of the 1st subcase should be repeated +C during the 2nd or later pass. This is useful for debugging. +C Controls of following request card: MAXKNT IOPCVP { Loop 8 times with +POCKET CALCULATOR VARIES PARAMETERS 8 1 { minimum printout +PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) + .0001 -1. 50. 50. + 1 1 1 0 1 2 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b +MHENRY = 1.2 - ( KNT - 1.0 ) * .05 { L =1.2, 1.15, 1.1, etc. +FREQUENCY = 25. * KNT { Here we vary the frequency for use on source card +BLANK card ends $PARAMETER definitions that are processed just b4 branch card + TRAN 1.0 1.E+5 3 + GEN TRAN 1.0MHENRY +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 70. FREQUENCY -1. +BLANK card terminating program source cards. + GEN TRAN +BLANK card ending program output-variable requests. + PRINTER PLOT + 19660. 0.0200. TRAN +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 2nd of 19 subcases that illustrate PCVP will superimpose the PCVP loop on +C top of a normal FREQUENCY SCAN loop. In fact, the latter controls the +C loop, so MAXKNT is arbitrary and can be ignored. By design, the +C solution to this 2nd subcase is the same as the first. As for subcase one, +C the .PL4 file will not be saved, although it could be. +FREQUENCY SCAN, 25., 25., 200., 0, { 25 < f < 200 Hz in 25-Hz increments +C Controls of following request card: MAXKNT IOPCVP { Ignore MAXKNT value +POCKET CALCULATOR VARIES PARAMETERS 0 1 { since also FS loop +PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) + .0001 -1. 50. 50. + 1 1 1 0 1 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b +MHENRY = 1.2 - ( KNT - 1.0 ) * .05 { L =1.2, 1.15, 1.1, .... +BLANK card ends $PARAMETER definitions that are processed just b4 branch card + TRAN 1.0 1.E+5 3 + GEN TRAN 1.0MHENRY +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 70. 25. -1. +BLANK card terminating program source cards. + GEN TRAN +BLANK card ending program output-variable requests. + PRINTER PLOT + 19660. 0.0200. TRAN +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 3rd of 19 subcases that illustrate PCVP will superimpose the PCVP loop on +C top of a normal HARMONIC FREQUENCY SCAN loop. In fact, the latter controls +C the loop, so MAXKNT is arbitrary and can be ignored. By design, +C solution to this 3rd subcase is the same as the preceding two. This is an +C abnormal, special case (normally, HFS data will skip harmonics, and source +C amplitudes will be varied). +C 27 December 2000, MODELS data is added. Nothing else is changed. This +C use of MODELS is imagined as an alternative to $POCKET as illustrated +C in the following (4th) subcase. Automatically, MODELS is called at the +C end of the phasor solution so the user need do nothing more than add his +C MODELS data in the usual place. 28 December 2000, add a switch, too, to +C prove that phasor switch current also can be viewed using MODELS. +POWER FREQUENCY, 25., ! Needed so minimum frequency is recognized as fundamental +HARMONIC FREQUENCY SCAN { Non-negative DELFFS in 25-32 means F in Hz (not log F) +C Controls of following request card: MAXKNT IOPCVP { Ignore MAXKNT value +POCKET CALCULATOR VARIES PARAMETERS 0 1 { since also HFS loop +PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) + .0001 -1. 50. 50. + 1 1 1 0 1 +MODELS +INPUT v_Re_tran {v(tran)} -- real part of voltage at node TRAN + v_Im_tran {imssv(tran)} -- imag part of voltage at node TRAN + i_Re_swit {i(swit)} -- real part of current at switch SWIT + i_Im_swit {imssi(swit)} -- imag part of current at switch SWIT +MODEL m1 ----------------------------------------------------------------------- + VAR pass, ampl, angle, sum -- list of all local variables + INPUT v_Re_1 {dflt: 0} -- to verify that real part of node voltage and the + v_Im_1 {dflt: 0} -- imaginary part of node voltage are accessed correctly + i_Re_1 {dflt: 0} -- to verify the real part of switch current + i_Im_1 {dflt: 0} -- to verify the imaginary part of switch current + INIT pass:=0 + sum:=0 ENDINIT + EXEC + pass:=pass+1 + ampl:= sqrt( v_Re_1 **2 + v_Im_1 ** 2 ) -- Magnitude of the node voltage + angle:= atan2( v_Im_1 , v_Re_1 ) -- Angle in radians of the node voltage + sum:= sum + ampl -- Sum the voltage amplitudes for different harmonics +write(' ************** In MODELS, pass = ', pass, ' freq = ', t ) +write(' Re{i1}, Im{i1} = ', i_Re_1, ', ', i_Im_1 ) +write(' Re{v1}, Im{v1} = ', v_Re_1, ', ', v_Im_1 ) +write(' Polar =', ampl, angle ) +write(' Sum of amplitudes =', sum ) + ENDEXEC +ENDMODEL +USE m1 as m1 + INPUT v_Re_1 := v_Re_tran + v_Im_1 := v_Im_tran + i_Re_1 := i_Re_swit + i_Im_1 := i_Im_swit +ENDUSE +ENDMODELS +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b +MHENRY = 1.2 - ( KNT - 1.0 ) * .05 { L =1.2, 1.15, 1.1, .... +BLANK card ends $PARAMETER definitions that are processed just b4 branch card + TRAN SWIT 1.0 1.E+5 3 + GEN TRAN 1.0MHENRY +BLANK card ending program branch cards. + SWIT -1. 1.0 1 +BLANK card terminating program switch cards (none, for this case) +14GEN 70. 25. 0.0 { Note comment and no negative T-start +14GEN 70. 50. 0.0 { This is normal HFS source data. It +14GEN 70. 75. 0.0 { is unusual, however, because of our +14GEN 70. 100. 0.0 { match preceding subases. Unusual +14GEN 70. 125. 0.0 { for HFS usage, the amplitude is not +14GEN 70. 150. 0.0 { varied and all harmonics are present +14GEN 70. 175. 0.0 +14GEN 70. 200. 0.0 +BLANK card ending source cards +BLANK card ending F-dependent series R-L-C branches (none, for this subcase) + GEN TRAN +BLANK card ending program output-variable requests. + CALCOMP PLOT { Use BOTH rather than DISK to see a bar chart of harmonics + 19680. 0.0200. TRAN SWIT { Units/inch of 5-7 are ingored since bar chart +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 4th of 19 subcases that illustrate PCVP comes from Gabor Furst. As can be +C read in the January, 1999, newsletter, this is "a semi-realistic HFS file +C .... This resistance is varied by the Funk-Hantel formula." Note harmonic +C source amplitudes vary (unlike 3rd subcase), and some harmonics are missing. +C 27 December 2000, modify existing case by addition of SAVE and $POCKET +C to illustrate operation of these new features. Included is WRITE( for +C output that is part of regular ATP output. +POWER FREQUENCY, 50., ! Needed so minimum frequency is recognized as fundamental +HARMONIC FREQUENCY SCAN { Non-negative DELFFS in 25-32 means F in Hz (not log F) +C Controls of following request card: MAXKNT IOPCVP { Ignore MAXKNT value +POCKET CALCULATOR VARIES PARAMETERS 0 1 { since also HFS loop +PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) + .0001 -1. 50. + 1 1 1 0 1 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would be +C SAVE ALL { Until selective SAVEs work better, make all variables non-volatile +C The preceding should be necessary. Instead, it results in KILL +C code because following resistance is not encoded right. So, +C remove temporarily. Amazingly, variable SUM seems believable. +RESISTANCE = 2.0 * (1.0 + 0.2 * (KNT - 1.0) ** 1.5)) { Funk-Hantel freq. dependence +BLANK card ends $PARAMETER definitions that are processed just b4 branch cards +$VINTAGE, 1, { Wide: RRRRRRRRRRRRRRRRLLL .... + TRAN RESISTANCE 3 +$VINTAGE, 0, { Narrow: RRRRRRLLLLLLCCCCCC + GEN TRAN 2.0 6.1 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) + POLAR OUTPUT VARIABLES { 2nd of 3 alternatives gives mag, angle (not mag only) +14GEN 10200. 1. 0.0 { Power freq source has harmonic # 1 +14GEN 500. 3. 0.0 { 3rd harmonic source has harmonic # 3 +14GEN 3000. 5. 0.0 +14GEN 2000. 7. 0.0 +14GEN 600. 11. 0.0 +BLANK card ending source cards +BLANK card ending F-dependent series R-L-C branches (none, for this subcase) + GEN TRAN +$POCKET { This will be serviced by CIMAGE just as any other $-card would be +C This is a comment within $POCKET. It shows they are ignored. Note placement +C of this block as part of node-voltage outputs, which ensures that it will be +C processed after each phasor solution. If one wanted to perform mathematics +C on the solution, this would be the appropriate location. About summation, +C values should be as follows: Harmonic # KNT: 1 3 5 7 11 +C Summation: 1 4 9 16 27 +STEP = KNT { Output is allowed only for local variables, so copy KNT to local +SUM = SUM + KNT { Preceding SAVE is for SUM, to accumulate sum of loop numbers +WRITE( ' -----------------------' ) { Statement can not go beyond column 40 +WRITE( ' >> Step KNT =', STEP ) $$ { The double dollar => continue on next +WRITE( ' Value of SUM =', SUM ) { No double dollar ===> end of out line +WRITE( ' =======================' ) { This separator line goes on new line +C There is a limit on the $$ concatensation. ATP output line length is +C 132 bytes maximum. It is the user's responsibility not to exceed this. +BLANK card ends $POCKET definitions to be processed after each phasor solution +BLANK card ending program output-variable requests. + CALCOMP PLOT { Use BOTH rather than DISK to see a bar chart of harmonics + 19680. 0.0600. TRAN { Units/inch of 5-7 are ingored since bar chart +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 5th of 19 subcases is unrelated to the preceding. Rather, it illustrates the +C first use of superposition of phasor solutions. Success occurred 8 Nov 98. +C This is a brand new type of frequency scan. It is fundamentally different +C because of the superposition and continuation into the time-step loop. In +C the following, only the minimum of 2 phasor solutions are superimposed. The +C more general case will involve repeated sets of excitation. Each such set +C begins with the request word SUPERIMPOSE MORE SOURCES: and ends with the +C blank card ending sources. The final set of excitation is the one actually +C used in the time-step loop, so must include all sources of the preceding +C sets, if this is what the user wants. Note T-start is not negative for +C such additions, however. Initial restriction: series R-L-C branches only. +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .01 2.0 + 1 1 1 1 1 -1 + 5 5 20 20 + GEN TRAN 1.0 10. 3 + TRAN 0.20E6 1 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 1.0 1.0 -1. +BLANK card terminating source cards for 1st phasor solution. + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GEN 1.0 1.0 { Source of 1st solution represented in dT-loop +14GEN -1.0 2.0 { 2nd harmonic source } -1. +BLANK card terminating source cards for 2nd (and final) phasor solution. + GEN TRAN +C First 3 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 Step Time GEN GEN TRAN GEN TRAN +C TRAN TRAN TERRA +C 0 0.0 -.278585723 0.0 .2785857231 -.280475835 -.280475835 +C 1 .01 -.258845905 .0059120271 .2647579317 -.272635821 -.272635821 +C 2 .02 -.228080808 .0235315402 .2516123484 -.253187511 -.253187511 +C 3 .03 -.187201315 .0525107648 .2397120803 -.222823212 -.222823212 +BLANK card ending program output-variable requests. +C 100 1.0 -.278768592 0.0 .2787685919 -.280624596 -.280624596 +C 120 1.2 .6447560418 1.118033989 .473277947 .6497741745 .6497741745 +C 140 1.4 -1.4363837 -1.11803399 .3183497136 -1.33496603 -1.33496603 +C 160 1.6 -.123925906 -1.11803399 -.994108083 -.293493889 -.293493889 +C 180 1.8 1.194323606 1.118033989 -.076289617 1.259311865 1.259311865 +C 200 2.0 -.278769536 0.0 .2787695359 -.280625592 -.280625592 +C Variable max: 1.417160223 1.124996567 .6739357604 1.408960752 1.408960752 +C Times of max: .75 .79 1.3 .76 .76 +C Variable min: -1.69650701 -2. -1.00648406 -1.68753921 -1.68753921 +C Times of min: 1.46 0.5 .62 1.47 1.47 + CALCOMP PLOT + 143 .2 0.0 2.0 GEN TRAN { Screen plot is perfectly smooth & periodic +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C 6th of 19 subcases extends uncoupled R-L-C elements of the 5th to multiphase +C lumped elements. Take a real 10-mile, 3-phase Pi-circuit from DC-3. But we +C want a longer line, so increase frequency to 400 Hz. Of course, single phase +C becomes 3-phase. Line section is terminated by balanced R-L load to ground. +C This solution was verified by comparing with one that omitted the second +C harmonic during the phasor solution, and extended T-max to 20 msec. The +C final 5 msec of voltages at TRANA, TRANB, and TRANC are very close to the +C solution about to be demonstrated. Also, in output below, compare step 0 +C with step 100. They are very close. But step 200 is even closer to 100 (a +C little settling has occurred; this is discretization error). +PRINTED NUMBER WIDTH, 10, 2, +.0000250 .005 3000. + 1 1 1 1 1 -1 + 5 5 20 20 + 1GENA TRANA 34.372457.68.15781 + 2GENB TRANB 35.735164.43-.031538.002451.79.16587 + 3GENC TRANC 35.735164.43-.031537.455151.72-.021938.002451.79.16587 + TRANA 20. 100. 1 + TRANB TRANA 1 + TRANC TRANA 1 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GENA 1.0 400. 0.0 -1. +14GENB 1.0 400. -120. -1. +14GENC 1.0 400. +120. -1. +BLANK card terminating source cards for 1st phasor solution. +C .... or "Q". The first solution frequency = 4.00000000E+02 Hertz. +C GENA 1.0 1.0 .00651879225127 .01786478870371 .00325939612564 -.3270376674E-4 +C 0.0 0.0 -.0166329799799 -68.5988095 .00831648998997 0.0061218 +C Total network loss P-loss by summing injections = 9.742991063317E-03 + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GENA 1.0 400. 0.0 +14GENB 1.0 400. -120. +14GENC 1.0 400. +120. +14GENA 1.0 800. 0.0 -1. +14GENB 1.0 800. -120. -1. +14GENC 1.0 800. +120. -1. +C .... or "Q". The first solution frequency = 8.00000000E+02 Hertz. +C GENA 1.0 1.0 .00179922688615 .00899257899614 .89961344308E-3 -.1832029691E-4 +C 0.0 0.0 -.008810746825 -78.4584168 .00440537341249 0.0031815 +C Total network loss P-loss by summing injections = 2.746542759765E-03 +BLANK card terminating source cards for 2nd (and final) phasor solution. + GENA TRANA GENB TRANB GENC TRANC +C First 6 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 3 output variables are branch currents (flowing from the upper node to the lower node); +C Step Time GENA TRANA GENB TRANB GENC TRANC TRANA TRANB TRANC +C TERRA TERRA TERRA +C 0 0.0 2.0 .6411782 -1. -.655262 -1. .0054071 .0081912 -.026551 .0176749 +C 1 .25E-4 1.990141 .6719155 -.832151 -.619565 -1.15799 -.059291 .0104086 -.02703 .0159628 +C 2 .5E-4 1.960698 .6964416 -.656435 -.578549 -1.30426 -.123014 .0125508 -.027293 .0141158 +C 3 .75E-4 1.912064 .7145235 -.47495 -.532698 -1.43711 -.18506 .014596 -.027337 .0121521 +C 4 .1E-3 1.84489 .7260175 -.289862 -.482595 -1.55503 -.244721 .0165237 -.027161 .0100912 +C 5 .125E-3 1.760074 .7308739 -.103383 -.428865 -1.65669 -.301342 .0183147 -.026767 .0079542 +C 10 .25E-3 1.118034 .6596136 .7736591 -.128596 -1.89169 -.520814 .0246777 -.021779 -.003099 +C 15 .375E-3 .2787683 .4603086 1.384884 .1607557 -1.66365 -.603533 .0259314 -.013004 -.012782 +C 20 .5E-3 -.5 .2026547 1.582676 .3644382 -1.08268 -.546164 .0222562 -.002864 -.018931 +BLANK card ending program output-variable requests. +C 40 1.E-3 -.5 -.185466 -.064602 .1746455 .5646021 .0137459 -.003673 .0129232 -.008711 +C 60 .0015 -.5 -.212516 .5646021 .2936534 -.064602 -.082817 -.005699 .0101002 -.004403 +C 80 .002 -.5 -.445559 -1.08268 -.178062 1.582676 .6101375 -.021055 .0063503 .0143909 +C 100 .0025 2.0 .6410699 -1. -.655039 -1. .0052891 .0081836 -.026535 .0176671 +C 120 .003 -.5 .2025812 1.582676 .3645686 -1.08268 -.546224 .0222515 -.002854 -.018936 +C 140 .0035 -.5 -.185508 -.064602 .1747437 .5646021 .0136866 -.003676 .0129291 -.008714 +C 160 .004 -.5 -.212542 .5646021 .2937035 -.064602 -.082843 -.005701 .0101038 -.004405 +C 180 .0045 -.5 -.445576 -1.08268 -.17803 1.582676 .6101221 -.021056 .0063525 .0143897 +C 200 .005 2.0 .6410606 -1. -.655017 -1. .005275 .008183 -.026534 .0176664 +C Variable maxima : 2.0 .7308739 1.582676 .443335 1.582676 .660407 .0261033 .013903 .0249648 +C Times of maxima : 0.0 .125E-3 .003 .00315 .002 .0021 .35E-3 .0034 .002275 +C Variable minima : -1.125 -.540524 -1.89169 -.742149 -1.89169 -.603595 -.021087 -.027337 -.020627 +C Times of minima : .004275 .00435 .00225 .002375 .25E-3 .002875 .004525 .75E-4 .0031 + CALCOMP PLOT + 144 .5 0.0 5.0 GENA GENB GENC + 144 .5 0.0 5.0 TRANA TRANB TRANC +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C 7th of 19 subcases demonstrates proper working of superposition of phasor +C solutions for a constant-parameter, distributed line. This first was made +C to work around sunset on 15 Nov 1998. Two hours later, the switch was made +C to work. It had a bug of its own. The switch is permanently closed, and +C does not change the answer. Its addition merely demonstrates that switches +C now are tolerated. Still no frequency-dependent line, however. That will +C be left for the expert (BPA's Dr. Tsu-huei Liu), if she thinks it is doable. +POWER FREQUENCY, 50., { Unnecessary. This merely avoids warning about 50 Hz. +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) + .000100 .040 50. 50. { Note L in ohms and C in micromhos at 50 Hz + 1 1 1 1 1 -1 + 2 -1 5 5 20 20 +-1SEND REC 0.3 0.4 12.6 80. { 1-phase distributed line + REC 68. 1.2E4 { Series L-C load on line } 1 + GEN 10. 40. { Series R-L at source has no effect } 1 +BLANK card ending branch cards + GEN SEND -1.0 1.0 { Permanently close switch at sending end } 1 +BLANK card ending switch cards +14GEN 1.0 50. { 60-Hz phasor solution } -1. +BLANK card ending source cards +C Total network loss P-loss by summing injections = 1.702540183859E-02 +C Output for steady-state phasor switch currents. +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GEN SEND 2.81684507E-02 -1.93583035E-02 3.41790218E-02 -34.4982 1.40842254E-02 9.67915173E-03 + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GEN 1.0 100. { 60-Hz phasor solution } -1. +14GEN 1.0 50. +BLANK card terminating source cards for 2nd (and final) phasor solution. +C Total network loss P-loss by summing injections = 1.174643954145E-03 +C Output for steady-state phasor switch currents. +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GEN SEND 8.10826370E-04 -4.85126124E-03 4.91855416E-03 -80.5114 4.05413185E-04 2.42563062E-03 +C Add these 2 I-real parts: i(0) = .0810826370E-02 + 2.81684507E-02 = 2.89792770 + GEN SEND REC { Names of nodes for voltage output +C First 3 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 3 output variables are branch currents (flowing from the upper node to the lower node); +C Step Time GEN SEND REC GEN REC GEN +C SEND TERRA TERRA +C *** Phasor I(0) = 2.8979277E-02 Switch "GEN " to "SEND " closed in the steady-state. +C 0 0.0 2.0 2.0 .3104733 .0289793 .0288367 .0074208 +C 1 .1E-3 1.997533 1.997533 .3168525 .0298758 .0298478 .0089264 +C 2 .2E-3 1.990141 1.990141 .321957 .03074 .0308254 .0104164 +C 3 .3E-3 1.977849 1.977849 .3249705 .0315697 .0317676 .0118871 +C 4 .4E-3 1.960698 1.960698 .3259235 .0323632 .0326723 .0133347 +BLANK card ending node voltage output requests +C 200 .02 2.0 2.0 .310624 .0289646 .0288207 .0074193 +C 220 .022 1.118034 1.118034 .1167098 .0390118 .0406289 .0307644 +C 240 .024 -.5 -.5 -.231522 .0292955 .0304699 .0301821 +C 260 .026 -1.11803 -1.11803 .0536218 .0061988 .0055751 .0120831 +C 280 .028 -.5 -.5 .7647646 -.015757 -.016915 -.002155 +C 300 .03 0.0 0.0 .9147476 -.027333 -.027043 -.004344 +C 320 .032 -.5 -.5 .1025405 -.029283 -.027728 -.006411 +C 340 .034 -1.11803 -1.11803 -.858517 -.02491 -.02427 -.018205 +C 360 .036 -.5 -.5 -.946641 -.013213 -.014641 -.029034 +C 380 .038 1.118034 1.118034 -.226641 .0070384 .0051148 -.020297 +C 400 .04 2.0 2.0 .310915 .0289599 .0288158 .007419 +C Variable maxima : 2.0 2.0 .9786535 .0390245 .0406716 .0337317 +C Times of maxima : 0.0 0.0 .0093 .002 .0021 .0029 +C Variable minima : -1.125 -1.125 -1.04338 -.029475 -.028141 -.029278 +C Times of minima : .0342 .0342 .0151 .0115 .0111 .0363 + CALCOMP PLOT + 144 4. 0.0 40. SEND REC { Show node voltages + 194 4. 0.0 40. GEN SEND REC { Show branch currents +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 8th of 19 subcases is like the preceding, only for 3-phase distributed line +C modeling rather than single phase. For simplicity, the line is taken to be +C continuously transposed. So, to make the solution interesting, the load +C is unbalanced slightly. The effect is substantial (must be near resonance). +C As phasor output is becoming substantial, illustrate lack of phasor output. +C About receiving end voltages (1st plot), phase b is dominated by fundamental +C whereas phase c is dominated by the harmonic. On the other hand, phase a +C is a nice mixture of the two. Done. WSM, 15 November 1998 19:52 +POWER FREQUENCY, 50., { Unnecessary. This merely avoids warning about 50 Hz. +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) + .000100 .040 50. 50. { Note L in ohms and C in micromhos at 50 Hz + 1 1 1 0 1 -1 + 2 -1 5 5 20 20 +-1SENDA RECA 0.3 0.4 12.6 80. { 3-phase distributed line +-2SENDB RECB .03 0.1 6.0 80. +-3SENDC RECC + RECA 68. 1.2E4 { Series L-C load on line } 1 + RECB 68. 1.3E4 { Series L-C: note imbalance } 1 + RECC 68. 1.4E4 { Series L-C: note imbalance } 1 +BLANK card ending branch cards + GENA SENDA -1.0 1.0 { Permanently closed switch at sending end + GENB SENDB -1.0 1.0 { Output 1 of 3 switch currents as demo } 1 + GENC SENDC -1.0 1.0 +BLANK card ending switch cards +14GENA 1.0 50. 0.0 -1. +14GENB 1.0 50. -120. -1. +14GENC 1.0 50. +120. -1. +BLANK card ending source cards +C Total network loss P-loss by summing injections = 2.454594157981E-01 + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GENA 1.0 100. 0.0 -1. +14GENB 1.0 100. -120. -1. +14GENC 1.0 100. +120. -1. +14GENA 1.0 50. 0.0 +14GENB 1.0 50. -120. +14GENC 1.0 50. +120. +BLANK card terminating source cards for 2nd (and final) phasor solution. +C Total network loss P-loss by summing injections = 2.573693449262E-04 + SENDA SENDB SENDC RECA RECB RECC { Names of nodes for voltage output +C First 6 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 3 output variables are branch currents (flowing from the upper node to the lower node); +C Step Time SENDA SENDB SENDC RECA RECB RECC RECA RECB RECC +C TERRA TERRA TERRA +C *** Phasor I(0) = 2.0781091E-04 Switch "GENA " to "SENDA " closed in the steady-state. +C *** Phasor I(0) = -6.9613077E-03 Switch "GENB " to "SENDB " closed in the steady-state. +C *** Phasor I(0) = 6.6629631E-03 Switch "GENC " to "SENDC " closed in the steady-state. +C 0 0.0 2.0 -1. -1. 2.598233 -3.43615 .2464431 .0129326 -.145918 .1293339 +C 1 .1E-3 1.997533 -.917186 -1.08035 2.60044 -3.42472 .2130051 .0099301 -.135491 .1225758 +C 2 .2E-3 1.990141 -.832151 -1.15799 2.597607 -3.40919 .1810295 .0069154 -.124907 .1156755 +C 3 .3E-3 1.977849 -.745147 -1.2327 2.589741 -3.38914 .1502293 .0038899 -.114174 .1086408 +BLANK card ending node voltage output requests +C 200 .02 2.0 -1. -1. 2.59818 -3.43448 .2454012 .0129912 -.146066 .1294166 +C 220 .022 1.118034 .7736591 -1.89169 1.768526 -2.57608 -.038877 -.047694 .0801733 -.023028 +C 240 .024 -.5 1.582676 -1.08268 .0182262 -1.32164 .5410584 -.098873 .2822797 -.164686 +C 260 .026 -1.11803 .8736191 .2444149 -1.03326 -.346505 .9871438 -.117453 .3728534 -.234779 +C 280 .028 -.5 -.064602 .5646021 -1.00872 .8642096 .2517987 -.08566 .3122374 -.211838 +C 300 .03 0.0 -.34E-13 .343E-13 -.883622 2.5928 -1.14969 -.01258 .1306647 -.114579 +C 320 .032 -.5 .5646021 -.064602 -1.27837 3.751705 -1.65974 .0650869 -.093075 .0189995 +C 340 .034 -1.11803 .2444149 .8736191 -1.43002 2.889824 -.686481 .1091923 -.274776 .1473061 +C 360 .036 -.5 -1.08268 1.582676 -.329683 .1403614 .6215637 .1064364 -.355258 .2280163 +C 380 .038 1.118034 -1.89169 .7736591 1.57799 -2.55929 .8876828 .0685474 -.308851 .2250084 +C 400 .04 2.0 -1. -1. 2.598318 -3.43353 .2448074 .0130473 -.146197 .1294821 + CALCOMP PLOT + 144 4. 0.0 40. RECA RECB RECC { Show node voltages + 194 4. 0.0 40. BRANCH + RECA RECB RECC { Show branch currents +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 9th of 19 subcases is like 7th, but with a frequency-dependent JMarti line +C replacing the constant-parameter, distributed line. BPA's Dr. Tsu-huei Liu +C and yours truly (WSM) force this through together on December 2, 1998. The +C switch has nothing to do with JMarti, but was corrected at this time (before +C the printed switch current before step 0 was wrong, not being superimposed). +POWER FREQUENCY, 50., { Unnecessary. This merely avoids warning about 50 Hz. +PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 8 output columns) + .000100 .040 50. 50. { Note L in ohms and C in micromhos at 50 Hz + 1 1 1 1 1 -1 + 2 -1 5 5 20 20 +-1SEND REC 1. 0.00 -2 1 + 25 4.6006674814379022100E+02 + -3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02 + 1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03 + 2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03 + 2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04 + 6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06 + 5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07 + 1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07 + 1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07 + 3.72990420164290220E+07 + 2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01 + 1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00 + 9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01 + 5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03 + 1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04 + 4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06 + 6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07 + 1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07 + 2.89354778735701405E+07 + 14 3.8038812947515656600E-04 + 1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00 + 9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01 + 2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04 + -3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10 + -7.93497908986478425E+10 3.97029750605793992E+10 + 3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02 + 1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03 + 1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04 + 5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05 + 1.82589462762646523E+05 1.82772052225409076E+05 + 0.58951119 + 0.00000000 + REC 68. 1.2E4 { Series L-C load on line } 1 +BLANK card ending branch cards + GEN SEND -1.0 1.0 { Permanent closed switch at sending end } 1 +BLANK card ending switch cards +14GEN 1.0 50. { 50-Hz phasor solution } -1. +C --------------+------------------------------ +C From bus name | Names of all adjacent busses. +C --------------+------------------------------ +C SEND |REC *GEN * +C REC |TERRA *SEND * +C GEN |SEND * +C TERRA |REC * +C --------------+------------------------------ +C 1st phasor flow: SEND 1.0 1.0 .00104886622776 .00412295654806 .52443311388E-3 .52443311388E-3 +C 0.0 0.0 -.0039873111659 -75.2621647 .00199365558294 0.0021265 +C Total network loss P-loss by summing injections = 5.244331138815E-04 +BLANK card ending source cards + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GEN 1.0 100. { 100-Hz phasor solution } -1. +14GEN 1.0 50. { 50-Hz but without phasor } +C 2nd phasor flow: SEND 1.0 1.0 .28023824311E-3 .00167781755608 .14011912156E-3 .14011912156E-3 +C 0.0 0.0 -.001654248554 -80.3850824 .82712427702E-3 0.0006777 +C Total network loss P-loss by summing injections = 1.401191215566E-04 +BLANK card terminating source cards for 2nd (and final) phasor solution. +C First 3 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 Step Time GEN SEND REC GEN REC +C SEND TERRA +C *** Phasor I(0) = 1.3291045E-03 Switch "GEN " to "SEND " closed in the steady-state. +C 0 0.0 2.0 2.0 .103985185 .001329104 .001331291 +C 1 .1E-3 1.99753329 1.99753329 .102188523 .001557337 .001566943 +C 2 .2E-3 1.99014143 1.99014143 .100246634 .001782631 .001799837 +C 3 .3E-3 1.97784922 1.97784922 .097684893 .002004723 .002029515 +C 4 .4E-3 1.96069786 1.96069786 .09455166 .002223115 .002255424 +C 5 .5E-3 1.93874486 1.93874486 .090914028 .002437299 .002477041 +C 10 .1E-2 1.76007351 1.76007351 .064262337 .003427592 .003502327 +C 15 .0015 1.47879178 1.47879178 .027694976 .00424585 .004348298 +C 20 .002 1.11803399 1.11803399 -.01526506 .004848781 .004971008 +C 40 .004 -.5 -.5 -.15428049 .004854478 .004951474 + GEN SEND REC { Names of nodes for voltage output +BLANK card ending node voltage output requests +C 360 .036 -.5 -.5 -.15242702 -.00466918 -.00477064 +C 380 .038 1.11803399 1.11803399 .017604346 -.00297868 -.00310013 +C 400 .04 2.0 2.0 .104051169 .001333593 .001335541 +C Variable max: 2.0 2.0 .231329745 .005331241 .005460826 +C Times of max: 0.0 0.0 .0297 .023 .023 +C Variable min: -1.1249966 -1.1249966 -.1823352 -.00466918 -.00477457 +C Times of min: .0342 .0342 .015 .036 .0361 + CALCOMP PLOT + 144 4. 0.0 40. SEND REC { Show node voltages + 194 4. 0.0 40. REC GEN SEND { Show branch currents +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 10th of 19 subcases is like 8th, but with a frequency-dependent JMarti line +C replacing the constant-parameter, distributed line. Without difficulty, +C this was demonstrated to be believable early in the morning of 4 Dec 1998. +C The JMarti branch cards are copied from DCNEW-27 (in turn, from DCNEW-4). +POWER FREQUENCY, 50., { Unnecessary. This merely avoids warning about 50 Hz. +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) + .000100 .040 50. 50. { Note L in ohms and C in micromhos at 50 Hz +C .000100 .0002 50. 50. { Note L in ohms and C in micromhos at 50 Hz + 1 1 1 1 1 -1 + 2 -1 5 5 20 20 +C PUNCHED CARD OUTPUT OF "JMARTI SETUP" WHICH BEGAN AT 14.32.51 09/14/82 +C 1.3636 .05215 4 1.602 -20.75 50. 50. +C 1.3636 .05215 4 1.602 -19.25 50. 50. +C 2.3636 .05215 4 1.602 - 0.75 77.5 77.5 +C 2.3636 .05215 4 1.602 0.75 77.5 77.5 +C 3.3636 .05215 4 1.602 19.25 50. 50. +C 3.3636 .05215 4 1.602 20.75 50. 50. +C 0.5 2.61 4 0.386 -12.9 98.5 98.5 +C 0.5 2.61 4 0.386 12.9 98.5 98.5 +C +C 100. 5000. 1 138. 1 3-2 +C 100. 60.00 1 138. 1 3 +C 100. .01 1 138. 1 9 10 3 +C +-1SENDA RECA 1. -2 3 1 + 18 0.47487144790587633736E+03 + -0.117967342612603676E+01 -0.316671345021416945E+01 0.255662317832426673E+02 + -0.601654631032706213E+02 -0.118825306534453652E+03 0.128169142098635430E+04 + 0.589610442538186533E+04 0.292891488173549769E+05 0.115665757459954655E+06 + 0.456394278488149546E+06 0.153681652821077613E+07 0.681848054466802243E+07 + 0.155477562610202006E+08 0.231789897204889809E+08 0.143827766961221036E+08 + 0.241524708800130989E+08 0.245049034784953832E+08 0.476808364145524129E+08 + 0.321645453133055891E+00 0.897616891499021768E+00 0.165632367032913183E+01 + 0.166669419493737986E+01 0.253073682990780913E+01 0.200014655383098283E+02 + 0.124257951002701473E+03 0.657345842950598183E+03 0.278623017775457413E+04 + 0.117145310573639467E+05 0.427410884649070877E+05 0.202553264108715757E+06 + 0.968789502295299084E+06 0.296908522006270138E+07 0.737968316276072920E+07 + 0.121009167337626640E+08 0.132891445947112725E+08 0.252097199908134253E+08 + 14 0.86849922355034181931E-03 + 0.580543547692016874E-01 0.269744993210885436E+00 0.358646519088149721E+00 + 0.775351125410302933E+00 0.512138885239557462E+01 0.838265420474258915E+01 + 0.700229782114550385E+02 0.393084113996385128E+03 0.138208543918962840E+04 + 0.706160326497018531E+03 0.527246586822120662E+04 0.528343251309079615E+05 + 0.178307801042203880E+08 -0.178914532139174500E+08 + 0.228360173022402533E+02 0.102750699821169636E+03 0.139222638455845377E+03 + 0.269608322018776320E+03 0.389220463867018722E+03 0.427793721751318500E+03 + 0.132148883486994549E+04 0.277178340638843952E+04 0.449341296469570784E+04 + 0.596603161966094024E+04 0.114714152641659587E+05 0.268314807408425590E+05 + 0.204027813729699328E+05 0.204231841543429041E+05 +-2SENDB RECB 1. -2 3 1 + 13 0.28580875773089406522E+03 + 0.381761604345214016E+04 -0.269529222495194318E+04 0.445018224900628539E+03 + 0.214658130219351907E+03 0.897325461203551775E+02 0.130480129629523390E+03 + 0.683037485206990933E+02 0.746626715384045969E+02 0.133938729387751110E+04 + 0.249954495902361043E+04 0.187764116219525090E+05 0.391031620802855759E+06 + 0.562668431190421223E+07 + 0.358668999787749571E+01 0.376679186833120605E+01 0.706439034439901437E+01 + 0.109721634064781350E+02 0.140506334140286961E+02 0.234292596911152251E+02 + 0.407686998675834982E+02 0.797598814775078289E+02 0.129106739969348078E+04 + 0.242138030939135859E+04 0.182362646041542116E+05 0.379814243183306855E+06 + 0.548808385606378759E+07 + 13 0.74123751236187229426E-03 + 0.192189595814863363E+02 0.210970481418782463E+01 0.383406546226866283E+02 + 0.307999912886148377E+02 0.634013410317139090E+03 -0.253243195247046863E+03 + 0.210540381301789816E+05 0.755533947306502050E+05 0.229576557459444489E+07 + 0.607777584135737231E+13 -0.156620677810598129E+13 0.153143851607186420E+13 + -0.604300997216751013E+13 + 0.325736126030671755E+04 0.366694913707514260E+03 0.675727697588602189E+04 + 0.516024976370228876E+04 0.249456113770459378E+05 0.115067772934717399E+06 + 0.105205466546493019E+06 0.240930799961836921E+06 0.633276538663731131E+06 + 0.975820481990592802E+06 0.976796302472583498E+06 0.974176721795814519E+06 + 0.975150898517610403E+06 +-3SENDC RECC 1. -2 3 1 + 13 0.27248404992352924125E+03 + 0.928429573805075023E+03 0.205756155237142632E+03 0.538106630554037608E+03 + 0.105199818376817531E+03 0.821382476778671116E+02 0.133356875032521778E+03 + 0.654774478488469374E+02 0.612546121562179922E+02 0.905544916173387175E+02 + 0.973244316091135204E+03 0.267481871266201051E+04 0.138566411938941123E+05 + 0.646075558955914967E+06 + 0.297277477955690506E+01 0.413177069940672514E+01 0.739410817322770197E+01 + 0.103072536754267949E+02 0.135874091964585515E+02 0.234792992956042887E+02 + 0.384075708508238121E+02 0.692539345057236435E+02 0.997101922255402791E+02 + 0.941460937860312015E+03 0.259295974388396348E+04 0.134538093446956691E+05 + 0.628668738677182439E+06 + 20 0.74871165651714274120E-03 + 0.249456855794677017E-01 0.516637614002139312E+01 0.435227595525510835E+01 + 0.696293545843379624E+01 0.113102519337405665E+02 0.107077756312168120E+02 + 0.136526319661978881E+02 0.103057494642840792E+03 0.647491597081196815E+03 + 0.221213037293261823E+04 0.164010494681972204E+05 0.517776934139508076E+05 + 0.393233870054795443E+06 -0.101893253559274247E+05 -0.173413454816459129E+06 + 0.524218185885897590E+06 0.711044562955687729E+12 -0.559640645251973480E+12 + 0.582892127161358521E+12 -0.734296849907947113E+12 + 0.887741958857318814E+01 0.183177145564733675E+04 0.151655674045953921E+04 + 0.245910936284398730E+04 0.383273362703089043E+04 0.399285585783856465E+04 + 0.496339281485561617E+04 0.253464455063045348E+05 0.272156609769823021E+05 + 0.475215596432249440E+05 0.149791320991727040E+06 0.333884409922466803E+06 + 0.541583000670874084E+06 0.670984155731562481E+06 0.569769858667429973E+06 + 0.132622344670840449E+07 0.440534234233555838E+07 0.440974768467789434E+07 + 0.436787151891064062E+07 0.437223939042955148E+07 + 0.57154434 0.70710678 -0.41761362 + 0.00000000 0.00000000 0.00000000 + 0.58879039 0.00000000 0.80696823 + 0.00000000 0.00000000 0.00000000 + 0.57154434 -0.70710678 -0.41761362 + 0.00000000 0.00000000 0.00000000 + RECA 68. 1.2E4 { Series L-C load on line } + RECB 68. 1.3E4 { Series L-C: note imbalance } + RECC 68. 1.4E4 { Series L-C: note imbalance } +BLANK card ending branch cards + GENA SENDA -1.0 1.0 { Permanently closed switch at sending end + GENB SENDB -1.0 1.0 { Output 1 of 3 switch currents as demo } 1 + GENC SENDC -1.0 1.0 +BLANK card ending switch cards +14GENA 1.0 50. 0.0 -1. +14GENB 1.0 50. -120. -1. +14GENC 1.0 50. +120. -1. +BLANK card ending source cards + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GENA 1.0 100. 0.0 -1. +14GENB 1.0 100. -120. -1. +14GENC 1.0 100. +120. -1. +14GENA 1.0 50. 0.0 +14GENB 1.0 50. -120. +14GENC 1.0 50. +120. +BLANK card terminating source cards for 2nd (and final) phasor solution. + SENDA SENDB SENDC RECA RECB RECC { Names of nodes for voltage output +BLANK card ending node voltage output requests + CALCOMP PLOT + 144 4. 0.0 40. RECA RECB RECC { Show node voltages + 194 4. 0.0 40. BRANCH + SENDA RECA SENDB RECB SENDC RECC { Show branch currents +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 11th of 19 subcases is like 9th, only with the single-phase JMarti line +C replaced by a Semlyen line. BPA's Dr. Tsu-huei Liu and WSM force this +C through together on December 8, 1998. The same line data was used to +C produce the Semlyen branch cards as was used to produce the Marti branch +C cards, according to Dr. Liu. +POWER FREQUENCY, 50., { Unnecessary. This merely avoids warning about 50 Hz. +PRINTED NUMBER WIDTH, 12, 2, { Request maximum precision (for 8 output columns) + .000100 .040 50. 50. { Note L in ohms and C in micromhos at 50 Hz + 1 1 1 1 1 -1 + 5 5 20 20 +C <++++++> Cards punched by support routine on 11-Nov-18 11.00.00 <++++++> +C SEMLYEN SETUP +C $ERASE +C TOLERANCES 10 5000. { Illustration only; value of FMED actually unchang +C 200 150 10 230 1 7 7777 0 { Semlyen miscellaneous da +C LINE CONSTANTS +C METRIC +C 1.3161 .1151 4 2.355 -6.5 22. 8.0 +C 0.5 4.189 4 0.914 4.5 27.3 15.7 +C BLANK card ending conductor cards within "LINE CONSTANTS" data +C 1500. 5000. 108.7 +C 1500. 60.00 108.7 { Phasor solution frequen +C 1500. 6.00 108.7 6 20 { log loopi +C BLANK card ending frequency cards of "LINE CONSTANTS" data +C BLANK card ending "LINE CONSTANTS" data +C L= 67.5 miles, rho= 1500.0, ss freq= 60.00, NSS=0, KFIT=10, KPS=2, KYC=30 +-1SEND REC 2.16583E-03 3.86702E-04 1 1 2 2 1 + 1.58360657E+01 9.71198911E+01-4.13400101E-05 2.87567721E-04 6.00000000E+01 + 0.00000E+00 4.76151E+04 8.46981E-01 0.00000E+00 3.44533E+03 1.53019E-01 + 0.00000E+00 5.15658E+04-3.35305E-04 0.00000E+00 1.93051E+02-6.36116E-04 + 1.00000E+00 0.00000E+00 + 1.00000E+00 0.00000E+00 + REC 68. 1.2E4 { Series L-C load on line } 1 +BLANK card ending branch cards + GEN SEND -1.0 1.0 { Permanent closed switch at sending end } 1 +BLANK card ending switch cards +14GEN 1.0 50. { 50-Hz phasor solution } -1. +BLANK card ending source cards +C Total network loss P-loss by summing injections = 1.734151791008E-03 +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GEN SEND 3.46830358E-03 -1.43381118E-02 1.47516298E-02 -76.4017 1.73415179E-03 7.16905591E-03 + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +C 14GEN 1.0 100. { 100-Hz phasor solution } -1. +14GEN 1.0 100. { 50-Hz but without phasor } -1. +14GEN 1.0 50. +BLANK card terminating source cards for 2nd (and final) phasor solution. +C Total network loss P-loss by summing injections = 1.065011729272E-04 +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GEN SEND 2.13002346E-04 -3.64646410E-03 3.65267990E-03 -86.6570 1.06501173E-04 1.82323205E-03 + GEN SEND REC { Names of nodes for voltage output +C First 3 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 Step Time GEN SEND REC GEN REC +C SEND TERRA +C *** Phasor I(0) = 3.6813059E-03 Switch "GEN " to "SEND " closed in the steady-state. +C 0 0.0 2.0 2.0 .153929313 .003681306 .003725237 +C 1 .1E-3 1.99753329 1.99753329 .1522527 .00435802 .004425728 +C 2 .2E-3 1.99014143 1.99014143 .149645253 .005027418 .005119878 +C 3 .3E-3 1.97784922 1.97784922 .145951055 .00568938 .005806382 +C 4 .4E-3 1.96069786 1.96069786 .141287156 .00634313 .006483935 +C 5 .5E-3 1.93874486 1.93874486 .136275468 .006989792 .007151414 +BLANK card ending node voltage output requests +C 10 .1E-2 1.76007351 1.76007351 .092891975 .010026854 .010291731 +C 15 .0015 1.47879178 1.47879178 .023317891 .012640295 .012991794 +C 20 .002 1.11803399 1.11803399 -.06174101 .014720636 .015118963 +C 40 .004 -.5 -.5 -.34443045 .016575341 .016838847 +C 60 .006 -1.118034 -1.118034 -.17058748 .010182945 .010031014 +C 80 .008 -.5 -.5 .362487911 .002326162 .002053321 +C 100 .01 0.0 0.0 .612908752 -.00296002 -.00293285 +C 120 .012 -.5 -.5 .248154876 -.00738972 -.00710881 +C 140 .014 -1.118034 -1.118034 -.30076592 -.01262818 -.01254386 +C 160 .016 -.5 -.5 -.40730556 -.01503915 -.01535767 +C 180 .018 1.11803399 1.11803399 -.06530854 -.00934344 -.00970902 +C 200 .02 2.0 2.0 .154986325 .00333684 .003383348 +C 220 .022 1.11803399 1.11803399 -.06853358 .01445717 .014861049 +C 240 .024 -.5 -.5 -.35906038 .016441167 .016706006 +C 260 .026 -1.118034 -1.118034 -.18786425 .010191471 .010041491 +C 280 .028 -.5 -.5 .346320129 .002467142 .002193761 +C 300 .03 0.0 0.0 .601636815 -.002726 -.00269921 +C 320 .032 -.5 -.5 .244116524 -.00711585 -.00683621 +C 340 .034 -1.118034 -1.118034 -.29695787 -.01237389 -.01229047 +C 360 .036 -.5 -.5 -.39659111 -.01485549 -.01517475 +C 380 .038 1.11803399 1.11803399 -.05015223 -.00926464 -.00963034 +C 400 .04 2.0 2.0 .171432693 .003300084 .003347 +C Variable maxima : 2.0 2.0 .616129386 .017112417 .017478995 +C Times of maxima : 0.0 0.0 .0098 .0033 .0033 +C Variable minima : -1.1249966 -1.1249966 -.43847067 -.01510408 -.01538971 +C Times of minima : .0342 .0342 .0153 .0157 .0158 + CALCOMP PLOT + 144 4. 0.0 40. SEND REC { Show node voltages + 194 4. 0.0 40. REC GEN SEND { Show branch currents +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 12th of 19 subcases is like 10th, only with the three-phase JMarti line +C replaced by a 3-phase Semlyen line. The Semlyen branch cards are copied +C from DC-31. Note XOPT and COPT have been set to 60 Hz to match the +C branch cards at the 1st of 2 frequencies. Unlike the high-order JMarti +C model, 2nd-order Semlyen involves substantial error, & is more biased +C toward a match at power frequency. So, the switch to 60 favors smoother +C Semlyen operation at the dominant frequency. Of course, the second freq +C of 120 Hz is not well fitted, but it does not dominate, so error is of +C less importance. The resulting plots are reasonably smooth. If really +C smooth plots are wanted, change the 2nd frequency to the 1st, so that +C the superposition has the effect of doubling power-frequency excitation. +C This is an independent, separate way to verify the superposition while +C avoiding error that is associated with the constant-inductance phasor +C representation of Semlyen modeling. +POWER FREQUENCY, 60., { Ensure that 60. two cards below does not produce warning +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) +C Steps/cycle Cycles --- Alternate to DELTAT and TMAX ensures full accuracy + 200 2 60. 60. { Note L in ohms and C in micromhos at 60 Hz + 1 1 1 1 1 -1 + 2 -1 5 5 20 20 +C L= 138.0MILES, RHO= 27.0, SS FREQ= 60.00, NSS=0, KFIT= 0, KPS=2, KYC=10 +C 1.3636 .05215 4 1.602 -20.75 50. 50. +C 1.3636 .05215 4 1.602 -19.25 50. 50. +C 2.3636 .05215 4 1.602 - 0.75 77.5 77.5 +C 2.3636 .05215 4 1.602 0.75 77.5 77.5 +C 3.3636 .05215 4 1.602 19.25 50. 50. +C 3.3636 .05215 4 1.602 20.75 50. 50. +C 0.5 2.61 4 0.386 -12.9 98.5 98.5 +C 0.5 2.61 4 0.386 12.9 98.5 98.5 +C +C 27. 5000. 1 138. +C 27. 60.00 1 138. +C 27. 6.00 1 138. 6 20 +C +-1SENDA RECA 0.55456E-02 0.77998E-03 1 1 2 2 3 1 + 0.14938602E+02 0.90516313E+02-0.16261865E-03 0.14500324E-02 0.60000000E+02 + 0.00000E+00 0.15161E+05 0.75119E+00 0.00000E+00 0.17105E+04 0.24881E+00 + 0.00000E+00 0.59584E+03-0.11954E-02 0.00000E+00 0.39933E+05-0.74162E-03 +-1SENDB RECB 0.74392E-02 0.74149E-03 2 2 2 2 3 1 + 0.83801231E+00 0.38634735E+02-0.44536501E-04 0.20637999E-02 0.60000000E+02 + 0.00000E+00 0.26608E+06 0.83767E+00 0.00000E+00 0.35627E+05 0.16233E+00 + 0.00000E+00 0.13630E+03-0.49171E-03 0.00000E+00 0.32940E+05-0.75349E-04 +-1SENDC RECC 0.42823E-02 0.74017E-03 3 3 2 2 3 1 + 0.26535168E+01 0.67894100E+02-0.45225808E-04 0.11682248E-02 0.60000000E+02 + 0.00000E+00 0.61698E+06 0.90903E+00 0.00000E+00 0.15239E+05 0.90969E-01 + 0.00000E+00 0.25667E+03-0.34465E-03 0.00000E+00 0.11831E+05-0.28233E-04 + 0.10000E+01 0.00000E+00 0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00 + 0.71114E+00 0.00000E+00-0.41983E-15 0.00000E+00 0.10000E+01 0.00000E+00 + 0.10000E+01 0.00000E+00-0.10000E+01 0.00000E+00-0.26698E+00 0.00000E+00 + 0.42027E+00 0.00000E+00 0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00 + 0.22471E+00 0.00000E+00-0.27373E-15 0.00000E+00 0.84054E+00 0.00000E+00 + 0.42027E+00 0.00000E+00-0.50000E+00 0.00000E+00-0.29908E+00 0.00000E+00 + RECA 68. 1.2E4 { Series L-C load on line } + RECB 68. 1.3E4 { Series L-C: note imbalance } + RECC 68. 1.4E4 { Series L-C: note imbalance } +BLANK card ending branch cards + GENA SENDA -1.0 1.0 { Permanently closed switch at sending end + GENB SENDB -1.0 1.0 { Output 1 of 3 switch currents as demo } 1 + GENC SENDC -1.0 1.0 +BLANK card ending switch cards +14GENA 1.0 60. 0.0 -1. +14GENB 1.0 60. -120. -1. +14GENC 1.0 60. +120. -1. +C Total network loss P-loss by summing injections = 7.920824547881E-04 +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GENA SENDA 2.13684848E-04 -1.46728245E-02 1.46743804E-02 -89.1656 1.06842424E-04 7.33641223E-03 +C GENB SENDB -1.32003123E-02 7.47886007E-03 1.51717367E-02 150.4655 6.16366730E-05 7.58561792E-03 +C GENC SENDC 1.11079349E-02 7.85331948E-03 1.36037070E-02 35.2603 6.23603358E-04 6.77320678E-03 +BLANK card ending source cards + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GENA 1.0 120. 0.0 -1. +14GENB 1.0 120. -120. -1. +14GENC 1.0 120. +120. -1. +14GENA 1.0 60. 0.0 +14GENB 1.0 60. -120. +14GENC 1.0 60. +120. +C Total network loss P-loss by summing injections = 2.436936258789E-05 +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GENA SENDA -9.06090319E-05 -2.82087453E-03 2.82232937E-03 -91.8398 -4.53045160E-05 1.41043726E-03 +C GENB SENDB -2.57770624E-03 1.49235343E-03 2.97853794E-03 149.9315 -1.78143046E-06 1.48926790E-03 +C GENC SENDC 2.32482625E-03 1.50725803E-03 2.77067570E-03 32.9566 7.14553090E-05 1.38349380E-03 +BLANK card terminating source cards for 2nd (and final) phasor solution. + RECA RECB RECC { Names of nodes for voltage output +C First 3 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 4 output variables are branch currents (flowing from the upper node to the lower node); +C Step Time RECA RECB RECC GENB SENDA SENDB SENDC +C SENDB RECA RECB RECC +C *** Phasor I(0) = 1.2307582E-04 Switch "GENA " to "SENDA " closed in the steady-state. +C *** Phasor I(0) = -1.5778019E-02 Switch "GENB " to "SENDB " closed in the steady-state. +C *** Phasor I(0) = 1.3432761E-02 Switch "GENC " to "SENDC " closed in the steady-state. +C 0 0.0 .1754598 -.149967 -.202628 -.015778 .128E-3 -.015779 .0134377 +C 1 .83333E-4 .19681 -.130782 -.20007 -.016119 .7992E-3 -.016119 .0131218 +C 2 .16667E-3 .1936652 -.115517 -.222603 -.016425 .0014416 -.016425 .0127552 +C 3 .25E-3 .1874392 -.096275 -.245921 -.016699 .0020782 -.016699 .0123656 +C 4 .33333E-3 .1801201 -.074648 -.267768 -.016946 .0027102 -.016946 .0119566 +C 5 .41667E-3 .1719231 -.051722 -.287869 -.017166 .0033372 -.017166 .0115298 +C 10 .83333E-3 .117421 .0693951 -.362436 -.017776 .0063154 -.017776 .0090996 +C 15 .00125 .0387023 .185327 -.393131 -.017831 .0090412 -.017831 .0064327 +C 20 .00166667 -.070868 .3086067 -.396284 -.017281 .0114347 -.017281 .0036444 +C 40 .00333333 -.388621 .3051482 .0228398 -.009941 .0156904 -.009941 -.006866 +C 60 .005 -.253038 -.204985 .4417365 -.778E-4 .0122998 -.778E-4 -.01189 +C 80 .00666667 .3044138 -.552418 .274188 .0068829 .0057587 .0068829 -.011409 +C 100 .00833333 .6360443 -.274231 -.265214 .0106411 -.3E-3 .0106411 -.008773 +BLANK card ending node voltage output requests +C 120 .01 .3181122 .3033606 -.463545 .0128573 -.006125 .0128573 -.005087 +C 140 .01166667 -.263527 .4887753 -.075016 .0123479 -.012279 .0123479 .0012662 +C 160 .01333333 -.410804 .1013287 .3706173 .0059884 -.015486 .0059884 .0098943 +C 180 .015 -.059755 -.30237 .2804666 -.005615 -.011145 -.005615 .0157608 +C 200 .01666667 .1817737 -.168327 -.19918 -.015787 .1853E-3 -.015787 .0134703 +C 220 .01833333 -.060762 .2780688 -.379241 -.017345 .0114781 -.017345 .0036718 +C 240 .02 -.398657 .32646 .0125865 -.009966 .0157061 -.009966 -.006836 +C 260 .02166667 -.250486 -.182702 .436643 -.129E-5 .0122619 -.129E-5 -.011919 +C 280 .02333333 .3159701 -.576173 .2854261 .0069158 .0057519 .0069158 -.011419 +C 300 .025 .6347384 -.288042 -.260816 .0105952 -.293E-3 .0105952 -.008748 +C 320 .02666667 .31108 .3297626 -.471975 .0128834 -.006175 .0128834 -.005106 +C 340 .02833333 -.257672 .4991878 -.073957 .0124327 -.012327 .0124327 .0012257 +C 360 .03 -.405636 .0768886 .3794062 .0059618 -.015473 .0059618 .0098998 +C 380 .03166667 -.0677 -.305613 .2785391 -.005706 -.011136 -.005706 .0157771 +C 400 .03333333 .1769868 -.142296 -.206482 -.015773 .1537E-3 -.015773 .0134395 +C variable max: .6370846 .5329181 .467254 .0132316 .0157061 .0132316 .016083 +C Times of max: .00825 .0111667 .0054167 .0105833 .02 .0105833 .0320833 +C Variable min: -.431775 -.577185 -.481248 -.017965 -.015486 -.017965 -.012246 +C Times of min: .0295 .0234167 .0264167 .01775 .0133333 .01775 .0056667 + CALCOMP PLOT + 142 .2 0.0 2.0 RECA RECB RECC { Show node voltages + 192 .2 0.0 2.0 BRANCH + SENDA RECA SENDB RECB SENDC RECC { Show branch currents +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 13th of 19 subcases is a generalization of the 5th. Superposition of phasor +C solutions is illustrated for simple series R-L-C elements. But here, we +C add a 3rd harmonic. Use 1, 2, and 3 Hz. The resulting solution has period +C 1 Hz, of course, so simulating and plotting to 2 sec shows 2 cycles. The +C periodicity is seen by comparing step 0 with 100 with 200 (very close). +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .01 2.0 + 1 1 1 1 1 -1 + 5 5 20 20 + GEN TRAN 1.0 10. 3 + TRAN 0.20E6 1 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 1.0 1.0 -1. +C TRAN .37942714659548 .64183619997635 .65053122009197 .80655515626147 -.752165257E-17 -.752165257E-17 +C -.5176762965662 -53.7607135 .47680221452676 36.2392865 -.2588381482831 -0.2588381 +C Total network loss P-loss by summing injections = 3.252656100460E-01 +BLANK card terminating source cards for 1st phasor solution. + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GEN -1.0 2.0 { 2nd harmonic source } -1. +C TRAN -.1008414235014 .38391636622821 -.931007055129 .96488706858836 -.545150405E-17 -.545150405E-17 +C .37043593719304 105.2282576 -.2534421401997 -164.7717424 -.1852179685965 -0.1852180 +C Total network loss P-loss by summing injections = 4.655035275645E-01 +BLANK card terminating source cards for 2nd (and final) phasor solution. + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GEN -1.0 3.0 { 2nd harmonic source } -1. +14GEN 1.0 1.0 { Source of 1st solution represented in dT-loop +14GEN -1.0 2.0 { Source of 2nd solution represented in dT-loop +C TRAN -.020242653071 .26448015717188 -.9941420093334 .99706670254975 -.575982404E-19 -.575982404E-19 +C .26370435820121 94.3895692 -.0763130042124 -175.6104308 -.1318521791006 -0.1318522 +C Total network loss P-loss by summing injections = 4.970710046667E-01 +BLANK card terminating source cards for 2nd (and final) phasor solution. + GEN TRAN +C First 3 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 Step Time GEN GEN TRAN GEN TRAN +C TRAN TRAN TERRA +C 0 0.0 -1.25834307 -1. .25834307 -1.27461784 -1.27461784 +C 1 .01 -1.17197835 -.976375224 .19560313 -1.23497976 -1.23497976 +C 2 .02 -1.04223034 -.906244946 .135985399 -1.14972948 -1.14972948 +C 3 .03 -.87351371 -.791817161 .081696549 -1.02182451 -1.02182451 +C 4 .04 -.671453154 -.636692146 .0347610081 -.855597126 -.855597126 +C 5 .05 -.442702206 -.44574573 -.003043524 -.656584162 -.656584162 +C 10 0.1 .8252797519 .8090169944 -.016262758 .5778901701 .5778901701 +C 15 .15 1.589742247 1.847858763 .2581165163 1.494230406 1.494230406 +C 20 0.2 1.283365647 1.927050983 .6436853361 1.410360926 1.410360926 +C 40 0.4 -1.48909279 -1.42705098 .0620418032 -1.57083779 -1.57083779 +C 60 0.6 -.676894279 -1.42705098 -.750156704 -.672056125 -.672056125 +C 80 0.8 2.141878198 1.927050983 -.214827215 2.10784098 2.10784098 +C 100 1.0 -1.2589336 -1. .2589336049 -1.27496681 -1.27496681 +BLANK card ending program output-variable requests. +C 120 1.2 1.28315489 1.927050983 .6438960933 1.410138424 1.410138424 +C 140 1.4 -1.4891661 -1.42705098 .0621151159 -1.57091519 -1.57091519 +C 160 1.6 -.676919781 -1.42705098 -.750131202 -.672083048 -.672083048 +C 180 1.8 2.141869327 1.927050983 -.214818344 2.107831615 2.107831615 +C 200 2.0 -1.25893669 -1. .2589366907 -1.27497007 -1.27497007 +C Variable max: 2.141878198 2.031786442 .8934909897 2.120254744 2.120254744 +C Times of max: 0.8 1.82 1.26 .81 .81 +C Variable min: -1.7217357 -1.43964823 -.935581739 -1.70523806 -1.70523806 +C Times of min: 1.36 1.61 .68 1.37 1.37 + CALCOMP PLOT + 143 .2 0.0 2.0 GEN TRAN { Screen plot is perfectly smooth & periodic +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C 14th of 19 subcases is a combination of the preceding. Use three single- +C phase lines: 1) Constant-parameter from 7th subcase; 2) JMARTI from 9th +C subcase; and 3) Semlyen from the 11th. Connect all 3 in series to verify +C that they all can coexist. Also, use 3 harmonics: 50, 75, and 100 Hz. The +C result is seen to be perfectly smooth. Terminal nodes SEND and REC are +C driven by voltage sources, so voltage here indicates nothing. But interior +C nodes SEM and REC provide a good measure of smoothness. The period is +C equal to 40 msec, note, because 75 Hz has a period of 13.3333 msec. Two +C cycles at 50 Hz, 3 cycles at 75 Hz, and 4 cycles at 100 Hz --- all are +C equal to 40 msec. So, step 400 agrees closely with step 0. Consider the +C two node voltages and 3 line currents: +C Step Time MARTI SEM GEN MARTI CAP +C SEND CAP SEM +C 0 0.0 2.746223 .7338817 -.004095 -.004261 -.006559 +C 400 .04 2.745004 .7255618 -.004087 -.004257 -.006555 +POWER FREQUENCY, 50., { Unnecessary. This merely avoids warning about 50 Hz. +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) + .000100 .040 50. 50. { Note L in ohms and C in micromhos at 50 Hz + 1 1 1 1 1 -1 + 2 -1 5 5 20 20 +-1SEND MARTI 0.3 0.4 12.6 80. { 1-phase distributed line + MARTI 2.0 7.E2 { Shunt compensation at end of 1st line +-1MARTI CAP 1. 0.00 -2 1 1 + 25 4.6006674814379022100E+02 + -3.53674337977459482E+03 4.74463047336912496E+03 4.19877687667209159E+02 + 1.20883875640704401E+03 1.37235832926770309E+03 1.76629378124055439E+03 + 2.57098916005431466E+03 3.78828017768501514E+03 1.42803850501766374E+03 + 2.65038870126655566E+04 5.80549633718908008E+04 2.89862514502709019E+04 + 6.76266799893490243E+04 1.75748067781355756E+05 1.15921506756057055E+06 + 5.68321189666407090E+06 9.19433740923774988E+06 1.70951407206966021E+07 + 1.81624149219296910E+07 1.63578536348816883E+07 2.07902744942919687E+07 + 1.63600113702911996E+07 3.05274574257793949E+07 2.82210286690547951E+07 + 3.72990420164290220E+07 + 2.77486902230361210E-01 2.68023805044265262E-01 4.15219939750085032E-01 + 1.00965215217368676E+00 2.08586147641298370E+00 4.31742025828644139E+00 + 9.38105998645288964E+00 2.16289190303449140E+01 4.71679444226597938E+01 + 5.00972850164923444E+02 2.19667801322956757E+03 2.86760537318643128E+03 + 1.14040516604910099E+04 2.97247489619667322E+04 9.91819942418309948E+04 + 4.96352587435292895E+05 1.62826284096999770E+06 3.10812825436382648E+06 + 6.60681297691485007E+06 1.23076395640831069E+07 1.54718126469801944E+07 + 1.14181900748053603E+07 2.33626050734163225E+07 2.05401226659539938E+07 + 2.89354778735701405E+07 + 14 3.8038812947515656600E-04 + 1.91758353605922701E-02 2.66649324474153904E-01 2.54923849247581425E+00 + 9.46636313120961682E+00 2.85381149799323453E+01 7.38766350382716724E+01 + 2.49452844472431571E+03 7.01441951080157833E+03 6.88095829884308187E+04 + -3.34847951818855899E+06 3.31843503484260664E+06 3.96467674493048249E+10 + -7.93497908986478425E+10 3.97029750605793992E+10 + 3.35859634392794338E+00 4.68146974610416323E+01 4.33602484888353672E+02 + 1.24257407602821491E+03 1.74575005233975117E+03 1.47765127586259360E+03 + 1.69531663286821313E+04 2.85574137491229049E+04 7.36636621071784758E+04 + 5.22749591784824327E+05 5.23272341376608878E+05 1.82407055706939399E+05 + 1.82589462762646523E+05 1.82772052225409076E+05 + 0.58951119 + 0.00000000 + CAP 2.0 7.E2 { Shunt compensation at end of 2nd line + CAP SEM 3.5 2.4E4 { Series capacitor after Marti } 1 +-1SEM REC 2.16583E-03 3.86702E-04 1 1 2 2 1 + 1.58360657E+01 9.71198911E+01-4.13400101E-05 2.87567721E-04 6.00000000E+01 + 0.00000E+00 4.76151E+04 8.46981E-01 0.00000E+00 3.44533E+03 1.53019E-01 + 0.00000E+00 5.15658E+04-3.35305E-04 0.00000E+00 1.93051E+02-6.36116E-04 + 1.00000E+00 0.00000E+00 + 1.00000E+00 0.00000E+00 +BLANK card ending branch cards + GEN SEND -1.0 1.0 { Permanently close switch at sending end } 1 +BLANK card ending switch cards +14GEN 1.0 50. 0.0 { 50-Hz phasor } -1. +14REC 1.0 50. 90. -1. +C Total network loss P-loss by summing injections = 9.724220340709E-04 +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GEN SEND -1.63570226E-03 -3.98415669E-03 4.30685807E-03 -112.3207 -8.17851130E-04 1.99207835E-03 +BLANK card ending source cards for 1st phasor solution + SUPERIMPOSE MORE SOURCES: { Following are sources for 2nd phasor solution +14GEN 1.0 75. 0.0 { 75-Hz phasor } -1. +14REC 1.0 75. 90. -1. +C Total network loss P-loss by summing injections = 4.417836206255E-04 +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GEN SEND -1.32752292E-03 -1.70816724E-03 2.16336599E-03 -127.8530 -6.63761459E-04 8.54083619E-04 +BLANK card ending source cards for 2nd phasor solution + SUPERIMPOSE MORE SOURCES: { Following are sources for 3rd phasor solution +14GEN 1.0 100. 0.0 { 100-Hz phasor } -1. +14REC 1.0 100. 90. -1. +$DISABLE { Begin block of data cards to be ignored +C Note about preceding, third, excitation set. Because there is no source +C having T-start zero or positive, ATP logic that began 15 December 1998 +C will internally copy all sources of all preceding excitation sets, changing +C T-start to zero in order that these sources become effective on the first +C time step. In effect, the following sources are added internally: +14GEN 1.0 50. 0.0 { 50-Hz phasor } 0.0 +14REC 1.0 50. 90. 0.0 +14GEN 1.0 75. 0.0 { 75-Hz phasor } 0.0 +14REC 1.0 75. 90. 0.0 +C But note carefully that this helpful feature relies upon lack of any source +C other than phasor excitation. If any source of the final data set has zero +C or positive T-start, there will be no such internal addition. If the user +C specifies one, he must specify all sources having T-start non-negative. +$ENABLE { End block of data cards to be ignored +C Total network loss P-loss by summing injections = 2.720584624880E-04 +C Node-K Node-M I-real I-imag I-magn Degrees Power Reactive +C GEN SEND -1.13140476E-03 -3.94389795E-04 1.19817363E-03 -160.7823 -5.65702382E-04 1.97194897E-04 +BLANK card terminating source cards for 3rd (and final) phasor solution. + SEND MARTI SEM REC { Names of nodes for voltage output +C First 4 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 3 output variables are branch currents (flowing from the upper node to the lower node); +C Step Time SEND MARTI SEM REC GEN MARTI CAP +C SEND CAP SEM +C *** Phasor I(0) = -4.0946299E-03 Switch "GEN " to "SEND " closed in the steady-state. +C 0 0.0 3.0 2.746223 .7338817 .184E-15 -.004095 -.004261 -.006559 +C 1 .1E-3 2.996423 2.7249 .6273332 -.141308 -.003855 -.003925 -.00623 +C 2 .2E-3 2.985703 2.697847 .5194808 -.282232 -.003612 -.003582 -.005889 +C 3 .3E-3 2.967873 2.664322 .410436 -.422391 -.00336 -.003232 -.005536 +C 4 .4E-3 2.942985 2.624337 .3004596 -.561404 -.003101 -.002874 -.005171 +C 5 .5E-3 2.911115 2.57779 .1894443 -.698897 -.002834 -.002509 -.004796 +BLANK card ending node voltage output requests +C 400 .04 3.0 2.745004 .7255618 -.17E-12 -.004087 -.004257 -.006555 +C Variable max: 3.0 2.766848 2.205032 2.750818 .0072415 .007581 .0081707 +C Times of max: 0.0 .0397 .0375 .0369 .0051 .0044 .0049 +C Variable min:-2.07381 -2.03357 -2.03363 -2.75082 -.006569 -.008275 -.009585 +C Times of min: .0339 .0059 .0037 .0031 .0379 .0378 .0382 + CALCOMP PLOT + 144 4. 0.0 40. -3.0 3.0SEND MARTI SEM REC { Show node voltages + 194 4. 0.0 40. BRANCH + GEN SEND MARTI CAP CAP SEM { Show branch currents +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 15th of 19 subcases has the same solution as the 5th, which has one series +C R-L-C circuit. Here, the L remains linear, but is represented in the +C form of a Type-93 nonlinear reactor. Since the solution is identical to +C that already shown, only 3 time steps will be taken. +C 18 Nov 98, expand by the addition of Type-98 pseudo-nonlinear element. +C Add a second, independent circuit in parallel with the first. Node +C name TRAN93 will be used for the first and TRAN98 for the second. +C The solutions are identical, of course. +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) + .01 .03 + 1 1 1 1 +C GEN TRAN 10. 3 +C The preceding 10-mHenry inductor can be replaced by the following equivalent +C nonlinear element. The phasor solution will be identical. Inductance 10 mH +C means that when the current is one amp, the flux is .01 volt-sec, so: +93GEN TRAN93 1.0 .01 3 + 0.0 0.0 { 1st point being origin is request to reflect + 1.0 .01 { i = 1 amp corresponds to flux = .01 volt-sec + 2.0 .02 { 2nd segment is, in fact, an extension of 1st + 9999 + TRAN93 1.0 0.20E6 1 +93GEN TRAN98 1.0 .01 3 + 1.0 .01 { i = 1 amp corresponds to flux = .01 volt-sec + 2.0 .02 { 2nd segment is, in fact, an extension of 1st + 9999 + TRAN98 1.0 0.20E6 1 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 1.0 1.0 -1. +C --------------+------------------------------ +C From bus name | Names of all adjacent busses. +C --------------+------------------------------ +C GEN |TRAN93*TRAN98* +C TRAN93 |TERRA *GEN * +C TRAN98 |TERRA *GEN * +C TERRA |TRAN93*TRAN98* +C --------------+------------------------------ +BLANK card terminating source cards for 1st phasor solution. +C TRAN93 1.0299583666875 1.0307690952352 .65053122009197 .80655515626147 .32526561004599 .32526561004599 +C -.0408740820394 -2.2726007 .47680221452676 36.2392865 -.2588381482831 -0.2588381 +C Total network loss P-loss by summing injections = 6.505312200920E-01 +C ---- Initial flux of coil "GEN " to "TRAN93" = 6.50531220E-03 +C ---- Initial flux of coil "GEN " to "TRAN98" = 6.50531220E-03 + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GEN -1.0 2.0 { 2nd harmonic source } -1. +14GEN 1.0 1.0 { Source of 1st solution represented in dT-loop +BLANK card terminating source cards for 2nd (and final) phasor solution. +C TRAN93 -1.03184847863 1.0384598361934 -.931007055129 .96488706858836 .46550352756452 .46550352756452 +C .11699379699335 173.5312736 -.2534421401997 -164.7717424 -.1852179685965 -0.1852180 +C Total network loss P-loss by summing injections = 9.310070551290E-01 +C ---- Initial flux of coil "GEN " to "TRAN93" = -9.31007055E-03 +C ---- Initial flux of coil "GEN " to "TRAN98" = -9.31007055E-03 + GEN TRAN93TRAN98 { Request for voltage output of these nodes +C First 5 output variables are electric-network voltage differences (upper voltage minus lower voltage); +C Next 4 output variables are branch currents (flowing from the upper node to the lower node); +C Step Time GEN GEN GEN TRAN93 TRAN98 GEN GEN TRAN93 TRAN98 +C TRAN93 TRAN98 TRAN93 TRAN98 TERRA TERRA +C 0 0.0 .0018901 .0018901 0.0 -.00189 -.00189 -.280476 -.280476 -.280476 -.280476 +C 1 .01 .0137899 .0137899 .005912 -.007878 -.007878 -.272636 -.272636 -.272636 -.272636 +C 2 .02 .0251067 .0251067 .0235315 -.001575 -.001575 -.253188 -.253188 -.253188 -.253188 +C 3 .03 .0356219 .0356219 .0525108 .0168889 .0168889 -.222823 -.222823 -.222823 -.222823 +BLANK card ending program output-variable requests. +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C 16th of 19 subcases illustrates KILL = 222 as mentioned in the January, +C 1999, newsletter. The R-L-C connection of the preceding subcase is +C modified to connect to a second generator rather than ground. The 2nd +C generator, GEN2, has a different frequency than the 1st. Prior to +C correction on November 17th, execution continued in spite of the mixing +C of frequencies through a nonlinear element. Now, the error is trapped. +C Note no use of superposition, which is an extraneous complication. A +C single frequency is all that is required to demonstrate the phenomenon. +PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns) + .01 .03 + 1 1 1 1 +93GEN TRAN 1.0 .01 3 + 0.0 0.0 { 1st point being origin is request to reflect + 1.0 .01 { i = 1 amp corresponds to flux = .01 volt-sec + 2.0 .02 { 2nd segment is, in fact, an extension of 1st + 9999 + TRAN GEN2 1.0 0.20E6 1 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 1.0 1.0 -1. +14GEN2 1.0 2.0 -1. +BLANK card terminating source cards for 1st phasor solution. +BLANK card terminating source cards for 2nd (and final) phasor solution. + GEN TRAN +BLANK card ending program output-variable requests. +BLANK card ending all plot cards +BEGIN NEW DATA CASE +C 17th of 19 subcases illustrates IF-THEN-ELSE-ENDIF within HFS data that +C is the same as 3rd subcase. There are 8 frequency steps, and for the +C first five, the solution is identical to that of 3rd subcase. But for +C steps 6 onward, the formula is changed (the ELSE alternative). This +C data is being added by WSM on 3 December 2000. +POWER FREQUENCY, 25., ! Needed so minimum frequency is recognized as fundamental +HARMONIC FREQUENCY SCAN { Non-negative DELFFS in 25-32 means F in Hz (not log F) +C Controls of following request card: MAXKNT IOPCVP { Ignore MAXKNT value +POCKET CALCULATOR VARIES PARAMETERS 0 1 { since also HFS loop +PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) + .0001 -1. 50. 50. + 1 1 1 0 1 2 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b +IF( KNT .LT. 6.0 ) THEN { If frequency (pass number) is low, use one formula: +C Note: the following formula is identical to that used in the 3rd subcase. +C The present comment card illustrates that these are tolerated within +C an IF-THEN-ELSE-ENDIF block. +MHENRY = 1.2 - ( KNT - 1.0 ) * .05 { L = 1.2, 1.15, 1.1, 1.05, and 1.0 +ELSE { Alternatively (if frequency is higher, use another formula: +C Note: the following formula is different from that used in the 3rd subcase. +C As a result, the solution to steps 6, 7, and 8 differ slightly. We +C use the ELSE alternative to illustrate saturation of the previously- +C linear inductance. Slope of L vs. f changes by a factor of 5: +MHENRY = 1.0 - ( KNT - 5.0 ) * .01 { L =.99, .98, .97 +ENDIF { Terminate 5-line block that provides two alternatives for MHENRY +C The following variable is not actually used, but it is added to illustrate +C that other variables can follow an IF-THEN-ELSE-ENDIF block. As for the +C formula, it can be found in the 1st subcase. In fact, this addition makes +C the solution comparable to that of the 1st subcase. They are identical +C through step 5. +FREQUENCY = 25. * KNT { This computes the frequency in Hz of the scan +BLANK card ends $PARAMETER definitions that are processed just b4 branch card + TRAN 1.0 1.E+5 3 + GEN TRAN 1.0MHENRY +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 70. 25. 0.0 { Note comment and no negative T-start +14GEN 70. 50. 0.0 { This is normal HFS source data. It +14GEN 70. 75. 0.0 { is unusual, however, because of our +14GEN 70. 100. 0.0 { match preceding subases. Unusual +14GEN 70. 125. 0.0 { for HFS usage, the amplitude is not +14GEN 70. 150. 0.0 { varied and all harmonics are present +14GEN 70. 175. 0.0 +14GEN 70. 200. 0.0 +BLANK card ending source cards +BLANK card ending F-dependent series R-L-C branches (none, for this subcase) + GEN TRAN +BLANK card ending program output-variable requests. + CALCOMP PLOT { Use BOTH rather than DISK to see a bar chart of harmonics + 19680. 0.0200. TRAN { Units/inch of 5-7 are ingored since bar chart +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 18th of 19 subcases illustrates ELSEIF use within an IF-THEN-ELSE-ENDIF +C block. It is related to the preceding, but for variety has replaced +C the HFS control of looping by simpler pocket calculator control as was +C seen in the 1st subcase. The first condition, giving 5 steps, has +C been left unchanged. But then the remaining 3 steps would need to be +C split among two formulas. Rather than 8 frequency steps, 9 steps will +C be taken in order that each formula is used at least twice. From top to +C bottom: a) below IF b) below ELSEIF c) below ELSE +C KNT = 1, 2, 3, 4, 5 steps 6 and 7 steps 8 and 9. +C This data is added by WSM on 5 December 2000. +C Controls of following request card: MAXKNT IOPCVP { Loop 9 times with +POCKET CALCULATOR VARIES PARAMETERS 9 1 { minimum printout +PRINTED NUMBER WIDTH, 11, 2, { Request maximum precision (for 8 output columns) + .0001 -1. 50. 50. + 1 1 1 0 1 2 +$PARAMETER { This will be serviced by CIMAGE just as any other $-card would b +IF( KNT .LT. 6.0 ) THEN { If frequency (pass number) is low, use one formula: +C Note: the following formula is identical to that used in the 3rd subcase. +C It is used for KNT < 6, so for values KNT = 1, 2, ... through 5: +MHENRY = 1.2 - ( KNT - 1.0 ) * .05 { L = 1.2, 1.15, 1.1, 1.05, and 1.0 +C Having evaluated this or any other formula of the block, control +C passes to the end of the block, of course. Control reaches the +C following statement if and only if the previous condition ( KNT < 6 ) +C was not satisfied. +ELSEIF( KNT .LT. 8.0 ) THEN { If frequency (pass number) satisfies KNT < 8 : +C Control reaches this point only for 5 < KNT < 8, so for KNT = 6 and 7. +C The slope has been reduced from .05 to .02 and the offset has been +C reduced to ensure continuity of the piecewise-linear segments: +MHENRY = 1.0 - ( KNT - 5.0 ) * .02 { Produces values L = .98 and .96 +C Having evaluated this formula, control passes to end of block (ENDIF) +C Etc. The preceding ELSEIF could be repeated an arbitrary number of +C times, providing different formulas for an arbitrary number of ranges +C of KNT. If the condition of the first ELSEIF was not satisfied, the +C second would be checked, etc. in order from top to bottom. If none is +C satisfied when ELSE finally was reached, the formula below ELSE +C would be used. +ELSE { Alternatively (if none of preceding conditions were satisfied), use this: +C Control reaches this point only if no preceding condition is satisfied. +C I.e., for 7 < KNT. The slope and offset will be changed once again, +C to provide a 3rd segment of the piecewise-linear curve: +MHENRY = .96 - ( KNT - 7.0 ) * .01 { L =.95, .94 +ENDIF { Terminate 26-line block that provides three alternatives for MHENRY +C Unlike preceding subcase, the following variable is required because HFS +C is not being used. Instead, the pocket calculator must itself provide +C the frequency for any source card (node GEN in this case). +FREQUENCY = 25. * KNT { This computes the frequency in Hz of the scan +BLANK card ends $PARAMETER definitions that are processed just b4 branch card + TRAN 1.0 1.E+5 3 + GEN TRAN 1.0MHENRY +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN 70. FREQUENCY -1. +BLANK card ending source cards + GEN TRAN +C Step F [Hz] TRAN GEN TRAN TRAN +C TERRA TERRA +C 1 25. 71.874164 70. 71.874164 3.5892245 +C New parameter values follow: 1) 1.15 2) 50. +C 2 50. 77.535296 70. 77.535296 7.7150503 +C New parameter values follow: 1) 1.1 2) 75. +C 3 75. 87.376164 70. 87.376164 12.96142 +C New parameter values follow: 1) 1.05 2) 100. +C 4 100. 101.32085 70. 101.32085 19.870653 +BLANK card ending program output-variable requests. +C New parameter values follow: 1) 1.0 2) 125. +C 5 125. 115.44696 70. 115.44696 28. +C New parameter values follow: 1) .98 2) 150. +C 6 150. 119.51423 70. 119.51423 34.342161 +C New parameter values follow: 1) .96 2) 175. +C 7 175. 102.75012 70. 102.75012 33.943547 +C New parameter values follow: 1) .95 2) 200. +C 8 200. 79.015248 70. 79.015248 29.345527 +C New parameter values follow: 1) .94 2) 225. +C 9 225. 60.191978 70. 60.191978 24.700662 + PRINTER PLOT + 19650. 0.0250. TRAN Plot limits: ( 0.000, 3.434 ) +BLANK card ending plot cards +BEGIN NEW DATA CASE +C 19th of 19 subcases is same as the 5th except that the voltage source used +C previously becomes a current source. Also, as first illustrated by the 14th +C subcase, the final excitation set will involve no preceding harmonics so it +C will be the job of ATP to add these automatically. Prior to correction on +C 13 March 2001, there was an enormous oscillation in dT loop as seen on plot. +C That is, the automatic addition of all preceding harmonics (here, just one) +C did not work properly for a current source. Orlando Hevia in Santa Fe, +C Argentina, first made WSM aware of the problem by E-mail. +PRINTED NUMBER WIDTH, 13, 2, { Request maximum precision (for 8 output columns) + .01 2.0 + 1 1 1 1 1 -1 + 5 5 20 20 + GEN TRAN 1.0 10. 3 + TRAN 0.20E6 1 +BLANK card ending program branch cards. +BLANK card terminating program switch cards (none, for this case) +14GEN -1 1.0 1.0 -1. +BLANK card terminating source cards for 1st phasor solution. + SUPERIMPOSE MORE SOURCES: { Request for superposition produced by more: +14GEN -1 -1.0 2.0 { 2nd harmonic source } -1. +BLANK card terminating source cards for 2nd (and final) phasor solution. + GEN TRAN +BLANK card ending program output-variable requests. + CALCOMP PLOT + 143 .2 0.0 2.0 GEN TRAN { Screen plot is perfectly smooth & periodic +BLANK card ending all plot cards +BEGIN NEW DATA CASE +BLANK -- cgit v1.2.3