A formulation is presented to perform crack propagation analyses in cohesive materials with the dual boundary element method (DBEM) using the tangential differential operator in the traction boundaryintegral equations. The cohesive law is introduced in the system of equations to directly compute the cohesive forces at each loading step. A single edge crack is analyzed with the linear function to describe the material softening law in the cohesive zone, and the results are compared with those from the literature. Keywords: cohesive model, crack analysis, dual boundary element model, plane problems, tangential differential operator.
INTRODUCTIONThe surfaces of cracks behind the (fictitious) crack tip are not completely separated in some materials, such as concrete, brittle polymers, fiber-reinforced composites, tough ceramics and some alloys. The tractions can be transferred across the crack line along a relatively long extension of the crack, which is commonly called the wake zone, the bridging zone, or the cohesive zone. The main assumption is the occurrence of material softening beyond the peak load in a narrow layer with a negligible volume behind the fictitious crack tip, in which the action can be replaced by cohesive forces (Fig. 1). Two types of constitutive laws for the material in the cohesive zone have been used in the literature. This study employs a traction-displacement relationship in the cohesive zone instead of using a material constitutive law that is defined in terms of stresses and strains accompanied by a layer thickening law. Barenblatt [1] introduced a cohesive model that employed a fictitious crack. Hillerborg et al.[2] proposed a function for a softening model related to an opening (mode I) crack (Fig. 2), which allowed finite element analyses of the problem to be performed, such as those by Petersson [3], Carpinteri [4], and Rots [5].The boundary element method (BEM) is superior to the finite element method (FEM) in crack propagation analysis because remeshing of the domain is not necessary when the crack grows, and new elements can be introduced without affecting the existing elements. The coincidence of two crack surfaces requires two different boundary integral equations (BIEs) for the solution. The dual boundary element method (DBEM) is a technique that is employed for crack propagation analysis. The position of the collocation point to solve the traction BIE and the strategy to treat improper integrals are the essential features of the formulation. The tangential differential operator (TDO) in conjunction with integration by parts is an interesting procedure for reducing the order of singularity in fundamental solution kernels of traction BIEs when Kelvin-type fundamental solutions are employed. Kupradze [7] first presented an application using the TDO, and Sladek [8] employed the TDO in a curved crack solution. Bonnet [9] presented regularized formulations for the BEM employing the TDO for gradients in potential problems and in stress BIEs for elasticity problems, includin...