Python for IPSA (Power flow analysis)
Dave Angel
davea at davea.name
Tue May 28 07:14:37 EDT 2013
On 05/28/2013 06:00 AM, Debbie wrote:
> Hi there,
> I am new to Python,
Welcome. Could you tell us a little about yourself, such as whether
you've experience in a few other languages, or if Python is your first
programming experience? Also, what version of Python (presumably 2.7 or
2.6) and what OS ?
> and wondering if you could help me with python based coding for the IPSA (Power system analysis software). I have a electrical distribution network with generators, buses and loads, on which I am performing the load flow analysis every 1 hour for a period of 1 year.
>
> The code to perform instantaneous load/power flow analysis is given below. I need to modify it such that I can perform this load flow analysis every 1 hour for a period of 1 year. Please help.
>
> from ipsa import *
I have no idea what functionality is in ipsa, so this whole message is a
guess. First question is whether the ReadFile() below and the
DoLoadFlow() read the whole year's data, or data for a particular hour.
And if the latter, how do you get the next set of data?
>
> ipsasys = IscInterface()
> net = ipsasys.ReadFile("refinery.iif")
> bok = net.DoLoadFlow();
> if bok:
> busbars = net.GetBusbars()
> print "Load Flow results:"
> print ""
> print "BusName Vmag(kV)"
> print "===================="
> for bus in busbars.itervalues():
> name = bus.GetName()
> vm = bus.GetVoltageMagnitudekV()
> res = "%-8s %10.5f" % (name, vm)
> print res
> else:
> print "Load Flow failed!"
>
> Regards,
> Debbie
>
First you want to make a function to print out a particular hour's data.
That might turn out to be something like:
def one_hour(net):
busbars = net.GetBusbars()
print "Load Flow results:"
print ""
print "BusName Vmag(kV)"
print "===================="
for bus in busbars.itervalues():
name = bus.GetName()
vm = bus.GetVoltageMagnitudekV()
res = "%-8s %10.5f" % (name, vm)
print res
(Just pasted from your code, I added in the probable parameter to the
definition)
Now your main function might be something like:
def main(filename):
ipsasys = IscInterface()
net = ipsasys.ReadFile(filename)
#if net gets you data for one hour, figure out how
#to get the data for the whole year, in the form
#of a list or an iterator called nets
bok = net.DoLoadFlow();
if bok:
for net in nets:
one_hour(net)
else:
print "Load Flow failed!"
#and your top-level code is:
main("refinery.iif")
As to ordering in the source file, put the import first:
from ipsa import *
then your function definitions one_hour() and main(), then
your top-level code.
--
DaveA
More information about the Python-list
mailing list