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

duym

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

MODEL breaker

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


0 e* l% Y- t  L( \! \( tThe model of this example  updates  the  present  state  of  the
circuit  breaker  by  taking   into   account   the   conditions
4 S: C- w; _; econtrolling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.

1 i: j$ g, M1 ?2 G7 d>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT

DATA
  \# c" E" P: |% K; p# w8 ?  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

2 K; x$ T' x! P: @INPUT
0 b4 \6 Z9 T; b  }+ p  v1   -- voltage across breaker terminals
  v2
  current  -- current through breaker
( W# @2 j3 \- |9 g
, j/ P+ ^' R/ _6 i' Y7 G$ u& D6 }VAR
/ [7 t- Y! t! {9 U- _2 M' c  V' J& v  state    -- 0= fully open
, S' U7 q: {" [9 \4 a0 M! h           -- 1= normal closed
8 r8 |# E, r2 a           -- 2= tentatively open
           -- 3= reignited
) [5 t$ K8 u3 t  first_opening  -- transition from 1 to 2
! X+ s& Y7 B1 Q7 ]- k' w  reigniting     -- transition from 2 to 3
$ R( S% n0 ^0 d  reopening      -- transition from 3 to 2
  fully_opening  -- transition from 2 to 0
! h; p% j- k1 T  reopen_clock   -- state variable of 3
  elapsed        -- time elapsed since tbeg
- b8 H/ S. q0 S- d- v0 ^  cbstatus       -- open/close signal to switch
2 W" l0 r; ~; v0 [  Q  deltav  
7 D) e4 V) e0 w+ _( O  M+ a  ds
5 M& b/ o, h- _7 S3 K  nds
  tdelta
  itopen
  ich2

" j, n6 ^, O. `) e7 Y7 d; d( YOUTPUT
  cbstatus ,ds,nds

HISTORY
$ `2 M+ x7 R  H0 _$ J( S  current {dflt: 0}
$ E9 ?& R) T7 f! p# I0 O
INIT
5 ~. v; _, x( m' e/ Y0 `' U  state    :=1       -- normal closed
  cbstatus :=closed
" u* w; E1 z& y% `7 \4 M  ds:=1
  nds:=-1
. J7 s+ o9 i: Q* F2 P- G9 xENDINIT
/ S6 C3 e- \! P1 e3 A
EXEC
6 U; P; M( J, p$ a) T deltav:=v1-v2
IF t>=topen THEN              -- otherwise do nothing
  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
    THEN ds:=(38.5*tdelta**1.5)*1000
. _5 q0 @$ }5 G; t  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
  ENDIF

  nds:=-ds
  itopen:=abs(current)
* R: S$ d- g; ?7 q  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))

1 u/ D& P+ d8 d* c7 ~4 W  first_opening :=false;  fully_opening :=false
  reopening     :=false;  reigniting    :=false
7 Z9 Q% B+ `& k, t2 o
6 m& h" Z( A! q8 s  `1 O8 w0 l6 i8 k  IF state=1                      -- was normal closed
# y" a  D/ k# C( n3 D   AND (abs(current) <ich1
        OR current*prevval(current)<0)
6 U; G& X7 t- i. q   THEN first_opening :=true
8 L4 z8 I/ _3 A8 v/ D! ?5 ~* [! x6 I        elapsed:=t-topen

  ELSIF state=2 THEN           -- was tentatively open   

   IF elapsed>reignit_w  THEN
$ e- P7 M; d- w) X      fully_opening :=true
   ELSIF  abs(deltav) > ds
: A. y* Q6 e" Q2 O# S& H      THEN reigniting:=true
   ENDIF
# x9 V7 H" |' w- X7 W
4 A! U& X# j4 Y& |5 R  ELSIF state=3 THEN                  -- was reignited
7 ^+ ~% t$ }* e7 s
9 f- Q* C$ ~3 V  w7 j% ^   reopen_clock :=reopen_clock +timestep
! ^$ K. f! }7 o6 J. G0 b7 c  O; a
; F: g- O& F  h! F1 S   IF reopen_clock >=reopen_delay
    AND (abs(current) <ich1
         OR current*prevval(current)<0)
    AND abs(deriv(current))<slope
        THEN reopening:=true
9 @% O' X- T4 L# v# u! j+ @* G   ENDIF

# M5 L0 G8 a  A. E  A  ELSIF state=0 THEN                 -- was fully open
7 i7 H* _0 B, [4 R, V  ENDIF

  IF first_opening OR reopening THEN
   state :=2               -- becomes tentatively open

- r2 @! m5 g2 j2 B: ^! ]  ELSIF reigniting THEN
5 y7 W, A( h0 M) j; S   reopen_clock :=0
7 ^5 ?' T# C6 c& n! m* U! F   state :=3                      -- becomes reignited
5 o& C% ?/ w, Y" Z9 y  ELSIF fully_opening THEN
   state :=0                     -- becomes fully open
8 b& `' S% J3 J  ENDIF

: K" r1 U+ |/ W  IF state = 2 OR state = 0             -- is now open
" ~6 |% [! ~2 H# P$ m  THEN cbstatus :=open
  ELSE cbstatus :=closed
5 l5 N9 n6 B( o) T  ENDIF
$ g* t2 X0 W7 k6 y4 Z/ d$ a4 |
9 e% B/ u* F! j3 n7 y3 V
ENDIF

ENDEXEC
6 [) u8 \0 P+ x/ H, v, r" sENDMODEL

. a  s* p6 o9 H
. H: e; `  \/ D

ziwaichengxiang

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

duym

@ZY2643181404