Equivalent Nodal Intensity-Based Contact Model for 3D Finite–Discrete Element Method in Rock Fracturing Analysis
摘要
In this paper, a novel equivalent nodal intensity-based contact model is proposed to enhance 3D FDEM (Combined Finite–Discrete-Element Method) in fracturing analysis. The normal contact force is uniformly distributed on the overlapped region between two contact bodies, with its magnitude being linearly proportional to the equivalent nodal intensity, which can be readily obtained based on the geometric characteristics of the overlapped region. In comparison with the existing contact models in 3D FDEM, the proposed contact model achieves significant simplification and enables efficient parallelization on General Purpose Graphic Processing Unit (GPGPU). By conducting a range of comparative tests, it has been confirmed that the proposed contact model correctly captures the interaction between discrete bodies and produces results consistent with those generated by the original contact model and laboratory experiment, thus validating the feasibility and applicability of the proposed contact model. Meanwhile, the efficiency analysis indicates that, for GPGPU parallel version, the proposed contact model achieves a speed-up ratio of 7.2 compared to the original model in 3D FDEM. The proposed contact model is physically sound and computationally efficient, which holds promise for driving advancements in 3D FDEM.