求助四机两区域动态仿真

actfate

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

1 `$ z  y" u3 O: w" d8 \+ m0 }' u
lofl
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
7 e2 U* b- s5 X6 I& m1 O/ G Executing activity read,file
read,file
$ }# k: [6 L3 u. G ENTER INPUT FILE NAME (0 TO EXIT, 1 FOR TERMINAL): 4g2a
0 B, R7 d1 c( W7 Z4 [  K ENTER IC, SBASE
ENTER TWO LINE HEADING
ENTER BUS DATA
$ J% b7 ~; c4 `7 a4 I ENTER LOAD DATA
ENTER GENERATOR DATA
4 ^8 L4 s1 Y, X ENTER NON-TRANSFORMER BRANCH DATA
ENTER TRANSFORMER DATA
  k( X! A$ F: _! O ENTER AREA INTERCHANGE DATA
ENTER TWO-TERMINAL DC LINE DATA
ENTER VSC DC LINE DATA
ENTER SWITCHED SHUNT DATA
ENTER TRANSFORMER IMPEDANCE CORRECTION DATA
0 ]" v* v6 r3 t# L ENTER MULTI-TERMINAL DC LINE DATA
ENTER MULTI-SECTION LINE DATA
ENTER ZONE NAME DATA
* ^9 O: x2 B' T, `" ~ ENTER INTER-AREA TRANSFER DATA
ENTER OWNER NAME DATA
5 p* E. Q4 V9 F ENTER FACTS DEVICE DATA
9 H7 ]. p6 V. Y3 T ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fnsl
fnsl
ORDERING NETWORK
DIAGONALS =     10  OFF-DIAGONALS =      9  MAX SIZE =     18
6 t- o: @; j1 n 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 P. |% ]$ ~   0        21.2968(      401   )    2.7013(      801   )    0.12761(      801   )   0.58266(      801   )
; \$ g5 y% _1 e! q! _: \   1         4.2106(      801   )    8.1793(      901   )    0.12606(      701   )   0.04352(        1   )
   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   )
* E7 i& \/ A3 z9 R* K! o8 c( a, Y   4         0.0000(      801   )    0.0000(      901   )    0.00000(      801   )   0.00000(      801   )
9 G% X1 s5 ]6 P. F# y   5         0.0000(      801   )    0.0000(      801   )
REACHED TOLERANCE IN   5 ITERATIONS
LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    801 [9           230.00]
SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
/ F* a! ]/ C; }+ P/ B/ ~5 G SWING BUS SUMMARY:
% n8 D, R# W5 ^, L- e   BUS# X-- NAME --X BASKV      PGEN     PMAX    PMIN      QGEN     QMAX    QMIN
4 x( d% `- k9 W5 o    301 4           20.000     700.1   9999.0 -9999.0     202.3   9999.0 -9999.0
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
: z4 U! S4 L, b, S; J( r. D$ j) T9 t Executing activity save,4g2a00
save,4g2a00
CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:06
4 |; Y7 e' g7 m4 n2 S ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
- ?+ u- g: w/ s( X& c' T Executing activity cong
cong
GENERATOR CONVERSION COMPLETED
% B1 H3 S% ~4 P2 ~6 ^ ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
! H9 S& \- v. b7 Z Executing activity conl all
+ x! V& G% K! K# [conl all
ENTER % CONSTANT I, % CONSTANT G FOR ACTIVE POWER:
ENTER % CONSTANT I, % CONSTANT B FOR REACTIVE POWER:
( M" m$ c6 ?7 z. p LOAD TO BE REPRESENTED AS:
* Z4 C( _0 K0 [7 c    ACTIVE  REACTIVE
! d) A+ a+ f6 y" q! w   100.00%   100.00% CONSTANT POWER
1 H3 D+ }; ?* P3 _( `     0.00%     0.00% CONSTANT CURRENT
0 |& j* ]6 ?+ o8 D! S  Y% E8 `     0.00%     0.00% CONSTANT ADMITTANCE
ENTER 1 IF O.K., 0 OTHERWISE: 1
: ~; G! T- h! e0 X  O! K5 j      2 OF      2 LOADS CONVERTED
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fact
# c6 I) m" k, D7 g+ j2 Z2 zfact
1 [0 Y& Q2 [0 G; L9 Q$ r ORDERING NETWORK
! `  Z# G# H9 E3 O. c2 [, m- a' Z DIAGONALS =     11  OFF-DIAGONALS =     10  MAX SIZE =     20
' I. Y- d( v% d9 ~( w     11 DIAGONAL AND     10 OFF-DIAGONAL ELEMENTS
1 l+ f- {4 a2 t; ]6 h! i ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
/ M; \% b/ ^9 @* g, n& o Executing activity tysl
tysl
% b, B5 r8 t  y( C ITER DELTAV/TOL X-------- AT BUS ---------X  REAL(DELTAV)  IMAG(DELTAV)
* b2 B) _% f, D7 L    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]
. O5 ~) k8 P" [/ j' H2 m$ R SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity save,4g2a00
save,4g2a00
8 ^: O, q7 V' t+ W$ J CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:07
, r* o3 c1 `1 l) h8 {+ l ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity RTRN
RTRN
ACTIVITY?
& {, z. f4 g0 e+ `" n6 B Executing activity DYRE
& x0 p9 E+ N5 y$ ^( W+ k OUT OF FILE DATA--SWITCH TO TERMINAL INPUT MODE
GENERATOR MODELS USE:
1 ]; w- _& c  q! p  CONS        1-    56
; V8 Q7 Z$ ]3 O" ?  STATES      1-    24

3 w" K- A8 q, x1 b6 q8 H8 O8 @& b1 j SUMMARY OF MODELS READ:
GENS:  GENROU
3 w  H% x7 k$ _; V* s           4
NEXT AVAILABLE ADDRESSES ARE:
    CON  STATE    VAR   ICON
/ X+ S9 d: a$ s" Q     57     25      1      1
NO MODEL CALLS IN CONNECTION ROUTINES--DYNAMICS SKELETON MAY BE USED
Executing activity CHAN
8 m- p8 C) {& C, j9 E3 [. O
      3 OUTPUT CHANNELS HAVE BEEN ADDED
NEXT AVAILABLE ADDRESSES ARE:
& A5 M1 B* Q' A3 }- i3 F$ \7 E2 b0 ^   CHANNEL   VAR   ICON
4 H: j4 U. J) Z3 y# i, [8 ?. [" Q        4      2      2
Executing activity STRT
GENROU 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
GENROU AT BUS    301 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
     PTI INTERACTIVE POWER SYSTEM SIMULATOR--PSS/E      SAT, JUN 02 2012  15:08
2 N4 O/ D0 ^5 K9 e 1
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
      1 1           20.000 1  1.0300 2.4652  700.00  185.25 0.9667  70.87 0.7104 0.3249
( d2 D, m! Q0 O) [6 `    101 2           20.000 1  1.0100 2.5768  700.00  234.84 0.9481  59.45 0.7504 0.3109
; H5 l- J0 m/ |. F' t    201 3           20.000 1  1.0300 2.4910  719.00  176.22 0.9713  44.98 0.7270 0.3303
  C3 p: x6 n+ _5 q" ?+ n, D: B9 ]1 g    301 4           20.000 1  1.0100 2.5189  700.10  202.33 0.9607  33.50 0.7359 0.3181
) ^) W8 v2 L  N8 v INITIAL CONDITIONS CHECK O.K.
! L& H6 _* n* G; B) g0 P' P# f SNAPSHOT STORED IN FILE D:\PSSE30\EXAMPLE\4g2a.snp AT TIME =  -0.020
CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out
: r, ^6 j' H+ ]5 F0 I) a+ I3 y NETWORK NOT CONVERGED AT TIME =  0.1700
  26   191.996    801 [9           230.00]     0.1697E-01     0.8974E-02
5 L6 R' p) W7 [( |( D& GDecimal number outside allowable range!
CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out


" I: Q* E- {" h7 z潮流数据预处理应该没问题，dyr文件只包括发电机模型，用PSASP6型模型仿真结果正常，到PSS/E就出问题了
) S; \5 I+ ]" a
- E+ B( ?: t1 B9 @% E. v最后两项饱和度是自己随便写的，其他数据来自kunder的书

apsal

问题很明显，你在执行STRT初始化操作时提示：
9 X5 B% ]5 A# ~. S8 n* E& GExecuting activity STRT
GENROU AT BUS      1 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
% D! a/ T* y3 U2 ]3 a2 b& kGENROU 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
GENROU AT BUS    301 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
& E, x" F, }5 j2 N4 P# S; O: b
% Z4 X9 w1 e! C这表明，你在潮流数据中的ZSORCE（CONG操作中将发电机等效为Norton电流源就是用的这个ZSORCE）与动态仿真中的发电机的次暂态电抗X"不相等。两个不相等，仿真就是错的。
5 M; ]* Q6 N4 d& b0 {这就是问题的所在。
8 N2 y8 C( q0 ^  s2 `6 S, B看起来你的X"可能是对的，ZSORCE用的的默认值1.0，那么将潮流数据中发电机的ZSORCE的电抗值改为对应X"就可以了。

actfate

厉害！！！光看INITIAL CONDITIONS CHECK O.K.去了
  y2 m+ c8 l$ s! d多谢

mingyu

学习一下。

wlm_28

看不懂呀

lxself

强啊强啊强啊强啊强啊

DUTshen

谢谢！谢谢！谢谢！

BPA_CN

同问，哪位大神指教一下？