Revision history [back]

I have simulated Alberta power system successfully. To do this, you have to run file 2014AESODynamicGNETR1.idv before import dyr file and reduce the time step of the simulation to 0.002s. Here is python code for this simulation: (you can change output channel and contigencies if you want)

psspy.cong(0)
psspy.conl(_i,_i,1,[0,_i],[_f,_f,_f,_f])
psspy.conl(1,1,2,[_i,_i],[ 20.0, 30.0, 20.0, 30.0])
psspy.conl(_i,_i,3,[_i,_i],[_f,_f,_f,_f])
psspy.ordr(0)
psspy.fact()
psspy.tysl(0)
psspy.dyre_new([1,1,1,1],r"""D:\Huy\Alberta power system\2014_AESO_Dynamic_R5\2014_AESO_Dynamic_V33_R1.dyr""","","","")
psspy.bsysinit(1)
psspy.bsys(1,0,[0.0,0.0],0,[],1,[19241],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[18307],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[18351],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[19351],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[18420],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[2674],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[2954],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[2358],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[2603],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[3227],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4205],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4227],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4735],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4740],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[536],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[17101],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[3043],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4043],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4060],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4997],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[1681],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[11013],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[19834],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[61683],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[61686],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[3085],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4127],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[25352],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4017],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[3272],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4269],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4690],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[2556],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4235],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4275],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[16834],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[17834],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[18834],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[80000],0,[],0,[])
psspy.bsysadd(1,0,[0.0,0.0],0,[],1,[4105],0,[],0,[])
psspy.gnet(1,0)
psspy.dynamics_solution_param_2([_i,_i,_i,_i,_i,_i,_i,_i],[_f,_f, 0.002,_f,_f,_f,_f,_f])
psspy.chsb(0,1,[1,3092,1247,1,7,0])
psspy.strt(0,r"""1234""")
psspy.run(0, 2.5,500,5,10)