模型大致是,在一个基板上有一个0.3mm熔覆层,想在熔覆层上表面加载移动的高斯热源模型。在设置时,我用到了设定了一个局部坐标系12,不知道是否有影响
命令流如下:
Finish
/filename,TA15 Ti2AlNb SINGLE
/prep7
et,1,solid70
et,2,surf152
keyopt,2,4,1 !solid 70单元没有中间节点
keyopt,2,5,0 !表面效应单元加载对流,0和1均可?
keyopt,2,8,4 !根据流体温度计算对流系数
!根据Jmat Pro 模拟数据输入材料参数
!2号材料为薄壁
MPREAD,'TA15','dat','C:\Users\laser\Desktop\'
MPREAD,'Ti-22Al-25Nb','dat','C:\Users\laser\Desktop\'
block,0,0.15,0,0.006,0,0.1
wpave,0.05,0.006,0.05
block,0,0.048,0,0.0003,0,0.005
WPAVE,0,0,0
CSYS,0
vglue,all
wpave,0.05,0.006,0.05
vsbw,all,,delete
wpave,0.05,0.006,0.055
vsbw,all,,delete
wprot,0,0,90
vsbw,all,,delete
wpave,0.098,0.006,0.055
vsbw,all,,delete
nummrg,all
numcmp,all
vsel,s,,,2
vatt,2,,1,0
vsel,s,,,1
vsel,a,,,3,10,1
vatt,1,,1,0
allsel
vsel,s,mat,,2
esize,0.001
MSHAPE,0,3d
MSHKEY,1
vmesh,all
vsel,s,mat,,1
esize,0.003
MSHAPE,0,3d
MSHKEY,1
vmesh,all !建模 划分网格
/solu
antype,trans !瞬态分析
OUTRES,,ALL !输出设置
toffst,273 !温度偏移量设置
nropt,full
timint,1,therm !打开瞬态热分析
tunif,25 !初始温度
autos,on !自动时间步长
solcontrol,on
kbc,0 !渐变载荷
pred,on !打开时间预测器
lnsrch,on !瞬态热分析设置
CNVTOL,HEAT, ,0.01,1,0.000001, !控制收敛准则
LUMPM,0 !不使用集中质量矩阵
asel,s,loc,x,0
sfa,all,,conv,30,25
asel,s,loc,x,0.15
sfa,all,,conv,30,25
asel,s,loc,z,0
sfa,all,,conv,30,25
asel,s,loc,z,0.1
sfa,all,,conv,30,25
asel,s,loc,y,0
sfa,all,,conv,30,25
asel,s,loc,y,0.006
asel,u,,,9
sfa,all,,conv,30,25
local,12,0,0.05,0.0063,0.05
csys,12
asel,s,loc,y,0
lsize=0.001
v=0.003
tinc=lsize/v
r=0.002
qmax=1.31e8
max_time=49
max_x=49
max_z=6
*dim,flux2,table,max_x,max_z,max_time,x,z,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*lsize
*enddo
*do,j,1,max_z,1
flux2(0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_z,1
xcenter=v*(k-1)*tinc
zcenter=0.0025
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-zcenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/exp(3*(distance**2)/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
*do,t,tinc,16,tinc
asel,s,loc,y,0
TSRES,ERASE
sfa,all,1,hflux,%flux2%
allsel,all
deltim,0.1,0.05,1/3
time,t
solve
*enddo
求指点,哪出错了
![求助,利用table加载高斯热源模型求解没用]()
模型.jpg |