逆变器自定义模型求助，编译不成功

bao1872

问题是conec文件编译后不能生成obj文件，请高手帮看看代码
6 S$ E7 `4 e$ N- q
! \7 R+ h  p0 T
5 O+ U' I- E. A7 N3 H5 U+ V# f需要自定义一个逆变器模型，数学模型如下：

数学模型

Pref (pu) 3-phase total real power reference, per- unit value
/ d8 ^: O. N. M, tvPCC (pu) voltage in PCC, per-unit line-line RMS value
Iref_p (pu) real power current reference, per-unit phase RMS value
- q, s' O& g5 I# d9 E' LIp (pu) real power current generated by inverter, per-unit phase RMS value


dyr文件内容如下：
11 'USRMDL' 1 'SLX_G1' 101 1 0 2 2 0  0.01 0.004166666666667   /


" V9 ?% E3 t+ y8 p/ L第一个con是指参考功率Pref，第二个con是指Tinv

& m# T  T( K& v: a9 K
首先在生成CONCE CONET的时候，没有出来call语句，psse33版好像都没有call语句了

, o# ?8 z- E9 @! t) a0 ]$ [# x
然后编辑CONEC 文件如下：
* _9 V8 P/ E. Y/ Y4 _
- g5 B  `# ^& ^, d) s3 F/ |    SUBROUTINE CONEC
6 O0 o( I+ A) v6 r$ t6 X" v* d* s+ VC
    INCLUDE 'COMON4.INS'
C
/ u) z3 k% m2 S. c3 I0 y$ RC
    RETURN
0 c  W5 [# E* @. A( s9 B) D    END
/ z6 w; y4 i! ^; l/ O, p6 b; O
SUBROUTINE SLX_G1(MC,ISLOT)
. n3 s( O2 T+ U" H% {! w& `3 P3 v
  @$ U1 d; l% L: G% m/ J" w; ]INTEGER MC,ISLOT
4 T6 S3 s  U$ P) OC
C MC =MACHINE ARRAY INDEX
C ISLOT=ARRAY ALLOCATION TQABLE INDEX
' y, Q$ M! v% }5 w& f( i4 CC J =STRTIN(2,ISLOT) [USES CON(J) THROUGH CON(J+1) ]
C K =STRTIN(1,ISLOT) [USES STATE(K) THROUGH STATE(K+1)]
* U4 e( O0 }0 n' |. B4 hC

$INSERT COMON4
C
& x9 g% j4 b" e; [3 v
INTRINSIC MAX

9 X/ z, t' T/ Y- uEXTERNAL BADMID
1 s! f' O6 }# C* ]  q! t8 DC
" P7 a* k* o& E) u& ]8 `
' Z# I8 V1 `: V5 wINTEGER IB, J, K
: C  I6 w7 E3 d
REAL IP
4 y! m. ^& y- b% Q  VC
C DO NOTHING IN ACTIVITIES ‘DOCU’ AND ‘DYDA’ FOR NOW
C
: j1 Z8 h2 T: X- P7 v; n
* y7 ?/ x# X& _% h& N4 aIF(MODE.EQ.8)

CON_DSCRPT(1)='Pref'

CON_DSCRPT(2)='Tinv'
6 z( H4 F1 v# q7 p5 m8 f
1 V% a4 \  w" n3 K. i8 u) ZRETURN
% |4 e+ P2 I% S* |2 M9 {5 L3 k
! E9 _7 x$ l& |! BIF (MODE.GT.4) RETURN
C
5 P0 u* q4 L: V  G2 kC BUS SEQUENCS NUMBER NEGATIVE IF MACHINE
C IS OFF LINE,SVS OR INDUCTION MACHINE
2 V; p5 R- [4 X% u' p, mC
8 I4 y# {2 _4 c
- P4 J8 Q# K/ q% R, X. LIB=NUMTRM(MC)

& r: z# j, H) a) k! D0 lIF (IB.LE.0) RETURN
3 L% Y& ?1 o+ F$ ?2 X; E1 e9 q$ r3 RC
: K) F+ h6 f' s
; M' x8 ^7 y. y9 ?- YIF (MIDTRM) GO TO 900
; Q$ M' L7 ^/ `5 Z9 g, E. VC
9 m2 e3 o. k2 [4 l% c; ~C GET STARTING ‘CON’ AND ‘STATE’ INDICES
C

J=STRTIN(1,ISLOT)
! Y5 o. Z5 ]$ ^: g
K=STRTIN(2,ISLOT)
* n8 m+ s. ~/ L+ f( EC
9 u1 O* T* ~0 V; t7 f3 V! m/ e
GO TO (100,200,300,400),MODE
C
C MODE = 1 – INITIALIZE
C
. [& P% R- `$ S! j4 }
% w- b0 \0 K/ I100 STATE(K)=PELEC(MC)/ETERM(MC)
. |, o  x7 G. P- l" F. ^) Y2 E
9 M0 D. O4 D) J" x/ Y; ~STATE(K+1)=PELEC(MC)/ETERM(MC)

RETURN
C
" y! k) S. V4 c1 Q, l$ BC MODE = 2 – CALCULATE DERIVATIVES
5 w7 p6 B+ x" E8 p' mC
5 m  h: Y1 F! n, {# L; `
) U$ w2 ~/ h, ~% m  I& [, H200 DSTATE(K)=CON(J)/ETERM(MC)

DSTATE(K+1)=(DSTATE(K)-STATE(K+1))/CON(J+1)

RETURN
C
C MODE = 3 – SET IP
C
3 A8 ^1 |/ R# t3 c9 ?0 M
' U, n: [: w. f9 M+ {6 F& V300 IP=STATE(K+1)
  _& v8 I5 m3 `* H$ |  F0 C5 A
RETURN
C
C MODE = 4 – SET NINTEG
5 M7 S) ?" T+ m# w, rC
' A0 B" {1 U2 x' Z$ k
+ Q/ S3 o) E2 K& }! a9 c400 NINTEG=MAX(NINTEG,K+1)

RETURN
C
C MODEL NOT IMPLEMENTED FOR MSTR/MRUN
C
& T0 ]9 P( S- Z, o
900 CALL BADMID(MC,IB,’SLX_G1’)

3 f  d1 C- _- ]; b8 Y) s& e5 }  {RETURN
( m6 ?4 C" h3 `2 [7 }C
% P/ G3 [2 d( D; T* M/ c! ]  k* _7 v, y
) v, v5 V1 T, A9 d8 m307 FORMAT(‘SLX_G1 AT BUS’,16, ‘MACHINE’,A,
3 b' c7 @+ h( u& r' i9 S
* , ‘INITIALIZED OUT OF LIMITS’)
9 I% u( R4 i( v- o, G2 B% W% B! R( x
8 F' T$ T" E8 @2 n* g; x1 ~3 _END
8 @  ]- \& m5 q* h/ G

3 N- i( _0 g8 W, @; D6 J2 P
还有一个问题想请教最后Ip这个参量在PSSE里好像没有内部变量，因为这个是输出变量，现在Ip表示的好像不对，又不知道该怎么解决

Khalifa

强烈支持楼主ing……

cxxgtxy

回复 1# bao1872
0 O# q* ?9 }4 ]
+ ?2 i  n$ K: n$ `
0 Y6 Y% f4 B) g: |, H    Simulink不是很擅长做这个么？为啥不用，然后调用simulink

redplum

有注释吗

redplum

的的确确看不懂

大萝卜

顶楼主

bao1872

回复 2# cxxgtxy
6 o2 i: f9 B0 b* C7 [

; i, N/ p/ L2 \0 v  我是33版的psse，已经不能调用simulink了，如果能装回psse30的话是否可以调用simulink来实现？能否推荐一个帮助看看，谢谢！

bao1872

回复 3# redplum

! L/ C+ N& }, {
& ?; A& Z+ N- e3 \   其实主要内容是mode1-mode4，注释都是用英文表示的，C后面都是注释

cxxgtxy

回复 6# bao1872
3 o7 U8 z7 a, p
3 f6 ?$ q3 g" y4 [
1 S. H5 O/ I0 F& \! h  `- @    http://wenku.baidu.com/view/be41c8c3d5bbfd0a79567382.html
- }; u1 k$ @5 t# w1 A( z不错的内容，希望能帮助你

bao1872

回复 8# cxxgtxy


   谢谢了，这个文档我已经下载过，没有具体讲simulink和psse怎么互相导入，而且上面很明确的写了不推荐用simulink。。。。因为容易产生未知错误，当时考虑过这个方法，但是自己用fortran写出来的话更容易让人信服

cxxgtxy

回复 9# bao1872
6 I# `, }' E. L! i

3 U8 v) x9 |6 t5 q) r7 f4 k1 l9 D    其实感觉，这种情况下升级一下simulink或许就可以解决了，你的simulink版本是？