Publication: Validating Theories for Brittle Damage

R.M. Brannon, J.M. Wells, and O.E. Strack

Realistic-looking, uneven damage zones in Brazilian simulations compare favorably with laboratory data for observable damage

Validating simulated predictions of internal damage within armor ceramics is preferable to simply assessing a models ability to predict penetration depth, especially if one hopes to perform subsequent ‘‘second strike’’ analyses. We present the results of a study in which crack networks are seeded by using a statistically perturbed strength, the median of which is inherited from a deterministic ‘‘smeared damage’’ model, with adjustments to reflect experimentally established size effects. This minor alteration of an otherwise conventional damage model noticeably mitigates mesh dependencies and, at virtually no computational cost, produces far more realistic cracking patterns that are well suited for validation against X-ray computed tomography (XCT) images of internal damage patterns. For Brazilian, spall, and indentation tests, simulations share qualitative features with externally visible damage. However, the need for more stringent quantitative validation, software quality testing, and subsurface XCT validation, is emphasized.

Continue reading

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

A. F. Fossum and R. M. Brannon

Rate Dependance

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:


Publication: Initial inclusion of thermodynamic considerations in Kayenta

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 [1] (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.

Available Online:

Publication: Application of Uintah-MPM to shaped charge jet penetration of aluminum

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.

Available Online:

Publication: The Use of Sphere Indentation Experiments to Characterize Ceramic Damage Models

R.B. Leavy; R.M. Brannon; O.E. Strack

Weibull modulus effect on radial cracking in boron carbide simulations impacted at 400 m/s.

Sphere impact experiments are used to calibrate and validate ceramic models that include statistical variability and/or scale effects in strength and toughness parameters. These dynamic experiments supplement traditional characterization experiments such as tension, triaxial compression, Brazilian, and plate impact, which are commonly used for ceramic model calibration.The fractured ceramic specimens are analyzed using sectioning, X-ray computed tomography, microscopy, and other techniques. These experimental observations indicate that a predictive material model must incorporate a standard deviation in strength that varies with the nature of the loading. Methods of using the spherical indentation data to calibrate a statistical damage model are presented in which it is assumed that variability in strength is tied to microscale stress concentrations associated with microscale heterogeneity.

Available Online:

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: