Publications

Abstract not provided.

Abstract not provided.

We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The matters of interest here are the elastic-plastic response and ductile failure behavior of 6061-T651 aluminum alloy under the conditions described above. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of a modular $J_2$ plasticity model with isotropic hardening as well as a ductile failure model. Both included modules that accounted for temperature and strain rate dependence. The models were coupled with an adiabatic heating module to calculate the temperature rise due to the conversion of plastic work to heat. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10^–4 and 10³ 1/s as well as four geometries of notched tension specimens and two tests on specimens with shear-dominated deformations. The test data collected allowed the calibration of both the plasticity and the ductile failure models. Most test specimens were extracted from a single piece of plate to maintain consistency. Notched tension tests came from a possibly different plate, but from the same lot. When using the model in structural finite element calculations, element formulations and sizes different from those used to model the test specimens in the calibration are likely to be used. A brief investigation demonstrated that the failure model can be particularly sensitive to the element selection and provided an initial guide to compensate in a specific example.

This is the second part of a two-part contribution on modeling of the anisotropic elastic-plastic response of aluminum 7079 from an extruded tube. Part I focused on calibrating a suite of yield and hardening functions from tension test data; Part II concentrates on evaluating those calibrations. Here, a rectangular validation specimen with a blind hole was designed to provide heterogeneous strain fields that exercise the material anisotropy, while at the same time avoiding strain concentrations near sample edges where Digital Image Correlation (DIC) measurements are difficult to make. Specimens were extracted from the tube in four different orientations and tested in tension with stereo-DIC measurements on both sides of the specimen. Corresponding Finite Element Analysis (FEA) with calibrated isotropic (von Mises) and anisotropic (Yld2004-18p) yield functions were also conducted, and both global force-extension curves as well as full-field strains were compared between the experiments and simulations. Specifically, quantitative full-field strain error maps were computed using the DIC-leveling approach proposed by Lava et al. The specimens experienced small deviations from ideal boundary conditions in the experiments, which had a first-order effect on the results. Therefore, the actual experimental boundary conditions had to be applied to the FEA in order to make valid comparisons. The predicted global force-extension curves agreed well with the measurements overall, but were sensitive to the boundary conditions in the nonlinear regime and could not differentiate between the two yield functions. Interrogation of the strain fields both qualitatively and quantitatively showed that the Yld2004-18p model was clearly able to better describe the strain fields on the surface of the specimen compared to the von Mises model. These results justify the increased complexity of the calibration process required for the Yld2004-18p model in applications where capturing the strain field evolution accurately is important, but not if only the global force-extension response of the elastic–plastic region is of interest.

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The material of interest here is 304L stainless steel, and the objective of this work is to calibrate thermal-mechanical models: one for the constitutive behavior and another for failure. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of the models. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10^–4 and 10³ 1/s. It also included notched tension and shear-dominated compression hat tests specifically designed to calibrate the failure model. All test specimens were extracted from a single piece of plate to maintain consistency. The constitutive model adopted was a modular $J_2$ plasticity model with isotropic hardening that included rate and temperature dependence. A criterion for failure initiation based on a critical value of equivalent plastic strain fitted the failure data appropriately and was adopted. Possible ranges of the values of the parameters of the models were determined partially on historical data from calibrations of the same alloy from other lots and are given here. The calibration of the parameters of the models were based on finite element simulations of the various material tests using relatively ne meshes and hexahedral elements. When using the model in structural finite element calculations, however, element formulations and sizes different from those in the calibration are likely to be used. A brief investigation demonstrated that the failure initiation predictions can be particularly sensitive to the element selection and provided an initial guide to compensate for the effect of element size in a specific example.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this work, we investigated microstructural features of elastomeric foam with the goal of identifying descriptors other than porosity that have a significant effect on the macroscale mechanical response. X-ray computed tomography (XCT) provided three-dimensional images of several flexible polyurethane foam samples prior to mechanical testing. The samples were then compressed to approximately 80% engineering strain. Stereo digital image correlation was used to measure the three-dimensional surface displacement data, from which strain was determined. The strain data, which were calculated with respect to the undeformed coordinates, were then overlaid on the corresponding surface generated from XCT. Heterogeneities in the strain-field were cross-correlated with topological quantities such as pore size distribution. A statistically significant correlation was identified between the distance transform of the pore phase and strain fluctuations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Sandia Fracture Challenges provide the mechanics community a forum for assessing its ability to predict ductile fracture through a blind, round-robin format where computationalists are asked to predict the deformation and failure of an arbitrary geometry given experimental calibration data. This presentation will cover the three Sandia Fracture Challenges, with emphasis on the third. The third Challenge, issued in 2017, consisted of an additively manufactured 316L stainless steel tensile bar with through holes and internal cavities that could not have been conventionally machined. The volunteer prediction teams were provided extensive materials data from tensile tests of specimens printed on the same build tray to electron backscatter diffraction microstructural maps and micro-computed tomography scans of the Challenge geometry. The teams were asked a variety of questions, including predictions of variability in the resulting fracture response, as the basis for assessment of their predictive capabilities. This presentation will describe the Challenges and compare the experimental results to the predictions, identifying gaps in capabilities, both experimentally and computationally, to inform future investments. The Sandia Fracture Challenge has evolved into the Structural Reliability Partnership, where researchers will create several blind challenges covering a wider variety of topics in structural reliability. This presentation will also describe this new venture.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This project targeted a full-field understanding of the conversion of plastic work into heat us- ing advanced diagnostics (digital image correlation, DIC, combined with infrared, IR, imaging). This understanding will act as a catalyst for reformulating the prevalent simplistic model, which will ultimately transform Sandia's ability to design for and predict thermomechanical behavior, impacting national security applications including nuclear weapon assessments of accident scenar- ios. Tensile 304L stainless steel dogbones are pulled in tension at quasi-static rates until failure and full-field deformation and temperature data are captured, while accounting for thermal losses. The IR temperature fields are mapped onto the DIC coordinate system (Lagrangian formulation). The resultant fields are used to calculate the Taylor-Quinney coefficient, p, at two strain rates rates (0.002 s -1 and 0.08 s -1 ) and two temperatures (room temperature, RT, and 250degC).

Modeling material and component behavior using finite element analysis (FEA) is critical for modern engineering. One key to a credible model is having an accurate material model, with calibrated model parameters, which describes the constitutive relationship between the deformation and the resulting stress in the material. As such, identifying material model parameters is critical to accurate and predictive FEA. Traditional calibration approaches use only global data (e.g. extensometers and resultant force) and simplified geometries to find the parameters. However, the utilization of rapidly maturing full-field characterization tech- niques (e.g. Digital Image Correlation (DIC)) with inverse techniques (e.g. the Virtual Feilds Method (VFM)) provide a new, novel and improved method for parameter identification. This LDRD tested that idea: in particular, whether more parameters could be identified per test when using full-field data. The research described in this report successfully proves this hypothesis by comparing the VFM results with traditional calibration methods. Important products of the research include: verified VFM codes for identifying model parameters, a new look at parameter covariance in material model parameter estimation, new validation tech- niques to better utilize full-field measurements, and an exploration of optimized specimen design for improved data richness.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Partial penetration laser welds join metal surfaces without additional filler material, providing hermetic seals for a variety of components. The crack-like geometry of a partial penetration weld is a local stress riser that may lead to failure of the component in the weld. Computational modeling of laser welds has shown that the model should include damage evolution to predict the large deformation and failure. We have performed interrupted tensile experiments both to characterize the damage evolution and failure in laser welds and to aid computational modeling of these welds. Several EDM-notched and laser-welded 304L stainless steel tensile coupons were pulled in tension, each one to a different load level, and then sectioned and imaged to show the evolution of damage in the laser weld and in the EDM-notched parent 304L material (having a similar geometry to the partial penetration laser-welded material). SEM imaging of these specimens revealed considerable cracking at the root of the laser welds and some visible micro-cracking in the root of the EDM notch even before peak load was achieved in these specimens. The images also showed deformation-induced damage in the root of the notch and laser weld prior to the appearance of the main crack, though the laser-welded specimens tended to have more extensive damage than the notched material. These experiments show that the local geometry alone is not the cause of the damage, but also microstructure of the laser weld, which requires additional investigation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A “good” speckle pattern enables DIC to make its full-field measurements, but oftentimes this artistic part of the DIC setup takes a considerable amount of time to develop and evaluate for a given optical configuration. A catalog of well-quantified speckle patterns for various fields of view would greatly decrease the time it would take to start making DIC measurements. The purpose of this speckle patterning study is to evaluate various speckling techniques we had readily available in our laboratories for fields of view from around 100 mm down to 5 mm that are common for laboratory-scale experiments. The list of speckling techniques is not exhaustive: spray painting, UV-printing of computer-designed speckle patterns, airbrushing, and particle dispersion. First, we quantified the resolution of our optical configurations for each of the fields of view to determine the smallest speckle we could resolve. Second, we imaged several speckle patterns at each field of view. Third, we quantified the average and standard deviation of the speckle size, speckle contrast, and density to characterize the quality of the speckle pattern. Finally, we performed computer-aided sub-pixel translation of the speckle patterns and ran correlations to examine how well DIC tracked the pattern translations. We discuss our metrics for a “good” speckle pattern and outline how others may perform similar studies for their desired optical configurations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Evermore sophisticated ductile plasticity and failure models demand experimental material characterization of shear behavior; yet, the mechanics community lacks a widely accepted, standard test method for shear-dominated deformation and failure of ductile metals. We investigated the use of the V-notched rail test, borrowed from the ASTM D7078 standard for shear testing of composites, for shear testing of Ti-6Al-4V titanium alloy sheet material, considering sheet rolling direction and quasi-static and transient load rates. In this paper, we discuss practical aspects of testing, modifications to the specimen geometry, and the experimental shear behavior of Ti-6Al-4V. Specimen installation, machine compliance, specimen-grip slip during testing, and specimen V-notched geometry all influenced the measured specimen behavior such that repeatable shear-dominated behavior was initially difficult to obtain. We will discuss the careful experimental procedure and set of measurements necessary to extract meaningful shear information for Ti-6Al-4V. We also evaluate the merits and deficiencies, including practicality of testing for engineering applications and quality of results, of the V-notched rail test for characterization of ductile shear behavior.

Abstract not provided.

Abstract not provided.

The 2014 WSEAT X-Prize is modeled as a double blind study to challenge the computational and material mechanics communities methodologies to develop better capabilities in modeling and experimentation to predict the failure in ductile metals. The challenge is presented as a distinct, yet relatively, simple geometry with all reported modeling predictions blind to each of the modeling teams. The experimental testing is validated by two independent test labs to confirm the experimentally observed behavior and results are unbiased and repeatable. The WSEAT X-Prize was issued to both external participants and internal participants as the Sandia Fracture Challenge 2 (SFC2) on May 30, 2014. A Challenge Supplemental Information Packet was sent to participants on August 13, 2014 to Prior years SFCs focused on the ability to predict failures under a quasi-static loading condition that focused on either a shear or tensile-dominated failure mode. This year’s challenge focuses on a geometry with a shear and/or tensile-dominated failure mode influenced by a moderate strain-rate ductile fracture in a metallic alloy.

The Virtual Fields Method (VFM) is an inverse method for constitutive model parameter identication that relies on full-eld experimental measurements of displacements. VFM is an alternative to standard approaches that require several experiments of simple geometries to calibrate a constitutive model. VFM is one of several techniques that use full-eld exper- imental data, including Finite Element Method Updating (FEMU) techniques, but VFM is computationally fast, not requiring iterative FEM analyses. This report describes the im- plementation and evaluation of VFM primarily for nite-deformation plasticity constitutive models. VFM was successfully implemented in MATLAB and evaluated using simulated FEM data that included representative experimental noise found in the Digital Image Cor- relation (DIC) optical technique that provides full-eld displacement measurements. VFM was able to identify constitutive model parameters for the BCJ plasticity model even in the presence of simulated DIC noise, demonstrating VFM as a viable alternative inverse method. Further research is required before VFM can be adopted as a standard method for constitu- tive model parameter identication, but this study is a foundation for ongoing research at Sandia for improving constitutive model calibration.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.