激光焊的命令流（请教高手指点）

srwszsc

finish
/clear
/UNITS,SI
/filname,work1,1           
/title,kdjkfkdsjg      
  
/CONFIG,NRES,10000
/PREP7
*set,length,0.04   
*set,width,0.04
*set,height,0.005  
*SET,Vol,0.02                                 
*SET,powe,80                                 
*SET,lSIZE,0.001                             
*SET,TINC,lSIZE/Vol                           
*SET,Rad,0.001                                
*SET,pi,3.1415926
*SET,MAX_TIME,LENGTH/lSIZE+1                  
*SET,MAX_X,length/lSIZE+1                     
*SET,MAX_Y,width/lSIZE+1                     
*SET,power_max,0.85*powe/(pi*Rad*Rad)         
*set,kk,1                                    
et,1,solid70
keyopt,1,2,1
et,2,surf152                              
r,2
keyopt,2,4,1
keyopt,2,5,0
keyopt,2,8,4
MPTEMP,,,,,,,,  
MPTEMP,1,20,100,200,300,400,500
mptemp,7,600,700,800,1000,1200   
MPDATA,KXX,1,1,1.38,1.53,1.65,1.74,1.86,1.91     
mpdata,kxx,1,7,2.02,2.09,2.15,2.29,2.43   
mpdata,c,1,1,600.0,660.4,836.0,900.0,973.9,1000.0   
mpdata,c,1,7,1061.7,1108.0,1086.8,1103.5,1128.6   
mpdata,ex,1,1,7.78e10,7.78e10,7.78e10,8.2e10,8.2e10,8.2e10   
mpdata,ex,1,7,8.2e10,8.2e10,8.52e10,8.52e10,8.52e10
mpdata,alpx,1,1,4.8e-10,5.1e-10,5.85e-10,6.27e-10,5.35e-10,6.12e-10   
mpdata,alpx,1,7,6.00e-10,5.71e-10,5.62e-10,5.56e-10,5.42e-10
mp,prxy,1,1,0.17                     
MPDATA,DENS,1,1,2200                  
block,0,length,0, width, 0,height     
esize,lsize                          
mshkey,1
vatt,1,1,1  
vmesh,all
type,2      
asel,s,,,1,6
amesh,all
eplot
Save
*dim,flux2,table,max_x,max_y,max_time,x,y,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_y,1   
   flux2(0,j,k)=(j-1)*lsize  
  *enddo

  *do,i,1,max_x,1  
    *do,j,1,max_y,1  
     xcenter=vol*(k-1)*tinc     
     ycenter=0            
     distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
    *if,distance,le,rad,then
        flux2(i,j,k)=power_max/exp(kk*distance**2/(rad**2))  
    *else
        flux2(i,j,k)=0
    *endif  
    *enddo
  *enddo
  flux2(0,0,k)=(k-1)*tinc
*enddo
save
FINISH
/SOL
antype,trans      
timint,on         
autots,off        
toffst,20        
tunif,20         
kbc,1            
esel,s,type,,2   
nsle,s
sf,all,conv,10,25
   
allsel  
*SET,tim,1e-12      
*do,I,1,MAX_TIME,1  
*SET,tim,tim+tinc   
time,tim
autots,off  
deltim,,,,1
kbc,1   
tsres,erase
Nsel,s,loc,z,height,height
Nsel,r,loc,x,0,length     
Nsel,r,loc,y,0,0.02      
Sf,all,hflux,%flux2%
allsel,all
asel,s,,,2            
sfa,2,,hflux,%flux2%  
Allsel,all  
Solve   
*enddo
save
FINISH
/POST1
set,last        
plnsol,temp
请教各位高手，看看我的命令流有什么错误的地方，我模拟的图总感觉不对，待会上图

srwszsc

怎么都不见人指教呀

srwszsc

这是我的模拟图片焊缝温度过高，而且等温曲线也不对，请教高手指点

zxw19851015

激光焊接采用椭球或锥体热源,高斯穿透能力有限

srwszsc

谢谢zxw19851015 ,再请问有没有加载锥体的命令流？热源模型我已经学习了不少，但是命令流还是不行

zxw19851015

好久没有来论坛了，不知道你还在不