R.M. Brannon, A.F. Fossum, and O.E. Strack
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.
J.M. Wells and R.M. Brannon
With the relatively recent introduction of quantitative and volumetric X-ray computedtomography (XCT) applied to ballistic impact damage diagnostics, signiﬁcant 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 inﬂuence that impact damage may play in the mitigation orprevention of ballistic penetration. Examples of current capabilities of XCT characterization,quantiﬁcation, and visualization of complex impact damage variants are demonstrated anddiscussed for impacted ceramic and metallic terminal ballistic target materials. Potential beneﬁtsof incorporating such impact damage diagnostics in future ballistic computational modeling arealso brieﬂy discussed.
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
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.
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
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.
Michael Braginski (postdoc, Mech. Engr., UofU)
Jeff Burghardt (PhD student, Mech. Engr., UofU)
Stephen Bauer (Manager, Sandia National Labs geomechanics testing lab)
David Bronowski (Sandia geomechanics lab technician)
Erik Strack (Manager, Sandia Labs Computational Physics)
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.
Krishna Kamojjala (PhD student, Mech. Engr., UofU)
Scot Swan (MS student, Mech. Engr., UofU)