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

duym

本来想把电抗器重燃模型的model模块model用txt文件贴出来，好像不能上传txt文件，故在下面贴出了model的内容。
' E- e$ ~  }) y
& p  \* A' P: }MODEL breaker

COMMENT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
8 o- B" `1 Q: l7 u* \* i! R
! `7 ~( n. R$ A7 Z2 w9 q
The model of this example  updates  the  present  state  of  the
circuit  breaker  by  taking   into   account   the   conditions
5 K: V) ?) \, I" b  Ycontrolling four represented state transitions:  first  opening,
reigniting, reopening, and full opening.
% u! U, x! w4 K
* I6 f9 F* V" w* o& X# @>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ENDCOMMENT

, }1 g8 c: H# {8 ?DATA
  topen                -- [s]   earliest elec opening
) U0 N6 m! B' _% D8 j+ J  ich1                 -- [A]   chopping level at power frequency
  slope                -- [A/s] current slope
  reopen_delay {dflt:2.5e-6}         -- [s]  reopening delay
6 R* l2 V" ^* V4 z. i1 Q  reignit_w {dflt:1e-2}              -- [s]  reignition window

. ^& U5 |9 i. h: W4 AINPUT
9 U& V# e) m1 u8 g2 U  v1   -- voltage across breaker terminals
  v2
  current  -- current through breaker

VAR
  state    -- 0= fully open
7 J: F5 a! i' N# [! g0 J           -- 1= normal closed
           -- 2= tentatively open
           -- 3= reignited
  first_opening  -- transition from 1 to 2
" J! \* h+ b3 k% O  N6 b) L  reigniting     -- transition from 2 to 3
  reopening      -- transition from 3 to 2
$ r0 ?, @+ Q; a3 W  fully_opening  -- transition from 2 to 0
& }* A4 t, m' C' C! O2 ], |# d  reopen_clock   -- state variable of 3
  elapsed        -- time elapsed since tbeg
  cbstatus       -- open/close signal to switch
  deltav  
  ds
  nds
1 S$ K$ c* w8 \2 }0 W  tdelta
  itopen
+ R) N5 e: {( e+ w' V& x# I  ich2
) p6 R' N' y/ K! G
( a* D+ [6 J* M6 X! w$ K" x  tOUTPUT
0 u0 w& n$ [- i  cbstatus ,ds,nds

% |, T7 ?6 E+ F2 HHISTORY
  current {dflt: 0}

' c6 _7 F' E5 r- s" ZINIT
  state    :=1       -- normal closed
  cbstatus :=closed
! A. U$ _; P5 e- s/ X  ds:=1
! u: A2 B+ N3 V+ i  nds:=-1
% c2 [7 E9 T" g; n4 t0 DENDINIT

! ]6 X1 y; W9 XEXEC
4 W) c+ B7 V. y2 M) q" M+ P deltav:=v1-v2
' q" e! X1 X) s% v IF t>=topen THEN              -- otherwise do nothing
  tdelta:=(t-topen)*1000
  IF tdelta<=0.44
    THEN ds:=(38.5*tdelta**1.5)*1000
+ ~0 Z  T3 Q6 J  T! i$ x  ELSIF tdelta>0.44 AND tdelta<=5.00
    THEN ds:=55.1*tdelta**0.55*1000
  ENDIF
7 ~- T( E- O8 N/ A
  nds:=-ds
, f* h3 \! J7 A) J  itopen:=abs(current)
' V# l1 J9 S) G/ K2 k  ich2:=(2*pi*50*14.3*itopen*6.3*exp(-16))**((1-14.3)**(-1))

  first_opening :=false;  fully_opening :=false
1 _, @- C5 a' v0 G1 b3 h- Z  reopening     :=false;  reigniting    :=false
. T* Q% e$ i( y
  IF state=1                      -- was normal closed
   AND (abs(current) <ich1
1 [/ c0 n3 E4 O9 J        OR current*prevval(current)<0)
' U9 M0 Z- g/ x' Z" s% ?% Y) ^   THEN first_opening :=true
. h# b/ [# P4 s        elapsed:=t-topen
3 B7 T0 t9 ~* [( c* K: I# J
  ELSIF state=2 THEN           -- was tentatively open   

   IF elapsed>reignit_w  THEN
0 x3 j- V$ @; n0 _8 O      fully_opening :=true
' |$ |1 R8 A0 R' ^1 U   ELSIF  abs(deltav) > ds
8 f" q! W2 t$ {9 U3 W" q! u      THEN reigniting:=true
   ENDIF

+ W4 B7 J! \4 A/ A' n8 @0 n7 u  ELSIF state=3 THEN                  -- was reignited

   reopen_clock :=reopen_clock +timestep
  ?- {  T, |* e( f7 h
   IF reopen_clock >=reopen_delay
- p. f; r& Z8 W: Z+ t6 ]    AND (abs(current) <ich1
         OR current*prevval(current)<0)
$ D' _/ j& y& e* p5 C    AND abs(deriv(current))<slope
6 L2 M9 `3 z( g! C5 x" y! Q        THEN reopening:=true
3 w. b! d: x7 O( }" M' {   ENDIF
; u5 P/ e2 m& M9 |4 X+ U0 o
- F3 i+ N. {. f# ]) q  ELSIF state=0 THEN                 -- was fully open
, @5 k- T6 n8 W- Q" n' K3 ^  ENDIF

" z0 h8 Y4 g+ ^1 O7 e  V  IF first_opening OR reopening THEN
  |* ^1 i. e- u7 d   state :=2               -- becomes tentatively open
: {- U$ l! K7 F) V1 `7 j9 [7 C% R
  ELSIF reigniting THEN
9 ^/ ~7 n. M% j% M2 z   reopen_clock :=0
! O) i8 r* a  C$ B% t/ J, v6 ~   state :=3                      -- becomes reignited
9 K- _, {# I" f6 S# I- @! m* _  ELSIF fully_opening THEN
+ o' Y: J. E5 S" P7 L+ [( w+ t   state :=0                     -- becomes fully open
  ENDIF
2 k3 l; p) x7 x, W) q  `
  IF state = 2 OR state = 0             -- is now open
  THEN cbstatus :=open
# Z$ I$ u5 R( l! Y/ R  ELSE cbstatus :=closed
" [3 p  K9 S* J7 Z* C( q! Q" L  ENDIF
) \1 c3 Y6 n4 p  _6 A/ h

ENDIF

8 B* A' N  ^* rENDEXEC
0 Z( Q! ~) F3 [1 GENDMODEL


1 S) S" x- n1 W+ q' x* `; I

ziwaichengxiang

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

duym

@ZY2643181404  
9 g; I2 ~2 Y2 n. N( J* U