/prep7

/units,si

mptemp,1,20,250,500,750,1000,1500

mptemp,7,1700,2500

mpdata,kxx,1,1,50,47,40,27,30,35

mpdata,kxx,1,7,140,142

mpdata,ex,1,1,2.05e+11,1.87e+11,1.50e+11,7.00e+10,2.00e+10,1.90e+10

mpdata,ex,1,7,1.80e+10,1.20e+10

mpdata,dens,1,1,7820,7700,7610,7550,7490,7350

mpdata,dens,1,7,7300,7090

mpdata,c,1,1,460,480,530,675,670,660

mpdata,c,1,7,780,820

mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.4,0.45

mpdata,prxy,1,7,0.48,0.5

mpdata,alpx,1,1,1.1e-05,1.22e-05,1.39e-05,1.48e-05,1.34e-05,1.33e-05,1.32e-05,1.31e-05

mpdata,alpx,1,7,1.32e-05,1.31e-05

mptemp,,,,,,,,

mptemp,1,20,1533,1590,1670

mpdata,enth,1,1,0,7.5e+009,9.6e+009,1.1e+010

rb=0.002                                        ！單元尺寸     

length=0.12                                    ！長方體幾何參數

width=0.06

height=0.006

et,1,solid70

block,0,length,0,width,0,height

Block,0,length,0,3*Rb,0,height

Block,0,length,0,6*Rb,0,height

Block,0,length,0,12*Rb,0,height

VOVLAP,ALL

Lesize,4,rb/2                        ！mapped劃分單元尺寸

lesize,5,rb/2

lesize,19,rb/2

lesize,14,rb/2

lesize,50,rb/2

lesize,10,rb/2

lesize,52,rb/2

lesize,23,rb/2

lesize,51,rb/2

lesize,9,rb/2

lesize,49,rb/2

lesize,24,rb/2

lesize,31,rb

lesize,26,rb

lesize,43,4*rb

lesize,38,4*rb

lesize,2,4*rb

lesize,7,4*rb

lesize,53,4*rb

lesize,54,4*rb

lesize,55,4*rb

lesize,56,4*rb

lesize,11,0.003

lesize,12,0.003

lesize,47,0.003

lesize,48,0.003

lesize,35,0.003

lesize,36,0.003

lesize,64,0.004

lesize,61,0.004

lesize,62,0.004

lesize,63,0.004

lesize,60,0.002

lesize,57,0.002

lesize,58,0.002

lesize,59,0.002

type,1

mshkey,1

mshape,0

vmesh,5

vmesh,6

type,1

mshkey,0

mshape,1

vmesh,7

vmesh,8

/solu                ！求解設置

antype,trans

timint,on,therm

tref,25

nropt,full,

kbc,1

autots,on

*do,i,1,60,1   ！高斯移動熱源加載

b=0.2*i            ！時間步

time,b

deltim,0.2,0.01,0.2

xcenter=i*0.002 ！加載中心點

*DEL,_FNCNAME   

*DEL,_FNCMTID   

*DEL,_FNCCSYS   

*SET,_FNCNAME,'flux2'   

*SET,_FNCCSYS,0 

! /INPUT,C:\Users\Administrator\Desktop\gaosi.func,,,1  

*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%  

!   

! Begin of equation: 89524655/exp(222222*({x}-0.01*{time})^2+({y})^2+   

! ({z}-0.006)^2)

*SET,%_FNCNAME%(0,0,1), 0.0, -999   

*SET,%_FNCNAME%(2,0,1), 0.0 

*SET,%_FNCNAME%(3,0,1), 0.0 

*SET,%_FNCNAME%(4,0,1), 0.0 

*SET,%_FNCNAME%(5,0,1), 0.0 

*SET,%_FNCNAME%(6,0,1), 0.0 

*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.01, 0, 0, 1   

*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 1 

*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 1, 2, 2, -2 

*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 2, 0, 0, -1 

*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 17, -2   

*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 222222, 0, 0, -3

*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, -1, 3, -3

*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, 3  

*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 3, 17, -1

*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -2, 1, -3   

*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 0.006, 0, 0, 4 

*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, 4, 2, -2

*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, -3

*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -3, 17, -2  

*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, -1, 1, -4   

*SET,%_FNCNAME%(0,16,1), 0.0, -1, 7, 1, -2, 0, 0

*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 89524655, 0, 0, -1 

*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -2, 4, -1   

*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -3, 0, 0

! End of equation: 89524655/exp(222222*({x}-0.01*{time})^2+({y})^2+({z}-0.006)^ 

! 2)

allsel,all                            ！選取加載中心點半徑6mm內的點施加熱流密度

nsel,s,loc,z,0.006

nsel,r,loc,x,xcenter-0.006,xcenter+0.006

nsel,r,loc,y,0,0.006

local,11,1,xcenter,0,0.006

csys,11

nsel,r,loc,x,0,0.006

csys,0

sf,all,hflux,%flux2%allsel,all            ！高度方向4個面施加對流載荷

nsel,s,loc,x,0

sf,all,conv,10,20

allsel,all

nsel,s,loc,x,0.12

sf,all,conv,10,20

allsel,all

nsel,s,loc,y,0.06

sf,all,conv,10,20

allsel,all

nsel,s,loc,y,0

sf,all,conv,10,20

allsel,all

solve

*enddo