amaity's profile - activity

Hi All!

Please correct/improve this QV curve generating script.

# PSS/E Saved case
CASE = r"C:\Program Files\PTI\PSSE32\EXAMPLE\12TH_PLAN_MAR_2017R1.sav"
psspy.psseinit(12000)
psspy.case(CASE)

def getdata(subsystem_data):
    ierr, values = subsystem_data
    return zip(*values)

#--------------------------------
ierr = psspy.bsys(sid=1, numzone=1, zones=[44])
buses = reduce(tuple.__add__, getdata(psspy.abusint(sid=1, string="NUMBER")))
buskv = reduce(tuple.__add__, getdata(psspy.abusreal(sid=1, string="BASE")))
#print zip(buses, buskv)
selected_buses = []
for bus, voltage in zip(buses, buskv):
    if voltage > 100.0:
        selected_buses.append((bus))
print selected_buses

busno = 44999 # Fictitious generator bus
genid = 1
status = 1
pgen = 0.0 # Fict gen P output
Qlimit = 9999.0 # Fict. gen Q limit
pmax = 0.0 # Fict gen P limit
ireg = 44106 # Regulated bus for which QV curve is to be determined

def add_machine():
    psspy.plant_data(busno, intgar1=ireg)
    psspy.machine_data_2(
        busno,
        str(genid),
        [int(status), 1, 0, 0, 0, 0],
        [pgen, 0.0,
        Qlimit,
        -Qlimit,
        pmax, 0.0, 100.0, 0.0, 1.0, 0.0, 0.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0])
    return

def get_mvar(i):
    psspy.plant_data(busno, realar1=i)
    psspy.solution_parameters_3(intgar2=60)
    ierr = psspy.fnsl()
    mvar = psspy.macdat(busno, str(genid), 'Q')
    return mvar[1]

psspy.bus_data_2(busno, intgar1=2, name='TEST')
psspy.branch_data(i=busno, j=ireg)
add_machine()

import numpy as np
import matplotlib.pyplot as plt
pu = [x for x in np.arange(.8, 1.2, .05)]
varlist = [get_mvar(v) for v in pu]
print varlist
plt.plot(pu, varlist, '-o')
plt.xlabel('PU')
plt.ylabel('MVar')
plt.grid()
plt.show()

Here is a back to basics question. Let me know if this seems irrelevant in which case I will withdraw it. How do you determine the line constants? Do you determine them using your own scripts or use the TLC module built into psse? Is the conductor database in TLC adequate for your needs or do you use your own database?

This link should be able to tell you all that you are looking for. I wouldn't be able to explain it better than this even if I tried.

Can you give us the transformer data? Perhaps I could plug it into my system and see why PSSE crashes.

Umm ..... I guess I have missed something in the psse manual. Hey folks, is there something in the manual that says that automatic contingencies (like SINGLE LINE FROM ...) cannot be mixed with single event contingencies (like TRIP BRANCH FROM ...) ???

I am having problems running a con file with the following contingencies in one grouping: REMOVE UNIT 1 FROM BUS 101, SINGLE BRANCH FROM BUS 101, END, END. The psse output says there is a problem at SINGLE. What am I doing wrong??

Hi @Eli Pack, can we mix individual trips/disconnects/remove with an automatic outage in one contingency group?

What are the names of the example files?

No, your link doesn't work.

A pyd extension means it's probably a dll file. So how are you using it ??? Can you share the documentation and the example files?

Not at all. I don't think you need to register to access the public folders in github. Folks correct me if I am wrong. How you can adapt the function for your needs has been explained in the function description. Just plug 'trn' or 'brn' in place of 'bus'. HTH.

Use the subsystem_info function described here. The full listing is available here. This little function saves you the bother of distinguishing lines from transformers. This is how I use it:

branchinfo = subsystem_info('brn', ['FROMNUMBER', 'TONUMBER', 'MVA', 'P'], sid=1)
#print branchinfo
trfinfo = subsystem_info('trn', ['FROMNUMBER', 'TONUMBER', 'MVA', 'P'], sid=1)
#print trfinfo

Hope that helps.

Which version of PSSE is this in and in what directory? I can't find it my PSSE v33.5 install directory.

What is the API for extracting the admittance matrix?

Actually there is a module called the Graphical Model Builder. Look here.

I think this is what you are looking for:

@rich how are you doing this? Have you written a script for this kind of dynamic study?