Revision history [back]

Hello,

I just started working extensively with the dynamic simulator in PSS/E 33, and I am currently running a flat (No Contingencies), 1 second, dynamic simulation of a simple 9 bus (3 gens, 3 loads) case. I have loaded all the .raw data and .dyr data into the case via python and the PSSE API (see applicable source code below).

Does anyone have any idea why the PU bus voltages blow up to infinity? I must be missing something somewhere - This exact case simulates perfectly in PowerWorld.

Here is some of the applicable source code:

p.case(static_case_file) #Loads a succesfully solved (with FNSL) case. (Verified against PowerWorld output)
##convert to dynamic case
p.conl(-1,1,1)  #setup
p.conl(-1,1,2,[0,0],[100,0,0,100]) #convert loads
p.conl(-1,1,3) #cleanup
p.cong() # convert Gens
p.ordr() # preserve sparsity of network matrices
p.fact() # factorize the network admittance matrix
p.tysl() # switching study network solutions
## Save dynamic case binary file
p.save(dynamic_case_file)
## load dynamic model definitions
p.dyre_new([1,1,1,1], dyr_infile) # loads dynamic parameters in the .dyr file overwrites existing dynamic data
## Save binary file for dynamic information
p.snap([-1,-1,-1,-1,-1], dynamic_snap_file) 
##LOAD Saved Files into working memory
p.case(dynamic_case_file)
p.rstr(dynamic_snap_file)
##Setup Dynamic Sim
p.chsb(0,1,[-1,-1,-1,1,13,0])
## strt(Option, Outfile)  execute before 'run'
p.strt(0, dynamic_channel_outfile)
## run(OPTION, TPause, numTimeStepsBTWPrint, numTimeStepsOUT, NumtimeStepPLT)
p.run(tpause = .1) # initial run to stabalize (suggested in PSSe)
p.run(tpause = 1)  # proceding runs...

Following are the generator models that I am using in the dyr_infile:

1, 'GENSAL', 1, 8.96, .0333, .0333, 23.64, 23.64, .146, .0969, .0608, .05, 0, 0, 0 /
2, 'GENSAL', 1, 6, .0333, .0333, 6.4, 6.4, .8958, .8645, .1198, .05, 0, 0, 0 /
3, 'GENSAL', 1, 5.89, .0333, .0333, 3.01, 3.01, 1.3125, 1.2578, .1813, .05, 0, 0, 0 /

Finally here is the channel_outfile graphed (and zoomed in...)

CANT UPLOAD THIS UNTIL I have more Karma.......

A few other notes:

1. I have already verified that complex machine impedance (Zx of ZSOURCE) matches the Xd" and Xq" value in the GENSAL model

2. For simplicity sake I am modeling 3 transformer branches as simple Non-Transformer Branches.

I am brand new to PSSE, and would really appreciate any thoughts on what might be causing this. I assume there is some parameter causing the differential equations in the dynamic simulation to blow up - or I am missing a command/step in allowing PSSE to completely understand the topology. Any thoughts on how to troubleshoot this beast would also be helpful. In the meantime, I will continue reading through the loads of documentation...

thanks

Hello,

I just started working extensively with the dynamic simulator in PSS/E 33, and I am currently running a flat (No Contingencies), 1 second, dynamic simulation of a simple 9 bus (3 gens, 3 loads) case. I have loaded all the .raw data and .dyr data into the case via python and the PSSE API (see applicable source code below).

Does anyone have any idea why the PU bus voltages blow up to infinity? I must be missing something somewhere - This exact case simulates perfectly in PowerWorld.

Here is some of the applicable source code:

p.case(static_case_file) #Loads a succesfully solved (with FNSL) case. (Verified against PowerWorld output)
##convert to dynamic case
p.conl(-1,1,1)  #setup
p.conl(-1,1,2,[0,0],[100,0,0,100]) #convert loads
p.conl(-1,1,3) #cleanup
p.cong() # convert Gens
p.ordr() # preserve sparsity of network matrices
p.fact() # factorize the network admittance matrix
p.tysl() # switching study network solutions
## Save dynamic case binary file
p.save(dynamic_case_file)
## load dynamic model definitions
p.dyre_new([1,1,1,1], dyr_infile) # loads dynamic parameters in the .dyr file overwrites existing dynamic data
## Save binary file for dynamic information
p.snap([-1,-1,-1,-1,-1], dynamic_snap_file) 
##LOAD Saved Files into working memory
p.case(dynamic_case_file)
p.rstr(dynamic_snap_file)
##Setup Dynamic Sim
p.chsb(0,1,[-1,-1,-1,1,13,0])
## strt(Option, Outfile)  execute before 'run'
p.strt(0, dynamic_channel_outfile)
## run(OPTION, TPause, numTimeStepsBTWPrint, numTimeStepsOUT, NumtimeStepPLT)
p.run(tpause = .1) # initial run to stabalize (suggested in PSSe)
p.run(tpause = 1)  # proceding runs...

Following are the generator models that I am using in the dyr_infile:

1, 'GENSAL', 1, 8.96, .0333, .0333, 23.64, 23.64, .146, .0969, .0608, .05, 0, 0, 0 /
2, 'GENSAL', 1, 6, .0333, .0333, 6.4, 6.4, .8958, .8645, .1198, .05, 0, 0, 0 /
3, 'GENSAL', 1, 5.89, .0333, .0333, 3.01, 3.01, 1.3125, 1.2578, .1813, .05, 0, 0, 0 /

Finally here is the channel_outfile graphed (and zoomed in...)

CANT UPLOAD THIS UNTIL I have more Karma.......

A few other notes:

1. I have already verified that complex machine impedance (Zx of ZSOURCE) matches the Xd" and Xq" value in the GENSAL model

2. For simplicity sake I am modeling 3 transformer branches as simple Non-Transformer Branches.

I am brand new to PSSE, and would really appreciate any thoughts on what might be causing this. I assume there is some parameter causing the differential equations in the dynamic simulation to blow up - or I am missing a command/step in allowing PSSE to completely understand the topology. Any thoughts on how to troubleshoot this beast would also be helpful. In the meantime, I will continue reading through the loads of documentation...

thanks