Multiscale semi-lagrangian reproducing kernel particle method for modeling damage evolution in geomaterials

Damage processes in geomaterials typically involve moving strong and weak discontinuities, multiscale phenomena, excessive deformation, and multibody contact that cannot be effectively modeled by a single-scale Lagrangian finite element formulation. In this work, we introduce a semi-Lagrangian Reproducing Kernel Particle Method (RKPM) which allows flexible adjustment of locality, continuity, polynomial reproducibility, and h- and p-adaptivity as the computational framework for modeling complex damage processes in geomaterials. Under this work, we consider damage in the continua as the homogenization of micro-cracks in the microstructures. Bridging between the cracked microstructure and the damaged continuum is facilitated by the equivalence of Helmholtz free energy between the two scales. As such, damage in the continua, represented by the degradation of continua, can be characterized from the Helmholtz free energy. Under this framework, a unified approach for numerical characterization of a class of damage evolution functions has been proposed. An implicit gradient operator is embedded in the reproduction kernel approximation as a regularization of ill-posedness in strain localization. Demonstration problems include numerical simulation of fragment-impact of concrete materials.