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

xuning0076

import globedef
) j" J  ~$ v. g, q  y; ]1 rimport fault
& j- m5 @3 j, P* M& j; t, Bpsspy.case('c:\sav.sav')
psspy.close_report()  # close REPORT windows
! V- n  a; }6 A# r# d$ C% ~! Loptions=(0,0,0,0,0,0,0)
9 S7 S- A4 `8 F+ n) N. i2 }5 A( s$ Spsspy.progress_output(6, "", options)
psspy.report_output(6, "", options)
; ?' \1 v5 c6 J' X' Ufault.three_phase(33600)
4 t. y7 k# M0 c' ^1 D& \: ^
[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
import globedef
$ @' B# e! ^2 z  {import psspy
5 u& d4 M% Q0 P7 G: |( k6 G& E
#a=(1,1,1)
#psspy.seqd(a)

& k3 B. [2 t' J4 y" E
def one_phase_ground(bus,r,x):
) \5 W0 e  D. A  j3 x        options=(1,1,1,0,0,0,0)
5 i% E1 r! s: e, M9 b7 h* B* @        rval=(r,x,0.0,0.0,0.0)
        cktarg=''
/ r- M' l; N# C  T8 Q3 }) q        psspy.scmu(1,options,rval,cktarg)
  {# J. M9 j7 K0 N7 A# _, H        options=(1,bus,1,0,0,0,0)
& R# e- g+ r4 w& u9 P        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
% m6 v, W0 x3 _- P        return 0
9 g, `) R1 W6 n6 }, o: s  e( T       
def two_phase(i):

        return 0
- ^2 Y( g/ g: q
- J. n8 U+ `# A0 L( ^2 S  k( X) }def three_phase(bus):
, A; u8 T2 F+ d% o        options=(1,1,1,0,0,0,0)
) a5 M, j( H4 m/ V6 \. G( D2 ~: p5 j        rval=(0.0,0.0,0.0,0.0,0.0)
/ J% a$ o% j; ]2 h$ m- m7 S! y        cktarg=''
3 S/ X4 ~# j/ a. c# G; Q        psspy.scmu(1,options,rval,cktarg)
3 a7 Y; ~0 O4 F) w9 y0 `        options=(7,bus,1,0,0,0,0)
1 k  y7 B9 V7 j$ J3 E$ o. y        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
        return 0

5 A+ M$ H! c" R: ]+ xdef line_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
& U7 l* ~) e+ Q1 |* k$ K: F        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,location)
$ O2 }! X, g0 l+ h5 t! Z+ @        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)
       
        return 0
" E# J7 y2 ^- K% _/ n) \
def line_LL_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
0 c- g# K' d% _3 {( F        options=(1,1,1,0,0,0,0)
- ~# c' F3 E7 \# p! Y        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
0 h2 w, ~8 ?7 |6 D$ W2 r  `8 p        psspy.scmu(1,options,rval,cktarg)
3 J2 W; ~( o  R& W        options=(9,bus_i,bus_j,2,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
) M7 B* Z9 t) n! L( p        return 0
- h/ c" D" g' y8 T' J! L0 J
def get_breaker_cline_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
. u- y$ V2 Z7 J+ ~5 s4 I* \7 N        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
8 o: T! i, j/ K        psspy.scmu(1,options,rval,cktarg)
2 o& E, E  @. }: A! Y8 [2 O        options=(9,bus_i,bus_j,1,1,0,0)
7 i; i$ Y8 S8 F5 B1 _3 M5 p1 H        psspy.scmu(2,options,rval,cktarg)
3 O( X- p  z  F$ x; Y        psspy.scmu(3,options,rval,cktarg)
. g! L& G7 v$ x$ ]7 A0 C       
! l0 q3 H# _$ u( {        return 0

Henryzhang

咋一个三相短路写这么复杂呢？
$ c! E) H( r. ~$ L; u8 p) x
高手详细对比一下python和IPLAN吧
各有哪些优势

meng823222

恩
# x! y$ C# P, z高手说说两者各自的优缺点吧

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
不知为什么
看起来官方的dist_branch_fault、dist_bus_fault等更简洁
8 o, L8 _+ E  g# [0 |' w$ {但没试过

Henryzhang

原先用IPLAN写branch trip：
  _( G; ~2 R) [' B
                PUSH        'ALTR'
; W2 U: O8 X# {0 [: ?, |5 D                PUSH        '0'
                PUSH        '1'        !CHANGE NETWORKDATA
  C$ s/ G% M& W                PUSH        '0'
& P  f8 W; ]# t                PUSH        '3'        !BRANCH DATA
                PUSH        '151,201,1'
2 p: w) f! I9 p                PUSH        '1'
                PUSH        '0,'        !ISOLATE THIS BRANCH
8 }! m4 ^! N" g6 B                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
. a4 M0 [2 s! Y, H1 m                PUSH        ' '                !NO CHANGE IN METERED END
                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
+ d, C, X! x+ B& H& T& Y, L4 b                PUSH  '0'   !EXIT ALTR ACTIVITY

: y' \5 z8 ^0 Y: j现在python：
  H% |! b0 z- d5 b! cdist_branch_trip(151,201,'1')
一句就可以了

2 o7 w8 f% {7 `- t& |' `确实比较简洁

loadflow

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

Henryzhang

用官方提供的应该可靠性高点
- i0 I, J) O1 ]" _& i就是不容易出错

要想玩PSSE，还是基于最基本命令单元，呵呵

向前向前

直接看不懂呀