Menu ▾ ▴

Error Encountered in Solve: Floating point overflow

Beginners
Kim Neumann
2019-10-14
2019-10-15

  • Kim Neumann

    Kim Neumann - 2019-10-14

    Hi
    I'm trying to simulate a distribution system of 10 loads with PV generation, i've set a daily curve for the load and the PV system and I want to see the variation of the results along the day.
    When I try to solve mode = daily or simply solve i get the "Error Encountered in Solve: Floating point overflow" error.
    If I try solve mode = direct or snap it runs, but it only gives me one point of the day.
    Any hint of where may be the floating point??
    Thank you for the help.

    clear

new object=circuit.TCC basekv=138 pu=1.00 MVAsc1=250.25 MVAsc3=511.11

//Metodos de solucao

set trapezoidal=true
set algorithm=Newton
set tolerance=0.00001
set maxiterations=1000

//Subestacao - Barra Infinita (fonte)

edit Vsource.source basekv=138 pu=1.00 MVAsc1=250.25 MVAsc3=511.11 bus1=10

//Tensao base

set voltagebases=[0.38 13.8 138]

calcVoltagebases

//Transformadores subestacao

New Transformer.Trafo1 Buses=[10, 11] Conns=[Delta Wye] kVs= [138 13.8] kVAs= [40000 50000] XHL=14.03

//Parametros dos alimentadores (cabos)

New Linecode.linha R1=0.0498 X1=0.3497 R0=0.258 X0=1.079 Units=km normamps=300

//Alimentadores Media Tensao

New Line.linha1112 Bus1=11 bus2=12 Linecode=linha Length=0.1 Units=km
New Line.linha1213 Bus1=12 bus2=13 Linecode=linha Length=0.2 Units=km
New Line.linha1314 Bus1=13 bus2=14 Linecode=linha Length=0.2 Units=km
New Line.linha1115 Bus1=11 bus2=15 Linecode=linha Length=0.1 Units=km
New Line.linha1116 Bus1=11 bus2=16 Linecode=linha Length=0.5 Units=km

New Line.linha1212c Bus1=12 bus2=12c Linecode=linha Length=0.02 Units=km
New Line.linha1313c Bus1=13 bus2=13c Linecode=linha Length=0.02 Units=km
New Line.linha1414c Bus1=14 bus2=14c Linecode=linha Length=0.02 Units=km

//Transformadores MT para BT

New Transformer.Trafo2 Buses=[16, 30] Conns=[Delta Wye] kVs= [13.8 0.380] kVAs= [0.225 0.250] XHL=7.03
New Transformer.Trafo3 Buses=[15, 20] Conns=[Delta Wye] kVs= [13.8 0.380] kVAs= [0.225 0.250] XHL=7.03

//Alimentadores Baixa Tensao

New Line.linha2021c Phases=3 Bus1=20 bus2=21c Linecode=linha Length=0.03 Units=km
New Line.linha2022 Phases=3 Bus1=20 bus2=22 Linecode=linha Length=0.05 Units=km
New Line.linha2223 Phases=3 Bus1=22 bus2=23 Linecode=linha Length=0.05 Units=km
New Line.linha2324 Phases=3 Bus1=23 bus2=24 Linecode=linha Length=0.05 Units=km

New Line.linha2222c Phases=3 Bus1=22 bus2=22c Linecode=linha Length=0.03 Units=km
New Line.linha2323c Phases=3 Bus1=23 bus2=23c Linecode=linha Length=0.03 Units=km
New Line.linha2424c Phases=3 Bus1=24 bus2=24c Linecode=linha Length=0.03 Units=km

New Line.linha3032c Phases=3 Bus1=30 bus2=32c Linecode=linha Length=0.05 Units=km
New Line.linha3031 Phases=3 Bus1=30 bus2=31 Linecode=linha Length=0.05 Units=km
New Line.linha3033 Phases=3 Bus1=30 bus2=33 Linecode=linha Length=0.05 Units=km

New Line.linha3131c Phases=3 Bus1=31 bus2=31c Linecode=linha Length=0.05 Units=km
New Line.linha3333c Phases=3 Bus1=33 bus2=33c Linecode=linha Length=0.05 Units=km


//Perfis de carga

New LoadShape.industrial1mw npts=24 interval=1 
~ mult=(0.1 0.09 0.08 0.08 0.075 0.077 0.075 0.075 0.087 0.11 0.13 0.12 0.12 0.11 0.125 0.15 0.14 0.125 0.17 0.18 0.24 0.175 0.14 0.125)
New LoadShape.industrial2 npts=24 interval=1 
~ mult=(5 5 5 5 5 5 5 12.5 110 80 75 110 25 115 120 60 75 60 10 9 8 7 6 5)
New LoadShape.comercial npts=24 interval=1 
~ mult=(0.31 0.3 0.29 0.27 0.25 0.28 0.3 0.425 0.9 1.09 1.1 0.95 0.75 0.99 1.1 1.09 1.05 1.15 1.02 1.01 0.9 0.6 0.45 0.38)
New LoadShape.residencial npts=24 interval=1 
~ mult=(0.32 0.25 0.24 0.22 0.24 0.28 0.4 0.44 0.3 0.39 0.38 0.4 0.35 0.4 0.43 0.45 0.41 0.92 1.1 0.82 1.4 0.75 0.55 0.4)

//Cargas 12 ~ 14

New Load.Carga1 Phases=3 Bus1=12c kV=13.8 kW=700 PF=0.98 Daily=industrial1mw
New Load.Carga2 Phases=3 Bus1=13c kV=13.8 kW=500 PF=0.95 Daily=industrial2
New Load.Carga3 Phases=3 Bus1=14c kV=13.8 kW=300 PF=0.96 Daily=comercial

//Cargas 21~24

New Load.Carga4 Phases=3 Bus1=21c kV=0.380 kW=60 PF=0.97 Daily=comercial
New Load.Carga5 Phases=3 Bus1=22c kV=0.380 kW=40 PF=0.95 Daily=comercial
New Load.Carga6 Phases=3 Bus1=23c kV=0.380 kW=30 PF=0.96 Daily=residencial
New Load.Carga7 Phases=3 Bus1=24c kV=0.380 kW=25 PF=0.96 Daily=residencial

//Cargas 31 ~ 33

New Load.Carga8 Phases=3 Bus=31c kV=0.380 kW=45 PF=0.95 Daily=comercial
New Load.Carga9 Phases=3 Bus1=32c kV=0.380 kW=36 PF=0.96 Daily=residencial
New Load.Carga10 Phases=3 Bus1=33c kV=0.380 kW=32 PF=0.96 Daily=residencial

//Sistemas Fotovoltaico

//Curva Pmpp x Temperatura
//Temperatura base de 25 graus, Modulo Jinko 245-265

New XYCurve.curvatemp npts=4 xarray=[0 25 50 75] yarray=[1.1 1.0 0.9 0.78]

//Curva de eficiencia (Potencia x eficiencia)

New XYCurve.efici npts=4 xarray=[.4 .6 .8 1.0] yarray=[.982 .98 .978 .976]

//Curva de irradiacao

New Loadshape.irradi npts=24 interval=1 
~ mult=[0 0 0 0 0 0 .06 .2 .38 .65 .88 .93 .85 .59 .38 .18 .05 0 0 0 0 0 0 0]

//Curva de temperatura

New Tshape.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]

//Gerador sistema PV de cada carga

New PVSystem.PV1 phases=3 bus1=12c kV=13.8 kVA=700 irrad=0.975 Pmpp=700
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp
New PVSystem.PV2 phases=3 bus1=13c kV=13.8 kVA=500 irrad=0.975 Pmpp=500
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp
New PVSystem.PV3 phases=3 bus1=14c kV=13.8 kVA=300 irrad=0.975 Pmpp=300
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp

New PVSystem.PV4 phases=3 bus1=21c kV=0.380 kVA=60 irrad=0.975 Pmpp=60
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp
New PVSystem.PV5 phases=3 bus1=22c kV=0.380 kVA=40 irrad=0.975 Pmpp=40
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp
New PVSystem.PV6 phases=3 bus1=23c kV=0.380 kVA=30 irrad=0.975 Pmpp=30
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp
New PVSystem.PV7 phases=3 bus1=24c kV=0.380 kVA=25 irrad=0.975 Pmpp=25
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp

New PVSystem.PV8 phases=3 bus1=31c kV=0.380 kVA=45 irrad=0.975 Pmpp=45
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp
New PVSystem.PV9 phases=3 bus1=32c kV=0.380 kVA=36 irrad=0.975 Pmpp=32
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp
New PVSystem.PV10 phases=3 bus1=33c kV=0.380 kVA=32 irrad=0.975 Pmpp=32
~ temperature=25 PF=1 effcurve=efici P-TCurve=curvatemp Daily=irradi TDaily=temp

//Monitores

New monitor.linha1112_power element=line.linha1112 terminal=1 mode=1 ppolar=no
New monitor.linha1112_voltage element=line.linha1112 terminal=1 mode=0 ppolar=no
New monitor.carga1_power element=load.carga1 terminal=1 mode=1 ppolar=no
New monitor.carga1_voltage element=load.carga1 terminal=1 mode=0 ppolar=no

set mode = daily
set stepsize = 1h
set number = 24
set hour = 0
set controlmode=time

make bus list
BusCoords coordenadas.csv
//solve mode=direct
//solve mode=snap
solve mode=daily

show powers
//show voltages
//show losses
//export voltages
//export currents
//export overloads

Plot monitor object= linha1112_power
Plot monitor object= carga1_power
     

  • Kim Neumann

    Kim Neumann - 2019-10-15

    Ok so I found one place where it gives me the floating point error:

    New LoadShape.residencial npts=24 interval=1 
~ mult=(0.32 0.25 0.24 0.22 0.24 0.28 0.4 0.44 0.3 0.39 0.38 0.4 0.35 0.4 0.43 0.45 0.41 0.92 1.1 0.82 1.4 0.75 0.55 0.4)

    If I change all values to 1 it gives no floating point error.
    Example:

    New LoadShape.residencial npts=24 interval=1 
~ mult=(1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1)
     

Log in to post a comment.

MongoDB Logo MongoDB