复合材料干涉螺栓连接的压入力提取？

请问各位老师 目前在做复合材料干涉螺栓压入的仿真 想要提取螺栓头的压入荷载 用的是动力显式的求解器 并且使用自己编写的VUMAT定义复合材料本构来考虑干涉产生的损伤动能与内能的比值应该是满足准静态要求的 

但是提取的压入荷载很不理想 请问出现这种情况是因为质量缩放造成的 还是编写的子程序错误 或是提取压入荷载的方式错误