能否按照混凝土CDP本构或者其他本构,写一个Vumat梁单元子程序,代入如下的三维梁单元模型中,分析成功后能导出CDP本构一样的荷载-位移曲线、应力-应变曲线、弯矩-转角曲线等,最重要的是,在写代码的过程中能把本构中相关的理论、公式等进行改进,有所创新,并像我发的参考文献里那样把创新的理论、公式写一下,因为光有代码肯定不行,要把代码涉及到的理论、公式等改进创新后展现出来,毕竟之后自己要写在论文中。(改进理论、公式肯定不能和参考文献里一样)