ANSYS固定TIG焊接的例子

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