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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
5 @2 R7 ~2 ~9 @6 ]* a
MODEL breaker

COMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

4 o! V3 k% Z, s+ _
. S4 d) o) j4 g8 w" u. mThe model of this example  updates  the  present  state  of  the
circuit  breaker  by  taking   into   account   the   conditions
% B6 {* Q7 g, \1 f/ }, r% ncontrolling four represented state transitions:  first  opening,
, B: e* t8 r) z1 |0 {reigniting, reopening, and full opening.

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

/ A- _, g# N3 E, [( d/ BDATA
  topen                -- [s]   earliest elec opening
7 }2 O" h/ k9 l/ n1 A6 i( C  ich1                 -- [A]   chopping level at power frequency
: a  w& Y& v& D5 V  slope                -- [A/s] current slope
  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
  reignit_w {dflt:1e-2}              -- [s]  reignition window
" T4 N; @$ a3 ?% O8 y9 j: W
3 n& {8 k' H  d7 lINPUT
  v1   -- voltage across breaker terminals
" Q8 h) X3 z' Q  v2
6 ~1 L8 p7 D5 V% f! H4 [% F0 Z6 r; E  current  -- current through breaker
  c8 ]8 H' h. H7 }( ~
VAR
& x- `$ d4 P3 _  state    -- 0= fully open
           -- 1= normal closed
. w! {8 E8 ?. j* z# I/ [           -- 2= tentatively open
% D* F1 H+ |2 s' x4 f  H, S           -- 3= reignited
  first_opening  -- transition from 1 to 2
  reigniting     -- transition from 2 to 3
6 l6 U. E7 v2 U0 @) r) ]9 q- L  reopening      -- transition from 3 to 2
! ]5 Y5 v( P( d- ^/ c/ I# L  fully_opening  -- transition from 2 to 0
0 H' |2 [: G3 ~4 b: Z& U# N+ M  reopen_clock   -- state variable of 3
  elapsed        -- time elapsed since tbeg
  cbstatus       -- open/close signal to switch
  deltav  
  r8 {7 o- ^4 ]8 p. H  ds
  nds
  tdelta
  itopen
, T" o/ y: E& R; R* F$ f/ o! W  ich2
( n# u& ?% N$ s- i5 p& t
) w5 b* c* }1 R2 DOUTPUT
1 v$ l9 I; H  T2 f  cbstatus ,ds,nds

! r" g- d6 M8 k. R' H1 CHISTORY
. X. Q# z, t4 O7 @8 r3 T6 a  current {dflt: 0}

INIT
  state    :=1       -- normal closed
  cbstatus :=closed
+ P2 s0 ^  l! N: k% J! E; X  ds:=1
3 R0 m  t2 d( f" q5 l  nds:=-1
ENDINIT

# O, E4 J7 E( ^) ~EXEC
deltav:=v1-v2
IF t>=topen THEN              -- otherwise do nothing
  tdelta:=(t-topen)*1000
: E# X5 l$ ^! O+ q) `  IF tdelta<=0.44
6 b4 D3 b/ o: W: h- c, o3 L  ?+ {    THEN ds:=(38.5*tdelta**1.5)*1000
  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
- B1 Z! B; a/ @8 C7 S( P  ENDIF
, g0 F( k' Q; Q8 B( P3 a* N
  nds:=-ds
' N2 Y* f3 s9 y% y1 U: D6 L- Y" N/ g  itopen:=abs(current)
  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))

1 R$ v) z5 k% _  first_opening :=false;  fully_opening :=false
( p9 s& A0 V2 v2 Y& T- |  reopening     :=false;  reigniting    :=false
6 n4 w" r* D3 K# p$ j2 ?" H
2 Y, N- q+ ]$ f% f# k  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
. K+ H, _; p0 J$ [        OR current*prevval(current)<0)
7 e9 W' s: F/ V& u7 V) i2 g8 h   THEN first_opening :=true
2 g7 N' `# c2 u% j        elapsed:=t-topen
+ l1 U, o. z6 Z; Q
1 q3 e8 g- x8 S  ELSIF state=2 THEN           -- was tentatively open   
8 r  l0 f; a& e# X
4 M: {' S4 Y+ L5 w  S   IF elapsed>reignit_w  THEN
  M; Q8 `$ [& r7 P3 j9 n4 w      fully_opening :=true
   ELSIF  abs(deltav) > ds
      THEN reigniting:=true
0 L- @: a6 P7 T& I   ENDIF

  ELSIF state=3 THEN                  -- was reignited

6 H, `- [: w9 f0 ?- Y) z+ H* O   reopen_clock :=reopen_clock +timestep

( w0 ^( ^0 l0 I  R# R& n   IF reopen_clock >=reopen_delay
    AND (abs(current) <ich1
         OR current*prevval(current)<0)
    AND abs(deriv(current))<slope
        THEN reopening:=true
+ m# a. t+ F/ V) t3 I6 y- q) m   ENDIF
- @; B4 v3 H) R$ g; B% |2 T
  d" C$ h3 ?3 h6 O1 A  ELSIF state=0 THEN                 -- was fully open
, F$ N* s( Q3 v1 U/ i  ENDIF

: @' {- x& ^# k; R# z  IF first_opening OR reopening THEN
: J0 e  W0 o) {9 r0 z   state :=2               -- becomes tentatively open
8 O; Z# S- @* x* S8 ^, h+ I: _- P
  ELSIF reigniting THEN
, H9 s' F. R* r, C2 I+ T   reopen_clock :=0
* e! }# Y# m+ N0 }# v% k  E8 F   state :=3                      -- becomes reignited
& I$ E7 A, q: T+ V( _  ELSIF fully_opening THEN
   state :=0                     -- becomes fully open
, I9 A# S" y: D% I: P0 l# z2 e& U* b! ~  ENDIF

! z/ |, i! @" H  IF state = 2 OR state = 0             -- is now open
  THEN cbstatus :=open
  ELSE cbstatus :=closed
  ENDIF
7 I6 M- V0 q& Y/ ~9 C, o
+ f" h3 E; h2 \# W" g! n6 K7 B
9 V5 V- M# I/ m% z- a* S4 m) m ENDIF

! T4 ]- R) Q& n# hENDEXEC
ENDMODEL
$ l$ a: P$ l! @2 s

, w& M7 X- d  Q7 e3 I9 b* R% L5 D

ziwaichengxiang

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

duym

@ZY2643181404