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