A free-Lagrange augmented Godunov method for the simulation of elastic-plastic solids

A Lagrangian finite-volume Godunov scheme is extended to simulate two-dimensional solids in planar geometry. The scheme employs an elastic-perfectly plastic material model, implemented using the method of radial return, and either the 'stiffened' gas or Osborne equation of state to describe the material. The problem of mesh entanglement, common to conventional two-dimensional Lagrangian schemes, is avoided by utilising the free-Lagrange Method. The Lagrangian formulation enables features convecting at the local velocity, such as material interfaces, to be resolved with minimal numerical dissipation. The governing equations are split into separate subproblems and solved sequentially in time using a time-operator split procedure. Local Riemann problems are solved using a two-shock approximate Riemann solver, and piecewise-linear data reconstruction is employed using a MUSCL-based approach to improve spatial accuracy. To illustrate the effectiveness of the technique, numerical simulations are presented and compared with results from commercial fixed-connectivity Lagrangian and smooth particle hydrodynamics solvers (AUTODYN-2D). The simulations comprise the low-velocity impact of an aluminium projectile on a semi-infinite target, the collapse of a thick-walled beryllium cylinder, and the high-velocity impact of cylindrical aluminium and steel projectiles on a thin aluminium target. The analytical solution for the collapse of a thick-walled cylinder is also presented for comparison.

[1]  W. Gust.  High impact deformation of metal cylinders at elevated temperatures , 1982 .

[2]  E. J. Caramana,et al.  Numerical Preservation of Symmetry Properties of Continuum Problems , 1998 .

[3]  T. Ting,et al.  Wave curves for the riemann problem of plane waves in isotropic elastic solids , 1987 .

[4]  Tai-Ping Liu,et al.  The Riemann problem for general systems of conservation laws , 1975 .

[5]  M. Wilkins.  Calculation of Elastic-Plastic Flow , 1963 .

[6]  Phillip Colella,et al.  A higher-order Godunov method for modeling finite deformation in elastic-plastic solids , 1991 .

[7]  G. J. Ball,et al.  Damping of mesh-induced errors in Free-Lagrange simulations of Richtmyer-Meshkov instability , 2000 .

[8]  D. Causon,et al.  A time-splitting approach to solving the Navier-Stokes equations , 1996 .

[9]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[10]  John K. Dukowicz,et al.  A general, non-iterative Riemann solver for Godunov's method☆ , 1985 .

[11]  N. N. Yanenko,et al.  The Method of Fractional Steps , 1971 .

[12]  A. Nádai.  Theory of flow and fracture of solids , 1950 .

[13]  D. Steinberg,et al.  A constitutive model for metals applicable at high-strain rate , 1980 .

[14]  M. B. Tyndall,et al.  Numerical modelling of shocks in solids with elastic-plastic conditions , 1993 .

[15]  Robin Sibson,et al.  Computing Dirichlet Tessellations in the Plane , 1978, Comput. J..

[16]  M. Shashkov,et al.  The Construction of Compatible Hydrodynamics Algorithms Utilizing Conservation of Total Energy , 1998 .

[17]  E. J. Caramana.  Timestep relaxation with symmetry preservation on high aspect-ratio angular or tapered grids , 2001 .

[18]  M. Shashkov,et al.  Elimination of Artificial Grid Distortion and Hourglass-Type Motions by Means of Lagrangian Subzonal Masses and Pressures , 1998 .

[19]  Michael Shearer,et al.  The Riemann problem for a class of conservation laws of mixed type , 1982 .

[20]  J. Monaghan.  On the problem of penetration in particle methods , 1989 .

[21]  John A. Trangenstein,et al.  The Riemann Problem for Longitudinal Motion in an Elastic-Plastic Bar , 1991, SIAM J. Sci. Comput..

[22]  E. Puckett,et al.  A High-Order Godunov Method for Multiple Condensed Phases , 1996 .

[23]  J. Ballmann,et al.  A Riemann solver and a second-order Godunov method for elastic-plastic wave propagation in solids , 1993 .

[24]  Timothy G. Leighton,et al.  Shock-induced collapse of a cylindrical air cavity in water: a Free-Lagrange simulation , 2000 .

[25]  S. Armfield,et al.  The Fractional-Step Method for the Navier-Stokes Equations on Staggered Grids , 1999 .

[26]  G. J. Ball,et al.  A Free-Lagrange method for unsteady compressible flow: simulation of a confined cylindrical blast wave , 1996 .

[27]  Richard Saurel,et al.  Treatment of interface problems with Godunov-type schemes , 1996 .

[28]  T. D. Riney.  CHAPTER V – NUMERICAL EVALUATION OF HYPERVELOCITY IMPACT PHENOMENA* , 1970 .

[29]  Phillip Colella,et al.  A high-order Eulerian Godunov method for elastic-plastic flow in solids , 2001 .

[30]  P. Roache.  Perspective: A Method for Uniform Reporting of Grid Refinement Studies , 1994 .

[31]  Fotis Sotiropoulos,et al.  A Second-Order Godunov Method for Wave Problems in Coupled Solid-Water-Gas Systems , 1999 .

[32]  Burton Wendroff,et al.  The Riemann problem for materials with nonconvex equations of state I: Isentropic flow☆ , 1972 .

[33]  Bradley J. Plohr,et al.  Shockless acceleration of thin plates modeled by a tracked random choice method , 1988 .

[34]  Harold Trease,et al.  The Free-Lagrange Method , 1985 .

[35]  Barbara Lee Keyfitz,et al.  A system of non-strictly hyperbolic conservation laws arising in elasticity theory , 1980 .