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

xuning0076

import globedef
import fault
6 U1 c: X1 `+ R+ h2 q4 bpsspy.case('c:\sav.sav')
" U" l* E5 b6 m* }2 H( k: ypsspy.close_report()  # close REPORT windows
: Q/ r7 p" |, n) m) qoptions=(0,0,0,0,0,0,0)
- F6 ~* }7 g  N4 V, m2 Apsspy.progress_output(6, "", options)
psspy.report_output(6, "", options)
; A: x& _- n& M! w- Kfault.three_phase(33600)
, L6 j: K! A  L6 j$ V0 v
[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
import globedef
. x! `& C3 r$ t. K0 G7 ?import psspy
7 M: z# ^; }. X( ]9 B* x8 U- }
#a=(1,1,1)
#psspy.seqd(a)
/ ^1 o: q: ?+ T/ }
9 P( m7 n8 M) c: B; U0 Z+ ^( t  }! ^
8 [  Q2 v" W8 D& c' G& ldef one_phase_ground(bus,r,x):
) v3 J1 {- M  [8 G8 q' `        options=(1,1,1,0,0,0,0)
6 I/ q' ~/ b% v/ f; q% X/ Q9 L* X  I        rval=(r,x,0.0,0.0,0.0)
1 a' F6 |0 E) V* V7 m        cktarg=''
        psspy.scmu(1,options,rval,cktarg)
% @2 D& P: a. {, m2 L4 `/ I        options=(1,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
  k9 J! c5 e6 w3 @/ d& z       
        return 0
; D9 x. W4 N, Y6 D4 x       
def two_phase(i):

- o4 U* f# _, I; H, B) K        return 0
9 ?4 p. I! F- B3 ~/ o/ i' ^/ ^. e
! J& |6 w$ s/ s" f) \0 Odef three_phase(bus):
        options=(1,1,1,0,0,0,0)
, b* r5 G7 K9 k4 m        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)
        psspy.scmu(3,options,rval,cktarg)
       
1 k' T! Q: Y1 u) Q9 a: [        return 0
: z% b+ X; G5 t8 A
def line_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
) ?+ [. h: b4 a) u; ^* W/ ~        rval=(0.0,0.0,0.0,0.0,location)
! J4 b- u" r- u' W6 Y: z9 m        cktarg=ckt
# G; Z1 c9 }1 \: C# @$ N$ Z        psspy.scmu(1,options,rval,cktarg)
0 B: g! ?+ N6 x        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
$ i" @: I- e" ?: m- U' z) I8 p        psspy.scmu(3,options,rval,cktarg)
# v& R# S) G% R2 b* T       
3 S* L' V' G8 \8 G& m        return 0

def line_LL_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
/ _) Q' y! n7 D- e8 `' U! t* z" X        options=(1,1,1,0,0,0,0)
- u. y2 p( N2 `& y: S( Q        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
+ Y# t3 l' X' x8 Q) d# F        options=(9,bus_i,bus_j,2,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
: N6 t! b. {  g! e; }- M) l' V        psspy.scmu(3,options,rval,cktarg)
$ f! q' m8 ^# x+ e/ m       
3 [! b  K% I9 {  {6 P. ^        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)
        rval=(0.0,0.0,0.0,0.0,location)
# y2 f& v& C9 k; r% K        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
$ r$ s% N) h) N7 |5 L8 Y        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
) k( r' z! q, |) D1 @; H4 J       
% K- E% `! U3 w. l/ }4 J& Q8 {        return 0

Henryzhang

咋一个三相短路写这么复杂呢？
9 a0 u" _1 \! x( Y3 l/ R
高手详细对比一下python和IPLAN吧
各有哪些优势

meng823222

恩
高手说说两者各自的优缺点吧

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
6 s7 k! H5 x" h9 K4 `不知为什么
2 G6 ?" e/ |  t% B看起来官方的dist_branch_fault、dist_bus_fault等更简洁
但没试过

Henryzhang

原先用IPLAN写branch trip：
1 R+ Q+ v8 z- x
2 ~" |8 i2 Y; }) }3 Y) c3 @! W8 W0 c                PUSH        'ALTR'
                PUSH        '0'
                PUSH        '1'        !CHANGE NETWORKDATA
9 s% Q2 s4 y6 ]1 ]                PUSH        '0'
$ C# M, a' a. |& Z8 q2 ^                PUSH        '3'        !BRANCH DATA
0 W2 E6 ?2 w# T% `                PUSH        '151,201,1'
                PUSH        '1'
                PUSH        '0,'        !ISOLATE THIS BRANCH
                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
+ w( D+ I) H- e5 \! @* ^4 z( j9 P, W                PUSH        ' '                !NO CHANGE IN METERED END
                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
                PUSH  '0'   !EXIT ALTR ACTIVITY

现在python：
dist_branch_trip(151,201,'1')
一句就可以了
9 s  g; a9 g0 a' o: z' ^! f
6 e9 w# i1 c% H! K$ O确实比较简洁

loadflow

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

Henryzhang

用官方提供的应该可靠性高点
就是不容易出错
. J# S3 H' z- Q
$ D9 g2 O) _) ^要想玩PSSE，还是基于最基本命令单元，呵呵

向前向前

直接看不懂呀