复合材料固化历程仿真 heTVal 子程序编写问题？

题主对一篇复合材料固化过程进行仿真复现， 完全找论文中给定的材料参数及固化动力学设定进行编写了子程序heTVal，可是计算所得的固化度历程完全与作者论文中给出的结果完全不同，纠结了很久，总是计算不出正确的结果，请问各位前辈可能是因为哪方面的问题呢？  我的错误计算方法所得的固化度会较早以很快的固化速率达到固化度为1的状态，对比如下图，图1为我计算所得的  图2为作者论文给出的