激光熔覆ANSYS仿真

!上表面上没有对流换热边界条件!单位制:米、秒、摄氏度!/CLEAR,START/FILNAME,temp,0/COM,ANSYSRELEASE10.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.5e9MPDATA,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.7MPDATA,KXX,1,,51.5MPTEMP,,MPTEMP,1,0MPTEMP,2,200MPTEMP,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!层数Height_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)/Velocity!载荷持续时间(扫描一层)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,,,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,,,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!基底左端面，对流换热边界条件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/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,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-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!激活单元的上表面，指定为对流换热边条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,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!激活单元的上表面，指定为对流换热边条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*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!载荷步内载荷随时间