ansys如何模拟多脉冲激光烧蚀？ 50
我想请教一下诸位大佬，ansys里在已经对材料进行了一次单脉冲激光烧蚀后，如何让材料自由散热一段时间后，进行第二次单脉冲的激光烧蚀？以此达到模拟多脉冲激光的效果。如果考虑到热应力的作用，是否可以让材料每受到一次激光作用，就分析一次热应力情况，在具有残余应力的情况下受到第二次激光作用？
附上命令流，望各位大佬不吝赐教。
fini
/clear,all
/filname,therma
/prep7
qhwd=550 !K
nengl=200e-6
banj=50e-6
shij=1e-9
SHASI=qhwd-279
jisuancs=10
et,1,55
mp,ex,1,25e9
mp,prxy,1,0.1
mp,alpx,1,7e-5
mp,c,1,1982
mp,kxx,1,0.8
mp,dens,1,1800
rect,0,0.1e-3,0,-0.1e-3
rect,0,banj*1.1,0,-banj/4
aovlap,all
esize,0.2e-6
asel,,loc,x,0,banj*1.1
asel,r,loc,y,0,-banj/5
lsel,,loc,y,-1e-9,-banj/5+1e-9
lesize,all,,,90
amesh,all
alls
esize,2e-6
amesh,all
!!!!!!!!!!
!下面是定义function，名称为reliu
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS
*SET,_FNCNAME,'reliu'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),nengl
*SET,_FNC_C2(1),shij
*SET,_FNC_C3(1),banj
*SET,_FNCCSYS,0
! /INPUT,aa.func,,,1
*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%
!
! Begin of equation: 2*nengl/shij/3.1415/banj/banj*exp(-2*{X}*{X}/banj/banj)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 17
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17
*SET,%_FNCNAME%(0,3,1), 0, -1, 0, 1, -2, 4, 18
*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 3.1415, 0, 0, -1
*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 4, -2
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -3, 4, 19
*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, -1, 4, 19
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 3, 2
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 1, -1, 3, 2
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -3, 4, 19
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, -1, 4, 19
*SET,%_FNCNAME%(0,17,1), 0.0, -1, 7, 1, -3, 0, 0
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 2*nengl/shij/3.1415/banj/banj*exp(-2*{X}*{X}/banj/banj)
!-->
alls
/solu
antype,trans
NROPT,FULL
ic,all,temp,20
*do,I,1,jisuancs
esel,,live
nsle,s
nsel,r,ext
nsel,r,loc,x,1e-9,banj
nsel,u,loc,y,-0.1e-3+1e-9,-0.1e-3-1e-9
cm,zaihen,node
esel,,live
nsle,s
nsel,r,ext
nsel,r,loc,x,1e-9,-1e-9
cm,nodeee,node
n11=node(0,0,0)
cmsel,,zaihen,node
nsel,a,node,,n11
cm,jiazai,node
sf,all,hflux,%reliu%
alls
time,shij/jisuancs*i
alls
solve
esel,,live
nsle,s
cmsel,r,jiazai,node
sfdele,all,all
esel,,live
cm,lins,elem
*get,nmax,elem,,count
*do,kk,1,nmax
cmsel,,lins,elem
*get,nmin,elem,,num,min
n1=nelem(nmin,1)
n2=nelem(nmin,2)
n3=nelem(nmin,3)
n4=nelem(nmin,4)
alls
temp1=temp(n1)
temp2=temp(n2)
temp3=temp(n3)
temp4=temp(n4)
nmm=(temp1+temp2+temp3+temp4)/4
*if,nmm,gt,SHASI,then
ekill,nmin
*endif
cmsel,,lins,elem
esel,u,elem,,nmin
cm,lins,elem
*enddo
! *endif
alls
esel,u,live
cm,ekill%i%,elem
*enddo
FINI
/PREP7
ETCHG,TTS
/SOLU
nropt,full
alls
ealive,all !激活所有单元
NSEL,,LOC,Y,-0.1e-3+1e-9,-0.1e-3-1e-9
D,ALL,ALL
*do,I,1,jisuancs
time,shij/jisuancs*i
*if,i,gt,1,then
cmsel,,ekill%i%,elem
ekill,all
*endif
alls
LDREAD,TEMP,,,shij/jisuancs*i,,therma,rth
solve
*enddo
大龙猫🐱 