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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
# ^: x; B: H0 k7 t$ A. [
4 M. _# D8 J2 q8 {5 f( C2 [* dMODEL breaker

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

" ?2 r" [$ V$ Q! x) [9 r
The model of this example  updates  the  present  state  of  the
, \5 y) T9 \- Q$ s# ~circuit  breaker  by  taking   into   account   the   conditions
( Y2 |# p+ H7 Mcontrolling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.

; z' V( S9 C7 m5 [: A2 }' [( J" K>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT
, i/ ?) G1 ?: A" y1 o
DATA
8 H: ^- m% v: P8 O% b1 a/ V" C  topen                -- [s]   earliest elec opening
# z0 j5 u1 D3 T9 N) r0 c/ p  ich1                 -- [A]   chopping level at power frequency
  slope                -- [A/s] current slope
( e$ Y0 ~* l# v5 M8 R) l  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
- r1 b( @3 v% o+ J3 {  reignit_w {dflt:1e-2}              -- [s]  reignition window
& F  `/ {  a, \& }% W4 l/ X
INPUT
  v1   -- voltage across breaker terminals
  v2
  current  -- current through breaker
2 Y7 L" Z( E) D4 R
+ r( U3 V5 a6 h: t' qVAR
  state    -- 0= fully open
           -- 1= normal closed
# H* O+ r- i  c, _5 J  t2 E           -- 2= tentatively open
# B: L: U* f4 B/ r           -- 3= reignited
$ e% h4 j7 I( H0 t0 d8 {  first_opening  -- transition from 1 to 2
  reigniting     -- transition from 2 to 3
  reopening      -- transition from 3 to 2
( `: E/ a% b( W0 y) Y' U# c* y" {  fully_opening  -- transition from 2 to 0
  reopen_clock   -- state variable of 3
  O' R) S) b7 s" i: R+ S2 U  elapsed        -- time elapsed since tbeg
  cbstatus       -- open/close signal to switch
  deltav  
  ds
  nds
  tdelta
$ y: D3 I( Y& R) [) h8 Z  itopen
6 R9 n. j! t5 e; v& W: O0 |  ich2
6 _8 T  z, O0 |" S7 k2 [
OUTPUT
  cbstatus ,ds,nds

HISTORY
# s  ?" V) f: a: D! r' `  current {dflt: 0}

! F4 r$ F/ J# f; jINIT
  state    :=1       -- normal closed
: G, B1 a1 `9 w- l9 `4 M" _* x  cbstatus :=closed
5 Y) a% g) n; W# q  ds:=1
  nds:=-1
ENDINIT
8 q! Z7 K$ o  \# l% T" o
EXEC
) u% l% M  g- q8 P4 y3 M deltav:=v1-v2
3 j4 Q* ]; n  ` IF t>=topen THEN              -- otherwise do nothing
* @$ g. V; V0 Z1 J0 s8 o  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
8 s4 Y3 p# L' _1 V* @% q' U! k0 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
  ENDIF

  nds:=-ds
- m4 X8 ^) @# h+ S7 B  O  itopen:=abs(current)
' j' _# c' o" M$ A  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
2 {! h4 |' X# D% v6 n+ ?
7 g# ^7 C$ s- }5 h3 d  first_opening :=false;  fully_opening :=false
  reopening     :=false;  reigniting    :=false
5 \: b7 Z* W: \
  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
. ?9 m2 U* s4 k9 e+ X1 m        OR current*prevval(current)<0)
   THEN first_opening :=true
0 ~5 e" ]6 d" v; u7 Z: ?        elapsed:=t-topen

  ELSIF state=2 THEN           -- was tentatively open   
8 W2 C) H/ `9 K# ^
   IF elapsed>reignit_w  THEN
8 N  n7 M* P( O# t; y( ~& m      fully_opening :=true
# N- D$ o; ^: c   ELSIF  abs(deltav) > ds
5 Q6 b% _4 G  @9 {  E' R      THEN reigniting:=true
6 b4 T( I. J+ V( K; y3 A; G   ENDIF
  [$ a8 z. j4 V) Y
  ELSIF state=3 THEN                  -- was reignited

* q, M9 d; t' b' c2 X9 q; s1 T   reopen_clock :=reopen_clock +timestep
4 r7 C+ j' h. V* Z5 @
4 P0 P5 i0 o2 b! @9 o5 N; d   IF reopen_clock >=reopen_delay
6 _2 h5 c1 ^! m6 T# V    AND (abs(current) <ich1
         OR current*prevval(current)<0)
    AND abs(deriv(current))<slope
: l6 x& C& y" x$ f8 d        THEN reopening:=true
   ENDIF
  ^- }& J% I$ Q3 i2 r
$ Y5 X2 S% P* J5 u/ M  ELSIF state=0 THEN                 -- was fully open
9 c" c' X# S7 Z- P  ENDIF
0 S3 r" F" L0 ?; t! Q! t( p4 o+ z
7 J9 V' w3 U8 c$ h  IF first_opening OR reopening THEN
   state :=2               -- becomes tentatively open
* R% ^. g8 W0 e. z
  ELSIF reigniting THEN
3 C# l5 t+ c* C* d& p! L   reopen_clock :=0
2 z- @3 C3 [: }) q2 P; ^   state :=3                      -- becomes reignited
  ELSIF fully_opening THEN
2 h/ S# u2 a2 X/ z2 T   state :=0                     -- becomes fully open
  ENDIF
& G; n' F7 l4 a) ]# F3 A; C
  IF state = 2 OR state = 0             -- is now open
5 s# ?' f6 |; G3 X+ Y  THEN cbstatus :=open
5 K1 R3 z8 g, J+ g' ^  ELSE cbstatus :=closed
  ENDIF
4 z6 T: U  y- Z2 n. ]6 x
, G, K0 Z; f. q- d, I+ Y
ENDIF
: Q! a- o. m8 D& a  @
ENDEXEC
1 X; C/ X; G& j$ B4 XENDMODEL

ziwaichengxiang

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

duym

@ZY2643181404