求助四机两区域动态仿真

actfate

潮流结果和PSASP的完全对得上，但在不设故障情况下在0.17s失稳，G1的Pelec如图
# `7 g8 S3 B" ?1 |/ y
; I" p: }: p3 s8 {" L2 D8 p
) ]& m! H% U9 A; E4 N2 Q* Llofl
+ }( W. ]' o& O  F ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity read,file
read,file
ENTER INPUT FILE NAME (0 TO EXIT, 1 FOR TERMINAL): 4g2a
8 C9 _& O6 E: G7 B$ g1 q$ q ENTER IC, SBASE
ENTER TWO LINE HEADING
ENTER BUS DATA
ENTER LOAD DATA
9 N) q7 a6 `5 A( _( A+ P ENTER GENERATOR DATA
ENTER NON-TRANSFORMER BRANCH DATA
9 [5 u, H% z7 [' Q7 Q" } ENTER TRANSFORMER DATA
3 T" z" Q3 X1 y8 c ENTER AREA INTERCHANGE DATA
% |% c* J* m: c2 ~( X5 h ENTER TWO-TERMINAL DC LINE DATA
0 X0 L: O' ~0 m  J' z& e4 u7 o ENTER VSC DC LINE DATA
5 M0 E# h! [. Q8 [  u3 I  A6 `' ? ENTER SWITCHED SHUNT DATA
" c  {5 T9 O0 q! R ENTER TRANSFORMER IMPEDANCE CORRECTION DATA
+ f) Y% X- N. I8 l2 c4 b ENTER MULTI-TERMINAL DC LINE DATA
ENTER MULTI-SECTION LINE DATA
ENTER ZONE NAME DATA
" N7 F8 d7 A  m ENTER INTER-AREA TRANSFER DATA
ENTER OWNER NAME DATA
% E: ?9 B5 }& p# _" {& L/ Y ENTER FACTS DEVICE DATA
" a4 }. {- q% t" F6 T ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fnsl
fnsl
2 {) O; Y! Q% R6 F" l ORDERING NETWORK
DIAGONALS =     10  OFF-DIAGONALS =      9  MAX SIZE =     18
6 F6 Y2 d  ~, | ENTER ITERATION NUMBER FOR VAR LIMITS
7 ?2 A8 e3 J2 d9 y  0 FOR IMMEDIATELY, -1 TO IGNORE COMPLETELY: 99
( a4 r/ t4 r& M% a0 Y. l  ITER       DELTAP      BUS         DELTAQ      BUS        DELTA/V/      BUS       DELTAANG      BUS
" W" V8 k" Y& u' f   0        21.2968(      401   )    2.7013(      801   )    0.12761(      801   )   0.58266(      801   )
# V7 c0 y( `6 I7 j& o' m   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 f: f4 X7 `1 H) P9 u7 U   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   )
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
' S7 L$ ~! i7 a SWING BUS SUMMARY:
   BUS# X-- NAME --X BASKV      PGEN     PMAX    PMIN      QGEN     QMAX    QMIN
0 p0 b# F, X3 y8 Z7 m( y    301 4           20.000     700.1   9999.0 -9999.0     202.3   9999.0 -9999.0
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity save,4g2a00
9 Q. Q8 b" x$ o1 d1 Msave,4g2a00
CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:06
# \# y2 m5 A' `6 K0 e. B ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
% Q- H" @* D+ T# C, `% f3 P* S Executing activity cong
+ N- y/ W" F1 x5 D8 |3 C7 Bcong
) B& k* C5 u$ f* Z# w GENERATOR CONVERSION COMPLETED
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
, r$ v3 r; _4 }2 n' S5 f Executing activity conl all
conl all
# [) T7 S; u; K0 d+ q/ s ENTER % CONSTANT I, % CONSTANT G FOR ACTIVE POWER:
ENTER % CONSTANT I, % CONSTANT B FOR REACTIVE POWER:
LOAD TO BE REPRESENTED AS:
    ACTIVE  REACTIVE
   100.00%   100.00% CONSTANT POWER
3 _2 q$ R, h* f, D. J! E3 F# Q     0.00%     0.00% CONSTANT CURRENT
     0.00%     0.00% CONSTANT ADMITTANCE
& ?7 O: M* {& l+ k' c ENTER 1 IF O.K., 0 OTHERWISE: 1
0 Y' F& V' r& [8 Y# M      2 OF      2 LOADS CONVERTED
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fact
fact
ORDERING NETWORK
DIAGONALS =     11  OFF-DIAGONALS =     10  MAX SIZE =     20
     11 DIAGONAL AND     10 OFF-DIAGONAL ELEMENTS
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity tysl
tysl
7 R, z- C' x+ T8 C, v! ]) H ITER DELTAV/TOL X-------- AT BUS ---------X  REAL(DELTAV)  IMAG(DELTAV)
$ I$ M' V/ y, i+ v3 o: H    1     0.003       1 [1           20.000]    0.3029E-07   -0.1089E-07
REACHED TOLERANCE IN   1 ITERATIONS
: L  o+ d* ]: @2 n' c) N- F; l8 q LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    601 [7           230.00]
! P7 G% f% n6 B' }. }5 m SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
7 T, w4 b  ?- y" D Executing activity save,4g2a00
save,4g2a00
CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:07
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity RTRN
RTRN
ACTIVITY?
Executing activity DYRE
+ _5 H: y. i; a! Q6 M4 h OUT OF FILE DATA--SWITCH TO TERMINAL INPUT MODE
GENERATOR MODELS USE:
  CONS        1-    56
  STATES      1-    24

2 C! k4 {+ v* Z1 M SUMMARY OF MODELS READ:
, k" A3 L& K# s4 D/ O GENS:  GENROU
5 I' \( ~; \: U# A( H) ^           4
9 q7 B3 q+ r$ ` NEXT AVAILABLE ADDRESSES ARE:
    CON  STATE    VAR   ICON
: }' V! v& M9 i3 M9 p: |     57     25      1      1
NO MODEL CALLS IN CONNECTION ROUTINES--DYNAMICS SKELETON MAY BE USED
Executing activity CHAN

      3 OUTPUT CHANNELS HAVE BEEN ADDED
NEXT AVAILABLE ADDRESSES ARE:
) }6 N/ \% a! |: D; s   CHANNEL   VAR   ICON
        4      2      2
: @9 l$ x7 v! g" V2 N& d Executing activity STRT
" H+ M" `# |  s+ b# {+ m" W2 h 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
1 Y6 P! A% |+ q" Z/ X GENROU AT BUS    201 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
: X1 [" D) h7 |3 p* R8 [ 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 h9 T- T" W- k 1
( @4 N; ]( T1 `" ]- Q0 J 2
$ ]6 M: }& {; t$ A. i 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
    101 2           20.000 1  1.0100 2.5768  700.00  234.84 0.9481  59.45 0.7504 0.3109
; h% R- V: v5 P  {. q8 v9 m    201 3           20.000 1  1.0300 2.4910  719.00  176.22 0.9713  44.98 0.7270 0.3303
8 }0 q& J8 D0 z    301 4           20.000 1  1.0100 2.5189  700.10  202.33 0.9607  33.50 0.7359 0.3181
INITIAL CONDITIONS CHECK O.K.
' `' e( w/ ]' d  H' i0 |* b SNAPSHOT STORED IN FILE D:\PSSE30\EXAMPLE\4g2a.snp AT TIME =  -0.020
& A1 i* l; V' K; r1 t$ b1 k CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out
NETWORK NOT CONVERGED AT TIME =  0.1700
' v/ [( J$ }; I$ C& X9 w  26   191.996    801 [9           230.00]     0.1697E-01     0.8974E-02
Decimal number outside allowable range!
CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out

" f" j( M# r* O% C0 e$ v3 D
潮流数据预处理应该没问题，dyr文件只包括发电机模型，用PSASP6型模型仿真结果正常，到PSS/E就出问题了

最后两项饱和度是自己随便写的，其他数据来自kunder的书

apsal

问题很明显，你在执行STRT初始化操作时提示：
Executing activity STRT
3 l, D3 ?& g1 Z+ U% _- mGENROU AT BUS      1 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
* B6 B% K# _3 }( P+ L, L  p; LGENROU AT BUS    101 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
$ N' @- ?' N6 `) M7 jGENROU AT BUS    201 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
8 x2 `3 e! P' k2 HGENROU AT BUS    301 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000

0 B, V& T. {" w7 d' [这表明，你在潮流数据中的ZSORCE（CONG操作中将发电机等效为Norton电流源就是用的这个ZSORCE）与动态仿真中的发电机的次暂态电抗X"不相等。两个不相等，仿真就是错的。
这就是问题的所在。
% r& w! C% F: k- |看起来你的X"可能是对的，ZSORCE用的的默认值1.0，那么将潮流数据中发电机的ZSORCE的电抗值改为对应X"就可以了。

actfate

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

mingyu

学习一下。

wlm_28

看不懂呀

lxself

强啊强啊强啊强啊强啊

DUTshen

谢谢！谢谢！谢谢！

BPA_CN

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