激光熔覆ANSYS仿真

1、!上表面上没有对流换热边界条件!单位制:米、秒、摄氏度!/CLEAR,START/FILNAME,temp,0/COM,ANSYS RELEASE 10.0 UP 20:15:52 07/14/2007/CONFIG, NRES, 5000/PREP7!* !=!指定单元 ET,1,SOLID70!*!*!=!材料属性MPTEMP,MPTEMP,1,0MPTEMP,2,1350 MPTEMP,3,1400MPTEMP,4,1420 MPTEMP,5,1440MPTEMP,6,1490 MPTEMP,7,2000 MPDATA,ENTH,1,0MPDATA,ENTH,1,6.5e9MPDATA,

2、ENTH,1,7.2e9MPDATA,ENTH,1,8e9MPDATA,ENTH,1,9e9MPDATA,ENTH,1,9.9e9 MPDATA,ENTH,1,12.84e9MPTEMP,MPTEMP,1,0MPTEMP,2,200MPTEMP,3,400MPTEMP,4,800MPTEMP,5,1440MPTEMP,6,2000 MPDATA,KXX,1,14.7MPDATA,KXX,1,18MPDATA,KXX,1,20.8MPDATA,KXX,1,26.3MPDATA,KXX,1,34.7 MPDATA,KXX,1,51.5MPTEMP,MPTEMP,1,0MPTEMP,2,200MPT

3、EMP,3,400MPTEMP,4,800MPTEMP,5,1200MPTEMP,6,2000 MPDATA,C,1,494 MPDATA,C,1,536 MPDATA,C,1,569 MPDATA,C,1,644 MPDATA,C,1,669 MPDATA,C,1,729MPTEMP,MPTEMP,1,0MPDATA,DENS,1,7700!=!定义常量Width_Base=0.025 !宽度Height_Base=0.02 !基底高度Length=0.05 !长度Width_Clad=0.0015 !宽度Height_Deposition=0.0025 !覆层高度Layer=10 !层数H

4、eight_Clad=Height_Deposition/layerdt=0.00001 !小量Temp=50 !环境温度Init_Temp=200 !初始温度Coff_Conv=30 !对流换热系数!=!定义常量Velocity=0.003 !扫描速度StepDis=0.001 !每个载荷步位移LaserPower=700 !激光功率Radius=0.0015 !激光光斑半径Area=3.*(Radius*2) !激光光斑面积Factor=0.25 !吸收因子StepTime=StepDis/Velocity !每个载荷步时间TotalTime=(Length+Radius*2)/Veloc

5、ity !载荷持续时间(扫描一层)StepNum=(Length+Radius*2)/StepDis !载荷步数!=!建模BLOCK,0,Length,0,Height_Deposition,0,Width_Clad, !覆层BLOCK,0,Length,0,-Height_Base,0,Width_Clad, !基底BLOCK,0,Length,0,-Height_Base,Width_Clad,Width_BaseVGLUE,ALLNUMCMP,ALL ！压缩对象编号!=!划分网格LESIZE,18,15,0.125 LESIZE,22,15,8 LESIZE,24,15,8 LESIZE

6、,26,15,8LESIZE,28,15,8LESIZE,20,15,8 LESIZE,21,15,8LESIZE,23,15,8 LESIZE,25,15,8LESIZE,27,15,8LSEL, S, LOC, Y, dt, Height_Deposition-dt, !覆层高度方向的单元数目LESIZE, ALL, , , Layer,LSEL, S, LOC, X, dt, Length-dt, !长度方向的单元数目LESIZE, ALL, , , Length/StepDis,LSEL, S, LOC, Z, dt, Width_Clad-dt, !宽度方向的单元数目LESIZE,

7、ALL, , , 5,VSEL,ALL !网格划分TYPE,1MAT,1MSHAPE,0,3D MSHKEY,1VMESH,ALL!=!基底边界条件、初始条件NSEL, S, LOC, Y, -Height_Base, 0 !基底初始温度IC,ALL,TEMP,Init_TempALLSEL,ALLNSEL, S, LOC, Y, -Height_Base, 0 !基底侧面，对流换热边界条件NSEL, R, LOC, Z, Width_BaseSF, ALL, CONV, Coff_Conv, TempALLSEL,ALLNSEL, S, LOC, Y, -Height_Base, 0 !基底

8、左端面，对流换热边界条件NSEL, R, LOC, X, 0SF, ALL, CONV, Coff_Conv, TempALLSEL,ALLNSEL, S, LOC, Y, -Height_Base, 0 !基底右端面，对流换热边界条件NSEL, R, LOC, X, LengthSF, ALL, CONV, Coff_Conv, TempALLSEL,ALLNSEL, S, LOC, Y, 0 !基底上表面，对流换热边界条件NSEL, R, LOC, Z, Width_Clad, Width_BaseSF, ALL, CONV, Coff_Conv, TempALLSEL,ALLFINISH

9、/SOLU!=!瞬态分析参数设置ANTYPE,4 !分析类型：瞬态!*TRNOPT,FULL !求解方法：完全的N-R方法LUMPM,0 !不使用集中质量矩阵!*!NSUBST, 1 !指定载荷子步数目DELTIM,0.01,0.001,0.05 !载荷子步(默认子步时间步长、最小、最大)载荷步为0.333CNVTOL,HEAT, ,0.01,1,0., !收敛准则：控制热流OUTRES,NSOL,LAST !结果输出：所有!=!杀死单元NSEL, S, LOC, Y, 0, Height_Deposition !杀死熔覆层单元ESLN, S, 1, ALLEKILL,ALLALLSEL,AL

10、LESEL,S,LIVEEPLOTESEL,S,LIVE !激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,HeightBaseNSEL,R,LOC,Z,0,RadiusSF,ALL,CONV,Coff_Conv,TempALLSEL,ALL!=!预热*DO, i, 1, 2m=mod(i,2)*IF,m,EQ,1,THEN !如果为奇数层，向右扫描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i-1)+StepTime*k !载荷步结束时间 KBC, 1 !载荷步内载荷随时间分布：常数 LeftX=StepDis*(k-1) Ri

11、ghtX=StepDis*k NSEL, S, LOC, Y, 0 ESEL, S, LIVE ESLN, R, 0 NSLE, S, 1 NSEL, R, LOC, x, RightX-2*Radius, RightX NSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度 ESLN, S, 1 SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , , ALLSEL, ALL SOLVE SFEDELE,ALL,4,HFLUX !删除热流密度载荷 ALLSEL, ALL ESEL,S,LIVE !激活单元的上表面，指定为

12、对流换热边条 NSLE,S,1 NSEL,R,LOC,Y,0 NSEL,R,LOC,x,LeftX-2*StepDis,RightX-2*StepDis SF,ALL,CONV,Coff_Conv,Temp ALLSEL,ALL ESEL,S,LIVE EPLOT *ENDDO *ELSE !如果为偶数层，向左扫描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i-1)+StepTime*k !载荷步结束时间 KBC, 1 !载荷步内载荷随时间分布：常数 LeftX=Length-StepDis*k RightX=Length-StepDis*(k-1) NSEL

13、, S, LOC, Y, 0 ESEL, S, LIVE ESLN, R, 0 NSLE, S, 1 NSEL, R, LOC, x, RightX-2*Radius, RightX NSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度 ESLN, S, 1 SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , , ALLSEL, ALL SOLVE SFEDELE,ALL,4,HFLUX !激活单元的上表面，删除载荷 ALLSEL, ALL ESEL,S,LIVE !激活单元的上表面，指定为对流换热边条 NSLE,S,

14、1 NSEL,R,LOC,Y,0 NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDis SF,ALL,CONV,Coff_Conv,Temp ALLSEL,ALL ESEL,S,LIVE EPLOT *ENDDO *ENDIF*ENDDO ESEL,S,LIVEEPLOT!=!熔覆*DO, i, 1, Layer, 1m=mod(i,2)*IF,m,EQ,1,THEN !如果为奇数层，向右扫描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i+1)+StepTime*k !载荷步结束时间 KBC, 1 !载荷步内载荷

15、随时间分布：常数 LeftX=StepDis*(k-1) RightX=StepDis*k NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i NSEL, R, LOC, x, LeftX, RightX ESLN, S, 1 EALIVE,ALL !激活单元 ALLSEL, ALL NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i NSEL, R, LOC, x, RightX-2*Radius, R

16、ightX NSEL, R, LOC, Z, 0, Radius ESLN, S, 1 SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , , !激活单元的上表面，加热流密度 ALLSEL, ALL ESEL,S,LIVE !激活单元的表面，如果包含左端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,0 NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i SF,ALL,CONV,Coff_Conv,Temp ALLSEL,ALL ESEL,S,

17、LIVE !激活单元的表面，如果包含右端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,Length NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i SF,ALL,CONV,Coff_Conv,Temp ALLSEL,ALL ESEL,S,LIVE !激活单元的侧面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,Z,Width_Clad NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer

18、*i SF,ALL,CONV,Coff_Conv,Temp ALLSEL,ALL SOLVE SFEDELE,ALL,4,HFLUX !删除热流密度载荷 ALLSEL, ALL ESEL,S,LIVE !激活单元的上表面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,Y,Height_Deposition/Layer*i NSEL,R,LOC,x,LeftX-2*StepDis,RightX-2*StepDis SF,ALL,CONV,Coff_Conv,Temp ALLSEL,ALL NSEL,S,LOC,Y,Height_Deposition/Layer*(i-1) !激活单

19、元的下表面，删除对流换热边条 ESEL, S, LIVE ESLN, R, 0 NSLE, S, 1 NSEL, R, LOC, x, LeftX-2*StepDis,RightX-2*StepDis NSEL, R, LOC, Z, 0, Radius NSEL, U, LOC, Y, Height_Deposition/Layer*iESLN, S, 1 SFDELE, ALL, CONV ALLSEL,ALL ESEL,S,LIVE EPLOT *ENDDO *ELSE !如果为偶数层，向左扫描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i+1)+Ste

20、pTime*k !载荷步结束时间 KBC, 1 !载荷步内载荷随时间分布：常数 LeftX=Length-StepDis*k RightX=Length-StepDis*(k-1) NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i NSEL, R, LOC, x, LeftX, RightX ESLN, S, 1 EALIVE,ALL ALLSEL, ALL NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Laye

21、r*i NSEL, R, LOC, x, LeftX, LeftX+2*Radius NSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度 ESLN, S, 1 SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , , ALLSEL, ALL ESEL,S,LIVE !激活单元的表面，如果包含左端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,0 NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i SF,ALL,CONV,Coff_Conv,Temp ALLSEL,ALL ESEL,S,LIVE !激活单元的表面，如果包含右端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,Length NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i SF,ALL,CONV,Coff_Conv,