求助四机两区域动态仿真

actfate

潮流结果和PSASP的完全对得上，但在不设故障情况下在0.17s失稳，G1的Pelec如图


lofl
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
) D) J4 y% E$ ?" i. Q1 @ Executing activity read,file
read,file
0 L4 v  K' h+ l: u+ x ENTER INPUT FILE NAME (0 TO EXIT, 1 FOR TERMINAL): 4g2a
ENTER IC, SBASE
9 V4 i' g' m8 S ENTER TWO LINE HEADING
' D) _, V! Q; r. Z ENTER BUS DATA
1 y: x+ _" x) ~- y. B4 P ENTER LOAD DATA
0 s, I" e/ T, u3 C2 Y# {! S5 n ENTER GENERATOR DATA
ENTER NON-TRANSFORMER BRANCH DATA
0 I  r- g  O8 a5 |) {3 p& i* ] ENTER TRANSFORMER DATA
ENTER AREA INTERCHANGE DATA
ENTER TWO-TERMINAL DC LINE DATA
0 ?/ `* O) T3 i$ J ENTER VSC DC LINE DATA
ENTER SWITCHED SHUNT DATA
$ A7 V9 ]8 C/ X3 ? ENTER TRANSFORMER IMPEDANCE CORRECTION DATA
  P7 R( O4 g/ l$ a' v) d. { ENTER MULTI-TERMINAL DC LINE DATA
ENTER MULTI-SECTION LINE DATA
( _( R/ z) F- } ENTER ZONE NAME DATA
' }" X$ F2 W; `( Z2 I" t- X ENTER INTER-AREA TRANSFER DATA
ENTER OWNER NAME DATA
ENTER FACTS DEVICE DATA
  Q( g4 Q' i# Q+ |: ^ ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
7 Z4 Q4 T: i7 @6 G1 f Executing activity fnsl
4 d; d4 z6 r4 N4 j2 Y! Yfnsl
* D' Y! _' B% U( M) X; J( K ORDERING NETWORK
DIAGONALS =     10  OFF-DIAGONALS =      9  MAX SIZE =     18
% u- `: n4 l9 x/ J/ U ENTER ITERATION NUMBER FOR VAR LIMITS
  0 FOR IMMEDIATELY, -1 TO IGNORE COMPLETELY: 99
  ITER       DELTAP      BUS         DELTAQ      BUS        DELTA/V/      BUS       DELTAANG      BUS
   0        21.2968(      401   )    2.7013(      801   )    0.12761(      801   )   0.58266(      801   )
   1         4.2106(      801   )    8.1793(      901   )    0.12606(      701   )   0.04352(        1   )
- Y8 f; D# i8 g2 n  |% e   2         0.2518(      901   )    0.7799(      901   )    0.02045(      701   )   0.00508(        1   )
   3         0.0122(      901   )    0.0124(      901   )    0.00048(      701   )   0.00029(      801   )
   4         0.0000(      801   )    0.0000(      901   )    0.00000(      801   )   0.00000(      801   )
   5         0.0000(      801   )    0.0000(      801   )
& ^* s# c" T' _; i& U: z& n7 n' G0 s REACHED TOLERANCE IN   5 ITERATIONS
LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    801 [9           230.00]
9 h" |0 p6 F+ Y! C0 d6 T% @ SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
SWING BUS SUMMARY:
( x; Z$ W% t6 g* B   BUS# X-- NAME --X BASKV      PGEN     PMAX    PMIN      QGEN     QMAX    QMIN
) O, G2 A# N4 ~- A    301 4           20.000     700.1   9999.0 -9999.0     202.3   9999.0 -9999.0
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
$ U. \4 i% g' g7 w3 q" ^# y7 e; l Executing activity save,4g2a00
save,4g2a00
CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:06
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
5 }+ F0 C" R  M8 r Executing activity cong
+ i+ c4 T) ~6 V# |$ n# w+ t! X$ |cong
GENERATOR CONVERSION COMPLETED
- ~# W/ G; ?# L$ k3 j$ U& V5 z% l ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
9 _& P: R/ x" U# J. f9 R$ b Executing activity conl all
conl all
7 e' R) g& p0 P% ~+ j4 u4 h ENTER % CONSTANT I, % CONSTANT G FOR ACTIVE POWER:
6 n9 _( R# Y* a6 _! T# c ENTER % CONSTANT I, % CONSTANT B FOR REACTIVE POWER:
LOAD TO BE REPRESENTED AS:
4 K+ w0 Y. j; a! {, D. ?0 g    ACTIVE  REACTIVE
   100.00%   100.00% CONSTANT POWER
     0.00%     0.00% CONSTANT CURRENT
2 j9 ~. R6 q3 j9 ^" u     0.00%     0.00% CONSTANT ADMITTANCE
ENTER 1 IF O.K., 0 OTHERWISE: 1
      2 OF      2 LOADS CONVERTED
. f" m; m6 f; G/ J/ H+ v$ W ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
( w+ O4 }' i2 T, A2 N. B Executing activity fact
fact
; D" u) R- K+ s. N* N ORDERING NETWORK
DIAGONALS =     11  OFF-DIAGONALS =     10  MAX SIZE =     20
9 W, s* ]8 Z( K3 \% H     11 DIAGONAL AND     10 OFF-DIAGONAL ELEMENTS
4 y. T  I" @$ W; X7 l, [ ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity tysl
tysl
# |' k( J  N, b6 @ ITER DELTAV/TOL X-------- AT BUS ---------X  REAL(DELTAV)  IMAG(DELTAV)
) G( z% n* d- z8 |    1     0.003       1 [1           20.000]    0.3029E-07   -0.1089E-07
REACHED TOLERANCE IN   1 ITERATIONS
LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    601 [7           230.00]
) {  F- p" t; W5 K SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
( O  m- R1 ]  A) I' T ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
6 L  K( e# t' Q' N$ t! b' q" D Executing activity save,4g2a00
  Y- c% b, k" d2 s2 o9 Lsave,4g2a00
. r3 W+ x( p  ?% j CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:07
# ?1 S) u! K: o( h& U; [6 w# Q7 Q& C ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity RTRN
6 W) X, _, Q+ r. M+ gRTRN
; _/ b$ `7 d4 h' o: h& w, L5 N ACTIVITY?
Executing activity DYRE
OUT OF FILE DATA--SWITCH TO TERMINAL INPUT MODE
GENERATOR MODELS USE:
. y, c* O0 Q9 k5 e' A+ e4 l2 _  CONS        1-    56
( h, I1 j% q# L4 J& ^$ k  STATES      1-    24

: d" _' _) x8 f& h, b- Q( R$ ^ SUMMARY OF MODELS READ:
GENS:  GENROU
: k! x' ]) b4 h; m: ~9 s6 g           4
NEXT AVAILABLE ADDRESSES ARE:
: {% G% q1 D  c( S, P8 P    CON  STATE    VAR   ICON
9 p3 [4 a3 N* @( Y     57     25      1      1
NO MODEL CALLS IN CONNECTION ROUTINES--DYNAMICS SKELETON MAY BE USED
( P: Y& b0 R' M* r6 C, G% ~: ` Executing activity CHAN

      3 OUTPUT CHANNELS HAVE BEEN ADDED
NEXT AVAILABLE ADDRESSES ARE:
1 _8 }/ {0 y+ o4 O: T; F   CHANNEL   VAR   ICON
        4      2      2
Executing activity STRT
2 k5 U" U) f- t4 m GENROU AT BUS      1 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
, n0 r  _- a5 z/ E8 }  p. ^ GENROU AT BUS    101 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
GENROU AT BUS    201 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
4 T6 g0 E$ V7 C7 i$ X1 [# \+ t3 N GENROU AT BUS    301 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
6 x$ X' ^: p0 A; n7 l4 q+ j6 I     PTI INTERACTIVE POWER SYSTEM SIMULATOR--PSS/E      SAT, JUN 02 2012  15:08
1
4 S/ n: B+ H7 D 2
INITIAL CONDITION LOAD FLOW USED  1 ITERATIONS
----------------------------- MACHINE INITIAL CONDITIONS -----------------------------
   BUS# X-- NAME --X BASKV ID  ETERM   EFD    POWER    VARS   P.F.  ANGLE   ID     IQ
: o9 `# C  e9 k3 e      1 1           20.000 1  1.0300 2.4652  700.00  185.25 0.9667  70.87 0.7104 0.3249
    101 2           20.000 1  1.0100 2.5768  700.00  234.84 0.9481  59.45 0.7504 0.3109
' C% A9 B; j! ?) v/ _7 W2 I    201 3           20.000 1  1.0300 2.4910  719.00  176.22 0.9713  44.98 0.7270 0.3303
" z$ b* L* y3 l* k    301 4           20.000 1  1.0100 2.5189  700.10  202.33 0.9607  33.50 0.7359 0.3181
4 N+ M3 ^$ C; ~/ ?' L INITIAL CONDITIONS CHECK O.K.
) x! ~& Q' K! e1 P SNAPSHOT STORED IN FILE D:\PSSE30\EXAMPLE\4g2a.snp AT TIME =  -0.020
/ H. l2 t- B3 n) X1 r, p CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out
NETWORK NOT CONVERGED AT TIME =  0.1700
+ Z$ E* w" f3 f/ ^8 }  26   191.996    801 [9           230.00]     0.1697E-01     0.8974E-02
# S* N' n; K/ _( gDecimal number outside allowable range!
CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out
6 t  O4 V% A3 w- p$ T/ Q

" ^5 A' s4 P  C, U/ h潮流数据预处理应该没问题，dyr文件只包括发电机模型，用PSASP6型模型仿真结果正常，到PSS/E就出问题了
9 Y& l, b5 }$ z8 a( ?& ?8 h" ^
) }9 v( `  m4 K! f% p6 @6 ~8 W最后两项饱和度是自己随便写的，其他数据来自kunder的书

apsal

问题很明显，你在执行STRT初始化操作时提示：
Executing activity STRT
+ F# O6 H# B$ n; qGENROU AT BUS      1 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
GENROU AT BUS    101 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
GENROU AT BUS    201 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
7 f3 q7 c9 S) ]0 D& `GENROU AT BUS    301 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000

$ o- I* {4 v3 E# m0 l7 `% f/ \这表明，你在潮流数据中的ZSORCE（CONG操作中将发电机等效为Norton电流源就是用的这个ZSORCE）与动态仿真中的发电机的次暂态电抗X"不相等。两个不相等，仿真就是错的。
这就是问题的所在。
( U9 D/ y3 d' s; {( d  G看起来你的X"可能是对的，ZSORCE用的的默认值1.0，那么将潮流数据中发电机的ZSORCE的电抗值改为对应X"就可以了。

actfate

厉害！！！光看INITIAL CONDITIONS CHECK O.K.去了
多谢

mingyu

学习一下。

wlm_28

看不懂呀

lxself

强啊强啊强啊强啊强啊

DUTshen

谢谢！谢谢！谢谢！

BPA_CN

同问，哪位大神指教一下？
/ S. P1 ]1 ]6 L- I9 T