Analysis of a rate-dependent cohesive model for dynamic crack propagation

The effect of including rate-dependence in the cohesive zone modeling of steady-state and transient dynamic crack propagation is analyzed. Spontaneous crack propagation simulations are performed using a spectral form of the elastodynamic boundary integral equations, while the solution to the steady-state problem is obtained by solving the governing Cauchy singular equation on the crack plane. The steady-state analysis shows that the existing techniques for solving the Cauchy singular integral equation are not suitable. A solution technique for the underlying Riemann-Hilbert problem for the chosen rate and damage-dependent cohesive law is presented. Under spontaneous propagation conditions, quasi-steady-state speeds slower than the theoretically predicted shear wave speed are possible. Results also show that, due to the dissipation of energy inside the cohesive zone, the energy required for crack propagation increases with the crack speed.     

Numerical analysis of dynamic crack propagation in rubber

In the present paper, dynamic crack propagation in rubber is analyzed numerically using the finite element method. The problem of a suddenly initiated crack at the center of stretched sheet is studied under plane stress conditions. A nonlinear finite element analysis using implicit time integration scheme is used. The bulk material behavior is described by finite-viscoelasticity theory and the fracture separation process is characterized using a cohesive zone model with a bilinear traction-separation law. Hence, the numerical model is able to model and predict the different contributions to the fracture toughness, i.e. the surface energy, viscoelastic dissipation, and inertia effects. The separation work per unit area and the strength of the cohesive zone have been parameterized, and their influence on the separation process has been investigated. A steadily propagating crack is obtained and the corresponding crack tip position and velocity history as well as the steady crack propagation velocity are evaluated and compared with experimental data. A minimum threshold stretch of 3.0 is required for crack propagation. The numerical model is able to predict the dynamic crack growth. It appears that the strength and the surface energy vary with the crack speed. Finally, the maximum principal stretch and stress distribution around steadily propagation crack tip suggest that crystallization and cavity formation may take place.     

A semi-analytical solution method for problems of cohesive fracture and some of its applications

Cohesive zone models are extensively used for the failure load estimates for structure elements with cracks. This paper focuses on some features of the models associated with the failure load and size of the cohesive zone predictions. For simplicity, considered is a mode I crack in an infinite plane under symmetrical tensile stresses. A traction–separation law is prescribed in the crack process zone. It is assumed by the problem statement that the crack faces close smoothly. This requirement is satisfied numerically by a formulation of the modified boundary conditions. The critical state of a plate with a cohesive crack is analyzed using singular integral equations. A numerical procedure is proposed to solve the obtained systems of integral equations and inequalities. The presented solution is in agreement with other published results for some limiting cases. Thus, an effective methodology is devised to solve crack mechanics problems within the framework of a cohesive zone model. Using this methodology, some problems are solved to illustrate the (i) influence of shape parameters of traction–separation law on the failure load, (ii) ability to account for contact stress for contacting crack faces, (iii) influence of getting rid of stress finiteness condition in the problem statement.     

Analysis of a rate-dependent cohesive model for dynamic crack propagation

The effect of including rate-dependence in the cohesive zone modeling of steady-state and transient dynamic crack propagation is analyzed. Spontaneous crack propagation simulations are performed using a spectral form of the elastodynamic boundary integral equations, while the solution to the steady-state problem is obtained by solving the governing Cauchy singular equation on the crack plane. The steady-state analysis shows that the existing techniques for solving the Cauchy singular integral equation are not suitable. A solution technique for the underlying Riemann–Hilbert problem for the chosen rate and damage-dependent cohesive law is presented. Under spontaneous propagation conditions, quasi-steady-state speeds slower than the theoretically predicted shear wave speed are possible. Results also show that, due to the dissipation of energy inside the cohesive zone, the energy required for crack propagation increases with the crack speed.     

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.     

Numerical simulations of a dynamically propagating crack with a nonlinear cohesive zone

Numerical solutions of a dynamic crack propagation problem are presented. Specifically, a mode III semi-infinite crack is assumed to be moving in an unbounded homogeneous linear elastic continuum while the crack tip consists of a nonlinear cohesive (or failure) zone. The numerical results are obtained via a novel semi-analytical technique based on complex variables and integral transforms. The relation between the properties of the failure zone and the resulting crack growth regime are investigated for several rate independent as well as rate dependent cohesive zone models. Based on obtained results, an hypothesis is formulated to explain the origin of the crack tip velocity periodic fluctuations that have been detected in recent dynamic crack propagation experiments. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

An extended finite element method with analytical enrichment for cohesive crack modeling

A recent approach to fracture modeling has combined the extended finite element method (XFEM) with cohesive zone models. Most studies have used simplified enrichment functions to represent the strong discontinuity but have lacked an analytical basis to represent the displacement gradients in the vicinity of the cohesive crack. In this study enrichment functions based upon an existing analytical investigation of the cohesive crack problem are proposed. These functions have the potential of representing displacement gradients in the vicinity of the cohesive crack and allow the crack to incrementally advance across each element. Key aspects of the corresponding numerical formulation and enrichment functions are discussed. A parameter study for a simple mode I model problem is presented to evaluate if quasi‐static crack propagation can be accurately followed with the proposed formulation. The effects of mesh refinement and mesh orientation are considered. Propagation of the cohesive zone tip and crack tip, time variation of the cohesive zone length, and crack profiles are examined. The analysis results indicate that the analytically based enrichment functions can accurately track the cohesive crack propagation of a mode I crack independent of mesh orientation. A mixed mode example further demonstrates the potential of the formulation. Copyright © 2008 John Wiley & Sons, Ltd.     

A micromechanical model for a viscoelastic cohesive zone

A micromechanical model for a viscoelastic cohesive zone is formulated herein. Care has been taken in the construction of a physically-based continuum mechanics model of the damaged region ahead of the crack tip. The homogenization of the cohesive forces encountered in this region results in a damage dependent traction-displacement law which is both single integral and internal variable-type. An incrementalized form of this traction-displacement law has been integrated numerically and placed within an implicit finite element program designed to predict crack propagation in viscoelastic media. This research concludes with several example problems on the response of this model for various displacement boundary conditions.     

Experimental validation of large-scale simulations of dynamic fracture along weak planes

A well-controlled and minimal experimental scheme for dynamic fracture along weak planes is specifically designed for the validation of large-scale simulations using cohesive finite elements. The role of the experiments in the integrated approach is two-fold. On the one hand, careful measurements provide accurate boundary conditions and material parameters for a complete setup of the simulations without free parameters. On the other hand, quantitative performance metrics are provided by the experiments, which are compared a posteriori with the results of the simulations. A modified Hopkinson bar setup in association with notch-face loading is used to obtain controlled loading of the fracture specimens. An inverse problem of cohesive zone modeling is performed to obtain accurate mode-I cohesive zone laws from experimentally measured deformation fields. The speckle interferometry technique is employed to obtain the experimentally measured deformation field. Dynamic photoelasticity in conjunction with high-speed photography is used to capture experimental records of crack propagation. The comparison shows that both the experiments and the numerical simulations result in very similar crack initiation times and produce crack tip velocities which differ by less than 6%. The results also confirm that the detailed shape of the non-linear cohesive zone law has no significant influence on the numerical results.     

Mode I crack in a soft ferromagnetic material

ABSTRACT In the existing magnetoelastic theories, stress is proportional to the square of magnetic intensity and the linear model developed is usually used to analyse magnetoelastic problems. For a crack problem, the perturbation of the magnetic field caused by deformation is not much less than the applied field. In this paper, complex potentials for a mode I crack with a nonlinear relation for magnetic intensity are developed. The boundary conditions on crack faces are represented in terms of the continuity of the magnetic field. A solution of the crack problem is obtained by solving the Riemann‐Hilbert problem. Making use of the solution, the effects of the boundary conditions on the crack faces on the magnetoelastic coupling are discussed.     

A numerical method for elastic and viscoelastic dynamic fracture problems in homogeneous and bimaterial systems

We present a summary of recent advances in the development of an efficient numerical scheme to be used in the investigation of a wide range of 2D and 3D dynamic fracture problems. The numerical scheme, which is based on a spectral representation of the boundary integral relations, can be applied to homogeneous and interfacial dynamic fracture problems involving planar cracks and faults of arbitrary shapes buried in elastic and viscoelastic media. Spontaneous propagation of the crack is achieved by combining the elastodynamic integral relations with a stress-based cohesive failure model. The objective of this paper is to present some of the major differences existing between the various formulations within the simpler 2D scalar framework of anti-plane shear (mode III) loading conditions. Examples are presented to illustrate some capabilities of the method.     

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.     

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.     

Modelling of crack propagation of gravity dams by scaled boundary polygons and cohesive crack model

Crack propagation in concrete gravity dams is investigated using scaled boundary polygons coupled with interface elements. The concrete bulk is assumed to be linear elastic and is modelled by the scaled boundary polygons. The interface elements model the fracture process zone between the crack faces. The cohesive tractions are modelled as side-face tractions in the scaled boundary polygons. The solution of the stress field around the crack tip is expressed semi-analytically as a power series. It reproduces the singular and higher-order terms in an asymptotic solution, such as the William’s eigenfunction expansion when the cohesive tractions vanish. Accurate results can be obtained without asymptotic enrichment or local mesh refinement. The stress intensity factors are obtained directly from their definition and provide a convenient and accurate means to assess the zero-K condition, which determines the stability of a cohesive crack. The direction of crack propagation is determined from the maximum circumferential stress criterion. To accommodate crack propagation, a local remeshing algorithm that is applicable to any polygon mesh is augmented by inserting cohesive interface elements between the crack surfaces as the cracks propagate. Three numerical benchmarks involving crack propagation in concrete gravity dams are modelled. The results are compared to the experimental and other numerical simulations reported in the literature.     

Computational model for predicting nonlinear viscoelastic damage evolution in materials subjected to dynamic loading

Many inelastic solids accumulate numerous cracks before failure due to impact loading, thus rendering any exact solution of the IBVP untenable. It is therefore useful to construct computational models that can accurately predict the evolution of damage during actual impact/dynamic events in order to develop design tools for assessing performance characteristics. This paper presents a computational model for predicting the evolution of cracking in structures subjected to dynamic loading. Fracture is modeled via a nonlinear viscoelastic cohesive zone model. Two example problems are shown: one for model validation through comparison with a one-dimensional analytical solution for dynamic viscoelastic debonding, and the other demonstrates the applicability of the approach to model dynamic fracture propagation in the double cantilever beam test with a viscoelastic cohesive zone.     

Lifetime evaluation of concrete structures under sustained post-peak loading

Experimental tests on crack propagation in concrete under constant post-peak loading are simulated using the finite element method and the cohesive crack model, in both Mode I and Mixed-mode conditions. The time-dependent behaviour of concrete in the process zone is due to the interaction and growth of microcracks, a phenomenon which, for high constant load levels, turns out to be predominant over linear viscoelastic creep in the bulk material. In mechanical systems based on this type of material behaviour (creep and strain-softening taking place simultaneously), the initial value problem is non-parabolic, i.e., the error at one time level is affected by the accumulation of errors introduced at earlier time levels. Despite these difficulties, the scatter in numerical failure lifetime vs. load level turns out to be negligible in Mode I conditions and practically acceptable in Mixed-mode conditions. Therefore the time-dependent behaviour of the process zone can be inferred solely from the results of direct tensile tests.     

Dynamic steady-state crack propagation in a transversely isotropic viscoelastic body

The problem considered herein is the dynamic, subsonic, steady-state propagation of a semi-infinite, generalized plane strain crack in an infinite, transversely isotropic, linear viscoelastic body. The corresponding boundary value problem is considered initially for a general anisotropic, linear viscoelastic body and reduced via transform methods to a matrix Riemann–Hilbert problem. The general problem does not readily yield explicit closed form solutions, so attention is addressed to the special case of a transversely isotropic viscoelastic body whose principal axis of material symmetry is parallel to the crack edge. For this special case, the out-of-plane shear (Mode III), in-plane shear (Mode II) and in-plane opening (Mode I) modes uncouple. Explicit expressions are then constructed for all three Stress Intensity Factors (SIF). The analysis is valid for quite general forms for the relevant viscoelastic relaxation functions subject only to the thermodynamic restriction that work done in closed cycles be non-negative. As a special case, an analytical solution of the Mode I problem for a general isotropic linear viscoelastic material is obtained without the usual assumption of a constant Poissons ratio or exponential decay of the bulk and shear relaxation functions. The Mode I SIF is then calculated for a generalized standard linear solid with unequal mean relaxation times in bulk and shear leading to a non-constant Poissons ratio. Numerical simulations are performed for both point loading on the crack faces and for a uniform traction applied to a compact portion of the crack faces. In both cases, it is observed that the SIF can vanish for crack speeds well below the glassy Rayleigh wave speed. This phenomenon is not seen for Mode I cracks in elastic material or for Mode III cracks in viscoelastic material.     

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.