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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
2 I, k8 ?- _6 Y( {- _3 r' j+ {
2 s  o0 M  k: q! R$ ~/ s/ V/ QMODEL breaker
0 l+ `' B& N  V) `% j! n
COMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
% ?; {7 d, U. p( _* {1 l

The model of this example  updates  the  present  state  of  the
  ?+ C- D% G- |( r% icircuit  breaker  by  taking   into   account   the   conditions
controlling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT
+ C# l+ c/ N6 s! e9 C
; @, P+ Q- B# I4 V/ E. s( dDATA
  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

INPUT
  v1   -- voltage across breaker terminals
/ m2 `) \6 }" `) W  v2
0 a; ]/ I3 f2 F" u8 L6 D. q  current  -- current through breaker

VAR
  state    -- 0= fully open
           -- 1= normal closed
  o% q# j# W0 {6 D& {5 j           -- 2= tentatively open
           -- 3= reignited
  first_opening  -- transition from 1 to 2
, T/ M: f( b9 f- }9 d  reigniting     -- transition from 2 to 3
5 B+ G+ _9 Y) Z! v) P' v% G  reopening      -- transition from 3 to 2
6 U1 ]7 o- j3 i- K1 e9 k, ?  fully_opening  -- transition from 2 to 0
. m& A0 R; `' K; Q; ?  reopen_clock   -- state variable of 3
- n, _2 E/ h) ]( A' u  elapsed        -- time elapsed since tbeg
" R% b/ ^- t+ z# j2 j9 E  c5 u  cbstatus       -- open/close signal to switch
) X, x, o- _- f9 k6 _  deltav  
" ~' z% _3 u0 g$ Z% ]  ds
; l9 f- V) P9 j+ [  nds
  tdelta
4 P/ L# G, I6 J  itopen
6 Y7 U) O+ b* n# l  ich2
. }, z* D) B8 V. k
0 b. W* }# _5 e) s9 i$ SOUTPUT
( i5 f8 P! _8 S5 s! W  cbstatus ,ds,nds
# o) J* q% A) X
4 c4 g0 p& S1 v' V  mHISTORY
7 |/ t' X1 [1 t/ d  current {dflt: 0}
1 h8 b4 P, T5 v, w* t0 q7 r
: F$ p5 [: h  j, R$ \* m$ s3 UINIT
' g+ X; O3 m9 U% n  state    :=1       -- normal closed
  cbstatus :=closed
2 ]2 ]' ~- L9 U" _& w% [, P  ds:=1
. n5 _$ X, D" v5 |  nds:=-1
ENDINIT
/ M5 x  o5 Y+ y! X. Q/ E4 R' {( v- z
& {. m+ g- D- [# fEXEC
8 n; x& X1 K0 ]7 M) j. S deltav:=v1-v2
* t: E9 D! g  [, \ IF t>=topen THEN              -- otherwise do nothing
  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
    THEN ds:=(38.5*tdelta**1.5)*1000
3 X1 v! a, _, [* I7 ?3 `) }9 n  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
1 E, K5 D0 ]/ b) c# K7 p# z7 P  ENDIF
/ H' X* v# M8 l, N' }1 y
  nds:=-ds
" m2 F1 c; N/ q/ @  itopen:=abs(current)
  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))
! Q% e* s: _6 @1 s) l/ Y  R) X; }/ v' C( c
5 b& B/ i, U- `1 |. X  first_opening :=false;  fully_opening :=false
6 w3 g( g0 t! B6 B# _  reopening     :=false;  reigniting    :=false

9 t6 q2 l+ K. i$ @; B2 @1 U! h( ^+ ~  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
        OR current*prevval(current)<0)
7 U. e* Q2 [, x   THEN first_opening :=true
' V7 F/ `  J6 o. n& Z- B' X        elapsed:=t-topen
- I0 J: V$ t+ I8 ]7 p" x  z+ e
  ELSIF state=2 THEN           -- was tentatively open   
5 o/ s9 m7 v: k9 i! T1 E
   IF elapsed>reignit_w  THEN
4 h+ M* x2 M( ~; F+ Z, G$ T      fully_opening :=true
   ELSIF  abs(deltav) > ds
' l$ W" I" Z1 J5 Z      THEN reigniting:=true
4 D+ ^# }4 }* L   ENDIF
0 t; b0 ~: m( M9 }" y; k
: Q( \, B  u8 Y9 }8 A1 n% s+ z  ELSIF state=3 THEN                  -- was reignited

   reopen_clock :=reopen_clock +timestep

   IF reopen_clock >=reopen_delay
    AND (abs(current) <ich1
6 g- Y8 I" u% |8 e" A2 m  k         OR current*prevval(current)<0)
    AND abs(deriv(current))<slope
) [) {" B8 q% @- g" j8 V        THEN reopening:=true
   ENDIF

  ELSIF state=0 THEN                 -- was fully open
  ENDIF
" \/ U* O$ x* \, N' k0 {
( y* e4 M+ z! P8 n" v( [" [  IF first_opening OR reopening THEN
   state :=2               -- becomes tentatively open
2 j4 ?4 Z4 U1 t7 w+ n/ {1 z
  ELSIF reigniting THEN
   reopen_clock :=0
3 F: m: F0 _( r! Q" }   state :=3                      -- becomes reignited
- U( p/ s+ C7 j6 J3 A- S6 u  ELSIF fully_opening THEN
) D" c, c+ F! R   state :=0                     -- becomes fully open
2 H& T7 J- A5 z! e  ENDIF

  IF state = 2 OR state = 0             -- is now open
+ x/ ]2 K* ^  v' t0 o0 n  THEN cbstatus :=open
# Y6 N# K! \# Q! G. N5 B  ELSE cbstatus :=closed
  ENDIF

' E5 E) @( \, W2 C3 W" _. C
ENDIF

ENDEXEC
/ w, u3 }; L  B; P& hENDMODEL
4 P) a: C( ~# T' j9 n& p

ziwaichengxiang

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

duym

@ZY2643181404