Plate puncture simulations are challenging computational tasks that require advanced material models including high strain rate and thermal-mechanical effects on both deformation and failure, plus finite element techniques capable of representing large deformations and material failure. The focus of this work is on the material issues, which require large sets of experiments, flexible material models and challenging calibration procedures. In this study, we consider the puncture of 12.7 mm thick, 7075-T651 aluminum alloy plates by a cylindrical punch with a hemispherical nose and diameter of 12.7 mm. The plasticity and ductile failure models were isotropic with calibration data obtained from uniaxial tension tests at different temperatures and strain rates plus quasi-static notched tension tests and shear-dominated tests described here. Sixteen puncture experiments were conducted to identify the threshold penetration energy, mode of puncture and punch acceleration during impact, The punch was mounted on a 139 kg mass and dropped on the plates with different impact speeds. Since the mass was the same in all tests, the quantity of interest was the impact speed. The axis and velocity of the punch were perpendicular to the plate surface. The mean threshold punch speed was 3.05 m/s, and the mode of failure was plugging by thermal-mechanical shear banding accompanied by scabbing fragments. Application of the material models in simulations of the tests yielded accurate estimates of the threshold puncture speed and of the mode of failure. Time histories of the punch acceleration compared well between simulation and test. Remarkably, the success of the simulations occurred in spite of even the smallest element used being larger than the width of the shear bands.
A mass property calculator has been developed to compute the moment of inertia properties of an assemblage of parts that make up a system. The calculator can take input from spreadsheets or Creo mass property files or it can be interfaced with Phoenix Integration Model Center. The input must include the centroidal moments of inertia of each part with respect to its local coordinates, the location of the centroid of each part in the system coordinates and the Euler angles needed to rotate from the part coordinates to the system coordinates. The output includes the system total mass, centroid and mass moment of inertia properties. The input/output capabilities allow the calculator to interface with external optimizers. In addition to describing the calculator, this document serves as its user's manual. The up-to-date version of the calculator can be found in the Git repository https://cee-gitlab.sandia.gov/cj?ete/mass-properties-calculator.
This memo’s objective is to report a calibration of the J2 plasticity model with the Wilkins ductile failure criterion for 17-4 PH H1150 stainless steel under slow loading at room temperature. The calibration of the hardening function was based on uniaxial tension tests, while that of the failure model included data from tension tests on notched specimens, a butterfly specimen shear test, and a set of interrupted compression tests on shear hat specimens. The procedure was that described in, minus the rate and temperature dependence.
The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.
The objective of this work is to create an accurate elastic-plastic J2 plasticity model calibration for the Inconel 718 material at room temperature for use in finite element models. This calibration was made using a power-law hardening model of form σ = σy + $Aε^{n}_{p}$ where A and n are empirically determined constants, and σy is the proportional limit.
The objective of this work is to extend the thermal-mechanical, elastic-plastic calibrations for 304L stainless steel [1] and and 6061-T651 aluminum alloy [2] to the regime between room temperature and -40 °C. The basis to extend the calibration consisted of new uniaxial tension tests conducted at -40 °C using the same plate material stocks, circular cylindrical specimen geometries and testing apparatus as previously, followed by attempts to fit power-law hardening functions to replicate the response observed in the specimens and then extend the yield, hardening constant, hardening exponent and rate constant functions in the calibrations to cover the new temperature regime.
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. 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 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.
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 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.
Impact problems of plate-like parts by punch-like objects with relatively large mass moving at slow speeds of a few feet per second constitute a subset of impact problems of interest at Sandia. This is in contrast to small objects moving in the range of hundreds or thousands of feet per second or higher. The objective of this work is to develop a simple formula that can be used to estimate a lower bound for the puncture energy of metal plates impacted by cylindrical, essentially rigid punches of circular cross-section and at nose. Such geometry is used as a basis in the design of puncture mitigation barriers or procedures. This was accomplished by deriving an expression using non-dimensional analysis and then calibrating it based on tests results in the range of speeds of interest. Lower bounds can then be determined based on confidence intervals or factors of safety.
This memo concerns calibration of an elastic-plastic J2 material model for Ti-6Al-4V (grade 5) alloy based on tensile uniaxial stress-strain data obtained in the laboratory. In addition, tension tests on notched specimens provided data to calibrate two ductile failure models: Johnson-Cook and Wellman's tearing parameter. The tests were conducted by Kim Haulen- beek and Dave Johnson (1528) in the Structural Mechanics Laboratory (SML) during late March and early April, 2017. The SML EWP number was 4162. The stock material was a TIMETALR® 6-4 Titanium billet with 9 in. by 9 in. square section and length of 137 in. The product description indicates that it was a forging delivered in annealed condition (2 hours @ 1300oF, AC at the mill). The tensile mechanical properties reported in the material certi cation are given in Table 1, where σo represents the 0.2% strain offset yield stress, σu the ultimate stress, εf the elongation at failure and R.A. the reduction in area.