Returning Dynamic simulation data.

asked 2018-05-19 18:36:41 -0600

Chris.R
1 ●1 ●1 ●2

updated 2018-05-19 18:50:51 -0600

Hello,

I am trying to run a dynamic simulation and return the data from each time step rather than outputting it to the .out file so that it can be used in Matlab. Is there a function or way that I would be able to do that?

dryfile, Outfile, and subsystem 2 are all defined earlier in the code and are working fine.

#initiate DYR    
psspy.dyre_new([1,1,1,1],dyrfile)
psspy.chsb(2,1,[-1,-1,-1,1,13,0]) #voltage at bus

#convert loads and generators    
psspy.cong(0) 
psspy.conl(-1,1,1)
psspy.conl(-1,1,2,[0,0],[0,100,0,100])
psspy.conl(-1,1,3)

psspy.ordr()
psspy.fact()
psspy.tysl()

psspy.strt(0,OutFile)
psspy.run(0,5,1,1,0)
psspy.dist_bus_fault(5,3, 0.0, [.1,0]); #SET Fault at bus 5
psspy.run(0,5.1,1,1,0)
psspy.dist_clear_fault(1)

t = 5.1

for TS in range(600): #runs to roughly 10 seconds
    t=t+0.0087
    psspy.run(tpause=t)
    #output voltage (or other CHSB data) as a variable for matlab use at each pause point. 
    pausePoint = input('pause here for Matlab to be run')

I want to set it up so that bus voltage is returned at each time step in the loop after the fault occurs so that some math can be done on with it before initiating the next time step.

edit retag flag offensive close merge delete

add a comment

3 answers

Sort by » oldest newest most voted

answered 2018-05-21 04:06:28 -0600

perolofl

3332 ●1 ●8 ●19

What kind of math are you going to apply to the bus voltage?

Do you want to manipulate PSSE during the simulation? If so, which parameter in PSSE will be updated based on the Matlab calculation?

edit flag offensive delete link

Comments

The objective is to run a single time step in PSSE, take the output values of the dynamic simulation and use them in Matlab to run Kalman Filter. Then apply the result of KF in PSSE simulation and run the next PSSE time step. In short I want to run KF in every time step of PSSE dynamic simulation.

Chris.R ( 2018-05-21 16:36:27 -0600 )edit

Hi Chris, I met the same problem and have no idea how to work out. Did you find a way to return the voltage value of each step into a variable?

ywan1246 ( 2021-04-18 19:44:35 -0600 )edit

Hi ywan1246 and Chris, you can use the ierr, rarray = abusreal(sid, flag, string) python API function to get different parameters at each time step during a simulation. You can look into detail in Ch 12 (subsystem data retrieval) of API documentation (PSS/E ver. 35)

powerengineer ( 2022-09-27 21:13:08 -0600 )edit

add a comment

answered 2018-05-20 04:27:51 -0600

AlexYang

1 ●6 ●7 ●8

updated 2018-05-20 04:34:27 -0600

You can use dyntool to export output file to excel file, then using matlab process excel file. Or catch dyntool objection, use matplotlib package in python to plot what you want.

You can refer to this script [search: Extract data with dyntools and plot with matplotlib]：

import dyntools
import matplotlib.pyplot as plt

outfile = r"""test.out"""
chnfobj = dyntools.CHNF(outfile)
short_title, chanid, chandata = chnfobj.get_data()

t = chandata['time']
v = chandata[1]              #based on the channel identifier set in the dynamic simulation

fig = plt.figure()
fig.patch.set_facecolor('0.8')
ax1.plot(t,v,linestyle='-', linewidth=1, color='green',label="Voltage")
plt.grid(linestyle='--', color='grey',linewidth=0.5)
plt.xlabel("Time")
plt.ylabel("Voltage")
plt.legend()
axes = plt.gca()
axes.set_facecolor('w')
axes.set_ylim([0.6,1.2])
axes.set_xlim([0,10])
plt.savefig(r"""figure.pdf""",dpi=fig.dpi,facecolor='0.8')
plt.show()

edit flag offensive delete link

Comments

Completely off topic!

perolofl ( 2018-05-22 02:10:20 -0600 )edit

add a comment

answered 2018-05-21 17:54:57 -0600

jconto
2701 ●10 ●30 ●41

AlexYang's code can run outside the PSSe GUI! I added code to export to XLS, change ax1.plot to plt.plot, and commented "axes.set_facecolor('w')" to make it work on my pc:

import dyntools
import matplotlib.pyplot as plt

outfile = r"""test.out"""
chnfobj = dyntools.CHNF(outfile)
short_title, chanid, chandata = chnfobj.get_data()

t = chandata['time']
v = chandata[1]              #based on the channel identifier set in the dynamic simulation
XLSfile = 'voltage.xls'
sheetname= 'V_chan1'
chnfobj.xlsout(channels=[1,],outfile=outfile,xlsfile=XLSfile,sheet=sheetname,show=False)

fig = plt.figure()
fig.patch.set_facecolor('0.8')
plt.plot(t,v,linestyle='-', linewidth=1, color='green',label="Voltage")
plt.grid(linestyle='--', color='grey',linewidth=0.5)
plt.xlabel("Time")
plt.ylabel("Voltage")
plt.legend()
axes = plt.gca()
#axes.set_facecolor('w')
axes.set_ylim([0.6,1.2])
axes.set_xlim([0,10])
plt.savefig(r"""figure.pdf""",dpi=fig.dpi,facecolor='0.8')
plt.show()

edit flag offensive delete link

add a comment

Returning Dynamic simulation data.

3 answers

Comments

Comments

Your Answer

Question Tools

Stats

Related questions