請問UMAT子程序中對于這種論文里所寫的載荷應該如何施加？

我嘗試對右兩個端面設置參考點，之后對參考點施加concentrated force，數值為論文中的平均強度乘板的截面面積25mm×2.5mm再乘應力水平的百分數，但是在第一個循環里板的90度鋪層就幾乎完全失效，0度和45度的鋪層沒有出現失效單元，請問可能是什么原因造成的？

此外小白剛剛開始接觸umat子程序，課題需要求復合材料層合板的剛度退化曲線，請問各位大神有沒有推薦的課程或者教材？最好能包括結果的正確性驗證，十分感謝！