A scanned copy is available here.
ABSTRACT (from OCR so has some typos): The influence of non-classical elastic-plastic constitutive features on dynamically moving discontinuities in stress, strain, and material velocity is investigated. Non-classical behavior here includes non-normality of the plastic strain increment to the yield surface, plastic compressibility, pressure sensitivity of yield, and dependence of the elastic moduli on plastic strain. DRUGAN and SHEN’s (1987) analysis of dynamically moving discontinuities with strain as well as stress jumps in classical materials is shown to be valid for a broad class of non-associative material models until deviation from normality exceeds a critical (non-infinitesimal) level. For these non-classical materials, an inequality that bounds the magnitude of the stress jump is derived, which is information not obtainable from a standard spectral analysis of a shock. For the special case of stress discontinuities with continuous strain, or for quasi-static deformations, this inequality is shown to rule out jumps in specific projections of the stress tensor unless the non-normality is sufficiently large. These results invalidate a … claim in the literature that an infinitesimal amount of non-normality permits moving surfaces of discontinuity in stress (with no strain jump) near the tip of a dynamically advancing crack tip.
Using a very general plastic constitutive law that subsumes most non-classical (and classical) descriptions currently in use, a complete closed-form solution is obtained for the plastic wave speeds and eigenvectors. A novel feature of the analysis is the clarity and completeness of the solutions. If the elastic part of the response is isotropic, one plastic wave speed equals the elastic shear wave speed, while the other two possible wave speeds depend in general on the stress and plastic strain within the shock transition layer. Concise necessary and sufficient conditions for real eigenvalues and for vanishing eigenvalues are derived. The real eigenvalues are classified by numerical sign and ordering relative to the elastic eigenvalues. The geometric multiplicity of plastic eigenvectors associated with elastic eigenvalues is shown to depend on the stress state within the shock transition layer. These solutions, several of which hold for arbitrary elastic anisotropy, are also applicable to acceleration waves and localization problems, and to materials with dependence of the elastic moduli on plastic strain. Such elastic–plastic coupling is shown to imply a non-self-adjoint fourth order tangent stiffness tensor even if the plastic constitutive law is associative.
A scanned copy is available here.
This 2007 Book Chapter on the basics of plasticity theory reviews the terminology and governing equations of plasticity, with emphasis on amending misconceptions, providing physical insights, and outlining computational algorithms. Plasticity theory is part of a larger class of material models in which a pronounced change in material response occurs when the stress (or strain) reaches a critical threshold level. If the stress state is subcritical, then the material is modeled by classical elasticity. The bound- ary of the subcritical (elastic) stress states is called the yield surface. Plasticity equations apply if continuing to apply elasticity theory would predict stress states that extend beyond this the yield surface. The onset of plasticity is typically characterized by a pronounced slope change in a stress–strain dia-gram, but load reversals in experiments are necessary to verify that the slope change is not merely nonlinear elasticity or reversible phase transformation.
The threshold yield surface can appear to be significantly affected by the loading rate, which has a dominant effect in shock physics applications.
In addition to providing a much-needed tutorial survey of the governing equations and their solution (defining Lode angle and other Lode invariants and addressing the surprisingly persistent myth that closest-point return satisfies the governing equations), this book chapter includes some distinctive contributions such as a simple 2d analog of plasticity that exhibits the same basic features of plasticity (such as existence of a “yield” surface with associative flow and vertex theory), an extended discussion of apparent nonassociativity, stability and uniqueness concerns about nonassociativity, and a summary of apparent plastic wave speeds in relation to elastic wave speeds (especially noting that non-associativity admits plastic waves that travel faster than elastic waves).
For the full manuscript with errata, click 2007 Book Chapter on the basics of plasticity theory.
R. Brannon, J.A. Burghardt, D. Bronowski, and S. Bauer
Common isotropic yield surfaces. Von Mises and Drucker-Prager models are often used for metals. Gurson’s function, and others like it, are used for porous media. Tresca and Mohr-Coulomb models approximate the yield threshold for brittle media. Fossum’s model, and others like it, combine these features to model realistic geological media.
This report investigates the validity of several key assumptions in classical plasticity theory regarding material response to changes in the loading direction. Three metals, two rock types, and one ceramic were subjected to non-standard loading directions, and the resulting strain response increments were displayed in Gudehus diagrams to illustrate the approximation error of classical plasticity theories. A rigorous mathematical framework for ﬁtting classical theories to the data,thus quantifying the error, is provided. Further data analysis techniques are presented that allow testing for the effect of changes in loading direction without having to use a new sample and for inferring the yield normal and ﬂow directions without having to measure the yield surface. Though the data are inconclusive, there is indication that classical, incrementally linear, plasticity theory may be inadequate over a certain range of loading directions. This range of loading directions also coincides with loading directions that are known to produce a physically inadmissible instability for any nonassociative plasticity model.
A plot of the frequency-dependent wave propagation velocity for the case study problem with an overlocal plasticity model, with the elastic and local hardening wave speeds shown for reference (left). Stress histories using an overlocal plasticity model with a nonlocal length scale of 1m and a mesh resolution of 0.125m (right)
The following series of three articles (with common authors J. Burghardt and R. Brannon of the University of Utah) describes a state of insufficient experimental validation of conventional formulations of nonassociative plasticity (AKA nonassociated and non-normality). This work provides a confirmation that such models theoretically admit negative net work in closed strain cycles, but this simple prediction has never been validated or disproved in the laboratory!
- An early (mostly failed) attempt at experimental investigation of unvalidated plasticity assumptions (click to view),
- A simple case study confirming that nonassociativity can cause non-unique and unstable solutions to wave motion problems (click to view),
- An extensive study showing that features like rate dependence, hardening, etc. do not eliminate the instability and also showing that it is NOT related to conventional localization (click to view).
A. F. Fossum and R. M. Brannon
This paper summarizes the results of a theoretical and experimental program at Sandia National Laboratories aimed at identifying and modeling key physical features of rocks and rock-like materials at the laboratory scale over a broad range of strain rates. The mathematical development of a constitutive model is discussed and model predictions versus experimental data are given for a suite of laboratory tests. Concurrent pore collapse and cracking at the microscale are seen as competitive micromechanisms that give rise to the well-known macroscale phenomenon of a transition from volumetric compaction to dilatation under quasistatic triaxial compression. For high-rate loading, this competition between pore collapse and microcracking also seems to account for recently identiﬁed differences in strain-rate sensitivity between uniaxial-strain ‘‘plate slap’’ data compared to uniaxial-stress Kolsky bar data. A description is given of how this work supports ongoing efforts to develop a predictive capability in simulating deformation and failure of natural geological materials, including those that contain structural features such as joints and other spatial heterogeneities.
Tough disk impacting brittle disk
Below are links to two simulations of disks colliding. The first is elastic and the second uses a fracture model with spatially variable strength based on a scale-dependent Weibull realization. Both take advantage of the automatic contact property of the MPM.
WeibConstMovie: disks colliding without fracture
WeibPerturbedGood: disks colliding with heterogeneous fracture
This basic capability to support statistically variable strength in a damage model has been extended to the Kayenta plasticity model in Uintah.
T.J. Fuller, R.M. Brannon, O.E. Strack, J.E. Bishop
Displacement profile for Thermo-Kayenta at the end of the simulation. the red dots represent the experimental profiles
A persistent challenge in simulating damage of natural geological materials, as well as rock-like engineered materials, is the development of efficient and accurate constitutive models.The common feature for these brittle and quasi-brittle materials are the presence of flaws such as porosity and network of microcracks. The desired models need to be able to predict the material responses over a wide range of porosities and strain rate. Kayenta  (formerly called the Sandia GeoModel) is a unifi ed general-purpose constitutive model that strikes a balance between rst-principles micromechanics and phenomenological or semi-empirical modeling strategies. However, despite its sophistication and ability to reduce to several classical plasticity theories, Kayenta is incapable of modeling deformation of ductile materials in which deformation is dominated by dislocation generation and movement which can lead to signi cant heating. This stems from Kayenta’s roots as a geological model, where heating due to inelastic deformation is often neglected or presumed to be incorporated implicitly through the elastic moduli.The sophistication of Kayenta and its large set of extensive features, however, make Kayenta an attractive candidate model to which thermal eff ects can be added. This report outlines the initial work in doing just that, extending the capabilities of Kayenta to include deformation of ductile materials, for which thermal e ffects cannot be neglected. Thermal e ffects are included based on an assumption of adiabatic loading by computing the bulk and thermal responses of the material with the Kerley Mie-Gruneisen equation of state and adjusting the yield surface according to the updated thermal state. This new version of Kayenta, referred to as Thermo-Kayenta throughout this report, is capable of reducing to classical Johnson-Cook plasticity in special case single element simulations and has been used to obtain reasonable results in more complicated Taylor impact simulations in LS-Dyna. Despite these successes, however, Thermo-Kayenta requires additional re nement for it to be consistent in the thermodynamic sense and for it to be considered superior to other, more mature thermoplastic models. The initial thermal development, results, and required refinements are all detailed in the following report.
J. Burghardt, B. Leavy, J. Guilkey, Z. Xue, R. Brannon
The capability of the generalized interpolation material point (GIMP) method in simulation of penetration events is investigated. A series of experiments was performed wherein a shaped charge jet penetrates into a stack of aluminum plates. Electronic switches were used to measure the penetration time history. Flash x-ray techniques were used to measure the density,length, radius and velocity of the shaped charge jet. Simulations of the penetration event were performed using the Uintah MPM/GIMP code with several different models of the shaped charge jet being used. The predicted penetration time history for each jet model is compared with the experimentally observed penetration history. It was found that the characteristics of the predicted penetration were dependent on the way that the jet data are translated to a discrete description. The discrete jet descriptions were modified such that the predicted penetration histories fell very close to the range of the experimental data. In comparing the various discrete jet descriptions it was found that the cumulative kinetic energy flux curve represents an important way of characterizing the penetration characteristics of the jet. The GIMP method was found to be well suited for simulation of high rate penetration events.
R.M. Brannon, A.F. Fossum, and O.E. Strack
Kayenta continuous yield surface. (a) three-dimensional view in principal stress space, (b) the meridional “side” view (thick line), and (c) the octahedral view
The physical foundations and domain of applicability of the Kayenta constitutive model are presented along with descriptions of the source code and user instructions. Kayenta, which is an outgrowth of the Sandia GeoModel, includes features and fitting functions appropriate to a broad class of materials including rocks, rock-like engineered materials (such as concretes and ceramics),and metals. Fundamentally, Kayenta is a computational framework for generalized plasticity models. As such, it includes a yield surface, but the term“yield” is generalized to include any form of inelastic material response including microcrack growth and pore collapse. Kayenta supports optional anisotropic elasticity associated with ubiquitous joint sets. Kayenta support optional deformation-induced anisotropy through kinematic hardening (inwhich the initially isotropic yield surface is permitted to translate in deviatoric stress space to model Bauschinger effects). The governing equations are otherwise isotropic. Because Kayenta is a unification and generalization of simple models, it can be run using as few as 2 parameters (for linear elasticity) to as many as 40 material and control parameters in the exceptionally rare case when all features are used. For high-strain-rate applications, Kayenta support rate dependence through an overstress model. Isotropic damage is model through loss of stiffness and strength.
R.M. Brannon and S. Leelavanichkul
RHT Model: Contour plots of damage: side, front, and back view of the target (top to bottom).
Four conventional damage plasticity models for concrete, the Karagozian and Case model (K&C),the Riedel-Hiermaier-Thoma model (RHT), the Brannon-Fossum model (BF1), and the Continuous Surface Cap Model (CSCM) are compared. The K&C and RHT models have been used in commercial finite element programs many years, whereas the BF1 and CSCM models are relatively new. All four models are essentially isotropic plasticity models for which plasticity is regarded as any form of inelasticity. All of the models support nonlinear elasticity, but with different formulations.All four models employ three shear strength surfaces. The yield surface bounds an evolving set of elastically obtainable stress states. The limit surface bounds stress states that can be reached by any means (elastic or plastic). To model softening, it is recognized that some stress states might be reached once, but, because of irreversible damage, might not be achievable again. In other words, softening is the process of collapse of the limit surface, ultimately down to a final residual surface for fully failed material. The four models being compared differ in their softening evolution equations, as well as in their equations used to degrade the elastic stiffness. For all four models, the strength surfaces are cast in stress space. For all four models, it is recognized that scale effects are important for softening, but the models differ significantly in their approaches. The K&C documentation, for example, mentions that a particular material parameter affecting the damage evolution rate must be set by the user according to the mesh size to preserve energy to failure. Similarly, the BF1 model presumes that all material parameters are set to values appropriate to the scale of the element, and automated assignment of scale-appropriate values is available only through an enhanced implementation of BF1 (called BFS) that regards scale effects to be coupled to statistical variability of material properties. The RHT model appears to similarly support optional uncertainty and automated settings for scale-dependent material parameters. The K&C, RHT, and CSCM models support rate dependence by allowing the strength to be a function of strain rate, whereas the BF1 model uses Duvaut-Lion viscoplasticity theory to give a smoother prediction of transient effects. During softening, all four models require a certain amount of strain to develop before allowing significant damage accumulation. For the K&C, RHT, and CSCM models, the strain-to-failure is tied to fracture energy release, whereas a similar effect is achieved indirectly in the BF1 model by a time-based criterion that is tied to crack propagation speed.