激光熔覆ANSYS仿真

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

2、TA,ENTH,1,7.2e9MP DATA,ENTH,1,8e9MP DATA,ENTH,1,9e9MP DATA,ENTH,1,9.9e9MP DATA,ENTH,1,12.84e9MP TEMP,MP TEMP,1,0MP TEMP,2,200MP TEMP,3,400MP TEMP,4,800MP TEMP,5,1440MP TEMP,6,2000MP DATA,KXX,1,14.7MP DATA,KXX,1,18MP DATA,KXX,1,20.8MP DATA,KXX,1,26.3MP DATA,KXX,1,34.7MP DATA,KXX,1,51.5MP TEMP,MP TEMP

3、,1,0MP TEMP,2,200MP TEMP,3,400MP TEMP,4,800MP TEMP,5,1200MP TEMP,6,2000MP DATA,C,1,494MP DATA,C,1,536MP DATA,C,1,569MP DATA,C,1,644MP DATA,C,1,669MP DATA,C,1,729MP TEMP,MP TEMP,1,0MP DATA,DENS,1,7700!1Width_Base=0.025!宽度Height_Base=0.02!基底高度Len gth=0.05!长度Width_Clad=0.0015!宽度Height_De positio n=0.00

4、25!覆层高度Layer=10!层数Height_Clad=Height_De positio n/layerdt=0.00001!小量Temp=50!环境温度Ini t_Te mp=200!初始温度Coff_Co nv=301 -!对流换热系数1Velocity=0.003!扫描速度Ste pDis=0.001!每个载荷步位移Laser Po wer=700!激光功率Radius=0.0015!激光光斑半径定义常量定义常量!激光光斑面积!吸收因子Area=3.14159265*(Radius*2)Factor=0.25TotalTime=(Le ngth+Radius*2)/Velocity

5、Ste pNum=(Le ngth+Radius*2)/Ste pDis=Ste pTime=Ste pDis/Velocity!每个载荷步时间!载荷持续时间（扫描一层）!载荷步数建模BLOCK,0,Le ngth,0,Height_De positio n,0,Width_Clad,BLOCK,0,Le ngth,0,-Height_Base,0,Width_Clad,BLOCK,0,Le ngth,0,-Height_Base,Width_Clad,Width_BaseVGLUE,ALLNUMCMP,ALL!覆层!基底!压缩对象编号划分网格LESIZE,18,15,0.125LESIZE,

6、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_De positio n-dt,LESIZE, ALL, , , Layer ,!覆层高度方向的单元数目LSEL, S, LOC, X, dt, Le ngth-dt,LESIZE, ALL, , , Le ngth/Ste pDis,!长度方向的单元数目LSEL, S, LOC,乙 dt, Width_

7、Clad-dt,LESIZE, ALL, , , 5,!宽度方向的单元数目VSEL,ALL!网格划分TYP E,1MAT,1MSHA PE,0,3DMSHKEY1VMESH,ALL11基底边界条件、初始条件NSEL, S, LOC, Y, -Height_Base, 0 IC,ALL,TEM P,lnit_Temp ALLSEL,ALL!基底初始温度NSEL, S, LOC, Y, -Height_Base, 0 NSEL, R, LOC,乙 Width_Base SF ALL, CONV, Coff_Co nv, Temp ALLSEL,ALL!基底侧面，对流换热边界条件NSEL, S, L

8、OC, Y, -Height_Base, 0NSEL, R, LOC, X, 0SF ALL, CONV, Coff_Co nv, Temp ALLSEL,ALL!基底左端面，对流换热边界条件NSEL, S, LOC, Y, -Height_Base, 0NSEL, R, LOC, X, Le ngthSF ALL, CONV, Coff_Co nv, Temp ALLSEL,ALL!基底右端面，对流换热边界条件NSEL, S, LOC, Y, 0NSEL, R, LOC,乙 Width_Clad, Width_Base SF ALL, CONV, Coff_Co nv, Temp ALLSE

9、L,ALL!基底上表面，对流换热边界条件FINISH/SOLU11瞬态分析参数设置ANT YP E,4!分析类型：瞬态!*TRN OP T,FULL!求解方法：完全的 N-R方法LUMP M,01*!不使用集中质量矩阵!NSUBST, 1!指定载荷子步数目DELTIM,0.01,0.001,0.05!载荷子步（默认子步时间步长、 最小、最大）载荷步为0.333CNVTOL,HEAT, ,0.01,1,0.000001,!收敛准则：控制热流OUTRES,NSOL,LAST!结果输岀：所有!:1杀死单元!杀死熔覆层单元NSEL, S, LOC, Y , 0, Height_DepositionES

10、LN, S, 1, ALLEKILL,ALLALLSEL,ALLESEL,S,LIVEE P LOTESEL,S,LIVENSLE,S,1NSEL,R,LOC,Y,HeightBaseNSEL,R,LOC, Z, 0,RadiusSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALL1!激活单元的上表面，指定为对流换热边条1预热*DO, i, 1, 2m=mod(i,2)*IF,m,EQ,1,THEN!如果为奇数层，向右扫描*DO, k, 1, Step Num, 1TIME,TotalTime*(i-1)+Ste pTime*kKBC, 1!载荷步结束时间!载荷步内载荷随

11、时间分布：常数LeftX=Ste pDis*(k-1)RightX=Ste pDis*kNSEL, S, LOC, Y , 0ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC,乙 0, Radius!激活单元的上表面，加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , Laser Po werDO, k, 1, Step Num, 1Factor/Area,ALLSEL, ALLSOLVE!删除热流密度载荷SFEDELE,ALL,4,HFLUXALLS

12、EL, ALLESEL,S,LIVE!激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,0NSEL,R,LOC,x,LeftX-2*Ste pDis,RightX-2*Ste pDisSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALLESEL,S,LIVE EP LOT*ENDDO*ELSE!如果为偶数层，向左扫描TIME,TotalTime*(i-1)+Ste pTime*kKBC, 1!载荷步结束时间!载荷步内载荷随时间分布：常数LeftX=Le ngth-Ste pDis*kRightX=Le ngth-Ste pDis*(k-1)NS

13、EL, S, LOC, Y , 0ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, RightX-2*Radius, RightX!激活单元的上表面，加热流密度NSEL, R, LOC,乙 0, RadiusESLN, S, 1SFE, ALL, 4, HFLUX, , Laser Po wer*Factor/Area,ALLSEL, ALLSOLVE!激活单元的上表面，删除载荷!激活单兀的上表面， 指定为对流换热边条SFEDELE,ALL,4,HFLUXALLSEL, ALLESEL,S,LIVENSLE,S,1NSEL,R,LOC,Y,0N

14、SEL, R, LOC, x, LeftX-2*Ste pDis, RightX-2*Ste pDisSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALLESEL,S,LIVEE P LOT*ENDDO*ENDIF*ENDDOeselsliveE P LOT熔覆*DO, i, 1, Layer , 1m=mod(i,2)TF,m,EQ,1,THEN!如果为奇数层，向右扫描TIME,TotalTime*(i+1)+Ste pTime*kKBC, 1!载荷步结束时间!载荷步内载荷随时间分布：常数LeftX=Ste pDis*(k-1)RightX=Ste pDis*kNSEL

15、, S, LOC, Y , Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i NSEL, R, LOC, x, LeftX, RightX!激活单元ESLN, S, 1 EALIVE,ALLALLSEL, ALLNSEL, S, LOC, Y , Height_Deposition/Layer*(i-1), Height_Deposition/Layer*iNSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC,乙 0, RadiusESLN, S, 1SFE, ALL, 4, H

16、FLUX, , Laser Po wer*Factor/Area, , , !激活单元的上表面，加热流密度ALLSEL, ALLESEL,S,LIVE!激活单元的表面，如果包含左端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,Height_De positio n/Layer*(i-1), Height_De positio n/Layer*iSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALLESEL,S,LIVE!激活单元的表面，如果包含右端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,Le ngthNS

17、EL,R,LOC,Y,Height_De positio n/Layer*(i-1), Height_De po sitio n/Layer*iSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALL!激活单元的侧面，指定为对流换热边条ESEL,S,LIVENSLE,S,1NSEL,R,LOC, Z, Width_CladNSEL,R,LOC,Y,Height_De positio n/Layer*(i-1), Height_De positio n/Layer*iSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALLSOLVESFEDELE,ALL,4,

18、HFLUX!删除热流密度载荷ALLSEL, ALLESEL,S,LIVE!激活单兀的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,Height_De positi on/Layer*iNSEL,R,LOC,x,LeftX-2*Ste pDis,RightX-2*Ste pDisSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALLNSEL,S,LOC,Y,Height_De po sitio n/Layer*(i-1)ESEL, S, LIVE!激活单元的下表面， 删除对流换热边条ESLN, R, 0NSLE, S, 1NSEL, R, LOC, x,

19、 LeftX-2*Ste pDis,RightX-2*Ste pDisNSEL, R, LOC,乙 0, RadiusNSEL, U, LOC, Y , Height_Deposition/Layer*iESLN, S, 1SFDELE, ALL, CONVALLSEL,ALLESEL,S,LIVEE P LOT*ENDDO*ELSE!如果为偶数层，向左扫描*DO, k, 1, Step Num, 1TIME,TotalTime*(i+1)+Ste pTime*kKBC, 1!载荷步结束时间!载荷步内载荷随时间分布：常数LeftX=Le ngth-Ste pDis*kRightX=Le ngt

20、h-Ste pDis*(k-1)NSEL, S, 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*iNSEL, R, LOC, x, LeftX, LeftX+2*Radius!激活单元的上表面，加热流密度NSEL, R, LOC,乙 0, Radi

21、usESLN, S, 1SFE, ALL, 4, HFLUX, , Laser Po wer*Factor/Area,ALLSEL, ALLESEL,S,LIVE!激活单元的表面，如果包含左端面，指定为对流换热边条NSLE,S,1NSELRLOCXONSEL,R,LOC,Y,Height_De positio n/Layer*(i-1), Height_De po sitio n/Layer*iSFALL,CONV,Coff_Co nv,Te mpALLSEL,ALLESELSLIVE!激活单元的表面，如果包含右端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,Le ngthNSEL,R,LOC,Y,Height_De positio n/Layer*(i-1), Height_De positio n/Layer*iSFALL,CONV,Coff_Co