workbench中如何使用生死单元？ 50

根据这个视频：http://www.jishulink.com/college/video/c10300

在Workbench中要实现生死单元功能，需要插入APDL代码来实现，可我在操作中有一些问题。

【模型介绍】这里创建一个简单的模型做练习之用。创建一个6*6*1mm的长方体，中间分割出一个直径1mm的圆柱，在圆柱的上表面设置HeatFlux边界，1e7W/m^2。其他表面设为Convection边界，10W/m^2 K。

【插入APDL测试】在B5右键Insert->Commands，然后在右侧文本中插入一段APDL，这段程序很简单，就是设置了两个载荷步，第1秒分为10个子步，第2秒分为5个子步

/solu
autots,off

time,1
nsubst,10
solve

time,2
nsubst,5
solve

最后得到结果也符合预期，可以看到仿真了2秒，分别是10步和5步。

【生死单元测试】从上图可以看到1秒时最高温度达到了159度，因此打算将超过120度的单元杀死再进行后面的加载。插入的APDL改成下面这样

/solu
autots,off

time,1
NSUBST,10
solve
save
finish

/post1
    etable,temper1,temp         
    esel,s,etab,temper1,120    
    finish
/solu
ANTYPE,,REST
ekill,all
esel,s,live
time,2
NSUBST,5
solve

但实际计算中并没有按照设想的那样执行，只计算了第1秒的10步，如下图所示。弄了半天也不知道为什么会这样，还请大家指点啊！！