【急】大侠帮我看看 焊接移动 问题出在哪里
我想用来模拟焊接移动，并求它的温度场，但是算到最后总是出问题，结果不理想，想请大侠帮我看看，问问出在哪里，非常感谢，以下是我的命令流：
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 condition25degrees
!
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 我也要学这种东西，请问楼主，关于焊接的有实例的书籍有那些噢，找不到合适的书看。
页:
[1]