激光熔覆ANSYS仿真

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

2、PDATA,ENTH,1,7.2e9MPDATA,ENTH,1,8e9MPDATA,ENTH,1,9e9MPDATA,ENTH,1,9.9e9MPDATA,ENTH,1,12.84e9MPTEMP,MPTEMP,1,0MPTEMP,2,200MPTEMP,3,400MPTEMP,4,800MPTEMP,5,1440MPTEMP,6,2000MPDATA,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,20

3、0MPTEMP,3,400MPTEMP,4,800MPTEMP,5,1200MPTEMP,6,2000MPDATA,C,1,494MPDATA,C,1,536MPDATA,C,1,569MPDATA,C,1,644MPDATA,C,1,669MPDATA,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!层数Heigh

4、t_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.14159265*(Radius*2)! 激光光斑面积Factor=0.25!吸收因子StepTime=StepDis/Velocity! 每个载荷步时间TotalTime=(Length+Radius*2)/Vel

5、ocity!载荷持续时间(扫描一层 )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.125LESIZE,22,15,8LESIZE,24,15,8LESIZE,26

6、,15,8LESIZE,28,15,8LESIZE,20,15,8LESIZE,21,15,8LESIZE,23,15,8LESIZE,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, ALL,

7、, , 5,VSEL,ALL! 网格划分.TYPE,1MAT,1MSHAPE,0,3DMSHKEY,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,ALL

9、FINISH/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.000001,!收敛准则：控制热流OUTRES,NSOL,LAST! 结果输出：所有!=!杀死单元NSEL, S, LOC, Y , 0, Height_Deposition!杀死熔覆层单元ESLN, S, 1, ALLEKILL,

10、ALLALLSEL,ALLESEL,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, 1TIME,TotalTime*(i-1)+StepTime*k! 载荷步结束时间KBC, 1!载荷步内载荷随时间分布：常数LeftX=StepDis*(k

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

12、NSLE,S,1NSEL,R,LOC,Y,0NSEL,R,LOC,x,LeftX-2*StepDis,RightX-2*StepDisSF,ALL,CONV ,Coff_Conv,TempALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ELSE! 如果为偶数层，向左扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i-1)+StepTime*k!载荷步结束时间KBC, 1!载荷步内载荷随时间分布：常数LeftX=Length-StepDis*kRightX=Length-StepDis*(k-1)NSEL, S, LOC, Y , 0ESEL, S

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

14、 LeftX-2*StepDis, RightX-2*StepDisSF,ALL,CONV ,Coff_Conv,TempALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ENDIF*ENDDOESEL,S,LIVEEPLOT!=!熔覆*DO, i, 1, Layer, 1m=mod(i,2)*IF ,m,EQ,1,THEN!如果为奇数层，向右扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i+1)+StepTime*k!载荷步结束时间KBC, 1!载荷步内载荷随时间分布：常数LeftX=StepDis*(k-1)RightX=StepDis*kN

15、SEL, S, LOC, Y , Height_Deposition/Layer*(i-1), Height_Deposition/Layer*iNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALL!激活单元ALLSEL, ALL.NSEL, S, LOC, Y , Height_Deposition/Layer*(i-1), Height_Deposition/Layer*iNSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC, Z, 0, RadiusESLN, S, 1SFE, ALL,

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

17、SEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i SF,ALL,CONV ,Coff_Conv,TempALLSEL,ALLESEL,S,LIVE! 激活单元的侧面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Z,Width_CladNSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*iSF,ALL,CONV ,Coff_Conv,TempALLSEL,ALLSOLVESFEDELE,ALL,4,HFLUX!删除

18、热流密度载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,Height_Deposition/Layer*iNSEL,R,LOC,x,LeftX-2*StepDis,RightX-2*StepDisSF,ALL,CONV ,Coff_Conv,TempALLSEL,ALLNSEL,S,LOC,Y,Height_Deposition/Layer*(i-1)!激活单元的下表面，删除对流换热边条ESEL, S, LIVEESLN, R, 0.NSLE, S, 1NSEL, R, LOC, x, LeftX-2*StepDi

19、s,RightX-2*StepDisNSEL, R, LOC, Z, 0, RadiusNSEL, U, LOC, Y , Height_Deposition/Layer*iESLN, S, 1SFDELE, ALL, CONVALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ELSE!如果为偶数层，向左扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i+1)+StepTime*k!载荷步结束时间KBC, 1!载荷步内载荷随时间分布：常数LeftX=Length-StepDis*kRightX=Length-StepDis*(k-1)NSEL, S

20、, LOC, Y , Height_Deposition/Layer*(i-1), Height_Deposition/Layer*iNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLNSEL, S, LOC, Y , Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i NSEL, R, LOC, x, LeftX, LeftX+2*RadiusNSEL, R, LOC, Z, 0, Radius ! 激活单元的上表面，加热流密度 ESLN, S, 1SFE,

21、ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLESEL,S,LIVE! 激活单元的表面，如果包含左端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*iSF,ALL,CONV ,Coff_Conv,TempALLSEL,ALL.ESEL,S,LIVE!激活单元的表面，如果包含右端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,LengthNSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*iSF,ALL,CON