求助，SLM仿真生死单元问题？ 26

请问各位大佬，我在模拟两层铺粉的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,900,873,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,31.95,63.78,63.78,63.78 !* mptemp mptemp,1,0,500,1000,1500 mpdata,dens,1,1,8000,8000,8000 !* mptemp mptemp,1,20,1380,1381 mpdata,dens,2,1,5013,7833,8000 block,0,0.005,0,0.003,0,0.001 block,0.0019,0.0031,0.0009,0.0021,0.001,0.00103 block,0.0019,0.0031,0.0009,0.0021,0.00103,0.00106 vglue,all K,17,0,0.0009,0, K,18,0.005,0.0009,0, K,19,0.005,0.0009,0.001, K,20,0,0.0009,0.001, K,25,0,0.0021,0, K,26,0.005,0.0021,0, K,27,0.005,0.0021,0.001, K,28,0,0.0021,0.001, K,29,0.0019,0,0, K,30,0.0019,0.003,0, K,31,0.0019,0.003,0.001, K,32,0.0019,0,0.001, K,33,0.0031,0,0, K,34,0.0031,0.003,0, K,35,0.0031,0.003,0.001, K,36,0.0031,0,0.001, A,17,18,19,20 A,25,26,27,28 A,29,30,31,32 A,33,34,35,36 ASEL,S,,,15,16,1 ASEL,A,,,2 ASEL,A,,,13 VSBA,5,ALL ALLSEL,ALL vglue,all lesize,17,0.00002 lesize,20,0.00002 lesize,21,0.000015 lesize,30,0.00002 lesize,29,0.00002 lesize,37,0.000015 type,1 mat,2 vsweep,2 type,1 mat,2 vsweep,4 !***********竖线尺寸************** LESIZE,9, , ,15,2, , , ,1 LESIZE,12, , ,15,2, , , ,1 LESIZE,28, , ,15,2, , , ,1 LESIZE,36, , ,15,2, , , ,1 LESIZE,44, , ,15,2, , , ,1 LESIZE,48, , ,15,2, , , ,1 LESIZE,49, , ,15,2, , , ,1 LESIZE,50, , ,15,2, , , ,1 LESIZE,52, , ,15,2, , , ,1 LESIZE,54, , ,15,2, , , ,1 LESIZE,10, , ,15,0.5, , , ,1 LESIZE,11, , ,15,0.5, , , ,1 LESIZE,26, , ,15,0.5, , , ,1 LESIZE,34, , ,15,0.5, , , ,1 LESIZE,42, , ,15,0.5, , , ,1 LESIZE,46, , ,15,0.5, , , ,1 !*********横线尺寸**************** LESIZE,59, , ,20,0.1, , , ,1 LESIZE,60, , ,20,0.1, , , ,1 LESIZE,62, , ,20,0.1, , , ,1 LESIZE,67, , ,20,0.1, , , ,1 LESIZE,68, , ,20,0.1, , , ,1 LESIZE,70, , ,20,0.1, , , ,1 LESIZE,75, , ,20,0.1, , , ,1 LESIZE,76, , ,20,0.1, , , ,1 LESIZE,78, , ,20,0.1, , , ,1 LESIZE,83, , ,20,0.1, , , ,1 LESIZE,84, , ,20,0.1, , , ,1 LESIZE,86, , ,20,0.1, , , ,1 LESIZE,59, , ,20,0.1, , , ,1 LESIZE,61, , ,20,10, , , ,1 LESIZE,69, , ,20,10, , , ,1 LESIZE,77, , ,20,10, , , ,1 LESIZE,85, , ,20,10, , , ,1 !************中线尺寸************* LESIZE,65, , ,10,0.1, , , ,1 LESIZE,66, , ,10,0.1, , , ,1 LESIZE,72, , ,10,0.1, , , ,1 LESIZE,73, , ,10,0.1, , , ,1 LESIZE,74, , ,10,0.1, , , ,1 LESIZE,80, , ,10,0.1, , , ,1 LESIZE,88, , ,10,0.1, , , ,1 LESIZE,64, , ,10,0.1, , , ,1 LESIZE,82, , ,10,0.1, , , ,1 LESIZE,58, , ,10,0.1, , , ,1 LESIZE,81, , ,10,0.1, , , ,1 LESIZE,57, , ,10,0.1, , , ,1 LESIZE,55, , ,10,10, , , ,1 LESIZE,63, , ,10,10, , , ,1 LESIZE,71, , ,10,10, , , ,1 LESIZE,79, , ,10,10, , , ,1 LESIZE,87, , ,10,10, , , ,1 LESIZE,63, , ,10,10, , , ,1 type,1 mat,1 VSEL,S,,,6,12,1 VSEL,A,,,1 VSEL,A,,,3 vsweep,ALL ALLSEL,ALL !****对流载荷加载(除加工面外其他外表面都选上)******** asel,s,,,9,12,1 asel,a,,,26,28,1 asel,a,,,31,34,1 asel,a,,,37,40,1 asel,a,,,43,46,1 asel,a,,,40 asel,a,,,49,62,1 SFA,all,,CONV,480,22 ALLSEL,ALL !**************************** TUNIF,22 !预热温度 FINISH /SOL !****激光参数设置**** P=275 !功率 w0=0.05e-3 !光斑半径 v=0.8 !激光扫描速度 d=0.06e-3 !激光扫描间距 f=0.07 pi=3.1415926 a=0.05E-3 b=0.05E-3 c=5 !********************** ! ****瞬态求解参数设置**** lstep=0.08e-3 !扫描步长为0.1mm time_step=lstep/v !载荷步时间大小 ANTYPE,4 !定义新分析 TRNOPT,FULL NROPT,FULL TOFFST,273, !定义摄氏度的温度偏移量 LNSRCH,1 TIMINT,1 !打开时间积分选项进行瞬态分析 TINTP,0.005, , ,1, , , !************************ OUTRES,ALL,ALL my_time=0 XX=1.9e-3 !XX,YY用于定义局部坐标的原点 YY=1.32e-3 ZZ=1.03e-3 LOCAL,11,0,XX,YY,ZZ !**********热源方程设置********** *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNCCSYS *SET,_FNCNAME,'hflux' *DIM,_FNC_C1,,3 *DIM,_FNC_C2,,3 *DIM,_FNC_C3,,3 *DIM,_FNC_C4,,3 *SET,_FNC_C1(2),f *SET,_FNC_C2(2),p *SET,_FNC_C3(2),pi *SET,_FNC_C4(2),w0 *SET,_FNC_C1(3),f *SET,_FNC_C2(3),p *SET,_FNC_C3(3),pi *SET,_FNC_C4(3),w0 *SET,_FNCCSYS,11 ! /INPUT,hflux.func,,,1 *DIM,%_FNCNAME%,TABLE,6,31,3,,,,%_FNCCSYS% ! ! Begin of equation: {z} *SET,%_FNCNAME%(0,0,1), 0.001, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), 0.0 *SET,%_FNCNAME%(4,0,1), 0.0 *SET,%_FNCNAME%(5,0,1), 0.0 *SET,%_FNCNAME%(6,0,1), 0.0 *SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1, 4, 0, 0 *SET,%_FNCNAME%(0,2,1), 0 *SET,%_FNCNAME%(0,3,1), 0 *SET,%_FNCNAME%(0,4,1), 0 *SET,%_FNCNAME%(0,5,1), 0 *SET,%_FNCNAME%(0,6,1), 0 *SET,%_FNCNAME%(0,7,1), 0 *SET,%_FNCNAME%(0,8,1), 0 *SET,%_FNCNAME%(0,9,1), 0 *SET,%_FNCNAME%(0,10,1), 0 *SET,%_FNCNAME%(0,11,1), 0 *SET,%_FNCNAME%(0,12,1), 0 *SET,%_FNCNAME%(0,13,1), 0 *SET,%_FNCNAME%(0,14,1), 0 *SET,%_FNCNAME%(0,15,1), 0 *SET,%_FNCNAME%(0,16,1), 0 *SET,%_FNCNAME%(0,17,1), 0 *SET,%_FNCNAME%(0,18,1), 0 *SET,%_FNCNAME%(0,19,1), 0 *SET,%_FNCNAME%(0,20,1), 0 *SET,%_FNCNAME%(0,21,1), 0 *SET,%_FNCNAME%(0,22,1), 0 *SET,%_FNCNAME%(0,23,1), 0 *SET,%_FNCNAME%(0,24,1), 0 *SET,%_FNCNAME%(0,25,1), 0 *SET,%_FNCNAME%(0,26,1), 0 *SET,%_FNCNAME%(0,27,1), 0 *SET,%_FNCNAME%(0,28,1), 0 *SET,%_FNCNAME%(0,29,1), 0 *SET,%_FNCNAME%(0,30,1), 0 *SET,%_FNCNAME%(0,31,1), 0 ! End of equation: {z} ! ! Begin of equation: (f*p/(pi*w0^2))*exp(-({x}^2+{y}^2)/w0^2) *SET,%_FNCNAME%(0,0,2), 0.0010001, -999 *SET,%_FNCNAME%(2,0,2), 0.0 *SET,%_FNCNAME%(3,0,2), %_FNC_C1(2)% *SET,%_FNCNAME%(4,0,2), %_FNC_C2(2)% *SET,%_FNCNAME%(5,0,2), %_FNC_C3(2)% *SET,%_FNCNAME%(6,0,2), %_FNC_C4(2)% *SET,%_FNCNAME%(0,1,2), 1.0, -1, 0, 1, 17, 3, 18 *SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 2, 0, 0, 20 *SET,%_FNCNAME%(0,3,2), 0, -3, 0, 1, 20, 17, -2 *SET,%_FNCNAME%(0,4,2), 0.0, -2, 0, 1, 19, 3, -3 *SET,%_FNCNAME%(0,5,2), 0.0, -3, 0, 1, -1, 4, -2 *SET,%_FNCNAME%(0,6,2), 0.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,7,2), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,8,2), 0.0, -4, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,9,2), 0.0, -1, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,10,2), 0.0, -2, 0, 1, 2, 17, -1 *SET,%_FNCNAME%(0,11,2), 0.0, -1, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,12,2), 0.0, -5, 0, 1, 3, 17, -1 *SET,%_FNCNAME%(0,13,2), 0.0, -1, 0, 1, -2, 1, -5 *SET,%_FNCNAME%(0,14,2), 0.0, -2, 0, 1, -4, 3, -1 *SET,%_FNCNAME%(0,15,2), 0.0, -1, 0, 2, 0, 0, 20 *SET,%_FNCNAME%(0,16,2), 0.0, -4, 0, 1, 20, 17, -1 *SET,%_FNCNAME%(0,17,2), 0.0, -1, 0, 1, -2, 4, -4 *SET,%_FNCNAME%(0,18,2), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,19,2), 0.0, -2, 0, 1, -3, 3, -1 *SET,%_FNCNAME%(0,20,2), 0.0, 99, 0, 1, -2, 0, 0 *SET,%_FNCNAME%(0,21,2), 0 *SET,%_FNCNAME%(0,22,2), 0 *SET,%_FNCNAME%(0,23,2), 0 *SET,%_FNCNAME%(0,24,2), 0 *SET,%_FNCNAME%(0,25,2), 0 *SET,%_FNCNAME%(0,26,2), 0 *SET,%_FNCNAME%(0,27,2), 0 *SET,%_FNCNAME%(0,28,2), 0 *SET,%_FNCNAME%(0,29,2), 0 *SET,%_FNCNAME%(0,30,2), 0 *SET,%_FNCNAME%(0,31,2), 0 ! End of equation: (f*p/(pi*w0^2))*exp(-({x}^2+{y}^2)/w0^2) ! ! Begin of equation: (f*p/(0.00003*pi*w0^2))*exp(-({x}^2+{y}^2)/w0^2)*exp(-abs ! ({z})/0.00003) *SET,%_FNCNAME%(0,0,3), 0.00103, -999 *SET,%_FNCNAME%(2,0,3), 0.0 *SET,%_FNCNAME%(3,0,3), %_FNC_C1(3)% *SET,%_FNCNAME%(4,0,3), %_FNC_C2(3)% *SET,%_FNCNAME%(5,0,3), %_FNC_C3(3)% *SET,%_FNCNAME%(6,0,3), %_FNC_C4(3)% *SET,%_FNCNAME%(0,1,3), 1.0, -1, 0, 1, 17, 3, 18 *SET,%_FNCNAME%(0,2,3), 0.0, -2, 0, 0.00003, 0, 0, 19 *SET,%_FNCNAME%(0,3,3), 0, -3, 0, 1, -2, 3, 19 *SET,%_FNCNAME%(0,4,3), 0.0, -2, 0, 2, 0, 0, 20 *SET,%_FNCNAME%(0,5,3), 0.0, -4, 0, 1, 20, 17, -2 *SET,%_FNCNAME%(0,6,3), 0.0, -2, 0, 1, -3, 3, -4 *SET,%_FNCNAME%(0,7,3), 0.0, -3, 0, 1, -1, 4, -2 *SET,%_FNCNAME%(0,8,3), 0.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,9,3), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,10,3), 0.0, -4, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,11,3), 0.0, -1, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,12,3), 0.0, -2, 0, 1, 2, 17, -1 *SET,%_FNCNAME%(0,13,3), 0.0, -1, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,14,3), 0.0, -5, 0, 1, 3, 17, -1 *SET,%_FNCNAME%(0,15,3), 0.0, -1, 0, 1, -2, 1, -5 *SET,%_FNCNAME%(0,16,3), 0.0, -2, 0, 1, -4, 3, -1 *SET,%_FNCNAME%(0,17,3), 0.0, -1, 0, 2, 0, 0, 20 *SET,%_FNCNAME%(0,18,3), 0.0, -4, 0, 1, 20, 17, -1 *SET,%_FNCNAME%(0,19,3), 0.0, -1, 0, 1, -2, 4, -4 *SET,%_FNCNAME%(0,20,3), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,21,3), 0.0, -2, 0, 1, -3, 3, -1 *SET,%_FNCNAME%(0,22,3), 0.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,23,3), 0.0, -3, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,24,3), 0.0, -4, 0, 1, -1, 2, -3 *SET,%_FNCNAME%(0,25,3), 0.0, -1, 15, 1, 4, 0, 0 *SET,%_FNCNAME%(0,26,3), 0.0, -3, 0, 1, -4, 3, -1 *SET,%_FNCNAME%(0,27,3), 0.0, -1, 0, 0.00003, 0, 0, -3 *SET,%_FNCNAME%(0,28,3), 0.0, -4, 0, 1, -3, 4, -1 *SET,%_FNCNAME%(0,29,3), 0.0, -1, 7, 1, -4, 0, 0 *SET,%_FNCNAME%(0,30,3), 0.0, -3, 0, 1, -2, 3, -1 *SET,%_FNCNAME%(0,31,3), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation: (f*p/(0.00003*pi*w0^2))*exp(-({x}^2+{y}^2)/w0^2)*exp(-abs ! ({z})/0.00003) !--> !*******************第一层**************************** Vsel,s,,,4 Eslv,s Ekill,all allsel，all !*******************第一道****************************** *DO,i,0,2,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 !******************************** !****加载**** NSEL,S,LOC,X,0.0019-XX,0.0031-XX NSEL,R,LOC,Y,0.0009-YY,0.0021-YY NSEL,R,LOC,Z,-0.00004,0 ESLN,S SFE,ALL,1,HFLUX,,%hflux% ALLSEL,ALL ESEL,S,LIVE !************* !****求解，完成后删除载荷**** !save SOLVE SFEDELE,ALL,,HFLUX !****修改超过熔点的单元的材料模型**** ESEL,S,LIVE ESEL,R,MAT,,2 NSLE,S *GET,SELNNUM,NODE,,COUNT *get,SLNMIN,node,,num,min *dim,node_temp1,,selnnum *dim,SELNODE1,,SELNNUM SELNODE1(1)=SLNMIN *do,x,2,SELNNUM,1 SELNODE1(x)=ndnext(SELNODE1(x-1)) *enddo *do,f,1,selnnum,1 *get,node_temp1(f),node,selnode1(f),temp *enddo *do,f,1,selnnum,1 *IF,NODE_TEMP1(F),GE,1381,THEN NSEL,S,NODE,,SELNODE1(F) ESLN,S MPCHG,1,ALL *ENDIF *ENDDO ALLSEL,ALL !**************************** XX=XX+lstep *ENDDO