我有一个关于三维梁单元Fortran方面的Vumat代码，想修改一下，然后代入Abaqus梁单元模型中，提交分析后，输出合理的荷载-位移曲线、应力应变曲线等，想让曲线基本达到输入CDP本构的曲线效果。