ansys.txt

焊接移动热源程序

!v=1mm/s,power=60w,r=100um
!time=5s,tinc=esize/v=0.1s
tm=5
tinc=1
esize=0.1e-3
v=0.001
R=0.1e-3
pi=3.1415926
!max_time=time/tinc+1=51
!max_x=x/esize+1=2e-3/0.1e-3+1=21
!max_y=y/esize+1=1e-2/0.1e-3+1=101
max_time=51
max_x=21
max_y=101
power=40
power_max=40/(pi*R*R)

/units,si
/prep7
et,1,solid70
mp,kxx,1,157  ! 定义si的材料属性
mp,dens,1,2.4e3 
mp,c,1,700  

mp,kxx,2,1.09  ! 定义glass(pyrex7740)的材料属性
mp,dens,2,2.23e3 
mp,c,2,800   !不知道对不对？

blc4,0,0,0.002,0.01,0.5e-3 !定义si实体模型
blc4,0.002,0,0.004,0.01,0.5e-3
blc4,-0.004,0,0.004,0.01,0.5e-3

block,0,0.002,0,0.01,0.5e-3,1e-3 !定义glass实体模型
block,0.002,0.006,0,0.01,0.5e-3,1e-3
block,-0.004,0,0,0.01,0.5e-3,1e-3

vglue,all
numcmp,all

esize,0.1e-3
mat,1
mshkey,1
vmesh,1

esize,0.1e-3
mat,2
mshkey,1
vmesh,2

esize,0.5e-3
mat,1
mshkey,1
vmesh,3,4

esize,0.5e-3
mat,2
mshkey,1
vmesh,5,6

*DIM,laser,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
   *DO,I,1,MAX_X,1
     laser(I,0,K)=(I-1)*ESIZE
   *ENDDO
   *DO,J,1,MAX_Y,1
     laser(0,J,K)=(J-1)*ESIZE
   *ENDDO
  *DO,I,1,MAX_X,1
  *DO,J,1,MAX_Y,1
     XCENTER=1e-3
     YCENTER=V*(K-1)*TINC+2e-3
     DISTANCE=SQRT(((I-1)*ESIZE-XCENTER)**2+((J-1)*ESIZE-YCENTER)**2)
       *IF,DISTANCE,LE,R,THEN
         laser(I,J,K)=power_max*EXP(-DISTANCE**2/(R**2)) ! 激光的表面热流
       *ELSE
         laser(I,J,K)=0
       *ENDIF
  *ENDDO
  *ENDDO
     laser(0,0,K)=(K-1)*TINC
*ENDDO
/solu
antype,trans
timint,on
autots,off
tunif,27
kbc,1 !step load
asel,s,,,1
asel,s,,,14
asel,s,,,18
sfa,all,,hflux,-200000  !在体的下表面加负的热流密度，模拟冷却
allsel
asel,s,,,23  
asel,s,,,27
asel,s,,,9
asel,s,,,3,4
asel,s,,,7,8
asel,s,,,10
asel,s,,,12
asel,s,,,16,17
asel,s,,,20,22,1
asel,s,,,24,26,1
asel,s,,,28,29
sfa,all,,conv,110,27 !在体的上表面和周边加对流
allsel

*do,t,1e-6,tm,tinc
time,t
autots,off
deltim,,,,1 !Use final time step from previous load step 
kbc,1
tsres,erase
sfa,2,,hflux,%laser% 
solve 
outres,all,all 
*enddo 
outres,all,all 
finish