Resultant Normal Contact Force-Based Contact Friction Model for the Combined Finite-Discrete Element Method and Its Validation

In the conventional FDEM (Combined Finite and Discrete Element Method), each contact pair might have multiple contact points where friction forces are applied, leading to non-unique friction force assignments and potentially introducing computational errors. This study introduces a new contact friction algorithm for FDEM based on the resultant normal contact force. This method necessitates determining the friction force at a unique equivalent contact point, thereby significantly simplifying the computational flow and reducing memory usage. A series of numerical tests are performed to validate the effectiveness of the proposed contact model. Using collision and block sliding tests, the proposed contact friction model is verified to be able to accurately capture the frictional effect between discrete bodies and circumvent the problematic kinetic energy dissipation issue associated with the original contact friction algorithm. For the Brazilian splitting and uniaxial compression tests, the simulated results closely align with those generated using the original contact friction algorithm and match the experimental measurements well, demonstrating the applicability of the proposed algorithm in fracturing analysis. Furthermore, by using the proposed contact friction algorithm, a computational efficiency enhancement of 8% in contact force evaluation can be achieved.