Theoretical and numerical comparison of hyperelastic and hypoelastic formulations for Eulerian non-linear elastoplasticity

The aim of this paper is to compare a hyperelastic with a hypoelastic model describing the Eulerian dynamics of solids in the context of non-linear elastoplastic deformations. Specifically, we consider the well-known hypoelastic Wilkins model, which is compared against a hyperelastic model based on the work of Godunov and Romenski. First, we discuss some general conceptual differences between the two approaches. Second, a detailed study of both models is proposed, where differences are made evident at the aid of deriving a hypoelastic-type model corresponding to the hyperelastic model and a particular equation of state used in this paper. Third, using the same high order ADER Finite Volume and Discontinuous Galerkin methods on fixed and moving unstructured meshes for both models, a wide range of numerical benchmark test problems has been solved. The numerical solutions obtained for the two different models are directly compared with each other. For small elastic deformations, the two models produce very similar solutions that are close to each other. However, if large elastic or elastoplastic deformations occur, the solutions present larger differences.Comment: 14 figure

Dynamic finite-strain modelling of the human left ventricle in health and disease using an immersed boundary-finite element method

Detailed models of the biomechanics of the heart are important both for developing improved interventions for patients with heart disease and also for patient risk stratification and treatment planning. For instance, stress distributions in the heart affect cardiac remodelling, but such distributions are not presently accessible in patients. Biomechanical models of the heart offer detailed three-dimensional deformation, stress and strain fields that can supplement conventional clinical data. In this work, we introduce dynamic computational models of the human left ventricle (LV) that are derived from clinical imaging data obtained from a healthy subject and from a patient with a myocardial infarction (MI). Both models incorporate a detailed invariant-based orthotropic description of the passive elasticity of the ventricular myocardium along with a detailed biophysical model of active tension generation in the ventricular muscle. These constitutive models are employed within a dynamic simulation framework that accounts for the inertia of the ventricular muscle and the blood that is based on an immersed boundary (IB) method with a finite element description of the structural mechanics. The geometry of the models is based on data obtained non-invasively by cardiac magnetic resonance (CMR). CMR imaging data are also used to estimate the parameters of the passive and active constitutive models, which are determined so that the simulated end-diastolic and end-systolic volumes agree with the corresponding volumes determined from the CMR imaging studies. Using these models, we simulate LV dynamics from end-diastole to end-systole. The results of our simulations are shown to be in good agreement with subject-specific CMR-derived strain measurements and also with earlier clinical studies on human LV strain distributions

An immersed peridynamics model of fluid-structure interaction accounting for material damage and failure

This paper develops and benchmarks an immersed peridynamics method to simulate the deformation, damage, and failure of hyperelastic materials within a fluid-structure interaction framework. The immersed peridynamics method describes an incompressible structure immersed in a viscous incompressible fluid. It expresses the momentum equation and incompressibility constraint in Eulerian form, and it describes the structural motion and resultant forces in Lagrangian form. Coupling between Eulerian and Lagrangian variables is achieved by integral transforms with Dirac delta function kernels, as in standard immersed boundary methods. The major difference between our approach and conventional immersed boundary methods is that we use peridynamics, instead of classical continuum mechanics, to determine the structural forces. We focus on non-ordinary state-based peridynamic material descriptions that allow us to use a constitutive correspondence framework that can leverage well characterized nonlinear constitutive models of soft materials. The convergence and accuracy of our approach are compared to both conventional and immersed finite element methods using widely used benchmark problems of nonlinear incompressible elasticity. We demonstrate that the immersed peridynamics method yields comparable accuracy with similar numbers of structural degrees of freedom for several choices of the size of the peridynamic horizon. We also demonstrate that the method can generate grid-converged simulations of fluid-driven material damage growth, crack formation and propagation, and rupture under large deformations

Validation of numerical codes for impact and explosion cratering: Impacts on strengthless and metal targets

Over the last few decades, rapid improvement of computer capabilities has allowed impact cratering to be modeled with increasing complexity and realism, and has paved the way for a new era of numerical modeling of the impact process, including full, three-dimensional (3D) simulations. When properly benchmarked and validated against observation, computer models offer a powerful tool for understanding the mechanics of impact crater formation. This work presents results from the first phase of a project to benchmark and validate shock codes. A variety of 2D and 3D codes were used in this study, from commercial products like AUTODYN, to codes developed within the scientific community like SOVA, SPH, ZEUS-MP, iSALE, and codes developed at U.S. National Laboratories like CTH, SAGE/RAGE, and ALE3D. Benchmark calculations of shock wave propagation in aluminum-on-aluminum impacts were performed to examine the agreement between codes for simple idealized problems. The benchmark simulations show that variability in code results is to be expected due to differences in the underlying solution algorithm of each code, artificial stability parameters, spatial and temporal resolution, and material models. Overall, the inter-code variability in peak shock pressure as a function of distance is around 10 to 20%. In general, if the impactor is resolved by at least 20 cells across its radius, the underestimation of peak shock pressure due to spatial resolution is less than 10%. In addition to the benchmark tests, three validation tests were performed to examine the ability of the codes to reproduce the time evolution of crater radius and depth observed in vertical laboratory impacts in water and two well-characterized aluminum alloys. Results from these calculations are in good agreement with experiments. There appears to be a general tendency of shock physics codes to underestimate the radius of the forming crater. Overall, the discrepancy between the model and experiment results is between 10 and 20%, similar to the inter-code variability.The Meteoritics & Planetary Science archives are made available by the Meteoritical Society and the University of Arizona Libraries. Contact [email protected] for further information.Migrated from OJS platform February 202

Studies on strain localization, ductile fracture and damage in structural metals

One of the most important limit states in structural metals is ductile fracture, and the prediction of ductile fracture is of great importance in many engineering applications. The overall objective of the research reported in this dissertation is to advance the understanding and modeling of ductile fracture in metals. This research addresses three main issues: micromechanical modeling of ductile fracture, the development of a micromechanics-based ductile fracture model and its numerical implementation, and a numerical investigation of geometry and damage induced strain localization based on a nonlocal formulation. It has long been recognized that stress triaxiality is a key parameter affecting initiation of ductile fracture. More recently, shear stress has been identified as another important parameter, in addition to stress triaxiality, that influences the process of ductile fracture. In this research, a micromechanics-based model is proposed for predicting initiation of ductile fracture that couples both stress triaxiality and shear stress. The new model is based on a combination of the existing Rice-Tracey and modified maximum shear stress models. The new model is applied to construct the fracture locus of different types of metal alloys and is used to predict fracture initiation by numerical tools. The predicted results are in good agreement with experimental data reported in literature that covers a wide range of triaxialities and shear stress. Another portion of this research, within the framework of micromechanics, investigated the effect of combined normal and shear stress components on micro-void evolution and material behavior. This work involved finite element modeling of a cubic unit cell associated with a spherical void. The results show that the void growth process and macroscopic stress-strain response is highly dependent on the shear stress component. At different ranges of triaxialities, and with different void growth and coalescence mechanisms, shear stress has an important effect on the ductile fracture process. Numerical modeling of strain localization in ductile metals based on standard continuum mechanics exhibits non-convergent mesh sensitivity. This issue is addressed in the final portion of this research. A one-dimensional model based on the nonlocal theory is proposed to analyze geometry-induced strain localization, i.e., necking in structural metals. A nonlocal continuum damage model using the same enhanced continuum law is developed to deal with the damage induced strain localization in metals. Both models provide encouraging performance in eliminating the non-convergent mesh sensitivity problem. Such improved strain localization modeling techniques show potential to be useful for further exploration of ductile fracture phenomena.Civil, Architectural, and Environmental Engineerin

The Material Point Method for Solid and Fluid Simulation

The Material Point Method (MPM) has shown its high potential for physics-based simulation in the area of computer graphics. In this dissertation, we introduce a couple of improvements to the traditional MPM for different applications and demonstrate the advantages of our methods over the previous methods.First, we present a generalized transfer scheme for the hybrid Eulerian/Lagrangian method: the Polynomial Particle-In-Cell Method (PolyPIC). PolyPIC improves kinetic energy conservation during transfers, which leads to better vorticity resolution in fluid simulations and less numerical damping in elastoplasticity simulations. Our transfers are designed to select particle-wise polynomial approximations to the grid velocity that are optimal in the local mass-weighted L2 norm. Indeed our notion of transfers reproduces the original Particle-In-Cell Method (PIC) and recent Affine Particle-In-Cell Method (APIC). Furthermore, we derive a polynomial basis that is mass orthogonal to facilitate the rapid solution of the optimality condition. Our method applies to both of the collocated and staggered grid.As the second contribution, we present a novel method for the simulation of thin shells with frictional contact using a combination of MPM and subdivision finite elements. The shell kinematics are assumed to follow a continuum shell model which is decomposed into a Kirchhoff-Love motion that rotates the mid-surface normals followed by shearing and compression/extension of the material along the mid-surface normal. We use this decomposition to design an elastoplastic constitutive model to resolve frictional contact by decoupling resistance to contact and shearing from the bending resistance components of stress. We show that by resolving frictional contact with a continuum approach, our hybrid Lagrangian/Eulerian approach is capable of simulating challenging shell contact scenarios with hundreds of thousands to millions of degrees of freedom. Without the need for collision detection or resolution, our method runs in a few minutes per frame in these high-resolution examples. Furthermore, we show that our technique naturally couples with other traditional MPM methods for simulating granular and related materials.In the third part, we present a new hybrid Lagrangian Material Point Method for simulating volumetric objects with frictional contact. The resolution of frictional contact in the thin shell simulation cannot be generalized to the case of volumetric materials directly. Also, even though MPM allows for the natural simulation of hyperelastic materials represented with Lagrangian meshes, it usually coarsens the degrees of freedom of the Lagrangian mesh and can lead to artifacts, e.g., numerical cohesion. We demonstrate that our hybrid method can efficiently resolve these issues. We show the efficacy of our technique with examples that involve elastic soft tissues coupled with kinematic skeletons, extreme deformation, and coupling with various elastoplastic materials. Our approach also naturally allows for two-way rigid body coupling

A Material Point Method for Simulating Frictional Contact with Diverse Materials

We present an extension to the Material Point Method (MPM) for simulating elastic objects with various co-dimensions like hair (1D), thin shells (2D), and volumetric objects (3D). We simulate thin shells with frictional contact using a combination of MPM and subdivision finite elements. The shell kinematics are assumed to follow a continuum shell model which is decomposed into a Kirchhoff-Love motion that rotates the mid-surface normals followed by shearing and compression/extension of the material along the mid-surface normal. We use this decomposition to design an elastoplastic constitutive model to resolve frictional contact by decoupling resistance to contact and shearing from the bending resistance components of stress. We show that by resolving frictional contact with a continuum approach, our hybrid Lagrangian/Eulerian approach is capable of simulating challenging shell contact scenarios with hundreds of thousands to millions of degrees of freedom. Furthermore our technique naturally couples with other traditional MPM methods for simulating granular materials. Without the need for collision detection or resolution, our method runs in a few minutes per frame in these high resolution examples. For the simulation of hair and volumetric elastic objects, we utilize a Lagrangian mesh for internal force computation and an Eulerian mesh for self collision as well as coupling with external materials. While the updated Lagrangian discretization where the Eulerian grid degrees of freedom are used to take variations of the potential energy is effective in simulating thin shells, its frictional contact response strategy does not generalize to volumetric objects. Therefore, we develop a hybrid approach that retains Lagrangian degrees of freedom while still allowing for natural coupling with other materials simulated with traditional MPM. We demonstrate the efficacy of our technique with examples that involve elastic soft tissues coupled with kinematic skeletons, extreme deformation, and coupling with multiple elastoplastic materials. Our approach also naturally allows for two-way rigid body coupling