ANSYS中被杀死单元为何还能传递热量? 100

下图是一个简单示例,热源只加载在中间的一个小区域上,但当热源周围的单元被杀死后,周围部分的温度依然在升高,并且被杀死,也就是被杀死的单元传递了热量。但这是我不想看到的,也是与实际不符合的。

请问这是为什么?我该怎么做来避免出现这种情况呢?

下图是用APDL实现的,相关程序附后。

1.gif

dt1=0.001
T=0
 
allsel
T=T+dt1
time,T
NSUBST,1 
solve                       !设置的时间先计算一下
save
finish
 
*do,i,1,10
    /post1
    etable,temper1,temp         !将温度存储在temperature表中
    esel,s,etab,temper1,1400    !选择所有温度大于等于1400的单元
    finish
    /solu
    *get,elemnum,elem,0,count
    *if,elemnum,gt,0,then
       ekill,all
       esel,s,live
    *else
       allsel
    *endif
    eplot
    resume                  
    ANTYPE,,REST            
    T=T+dt1
    time,T
    nropt,full
    NSUBST,1 
    solve
    save
    finish
*enddo


************************************************************************************************************

可直接运行的程序:

ANSYS中被杀死单元为何还能传递热量?的图2二维轴对称模型激光打孔命令流.txt


应答将前几步的结果整理放在下面:


1.模型厚度1mm,试件半径3mm,热源半径50um,下图红色箭头所示即为热源,为固定热源

1.png

2.进行一步求解后(0.1s),未杀死单元时,得到温度分布,最高已经达到5万多度

2.png

3.选中大于1400度的单元进行杀死,剩下的活单元如图所示

3.png

4.此时的温度场如图所示,此时热源已经完全施加在被杀死单元上,后面应该没有影响才对

4.png

5.可进行下一步,结束时间0.2s的计算后,周围温度依然在升高,一直计算下去如上动画所示。

5.png


邀请回答 举报

共5个回答

默认排序
  • 孙冬
    2018年3月2日
    使用allsel,all选择了所有的单元,包括被杀死的单元,发现被杀死的单元温度也是很高的,上万摄氏度,已经被杀死的单元为什么还有这么高的温度?
    0评论 采纳答案 举报

  • frishing
    2018年3月21日
    杀死之前要用D命令限定温度,不然被杀死的单元自由度可能出现漂移
    1评论 采纳答案 举报

  • ANSYS中被杀死单元为何还能传递热量?的相关案例教程

    Spring comes 小伙们大家好!之前有一期讲了ABAQUS进行反应谱分析时怎么进行静动叠加,利用工况叠加的原理。这一期主要讲解了ANSYS反应谱分析时怎么进行静动叠加。其实,同样是采用了工况叠加。 我接下来以一根柱子来做这个案例。反应谱采用水工抗震规范[1]的反应谱。柱子的尺寸是1×2×5m,弹性模量为2.1E8Pa,泊松比0.2,密度2400kg/m3。 假设本案例地面最大加速度为 a=
    为什么要导出单元刚度矩阵 在学习有限元方法时,我们会需要编写程序计算结构的单元刚度矩阵。此外,当我们需要做有限元软件二次开发时,我们也需要验证所做的开发是否正确。为了验证程序正确性,我们可以从商业有限元软件中导出单元刚度矩阵来验证程序的计算结果。下面简单介绍从ansys软件中导出平面四边形四节点单元的单元刚度矩阵。 平面四边形四节点单元示例 如图所示,计算这两个单元组成单元刚度矩阵,并组装成整体刚
    点击蓝字 关注我们 好久没有更新算例了,这期更新一个小算例,用强度折减法计算边坡的稳定性,用ANSYS软件计算,ABAQUS的算例我前面写了两篇,感兴趣的可以看看。 边坡的尺寸及参数与之前的ABAQUS的边坡算例相同。模型如下。 接下来讲解一下建模的过程。先定义一下单元的类型。这里采用平面应变单元,PLANE42单元。 !定义单元类型ET,1,PLANE82     !定义边坡围岩单元KEYOPT
    作者:水哥ANSYS 来源:本文源于ANSYS结构院,上海安世亚太授权转载 随机分布在材料微观力学分析中扮演着重要角色,例如混凝土骨料力学、新型材料纤维力学分析等内容,提及随机分布,更多的同学可能会联想到采用第三方软件如Matlab来生成,并导入ANSYS计算,其实ANSYS本身自带随机分布功能,只是功能略有限制。 ANSYS中产生随机分布的一个重要函数是 *VFILL,该函数主要的作用是对数组进
    该命令流为计算单转子-支承系统在加速运行过程中,受质量不平衡激励下的瞬态动力学响应。可以准确计算出在共振转速下的峰值及弯曲应变能情况。给出了详细的表加载转速和不平衡力的方法,可供参考。 /prep7 ! ** parameters length = 0.4 ro_shaft = 0.01 ro_disk = 0.15 md = 16.47 id = 9.427e-2 ip = 0.1861 kxx