Y-Mat: an improved hybrid finite-discrete element code for addressing geotechnical and geological engineering problems

Abstract

PurposeContinua and discontinua coexist in natural rock materials. This paper aims to present an improved approach for addressing the mechanical response of rock masses based on the combined finite-discrete element method (FDEM) proposed by Munjiza.Design/methodology/approachSeveral algorithms have been programmed in the new approach. The algorithms include (1) a simpler and more efficient algorithm to calculate the contact force; (2) An algorithm for tangential contact force closer to the actual physical process; (3) a plastic yielding criterion (e.g. Mohr-Coulomb) to modify the elastic stress for fitting the mechanical behavior of elastoplastic materials; and (4) a complete code for the mechanical calculation to be implemented in Matrix Laboratory (MATLAB).FindingsThree case studies, including two standard laboratory experiments (uniaxial compression and Brazilian split test) and one engineering-scale anti-dip slop model, are presented to illustrate the feasibility of the Y-Mat code and its ability to deal with multi-scale rock mechanics problems. The results, including the progressive failure process, failure mode and trajectory of each case, are acceptable compared to other corresponding studies. It is shown that, the code is capable of modeling geotechnical and geological engineering problems.Originality/valueThis article gives an improved FDEM-based numerical calculation code. And, feasibility of the code is verified through three cases. It can effectively solve the geotechnical and geological engineering problems.