FRIC子程序结果不对，求帮助

用fortran编写了一个FRIC子程序，想要描述的是随着接触面间相对位移的增大，剪应力逐渐增大。可是放到abaqus中分析，发现增大的比例比理论的小很多（直线的斜率比理论的小），请知道的老师同学指点一下，谢谢。