非常着急!我用的do循环语句有问题吗?热源就是不动,请高手们指点!谢谢了。下面是命令流,请帮忙看看哪里有问题。
/prep7
/units,si
length=0.06
v=0.01
step=0.002
count=length/step
ET,1,SOLID70
!材料参数
!建立模型 划分网格
!位移约束
!-->!定义高斯函数
*DEL,_FNCNAME
*DEL,_FNCMTID
*SET,_FNCNAME,'gaosi3'
! /INPUT,gaosi3.func
*DIM,%_FNCNAME%,TABLE,6,26,1
!
! Begin of equation: 2.1486*10^7*exp(-8.33*10^4*({X}^2+{Y}^2+{Z}^2))
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 10, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 7, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 2.1486, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 8.33, 0, 0, -4
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 10, 0, 0, 0
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 4, 0, 0, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -1, 17, -4
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 1, -3, 3, -5
*SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, 2, 17, -3
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,18,1), 0.0, -5, 0, 1, 3, 17, -3
*SET,%_FNCNAME%(0,19,1), 0.0, -3, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,20,1), 0.0, -4, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,21,1), 0.0, -5, 0, 1, 4, 17, -4
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, -3, 1, -5
*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 7, 1, -3, 0, 0
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,26,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 2.1486*10^7*exp(-8.33*10^4*({X}^2+{Y}^2+{Z}^2))
!-->
fini
/solu
!
antype,trans!设置分析选项
trnopt,full
nropt,auto, ,
eqslv,
!
solcontrol,on
autots,on
!
kbc,1 !step
outres,all,5,
!
!初始条件!
time,0.001
deltim,0.001,0.001,0.001
timint,off! time integration=off
tunif,25,
sfa,1,,conv,18,25!空气自然流动
sfa,2,,conv,18,25!空气自然流动
sfa,7,,conv,18,25!空气自然流动
sfa,8,,conv,18,25!空气自然流动
sfa,16,,conv,18,25!空气自然流动
sfa,17,,conv,18,25!空气自然流动
sfa,4,,conv,18,25!空气自然流动
sfa,14,,conv,18,25!空气自然流动
solve
!
WPSTYLE,,,,,,,,0
*do,I,0,count !循环加载热源
disx=I*step
T=(disx+step)/v
Time,T
autots,1
antype,4
kbc,1
outres,basic,last
trnopt,full
nropt,full,,on
nsubst,5,,,1
lnsrch,1
tsres,erase
timint,1,therm
local,11,1,disx,0,0!定义局部坐标系,选节点
nsel,s,loc,y,-0.003,0.003
nsel,r,loc,x,0,0.002
nsel,r,loc,z,0
nplot
local,12,0,disx,0,0!定义局部坐标系,加载荷
sf,all,hflux,%gaosi3%
allsel,all
solve
sfdele,all,hflux
*enddo
time,30 !冷却时间
nsubst,5,,,1
solve
outres,all,5 !写入所有数据
finish |