def run_savnw_simulation(datapath, outfile1, outfile2, outfile3, prgfile):
import psspy
psspy.psseinit()
savfile = 'Converted_NETS-NYPS 68 Bus System_C.sav'
snpfile = 'NETS-NYPS 68 Bus System.snp'
if datapath:
savfile = os.path.join(datapath, savfile)
snpfile = os.path.join(datapath,
snpfile) #why produce these two kinds of files?
psspy.lines_per_page_one_device(1, 90)
psspy.progress_output(2, prgfile, [0, 0])
ierr = psspy.case(savfile)
if ierr:
psspy.progress_output(1, "", [0, 0])
print(" psspy.case Error")
return
ierr = psspy.rstr(snpfile)
if ierr:
psspy.progress_output(1, "", [0, 0])
print(" psspy.rstr Error")
return
# fault + line trip
psspy.strt(0, outfile1)
psspy.run(0, 1.0, 1000, 1, 0)
psspy.dist_bus_fault(52, 1, 138.0, [0.0, -0.2E+10])
psspy.run(0, 1.1, 1000, 1, 0)
psspy.dist_clear_fault(1)
psspy.dist_branch_trip(52, 55, '1')
psspy.run(0, 1.2, 1000, 1, 0)
psspy.dist_machine_trip(1, '1')
psspy.run(0, 5.0, 1000, 1, 0)
# line trip (with faults) + generator trip
psspy.case(savfile)
psspy.rstr(snpfile)
psspy.strt(0, outfile2)
psspy.run(0, 1.0, 1000, 1, 0)
psspy.dist_bus_fault(52, 1, 138.0, [0.0, -0.2E+10])
psspy.run(0, 1.1, 1000, 1, 0)
psspy.dist_clear_fault(1)
psspy.run(0, 1.2, 1000, 1, 0)
psspy.dist_machine_trip(8, '1')
psspy.run(0, 5.0, 1000, 1, 0)
psspy.case(savfile)
psspy.rstr(snpfile)
psspy.strt(0, outfile3)
psspy.run(0, 1.0, 1000, 1, 0)
psspy.dist_branch_trip(32, 33, '1')
psspy.run(0, 5.0, 1000, 1, 0)
psspy.lines_per_page_one_device(2, 10000000)
psspy.progress_output(1, "", [0, 0])
def run_savnw_simulation(datapath, outfile1, outfile2, outfile3, prgfile):
import psspy
psspy.psseinit()
savfile = 'savcnv.sav'
snpfile = 'savnw.snp'
if datapath:
savfile = os.path.join(datapath, savfile)
snpfile = os.path.join(datapath, snpfile)
psspy.lines_per_page_one_device(1, 90)
psspy.progress_output(2, prgfile, [0, 0])
ierr = psspy.case(savfile)
if ierr:
psspy.progress_output(1, "", [0, 0])
print(" psspy.case Error")
return
ierr = psspy.rstr(snpfile)
if ierr:
psspy.progress_output(1, "", [0, 0])
print(" psspy.rstr Error")
return
psspy.strt(0, outfile1)
psspy.run(0, 1.0, 1000, 1, 0)
psspy.dist_bus_fault(154, 1, 230.0, [0.0, -0.2E+10])
psspy.run(0, 1.05, 1000, 1, 0)
psspy.dist_clear_fault(1)
psspy.run(0, 5.0, 1000, 1, 0)
psspy.case(savfile)
psspy.rstr(snpfile)
psspy.strt(0, outfile2)
psspy.run(0, 1.0, 1000, 1, 0)
psspy.dist_machine_trip(3018, '1')
psspy.run(0, 5.0, 1000, 1, 0)
psspy.case(savfile)
psspy.rstr(snpfile)
psspy.strt(0, outfile3)
psspy.run(0, 1.0, 1000, 1, 0)
psspy.dist_branch_trip(3005, 3007, '1')
psspy.run(0, 5.0, 1000, 1, 0)
psspy.lines_per_page_one_device(2, 10000000)
psspy.progress_output(1, "", [0, 0])
def Islanding(Time_Trip, Time_Reconnect, End_Time, Out_File): psspy.strt(0,Out_File) #Start our case and specify our output file psspy.run(0,0.0,1,1,0) #Run until 0 seconds psspy.run(0, 10,1,1,0) #Run until 50 seconds psspy.dist_branch_trip(323,325,r"""1""") #Open branch between main grid and microgrid psspy.change_channel_out_file(Out_File) #Resume our output file psspy.run(0, 10,1,1,0) #Run until 100 seconds psspy.dist_branch_trip(223,318,r"""1""") #Open branch between main grid and microgrid psspy.change_channel_out_file(Out_File) #Resume our output file psspy.run(0, 10,1,1,0) #Run until 100 seconds psspy.dist_branch_trip(121,325,r"""1""") #Open branch between main grid and microgrid psspy.change_channel_out_file(Out_File) #Resume our output file psspy.run(0, 10,1,1,0) #Run until 100 seconds
# get the nominal voltages as well as the fault impedance in ohms
FaultBusNomVolt = float(BusDataDict[str(FaultBus)].NominalVolt)
Zbase = FaultBusNomVolt**2/Sbase # float since Sbase is a float
Rohm = FaultRpu*Zbase # fault impedance in ohms
# run simulation till just before the fault
output = StringIO.StringIO()
with silence(output):
# load the sav and snp file
psspy.case(savFile)
psspy.rstr(snpFile)
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.strt(0,out_file)
ierr = psspy.run(0,0.1,1,1,1)
ierr = psspy.dist_branch_trip(L1Bus1, L1Bus2, L1cktID)
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0,0.2,1,1,1) #fault on time
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print('For ' + event + ':')
print('Network did not converge between branch 1 trip and fault application, skipping...')
continue
#######
# check for convergence during fault
#output = StringIO.StringIO()
with silence(output):
#
# Save snapshot
psspy.snap(sfile=r'{0}\PythonDynTest.snp'.format(example_path))
# Initialize and run the dynamic scenario
psspy.strt(option=0, outfile=out_file)
psspy.run(0, 1, 0, 0, 0)
# 3-phase fault on bus 151 (default bus fault is a 3phase and there is no bus 151)
psspy.dist_bus_fault(ibus=151)
# Run to 3 cycles
time = 3.0 / 60.0
psspy.run(0, 1 + time, 0, 0, 0)
# Clear fault (assuming only part of bus faults)
psspy.dist_clear_fault()
psspy.dist_branch_trip(ibus=151, jbus=201, id='1')
# Run to 10 seconds
time = 10
psspy.run(0, time, 0, 0, 0)
# Export channel data to Excel
dyntools.CHNF.xlsout(dyntools.CHNF(out_file),
channels='',
show='True',
xlsfile='out.xls',
sheet='',
overwritesheet=True)
def main():
try:
''' Drives a PSS/E Dynamic simulation and returns values '''
##### Get everything set up on the PSSE side
redirect.psse2py()
#output = StringIO.StringIO()
with silence():
psspy.psseinit(buses=80000)
_i = psspy.getdefaultint()
_f = psspy.getdefaultreal()
_s = psspy.getdefaultchar()
"""
# Redirect any psse outputs to psse_log
psspy.report_output(2,psse_log,[0,0])
psspy.progress_output(2,psse_log,[0,0]) #ignored
psspy.alert_output(2,psse_log,[0,0]) #ignored
psspy.prompt_output(2,psse_log,[0,0]) #ignored
"""
k = 1
for rawFile in RawFileList:
# get the percentage loading from the raw file name
if rawFile == 'savnw_conp.raw':
PL = '100'
else:
rawFileName = rawFile.replace('.raw', '')
PL = rawFileName[-3:]
#Parameters. CONFIGURE THIS
settings = {
# use the same raw data in PSS/E and TS3ph #####################################
'filename':
rawFile, #use the same raw data in PSS/E and TS3ph
################################################################################
'dyr_file':
dyrFile,
'out_file':
'output2.out',
'pf_options': [
0, #disable taps
0, #disable area exchange
0, #disable phase-shift
0, #disable dc-tap
0, #disable switched shunts
0, #do not flat start
0, #apply var limits immediately
0, #disable non-div solution
]
}
##### Load Raw Datafile and do power flow
print "\n Reading raw file:", settings['filename']
# " Reading raw file: {0:s}".format('text')
FaultRpu = 1e-06
Sbase = 100.0
#FaultBusNomVolt = float(BusDataDict[FaultBus].NominalVolt)
#Zbase = FaultBusNomVolt**2/Sbase # float since Sbase is a float
#Rohm = FaultRpu*Zbase # fault impedance in ohms
##########################
# run nested loops to see if there are any abnormal low voltages
simCount = 0 # to keep track of how many simulations are already done
croppedHVLineSet = list(HVLineSet)
for line1 in croppedHVLineSet:
for line2 in croppedHVLineSet:
# stability_indicator = 1
# Bus_issues = [] # list of buses where issues (low voltage or high dv_dt) are reported
# the lines cannot be the same
if line1 == line2:
continue
# part to ensure there is no duplication of events
currentSet = line1 + ';' + line2
currentSetReverse = line2 + ';' + line1
# if case causes topology inconsistencies, continue
if currentSet in topology_inconsistent_set or currentSetReverse in topology_inconsistent_set:
continue
line1Elements = line1.split(',')
line2Elements = line2.split(',')
# Line 1 params
L1Bus1 = int(line1Elements[0])
L1Bus2 = int(line1Elements[1])
L1cktID = line1Elements[2].strip("'").strip()
# Line 2 params
L2Bus1 = int(line2Elements[0])
L2Bus2 = int(line2Elements[1])
L2cktID = line2Elements[2].strip("'").strip()
FaultBusList = [L2Bus1,
L2Bus2] # apply faults at both buses
for FaultBus in FaultBusList:
output = StringIO.StringIO()
with silence():
ierr = psspy.read(0, settings['filename'])
#This is for the power flow. I'll use the solved case instead
ierr = psspy.fnsl(settings['pf_options'])
##### Prepare case for dynamic simulation
# Load conversion (multiple-step)
psspy.conl(_i, _i, 1, [0, _i], [_f, _f, _f, _f])
# all constant power load to constant current, constant reactive power load to constant admittance
# standard practice for dynamic simulations, constant MVA load is not acceptable
psspy.conl(1, 1, 2, [_i, _i],
[100.0, 0.0, 0.0, 100.0])
psspy.conl(_i, _i, 3, [_i, _i], [_f, _f, _f, _f])
ierr = psspy.cong(0) #converting generators
ierr = psspy.ordr(
0
) #order the network nodes to maintain sparsity
ierr = psspy.fact(
) #factorise the network admittance matrix
ierr = psspy.tysl(0) #solving the converted case
ierr = psspy.dynamicsmode(0) #enter dynamics mode
print "\n Reading dyr file:", settings['dyr_file']
ierr = psspy.dyre_new([1, 1, 1, 1],
settings['dyr_file'])
ierr = psspy.docu(0, 1, [
0, 3, 1
]) #print the starting point of state variables
# select time step ##############################################################
ierr = psspy.dynamics_solution_params(
[_i, _i, _i, _i, _i, _i, _i, _i], [
_f, _f, 0.00833333333333333, _f, _f, _f,
_f, _f
],
'out_file') # the number here is the time step
################################################################################
##### select channels
ierr = psspy.delete_all_plot_channels(
) # clear channels
# get all the bus voltages, angles and frequencies
for bus in BusDataDict:
bus = int(bus)
ierr = psspy.voltage_and_angle_channel(
[-1, -1, -1, bus])
ierr = psspy.bus_frequency_channel([-1, bus])
eventStr = PL + '/' + line1 + ';' + line2 + '/F' + str(
FaultBus)
print 'Event: {}'.format(eventStr)
# get the nominal voltages as well as the fault impedance in ohms
FaultBusNomVolt = float(
BusDataDict[str(FaultBus)].NominalVolt)
Zbase = FaultBusNomVolt**2 / Sbase # float since Sbase is a float
Rohm = FaultRpu * Zbase # fault impedance in ohms
# run simulation till just before the fault
ResultsDict = {}
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.strt(0, settings['out_file'])
ierr = psspy.run(0, 0.1, 1, 1, 1)
ierr = psspy.dist_branch_trip(
L1Bus1, L1Bus2, L1cktID)
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0, 0.2, 1, 1, 1) #fault on time
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + eventStr + ':'
print 'Network did not converge between branch 1 trip and fault application, skipping...'
continue
#######
# check for convergence during fault
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.dist_bus_fault(
FaultBus, 3, 0.0, [Rohm, 0.0])
ierr = psspy.run(0, 0.3, 1, 1, 1) #fault off time
ierr = psspy.dist_clear_fault(1)
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + eventStr + ':'
print 'Network did not converge during fault, skipping...'
continue
# check for convergence between fault clearance and second branch trip
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0, 0.31, 1, 1, 1) #fault off time
ierr = psspy.dist_branch_trip(
L2Bus1, L2Bus2, L2cktID)
ierr = psspy.run(0, 0.35, 1, 1, 1) #fault off time
# check for non-convergence
#output = StringIO.StringIO()
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + eventStr + ':'
print 'Network did not converge between fault clearance and branch 2 trip, skipping...'
continue
# select run time ##############################################################
output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(
0, 10.0, 1, 1, 1
) #exit time (second argument is the end time)
################################################################################
# check for non-convergence
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + eventStr + ':'
print 'Network did not converge sometime after 2nd branch trip, skipping...'
continue
outputData = dyntools.CHNF(settings['out_file'])
data = outputData.get_data()
channelDict = data[
1] # dictionary where the value is the channel description
valueDict = data[
2] # dictionary where the values are the signal values, keys match that of channelDict
tme = valueDict['time'] # get time
ResultsDict['time'] = tme
for key in channelDict:
if key == 'time':
continue
signalDescr = channelDict[key]
words = signalDescr.split()
signalType = words[0].strip()
bus = words[1].strip()
#print Bus + ' ' + signalType
if bus not in ResultsDict:
ResultsDict[bus] = Results()
if signalType == 'VOLT':
ResultsDict[bus].volt = valueDict[key]
elif signalType == 'ANGL':
ResultsDict[bus].angle = valueDict[key]
elif signalType == 'FREQ':
ResultsDict[bus].freq = valueDict[key]
EventsDict[eventStr] = ResultsDict
simCount += 1
print 'Simulation ' + str(simCount) + ' out of ' + str(
totalSims)
# Uncomment next two lines if you want to see the output
#with open('output'+str(k) + '.txt','w') as f:
# f.write(outputStr)
k += 1
save_obj(EventsDict, 'EventData')
except Exception:
traceback.print_exc(file=logfile)
sys.exit(0)
if Branch_kV[i] == 330 or Branch_kV[i] == 275:
fault_time = 0.250
if Branch_kV[i] == 132:
fault_time = 0.720
if fault_type == 4: # phase to phase to ground #line to line to ground, line to line
psspy.dist_spcb_fault(Branch_Outage_List_Start[i], Branch_Outage_List_End[i], '1', [3, 0, 2, 1, 0, 0],
[0.5, 0.0, 0.000001, 0.0, 0.000001])
fault_name = 'TwoPhaseGround'
if Branch_kV[i] == 330 or Branch_kV[i] == 275:
fault_time = 0.250
if Branch_kV[i] == 132:
fault_time = 0.720
psspy.run(0, 5.0 + fault_time, 1000, 1, 0)
psspy.dist_branch_trip(Branch_Outage_List_Start[i], Branch_Outage_List_End[i], '1')
psspy.run(0, 15.000, 1000, 1, 0)
# start draw curves
# new folder if necessary
GraphPath = FigurePath + ClauseName + '/' + LoadScenario + '/' + fault_name
if not os.path.exists(GraphPath):
os.makedirs(GraphPath)
# read data curves
chnfobj = dyntools.CHNF(OutputFilePath)
short_title, chanid, chandata = chnfobj.get_data()
freq_data = numpy.array(chandata[1])
p_data_101 = numpy.array(chandata[7])
q_data_101 = numpy.array(chandata[8])
p_data_102 = numpy.array(chandata[9])
bus3 = tfElements[2].strip()
cktID = tfElements[3].strip("'").strip()
eventStr = 'R{}/T{}'.format(PL, tf)
#cktID = GenDict[genbus]
output = StringIO.StringIO()
with silence(output):
# load the sav and snp file
psspy.case(savFile)
psspy.rstr(snpFile)
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.strt(0, settings['out_file'])
ierr = psspy.run(0, 0.1, 1, 1, 1)
ierr = psspy.dist_branch_trip(
int(bus1), int(bus2),
cktID) # this line can be used for 2 winding tf as well
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0, totSimTime, 1, 1, 1) #simulation end time
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print('For ' + eventStr + ':')
print('Network did not converge after tf trip, skipping...')
continue
outputData = dyntools.CHNF(settings['out_file'])
data = outputData.get_data()
def run_savnw_simulation(datapath, outfile1, outfile2, outfile3, prgfile):
import psspy
psspy.psseinit()
savfile = 'savcnv.sav'
snpfile = 'savnw.snp'
if datapath:
savfile = os.path.join(datapath, savfile)
snpfile = os.path.join(datapath, snpfile)
psspy.lines_per_page_one_device(1,90)
psspy.progress_output(2,prgfile,[0,0]) # directly output to file
ierr = psspy.case(savfile)
if ierr:
psspy.progress_output(1,"",[0,0])
print(" psspy.case Error")
return
ierr = psspy.rstr(snpfile)
if ierr:
psspy.progress_output(1,"",[0,0])
print(" psspy.rstr Error")
return
# fault + line trip
psspy.strt(0,outfile1)
psspy.chsb(0,1, [-1,-1,-1,1,13,0])
psspy.run(0, 1.0,1000,1,0)# start from 1 second, 1000 steps, and 1 writing for 1 output step
psspy.dist_bus_fault(154,1, 230.0,[0.0,-0.2E+10]) # ibus, units, voltage kv
psspy.run(0, 1.1,1000,1,0)# start from 1.1 second, 1000 steps, and 1 writing for 1 output step
psspy.dist_clear_fault(1)
psspy.dist_branch_trip(3005,3007,'1')
psspy.run(0,1.2,1000,1,0)
psspy.dist_machine_trip(3018,'1')
psspy.run(0, 5.0,1000,1,0)
# line trip (with faults) + generator trip
psspy.case(savfile)
psspy.rstr(snpfile)
psspy.strt(0,outfile2)
psspy.chsb(0,1, [-1,-1,-1,1,13,0])
psspy.run(0, 1.0,1000,1,0)
psspy.dist_bus_fault(3005,1,230.0,[0.0,-0.2E+10])
psspy.run(0,1.1,1000,1,0)
psspy.dist_clear_fault(1)
psspy.run(0,1.2,1000,1,0)
psspy.dist_machine_trip(3018,'1')
psspy.run(0, 5.0,1000,1,0)
psspy.case(savfile)
psspy.rstr(snpfile)
psspy.strt(0,outfile3)
psspy.chsb(0,1, [-1,-1,-1,1,13,0])
psspy.run(0, 1.0,1000,1,0)
psspy.dist_branch_trip(3005,3007,'1')
psspy.run(0, 5.0,1000,1,0)
psspy.lines_per_page_one_device(2,10000000)
psspy.progress_output(1,"",[0,0])
def applyDisturbance(self):
psspy.dist_branch_trip(self._ibus,self._jbus,self._branch_id)
def runPSSESimBatches(simList, dyrFile, objName):
import sys, os
# add psspy to the system path
sys.path.append(r"C:\Program Files (x86)\PTI\PSSE33\PSSBIN")
os.environ['PATH'] = (r"C:\Program Files (x86)\PTI\PSSE33\PSSBIN;" +
os.environ['PATH'])
from contextlib import contextmanager
import StringIO
from getBusDataFn import getBusData
@contextmanager
def silence(file_object=None):
#Discard stdout (i.e. write to null device) or
#optionally write to given file-like object.
if file_object is None:
file_object = open(os.devnull, 'w')
old_stdout = sys.stdout
try:
sys.stdout = file_object
yield
finally:
sys.stdout = old_stdout
if file_object is None:
file_object.close()
# Local imports
import redirect
import psspy
import dyntools
# getting the raw file
##### Get everything set up on the PSSE side
redirect.psse2py()
#output = StringIO.StringIO()
with silence():
psspy.psseinit(buses=80000)
_i = psspy.getdefaultint()
_f = psspy.getdefaultreal()
_s = psspy.getdefaultchar()
# some important parameters
FaultRpu = 1e-06
Sbase = 100.0
EventsDict = {}
for event in simList:
eventWords = event.split('/')
RawFileIndicator = eventWords[0].strip()
linesOutage = eventWords[1].strip()
FaultBus = eventWords[2].strip()[
1:] # exclude the 'F' at the beginning
# get the raw file
if RawFileIndicator == '100':
rawFile = 'savnw_conp.raw'
else:
rawFile = 'savnw_conp{}.raw'.format(RawFileIndicator)
#Parameters. CONFIGURE THIS
settings = {
# use the same raw data in PSS/E and TS3ph #####################################
'filename':
rawFile, #use the same raw data in PSS/E and TS3ph
################################################################################
'dyr_file':
dyrFile,
'out_file':
'output2.out',
'pf_options': [
0, #disable taps
0, #disable area exchange
0, #disable phase-shift
0, #disable dc-tap
0, #disable switched shunts
0, #do not flat start
0, #apply var limits immediately
0, #disable non-div solution
]
}
output = StringIO.StringIO()
with silence(output):
ierr = psspy.read(0, settings['filename'])
#This is for the power flow. I'll use the solved case instead
ierr = psspy.fnsl(settings['pf_options'])
##### Prepare case for dynamic simulation
# Load conversion (multiple-step)
psspy.conl(_i, _i, 1, [0, _i], [_f, _f, _f, _f])
# all constant power load to constant current, constant reactive power load to constant admittance
# standard practice for dynamic simulations, constant MVA load is not acceptable
psspy.conl(1, 1, 2, [_i, _i], [100.0, 0.0, 0.0, 100.0])
psspy.conl(_i, _i, 3, [_i, _i], [_f, _f, _f, _f])
ierr = psspy.cong(0) #converting generators
ierr = psspy.ordr(0) #order the network nodes to maintain sparsity
ierr = psspy.fact() #factorise the network admittance matrix
ierr = psspy.tysl(0) #solving the converted case
ierr = psspy.dynamicsmode(0) #enter dynamics mode
print "\n Reading dyr file:", settings['dyr_file']
ierr = psspy.dyre_new([1, 1, 1, 1], settings['dyr_file'])
ierr = psspy.docu(
0, 1, [0, 3, 1]) #print the starting point of state variables
# select time step ##############################################################
ierr = psspy.dynamics_solution_params(
[_i, _i, _i, _i, _i, _i, _i, _i],
[_f, _f, 0.00833333333333333, _f, _f, _f, _f, _f],
'out_file') # the number here is the time step
################################################################################
##### select channels
ierr = psspy.delete_all_plot_channels() # clear channels
BusDataDict = getBusData(rawFile)
# get all the bus voltages, angles and frequencies
for bus in BusDataDict:
bus = int(bus)
ierr = psspy.voltage_and_angle_channel([-1, -1, -1, bus])
ierr = psspy.bus_frequency_channel([-1, bus])
print 'Event: {}'.format(event)
# get the nominal voltages as well as the fault impedance in ohms
FaultBusNomVolt = float(BusDataDict[str(FaultBus)].NominalVolt)
Zbase = FaultBusNomVolt**2 / Sbase # float since Sbase is a float
Rohm = FaultRpu * Zbase # fault impedance in ohms
# run simulation till just before the fault
ResultsDict = {}
# get the line params
line1Elements = linesOutage.split(';')[0].strip()
line2Elements = linesOutage.split(';')[1].strip()
# Line 1 params
line1 = line1Elements.split(',')
L1Bus1 = int(line1[0].strip())
L1Bus2 = int(line1[1].strip())
L1cktID = line1[2].strip("'").strip()
#print L1Bus1
#print L1Bus2
#print L1cktID
# Line 2 params
line2 = line2Elements.split(',')
L2Bus1 = int(line2[0].strip())
L2Bus2 = int(line2[1].strip())
L2cktID = line2[2].strip("'").strip()
#print L2Bus1
#print L2Bus2
#print L2cktID
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.strt(0, settings['out_file'])
ierr = psspy.run(0, 0.1, 1, 1, 1)
ierr = psspy.dist_branch_trip(L1Bus1, L1Bus2, L1cktID)
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0, 0.2, 1, 1, 1) #fault on time
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + event + ':'
print 'Network did not converge between branch 1 trip and fault application, skipping...'
continue
#######
# check for convergence during fault
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.dist_bus_fault(int(FaultBus), 3, 0.0, [Rohm, 0.0])
ierr = psspy.run(0, 0.3, 1, 1, 1) #fault off time
ierr = psspy.dist_clear_fault(1)
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + event + ':'
print 'Network did not converge during fault, skipping...'
continue
# check for convergence between fault clearance and second branch trip
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0, 0.31, 1, 1, 1) #fault off time
ierr = psspy.dist_branch_trip(L2Bus1, L2Bus2, L2cktID)
ierr = psspy.run(0, 0.35, 1, 1, 1) #fault off time
# check for non-convergence
#output = StringIO.StringIO()
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + event + ':'
print 'Network did not converge between fault clearance and branch 2 trip, skipping...'
continue
# select run time ##############################################################
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0, 10.0, 1, 1,
1) #exit time (second argument is the end time)
################################################################################
# check for non-convergence
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print 'For ' + event + ':'
print 'Network did not converge sometime after 2nd branch trip, skipping...'
continue
# write to output file
#with open('outputTmp.txt','w') as f:
# f.write(outputStr)
outputData = dyntools.CHNF(settings['out_file'])
data = outputData.get_data()
channelDict = data[
1] # dictionary where the value is the channel description
valueDict = data[
2] # dictionary where the values are the signal values, keys match that of channelDict
tme = valueDict['time'] # get time
ResultsDict['time'] = tme
for key in channelDict:
if key == 'time':
continue
signalDescr = channelDict[key]
words = signalDescr.split()
signalType = words[0].strip()
bus = words[1].strip()
#print Bus + ' ' + signalType
if bus not in ResultsDict:
ResultsDict[bus] = Results()
if signalType == 'VOLT':
ResultsDict[bus].volt = valueDict[key]
elif signalType == 'ANGL':
ResultsDict[bus].angle = valueDict[key]
elif signalType == 'FREQ':
ResultsDict[bus].freq = valueDict[key]
EventsDict[event] = ResultsDict
return EventsDict
# Run simulation DELTA = 0.005 psspy.dynamics_solution_param_2(realar3=DELTA) # set sample time, DELTA sec. # Simulation Procedure ##1. System initialization and run simulation to t = 1 s. psspy.run(0, 1, 0, 0, 0) ##2. Apply 3-phase fault at line B6 – B7 at t = 1s. psspy.dist_branch_fault(6, 7, '1', 3, 0.0, [75.625e-5, 0]) ##3. Continue simulation to t = 1.1s. psspy.run(0, 1.1, 0, 0, 0) ##4. Trip line B6 – B7 and clear fault at t = 1.1s. psspy.dist_clear_fault() psspy.dist_branch_trip(6, 7, '1') ##5. Continue simulation to t = 6.1s. psspy.run(0, 6.1, 0, 0, 0) ##6. Reschedule generator with run-up M1 from 250 MW to 500 MW and run-down M2 from about 573 MW to 300 MW, and trip load at B7 at t = 6.1s at t = 6.1s psspy.increment_gref(30, '1', 0.25) psspy.increment_gref(31, '1', -0.273) psspy.dist_branch_trip(7, 71, '1') ##7. Continue simulation to t = 30s. psspy.run(0, 30, 0, 0, 0) ### Reconnect Load ##psspy.busdata2(jbus,INTGAR1=1); ##psspy.DIST_BRANCH_CLOSE(ibus, jbus, id); # Data Retrieval achnf = dyntools.CHNF(out_file)
bus2 = lineElements[1].strip()
cktID = lineElements[2].strip()
eventStr = 'R{}/L{}'.format(PL,line)
#cktID = GenDict[genbus]
output = StringIO.StringIO()
with silence(output):
# load the sav and snp file
psspy.case(savFile)
psspy.rstr(snpFile)
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.strt(0,settings['out_file'])
ierr = psspy.run(0,0.1,1,1,1)
ierr = psspy.dist_branch_trip(int(bus1), int(bus2),cktID)
#output = StringIO.StringIO()
with silence(output):
ierr = psspy.run(0,totSimTime,1,1,1) #simulation end time
outputStr = output.getvalue()
if "Network not converged" in outputStr:
print('For ' + eventStr + ':')
print('Network did not converge after line trip, skipping...')
continue
outputData = dyntools.CHNF(settings['out_file'])
data = outputData.get_data()
psspy.conl(1, 0, 2, [0, 0], [100.0, 0.0, 0.0, 100.0])
psspy.conl(0, 1, 2, [0, 0], [100.0, 0.0, -306.02, 303.0])
psspy.conl(0, 1, 3, [0, 0], [100.0, 0.0, -306.02, 303.0])
psspy.ordr(0)
psspy.fact()
psspy.tysl(0)
psspy.bsys(0, 0, [0.4, 500.], 0, [], 0, [], 0, [], 0, [])
# start simulation
psspy.strt_2([0,0], OutputFilePath)
psspy.run(0, 1, 1000, 1, 0)
psspy.run(0, 5, 1000, 1, 0)
psspy.dist_branch_fault(46523, 46660, '1', 1, 275.0, [0.00, -0.2E+10])
fault_name = 'ThreePhase'
psspy.run(0, 5.12, 1000, 1, 0)
psspy.dist_branch_trip(46523, 46660, '1')
psspy.run(0, 20.12, 1000, 1, 0)
psspy.dist_branch_close(46523,46660,r"""1""")
psspy.run(0, 20.22, 1000, 1, 0)
psspy.dist_branch_trip(46523, 46660, '1')
psspy.run(0, 30, 1000, 1, 0)
# start draw curves
# new folder if necessary
GraphPath = FigurePath + ClauseName + '/' + LoadScenario + '/' + fault_name
if not os.path.exists(GraphPath):
os.makedirs(GraphPath)
# read data curves
chnfobj = dyntools.CHNF(OutputFilePath)
short_title, chanid, chandata = chnfobj.get_data()