关于结冰问题，希望大家帮忙看看

TI99

首先，目前我想仿真的是水结冰对铜管壁造成的应力与应变。
我将水结冰产生的体积变化用热应力分析中的热膨胀系数来代替，用热失配造成的相互作用力来仿真结冰造成的膨胀效果。
仿真结果铜管的应力已经超过了弹性限，已经发生了塑性变形。
下面我想仿真经过10次的温度循环之后，经过10次的结冰与融化过程后，铜管的变形有多大
用多载荷步进行加载的时候，第一步载荷将温度加到1度（因为热膨胀系数已经设成上升1度，水的体积变成冰的体积）
第二步将温度delete掉
第三步再将温度加到1度
计算载荷步1到载荷步3
结果是载荷步2的时候塑性形变比载荷1的塑性形变小。
请问我这样加载有问题吗？
我将载荷去掉了，不就相当于冰化了吗？为什么塑性形变会变小呢？
大家能不能给点建议，因为这个问题已经搞了很久了，谢谢！