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

xuning0076

import globedef
5 X( F6 T; M' ximport fault
psspy.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)
# U, a  I3 x- n: [fault.three_phase(33600)

[ 本帖最后由 dsyljh 于 2008-8-8 08:40 编辑 ]-%-96529-%-

xuning0076

这是FAULT.PY的内容
, w4 L4 D) V2 J2 n) F+ himport globedef
import psspy
% ~$ \  C* m# f- O4 j% H
; ?# A& ~4 e0 d. e: P3 y#a=(1,1,1)
% o0 V  @4 H, l7 b, n. K; \#psspy.seqd(a)

0 R( Q  N$ ~6 ~- L
- O+ I% z4 z7 c: A. p8 W& zdef one_phase_ground(bus,r,x):
        options=(1,1,1,0,0,0,0)
! i  _8 K7 g0 U        rval=(r,x,0.0,0.0,0.0)
. o9 \$ F4 n7 n9 S        cktarg=''
1 x* [! U" X# g, [        psspy.scmu(1,options,rval,cktarg)
        options=(1,bus,1,0,0,0,0)
        psspy.scmu(2,options,rval,cktarg)
& b0 M" p6 m# M0 O        psspy.scmu(3,options,rval,cktarg)
       
        return 0
% z/ Q9 j( s# ^2 k: S. J       
def two_phase(i):
7 z% f( f4 \% y6 l% V# B
        return 0

def three_phase(bus):
  l; w/ G) c  _9 v. R        options=(1,1,1,0,0,0,0)
! h# {* w% o6 D- R1 P8 y+ p+ n        rval=(0.0,0.0,0.0,0.0,0.0)
! A. g2 n  c2 s& G2 Q        cktarg=''
        psspy.scmu(1,options,rval,cktarg)
& j9 |+ g( J  S8 }: |        options=(7,bus,1,0,0,0,0)
* b! u, g( C: g$ G- }  @        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
3 p; O2 A. s  k8 B& U7 @# _        return 0
! W: c1 I$ q* @. }
- ]# w. n* ^# R( e1 hdef line_LG_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
4 r4 T" _: |- x        options=(1,1,1,0,0,0,0)
        rval=(0.0,0.0,0.0,0.0,location)
, I# M) {6 G3 j5 V( e6 v        cktarg=ckt
        psspy.scmu(1,options,rval,cktarg)
& x0 j1 m0 O1 l        options=(9,bus_i,bus_j,1,1,0,0)
0 @# L0 J$ q5 j7 a( V3 M        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
1 g# V: d; E0 |  \) ~. m       
( B+ U+ L% q" n# T        return 0
  s! m( O- m: v7 e/ t
def line_LL_fault_both_ends_close(bus_i,bus_j,ckt,location,r,x):
        options=(1,1,1,0,0,0,0)
' |/ l1 }- [$ e, L% m% N% |6 d# ^        rval=(0.0,0.0,0.0,0.0,location)
" J0 f( t6 r( |- d; g6 m7 J# b( a, {        cktarg=ckt
$ `1 v1 H! R/ ~& W6 a        psspy.scmu(1,options,rval,cktarg)
1 O6 ]0 a: g1 T+ t4 R- f        options=(9,bus_i,bus_j,2,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
        psspy.scmu(3,options,rval,cktarg)
       
7 V/ N3 O" F) W4 ?1 T7 D        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)
& k: p( \& e& r        cktarg=ckt
+ R' T, ]& L+ I" U$ q6 N4 p/ G        psspy.scmu(1,options,rval,cktarg)
        options=(9,bus_i,bus_j,1,1,0,0)
        psspy.scmu(2,options,rval,cktarg)
- @" F8 P! A* ^$ i        psspy.scmu(3,options,rval,cktarg)
' y$ [% {  a# Z# x       
        return 0

Henryzhang

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

/ R* q" S1 I5 m  b( p1 I' w高手详细对比一下python和IPLAN吧
各有哪些优势

meng823222

恩
* |5 L8 I3 D: E6 S0 }高手说说两者各自的优缺点吧

Henryzhang

楼主上述故障实现没有直接用官方提供的disturbance APIs
& d2 L* o: A+ c. h不知为什么
& W% E4 x  A) m  Q看起来官方的dist_branch_fault、dist_bus_fault等更简洁
但没试过

Henryzhang

原先用IPLAN写branch trip：

                PUSH        'ALTR'
                PUSH        '0'
# i; C0 B7 W/ G                PUSH        '1'        !CHANGE NETWORKDATA
2 k* C! x9 C5 }& y0 T                PUSH        '0'
                PUSH        '3'        !BRANCH DATA
7 F+ Q1 B) u/ _5 @/ c                PUSH        '151,201,1'
                PUSH        '1'
6 v3 n/ T3 c+ e7 q" j- Z                PUSH        '0,'        !ISOLATE THIS BRANCH
                PUSH        '0'          !NO MORE CHANGES TO THIS BRANCH
                PUSH        ' '                !NO CHANGE IN METERED END
6 P: N1 n* S+ p5 E6 H; G                PUSH        ' '          !NO MORE BRANCHES TO BE CHANGED
                PUSH  '0'   !EXIT ALTR ACTIVITY

/ i5 O& f- C! S; l现在python：
dist_branch_trip(151,201,'1')
" T' x/ L1 c0 I/ ?- m/ q% x: |一句就可以了
$ U4 R, B; b! \# Y
确实比较简洁

loadflow

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

Henryzhang

用官方提供的应该可靠性高点
就是不容易出错

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

向前向前

直接看不懂呀