高斯热源 & 生死技术 & APDL

huangfuminkyjy

各位同仁：我用生死网格技术和高斯热源建立一个简单模型，是APDL的，麻烦大家大家帮忙看一下，谢谢了。


/prep7                                             !进入前处理
/unit,SI                                           !指定单位制
et,1,plane13                                       !单元类型选择

               %几何尺寸
length=0.01                                        !长度
width=0.01                                         !宽度
blc4,0,0,length,width,0                            !矩形框


             %材料参数
mat,1                                             
mp,dens,1,2520                                     !密度
mp,kxx,1,1.03                                      !热导率
mp,C,1,800                                         !热熔

             %网格划分
rb=0.0001                                          !单元尺寸值        
esize,rb,0                                         !单元尺寸设定
eshape,0,2D                                        !单元形状
mshkey,0                                           !自由划分
allsel                                             !选择所有实体
amesh,all                                          !划分网格


            %求解
/solu

                       %高斯光束=======开始======                              

              %激光参数
p=40                                              !激光能量（W）
r=0.0035                                          !光斑直径（m）
qmax=p/(3.14159*r**2)                             !峰值功率（W/m2）
A=0.03                                            !吸收率
              %光斑中心
xcenter=length                                    !中心点的x坐标
ycenter=width/2                                   !中心点的x坐标
              %表维数
max_x=length/rb+1                                 !x向最大维数
max_y=width/rb+1                                  !y向最大维数                     
*dim,flux2,table,max_x,max_y,,x,y                 !定义表
              %表头
*do,i,1,max_x,1                                   !x表头
flux2(i,0)=(i-1)*rb
*enddo
*do,j,1,max_y,1                                   !y表头
*flux2(0,j)=(j-1)*rb
*enddo
              %表元素
*do,i,1,max_x,1
*do,j,1,max_y,1  
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2)   !各个单元距离中心的距离
*if,distance,le,r,then                                       !半径内
flux2(i,j)=A*qmax*exp(-distance**2/r**2)                     !r为光强为1/e时的半径
*else                                                        !半径外
flux2(i,j)=0
*endif
*enddo
*enddo
                        %高斯光束=======结束======

        %时间总时间及时间步长
tall=0.5                                            !总时间
nt=30                                               !载荷步数
tsize=tall/nt                                       !时间子步

       %输出及求解设置                           
antype,4                                            !瞬态分析
autos,on                                            !自动时间步长跟踪
outres,,all                                         !输出控制
nropt,full                                          !牛顿-拉夫逊
kbc,1                                               !阶跃加载(貌似没有用)
        %初始温度和边界条件加载
tunif,20                                            !初始温度
nsel,s,loc,y,width                                  !选择上边界
sf,all,hflux2,%flux2%                               !热流加载

                              %循环求解及单元杀死=====开始=====
*do k,1,nt,1
time,k*tsize                                       !每次循环的时间                                                
nsel,all
solve                                               !求解

                               %杀死单元==开始==
Tr=3535                                             !材料去除温度
          %后处理-数据读取及单元选择
/post1
set,last                                            !读取最后一子步的数据
etable,melty,temp                                   !读取其温度值，并存在表melty中
esel,s,etab,melty,Tr                                !选择了温度大于Tr的单元格
finish

        %返回求解
/solu
ekill,all                                           !单元杀死
                               %杀死单元==结束==

antypes,,rest                                        !重复之前的设置/前一次分析的重新启动                                                         
*enddo
                             %循环求解及单元杀死=====结束=====


%