vumat子程序事计算发散问题

根据2020Jiang的聚合物本构文章写了应用于显式分析计算的vumat子程序。单胞计算正常，但在计算多单元模型时，在前面弹性阶段就有严重的发散问题，导致计算中途崩溃。希望能够进行解答。