Hello Dear,
I am new to OpenDss. I have a question. I am trying to add PV (on the house rooftop) on the power distribution system. How can I simulate a 24 hour with load? I saw that it was giving me a snap short after simulation, but I want it over the time if load change. What COM interface use for?
Here is the details about the connection.
The 3-phase transformer is connected to the substation transformer to step down from 66 KV to 12.47 KV voltage. Than, A single-phase transformer is connected from the 3-phase transformer to step down from 12.47 to 240/120 V. The output of single-phase transformer goes into the house panel which is connect PV invertor wire. Thank you.
Here is the code that I am simulating on OpenDss:
Clear
New Circuit.AEC
// ................................Creates voltage source (Vsource.Source)..............................
Edit Vsource.Source BasekV=66 pu=1.05 ISC3=3000 ISC1=2500
//Define base v,ISC3(3 phase source impedance), and 1 phase Z
! Define Transformers for each location (T1.substation)
New Transformer.T1 Buses=[SourceBus.1.2.3, B1.1.2.3] Conns=[Delta Wye] kVs= [66 12.47]
~ kVAs=[20000 20000] XHL=10
// Define TR1 propertics xhl means Percent reactance high to low (winding 1 to winding 2).
! Define Transformers for each location (T2.T1(split_phase)) A is primary and B is secondary
New Transformer.T2 phases=1 Windings=3
~ Xhl=2.04 Xht=2.04 Xlt=1.36 %noloadloss=.2
~ Buses=[A2.1.2 B2.1.0 B2.0.2] !!! Note polarity
~ kVs=[12.47 .12 .12] ! ratings of windings
~ kVAs=[25 25 25]
~ %Rs = [0.6 1.2 1.2]
~ conns=[delta wye wye]
//Note that the winding voltage
//............................................WIRES....................................................//
New Wiredata.Dove_T1 Diam=0.184 Rdc=0.0000306 Rac=0.0000375
~ NormAmps=726
~ RUnits=ft !Length units for resistance: ohms per {mi|kft|km|m|Ft|in|cm } Default=none.
~ RadUnits=in !Units for outside radius: {mi|kft|km|m|Ft|in|cm } Default=none.
~ GMRUnits=ft !Units for GMR: {mi|kft|km|m|Ft|in|cm } Default=none.
New Wiredata.Dove_T2 Diam=0.708 Rdc=0.0000306 Rac=0.0000375
~ NormAmps=726
~ RUnits=ft !Length units for resistance: ohms per {mi|kft|km|m|Ft|in|cm } Default=none.
~ RadUnits=in !Units for outside radius: {mi|kft|km|m|Ft|in|cm } Default=none.
~ GMRUnits=ft !Units for GMR: {mi|kft|km|m|Ft|in|cm } Default=none.
New PVSystem.PV1 phases=2 bus= PvBus kV=0.240 kVA=50 irrad=0.8 Pmpp=50
~ temperature=25 PF=1 effcurve=My_Eff P-TCurve=MyPv_T
~ Daily=My_Irrad TDaily=My_Temp
~ %cutin=0.1 %cutout=0.1
// Define Pv's parameter (% cut in power -- % of kVA rating of inverter), (% cut out power -- % of kVA rating of inverter)
Set VoltageBases=[66, 12.47, 0.12]
CalcVoltageBases
Study the example in Example Circuit 1 in OpenDSSManual.pdf in the Doc folder. It will show you the basics of defining a Loadshape object and assigning it to a load. In your case you will define irradiance shapes (as Loadshape objects) and temperature shapes as well as efficiency curves.
Yes. You defined the circuit the same way as the MV system (primary distribution). In fact you can put them both in the same model. See examples such as the IEEE 8500 node or any of the EPRI test circuits.
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Hello Dear,
I am new to OpenDss. I have a question. I am trying to add PV (on the house rooftop) on the power distribution system. How can I simulate a 24 hour with load? I saw that it was giving me a snap short after simulation, but I want it over the time if load change. What COM interface use for?
Here is the details about the connection.
The 3-phase transformer is connected to the substation transformer to step down from 66 KV to 12.47 KV voltage. Than, A single-phase transformer is connected from the 3-phase transformer to step down from 12.47 to 240/120 V. The output of single-phase transformer goes into the house panel which is connect PV invertor wire. Thank you.
Here is the code that I am simulating on OpenDss:
Clear
New Circuit.AEC
// ................................Creates voltage source (Vsource.Source)..............................
Edit Vsource.Source BasekV=66 pu=1.05 ISC3=3000 ISC1=2500
//Define base v,ISC3(3 phase source impedance), and 1 phase Z
//.....................................Creates Transformer..............................................................
! Define Transformers for each location (T1.substation)
New Transformer.T1 Buses=[SourceBus.1.2.3, B1.1.2.3] Conns=[Delta Wye] kVs= [66 12.47]
~ kVAs=[20000 20000] XHL=10
// Define TR1 propertics xhl means Percent reactance high to low (winding 1 to winding 2).
! Define Transformers for each location (T2.T1(split_phase)) A is primary and B is secondary
New Transformer.T2 phases=1 Windings=3
~ Xhl=2.04 Xht=2.04 Xlt=1.36 %noloadloss=.2
~ Buses=[A2.1.2 B2.1.0 B2.0.2] !!! Note polarity
~ kVs=[12.47 .12 .12] ! ratings of windings
~ kVAs=[25 25 25]
~ %Rs = [0.6 1.2 1.2]
~ conns=[delta wye wye]
//Note that the winding voltage
//............................................WIRES....................................................//
New Wiredata.Dove_T1 Diam=0.184 Rdc=0.0000306 Rac=0.0000375
~ NormAmps=726
~ RUnits=ft !Length units for resistance: ohms per {mi|kft|km|m|Ft|in|cm } Default=none.
~ RadUnits=in !Units for outside radius: {mi|kft|km|m|Ft|in|cm } Default=none.
~ GMRUnits=ft !Units for GMR: {mi|kft|km|m|Ft|in|cm } Default=none.
New Wiredata.Dove_T2 Diam=0.708 Rdc=0.0000306 Rac=0.0000375
~ NormAmps=726
~ RUnits=ft !Length units for resistance: ohms per {mi|kft|km|m|Ft|in|cm } Default=none.
~ RadUnits=in !Units for outside radius: {mi|kft|km|m|Ft|in|cm } Default=none.
~ GMRUnits=ft !Units for GMR: {mi|kft|km|m|Ft|in|cm } Default=none.
New Wiredata.ACSR1/0 GMR=0.0056 DIAM=0.3980000 RAC=1.120000
~ NormAmps=230.0000
~ Runits=mi radunits=in gmrunits=ft
// ...........................................................creates linecode....................
New Linecode.T1_LC R1=0.058 X1=.1206 R0=.1784 X0=.4047 C1=3.4 C0=1.6 Units=kft // default it has 3 phase
New Linecode.T2_LC R1=0.058 X1=.1206 R0=.1784 X0=.4047 C1=3.4 C0=1.6 Units=kft nphases=2
//........................................................GEOMETRY...................................................//
New LineGeometry.T1_GEO nconds=4 nphases=3
~ cond=1 Wire=Dove_T1 x=-1.3 h=15 units=ft
~ cond=2 Wire=Dove_T1 x=-1.5 h=15 units=ft
~ cond=3 Wire=Dove_T1 x=1.5 h=15 units=ft
~ cond=4 Wire=ACSR1/0 x=0 h=15 units=ft !neutral
New LineGeometry.T2_GEO nconds=2 nphases=1
~ cond=1 Wire=Dove_T2 x=-.1524 h=8.83 units=m
~ cond=2 Wire=ACSR1/0 x=0 h=7.3152 units=m !neutral
// Define Zero Sequence resistance(RO), reactance XO, capacitance CO, and psositive R1,X1,C1
// .............................................connecting line with load and T1,T2................
New Line.Line1 Bus1=B1.1.2.3 Bus2=A2.1.2 Geometry=T1_GEO linecode=T1_LC Length=10 units=ft
New Line.Line2 Bus1=B2.1.0 Bus2=PvBus.1.0 Geometry=T2_GEO linecode=T2_LC Length=10 units=ft
New Line.Line3 Bus1=B2.0.2 Bus2=PvBus.0.2 Geometry=T2_GEO linecode=T2_LC Length=10 units=ft
// ............... ...................Load .........................................................................
New Load.LOAD1 phases=1 Bus1= PvBus.0.2 Model=1 kV=0.120 kVA=5 PF=0.98
New Load.LOAD2 phases=1 Bus1= PvBus.1.0 Model=1 kV=0.120 kVA=5 PF=0.98
// ............................... connecting line with Pv .....................................................
// P-T curve is per unit of rated Pmpp(max power poin in pannel) vs temperature
// This one is for a Pmpp stated at 25 deg
New XYCurve.MyPv_T 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.My_Eff npts=4 xarray=[.1 .2 .4 1.0] yarray=[.86 .9 .93 .97]
// per unit irradiance curve (per unit if "irradiance" property)
New Loadshape.My_Irrad npts=24 interval=1 mult=[0 0 0 0 0 0 .1 .2 .3 .5 .8 .9 1.0 1.0 .99 .9 .7 .4 .1 0 0 0 0 0]
//plot loadshape object=My_Irrad
// 24-hr temp shape curve
New Tshape.My_Temp 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]
New PVSystem.PV1 phases=2 bus= PvBus kV=0.240 kVA=50 irrad=0.8 Pmpp=50
~ temperature=25 PF=1 effcurve=My_Eff P-TCurve=MyPv_T
~ Daily=My_Irrad TDaily=My_Temp
~ %cutin=0.1 %cutout=0.1
// Define Pv's parameter (% cut in power -- % of kVA rating of inverter), (% cut out power -- % of kVA rating of inverter)
Set VoltageBases=[66, 12.47, 0.12]
CalcVoltageBases
// ....................................solve.........................
Solve
// ....................................solve with monitoring .........................
// ............................... display results........................
Show Voltages LN Nodes
Show Currents Element
Study the example in Example Circuit 1 in OpenDSSManual.pdf in the Doc folder. It will show you the basics of defining a Loadshape object and assigning it to a load. In your case you will define irradiance shapes (as Loadshape objects) and temperature shapes as well as efficiency curves.
See the attached example of a PVSystem definition
hlo
can we design secondry distribution system in open dss with specified data
Hello,
Yes you can.
Best regards
Davis
Yes. You defined the circuit the same way as the MV system (primary distribution). In fact you can put them both in the same model. See examples such as the IEEE 8500 node or any of the EPRI test circuits.