A cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.
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.
The method of thin-wall tube torsion to characterize metal’s shear response is well-known. Unfortunately, the thin wall tube specimen tends to buckle before reaching large shear deformation and failure. An alternative technique, which has rarely been considered, is Nadai’s surface stress method (Nadai, Theory of Flow and Fracture of Solids. McGraw-Hill, New York, 1950). It derives shear stress-strain curve from the torque-twist relationship of a solid bar. Although the analysis is more complex due to nonlinear shear stress distribution along the radius, the deformation is stable through large shear deformation to failure. Solid bar torsion experiments were conducted to study large shear deformation of Al6061-T6. Two experiments were described in this study. Since few tests were available in the literature, these experiments were to explore the large deformation behaviors of an engineering alloy and the application of modern measurement techniques, such as 3D DIC method, under torsion. Results show during twisting, the surface shear strain distribution was uniform initially and then localized on a narrow band; eventually, the specimen was cracked and failed within the band. Depending on the specimen size, the twist could be greater than 360°. Details are discussed.
Limited fatigue data exists for small-volume welded austenitic stainless steel components typically employed in hydrogen infrastructure due to the difficulty of testing these components with conventional specimen designs. To assess the fatigue performance of orbital tube welds of austenitic stainless steels, a hole-drilled tubular specimen was designed to produce a stress concentration in the center of the orbital weld. Fatigue life testing was performed on welded and non-welded 316L stainless steel hole-drilled tubular specimens, and the effects of hydrogen were evaluated by testing specimens with no added hydrogen and with internal hydrogen introduced through gaseous precharging. When accounting for the differences in flow stress caused by microstructural variations and the presence of internal hydrogen, the total fatigue life and fatigue crack initiation life of the welded and non-welded tubes were comparable and were reduced by internal hydrogen. In addition, the fatigue life results produced with the hole drilled tubular specimens were consistent with fatigue life data from circumferentially notched stainless steel specimens that have a similar elastic stress concentration factor. To better understand the mechanics of this specimen geometry, mechanics modeling was performed to compare the stress and strain distributions that develop at the stress concentration in the hole-drilled tubular and circumferentially notched specimens during fatigue cycling.
Crystal plasticity-finite element method (CP-FEM) is now widely used to understand the mechanical response of polycrystalline materials. However, quantitative mesh convergence tests and verification of the necessary size of polycrystalline representative volume elements (RVE) are often overlooked in CP-FEM simulations. Mesh convergence studies in CP-FEM models are more challenging compared to conventional finite element analysis (FEA) as they are not only computationally expensive but also require explicit discretization of individual grains using many finite elements. Resolving each grains within a polycrystalline domain complicates mesh convergence study since mesh convergence is strongly affected by the initial crystal orientations of grains and local loading conditions. In this work, large-scale CP-FEM simulations of single crystals and polycrystals are conducted to study mesh sensitivity in CP-FEM models. Various factors that may affect the mesh convergence in CP-FEM simulations, such as initial textures, hardening models and boundary conditions are investigated. In addition, the total number of grains required to obtain adequate RVE is investigated. This work provides a list of guidelines for mesh convergence and RVE generation in CP-FEM modeling.
The third Sandia Fracture Challenge highlighted the geometric and material uncertainties introduced by modern additive manufacturing techniques. Tasked with the challenge of predicting failure of a complex additively-manufactured geometry made of 316L stainless steel, we combined a rigorous material calibration scheme with a number of statistical assessments of problem uncertainties. Specifically, we used optimization techniques to calibrate a rate-dependent and anisotropic Hill plasticity model to represent material deformation coupled with a damage model driven by void growth and nucleation. Through targeted simulation studies we assessed the influence of internal voids and surface flaws on the specimens of interest in the challenge which guided our material modeling choices. Employing the Kolmogorov–Smirnov test statistic, we developed a representative suite of simulations to account for the geometric variability of test specimens and the variability introduced by material parameter uncertainty. This approach allowed the team to successfully predict the failure mode of the experimental test population as well as the global response with a high degree of accuracy.
The heterogeneity in mechanical fields introduced by microstructure plays a critical role in the localization of deformation. To resolve this incipient stage of failure, it is therefore necessary to incorporate microstructure with sufficient resolution. On the other hand, computational limitations make it infeasible to represent the microstructure in the entire domain at the component scale. In this study, the authors demonstrate the use of concurrent multiscale modeling to incorporate explicit, finely resolved microstructure in a critical region while resolving the smoother mechanical fields outside this region with a coarser discretization to limit computational cost. The microstructural physics is modeled with a high-fidelity model that incorporates anisotropic crystal elasticity and rate-dependent crystal plasticity to simulate the behavior of a stainless steel alloy. The component-scale material behavior is treated with a lower fidelity model incorporating isotropic linear elasticity and rate-independent J2 plasticity. The microstructural and component scale subdomains are modeled concurrently, with coupling via the Schwarz alternating method, which solves boundary-value problems in each subdomain separately and transfers solution information between subdomains via Dirichlet boundary conditions. In this study, the framework is applied to model incipient localization in tensile specimens during necking.
The heterogeneity in mechanical fields introduced by microstructure plays a critical role in the localization of deformation. To resolve this incipient stage of failure, it is therefore necessary to incorporate microstructure with sufficient resolution. On the other hand, computational limitations make it infeasible to represent the microstructure in the entire domain at the component scale. In this study, the authors demonstrate the use of concurrent multi- scale modeling to incorporate explicit, finely resolved microstructure in a critical region while resolving the smoother mechanical fields outside this region with a coarser discretization to limit computational cost. The microstructural physics is modeled with a high-fidelity model that incorporates anisotropic crystal elasticity and rate-dependent crystal plasticity to simulate the behavior of a stainless steel alloy. The component-scale material behavior is treated with a lower fidelity model incorporating isotropic linear elasticity and rate-independent J 2 plas- ticity. The microstructural and component scale subdomains are modeled concurrently, with coupling via the Schwarz alternating method, which solves boundary-value problems in each subdomain separately and transfers solution information between subdomains via Dirichlet boundary conditions. Beyond cases studies in concurrent multiscale, we explore progress in crystal plastic- ity through modular designs, solution methodologies, model verification, and extensions to Sierra/SM and manycore applications. Advances in conformal microstructures having both hexahedral and tetrahedral workflows in Sculpt and Cubit are highlighted. A structure-property case study in two-phase metallic composites applies the Materials Knowledge System to local metrics for void evolution. Discussion includes lessons learned, future work, and a summary of funded efforts and proposed work. Finally, an appendix illustrates the need for two-way coupling through a single degree of freedom.
Circumferentially notched specimens of several austenitic stainless steel alloys subjected to positive load ratio, load-controlled fatigue have been cycled to failure in high pressure hydrogen gas. The number of cycles to failure for a given applied stress amplitude varies among the alloys tested indicating that control of fatigue life in hydrogen environments may be attained through informed alloy selection. The number of cycles to initiate a crack does not vary significantly among the alloys tested, however the total life to failure varied by over an order of magnitude. This difference in life is attributed to variations of the stress and strain fields ahead of the blunt notch. These fields are influenced by the strength and strain hardening characteristics of each alloy and they dictate the driving force for fatigue crack growth while the crack is small.
The second Sandia Fracture Challenge illustrates that predicting the ductile fracture of Ti-6Al-4V subjected to moderate and elevated rates of loading requires thermomechanical coupling, elasto-thermo-poro-viscoplastic constitutive models with the physics of anisotropy and regularized numerical methods for crack initiation and propagation. We detail our initial approach with an emphasis on iterative calibration and systematically increasing complexity to accommodate anisotropy in the context of an isotropic material model. Blind predictions illustrate strengths and weaknesses of our initial approach. We then revisit our findings to illustrate the importance of including anisotropy in the failure process. Mesh-independent solutions of continuum damage models having both isotropic and anisotropic yields surfaces are obtained through nonlocality and localization elements.
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.
Laser welds are prevalent in complex engineering systems and they frequently govern failure. The weld process often results in partial penetration of the base metals, leaving sharp crack-like features with a high degree of variability in the geometry and material properties of the welded structure. Furthermore, accurate finite element predictions of the structural reliability of components containing laser welds requires the analysis of a large number of finite element meshes with very fine spatial resolution, where each mesh has different geometry and/or material properties in the welded region to address variability. We found that traditional modeling approaches could not be efficiently employed. Consequently, a method is presented for constructing a surrogate model, based on stochastic reduced-order models, and is proposed to represent the laser welds within the component. Here, the uncertainty in weld microstructure and geometry is captured by calibrating plasticity parameters to experimental observations of necking as, because of the ductility of the welds, necking – and thus peak load – plays the pivotal role in structural failure. The proposed method is exercised for a simplified verification problem and compared with the traditional Monte Carlo simulation with rather remarkable results.
Team Sandia California (Team H) used the Sandia code SIERRA Solid Mechanics: Implicit (SIERRA SM) to model the SFC2 challenge problem. SIERRA SM is a Lagrangian, three-dimensional, implicit code for the analysis of solids and structures. It contains a versatile library of continuum and structural elements, and an extensive library of material models. For all SFC2 related simulations, our team used Q1P0, 8 node hexahedral elements with element side lengths on the order 0.175 mm in failure regions. To model crack initiation and failure, element death removed elements from the simulation according to a continuum damage model. SIERRA SM’s implicit dynamics, implemented with an HHT time integration scheme for numerical damping [1], was used to model the unstable failure modes of the models. We chose SIERRA SM’s isotropic Elasto Viscoplastic material model for our simulations because it contains most of the physics required to accurately model the SFC2 challenge problem such as the flexibility to include temperature and rate dependence for a material.
Fracture or tearing of ductile metals is a pervasive engineering concern, yet accurate prediction of the critical conditions of fracture remains elusive. Sandia National Laboratories has been developing and implementing several new modeling methodologies to address problems in fracture, including both new physical models and new numerical schemes. The present study provides a double-blind quantitative assessment of several computational capabilities including tearing parameters embedded in a conventional finite element code, localization elements, extended finite elements (XFEM), and peridynamics. For this assessment, each of four teams reported blind predictions for three challenge problems spanning crack initiation and crack propagation. After predictions had been reported, the predictions were compared to experimentally observed behavior. The metal alloys for these three problems were aluminum alloy 2024-T3 and precipitation hardened stainless steel PH13-8Mo H950. The predictive accuracies of the various methods are demonstrated, and the potential sources of error are discussed.
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.
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.