powerio-dist 0.5.0

Multiconductor distribution network model and lossless converters for OpenDSS, PMD JSON, and BMOPF JSON.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49

"""Solve a .dss case with the OpenDSS engine and print node voltages as JSON.

Usage: <python-with-opendssdirect> solve_dss.py case.dss

Output: {"converged": bool, "voltages": {"<bus>.<node>": [re, im]}, ...} with
voltages in volts. Run it under an interpreter that has opendssdirect
installed (the repo's .venv works). Nothing in the test suite shells out to
this script; it exists for hand checks and for regenerating the engine bus
maps pinned in tests/dss_reader.rs (print AllBusNames + Bus.Nodes per bus
with the same staging as physics_check.py).
"""

import json
import sys


def solve(path):
    import opendssdirect as dss

    dss.Text.Command("Clear")
    dss.Text.Command(f'Redirect "{path}"')
    dss.Text.Command("Solve")

    volts = {}
    for bus in dss.Circuit.AllBusNames():
        dss.Circuit.SetActiveBus(bus)
        nodes = dss.Bus.Nodes()
        raw = dss.Bus.Voltages()  # interleaved re, im per node
        for k, node in enumerate(nodes):
            volts[f"{bus}.{node}"] = [raw[2 * k], raw[2 * k + 1]]

    return {
        "case": path,
        "converged": bool(dss.Solution.Converged()),
        "iterations": dss.Solution.Iterations(),
        "voltages": volts,
    }


def main():
    if len(sys.argv) != 2:
        print(__doc__, file=sys.stderr)
        return 2
    print(json.dumps(solve(sys.argv[1]), indent=1, sort_keys=True))
    return 0


if __name__ == "__main__":
    sys.exit(main())