ikuangke 发表于 2011-4-15 17:28:21

求助:ansys高斯热源加载

我做了一个三维高斯热源加载的命令流,定义了一个四维表格,但运行的时候就老出现如图所示的错误
应该是其中第三个变量超出了表格范围,请教高手,到底是哪个步骤出了问题,该如何解决?
命令流如下:
clear
/filname,seal
/title,temp seal
/units,si
/prep7
lx=0.1
ly=0.04
lz=0.01
lsize=0.0025
pi=3.1415926
r=0.002
v=0.005
p=1200
absor=0.8
tm=lx/v
tinc=lsize/v
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=ly/lsize
max_z=lz/lsize
qm=3*absor*p/(pi*r**2)
et,1,solid70
mp,dens,1,7840
mptemp,1,20,100,300,500,700,900,1100
mpdata,kxx,1,1,47,43,38,34,28,25,24
mpdata,c,1,1,472,480,524,615,854,637,595
/view,1,1,1,1
block,0,lx,0,ly,0,lz
esize,lsize
mshkey,1
vatt,1,1,1
vmesh,all
*dim,flux2,tab4,max_x,max_y,max_z,max_time,x,y,z,time
*do,k,1,max_time,1
*do,h,1,max_x,1
flux2(h,0,0,k)=(h-1)*lsize
*enddo
*do,i,1,max_y,1
flux2(0,i,0,l)=(i-1)*lsize
*enddo
*do,j,1,max_z,1
flux2(0,0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
*do,k,1,max_z,1
xcenter=v*(k-1)*tinc
ycenter=0
zcenter=0
distance=sqrt(((h-1)*lsize-xcenter)**2+((i-1)*lsize-ycenter)**2+((j-1)*lsize-zcenter)**2)
*if,distance,le,r,then
flux2(h,i,j,k)=qm*exp(-3*distance**2/(r**2))
*else
flux2(h,i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,0,k)=(k-1)*tinc
*enddo
/solu
antype,trans
tunif,20
kbc,1,timint,on
autots,off
esel,s,type,,2
nesl,s
sf,all,conv,12.5,20
allsel
*do,t,1e-6,tm,tinc
time,t
autots,off
deltim,,,,1
kbc,0
outres,all,all
tsres,erase
sf,all,hflux,%flux2%
solve
*enddo
finish

三毛 发表于 2011-11-15 10:00:02

哥们你的的问题解决了吗?我也遇到相同的问题。

qiminliu 发表于 2011-11-16 16:45:14

我拷回去看看

poetal 发表于 2012-5-15 20:23:12

也遇到这个问题,楼上兄弟们解决了吗
页: [1]
查看完整版本: 求助:ansys高斯热源加载