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

bao1872

问题是conec文件编译后不能生成obj文件，请高手帮看看代码


. e! y. v- H3 Y需要自定义一个逆变器模型，数学模型如下：

数学模型

0 F% R8 R& a; l: VPref (pu) 3-phase total real power reference, per- unit value
1 ~) g' j5 [: ~/ e$ uvPCC (pu) voltage in PCC, per-unit line-line RMS value
Iref_p (pu) real power current reference, per-unit phase RMS value
Ip (pu) real power current generated by inverter, per-unit phase RMS value
) n# k6 p* ^! O
- W6 d3 f4 O& p5 V
dyr文件内容如下：
( m  h9 r5 f6 Y4 R, N11 'USRMDL' 1 'SLX_G1' 101 1 0 2 2 0  0.01 0.004166666666667   /


, R6 ]1 B8 C  c* G第一个con是指参考功率Pref，第二个con是指Tinv

. p6 ?) ^; b8 k% V/ N# b
首先在生成CONCE CONET的时候，没有出来call语句，psse33版好像都没有call语句了
9 ]* }+ A  Q, x; E0 e
7 k1 I) z) A" G2 n) \' h* H. {
7 L3 Q6 x' `) x然后编辑CONEC 文件如下：
+ f$ O: I) `' {" O1 c- G; n
5 Z3 ]# Y& I, R+ U4 m    SUBROUTINE CONEC
C
3 R/ s1 h/ o+ F    INCLUDE 'COMON4.INS'
C
C
    RETURN
    END
6 U0 I- v7 s" Y# u5 O/ i' @5 d9 J1 P; @
SUBROUTINE SLX_G1(MC,ISLOT)

INTEGER MC,ISLOT
C
% T: w8 d& r( w1 Z! e8 G' C( M, [C MC =MACHINE ARRAY INDEX
7 r2 i& [& `0 g/ \$ t$ m% tC ISLOT=ARRAY ALLOCATION TQABLE INDEX
) s5 |$ n) R0 ]C J =STRTIN(2,ISLOT) [USES CON(J) THROUGH CON(J+1) ]
/ W- Y1 x7 L2 k0 c' i: AC K =STRTIN(1,ISLOT) [USES STATE(K) THROUGH STATE(K+1)]
C

# v! F7 k/ Z( n$INSERT COMON4
" C9 \5 o& o* y5 @( lC
+ ^( E7 S8 [# }, U
INTRINSIC MAX

# t" @% l2 S- G. `7 U6 ~EXTERNAL BADMID
1 S( T4 U# ?- \5 ^- d1 W  }C
; W; r3 Q2 V5 N! `+ e
INTEGER IB, J, K

REAL IP
& \" r- ~) y2 g4 N& s% w. WC
C DO NOTHING IN ACTIVITIES ‘DOCU’ AND ‘DYDA’ FOR NOW
* n, Q) O* t( _0 Z3 g9 |0 G' |C

IF(MODE.EQ.8)

CON_DSCRPT(1)='Pref'

4 X  E- z) \  J9 W8 C3 @4 pCON_DSCRPT(2)='Tinv'

RETURN
- d7 A( R% q0 I6 g) d. l" d
IF (MODE.GT.4) RETURN
C
9 F: t: Y7 a* _7 _0 Q  `C BUS SEQUENCS NUMBER NEGATIVE IF MACHINE
; ^0 @. G9 Y5 V7 UC IS OFF LINE,SVS OR INDUCTION MACHINE
C
* R4 H4 ~4 K" H) V/ t5 ?
IB=NUMTRM(MC)
- a) x. }  }% i$ P
# m. Z. M% A7 p8 S4 xIF (IB.LE.0) RETURN
4 `. r/ ^* F. F- V+ BC
- A$ ?+ g4 X, n  @+ S3 l8 K+ r& ^
- B; u+ ?7 t2 D+ V8 H7 E$ \  |; }) FIF (MIDTRM) GO TO 900
6 M, ]! q0 t, n- MC
C GET STARTING ‘CON’ AND ‘STATE’ INDICES
C
) }0 S& B2 H' W7 ~; ?
J=STRTIN(1,ISLOT)
4 i$ t4 U2 A; a
4 K! Z; p; M0 a/ g( y" x( i- `K=STRTIN(2,ISLOT)
$ A3 ^2 G8 l3 OC
% w! i- f7 v9 {. p6 d
GO TO (100,200,300,400),MODE
C
( u0 ^! L3 ~9 `& BC MODE = 1 – INITIALIZE
; k9 c# w; T' a; P) Q% i$ L/ R* p3 \C
4 a4 v; _- k/ X7 e
100 STATE(K)=PELEC(MC)/ETERM(MC)

STATE(K+1)=PELEC(MC)/ETERM(MC)
; ^2 \# G, f& {5 L  E
RETURN
8 _  U% C6 |3 l: c$ I( OC
6 O) a. y+ K( o9 _C MODE = 2 – CALCULATE DERIVATIVES
C

- S7 s) j. c& r  N200 DSTATE(K)=CON(J)/ETERM(MC)

DSTATE(K+1)=(DSTATE(K)-STATE(K+1))/CON(J+1)
% N. [' [# }3 U( b2 K  b6 p1 m/ P
# C0 r) u. W. p/ H! ~RETURN
C
C MODE = 3 – SET IP
% X6 x. p0 h( ?! HC
# V, ^  s) q; s7 G
300 IP=STATE(K+1)

RETURN
5 ]4 \3 F* Y5 z( B# SC
C MODE = 4 – SET NINTEG
C
) v7 X. e  ]( _: }+ j9 n
% M* x4 b" V6 I6 M400 NINTEG=MAX(NINTEG,K+1)

" }" W( Y8 W; MRETURN
C
* W! J0 P3 A! Z5 RC MODEL NOT IMPLEMENTED FOR MSTR/MRUN
& f6 B) g: R2 q2 UC
) m% i3 `+ U* N* o
9 e7 F8 C7 \, h* X$ Z0 T! _900 CALL BADMID(MC,IB,’SLX_G1’)
$ Y8 u$ ^  F+ Y! v- F6 j" \6 |
. z) J# {) U& A7 ~RETURN
C

307 FORMAT(‘SLX_G1 AT BUS’,16, ‘MACHINE’,A,
! ~: G. ]+ o0 C3 S9 F0 s: z) D
0 L; a6 Z2 A, q  p) R) P* , ‘INITIALIZED OUT OF LIMITS’)

/ z5 B" W% H* QEND

  H; l$ j$ o6 o9 I( w

还有一个问题想请教最后Ip这个参量在PSSE里好像没有内部变量，因为这个是输出变量，现在Ip表示的好像不对，又不知道该怎么解决

Khalifa

强烈支持楼主ing……

cxxgtxy

回复 1# bao1872
1 t1 A8 [, @- M2 @0 W; H/ b9 T1 G! @5 S1 u

    Simulink不是很擅长做这个么？为啥不用，然后调用simulink

redplum

有注释吗

redplum

的的确确看不懂

大萝卜

顶楼主

bao1872

回复 2# cxxgtxy

4 x$ M% `" }9 G
% j9 b8 c- }2 P" u. q  我是33版的psse，已经不能调用simulink了，如果能装回psse30的话是否可以调用simulink来实现？能否推荐一个帮助看看，谢谢！

bao1872

回复 3# redplum
2 w' v# [* j& s' T  X. g
# K  y" a9 x0 ~# E( \  _
   其实主要内容是mode1-mode4，注释都是用英文表示的，C后面都是注释

cxxgtxy

回复 6# bao1872


# \0 b8 ]3 |2 l0 m; e/ d    http://wenku.baidu.com/view/be41c8c3d5bbfd0a79567382.html
不错的内容，希望能帮助你

bao1872

回复 8# cxxgtxy

0 u$ H4 o7 b) Y6 h6 P  C' Z7 U
: q  x& `& X2 s4 V) }; n   谢谢了，这个文档我已经下载过，没有具体讲simulink和psse怎么互相导入，而且上面很明确的写了不推荐用simulink。。。。因为容易产生未知错误，当时考虑过这个方法，但是自己用fortran写出来的话更容易让人信服

cxxgtxy

回复 9# bao1872

* `1 [0 \% D$ R1 H# E4 g, x' O
    其实感觉，这种情况下升级一下simulink或许就可以解决了，你的simulink版本是？