No Cover Image

Journal article 1389 views 450 downloads

A first-order hyperbolic framework for large strain computational solid dynamics: An upwind cell centred Total Lagrangian scheme

Jibran Haider, Chun Hean Lee Orcid Logo, Antonio Gil Orcid Logo, Javier Bonet

International Journal for Numerical Methods in Engineering, Volume: 109, Issue: 3, Pages: 407 - 456

Swansea University Authors: Jibran Haider, Chun Hean Lee Orcid Logo, Antonio Gil Orcid Logo

  • An_upwind_cell_centred_Total_Lagrangian_Scheme.pdf

    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)

Check full text

DOI (Published version): 10.1002/nme.5293

Abstract

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...

Full description

Published in: International Journal for Numerical Methods in Engineering
ISSN: 0029-5981
Published: 2017
Online Access: Check full text

URI: https://cronfa.swan.ac.uk/Record/cronfa27826
Tags: Add Tag
No Tags, Be the first to tag this record!
Abstract: 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.
College: Faculty of Science and Engineering
Issue: 3
Start Page: 407
End Page: 456