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

摘要

Purpose Continua 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/approach Several 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). Findings Three 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/value This 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.