A Method for Combining Experimentation and Molecular Dynamics Simulation to Improve Cohesive Zone Models for Metallic Microstructures

Fracture processes within a material begin at the nanometer length scale at which the formation, propagation, and interaction of fundamental damage mechanisms occur. Physics-based modeling of these atomic processes quickly becomes computationally intractable as the system size increases. Thus, a multiscale modeling method, based on the aggregation of fundamental damage processes occurring at the nanoscale within a cohesive zone model, is under development and will enable computationally feasible and physically meaningful microscale fracture simulation in polycrystalline metals. This method employs atomistic simulation to provide an optimization loop with an initial prediction of a cohesive zone model (CZM). This initial CZM is then applied at the crack front region within a finite element model. The optimization procedure iterates upon the CZM until the finite element model acceptably reproduces the near-crack-front displacement fields obtained from experimental observation. With this approach, a comparison can be made between the original CZM predicted by atomistic simulation and the converged CZM that is based on experimental observation. Comparison of the two CZMs gives insight into how atomistic simulation scales

Discrete crack growth analysis methodology for through cracks in pressurized fuselage structures

A methodology for simulating the growth of long through cracks in the skin of pressurized aircraft fuselage structures is described. Crack trajectories are allowed to be arbitrary and are computed as part of the simulation. The interaction between the mechanical loads acting on the superstructure and the local structural response near the crack tips is accounted for by employing a hierarchical modeling strategy. The structural response for each cracked configuration is obtained using a geometrically nonlinear shell finite element analysis procedure. Four stress intensity factors, two for membrane behavior and two for bending using Kirchhoff plate theory, are computed using an extension of the modified crack closure integral method. Crack trajectories are determined by applying the maximum tangential stress criterion. Crack growth results in localized mesh deletion, and the deletion regions are remeshed automatically using a newly developed all-quadrilateral meshing algorithm. The effectiveness of the methodology and its applicability to performing practical analyses of realistic structures is demonstrated by simulating curvilinear crack growth in a fuselage panel that is representative of a typical narrow-body aircraft. The predicted crack trajectory and fatigue life compare well with measurements of these same quantities from a full-scale pressurized panel test

Residual Strength Prediction of Fuselage Structures with Multiple Site Damage

This paper summarizes recent results on simulating full-scale pressure tests of wide body, lap-jointed fuselage panels with multiple site damage (MSD). The crack tip opening angle (CTOA) fracture criterion and the FRANC3D/STAGS software program were used to analyze stable crack growth under conditions of general yielding. The link-up of multiple cracks and residual strength of damaged structures were predicted. Elastic-plastic finite element analysis based on the von Mises yield criterion and incremental flow theory with small strain assumption was used. A global-local modeling procedure was employed in the numerical analyses. Stress distributions from the numerical simulations are compared with strain gage measurements. Analysis results show that accurate representation of the load transfer through the rivets is crucial for the model to predict the stress distribution accurately. Predicted crack growth and residual strength are compared with test data. Observed and predicted results both indicate that the occurrence of small MSD cracks substantially reduces the residual strength. Modeling fatigue closure is essential to capture the fracture behavior during the early stable crack growth. Breakage of a tear strap can have a major influence on residual strength prediction

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework

Fracture mechanics life analytical methods verification testing

The objective was to evaluate NASCRAC (trademark) version 2.0, a second generation fracture analysis code, for verification and validity. NASCRAC was evaluated using a combination of comparisons to the literature, closed-form solutions, numerical analyses, and tests. Several limitations and minor errors were detected. Additionally, a number of major flaws were discovered. These major flaws were generally due to application of a specific method or theory, not due to programming logic. Results are presented for the following program capabilities: K versus a, J versus a, crack opening area, life calculation due to fatigue crack growth, tolerable crack size, proof test logic, tearing instability, creep crack growth, crack transitioning, crack retardation due to overloads, and elastic-plastic stress redistribution. It is concluded that the code is an acceptable fracture tool for K solutions of simplified geometries, for a limited number of J and crack opening area solutions, and for fatigue crack propagation with the Paris equation and constant amplitude loads when the Paris equation is applicable

Failure patterns caused by localized rise in pore-fluid overpressure and effective strength of rocks

In order to better understand the interaction between pore-fluid overpressure and failure patterns in rocks we consider a porous elasto-plastic medium in which a laterally localized overpressure line source is imposed at depth below the free surface. We solve numerically the fluid filtration equation coupled to the gravitational force balance and poro-elasto-plastic rheology equations. Systematic numerical simulations, varying initial stress, intrinsic material properties and geometry, show the existence of five distinct failure patterns caused by either shear banding or tensile fracturing. The value of the critical pore-fluid overpressure at the onset of failure is derived from an analytical solution that is in excellent agreement with numerical simulations. Finally, we construct a phase-diagram that predicts the domains of the different failure patterns and at the onset of failure

A Geometric Approach to Modeling Microstructurally Small Fatigue Crack Formation

The objective of this paper is to develop further a framework for computationally modeling microstructurally small fatigue crack growth in AA 7075-T651 [1]. The focus is on the nucleation event, when a crack extends from within a second-phase particle into a surrounding grain, since this has been observed to be an initiating mechanism for fatigue crack growth in this alloy. It is hypothesized that nucleation can be predicted by computing a non-local nucleation metric near the crack front. The hypothesis is tested by employing a combination of experimentation and nite element modeling in which various slip-based and energy-based nucleation metrics are tested for validity, where each metric is derived from a continuum crystal plasticity formulation. To investigate each metric, a non-local procedure is developed for the calculation of nucleation metrics in the neighborhood of a crack front. Initially, an idealized baseline model consisting of a single grain containing a semi-ellipsoidal surface particle is studied to investigate the dependence of each nucleation metric on lattice orientation, number of load cycles, and non-local regularization method. This is followed by a comparison of experimental observations and computational results for microstructural models constructed by replicating the observed microstructural geometry near second-phase particles in fatigue specimens. It is found that orientation strongly influences the direction of slip localization and, as a result, in uences the nucleation mechanism. Also, the baseline models, replication models, and past experimental observation consistently suggest that a set of particular grain orientations is most likely to nucleate fatigue cracks. It is found that a continuum crystal plasticity model and a non-local nucleation metric can be used to predict the nucleation event in AA 7075-T651. However, nucleation metric threshold values that correspond to various nucleation governing mechanisms must be calibrated

The Sandia Fracture Challenge: blind round robin predictions of ductile tearing

Existing and emerging methods in computational mechanics are rarely validated against problems with an unknown outcome. For this reason, Sandia National Laboratories, in partnership with US National Science Foundation and Naval Surface Warfare Center Carderock Division, launched a computational challenge in mid-summer, 2012. Researchers and engineers were invited to predict crack initiation and propagation in a simple but novel geometry fabricated from a common off-the-shelf commercial engineering alloy. The goal of this international Sandia Fracture Challenge was to benchmark the capabilities for the prediction of deformation and damage evolution associated with ductile tearing in structural metals, including physics models, computational methods, and numerical implementations currently available in the computational fracture community. Thirteen teams participated, reporting blind predictions for the outcome of the Challenge. The simulations and experiments were performed independently and kept confidential. The methods for fracture prediction taken by the thirteen teams ranged from very simple engineering calculations to complicated multiscale simulations. The wide variation in modeling results showed a striking lack of consistency across research groups in addressing problems of ductile fracture. While some methods were more successful than others, it is clear that the problem of ductile fracture prediction continues to be challenging. Specific areas of deficiency have been identified through this effort. Also, the effort has underscored the need for additional blind prediction-based assessments

A Criterion for Brittle Failure of Rocks Using the Theory of Critical Distances

This paper presents a new analytical criterion for brittle failure of rocks and heavily overconsolidated soils. Griffith’s model of a randomly oriented defect under a biaxial stress state is used to keep the criterion simple. The Griffith’s criterion is improved because the maximum tensile strength is not evaluated at the boundary of the defect but at a certain distance from the boundary, known as the critical distance. This fracture criterion is known as the Point Method, and is part of the Theory of Critical Distances, which is utilized in fracture mechanics. The proposed failure criterion has two parameters: the inherent tensile strength, ó0, and the ratio of the half-length of the initial crack/flaw to the critical distance, a/L. These parameters are difficult to measure but they may be correlated with the uniaxial compressive and tensile strengths, óc and ót. The proposed criterion is able to reproduce the common range of strength ratios for rocks and heavily overconsolidated soils (óc/ót=3-50) and the influence of several microstructural rock properties, such as texture and porosity. Good agreement with laboratory tests reported in the literature is found for tensile and low confining stresses.The work presented was initiated during a research project on “Structural integrity assessments of notch-type defects", for the Spanish Ministry of Science and Innovation (Ref.: MAT2010-15721)