【编程实例PSS/E】自己用python写的一个三相短路故障的例子

xuning0076

import globedef
  J, P* w6 z: W' x1 [$ g7 r! qimport fault
! F5 ^2 T5 [( Y) @& z+ f$ D6 C( P+ n  gpsspy.case('c:\sav.sav')
psspy.close_report()  # close REPORT windows
options=(0,0,0,0,0,0,0)
* }/ g" w4 P, c" K, `psspy.progress_output(6, "", options)
- N& l' C& [" ], Vpsspy.report_output(6, "", options)
8 }8 @% w; d3 M6 b2 C9 tfault.three_phase(33600)
" N# M7 L. M  s' _' s
9 W  F/ R$ v+ g% b& Z; @4 G- f5 M[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
import globedef
import psspy
8 m8 b% q/ {; j% T
#a=(1,1,1)
+ @8 V! }" G7 o9 k2 Y, S- n4 G#psspy.seqd(a)

$ W8 b/ _' T5 ]& ^7 X5 K1 s' E' f8 K
6 j& S) h1 ?9 f7 C3 i/ K7 V) w0 Mdef one_phase_ground(bus,r,x):
        options=(1,1,1,0,0,0,0)
        rval=(r,x,0.0,0.0,0.0)
+ V' S1 B, q+ m( C2 f! Z0 P& F0 u0 x        cktarg=''
        psspy.scmu(1,options,rval,cktarg)
0 ?( E. c4 l7 H3 W( Z        options=(1,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
! |& i! K" Y, ^7 G& T; v  T! ]        psspy.scmu(3,options,rval,cktarg)
& r3 X. ]: j- `8 V: T/ l       
, v+ g% _5 L( D# G/ ~/ y+ w% c' K        return 0
       
& K8 c% k8 ^! K% P0 Adef two_phase(i):

        return 0
' K  L3 n# m' U, o; l
) Z, _! g; k' x: [$ rdef three_phase(bus):
        options=(1,1,1,0,0,0,0)
& [& @' j5 C* \" ]- i5 C" Q        rval=(0.0,0.0,0.0,0.0,0.0)
        cktarg=''
3 i, ?+ f/ d: U7 _        psspy.scmu(1,options,rval,cktarg)
        options=(7,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
6 R6 [  Y5 }9 [. Y2 u        return 0

2 r3 \; q, ?* \& [) ]4 E0 Ydef line_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
9 c) Z, M0 A( L, s1 Q" x: g        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,location)
, t0 Z8 L/ ^5 S+ k( F  `! k        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
" V5 y* m5 G9 \6 _' x4 j        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
( u" H# |8 u7 c! ?0 o$ D        return 0

! p0 l1 e2 J; \def line_LL_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
) j3 \" z  A5 R( F& [( n3 t        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,location)
. Q# t9 _8 H1 e/ W2 i: D3 h9 ?        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,2,1,0,0)
* d) V' k* q# I0 T0 x7 N        psspy.scmu(2,options,rval,cktarg)
) ^6 N% R% w5 H9 t- c& P& Z& B        psspy.scmu(3,options,rval,cktarg)
, ^# T5 c; S1 I( \       
        return 0

( U3 g! j. D2 udef get_breaker_cline_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
1 `/ _$ w8 \% r2 H1 I8 C9 f4 N/ v7 k# \        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
5 X& ]8 y- g3 H/ Y; E        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
: k( T6 Z9 ]& _        psspy.scmu(3,options,rval,cktarg)
# I* D6 |& p* ~6 J! i       
! K$ Y* M: h: G; Q& }4 Z        return 0

Henryzhang

咋一个三相短路写这么复杂呢？

- o% }$ Q* o7 e. ?2 _( z6 R; }. `$ n! W高手详细对比一下python和IPLAN吧
1 g3 z: \9 k4 J: S7 {! C各有哪些优势

meng823222

恩
$ F4 ?5 R0 X  j5 {0 |; F4 x! S高手说说两者各自的优缺点吧

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
& g; r  X+ R8 ]- X; n不知为什么
2 h, N" G2 B2 @/ y/ O看起来官方的dist_branch_fault、dist_bus_fault等更简洁
但没试过

Henryzhang

原先用IPLAN写branch trip：

                PUSH        'ALTR'
                PUSH        '0'
' n0 D4 ?( K5 ^$ B                PUSH        '1'        !CHANGE NETWORKDATA
2 h  i; w+ R( k1 y$ ?( g/ B6 E                PUSH        '0'
( `. I+ b# P3 c7 S9 i6 ?" L                PUSH        '3'        !BRANCH DATA
/ h% Z9 {! E) J3 x( ]) A9 Q$ v                PUSH        '151,201,1'
6 ?1 g7 q4 G; |) P$ }                PUSH        '1'
+ t, D. ?( u5 x9 k                PUSH        '0,'        !ISOLATE THIS BRANCH
                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
" c0 C9 e6 t8 f8 h1 }, u! |1 }                PUSH        ' '                !NO CHANGE IN METERED END
                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
, J, y' Y& R1 L% P8 q                PUSH  '0'   !EXIT ALTR ACTIVITY

现在python：
dist_branch_trip(151,201,'1')
9 k0 M2 \0 ^9 |  `1 U一句就可以了
' d& v1 G$ y" }, H
8 a5 y/ t: a$ |7 ]3 @! x; S' ~确实比较简洁

loadflow

亨利版主是ＩＰＬＡＮ的行家。觉得ｐｙｔｈｏｎ很多东西藏的太深了。不像直接调用ａｃｔｉｖｉｔｉｅｓ这样明显。

Henryzhang

用官方提供的应该可靠性高点
就是不容易出错

# }0 O% ]8 i2 p3 n# R- V( b: K要想玩PSSE，还是基于最基本命令单元，呵呵

向前向前

直接看不懂呀