Cohesive zone laws for fatigue crack growth: Numerical field projection of the micromechanical damage process in an elasto-plastic medium

Cohesive zone failure models are widely used to simulate fatigue crack propagation under cyclic loading, but the model parameters are phenomenological and are not closely tied to the underlying micromechanics of the problem. In this paper, we will inversely extract the cohesive zone laws for fatigue crack growth in an elasto-plastic ductile solid using a field projection method (FPM), which projects the equivalent tractions and separations at the cohesive crack-tip from field information outside the process zone. In our small-scale yielding model, a single row of discrete voids is deployed directly ahead of a crack in an elasto-plastic medium subjected to cyclic mode I K-field loading. Damage accumulation under cyclic loading is captured by the growth of voids within the micro-voiding zone ahead of the crack, while the evolution of the cohesive zone law representing the micro-voiding zone is inversely extracted via the FPM. We show that the field-projected cohesive zone law captures the essential micromechanisms of fatigue crack growth in the ductile medium: from loading and unloading hysteresis caused by void growth and plastic hardening, to the softening damage locus associated with crack propagation via a void by void growth mechanism. The results demonstrate the effectiveness of the FPM in obtaining a micromechanics-based cohesive zone law in-place of phenomenological models, which opens the way for a unified treatment of fatigue crack problems.     

Plastic strip model of circular crack under cyclic loading with damage accumulation

A cylindrical solid containing a penny-shaped crack in its mid plane is subjected to a remote tensile stress field. In the plastic region near the crack border, damage due to continuous deterioration of the material accumulates, and may lead to crack instability and crack growth. A damage model relating the crack opening displacement in the plastic zone to the fraction of the cross-sectional area occupied by voids is used to describe the conditions necessary for the onset of crack instability, fatigue crack propagation by cyclic loading, and rates of crack growth for time-dependent environmental effects.     

Analysis of crack growth and crack-tip plasticity in ductile materials using cohesive zone models

Crack initiation and crack growth resistance in elastic plastic materials, dominated by crack-tip plasticity are analyzed with the crack modeled as a cohesive zone. Two different types (exponential and bilinear) of cohesive zone models (CZMs) have been used to represent the mechanical behavior of the cohesive zones. In this work, it is suggested that different forms of CZMs (e.g., exponential, bilinear) are the manifestations of different micromechanisms-based inelastic processes that participate in dissipating energy during the fracture process and each form is specific to each material system. It is postulated that the total energy release rate comprises the plastic dissipation rate in the bounding material and the separation energy rate within the fracture process zone, the latter is determined by CZMs. The total energy release rate then becomes a function of the material properties (e.g., yield strength, strain hardening exponent) and cohesive properties of the fracture process zone (e.g., cohesive strength and cohesive energy), and the form of cohesive zone model (CZM) that determines the rate of energy dissipation in the forward and wake regions of the crack. The effects of material parameters, cohesive zone parameters as well as the form/shape of CZMs in predicting the crack growth resistance and the size of plastic zone (SPZ) surrounding the crack tip are systematically examined. It is found that in addition to the cohesive strength and cohesive energy, the form (shape) of the traction–separation law of CZM plays a very critical role in determining the crack growth resistance (R-curve) of a given material. It is further observed that the shape of the CZM corresponds to inelastic processes active in the forward and wake regions of the crack, and has a profound influence on the R-curve and SPZ.     

A comparison of cohesive zone modeling and classical fracture mechanics based on near tip stress field

A mode III crack with a cohesive zone in a power-law hardening material is studied under small scale yielding conditions. The cohesive law follows a softening path with the peak traction at the start of separation process. The stress and strain fields in the plastic zone, and the cohesive traction and separation displacement in the cohesive zone are obtained. The results show that for a modest hardening material (with a hardening exponent N = 0.3), the stress distribution in a large portion of the plastic zone is significantly altered with the introduction of the cohesive zone if the peak cohesive traction is less than two times yield stress, which implies the disparity in terms of the fracture prediction between the classical approach of elastic–plastic fracture mechanics and the cohesive zone approach. The stress distributions with and without the cohesive zone converge when the peak cohesive traction becomes infinitely large. A qualitative study on the equivalency between the cohesive zone approach and the classical linear elastic fracture mechanics indicates that smaller cracks require a higher peak cohesive traction than that for longer cracks if similar fracture initiations are to be predicted by the two approaches.     

STOCHASTIC ELASTO-PLASTIC FRACTURE ANALYSIS OF ALUMINUM FOAMS

Model I quasi-static nonlinear fracture of aluminum foams is analyzed by considering the effect of microscopic heterogeneity. Firstly, a continuum constitutive model is adopted to account for the plastic compressibility of the metallic foams. The yield strain modeled by a two- parameter Weibull-type function is adopted in the constitutive model. Then, a modified cohesive zone model is established to characterize the fracture behavior of aluminum foams with a cohesive zone ahead of the initial crack. The tensile traction versus local crack opening displacement relation is employed to describe the softening characteristics of the material. And a Weibull statistical model for peak bridging stress within the fracture process zone is used for considering microscopic heterogeneity of aluminum foams. Lastly, the influence of stochastic parameters on the curve of stress-strain is given. Numerical examples are given to illustrate the numerical model presented in this paper and the effects of Weibull parameters and material properties on J-integral are discussed.     

Finite deformation analysis of crack tip fields in plastically compressible hardening–softening–hardening solids

Crack tip fields are calculated under plane strain small scale yielding conditions. The material is characterized by a finite strain elastic–viscoplastic constitutive relation with various hardening–softening–hardening hardness functions. Both plastically compressible and plastically incompressible solids are considered. Displacements corresponding to the isotropic linear elastic mode I crack field are prescribed on a remote boundary. The initial crack is taken to be a semi-circular notch and symmetry about the crack plane is imposed. Plastic compressibility is found to give an increased crack opening displacement for a given value of the applied loading. The plastic zone size and shape are found to depend on the plastic compressibility, but not much on whether material softening occurs near the crack tip.On the other hand, the near crack tip stress and deformation fields depend sensitively on whether or not material softening occurs. The combination of plastic compressibility and softening(or softening–hardening) has a particularly strong effect on the near crack tip stress and deformation fields.     

An isolated cohesive crack in tension

The shape of the back-calculated stress-separation curves obtained from the in-situ fracture of first-year sea ice in the Arctic and Antarctic is convex, and radically different from those for concrete. In these tests, the process zone size changes with crack growth, but the nature of this change differs with the test conditions, load-control or displacement control. These results prompted a closer examination of the cohesive crack model using the simplest cracked configuration, a finite cohesive crack in an infinite elastic medium loaded in tension by a uniform stress at infinity. Different types of strain softening are examined: rectangular softening, linear softening, prescribed cohesive stresses, and prescribed cohesive crack-opening displacements. For each of these cases, crack nucleation is examined; close attention is paid to test control conditions, be they load-control or fixed-grip. The test control conditions alter the fracture, crack nucleation, crack growth, and process zone size behavior significantly. Accurate approximate solutions to linear softening are presented and examined.     

Cohesive zone laws for void growth — II. Numerical field projection of elasto-plastic fracture processes with vapor pressure

Modeling ductile fracture processes using Gurson-type cell elements has achieved considerable success in recent years. However, incorporating the full mechanisms of void growth and coalescence in cohesive zone laws for ductile fracture still remains an open challenge. In this work, a planar field projection method, combined with equilibrium field regularization, is used to extract crack-tip cohesive zone laws of void growth in an elastic-plastic solid. To this end, a single row of void-containing cell elements is deployed directly ahead of a crack in an elastic-plastic medium subjected to a remote K-field loading; the macroscopic behavior of each cell element is governed by the Gurson porous material relation, extended to incorporate vapor pressure effects. A thin elastic strip surrounding this fracture process zone is introduced, from which the cohesive zone variables can be extracted via the planar field projection method. We show that the material's initial porosity induces a highly convex traction-separation relationship — the cohesive traction reaches the peak almost instantaneously and decreases gradually with void growth, before succumbing to rapid softening during coalescence. The profile of this numerically extracted cohesive zone law is consistent with experimentally determined cohesive zone law in Part I for multiple micro-crazing in HIPS. In the presence of vapor pressure, both the cohesive traction and energy are dramatically lowered; the shape of the cohesive zone law, however, remains highly convex, which suggests that diffusive damage is still the governing failure mechanism.     

Delamination of an elastic film from an elastic–plastic substrate during adhesive contact loading and unloading

Adhesive contact between a rigid sphere and an elastic film on an elastic–perfectly plastic substrate was examined in the context of finite element simulation results. Surface adhesion was modeled by nonlinear springs obeying a force-displacement relationship governed by the Lennard–Jones potential. A bilinear cohesive zone law with prescribed cohesive strength and work of adhesion was used to simulate crack initiation and growth at the film/substrate interface. It is shown that the unloading response consists of five sequential stages: elastic recovery, interface damage (crack) initiation, damage evolution (delamination), film elastic bending, and abrupt surface separation (jump-out), with plastic deformation in the substrate occurring only during damage initiation. Substrate plasticity produces partial closure of the cohesive zone upon full unloading (jump-out), residual tensile stresses at the front of the crack tip, and irreversible downward bending of the elastic film. Finite element simulations illustrate the effects of minimum surface separation (i.e., maximum compressive surface force), work of adhesion and cohesive strength of the film/substrate interface, substrate yield strength, and initial crack size on the evolution of the surface force, residual deflection of the elastic film, film-substrate separation (debonding), crack-tip opening displacement, and contact instabilities (jump-in and jump-out) during a full load–unload cycle. The results of this study provide insight into the interdependence of contact instabilities and interfacial damage (cracking) encountered in layered media during adhesive contact loading and unloading.     

Simulation of fatigue crack growth at plastically mismatched bi-material interfaces

The present paper reports the results of a study on fatigue crack growth in a system with an interface between two elastic–plastic solids of different yield strength. Fatigue crack growth analysis is conducted by the use of a cohesive zone model. Irreversibility of material separation processes is introduced through the load history dependent degradation of the cohesive strength. Fatigue crack growth is considered to occur along the direction perpendicular to the interface and along the interface. Crack growth rate acceleration, deceleration or arrest, as well as crack bifurcation at the interface are predicted in dependence of the plastic property mismatch of the two solids and the interface properties. The outcome of the simulations is in very good agreement with trends of published experimental data.     

Damage accumulation for growth of anti-plane crack in work-hardening material

Elastic–plastic solutions of an anti-plane crack in an infinite body are used in conjunction with a continuum damage model to describe the conditions necessary for the onset of crack instability, fatigue crack propagation due to cyclic loading, and rates of crack growth due to time dependent events. A power law relates the stress to the strain of the material. The damage, which invokes nucleation, growth and coalescence of microvoids due to elevated strain, is confined to the plastic zone surrounding the crack tip. For applied loading below the yield stress, the small-scale and large-scale yielding solutions are used to determine the influence of strain hardening on crack instability and failure. Crack growth due to cyclic loading and time-dependent deformations are studied using the small-scale yielding solution of the deformation theory of plasticity.     

Simulation of local instabilities during crack propagation in the ductile–brittle transition region

The crack propagation for a cohesive zone model within an elastic-plastic material under small-scale yielding conditions is simulated numerically. The resulting crack growth resistance curves show local instabilities, so-called pop-ins even for homogeneous cohesive properties if the cohesive strength lies sufficiently close to the maximum stress of the corresponding blunting solution. The formation of secondary cracks and unloading zones in front of the actual crack tip is identified as the underlying mechanism. It is found that the shape of the cohesive law has a considerable influence on the crack arrest behavior. Furthermore, requirements to the mesh resolution are derived.     

CRACK PROPAGATION IN POLYCRYSTALLINE ELASTIC-VISCOPLASTIC MATERIALS USING COHESIVE ZONE MODELS

Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior's material properties (e.g. strain rate sensitivity) characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.     

Intergranular strain evolution near fatigue crack tips in polycrystalline metals

The deformation field near a steady fatigue crack includes a plastic zone in front of the crack tip and a plastic wake behind it, and the magnitude, distribution, and history of the residual strain along the crack path depend on the stress multiaxiality, material properties, and history of stress intensity factor and crack growth rate. An in situ, full-field, non-destructive measurement of lattice strain (which relies on the intergranular interactions of the inhomogeneous deformation fields in neighboring grains) by neutron diffraction techniques has been performed for the fatigue test of a Ni-based superalloy compact tension specimen. These microscopic grain level measurements provided unprecedented information on the fatigue growth mechanisms. A two-scale model is developed to predict the lattice strain evolution near fatigue crack tips in polycrystalline materials. An irreversible, hysteretic cohesive interface model is adopted to simulate a steady fatigue crack, which allows us to generate the stress/strain distribution and history near the fatigue crack tip. The continuum deformation history is used as inputs for the micromechanical analysis of lattice strain evolution using the slip-based crystal plasticity model, thus making a mechanistic connection between macro- and micro-strains. Predictions from perfect grain-boundary simulations exhibit the same lattice strain distributions as in neutron diffraction measurements, except for discrepancies near the crack tip within about one-tenth of the plastic zone size. By considering the intergranular damage, which leads to vanishing intergranular strains as damage proceeds, we find a significantly improved agreement between predicted and measured lattice strains inside the fatigue process zone. Consequently, the intergranular damage near fatigue crack tip is concluded to be responsible for fatigue crack growth.     

Crack growth resistance of shape memory alloys by means of a cohesive zone model

Crack growth resistance of shape memory alloys (SMAs) is dominated by the transformation zone in the vicinity of the crack tip. In this study, the transformation toughening behavior of a slowly propagating crack in an SMA under plane strain conditions and mode I deformation is numerically investigated. A small-scale transformation zone is assumed. A cohesive zone model is implemented to simulate crack growth within a finite element scheme. Resistance curves are obtained for a range of parameters that specify the cohesive traction-separation constitutive law. It is found that the choice of the cohesive strength t₀ has a great influence on the toughening behavior of the material. Moreover, the reversibility of the transformation can significantly reduce the toughening of the alloy. The shape of the initial transformation zone, as well as that of a growing crack is determined. The effect of the Young's moduli ratio of the martensite and austenite phases is examined.     

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

Theoretical development and experimental validation of a thermally dissipative cohesive zone model for dynamic fracture of amorphous polymers

A thermally dissipative cohesive zone model is developed for predicting the temperature increase at the tip of a crack propagating dynamically in a nominally brittle material exhibiting a cohesive-type failure such as crazing. The model assumes that fracture energy supplied to the crack tip region that is in excess of that needed for the creation of new free surfaces during crack advance is converted to heat within the cohesive zone. Bulk dissipation mechanisms, such as plasticity, are not accounted for. Several cohesive traction laws are examined, and the model is then used to make predictions of crack tip heating at various crack propagation speeds in the nominally brittle amorphous polymer PMMA, observed to fail by a crazing-type mechanism. The heating predictions are compared to experimental data where the temperature field surrounding a high speed crack in PMMA was measured. Measurements are made in real time using a multi-point high speed HgCdTe infrared radiation detector array. At the same time as temperature, simultaneous measurement of fracture energy is made by a strain gauge technique, and crack tip speed is monitored through a resistance ladder method. Material strength can be estimated through uniaxial tension tests, thus minimizing the need for parameter fitting in the stress-opening traction law. Excellent agreement between experiments and theory is found for two of the cohesive traction law temperature predictions, but only for the case where a single craze is active during the dynamic fracture of PMMA, i.e. crack tip speed up to approximately 0.2c_R. For higher speed fracture where subsurface damage becomes prominent, the line dissipation model of a cohesive zone is inadequate, and a distributed damage model is needed.     

Mode-3 spontaneous crack propagation along functionally graded bimaterial interfaces

The effects of combining functionally graded materials (FGMs) of different inhomogeneous property gradients on the mode-3 propagation characteristics of an interfacial crack are numerically investigated. Spontaneous interfacial crack propagation simulations were performed using the newly developed spectral scheme. The numerical scheme derived and implemented in the present work can efficiently simulate planar crack propagation along functionally graded bimaterial interfaces. The material property inhomogeneity was assumed to be in the direction normal to the interface. Various bimaterial combinations were simulated by varying the material property inhomogeneity length scale. Our parametric study showed that the inclusion of a softening type FGM in the bimaterial system leads to a reduction in the fracture resistance indicated by the increase in crack propagation velocity and power absorbed. An opposite trend of increased fracture resistance was predicted when a hardening material was included in the bimaterial system. The cohesive tractions and crack opening displacements were altered due to the material property inhomogeneity, but the stresses ahead of the cohesive zone remained unaffected.