A two‐step Taylor‐Galerkin formulation for fast dynamics. Issue 3 (2014)

• Record Type: Journal Article
• Authors: Izian Abd. Karim; Chun Hean Lee; Antonio J. Gil; Javier Bonet
• Is Part Of: Engineering computations. Volume 31:Issue 3(2014)
• Page Start: 366
• Page End: 387
• Publication Date: 2014
• Subjects: Computer-aided engineering -- Periodicals; Computer graphics -- Periodicals; 620.00285
• DOI: 10.1108/EC-12-2012-0319
• Languages: English
• ISSNs: 0264-4401

Abstract

<abstract> <title> <x xml:space="preserve"> Abstract </x> </title> <p> <bold>Purpose</bold> – The purpose of this paper is to present a new stabilised low‐order finite element methodology for large strain fast dynamics. <bold>Design/methodology/approach</bold> – The numerical technique describing the motion is formulated upon the mixed set of first‐order hyperbolic conservation laws already presented by Lee <italic>et al.</italic> (2013) where the main variables are the linear momentum, the deformation gradient tensor and the total energy. The mixed formulation is discretised using the standard explicit two‐step Taylor‐Galerkin (2TG) approach, which has been successfully employed in computational fluid dynamics (CFD). Unfortunately, the results display non‐physical spurious (or hourglassing) modes, leading to the breakdown of the numerical scheme. For this reason, the 2TG methodology is further improved by means of two ingredients, namely a curl‐free projection of the deformation gradient tensor and the inclusion of an additional stiffness stabilisation term. <bold>Findings</bold> – A series of numerical examples are carried out drawing key comparisons between the proposed formulation and some other recently published numerical techniques. <bold>Originality/value</bold> – Both velocities (or displacements) and stresses display the same rate of convergence, which proves ideal in the case of industrial applications where low‐order discretisations tend to be preferred. The enhancements introduced in this paper enable the use of linear triangular (or bilinear quadrilateral) elements in two dimensional nearly incompressible dynamics applications without locking difficulties. In addition, an artificial viscosity term has been added into the formulation to eliminate the appearance of spurious oscillations in the vicinity of sharp spatial gradients induced by shocks.</p> <ack> <title> <x xml:space="preserve"> Acknowledgements </x> </title> <p>The first author acknowledges the financial support received from the Ministry of Higher Education (Malaysia) and the University of Putra Malaysia (Serdang). The third author acknowledges the financial support received through "The Leverhulme Prize" awarded by The Leverhulme Trust, United Kingdom.</p> </ack> </abstract>