24小时热门版块排行榜    

查看: 1979  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

herozhu

木虫 (正式写手)

[求助] ansys 生死单元模拟融化过程 已有1人参与

刚接触ansys,需要用到生死单元,网上找到一个例子,阿尔伯特大学较早的一个例子,原址如下:http://www.mece.ualberta.ca/tuto ... ath/BirthDeath.html。按他的讲解一步一步做,最后老是不对,报错如下图。原文apdl代码如下,请高手指教如何能调试正确,多谢!


finish
/clear

/title, Convection Example       
/prep7                                ! Enter the preprocessor

! define geometry

k,1,0,0                                ! Define keypoints
k,2,0.03,0
k,3,0.03,0.03
k,4,0,0.03
a,1,2,3,4                        ! Connect the keypoints to form area


! mesh 2D areas

ET,1,plane55                        ! Element type

MP,Dens,1,920                        ! Define density
mp,c,1,2040                        ! Define specific heat
mp,kxx,1,1.8                        ! Define heat transfer coefficient

esize,0.0005                        ! Mesh size
amesh,all                        ! Mesh area

finish
/solu                                ! Enter solution phase

antype,4                        ! Transient analysis

time,60                                ! Time at end of analysis

nropt,full                        ! Newton Raphson - full
lumpm,0                                ! Lumped mass off
nsubst,20                        ! Number of substeps, 20
neqit,100                        ! Max no. of iterations
autots,off                        ! Auto time search off
lnsrch,on                        ! Line search on
outres,all,all                        ! Output data for all substeps
kbc,1                                ! Load applied in steps, not ramped

IC,all,temp,268                        ! Initial conditions, temp = 268

nsel,s,ext                        ! Node select all exterior nodes
sf,all,conv,10,368                ! Apply a convection BC
nsel,all                        ! Reselect all nodes
/gst,off                        ! Turn off graphical convergence monitor

solve
finish

/post1                                ! Enter postprocessor
set,last                        ! Read in last subset of data
etable,melty,temp,                ! Create an element table
esel,s,etab,melty,273                ! Select all elements from table above 273
finish

/solu                                ! Re-enter solution phase
antype,,rest                        ! Restart analysis
ekill,all                        ! Kill all selected elements
esel,all                        ! Re-select all elements

finish

/post1                                ! Re-enter postprocessor
set,last                        ! Read in last subset of data
esel,s,live                        ! Select all live elements
plnsol,temp                        ! Plot the temp contour of the live elements

ansys 生死单元模拟融化过程
最后一步报错.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

最怕的即是让爱自己的人失望
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

问题可能是出在如下命令:
       esel,s,etab,melty,273                ! Select all elements from table above 273
你可以先执行到这个命令,然后列表单元表 melty 和绘制单元温度云图,看看结果是否正确。
4楼2015-09-19 08:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
问题出在第一次求解之后,见如下命令流部分,有问题的地方以 !*** 作为标记:

/post1                                ! Enter postprocessor
set,last                        ! Read in last subset of data
etable,melty,temp,                ! Create an element table
esel,s,etab,melty,273                ! Select all elements from table above 273
finish                          !***  这一次进入 /post1 是要做什么?只是定义第一部吗?

/solu                                ! Re-enter solution phase
antype,,rest                        ! Restart analysis         !*** 此处表示执行重启动
ekill,all                        ! Kill all selected elements     !*** 此处杀死了所有的单元
esel,all                        ! Re-select all elements
         !***    这里没有求解:   solve,为什么?
finish

/post1                                ! Re-enter postprocessor
    /graphics, full        !***  针对你遇到的警告,此处需要关闭 Power Graphics
set,last                        ! Read in last subset of data
esel,s,live                        ! Select all live elements
           !***   经检查,不存在 live 单元 - 前面杀死了所有单元
plnsol,temp       ! Plot the temp contour of the live elements
           !***  没有 live 单元,无法显示结果;如果先执行 allsel,all,则可以显示结果

» 本帖已获得的红花(最新10朵)

2楼2015-09-18 09:31:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herozhu

木虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by shgao20 at 2015-09-18 09:31:30
问题出在第一次求解之后,见如下命令流部分,有问题的地方以 !*** 作为标记:

/post1                                ! Enter postprocessor
set,last                        ! Read in last subset of data
...

多谢指教。如你所述,最终没有结果显示是因为杀死步骤错误,杀死了所有单元,如何调整才能做到只杀死温度大于273的单元?etable 定义位置出错?多谢!
最怕的即是让爱自己的人失望
3楼2015-09-18 12:26:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
herozhu: 金币+30, ★★★★★最佳答案, 全给你了,多谢了 2015-09-22 09:58:26
上述回答有误。你可以将 solve $ fini 之后的命令流去掉,改为如下命令流试试:
  /post1                                ! Enter postprocessor
  set,last                               ! Read in last subset of data
  allsel,all
  etable,melty,temp,                ! Create an element table
  esel,s,etab,melty,273             ! Select all elements from table above 273

  cm, e_273, ELEM         !***  温度超过 273 度的单元,可以杀死。本句也可以不要。
  esel, inve                    !***  颠倒单元选择,即选择温度小于 273 的单元
  /graphics, full               !***  针对你遇到的警告,此处需要关闭 Power Graphics

  plnsol,temp       ! Plot the temp contour of the live elements
5楼2015-09-19 09:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +10 Tide man 2026-08-11 10/500 2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 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
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 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 医学老男孩
信息提示
请填处理意见