Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Steel grades such as A572 and AISI 4140 are often used for applications where high rate or impact loading may occur. A572 is a hot-rolled carbon steel that is used where a high strength to weight ratio is desired. A grade such as AISI 4140 offers decent corrosion resistance due to higher chromium and molybdenum content and is commonly used in firearm parts, pressurized gas tubes, and structural tubing for roll cages. In these scenarios, the material may undergo high rate loading. Thus, material properties including failure and fracture response at relevant loading rates must be understood so that numerical simulations of impact events accurately capture the deformation and failure/fracture behavior of the involved materials. In this study, the high strain rate tensile response of A572 and 4140 steel are investigated. An increase in yield strength of approximately 28% was observed for 4140 steel when comparing 0.001 s-1 strain rate to 3000 s-1 experiments. A572 showed an increase in yield strength of approximately 52% when the strain rate increased from quasi-static to 2750 s-1. Effects on true stress and strain at failure for the two materials are also discussed.

Abstract not provided.

Abstract not provided.

The exploration of large parameter spaces in search of problem solution and uncertainty quantifcation produces very large ensembles of data. Processing ensemble data will continue to require more resources as simulation complexity and HPC platform throughput increase. More tools are needed to help provide rapid insight into these data sets to decrease manual processing time by the analyst and to increase knowledge the data can provide. One such tool is Tecplot Chorus, whose strengths are visualizing ensemble metadata and linked images. This report contains the analysis and conclusions from evaluating Tecplot Chorus with an example problem that is relevant to Sandia National Laboratories. This report documents a preliminary evaluation of Tecplot Chorus for analyzing ensemble data from CTH simulations. The project that funded this report and evaluation is also evaluating and guiding development with SNL’s Slycat. Slycat and Tecplot Chorus each have their strengths, weaknesses, and overlapping capabilities. It is quite likely that, as the scale of ensemble data increases, both of these tools (and possibly others) will be needed for different processing goals. This report will focus on Tecplot Chorus and its application to an example ensemble of data supplied by David J. Peterson and John P. Korbin; this example is of a flyer plate impact and weld study henceforth referred to as CTH Impact Example. This evaluation also defines a workflow for analysts that can help reduce the time and resources for processing ensemble data.

Abstract not provided.

Blast waves from an explosion in air can cause significant structural damage. As an example, cylindrically-shaped charges have been used for over a century as dynamite sticks for mining, excavation, and demolition. Near the charge, the effects of geometry, standoff from the ground, the proximity to other objects, confinement (tamping), and location of the detonator can significantly affect blast wave characteristics. Furthermore, nonuniformity in the surface characteristics and the density of the charge can affect fireball and shockwave structure. Currently, the best method for predicting the shock structure near a charge and the dynamic loading on nearby structures is to use a multidimensional, multimaterial shock physics code. However, no single numerical technique currently exists for predicting secondary combustion, especially when particulates from the charge are propelled through the fireball and ahead of the leading shock lens. Furthermore, the air within the thin shocked layer can dissociate and ionize. Hence, an appropriate equation of state for air is needed in these extreme environments. As a step towards predicting this complex phenomenon, a technique was developed to provide the equilibrium species composition at every computational cell in an air blast simulation as an initial condition for hand-off to other analysis codes for combustion fluid dynamics or radiation transport. Here, a bare cylindrical charge of TNT detonated in air is simulated using CTH, an Eulerian, finite volume, shock propagation code developed and maintained at Sandia National Laboratories. The shock front propagation is computed at early times, including the detonation wave structure in the explosive and the subsequent air shock up to 100 microseconds, where ambient air entrainment is not significant. At each computational cell, which could have TNT detonation products, air, or both TNT and air, the equilibrium species concentration at the density-energy state is computed using the JCZS2i database in the thermochemical code TIGER. This extensive database of 1267 gas (including 189 ionized species) and 490 condensed species can predict thermodynamic states up to 20,000 K. The results of these calculations provide the detailed three-dimensional structure of a thin shock front, and spatial species concentrations including free radicals and ions. Furthermore, air shock predictions are compared with experimental pressure gage data from a right circular cylinder of pressed TNT, detonated at one end. These complimentary predictions show excellent agreement with the data for the primary wave structure.

The 9/30/2008 ASC Level 2 Post-Processing V&V Milestone (Milestone 2843) contains functionality required by the user community for certain verification and validation tasks. These capabilities include fragment detection from CTH simulation data, fragment characterization and analysis, and fragment sorting and display operations. The capabilities were tested extensively both on sample and actual simulations. In addition, a number of stretch criteria were met including a comparison between simulated and test data, and the ability to output each fragment as an individual geometric file.

Abstract not provided.

In 2004, at the request of the Department of Energy, Sandia National Laboratories (Sandia) prepared a report, ''Guidance on the Risk and Safety Analysis of Large Liquefied Natural Gas (LNG) Spills Over Water''. That report provided framework for assessing hazards and identifying approaches to minimize the consequences to people and property from an LNG spill over water. The report also presented the general scale of possible hazards from a spill from 125,000 m3 o 150,000 m3 class LNG carriers, at the time the most common LNG carrier capacity.

In order to predict blast damage on structures, it is current industry practice to decouple shock calculations from computational structural dynamics calculations. Pressure-time histories from experimental tests were used to assess computational models developed using a shock physics code (CTH) and a structural dynamics code (PRONTO3D). CTH was shown to be able to reproduce three independent characteristics of a blast wave: arrival time, peak overpressure, and decay time. Excellent agreement was achieved for early times, where the rigid wall assumptions used in the model analysis were valid. A one-way coupling was performed for this blast-structure interaction problem by taking the pressure-time history from the shock physics simulation and applying it to the structure at the corresponding locations in the PRONTO3D simulation to capture the structural deformation. In general, the one-way coupling was shown to be a cost-effective means of predicting the structural response when the time duration of the load was less than the response time of the structure. Therefore, the computational models were successfully evaluated for the internal blast problems studied herein.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report provides soil evaluation and characterization testing for the submarine bases at Kings Bay, Georgia, and Bangor, Washington, using triaxial testing at high confining pressures with different moisture contents. In general, the samples from the Bangor and Kings Bay sites appeared to be stronger than a previously used reference soil. Assuming the samples of the material were representative of the material found at the sites, they should be adequate for use in the planned construction. Since soils can vary greatly over even a small site, a soil specification for the construction contractor would be needed to insure that soil variations found at the site would meet or exceed the requirements. A suggested specification for the Bangor and Kings Bay soils was presented based on information gathered from references plus data obtained from this study, which could be used as a basis for design by the construction contractor.

A new capability for modeling thin-shell structures within the coupled Euler-Lagrange code, Zapotec, is under development. The new algorithm creates an artificial material interface for the Eulerian portion of the problem by expanding a Lagrangian shell element such that it has an effective thickness that spans one or more Eulerian cells. The algorithm implementation is discussed along with several examples involving blast loading on plates.

This report outlines the application of finite element methodology to large deformation solid mechanics problems, detailing also some of the key technological issues that effective finite element formulations must address. The presentation is organized into three major portions: first, a discussion of finite element discretization from the global point of view, emphasizing the relationship between a virtual work principle and the associated fully discrete system, second, a discussion of finite element technology, emphasizing the important theoretical and practical features associated with an individual finite element; and third, detailed description of specific elements that enjoy widespread use, providing some examples of the theoretical ideas already described. Descriptions of problem formulation in nonlinear solid mechanics, nonlinear continuum mechanics, and constitutive modeling are given in three companion reports.

This report provides an updated set of users` instructions for PRONTO3D. PRONTO3D is a three-dimensional, transient, solid dynamics code for analyzing large deformations of highly nonlinear materials subjected to extremely high strain rates. This Lagrangian finite element program uses an explicit time integration operator to integrate the equations of motion. Eight-node, uniform strain, hexahedral elements and four-node, quadrilateral, uniform strain shells are used in the finite element formulation. An adaptive time step control algorithm is used to improve stability and performance in plasticity problems. Hourglass distortions can be eliminated without disturbing the finite element solution using either the Flanagan-Belytschko hourglass control scheme or an assumed strain hourglass control scheme. All constitutive models in PRONTO3D are cast in an unrotated configuration defined using the rotation determined from the polar decomposition of the deformation gradient. A robust contact algorithm allows for the impact and interaction of deforming contact surfaces of quite general geometry. The Smooth Particle Hydrodynamics method has been embedded into PRONTO3D using the contact algorithm to couple it with the finite element method.

An efficient, scalable, parallel algorithm for treating contacts in solid mechanics has been applied to interactions between particles in smooth particle hydrodynamics (SPH). The algorithm uses three different decompositions within a single timestep: (1) a static FE-decomposition of mesh elements; (2) a dynamic SPH-decomposition of SPH particles; (3) and a dynamic contact-decomposition of contact nodes and SPH particles. The overhead cost of such a scheme is the cost of moving mesh and particle data between the decompositions. This cost turns out to be small in practice, leading to a highly load-balanced decomposition in which to perform each of the three major computational states within a timestep.

A gridless technique called smooth particle hydrodynamics (SPH) has been coupled with the transient dynamics finite element code pronto. In this paper, a new weighted residual derivation for the SPH method will be presented, and the methods used to embed SPH within pronto will be outlined. Example SPH pronto calculations will also be presented. One major difficulty associated with the Lagrangian finite element method is modeling materials with no shear strength; for example, gases, fluids and explosive biproducts. Typically, these materials can be modeled for only a short time with a Lagrangian finite element code. Large distortions cause tangling of the mesh, which will eventually lead to numerical difficulties, such as negative element area or "bow tie" elements. Remeshing will allow the problem to continue for a short while, but the large distortions can prevent a complete analysis. SPH is a gridless Lagrangian technique. Requiring no mesh, SPH has the potential to model material fracture, large shear flows and penetration. SPH computes the strain rate and the stress divergence based on the nearest neighbors of a particle, which are determined using an efficient particle-sorting technique. Embedding the SPH method within pronto allows part of the problem to be modeled with quadrilateral finite elements, while other parts are modeled with the gridless SPH method. SPH elements are coupled to the quadrilateral elements through a contact-like algorithm. © 1994.

A gridless numerical technique called smooth particle hydrodynamics (SPH) has been coupled the transient dynamics finite element code, PRONTO. In this paper, a new weighted residual derivation for the SPH method will be presented, and the methods used to embed SPH within PRONTO will be outlined. Example SPH-PRONTO calculations will also be presented. Smooth particle hydrodynamics is a gridless Lagrangian technique. Requiring no mesh, SPH has the potential to model material fracture, large shear flows, and penetration. SPH computes the strain rate and the stress divergence based on the nearest neighbors of a particle, which are determined using an efficient particle sorting technique. Embedding the SPH method within PRONTO allows part of the problem to be modeled with quadrilateral finite elements while other parts are modeled with the gridless SPH method. SPH elements are coupled to the quadrilateral elements through a contact like algorithm.

This report will present a brief overview of the transient dynamics capabilities at Sandia National Laboratories, with an emphasis on recent new developments and current research. In addition, the Sandia National Laboratories (SNL) Engineering Analysis Code Access System (SEACAS), which is a collection of structural and thermal codes and utilities used by analysts at SNL, will be described. The SEACAS system includes pre- and post-processing codes, analysis codes, database translation codes, support libraries, Unix shell scripts for execution, and an installation system. SEACAS is used at SNL on a daily basis as a production, research, and development system for the engineering analysts and code developers. Over the past year, approximately 190 days of CPU time have been used by SEACAS codes on jobs running from a few seconds up to two and one-half days of CPU time. SEACAS is running on several different systems at SNL including Cray Unicos, Hewlett Packard HP-UX, Digital Equipment Ultrix, and Sun SunOS. An overview of SEACAS, including a short description of the codes in the system, will be presented. Abstracts and references for the codes are listed at the end of the report.

PRONTO 2D and PRONTO 3D are two- and three-dimensional transient solid dynamics codes for analyzing large deformations of highly nonlinear materials subjected to high strain rates. This newsletter is issued to document changes to these codes. As of this writing, the latest version of PRONTO 2D is Version 4.5.6, and the latest version of PRONTO 3D is Version 4.5.6. This update of the two codes discusses two major modifications to the numerical formulations, three new constitutive models, and the additions and improvements of contact surfaces. Changes in file formats, other miscellaneous revisions, and the availability of PRONTO 2D and PRONTO 3D are also discussed. In addition, updated commands for PRONTO 2D are provided in Appendix A of this newsletter. 29 refs., 12 figs., 2 tabs.