A two-dimensional in situ fatigue cohesive zone model for crack propagation in composites under cyclic loading

In this paper a two-dimensional fatigue cohesive zone model (CZM) for crack propagation in composites under cyclic loading has been formulated and validated through successful predictions of fatigue crack growth under pure and mixed mode conditions for several different composites. The proposed fatigue CZM assumes simple power-law functions for fatigue damage accumulation of which the damage parameters can be calibrated from simple fatigue tests under pure mode I and mode II conditions. The model relies solely on the in situ cohesive responses for fatigue damage rate calculation, enabling the differentiation of the local elemental load history from the global load history. An effective cycle jump strategy for high-cycle fatigue has also been proposed. It has been demonstrated that once calibrated, the fatigue CZM can predict the Paris laws for the pure modes. Furthermore, it can predict the Paris laws of any mixed-mode conditions without the need of additional empirical parameters. This is of significant practical importance because it leads to greatly reduced experimental needs for mixed mode crack propagation widely observed in composites under cyclic loads.     

Fatigue crack growth simulation in adhesively bonded composite joints

A numerical model has been developed for simulating fatigue crack growth (debonding) in adhesively bonded composite joints subjected to mode‐I, mode‐II, and mixed‐mode I + II loading conditions. The model employs a cohesive zone model described by a modified bilinear traction‐separation law. Fatigue damage in the composite adherends is not considered in the model. To account for crack divergence and reduce sensitivity of numerical results on mesh density, a crack front detection algorithm based on the effective element's length was employed. The model is implemented as a user‐defined subroutine (UMAT) in the commercial FE code LS‐DYNA. The model's input parameters, in the form of a modified Paris law, and the validation data were obtained from experimental tests conducted by the authors. It was found that the model is able to successfully simulate crack growth in the regime of the experimental data.     

Fatigue crack growth of a railway wheel

Typically, fatigue crack propagation in railway wheels is initiated at some subsurface defect and occurs under mixed mode (I–II) conditions. For a Spanish AVE train wheel, fatigue crack growth characterization of the steel in mode I, mixed mode I–II, and evaluation of crack path starting from an assumed flaw are presented and discussed.Mode I fatigue crack growth rate measurement were performed in compact tension C(T) specimens according to the ASTM E647 standard. Three different load ratios were used, and fatigue crack growth thresholds were determined according to two different procedures. Load shedding and constant maximum stress intensity factor with increasing load ratio R were used for evaluation of fatigue crack growth threshold.To model a crack growth scenario in a railway wheel, mixed mode I–II fatigue crack growth tests were performed using CTS specimens. Fatigue crack growth rates and propagation direction of a crack subjected to mixed mode loading were measured. A finite element analysis was performed in order to obtain the K_I and K_II values for the tested loading angles. The crack propagation direction for the tested mixed mode loading conditions was experimentally measured and numerically calculated, and the obtained results were then compared in order to validate the used numerical techniques.The modelled crack growth, up to final fracture in the wheel, is consistent with the expectation for the type of initial damage considered.     

Computational modeling of mixed-mode fatigue crack growth using extended finite element methods

The extended finite element method (XFEM) combined with a cyclic cohesive zone model (CCZM) is discussed and implemented for analysis of fatigue crack propagation under mixed-mode loading conditions. Fatigue damage in elastic-plastic materials is described by a damage evolution equation in the cohesive zone model. Both the computational implementation and the CCZM are investigated based on the modified boundary layer formulation under mixed-mode loading conditions. Computational results confirm that the maximum principal stress criterion gives accurate predictions of crack direction in comparison with known experiments. Further popular multi-axial fatigue criteria are compared and discussed. Computations show that the Findley criterion agrees with tensile stress dominant failure and deviates from experiments for shear failure. Furthermore, the crack propagation rate under mixed mode loading has been investigated systematically. It is confirmed that the CCZM can agree with experiments.     

A cohesive crack model coupled with damage for interface fatigue problems

An semi-analytical formulation based on the cohesive crack model is proposed to describe the phenomenon of fatigue crack growth along an interface. Since the process of material separation under cyclic loading is physically governed by cumulative damage, the material deterioration due to fatigue is taken into account in terms of interfacial cohesive properties degradation. More specifically, the damage increment is determined by the current separation and a history variable. The damage variable is introduced into the constitutive cohesive crack law in order to capture the history-dependent property of fatigue. Parametric studies are presented to understand the influences of the two parameters entering the damage evolution law. An application to a pre-cracked double-cantilever beam is discussed. The model is validated by experimental data. Finally, the effect of using different shapes of the cohesive crack law is illustrated     

Experimental characterization of dynamic crack growth behavior in CFRP adhesive interface

The dynamic crack growth behavior of adhesively bonded joints under mode I and mixed mode (I + II) loading were investigated. The split Hopkinson pressure bar (SHPB) apparatus and the digital image correlation (DIC) technique were employed to determine the mode I fracture toughness of the adhesively bonded joints during crack propagation under impact loading. The dynamic crack growth behavior for carbon fiber reinforced plastics (CFRP) adhesively bonded joints under mode I loading was studied using this method. In order to verify the proposed method, the dynamic crack growth behavior of titanium alloy adhesively bonded joints was also studied. Moreover, the crack growth behavior of CFRP adhesively bonded joints under mixed mode loading was studied using the SHPB technique. For the considered CFRP adhesively bonded joints, the fracture toughness decreased under both mode I and mixed mode loading as the loading rate increased. Microscope observation showed that a shift in the crack location occurred in the high loading tests.     

An automatic algorithm for mixed mode crack growth rate based on drop potential method

A new automatic algorithm for the assessment of mixed mode crack growth rate characteristics is presented based on the concept of an equivalent crack. The residual ligament size approach is introduced to implementation this algorithm for identifying the crack tip position on a curved path with respect to the drop potential signal. The automatic algorithm accounting for the curvilinear crack trajectory and employing an electrical potential difference was calibrated with respect to the optical measurements for the growing crack under cyclic mixed mode loading conditions. The effectiveness of the proposed algorithm is confirmed by fatigue tests performed on ST3 steel compact tension–shear specimens in the full range of mode mixities from pure mode I to pure mode II.     

History effect in fatigue crack growth under mixed-mode loading conditions

Plastic deformation within the crack tip region introduces internal stresses that modify subsequent behaviour of the crack and are at the origin of history effects in fatigue crack growth. Consequently, fatigue crack growth models should include plasticity-induced history effects. A model was developed and validated for mode I fatigue crack growth under variable amplitude loading conditions. The purpose of this study was to extend this model to mixed-mode loading conditions. Finite element analyses are commonly employed to model crack tip plasticity and were shown to give very satisfactory results. However, if millions of cycles need to be modelled to predict the fatigue behaviour of an industrial component, the finite element method becomes computationally too expensive. By employing a multiscale approach, the local results of FE computations can be brought to the global scale. This approach consists of partitioning the velocity field at the crack tip into plastic and elastic parts. Each part is partitioned into mode I and mode II components, and finally each component is the product of a reference spatial field and an intensity factor. The intensity factor of the mode I and mode II plastic parts of the velocity fields, denoted by dρ_I/dt and dρ_II/dt, allow measuring mixed-mode plasticity in the crack tip region at the global scale. Evolutions of dρ_I/dt and dρ_II/dt, generated using the FE method for various loading histories, enable the identification of an empirical cyclic elastic–plastic constitutive model for the crack tip region at the global scale. Once identified, this empirical model can be employed, with no need of additional FE computations, resulting in faster computations. With the additional hypothesis that the fatigue crack growth rate and direction can be determined from mixed-mode crack tip plasticity (dρ_I/dt and dρ_II/dt), it becomes possible to predict fatigue crack growth under I/II mixed-mode and variable amplitude loading conditions. To compare the predictions of this model with experiments, an asymmetric four point bend test system was setup. It allows applying any mixed-mode loading case from a pure mode I condition to a pure mode II. Initial experimental results showed an increase of the mode I fatigue crack growth rate after the application of a set of mode II overload cycles.     

CRACK PROPAGATION UNDER MIXED MODE (I + III) LOADING

Abstract— In this paper are presented the results of fatigue crack propagation tests on angled-slit, three point bend mixed-mode (I + III) specimens manufactured from a low pressure steam turbine rotor forging. The path of crack propagation has been studied for two mixed mode (I + III) loading conditions. It has been observed that crack growth occurs by a mode I mechanism and a model has been developed to correlate crack growth rates in mixed mode (I + III) specimens with data from pure mode I fatigue tests.     

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.     

Experimental and numerical investigations on the influence of the loading direction on the fatigue crack growth

During a service loading fatigue cracks can be subjected to a mixed mode loading if, due to the alteration of the loading direction, the basic crack modes (Modes I, II and III) are combined. An alteration of the loading direction, e.g. can occur either occasionally paired with an overload (mixed mode overload) or permanently in terms of a mixed mode block loading as a combination of normal and shear stresses.Within the scope of this paper, experimental investigations on both mixed mode overloads, which are interspersed into a Mode I baseline level loading, and mixed mode block loadings are presented. The experimental investigations show that the retardation effect decreases with an increasing amount of Mode II of the overload. Due to the block loading, the fatigue crack growth rate is retarded as well, and the crack is also deflected. The kinking angle depends on the fraction of shear stresses. Furthermore, a detailed elastic–plastic finite element analysis of the fatigue crack growth after mixed mode overloads is presented in order to understand the mechanism of the load interaction effects. By such numerical simulations, it can be shown that, due to mixed mode overloads, plastic deformations occur, which on the one hand reduce the near-tip closure and on the other hand cause a far-field closure. Also the stress distribution before and after the crack tip changes. A mixed mode overload causes lower closure and the crack tip deformations become asymmetrical, which is a reason for the smaller retardation effect of a mixed mode overload.     

The effect of mode ratio and bond interface on the fatigue behavior of a highly-toughened epoxy

The fatigue threshold and the cyclic crack growth of a highly-toughened epoxy adhesive were studied under mode I and several mixed-mode loading cases and compared with the quasi-static critical fracture energies. Four different adhesive systems were examined using steel and aluminum substrates having different surface roughness, and surface treatment. The effect of increasing the amount of mode II (increasing the phase angle) on the fatigue threshold strain energy release rate and the cyclic crack growth rate was found to be insignificant at low phase angles. However, a significant increase in the fatigue threshold and decrease in the cyclic crack growth rate was observed at higher phase angles. These trends were similar to that seen in adhesive joint fracture. Adherend surface roughness and surface preparation affected the fatigue behavior significantly, particularly at low crack speeds and high phase angles. The fatigue properties were essentially the same for both steel and aluminum adherends provided that the crack paths were cohesive. A general observation was that the fatigue crack path moved progressively closer to the more highly strained adherend under mixed-mode loading as the applied strain energy release rate and hence the crack speed, decreased. This caused mixed-mode cracks to be nearly interfacial in the threshold region.     

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^®.     

A computationally efficient cohesive zone model for fatigue

A cohesive zone model has been developed for the simulation of both high and low cycle fatigue crack growth. The developed model provides an alternative approach that reflects the computational efficiency of the well‐established envelop‐load damage model yet can deliver the accuracy of the equally well‐established loading‐unloading hysteresis damage model. A feature included in the new cohesive zone model is a damage mechanism that accumulates as a result of cyclic plastic separation and material deterioration to capture a finite fatigue life. The accumulation of damage is reflected in the loading‐unloading hysteresis curve, but additionally, the model incorporates a fast‐track feature. This is achieved by “freezing in” a particular damage state for one loading cycle over a predefined number of cycles. The new model is used to simulate mode I fatigue crack growth in austenitic stainless steel 304 at significant reduction in the computational cost.     

Fatigue crack growth in adhesively bonded composite-metal double-lap joints

This paper presents experimental and numerical investigations of the fatigue crack initiation and growth mechanism in metal-to-composite bonded double-lap joints. Fatigue tests were conducted under tension dominated loading, with crack lengths being measured optically. Examination of the fracture surface using scanning electron microscope revealed that fatigue cracks were near the interface between the co-cured adhesive and the first ply of the composite adherend. The finite element method has been used to determine the strain-energy release rate of a fatigue crack growing along the first ply of the composite. The effects of spew fillet size and crack initiation modes have also been studied by the finite element method. Comparison of the present experimental crack growth results with those measured using double-overlap joints, where the fatigue cracks were driven by pure mode II loading, indicate that the tensile mode loading has a overwhelming effect on the fatigue crack growth rates. The present results suggest that fatigue failure of metal-composite double-lap joints is mainly driven by tensile mode loading due to the peel stress.     

A cohesive XFEM model for simulating fatigue crack growth under mixed-mode loading and overloading

Structures are subjected to cyclic loads that can vary in direction and magnitude, causing constant amplitude mode I simulations to be too simplistic. This study presents a new approach for fatigue crack propagation in ductile materials that can capture mixed-mode loading and overloading. The extended finite element method is used to deal with arbitrary crack paths. Furthermore, adaptive meshing is applied to minimize computation time. A fracture process zone ahead of the physical crack tip is represented by means of cohesive tractions from which the energy release rate, and thus the stress intensity factor can be extracted for an elastic-plastic material. The approach is therefore compatible with the Paris equation, which is an empirical relation to compute the fatigue crack growth rate. Two different models to compute the cohesive tractions are compared. First, a cohesive zone model with a static cohesive law is used. The second model is based on the interfacial thick level set method in which tractions follow from a given damage profile. Both models show good agreement with a mode I analytical relation and a mixed-mode experiment. Furthermore, it is shown that the presented models can capture crack growth retardation as a result of an overload.     

Numerical analysis of the ENF and ELS tests applied to mode II fracture characterization of cortical bone tissue

The objective of this work is to verify numerically the adequacy of the ENF and the ELS tests to determine the fracture toughness under mode II loading of cortical bovine bone tissue. A data‐reduction scheme based on the specimen compliance and the equivalent crack concept is proposed to overcome the difficulties inherent to crack monitoring during its growth. A cohesive damage model was used to simulate damage initiation and growth, thus assessing the efficacy of the proposed data‐reduction scheme. The influences of the initial crack length, local strength and toughness on the measured fracture energy were analysed, taking into account the specimen length restriction. Some limitations related to spurious influence on the fracture process zone of the central loading in the ENF test, and clamping conditions in the ELS test were identified. However, it was verified that a judicious selection of the geometry allows, in both cases, a rigorous estimation of bone toughness in mode II.     

Fatigue crack growth acceleration due to intermittent overstressing in adhesively bonded CFRP joints

Double cantilever beam joints were used to investigate cohesive and interlaminar crack growth in bonded composite joints under constant and variable amplitude (VA) loading. Numerical crack growth integration was used to predict the VA fatigue life using constant amplitude data. This underestimated the fatigue crack growth rate for interlaminar cracks, indicating crack growth acceleration due to load interactions. This was also the case for cohesive cracks subjected to a moderate initial strain energy release rate (G_max). An unstable crack growth regime was also identified for the case of high initial G_max cohesive crack propagation. This behaviour is attributed to the development of a damage zone ahead of the crack tip.     

MIXED MODE SMALL CRACK GROWTH

Abstract— The fatigue crack growth behavior of small part-through cracks in 1045 steel and Inconel 718 subjected to biaxial loading has been investigated. Experiments were performed on thin-wall tubular specimens loaded in tension, torsion and combined tension torsion. Crack sizes analyzed ranged from 20 μm to 1 mm and growth rates ranged from 10^-7 to 10^-4 mm/cycle for 1045 steel and from 10^-5 to 10^-2 mm/cycle for Inconel. Nucleation and the early growth of cracks occurs on planes of maximum shear strain amplitude for both of these materials even in tensile loading. An equivalent strain based intensity factor was employed to correlate the crack growth rate under mixed mode loading conditions In loading conditions other than torsion, a transition from mode II to mode I was observed for 1045 steel. Principal strains were used to analyze mode I cracks. Cracks in Inconel 718 grow in mode II for the majority of the fatigue life. The maximum shear strain amplitude and the tensile strain normal to the maximum shear strain amplitude plane were used to calculate the strain based intensity factor for mixed mode loading.     

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.