我正在尝试运行我在PSSE上制作的自动化程序。我几乎从excel表访问数据并将其传输到PSSE。代码如下所示。有没有办法解决这个问题，更改程序设置或任何解决方法？任何帮助，将不胜感激。我使用的函数是psspy.fdns，我需要使用固定斜率解耦Newton Raphson方法。并且在代码中多次显示。
from __future__ import print_function
import os.path
import win32com.client
#---------------------------------------------------------------------------------------------------------------------
xlApp = win32com.client.DispatchEx('Excel.Application') # Running Excel
xlsPath = os.path.expanduser('I:\TPGeneral\Ten Year Load Forecasts 2017-2026.xlsm')# Reading xlsm file
wb = xlApp.Workbooks.Open(Filename=xlsPath) # Opening file
xlApp.Run('csvfile4')# Running macro---- csvfile2 is the macro name. It is under the "csv" module in the VBA editor
wb.Save()
xlApp.Quit()
#---------------------------------------------------------------------------------------------------------------------
import os, sys
import psspy
import redirect
import csv
psspy.throwPsseExceptions = True
from Tkinter import * # File Dialog
import tkFileDialog
import tkSimpleDialog
import tkMessageBox
import Tkinter as tk
LOAD_GEN_DATAFILE = 'I:\TPGeneral\Load_Forecasts.csv' # CSV File to Read
#-------------------------------------------------------------------------------
# read the entire CSV into Python.
# CSV has columns starting with Year,busnum,busname,scaled_power,tla,location
data = list(csv.reader(open(LOAD_GEN_DATAFILE)))
mydict = {}
for row in data:
Year,busnum,busname,scaled_power,tla,data_location,empty,year_link,from_,to,digit,name2,tla_2,min_value,max_value = row[0:15]
#If this is a year not seen before, add it to the dictionary
if Year not in mydict:
mydict[Year] = {}
busses_in_year = mydict[Year]
if data_location not in busses_in_year:
busses_in_year[data_location] = []
#Add the bus to the list of busses that stop at this location
busses_in_year[data_location].append((busnum,busname,scaled_power))
#-------------------------------------------------------------------------------
#-------------------------------------------------------------------------------------------------------------------------------------
#User Input Statement
year = raw_input("Please Select Year of Study: ")
print("\n")
location=raw_input(" \n The list above show the TLA Pockets as well as the ID numbers assigned to them ()\n\n Please enter the ID #: ")
print("\n")
Year=year
data_location=location
#---------------------------------------------------------------------------------------------------------------------------------------
#-------------------------------------------------LOAD FORECASTS--------------------------------------------------------------------
outputs=[]
if Year in mydict and data_location in mydict[Year]:
busses_in_year = mydict[Year]
#print("Here are all the busses at that location for that year and the new LOAD TOTAL: ")
#print("\n")
#Busnum, busname,scaled_power read from excel sheet matching year and location
for busnum,busname,scaled_power in busses_in_year[data_location]:
scaled_power= float(scaled_power)
busnum = int(busnum)
output='Bus #: {}\t Area Station: {}\t New Load Total: {} MW\t'
formatted = output.format(busnum, busname, scaled_power)
outputs.append(formatted)
psspy.bsys(1,0,[0.0,0.0],0,[],1,[busnum],0,[],0,[])
psspy.scal_2(1,0,1,[0,0,0,0,0],[0.0,0.0,0.0,0.0,0.0,0.0,0.0])
psspy.scal_2(0,1,2,[0,1,0,1,0],[scaled_power,0.0,0,-.0,0.0,-.0,0])
#psspy.fdns([0,0,0,0,0,0,0,0])
#psspy.fdns([1,0,0,0,0,0,0,0])
#psspy.fdns([1,0,1,0,0,0,0,0])
else:
exit
psspy.fdns([0,0,0,0,0,0,0,0])
psspy.fdns([1,0,0,0,0,0,0,0])
psspy.fdns([1,0,1,0,0,0,0,0])
#-------------------------------------------------LOAD FORECASTS--------------------------------------------------------------------
#------------------------------------------------VOLTAGE ENVLOPES---------------------------------------------------------------------------------------------------------------------------------
outputs_3=[]
for row in data:
year_link, from_,to,digit,name2,tla_2,min_value,max_value = row[7:15]
output = 'From Bus #: {}\tTo Bus #: {}\t Area Station: {}\t VMIN: {} pu\tVMAX: {} pu\t'
if year_link == year and tla_2==location:
from_=int(from_)
to=int(to)
min_value=float(min_value)
max_value=float(max_value)
digit=int(digit)
output = 'From Bus #: {}\tTo Bus #: {}\t Area Station: {}\t VMIN: {} pu\tVMAX: {} pu\t'
formatted_3 = output.format(from_, to,name2, min_value, max_value)
outputs_3.append(formatted_3)
_i=psspy.getdefaultint()
_f=psspy.getdefaultreal()
psspy.two_winding_chng_4(from_,to,'%d'% digit,[_i,_i,_i,_i,_i,_i,_i,_i,_i,_i,_i,_i,_i,_i,_i],[_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f,_f, max_value, min_value,_f,_f,_f],[])
psspy.fdns([0,0,0,0,0,0,0,0])
psspy.fdns([1,0,0,0,0,0,0,0])
psspy.fdns([1,0,1,0,0,0,0,0])
else:
exit
#psspy.fdns([0,0,1,1,0,0,0,0])
outputs_2=[]
for row in data:
year_link, from_,to,digit,name2,tla_2,min_value,max_value,last_bus = row[7:16]
if isinstance(last_bus, str) and year_link==year and tla_2==location and last_bus.isdigit() is True:
min_value=float(min_value)
max_value=float(max_value)
from_=int(from_)
to=int(to)
last_bus=int(last_bus)
output = '3 WINDING VOLTAGE ENVELOPE\tFrom Bus #: {}\tTo Bus #: {}\tLast Bus #: {}\t Area Station: {}\t VMIN: {} pu\tVMAX: {} pu\t'
formatted_2 = output.format(from_, to,last_bus,name2, min_value, max_value)
outputs_2.append(formatted_2)
_i=psspy.getdefaultint()
_f=psspy.getdefaultreal()
psspy.three_wnd_winding_data_3(from_,to, last_bus,r"""1""",1,[_i,_i,_i,_i,_i],[_f,_f,_f,_f,_f,_f,_f,_f, max_value, min_value,_f,_f,_f])
psspy.fdns([0,0,0,0,0,0,0,0])
psspy.fdns([1,0,0,0,0,0,0,0])
psspy.fdns([1,0,1,0,0,0,0,0])
else:
exit
#------------------------------------------------VOLTAGE ENVLOPES----------------------------------------------------------------------------------------------------------------------------------
print('\n','---------------LOAD FORECASTS---------------','\n')
s='\n'*2
print(s.join(outputs))
print('\n','---------------2 WINDING VOLTAGE ENVELOPES---------------','\n')
print(s.join(outputs_3))
print('\n')
print(s.join(outputs_2))
os.remove(LOAD_GEN_DATAFILE)