European LV distribution grid example #96

doc/elements/bus_res_3ph.csv

@@ -1,10 +1,10 @@
 **Parameter**;**Datatype**;**Explanation**
-vm_A_pu;float;voltage magnitude:Phase A [p.u]
-va_A_degree;float;voltage angle:Phase A [degree]
-vm_B_pu;float;voltage magnitude:Phase B [p.u]
-va_B_degree;float;voltage angle:Phase B [degree]
-vm_C_pu;float;voltage magnitude:Phase C [p.u]
-va_C_degree;float;voltage angle:Phase C [degree]
+vmA_pu;float;voltage magnitude:Phase A [p.u]
+vaA_degree;float;voltage angle:Phase A [degree]
+vmB_pu;float;voltage magnitude:Phase B [p.u]
+vaB_degree;float;voltage angle:Phase B [degree]
+vmC_pu;float;voltage magnitude:Phase C [p.u]
+vaC_degree;float;voltage angle:Phase C [degree]
 p_A_mw;float;resulting active power demand:Phase A [MW]
 q_A_mvar;float;resulting reactive power demand:Phase A [Mvar]
 p_B_mw;float;resulting active power demand:Phase B [MW]

pandapower/networks/__init__.py

@@ -7,4 +7,5 @@
 from pandapower.networks.mv_oberrhein import *
 from pandapower.networks.power_system_test_cases import *
 from pandapower.networks.simple_pandapower_test_networks import *
+from pandapower.networks.ieee_europen_lv_asymmetric import *
 

pandapower/networks/ieee_europen_lv_asymmetric.py

@@ -0,0 +1,52 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Wed Aug 28 12:15:13 2019
+
+@author: uk067483
+"""
+
+import os
+import pandapower as pp
+from pandapower.pf.runpp_3ph import runpp_3ph
+from pandapower import pp_dir
+
+
+def ieee_european_lv_asymmetric(scenario="on_mid"):
+    """
+    Loads the IEEE European LV network, a generic 0.416 kV network serviced by one 0.8 MVA MV/LV transformer
+    station. The network supplies 906 LV buses and 55 1-PH loads
+    The network layout is mostly radial
+
+    The network can be loaded with three different scenarios for On-Peak and Off-Peak load
+    which are defined by scaling factors for loads / generators . 
+    These scenarios are a good starting point for working with this
+    network, but you are of course free to parametrize the network for your use case.
+
+    The network also includes geographical information of lines and buses for plotting.
+
+    OPTIONAL:
+        **scenario** - (str, "on_mid"): defines the scaling for load and generation
+
+                - "on_mid": high load scenario
+                - "off_start": low load scenario at 12:01 AM
+                - "off_end": low load scenario at mignight
+
+    OUTPUT:
+         **net** - pandapower network
+
+    EXAMPLE:
+
+    import pandapower.networks
+    net = pandapower.networks.ieee_european_lv_asymmetric("off_start")
+    """
+    if scenario == "on_mid":
+        net = pp.from_json(os.path.join(pp_dir, "networks", "IEEE_European_LV_On_Peak_mid.json"))
+    elif scenario == "off_start":
+        net = pp.from_json(os.path.join(pp_dir, "networks", "IEEE_European_LV_Off_Peak_start.json"))
+    elif scenario == "off_end":
+        net = pp.from_json(os.path.join(pp_dir, "networks", "IEEE_European_LV_Off_Peak_end.json")) 
+    else:
+        raise ValueError("Unknown scenario %s - chose 'on_mid' or 'off_start' or 'off_start'" % scenario)
+
+    runpp_3ph(net)
+    return net

pandapower/pf/runpp_3ph.py

@@ -137,7 +137,7 @@ def _load_mapping(net, ppci1):
 
     # last return varaible left for constant impedance loads
     return np.vstack([params['S'+phase+'delta']] for phase in phases),\
-        np.vstack([params['S'+phase+'wye']] for phase in phases), {}
+        np.vstack([params['S'+phase+'wye']] for phase in phases)
 
 
 # =============================================================================
@@ -418,7 +418,7 @@ def runpp_3ph(net, calculate_voltage_angles="auto", init="auto",
 #     P Q values aggragated and summed up for each bus to make s_abc matrix
 #     s_abc for wye connections ; s_abc_delta for delta connection
 # =============================================================================
-    s_abc_delta, s_abc, _ = _load_mapping(net, ppci1)
+    s_abc_delta, s_abc = _load_mapping(net, ppci1)
     # =========================================================================
     # Construct Sequence Frame Bus admittance matrices Ybus
     # =========================================================================

tutorials/create_advanced.ipynb

@@ -5177,7 +5177,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.6.5"
+   "version": "3.7.1"
   }
  },
  "nbformat": 4,