24小时热门版块排行榜    

查看: 751  |  回复: 0

grinmtyj

新虫 (小有名气)

[交流] 求助,利用table加载高斯热源模型求解没用

模型大致是,在一个基板上有一个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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 grinmtyj 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见