Abstract
Mesh-free methods, such as the Reproducing Kernel Particle Method, have been proven advantageous in modeling excessive deformation problems involving material separation, fracture, impact, etc. However, the domain integration in the RKPM remains challenging due to instability and sub-optimal convergence for high strain rate events. Although some novel developments alleviate the above issue, they are either computationally expensive or require evaluating the contour integral, which is not straightforward to obtain in contact and material separation problems using meshfree discretization. This work develops a simple and stable integration method based on the extension of modified Simpson’s rule. The method is free from conforming subdomains and can straightforwardly be applied to the meshfree formulation with updated configuration. To model penetration into the earth, a standard viscous boundary is introduced to address the issue of reflecting waves from the truncated computational domain for the ground target. proposed semi-Lagrangian RK framework to mimic the earth domain. The numerical results are validated with experimental data for various geo-materials and experimental setups.