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