24小时热门版块排行榜

返回列表

huangfumin

木虫 (著名写手)

应助: 0 (幼儿园)
金币: 4219.2
散金: 5
红花: 7
帖子: 1038
在线: 306.1小时
虫号: 909183
注册: 2009-11-21
专业: 制造系统与自动化

[求助] ansys 生死技术 & 继承性循环 & 温度漂移 & APDL

各位ansys同仁：

  小弟在用ansys做简单生死单元技术时，遇到了一些问题，还望各位帮忙解答一下。
模型的APDL程序如下：

finish
/clear

/prep7
k,1,0,0
k,2,0.03,0
k,3,0.03,0.03
k,4,0,0.03
a,1,2,3,4

nt=20
tall=60
tsize=tall/nt

ET,1,Plane13,2
MP,Dens,1,920
mp,c,1,2040
mp,kxx,1,1.8

esize,0.0005
amesh,all
finish

/solu
tunif,300
nsel,s,ext
sf,all,conv,10,2000

antype,4
nropt,full
lumpm,0
neqit,100
autots,on
outres,all,all

time,10*tsize
nsel,all
solve
finish

*do,i,2,4
/solu
   resume
antype,,rest
time,i*tsize*10
nsel,s,ext
sfdele,all,conv

nropt,full
allsel
solve
save
finish

/post1
set,last
etable,melty,temp
esel,s,etab,melty,340
finish

/solu
ekill,all
esel,s,live
nsle,s
nsel,inve
D,all,all,0
finish
*enddo

不过程序的结果有问题，这个程序参考的是University of Alberta一个简单程序（如上或者附件所示），几何边界条件为对流。我单独用了一个tsize步长求解，为了循环语句中的resume和antype,,rest。同时为了研究生死单元技术，我在循环中就将载荷删除施加的载荷，为的是检测是否生死单元技术能正确解决问题。
现在主要问题如下：
1. 得到的时间段有问题，理论上求解时间应该为30(最开始的单独语句-非*do循环内)、60、90和120[s]，但实际得到的是30、31、32和33[s]。

2. 在循环语句中，我已经将载荷删除，且对需要杀死单元的节点进行约束，按照理论上说被杀死的单元或者节点在后续求解时间内就不应该有温度变化才对啊？（是不是我对有限元或者生死单元技术理解有问题？）。部分图片如附件所示。
图片：31秒的所有单元和仅仅为活单元的温度场

32秒的所有单元和仅仅为活单元的温度场

33秒的所有单元和仅仅为活单元的温度场

不知道有人遇到相同问题没有？又是如何解决的呢？谢谢您的指点。

回复此楼