UMAT子程序计算复合材料疲劳问题 相关求助？ 100

请问UMAT子程序中对于这种论文里所写的载荷应该如何施加？

我尝试对右两个端面设置参考点，之后对参考点施加concentrated force，数值为论文中的平均强度乘板的截面面积25mm×2.5mm再乘应力水平的百分数，但是在第一个循环里板的90度铺层就几乎完全失效，0度和45度的铺层没有出现失效单元，请问可能是什么原因造成的？

此外小白刚刚开始接触umat子程序，课题需要求复合材料层合板的刚度退化曲线，请问各位大神有没有推荐的课程或者教材？最好能包括结果的正确性验证，十分感谢！