论文部分内容阅读
在弹塑性模型数值实现的过程中,需要进行应力更新的回映算法。针对三维应力空间回映算法在奇异点收敛性方面的不足,提出主应力空间的回映算法,讨论了算法实现过程的应力空间转化问题,分析了应力更新过程中确定回映区域的方法,建立了相应的一致性刚度矩阵。基于大型有限元软件ABAQUS提供的用户材料子程序UMAT接口,利用Fortran编程语言,实现了Hoek-Brown准则主应力空间的隐式积分算法,利用开发的模型,进行了岩石常规三轴压缩试验的数值模拟,通过与ABAQUS内部的Mohr-Coulomb准则计算结果的对比,验证了模型和程序的可行性和准确性。
In the process of numerical implementation of elasto-plastic model, the mapping algorithm of stress update is needed. In order to solve the problem of the singularity convergence in the three-dimensional stress space mapping algorithm, the mapping method of the main stress space is proposed. The transformation of stress space in the process of algorithm implementation is discussed, and the method of determining the mapping area in stress updating is analyzed. The corresponding consistent stiffness matrix. Based on the user material subroutine UMAT interface provided by the large-scale finite element software ABAQUS, the implicit integration algorithm of the Hoek-Brown principle principal stress space is implemented by using the Fortran programming language. The numerical values of conventional rock triaxial compression tests The simulation and the comparison with the calculation results of Mohr-Coulomb criterion in ABAQUS validate the feasibility and accuracy of the model and program.