激光熔覆ANSYS仿真

1、!上表面上沒(méi)有對(duì)流換熱邊界條件!單位制:米、秒、攝氏度!/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 !長(zhǎng)度Width_Clad=0.0015 !寬度Height_Deposition=0.0025 !覆層高度Layer=10 !層數(shù)H

4、eight_Clad=Height_Deposition/layerdt=0.00001 !小量Temp=50 !環(huán)境溫度Init_Temp=200 !初始溫度Coff_Conv=30 !對(duì)流換熱系數(shù)!=!定義常量Velocity=0.003 !掃描速度StepDis=0.001 !每個(gè)載荷步位移LaserPower=700 !激光功率Radius=0.0015 !激光光斑半徑Area=3.*(Radius*2) !激光光斑面積Factor=0.25 !吸收因子StepTime=StepDis/Velocity !每個(gè)載荷步時(shí)間TotalTime=(Length+Radius*2)/Veloc

5、ity !載荷持續(xù)時(shí)間(掃描一層)StepNum=(Length+Radius*2)/StepDis !載荷步數(shù)!=!建模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 ！壓縮對(duì)象編號(hào)!=!劃分網(wǎng)格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, !覆層高度方向的單元數(shù)目LESIZE, ALL, , , Layer,LSEL, S, LOC, X, dt, Length-dt, !長(zhǎng)度方向的單元數(shù)目LESIZE, ALL, , , Length/StepDis,LSEL, S, LOC, Z, dt, Width_Clad-dt, !寬度方向的單元數(shù)目LESIZE,

7、ALL, , , 5,VSEL,ALL !網(wǎng)格劃分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 !基底側(cè)面，對(duì)流換熱邊界條件NSEL, R, LOC, Z, Width_BaseSF, ALL, CONV, Coff_Conv, TempALLSEL,ALLNSEL, S, LOC, Y, -Height_Base, 0 !基底

8、左端面，對(duì)流換熱邊界條件NSEL, R, LOC, X, 0SF, ALL, CONV, Coff_Conv, TempALLSEL,ALLNSEL, S, LOC, Y, -Height_Base, 0 !基底右端面，對(duì)流換熱邊界條件NSEL, R, LOC, X, LengthSF, ALL, CONV, Coff_Conv, TempALLSEL,ALLNSEL, S, LOC, Y, 0 !基底上表面，對(duì)流換熱邊界條件NSEL, R, LOC, Z, Width_Clad, Width_BaseSF, ALL, CONV, Coff_Conv, TempALLSEL,ALLFINISH

9、/SOLU!=!瞬態(tài)分析參數(shù)設(shè)置ANTYPE,4 !分析類型：瞬態(tài)!*TRNOPT,FULL !求解方法：完全的N-R方法LUMPM,0 !不使用集中質(zhì)量矩陣!*!NSUBST, 1 !指定載荷子步數(shù)目DELTIM,0.01,0.001,0.05 !載荷子步(默認(rèn)子步時(shí)間步長(zhǎng)、最小、最大)載荷步為0.333CNVTOL,HEAT, ,0.01,1,0., !收斂準(zhǔn)則：控制熱流OUTRES,NSOL,LAST !結(jié)果輸出：所有!=!殺死單元NSEL, S, LOC, Y, 0, Height_Deposition !殺死熔覆層單元ESLN, S, 1, ALLEKILL,ALLALLSEL,AL

10、LESEL,S,LIVEEPLOTESEL,S,LIVE !激活單元的上表面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,Y,HeightBaseNSEL,R,LOC,Z,0,RadiusSF,ALL,CONV,Coff_Conv,TempALLSEL,ALL!=!預(yù)熱*DO, i, 1, 2m=mod(i,2)*IF,m,EQ,1,THEN !如果為奇數(shù)層，向右掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i-1)+StepTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷隨時(shí)間分布：常數(shù) 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、對(duì)流換熱邊條 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 !如果為偶數(shù)層，向左掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i-1)+StepTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷隨時(shí)間分布：常數(shù) 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 !激活單元的上表面，指定為對(duì)流換熱邊條 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 !如果為奇數(shù)層，向右掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i+1)+StepTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷

15、隨時(shí)間分布：常數(shù) 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 !激活單元的表面，如果包含左端面，指定為對(duì)流換熱邊條 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 !激活單元的表面，如果包含右端面，指定為對(duì)流換熱邊條 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 !激活單元的側(cè)面，指定為對(duì)流換熱邊條 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 !激活單元的上表面，指定為對(duì)流換熱邊條 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、元的下表面，刪除對(duì)流換熱邊條 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 !如果為偶數(shù)層，向左掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i+1)+Ste

20、pTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷隨時(shí)間分布：常數(shù) 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 !激活單元的表面，如果包含左端面，指定為對(duì)流換熱邊條 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 !激活單元的表面，如果包含右端面，指定為對(duì)流換熱邊條 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,