【急】大侠帮我看看 焊接移动 问题出在哪里

yjh3308

我想用来模拟焊接移动，并求它的温度场，但是算到最后总是出问题，结果不理想，想请大侠帮我看看，问问出在哪里，非常感谢，
以下是我的命令流：
finish
/clear
/title,3D moving heat model
/unit,si
/prep7
et,1,plane55
et,2,solid90
mptemp,1,20,400,800,1200,1450,1500
mpdata,kxx,1,1,55,46,25,30,100,100             !W/m/C
mpdata,c,1,1,400,550,850,795,120,100           !J/kg/c
mpdata,dens,1,1,7860,7860,7860,7860,7860,7800  !kg/m**3  
!建立模型
k,1,0,0,0
k,2,0.05,0,0
k,3,0.1,0,0
k,4,0,0.01,0
k,5,0.05,0.01,0
k,6,0.06,0.01,0
k,7,0.1,0.01,0
k,8,0.055,0.012,0
a,1,2,5,4
a,2,3,7,6
larc,5,6,8
al,2,8,9
!划分网格
type,1
mat,1
esize,0.001
amesh,1
amesh,2
mshkey,0
amesh,3
eplot
finish
!开始3D生成
/PREP7  
ESIZE, ,50,
TYPE,   2   
EXTOPT,ESIZE,0,0,   !定义扫掠方向的单元尺寸
EXTOPT,ACLEAR,1   !指定在体扫掠完成后是否删除源面上的网格
!*  
EXTOPT,ATTR,0,0,0   !定义扫掠体单元属性
MAT,1   
REAL,_Z4
ESYS,0  
!*  
!*  
VOFFST,1,0.1, ,
!*  
VOFFST,3,-0.1, ,
!*  
VOFFST,2,0.1, ,
allsel,all
/view,1,1,1,1

finish

!开始加载
lz=0.1
v=1/300
t1=lz/v
tinc1=0.002/v
z1=-0.002
/solu
antype,trans
trnopt,full
nropt,auto, ,
!
solcontrol,on
autots,on
outres,all,all
!load step 1 ,initial condition  25  degrees
!
time,0.001
deltim,0.001,0.001,0.001
timint,off! time integration=off
tunif,20
! 施加对流系数
sfa,1,,conv,110,20
sfa,2,,conv,110,20
sfa,4,,conv,110,20
sfa,7,,conv,110,20
sfa,8,,conv,110,20
sfa,13,,conv,110,20
sfa,15,,conv,110,20
sfa,16,,conv,110,20
sf,all,conv,110,20
solve
!单元生死技术load step 2,  apply moving heat source
!杀死所有单元
  vsel,s,volu,,2,,,1
  eslv,s
  ekill,all
  esel,s,live
  eplot
  timint,on        ! time integration=on
  timint,0,struct
  timint,1,therm
  tintp,0.005,,,1,0.5,0.2

*do,t,tinc1,t1,tinc1
  time,t
  autots,0
  deltim,tinc1,0.2,2
  kbc,1   
  outres,all,all
  z1=z1+0.002
  tsres,erase
  vsel,s,volu,,2,,,1      !定位
  nslv,r,1
  nsel,r,loc,z,z1
  esln,s,0,all
  ealive,all
  d,all,temp,1500        !加载
  esel,s,live
  eplot
  allsel,all
  solve
*enddo
time,40
nsubst,20
solve
outers,all,all
finish

wangyan43

我也要学这种东西，请问楼主，关于焊接的有实例的书籍有那些噢，找不到合适的书看。