Material Testing 2.0 (MT2.0) is a paradigm that advocates for the use of rich, full-field data, such as from digital image correlation and infrared thermography, for material identification. By employing heterogeneous, multi-axial data in conjunction with sophisticated inverse calibration techniques such as finite element model updating and the virtual fields method, MT2.0 aims to reduce the number of specimens needed for material identification and to increase confidence in the calibration results. To support continued development, improvement, and validation of such inverse methods—specifically for rate-dependent, temperature-dependent, and anisotropic metal plasticity models—we provide here a thorough experimental data set for 304L stainless steel sheet metal. The data set includes full-field displacement, strain, and temperature data for seven unique specimen geometries tested at different strain rates and in different material orientations. Commensurate extensometer strain data from tensile dog bones is provided as well for comparison. We believe this complete data set will be a valuable contribution to the experimental and computational mechanics communities, supporting continued advances in material identification methods.
Additive manufacturing (AM) maintains a wide process window that enables complex designs otherwise unattainable via conventional production technologies. However, the lack of confidence in qualifying AM parts that leverage AM process–structure–property–performance (PSPP) relationships stymies design optimization and adoption of AM. While continuing efforts to map fundamental PSPP relationships that cover the potential design space, we first need pragmatic and then long-term solutions that overcome challenges associated with qualifying AM-designed parts. Two pragmatic solutions include: (1) AM material specifications to substantiate process reproducibility, and (2) component risk categorization to associate system risk relative to part performance and required part quality. A novel qualification paradigm under development involves efficient prediction of part performance over wide-ranging PSPP relationships through targeted testing and computational simulation. This paper describes projects at Sandia National Laboratories on PSPP relationship discovery, these pragmatic approaches, and the novel qualification approach.
Traditional electronics assemblies are typically packaged using physically or chemically blown potted foams to reduce the effects of shock and vibration. These potting materials have several drawbacks including manufacturing reliability, lack of internal preload control, and poor serviceability. A modular foam encapsulation approach combined with additively manufactured (AM) silicone lattice compression structures can address these issues for packaged electronics. These preloaded silicone lattice structures, known as foam replacement structures (FRSs), are an integral part of the encapsulation approach and must be properly characterized to model the assembly stresses and dynamics. In this study, dynamic test data is used to validate finite element models of an electronics assembly with modular encapsulation and a direct ink write (DIW) AM silicone FRS. A variety of DIW compression architectures are characterized, and their nominal stress-strain behavior is represented with hyperfoam constitutive model parameterizations. Modeling is conducted with Sierra finite element software, specifically with a handoff from assembly preloading and uniaxial compression in Sierra/Solid Mechanics to linear modal and vibration analysis in Sierra/Structural Dynamics. This work demonstrates the application of this advanced modeling workflow, and results show good agreement with test data for both static and dynamic quantities of interest, including preload, modal, and vibration response.
The Sandia Mechanics Challenge (SMC) provides the solid-mechanics community a forum for assessing its ability to predict mechanical behavior in structures and materials through a blind, round-robin format. Computationalists are asked to predict the behavior of an unfamiliar geometry given experimental calibration data, their predictions are compared to experimental measurements of the SMC scenario, and then the participants assess and compare their approaches, documenting their findings. The SMC broadens the scope of Sandia-hosted benchmarking problems that previously focused on ductile failure through the Sandia Fracture Challenges, enabling an enduring, community-wide self-assessment of predictive capabilities for a variety of mechanics topics. The SMC is part of the Structural Reliability Partnership, which offers other benchmarking challenges hosted by several participating institutions.
Here this paper introduces a publicly available PyTorch-ABAQUS deep-learning framework of a family of plasticity models where the yield surface is implicitly represented by a scalar-valued function. In particular, our focus is to introduce a practical framework that can be deployed for engineering analysis that employs a user-defined material subroutine (UMAT/VUMAT) for ABAQUS, which is written in FORTRAN. To accomplish this task while leveraging the back-propagation learning algorithm to speed up the neural-network training, we introduce an interface code where the weights and biases of the trained neural networks obtained via the PyTorch library can be automatically converted into a generic FORTRAN code that can be a part of the UMAT/VUMAT algorithm. To enable third-party validation, we purposely make all the data sets, source code used to train the neural-network-based constitutive models, and the trained models available in a public repository. Furthermore, the practicality of the workflow is then further tested on a dataset for anisotropic yield function to showcase the extensibility of the proposed framework. A number of representative numerical experiments are used to examine the accuracy, robustness and reproducibility of the results generated by the neural network models.
The mechanical behavior of partial-penetration laser welds exhibits significant variability in engineering quantities such as strength and apparent ductility. Understanding the root cause of this variability is important when using such welds in engineering designs. In Part II of this work, we develop finite element simulations with geometry derived from micro-computed tomography (μCT) scans of partial-penetration 304L stainless steel laser welds that were analyzed in Part I. We use these models to study the effects of the welds’ small-scale geometry, including porosity and weld depth variability, on the structural performance metrics of weld ductility and strength under quasi-static tensile loading. We show that this small-scale geometry is the primary cause of the observed variability for these mechanical response quantities. Additionally, we explore the sensitivity of model results to the conversion of the μCT data to discretized model geometry using different segmentation algorithms, and to the effect of small-scale geometry simplifications for pore shape and weld root texture. The modeling approach outlined and results of this work may be applicable to other material systems with small-scale geometric features and defects, such as additively manufactured materials.
The calibration of solid constitutive models with full-field experimental data is a long-standing challenge, especially in materials that undergo large deformations. In this paper, we propose a physics-informed deep-learning framework for the discovery of hyperelastic constitutive model parameterizations given full-field surface displacement data and global force-displacement data. Contrary to the majority of recent literature in this field, we work with the weak form of the governing equations rather than the strong form to impose physical constraints upon the neural network predictions. The approach presented in this paper is computationally efficient, suitable for irregular geometric domains, and readily ingests displacement data without the need for interpolation onto a computational grid. A selection of canonical hyperelastic material models suitable for different material classes is considered including the Neo–Hookean, Gent, and Blatz–Ko constitutive models as exemplars for general non-linear elastic behaviour, elastomer behaviour with finite strain lock-up, and compressible foam behaviour, respectively. We demonstrate that physics informed machine learning is an enabling technology and may shift the paradigm of how full-field experimental data are utilized to calibrate constitutive models under finite deformations.
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.
The choice of model form used to represent the anisotropic yield response of metals can depend strongly on the type and amount of data available for calibration. This two-part contribution considers the calibration (part I) of three yield functions: von Mises, Hill-48 and Yld2004-18p by Barlat and co-workers. This is followed by model verification exercises (part II). The material used was a 7079 aluminum alloy extruded tube. The calibration data were measurements of yield stress and Lankford ratio from uniaxial tension specimens cut along 12 orientations. Given that the tube was relatively thick-walled, some of the orientations included through-thickness components. This allowed the calibrations to be based exclusively on test data, without the need for parameter assumptions or supplemental crystal plasticity calculations. The Yld2004-18p function provided the best fit to the data available due to its 18 anisotropy parameters plus an unspecified exponent, compared to the quadratic Hill function with 6 anisotropy parameters and to the isotropic von Mises function. Whereas the Yld2004-18p function did not warrant further exploration due to the excellent fit it provided, the results showed that care must be taken when using Hill’s function. Finally, due to its parametrization with only 6 anisotropy parameters, it can significantly misrepresent the yield behavior depending on the calibration data used, possibly rendering it less desirable than a simple isotropic function in some applications.
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 103 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.
Background: Using a thin-walled tube torsion test to characterize a material’s shear response is a well-known technique; however, the thin walled specimen tends to buckle before reaching large shear deformation and failure. An alternative technique is the surface stress method (Nadai 1950; Wu et al. J Test Eval 20:396–402, 1992), which derives a shear stress-strain curve from the torque-angular displacement relationship of a solid cylindrical bar. The solid bar torsion test uniquely stabilizes the deformation which allows us to control and explore very large shear deformation up to failure. However, this method has rarely been considered in the literature, possibly due to the complexity of the analysis and experimental issues such as twist measurement and specimen uniformity. Objective: In this investigation, we develop a method to measure the large angular displacement in the solid bar torsion experiments to study the large shear deformation of two common engineering materials, Al6061-T6 and SS304L, which have distinctive hardening behaviors. Methods: Modern stereo-DIC methods were applied to make deformation measurements. The large angular displacement of the specimen posed challenges for the DIC analysis. An analysis method using multiple reference configurations and transformation of deformation gradient is developed to make the large shear deformation measurement successful. Results: We successfully applied the solid bar torsion experiment and the new analysis method to measure the large shear deformation of Al6061-T6 and SS304L till specimen failure. The engineering shear strains at failure are on the order of 2–3 for Al6061-T6 and 3–4 for SS304L. Shear stress-strain curves of Al6061-T6 and SS304L are also obtained. Conclusions: Solid bar torsion experiments coupled with 3D-DIC technique and the new analysis method of deformation gradient transformation enable measurement of very large shear deformation up to specimen failure.
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 103 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.
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.
This report describes the efforts to characterize and model General Plastics TF6070 and EF4000 flexible polyurethane foams under room temperature, large deformation quasi-static cyclic mechanical loading conditions. Densities from three to fifteen pound per cubic foot (PCF) are examined, which correspond to relative densities of approximately 4 to 20%. These foams are open cell structured and flexible at room temperature with a glass transition transition less than -30°C, and they fully recover their original shape when unloaded (at room temperature). Uniaxial compression tests were conducted with accompanying lateral image series for Digital Image Correlation (DIC) analysis with the goal of extracting transverse strain responses. Due to difficulties with DIC analysis at large strains, lateral strains were instead extracted for each test via edge tracking. The experimental results exhibit a nonlinear elastic response and anisotropic material behavior (particularly for the lower densities). Some hysteresis is observed that is different between the first and subsequent cycles of deformation indicating both a small degree of permanent damage (reduced stiffness during reloading) and viscoelasticity. These inelastic mechanisms are not considered in the modeling and calibration in this report. This work considers only the rate independent, room temperature foam behavior. Individual foam densities were calibrated for loading in two directions, parallel and perpendicular to the foam bubble rise direction, since the mechanical behavior is different in these two directions. The Flex Foam constitutive model was used for all parameterizations despite the fact that the model is isotropic. A review of the constitutive model is given as well as necessary data reduction procedures to parameterize it for each foam density and orientation are discussed. Finally, two different parameterizations are developed that take the undeformed foam density as an input that span all densities considered. These two parameterized models represent foams loaded in the rise and transverse directions respectively. We summarize the assumptions and limitations of the parameterizations provided in this report to guide analysis choices with them. All parameterizations presented herein have the following traits, room temperature, rate independent, damage-free, and non-dissipative . Isotropy (even if they are representing anisotropic data). Supplied Sierra Solid Mechanics Flex Foam Model Inputs are in units: pounds, inches, Celsius, and seconds
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.
The classic models for ductile fracture of metals were based on experimental observations dating back to the 1950’s. Using advanced microscopy techniques and modeling algorithms that have been developed over the past several decades, it is possible now to examine the micro- and nano-scale mechanisms of ductile rupture in more detail. This new information enables a revised understanding of the ductile rupture process under quasi-static room temperature conditions in ductile pure metals and alloys containing hard particles. While ductile rupture has traditionally been viewed through the lens of nucleation-growth-and-coalescence, a new taxonomy is proposed involving the competition or cooperation of up to seven distinct rupture mechanisms. Generally, void nucleation via vacancy condensation is not rate limiting, but is extensive within localized shear bands of intense deformation. Instead, the controlling process appears to be the development of intense local dislocation activity which enables void growth via dislocation absorption.
An initial foray into the design of specimens that can be used to provide data about the quasistatic ductile failure of metals when subjected to shear-dominated (low triaxiality) states of stress was undertaken. Four specimen geometries made from two materials with different ductility (Al 7075, lower ductility and steel A286, higher ductility) were considered as candidates. Based on results from analysis and experimentation, it seems that two show promise for further consideration. Whereas preliminary results indicate that the Johnson-Cook model fit the failure data for Al 7075 well, it did not fit the data for steel A286. Further work is needed to consolidate the results and evaluate other failure models that may fit the steel data better, as well as to extend the results of this work to the dynamic loading regime.
The third Sandia Fracture Challenge (SFC3) was a benchmark problem for comparing experimental and simulated ductile deformation and failure in an additively manufactured (AM) 316L stainless steel structure. One surprising observation from the SFC3 was the Challenge-geometry specimens had low variability in global load versus displacement behavior, attributed to the large stress-concentrating geometric features dominating the global behavior, rather than the AM voids that tend to significantly influence geometries with uniform cross-sections. This current study reinvestigates the damage and failure evolution of the Challenge-geometry specimens, utilizing interrupted tensile testing with micro-computed tomography (micro-CT) scans to monitor AM void and crack growth from a virgin state through complete failure. This study did not find a correlation between global load versus displacement behavior and AM void attributes, such as void volume, location, quantity, and relative size, which incidentally corroborates the observation from the SFC3. However, this study does show that the voids affect the local behavior of damage and failure. Surface defects (i.e. large voids located on the surface, far exceeding the nominal surface roughness) that were near the primary stress concentration affected the location of crack initiation in some cases, but they did not noticeably affect the global response. The fracture surfaces were a combination of classic ductile dimples and crack deviation from a more direct path favoring intersection with AM voids. Even though the AM voids promoted crack deviation, pre-test micro-CT scan statistics of the voids did not allow for conclusive predictions of preferred crack paths. This study is a first step towards investigating the importance of voids on the ductile failure of AM structures with stress concentrations.
The mounting reliance on computational simulations to predict all aspects of the lifecycle of a mechanical system, from fabrication to failure, has prompted the mechanics community to selfassess its abilities to perform those predictions. Benchmark problems in mechanics that compare simulations that use different computational approaches with experiments have sprung up lately, including the NIST AM-Bench looking at additively manufactured (AM) materials (https://www.nist.gov/ambench),the Contact-Mechanics Challenge (Miiser, 2017) considering adhesion between two nominally flat surfaces, Numisheet providing semiannual benchmarking activities in sheet metal forming (http://numisheet2018.org),and the Sandia Fracture Challenge (SFC) (Boyce, 2014 and Boyce, 2016) investigating ductile failure. The previous SFCs have shown that progress has been made in computations of ductile failure, but improvements still can be made, hence the third Sandia Fracture Challenge (SFC3), the subject of this Special Volume. The most recent installment of SFC is building on previous successes and tackling the difficult problem of fracture in an AM 316L stainless steel structure.
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.
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 techniques (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 techniques to better utilize full-field measurements, and an exploration of optimized specimen design for improved data richness.
This project targeted a full-field understanding of the conversion of plastic work into heat using 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 scenarios. 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, β, at two strain rates rates (0.002 s-1 and 0.08 s-1) and two temperatures (room temperature, RT, and 250°C).
Traditionally, material identification is performed using global load and displacement data from simple boundary-value problems such as uni-axial tensile and simple shear tests. More recently, however, inverse techniques such as the Virtual Fields Method (VFM) that capitalize on heterogeneous, full-field deformation data have gained popularity. In this work, we have written a VFM code in a finite-deformation framework for calibration of a viscoplastic (i.e. strain-rate dependent) material model for 304L stainless steel. Using simulated experimental data generated via finite-element analysis (FEA), we verified our VFM code and compared the identified parameters with the reference parameters input into the FEA. The identified material model parameters had surprisingly large error compared to the reference parameters, which was traced to parameter covariance and the existence of many essentially equivalent parameter sets. This parameter non-uniqueness and its implications for FEA predictions is discussed in detail. Lastly, we present two strategies to reduce parameter covariance – reduced parametrization of the material model and increased richness of the calibration data – which allow for the recovery of a unique solution.
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.
This memo addresses the calibration of anisotropic yield functions based on data obtained from a series of uniaxial tension specimens extracted from a tubular Al 7079 circular cylindrical extrusion. Achieving the calibrations completed an important step in the Experimentally Enhanced Computations (ExEC) project. The focus of the project is on novel calibration approaches that will be based on advanced diagnostics and numerical simulations with the intention of reducing the overall calibration effort. The test data used here resulted from traditional tensile tests on specimens cut at 12 orientations within the extrusion. Two anisotropic yield surfaces — Hill’s (1948) and Barlat’s (2005) — were calibrated based on the test data. The methods used to conduct the calibrations are described, and the results show that the material exhibited significant yield anisotropy. The larger number of parameters in Barlat’s yield function allowed it to fit the test data more accurately than Hill’s. Although work remains to assess the sensitivity of the calibrated model parameters to various factors, the methods implemented and the results obtained here provide bases for further work and useful benchmarks for future calibrations to be conducted using the novel approach.
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.
Ductile failure of structural metals is relevant to a wide range of engineering scenarios. Computational methods are employed to anticipate the critical conditions of failure, yet they sometimes provide inaccurate and misleading predictions. Challenge scenarios, such as the one presented in the current work, provide an opportunity to assess the blind, quantitative predictive ability of simulation methods against a previously unseen failure problem. Rather than evaluate the predictions of a single simulation approach, the Sandia Fracture Challenge relies on numerous volunteer teams with expertise in computational mechanics to apply a broad range of computational methods, numerical algorithms, and constitutive models to the challenge. This exercise is intended to evaluate the state of health of technologies available for failure prediction. In the first Sandia Fracture Challenge, a wide range of issues were raised in ductile failure modeling, including a lack of consistency in failure models, the importance of shear calibration data, and difficulties in quantifying the uncertainty of prediction [see Boyce et al. (Int J Fract 186:5–68, 2014) for details of these observations]. This second Sandia Fracture Challenge investigated the ductile rupture of a Ti–6Al–4V sheet under both quasi-static and modest-rate dynamic loading (failure in (Formula presented.) 0.1 s). Like the previous challenge, the sheet had an unusual arrangement of notches and holes that added geometric complexity and fostered a competition between tensile- and shear-dominated failure modes. The teams were asked to predict the fracture path and quantitative far-field failure metrics such as the peak force and displacement to cause crack initiation. Fourteen teams contributed blind predictions, and the experimental outcomes were quantified in three independent test labs. Additional shortcomings were revealed in this second challenge such as inconsistency in the application of appropriate boundary conditions, need for a thermomechanical treatment of the heat generation in the dynamic loading condition, and further difficulties in model calibration based on limited real-world engineering data. As with the prior challenge, this work not only documents the ‘state-of-the-art’ in computational failure prediction of ductile tearing scenarios, but also provides a detailed dataset for non-blind assessment of alternative methods.
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.
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.
In this study, ductile failure of structural metals is a pervasive issue for applications such as automotive manufacturing, transportation infrastructures, munitions and armor, and energy generation. Experimental investigation of all relevant failure scenarios is intractable, requiring reliance on computation models. Our confidence in model predictions rests on unbiased assessments of the entire predictive capability, including the mathematical formulation, numerical implementation, calibration, and execution.