-
Notifications
You must be signed in to change notification settings - Fork 2
/
demo.py
66 lines (44 loc) · 1.89 KB
/
demo.py
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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
# test
import sys
sys.path.append('.\src')
from PSATool import PSAT
raw_file = '.\data\IEEE9-2阶.raw'
seq_file = '.\data\IEEE9.seq'
dyr_file = '.\data\IEEE9-2阶.dyr'
simulator = PSAT()
simulator.load_simulation_data(raw_file, 'powerflow')
buses = simulator.get_all_devices('BUS')
print(buses)
loads = simulator.get_all_devices('LOAD')
print(loads)
generators = simulator.get_all_devices('GENERATOR')
print(generators)
simulator.build_network_Y_matrix('basic')
simulator.save_network_Y_matrix('./test_result/IEEE9_Y_mat.csv', 'basic')
simulator.solve_powerflow('NR')
simulator.save_powerflow_result('./test_result/IEEE9_powerflow_result.csv')
simulator.load_simulation_data(dyr_file, 'dynamic')
simulator.build_network_Y_matrix('dynamic')
simulator.save_network_Y_matrix('./test_result/IEEE9_Y_dynamic_mat.csv', 'dynamic')
#for bus in buses:
# simulator.prepare_dynamic_output_meter(bus, 'BUS', 'VM')
#for generator in generators:
# simulator.prepare_dynamic_output_meter(generator, 'GEN', 'ROTOR ANGLE IN DEG')
simulator.set_simulator_parameter('dynamic', 'output_file', './test_result/IEEE9_dynamic_result.csv')
simulator.start_dynamic_simulation()
simulator.run_dynamic_simulation_to_time(0.1)
bus = 7
Yf = -1j*1e6
simulator.set_bus_fault(bus, Yf)
simulator.run_dynamic_simulation_to_time(0.2)
simulator.clear_bus_fault(bus, Yf)
simulator.run_dynamic_simulation_to_time(2.0)
'''
simulator.load_simulation_data(seq_file, 'sequence')
simulator.build_network_Y_matrix('positive')
simulator.save_network_Y_matrix('./test_result/IEEE9_Y_positive_mat.csv', 'positive')
simulator.build_network_Y_matrix('negative')
simulator.build_network_Y_matrix('zero')
simulator.solve_bus_asymmetry_fault(3, 'single phase grounding', 0.0)
simulator.save_fault_analysis_result('./test_result/IEEE9_fault_analysis.csv')
'''