电抗器重燃-model-tacsres例子分享

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
4 y; @0 p1 r( d
MODEL breaker
+ Y5 j. \, l  {6 e4 g  ]9 i4 `
" k5 r1 P" f0 {6 z' D7 E( qCOMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
, M8 M; e4 \' C

( Q8 W& z3 D& WThe model of this example  updates  the  present  state  of  the
% u" m$ `* U$ D4 E: E$ U6 ycircuit  breaker  by  taking   into   account   the   conditions
( K5 ?0 o7 a! C$ l# Pcontrolling four represented state transitions:  first  opening,
% Z! @+ w& {( c# T5 N; r1 ~reigniting, reopening, and full opening.

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT
0 f3 g+ F4 F0 v5 ?
2 b7 }+ M* E1 |( t/ l: r: u0 gDATA
  topen                -- [s]   earliest elec opening
. b3 u; A; h+ v% h  ich1                 -- [A]   chopping level at power frequency
1 }4 s0 U! r1 s6 @: Y# Q" T  slope                -- [A/s] current slope
9 ?$ x1 X, \' B2 `$ j$ V! f1 C  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
+ m% o) C: ?1 @# q, K/ u# e  reignit_w {dflt:1e-2}              -- [s]  reignition window

INPUT
2 X6 }  K# W7 w# t  v1   -- voltage across breaker terminals
4 [+ Q7 _+ K3 [9 t! \  v2
1 [. }4 X; y, o9 d" y  current  -- current through breaker
- `! l& S3 W. f0 x
VAR
4 H! x% _! T, m. M: a7 U4 w  K: X9 [  state    -- 0= fully open
/ P, o( b, j' }4 o# Y           -- 1= normal closed
           -- 2= tentatively open
4 G3 h% g9 [; T$ n1 D           -- 3= reignited
  first_opening  -- transition from 1 to 2
/ a; X, D) q9 M, ^8 B6 l  reigniting     -- transition from 2 to 3
9 L6 z+ i; `; Y- d, `  reopening      -- transition from 3 to 2
! I4 l. [  X* R* e, }1 V  fully_opening  -- transition from 2 to 0
  reopen_clock   -- state variable of 3
2 }2 T- r! b, f! z6 L  p/ R7 |  elapsed        -- time elapsed since tbeg
  cbstatus       -- open/close signal to switch
  deltav  
  ds
  nds
1 O) T4 Y8 f$ j# o$ Y' ?3 Q5 x+ ?  tdelta
  itopen
  ich2

OUTPUT
  cbstatus ,ds,nds

HISTORY
  current {dflt: 0}
( A# e3 z9 `9 L  n2 z  o
INIT
  state    :=1       -- normal closed
7 w: A9 n# N3 O( s; r: S; x  cbstatus :=closed
* F  K- v% G/ X  ds:=1
  nds:=-1
; b, D) C; m8 W2 m. e3 {: e# fENDINIT

; }! e& j" g* K! e6 X% k! w0 n, zEXEC
5 Y6 _2 z% t% X* p: [2 J" R% e deltav:=v1-v2
IF t>=topen THEN              -- otherwise do nothing
  tdelta:=(t-topen)*1000
6 q4 k4 q) E3 C  IF tdelta<=0.44
    THEN ds:=(38.5*tdelta**1.5)*1000
  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
/ C( {4 J5 M/ J2 w  M6 @7 q  ENDIF

  nds:=-ds
  E' m' a( j6 T2 R  itopen:=abs(current)
. E  v0 k* Z) E5 D' o. w# o$ ]  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))

; {. G+ l& W( U$ U; j0 D. r: O; z  first_opening :=false;  fully_opening :=false
  reopening     :=false;  reigniting    :=false

  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
. r4 k# W8 P4 a+ m        OR current*prevval(current)<0)
- v: Y! g! i4 ^" d9 k: S   THEN first_opening :=true
        elapsed:=t-topen
* B5 q: V/ L0 I2 p8 u& U8 k" K7 Q
1 A% S$ F- M% N, b; E. ~  ELSIF state=2 THEN           -- was tentatively open   

   IF elapsed>reignit_w  THEN
      fully_opening :=true
   ELSIF  abs(deltav) > ds
# e: }0 _0 T) c      THEN reigniting:=true
   ENDIF

  ELSIF state=3 THEN                  -- was reignited

6 B8 {5 [" K9 t6 Z! W6 }* y, K( J   reopen_clock :=reopen_clock +timestep

! P" ?- W6 s" T  Y$ q" J) s2 X  \5 C   IF reopen_clock >=reopen_delay
/ R: m. X+ g8 E- ~    AND (abs(current) <ich1
         OR current*prevval(current)<0)
    AND abs(deriv(current))<slope
: u5 m$ b* [  i( y        THEN reopening:=true
   ENDIF

  ELSIF state=0 THEN                 -- was fully open
  ENDIF

  IF first_opening OR reopening THEN
   state :=2               -- becomes tentatively open

  E1 [7 ?/ ?4 l5 l2 @6 X  ELSIF reigniting THEN
   reopen_clock :=0
   state :=3                      -- becomes reignited
  ELSIF fully_opening THEN
4 `' C* n4 K, W( C9 x3 b$ u   state :=0                     -- becomes fully open
/ P. g& @/ p) u) d0 ?8 ?; b  ENDIF
/ I. }/ B  h! ]: _
9 C8 g( A$ _- j/ A  IF state = 2 OR state = 0             -- is now open
  THEN cbstatus :=open
: `7 x- m. M% g" k: n( }' M0 y. A  ELSE cbstatus :=closed
9 k7 I' a5 i. q9 P8 m- f7 b  ENDIF


ENDIF

ENDEXEC
ENDMODEL

- e% ?5 B9 [% |7 ]
; V& D$ B- `! ^+ e0 x# u* x+ o

ziwaichengxiang

这个帖子不错，大家快来顶起来！

duym

@ZY2643181404