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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
. y; z, g4 d) r# g+ r+ a/ y( g
; q+ d: i- P4 W; u- B, ?& O$ bMODEL breaker

  N4 Z; |& R, J) g$ h' JCOMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
) D! W) Y' M% Z( J

The model of this example  updates  the  present  state  of  the
7 ~: ^" U2 a4 h  Z2 @circuit  breaker  by  taking   into   account   the   conditions
controlling four represented state transitions:  first  opening,
# w. x2 E, y& C; X6 ]# B6 Yreigniting, reopening, and full opening.
  {$ N9 ]* g5 w$ g; e
8 y6 Z0 ]2 q5 }6 O) [>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT

# R6 E: \' o" T, N# dDATA
3 u6 g5 l$ ]& ~0 k$ F  topen                -- [s]   earliest elec opening
- `7 w2 ]7 Z2 U1 _  ich1                 -- [A]   chopping level at power frequency
  slope                -- [A/s] current slope
. {( E0 j0 H& E  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
4 E. Y, t- q+ @6 D/ s/ d3 Y; ^  reignit_w {dflt:1e-2}              -- [s]  reignition window

INPUT
  v1   -- voltage across breaker terminals
  v2
, ?6 O2 [1 J$ _  current  -- current through breaker

VAR
  state    -- 0= fully open
           -- 1= normal closed
           -- 2= tentatively open
- Q* I- L+ H1 {+ D7 n           -- 3= reignited
  first_opening  -- transition from 1 to 2
8 J) [/ T$ W5 D, v7 e$ O) I' O  reigniting     -- transition from 2 to 3
  reopening      -- transition from 3 to 2
1 N% b& o/ W# J6 r8 L  fully_opening  -- transition from 2 to 0
  reopen_clock   -- state variable of 3
  elapsed        -- time elapsed since tbeg
. \7 `: u% C1 T% N9 z  ~  cbstatus       -- open/close signal to switch
) @& W: B% G+ E  deltav  
  ds
  nds
  tdelta
  itopen
' \- O! c8 Y1 H' I5 F( X  ich2
+ ]! p% `1 A( u. z! Z( K
9 L: `8 m# Z  k4 v7 G+ qOUTPUT
  cbstatus ,ds,nds

HISTORY
  current {dflt: 0}

5 B, k  F# [; q9 R: h& @8 ~$ j- N# zINIT
  state    :=1       -- normal closed
  cbstatus :=closed
  ds:=1
- M: @; }  _! ]# i/ @  nds:=-1
ENDINIT
) ?2 i) S! c6 b9 @9 u" e
9 t& X9 w0 ]7 i+ L4 k9 ]% jEXEC
deltav:=v1-v2
. X* ~, d- a5 | IF t>=topen THEN              -- otherwise do nothing
5 {) W' ?% y9 s  Y  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
; M$ A* |* {1 z( V5 b    THEN ds:=(38.5*tdelta**1.5)*1000
  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
* A  J5 v/ K; [. ^  g  ENDIF

  nds:=-ds
/ |# A, P/ N  D4 E5 U4 G* p  itopen:=abs(current)
' r4 O3 }. H) a  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))

3 X* x! [9 {/ i# y' \  first_opening :=false;  fully_opening :=false
  reopening     :=false;  reigniting    :=false

  IF state=1                      -- was normal closed
4 j7 z6 m' y4 H+ ^   AND (abs(current) <ich1
        OR current*prevval(current)<0)
$ Y; X" j6 h$ {8 ?3 _# O   THEN first_opening :=true
        elapsed:=t-topen

& p4 v: z* a& U( ^; A! _  ELSIF state=2 THEN           -- was tentatively open   

$ V8 R1 V! [: x) E- ]5 a2 G) Y   IF elapsed>reignit_w  THEN
      fully_opening :=true
1 T5 x( Z' N6 d4 E  g  E   ELSIF  abs(deltav) > ds
      THEN reigniting:=true
   ENDIF

( N8 O) g+ u" B  ]7 l1 [  ELSIF state=3 THEN                  -- was reignited
' j' W6 r+ C  |
+ ?& @4 M0 W1 n; E   reopen_clock :=reopen_clock +timestep

   IF reopen_clock >=reopen_delay
" v0 v; y- `% ?) r    AND (abs(current) <ich1
         OR current*prevval(current)<0)
' t& ?1 g% }( R" S, B7 A! Z$ m7 I8 b/ W    AND abs(deriv(current))<slope
        THEN reopening:=true
   ENDIF

  ELSIF state=0 THEN                 -- was fully open
* U5 `" W7 f+ l0 b  ENDIF
, Z/ X& W0 N5 [
  IF first_opening OR reopening THEN
' I) v$ l; ^8 w; p6 D   state :=2               -- becomes tentatively open
5 ~1 e: o! _3 m4 }- h
  ELSIF reigniting THEN
   reopen_clock :=0
   state :=3                      -- becomes reignited
- F& S  `# L) |9 [  ELSIF fully_opening THEN
$ B& s- o( u. `- d7 ]* h   state :=0                     -- becomes fully open
9 m* [* j' S) R" E7 c5 R6 H- L  ENDIF
. X2 k3 q/ O& Y. d
  IF state = 2 OR state = 0             -- is now open
  THEN cbstatus :=open
  A/ ~, y, M- O3 F% a2 z0 R  ELSE cbstatus :=closed
  ENDIF


ENDIF
; b9 G1 m  N2 U* ]
ENDEXEC
ENDMODEL


, J# ], X8 f' v4 G

ziwaichengxiang

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

duym

@ZY2643181404