Journal article 1136 views 297 downloads
A first-order hyperbolic framework for large strain computational solid dynamics: An upwind cell centred Total Lagrangian scheme
International Journal for Numerical Methods in Engineering, Volume: 109, Issue: 3, Pages: 407 - 456
PDF | Accepted Manuscript
This is the peer reviewed version of the article which has been published in final form. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions.Download (19.05MB)
This paper builds on recent work developed by the authors for the numerical analysis of large strain solid dynamics, by introducing an upwind cell centred hexahedral finite volume framework implemented within the open source code OpenFOAM [http://www.openfoam.com/]. In Lee, Gil and Bonet (2013), a f...
|Published in:||International Journal for Numerical Methods in Engineering|
Check full text
No Tags, Be the first to tag this record!
This paper builds on recent work developed by the authors for the numerical analysis of large strain solid dynamics, by introducing an upwind cell centred hexahedral finite volume framework implemented within the open source code OpenFOAM [http://www.openfoam.com/]. In Lee, Gil and Bonet (2013), a first-order hyperbolic system of conservation laws was introduced in terms of the linear momentum and the deformation gradient tensor of the system, leading to excellent behaviour in two-dimensional bending dominated nearly incompressible scenarios. The main aim of this paper is the extension of this algorithm into three dimensions, its tailor-made implementation into OpenFOAM and the enhancement of the formulation with three key novelties. First, the introduction of two different strategies in order to ensure the satisfaction of the underlying involutions of the system, that is, that the deformation gradient tensor must be curl-free throughout the deformation process. Second, the use of a discrete angular momentum projection algorithm and a monolithic Total Variation Diminishing Runge-Kutta time integrator combined in order to guarantee the conservation of angular momentum. Third, and for comparison purposes, an adapted Total Lagrangian version of the hyperelastic-GLACE nodal scheme of Kluth and Després (2010) is presented. A series of challenging numerical examples are examined in order to assess the robustness and accuracy of the proposed algorithm, benchmarking it against an ample spectrum of alternative numerical strategies developed by the authors in recent publications.
Faculty of Science and Engineering