求助四机两区域动态仿真

actfate

潮流结果和PSASP的完全对得上，但在不设故障情况下在0.17s失稳，G1的Pelec如图
9 O. n+ u- o9 H* i) U* x

lofl
+ }1 @* Y) F% y# M ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
+ J% d# l; V5 A/ C Executing activity read,file
read,file
ENTER INPUT FILE NAME (0 TO EXIT, 1 FOR TERMINAL): 4g2a
ENTER IC, SBASE
ENTER TWO LINE HEADING
ENTER BUS DATA
1 i) `: n" g7 t  }7 X: e ENTER LOAD DATA
ENTER GENERATOR DATA
ENTER NON-TRANSFORMER BRANCH DATA
1 w4 ?, w9 ~- Q8 d& | ENTER TRANSFORMER DATA
" v% E, I9 V$ E2 o% G# } ENTER AREA INTERCHANGE DATA
ENTER TWO-TERMINAL DC LINE DATA
ENTER VSC DC LINE DATA
2 h' g5 j9 T: S( G3 j, ^ ENTER SWITCHED SHUNT DATA
( N0 z( t# x# j' l4 S ENTER TRANSFORMER IMPEDANCE CORRECTION DATA
ENTER MULTI-TERMINAL DC LINE DATA
ENTER MULTI-SECTION LINE DATA
8 y: N' E1 R8 a' m2 i ENTER ZONE NAME DATA
" N. `) ]+ a' P$ D" g ENTER INTER-AREA TRANSFER DATA
ENTER OWNER NAME DATA
  x9 l& W( U& V" U4 }# Z ENTER FACTS DEVICE DATA
' m! \% Z# a3 O2 U, {; A ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fnsl
+ j- ]0 U: T& A! _3 l' qfnsl
ORDERING NETWORK
DIAGONALS =     10  OFF-DIAGONALS =      9  MAX SIZE =     18
$ \3 W' l; @$ a! R4 _8 C. V0 J ENTER ITERATION NUMBER FOR VAR LIMITS
  0 FOR IMMEDIATELY, -1 TO IGNORE COMPLETELY: 99
0 j8 S% H4 M, }' E$ C; N" j  ITER       DELTAP      BUS         DELTAQ      BUS        DELTA/V/      BUS       DELTAANG      BUS
" g2 }% j$ ^( X! V5 F! r   0        21.2968(      401   )    2.7013(      801   )    0.12761(      801   )   0.58266(      801   )
7 m, s. A2 L; D$ X3 S   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   )
/ p* @" ~  W$ r7 z5 Z* k+ Y   3         0.0122(      901   )    0.0124(      901   )    0.00048(      701   )   0.00029(      801   )
. a. q, R3 b3 ^" l1 _   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
  M5 [1 T: }* ~7 x6 ]# E. _ LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    801 [9           230.00]
2 @7 m2 p  \) I9 b8 E SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
SWING BUS SUMMARY:
   BUS# X-- NAME --X BASKV      PGEN     PMAX    PMIN      QGEN     QMAX    QMIN
    301 4           20.000     700.1   9999.0 -9999.0     202.3   9999.0 -9999.0
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
  g! k7 a, v  y+ C& [( c$ t Executing activity save,4g2a00
* C8 ?6 J3 E; m, }* c, w; B+ Usave,4g2a00
CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:06
) E! y! e: F* }+ b+ k ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity cong
- ]2 O6 F  L8 i. j8 c& i6 @cong
GENERATOR CONVERSION COMPLETED
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity conl all
0 X$ ^9 o9 U9 O+ V7 ]; qconl all
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
     0.00%     0.00% CONSTANT CURRENT
) q$ @* S& \( h     0.00%     0.00% CONSTANT ADMITTANCE
' c* Y5 C5 o; z, K ENTER 1 IF O.K., 0 OTHERWISE: 1
      2 OF      2 LOADS CONVERTED
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fact
fact
/ Y6 ^$ M* ~6 Z. v5 R! U' E6 b: q  R ORDERING NETWORK
+ p& Z4 `: O9 j2 H1 T5 ? DIAGONALS =     11  OFF-DIAGONALS =     10  MAX SIZE =     20
! r+ Q- f1 p1 F6 E% \     11 DIAGONAL AND     10 OFF-DIAGONAL ELEMENTS
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity tysl
7 G& t# i. H/ y+ otysl
ITER DELTAV/TOL X-------- AT BUS ---------X  REAL(DELTAV)  IMAG(DELTAV)
& N' ]! F' q% O+ m/ @7 x, M    1     0.003       1 [1           20.000]    0.3029E-07   -0.1089E-07
REACHED TOLERANCE IN   1 ITERATIONS
# Q  Q$ V, _5 N7 v+ Q: h2 x LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    601 [7           230.00]
SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
; r, s1 \  ~. w* J( o# `# m" q ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity save,4g2a00
, g2 f2 f3 |" \$ t- ]1 ?save,4g2a00
CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:07
2 E% a4 C3 n* ]4 z ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity RTRN
RTRN
, r; i8 q! Z3 r ACTIVITY?
( n' q( r! Q. g8 n Executing activity DYRE
OUT OF FILE DATA--SWITCH TO TERMINAL INPUT MODE
. N/ w1 O# q% q6 N7 H0 r" C; I3 x GENERATOR MODELS USE:
  CONS        1-    56
  STATES      1-    24

1 H: ?" q& e1 f0 h# H6 Z+ E6 | SUMMARY OF MODELS READ:
' |2 `$ V4 @3 t" M# A1 o GENS:  GENROU
           4
" b$ m! r9 c- q+ {! _, R NEXT AVAILABLE ADDRESSES ARE:
1 _6 X. s# j  T) {    CON  STATE    VAR   ICON
     57     25      1      1
8 M/ g5 }+ m! d- `4 s NO MODEL CALLS IN CONNECTION ROUTINES--DYNAMICS SKELETON MAY BE USED
Executing activity CHAN
' `0 a; G# A' y4 V* V, G6 f1 J
      3 OUTPUT CHANNELS HAVE BEEN ADDED
8 w9 |0 U; w9 l' r6 h- C- L* D NEXT AVAILABLE ADDRESSES ARE:
   CHANNEL   VAR   ICON
        4      2      2
Executing activity STRT
  ^; t( o4 r( e/ |9 }& L GENROU AT BUS      1 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
. s4 q4 O, h* L& Z( @ GENROU AT BUS    101 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
: _8 R0 L) N+ [ 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
# }$ t- t& I* g! Q+ F, ~/ p     PTI INTERACTIVE POWER SYSTEM SIMULATOR--PSS/E      SAT, JUN 02 2012  15:08
1
2
& N; @) c, k% t' G INITIAL CONDITION LOAD FLOW USED  1 ITERATIONS
( ?" t% R" ^5 n! l ----------------------------- 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
5 X0 r( f1 @, g# W! \    101 2           20.000 1  1.0100 2.5768  700.00  234.84 0.9481  59.45 0.7504 0.3109
- @( ^2 z1 E: \1 u) @" }4 `+ K    201 3           20.000 1  1.0300 2.4910  719.00  176.22 0.9713  44.98 0.7270 0.3303
    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.
SNAPSHOT STORED IN FILE D:\PSSE30\EXAMPLE\4g2a.snp AT TIME =  -0.020
8 j0 w  c/ X  Q5 O* F, B CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out
, J: o6 m" {; A$ Q  d NETWORK NOT CONVERGED AT TIME =  0.1700
  26   191.996    801 [9           230.00]     0.1697E-01     0.8974E-02
Decimal number outside allowable range!
) f- X) o. X9 Z' K2 |' E, j' ^ CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out


潮流数据预处理应该没问题，dyr文件只包括发电机模型，用PSASP6型模型仿真结果正常，到PSS/E就出问题了
6 z, }4 z! k  b: g. ^5 y
最后两项饱和度是自己随便写的，其他数据来自kunder的书

apsal

问题很明显，你在执行STRT初始化操作时提示：
Executing activity STRT
3 Q9 ?) {' }5 B! V  r$ fGENROU AT BUS      1 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
' ~3 @8 I( H5 F3 D# eGENROU 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
" b" h. g9 Z- z+ C2 H# s
这表明，你在潮流数据中的ZSORCE（CONG操作中将发电机等效为Norton电流源就是用的这个ZSORCE）与动态仿真中的发电机的次暂态电抗X"不相等。两个不相等，仿真就是错的。
这就是问题的所在。
看起来你的X"可能是对的，ZSORCE用的的默认值1.0，那么将潮流数据中发电机的ZSORCE的电抗值改为对应X"就可以了。

actfate

厉害！！！光看INITIAL CONDITIONS CHECK O.K.去了
' Y9 x& }2 I2 Q9 t5 V多谢

mingyu

学习一下。

wlm_28

看不懂呀

lxself

强啊强啊强啊强啊强啊

DUTshen

谢谢！谢谢！谢谢！

BPA_CN

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