Publication: Survey of Four Damage Models for Concrete

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.

Available Online:

Publication: A multi-stage return algorithm for solving the classical damage component of constitutive models for rocks, ceramics, and other rock-like media

R. M. Brannon and S. Leelavanichkul

Octahedral isosurfaces for a) the unacceptable, b) the admissible, and c) the admissible

Classical plasticity and damage models for porous quasi-brittle media usually suffer from mathematical defects such as non-convergence and nonuniqueness.Yield or damage functions for porous quasi-brittle media often have yield functions with contours so distorted that following those contours to the yield surface in a return algorithm can take the solution to a false elastic domain. A steepest-descent return algorithm must include iterative corrections; otherwise,the solution is non-unique because contours of any yield function are non-unique. A multi-stage algorithm has been developed to address both spurious convergence and non-uniqueness, as well as to improve efficiency. The region of pathological isosurfaces is masked by first returning the stress state to the Drucker–Prager surface circumscribing the actual yield surface. From there, steepest-descent is used to locate a point on the yield surface. This first-stage solution,which is extremely efficient because it is applied in a 2D subspace, is generally not the correct solution,but it is used to estimate the correct return direction.The first-stage solution is projected onto the estimated correct return direction in 6D stress space. Third invariant dependence and anisotropy are accommodated in this second-stage correction. The projection operation introduces errors associated with yield surface curvature,so the two-stage iteration is applied repeatedly to converge. Regions of extremely high curvature are detected and handled separately using an approximation to vertex theory. The multi-stage return is applied holding internal variables constant to produce a non-hardening solution. To account for hardening from pore collapse (or softening from damage), geometrical arguments are used to clearly illustrate the appropriate scaling of the non-hardening solution needed to obtain the hardening (or softening) solution.

For errata (transcription errors in two of the verification solutions), please see:

Available Online:

Publication: Decomposition and Visualization of Fourth-Order Elastic-Plastic Tensors

A.G. Neeman; R.M. Brannon; B. Jeremic; A. Van Gelderand;  A. Pang

Top view (Z from above) of eigentensors for Drucker-Prager material, time step 124, colored by minimum stretch eigenvalue.

Visualization of fourth-order tensors from solid mechanics has not been explored in depth previously. Challenges include the large number of components (3x3x3x3 for 3D), loss of major symmetry and loss of positive definiteness(with possibly zero or negative eigenvalues). This paper presents a decomposition of fourth-order tensors that facilitates their visualization and understanding. Fourth-order tensors are used to represent a solid’s stiffness.The stiffness tensor represents the relationship between increments of stress and increments of strain. Visualizing stiffness is important to understand the changing state of solids during plastification and failure. In this work,we present a method to reduce the number of stiffness components to second-order 3×3 tensors for visualization.The reduction is based on polar decomposition, followed by eigen-decomposition on the polar “stretch”. If any resulting eigenvalue is significantly lower than the others, the material has softened in that eigen-direction. The associated second-order eigentensor represents the mode of stress (such as compression, tension, shear, or some combination of these) to which the material becomes vulnerable. Thus we can visualize the physical meaning of plastification with techniques for visualizing second-order symmetric tensors.

Available Online:



Research: Instability of *ANY* nonassociative plasticity model

The CSM group has independently confirmed  a case study demonstrating the truth of a claim in the literature that any non-associative rate-independent model admits a non-physical dynamic achronistity instability. By stimulating a non-associative material in the “Sandler-Rubin wedge” (above yield but below the flow surface), plastic waves are generated that travel faster than elastic waves, thus introducing a negative net work in a closed strain cycle that essentially feeds energy into a propagating wave to produce unbounded increases in displacement with time.

Sandler-Rubin instability: an infinitesimal pulse grows as it propagates

Continue reading

Nonclassical plasticity validation

Analysis and computations have been performed by the Utah CSM group to support experimental investigations of unvalidated assumptions in plasticity theory. The primary untested assumption is that of a regular flow rule in which it is often assumed that the direction of the inelastic strain increment is unaffected by the total strain increment itself. To support laboratory testing of this hypothesis, the general equations of classical plasticity theory were simplified for the case of axisymmetric loading to provide experimentalists with two-parameter control of the axial and lateral stress increments corresponding to a specified loading trajectory in stress space. Loading programs involving changes in loading directions were designed. New methods for analyzing the data via a moving least squares fit to tensor-valued input-output data were used to quantitatively infer the apparent plastic tangent modulus matrix and thereby detect violations of the regular flow rule. Loading programs were designed for validating isotropic cap hardening models by directly measuring the effect of shear loading on the hydrostatic elastic limit.

UofU Contributors/collaborators:
Michael Braginski (postdoc, Mech. Engr., UofU)
Jeff Burghardt (PhD student, Mech. Engr., UofU)

External collaborators/mentors:
Stephen Bauer (Manager, Sandia National Labs geomechanics testing lab)
David Bronowski (Sandia geomechanics lab technician)
Erik Strack (Manager, Sandia Labs Computational Physics)

Publication: On a viscoplastic model for rocks with mechanism-dependent characteristic times

A.F. Fossum and R.M. Brannon (2006)

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 identified 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.

Available online: