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:
https://csmbrannon.net/2015/07/12/errata-for-two-verification-publications/

Available Online:
http://www.mech.utah.edu/~brannon/pubs/7-2009BrannonLeelavanichkul-IJF.pdf
http://dx.doi.org/10.1007/s10704-009-9398-4

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

Powder metal jet penetration into stressed rock

The Uintah computational framework (UCF) has been adopted for simulation of shaped charge jet penetration and subsequent damage to geological formations.  The Kayenta geomechanics model, as well as a simplified model for shakedown simulations has been  incorporated within the UCF and is undergoing extensive development to enhance it to account for fluid in pore space.

A generic penetration simulation using Uintah

The host code (Uintah) itself has been enhanced to accommodate  material variability and scale effects. Simulations have been performed that import flash X-ray data for the velocity and geometry of a particulated metallic jet so that uncertainty about the jet can be reduced to develop predictive models for target response.  Uintah’s analytical polar decomposition has been replaced with an iterative algorithm to dramatically improve accuracy under large deformations. 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:

http://dx.doi.org/10.1007/s11440-006-0010-z
http://www.mech.utah.edu/~brannon/pubs/7-2006FossumBrannonMechanismDependentViscoplasticity.pdf

Publication: Experimental Assessment of Unvalidated Assumptions in Classical Plasticity Theory

Abstract: 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 fitting 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 flow 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.

You may download the full report here.