我要实现的是一个可以考虑高温情况下的混凝土弹塑性损伤本构模型(本构原理通过gpt解读或者查看abaqusCDP本构手册,就是在abaqus中CDP本构基础上升级的一个本构),目前的本构模型已经实现了一个单元在单轴拉伸、单轴压缩、纯剪切或者是循环加卸载情况下的正常计算,以及升温时相应的温度损伤变化等工况。总而言之,一个单元情况下计算结果是理想且正确的。
但是当将该本构模型应用复杂模型中进行计算时(说是复杂模型,其实也就是单元数量较多的模型,例如常温情况下的简支梁模型、高温情况下的盾构管片试验模型等),出现了计算无法收敛的情况。而根据这些工况的计算结果来看,其实只是受力较大的一些区域先出现了损伤(可能损伤单元仅占单元总数量的1%),然后计算增量步就不断变小、随后就无法收敛报错了。重要的是,我的本构基本上是在abaqus中自带的CDP本构上升级得到,相同的计算工况CDP本构可以完整计算。最终的结果就是,我的本构计算停止时仅有局部区域发生了损伤(甚至损伤最大值都没有达到损伤设置的上限0.9),而CDP本构的计算结果则是结构发生了大面的损伤。