24小时热门版块排行榜

返回列表

kanrunzhe

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 956
散金: 59
帖子: 338
在线: 44.5小时
虫号: 5305243
注册: 2016-12-01
性别: GG
专业: 爆炸与冲击动力学

[求助] 求助切割命令流问题出在哪

/SOLU

temp_dele=3000 !杀死单元的温度

LSIZE=0.001

V=0.0025 !速度

TINC=LSIZE/V !单步时间长度

KK=1 !能量集中系数

Qmax=0.3e8

R=0.01

MAX_TIME=101

MAX_X=101

MAX_Y=101

*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME !定义一个表数组

*DO,K,1,MAX_TIME,1

*DO,I,1,MAX_X,1

FLUX2(I,0,K)=(I-1)*LSIZE

*ENDDO

*DO,J,1,MAX_Y,1

FLUX2(0,J,K)=(J-1)*LSIZE

*ENDDO

*DO,I,1,MAX_X,1

*DO,J,1,MAX_Y,1

XCENTER=V*(K-1)*TINC !热源X位置

YCENTER=0.05    !热源Y位置

DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)

*IF,DISTANCE,LE,R,THEN

FLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(3*r**2))

*ELSE

FLUX2(I,J,K)=0

*ENDIF

*ENDDO

*ENDDO

FLUX2(0,0,K)=(K-1)*TINC

*ENDDO

!开始计算

ANTYPE,TRANS

TUNIF,20 !初始温度20度

KBC,1

*SET,DELT,1

TIMINT,ON

AUTOTS,OFF

TT=5

*DO,TM,1E-6,TT,TINC
/solu

TIME,TM

AUTOTS,0

DELTIM, , , ,1

KBC,0

TSRES,ERASE

asel,s,loc,z,0

SFA,all,1,HFLUX,%FLUX2%    !用的表载荷

nsubst,40,80,20

SOLVE

allsel,all

OUTRES,ALL,ALL,

/post1
ETABLE,temp_over,TEMP,    !建立温度结果单元表
ESEL,S,ETAB,TEMP_OVE,temp_dele,, ,0 !选择温度超过temp_dele的单元
/solu
ekill,all    !杀死这些超温的单元
allsel,all

*ENDDO

OUTRES,ALL,ALL,

FINISH

计算出来只能杀死第一层（0.1mm）的单元，其他都没法杀死，请问问题出在哪？

回复此楼