Crack growth resistance behavior of a functionally graded material: computational studies

This paper describes crack growth resistance simulation in a ceramic/metal functionally graded material (FGM) using a cohesive zone ahead of the crack front. The plasticity in the background (bulk) material follows J₂ flow theory with the flow properties determined by a volume fraction based, elastic-plastic model (extension of the original Tamura-Tomota-Ozawa model). A phenomenological, cohesive zone model with six material-dependent parameters (the cohesive energy densities and the peak cohesive tractions of the ceramic and metal phases, respectively, and two cohesive gradation parameters) describes the constitutive response of the cohesive zone. Crack growth occurs when the complete separation of the cohesive surfaces takes place. The crack growth resistance of the FGM is characterized by a rising J-integral with crack extension (averaged over the specimen thickness) computed using a domain integral (DI) formulation. The 3-D analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates solid elements with graded elastic and plastic properties and interface-cohesive elements coupled with the functionally graded cohesive zone model. The paper describes applications of the cohesive zone model and the DI method to compute the J resistance curves for both single-edge notch bend, SE(B), and single-edge notch tension, SE(T), specimens having properties of a TiB/Ti FGM. The numerical results show that the TiB/Ti FGM exhibits significant crack growth resistance behavior when the crack grows from the ceramic-rich region into the metal-rich region. Under these conditions, the J-integral is generally higher than the cohesive energy density at the crack tip even when the background material response remains linearly elastic, which contrasts with the case for homogeneous materials wherein the J-integral equals the cohesive energy density for a quasi-statically growing crack.     

J resistance behavior in functionally graded materials using cohesive zone and modified boundary layer models

This paper describes elastic–plastic crack growth resistance simulation in a ceramic/metal functionally graded material (FGM) under mode I loading conditions using cohesive zone and modified boundary layer (MBL) models. For this purpose, we first explore the applicability of two existing, phenomenological cohesive zone models for FGMs. Based on these investigations, we propose a new cohesive zone model. Then, we perform crack growth simulations for TiB/Ti FGM SE(B) and SE(T) specimens using the three cohesive zone models mentioned above. The crack growth resistance of the FGM is characterized by the J-integral. These results show that the two existing cohesive zone models overestimate the actual J value, whereas the model proposed in the present study closely captures the actual fracture and crack growth behaviors of the FGM. Finally, the cohesive zone models are employed in conjunction with the MBL model. The two existing cohesive zone models fail to produce the desired K–T stress field for the MBL model. On the other hand, the proposed cohesive zone model yields the desired K–T stress field for the MBL model, and thus yields J _R curves that match the ones obtained from the SE(B) and SE(T) specimens. These results verify the application of the MBL model to simulate crack growth resistance in FGMs.     

A boundary layer framework considering material gradation effects

This paper describes the development and application of a novel modified boundary layer (MBL) model for graded nonhomogeneous materials, e.g. functionally graded materials (FGMs). The proposed model is based on a middle-crack tension, M(T), specimen with traction boundary conditions applied to the top and lateral edges of the model. Finite element analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates elements with graded elastic and plastic properties. Elastic crack-tip fields obtained from the proposed MBL model show excellent agreement with those obtained from full models of the cracked component for homogeneous and graded nonhomogeneous materials. The K-T dominance of FGMs is investigated by comparing the actual stress fields with the asymptotic stress fields (the Williams’ solution). The examples investigated in the present study consider a crack parallel to the material gradient. Results of the present study provide insight into the K-T dominance of FGMs and also show the range of applicability of the proposed MBL model. The MBL model is applied to analyze the elastic-plastic crack-tip response of a Ti/TiB FGM SE(T) specimen. The numerical results demonstrate that the proposed MBL model captures the effect of T-stress on plastic zone size and shape, constraint effects, etc. for such configurations.     

Cohesive modeling of dynamic fracture in functionally graded materials

The dynamic fracture of functionally graded materials (FGMs) is modeled using an explicit cohesive volumetric finite element scheme that incorporates spatially varying constitutive and failure properties. The cohesive element response is described by a rate-independent bilinear cohesive failure model between the cohesive traction acting along the cohesive zone and the associated crack opening displacement. A detailed convergence analysis is conducted to quantify the effect of the material gradient on the ability of the numerical scheme to capture elastodynamic wave propagation. To validate the numerical scheme, we simulate dynamic fracture experiments performed on model FGM compact tension specimens made of a polyester resin with varying amounts of plasticizer. The cohesive finite element scheme is then used in a parametric study of mode I dynamic failure of a Ti/TiB FGM, with special emphasis on the effect of the material gradient on the initiation, propagation and arrest of the crack.     

Effect of material gradation on K-dominance of fracture specimens

This work describes the effect of material gradation (parallel to the crack plane) on stress intensity factors and K-dominance, i.e. the dominance of the singular region, of fracture specimens; SE(T), SE(B) and C(T). The extent of K-dominance is investigated by comparing the actual stress field with the Williams’ asymptotic stress field. Linear-elastic finite element analyses are performed using graded elements which incorporate graded material properties at the element level. Material gradation and crack geometry are systematically varied to perform the parametric study. Results reveal that the effect of material gradation on K_I is most pronounced when a short crack is located on the stiffer side of the fracture specimen. For a given specimen and crack geometry, the extent of K-dominance yields a curve with a peak point at a certain material gradation. Results of the present study provide valuable insight into the K-dominance of FGMs.     

A cohesive zone model for fatigue and creep–fatigue crack growth in single crystal superalloys

A numerical analysis using cohesive zone model under cyclic loading is proposed to develop a coupled predictive approach of crack growth in single crystal. The process of material damage during fatigue crack growth is described using an irreversible cohesive zone model, which governs the separation of the crack flanks and eventually leads to the formation of free surfaces. The cohesive zone element is modeled to accumulate fatigue damage during loadings and no damage during unloadings. This paper presents the damage model and its application in the study of the crack growth for precracked specimens. The use of cohesive zone approach is validated through a convergence study. Then, a general procedure of parameters calibration is presented in pure fatigue crack growth. In the last section, an extension of the cohesive zone model is presented in the case of creep–fatigue regime at high temperature. The model showed its capability to predict with a good agreement the crack growth in the case of complex loading and complex specimen geometries.     

On the determination of the cohesive zone parameters for the modeling of micro-ductile crack growth in thick specimens

The paper deals with the determination of the cohesive zone parameters (separation energy, , and cohesive strength, T _max) for the 3D finite element modeling of the micro-ductile crack growth in thick, smooth-sided compact tension specimens made of a low-strength steel. Since the cohesive zone parameters depend, in general, on the local constraint conditions around the crack tip, their values will vary along the crack front and with crack extension. The experimental determination of the separation energy via automated fracture surface analysis is not accurate enough. The basic idea is, therefore, to estimate the cohesive zone parameters, and T _max, by fitting the simulated distribution of the local crack extension values along the crack front to the experimental data of a multi-specimen J _IC-test. Furthermore, the influence of the cohesive zone parameters on the crack growth behavior is investigated. The point of crack growth initiation is determined only by the magnitude of . Both and T _max affect the crack growth rate (or the crack growth resistance), but the influence of the cohesive strength is much stronger than that of the separation energy. It turns out that T _max as well as vary along the crack front. In the center of the specimen, where plane strain conditions prevail, the separation energy is lower and the cohesive strength is higher than at the side-surface.     

Al_2O_3-Ti系梯度功能材料残余热应力的有限元分析

采用有限元方法（ＦｉｎｉｔｅＥｌｅｍｅｎｔＭｅｔｈｏｄ）对Ａｌ２Ｏ３－Ｔｉ系梯度功能材料在制备过程中产生的残余热应力进行了线弹性分析。详细讨论了梯度层数目、梯度层厚度和成分梯度指数对应力大小和分布的影响，确定了各项最佳参数。非梯度功能材料（ＮＦＧＭ）与优化后的梯度功能材料的残余热应力对比结果显示：梯度功能材料缓和热应力的效果十分显著。     

Three-dimensional modeling of ductile crack growth: Cohesive zone parameters and crack tip triaxiality

For 10 mm thick smooth-sided compact tension specimens made of a pressure vessel steel 20MnMoNi55, the interrelations between the cohesive zone parameters (the cohesive strength, T_max, and the separation energy, Γ) and the crack tip triaxiality are investigated. The slant shear-lip fracture near the side-surfaces is modeled as a normal fracture along the symmetry plane of the specimen. The cohesive zone parameters are determined by fitting the simulated crack extensions to the experimental data of a multi-specimen test. It is found that for constant cohesive zone parameters, the simulated crack extension curves show a strong tunneling effect. For a good fit between simulated and experimental crack growth, both the cohesive strength and the separation energy near the side-surface should be considerably lower than near the midsection. When the same cohesive zone parameters are applied to the 3D model and a plane strain model, the stress triaxiality in the midsection of the 3D model is much lower, the von-Mises equivalent stress is distinctly higher, and the crack growth rate is significantly lower than in the plane strain model. Therefore, the specimen must be considered as a thin specimen. The stress triaxiality varies dramatically during the initial stages of crack growth, but varies only smoothly during the subsequent stable crack growth. In the midsection region, the decrease of the cohesive strength results in a decrease of the stress triaxiality, while the decrease of the separation energy results in an increase of the triaxiality.     

THREE-DIMENSIONAL FINITE ELEMENT SIMULATION OF CRACK EXTENSION IN ALUMINIUM ALLOY 2024FC

A simulation of 3D crack extension using a cohesive zone model (CZM) has been carried out for a side-grooved compact tension specimen and a surface-crack tension specimen of aluminium 2024FC. Detailed finite element calculations were conducted by assuming crack extension only along the crack plane (mode I). For comparison, a 2D plane strain simulation is also presented. Load, displacement and crack extension histories are predicted and compared with the experiment. It is shown that the 2D approximation appears to agree reasonably well with experimental results, and that the 3D calculation gives very good agreement with test data. The determination of the CZM parameters is also discussed. Numerical results show that the CZM is a workable computational model which involves only a few microstructurally motivated phenomenological parameters for crack extension simulation.     

Simulation of ductile crack growth in thin aluminum panels using 3-D surface cohesive elements

This work describes the formulation and application of a 3-D, interface-cohesive finite element model to predict quasi-static, ductile crack extension in thin aluminum panels for mode I loading and growth. The fracture model comprises an initially zero thickness, interface element with constitutive response described by a nonlinear traction-separation relationship. Conventional volumetric finite elements model the nonlinear (elastic-plastic) response of background (bulk) material. The interface-cohesive elements undergo gradual decohesion between faces of the volumetric elements to create new traction free crack faces. The paper describes applications of the computational model to simulate crack extension in C(T) and M(T) panels made of a 2.3 mm thick, Al 2024-T3 alloy tested as part of the NASA-Langley Aging Aircraft program. Parameters of the cohesive fracture model (peak opening traction and local work of separation) are calibrated using measured load vs. outside surface crack extensions of high constraint (T-stress > 0) C(T) specimens. Analyses of low constraint M(T) specimens, having widths of 300 and 600 mm and various a/W ratios, demonstrate the capabilities of the calibrated model to predict measured loads and outside surface crack extensions. The models capture accurately the strong 3-D effects leading to various degrees of crack front tunneling in the C(T) and M(T) specimens. The predicted crack growth response shows rapid convergence with through-thickness mesh refinement. Adaptive load increment procedures to control the rate of decohesion in the interface elements leads to stable, rapidly converging iterations in the globally implicit solution procedures.     

Fatigue simulation for titanium/CFRP hybrid laminates using cohesive elements

This paper presents a new numerical approach for predicting fatigue crack growth in fiber-metal laminate (FML). Cohesive elements are used to express the complicated damage consisting of transverse cracking, splitting, and interlaminar delamination. The damage growth in the cohesive elements due to cyclic loading is represented by the conventional damage-mechanics model. The simulation was applied to notched Ti/CFRP hybrid laminates of two stacking configurations. In both cases, the crack growth rate in the titanium layer and the delamination shape agreed well with experiments reported in the literature. Complementary analysis for crack extension in the metal sheet is performed out of consideration of the damage in internal FRP layers. The numerical results demonstrated that the underlying damage modes in the FRP layer must be taken into account to predict the fatigue crack growth at the metal layer in FMLs.     

Fatigue crack growth simulations of interfacial cracks in bi-layered FGMs using XFEM

An investigation of fatigue crack growth of interfacial cracks in bi-layered materials using the extended finite element method is presented. The bi-material consists of two layers of dissimilar materials. The bottom layer is made of aluminium alloy while the upper one is made of functionally graded material (FGM). The FGM layer consists of 100 % aluminium alloy on the left side and 100 % ceramic (alumina) on the right side. The gradation in material property of the FGM layer is assumed to be exponential from the alloy side to the ceramic side. The domain based interaction integral approach is extended to obtain the stress intensity factors for an interfacial crack under thermo-mechanical load. The edge and centre cracks are taken at the interface of bi-layered material. The fatigue life of the interface crack plate is obtained using the Paris law of fatigue crack growth under cyclic mode-I, mixed-mode and thermal loads. This study reveals that the crack propagates into the FGM layer under all types of loads.     

Experimental Investigations on Deformation and Fracture Behavior of Glass Sphere Filled Epoxy Functionally Graded Materials

In this paper, deformation and fracture behavior of glass sphere filled epoxy functionally graded materials (FGM) are numerically evaluated and experimentally studied. The fabrication of the FGM is described in detail, and the spatial gradation of elastic modulus and the microscopic structure in FGM are measured and analyzed. The deformation and fracture characterization of the FGM specimen with a crack oriented along the direction of the elastic gradient under three point bend are studied by the experimental and the finite element method. The influences of crack location at both the stiff and the compliant sides of the FGM specimen on crack initiation, deformation field and stress intensity factor are analyzed. The results are: (a) The neutral-axis in the FGM specimen under three-point-bending will shift toward the stiffer side; (b) The initial fracture load increases with the increase of elastic modulus at the crack tip; (c) The elastic gradients shield a crack on the compliant side and lower the stress intensity factor when compared to the one with crack on the stiff side. These results will be useful for better design and reliable evaluation of FGM.     

A cohesive finite element formulation for modelling fracture and delamination in solids

In recent years, cohesive zone models have been employed to simulate fracture and delamination in solids. This paper presents in detail the formulation for incorporating cohesive zone models within the framework of a large deformation finite element procedure. A special Ritz-finite element technique is employed to control nodal instabilities that may arise when the cohesive elements experience material softening and lose their stress carrying capacity. A few simple problems are presented to validate the implementation of the cohesive element formulation and to demonstrate the robustness of the Ritz solution method. Finally, quasi-static crack growth along the interface in an adhesively bonded system is simulated employing the cohesive zone model. The crack growth resistance curves obtained from the simulations show trends similar to those observed in experimental studies     

Cohesive zone simulations of crack growth along a rough interface between two elastic-plastic solids

It has been known for many years that crack propagation along interfaces is influenced by interface topography or roughness profile. This has given rise to a small body of literature in which interface toughening with stochastic surface finishes, produced by grinding, rolling, or grit blasting, has been the primary focus. However, there is very little information currently available on the effect of patterned interfaces that are characterized by a minimal number of geometric parameters. In the present article, roughness-enhanced toughening of a cohesive interface between two identical materials is explored with a pure sinusoidal interface morphology that is characterized by its aspect ratio or ratio of amplitude to wavelength. Sixteen finite element meshes, each with a different aspect ratio, were constructed to study initiation and growth of a semi-infinite interface crack due to remote mode-I loading. The cohesive interface was modeled with a viscosity-modified Xu-Needleman cohesive zone law, and the solids were characterized with continuum elastic and elastic-plastic constitutive models. Predicted relationships between the aspect ratio and the macroscopic toughness point to key differences in the material models. A set of critical parameters which include the aspect ratio, material and cohesive properties is predicted such that catastrophic crack growth is inhibited due to crack blunting. A clear boundary between brittle and ductile fracture behavior is thus identified. The results suggest some guidelines for practical design of failure resistant interfaces through appropriate choice of geometric, material, and cohesive parameters.     

Comparison of cohesive zone model and linear elastic fracture mechanics for a mode I crack near a compliant/stiff interface

Cohesive zone model has been widely applied to simulate crack growth along interfaces, but its application to crack growth perpendicularly across the interface is rare. In this paper, the cohesive zone model is applied to a crack perpendicularly approaching a compliant/stiff interface in a layered material model. One aim is to understand the differences between the cohesive zone model and linear elastic fracture mechanics in simulating mode I crack growth near a compliant/stiff interface. Another aim is to understand the effects of elastic modulus mismatch and cohesive strength of the stiff layer on the crack behavior near the interface. To simulate crack growth approaching an interface, the cohesive zone model which incorporates both the energy criterion and the strength criterion is an effective method.     

A level set model for simulating fatigue-driven delamination in composites

This paper proposes a level set model for simulating delamination propagation in composites under high-cycle fatigue loading. For quasi-static loading conditions, interface elements with a cohesive law are widely used for the simulation of delamination. However, basic concepts from fatigue analysis such as the notion that the crack growth rate is a function of energy release rate cannot be embedded in existing cohesive laws. Therefore, we propose a model in which the cohesive zone is eliminated from the computation while maintaining the flexibility that the crack shape is not bound to element edges. The model is able to predict the delamination growth rate and its front shape accurately. To demonstrate the validity of the model, several tests under different fracture modes are conducted and the results are compared with experimental data, analytical solutions and results from cohesive zone analysis.     

A combined experimental-numerical investigation of crack growth in a carbon–carbon composite

A combined experimental–numerical investigation of crack growth in a carbon–carbon composite is reported. In this material, both matrix fracture and fibre bridging contribute significantly to toughness. Crack growth experiments were performed using side‐notched DCB specimens with doublers. A digital image correlation method was used to measure displacements fields on the specimen surfaces, crack extension and crack opening profiles. An effective cohesive zone law was determined from the experimental data. The effective cohesive zone law is subsequently separated into the individual contributions from matrix cracking and fibre bridging. Numerical simulation of crack growth based on this cohesive zone law and experimental data are in good agreement. Special focus of the numerical study is on the investigation of the discontinuous nature of crack growth.     

Discrete cohesive zone model for mixed-mode fracture using finite element analysis

The discrete cohesive zone model (DCZM) is implemented using the finite element (FE) method to simulate fracture initiation and subsequent growth when material non-linear effects are significant. Different from the widely used continuum cohesive zone model (CCZM) where the cohesive zone model is implemented within continuum type elements and the cohesive law is applied at each integral point, DCZM uses rod type elements and applies the cohesive law as the rod internal force vs. nodal separation (or rod elongation). These rod elements have the provision of being represented as spring type elements and this is what is considered in the present paper. A series of 1D interface elements was placed between node pairs along the intended fracture path to simulate fracture initiation and growth. Dummy nodes were introduced within the interface element to extract information regarding the mesh size and the crack path orientation. To illustrate the DCZM, three popular fracture test configurations were examined. For pure mode I, the double cantilever beam configuration, using both uniform and biased meshes were analyzed and the results show that the DCZM is not sensitive to the mesh size. Results also show that DCZM is not sensitive to the loading increment, either. Next, the end notched flexure for pure mode II and, the mixed-mode bending were studied to further investigate the approach. No convergence difficulty was encountered during the crack growth analyses. Therefore, the proposed DCZM approach is a simple but promising tool in analyzing very general two-dimensional crack growth problems. This approach has been implemented in the commercial FEA software ABAQUS^® using a user defined subroutine and should be very useful in performing structural integrity analysis of cracked structures by engineers using ABAQUS^®.