Hello Dear,
I am new to OpenDss. I am trying to add PV (on the house rooftop) on the power distribution system. How can I sure that my simuation is correct? I am using the ETAP to validate my simuation, but it is not matching with my my OpenDss simuation results. Thank you
here is code:
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
Hello Dear,
I am new to OpenDss. I am trying to add PV (on the house rooftop) on the power distribution system. How can I sure that my simuation is correct? I am using the ETAP to validate my simuation, but it is not matching with my my OpenDss simuation results. Thank you
here is code:
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
Show Powers kVA Elements