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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
" ]' m0 R( y6 u: N/ V0 ~
MODEL breaker
; b# \4 `) g6 B& }8 r" _+ L6 K
COMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
' R  D* A) N% l+ [, M+ f

  A" L7 b* W0 E* L( k, Q/ u/ ZThe 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,
5 U/ u: ~0 C1 D% ]' z6 b, x  Hreigniting, reopening, and full opening.

' c- e. n) Q* c2 V" w+ [8 w6 J4 r>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT
! d& S9 v3 ?4 O/ I
6 \0 J( _' N. sDATA
  topen                -- [s]   earliest elec opening
7 g% p# p1 Y& |; t, R# |9 Q  ich1                 -- [A]   chopping level at power frequency
/ h; O6 F( p/ b! M  slope                -- [A/s] current slope
  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
  reignit_w {dflt:1e-2}              -- [s]  reignition window

INPUT
  v1   -- voltage across breaker terminals
6 P- c/ F- r" ^  O  v2
  current  -- current through breaker
9 Q% T' F7 J) j1 x8 K: \1 y
VAR
! j" C0 j1 K0 u% `  state    -- 0= fully open
           -- 1= normal closed
9 e* @5 Z5 S# ^( Z8 \& D1 A+ h8 J+ ^# C           -- 2= tentatively open
           -- 3= reignited
  first_opening  -- transition from 1 to 2
# R1 p" J- C5 E* [- _+ p: ?" H( P( Q1 [  reigniting     -- transition from 2 to 3
4 J! b) j/ O  H, G9 r0 N8 L+ i9 v0 U  reopening      -- transition from 3 to 2
( W, I# m. Q3 k  fully_opening  -- transition from 2 to 0
4 S  C+ x2 o& l7 \+ r6 x. y  reopen_clock   -- state variable of 3
  elapsed        -- time elapsed since tbeg
5 {/ G4 f3 A. a. s  cbstatus       -- open/close signal to switch
  deltav  
  ds
( k. M6 q- z$ p4 G5 v  nds
) Y( b9 V. Q: q6 u  tdelta
  itopen
  ich2
5 E, ~6 i+ o6 C
4 R, e' h0 D) U" U" vOUTPUT
  cbstatus ,ds,nds
: e# J) V( g- |' s+ ^  W
" S" a4 B6 j* Q: Y1 y5 V% qHISTORY
; I& {) v; a6 x; p7 U! ?0 ~  current {dflt: 0}

INIT
  state    :=1       -- normal closed
  cbstatus :=closed
8 X6 H3 @  [7 s0 t4 {9 I9 [  ds:=1
  nds:=-1
" |, P. ?$ w4 p5 O7 |ENDINIT

; `- m' W  j3 _% MEXEC
6 D1 `+ C  x5 ~ deltav:=v1-v2
IF t>=topen THEN              -- otherwise do nothing
, N' M  b$ C  E8 N) g$ o  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
    THEN ds:=(38.5*tdelta**1.5)*1000
8 r* Y' |7 ^/ @4 q  ELSIF tdelta>0.44 AND tdelta<=5.00
1 @0 h. }, G: x. O    THEN ds:=55.1*tdelta**0.55*1000
: g" W# D; J3 D- y  ENDIF
) ^: X6 W% [' z9 n. }. V. E! Z1 i
1 F  G& c1 C$ |9 O  nds:=-ds
; m6 x: }" U- d# z/ g6 N& s+ w" X# A% A  itopen:=abs(current)
  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
; c' n4 R8 T+ L1 R
  first_opening :=false;  fully_opening :=false
" ^! r& V, X8 v1 x/ j9 e9 ^  reopening     :=false;  reigniting    :=false

  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
        OR current*prevval(current)<0)
   THEN first_opening :=true
        elapsed:=t-topen
' H! n/ |# X2 o4 |. e
) g, c4 u3 s$ K  `5 u  ELSIF state=2 THEN           -- was tentatively open   
% W, \  N& c: ~/ s7 Y6 G7 {( Z+ g
   IF elapsed>reignit_w  THEN
, r( M) }7 c, \      fully_opening :=true
% |9 C, Q! |9 J* B2 q5 o   ELSIF  abs(deltav) > ds
0 `$ M4 L7 Q# y" ~+ Z( o      THEN reigniting:=true
   ENDIF

  ELSIF state=3 THEN                  -- was reignited
4 q  Z( o1 P& ^- L0 N
   reopen_clock :=reopen_clock +timestep

   IF reopen_clock >=reopen_delay
; `/ n( k* }, {1 W* ~    AND (abs(current) <ich1
' v! `% O. m% b7 d6 H  f2 p; X         OR current*prevval(current)<0)
2 D. D1 t5 c0 b% h    AND abs(deriv(current))<slope
        THEN reopening:=true
   ENDIF
4 U" H; g2 ]" C, i+ \) d
  ELSIF state=0 THEN                 -- was fully open
8 q9 e$ z- l) Q0 I; ?  w  ENDIF
) _+ W& d& p8 x! l
  IF first_opening OR reopening THEN
   state :=2               -- becomes tentatively open

  ELSIF reigniting THEN
/ M; _  _9 [" N" e7 m& O   reopen_clock :=0
+ Z1 v+ ]' F5 Y   state :=3                      -- becomes reignited
  ELSIF fully_opening THEN
   state :=0                     -- becomes fully open
  ENDIF

  IF state = 2 OR state = 0             -- is now open
  THEN cbstatus :=open
* `  p3 I+ ~1 T2 U+ a$ k  ELSE cbstatus :=closed
! }; B( N3 V0 u  ENDIF


" I" }5 r2 w, {* t; _1 v ENDIF
" p6 j: _" b; P' k4 m
2 _! X4 ~7 A; P! M: OENDEXEC
/ B5 z2 K& c. VENDMODEL

$ Y3 q, q) |* |. a
; J* C% r4 K% O+ c+ `+ N) [

ziwaichengxiang

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

duym

@ZY2643181404