Publications

Abstract not provided.

Recent work at Sandia National Laboratories has focused on preparing strong predictive models for the simulation of ductile failure in metals. The focus of this talk is on the development of engineering-ready models that use a phenomenological approach to represent the ductile fracture processes. As such, an empirical tearing parameter that accounts for mean stress effects along the crack front is presented. A critical value of the tearing parameter is used in finite element calculations as the criterion for crack growth. Regularization is achieved with three different methods and the results are compared. In the first method, upon reaching the critical tearing, the stress within a solid element is decayed by uniformly shrinking the yield surface over a user specified amount of strain. This yields mesh-size dependent results. As a second method for regularization, cohesive surface elements are inserted using an automatic remeshing technique. In the third method, strain-localization elements are inserted with the automated remeshing.

Designing reliable MEMS structures presents numerous challenges. Polycrystalline silicon fractures in a brittle manner with considerable variability in measured strength. Furthermore, it is not clear how to use a measured tensile strength distribution to predict the strength of a complex MEMS structure. To address such issues, two recently developed high throughput MEMS tensile test techniques have been used to measure strength distribution tails. The measured tensile strength distributions enable the definition of a threshold strength as well as an inferred maximum flaw size. The nature of strength-controlling flaws has been identified and sources of the observed variation in strength investigated. A double edge-notched specimen geometry was also tested to study the effect of a severe, micron-scale stress concentration on the measured strength distribution. Strength-based, Weibull-based, and fracture mechanics-based failure analyses were performed and compared with the experimental results.

This report describes a Laboratory Directed Research and Development (LDRD) project to use of synchrotron-radiation computed tomography (SRCT) data to determine the conditions and mechanisms that lead to void nucleation in rolled alloys. The Advanced Light Source (ALS) at Lawrence Berkeley National Laboratory (LBNL) has provided SRCT data of a few specimens of 7075-T7351 aluminum plate (widely used for aerospace applications) stretched to failure, loaded in directions perpendicular and parallel to the rolling direction. The resolution of SRCT data is 900nm, which allows elucidation of the mechanisms governing void growth and coalescence. This resolution is not fine enough, however, for nucleation. We propose the use statistics and image processing techniques to obtain sub-resolution scale information from these data, and thus determine where in the specimen and when during the loading program nucleation occurs and the mechanisms that lead to it. Quantitative analysis of the tomography data, however, leads to the conclusion that the reconstruction process compromises the information obtained from the scans. Alternate, more powerful reconstruction algorithms are needed to address this problem, but those fall beyond the scope of this project.

Threshold stress intensity factors were measured in high-pressure hydrogen gas for a variety of low alloy ferritic steels using both constant crack opening displacement and rising crack opening displacement procedures. The sustained load cracking procedures are generally consistent with those in ASME Article KD-10 of Section VIII Division 3 of the Boiler and Pressure Vessel Code, which was recently published to guide design of high-pressure hydrogen vessels. Three definitions of threshold were established for the two test methods: K{sub THi}* is the maximum applied stress intensity factor for which no crack extension was observed under constant displacement; K{sub THa} is the stress intensity factor at the arrest position for a crack that extended under constant displacement; and K{sub JH} is the stress intensity factor at the onset of crack extension under rising displacement. The apparent crack initiation threshold under constant displacement, K{sub THi}*, and the crack arrest threshold, K{sub THa}, were both found to be non-conservative due to the hydrogen exposure and crack-tip deformation histories associated with typical procedures for sustained-load cracking tests under constant displacement. In contrast, K{sub JH}, which is measured under concurrent rising displacement and hydrogen gas exposure, provides a more conservative hydrogen-assisted fracture threshold that is relevant to structural components in which sub-critical crack extension is driven by internal hydrogen gas pressure.

Fatigue cracking in metals has been and is an area of great importance to the science and technology of structural materials for quite some time. The earliest stages of fatigue crack nucleation and growth are dominated by the microstructure and yet few models are able to predict the fatigue behavior during these stages because of a lack of microstructural physics in the models. This program has developed several new simulation tools to increase the microstructural physics available for fatigue prediction. In addition, this program has extended and developed microscale experimental methods to allow the validation of new microstructural models for deformation in metals. We have applied these developments to fatigue experiments in metals where the microstructure has been intentionally varied.

The mechanical properties of some materials (Cu, Ni, Ag, etc.) have been shown to develop strong dependence on the geometric dimensions, resulting in a size effect. Several theories have been proposed to model size effects, but have been based on very few experiments conducted at appropriate scales. Some experimental results implied that size effects are caused by increasing strain gradients and have been used to confirm many strain gradient theories. On the other hand, some recent experiments show that a size effect exists in the absence of strain gradients. This report describes a brief analytical and experimental study trying to clarify the material and experimental issues surrounding the most influential size-effect experiments by Fleck et al (1994). This effort is to understand size effects intended to further develop predictive models.

Abstract not provided.

Abstract not provided.

A micromechanical model is developed for grain bridging in monolithic ceramics. Specifically, bridge formation of a single, non-equiaxed grain spanning adjacent grains is addressed. A cohesive zone framework enables crack initiation and propagation along grain boundaries. The evolution of the bridge is investigated through a variance in both grain angle and aspect ratio. We propose that the bridging process can be partitioned into five distinct regimes of resistance: propagate, kink, arrest, stall, and bridge. Although crack propagation and kinking are well understood, crack arrest and subsequent "stall" have been largely overlooked. Resistance during the stall regime exposes large volumes of microstructure to stresses well in excess of the grain boundary strength. Bridging can occur through continued propagation or reinitiation ahead of the stalled crack tip. The driving force required to reinitiate is substantially greater than the driving force required to kink. In addition, the critical driving force to reinitiate is sensitive to grain aspect ratio but relatively insensitive to grain angle. The marked increase in crack resistance occurs prior to bridge formation and provides an interpretation for the rapidly rising resistance curves which govern the strength of many brittle materials at realistically small flaw sizes. © 2006 Elsevier Ltd. All rights reserved.

Polymers and fiber-reinforced polymer matrix composites play an important role in many Defense Program applications. Recently an advanced nonlinear viscoelastic model for polymers has been developed and incorporated into ADAGIO, Sandia's SIERRA-based quasi-static analysis code. Standard linear elastic shell and continuum models for fiber-reinforced polymer-matrix composites have also been added to ADAGIO. This report details the use of these models for advanced adhesive joint and composites simulations carried out as part of an Advanced Simulation and Computing Advanced Deployment (ASC AD) project. More specifically, the thermo-mechanical response of an adhesive joint when loaded during repeated thermal cycling is simulated, the response of some composite rings under internal pressurization is calculated, and the performance of a composite container subjected to internal pressurization, thermal loading, and distributed mechanical loading is determined. Finally, general comparisons between the continuum and shell element approaches for modeling composites using ADAGIO are given.

This report documents activities related to the ASCI AD Resistance Weld Process Modeling Project AD2003-15. Activities up to and including FY2004 are discussed. This was the third year for this multi year project, the objective of which is to position the SIERRA computational tools for the solution of resistance welding problems. The process of interest is a three-way coupled problem involving current flow, temperature buildup and large plastic deformation. The DSW application is the reclamation stem weld used in the manufacture of high pressure gas bottles. This is the first year the CALAGIO suite of codes (eCALORE, CALORE, and ADAGIO) was used to successfully solve a three-way coupled problem in SIERRA. This report discusses the application of CALAGIO to the tapered bar acceptance problem and a similar but independent tapered bar simulation of a companion C6 experiment. New additions to the EMMI constitutive model and issues related to CALAGIO performance are also discussed.