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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。

; w3 `1 x$ ^( ]" y6 X7 ^7 V- |MODEL breaker
9 A% m& l, y( Y9 D+ Z% G/ a
6 B; _7 ?) e% Q) D% l6 QCOMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

! q' j1 {6 [" \2 K$ i/ k8 i) y) v' D
0 Q/ G) I+ t! K; }2 C" tThe model of this example  updates  the  present  state  of  the
0 i( P9 z. j6 U2 v" Y+ s5 x# ]circuit  breaker  by  taking   into   account   the   conditions
controlling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.
$ N% X% w8 `/ R. [3 o$ l
7 G5 \8 l, t$ p/ G2 z>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT
: M+ E% M# |! @; Z3 g+ y
DATA
  topen                -- [s]   earliest elec opening
  ich1                 -- [A]   chopping level at power frequency
; g  O1 Q- {: e5 }, {* J; ]  o  slope                -- [A/s] current slope
+ l; P+ ^  [5 |1 Q' B  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
( y& x$ y% q: X9 K  reignit_w {dflt:1e-2}              -- [s]  reignition window

INPUT
4 j3 f* s6 u9 c8 @9 a+ [& B  v1   -- voltage across breaker terminals
  v2
  current  -- current through breaker
2 g$ L! ?+ h: H2 q1 v
VAR
8 ~1 E( J: \! Y  state    -- 0= fully open
           -- 1= normal closed
           -- 2= tentatively open
           -- 3= reignited
2 m; E) ^" S2 O- m1 J/ ]  first_opening  -- transition from 1 to 2
2 l1 [7 w% U' s: Z, ]2 x. a  reigniting     -- transition from 2 to 3
  reopening      -- transition from 3 to 2
  fully_opening  -- transition from 2 to 0
  reopen_clock   -- state variable of 3
" w0 o" T' u) s1 D' H  elapsed        -- time elapsed since tbeg
  cbstatus       -- open/close signal to switch
; z" o' u( u5 s/ Q3 N% u  deltav  
  ds
  nds
  tdelta
; C9 M9 _6 ^: ~6 Z$ H: X5 G  itopen
4 ^/ l3 i9 f- f  ich2

OUTPUT
  cbstatus ,ds,nds

; Z2 H- `# u+ Q% y- DHISTORY
8 @" J4 r8 _" o7 L  current {dflt: 0}
! |8 Y4 c2 }7 q9 M- R4 v
* l- x. v) f0 A2 F6 L* vINIT
# w% a+ h1 [8 p  \  state    :=1       -- normal closed
  cbstatus :=closed
  ds:=1
0 S) d4 T5 G! x2 g6 k$ r  nds:=-1
ENDINIT

! a/ Y- f: r. y* AEXEC
deltav:=v1-v2
IF t>=topen THEN              -- otherwise do nothing
  tdelta:=(t-topen)*1000
  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
' @* j) Q) d5 o* A  ENDIF

  nds:=-ds
( ~' ~5 N- ]! X1 z7 T8 O  itopen:=abs(current)
2 O% E, L; T: `* u% ]3 v  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
, _; u% ~* c1 G9 H3 k
  i1 H5 J1 v  f4 B: t* N  first_opening :=false;  fully_opening :=false
+ t7 R& R' f' i6 a& N  reopening     :=false;  reigniting    :=false

  IF state=1                      -- was normal closed
8 z% g( D* E  W- j* f! q   AND (abs(current) <ich1
0 E: e" \% l* |( T" `5 Z- G        OR current*prevval(current)<0)
$ x# d& |9 a% S5 f+ {( x( W( n+ o   THEN first_opening :=true
        elapsed:=t-topen

  ELSIF state=2 THEN           -- was tentatively open   

   IF elapsed>reignit_w  THEN
' Z  s0 ]; ?0 z9 }$ u( r( f3 Z      fully_opening :=true
! K4 V( X% j7 N   ELSIF  abs(deltav) > ds
4 `2 P0 X5 m  O* E      THEN reigniting:=true
   ENDIF

4 @, g- k) d( Z0 v5 J4 ?7 H  ELSIF state=3 THEN                  -- was reignited

& R, A2 ?$ o0 K0 F0 i( q   reopen_clock :=reopen_clock +timestep
  j# }; a, P2 O9 k. b/ z/ N
% q5 k9 O5 O! t; u5 |   IF reopen_clock >=reopen_delay
" [, Z; A6 H- }9 e' |, O    AND (abs(current) <ich1
         OR current*prevval(current)<0)
$ g- y; d; }. a# W: S8 b4 R    AND abs(deriv(current))<slope
5 ~. i3 n, V& v; t8 [' ~6 O; H2 d        THEN reopening:=true
. U1 s2 \0 k* t5 U& }% L+ e   ENDIF
/ W, Z& O, E: i, e2 |
  ELSIF state=0 THEN                 -- was fully open
  ENDIF
, ^" F2 _( n, M# D* Q( e5 B
% ^  y1 N# Z) k5 c6 b  IF first_opening OR reopening THEN
3 r9 J- d# v4 t   state :=2               -- becomes tentatively open

3 C' A1 s# t- X& O, o* Z% N  ELSIF reigniting THEN
) W/ c/ o* ?4 k; H   reopen_clock :=0
0 ]  r! T5 D% F- _' M( X4 O* x2 O  [   state :=3                      -- becomes reignited
! r0 q2 K+ V& h9 `8 `; Z) _  ELSIF fully_opening THEN
   state :=0                     -- becomes fully open
  ENDIF
) w! y4 s+ o) ~) T. u) B& Q
2 t+ ]4 H- S- L2 W& O  IF state = 2 OR state = 0             -- is now open
  THEN cbstatus :=open
( o5 \# \* C, m2 R  ELSE cbstatus :=closed
& @/ ^9 I4 m+ E' r" Y& E! W  ENDIF
) a! Q- }8 p% F  Q3 R
# n2 \1 i! c, o2 ~: k
ENDIF
4 e# N8 t* e5 w* y; p, F
ENDEXEC
ENDMODEL
6 q' X0 e5 |6 S

ziwaichengxiang

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

duym

@ZY2643181404