Publication: Advances in X-ray Computed Tomography Diagnostics of Ballistic Impact Damage

J.M. Wells and R.M. Brannon

Dynamic indentation of SiC-N ceramic by a tungsten carbide sphere. Left: experimentally observed impact crater and radial cracking (both highlighted for clarity). Middle: BFS model prediction of externally visible damage. Right: prediction of internal damage (suitable for validation against XCT data).

With the relatively recent introduction of quantitative and volumetric X-ray computedtomography (XCT) applied to ballistic impact damage diagnostics, significant inroads have beenmade in expanding our knowledge base of the morphological variants of physical impactdamage. Yet, the current state of the art in computational and simulation modeling of terminalballistic performance remains predominantly focused on the penetration phenomenon, withoutdetailed consideration of the physical characteristics of actual impact damage. Similarly, armorceramic material improvements appear more focused on penetration resistance than on improved intrinsic damage tolerance and damage resistance. Basically, these approaches minimizeour understanding of the potential influence that impact damage may play in the mitigation orprevention of ballistic penetration. Examples of current capabilities of XCT characterization,quantification, and visualization of complex impact damage variants are demonstrated anddiscussed for impacted ceramic and metallic terminal ballistic target materials. Potential benefitsof incorporating such impact damage diagnostics in future ballistic computational modeling arealso briefly discussed.

Available Online:

http://dx.doi.org/10.1007/s11661-007-9304-5
http://www.mech.utah.edu/~brannon/pubs/7-2007WellsBrannonAdvancesInXrayComputedTomographyDiagnosticsOfBallisticDamage.pdf

CPDI shape functions for the Material Point Method

In a conventional MPM formulation, the shape functions on the grid are the same as in a traditional FEM solution. In the CPDI, the shape functions on the grid are replaced by alternative (and still linearly complete*) shape functions, given by piecewise linear interpolations of the traditional FEM shape functions to the boundaries of the particles.  This change provides FEM-level accuracy in moderately deforming regions while retaining the attractive feature of MPM that particles can move arbitrarily relative to one another in massively deforming regions (provided, of course, that the deformation is updated in a manner compatible with the constitutive model).

In the images below, the shaded regions are the traditional FEM “tent” linear shape functions in 1-D, and the solid lines are the CPDI interpolated shape functions, which clearly change based on particle position relative to the grid.  Both the traditional FEM tent functions and these new CPDI functions are linearly complete (i.e., they can exactly fit any affine function). The tremendous advantage of CPDI is that the basis functions are extraordinarily simple over a particle domain, thus facilitating exact and efficient evaluation of integrals over particle domains.

Continue reading

Research: Radial cracking as a means to infer aleatory uncertainty parameters

Aleatory uncertainty in constitutive modeling refers to the intrinsic variability in material properties caused by differences in micromorphology (e.g., grain orientation or size, microcracks, inclusions, etc.) from sample to sample. Accordingly, a numerical simulation of a nominally axisymmetric problem must be run in full 3D (non-axisymmetric) mode if there is any possibility of a bifurcation from stability.

Dynamic indentation experiments, in which a spherical ball impacts to top free surface of a cylindrical specimen, nicely illustrate that fracture properties must have spatial variability — in fact, the intrinsic instability that leads to radial cracking is regarded by the Utah CSM group as a potential inexpensive means of inferring the spatial frequency of natural variations in material properties.

Radial cracking in dynamic indentation experiments.

Continue reading

Research: Thermodynamic Consistency, and Strain-Based Failure

The Kayenta geological material model has been enhanced to span a broader range of pressures and loading rates. Temperature dependence of yield strength has been added along with nonlinear thermoelasticity that can accommodate pressure dependence of the shear modulus and entropy dependence of the bulk modulus in a thermodynamically consistent manner.   Continue reading

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

Merits and shortcomings of conventional smeared damage models

Initial teardrop yield and third-invariant limit surfaces in the Kayenta model

Four classical damage models for concrete (three of which are available in commercial codes) have been compared and critiqued, showing that they all share the notions of a “teardrop” yield surface that can harden and (for some models) translate until reaching a three-invariant fracture limit surface that then collapses to account for softening (i.e., permanent loss of strength).   Practical engineering models for rock and ceramics are similar.  The common drawbacks of these models (primarily severe mesh dependence) can be mitigated, though not eliminated, by seeding their material properties in the simulation with spatial variability (aleatory uncertainty) and by using appropriate scale effects for the strength and failure progression properties. 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

Accelerated hip implant wear testing

3D model Experimental Setup

Rim cracking of polyethylene acetabular liners and squeaking in ceramic components are two important potential failure modes of hip implants, but the loads and stresses that cause such failures are not well understood. Contact stresses in hip implants are analyzed under worst case load conditions to develop new wear testing methods to improve the pre-clinical evaluation of next-generation hip implants and their materials. Complicated full-scale hip implant simulator tests are expensive and take months to complete. A primary goal of this work is to find inexpensive surrogate specimen shapes and loading modes that can, in inexpensive lab tests taking only a few hours, produce the same wear patterns as seen in full-scale prototype testing. Continue reading

Computational approaches for dynamically loaded low-ductility metals

A generic Charpy simulation showing fracture at locations not observed in the lab

Eulerian simulations of un-notched Charpy impact specimens, provide unsatisfactory results in that experimentally observed bend angle, absorbed energy, and fracture mode are not reproduced. The Utah CSM group is independently confirming poor simulation fidelity using conventional constitutive models. From there, we aim to identify the cause, and investigate solutions using capabilities in the Kayenta material framework.

UofU Contributors/collaborators:
Krishna Kamojjala (PhD student, Mech. Engr., UofU)
Scot Swan (MS student, Mech. Engr., UofU)

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