我使用psspy(PSS / E，工程电力系统模拟器)来模拟电力系统。

我在公共汽车上应用了一个三相故障并且排队。现在我想绘制动态模拟过程中的电压大小，我的python代码如下：

import sys

import os

sys_path_psse=r'C:\Program Files (x86)\PTI\PSSEXplore34\PSSPY27'

sys.path.append(sys_path_psse)

os_path_psse=r'C:\Program Files (x86)\PTI\PSSEXplore34\PSSBIN'

os.environ['PATH'] += ';' + os_path_psse

os.environ['PATH'] += ';' + sys_path_psse

import psspy

import redirect

# Redirect output from PSSE to Python:

redirect.psse2py()

# Last case:

CASE = r"C:\Program Files (x86)\PTI\PSSEXplore34\EXAMPLE\savnw.sav"

psspy.psseinit(12000)

psspy.case(CASE)

# Convert loads (3 step process):

psspy.conl(-1,1,1)

psspy.conl(-1,1,2,[0,0],[100,0,0,100])

psspy.conl(-1,1,3)

# Convert generators:

psspy.cong()

# Solve for dynamics

psspy.ordr()

psspy.fact()

psspy.tysl()

# Save converted case

case_root = os.path.splitext(CASE)[0]

psspy.save(case_root + "_C.sav")

# Add dynamics data

psspy.dyre_new(dyrefile="C:\Program Files (x86)\PTI\PSSEXplore34\EXAMPLE\savnw.dyr")

# Add channels by subsystem

# BUS VOLTAGE

psspy.chsb(sid=0,all=1, status=[-1,-1,-1,1,13,0])

# Save snapshot

psspy.snap(sfile="C:\Program Files (x86)\PTI\PSSEXplore34\EXAMPLE\python_test.snp")

# Initialize

psspy.strt(outfile="C:\Program Files (x86)\PTI\PSSEXplore34\EXAMPLE\python_test_1.out")

psspy.run(tpause=0)

# 3-phase fault on bus 201 (default bus fault is a 3phase and there is no bus 200)

psspy.dist_bus_fault(ibus=201)

# Run to 3 cycles

time = 3.0/60.0

psspy.run(tpause=time)

# Clear fault (assuming only part of bus faults)

psspy.dist_clear_fault()

psspy.dist_branch_trip(ibus=201, jbus=151, id='1')

# Run to 20 seconds

time = 1

psspy.run(tpause=time)