我寫了一個混凝土水化熱過程中彈性模量隨時間變化的UMAT子程序，程序可以被調用，但結果卻與實際情況不符，正常情況下，混凝土放熱過程中應該是外面受拉，內部受壓，但子程序的結果卻是剛開始符合實際，后面應力卻發生突變，變成內部受拉，外部受壓，我試了下，如果不調用子程序，而是設定一個不變的彈性模量，應力的變化趨勢與實際情況相符，請問這種情況是什么原因呢？

這是彈性模量隨時間變化情況

AA=0.28
BB=0.52

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

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