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

xuning0076

import globedef
import fault
& X/ T8 |& y% J+ @: a. tpsspy.case('c:\sav.sav')
8 }' z8 w: @" F, W1 gpsspy.close_report()  # close REPORT windows
) j6 e) q& W1 Toptions=(0,0,0,0,0,0,0)
psspy.progress_output(6, "", options)
psspy.report_output(6, "", options)
* J3 c( l4 ?$ O' p9 j7 Tfault.three_phase(33600)

1 ]* N; C0 ?. T1 z5 J[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
import globedef
4 _" A5 a4 w! o) ^7 d# @) Kimport psspy

2 d, f! B4 c7 a. G+ p( h#a=(1,1,1)
#psspy.seqd(a)
# |: E' U" [) |. w1 V
0 P: b# I# P+ G4 ^# ^4 D! h, B
3 k1 S, o, I/ ~' Vdef one_phase_ground(bus,r,x):
        options=(1,1,1,0,0,0,0)
# w9 c# H7 ^: n% T        rval=(r,x,0.0,0.0,0.0)
        cktarg=''
( V3 P! }& u% G) q4 N' S5 V* _        psspy.scmu(1,options,rval,cktarg)
        options=(1,bus,1,0,0,0,0)
% H" Z+ M! e* v  C4 t4 O- x6 k        psspy.scmu(2,options,rval,cktarg)
; ]* X" i7 U( z  o1 D- ?        psspy.scmu(3,options,rval,cktarg)
. A/ q( F8 }( B( p       
9 C- B  u+ Y! ?' b7 W+ K' f' y        return 0
       
def two_phase(i):
* H2 p5 ]" |5 ?" O+ s( _3 t
; x$ V5 B) r* }1 }/ r' e/ d        return 0
8 [6 V6 Y, ~% ~$ t
* c9 b$ B# b& d9 d# G3 b& S: V8 zdef three_phase(bus):
        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,0.0)
        cktarg=''
        psspy.scmu(1,options,rval,cktarg)
        options=(7,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
8 F4 M. R# q) O        psspy.scmu(3,options,rval,cktarg)
9 S$ c# ?0 S( }) m% Z! G       
; l1 B# c" |0 ~2 w% O1 g        return 0

def line_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
& R. q; k2 B9 b" i( K        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
/ q3 |% y3 ~4 m2 ~       
        return 0

def line_LL_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
' }" b2 o! h, f) k; L        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
$ y8 e: g% `/ K4 c* U$ v2 V' U. D. K        psspy.scmu(1,options,rval,cktarg)
1 S8 L' F- |/ q" V7 }! {5 y. O        options=(9,bus_i,bus_j,2,1,0,0)
+ \4 ]" |4 p; [( l7 Z; |: a        psspy.scmu(2,options,rval,cktarg)
- w5 U  T* Q7 z        psspy.scmu(3,options,rval,cktarg)
       
        return 0

def get_breaker_cline_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
  D# }: p5 [1 x. s, F+ `* s        rval=(0.0,0.0,0.0,0.0,location)
0 B4 Z/ Q7 H1 i2 X1 ~- N3 d; V9 X        cktarg=ckt
' u4 |/ H" c: M9 n( L        psspy.scmu(1,options,rval,cktarg)
: }$ v6 [1 N4 E6 Q        options=(9,bus_i,bus_j,1,1,0,0)
( H6 ~* G# K# d. k. w        psspy.scmu(2,options,rval,cktarg)
6 S3 d. R2 {6 H3 l4 P2 }) H9 [4 x        psspy.scmu(3,options,rval,cktarg)
       
        return 0

Henryzhang

咋一个三相短路写这么复杂呢？
7 {- N1 X6 B! V* S
高手详细对比一下python和IPLAN吧
各有哪些优势

meng823222

恩
5 o. b/ w: ?' k2 v  d4 m高手说说两者各自的优缺点吧

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
不知为什么
看起来官方的dist_branch_fault、dist_bus_fault等更简洁
7 _" w- q  b2 o" t5 N9 \, O6 x3 P但没试过

Henryzhang

原先用IPLAN写branch trip：

                PUSH        'ALTR'
                PUSH        '0'
                PUSH        '1'        !CHANGE NETWORKDATA
                PUSH        '0'
                PUSH        '3'        !BRANCH DATA
& u4 O: g8 _) }( z                PUSH        '151,201,1'
                PUSH        '1'
                PUSH        '0,'        !ISOLATE THIS BRANCH
                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
                PUSH        ' '                !NO CHANGE IN METERED END
                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
/ p( ^4 B" X' _# }                PUSH  '0'   !EXIT ALTR ACTIVITY

: G4 G  G& t3 ^6 D9 ^) Q现在python：
, N/ S2 H, z7 L$ A4 ndist_branch_trip(151,201,'1')
& `1 x$ v4 S2 P1 b" A$ p一句就可以了

* L8 ], c! n* {4 C. y5 `确实比较简洁

loadflow

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

Henryzhang

用官方提供的应该可靠性高点
$ k# S/ K( e" f- u0 I1 _+ `' g' g  K- c就是不容易出错
/ t9 j+ R9 t  z) N7 M
要想玩PSSE，还是基于最基本命令单元，呵呵

向前向前

直接看不懂呀