BEGIN NEW DATA CASE C -------------------------------------------------------------------- C I *** IEEE POWER ENGINEERING SOCIETY HARMONICS WORKING GROUP *** I C I TASK FORCE ON HARMONICS MODELING AND SIMULATION I C I I C I HARMONICS TEST SYSTEMS PROJECT C I IEEE 13 BUS UNBALANCED DISTRIBUTION MODEL I C I Version 1.0 I C I *** SATISH J RANADE & JITENDRA S KANETKAR *** I C I I C ------------------------------------------------------------------------ C JUNE 1995 FIRST MODEL BY KANETKAR FOR PGE PROJECT C DECEMBER 1996 HARMONIC ANALYSIS USING 'MODELS' BY RANADE C --------------------------------------------------------------------- C --------------------------------------------------------------------- C C -------------------------------------------------------------------- C FOLLWING ARE USED FOR LOAD FLOW C C FIX SOURCE C -------------------------------------------------------------------- C -------------------------------------------------------------------- C -------------------------------------------------------------------- C SPECIAL REQUEST FOR FREQUENCY, COL. 33-40 C POWER FREQUENCY 60.0 C C -------------------------------------------------------------------- C -------------------------------------------------------------------- C TI,E STEP ETC .0001 0.046 060. 1 1 1 1 C 00000000000000000000000000000000000000000000000000000000000000000000000 C C THIS IMPLEMENTATION USES IMPLEMENTS A CURRENT SOURCE HARMONICS MODEL USING C 'MODELS' FOR HARMONIC PROPAGATION STUDIES. THE MODEL IS USED FOR DOCUMENTING C UNBALANCED HARMONIC TEST SYSTEMS IN THE PROOSED TASK FORCE PAPER C Acknowledgement: I want to recognize Mack Grady who originally suggested c and has developed 'MODELS' based harmonic simulation C THE NETWORK ( IEEE STANDARD 13 BUS DISTRIBUTION SYSTEM) IS MODELED IN THE C NETWORK PART. C C LOADS ARE SPECIFIED AT RELEVANT BUSES AND ARE ASSUMED TO CONTAIN BOTH LINEAR C AND HARMONIC PRODUCING COMPONENTS. SPECIFIC COMPNENTS ARE PASSIVE LINEAR, C MOTORS, AND HARMONIC PRODUCING FLUORESCENT LIGHTS, ADJUSTABLE SPEED DRIVES C AND OTHER(MIXED) COMPONENTS FOR WHICH A CURRENT SPECTRUM IS KNOWN C THESE COMPONENTS ARE SPECIFIED AS PERCENT OF TOTAL LOAD C C THE SOLUTION PROCESS IS AS FOLLOWS C C PERFORM EMTP FUNDAMENTAL FREQUENCY STEADY STATE SOLUTION C WITH ALL LOAD AS CONSTANT IMPEDANCE C C PERFORM SEVERAL CYCLES OF TIME SIMULATION WITH HARMONIC PRODUCING LOADS C MODELED AS CURRENT SOURCES C C CALCULATE SPECTRA C C 'MODELS' IS USED TO REPRESENT HARMONIC CURRENT SOURCES IN THE TIME SIMULATION. C THE SPECIFIC MODEL IS H_INJ C C INPUTS TO H_INJ ARE : C BUS VOLTAGES C IMAGINARY COMPONENT OF STEADY STATE (T=0) VOLTAGE C LOAD CURRENT C IMAGINARY COMPONENT OF STEADY STATE (T=0) LOAD CURRENT C PERCENTAGE OF EACH TYPE OF HARMONIC PRODUCING LOAD C AT EACH BUS C ///USER MUST SPECIFY BUS VOLTAGE NAMES IN THE MODELS INPUT SECTION// C INTENAL VARIABLES C NO OF BUSES AND HARMONIC SOURCE TYPES C MAGNITUDE ( % OF FUNDAMENTAL) AND PHASE WRT FUNADMENTAL C VOLTAGE FOR EACH SOURCE TYPE C ///USER MUST SET UP THIS DATA/// C PROCESS C INITIALIZATION C FROM THE EMTP STEADY STATE (T=0) SOLUTION CALCULATE C FUNDAMENTAL PHASOR VOLTAGES C C NET HARMONIC CURRENTS TO BE INJECTED C USING THE PERCENTAGE OF LOAD OF EACH HARMONIC LOADTYPE C AND THE STANDARD SPECTRUM CALCULATE THE TIME DOMAIN C CURRENT TO BE INJECTED. THE CURRENT INCLUDES THE C FUNDAMENTAL COMPONENT AND ALL COMPONENTS ARE ADJUSTED C FOR THE PHASE ANGLE OF ACTUAL BUS VOLTAGE AT FROM THE C STEADY STATE T<0 SOLUTIN C C COMPENSATING CURRENT C SINCE THE LOAD IMPEDANCE IN THE NETWORK MODEL REPRESENTS TOTAL C LOAD IT IS NECESSARY TO COMPENSATE THE PART OF THIS IMPEDANCE C THAT REPRESENTS THE HARMONIC PRODUCING LOAD. C THIS IS DONE BY CALCULATING A COMPENSATING CURRENT INJECTION. C FIRST WE ASSUME THAT CAPACITORS ARE REPRESENTED EXPLICITLY C THUS LOAD IS INDUCTIVE. WE CALCULATE THE PARALLEL R-L C REPRESENTATION OF THE PART OF THE BUS LOAD THAT CORRESPONDS C TO EACH HARMONIC PRODUCING LOAD TYPE. THE REQUIRED COMPENSATING C CURRENT IS THEN C V/R + (1/L) INTEGRAL(V) C THE INTEGRAL IS INITIALIZED FOR FUNDAMENTAL FREUENCY USING THE C PHASOR VALUE BASED ON STEADY STATE VOLTAGE PRIOR TO SIMULATION C C OUTPUT THE NET CURRENT SOURCE IS INJECTED INTO THE APPROPRIATE C BUS USING THE TYPE 60 SOURCE C C INTERFACE THE REQUIRED BUS VOLTAGES AND LOAD CURRENTS ARE MEASURED USING C C IDEAL 1:1 TRANSFORMER ACROSS LOAD C IDEAL 1:1 CT WITH 0.001 OHM BURDEN IN SERIES WITH LOAD C THE INJECTED TYPE 60 CURRENT SOURCE IS CONNECTED TO THE C IDEAL TRANSFORMER ACROSS THE LOAD C THIS ALLOWS PHASE-GROUND AND PHASE-PHASE CONNECTED LOAD C C MODELS C ------------------------------------------------------>> INPUT SECTION C INPUT VOLTAGE, SS IM(VOLTAGE) C CURRENT SS IM(CURRENT) C FOR EACH BUS WITH HARMONIC SOURCE C C RECALL FROM OPENING REMARKS THAT VOLTAGE IS MEASURED C THROUGH AN IDEAL TRANSFORMER ACROSS LOAD AND CURRENT C IS MEASURED VIA A SERIES IT WITH A 0.001 OHM BURDEN C MEASUREMENT BUSES WHICH ARE THE SECONDARIES OF THESE C TRANSFORMERS ARE DENOTED , EG FOR BUS #i(NOT ACTUSL EMTP NAME), C BY BUSNAME husi FOR VOLTAGE AND CUSi FOR CURRENT. C C BUSVRi{ v(HUSi)} VOLTAGE AT BUS i C BUSVIi{ imssv(CUSi)} IMAGINARY PART OF STEADY STATE VOLTAGE AT BUS i C C HUSIRi{ v(cUSi )} CURRENT AT BUS i C HUSIIi{ imssv(CUSi )} IMAGINARY PART OF STEADY STATE CURRENT AT BUS i C C INPUT BUSVR1{v(HUS1 )} BUSVI1{imssv(HUS1 )} LODCR1{v(CUS1 )} LODCI1{imssv(CUS1 )} BUSVR2{v(HUS2 )} BUSVI2{imssv(HUS2 )} LODCR2{v(CUS2 )} LODCI2{imssv(CUS2 )} BUSVR3{v(HUS3 )} BUSVI3{imssv(HUS3 )} LODCR3{v(CUS3 )} LODCI3{imssv(CUS3 )} BUSVR4{v(HUS4 )} BUSVI4{imssv(HUS4 )} LODCR4{v(CUS4 )} LODCI4{imssv(CUS4 )} BUSVR5{v(HUS5 )} BUSVI5{imssv(HUS5 )} LODCR5{v(CUS5 )} LODCI5{imssv(CUS5 )} BUSVR6{v(HUS6 )} BUSVI6{imssv(HUS6 )} LODCR6{v(CUS6 )} LODCI6{imssv(CUS6 )} BUSVR7{v(HUS7 )} BUSVI7{imssv(HUS7 )} LODCR7{v(CUS7 )} LODCI7{imssv(CUS7 )} BUSVR8{v(HUS8 )} BUSVI8{imssv(HUS8 )} LODCR8{v(CUS8 )} LODCI8{imssv(CUS8 )} BUSVR9{v(HUS9 )} BUSVI9{imssv(HUS9 )} LODCR9{v(CUS9 )} LODCI9{imssv(CUS9 )} BUSVR10{v(HUS10 )} BUSVI10{imssv(HUS10 )} LODCR10{v(CUS10 )} LODCI10{imssv(CUS10 )} BUSVR11{v(HUS11 )} BUSVI11{imssv(HUS11 )} LODCR11{v(CUS11 )} LODCI11{imssv(CUS11 )} C C ------------------------------------------------------>> OUTPUT SECTION C THE OUTPUTS ARE THE DISTORTED AND COMPENSATING CURRENT INJECTIONS INTO EACH C HARMONIC LOAD BUS. SUCH BUSES ARE DENOTED BY HUSI C AND ARE AT THE SECONDARY OF THE VOLTAGE MEASURING IDEAL C TRANSFORMER. A TYPE 60 SOURCE IS USED IN THE NETWORK C OUTPUT HUS1 ,HUS2,HUS3 ,HUS4,HUS5 ,HUS6 ,HUS7,HUS8,HUS9,HUS10, HUS11 C ,HUS12,HUS13,HUS14,HUS15,HUS16,HUS17,HUS18,HUS19,HUS20 C1,C2,P1,P2 C ------------------------------------------------------>> OUTPUT SECTION C C ------------------------------------------------------>> MODEL H_INJ MODEL H_INJ DATA tstart {dflt:0.0} C C NHBEG=1 IS FIRST HARMONIC CONSIDERIED; NHEND IS THE HIGHEST C NHS IS THE NUMBER OF HARMONIC SOURCE TYPES, NBUS IS THE # OF BUSES C CONST NHBEG{val : 1},nhend{val : 20},NHS{val:3} ,NBUS{val:11} C C THE INPUTS FROM THE MOIDWLS SECTION: C BUSVR IS THE INSTANTANEOUS VOLTAGE. AT T=0 IT IS THE REAL PART OF THE STEADY C STATE VOKTAGE WITHOUT HARMONICS C BUSVIR IS THE IMAGINARY PART OF STEADY STATE VOLTAGE (T<0) WITHOUT HARMONICS C LODCR IS THE INSTANTANEOUS CURRENT. AT T=0 IT IS THE REAL PART OF THE STEADY C STATE CURRENT WITHOUT HARMONICS C LODCI IS THE IMAGINARY PART OF STEADY STATE CURRENT (T<0) WITHOUT HARMONICS C PERi IS THE PERCENTAGE OF TOTAL LOAD CORRESPONDING TO HARMNIC SOURCE TYPE i INPUT BUSVR[1..NBUS],BUSVI[1..NBUS],LODCR[1..NBUS],LODCI[1..NBUS], PER1[1..NBUS],PER2[1..NBUS],PER3[1..NBUS] C C OMEGA IS SYSTEM FREQUENCY IN RAD/S C HSMi IS MAG IN PERCENT OF FUND OF HARMONIC SOURCE TYPE i C HSAi IS ANGLE WITH RESPECT TO FUNDAMENTAL VOLTAGE C G IS THE CONDUCTIONS (1/R) IN THE PARALLEL R-L REPRESENTATION OF C THE HARMONIC PRODUCING LOAD AT A BUS C LINV =1/L WHERE L IS THE INDUCTANCE IN THE PARALLEL R-L REPRESENTATION OF C THE HARMONIC PRODUCING LOAD AT A BUS C VMAG/PHASEV ID THE SS PHASOR VOLTAGE PRIOR TO TIME SIMULATION C CMAG IS THE SS PHASOR CURRENT MAGNITUDE PRIOR TO TIME SIMULATION C CINJT IS THE INJECTION THAT COMPENSATES FOR THE PARALLEL R-L PART OF LOAD IN C THE NETWORK MODEL THAT CORRESPONDS TO HARMONIC PRODUCING LOAD C HS IS THE TOTAL INJECTED CURRENT VAR OMEGA,I,J, HSM1[1..nhend], HSA1[1..nhend], HSM2[1..nhend], HSA2[1..nhend], HSM3[1..nhend], HSA3[1..nhend], Y[1..NBUS],G[1..NBUS], LINV[1..NBUS], VMAG[1..NBUS],CMAG[1..NBUS],PHASEV[1..NBUS] , CINJ[1..NBUS] ,CINJT[1..NBUS],HS[1..NBUS],HX , HISTORY BUSVR[1..NBUS]{dflt:0} C OUTPUT HS , PHASEV, CINJT INIT C spectrum of asd 1 C HSM1[1..nhend]:=[1 , 0 ,0.542, 0, 0.152, 0 , 0.069, 0,0.043, 0, 0.036, 0 ,0.029, 0, 0.025, 0.000,0.018,0.000,0.014,0 ] HSA1[1..nhend]:=[ -1.45, 0 ,0.66, 0,110.77 , 0,151.87, 0,-95.02, 0., -13.91, 0 ,95.2, 0,-182.7 , 0 , -91.59 , 0, 10.52, 0.] HSA1[1..nhend] := HSA1[1..nhend]*PI/180 C C spectrum of fluorescent light bank HSM2[1..nhend]:=[1 , 0 ,0.192, 0, 0.107, 0 , 0.021, 0,0.014, 0, 0.009, 0 ,0.006, 0, 0.005, 0.000,0.000,0.000,0.000,0 ] HSA2[1..nhend]:=[ -41.2, 0 ,273.40, 0,339.00 , 0,137.67, 0,263.20, 0., 39.80, 0 ,182.4 , 0,287.00 , 0 , 0 , 0, 0 , 0.] HSA2[1..nhend] := HSA2[1..nhend]*PI/180 C spectrum of general single phase load HSM3[1..nhend]:=[1 ,0,0.0074,0.0954 , 0.0017, 0.0832, 0.00, 0.0049,0, 0, 0.0, 0 ,0.0, 0, 0.0, 0.000,0.0,0.000,0.0,0 ] HSA3[1..nhend]:=[ -35.0, 0 ,-105.8,-167.4,-275.5 , -42.6, 0, -247.8,0, 0., 0, 0 ,0, 0,0 , 0 , -91.59 , 0, 10.52, 0.] HSA3[1..nhend] := HSA3[1..nhend]*PI/180 C OMEGA:=2*PI*60 C ENDINIT C Begin Routine EXEC IF (t=tstart) THEN C FOR TIME STEP ZERO WE CALCULATE THE G AND LINV FOR LOAD AND PHASE C OF FUNDAMENTAL VOLTAGE I:=0 WHILE I>> END MODEL HINJA C C USE H_INJ AS H_INJA C C SPECIFY PERCENTAGES C INPUT BUSVR[1] := BUSVR1 BUSVI[1] := BUSVI1 LODCR[1] := LODCR1 LODCI[1] := LODCI1 PER1[1] := 0.0 PER2[1] := 0.2 PER3[1] := 0.2 BUSVR[2] := BUSVR2 BUSVI[2] := BUSVI2 LODCR[2] := LODCR2 LODCI[2] := LODCI2 PER1[2] := 0.0 PER2[2] := 0.2 PER3[2] := 0.2 BUSVR[3] := BUSVR3 BUSVI[3] := BUSVI3 LODCR[3] := LODCR3 LODCI[3] := LODCI3 PER1[3] := 0.1 PER2[3] := 0.1 PER3[3] := 0. BUSVR[4] := BUSVR4 BUSVI[4] := BUSVI4 LODCR[4] := LODCR4 LODCI[4] := LODCI4 PER1[4] := 0.0 PER2[4] := 0.3 PER3[4] := 0. BUSVR[5] := BUSVR5 BUSVI[5] := BUSVI5 LODCR[5] := LODCR5 LODCI[5] := LODCI5 PER1[5] := 0.0 PER2[5] := 0.3 PER3[5] := 0. BUSVR[6] := BUSVR6 BUSVI[6] := BUSVI6 LODCR[6] := LODCR6 LODCI[6] := LODCI6 PER1[6] := 0.0 PER2[6] := 0.3 PER3[6] := 0. BUSVR[7] := BUSVR7 BUSVI[7] := BUSVI7 LODCR[7] := LODCR7 LODCI[7] := LODCI7 PER1[7] := 0.0 PER2[7] := 0.15 PER3[7] := 0.2 BUSVR[8] := BUSVR8 BUSVI[8] := BUSVI8 LODCR[8] := LODCR8 LODCI[8] := LODCI8 PER1[8] := 0.0 PER2[8] := 0.15 PER3[8] := 0.2 BUSVR[9] := BUSVR9 BUSVI[9] := BUSVI9 LODCR[9] := LODCR9 LODCI[9] := LODCI9 PER1[9] := 0.0 PER2[9] := 0.15 PER3[9] := 0.2 BUSVR[10] := BUSVR10 BUSVI[10] := BUSVI10 LODCR[10] := LODCR10 LODCI[10] := LODCI10 PER1[10] := 0.0 PER2[10] := 0.15 PER3[10] := 0.2 BUSVR[11] := BUSVR11 BUSVI[11] := BUSVI11 LODCR[11] := LODCR11 LODCI[11] := LODCI11 PER1[11] := 0.0 PER2[11] := 0.15 PER3[11] := 0.2 C OUTPUT HUS1 := HS[1] HUS2 := HS[2] HUS3 := HS[3] HUS4 := HS[4] HUS5 := HS[5] HUS6 := HS[6] HUS7 := HS[7] HUS8 := HS[8] HUS9 := HS[9] HUS10 := HS[10] HUS11 := HS[11] C ENDUSE RECORD H_INJA.HS[1] AS HS1 H_INJA.HS[2] AS HS2 C ------------------------------------------------------>> OUTPUT SECTION C C C ------------------------------------------------------>> CURRENT INJECTION C endmodels