应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ansys 生死单元模拟融化过程
61/1返回列表
查看: 1769  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
最怕的即是让爱自己的人失望
1楼2015-09-17 19:45:16
已阅   回复此楼   关注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，则可以显示结果
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

herozhu
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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭