import pandas as pd
pd.options.display.float_format = '{:,.5f}'.format
import igraph
import pandapower as pp
import pandapower.plotting
from pandapower.test.shortcircuit.test_1ph import *
LightSimBackend import error: No module named 'grid2op' PhysicalLawChecker import error: No module named 'grid2op' TimeSerie import error: cannot import name 'TimeSerie' from 'lightsim2grid.timeSerie' (C:\ProgramData\miniconda3\envs\py3.11\Lib\site-packages\lightsim2grid\timeSerie.py) ContingencyAnalysis import error: cannot import name 'ContingencyAnalysis' from 'lightsim2grid.contingencyAnalysis' (C:\ProgramData\miniconda3\envs\py3.11\Lib\site-packages\lightsim2grid\contingencyAnalysis.py) rewards import error: No module named 'grid2op'
# beauftragt: Y0y0d5, YN0y0d5, Y0yn0d5, YN0yn0d5, Y0y0y0, Y0d5d5, YN0d5d5, Y0d5y0, Y0y0d11 und D0d0d0
# implementiert zusätzlich: YNdy, Ydyn, YNdyn, YNyy
vg = "yyy"
net = single_3w_trafo_grid(vg)
pp.plotting.simple_plot(net, ext_grid_size=5)
1 in collections is of unknown type. Skipping
<Axes: >
%%capture
net = single_3w_trafo_grid("Yyy")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("Ydd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("Ddd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("Yynd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
1 | 3.33940 | 0.10727 | 18.19496 | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("Ydyn") # Extra vector group
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 8.83645 | 0.00670 | 1.13718 | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNynd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.78326 | 5.12452 | 88.22838 | 15.80517 | 158.05171 |
1 | 3.49934 | 0.21021 | 15.31723 | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNdyn") # Extra vector group
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.79376 | 4.86639 | 85.86380 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 9.04239 | 0.01219 | 0.98955 | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNyd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.78326 | 5.12452 | 88.22838 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNdy") # Extra vector group
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.79376 | 4.86639 | 85.86380 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNdd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.84354 | 3.76841 | 75.01993 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNyy") # Extra vector group
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc) # warum dasselbe Ergebnis wie Yyy?
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
all-bus fault vs single-bus fault vs LU factorization
%%capture
net = single_3w_trafo_grid("Ynynd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.78326 | 5.12452 | 88.22838 | 15.80517 | 158.05171 |
1 | 3.49934 | 0.21021 | 15.31723 | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("Ynynd")
sc.calc_sc(net, fault="1ph", case="max", bus=0)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.78326 | 5.12452 | 88.22838 | 15.80517 | 158.05171 |
%%capture
# now with LU factorization instead of Y inversion
net = single_3w_trafo_grid("Ynynd")
sc.calc_sc(net, fault="1ph", case="max", bus=0, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.78326 | 5.12452 | 88.22838 | 15.80517 | 158.05171 |
%%capture
# now with LU factorization instead of Y inversion
net = single_3w_trafo_grid("Ynynd")
sc.calc_sc(net, fault="1ph", case="max", bus=1, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
1 | 3.49934 | 0.21021 | 15.31723 | 1.52174 | 22.24258 |
%%capture
# now with LU factorization instead of Y inversion
net = single_3w_trafo_grid("Ynynd")
sc.calc_sc(net, fault="1ph", case="max", bus=2, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("Ydd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("Ydd")
sc.calc_sc(net, fault="1ph", case="max", bus=0, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.51934 | 15.80517 | 158.05171 | 15.80517 | 158.05171 |
%%capture
net = single_3w_trafo_grid("Ydd")
sc.calc_sc(net, fault="1ph", case="max", bus=1, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
%%capture
net = single_3w_trafo_grid("Ydd")
sc.calc_sc(net, fault="1ph", case="max", bus=2, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNyd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.78326 | 5.12452 | 88.22838 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNyd")
sc.calc_sc(net, fault="1ph", case="max", bus=0, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.78326 | 5.12452 | 88.22838 | 15.80517 | 158.05171 |
%%capture
net = single_3w_trafo_grid("YNyd")
sc.calc_sc(net, fault="1ph", case="max", bus=1, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
%%capture
net = single_3w_trafo_grid("YNyd")
sc.calc_sc(net, fault="1ph", case="max", bus=2, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNdd")
sc.calc_sc(net, fault="1ph", case="max")
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.84354 | 3.76841 | 75.01993 | 15.80517 | 158.05171 |
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
%%capture
net = single_3w_trafo_grid("YNdd")
sc.calc_sc(net, fault="1ph", case="max", bus=0, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
0 | 1.84354 | 3.76841 | 75.01993 | 15.80517 | 158.05171 |
%%capture
net = single_3w_trafo_grid("YNdd")
sc.calc_sc(net, fault="1ph", case="max", bus=1, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
1 | 0.00000 | inf | inf | 1.52174 | 22.24258 |
%%capture
net = single_3w_trafo_grid("YNdd")
sc.calc_sc(net, fault="1ph", case="max", bus=2, inverse_y=False)
display(net.res_bus_sc)
ikss_ka | rk0_ohm | xk0_ohm | rk_ohm | xk_ohm | |
---|---|---|---|---|---|
2 | 0.00000 | inf | inf | 0.11730 | 2.66335 |
For further vector groups, we included the check in the test, all the vector groups passed.