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

xuning0076

import globedef
. {! H* v! F, ^# E( D& vimport fault
, L( Y; k8 n' |* X1 X! q/ n8 \8 ^psspy.case('c:\sav.sav')
psspy.close_report()  # close REPORT windows
options=(0,0,0,0,0,0,0)
. b, n4 e/ g# L( epsspy.progress_output(6, "", options)
: v& P) n3 C. _" Kpsspy.report_output(6, "", options)
fault.three_phase(33600)
* r0 Z& j: G. v- |
[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
import globedef
  c* u, }9 V& C7 x- gimport psspy

#a=(1,1,1)
#psspy.seqd(a)


- H/ [5 g/ k3 Pdef one_phase_ground(bus,r,x):
        options=(1,1,1,0,0,0,0)
$ `9 i! X* T5 m( C        rval=(r,x,0.0,0.0,0.0)
        cktarg=''
! O& z2 B- e/ u  h, c: T        psspy.scmu(1,options,rval,cktarg)
        options=(1,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
1 n  ~# v- f8 M! \# Z4 p8 L  H0 \       
' |" J. @) k. [# }. F        return 0
: \- V  ]. |' _, d  a       
def two_phase(i):
* U% |! |1 q3 w
4 R) Q; F2 v; q7 U        return 0
2 F: m0 Z" E- h3 [# q, y+ h+ o8 g
3 v, x2 M2 s3 o8 f% r) O1 U/ gdef three_phase(bus):
        options=(1,1,1,0,0,0,0)
" l; V- D0 x; P# @- f5 _' s        rval=(0.0,0.0,0.0,0.0,0.0)
! c8 J; D! A6 P) M        cktarg=''
  T4 z5 x' a2 o) U0 l        psspy.scmu(1,options,rval,cktarg)
        options=(7,bus,1,0,0,0,0)
" r( }  A7 ^# b* z        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
        return 0
* l: b2 M( X- |5 u: ~
- d1 _4 p6 @8 D$ \0 fdef line_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)
        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,1,1,0,0)
/ M; Z  c9 B" v( m, D& i        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
! Q0 ~' H0 a9 t. t% n/ M' S       
        return 0

0 b) B9 I7 i" c$ t  i) x% ~2 Qdef line_LL_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
/ L- s& w" m! \2 A: _% h        rval=(0.0,0.0,0.0,0.0,location)
5 n% Y4 B+ A' `5 `9 g! c        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,2,1,0,0)
' {! L; ]$ k+ s( l( J        psspy.scmu(2,options,rval,cktarg)
2 V6 R0 E; p& P6 a        psspy.scmu(3,options,rval,cktarg)
       
        return 0
2 C! d1 x# B. |6 @: `
- ]# M5 N  n) a; Jdef get_breaker_cline_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
1 F: K- m5 O% [! T        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
9 B! W6 W: v/ D: S, s        psspy.scmu(1,options,rval,cktarg)
0 m3 w5 e* W) N  J, p8 q. ^        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
7 U" A- }1 o  N4 |- g9 @7 R        psspy.scmu(3,options,rval,cktarg)
       
5 b& K8 y& j" V        return 0

Henryzhang

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

, a: p1 Q/ ]% ]2 K, j* W2 E+ d高手详细对比一下python和IPLAN吧
各有哪些优势

meng823222

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

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
2 D8 ~) r' b3 A4 d+ N* J不知为什么
2 K$ M! {" \" |4 r; {看起来官方的dist_branch_fault、dist_bus_fault等更简洁
但没试过

Henryzhang

原先用IPLAN写branch trip：

                PUSH        'ALTR'
" C& ^$ [0 ~8 d: {8 Y0 Q- P                PUSH        '0'
. H2 {* t/ n. A) @3 Y: a  y6 s                PUSH        '1'        !CHANGE NETWORKDATA
                PUSH        '0'
& G  V6 u1 S. Y7 K- d2 I* O/ z9 x0 H                PUSH        '3'        !BRANCH DATA
                PUSH        '151,201,1'
: Q& \) A  R- R9 q5 [" W                PUSH        '1'
. N% x) Q, }& d: F                PUSH        '0,'        !ISOLATE THIS BRANCH
                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
                PUSH        ' '                !NO CHANGE IN METERED END
5 S  L4 E" D- J# }4 w2 o. b                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
                PUSH  '0'   !EXIT ALTR ACTIVITY
& d$ t$ Y. U9 R: B7 t4 D
! `. C4 x$ y( F8 X% L) Q现在python：
' {' {) b0 @8 p6 D) odist_branch_trip(151,201,'1')
一句就可以了

确实比较简洁

loadflow

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

Henryzhang

用官方提供的应该可靠性高点
就是不容易出错
4 Q$ c) f0 v# z: [" Q  R& p) B% `
要想玩PSSE，还是基于最基本命令单元，呵呵

向前向前

直接看不懂呀