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

xuning0076

import globedef
- @7 A. s9 g3 V! }; n' ?- f$ timport fault
2 R3 Z6 x- r) o( p6 U/ E% S# z: Ypsspy.case('c:\sav.sav')
psspy.close_report()  # close REPORT windows
options=(0,0,0,0,0,0,0)
psspy.progress_output(6, "", options)
psspy.report_output(6, "", options)
fault.three_phase(33600)
% m; B1 Y! x. q$ T* t
[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
import globedef
/ J/ U, H5 \% _" B7 _import psspy

#a=(1,1,1)
8 ~. F! Z0 V+ _#psspy.seqd(a)
# g2 P, m  `3 ]5 ^& _7 q4 H5 _3 S
7 Y( j5 i1 j2 @
def one_phase_ground(bus,r,x):
        options=(1,1,1,0,0,0,0)
. S/ J- u9 P. w( g8 C& ^        rval=(r,x,0.0,0.0,0.0)
        cktarg=''
* q8 A' p0 r: Z2 n        psspy.scmu(1,options,rval,cktarg)
        options=(1,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
% c# L0 X- B: l! t4 m8 c        psspy.scmu(3,options,rval,cktarg)
       
. r9 _( X+ n! l0 L6 H9 E        return 0
- Q" x% X( ?0 b2 s' j, w( z       
6 x1 [) J# i. [3 A: w" G& V! Xdef two_phase(i):
" ]7 d+ s0 L/ B
        return 0
+ K/ h. t' [7 {0 Y, X
9 L* E9 }% \5 I2 [: X8 V6 Xdef three_phase(bus):
        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,0.0)
6 \. |; d# P; z0 B# w        cktarg=''
        psspy.scmu(1,options,rval,cktarg)
        options=(7,bus,1,0,0,0,0)
( @; z$ d. Z# S  V        psspy.scmu(2,options,rval,cktarg)
- A  r; S# x' |0 A( ^& \( w2 P. _        psspy.scmu(3,options,rval,cktarg)
5 h6 a# T+ B. E; ?  A9 G- @       
0 g& V& j. t$ N' r6 L9 N        return 0
& m5 d5 Z! K9 _* }
def line_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
3 K( p# g5 M1 }2 \) ^        rval=(0.0,0.0,0.0,0.0,location)
3 M% I. S1 A5 _5 o        cktarg=ckt
& e/ U5 i: S+ Q( g        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
' }8 O7 e* N: F  `& z# x        psspy.scmu(3,options,rval,cktarg)
       
        return 0

def line_LL_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
, S% _, W3 c7 u  F2 k( C1 \  e        options=(1,1,1,0,0,0,0)
0 [! p' j" [2 g! ]        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
5 ?- K0 w& Q! {& N+ ~        psspy.scmu(1,options,rval,cktarg)
/ S) T# H8 r3 ^- b  k        options=(9,bus_i,bus_j,2,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
3 s/ ~" h0 v9 A% l  P        psspy.scmu(3,options,rval,cktarg)
. g) C) V- Y' a- s/ Q       
        return 0
4 i5 ^9 Z3 [' @5 C7 l! i: j; V0 Q# |
def get_breaker_cline_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
2 s8 {) k! V- I7 c0 Y        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,location)
        cktarg=ckt
6 l5 T2 {/ t# U" ~        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)
) u7 y  L! _" V+ v7 B, h       
) d+ V% Q- {& p! y0 z  I; ]% E        return 0

Henryzhang

咋一个三相短路写这么复杂呢？
9 x. j! M6 d3 |& j$ m$ K6 r3 h
( S5 F) v, O  l9 e8 Y2 \! ?高手详细对比一下python和IPLAN吧
# G  N9 n7 U( p& w0 F8 T: y各有哪些优势

meng823222

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

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
不知为什么
5 m' Q- H! Z" k: D1 G  p看起来官方的dist_branch_fault、dist_bus_fault等更简洁
但没试过

Henryzhang

原先用IPLAN写branch trip：

, c" B  i: o" E* d( S                PUSH        'ALTR'
$ a" v) u4 W# a" F% X) S                PUSH        '0'
; D6 x. l: W$ A( _9 T) F$ v1 |                PUSH        '1'        !CHANGE NETWORKDATA
0 k1 b! Y- P! m2 X6 t( e                PUSH        '0'
& F! [1 k! F) |9 ^' v, W                PUSH        '3'        !BRANCH DATA
: @, d- }/ M# S" g2 L, V                PUSH        '151,201,1'
! b4 y1 P; G/ Y; X' W/ A                PUSH        '1'
                PUSH        '0,'        !ISOLATE THIS BRANCH
                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
% s  E8 x+ D3 K- g6 q: Y0 e1 n                PUSH        ' '                !NO CHANGE IN METERED END
$ t! N& G% p4 i) a                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
                PUSH  '0'   !EXIT ALTR ACTIVITY
) s$ M) a+ y: Q. R) g& }
现在python：
dist_branch_trip(151,201,'1')
一句就可以了

5 u0 G! Z+ r1 |! K1 M确实比较简洁

loadflow

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

Henryzhang

用官方提供的应该可靠性高点
% e2 ^9 L# J/ y5 f就是不容易出错
2 ^0 @& {. N3 @* Q. v& R
2 @# a) N$ \8 ^& y( Q- x$ b要想玩PSSE，还是基于最基本命令单元，呵呵

向前向前

直接看不懂呀