A Total Lagrangian based method for recovering the un-deformed configuration in finite elasticity

Abstract The problem of finding the un-deformed configuration of an elastic body, when the deformed configuration and the loads are known, occurs in many engineering applications. Standard solution methods for such problems include conservation laws based on Eshelby’s energy–momentum tensor and re-parameterization of the standard equilibrium equations. In this paper we present a different method for solving such problems, based on a re-parameterization of the nodal forces using the Total Lagrangian formulation. The obtained nonlinear system of equations describing equilibrium can be solved using either Newton–Raphson or an explicit dynamic relaxation algorithm. The solution method requires only minor modifications to similar algorithms designed for forward motion calculations. Several examples involving large deformations and different boundary conditions and loads are presented.

[1]  Ted Belytschko,et al.  A survey of numerical methods and computer programs for dynamic structural analysis , 1976 .

[2]  Karol Miller,et al.  Computation of intra-operative brain shift using dynamic relaxation. , 2009, Computer methods in applied mechanics and engineering.

[3]  M. L. Raghavan,et al.  Computational method of inverse elastostatics for anisotropic hyperelastic solids , 2007 .

[4]  Richard T. Schield.  Inverse deformation results in finite elasticity , 1967 .

[5]  J. D. Eshelby.  The elastic energy-momentum tensor , 1975 .

[6]  Karol Miller,et al.  Real-Time Nonlinear Finite Element Computations on GPU - Application to Neurosurgical Simulation. , 2010, Computer methods in applied mechanics and engineering.

[7]  K. Bathe.  Finite Element Procedures , 1995 .

[8]  P. Chadwick,et al.  Applications of an energy-momentum tensor in non-linear elastostatics , 1975 .

[9]  S. Govindjee,et al.  Computational methods for inverse finite elastostatics , 1996 .

[10]  Karol Miller,et al.  Real-Time Nonlinear Finite Element Computations on GPU: Handling of Different Element Types , 2011 .

[11]  M. L. Raghavan,et al.  Inverse elastostatic stress analysis in pre-deformed biological structures: Demonstration using abdominal aortic aneurysms. , 2007, Journal of biomechanics.

[12]  J. Rodríguez,et al.  A Pull-Back Algorithm to Determine the Unloaded Vascular Geometry in Anisotropic Hyperelastic AAA Passive Mechanics , 2013, Annals of Biomedical Engineering.

[13]  K. Miller,et al.  Total Lagrangian explicit dynamics finite element algorithm for computing soft tissue deformation , 2006 .

[14]  Karol Miller,et al.  On the prospect of patient-specific biomechanics without patient-specific properties of tissues. , 2013, Journal of the mechanical behavior of biomedical materials.

[15]  Karol Miller,et al.  An adaptive dynamic relaxation method for solving nonlinear finite element problems. Application to brain shift estimation , 2011, International journal for numerical methods in biomedical engineering.

[16]  M. Nash,et al.  Determining the finite elasticity reference state from a loaded configuration , 2007 .

[17]  Kathrin Abendroth,et al.  Nonlinear Finite Elements For Continua And Structures , 2016 .

[18]  S. Govindjee,et al.  Computational methods for inverse de-formations in quasi-incompressible nite elasticity , 1998 .

[19]  Pierre Badel,et al.  In vitro analysis of localized aneurysm rupture. , 2014, Journal of biomechanics.

[20]  Karol Miller,et al.  Non-locking Tetrahedral Finite Element for Surgical Simulation. , 2009, Communications in numerical methods in engineering.