Constitutive models for cohesive zones in mixed-mode fracture of plain concrete

Discrete mixed-mode fracture (modes I and II) of plain concrete is investigated using a coupled and an uncoupled cohesive zone constitutive model in a finite element context. Fracture surfaces are confined to inter-element boundaries that are not necessarily coincident with the actual fracture surfaces. For this reason, traction components on the cohesive zone do not correspond to actual values either. In this work is demonstrated that only the coupled model is able to cope with these spurious traction components, that must decrease with crack opening. It is shown also that, in this regard, the key variable is the plastic potential adopted in the integration of tractions. Three mixed-mode fracture examples were tested in this work: a three-point single-edge notched beam, double-edge notched plates under variable lateral and normal deformation and four-point double-edge notched beams. A good fitting with experiments was obtained only for the coupled model. Mode II parameters can change in a large range without noticeable change in results, at least in the tested examples.     

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.     

Microcrack nucleation in thermal barrier coating systems

Crack nucleation in thermal-barrier coating (TBC) systems subjected to a monotonic cooling process is studied. The TBC system is modeled using the finite element method, where cracks are represented as discrete discontinuities across continuum elements using the partition-of-unity method. The numerical implementation used for crack nucleation is based on an algorithm where, at insertion of a discontinuity, the traction response is derived from a cohesive zone model that has been modified to (i) behave like an initially rigid cohesive model, and to (ii) ensure smoothness of the traction-separation law at zero crack opening. Accordingly, an adequate convergence behavior of the numerical formulation can be warranted in boundary value problems of systems with relatively complex geometries. In the present numerical study, a comparison is made between TBC systems composed of different constitutive models. The fracture patterns and evolutions of the overall crack growth resulting from the simulations clearly illustrate the importance of accounting for the effects of plasticity in the bond coating and anisotropy in the top coating. The computed fracture profile is in good correspondence with experimental observations reported in the literature.     

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.     

Applications of normal stress dominated cohesive zone models for mixed-mode crack simulation based on extended finite element methods

In conventional cohesive zone models the traction-separation law starts from zero load, so that the model cannot be applied to predict mixed-mode cracking. In the present work the cohesive zone model with a threshold is introduced and applied for simulating different mixed-mode cracks in combining with the extended finite element method. Computational results of cracked specimens show that the crack initiation and propagation under mixed-mode loading conditions can be characterized by the cohesive zone model for normal stress failure. The contribution of the shear stress is negligible. The maximum principal stress predicts crack direction accurately. Computations based on XFEM agree with known experiments very well. The shear stress becomes, however, important for uncracked specimens to catch the correct crack initiation angle. To study mixed-mode cracks one has to introduce a threshold into the cohesive law and to implement the new cohesive zone based on the fracture criterion. In monotonic loading cases it can be easily realized in the extended finite element formulation. For cyclic loading cases convergence of the inelastic computations can be critical.     

An embedded atom hyperelastic constitutive model and multiscale cohesive finite element method

Based on embedded atom method (EAM), an embedded atom hyperelastic (EAH) constitutive model is developed. The proposed EAH constitutive model provides a multiscale formalism to determine mesoscale or macroscale material behavior by atomistic information. By combining the EAH with cohesive zone model (CZM), a multiscale embedded atom cohesive finite element model (EA-cohesive FEM) is developed for simulating failure of materials at mesoscale and macroscale, e.g. fracture and crack propagation etc. Based on EAH, the EA-cohesive FEM applies the Cauchy-Born rule to calculate mesoscale or macroscale material response for bulk elements. Within the cohesive zone, a generalized Cauchy-Born rule is applied to find the effective normal and tangential traction-separation cohesive laws of EAH material. Since the EAM is a realistic semi-empirical interatomic potential formalism, the EAH constitutive model and the EA-cohesive FEM are physically meaningful when it is compared with experimental data. The proposed EA-cohesive FEM is validated by comparing the simulation results with the results of large scale molecular dynamics simulation. Simulation result of dynamic crack propagation is presented to demonstrate the capacity of EA-cohesive FEM in capturing the dynamic fracture.     

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.     

Finite‐deformation irreversible cohesive elements for three‐dimensional crack‐propagation analysis

We develop a three‐dimensional finite‐deformation cohesive element and a class of irreversible cohesive laws which enable the accurate and efficient tracking of dynamically growing cracks. The cohesive element governs the separation of the crack flanks in accordance with an irreversible cohesive law, eventually leading to the formation of free surfaces, and is compatible with a conventional finite element discretization of the bulk material. The versatility and predictive ability of the method is demonstrated through the simulation of a drop‐weight dynamic fracture test similar to those reported by Zehnder and Rosakis. The ability of the method to approximate the experimentally observed crack‐tip trajectory is particularly noteworthy. © 1999 John Wiley & Sons, Ltd.     

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.     

Dynamic ductile fracture in aluminum round bars: experiments and simulations

The application of rate-dependent cohesive elements is validated in simulation of ductile fracture in aluminum round bars under dynamic loading conditions. Smooth and notched round bars made of AA6060-T6 are tested and simulated under quasi-static and dynamic loadings. The smooth round bar is modeled using finite elements that obey Gurson–Tvergaard–Needleman (GTN) formulation as the constitutive equation. Comparing with experimental results, corresponding GTN parameters and rate-dependent plasticity of the alloy are obtained. A single strain rate-dependent GTN element with the obtained parameters is examined under different values of stress triaxiality and loading rates. The resulting stress-elongation curves represent the traction separation law (TSL) for cohesive elements and the variations of the maximum traction and the energy absorbed are investigated. The notched round bars are modeled by axisymmetric continuum and cohesive elements. The undamaged bulk material is elastic-visco plastic and the cohesive elements obey the TSL defined from the single element calculations. The experiments are simulated by these models in which the cohesive elements are rate sensitive and automatically obtain the values of the total strain rate from their adjacent continuum elements to update the values of the cohesive strength during the analysis. The results of the analysis, including maximum load, time of failure and diameter reduction are validated with the experimental results. The effects of element size, rate-dependent plasticity of the material and stress triaxiality are also discussed.     

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 cohesive finite element formulation for modelling fracture and delamination in solids

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

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

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

Modeling of hydrogen-assisted ductile crack propagation in metals and alloys

This paper presents a finite element study of the hydrogen effect on ductile crack propagation in metals and alloys by linking effects at the microstructural level (i.e., void growth and coalescence) to effects at the macro-level (i.e., bulk material deformation around a macroscopic crack). The purpose is to devise a mechanics methodology to simulate the conditions under which hydrogen enhanced plasticity induces fracture that macroscopically appears to be brittle. The hydrogen effect on enhanced dislocation mobility is described by a phenomenological constitutive relation in which the local flow stress is taken as a decreasing function of the hydrogen concentration which is determined in equilibrium with local stress and plastic strain. Crack propagation is modeled by cohesive elements whose traction separation law is determined through void cell calculations that address the hydrogen effect on void growth and coalescence. Numerical results for the A533B pressure vessel steel indicate that hydrogen, by accelerating void growth and coalescence, promotes crack propagation by linking simultaneously a finite number of voids with the crack tip. This “multiple-void” fracture mechanism knocks down the initiation fracture toughness of the material and diminishes the tearing resistance to crack propagation.     

An irreversible cohesive zone model for interface fatigue crack growth simulation

Fatigue crack growth (FCG) along an interface is studied. Instead of using the Paris equation, the actual process of material separation during FCG is described by the use of an irreversible constitutive equation for the cyclic interface traction-separation behavior within the cohesive zone model (CZM) approach. In contrast to past development of CZMs, the traction-separation behavior does not follows a predefined path. The model definition, its predicted cyclic material separation behavior and application to a numerical study of interface FCG in double-cantilever beam, end-loaded split and mixed-mode beam specimens are reported.     

Formulation of a cohesive zone model for a Mode III crack

Cohesive zone models have been proven effective in modeling crack initiation and propagation phenomena. In this work, a possible form for a Mode III cohesive zone model is formulated from elastic stress and displacement fields around a crack with a cohesive zone ahead of the crack tip. A traction-separation relation for the model is derived as a direct consequence of the formulation, which establishes some intrinsic connections between properties of the cohesive zone and those of the bulk material. Interestingly, this model states that the von Mises effective stress in the cohesive zone is constant, which may be related to the bulk material’s yield stress and is consistent with the assumption made in conventional strip-yield elastic-plastic solutions.     

Cohesive fracture modeling of elastic-plastic crack growth in functionally graded materials

This work investigates elastic-plastic crack growth in ceramic/metal functionally graded materials (FGMs). The study employs a phenomenological, cohesive zone model proposed by the authors and simulates crack growth by the gradual degradation of cohesive surfaces ahead of the crack front. The cohesive zone model uses six material-dependent parameters (the cohesive energy densities and the peak cohesive tractions of the ceramic and metal phases, respectively, and two cohesive gradation parameters) to describe the constitutive response of the material in the cohesive zone. A volume fraction based, elastic-plastic model (extension of the original Tamura-Tomota-Ozawa model) describes the elastic-plastic response of the bulk background material. The numerical analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates solid elements with graded elastic and plastic properties and interface-cohesive elements coupled with the functionally graded cohesive zone model. Numerical values of volume fractions for the constituents specified at nodes of the finite element model set the spatial gradation of material properties with isoparametric interpolations inside interface elements and background solid elements to define pointwise material property values. The paper describes applications of the cohesive zone model and the computational scheme to analyze crack growth in a single-edge notch bend, SE(B), specimen made of a TiB/Ti FGM. Cohesive parameters are calibrated using the experimentally measured load versus average crack extension (across the thickness) responses of both Ti metal and TiB/Ti FGM SE(B) specimens. The numerical results show that with the calibrated cohesive gradation parameters for the TiB/Ti system, the load to cause crack extension in the FGM is much smaller than that for the metal. However, the crack initiation load for the TiB/Ti FGM with reduced cohesive gradation parameters (which may be achieved under different manufacturing conditions) could compare to that for the metal. Crack growth responses vary strongly with values of the exponent describing the volume fraction profile for the metal. The investigation also shows significant crack tunneling in the Ti metal SE(B) specimen. For the TiB/Ti FGM system, however, crack tunneling is pronounced only for a metal-rich specimen with relatively smaller cohesive gradation parameter for the metal.     

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.     

Bounds for element size in a variable stiffness cohesive finite element model

The cohesive finite element method (CFEM) allows explicit modelling of fracture processes. One form of CFEM models integrates cohesive surfaces along all finite element boundaries, facilitating the explicit resolution of arbitrary fracture paths and fracture patterns. This framework also permits explicit account of arbitrary microstructures with multiple length scales, allowing the effects of material heterogeneity, phase morphology, phase size and phase distribution to be quantified. However, use of this form of CFEM with cohesive traction–separation laws with finite initial stiffness imposes two competing requirements on the finite element size. On one hand, an upper bound is needed to ensure that fields within crack‐tip cohesive zones are accurately described. On the other hand, a lower bound is also required to ensure that the discrete model closely approximates the physical problem at hand. Both issues are analysed in this paper within the context of fracture in multi‐phase composite microstructures and a variable stiffness bilinear cohesive model. The resulting criterion for solution convergence is given for meshes with uniform, cross‐triangle elements. A series of calculations is carried out to illustrate the issues discussed and to verify the criterion given. These simulations concern dynamic crack growth in an Al₂O₃ ceramic and in an Al₂O₃/TiB₂ ceramic composite whose phases are modelled as being hyperelastic in constitutive behaviour. Copyright © 2004 John Wiley & Sons, Ltd.     

含剪胀效应界面张力-位移关系的一种构造方法

为了研究异种材料界面的开裂过程,在Srensen等人工作的基础上,给出一种含剪胀效应张力-位移关系的构造方法。在界面受拉状态下,通过预先给定的剪胀函数及法向张力-位移关系导出切向张力-位移关系;在界面受压状态下,将切向张力分解为粘结力和摩擦力,摩擦力的大小与法向压力和粘结界面的破坏程度相关。该方法的结果解释了Srensen模型中切向张力-位移关系不连续及其不符合一致关联准则的原因。为便于进行数值计算,给出了用于三维有限元模型的界面刚度矩阵计算方法。选取了一种特定形式的剪胀函数,并将法向张力-位移关系假定为分段线性形式和指数形式,分别求得对应的切向张力-位移关系。最后给出了两个工程应用的例子,数值模拟结果与试验数据吻合较好。