Implementing a simple 2D constitutive model for rocks into finite element method

Our research group has developed a 2D constitutive model for rocks that incorporates a strain-dependent elastic modulus. This model efficiently represents nonlinear stress–strain curves using just three equations and four parameters. In this study, we have implemented the simple 2D model into finite element codes to tackle engineering problems. We have calculated and compared stresses and displacements around a pressurized thick-walled hollow cylinder, a lined non-circular tunnel, and a rock slope using the elastic, simple 2D, and elasto-plastic models. The results confirmed that the simple 2D model closely matches the stress distribution obtained from the elastic solution at low-stress levels. Additionally, the simple 2D model produces a plastic zone with a larger inward displacement than the elasto-plastic model at higher stress levels. Only the simple 2D model captured sidewall convergence, roof sag, and floor heave in non-circular tunnels, and the toppling failure for 90° slope and toe translational slide for 60° slope in rock slopes with smooth critical slip surface. We did not encounter any convergence difficulties while solving the simple 2D model.