马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
!这是一个加表面效应单元的例子,你看一下有没有帮助。
FINISH
/CLEAR
!!!gauss加载
/title, the simulate of laser cladding process
*SET,P,2000, !the power of the laser
*SET,V,0.008 !scanning speed
!*SET, RADIUS,0.003 !radius of the laser beam
*SET, R,0.003
*set,k,0.6 !the coefficient of the absorb
*SET,L,0.05 !the length of the block
/prep7
!熔覆层高5mm=0.005m
!基体高10mm=0.01m,长x=50mm=0.05,宽=20mm=0.02m
BLOCK,0,0.05,0, -0.0105,-0.01,0,
Vsel,all
vplot
!
/USER, 1
/VIEW, 1, -0.198026943995E-01, 0.382677277552 , 0.923669829831
/ANG, 1, -0.107453019382
/REPLO
!
ET,1,SOLID70
ET,2,Surf152
DOF,TEMP
KEYOPT,2,3,0
KEYOPT,2,4,1
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,4 !根据固体表面温度计算对流系数
KEYOPT,2,9,0
!用温度表方式定义材料性能
Mp,kxx,1,82.9
Mp,dens,1,8900
Mp,c,1,471
mptemp,1,20,300,600,750,900,1200,1500,1800,2000
mpdata,hf,1,1,6,50,120,180,200,250,378,700,850
!
!
Mshkey,0 !自由网格化
Mshape,1,3d
type,1
esize,0.002,
Vmesh,all
eplot
!
!
asel,s,loc,y,0 !对流边界条件要施加在5个面的周围
asel,a,loc,z,-0.01
asel,a,loc,x,0
asel,a,loc,x,0.05
type,2 !定义属性
amesh,all
esel,s,type,,2 !选择表面效应单元
sf,all,conv,-1,25 !确定对流边界条件,表格加载
!sf,all,conv,200,25 !确定对流边界条件
!!底面也施加轻微的对流
allsel,all
asel,s,loc,y, -0.0105
type,2 !定义属性
amesh,all
esel,s,type,,2
sf,all,conv,3,25
allsel,all
!
nsel,all
/psf,hflux,,2,0 !以箭头方式显示载荷
!
finish
/SOLU
!!!!!
antype,transient,new
TRNOPT,FULL
Nropt,,,off !关闭自适应下降
LUMPM,0
Solcontrol,on
Autots,on
Kbc,1
Outres,all,last,
Lnsrch,on
neqit,50 !每一子步中方程的迭代次数限值
!load step1,initial conditions 25
time,0.001
deltim,0.001,0.001,0.001
allsel,all
timint,off
tunif,25 !均匀的节点初始温度25,只对分析的第一个子步有效
solve
EPLOT
!
! load step2………,apply moving heat flux
allsel
nsel,all
t=0
dt=0.5
local,11,1, 0.001+0.005,0,0,0,-90,0
CSYS,11
esel,s,type,,1
nsel,r,loc,x, 0,0.003
nsel,r,loc,y,0,180
nsel,r,loc,z,0
sf,,hflux,7000000
!sf,,hflux,%flux2%
csys,0
t=t+dt
time,t !每一载荷步结束时的时间
deltim,0.1,0.05,0.1,on
kbc,1 !设定为阶越载荷
!neqit,25 !每个子步的默认次数
autots,on !打开自动时间步长
timint,on,therm !时间积分效果
TINTP,,,,1
Lnsrch,on
!outres,all,all
outres,all,last
!outres,all,10
allsel
solve
save |