查看: 1102  |  回复: 1

kanrunzhe

银虫 (正式写手)

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

/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)的单元,其他都没法杀死,请问问题出在哪?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kanrunzhe

银虫 (正式写手)

2楼2017-07-24 15:18:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kanrunzhe 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +3 Tide man 2026-08-10 3/150 2026-08-10 21:48 by GZUZM
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见