!非线性材料求解成功GAUSS分布Z向扫描!
!考虑了和基板对流情况下加载热流密度的热分析rb=0.15mm
!粉床7.8X3.3mm,烧结区域6X1.5mm
!考虑粉床表面的辐射、对流
!考虑粉末状、熔化状(液态)、凝固态的边界条件变化
!考虑金属蒸发会带走一部分能量,故效率减半eff=0.7
!r=0.2mm
FINISH
/CLEAR,START
/BATCH !设置程序模式为“批处理”
/CWD,'d:\Ansyswork\EBM FEM Simulation for Metals\3D Temperture distribution' !更改当前工作目录
/FILNAME,EBM_db,0
/TITLE,3D T emperture distrbution
!选择分析问题的类型:热
/NOPR
/PMETH,OFF,0
KEYW,PR_SET,1
KEYW,PR_STRUC,0
KEYW,PR_THERM,1
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
KEYW,PR_CFD,0
/GO
/UNITS,SI !Standard Units
!************************************************************************************** ***********!
!开始,进入前处理器,定义下列属性-
!0参数定义
!1定义单元类型
!2实常数
!3材料性能
!4几何模型
!5网格划分控制
!6其他
!************************************************************************************** ***********!
/PREP7
!0参数定义
!模型参数
!几何模型尺寸me-3=mm
*SET,Block_x,7.8e-3
*SET,Block_x1,0.9e-3
*SET,Block_x2,6.9e-3
*SET,Block_y,3.3e-3
*SET,Block_y1,0.9e-3
*SET,Block_y2,2.4e-3
*SET,Block_z,0.3e-3
*SET,Block_z1,0.9e-3
!划分网格的大小
*SET,LMeshx,3e-4 !line1
*SET,LMeshy,3e-4 !line4
*SET,LMeshz,3e-4 !line9
!边界条件参数
*SET,Tinit,30 !初始温度!1定义单元类型
ET,1,SOLID70
ET,2,surf152 !表面效应单元r,2
KEYOPT,2,4,1 !表面效应单元设置KEYOPT,2,8,4
ET,3,SOLID87
!2实常数
!Nothing
!3材料性能
MPTEMP
MPTEMP,1,25,400,1200,1600,2000,2500
MPDATA,DENS,1,1,3983,4688.4,5406.8,5948.8,6908,6331 !密度MPDATA,KXX,1,1,1.6278,2.3812,3.6396,29.8,30.6,32 !热导率MPDATA,c,1,1,472,535,679,859,945,1034 !比热
MPTEMP,,,,,,,,
MPTEMP,1,0,200,400,800,1200
MPTEMP,6,1400,1454,1800,2000
MPDATA,ENTH,1,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9 MPDATA,ENTH,1,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10 !热焓
MPTEMP,,,,,,,,
MPTEMP,1,20,300,600,900,1200,1500,1800,2000,3000,4000 MPDATA,HF,1,1,6,50,120,200,250,378,700,850,1746,3799 !对流系数
MPTEMP,,,,,,,,
MPTEMP,1,20,200,400,600,800,900,1000,1100,1200,1420,1460
MPDATA,DENS,2,1,7966,7893,7814,7724,7630,7583,7535,7486,7436,7320,7320
MPDATA,KXX,2,1,13.31,16.33,19.47,22.38,25.07,26.33,27.53,28.67,29.76,31.95,32 !热导率
mpdata,c,2,1,470,508,550,592,634,655,676,698,719,765,765
MPTEMP,,,,,,,,
MPTEMP,1,0,200,400,800,1200
MPTEMP,6,1400,1454,1800,2000
MPDATA,ENTH,2,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9
MPDATA,ENTH,2,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10
!4几何模型
BLOCK,0,Block_x,0,Block_y,0,Block_z,
BLOCK,Block_x1,Block_x2,Block_y1,Block_y2,0,Block_z,
vsel,all
VOVLAP,all
BLOCK,0,Block_x,0,Block_y,0,-Block_z1,
vglue,all
numcmp,all
!划分网格
lesize,8,LMeshx lesize,5,LMeshy lesize,13,LMeshz vsel,s,,,1 type,1
mat,1 mshkey,1 vmesh,all
allsel
esize,0.3e-3 mshkey,0 mshape,1,3d vsel,s,,,2 type,3
mat,2 vmesh,all
allsel
esize,0.3e-3 mshkey,0 mshape,1,3d