我的热源为什么不能移动 请大家指点下

复合材料22

我的模型是一个体积热源与面热源同时施加在体上移动
两个热源共圆心  为什么计算后并不能显示热源移动
而且文章的结果是同激光那个移动差不多
而我的结果怎么相差那么大？

/clear
finish
/CONFIG,NRES,10000
/prep7

et,1,70    !定义材料常数
et,2,surf152
keyopt,2,4,1
keyopt,2,5,0
keyopt,2,8,1       !Include heat flux, ignore convection

mp,dens,1,2710

mptemp,,25,100,200,300,400

mpdata,kxx,1,,217,230,247.2,268.2,272.4

mpdata,c,1,,930,946.4,962.2,999.2,994.4
/VIEW,1,1,1,1
XA=0.3     !模型长度
YA=0.1     !模型宽度
ZA=0.003    !模型厚度
XNUM=40     !X方向上单元个数
YNUM=25     !y方向上单元个数
ZNUM=5      !z方向上单元个数

BLOCK,0,XA,0,YA,0,ZA
SAVE
LESIZE,4,,,XNUM
LESIZE,1,,,YNUM
LESIZE,9,,,ZNUM
MSHKEY,1
VATT,1,1,1
VMESH,1

EPLOT

TYPE,2  ! 第2个面用surf152划分
asel,s,,,2
AMESH,ALL

ALLSEL,ALL
LSIZE=XA/XNUM      !x方向单元长度
LSIZE2=YA/YNUM     !Y方向单元长度
V=0.002         !焊接速度
TINC=LSIZE/V    !经过每两个节点间的时间
TT=150

RA=0.004    !焊接头半径

RB=0.012  !轴肩半径

/SOLU

ANTYPE,4

TUNIF,25
KBC,1

TIMINT,ON

asel,s,,,1   !施加conv
nsla,s,0
sf,all,conv,200,23
allsel,all

asel,s,,,2
nsla,s,0
sf,all,conv,40,23
ALLSEL,ALL

*DO,TM,8,TT,TINC  ! 开始循环

/VIEW,1,1,1,1

TIME,TM
NSUBST,5,20,5
AUTOTS,ON

KBC,0
!*
TSRES,ERASE        
XCENTERA=V*TM        !建立移动的局部坐标系为柱坐标系
YCENTERA=0
NUM=NODE(XCENTERA,0,0)
NWPAVE,NUM
CSWPLA,11,1,1,1,

NSEL,S,LOC,X,0,RA           !施加体积热源
BF,ALL,HGEN,0.5e6

nsel,s,loc,z,za             !施加表面热源
nsel,r,loc,x,ra,rb
esln
sfe,all,,hflux,,3e5

ALLSEL,ALL
CSDELE,11, ,1,  
WPAVE,0,0,0
OUTRES,ALL,ALL,
SOLVE

*ENDDO

FINISH

sunclark

我认真给你看了一遍你的程序，问题还是有点多！！修改好的给你贴在下面：
     需要说明的是：1、/CONFIG命令不是前处理命令，所以不能写在APDL程序中
     2、既然你用的单位都是统一的，你可以定义单位为国际单位，/units,si
     3、sf,all,conv,200,23和sf,all,conv,40,23用错了，第四个应该是值，后面一个才是温度值，反了
     4、你的移动坐标系应该和热源一起加载移动，才不会出现你的热源不能移动的问题，可在坐标系定义一个变量，然后让变量增大，可实现热源移动
     5、热源加载后随之删除，你没有删除
     6、时间TIME不是固定的，而应该是随着你的坐标系移动而不断增长；再者，时间一般放在加载后，卸载前
     你的程序有点乱，可参照比较成功易董的例子来写，就象C语言一样，一定要易懂