关于ansys模拟增材制造的问题？ 20

请教大佬们，我是做SLM仿真的，在模拟过程中我发现我得到的温度场始终是这样的圆形，而不是文献中这种类似慧星状的结果，是不是因为在循环加载过程中前一步计算的温度场结果没有被后一步继承，应该怎么解决，附上命令流求大佬们指导一下 finish /clear /prep7 /unit,si et,1,solid70 !****材料属性设置**** et,1,solid70 mptemp,1,20,200,400,600,800,1000 mptemp,7,1200,1400,1500,1600,1700 mpdata,c,1,1,470,508,550,592,634,676 mpdata,c,1,7,719,765,900,873,873 mpdata,kxx,1,1,13.31,16.33,19.47,22.38,25.07,27.53 mpdata,kxx,1,7,29.76,31.95,63.78,63.78,63.78 mpdata,c,2,1,470,508,550,592,634,676 mpdata,c,2,7,719,765,900,873,873 mpdata,kxx,2,1,0.15,0.23,0.32,0.41,0.50,0.59 mpdata,kxx,2,7,0.68,0.77,0.85,0.91,1.04 !* mptemp mptemp,1,0,500,1000,1500 mpdata,enth,1,1,1.25E+5,3.375E+5,6E+5,8.375E+5 mpdata,enth,2,1,1.25E+5,3.375E+5,6E+5,8.375E+5 mpdata,dens,1,1,8000,8000,8000 !* mptemp mptemp,1,20,1380,1381 mpdata,dens,2,1,5013,7833,8000 block,0,0.0012,0,0.0008,0,0.0006 block,0,0.0012,0,0.0008,0.0006,0.00063 vglue,all lesize,17,0.00001 lesize,18,0.00001 lesize,26,0.000015 type,1 mat,2 vsweep,3 !****粉床网格划分**** type,1 mat,1 Extopt,esize,15,0.2 Vsweep,1,1,2 !****基板网格划分*** !************************* ALLSEL,ALL !****对流载荷加载(除加工面外其他外表面都选上)******** SFA,5,,CONV,480,22 SFA,3,,CONV,480,22 SFA,6,,CONV,480,22 SFA,4,,CONV,480,22 SFA,16,,CONV,480,22 SFA,15,,CONV,480,22 SFA,14,,CONV,480,22 SFA,13,,CONV,480,22 SFA,1,,CONV,480,22 ALLSEL,ALL /REP,FAST !**************************** TUNIF,22 !预热温度 FINISH /SOL !****激光参数设置**** P=50 !激光功率 v=1 !激光扫描速度 d=0.12e-3 !激光扫描间距 f=0.35e-7 pi=3.1415926 a=0.00009 b=0.00009 c=0.00002 !********************** ! ****瞬态求解参数设置**** lstep=0.04e-3 !扫描步长 time_step=lstep/v !载荷步时间大小 ANTYPE,4 !定义新分析 TRNOPT,FULL TOFFST,273, !定义摄氏度的温度偏移量 !* TIMINT,1 !打开时间积分选项进行瞬态分析 !************************ OUTRES,ALL,ALL my_time=0 XX=0.2e-3 !XX,YY用于定义局部坐标的原点 YY=0.2e-3 ZZ=0.63e-3 LOCAL,11,0,XX,YY,ZZ !***************定义热源函数***** *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNC_C6 *DEL,_FNCCSYS *SET,_FNCNAME,'hflux' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *DIM,_FNC_C6,,1 *SET,_FNC_C1(1),p *SET,_FNC_C2(1),f *SET,_FNC_C3(1),pi *SET,_FNC_C4(1),a *SET,_FNC_C5(1),b *SET,_FNC_C6(1),c *SET,_FNCCSYS,11 ! /INPUT,hflux.func,,,1 *DIM,%_FNCNAME%,TABLE,8,44,1,,,,%_FNCCSYS% ! ! Begin of equation: ((6*(3^0.5)*p*f)/(pi*(pi^0.5)*a*b*c))*exp((-3* ! ({X}^2)/(a^2))-(3*({Y}^2)/(b^2))-(3*({Z}^2)/(c^2))) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)% *SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)% *SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)% *SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)% *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 3, 0, 0, 0 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 0.5, 0, 0, -1 *SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 6, 0, 0, -3 *SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, -3 *SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 1, -2, 3, 17 *SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, -1, 3, 18 *SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 0.5, 0, 0, 19 *SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 19, 17, -1 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, 19, 3, -3 *SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, -1, 3, 20 *SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -3, 3, 21 *SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, -1, 3, 22 *SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 1, -2, 4, -3 *SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -4, 0, 1, -2, 2, -3 *SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 3, 0, 0, -4 *SET,%_FNCNAME%(0,19,1), 0.0, -3, 0, 1, -4, 3, -2 *SET,%_FNCNAME%(0,20,1), 0.0, -2, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,21,1), 0.0, -4, 0, 1, 2, 17, -2 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, -3, 3, -4 *SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 2, 0, 0, 20 *SET,%_FNCNAME%(0,24,1), 0.0, -4, 0, 1, 20, 17, -3 *SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -2, 4, -4 *SET,%_FNCNAME%(0,26,1), 0.0, -2, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,27,1), 0.0, -4, 0, 1, 3, 17, -2 *SET,%_FNCNAME%(0,28,1), 0.0, -2, 0, 3, 0, 0, -4 *SET,%_FNCNAME%(0,29,1), 0.0, -5, 0, 1, -2, 3, -4 *SET,%_FNCNAME%(0,30,1), 0.0, -2, 0, 2, 0, 0, 21 *SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, 21, 17, -2 *SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -5, 4, -4 *SET,%_FNCNAME%(0,33,1), 0.0, -4, 0, 1, -3, 2, -2 *SET,%_FNCNAME%(0,34,1), 0.0, -2, 0, 2, 0, 0, 4 *SET,%_FNCNAME%(0,35,1), 0.0, -3, 0, 1, 4, 17, -2 *SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 3, 0, 0, -3 *SET,%_FNCNAME%(0,37,1), 0.0, -5, 0, 1, -2, 3, -3 *SET,%_FNCNAME%(0,38,1), 0.0, -2, 0, 2, 0, 0, 22 *SET,%_FNCNAME%(0,39,1), 0.0, -3, 0, 1, 22, 17, -2 *SET,%_FNCNAME%(0,40,1), 0.0, -2, 0, 1, -5, 4, -3 *SET,%_FNCNAME%(0,41,1), 0.0, -3, 0, 1, -4, 2, -2 *SET,%_FNCNAME%(0,42,1), 0.0, -2, 7, 1, -3, 0, 0 *SET,%_FNCNAME%(0,43,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,44,1), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation: ((6*(3^0.5)*p*f)/(pi*(pi^0.5)*a*b*c))*exp((-3*({X}^2)/(a^2)) ! -(3*({Y}^2)/(b^2))-(3*({Z}^2)/(c^2))) !--> !************************************************* *DO,i,0,5,1 LOCAL,11,0,XX,YY,ZZ !****每个载荷步的求解参数设置**** my_time=my_time+time_step TIME,my_time AUTOTS,1 DELTIM,lstep/v,lstep/v,lstep/v, KBC,1 !******************************** !****加载**** SFA,8,,HFLUX, %HFLUX% !************* !****求解，完成后删除载荷**** ALLSEL,ALL !save SOLVE SFADELE,8,,HFLUX !**************************** XX=XX+lstep *ENDDO !****显示结果**** FINISH /POST1 !* /EFACET,1 PLNSOL, TEMP,, 0 !****************