notebook.community

Edit and run

Tutorial 1: Transient simulation and figures

Import python powerfactory module





In [1]:



    


import sys
sys.path.append(r"C:\DigSILENT15p1p7\python") # set the path to python folder inside your 
                                              # digsilent instalation folder

# import  PowerFactory  module
import powerfactory

Opening digsilent and project and calling DigSILENT commands





In [2]:



    


# start PowerFactory  in engine  mode
app = powerfactory.GetApplication('jmmauricio','') # change 'jmmauricio' by 
                                                   # your user name and '' by your password

# activate project
project = app.ActivateProject("Nine Bus System")  # change "Nine Bus System" by 
                                                  # the name of your project
    
prj = app.GetActiveProject()    # active project instance

print('prj: {:s}'.format(prj)) # to check if the project is opened properly

ldf = app.GetFromStudyCase("ComLdf")    # Load flow
ini = app.GetFromStudyCase('ComInc')    # Dynamic initialization
sim = app.GetFromStudyCase('ComSim')    # Transient simulations



    


















    






prj: <l3>\jmmauricio.IntUser\Nine Bus System.IntPrj</l3>

Getting elements and defining ouput channels





In [3]:



    


buses = app.GetCalcRelevantObjects("*.ElmTerm")
syms = app.GetCalcRelevantObjects("*.ElmSym")
loads = app.GetCalcRelevantObjects("*.ElmLod")

elmres = app.GetFromStudyCase('Results.ElmRes')

# channels for buses
for bus in buses:
    elmres.AddVars(bus,
                   'm:u',         # voltage ('m:u')
                   'm:phiu',      # angle ('m:phi')
                   'm:fehz',      # frequency ('m:fehz')
                   )   
    
# channels for synchronous machines
for sym in syms:
    elmres.AddVars(sym, 
                   's:ve',    # p.u  Excitation Voltage
                   's:pt',    # p.u.   IN    Turbine Power
                   's:ut',    # p.u.   OUT   Terminal Voltage
                   's:ie',    # p.u.   OUT   Excitation Current
                   's:xspeed',# p.u.   OUT   Speed
                   's:xme',   # p.u.   OUT   Electrical Torque
                   's:xmt',   # p.u.   OUT   Mechanical Torque
                   's:cur1',  # p.u.   OUT   Positive-Sequence Current, Magnitude
                   's:P1',    # MW     OUT   Positive-Sequence, Active Power
                   's:Q1',    # Mvar   OUT   Positive-Sequence, Reactive Power
                   'c:firel', # deg    Rotor angle with reference to reference machine angle 
                   )

# channels for loads
for load in loads:                        # creating channels for:
    elmres.AddVars(load, 'n:u1:bus1',     # voltage ('m:u')
                         'm:I1:bus1',     # current ('s:P1')
                         'm:Psum:bus1',   # active power ('s:P1')
                         'm:Qsum:bus1')   # reactive power ('s:Q1')

Generation of events





In [4]:



    


events_folder = app.GetFromStudyCase('IntEvt');  # to get events folder

evt1 = events_folder.CreateObject('EvtSwitch', 'gen 1 trip event');  # generator switch 
                                                                     # event is created here
    
evt1[0].p_target = syms[2]   # event target is the second synchronous generator
evt1[0].time = 1.0

#evt1.clear()
#evt2.clear()



    











In [5]:



    


#print(evt1[0].time)
#print(evt2)

Power flow settings and execution





In [6]:



    


ldf.iopt_net = 0
ldf.Execute()



    


















    
Out[6]:








0

Dynamic initialization





In [7]:



    


ini.iopt_dtgrd = 0.001  # step size 
ini.Execute()



    


















    
Out[7]:








0

Transient simulation





In [8]:



    


sim.tstop = 15.0
sim.Execute()



    


















    
Out[8]:








0

Results exportation





In [1]:



    


comres = app.GetFromStudyCase('ComRes'); 
comres.iopt_csel = 0
comres.iopt_tsel = 0
comres.iopt_locn = 2
comres.ciopt_head = 1
comres.pResult=elmres
comres.f_name = r'C:\Users\jmmauricio\hola.txt'  # change 'C:\Users\jmmauricio\hola.txt' 
                                                 # with your path
comres.iopt_exp=4
comres.Execute()



    


















    






---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-1-1b0285a1296b> in <module>()
----> 1 comres = app.GetFromStudyCase('ComRes');
      2 comres.iopt_csel = 0
      3 comres.iopt_tsel = 0
      4 comres.iopt_locn = 2
      5 comres.ciopt_head = 1

NameError: name 'app' is not defined

Modifications removal





In [10]:



    


evt1.clear()

Import results with pypstools





In [2]:



    


import pypstools.digsilent_simulation as ds
res = ds.ds_2_dict(r'C:\Users\jmmauricio\hola.txt')
print(res['sys'])



    


















    






---------------------------------------------------------------------------
ImportError                               Traceback (most recent call last)
<ipython-input-2-485c23e1256f> in <module>()
----> 1 import pypstools.digsilent_simulation as ds
      2 res = ds.ds_2_dict(r'C:\Users\jmmauricio\hola.txt')
      3 print(res['sys'])

ImportError: No module named 'pypstools'

Plot results with matplolib

Figure with synchronous generators speeds and powers





In [12]:



    


%matplotlib inline
import matplotlib.pyplot as plt

# create figure and axis instances
fig, (ax0, ax1) = plt.subplots(nrows=2, figsize=(10, 10))

t = res['sys']['time']  # simulation time vector

# axe 1
ax0.plot(t,res['sym']['G1']['speed']['data'], label='G1') 
ax0.plot(t,res['sym']['G2']['speed']['data'], label='G2')
ax0.plot(t,res['sym']['G3']['speed']['data'], label='G3')
ax0.legend(loc='best')
ax0.grid(True)
ax0.set_ylabel('Speed (p.u.)')
ax0.set_xlim((0,t[-1]))

# ax 2
ax1.plot(t,res['sym']['G1']['p']['data'], label='G1')
ax1.plot(t,res['sym']['G2']['p']['data'], label='G2')
ax1.plot(t,res['sym']['G3']['p']['data'], label='G3')
ax1.legend(loc='best')
ax1.grid(True)
ax1.set_ylabel('Ative Power (MW)')
ax1.set_xlim((0,t[-1]))

ax1.set_xlabel('Time (s)')



    


















    
Out[12]:








<matplotlib.text.Text at 0xebce650>

Figure with load voltages and powers





In [13]:



    


# create figure and axis instances
fig, (ax0, ax1) = plt.subplots(nrows=2, figsize=(10, 10))

t = res['sys']['time']  # simulation time vector

# axe 1
for load in res['sys']['loads']:
    ax0.plot(t,res['load'][load]['n:u1:bus1 in p.u.']['data'], label=load) 
ax0.legend(loc='best')
ax0.grid(True)
ax0.set_ylabel('Volatage (p.u.)')
ax0.set_xlim((0,t[-1]))

# ax 2
for load in res['sys']['loads']:
    ax1.plot(t,res['load'][load]['p']['data'], label=load)
ax1.legend(loc='best')
ax1.grid(True)
ax1.set_ylabel('Ative Power (MW)')
ax1.set_xlim((0,t[-1]))

ax1.set_xlabel('Time (s)')



    


















    
Out[13]:








<matplotlib.text.Text at 0xeddaef0>

Voltages at selected buses





In [14]:



    


# create figure and axis instances
fig, (ax0) = plt.subplots(nrows=1, figsize=(10, 8))

t = res['sys']['time']  # simulation time vector

buses = ['Fault',
         'Bus 1',
         'Bus 2',
         'Bus 3',
         'Bus 7',
         'Bus 8',
         'Bus 9']
# axe 1
for bus in buses:
    ax0.plot(t,res['bus'][bus]['m:u in p.u.']['data'], label=bus) 
ax0.legend(loc='best')
ax0.grid(True)
ax0.set_ylabel('Volatage (p.u.)')
ax0.set_xlim((0,t[-1]))

ax0.set_xlabel('Time (s)')



    


















    
Out[14]:








<matplotlib.text.Text at 0xe2e4db0>










    
























In [ ]:

Content source: jmmauricio/pypstools

Similar notebooks:

notebook.community | gallery | about