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

xuning0076

import globedef
import fault
4 m  Q4 ^* H* ~4 e7 e3 Qpsspy.case('c:\sav.sav')
psspy.close_report()  # close REPORT windows
9 P7 x5 Z! e2 h" C# A1 {; Joptions=(0,0,0,0,0,0,0)
psspy.progress_output(6, "", options)
! I4 h( y! m! `1 E* M7 L8 cpsspy.report_output(6, "", options)
, L/ E5 ~0 n* |" l. B1 pfault.three_phase(33600)
3 O1 f5 v# r+ ?. D5 ?, X
; W! z3 p% R; g# w[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
import globedef
import psspy

0 J( w- t" M& y$ }/ O#a=(1,1,1)
9 K  g% ~* Q2 D6 ?$ g+ k4 g1 U# S#psspy.seqd(a)

# M( {7 x8 D2 h1 k5 v6 G
) T# X0 ^6 f( G" j/ odef one_phase_ground(bus,r,x):
9 w7 F' g. L' l4 ]* P        options=(1,1,1,0,0,0,0)
( N4 z9 l! A+ \9 A' P) b        rval=(r,x,0.0,0.0,0.0)
4 c8 c% ^( ]" D. [' M1 K        cktarg=''
        psspy.scmu(1,options,rval,cktarg)
5 J. [, y" |# Z5 k        options=(1,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
9 g+ O  `8 Y/ a0 N+ c, F( f        psspy.scmu(3,options,rval,cktarg)
! B6 G  i9 b5 s7 U4 _  I9 k! n0 G       
  y' A' a. M( E4 b/ d' O3 c! v( s( v        return 0
       
% ]. j- u3 T' ~def two_phase(i):

- q  ?/ j; o; T1 F0 D" I        return 0
% J: E2 o1 `& B$ ]& b  q) e
def three_phase(bus):
        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,0.0)
        cktarg=''
4 }$ S! A. T0 p4 Z        psspy.scmu(1,options,rval,cktarg)
        options=(7,bus,1,0,0,0,0)
" b* O. R$ x! n" t1 J3 n" O        psspy.scmu(2,options,rval,cktarg)
2 e. n' P0 b9 h: @2 ^        psspy.scmu(3,options,rval,cktarg)
       
        return 0

( C5 X7 x7 c2 Z! I+ n- T2 gdef line_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
, V- F& O3 @; Y7 c7 d/ s9 x        options=(1,1,1,0,0,0,0)
6 _$ g" {* d" d        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
3 K! J6 Z  O: h. @* k) d1 T" o2 F6 R        psspy.scmu(1,options,rval,cktarg)
" X/ y. z7 x( C        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
, {- @, `% g! k        psspy.scmu(3,options,rval,cktarg)
' |  l( Y, l6 T0 j       
& Q: V+ z9 C- k& O$ B7 z( R% H. o        return 0
4 E' h; Z6 T. ?1 _7 W
def line_LL_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)
        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,2,1,0,0)
* V( M& I8 [9 A) r7 x        psspy.scmu(2,options,rval,cktarg)
5 f: D4 |0 L2 c$ i. z! d        psspy.scmu(3,options,rval,cktarg)
+ h8 u. j2 N$ m3 ?       
        return 0
! o; N6 h  z) {1 i
def get_breaker_cline_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
3 L! f1 Q) V5 p# l. [        options=(1,1,1,0,0,0,0)
* q6 K; K) q( A- M        rval=(0.0,0.0,0.0,0.0,location)
; g; J, c4 ^" M* Q8 m# O2 \        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
: Z0 C; c1 n7 Y- ^        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
( }' Y* J2 [+ U  W2 v9 W        psspy.scmu(3,options,rval,cktarg)
       
        return 0

Henryzhang

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

4 P* ]. I4 C, H( M9 N  A3 b5 p4 B高手详细对比一下python和IPLAN吧
/ u0 Q# t9 ~: d各有哪些优势

meng823222

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

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
& ~! n* f) o; G' w% _! m. R5 X不知为什么
  ~3 k+ J2 b' T8 L( o1 G& e看起来官方的dist_branch_fault、dist_bus_fault等更简洁
: U. w- |+ Y9 p1 R  B1 o! w但没试过

Henryzhang

原先用IPLAN写branch trip：
" r5 q! x9 S$ c- f5 e, U6 ]  B
                PUSH        'ALTR'
                PUSH        '0'
                PUSH        '1'        !CHANGE NETWORKDATA
- w# u6 L0 E( g" Q' A                PUSH        '0'
% e! [& h! ?0 _( L  Y                PUSH        '3'        !BRANCH DATA
                PUSH        '151,201,1'
% [. [' G- x8 O& w                PUSH        '1'
" Y. a2 X$ a" u2 N# L2 }                PUSH        '0,'        !ISOLATE THIS BRANCH
4 i* m4 O) G1 Y                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
5 \% ^2 V9 T. D( E+ t1 Q                PUSH        ' '                !NO CHANGE IN METERED END
0 L! y. L' R2 c1 N                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
                PUSH  '0'   !EXIT ALTR ACTIVITY
1 ~8 L2 h8 v- }+ l- k( |: t* s
$ X7 z  b6 B) ~$ w" c: u% g; C现在python：
1 F4 I3 w, m7 u: P/ V/ ~dist_branch_trip(151,201,'1')
8 C% G! Q7 E! O* O8 q& x: E3 W- Q; |一句就可以了
" C$ P5 C$ N2 ]5 L6 f& A; p
8 X$ `& g( ~4 _. W确实比较简洁

loadflow

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

Henryzhang

用官方提供的应该可靠性高点
0 U, r# k3 V0 H5 A$ E就是不容易出错

- a3 z$ H: J5 b9 s要想玩PSSE，还是基于最基本命令单元，呵呵

向前向前

直接看不懂呀