求助四机两区域动态仿真

actfate

潮流结果和PSASP的完全对得上，但在不设故障情况下在0.17s失稳，G1的Pelec如图
3 q: g& G/ u' ^% Y2 \- e2 Z
8 A2 ~: Q& e8 ^7 j& ~
lofl
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity read,file
read,file
ENTER INPUT FILE NAME (0 TO EXIT, 1 FOR TERMINAL): 4g2a
. t: e$ |3 [" g# Q8 |- ^- x  v ENTER IC, SBASE
/ E& H' S# C8 G ENTER TWO LINE HEADING
+ C% k3 c0 _$ [* e# \/ v ENTER BUS DATA
ENTER LOAD DATA
7 N/ ]/ z/ ^, |; O: t1 q ENTER GENERATOR DATA
ENTER NON-TRANSFORMER BRANCH DATA
+ H1 c2 ^7 V) |5 o  t7 y ENTER TRANSFORMER DATA
. l; ]! i$ Y1 I: [1 Y ENTER AREA INTERCHANGE DATA
ENTER TWO-TERMINAL DC LINE DATA
+ p0 m6 }; F% X7 D& v+ \ ENTER VSC DC LINE DATA
ENTER SWITCHED SHUNT DATA
6 j3 J$ r; F& p& y' I# p& D1 ^ ENTER TRANSFORMER IMPEDANCE CORRECTION DATA
7 |; N; ?8 V$ K! l8 q, x, Y7 a ENTER MULTI-TERMINAL DC LINE DATA
$ q, K; l/ k: c# i) i/ `9 s+ W ENTER MULTI-SECTION LINE DATA
/ u" R- O1 y- N4 q0 S ENTER ZONE NAME DATA
ENTER INTER-AREA TRANSFER DATA
& G7 x9 l5 O$ g' j2 y7 D; R ENTER OWNER NAME DATA
ENTER FACTS DEVICE DATA
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fnsl
fnsl
ORDERING NETWORK
: `* u+ J# A. C! f DIAGONALS =     10  OFF-DIAGONALS =      9  MAX SIZE =     18
ENTER ITERATION NUMBER FOR VAR LIMITS
% I. ?8 M, e  j# k6 [  0 FOR IMMEDIATELY, -1 TO IGNORE COMPLETELY: 99
  ITER       DELTAP      BUS         DELTAQ      BUS        DELTA/V/      BUS       DELTAANG      BUS
$ J6 V( y& C* }5 }   0        21.2968(      401   )    2.7013(      801   )    0.12761(      801   )   0.58266(      801   )
# ~+ X8 u. Q7 m- F0 J# W, a% v: i4 V   1         4.2106(      801   )    8.1793(      901   )    0.12606(      701   )   0.04352(        1   )
0 {( b. ?. b' N: {   2         0.2518(      901   )    0.7799(      901   )    0.02045(      701   )   0.00508(        1   )
, B% e/ w9 Q$ `   3         0.0122(      901   )    0.0124(      901   )    0.00048(      701   )   0.00029(      801   )
0 G2 d/ [1 z' N$ r8 L& M2 b   4         0.0000(      801   )    0.0000(      901   )    0.00000(      801   )   0.00000(      801   )
  T# M9 ?, ]8 g: p6 F) _% k   5         0.0000(      801   )    0.0000(      801   )
+ g- \/ I7 A4 O7 n" j, B7 k9 u2 f REACHED TOLERANCE IN   5 ITERATIONS
LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    801 [9           230.00]
. L# Z0 J# F7 b: X! f. @9 B! A SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
# y$ z. Y$ i6 z 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?
9 g; K: k, p7 W$ G- m; D- Z% I% ^$ N Executing activity save,4g2a00
save,4g2a00
CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:06
" u" ^$ L3 T* T; b  D3 ^% p ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity cong
0 S: P2 B1 Z$ I+ e$ e8 n. Mcong
$ J2 {- Z. V" e3 u7 R GENERATOR CONVERSION COMPLETED
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
- G8 b! m& L6 R6 u) F. S Executing activity conl all
conl all
ENTER % CONSTANT I, % CONSTANT G FOR ACTIVE POWER:
ENTER % CONSTANT I, % CONSTANT B FOR REACTIVE POWER:
LOAD TO BE REPRESENTED AS:
( |0 M6 W3 C. G$ ^0 G  `    ACTIVE  REACTIVE
   100.00%   100.00% CONSTANT POWER
     0.00%     0.00% CONSTANT CURRENT
) F  T4 \5 F. t3 _- [     0.00%     0.00% CONSTANT ADMITTANCE
ENTER 1 IF O.K., 0 OTHERWISE: 1
& U0 |2 h5 y3 S# P; h( N  V      2 OF      2 LOADS CONVERTED
$ Q* Y4 R3 z6 F7 B ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity fact
fact
ORDERING NETWORK
; E; t1 H: t9 g' h DIAGONALS =     11  OFF-DIAGONALS =     10  MAX SIZE =     20
3 v" q1 T' m  }; H0 E     11 DIAGONAL AND     10 OFF-DIAGONAL ELEMENTS
ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
: T( i9 l& h, E8 d3 H- ^4 ^# F, S Executing activity tysl
tysl
8 c8 C! ~* D  H7 M6 b ITER DELTAV/TOL X-------- AT BUS ---------X  REAL(DELTAV)  IMAG(DELTAV)
: Z4 X. z' J' |, n9 P' p) E    1     0.003       1 [1           20.000]    0.3029E-07   -0.1089E-07
+ P# ^, |8 b) R$ _7 R5 K" a1 w; K REACHED TOLERANCE IN   1 ITERATIONS
; o+ d- s( E' z, P  d LARGEST MISMATCH:    0.00 MW    0.00 MVAR     0.00 MVA AT BUS    601 [7           230.00]
' T9 |, ~: l' o, A- F SYSTEM TOTAL ABSOLUTE MISMATCH:               0.00 MVA
2 Q( F& l; b- k; J/ H, u ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
Executing activity save,4g2a00
+ Z  i2 X% D$ tsave,4g2a00
  }0 ?# h5 v% p7 Z! v3 A# \ CASE SAVED IN FILE 4g2a00.sav ON SAT, JUN 02 2012  15:07
% `$ \; J+ K4 a! g0 v ACTIVITY RTRN TO RETURN TO DYNAMICS--ACTIVITY?
) z5 w7 F5 o7 |1 t* T Executing activity RTRN
RTRN
) l5 b, B* {; ~# J! L ACTIVITY?
. J5 H: \: n. z8 o, r9 i: R Executing activity DYRE
OUT OF FILE DATA--SWITCH TO TERMINAL INPUT MODE
9 O% Z: k, Q5 T+ C7 t+ G% h GENERATOR MODELS USE:
  CONS        1-    56
  STATES      1-    24
+ ?3 P! M4 `2 F" G/ b
SUMMARY OF MODELS READ:
7 [% l. M/ {  z; R+ i. ^/ k9 S GENS:  GENROU
! S6 U" V. ?7 F; a7 m           4
, E  c$ {0 E" U! S% a! Y' k+ s NEXT AVAILABLE ADDRESSES ARE:
    CON  STATE    VAR   ICON
     57     25      1      1
' y. o4 D. n6 ? NO MODEL CALLS IN CONNECTION ROUTINES--DYNAMICS SKELETON MAY BE USED
& H) }) c: [/ {7 W, }) B' j) v) b1 k# s Executing activity CHAN
  A3 }8 ?+ e) S) ]8 M3 [
      3 OUTPUT CHANNELS HAVE BEEN ADDED
& `- S9 P  t7 c+ x2 b+ x NEXT AVAILABLE ADDRESSES ARE:
4 r) a/ H0 c9 p4 U) i# B   CHANNEL   VAR   ICON
        4      2      2
3 F3 ~( V3 O# N" C# Z Executing activity STRT
8 ?$ E6 C& s# Y/ z1 p 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
5 t9 g) p9 f, k( c& w" ?/ f2 E. e 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
3 s  G" Q3 x2 `! M5 X     PTI INTERACTIVE POWER SYSTEM SIMULATOR--PSS/E      SAT, JUN 02 2012  15:08
1
1 A9 _7 L* _! X! Y* v  G4 ?, e 2
INITIAL CONDITION LOAD FLOW USED  1 ITERATIONS
----------------------------- MACHINE INITIAL CONDITIONS -----------------------------
, P: E7 ~6 o) p8 f" U   BUS# X-- NAME --X BASKV ID  ETERM   EFD    POWER    VARS   P.F.  ANGLE   ID     IQ
! k* ]( `! z( v2 G      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
$ p, R0 ^# v4 ?( b    201 3           20.000 1  1.0300 2.4910  719.00  176.22 0.9713  44.98 0.7270 0.3303
; b+ ]1 L7 e$ R" C7 X3 d. v5 O    301 4           20.000 1  1.0100 2.5189  700.10  202.33 0.9607  33.50 0.7359 0.3181
" U# f- z' J0 c" A4 O3 M* B& h6 Z INITIAL CONDITIONS CHECK O.K.
SNAPSHOT STORED IN FILE D:\PSSE30\EXAMPLE\4g2a.snp AT TIME =  -0.020
( U! r) ^) f/ ~. I: w" u CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out
* X: k( l" y& J4 l NETWORK NOT CONVERGED AT TIME =  0.1700
" }$ f+ J  |' ^  26   191.996    801 [9           230.00]     0.1697E-01     0.8974E-02
4 `. |: m4 y3 qDecimal number outside allowable range!
CHANNEL OUTPUT FILE IS D:\PSSE30\EXAMPLE\4g2a.out

( v/ g  E% ^6 l+ ~, M
+ E4 B- k8 Q" J! P' r潮流数据预处理应该没问题，dyr文件只包括发电机模型，用PSASP6型模型仿真结果正常，到PSS/E就出问题了

: A* m3 m) h2 m5 W/ P7 o/ E最后两项饱和度是自己随便写的，其他数据来自kunder的书

apsal

问题很明显，你在执行STRT初始化操作时提示：
Executing activity STRT
8 i7 f" T3 N2 o6 ^4 lGENROU AT BUS      1 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
# F: f$ r+ z+ y2 sGENROU AT BUS    101 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
% q* u0 V0 P0 t" U2 {& ?* ?  uGENROU AT BUS    201 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
, d, y  ^. P+ |7 fGENROU AT BUS    301 MACHINE 1 : X'' =0.25000 IM(ZSORCE) =1.00000
  _  ?+ W" J8 P' q, J7 A
这表明，你在潮流数据中的ZSORCE（CONG操作中将发电机等效为Norton电流源就是用的这个ZSORCE）与动态仿真中的发电机的次暂态电抗X"不相等。两个不相等，仿真就是错的。
# `) `# I2 L5 ]" i这就是问题的所在。
看起来你的X"可能是对的，ZSORCE用的的默认值1.0，那么将潮流数据中发电机的ZSORCE的电抗值改为对应X"就可以了。

actfate

厉害！！！光看INITIAL CONDITIONS CHECK O.K.去了
6 r+ o  P. i$ ~% V" k多谢

mingyu

学习一下。

wlm_28

看不懂呀

lxself

强啊强啊强啊强啊强啊

DUTshen

谢谢！谢谢！谢谢！

BPA_CN

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