24小时热门版块排行榜    

查看: 2024  |  回复: 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

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +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

专家顾问 (职业作家)

【答案】应助回帖

问题可能是出在如下命令:
       esel,s,etab,melty,273                ! Select all elements from table above 273
你可以先执行到这个命令,然后列表单元表 melty 和绘制单元温度云图,看看结果是否正确。
4楼2015-09-19 08:56:40
已阅   回复此楼   关注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的回帖

你要哪样子

金虫 (初入文坛)

楼主,哪里还可以找到这种有GUI操作的生死单元的例子,急求
6楼2017-01-04 10:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 herozhu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 7/350 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +4 SCSIOyxguo 2026-09-03 7/350 2026-09-05 10:39 by Honglianlian
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见