In [3]:
import numpy as np
from pydgrid.pydgrid import grid

3 bus 4 wire system with transformer

Data

In [13]:
data = {
        "buses":[
                 {"bus": "B1",  "pos_x":   0, "pos_y":   0, "units": "m", "U_kV":20.0},
                 {"bus": "B2",  "pos_x":  10, "pos_y":   0, "units": "m", "U_kV":0.4},
                 {"bus": "B3",  "pos_x": 100, "pos_y":  0, "units": "m", "U_kV":0.4}
                ],
        "grid_formers":[
                        {"bus": "B1",
                        "bus_nodes": [1, 2, 3], "deg": [0, -120, -240],
                        "kV": [11.547, 11.547, 11.547]}
                       ],
        "transformers":[
                        {"bus_j": "B1",  "bus_k": "B2",  "S_n_kVA": 1000.0, "U_j_kV":20, "U_k_kV":0.42,
                         "R_cc_pu": 0.01, "X_cc_pu":0.04, "connection": "Dyn11",   "conductors_j": 3, "conductors_k": 4},
                       ],
        "lines":[
                 {"bus_j": "B2",  "bus_k": "B3",  "code": "lv_cu_150", "m": 100.0},
                ],
        "loads":[{"bus": "B3" , "kVA": 300.0/3, "pf": 0.85,"type":"3P+N"},
                 {"bus": "B3" , "kVA": 300.0/3, "pf": 0.85,"type":"1P+N", 'bus_nodes':[1,4]}
                ],
        "shunts":[
                 {"bus": "B2" , "R": 0.001, "X": 0.0, "bus_nodes": [4,0]}
                 ],
        "line_codes":
            {"lv_cu_150":  {"Rph":0.167,"Xph":0.08, "Rn":0.167, "Xn": 0.08}
            }
       }

Execute power flow

In [14]:
grid_1 = grid()
grid_1.read(data)  # Load data

grid_1.pf()  # solve power flow

Graph with obtained results

In [15]:
from pydgrid.plot_bokeh import plot_results
plot_results(grid_1)
Out[15]:
Figure(
id = '1309', …)
In [ ]:

In [ ]: