Revision history [back]

I am kind of new in PSS/E but here is my question:

I have a single machine infinite-bus, how can I modulate the governor or exciter of machine with a sin wave when I am doing dynamic simulation?

I found out that change_gref can change the value of the element of the GREF (governor reference). But I am not sure if I can use it to modulate the governor with sin wave.

Here is my python script:

psspy.newcase_2([0,1], 100.0, 60.0,"","")

psspy.busdata3(1,[2,i,i,i],[f,f,f,f,f,f,f],_s)

psspy.busdata3(101,[3,i,i,i],[f,f,f,f,f,f,f],_s)

psspy.buschng3(1,[i,i,i,i],[ 100.0,f,f,f,f,f,f],_s)

psspy.buschng3(101,[i,i,i,i],[ 100.0,f,f,f,f,f,f],_s)

psspy.bus_number(101,2)

psspy.branchdata(1,2,r"""1""",[i,i,i,i,i,i],[f, 0.2,f,f,f,f,f,f,f,f,f,f,f,f,_f])

psspy.branchdata(1,2,r"""2""",[i,i,i,i,i,i],[f, 0.2,f,f,f,f,f,f,f,f,f,f,f,f,_f])

psspy.plant_data(1,0,[ 1.0, 100.0])

psspy.machinedata2(1,r"""1""",[i,i,i,i,i,i],[ 50.0, 10.0,f,f,f,f,f,f,f,f, 0.2,f,f,f,f,f,f])

psspy.plant_data(2,0,[ 1.0, 100.0])

psspy.machinedata2(2,r"""1""",[i,i,i,i,i,i],[ 50.0, 10.0,f,f,f,f,f,f,f,f,f,f,f,f,f,f,_f])

psspy.loaddata4(2,r"""1""",[i,i,i,i,i,i],[ 100.0, 50.0,f,f,f,f])

psspy.save(r"""C:\Users\m\Desktop\Test\SMIB.sav""")

psspy.fdns([0,0,0,1,1,0,99,0])

psspy.fdns([0,0,0,1,1,0,99,0])

psspy.cong(0)

psspy.conl(0,1,1,[0,0],[ 100.0,0.0,0.0, 100.0])

psspy.conl(0,1,2,[0,0],[ 100.0,0.0,0.0, 100.0])

psspy.conl(0,1,3,[0,0],[ 100.0,0.0,0.0, 100.0])

psspy.fact()

psspy.tysl(1)

psspy.save(r"""C:\Users\m\Desktop\Test\SIMBconv.sav""")

psspy.dyre_new([1,1,1,1],r"""C:\Users\mg\Desktop\Test\SMIB.dyr""",r"""CC1""",r"""CT1""",r"""N1,N2,N3""")

psspy.machinearraychannel([1,1,1],r"""1""","")# Angle M1

psspy.machinearraychannel([2,7,1],r"""1""","")#speed M1

psspy.machinearraychannel([3,1,2],r"""1""","")# Angle M2

psspy.machinearraychannel([4,7,2],r"""1""","")#speed M2

psspy.voltageandangle_channel([7,4,2,1],["",""])# Voltage and Angle Bus 1

psspy.voltageandangle_channel([9,6,3,2],["",""])# Voltage and Angle Bus 2

psspy.strt(0,r"""C:\Users\m\Desktop\Test\GOP1.out""")

psspy.snap([11,6,1,0,4],r"""C:\Users\m\Desktop\Test\SMIB.snp""")

psspy.run(0,0.0,1,1,0)

psspy.run(0, 5.0,1,1,0)

Thanks.

I am kind of new in PSS/E but here is my question:

I have a single machine infinite-bus, how can I modulate the governor or exciter of machine with a sin wave when I am doing dynamic simulation?

I found out that change_gref can change the value of the element of the GREF (governor reference). But I am not sure if I can use it to modulate the governor with sin wave.

Here is my python script:

psspy.newcase_2([0,1], 100.0, 60.0,"","")

psspy.busdata3(1,[2,i,i,i],[f,f,f,f,f,f,f],_s)

psspy.busdata3(101,[3,i,i,i],[f,f,f,f,f,f,f],_s)

psspy.buschng3(1,[i,i,i,i],[ 100.0,f,f,f,f,f,f],_s)

psspy.buschng3(101,[i,i,i,i],[ 100.0,f,f,f,f,f,f],_s)

psspy.bus_number(101,2)

psspy.branchdata(1,2,r"""1""",[i,i,i,i,i,i],[f, 0.2,f,f,f,f,f,f,f,f,f,f,f,f,_f])

psspy.branchdata(1,2,r"""2""",[i,i,i,i,i,i],[f, 0.2,f,f,f,f,f,f,f,f,f,f,f,f,_f])

psspy.plant_data(1,0,[ 1.0, 100.0])

psspy.machinedata2(1,r"""1""",[i,i,i,i,i,i],[ 50.0, 10.0,f,f,f,f,f,f,f,f, 0.2,f,f,f,f,f,f])

psspy.plant_data(2,0,[ 1.0, 100.0])

psspy.machinedata2(2,r"""1""",[i,i,i,i,i,i],[ 50.0, 10.0,f,f,f,f,f,f,f,f,f,f,f,f,f,f,_f])

psspy.loaddata4(2,r"""1""",[i,i,i,i,i,i],[ 100.0, 50.0,f,f,f,f])

psspy.save(r"""C:\Users\m\Desktop\Test\SMIB.sav""")

psspy.fdns([0,0,0,1,1,0,99,0])

psspy.fdns([0,0,0,1,1,0,99,0])

psspy.cong(0)

psspy.conl(0,1,1,[0,0],[ 100.0,0.0,0.0, 100.0])

psspy.conl(0,1,2,[0,0],[ 100.0,0.0,0.0, 100.0])

psspy.conl(0,1,3,[0,0],[ 100.0,0.0,0.0, 100.0])

psspy.fact()

psspy.tysl(1)

psspy.save(r"""C:\Users\m\Desktop\Test\SIMBconv.sav""")

psspy.dyre_new([1,1,1,1],r"""C:\Users\mg\Desktop\Test\SMIB.dyr""",r"""CC1""",r"""CT1""",r"""N1,N2,N3""")

psspy.machinearraychannel([1,1,1],r"""1""","")# Angle M1

psspy.machinearraychannel([2,7,1],r"""1""","")#speed M1

psspy.machinearraychannel([3,1,2],r"""1""","")# Angle M2

psspy.machinearraychannel([4,7,2],r"""1""","")#speed M2

psspy.voltageandangle_channel([7,4,2,1],["",""])# Voltage and Angle Bus 1

psspy.voltageandangle_channel([9,6,3,2],["",""])# Voltage and Angle Bus 2

psspy.strt(0,r"""C:\Users\m\Desktop\Test\GOP1.out""")

psspy.snap([11,6,1,0,4],r"""C:\Users\m\Desktop\Test\SMIB.snp""")

psspy.run(0,0.0,1,1,0)

psspy.run(0, 5.0,1,1,0)

Thanks.