关于计算线路感应电压模型讨论

不学无术

用MODELS计算线路感应电压：程序如下：
MODEL indloss
CONST
Tmax {VAL:500}
n {VAL:2}
c {VAL:3.e8}
I0 {VAL:1}
eps0 {val:8.85e-12}
INPUT UAP[0...n],UBP[0...n]
DATA Y,XA,XB
z[1...n] {DFLT:10}
Im       {DFLT:3e4}
tc       {DFLT:5e-5}
th       {DFLT:2e-6}
v        {DFLT:1.5e8}
sigma    {DFLT:0.001}
epsr     {DFLT:10}
OUTPUT UrA[0...n],UrB[0...n]
VAR UindA0[0...Tmax],UindB0[0...Tmax],g1[0...Tmax],g2[0...Tmax],
UindAD[0...Tmax],UindBD[0...Tmax],Ui0,UiD,
UrA[0...n],UrB[0...n],Tr,k,i,AB,dt,uk[1..4],
b,L,x,tau,a1,b1,c1,c2,c3
FUNCTION F0(x,tr):= (c*tr-x)/(y*y+(b*(c*tr-x))**2)
FUNCTION F1(x,tr):= (x+b*b*(c*tr-x))/sqrt((v*tr)**2+(1-b*b)*(x*x+y*y))
FUNCTION F2(x,tr):= x+b*b*(c*tr-x)+sqrt((v*tr)**2+(1-b*b)*(x*x+y*y))
FUNCTION F3(x,tr):= (v*tr+sqrt((v*tr)**2+(1-b*b)*(x*x+y*y)))/sqrt(x*x+y*y)
FUNCTION t0(x):= sqrt(x*x+y*y)/c
HISTORY
UrA[1...n] {dflt:0}
UrB[1...n] {dflt:0}
UAP[1...n] {dflt:0}
UBP[1...n] {dflt:0}
DELAY CELLS DFLT: (XA-XB)/(c*timestep)+1
INIT
k:=sqrt(eps0*epsr/(PI*sigma*dt))
dt:= timestep b:=v/c L:=XA-XB tau:=L/c
a1:=Im/(I0*th)
b1:=0.5*th/(tc-th)+1
c1:=a1*k*dt/sqrt(epsr)
c2:=-1.073*k+0.2153*(k**3)+4/3
c3:=-0.2153*(k**3)+1/6
FOR AB:=1 TO 2 DO
if AB=1 then
x:=XA else x:=-XB
endif
FOR i:=0 TO Tmax DO
Tr:=i*dt
if Tr>t0(x)
then
if Tr>t0(x-L)+tau
then
Ui0:=f0(x,Tr)*(f1(x,Tr)-f1(x-L,Tr-tau))
UiD:=ln( f2(x,Tr)/f2(x-L,Tr-tau) )-1/b*ln( f3(x,Tr)/f3(x-L,Tr-tau) )
else
Ui0:=f0(x,Tr)*(f1(x,Tr)+1)
UiD:=ln( f2(x,Tr)*f0(x,Tr) )-1/b*ln( f3(x,Tr)/(1+b) )
endif
else
Ui0:=0, UiD:=0
endif
if AB=1 then
UindA0[i]:=Ui0
UindAD[i]:=UiD
g1[i]:=a1*dt
if i>0
then
g2[i]:=a1*sqrt(eps0/(PI*sigma*Tr))*dt
endif
else
UindB0[i]:=Ui0
UindBD[i]:=UiD
endif
ENDFOR
ENDFOR
i:=Tmax
WHILE i>1 DO
uk[1...4]:=0
FOR Tr:=1 TO i-1 DO
uk[1]:=uk[1]+UindA0[Tr]*g1[i-Tr]
uk[2]:=uk[2]+UindB0[Tr]*g1[i-Tr]

uk[3]:=uk[3]+UindAD[Tr]*g2[i-Tr]

uk[4]:=uk[4]+UindBD[Tr]*g2[i-Tr]

ENDFOR
UindA0[i]:=uk[1]+0.5*UindA0[i]*g1[0]
UindB0[i]:=uk[2]+0.5*UindB0[i]*g1[0]
UindAD[i]:=uk[3]+(UindAD[i]*c2+UindAD[i-1]*c3)*c1
UindBD[i]:=uk[4]+(UindBD[i]*c2+UindBD[i-1]*c3)*c1
i:=i-1
ENDWHILE
Tr:=trunc(th/dt)

FOR i:=Tmax TO Tr BY -1 DO
UindA0[i]:=UindA0[i]-b1*UindA0[i-Tr]
UindB0[i]:=UindB0[i]-b1*UindB0[i-Tr]
UindAD[i]:=UindAD[i]-b1*UindAD[i-Tr]
UindBD[i]:=UindBD[i]-b1*UindBD[i-Tr]
endif
ENDFOR
ENDINIT
EXEC
FOR i:=1 to n DO
UrA[i]:=60*I0*b*(z[i]*UindA0[t/dt] - UindAD[t/dt]) +2*delay(UBP[i],tau-dt,1)-delay(UrB[i],tau,1)
UrB[i]:=60*I0*b*(z[i]*UindB0[t/dt] - UindBD[t/dt]) +2*delay(UAP[i],tau-dt,1)-delay(UrA[i],tau,1)
ENDFOR
ENDEXEC
ENDMODEL

程序报错如下
KILL = 422.  The present statement is interpreted as a value assignment statement, and the preceding name is interpreted as a
variable to which a value is being assigned.  The identification of this variable is interpreted as complete, and should now be
followed by ":=", in the following format:
     variable_identification := expression

本人已经仔细检查了每个等式都加：=了，不知道为什么还是不行，有兴趣的可以仿真一下，大家交流交流。

twjcly

Exa 15就是原模型 自己对照看下

不学无术

我看到你发的帖子了，但是已修改里面的参数，整个模型就不一样了，如果想修改sigma参数，你知道怎么改吗，改完之后模型不变

twjcly

对，要重新修改MODEL的界面呢，在Edit definitions里面，修改Data和NOdes界面，增加变量，更改变量的节点位置。另外你点回复帖子，这样你回复了我都看不到。

twjcly

不学无术 发表于 2017-10-13 14:02
我看到你发的帖子了，但是已修改里面的参数，整个模型就不一样了，如果想修改sigma参数，你知道怎么改吗， ...

每层楼下面有个回复，不是大的回复。