求助：热力耦合法进行热应力模拟时使用生死单元疑问

elflxg

我做金属凝固过程中的热应力模拟，在间接耦合中采用生死单元处理液相，
根据温度载荷作为"杀死"和激活的条件，看了坛子中的许多相关例子，但是还是有一个问题不是很明白：
       当被杀死的单元重新”激活“时，在计算热应变时 a*（T-T0），其中T0究竟是上一步载荷中相应的节点温度，还是起初设定的参考温度TREF呢？看参考资料中好像说是TREF，如果是这样的话，岂不是与实际情况不符合了，因为重新激活的单元在上一个载荷步中实际温度应该是低于TREF的。这样会不会造成刚激活的单元局部应力偏大呢？

JCHLI2000

不用生死单元处理液相，而通过设置液相的物性参数来处理这个问题不知是否可以？

elflxg

我以前也是这样做的，通过设置材料属性来处理液相，也成功的做过一些算例。但还是觉得，那样做有些问题：一是如果材料参数设置不合适，计算不容易收敛。二是觉得那样做不符合实际，想进一步提高模拟的精度。