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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
) B2 n& @  x* o$ v+ Q8 @) K
; G( ?! K3 y% Y  K8 zMODEL breaker
  D! M' x! [& t$ j, T
# a* _2 g' J, _$ W+ t/ w6 G6 V* KCOMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

- H! L2 Q( Q3 w& n
! N6 K8 V* F6 y, M4 G# p. n5 I% BThe 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.
5 o  Y  G$ B+ v2 p- t. J& L
: L5 _2 T$ L+ R! ~- Z7 U& N% J+ [>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT

DATA
+ P8 D& _+ k4 U" l  topen                -- [s]   earliest elec opening
  ich1                 -- [A]   chopping level at power frequency
1 S3 G. t8 a  f9 f  slope                -- [A/s] current slope
- t- X" A/ Q, X/ G3 y9 J! y  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
  reignit_w {dflt:1e-2}              -- [s]  reignition window

* q: y$ N, A$ g( R0 c& h8 rINPUT
  v1   -- voltage across breaker terminals
  v2
  current  -- current through breaker
! M* C  V! z* O! S2 T6 a0 k6 ]) g
VAR
  state    -- 0= fully open
           -- 1= normal closed
           -- 2= tentatively open
" B# Q$ P; S( i' Q3 H1 c5 K8 i           -- 3= reignited
& k4 X1 x) F* ?0 r4 R$ x; k  first_opening  -- transition from 1 to 2
4 C' b% Q1 K  h  ?0 ^7 V' s+ V  reigniting     -- transition from 2 to 3
% l9 |/ M0 `& p) R$ ^  reopening      -- transition from 3 to 2
  fully_opening  -- transition from 2 to 0
  reopen_clock   -- state variable of 3
+ C2 u0 o7 T# H! s  elapsed        -- time elapsed since tbeg
' i3 o! n5 ~7 [  cbstatus       -- open/close signal to switch
  deltav  
( {6 |( }+ r9 j+ o  x: W3 o% p/ ]  ds
( O# ~0 b, p: k5 g5 G3 j  nds
  tdelta
  itopen
  ich2
8 L) f  Y/ V3 r2 H8 c
- y6 g4 e/ F1 `7 xOUTPUT
  cbstatus ,ds,nds
5 x2 L0 p+ u$ ~7 g4 T
HISTORY
3 W# K* }+ ^# l2 q) Z$ ^' z  current {dflt: 0}

INIT
9 n1 w% n7 ~/ L( k' S  state    :=1       -- normal closed
  cbstatus :=closed
( h' A7 m) v; h" ?  ds:=1
- ~; M4 d9 R4 P7 b' b( w1 K  nds:=-1
ENDINIT

; D/ A1 x1 \0 ?) uEXEC
+ d7 I" M9 F9 D; K, K/ D* h deltav:=v1-v2
7 i8 n) w( a# M% Q8 G* N. ]2 i IF t>=topen THEN              -- otherwise do nothing
) U  R, C/ U* o1 {) C: \8 y2 L  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
  ENDIF

  nds:=-ds
  itopen:=abs(current)
  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
4 Q! ?0 G- ?* j% t9 J
  first_opening :=false;  fully_opening :=false
  reopening     :=false;  reigniting    :=false

  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
! C1 {+ H$ a# L+ ]' D( t        OR current*prevval(current)<0)
   THEN first_opening :=true
" c, ^' |) d- [" |9 @        elapsed:=t-topen
) e/ Z" {  B9 Z3 a# _( W
  ELSIF state=2 THEN           -- was tentatively open   
8 C1 N9 V$ L5 L& B* |$ E
   IF elapsed>reignit_w  THEN
      fully_opening :=true
   ELSIF  abs(deltav) > ds
" ^) |, y3 N9 i" T( ~5 x      THEN reigniting:=true
   ENDIF

: f* O3 Y$ a, e  ELSIF state=3 THEN                  -- was reignited

   reopen_clock :=reopen_clock +timestep

   IF reopen_clock >=reopen_delay
+ J7 ]; y& q: m- V# l# |) n    AND (abs(current) <ich1
( o3 a4 D0 E* V9 u5 V$ m- U         OR current*prevval(current)<0)
! Q* X& N0 K; f    AND abs(deriv(current))<slope
        THEN reopening:=true
8 Z  K+ M, @: ]& r# C   ENDIF

  ELSIF state=0 THEN                 -- was fully open
, O2 Z' e5 l7 K4 |- I/ V  ENDIF

  IF first_opening OR reopening THEN
% ?, x! p  I% y   state :=2               -- becomes tentatively open

  }+ ^) A4 |6 W3 b  ELSIF reigniting THEN
   reopen_clock :=0
   state :=3                      -- becomes reignited
- h) r" {' P: u  ELSIF fully_opening THEN
" J3 o4 `' D' _  W* p   state :=0                     -- becomes fully open
, M9 t% G  M  |+ W* B% m& ^' [8 H  ENDIF
2 g# D" I$ `& Y, [( W) ]! Z4 k9 B
& t, W/ \9 F) n# w; ^& u2 s$ @  IF state = 2 OR state = 0             -- is now open
! L" M% h- _' e8 @  THEN cbstatus :=open
: }. [! x9 g* R# v4 _. j  ELSE cbstatus :=closed
  ENDIF
9 q( ]* |" p9 }
8 C; q/ N7 H8 l8 e1 H5 j
) ], q) C. s" ?; d ENDIF

: b2 T. O) e; M7 UENDEXEC
3 H! k) i- _: f) N8 E- bENDMODEL


$ O) ?3 r9 D$ k2 G

ziwaichengxiang

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

duym

@ZY2643181404  
7 [# t3 n& Q1 D9 m. M$ A% O7 h+ W