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

duym

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

MODEL breaker
" t. O; P) c* {# [" g$ l
3 j9 _( i, }& E& @. t, UCOMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

# D7 H' [) e; |
The model of this example  updates  the  present  state  of  the
circuit  breaker  by  taking   into   account   the   conditions
controlling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.
2 N+ K8 v( L8 \& |7 E5 Q  u
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT

+ L4 G4 |  z) ], ZDATA
  topen                -- [s]   earliest elec opening
( X1 ]& a$ J& }, C; _  ich1                 -- [A]   chopping level at power frequency
% w: o8 s' d! ?  slope                -- [A/s] current slope
0 s. M9 q' M' B  l9 b$ n) Z  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
  reignit_w {dflt:1e-2}              -- [s]  reignition window
" ?; V) n7 N1 G8 E
0 P7 H/ `8 t  y, ~* w; T* e: zINPUT
  v1   -- voltage across breaker terminals
; Q4 s# f- N; v5 C- b4 _2 t  v2
* g9 U0 O1 ]8 Z) L: K! Z  current  -- current through breaker
. J0 V# n( a0 B+ M; F( y. j
VAR
  state    -- 0= fully open
$ D5 W  t  t$ K" O: p6 A           -- 1= normal closed
$ C, W+ o' I" R4 a           -- 2= tentatively open
           -- 3= reignited
  first_opening  -- transition from 1 to 2
; K$ N" q2 b7 E: Z  N% L4 k, @' M  reigniting     -- transition from 2 to 3
  reopening      -- transition from 3 to 2
( s+ {, D, Y7 j: n, _9 U2 J! C$ B  fully_opening  -- transition from 2 to 0
# ~6 ]! j) ?. @5 f$ i- t  reopen_clock   -- state variable of 3
1 Q+ g. T/ J5 z0 P5 V2 k1 w  elapsed        -- time elapsed since tbeg
  cbstatus       -- open/close signal to switch
  deltav  
  ds
3 ?9 N" N( N6 A! p, S/ o+ c  nds
  tdelta
% o8 K+ @) c) G. j) m& _  itopen
  ich2
, h  |! M; V! R: N" w
" C) F# N- t% ]4 R  ZOUTPUT
3 K5 b5 m1 m4 }  cbstatus ,ds,nds

$ p4 t' v, q) B% V: r# u% P% \$ YHISTORY
' P& }4 Y5 }# c- m  k, T  current {dflt: 0}
4 @7 B9 b9 |; M3 o/ u5 h
; L. B, J; D" Q# e; x1 B2 L" m2 PINIT
  state    :=1       -- normal closed
  cbstatus :=closed
  ds:=1
, |( L9 \6 f- s* h1 ^  nds:=-1
ENDINIT
8 \$ v6 p  u9 n" ?
/ Q6 G9 c& [6 R/ }EXEC
: J) D8 b& E* Q2 | deltav:=v1-v2
" w' w, U% R" l IF t>=topen THEN              -- otherwise do nothing
  tdelta:=(t-topen)*1000
: E. z- M4 Y5 Z  IF tdelta<=0.44
    THEN ds:=(38.5*tdelta**1.5)*1000
& F- G$ Y% t9 J  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
! k% ~" w8 U) N3 a0 ~) _  F  ENDIF

' `7 D6 O- |7 |2 x6 O  nds:=-ds
  itopen:=abs(current)
: e2 y% w, t, f: k: G" {  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
" e  d- d, y( }6 l7 ?  f0 O+ a
6 Q4 E7 C$ O* G! `+ }- `  first_opening :=false;  fully_opening :=false
" y+ i# W8 S. v7 i# p  reopening     :=false;  reigniting    :=false
- ?& ?3 H. ]# W# z' r& x# I
3 y! n+ h1 P5 X1 L2 p& i; G  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
        OR current*prevval(current)<0)
4 u( g( e2 W% ^& a% G1 \   THEN first_opening :=true
        elapsed:=t-topen

: Q$ v+ H  |: r: B: N( k1 P/ ~  ELSIF state=2 THEN           -- was tentatively open   

& `6 W( U8 D" O4 N% V5 D" I; M   IF elapsed>reignit_w  THEN
      fully_opening :=true
   ELSIF  abs(deltav) > ds
      THEN reigniting:=true
% f& i3 o3 \1 V7 Q/ L5 j. F# |. u   ENDIF

# N; @( d$ n/ u  ELSIF state=3 THEN                  -- was reignited
+ ^! f3 ?2 B" P* m3 ?/ o2 F, o7 ~% M
5 C& Y% B! L4 D0 n8 z5 ^   reopen_clock :=reopen_clock +timestep
& T/ o5 f2 w. G$ T% R. }  I4 g9 [$ K
   IF reopen_clock >=reopen_delay
" S) s# [( d* Y# }' m    AND (abs(current) <ich1
9 n) j+ T8 W8 K         OR current*prevval(current)<0)
    AND abs(deriv(current))<slope
% Q$ w; I/ s: y$ T        THEN reopening:=true
   ENDIF
* g7 A3 B. z7 l- d" f" ?6 X
  ELSIF state=0 THEN                 -- was fully open
4 T$ _4 j2 K" z9 `( L  ENDIF
( ]- d/ z: d5 h# U
/ B; m" ^2 }* _! G" S  H8 u% c  IF first_opening OR reopening THEN
   state :=2               -- becomes tentatively open

  ELSIF reigniting THEN
   reopen_clock :=0
3 J+ M% @. t5 g  X% {5 l   state :=3                      -- becomes reignited
# H. s) r2 ~$ y2 f  ELSIF fully_opening THEN
   state :=0                     -- becomes fully open
  ENDIF
2 Y$ h; a, w2 V  Q
  IF state = 2 OR state = 0             -- is now open
* g5 G9 r% H3 V# }  THEN cbstatus :=open
( u, J2 r) O' X: @- u& L  ELSE cbstatus :=closed
3 q' m5 x$ a: ?8 d9 k9 s) N  ENDIF

6 S) `6 y2 ~. w: C8 K: a7 Y/ e% b% K
ENDIF
6 f& N. I/ \. N* V; ]$ i3 O3 m+ i
ENDEXEC
ENDMODEL

ziwaichengxiang

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

duym

@ZY2643181404