马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
我最近在用Python调用PSS/E动态模拟计算.& N9 Y6 i6 u# ?0 j) l- U5 ]
我的版本是:! T3 {' y7 V7 L
Python 2.5# ]2 |" [6 f' n" G( u8 R
PSS/E 31
# S7 m" R" z0 L! r7 c$ w在31之前的版本是不支持Python的动态模拟调用的.
0 [- R( _) C' q9 Y3 b# x3 r6 \* I; C$ r
我的一个简单算例:
0 W% q2 R# Q" S/ i) i 0 L, `3 `1 i, ~4 s
import os,sys 1 [3 r" D2 d2 V' v
sys.path.append('C:\\Program Files\\PTI\\PSSE31\\PSSBIN')
6 s& S/ G6 D. u3 w; C- v; r' M os.environ['PATH'] += ';C:\\Program Files\\PTI\\PSSE31\\PSSBIN'
; u' H1 v" Y( C/ o$ L( w
% S' z1 M2 ?) |/ K4 ~ import redirect$ q6 P d' X [7 }, j8 e
redirect.psse2py()
% F ]! b: z/ y0 {2 a8 e2 w import psspy1 A1 N/ ?8 V- {
psspy.psseinit(80000)
6 _% \5 G! E5 ~ psspy.progress_output(islct=6) ( s- j6 D3 \9 `
psspy.prompt_output(islct=6)
9 e) ?* b8 K) L, q% _9 z0 K$ m& C) a3 e) b' |8 F, J
##open power flow data, convert and solve it
8 h# ^8 G+ f0 d* Z os.chdir('C:\\Python25\\pp_test\\martin\\case1')
7 w- ?/ O4 b% s6 k8 Y6 ` F" M psspy.case(r"""savcnv.sav""") ##this case is pre-conveted case
# n0 P K# ^) S$ Q5 F- C0 v" j4 J9 N; g' F. Y' e1 x& M
##load snap file ! Q% \+ I) c; P7 d
psspy.dynamicsmode(1)
9 K5 F% L: G1 `+ E7 v4 C psspy.rstr(r"""savnw.snp""")# d# C/ f* ~' X7 `0 k: }
' k' S9 a$ |' g' V' d$ |' I6 g ierr = psspy.strt(0, 'test.out')) U/ B+ l& x: S2 [. w2 d4 B$ I% R7 F
ierr = psspy.run(0, 10)$ P0 ^; p$ ^0 s& K% Q( t
5 g( f& V5 B1 z! R$ K5 _/ @6 n" }
如果要加故障操作,可以用:( D6 e9 R$ d) g. J. N
ierr = psspy.dist_bus_fault(faultbus, 1, 0, [0,-20000000000])+ s% @; l" Q/ q! \* h& m7 j E
ierr = psspy.dist_clear_fault(1)# J4 R. Y3 s$ Q
ierr = psspy.dist_branch_trip(faultbus, nb_bus,'1')0 \0 x7 a0 g& z7 P
/ H5 ^. n2 h6 K% y7 }
[ 本帖最后由 dsyljh 于 2008-7-28 18:21 编辑 ]-%-89529-%--%-97274-%--%-89538-%- |