Smoothed particle hydrodynamics (SPH) is a meshless method gaining popularity recently in geotechnical modeling. It is suitable to solve problems involving large deformation, free-surface, cracking and fragmentation. To promote the research and application of SPH in geotechnical engineering, we present LOQUAT, an open-source three-dimensional GPU accelerated SPH solver. LOQUAT employs the standard SPH formulations for solids with two geomechnical constitutive models which are the Drucker–Prager model and a hypoplastic model. Three stabilization techniques, namely, artificial viscosity, artificial pressure and stress regularization are included. A generalized boundary particle method is presented to model static and moving boundaries with arbitrary geometry. LOQUAT employs GPU acceleration technique to greatly increase the computational efficiency. Numerical examples show that the solver is convergent, stable and highly efficient. With a mainstream GPU, it can simulate large scale problems with tens of millions of particles, and easily performs more than one thousand times faster than serial CPU code
