Howdy!
I am trying to run the IEEE13 test case and i am wondering why a new 3 sets of single phase transformers were defined for the voltage regulator control. Can we just define one regulator control connected to the substation transformer of bus 650? I couldn't run the two scenarios to see the impact because my circuit is not running well. All the voltages are going to zero and i am still struggling to find out why?
Sophia.
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Hi,
I am a newbie but here is my 2cents.
So your two senarios are three phase regulator vs 3xsingle phase regulators?
Now the regulators control the taps to keep voltage at specified value of local or remote bus. Since the voltages of the three phases are different in typical distribution systemsit is not possible to have the same taps (one + seq control) on all three transformer windings to control the voltage. I believe that is why we use single phase regulators and independednt control for each winding. If you can post your code, I can have a look as to why you are getting 0 voltage.
Hasala
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Hi Hasala,
Paulo Ricardo Radatz de Freitas has responded and corrected the code i posted in the open forum. I have run it and i'm getting a solution. I am yet to study the solution to see whether it matches the posted solution by IEEE. Paulo's response to me is as follows:
Hello Sophia,
The issue is related in the way you defined the matrices of the linecodes when nphases is 1 or 2.
For example:
You did
New Linecode.605 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[0|0 0|0 0 1.3292] ! in ohms per mile
~ Xmatrix=[0|0 0|0 0 1.3475] ! in ohms per mile
~ Cmatrix=[0|0 0|0 0 11.988] ! in nF per mile.
It should be
New Linecode.605 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[1.3292] ! in ohms per mile
~ Xmatrix=[1.3475] ! in ohms per mile
~ Cmatrix=[11.988] ! in nF per mile.
You can't define those matrices with 3x3 dimension when nphases is not equal to 3.
The working code is below:
clear
New Circuit.IEEE13bus basekV = 115 bus1=SourceBus pu=1.0001 MVAsc3=20000 MVAsc1=21000 Angle=0 Frequency=60
// ------------------------------------Line Codes Definitions---------------------------------------------------------//
New Linecode.601 Nphases=3 basefreq=60 units=mi
~ Rmatrix=[0.3465|0.1560 0.3375|0.1580 0.1535 0.3414] ! in ohms per mile
~ Xmatrix=[1.0179|0.5017 1.0478|0.4236 0.3849 1.0348] ! in ohms per mile
~ Cmatrix=[16.711|-5.294 15.809|-3.341 -1.967 14.957] ! in nF per mile.
New Linecode.602 Nphases=3 basefreq=60 units=mi
~ Rmatrix=[0.7526|0.1580 0.7475|0.1560 0.1535 0.7436] ! in ohms per mile
~ Xmatrix=[1.1814|0.4236 1.1983|0.5017 0.3849 1.2112] ! in ohms per mile
~ Cmatrix=[5.6991|-1.0817 5.1795|-1.6905 -0.6588 5.4246] ! in nF per mile.
New Linecode.603 Nphases=2 basefreq=60 units=mi
~ Rmatrix=[1.3294|0.2066 1.3238] ! in ohms per mile
~ Xmatrix=[1.3471|0.4591 1.3569] ! in ohms per mile
~ Cmatrix=[12.493|-2.387 12.376] ! in nF per mile.
New Linecode.604 Nphases=2 basefreq=60 units=mi
~ Rmatrix=[1.3238|0.2066 1.3294] ! in ohms per mile
~ Xmatrix=[1.3569|0.4591 1.3471] ! in ohms per mile
~ Cmatrix=[12.376|-2.387 12.493] ! in nF per mile.
New Linecode.605 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[1.3292] ! in ohms per mile
~ Xmatrix=[1.3475] ! in ohms per mile
~ Cmatrix=[11.988] ! in nF per mile.
New Linecode.606 Nphases=3 basefreq=60 units=mi
~ Rmatrix=[0.7982|0.3192 0.7891|0.2849 0.3192 0.7982] ! in ohms per mile
~ Xmatrix=[0.4463|0.0328 0.4041|-0.0143 0.0328 0.4463] ! in ohms per mile
~ Cmatrix=[257.01|0 257.01|0 0 257.01] ! in nF per mile.
New Linecode.607 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[1.3425] ! in ohms per mile
~ Xmatrix=[0.5124] ! in ohms per mile
~ Cmatrix=[236.06] ! in nF per mile.
// --------------------Transformer Definitions---------------------------------------------------------------//
New Transformer.Sub phases=3 Windings=2 %r=0.001 XHL=0.008 %noloadloss=0 %imag=0 !This defines a substation transformer connected to bus 650
~ wdg=1 bus=SourceBus kV=115 kVA=5000 conn=delta
~ wdg=2 bus=650 kV=4.16 kVA=5000 conn=wye
New Transformer.XFM1 phases=3 windings=2 xhl=2 %r=1.1 %noloadloss=0 %imag=0 !This defines the transformer connected between node 633 and node 634
~ wdg=1 bus=633 kV=4.16 kVA=500 conn=wye
~ wdg=2 bus=634 kV=0.480 kVA=500 conn=wye
//--------------------Regulator Definitions-------------------------------------------------------------------------------//
New Transformer.Reg1 phases=1 bank=reg1 xhl=0.01 %loadloss=0.01
~ wdg=1 bus=650.1 conn=wye kVA=1666.67 kV=2.402
~ wdg=2 bus=REG.1 conn=wye kVA=1666.67 kV=2.402
New RegControl.regc1 Transformer=Reg1 Winding=2 Vreg=122 band=2 ptratio=20 CTprim=700 R=3 X=9
New Transformer.Reg2 phases=1 bank=reg1 xhl=0.01 %loadloss=0.01
~ wdg=1 bus=650.2 conn=wye kVA=1666.67 kV=2.402
~ wdg=2 bus=REG.2 conn=wye kVA=1666.67 kV=2.402
New RegControl.regc2 Transformer=Reg2 Winding=2 Vreg=122 band=2 ptratio=20 CTprim=700 R=3 X=9
New Transformer.Reg3 phases=1 bank=reg1 xhl=0.01 %loadloss=0.01
~ wdg=1 bus=650.3 conn=wye kVA=1666.67 kV=2.402
~ wdg=2 bus=REG.3 conn=wye kVA=1666.67 kV=2.402
New RegControl.regc3 Transformer=Reg3 Winding=2 Vreg=122 band=2 ptratio=20 CTprim=700 R=3 X=9
// --------------------Line Definitions ............................................................................................................................//
New Line.650632 phases=3 bus1=REG.2.1.3 bus2=632.2.1.3 length=2000 units=ft Linecode=601
New Line.646645 phases=2 bus1=645.3.2 bus2=646.3.2 length=300 units=ft Linecode=603
New Line.645632 phases=2 bus1=632.3.2 bus2=645.3.2 length=500 units=ft Linecode=603
New Line.632633 phases=3 bus1=632.3.1.2 bus2=633.3.1.2 length=500 units=ft Linecode=602
New Line.632671 phases=3 bus1=632.2.1.3 bus2=671.2.1.3 length=2000 units=ft Linecode=601
New Line.611684 phases=1 bus1=684.3 bus2=611.3 length=300 units=ft Linecode=605
New Line.684652 phases=1 bus1=684.1 bus2=652.1 length=800 units=ft Linecode=607
New Line.684671 phases=2 bus1=671.1.3 bus2=684.1.3 length=300 units=ft Linecode=604
New Line.671680 phases=3 bus1=671.2.1.3 bus2=680.2.1.3 length=1000 units=ft Linecode=601
New Line.692675 phases=3 bus1=692.1.2.3 bus2=675.1.2.3 length=500 units=ft Linecode=606
// --------------------Switch Definition---------------------------------------------------------------------//
New Line.671692 phases=3 bus1=671 bus2=692 Switch=Yes r1=1e-4 r0=1e-4 x1=0.000 x0=0.000 c1=0.000 c0=0.000
// --------------------Spot Load Definitions----------------------------------------------------------------------//
New Load.Load634A phases=1 bus1=634.1 kV=0.277 kW=160 kVAr=110 model=1 conn=wye
New Load.Load634B phases=1 bus1=634.2 kV=0.277 kW=120 kVAr=90 model=1 conn=wye
New Load.Load634C phases=1 bus1=634.3 kV=0.277 kW=120 kVAr=90 model=1 conn=wye
New Load.Load645 phases=1 bus1=645.2 kV=2.402 kW=170 kVAr=125 model=1 conn=wye
New Load.Load646 phases=1 bus1=646.2.3 kV=4.16 kW=230 kVAr=132 model=2 conn=delta
New Load.Load652 phases=1 bus1=652.1 kV=2.402 kW=128 kVAr=86 model=2 conn=wye
New Load.Load671 phases=3 bus1=671 kV=4.16 kW=1155 kVAr=660 model=1 conn=delta
New Load.Load675A phases=1 bus1=675.1 kV=2.402 kW=485 kVAr=190 model=1 conn=wye
New Load.Load675B phases=1 bus1=675.2 kV=2.402 kW=68 kVAr=60 model=1 conn=wye
New Load.Load675C phases=1 bus1=675.3 kV=2.402 kW=290 kVAr=212 model=1 conn=wye
New Load.Load692 phases=1 bus1=692.3.1 kV=4.16 kW=170 kVAr=151 model=5 conn=delta
New Load.Load611 phases=1 bus1=611.3 kV=2.402 kW=170 kVAr=80 model=5 conn=wye
// -------------------Distributed Load Definitions----------------------------------------------------------------//
! The load is split into two with half connected to one end of the line and the other half connected to the other end of the line
New Load.Load632A phases=1 bus1=632.1 kV=2.402 kW=8.5 kVAr=5 model=1 conn=wye
New Load.Load671A phases=1 bus1=671.1 kV=2.402 kW=8.5 kVAr=5 model=1 conn=wye
New Load.Load632B phases=1 bus1=632.2 kV=2.402 kW=33 kVAr=19 model=1 conn=wye
New Load.Load671B phases=1 bus1=671.2 kV=2.402 kW=33 kVAr=19 model=1 conn=wye
New Load.Load632C phases=1 bus1=632.3 kV=2.402 kW=58.5 kVAr=34 model=1 conn=wye
New Load.Load671C phases=1 bus1=671.3 kV=2.402 kW=58.5 kVAr=34 model=1 conn=wye
//---------------------Capacitor Definitions---------------------------------------------------------------------//
New Capacitor.Cap675 phases=3 bus1=675 bus2=675.0 kVAr=600 kV=4.16
New Capacitor.Cap611 phases=1 bus1=611.3 kVAr=100 kV=2.402
// ------------------------This section calculates the voltage bases and sets the solution mode-------------//
Set Voltagebases=[115 4.16 0.48]
CalcVoltageBases
Set mode=snapshot
New Energymeter.M1 Line.650632 ! This adds energy meter to monitor flow on the line
// ------------------------------- This section solves the circuit and show some results-----------------BusCoords C:\Users\satij\Documents\OpenDSS\685\Assumed_Cordinates.dss !This points to the file location of the assumed bus coordinates
solve
Plot profile phases=all
Show Voltage LN Nodes
Show Currents Elements
Show Powers kva Elements
Show taps ! shows the present regulator taps
Best Regards.
Paulo Radatz
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Howdy!
I am trying to run the IEEE13 test case and i am wondering why a new 3 sets of single phase transformers were defined for the voltage regulator control. Can we just define one regulator control connected to the substation transformer of bus 650? I couldn't run the two scenarios to see the impact because my circuit is not running well. All the voltages are going to zero and i am still struggling to find out why?
Sophia.
Hi,
I am a newbie but here is my 2cents.
So your two senarios are three phase regulator vs 3xsingle phase regulators?
Now the regulators control the taps to keep voltage at specified value of local or remote bus. Since the voltages of the three phases are different in typical distribution systemsit is not possible to have the same taps (one + seq control) on all three transformer windings to control the voltage. I believe that is why we use single phase regulators and independednt control for each winding. If you can post your code, I can have a look as to why you are getting 0 voltage.
Hasala
Hi Hasala,
Paulo Ricardo Radatz de Freitas has responded and corrected the code i posted in the open forum. I have run it and i'm getting a solution. I am yet to study the solution to see whether it matches the posted solution by IEEE. Paulo's response to me is as follows:
Hello Sophia,
The issue is related in the way you defined the matrices of the linecodes when nphases is 1 or 2.
For example:
You did
New Linecode.605 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[0|0 0|0 0 1.3292] ! in ohms per mile
~ Xmatrix=[0|0 0|0 0 1.3475] ! in ohms per mile
~ Cmatrix=[0|0 0|0 0 11.988] ! in nF per mile.
It should be
New Linecode.605 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[1.3292] ! in ohms per mile
~ Xmatrix=[1.3475] ! in ohms per mile
~ Cmatrix=[11.988] ! in nF per mile.
You can't define those matrices with 3x3 dimension when nphases is not equal to 3.
The working code is below:
clear
New Circuit.IEEE13bus basekV = 115 bus1=SourceBus pu=1.0001 MVAsc3=20000 MVAsc1=21000 Angle=0 Frequency=60
// ------------------------------------Line Codes Definitions---------------------------------------------------------//
New Linecode.601 Nphases=3 basefreq=60 units=mi
~ Rmatrix=[0.3465|0.1560 0.3375|0.1580 0.1535 0.3414] ! in ohms per mile
~ Xmatrix=[1.0179|0.5017 1.0478|0.4236 0.3849 1.0348] ! in ohms per mile
~ Cmatrix=[16.711|-5.294 15.809|-3.341 -1.967 14.957] ! in nF per mile.
New Linecode.602 Nphases=3 basefreq=60 units=mi
~ Rmatrix=[0.7526|0.1580 0.7475|0.1560 0.1535 0.7436] ! in ohms per mile
~ Xmatrix=[1.1814|0.4236 1.1983|0.5017 0.3849 1.2112] ! in ohms per mile
~ Cmatrix=[5.6991|-1.0817 5.1795|-1.6905 -0.6588 5.4246] ! in nF per mile.
New Linecode.603 Nphases=2 basefreq=60 units=mi
~ Rmatrix=[1.3294|0.2066 1.3238] ! in ohms per mile
~ Xmatrix=[1.3471|0.4591 1.3569] ! in ohms per mile
~ Cmatrix=[12.493|-2.387 12.376] ! in nF per mile.
New Linecode.604 Nphases=2 basefreq=60 units=mi
~ Rmatrix=[1.3238|0.2066 1.3294] ! in ohms per mile
~ Xmatrix=[1.3569|0.4591 1.3471] ! in ohms per mile
~ Cmatrix=[12.376|-2.387 12.493] ! in nF per mile.
New Linecode.605 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[1.3292] ! in ohms per mile
~ Xmatrix=[1.3475] ! in ohms per mile
~ Cmatrix=[11.988] ! in nF per mile.
New Linecode.606 Nphases=3 basefreq=60 units=mi
~ Rmatrix=[0.7982|0.3192 0.7891|0.2849 0.3192 0.7982] ! in ohms per mile
~ Xmatrix=[0.4463|0.0328 0.4041|-0.0143 0.0328 0.4463] ! in ohms per mile
~ Cmatrix=[257.01|0 257.01|0 0 257.01] ! in nF per mile.
New Linecode.607 Nphases=1 basefreq=60 units=mi
~ Rmatrix=[1.3425] ! in ohms per mile
~ Xmatrix=[0.5124] ! in ohms per mile
~ Cmatrix=[236.06] ! in nF per mile.
// --------------------Transformer Definitions---------------------------------------------------------------//
New Transformer.Sub phases=3 Windings=2 %r=0.001 XHL=0.008 %noloadloss=0 %imag=0 !This defines a substation transformer connected to bus 650
~ wdg=1 bus=SourceBus kV=115 kVA=5000 conn=delta
~ wdg=2 bus=650 kV=4.16 kVA=5000 conn=wye
New Transformer.XFM1 phases=3 windings=2 xhl=2 %r=1.1 %noloadloss=0 %imag=0 !This defines the transformer connected between node 633 and node 634
~ wdg=1 bus=633 kV=4.16 kVA=500 conn=wye
~ wdg=2 bus=634 kV=0.480 kVA=500 conn=wye
//--------------------Regulator Definitions-------------------------------------------------------------------------------//
New Transformer.Reg1 phases=1 bank=reg1 xhl=0.01 %loadloss=0.01
~ wdg=1 bus=650.1 conn=wye kVA=1666.67 kV=2.402
~ wdg=2 bus=REG.1 conn=wye kVA=1666.67 kV=2.402
New RegControl.regc1 Transformer=Reg1 Winding=2 Vreg=122 band=2 ptratio=20 CTprim=700 R=3 X=9
New Transformer.Reg2 phases=1 bank=reg1 xhl=0.01 %loadloss=0.01
~ wdg=1 bus=650.2 conn=wye kVA=1666.67 kV=2.402
~ wdg=2 bus=REG.2 conn=wye kVA=1666.67 kV=2.402
New RegControl.regc2 Transformer=Reg2 Winding=2 Vreg=122 band=2 ptratio=20 CTprim=700 R=3 X=9
New Transformer.Reg3 phases=1 bank=reg1 xhl=0.01 %loadloss=0.01
~ wdg=1 bus=650.3 conn=wye kVA=1666.67 kV=2.402
~ wdg=2 bus=REG.3 conn=wye kVA=1666.67 kV=2.402
New RegControl.regc3 Transformer=Reg3 Winding=2 Vreg=122 band=2 ptratio=20 CTprim=700 R=3 X=9
// --------------------Line Definitions ............................................................................................................................//
New Line.650632 phases=3 bus1=REG.2.1.3 bus2=632.2.1.3 length=2000 units=ft Linecode=601
New Line.646645 phases=2 bus1=645.3.2 bus2=646.3.2 length=300 units=ft Linecode=603
New Line.645632 phases=2 bus1=632.3.2 bus2=645.3.2 length=500 units=ft Linecode=603
New Line.632633 phases=3 bus1=632.3.1.2 bus2=633.3.1.2 length=500 units=ft Linecode=602
New Line.632671 phases=3 bus1=632.2.1.3 bus2=671.2.1.3 length=2000 units=ft Linecode=601
New Line.611684 phases=1 bus1=684.3 bus2=611.3 length=300 units=ft Linecode=605
New Line.684652 phases=1 bus1=684.1 bus2=652.1 length=800 units=ft Linecode=607
New Line.684671 phases=2 bus1=671.1.3 bus2=684.1.3 length=300 units=ft Linecode=604
New Line.671680 phases=3 bus1=671.2.1.3 bus2=680.2.1.3 length=1000 units=ft Linecode=601
New Line.692675 phases=3 bus1=692.1.2.3 bus2=675.1.2.3 length=500 units=ft Linecode=606
// --------------------Switch Definition---------------------------------------------------------------------//
New Line.671692 phases=3 bus1=671 bus2=692 Switch=Yes r1=1e-4 r0=1e-4 x1=0.000 x0=0.000 c1=0.000 c0=0.000
// --------------------Spot Load Definitions----------------------------------------------------------------------//
New Load.Load634A phases=1 bus1=634.1 kV=0.277 kW=160 kVAr=110 model=1 conn=wye
New Load.Load634B phases=1 bus1=634.2 kV=0.277 kW=120 kVAr=90 model=1 conn=wye
New Load.Load634C phases=1 bus1=634.3 kV=0.277 kW=120 kVAr=90 model=1 conn=wye
New Load.Load645 phases=1 bus1=645.2 kV=2.402 kW=170 kVAr=125 model=1 conn=wye
New Load.Load646 phases=1 bus1=646.2.3 kV=4.16 kW=230 kVAr=132 model=2 conn=delta
New Load.Load652 phases=1 bus1=652.1 kV=2.402 kW=128 kVAr=86 model=2 conn=wye
New Load.Load671 phases=3 bus1=671 kV=4.16 kW=1155 kVAr=660 model=1 conn=delta
New Load.Load675A phases=1 bus1=675.1 kV=2.402 kW=485 kVAr=190 model=1 conn=wye
New Load.Load675B phases=1 bus1=675.2 kV=2.402 kW=68 kVAr=60 model=1 conn=wye
New Load.Load675C phases=1 bus1=675.3 kV=2.402 kW=290 kVAr=212 model=1 conn=wye
New Load.Load692 phases=1 bus1=692.3.1 kV=4.16 kW=170 kVAr=151 model=5 conn=delta
New Load.Load611 phases=1 bus1=611.3 kV=2.402 kW=170 kVAr=80 model=5 conn=wye
// -------------------Distributed Load Definitions----------------------------------------------------------------//
! The load is split into two with half connected to one end of the line and the other half connected to the other end of the line
New Load.Load632A phases=1 bus1=632.1 kV=2.402 kW=8.5 kVAr=5 model=1 conn=wye
New Load.Load671A phases=1 bus1=671.1 kV=2.402 kW=8.5 kVAr=5 model=1 conn=wye
New Load.Load632B phases=1 bus1=632.2 kV=2.402 kW=33 kVAr=19 model=1 conn=wye
New Load.Load671B phases=1 bus1=671.2 kV=2.402 kW=33 kVAr=19 model=1 conn=wye
New Load.Load632C phases=1 bus1=632.3 kV=2.402 kW=58.5 kVAr=34 model=1 conn=wye
New Load.Load671C phases=1 bus1=671.3 kV=2.402 kW=58.5 kVAr=34 model=1 conn=wye
//---------------------Capacitor Definitions---------------------------------------------------------------------//
New Capacitor.Cap675 phases=3 bus1=675 bus2=675.0 kVAr=600 kV=4.16
New Capacitor.Cap611 phases=1 bus1=611.3 kVAr=100 kV=2.402
// ------------------------This section calculates the voltage bases and sets the solution mode-------------//
Set Voltagebases=[115 4.16 0.48]
CalcVoltageBases
Set mode=snapshot
New Energymeter.M1 Line.650632 ! This adds energy meter to monitor flow on the line
// ------------------------------- This section solves the circuit and show some results-----------------BusCoords C:\Users\satij\Documents\OpenDSS\685\Assumed_Cordinates.dss !This points to the file location of the assumed bus coordinates
solve
Plot profile phases=all
Show Voltage LN Nodes
Show Currents Elements
Show Powers kva Elements
Show taps ! shows the present regulator taps
Best Regards.
Paulo Radatz
That is fantastic!