A moving cohesive interface model for fracture in creeping materials

We present a new cohesive interface model for quasi-static creep crack growth that is implemented within a moving-grid finite element model. A pseudo crack tip separates the cohesive process zone from the free surfaces of the crack. The moving-grid formulation models continuous crack advance by describing relative motion between the pseudo crack tip and the material. This eliminates the need for extensive mesh refinement away from the current crack-tip location and supports both transient and direct steady-state solutions. A traction-separation law determines the energetic properties of the decohesion process and generates a simple criterion for crack advance. The new formulation remedies a problem in earlier models which permit a crack to heal on unloading. Numerical examples demonstrate the moving cohesive interface model in studies of steady-state crack growth. Adaptive grid refinement is used to control the accuracy of the solution.     

A single-domain dual-boundary-element formulation incorporating a cohesive zone model for elastostatic cracks

A cracked elastostatic structure is artificially divided into subdomains of simpler topology such that the well-developed classic dual integral equations can be applied appropriately to each domain. Applying the continuity and equilibrium conditions along artificial boundaries and properties of the integral kernels a single-domain dual-boundary-integral equation formulation is derived for a cracked elastic structure. A cohesive zone model is used to model the crack tip processes and is coupled with the single-domain dual-boundary-integral equation formulation; the resulting nonlinear equations are solved using the iterative method of successive-over-relaxation. The constitutive law used for a crack includes three parts: a law relating cohesive force to crack displacement difference when a crack is opening, a characterization of tangential interaction between crack surfaces when the crack surfaces are in contact, and a maximum principal stress criterion of crack advance. Incorporation of local unloading effect of the cohesive zone material has enabled a simulation of fracture with initial damage, partial development of the failure process zone at structural instability and multiple crack interaction. Some of the features of the method are demonstrated by considering three examples. The first problem is a single-edge-cracked specimen that exhibits a snap-back instability. The second example is the development of wing cracks from an angled crack under compression. The last example demonstrates the capability to consider mixed-mode crack growth and interaction of cracks. Thus, the problem of crack growth has been reduced to the determination of the cohesive model for the fracture process. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

A modified cohesive zone model for a high‐speed expanding crack

The present work studies a self‐similar high‐speed expanding crack of mode‐I in a ductile material with a modified cohesive zone model. Compared with existing Dugdale models for moving crack, the new features of the present model are that the normal stress parallel to crack faces is included in the yielding condition in the cohesive zone and traction force in the cohesive zone can be non‐uniform. For a ductile material defined by von Mises criterion without hardening, the present model confirms that the normal stress parallel to crack face increases with increasing crack speed and can be even larger than the normal traction in the cohesive zone, which justifies the necessity of including the normal stress parallel to the crack faces in the yielding condition at high crack speed. In addition, strain hardening effect is examined based on a non‐uniform traction distribution in the cohesive zone.     

A two-field modified Lagrangian formulation for robust simulations of extrinsic cohesive zone models

This paper presents the robust implementation of a cohesive zone model based on extrinsic cohesive laws (i.e. laws involving an infinite initial stiffness). To this end, a two-field Lagrangian weak formulation in which cohesive tractions are chosen as the field variables along the crack’s path is presented. Unfortunately, this formulation cannot model the infinite compliance of the broken elements accurately, and no simple criterion can be defined to determine the loading–unloading change of state at the integration points of the cohesive elements. Therefore, a modified Lagrangian formulation using a fictitious cohesive traction instead of the classical cohesive traction as the field variable is proposed. Thanks to this change of variable, the cohesive law becomes an increasing function of the equivalent displacement jump, which eliminates the problems mentioned previously. The ability of the proposed formulations to simulate fracture accurately and without field oscillations is investigated through three numerical test examples.     

Cohesive models for damage evolution in laminated composites

A trend in the last decade towards models in which nonlinear crack tip processes are represented explicitly, rather than being assigned to a point process at the crack tip (as in linear elastic fracture mechanics), is reviewed by a survey of the literature. A good compromise between computational efficiency and physical reality seems to be the cohesive zone formulation, which collapses the effect of the nonlinear crack process zone onto a surface of displacement discontinuity (generalized crack). Damage mechanisms that can be represented by cohesive models include delamination of plies, large splitting (shear) cracks within plies, multiple matrix cracking within plies, fiber rupture or microbuckling (kink band formation), friction acting between delaminated plies, process zones at crack tips representing crazing or other nonlinearity, and large scale bridging by through-thickness reinforcement or oblique crack-bridging fibers. The power of the technique is illustrated here for delamination and splitting cracks in laminates. A cohesive element is presented for simulating three-dimensional, mode-dependent process zones. An essential feature of the formulation is that the delamination crack shape can follow its natural evolution, according to the evolving mode conditions calculated within the simulation. But in numerical work, care must be taken that element sizes are defined consistently with the characteristic lengths of cohesive zones that are implied by the chosen cohesive laws. Qualitatively successful applications are reported to some practical problems in composite engineering, which cannot be adequately analyzed by conventional tools such as linear elastic fracture mechanics and the virtual crack closure technique. The simulations successfully reproduce experimentally measured crack shapes that have been reported in the literature over a decade ago, but have not been reproduced by prior models.     

On damage accumulations in the cyclic cohesive zone model for XFEM analysis of mixed-mode fatigue crack growth

Predicting mixed-mode fatigue crack propagation is an important and troublesome issue in structure assessment for decades. In the present paper an extended finite element method (XFEM) combined with a new cyclic cohesive zone model (CCZM) is introduced for simulating fatigue crack propagation under mixed-mode loading conditions, which has been implemented in the commercial general purpose software ABAQUS. The algorithm allows introducing a new crack surface at arbitrary locations and directions in a finite element mesh, without re-meshing. The cyclic cohesive zone model is based on the known S–N curves and Goodman diagram for metallic materials and validated by uniaxial tension results. Furthermore, the sensitivity of the model parameter is investigated for mixed-mode fatigue. The virtual crack closure technique has been extended to the cohesive zone model and proposed to calculate the energy release rate for the generalized Paris’ law. Finally, the crack propagation rate and direction under mixed-mode fatigue loading conditions are studied.     

A new nonlocal cohesive stress law and its applicable range

In this paper, we attempt to provide a new analytical method to determine the cohesive law in the framework of nonlocal continuum mechanics. Firstly, the equivalence between the cohesive stress and the surface-induced traction (nonlocal surface residual) is established on the basis of the nonlocal stress boundary condition. Then a new cohesive stress law is derived logically from the perspective of rational mechanics, which characterizes the dependence of the cohesive stress on the crack opening displacement (COD) within the cohesive zone. Finally, we apply this new cohesive crack model to two fracture examples with different cohesive zone sizes, and investigate the stress field ahead of the crack tip in detail. The results show that the stress singularity at the crack tip is removed, and the maximum stress occurs within the cohesive zone away from the crack tip. Moreover, the stress in the large-scale cohesive zone drops rapidly to a constant approaching zero, exhibiting a stronger softening behavior.     

Damage Dependent Constitutive Behavior and Energy Release Rate for a Cohesive Zone in a Thermoviscoelastic Solid

In this paper a cohesive zone is introduced ahead of a crack tip in order to avoid the singularity at the crack tip. By applying thermodynamics to the cohesive zone and the surrounding body, a fracture criterion will be established so that the inelastic energy dissipation both in the cohesive zone and the surrounding bulk material can be distinguished from the energy released by fracture, and the propagation of crack can be predicted. In addition, the cohesive zone constitutive equation is constructed utilizing the Helmholtz free energy in the form of a single hereditary integral for a nonlinear viscoelastic material. The resulting constitutive model for the cohesive zone contains an internal state variable which represents the damage state within the cohesive zone. When the cohesive zone opening displacement is known, the energy release rate is thus history dependent, which is expressed in terms of the damage state, the length of separation in the cohesive zone and the geometric configuration of the cohesive zone opening displacement. Example results contained herein demonstrate this effect. This revised version was published online in July 2006 with corrections to the Cover Date.     

A numerical analysis of constraint effects in fatigue crack growth by use of an irreversible cohesive zone model

The analysis of constraint effects in fatigue crack growth in multi-layer structures is discussed. The process of material separation under cyclic loading is described by a cohesive zone model (CZM) with an irreversible constitutive relationship. The traction–separation behavior does not follow a predefined path, but is dependent on the evolution of the damage dependent cohesive zone properties. A modified boundary layer model is used in simulations of fatigue crack growth along the centerline crack of the metal layer sandwiched between two elastic substrates. Fatigue crack growth is computed for a series of values of metal layer thickness under constant and variable amplitude loading conditions. The results of the computations demonstrate that certain combinations of load magnitude, layer thickness and material properties results in significant constrain effects in fatigue crack growth. The influence of these constraint effects on fatigue crack growth rates and on crack closure processes is determined. The evolutions of the traction–separation law, the accumulated and current plastic zones, as well as the stress fields during the crack propagation are discussed.     

On the effect of the residual stresses on the crack opening displacement in a cracked sheet

The residual stress and displacement fields caused by localized plastic flow near a mode I crack tip in a sheet under plane stress conditions are investigated. The present study founds on the classical Dugdale scheme of the plastic flow localization. The residual stress field is considered to be induced by reversed plastic flow near the crack tip caused by an unloading. As it is found the residual stresses around the crack compress the crack tip, while the residual tensile stresses in a distant from the crack tip zone occur. It is shown the maximum residual tensile stresses can reach the significant value of the one third of the yield limit. The length of the compressed plastic zone and the residual displacement distributions are obtained. The exact formula for the residual crack opening displacement to estimate the crack closure is found. Then the next loading of the cracked plate is considered. It is shown that the second loading causes the origin of two plastic zones localized near the crack tip and at the point, where the maximum residual tensile stresses are concentrated. Again, according to the Dugdale scheme of the plastic localization, both the plastic flow zones are modelled as narrow stripes on the line extending the crack. To determine three non-dimensional parameters, characterizing the position of the segment-like plastic flow zones, a non-linear system of equations is obtained and analyzed. The exact formula for the crack opening displacement after a loading–unloading cycle is obtained. An asymptotic analysis (as the linear dimension of the distant plastic flow zone compared with the actual crack length is small) is given. It shows that the effect of the distant plastic flow zone appears as some complementary crack closure.     

Effect of overload on crack closure in thick and thin specimens via digital image correlation

The displacement field of compact tension (CT) specimens have been mapped by digital image correlation (DIC) local to growing fatigue cracks to study overload effects for plane stress and plane strain. We have extracted crack opening displacement (ΔCOD) and stress intensity (K) determined by a Muskhelishvili fit to the crack tip displacement field to infer the closure load. In both cases a classical knee was observed upon unloading consistent with closure which disappeared during the accelerated growth following OL, before increasing during retardation. In both cases following OL the crack growth rate is perturbed for a distance similar to the plastic zone.     

The structure in the vicinity of a crack tip: A general theory based on the cohesive zone model

Consideration is given to the development of a general theory for the determination of the structure in the vicinity of the tip of a crack within an isotropic elastic solid. The basis of the model is the presence of a Barenblatt-type zone, in which cohesive stresses act, thereby allowing non-linear behaviour to he taken into account. The theory enables the stress and displacements within the cohesive zone, the zone size and critical crack extension stress to be determined, for a wide range of stress-displacement laws describing the non-linear behaviour within the cohesive zone, provided its size is small compared with that of the crack.     

Cohesive Fracture Model Based on Necking

A linear hardening model together with a linear elastic background material is first used to discuss some aspects of the mathematical and physical limitations and constraints on cohesive laws. Using an integral equation approach together with the cohesive crack assumption, it is found that in order to remove the stress singularity at the tip of the cohesive zone, the cohesive law must have a nonzero traction at the initial zero opening displacement. A cohesive zone model for ductile metals is then derived based on necking in thin cracked sheets. With this model, the cohesive behavior including peak cohesive traction, cohesive energy density and shape of the cohesive traction–separation curve is discussed. The peak cohesive traction is found to vary from 1.15 times the yield stress for perfectly plastic materials to about 2.5 times the yield stress for modest hardening materials (power hardening exponent of 0.2). The cohesive energy density depends on the critical relative plate thickness reduction at the root of the neck at crack initiation, which needs to be determined by experiments. Finally, an elastic background medium with a center crack is employed to re-examine the shape effect of cohesive traction–separation curve, and the relation between the linear elastic fracture mechanics (LEFM) and cohesive zone models by considering the cohesive zone development and crack growth in the infinite elastic medium. It is shown that the shape of the cohesive curve does affect the cohesive zone size and the apparent energy release rate of LEFM for the crack growth in the elastic background material. The apparent energy release rate of LEFM approaches the cohesive energy density when the crack extends significantly longer than the characteristic length of the cohesive zone.     

Fatigue crack closure study based on whole-field displacements

This study is concerned with crack tip strain field fluctuations at loads below the point of crack closure in fatigue cycling. Moiré interferometry was used to investigate crack tip fields in compact tension specimens, cracked under constant stress intensity range and fixed R-ratio conditions. An elastic-plastic finite element model of simulated closure was developed to provide a theoretical cross-reference for the moiré studies. The ‘stretched zone’, which is believed to be the most significant source of closure effects, was simulated by inserting a constant thickness strip of elements into the crack before unloading from the maximum load point. Analysis of the crack tip fields in the experimental and theoretical cases was made in terms of crack face opening profiles, compliance changes and elastic stress intensity parameters. The latter were inferred through stress and displacement measurements made along circular and radial paths relative to the crack tip. Closure on the stretched zone was found to generate non-proportional loading in the crack tip field, so that the resulting stress changes were not well characterized by the asymptotic elastic equations. It is concluded firstly, that significant strain fluctuations occur below the point of closure load and that these should not be ignored in crack propagation studies. Secondly, the effective stress intensity range in fatigue cycling is not simply related to the open-crack stress intensity range and the need therefore remains for R-ratio and geometry effects to be treated as variables in crack propagation data collection programmes.     

A continuum mechanics approach to some problems in subcritical crack propagation

The results of the so-called energetic approach to fracture for the cases of a sharp crack without and with a cohesive zone are briefly reviewed with particular attention to the crack tip singularity analysis and to the issue of energy dissipation due to crack propagation. The case of a crack with a cohesive zone removing all thermomechanical singularities is then further analyzed, focusing the attention on the question of the thermodynamic admissibility of subcritical crack growth, and on some of the hypotheses that lead to the derivation of subcritical crack growth laws. A two-phase cohesive zone model for discontinuous crack growth is presented and its thermodynamics analyzed, followed by an example of its possible application.     

Computations of fatigue crack growth with strain gradient plasticity and an irreversible cohesive zone model

Computations of fatigue crack growth with a first-order strain gradient plasticity (SGP) model and an irreversible cohesive zone model are reported. SGP plays a significant role in the model predictions and leads to increased fatigue crack growth rates relative to predictions with classical plasticity. Increased magnitudes of tractions and material separation at the crack tip together with reduced crack closure appear as the cause for accelerated crack growth in SGP. Under plane strain conditions SGP appears as an essential feature of the development of the crack closure zone. Size effects are explored relative to changes in internal material length scale as well as to structural length scales.     

Cohesive Crack with Rate-Dependent Opening and Viscoelasticity: I. Mathematical Model and Scaling

The time dependence of fracture has two sources: (1) the viscoelasticity of material behavior in the bulk of the structure, and (2) the rate process of the breakage of bonds in the fracture process zone which causes the softening law for the crack opening to be rate-dependent. The objective of this study is to clarify the differences between these two influences and their role in the size effect on the nominal strength of stucture. Previously developed theories of time-dependent cohesive crack growth in a viscoelastic material with or without aging are extended to a general compliance formulation of the cohesive crack model applicable to structures such as concrete structures, in which the fracture process zone (cohesive zone) is large, i.e., cannot be neglected in comparison to the structure dimensions. To deal with a large process zone interacting with the structure boundaries, a boundary integral formulation of the cohesive crack model in terms of the compliance functions for loads applied anywhere on the crack surfaces is introduced. Since an unopened cohesive crack (crack of zero width) transmits stresses and is equivalent to no crack at all, it is assumed that at the outset there exists such a crack, extending along the entire future crack path (which must be known). Thus it is unnecessary to deal mathematically with a moving crack tip, which keeps the formulation simple because the compliance functions for the surface points of such an imagined preexisting unopened crack do not change as the actual front of the opened part of the cohesive crack advances. First the compliance formulation of the cohesive crack model is generalized for aging viscoelastic material behavior, using the elastic-viscoelastic analog (correspondence principle). The formulation is then enriched by a rate-dependent softening law based on the activation energy theory for the rate process of bond ruptures on the atomic level, which was recently proposed and validated for concrete but is also applicable to polymers, rocks and ceramics, and can be applied to ice if the nonlinear creep of ice is approximated by linear viscoelasticity. Some implications for the characteristic length, scaling and size effect are also discussed. The problems of numerical algorithm, size effect, roles of the different sources of time dependence and rate effect, and experimental verification are left for a subsequent companion paper. This revised version was published online in July 2006 with corrections to the Cover Date.