Axisymmetric adaptive lower bound finite element method based on Mohr-Coulomb yield criterion and second-order cone programming
SUN Rui,ZHANG Jian,YANG Junsheng,YANG Feng
DOI: https://doi.org/10.11779/cjge20220781
2023-01-01
Abstract:The axisymmetric Mohr Coulomb yield surfaces have edges and corners in the three-dimensional stress space, which leads to difficulties in numerical calculation. Therefore, how to deal with axisymmetric Mohr Coulomb criteria efficiently has always been an important research content of the lower bound finite element limit analysis（LB-FELA） method. Firstly, the axisymmetric Mohr Coulomb criterion is transformed into a set of inequality constraints and three linear equality constraints by introducing the complete plasticity assumption; Then, the inequality constraint is directly transformed into a second-order cone constraint, which avoids the smooth approximation of numerical singularities; Finally, the axisymmetric LB-FELA model based on Mohr Coulomb criterion is transformed into a mathematical optimization calculation model of second-order cone programming. The linear stress element adopted by the axisymmetric LB-FELA cannot accurately simulate the stress change in the failure region. Therefore, the mesh distribution form has a great impact on the calculation accuracy of the LB-FELA. To solve this problem, an adaptive mesh refinement strategy based on element yield residual is proposed. By judging the degree of the node stress in the element close to the yield, the elements can be automatically identified and refined in the failure area. By analyzing the typical axisymmetric geotechnical engineering stability problems such as the bearing capacity of circular foundation and the ultimate uplift capacity of vertical anchor plate, it shows that the proposed method has high calculation efficiency and accuracy, and has certain theoretical value and application prospect.