Journal
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING
Volume 340, Issue -, Pages 684-727Publisher
ELSEVIER SCIENCE SA
DOI: 10.1016/j.cma.2018.06.010
Keywords
First order conservation laws; Large strain solid dynamics; Finite volume method; Riemann solver; OpenFOAM
Funding
- The Erasmus Mundus Joint Doctorate SEED programme
- European Regional Development Fund (ERDF)
- Ser Cymru National Research Network for Advanced Engineering and Materials, United Kingdom
Ask authors/readers for more resources
The paper presents a new computational framework for the numerical simulation of fast large strain solid dynamics, with particular emphasis on the treatment of near incompressibility. A complete set of first order hyperbolic conservation equations expressed in terms of the linear momentum and the minors of the deformation (namely the deformation gradient, its co-factor and its Jacobian), in conjunction with a polyconvex nearly incompressible constitutive law, is presented. Taking advantage of this elegant formalism, alternative implementations in terms of entropy-conjugate variables are also possible, through suitable symmetrisation of the original system of conservation variables. From the spatial discretisation standpoint, modern Computational Fluid Dynamics code OpenFOAM [http://www.openfoam.com/] is here adapted to the field of solid mechanics, with the aim to bridge the gap between computational fluid and solid dynamics. A cell centred finite volume algorithm is employed and suitably adapted. Naturally, discontinuity of the conservation variables across control volume interfaces leads to a Riemann problem, whose resolution requires special attention when attempting to model materials with predominant nearly incompressible behaviour (kappa/mu >= 500). For this reason, an acoustic Riemann solver combined with a preconditioning procedure is introduced. In addition, a global a posteriori angular momentum projection procedure proposed in Haider et al. (2017) is also presented and adapted to a Total Lagrangian version of the nodal scheme of Kluth and Despres (2010) used in this paper for comparison purposes. Finally, a series of challenging numerical examples is examined in order to assess the robustness and applicability of the proposed methodology with an eye on large scale simulation in future works. (C) 2018 Elsevier B.Y. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available