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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
) e8 a7 t6 I0 Y9 ]
, U! e& ]: ]4 K8 BMODEL breaker

COMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
8 [8 Z2 }) m  w- p. ]

The model of this example  updates  the  present  state  of  the
2 ]: \, K0 D$ n4 ^2 r* ~circuit  breaker  by  taking   into   account   the   conditions
  G  x" _+ |% O, l) c4 C5 Xcontrolling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.
3 e# u1 b+ t8 s: u
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT

& ]: a! v, D- DDATA
3 ?7 a+ ]& w  F% V  w) \' q4 _  topen                -- [s]   earliest elec opening
  ich1                 -- [A]   chopping level at power frequency
. M; D# j& L2 k' s. {: E  slope                -- [A/s] current slope
  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
" h7 @* H' ~2 m2 F  reignit_w {dflt:1e-2}              -- [s]  reignition window
7 S. \( ?& }0 ~& p+ A& x# M) F$ o
INPUT
+ ^3 i! p$ v, N: Y( t3 u  h  v1   -- voltage across breaker terminals
  v2
6 Z. Q3 b1 c/ O  K! Z/ |8 Q  current  -- current through breaker

# V) Z6 {/ K' ~- f3 I, s, {VAR
+ Z0 ^' P6 i$ |# M+ @  state    -- 0= fully open
           -- 1= normal closed
           -- 2= tentatively open
7 y) D. G' U/ Z           -- 3= reignited
  first_opening  -- transition from 1 to 2
  reigniting     -- transition from 2 to 3
  reopening      -- transition from 3 to 2
  fully_opening  -- transition from 2 to 0
! B: J. ?5 x3 |" Q3 ^* E  _; x  reopen_clock   -- state variable of 3
- k" f# n7 c, O, [- l: [2 B% C  elapsed        -- time elapsed since tbeg
7 I5 B& Y* s0 E5 R8 \5 n% k  cbstatus       -- open/close signal to switch
3 k) X! o3 D' t1 t' K8 P  deltav  
  ds
% k/ V: X, I  }  nds
8 ~1 w# _; C" j1 i5 ^  tdelta
0 [. \7 N4 t6 E) Q; U. j  itopen
  ich2
' z# a3 K# `% j% |: c) w
OUTPUT
6 s( Q, _  h# C, ?/ |# |% c" {  cbstatus ,ds,nds
& j& B. B1 D7 G, `  d6 v
HISTORY
0 \& p) n# o- l. Z0 L, S  current {dflt: 0}

: e3 t6 v& M4 e5 b5 g, [INIT
! P( P% Q+ }# }3 h5 I. c  state    :=1       -- normal closed
  cbstatus :=closed
4 K# p$ u4 L1 L9 e. L8 S( K4 ~  ds:=1
  nds:=-1
ENDINIT
/ h  G* z  `4 d( v' l, }
EXEC
  {. q/ l" x4 i" Q3 ` deltav:=v1-v2
IF t>=topen THEN              -- otherwise do nothing
# o' Z2 _  e' {/ M  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
+ L' X& r1 ^; Q( ^    THEN ds:=(38.5*tdelta**1.5)*1000
  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
  ENDIF

# C" ^( k3 r, z9 Q  nds:=-ds
  itopen:=abs(current)
  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
+ r) k4 [6 Q1 L. J. Z
5 v1 d- F, P7 L  first_opening :=false;  fully_opening :=false
  reopening     :=false;  reigniting    :=false

  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
        OR current*prevval(current)<0)
   THEN first_opening :=true
' o7 u+ j" t' M8 T. ~        elapsed:=t-topen
) d" V! |% M) \7 h; T
$ l, j) ?* z* G0 P% S' O  ELSIF state=2 THEN           -- was tentatively open   
8 w% ~' n" P8 o5 Y0 }* G, v$ T
# d" i5 O) c$ z5 |3 W" w; J& z   IF elapsed>reignit_w  THEN
1 k5 x0 S% w& Q" U      fully_opening :=true
& D& |; {6 Y6 u   ELSIF  abs(deltav) > ds
      THEN reigniting:=true
   ENDIF

  ELSIF state=3 THEN                  -- was reignited

6 J; }- y1 E. }# U6 k   reopen_clock :=reopen_clock +timestep
% A- w* q; @' H/ Q$ x$ q/ k
: J/ S+ Z5 |) c# o; g$ J   IF reopen_clock >=reopen_delay
1 w# }& B3 a: F    AND (abs(current) <ich1
0 F! F6 \9 p0 ?/ y/ x         OR current*prevval(current)<0)
    AND abs(deriv(current))<slope
        THEN reopening:=true
   ENDIF

. x, I# e' C) Y9 v  ELSIF state=0 THEN                 -- was fully open
  ENDIF

. P' q3 B$ A1 P* d/ G  IF first_opening OR reopening THEN
4 i- a3 {$ H% V+ m# R1 q7 ]/ Z   state :=2               -- becomes tentatively open
" k9 n6 a0 C# I* ?& x6 u
4 @1 d  ^$ t+ p' g/ Y' E! w  ELSIF reigniting THEN
   reopen_clock :=0
   state :=3                      -- becomes reignited
/ `% m- a, F) S8 W  ELSIF fully_opening THEN
# Y0 S2 J7 U* L- Z   state :=0                     -- becomes fully open
# G4 G, A) G, A% b  ENDIF
" M8 V: E' {9 {  h9 a. p* ^# ^
, C; [) }1 l1 H7 {9 r+ I$ S1 S  IF state = 2 OR state = 0             -- is now open
  THEN cbstatus :=open
# r( I, c$ |$ e1 n* W  ELSE cbstatus :=closed
  ENDIF
$ N7 \, z* }- c' d  t  ^
( o$ h* ]7 p! p  E
ENDIF
$ d! m2 r6 b# b* k
) b$ I, e$ ]- S8 R/ T$ ~ENDEXEC
4 @* t' `% }0 D& }) H& yENDMODEL
" @: Z1 {$ S6 {- d
" I$ K3 H+ ~, `
( Q1 C+ m; Z' f* N3 S6 L$ Q

ziwaichengxiang

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

duym

@ZY2643181404  
; e: ?+ B- m9 k% F- _1 @. `