混凝土弹性模量UMAT结果有误？ 20

我写了一个混凝土水化热过程中弹性模量随时间变化的UMAT子程序，程序可以被调用，但结果却与实际情况不符，正常情况下，混凝土放热过程中应该是外面受拉，内部受压，但子程序的结果却是刚开始符合实际，后面应力却发生突变，变成内部受拉，外部受压，我试了下，如果不调用子程序，而是设定一个不变的弹性模量，应力的变化趋势与实际情况相符，请问这种情况是什么原因呢？

这是弹性模量随时间变化情况

AA=0.28
BB=0.52

AAB=-AA*(TIME(2)**BB)

EMOD=PROPS(1)*(1.0D0-DEXP(AAB))