monitor the IEEE 30 Bus for each element

[Cheng Wei Lin]

Dear all
I am trying to run the IEEE 30 bus and want to connect the PV system on the bus 30.Now I think I have already completed my script.Then I want to analyze the condition for each element.But when I monitor the generator.G2,In principle it have to produce 40000kW,but the result show it produce 960000kW.It is totally exceed the limited value that I desired.So I think I may did something wrong in my script.As follow is my script.

---

Clear
New Circuit.IEEE_30bus BasekV=132 Bus1=B1.1.2.3 pu=1.06 MVASC3=5000000 MVASC1=5000000

! Line definitions
New line.1-2 Bus1=B1.1.2.3 Bus2=B2.1.2.3 R1=3.345408 X1=10.0188 R0=10.036224 X0=30.0564 C1=803.793907242183 C0=267.931302414061 Length=1
New line.1-3 Bus1=B1.1.2.3 Bus2=B3.1.2.3 R1=7.875648 X1=28.784448 R0=23.626944 X0=86.353344 C1=621.113473778051 C0=207.037824592684 Length=1
New line.2-4 Bus1=B2.1.2.3 Bus2=B4.1.2.3 R1=9.93168 X1=30.265488 R0=29.79504 X0=90.796464 C1=560.219995956673 C0=186.739998652224 Length=1
New line.3-4 Bus1=B3.1.2.3 Bus2=B4.1.2.3 R1=2.299968 X1=6.603696 R0=6.899904 X0=19.811088 C1=127.876303424893 C0=42.6254344749642 Length=1
New line.2-5 Bus1=B2.1.2.3 Bus2=B5.1.2.3 R1=8.224128 X1=34.551792 R0=24.672384 X0=103.655376 C1=636.336843233395 C0=212.112281077798 Length=1
New line.2-6 Bus1=B2.1.2.3 Bus2=B6.1.2.3 R1=10.123344 X1=30.718512 R0=30.370032 X0=92.155536 C1=569.35401762988 C0=189.784672543293 Length=1
New line.4-6 Bus1=B4.1.2.3 Bus2=B6.1.2.3 R1=2.073456 X1=7.213536 R0=6.220368 X0=21.640608 C1=137.010325098099 C0=45.6701083660331 Length=1
New line.5-7 Bus1=B5.1.2.3 Bus2=B7.1.2.3 R1=8.01504 X1=20.21184 R0=24.04512 X0=60.63552 C1=310.556736889025 C0=103.518912296342 Length=1
New line.6-7 Bus1=B6.1.2.3 Bus2=B7.1.2.3 R1=4.652208 X1=14.28768 R0=13.956624 X0=42.86304 C1=258.797280740854 C0=86.2657602469515 Length=1
New line.6-8 Bus1=B6.1.2.3 Bus2=B8.1.2.3 R1=2.09088 X1=7.31808 R0=6.27264 X0=21.95424 C1=137.010325098099 C0=45.6701083660331 Length=1
New line.12-14 Bus1=B12.1.2.3 Bus2=B14.1.2.3 R1=1.340559 X1=2.786751 R0=4.021677 X0=8.360253 C1=0 C0=0 Length=1
New line.12-15 Bus1=B12.1.2.3 Bus2=B15.1.2.3 R1=0.720918 X1=1.420056 R0=2.162754 X0=4.260168 C1=0 C0=0 Length=1
New line.12-16 Bus1=B12.1.2.3 Bus2=B16.1.2.3 R1=1.029105 X1=2.163843 R0=3.087315 X0=6.491529 C1=0 C0=0 Length=1
New line.14-15 Bus1=B14.1.2.3 Bus2=B15.1.2.3 R1=2.40669 X1=2.174733 R0=7.22007 X0=6.524199 C1=0 C0=0 Length=1
New line.16-17 Bus1=B16.1.2.3 Bus2=B17.1.2.3 R1=0.570636 X1=2.094147 R0=1.711908 X0=6.282441 C1=0 C0=0 Length=1
New line.15-18 Bus1=B15.1.2.3 Bus2=B18.1.2.3 R1=1.168497 X1=2.379465 R0=3.505491 X0=7.138395 C1=0 C0=0 Length=1
New line.18-19 Bus1=B18.1.2.3 Bus2=B19.1.2.3 R1=0.695871 X1=1.406988 R0=2.087613 X0=4.220964 C1=0 C0=0 Length=1
New line.19-20 Bus1=B19.1.2.3 Bus2=B20.1.2.3 R1=0.37026 X1=0.74052 R0=1.11078 X0=2.22156 C1=0 C0=0 Length=1
New line.10-20 Bus1=B10.1.2.3 Bus2=B20.1.2.3 R1=1.019304 X1=2.27601 R0=3.057912 X0=6.82803 C1=0 C0=0 Length=1
New line.10-17 Bus1=B10.1.2.3 Bus2=B17.1.2.3 R1=0.352836 X1=0.920205 R0=1.058508 X0=2.760615 C1=0 C0=0 Length=1
New line.10-21 Bus1=B10.1.2.3 Bus2=B21.1.2.3 R1=0.378972 X1=0.815661 R0=1.136916 X0=2.446983 C1=0 C0=0 Length=1
New line.10-22 Bus1=B10.1.2.3 Bus2=B22.1.2.3 R1=0.791703 X1=1.632411 R0=2.375109 X0=4.897233 C1=0 C0=0 Length=1
New line.21-22 Bus1=B21.1.2.3 Bus2=B22.1.2.3 R1=0.126324 X1=0.257004 R0=0.378972 X0=0.771012 C1=0 C0=0 Length=1
New line.15-23 Bus1=B15.1.2.3 Bus2=B23.1.2.3 R1=1.089 X1=2.19978 R0=3.267 X0=6.59934 C1=0 C0=0 Length=1
New line.22-24 Bus1=B22.1.2.3 Bus2=B24.1.2.3 R1=1.25235 X1=1.94931 R0=3.75705 X0=5.84793 C1=0 C0=0 Length=1
New line.23-24 Bus1=B23.1.2.3 Bus2=B24.1.2.3 R1=1.43748 X1=2.9403 R0=4.31244 X0=8.8209 C1=0 C0=0 Length=1
New line.24-25 Bus1=B24.1.2.3 Bus2=B25.1.2.3 R1=2.052765 X1=3.584988 R0=6.158295 X0=10.754964 C1=0 C0=0 Length=1
New line.25-26 Bus1=B25.1.2.3 Bus2=B26.1.2.3 R1=2.770416 X1=4.1382 R0=8.311248 X0=12.4146 C1=0 C0=0 Length=1
New line.25-27 Bus1=B25.1.2.3 Bus2=B27.1.2.3 R1=1.190277 X1=2.272743 R0=3.570831 X0=6.818229 C1=0 C0=0 Length=1
New line.27-29 Bus1=B27.1.2.3 Bus2=B29.1.2.3 R1=2.393622 X1=4.522617 R0=7.180866 X0=13.567851 C1=0 C0=0 Length=1
New line.27-30 Bus1=B27.1.2.3 Bus2=B30.1.2.3 R1=3.486978 X1=6.563403 R0=10.460934 X0=19.690209 C1=0 C0=0 Length=1
New line.29-30 Bus1=B29.1.2.3 Bus2=B30.1.2.3 R1=2.612511 X1=4.936437 R0=7.837533 X0=14.809311 C1=0 C0=0 Length=1
New line.8-28 Bus1=B8.1.2.3 Bus2=B28.1.2.3 R1=11.081664 X1=34.848 R0=33.244992 X0=104.544 C1=651.560212688739 C0=217.186737562913 Length=1
New line.6-28 Bus1=B6.1.2.3 Bus2=B28.1.2.3 R1=2.944656 X1=10.436976 R0=8.833968 X0=31.310928 C1=197.903802919477 C0=65.9679343064923 Length=1

! Transformer definitions
New Transformer.6-9 kVAs=[100000 100000] XHL=20.8 PPM=0
~ Wdg=1 R=0 kV=132 Bus=B6.1.2.3 Tap=0.978
~ Wdg=2 R=0 kV=1 Bus=B9.1.2.3
~ %loadloss=0
New Transformer.6-10 kVAs=[100000 100000] XHL=55.6 PPM=0
~ Wdg=1 R=0 kV=132 Bus=B6.1.2.3 Tap=0.969
~ Wdg=2 R=0 kV=33 Bus=B10.1.2.3
~ %loadloss=0
New Transformer.9-11 kVAs=[100000 100000] XHL=20.8 PPM=0
~ Wdg=1 R=0 kV=1 Bus=B9.1.2.3 Tap=1
~ Wdg=2 R=0 kV=11 Bus=B11.1.2.3
~ %loadloss=0
New Transformer.9-10 kVAs=[100000 100000] XHL=11 PPM=0
~ Wdg=1 R=0 kV=1 Bus=B9.1.2.3 Tap=1
~ Wdg=2 R=0 kV=33 Bus=B10.1.2.3
~ %loadloss=0
New Transformer.4-12 kVAs=[100000 100000] XHL=25.6 PPM=0
~ Wdg=1 R=0 kV=132 Bus=B4.1.2.3 Tap=0.932
~ Wdg=2 R=0 kV=33 Bus=B12.1.2.3
~ %loadloss=0
New Transformer.12-13 kVAs=[100000 100000] XHL=14 PPM=0
~ Wdg=1 R=0 kV=33 Bus=B12.1.2.3 Tap=1
~ Wdg=2 R=0 kV=11 Bus=B13.1.2.3
~ %loadloss=0
New Transformer.28-27 kVAs=[100000 100000] XHL=39.6 PPM=0
~ Wdg=1 R=0 kV=132 Bus=B28.1.2.3 Tap=0.968
~ Wdg=2 R=0 kV=33 Bus=B27.1.2.3
~ %loadloss=0

! Load definition
New Loadshape.MyLoad npts=24 interval=1 mult=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1]
New Load.L2 Bus1=B2.1.2.3 kV=132 kW=21700 kvar=12700 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L3 Bus1=B3.1.2.3 kV=132 kW=2400 kvar=1200 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L4 Bus1=B4.1.2.3 kV=132 kW=7600 kvar=1600 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L5 Bus1=B5.1.2.3 kV=132 kW=94200 kvar=19000 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L7 Bus1=B7.1.2.3 kV=132 kW=22800 kvar=10900 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L8 Bus1=B8.1.2.3 kV=132 kW=30000 kvar=30000 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L10 Bus1=B10.1.2.3 kV=33 kW=5800 kvar=2000 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L12 Bus1=B12.1.2.3 kV=33 kW=11200 kvar=7500 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L14 Bus1=B14.1.2.3 kV=33 kW=6200 kvar=1600 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L15 Bus1=B15.1.2.3 kV=33 kW=8200 kvar=2500 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L16 Bus1=B16.1.2.3 kV=33 kW=3500 kvar=1800 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L17 Bus1=B17.1.2.3 kV=33 kW=9000 kvar=5800 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L18 Bus1=B18.1.2.3 kV=33 kW=3200 kvar=900 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L19 Bus1=B19.1.2.3 kV=33 kW=9500 kvar=3400 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L20 Bus1=B20.1.2.3 kV=33 kW=2200 kvar=700 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L21 Bus1=B21.1.2.3 kV=33 kW=17500 kvar=11200 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L23 Bus1=B23.1.2.3 kV=33 kW=3200 kvar=1600 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L24 Bus1=B24.1.2.3 kV=33 kW=8700 kvar=6700 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L26 Bus1=B26.1.2.3 kV=33 kW=3500 kvar=2300 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L29 Bus1=B29.1.2.3 kV=33 kW=2400 kvar=900 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad
New Load.L30 Bus1=B30.1.2.3 kV=33 kW=10600 kvar=1900 vminpu=0.9 vmaxpu=1.10 Daily=MyLoad

! Shunt Capacitor definitions
New Capacitor.C10 Bus1=B10.1.2.3 kV= 33 kvar=19000
New Capacitor.C24 Bus1=B24.1.2.3 kV= 33 kvar=4300

! Generator Definitions
New Generator.G2 Bus1=B2.1.2.3 kV= 132 kW=40000 Model=3 Vpu=1.045 Maxkvar=50000 Minkvar=-40000
New Generator.G5 Bus1=B5.1.2.3 kV= 132 kW=1 Model=3 Vpu=1.01 Maxkvar=40000 Minkvar=-40000
New Generator.G8 Bus1=B8.1.2.3 kV= 132 kW=1 Model=3 Vpu=1.01 Maxkvar=40000 Minkvar=-10000
New Generator.G11 Bus1=B11.1.2.3 kV= 11 kW=1 Model=3 Vpu=1.082 Maxkvar=24000 Minkvar=-6000
New Generator.G13 Bus1=B13.1.2.3 kV= 11 kW=1 Model=3 Vpu=1.071 Maxkvar=24000 Minkvar=-6000

! PV system definitions
// P-T curve is per unit of rated Pmpp vs temperature
// This one is for a Pmpp stated at 25 deg
New XYCurve.MyPvsT npts=4 xarray=[0 25 75 100] yarray=[1.2 1.0 0.8 0.6]
// efficiency curve is per unit eff vs per unit power
New XYCurve.MyEff npts=4 xarray=[.1 .2 .4 1.0] yarray=[.86 .9 .93 .97]
// per unit irradiance curve (per unit if "irradiance" property)
New Loadshape.MyIrrad npts=24 interval=1 mult=[0 0 0 0 0 0 0.1 0.2 0.3 0.5 0.8 0.9 1.0 1.0 0.99 0.9 0.7 0.4 0.1 0 0 0 0 0]
// 24-hr temp shape curve
New Tshape.MyTemp npts=24 interval=1 temp=[25 25 25 25 25 25 25 25 35 40 45 50 60 60 55 40 35 30 25 25 25 25 25 25]
// **** plot tshape object=mytemp
// take the default line
// pv definition
New PVSystem.PV phases=3 bus1=B30.1.2.3 kV=33 kVA=500 irrad=0.8 Pmpp=500
~ temperature=25 PF=1 effcurve=Myeff P-TCurve=MyPvsT
~ Daily=MyIrrad TDaily=MyTemp

! Let the openDSS estimate the voltage bases
Set Voltagebases=[132, 33, 11, 1]
Calcvoltagebases
Set maxiterations=100
solve

! Moniter definitions
new monitor.m1 element=PVSystem.PV terminal=1 mode=1 ppolar=no
new monitor.m2 element=PVSystem.PV terminal=1 mode=0
new monitor.m3 element=PVSystem.PV terminal=1 mode=3
new monitor.m4 element=Load.L2 terminal=1 mode=1 ppolar=no
new monitor.m5 element=Load.L3 terminal=1 mode=1 ppolar=no
new monitor.m6 element=Load.L4 terminal=1 mode=1 ppolar=no
new monitor.m7 element=Load.L5 terminal=1 mode=1 ppolar=no
new monitor.m8 element=Load.L7 terminal=1 mode=1 ppolar=no
new monitor.m9 element=Load.L8 terminal=1 mode=1 ppolar=no
new monitor.m10 element=Load.L10 terminal=1 mode=1 ppolar=no
new monitor.m11 element=Load.L12 terminal=1 mode=1 ppolar=no
new monitor.m12 element=Load.L14 terminal=1 mode=1 ppolar=no
new monitor.m13 element=Load.L15 terminal=1 mode=1 ppolar=no
new monitor.m14 element=Load.L16 terminal=1 mode=1 ppolar=no
new monitor.m15 element=Load.L17 terminal=1 mode=1 ppolar=no
new monitor.m16 element=Load.L18 terminal=1 mode=1 ppolar=no
new monitor.m17 element=Load.L19 terminal=1 mode=1 ppolar=no
new monitor.m18 element=Load.L20 terminal=1 mode=1 ppolar=no
new monitor.m19 element=Load.L21 terminal=1 mode=1 ppolar=no
new monitor.m20 element=Load.L23 terminal=1 mode=1 ppolar=no
new monitor.m21 element=Load.L24 terminal=1 mode=1 ppolar=no
new monitor.m22 element=Load.L26 terminal=1 mode=1 ppolar=no
new monitor.m23 element=Load.L29 terminal=1 mode=1 ppolar=no
new monitor.m24 element=Load.L30 terminal=1 mode=1 ppolar=no
new monitor.m25 element=Line.1-2 terminal=1 mode=1 ppolar=no
new monitor.m26 element=Generator.G2 terminal=1 mode=1 ppolar=no

solve
solve mode=daily

show mon m1
show mon m2
show mon m3
show mon m4
show mon m5
show mon m6
show mon m7
show mon m8
show mon m9
show mon m10
show mon m11
show mon m12
show mon m13
show mon m14
show mon m15
show mon m16
show mon m17
show mon m18
show mon m19
show mon m20
show mon m21
show mon m22
show mon m23
show mon m24
show mon m25
show mon m26
Export monitors m1
Export monitors m2
Export monitors m3
Export monitors m4
Export monitors m5
Export monitors m6
Export monitors m7
Export monitors m8
Export monitors m9
Export monitors m10
Export monitors m11
Export monitors m12
Export monitors m13
Export monitors m14
Export monitors m15
Export monitors m16
Export monitors m17
Export monitors m18
Export monitors m19
Export monitors m20
Export monitors m21
Export monitors m22
Export monitors m23
Export monitors m24
Export monitors m25
Export monitors m25

Plot monitor object= m1 channels=(1,3,5 )
Plot monitor object= m1 channels=(2,4,6 )
Plot monitor object= m2 channels=(1,3,5 )
Plot monitor object= m2 channels=(9,11,13 )
Plot monitor object= m3 channels=(1 )
Plot monitor object= m3 channels=(2 )
Plot monitor object= m3 channels=(3 )
Plot monitor object= m4 channels=(1,3,5 )
Plot monitor object= m4 channels=(2,4,6 )
Plot monitor object= m5 channels=(1,3,5 )
Plot monitor object= m5 channels=(2,4,6 )
Plot monitor object= m6 channels=(1,3,5 )
Plot monitor object= m6 channels=(2,4,6 )
Plot monitor object= m7 channels=(1,3,5 )
Plot monitor object= m7 channels=(2,4,6 )
Plot monitor object= m8 channels=(1,3,5 )
Plot monitor object= m8 channels=(2,4,6 )
Plot monitor object= m9 channels=(1,3,5 )
Plot monitor object= m9 channels=(2,4,6 )
Plot monitor object= m10 channels=(1,3,5 )
Plot monitor object= m10 channels=(2,4,6 )
Plot monitor object= m11 channels=(1,3,5 )
Plot monitor object= m11 channels=(2,4,6 )
Plot monitor object= m12 channels=(1,3,5 )
Plot monitor object= m12 channels=(2,4,6 )
Plot monitor object= m13 channels=(1,3,5 )
Plot monitor object= m13 channels=(2,4,6 )
Plot monitor object= m14 channels=(1,3,5 )
Plot monitor object= m14 channels=(2,4,6 )
Plot monitor object= m15 channels=(1,3,5 )
Plot monitor object= m15 channels=(2,4,6 )
Plot monitor object= m16 channels=(1,3,5 )
Plot monitor object= m16 channels=(2,4,6 )
Plot monitor object= m17 channels=(1,3,5 )
Plot monitor object= m17 channels=(2,4,6 )
Plot monitor object= m18 channels=(1,3,5 )
Plot monitor object= m18 channels=(2,4,6 )
Plot monitor object= m19 channels=(1,3,5 )
Plot monitor object= m19 channels=(2,4,6 )
Plot monitor object= m20 channels=(1,3,5 )
Plot monitor object= m20 channels=(2,4,6 )
Plot monitor object= m21 channels=(1,3,5 )
Plot monitor object= m21 channels=(2,4,6 )
Plot monitor object= m22 channels=(1,3,5 )
Plot monitor object= m22 channels=(2,4,6 )
Plot monitor object= m23 channels=(1,3,5 )
Plot monitor object= m23 channels=(2,4,6 )
Plot monitor object= m24 channels=(1,3,5 )
Plot monitor object= m24 channels=(2,4,6 )
Plot monitor object= m25 channels=(1,3,5 )
Plot monitor object= m25 channels=(2,4,6 )
Plot monitor object= m26 channels=(1,3,5 )
Plot monitor object= m26 channels=(2,4,6 )

Thanks

C.W Lin

[wsunderm1]

C.W Lin,

Thank you for posting your script. I copied it into OpenDSS and 'ran' it down to the solve mode=daily, and then I exported the monitor that you have monitoring only G2.

The maximum active generation I see it producing is 40,002.6 kW, at one of the time steps. I also believe it is also obeying your reactive power limits as well.

I don't know if you intentionally left it out, but I noticed you are missing an export monitors m26.

I also plotted the output using your next to the last statement (channels 1,3 and 5 which are active power by phase), and I see each phase of the generator is exporting about 13,333 kW, which adds up to about 40,000 kW.

Here is the direct output of the monitor and an excel file where I summed up the three phases for P and Q for the monitor from the export file I generated from the monitor. The values are in kW and kvar.

I hope this is helpful.

[Cheng Wei Lin]

Dear wsunderm1
In fact I have solved the problem three hours ago,but I have to thanks for your answer.
But there is a small question,in my original opinion,I think it should produce total kW=40000 for 24hr,but now it has produce 40000kW in just one hour,so the total kW that G2 produce is 40000x12=480000kW.I think G2 produce too much kW.Can I have some method to produce 40000kW in 24hr?

Another question is I want to monitor the voltage for each bus.So I specified the command "new monitor.m27 element=bus.B1 terminal=1 mode=0".But the opendss can not run,it said that it can't find the bus.b1.So how can I detect the pu value for each bus?

Thanks

C.W Lin

[wsunderm1]

C. W Lin,

I think there is some difference in terminology. The generator G2 is set to produce 40,000 kW instantaneously. This is active power.

What I understand you asking for in your latest message (if I understand correctly), is you want to generate 40,000 kW-hours over 24 hours, which is energy.

In that case, just divide 40,000 kW by 24, and that is the per-hour output you want to set the generator to produce (you may want to proportionally reduce your reactive power, but that is up to you), if you want to generate approximately 40,000 kW-hours in 24 hours.

To your second question,a bus is not a power delivery or power conversion bus. When it is asking for an element, it is asking for either a power delivery or power conversion bus.

To monitor the voltage at bus1 you would set element=line.1-2 and terminal=1 (this is the end of the line that has bus1 on it).

To see a complete list of which elements are power delivery and which elements are power conversion, see the help. It lists all of these and groups them by type (in the first window on the left: PD=power delivery, PC=power conversion).

Thanks,

Wes

[Cheng Wei Lin]

Dear wsunderm1
Thanks a lot .I think it runs successfully now.But I'm not sure the result is correct.I want to check the impact of the relationship between the distance and the change of the voltage .For eaxample ,I set the monitor m27 that close to the PV system and set the monitor m28 that far from the PV system.But the voltage at m28 is much higher than the m27.In fact ,I know that the voltage close to the substation must higher than the end of the feeder.But I'm not sure the result I get is correct.Can you check it for me?And as follow is my new version for my script.

Thanks

C.W Lin

[Bandrea87]

Hi Wei

The IEEE 30 network is a little tricky case because it may requires more iterations to converge.

The script you just attached indeed does not converge (you can see it in Show-summary) and consequently the problem you are seeing are perhaps correlated with that. This happens because in the script you just posted the command Set maxiterations=100 is not there anymore (whilst it was in your first script in the first post).

Regards
Andrea

[Cheng Wei Lin]

Hi Andrea
Thanks a lot,I returned the command "maxiterations=100" to my script and it seem run successfully now.

Then I noticed that when I do "show summery" and the information there show the "max pu" and "min pu".Is there have any method to find where the nodes are.In fact,there is provision that after connecting the PV system,the rate of change for voltage can't over the 2.5%.So I want to know the pu value for each bus.No matter I connected the PV system or not,the max pu and min pu is the same.It is too strange.By the way,I am doing the time series,so I think "show" and "export" may not be suitable me because it always display the result for the last iteration point.

Thanks

C.W Lin

[Roger Dugan]

I would simply do "Show Voltage LN Node" and sort on the per unit column.

If your Editor won't sort for you, copy-and-paste into Excel and sort on the per unit column.

[Cheng Wei Lin]

Hi Roger
But when I do "Show Voltage LN Node" and I remember it just show the last point.For example,in my script,it show last point,the 24h pu value.If I want to see the pu value for each hour,for example,1h,2h,3h........,22h,23h,24h and what can I do?

[Bandrea87]

Hi Cheng

You might still use OpenDSS alone and producing a "Show Voltage LN Node" for every instant by defining a solution step of 1 (i.e., set Number=1) and repeating the command Solve - Show Voltage LN Node for any instant (e.g., 24 for a daily with 1 hour resolution).

' In the solution statement
set mode=daily
set number=1
set stepsize=1h

solve
Show Voltage LN Node

solve
Show Voltage LN Node

solve
Show Voltage LN Node

solve
Show Voltage LN Node

...

However, this solution might not be useful for you if you are going further. In this case you might need to pass at the COM interface. In this case you might consider the following script (MATLAB but you can do in Python, VBA... see Primer Manual on how to do it) for every instant:

....
DSSSolution.Number=1
....

for iInstant=1:24
DSSSolution.Solve()
% Collect the results in the A matrix
A(iInstant,:)=DSSCircuit.AllBusVmag
% Do something with A (find max-min..)
end

Regards
Andrea

[Paulo Radatz]

Hello Cheng,
I think I can help you. I know two ways to get the Voltage LN node for each iteration.

The first one is ok if you are not using another software to control OpenDSS via COM interface.
For example, you want to show Voltage LN node for hour 5.

set mode=daily stepsize=1h number=24
solve number = 5 ! Change this hour to get the voltages for that time
show Voltage LN

So, you can change the number to any hour you want to get the voltages.

The second one, you need to use (Python, Matlab, etc..) to control OpenDSS. In my case, I use Python.
Attached are two files:

1- Cheng_test.py: The Python file controls OpenDSS to do what you want automatically (show voltage ln node). You can open it using notepad++
2- TCC-Rede13-daily.zip: The folder has all the .dss files that describe the example I ran. (Realize that there is not any solve command in any .dss file. The solve command is in Cheng_test.py file.)

If you don't use any software to control OpenDSS, you can check the Python docs out.
https://sourceforge.net/p/electricdss/code/HEAD/tree/trunk/Distrib/Examples/Python/

I hope it helps.

Paulo Radatz

[Roger Dugan]

Several good suggestions. Let expand a little on Andrea's posting...

If you change the CaseName for each solution, you will get a different output file name, so you will have the solution at each time step.

! In the solution statement
set mode=daily
set number=1
set stepsize=1h

Set Casename=Hour01
solve
Show Voltage LN Node

Set Casename=Hour02
solve
Show Voltage LN Node

Set Casename=Hour03
solve
Show Voltage LN Node

Set Casename=Hour04
solve
Show Voltage LN Node

[Cheng Wei Lin]

Thanks for Andrea,Paulo and Roger 's answer
I think i still have a lot to learn.
Then I have a question,when I canceled all the shunt capacitors , generators and the PV system but why the pu value at the feeder(for example,at hour02 the B12) is higher than the B1(close to substation)?It is unreasonable because I think B1 should have the highest pu value.

[Bandrea87]

Hi Cheng

I am in a bus with a small phone   and slow connection hence I cannot test your circuit. However I can tell you what I would do at your place.

First I would run a power flow in only one instant as Paulo explained (at the 2nd hour for instance). Then I would plot the voltage profile of the whole network to see clearly what's going on. To do that you need to install an energy meter at the head of the network (line
1-2) and then go in plot-profile in the OpenDss menu (if I well remember).

Once you have the plot ask yourself if they are strange. If so look also at the absplute value. OpenDss changes load model in attempt to converge when the voltage drop below 0.95pu. However, if the circuit is well designed this should not happen.

Do the voltage have weird values? If so the problem is with the network. The network is given by topology+ loads. The topology, if you did not change it from the original, should be ok. Hence perhaps is the load too big? In this case I woulf try to reduce the load to see if the situation improves (you can use batch edit property to do it for all of them with on line of code). If you did change also the topology I would go for dwbugging function in OpenDsa as show loop-kvbasemismatch (in the show menu of Opendss).

Hope this helps.

Regards
Andrea

[Rakib Hossain]

How can I add a PV in a bus using python command?

[Davis Montenegro]

Hello,

By using the text interface, send the PV definition the same way you write it on a script:

DSSText.Command='New PVSystem.myPV phases=3 bus1=A Pmpp=1000 kV=13.8 kVA...'

Best regards

Davis

[Anshuman Srivastava]

Hey all,
I've been working on integrating a PV system into the IEEE 13-bus test feeder using OpenDSS. I've successfully set up the simulation and plotted the per-unit voltage (Vpu) for a 12-hour period based on the irradiance values from a CSV file. The Vpu plot shows how the voltage varies over time.
Now, I want to extract the Vpu at a specific time of the day (e.g., 6 hours into the simulation). I understand that OpenDSS can run simulations in different modes, and I believe I should be able to get the Vpu at a particular time as per the generated plot. But I'm not sure how to set it up to get the Vpu for a particular timestamp.

I tried adding these lines but it did not help. The Vpu at the bus remained constant for all the sec/hour.
//set mode=daily stepsize=1s number=43300
//solve number = 20000 ! Change this hour to get the voltages for that time
//show Voltage LN

I have attached my current script below.
Kindly

Clear

! Initial Setup for IEEE 13-Node Test Feeder
new circuit.IEEE13Nodeckt
~ basekv=115 pu=1.0001 phases=3 bus1=SourceBus
~ Angle=30
~ MVAsc3=20000 MVASC1=21000

! Sub Transformer Definition
New Transformer.Sub Phases=3 Windings=2 XHL=(8 1000 /)
~ wdg=1 bus=SourceBus conn=delta kv=115 kva=5000 %r=(.5 1000 /) XHT=4
~ wdg=2 bus=650 conn=wye kv=4.16 kva=5000 %r=(.5 1000 /) XLT=4

! Feeder 1-Phase Voltage Regulators
New Transformer.Reg1 phases=1 XHL=0.01 kVAs=[1666 1666]
~ Buses=[650.1 RG60.1] kVs=[2.4 2.4] %LoadLoss=0.01
new regcontrol.Reg1 transformer=Reg1 winding=2 vreg=122 band=2 ptratio=20 ctprim=700 R=3 X=9

New Transformer.Reg2 phases=1 XHL=0.01 kVAs=[1666 1666]
~ Buses=[650.2 RG60.2] kVs=[2.4 2.4] %LoadLoss=0.01
new regcontrol.Reg2 transformer=Reg2 winding=2 vreg=122 band=2 ptratio=20 ctprim=700 R=3 X=9

New Transformer.Reg3 phases=1 XHL=0.01 kVAs=[1666 1666]
~ Buses=[650.3 RG60.3] kVs=[2.4 2.4] %LoadLoss=0.01
new regcontrol.Reg3 transformer=Reg3 winding=2 vreg=122 band=2 ptratio=20 ctprim=700 R=3 X=9

! Transformer Definition
New Transformer.XFM1 Phases=3 Windings=2 XHL=2
~ wdg=1 bus=633 conn=Wye kv=4.16 kva=500 %r=.55 XHT=1
~ wdg=2 bus=634 conn=Wye kv=0.480 kva=500 %r=.55 XLT=1

! Line Codes
redirect IEEELineCodes.dss

New linecode.mtx601 nphases=3 BaseFreq=60
~ rmatrix = (0.3465 | 0.1560 0.3375 | 0.1580 0.1535 0.3414 )
~ xmatrix = (1.0179 | 0.5017 1.0478 | 0.4236 0.3849 1.0348 )
~ units=mi
New linecode.mtx602 nphases=3 BaseFreq=60
~ rmatrix = (0.7526 | 0.1580 0.7475 | 0.1560 0.1535 0.7436 )
~ xmatrix = (1.1814 | 0.4236 1.1983 | 0.5017 0.3849 1.2112 )
~ units=mi
New linecode.mtx603 nphases=2 BaseFreq=60
~ rmatrix = (1.3238 | 0.2066 1.3294 )
~ xmatrix = (1.3569 | 0.4591 1.3471 )
~ units=mi
New linecode.mtx604 nphases=2 BaseFreq=60
~ rmatrix = (1.3238 | 0.2066 1.3294 )
~ xmatrix = (1.3569 | 0.4591 1.3471 )
~ units=mi
New linecode.mtx605 nphases=1 BaseFreq=60
~ rmatrix = (1.3292 )
~ xmatrix = (1.3475 )
~ units=mi

New Linecode.mtx606 nphases=3 Units=mi
~ Rmatrix=[0.791721 |0.318476 0.781649 |0.28345 0.318476 0.791721 ]
~ Xmatrix=[0.438352 |0.0276838 0.396697 |-0.0184204 0.0276838 0.438352 ]
~ Cmatrix=[383.948 |0 383.948 |0 0 383.948 ]
New linecode.mtx607 nphases=1 BaseFreq=60
~ rmatrix = (1.3425 )
~ xmatrix = (0.5124 )
~ cmatrix = [236]
~ units=mi

! Load Definitions
New Load.671 Bus1=671.1.2.3 Phases=3 Conn=Delta Model=1 kV=4.16 kW=1155 kvar=660
New Load.634a Bus1=634.1 Phases=1 Conn=Wye Model=1 kV=0.277 kW=160 kvar=110
New Load.634b Bus1=634.2 Phases=1 Conn=Wye Model=1 kV=0.277 kW=120 kvar=90
New Load.634c Bus1=634.3 Phases=1 Conn=Wye Model=1 kV=0.277 kW=120 kvar=90
New Load.645 Bus1=645.2 Phases=1 Conn=Wye Model=1 kV=2.4 kW=170 kvar=125
New Load.646 Bus1=646.2.3 Phases=1 Conn=Delta Model=2 kV=4.16 kW=230 kvar=132
New Load.692 Bus1=692.3.1 Phases=1 Conn=Delta Model=5 kV=4.16 kW=170 kvar=151
New Load.675a Bus1=675.1 Phases=1 Conn=Wye Model=1 kV=2.4 kW=485 kvar=190
New Load.675b Bus1=675.2 Phases=1 Conn=Wye Model=1 kV=2.4 kW=68 kvar=60
New Load.675c Bus1=675.3 Phases=1 Conn=Wye Model=1 kV=2.4 kW=290 kvar=212
New Load.611 Bus1=611.3 Phases=1 Conn=Wye Model=5 kV=2.4 kW=170 kvar=80
New Load.652 Bus1=652.1 Phases=1 Conn=Wye Model=2 kV=2.4 kW=128 kvar=86
New Load.670a Bus1=670.1 Phases=1 Conn=Wye Model=1 kV=2.4 kW=17 kvar=10
New Load.670b Bus1=670.2 Phases=1 Conn=Wye Model=1 kV=2.4 kW=66 kvar=38
New Load.670c Bus1=670.3 Phases=1 Conn=Wye Model=1 kV=2.4 kW=117 kvar=68

! Capacitor Definitions
New Capacitor.Cap1 Bus1=675 phases=3 kVAR=600 kV=4.16
New Capacitor.Cap2 Bus1=611.3 phases=1 kVAR=100 kV=2.4

! Line Definitions
New Line.650632 Phases=3 Bus1=RG60.1.2.3 Bus2=632.1.2.3 LineCode=mtx601 Length=2000 units=ft
New Line.632670 Phases=3 Bus1=632.1.2.3 Bus2=670.1.2.3 LineCode=mtx601 Length=667 units=ft
New Line.670671 Phases=3 Bus1=670.1.2.3 Bus2=671.1.2.3 LineCode=mtx601 Length=1333 units=ft
New Line.671680 Phases=3 Bus1=671.1.2.3 Bus2=680.1.2.3 LineCode=mtx601 Length=1000 units=ft
New Line.632633 Phases=3 Bus1=632.1.2.3 Bus2=633.1.2.3 LineCode=mtx602 Length=500 units=ft
New Line.632645 Phases=2 Bus1=632.3.2 Bus2=645.3.2 LineCode=mtx603 Length=500 units=ft
New Line.645646 Phases=2 Bus1=645.3.2 Bus2=646.3.2 LineCode=mtx603 Length=300 units=ft
New Line.692675 Phases=3 Bus1=692.1.2.3 Bus2=675.1.2.3 LineCode=mtx606 Length=500 units=ft
New Line.671684 Phases=2 Bus1=671.1.3 Bus2=684.1.3 LineCode=mtx604 Length=300 units=ft
New Line.684611 Phases=1 Bus1=684.3 Bus2=611.3 LineCode=mtx605 Length=300 units=ft
New Line.684652 Phases=1 Bus1=684.1 Bus2=652.1 LineCode=mtx607 Length=800 units=ft

! Switch Definitions
New Line.671692 Phases=3 Bus1=671 Bus2=692 Switch=y r1=1e-4 r0=1e-4 x1=0.000 x0=0.000 c1=0.000 c0=0.000

Set Voltagebases=[115, 4.16, .48]
calcv
Solve
BusCoords IEEE13Node_BusXY.csv

! Show some results (commented out)
!Show Voltages LN Nodes
! Show Currents Elem
! Show Powers kVA Elem
! Show Losses
! Show Taps

//Alternate Solution Script
//To force the taps to be same as published results, set the transformer taps manually and disable the controls
//Transformer.Reg1.Taps=[1.0 1.0625]
//Transformer.Reg2.Taps=[1.0 1.0500]
//Transformer.Reg3.Taps=[1.0 1.06875]
//Set Controlmode=OFF
//Solve

! Define the LoadShape for the PV system
New LoadShape.PV npts=100000 sinterval=1.0 mult=(File=1-day-1-s-Solar-2.csv) UseActual=Yes

! Define the PV Generator
New Generator.Solar Bus=611.3 kV=2.4 kW=2000 PF=1 Duty=PV Phases=1

! Define and set up the monitor for the duty cycle simulation
New Monitor.Mon611 Element=Line.684611 Terminal=1
Solve mode=dutycycle number=43300 stepsize=1s

//set mode=daily stepsize=1s number=43300
//solve number = 20000 ! Change this hour to get the voltages for that time
//show Voltage LN

! Plot the monitor data
Plot Monitor object=Mon611 Channels=(1) Bases=[2400,2400,2400]

Please help!

Best Regards,
Anshuman

[Paulo Radatz]

Hi Anshuman,

Set mode=Daily
Set stepsize=1s
Set number=1
set time=(10, 1) ! OpenDSS simulates hour 10 second 2. (time + stepsize)

solve

If you want to watch a video about daily mode, check this out: https://www.youtube.com/watch?v=x5-f6e0UreI&list=PLhdRxvt3nJ8yZH-xVuE-4mifMHRlJK2n3&index=12

Best regards,

Paulo Radatz

[Anshuman Srivastava]

Thank you so much.
I will look into it.

[malaka hanhan]

Hi everyone
I'm new here and new in opendss ; I run the same code above and then add the pv system but that dosn't differ in the results ; and the output power of the pv is zero , what could be the problem?

also I wonder how can I get a power reduction!! .. the load power results always remains the same at all load flow scenarios even though I change the load model to 4 for example.
please help!!
best regards
malaka

[wsunderm1]

The output of the PV is not zero in a snapshot solution. However, it is zero at the beginning and the end of the daily solution, due to the irradiance curve that is defined.

There was at least 1 issue in the code where there are duplicate set voltagebases and calcv. Only 1 set is needed. However, I don't think this causes problems with the solution.

I noticed that the kVA of the PVSystem was defined to be 50, while the Pmpp is equal to 500. Was the PVSystem kVA also supposed to be equal to 500?

In any event, when running the PVSystem with kVA equal to 50 , and placing a monitor directly on the PVSystem (named PV) and plot the output, it looks like the attached.

Each phase is outputting approximately -12 kW at peak output due to irrad=0.8 times irradiance curve maximum value of 1.0 times 50 kVA divided by three (divide by three because we're plotting each phase separately).