The influence of the cohesive process zone in hydraulic fracturing modelling

This paper studies the importance of the cohesive zone in the modelling of a fluid driven fracture under plain strain conditions. The fracture is driven by pumping of an incompressible viscous fluid at the fracture inlet. Rock deformation is modeled for linear elastic and poroelastic solids. Fluid flow in the fracture is modeled by lubrication theory. The cohesive zone approach is used as the fracture propagation criterion. Finite element analysis was used to compute the solution for the crack length, the fracture opening and propagation pressure as a function of the time and distance from the wellbore. It is demonstrated that the crack profiles and the propagation pressures are larger in the case of elastic-softening cohesive model compared to the results of the rigid-softening cohesive model for both elastic and poroelastic cohesive solids. It is found that the results are affected by the slope of the loading branch of the cohesive model and they are nearly unaffected from the exact form of the softening branch. Furthermore, the size of the process zone, the fracture geometry and the propagation pressure increase with increasing confining stresses. These results may explain partially the discrepancies in net-pressures between field measurements and conventional model predictions.     

黏聚裂纹阻抗的弯曲梁承载力

在混凝土类软化材料断裂研究中,裂纹端部损伤区被简化为具有黏聚应力分布的非线性裂纹,该黏 聚力对裂纹扩展有阻抗作用。裂纹体的应力强度因子是断裂力学标志载荷作用与几何构型因素的量化表达指标; 黏聚力形成的阻抗强度因子数值,与黏聚裂纹长度和材料极值拉伸应力存在数量关系。通过双K断裂判据,以 带切口的三点弯曲梁为断裂力学模型,分析了裂纹黏聚阻力对断裂过程的影响规律,计算该弯曲梁结构断裂试 样的最大承担载荷;其结果比较符合实验数据。     

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.     

Numerical study of geometric constraint and cohesive parameters in steady-state viscoelastic crack growth

This paper presents a finite element study of cohesive crack growth in a thin infinite viscoelastic strip to investigate the effects of viscoelastic properties, strip height, and cohesive model parameters on the crack growth resistance. The results of the study show that the dependence of the fracture energy on the viscoelastic properties for the strip problem is similar to that obtained for the infinite body problem even when the cohesive zone length is large compared to the height of the strip. The fracture energy also depends on the crack speed v through the dimensionless parameter v τ/L _∞ where L _∞ is the characteristic length of the cohesive zone and τ is the characteristic relaxation time of the bulk material. This relationship confirms that at least two properties of the fracture process must be prescribed accurately to model viscoelastic crack growth. In contrast, the fracture energy and crack speed are insensitive to the strip height even in situations where the growth of the dissipation zone is severely constrained by the strip boundaries. We observe that at high speeds, where the fracture energy asymptotically approaches the maximum value, the material surrounding the cohesive zone is in the rubbery (equilibrium) state and not the glassy state.     

A Viscoelastic Fictitious Crack Model for the Fracture of Sea Ice

The fracture of sea ice is modeled using a viscoelastic fictitious crack(cohesive zone) model. The sea ice is modeled as a linear viscoelastic material. The fictitious crack model is implemented via the weight function method. The associated stress-separation curve can be rate dependent. The impact of assuming viscoelastic behavior in the bulk as opposed to elastic behavior is studied. Results from the model are compared to the available exact results for various test cases. The model is applied to a large scale in situ sea ice fracture test. Various implications of such applications are pointed out. This viscoelastic fictitious crack model is found to be a promising tool in investigations pertaining to the fracture of sea ice. This revised version was published online in July 2006 with corrections to the Cover Date.     

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 boundary integral method for a dynamic, transient mode I crack problem with viscoelastic cohesive zone

We consider the problem of the dynamic, transient propagation of a semi-infinite, mode I crack in an infinite elastic body with a nonlinear, viscoelastic cohesize zone. Our problem formulation includes boundary conditions that preclude crack face interpenetration, in contrast to the usual mode I boundary conditions that assume all unloaded crack faces are stress-free. The nonlinear viscoelastic cohesive zone behavior is motivated by dynamic fracture in brittle polymers in which crack propagation is preceeded by significant crazing in a thin region surrounding the crack tip. We present a combined analytical/numerical solution method that involves reducing the problem to a Dirichlet-to-Neumann map along the crack face plane, resulting in a differo-integral equation relating the displacement and stress along the crack faces and within the cohesive zone.     

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.     

Modeling fracture in nanomaterials via a virtual internal bond method

The recent surging interest in nanotechnology is providing a strong impetus to understanding fracture processes in nanoscale materials. There are open challenges because many classical concepts of fracture mechanics are no longer applicable as the characteristic dimension of a structure becomes comparable to or smaller than the size of the cohesive zone near a crack tip. In materials with a characteristic size on the nanometer scale, the fracture process is often strongly dominated by the surface energy and nonlinear material properties. In this paper, we apply a recently developed virtual-internal-bond (VIB) method to investigating fracture of such nanomaterials. In the VIB method, a cohesive interactive law is directly incorporated into the constitutive model so that separate fracture criteria are no longer necessary. We demonstrate that, at a critical length scale typically on the order of nanometer scale, the fracture mechanism changes from the classical Griffith fracture to one of homogeneous failure near the theoretical strength of solids; when this transition occurs, the classical singular deformation field near a crack tip disappears and is replaced by a uniform stress distribution with no stress concentration near the crack tip.     

The effects of inter-surface cohesive tractions on linear and penny-shaped cracks

A mesoscopic fracture model of equilibrium slit cracks in brittle solids, including inter-surface cohesive tractions acting near the crack tip, is analyzed and the effects of the cohesive tractions on the in-plane stress fields, crack-opening displacement profiles, and crack driving forces examined quantitatively for linear and penny-shaped cracks. The (numerical) analysis method is described in detail, along with results for four different cohesive forces. The assumed distribution of cohesive tractions were found to suppress the in-plane stress field adjacent to cracks in a homogeneous, isotropic medium when uniformly loaded in mode-I, and the suppression was a function of crack length. The crack-opening displacement profile was also perturbed and a new regime identified between the near-field Barenblatt zone and the far-field continuum zone. The extent of this `cohesive zone' was quantified by use of an interpolating function fit to the calculated profiles and found to be independent of crack size for a given cohesive tractions distribution. Furthermore, the crack-opening displacement at the edge of the cohesive zone was found to be independent of crack size, implying that despite significant perturbations to the stress field, the crack driving force at unstable equilibrium remains unchanged with crack size.     

The cohesive zone model in a problem of delayed hydride cracking of zirconium alloys

The crack tip model with the cohesive zone ahead of a finite crack tip has been presented. The estimation of the length of the cohesive zone and the crack tip opening displacement is based on the comparison of the local stress concentration, according to Westergaard's theory, with the cohesive stress. To calculate the cohesive stress, von Mises yield condition at the boundary of the cohesive zone is employed for plane strain and plane stress. The model of the stress distribution with the maximum stress within the cohesive zone is discussed. Local criterion of brittle fracture and modelling of the fracture process zone by cohesive zone were used to describe fracture initiation at the hydride platelet in the process zone ahead of the crack tip. It was shown that the theoretical K _IH-estimation applied to the case of mixed plane condition within the process zone is qualitatively consistent with experimental data for unirradiated Zr-2.5Nb alloy. In the framework of the proposed model, the theoretical value of K ^H _IC for a single hydride platelet at the crack tip has been also estimated.     

The cohesive zone model in a problem of delayed hydride cracking of zirconium alloys

The crack tip model with the cohesive zone ahead of a finite crack tip has been presented. The estimation of the length of the cohesive zone and the crack tip opening displacement is based on the comparison of the local stress concentration, according to Westergaard's theory, with the cohesive stress. To calculate the cohesive stress, von Mises yield condition at the boundary of the cohesive zone is employed for plane strain and plane stress. The model of the stress distribution with the maximum stress within the cohesive zone is discussed. Local criterion of brittle fracture and modelling of the fracture process zone by cohesive zone were used to describe fracture initiation at the hydride platelet in the process zone ahead of the crack tip. It was shown that the theoretical K _IH-estimation applied to the case of mixed plane condition within the process zone is qualitatively consistent with experimental data for unirradiated Zr-2.5Nb alloy. In the framework of the proposed model, the theoretical value of K ^H _IC for a single hydride platelet at the crack tip has been also estimated.     

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.     

Craze and yield zones in the fracture of rigid PVC

Two techniques have been used to examine the reported discrepancy between the measured fracture toughness of rigid PVC and that calculated from its crack tip craze. The use of a generalized model of a cohesive zone produced no better agreement between calculation and experiment than the Dugdale model used previously. The discrepancy was probably caused by the coexistence with the craze of a yield zone which was observed by optical and scanning electron microscopy.     

A combined dislocation—cohesive zone model for fracture in a confined ductile layer

The collective dislocation behavior near a crack tip in a ductile layer sandwiched between two brittle solids is analyzed via two-dimensional dislocation dynamics (DD) simulations that incorporate a cohesive zone (CZ) model. The cohesive crack tip is treated as part of a much larger finite crack confined in the ductile layer. The underlying boundary value problem is formulated with a set of boundary integral equations and numerically evaluated with a collocation method. The fracture energy of the layered composite material is shown to be strongly correlated with the layer thickness and is directly influenced by the cohesive strength of the ductile layer (Hsia KJ et al. (1994) J Mech Phys Solids 6 877–896).     

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.     

A bilinear cohesive zone model tailored for fracture of asphalt concrete considering viscoelastic bulk material

A bilinear cohesive zone model (CZM) is employed in conjunction with a viscoelastic bulk (background) material to investigate fracture behavior of asphalt concrete. An attractive feature of the bilinear CZM is a potential reduction of artificial compliance inherent in the intrinsic CZM. In this study, finite material strength and cohesive fracture energy, which are cohesive parameters, are obtained from laboratory experiments. Finite element implementation of the CZM is accomplished by means of a user-subroutine which is employed in a commercial finite element code (e.g., UEL in ABAQUS). The cohesive parameters are calibrated by simulation of mode I disk-shaped compact tension results. The ability to simulate mixed-mode fracture is demonstrated. The single-edge notched beam test is simulated where cohesive elements are inserted over an area to allow cracks to propagate in any general direction. The predicted mixed-mode crack trajectory is found to be in close agreement with experimental results. Furthermore, various aspects of CZMs and fracture behavior in asphalt concrete are discussed including: compliance, convergence, and energy balance.     

A nonlinear viscoelastic fracture analysis of concrete/FRP delamination in aggressive environments

In this paper, the durability of the bondline between concrete and FRP reinforcement was characterized at various temperature and humidity levels. The linear and nonlinear viscoelastic constitutive behavior of the epoxy bondline was characterized and used for a nonlinear viscoelastic fracture analysis of delamination. A hygrothermal nonlinear viscoelastic pseudo-stress model was developed and calibrated in order to compute a generalized J integral. Driven wedge tests were conducted for examining the fracture behavior of the interface. A finite element analysis was developed for determining the cohesive zone size and the generalized J integral at various temperature and humidity levels. The fracture energy obtained from these parameters greatly depended upon crack growth rate, temperature and humidity.     

Time-dependent cohesive zone modelling for discrete fracture simulation

The cohesive crack tip model became very popular in fracture and failure mechanics, starting with the original publications of (Dugdale, 1960) [1] and (Barenblatt, 1962) [2] and first practical applications in the early 1980s. A compact representation of the fracture mechanical basis, kinematic and constitutive issues as well as some special characteristics of the related finite element formulation are given in the first part of this contribution. The main part is dedicated to the presentation and discussion of recent developments of cohesive finite element methods for time-dependent fracture. On the basis of a rheological model assumption, a novel viscoelastic extension for cohesive traction separation laws is presented and the resultant characteristic behaviour is depicted and compared for different loading conditions. Adopting an industrial application of a peel foil specimen, the time-dependent characteristics as well as some aspects of parameter identification and application of the material model are shown.