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

duym

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

MODEL breaker

1 U" V2 Y9 y+ [) W! SCOMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>


: M7 k% F/ k% B5 [) x) sThe model of this example  updates  the  present  state  of  the
circuit  breaker  by  taking   into   account   the   conditions
7 A0 p1 M; A% Dcontrolling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT

/ N+ N& P' |) R7 xDATA
  topen                -- [s]   earliest elec opening
  ich1                 -- [A]   chopping level at power frequency
  slope                -- [A/s] current slope
  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
  reignit_w {dflt:1e-2}              -- [s]  reignition window

, l  Z. z2 @) ?2 D2 W: j4 }INPUT
  v1   -- voltage across breaker terminals
  v2
  current  -- current through breaker

. n$ o! @/ `* k, M6 M. xVAR
  state    -- 0= fully open
5 I* {0 I# f! _+ Z4 a           -- 1= normal closed
           -- 2= tentatively open
           -- 3= reignited
: U/ L4 w" ]5 V8 R  first_opening  -- transition from 1 to 2
  reigniting     -- transition from 2 to 3
2 c8 H0 E8 r, O: }, M/ U% I- L  reopening      -- transition from 3 to 2
  fully_opening  -- transition from 2 to 0
  reopen_clock   -- state variable of 3
. a: ~2 w; V" N  elapsed        -- time elapsed since tbeg
3 b$ G1 {8 V2 G1 k% S+ O7 A  cbstatus       -- open/close signal to switch
+ V4 I% j7 B# u  q4 H6 u  deltav  
% _+ E- M& o5 z  x% K- q, ]. B  ds
  nds
  tdelta
- `/ U. M5 N5 {$ n/ D/ m  itopen
  ich2

OUTPUT
  cbstatus ,ds,nds

9 O- J/ M( w  r+ dHISTORY
$ h  K4 k# ~  B8 f/ B7 q5 ?5 U* w  current {dflt: 0}
. |3 }& s. K# R3 J, Z9 f& H: L
" x5 }6 ^% r5 M: bINIT
  state    :=1       -- normal closed
: I# J3 L6 r0 R* W# ?  cbstatus :=closed
" i# c  p4 \. T% I( k- z3 i: U  ds:=1
  nds:=-1
+ N% Q" ^5 _) B1 Z: [9 M3 w: {7 W6 qENDINIT

EXEC
deltav:=v1-v2
IF t>=topen THEN              -- otherwise do nothing
% f6 x+ Y! _. h' E0 ?6 M) E  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
    THEN ds:=(38.5*tdelta**1.5)*1000
$ |0 y6 ?2 d# y3 |  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
  ENDIF

  nds:=-ds
& i! t- `2 n- S* k  itopen:=abs(current)
  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
3 J4 ]6 ?$ t8 V2 y7 z
  first_opening :=false;  fully_opening :=false
1 P8 j3 o* R0 N' m  reopening     :=false;  reigniting    :=false

/ m: F2 J! ~( t9 X# K  ^/ ]  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
        OR current*prevval(current)<0)
! L% k7 Y6 y* a6 h' \   THEN first_opening :=true
        elapsed:=t-topen
. D4 s" \3 A+ |5 p2 z& d
, l% S: ~( x( U8 A  ELSIF state=2 THEN           -- was tentatively open   
; `! D, @# ]$ ~* L! U
- h' A/ J6 ^0 [   IF elapsed>reignit_w  THEN
9 \2 ?' k# f" a# ~0 ]2 d' |1 F6 C/ c7 Q      fully_opening :=true
& ]5 \: [. ~2 L3 k  t# x: B( _" n   ELSIF  abs(deltav) > ds
0 N8 ]- j" c6 k( F      THEN reigniting:=true
   ENDIF
0 C) L" j& j! f! H
  ELSIF state=3 THEN                  -- was reignited

   reopen_clock :=reopen_clock +timestep
9 p) x5 J; A) Q
   IF reopen_clock >=reopen_delay
    AND (abs(current) <ich1
         OR current*prevval(current)<0)
& P: y1 l& I5 W0 N0 k/ i. I    AND abs(deriv(current))<slope
        THEN reopening:=true
   ENDIF

  ELSIF state=0 THEN                 -- was fully open
  ENDIF

; q2 y$ I1 p6 r! g1 }  IF first_opening OR reopening THEN
   state :=2               -- becomes tentatively open

  ELSIF reigniting THEN
5 P2 @8 l) N  G* S( A+ j+ C& r   reopen_clock :=0
, j6 ]3 h; y( d3 ]' a" Q( w   state :=3                      -- becomes reignited
  ELSIF fully_opening THEN
+ A/ U& }3 D& X/ i8 A  r   state :=0                     -- becomes fully open
  ENDIF
( R4 d3 E, l7 m
  IF state = 2 OR state = 0             -- is now open
& z3 k) w) q7 b  THEN cbstatus :=open
& M6 A2 v6 e* N  ELSE cbstatus :=closed
  ENDIF
* L; w6 p: k9 Z5 B' ^; w+ e  ]

ENDIF
( n, l: Y6 Z4 P0 E1 b
ENDEXEC
5 J6 ]# r" K# l, i) {ENDMODEL

ziwaichengxiang

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

duym

@ZY2643181404