ANSYS固定TIG焊接的例子
!采用固定高斯热源TIG焊接
!最高温度4500°
!数组加载,热源加载在surf152上面
finish
/clear
/prep7
/UNITS,si !International system (m, kg, s, K)
et,1,solid70
ET,2,SURF152
KEYOPT,2,8,1
MPTEMP,1,27,127,327,527,727 !定义材料1 性能的温度表
MPDATA,DENS,MNUM,1,7920,7890,7810,7730,7640 !KG/m3 定义各温度下的密度
MPDATA,KXX,MNUM,1,16.0,16.5,19.0,22.5,25.7 !W/m·K定义各温度下的X上热传导系数
MPDATA,Kyy,MNUM,1,16.0,16.5,19.0,22.5,25.7 !定义各温度下的Y上热传导系数
MPDATA,Kzz,MNUM,1,16.0,16.5,19.0,22.5,25.7 !定义各温度下的Z上热传导系数
MPTEMP,1,293,773,873,973,4500
MPDATA,C,mnum,1,496.964,592.004,653,628,903.95 !J/Kg·K热容
MPTEMP,1,200,700,800,1000,1200
MPTEMP,6,1400,1454,1800,2000
MPDATA,HF,1,1,30,41.86,46.56,60.9,90.72
MPDATA,HF,1,6,129.09,141.06,235.86,306.44 !W/m2·K等效对流系数
!********建模********
rh=3e-3
alength=50e-3
bwidth=25e-3
hheight=4e-3
block,0,alength,0,rh,0,hheight
block,0,alength,rh,bwidth,0,hheight
VGLUE,1,2
esizel=0.001
b1num=rh/esizel
b2num=bwidth/(3*esizel)
a1num=alength/esizel
hnum=hheight/esizel
a1rito=1
b2rito=6
hrito=1
lesize,7,,,a1num,-1/a1rito
lesize,5,,,a1num,-1/a1rito
lesize,2,,,a1num,-1/a1rito
lesize,4,,,a1num,-1/a1rito
lesize,14,,,a1num,-1/a1rito
lesize,19,,,a1num,-1/a1rito
lesize,3,,,b1num
lesize,6,,,b1num
lesize,8,,,b1num
lesize,1,,,b1num
lesize,25,,,b2num,b2rito
lesize,26,,,b2num,b2rito
lesize,27,,,b2num,b2rito
lesize,28,,,b2num,b2rito
lesize,11,,,hnum
lesize,9,,,hnum
lesize,12,,,hnum
lesize,10,,,hnum
lesize,23,,,hnum
lesize,24,,,hnum
vmesh,all
ASEL,S,LOC,Z,hheight
NSLA,S,1
NSEL,R,LOC,Y,0,rh
type,2
ESURF
ALLS
NUMCMP,AREA
!*******固定高斯热源TABLE*********
pi=acos(-1)
U=15
I=60
ef=0.85
Qm=U*I*ef
qmax=3*Qm/(pi*rh**2)
max_x=alength/esizel+1
max_y=bwidth/esizel+1
*dim,flux2,table,max_x,max_y,,x,y
*do,i,1,max_x,1
flux2(i,0)=(i-1)*esizel !数组赋值方式:固定列按行赋值
*enddo !对第0列,第i行赋值
*do,j,1,max_y,1
flux2(0,j)=(j-1)*esizel
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=alength/2
ycenter=0
a=(i-1)*esizel-xcenter
b=(j-1)*esizel-ycenter
1