A Total Lagrangian upwind Smooth Particle Hydrodynamics algorithm for large strain explicit solid dynamics

Lee, C. H. , Gil, A. J., Ghavamian, A. and Bonet, J. (2019) A Total Lagrangian upwind Smooth Particle Hydrodynamics algorithm for large strain explicit solid dynamics. Computer Methods in Applied Mechanics and Engineering, 344, pp. 209-250. (doi: 10.1016/j.cma.2018.09.033)

169457.pdf - Accepted Version
Available under License Creative Commons Attribution Non-commercial No Derivatives.



In previous work, Lee et al. introduced a new Smooth Particle Hydrodynamics (SPH) computational framework for large strain explicit solid dynamics with special emphasis on the treatment of near incompressibility. A first order system of hyperbolic equations was presented expressed in terms of the linear momentum and the minors of the deformation, namely the deformation gradient, its co-factor and its Jacobian. Taking advantage of this representation, the suppression of numerical deficiencies (e.g. spurious pressure, long term instability and/or consistency issues) was addressed through well-established stabilisation procedures. In Reference [1], the adaptation of the very efficient Jameson-Schmidt-Turkel algorithm was presented. Reference [2] introduced an adapted variationally consistent Streamline Upwind Petrov Galerkin methodology. In this paper, we now introduce a third alternative stabilisation strategy, extremely competitive, and which does not require the selection of any user-defined artificial stabilisation parameter. Specifically, a characteristic-based Riemann solver in conjunction with a linear reconstruction procedure is used, with the aim to guarantee both consistency and conservation of the overall algorithm. We show that the proposed SPH formulation is very similar in nature to that of the upwind vertex centred Finite Volume Method presented in [3]. In order to extend the application range towards the incompressibility limit, an artificial compressibility algorithm is also developed. Finally, an extensive set of challenging numerical examples is analysed. The new SPH algorithm shows excellent behaviour in compressible, nearly incompressible and truly incompressible scenarios, yielding second order of convergence for velocities, deviatoric and volumetric components of the stress.

Item Type:Articles
Glasgow Author(s) Enlighten ID:Lee, Dr Chun Hean
Authors: Lee, C. H., Gil, A. J., Ghavamian, A., and Bonet, J.
College/School:College of Science and Engineering > School of Engineering > Infrastructure and Environment
Journal Name:Computer Methods in Applied Mechanics and Engineering
Journal Abbr.:CMAME
ISSN (Online):1879-2138
Published Online:13 October 2018
Copyright Holders:Copyright © 2018 Elsevier B.V.
First Published:First published in Computer Methods in Applied Mechanics and Engineering 344: 209-250
Publisher Policy:Reproduced in accordance with the publisher copyright policy

University Staff: Request a correction | Enlighten Editors: Update this record