Adding course examples

Author: PauloRadatz

1-buses_min_max_v.py

@@ -0,0 +1,43 @@
+# -*- coding: utf-8 -*-
+# @Author  : Paulo Radatz
+# @Email   : [email protected]
+# @File    : 1-buses_min_max_v.py
+# @Software: PyCharm
+
+import py_dss_interface
+import os
+import pathlib
+import numpy as np
+
+script_path = os.path.dirname(os.path.abspath(__file__))
+
+dss_file = pathlib.Path(script_path).joinpath("../feeders", "123Bus", "IEEE123Master.dss")
+
+dss = py_dss_interface.DSS()  # using OpenDSS provided in the package
+dss.text(f"compile [{dss_file}]")
+dss.text("buscoords buscoords.dat")
+dss.text("solve")
+
+nodes = dss.circuit.nodes_names
+voltages = dss.circuit.buses_vmag_pu
+max_voltage = max(voltages)
+max_voltage_index = voltages.index(max_voltage)
+bus_max_voltages = nodes[max_voltage_index].split(".")[0]
+
+min_voltage = min(voltages)
+min_voltage_index = voltages.index(min_voltage)
+bus_min_voltages = nodes[min_voltage_index].split(".")[0]
+# bus_min = dss.circuit.nodes_names[dss.circuit.buses_vmag_pu.index(min(dss.circuit.buses_vmag_pu))].split(".")[0]
+
+dss.text("set markcapacitor=yes")
+dss.text("set markregulator=yes")
+dss.text(f"AddBusMarker Bus={bus_min_voltages} code=7 color=red size=10")
+dss.text(f"AddBusMarker Bus={bus_max_voltages} code=7 color=green size=10")
+dss.text("plot circuit Power max=2000 n n C1=$00FF0000")
+
+
+voltages_array = np.array(voltages)
+nodes_array = np.array(nodes)
+
+# Bus.Node with voltage less than 0.98 pu
+buses = nodes_array[voltages_array < 0.98]

2-futher_3ph_bus.py

@@ -0,0 +1,34 @@
+# -*- coding: utf-8 -*-
+# @Author  : Paulo Radatz
+# @Email   : [email protected]
+# @File    : 2-futher_3ph_bus.py
+# @Software: PyCharm
+
+import py_dss_interface
+import os
+import pathlib
+
+script_path = os.path.dirname(os.path.abspath(__file__))
+
+dss_file = pathlib.Path(script_path).joinpath("../feeders", "123Bus", "IEEE123Master.dss")
+
+dss = py_dss_interface.DSS()  # using OpenDSS provided in the package
+dss.text(f"compile [{dss_file}]")
+dss.text("New EnergyMeter.Feeder Line.L115 1")
+dss.text("buscoord buscoords.dat")
+dss.text("solve")
+
+bus_dist_actual = 0
+bus_name = None
+for bus in dss.circuit.buses_names:
+    dss.circuit.set_active_bus(bus)
+    if len(dss.bus.nodes) >= 3:
+        bus_dist = dss.bus.distance
+        if bus_dist > bus_dist_actual:
+            bus_dist_actual = bus_dist
+            bus_name = bus
+
+print(f"Bus: {bus_name} with distance: {bus_dist_actual} is the furthest Three-phase bus")
+
+dss.text(f"AddBusMarker Bus={bus_name} code=7 color=Red size=10")
+dss.text("plot circuit Power max=2000 n n C1=$00FF0000")

3-max_load_mult.py

@@ -0,0 +1,34 @@
+# -*- coding: utf-8 -*-
+# @Author  : Paulo Radatz
+# @Email   : [email protected]
+# @File    : 3-max_load_mult.py
+# @Software: PyCharm
+
+import py_dss_interface
+import os
+import pathlib
+
+script_path = os.path.dirname(os.path.abspath(__file__))
+
+dss_file = pathlib.Path(script_path).joinpath("../feeders", "123Bus", "IEEE123Master.dss")
+
+dss = py_dss_interface.DSS()  # using OpenDSS provided in the package
+
+under_voltage = False
+i = 0
+while not under_voltage and i < 100:
+    i = i + 1
+    load_mult = 1 + i / 100
+    dss.text(f"compile [{dss_file}]")
+    dss.text("New EnergyMeter.Feeder Line.L115 1")
+    dss.text(f"set loadmult={load_mult}")
+    dss.text("solve")
+
+    if min(dss.circuit.buses_vmag_pu) < 0.95:
+        under_voltage = True
+        load_mult = 1 + (i - 1) / 100
+
+if under_voltage:
+    print(f"The max loadmult: {load_mult}")
+else:
+    print(f"No problem with max loadmult of: {load_mult}")

4-line_dataframe.py

@@ -0,0 +1,43 @@
+# -*- coding: utf-8 -*-
+# @Author  : Paulo Radatz
+# @Email   : [email protected]
+# @File    : 4-line_dataframe.py
+# @Software: PyCharm
+
+import py_dss_interface
+import os
+import pathlib
+import pandas as pd
+
+script_path = os.path.dirname(os.path.abspath(__file__))
+
+dss_file = pathlib.Path(script_path).joinpath("../feeders", "123Bus", "IEEE123Master.dss")
+
+dss = py_dss_interface.DSS()
+dss.text(f"compile [{dss_file}]")
+
+dss.lines.first()
+line_properties = dss.cktelement.property_names
+
+dict_to_df = dict()
+
+line_name_list = list()
+dss.lines.first()
+for _ in range(dss.lines.count):
+    line_name_list.append(dss.lines.name)
+    dss.lines.next()
+dict_to_df["name"] = line_name_list
+
+for line_property in line_properties:
+    property_list = list()
+
+    dss.lines.first()
+    for _ in range(dss.lines.count):
+        property_list.append(dss.dssproperties.value_read(str(dss.cktelement.property_names.index(line_property) + 1)))
+        dss.lines.next()
+
+    dict_to_df[line_property] = property_list
+
+
+df = pd.DataFrame().from_dict(dict_to_df).set_index("name")
+print("This is cool!")