ansys瞬变电流激励下，计算中在某节点值过大和没有出趋肤效应的原因是什么？

woods12

　　正在做的一个瞬态分析，结果总在某步出现节点过大，超出计算范围程序中断的问题，而前面的计算结果显示没有出现趋肤效应。请大家帮忙看一下。
　　模型很简单，就是几根直导线置于一个导体上（为了简化，这个导体设为了非铁磁材料）。当导线中通有脉冲电压时，求导体内的涡流。（其实本意想得到导体内所受的洛仑兹力，但先求出涡流再说。）
　　文末是命令流。运行的中间会出现这样的错误： The value of VOLT at node 3387 is 7.389323872E+25.  It is greater than  the current limit of 1.E+25.  This generally indicates rigid body motion as a result of an unconstrained model.  Verify that your model is properly constrained. 但我想这个模型是由于受到洛仑兹力，而应该有小的变形，但还不至于运动量大到发生剧烈的运动。 所以应该有别的错误在里面。请教大家了。附件1中是线圈和导体的模型，2附件是某步中的磁力线图，结果太不正常了。
　　
finish  
/clear     
/CWD,'F:\YJT\ansys\example\322'
/FILNAME,example-eddycurrent323        
/prep7
pi=4*atan(1)
n=30
ra=0.03        !单位米
rb=0.001
l=0.2
h=0.001
radius_air=1.2*l
t=0.05
!f=1e5
!v=24
to=1e-6
tp=1e-2
tim=0
v0=200
f0=1/to
w=2*pi*f0
sw=0.8e-6                !单根导线面积
!*******************趋肤深度*******************************
cond=1e7/1.5
muzero=4e-7*pi
mur=1
skind=sqrt(1/(pi*f0*cond*muzero*mur))
/pnum,area,1
RECTNG,-ra/2,ra/2,h,h+rb,
RECTNG,-l/2,l/2,-t,0,
RECTNG,-l/2,l/2,-skind*2,0,
RECTNG,-l/2,l/2,0,2*h+rb,
CYL4,0,0,radius_air
aovlp,all
numcmp,area
/replot
et,1,53,2, ,  !线圈
et,2,53,1, ,  !被测导体
et,3,53, , ,  !空气层
emunit,mks
mp,murx,1,1       !线圈
mp,rsvx,1,1e-7/5.7
mp,murx,2,1.03    !被测导体
mp,rsvx,2,1.5e-7
mp,murx,3,1       !！空气层

asel,s,area, ,1    !！线圈
aatt,1,1,1
asum
*get,a,area, ,area
asel,s,area, ,2,3,1     !！被测导体
aatt,2,1,2

asel,s,area, ,4,5,1     !空气层
aatt,3,1,3
asel,all
lsel,s,, ,1,3,2
LESIZE,all, , ,n, , , , ,1  !！线圈
lsel,s,, ,23,24,1
LESIZE,all, , ,2, , , , ,1
lsel,s,, ,7
LESIZE,all, , ,150, , , , ,1
asel,s,,,1,3,1
mshape,0,2d
mshkey,1
amesh,all
smrt,1            ! 空气隙
asel,s,area, ,4
mshape,0,2d
mshkey,0
amesh,4
smrt,2
asel,s,area, ,5
mshape,0,2d      !！空气层
mshkey,0
amesh,all

r,1,a,n, 0.04,　,sw/a
n1=node(-ra/2,h,0)
esel,s,mat, ,1
nsle,s
cp,1,curr,all
allsel,all
finish
/solu
csys,1
nsel,s,loc,x,radius_air
d,all,az,0
allsel,all
to=1e-6
tim=0
v0=500
antype, trans
tstep=to/20
timerest=2*to
i=0
timint,on
!******************************脉冲信号循环加载*************************
*do,i,0,2*to,tstep
esel,,mat,,1
v=v0*(1-COS(w*i/n))*COS(w*i)　　　
bfe, all,vltg,1,v　　　　　　　!BFE,P51X,VLTG,1,v0*(1-COS(w*i/n))*COS(w*i), ,  
alls
tim=tim+tstep  
time,tim            !设定瞬态分析时间--载荷步的结束时间
autots,on                !打开自动时间步长
NSUBST,2  
!time,to/10
!deltim, to/10
outres,all,all
solve
*enddo

finish

[ 本帖最后由 woods12 于 2007-3-23 10:15 编辑 ]

perch88

我没有运行你的程序，但是看得出：

首先，空气隙似乎不够大
aovlp这个命令不对，导致overlap没有真正发生