在柱坐标系下加载table，实现移动热源功能

ifly7

求大家帮忙看一看，到底出了什么问题，不能实现移动热源功能，并且感觉table加载不上。
贴下命令流！
/prep7  
/UNITS,si
/view,1,1,1,1

et,1,solid70
!材料属性
mp,c,1,475
mp,kxx,1,45
mp,dens,1,7830

r1=0.02         !内圈半径
r2=0.023485  !外圈半径
lz=0.015        !圆柱厚度
n1=50           !周向1/4划分网格份数
n2=20           !轴向划分网格份数
n3=5             !径向划分网格份数
v=0.002         !热源速度
pi=3.1415926  !圆周率
lc=0.002         !热源长度
lsize1=2*pi*r2/4/n1   !外圈周向单元长度
lsize2=lz/n2    !外圈轴向单元长度
tinc=lsize1/v   !热源经过一个单元所用时间
tm=2*pi*r2/v !热源走过外圈一周所用时间
max_y=4*n1+1 !y坐标值个数
max_z=n2+1  !z坐标值个数
max_time=1+tm/tinc !time坐标值个数

cylind, r1, r2,0,lz ！建模
mshape,0      
mshkey,1
wprota,, 90
vsbw,all
wprota,,,90
vsbw,all
wpcsys,-1,0
type,1
lesize,1,,,n1
lesize,2,,,n1
lesize,3,,,n1
lesize,4,,,n1
lesize,5,,,n1
lesize,6,,,n1
lesize,7,,,n1
lesize,8,,,n1
lesize,17,,,n2
lesize,21,,,n3
vmesh,all       !网格划分
finish
CSWPLA,11,1,1,1,  !建立局部坐标系，与整体柱坐标系重合

*dim,heatflux,table,max_y,max_z,max_time,y,z,time!定义数组和它的尺寸
*do,k,1,max_time,1   !time坐标

*do,i,1,max_y,1        !y坐标，我认为是柱坐标系下y坐标，并且是角度，不是弧度
heatflux(i,0,k)=-180+(i-1)*lsize1/pi/r2*180  !坐标值范围从-180-180
*enddo

*do,j,1,max_z,1 !z坐标，也就是柱坐标系下z坐标
heatflux(0,j,k)=(j-1)*lsize2 !经过轴向的距离
*enddo

*do,i,1,max_y,1
*do,j,1,max_z,1
l_y=v*(k-1)*tinc   !热源所在位置
distance= (i-1)*lsize1- l_y
*if,distance,gt,0,then
heatflux(i,j,k)=0
*elseif,distance,ge,-1*lc,then
heatflux(i,j,k)=4500000* (1- abs(distance) /lc)
*else
heatflux(i,j,k)=0
*endif
*enddo
*enddo
heatflux(0,0,k)=(k-1)*tinc  !time坐标
*enddo
finish
、
/solu
antype,trans
tunif,20
*do,t,1e-6,tm,tinc
time,t
autots,0
deltim,,,,1
kbc,0
outres,all,all
tsres,erase
allsel,all
nsel,r,loc,x,r2
!NROTAT,ALL
sf,all,hflux,%heatflux%
allsel,all
solve
allsel,all
nsel,r,loc,x,r2
sfdele,all,hflux
*enddo
Finish

cg_hero

你好，请问你的问题解决了吗？

舒凤远

我也有这样的问题啊

湖之子

不知道解决没，我的圆柱面加载温度场，也遇到问题了，变成了一个均布的温度场了