Below is my updated DSS script, I am wondering whether my script logic is correct, I have been consistently updating the script to ascertain whether I have done it correctly.
Kind Regards,
Georgin
clear
set defaultbasefreq = 50
New object=circuit.Kaiwaka.
basekv=11 pu =1.0 bus1=sourcebus
I believe the problem is in the WireData definitions. By running your script and using the command FormEdit "WireData.11KVOH" for example, you can see that there are no default values for Rac, GMRac and the same applies for the units (Runits, GMRunits, radunits). So you have to specify them.
In the following video there is a complete example about how to define a line object based on its geometry and wiredata.
Thankyou Celso for your time and effort, I have made the adjustments as per your instructions. I am now getting very low losses, less than 1%. I believe I have defined everything correctly now and that this could possibly be fine?
Best regards,
Georgin
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
I have noticed that the LV lines are being neglected when I draw the voltage profile, I seem to just get dashed lines, is it possible to display the LV side of the network on a voltage plot?
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
I noticed that there are several loops in your LV network. If you use the command "Show Loops", you can verify that. All loops start right after line L7. For example, there is one loop between the following elements
New Transformer.TR0 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB8.1.2.3 conn=Wye kv=0.410 kva=150
New Transformer.TR1 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB9.1.2.3 conn=Wye kv=0.410 kva=150
New line.L9 LB8.1.2.3 LB9.1.2.3 geometry=lineC length=40 Units = m !LV
Is that what you really want to do? That's why you can barely see the LV lines in the voltage profile. They are all concentrated after line L7 (the last 3-phase line). The energy meter is getting confused with all those loops.
For example, try to comment all transformers except the first one(TR0). You will see the difference in the voltage profile and in the circuit losses!
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Hey Celso, this is what I have been stuck on for a while now :( I just wanted to define one transformer to deliver 230V to all my loads, but when I do this I get very high losses(Commenting out all transformers but TR0). Thankyou so much for your responses.
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
I have updated the script according to our suggestion Celso, I have already tried this before and got quite high losses, around 12% when only using 1 transfomer, but perhaps this is okay?
Last edit: Georgin Raju 2017-08-28
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Hey everyone,
Below is my updated DSS script, I am wondering whether my script logic is correct, I have been consistently updating the script to ascertain whether I have done it correctly.
Kind Regards,
Georgin
clear
set defaultbasefreq = 50
New object=circuit.Kaiwaka.
basekv=11 pu =1.0 bus1=sourcebus
New loadshape.day 12 2.0
~ mult=(.44 .38 .36 .50 .88 .67 .56 .49 .55 .93 .89 .66 .42)
New loadshape.EV 24 1.0
~ mult=(.3 .3 .3 .35 .36 .39 .41 .48 .52 .59 .62 .94 .87 .91 .95 .95 1.0 .98 .94 .92 .61 .60 .51 .44)
//Substation transformer connections
New Transformer.TR0 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB8.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR1 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB9.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR2 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB10.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR3 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB11.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR4 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB12.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR5 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB13.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR6 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB14.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR7 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB15.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR8 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB16.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR9 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB17.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR10 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB18.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR11 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB19.1.2.3 conn=Wye kv=0.410 kva=150
//New Transformer.TR12 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB20.1.2.3 conn=Wye kv=0.410 kva=150
//Wiredata
New Wiredata.11kVOH Radius=0.181 RDC= 0.524 GMR=0.144619
~ Runits=km radunits=m gmrunits=m
New Wiredata.400VOH Radius =0.071 RDC= 0.234 GMR=0.056729
~ Runits=km radunits=m gmrunits=m
New Wiredata.Neutral Radius =0.1475 RDC= 0.524 GMR=0.117853
~ Runits=km radunits=m gmrunits=m
//FormEdit "WireData.11KVOH"
//Linegeometry
New Linegeometry.LineB nconds =4 nphases=3 reduce = Yes
~ cond=1 Wire=11kVOH x=-1.5 h=9.2 units=m
~ cond=2 Wire=11kVOH x=1.5 h=9.2 units=m
~ cond=3 Wire=11kVOH x=0.6 h=9.2 units=m
~ cond=4 Wire=Neutral x=-0.6 h=9.2 units=m
New Linegeometry.LineC nconds =4 nphases=3 reduce = Yes
~ cond=1 Wire=400VOH x=-1.5 h=8.2 units=m
~ cond=2 Wire=400VOH x=1.5 h=8.2 units=m
~ cond=3 Wire=400VOH x=0.6 h=8.2 units=m
~ cond=4 Wire=Neutral x=-0.6 h=8.2 units=m
New LineCode.lineB nphases=3 R1=0.5755 X1=0.19745 R0=0.723544 X0=2.50369 C1=0 C0=0 Units=km !11kV
New LineCode.lineD nphases=3 R1=0.361 X1=0.12798 R0=0.509044 X0=2.57317 C1=0 C0=0 Units=km !ServiceCable
//Connect to the LV side of the transformer
// Define the lines, line lengths based off to scale diagram
New line.L1 bus1=sourcebus bus2=LB1.1.2.3 linecode= lineB length =35 Units = m
New line.L2 LB1.1.2.3 LB2.1.2.3 geometry=lineB length=67 Units = m !MV
New line.L3 LB2.1.2.3 LB3.1.2.3 geometry=lineB length=43 Units = m !MV
New line.L4 LB3.1.2.3 LB4.1.2.3 geometry=lineB length=37 Units = m !MV
New line.L5 LB4.1.2.3 LB5.1.2.3 geometry=lineB length=51 Units = m !MV
New line.L6 LB5.1.2.3 LB6.1.2.3 geometry=lineB length=51 Units = m !MV
New line.L7 LB6.1.2.3 LB7.1.2.3 geometry=lineB length=26 Units = m !MV
//MV feeder continues
New line.L9 LB8.1.2.3 LB9.1.2.3 geometry=lineC length=40 Units = m !LV
New line.L10 LB9.1.2.3 LB10.1.2.3 geometry=lineC length=37 Units = m !LV
New line.L11 LB10.1.2.3 LB11.1.2.3 linecode=lineD length=12 Units = m !LV
New line.L12 LB11.1.2.3 LB12.1.2.3 geometry=lineC length=28 Units = m !LV
New line.L13 LB12.1.2.3 LB13.1.2.3 geometry=lineC length=30 Units = m !LV
New line.L14 LB13.1.2.3 LB14.1.2.3 geometry=lineC length=26 Units = m !LV
New line.L15 LB8.1.2.3 LB15.1.2.3 geometry=lineC length=26 Units = m !LV
New line.L16 LB15.1.2.3 LB16.1.2.3 geometry=lineC length=51 Units = m !LV
New line.L17 LB16.1.2.3 LB17.1.2.3 geometry=lineC length=51 Units = m !LV
New line.L18 LB17.1.2.3 LB18.1.2.3 geometry=lineC length=37 Units = m !LV
New line.L19 LB18.1.2.3 LB19.1.2.3 geometry=lineC length=43 Units = m !LV
New line.L20 LB19.1.2.3 LB20.1.2.3 geometry=lineC length=67 Units = m !LV
// Define the loads, all are residential loads taking 1 phase and at 230V
New load.load37 bus1=LB20.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load39 bus1=LB20.2 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load36 bus1=LB19.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load26 bus1=LB18.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load28 bus1=LB18.2 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load25 bus1=LB17.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load22 bus1=LB17.2 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load24 bus1=LB17.3 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load20 bus1=LB16.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load16 bus1=LB15.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load18 bus1=LB15.2 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load7 bus1=LB8.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load3 bus1=LB9.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load5 bus1=LB9.2 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load6 bus1=LB9.3 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load10 bus1=LB10.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load2 bus1=LB10.2 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load1 bus1=LB12.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load8 bus1=LB13.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New load.load11 bus1=LB14.1 phases=1 kv=0.237 kw=23.23 pf=0.99 model=1 daily = day
New energymeter.m1 LINE.L2 1
//Define EV load with it's transformer here
New load.EVload bus1=LB11.1.2.3 phases=3 kv=0.410 kw=50 pf=0.98 model=1 daily =EV
set mode=daily
set stepsize = 1h
set number = 24
set voltagebases=[11.0 0.410] !required to arrive at per-unit values
calcvoltage bases
//Solve circuit
Solve
Show Currents [residual=yes] [Elements]
Show Powers kVA elements
Show Voltages [LN] [Nodes]
Show Losses
plot profile ph=all
Last edit: Georgin Raju 2017-08-28
The losses are all from L2-L7, I wonder if I am doing something wrong when defining them.
It seems to be the way I am defining the line's, I think perhaps I have used the linegeometry function incorrectly?
Hi Georgin,
I believe the problem is in the WireData definitions. By running your script and using the command FormEdit "WireData.11KVOH" for example, you can see that there are no default values for Rac, GMRac and the same applies for the units (Runits, GMRunits, radunits). So you have to specify them.
In the following video there is a complete example about how to define a line object based on its geometry and wiredata.
https://www.youtube.com/watch?v=s8c9xq8IxUQ&list=PLcOap2oqW_gEMEVH9dg2HoXJ4NvydfsZM&index=9
I hope it helps!
Thankyou Celso for your time and effort, I have made the adjustments as per your instructions. I am now getting very low losses, less than 1%. I believe I have defined everything correctly now and that this could possibly be fine?
Best regards,
Georgin
Note that I have adjusted the script above to show the updated script
I have noticed that the LV lines are being neglected when I draw the voltage profile, I seem to just get dashed lines, is it possible to display the LV side of the network on a voltage plot?
Hi Georgin,
I noticed that there are several loops in your LV network. If you use the command "Show Loops", you can verify that. All loops start right after line L7. For example, there is one loop between the following elements
New Transformer.TR0 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB8.1.2.3 conn=Wye kv=0.410 kva=150
New Transformer.TR1 Phases=3 Windings=2 XHL=4 BaseFreq=50
~ wdg=1 bus=LB7.1.2.3 conn=Delta kv=11 kva =150
~ wdg=2 bus=LB9.1.2.3 conn=Wye kv=0.410 kva=150
New line.L9 LB8.1.2.3 LB9.1.2.3 geometry=lineC length=40 Units = m !LV
Is that what you really want to do? That's why you can barely see the LV lines in the voltage profile. They are all concentrated after line L7 (the last 3-phase line). The energy meter is getting confused with all those loops.
For example, try to comment all transformers except the first one(TR0). You will see the difference in the voltage profile and in the circuit losses!
Hey Celso, this is what I have been stuck on for a while now :( I just wanted to define one transformer to deliver 230V to all my loads, but when I do this I get very high losses(Commenting out all transformers but TR0). Thankyou so much for your responses.
I have updated the script according to our suggestion Celso, I have already tried this before and got quite high losses, around 12% when only using 1 transfomer, but perhaps this is okay?
Last edit: Georgin Raju 2017-08-28