Modelling of growth of three-dimensional cracks by a continuous distribution of dislocation loops

This paper is concerned with a numerical simulation of growth of fatigue cracks in a three-dimensional geometry. A continuous distribution of infinitesimal dislocation loops is employed to model the crack faces, so that the crack problem can be formulated as a set of singular integral equations. A numerical procedure based on an analytical treatment of the associated finite part integral is developed to solve the singular integral equations. The Paris law is then used to predict the rate of crack growth, so that the evolution of the crack shape under fatigue can be traced by a step-by-step algorithm. Various crack growth problems, e.g., the growth of a subsurface crack in a surface treated specimen, are analyzed using the technique, providing new data for several cracks of practical interest.     

A distributed dislocation stress analysis for crazes and plastic zones at crack tips

A distributed dislocation method is developed to obtain analytically the applied stress as well as the surface stress profile along narrow plastic zones at the tip of a crack in a homogeneous tensile stress field. Replacing the plastic zone by a continuous array of mathematical dislocations, the stress field solution of this mixed boundary value problem (the displacement profile of the plastic zone is fixed while the tensile stresses are zero across the crack) can be solved. A computer program based on this stress field solution has been constructed and tested using the analytical results of the Dugdale model. The method is then applied to determining the surface stress profiles of crazes and plane-stress plastic deformation zones grown from electron microprobe cracks in polystyrene and polycarbonate respectively. The necessary craze and zone surface displacement profiles are determined by quantitative analysis of transmission electron micrographs. The surface stress profiles, which show small stress concentrations at the craze or zone tip falling to an approximately constant value which is maintained to the crack tip, are compared with those previously computed using an approximate Fourier transform method involving estimation of the displacement profile in the crack. The agreement between the approximate method and the exact distributed dislocation method is satisfactory.     

Application of the distributed dislocation technique for calculating cyclic crack tip plasticity effects

This paper describes a method for modelling cyclic crack tip plasticity effects based on the distributed dislocation technique (DDT). A strip‐yield model is utilised to allow for the determination of the crack opening displacement, size of the plastic zones and in the case of a fatigue crack, the wake of plasticity. The DDT can be easily implemented for a wide range of cracked geometries with reliable control over the accuracy and convergence. Thickness effects can also be incorporated through a recently obtained solution for an edge dislocation in an infinite plate of finite thickness. Results for finite length cracks that have had limited growth, such that there is no plastic wake, are presented for a range of applied loads and R‐ratios. Further results are provided for a steady‐state fatigue crack in a plate of finite thickness. The present results are compared with analytical solutions and they show an excellent agreement.     

Numerical analysis of 2-D crack propagation problems using the numerical manifold method

The numerical manifold method is a cover-based method using mathematical covers that are independent of the physical domain. As the unknowns are defined on individual physical covers, the numerical manifold method is very suitable for modeling discontinuities. This paper focuses on modeling complex crack propagation problems containing multiple or branched cracks. The displacement discontinuity across crack surface is modeled by independent cover functions over different physical covers, while additional functions, extracted from the asymptotic near tip field, are incorporated into cover functions of singular physical covers to reflect the stress singularity around the crack tips. In evaluating the element matrices, Gaussian quadrature is used over the sub-triangles of the element, replacing the simplex integration over the whole element. First, the method is validated by evaluating the fracture parameters in two examples involving stationary cracks. The results show good agreement with the reference solutions available. Next, three crack propagation problems involving multiple and branched cracks are simulated. It is found that when the crack growth increment is taken to be 0.5h≤da≤0.75h, the crack growth paths converge consistently and are satisfactory.     

Extended finite element method for cohesive crack growth

The extended finite element method allows one to model displacement discontinuities which do not conform to interelement surfaces. This method is applied to modeling growth of arbitrary cohesive cracks. The growth of the cohesive zone is governed by requiring the stress intensity factors at the tip of the cohesive zone to vanish. This energetic approach avoids the evaluation of stresses at the mathematical tip of the crack. The effectiveness of the proposed approach is demonstrated by simulations of cohesive crack growth in concrete.     

Shielding and amplification of a penny-shape crack due to the presence of dislocations

The interaction between a penny-shape crack and a dislocation in crystalline materials is investigated within the framework of dislocation dynamics. The long-range and singular stress field resulting from the crack is determined by modeling the crack as continuous distribution of dislocation loops. This distribution is determined by satisfying the traction boundary condition at the crack face, resulting into a singular integral equation of the first kind that is solved numerically. This crack model is integrated with the dislocation dynamics simulation technique to yield the stress field of the combine system of crack and different types of dislocations situated at different positions in a three dimensional space. The integrated system is then used to investigate the dislocation behavior and its influence on the crack opening displacement and the characteristic of the stress field near the crack tip. It is shown that, depending on the relative position of the dislocation and its character, the dislocation may result in reduction in the stress amplitude at the crack tip and in some cases in closure of the crack tip. These analyses yield shielding and amplification zones near the crack providing an insight of the dislocation influence on the crack. The full dislocation dynamic analysis reveals the nature of the crack dislocation interaction and the manner in which the dislocation morphology changes as it is attracted to the crack surfaces, as well as the changes it causes to the crack profile.     

Acoustic emission during jumps in subcritical growth of cracks in three-dimensional bodies

A model for jumps in subcritical growth of plane cracks (of arbitrary shape) is proposed for three-dimensional bodies and the accompanying acoustic emission (AE). An approximation method for determining the displacement field resulting from one jump of the crack, based on the fact that the required solution is found in the comparison of displacement fields obtained from currently known problems in crack theory, is proposed. A series of new analytical relationships between parameters of plane cracks in three-dimensional bodies and parameters of AE is given.     

On the computation of two-dimensional stress intensity factors using the boundary element method

General two-dimensional linear elastic fracture problems are investigated using the boundary element method. The √r displacement and 1/√r traction behaviour near a crack tip are incorporated in special crack elements. Stress intensity factors of both modes I and II are obtained directly from crack-tip nodal values for a variety of crack problems, including straight and curved cracks in finite and infinite bodies. A multidomain approach is adopted to treat cracks in an infinite body. The body is subdivided into two regions: an infinite part with a finite hole and a finite inclusion. Numerical results, compared with exact solution whenever possible, are accurate even with a coarse discretization.     

Modelling cracks in finite bodies by distributed dislocation dipoles

The method of continuous distribution of dislocations is extended here to model cracks in finite geometries. The cracks themselves are still modelled by distributed dislocations, whereas the finite boundaries are represented by a continuous distribution of dislocation dipoles. The use of dislocation dipoles, instead of dislocations, provides a unified formulation to treat both simple and arbitrary boundaries in a numerical solution. The method gives a set of singular integral equations with Cauchy kernels, which can be readily solved using Gauss–Chebyshev quadratures for finite bodies of simple shapes. When applied to arbitrary geometries, the continuous distribution of infinitesimal dislocation dipoles is approximated by a discrete distribution of finite dislocation dipoles. Both the stress intensity factor and the T -stress are evaluated for some well-known crack problems, in an attempt to assess the performance of the methods and to provide some new engineering data.     

THE TIME-DOMAIN DBEM FOR RAPIDLY GROWING CRACKS

The time-domain Dual Boundary Element Method (DBEM) is developed to analyse rapidly growing cracks in structures subjected to dynamic loads. Two-dimensional problems, where the velocity of the crack growth is constant and the path is not predefined are studied. The present method uses the dual boundary formulation, i.e. the displacement and the traction boundary integral equations to obtain the solution by discretizing the boundary of the body and the crack surfaces only. The crack growth is modelled by adding new elements ahead of the crack tip. It is assumed that the direction of the increment is perpendicular to the direction of maximum circumferential stress. The method is used to analyse growing cracks in infinite sheet and finite plates, and the results are compared with other reported solutions, showing good agreement. © 1997 by John Wiley & Sons, Ltd.     

A two-scale approach for the analysis of propagating three-dimensional fractures

This paper presents a generalized finite element method (GFEM) for crack growth simulations based on a two-scale decomposition of the solution—a smooth coarse-scale component and a singular fine-scale component. The smooth component is approximated by discretizations defined on coarse finite element meshes. The fine-scale component is approximated by the solution of local problems defined in neighborhoods of cracks. Boundary conditions for the local problems are provided by the available solution at a crack growth step. The methodology enables accurate modeling of 3-D propagating cracks on meshes with elements that are orders of magnitude larger than those required by the FEM. The coarse-scale mesh remains unchanged during the simulation. This, combined with the hierarchical nature of GFEM shape functions, allows the recycling of the factorization of the global stiffness matrix during a crack growth simulation. Numerical examples demonstrating the approximating properties of the proposed enrichment functions and the computational performance of the methodology are presented.     

A study of the influence of grain boundaries on short crack growth during varying load using a dislocation technique

The propagation of short cracks in the neighbourhood of grain boundaries have been investigated using a technique were the crack is modelled by distributed dislocation dipoles and the plastic deformation is represented by discrete dislocations. Discrete dislocations are emitted from the crack tip as the crack grows. Dislocations can also nucleate at the grain boundaries. The influence on crack growth characteristics of the distance between the initial crack tip and the grain boundary has been studied. It was found that crack growth rate is strongly correlated to the dislocation pile-ups at the grain boundaries.     

基于X-SBFEM的裂纹体非网格重剖分耦合模型研究

扩展有限元法利用了非网格重剖分技术,但需要基于裂尖解析解构造复杂的插值基函数,计算精度受网格疏密和插值基函数等因素影响。比例边界有限元法则在求解无限域和裂尖奇异性问题优势明显,两者衔接于有限元法理论内,可建立一种结合二者优势的断裂耦合数值模型。该文从虚功原理出发,利用位移协调与力平衡机制,提出了一种断裂计算的新方法X-SBFEM,达到了扩展有限元模拟裂纹主体、比例边界有限元模拟裂尖的目的。在数值算例中,通过边裂纹和混合型裂纹的应力强度因子计算,并与理论解对比,验证了该方法的准确性和有效性。     

Crack tip fields in ductile crystals

Results on the asymptotic analysis of crack tip fields in elastic-plastic single crystals are presented and some preliminary results of finite element solutions for cracked solids of this type are summarized. In the cases studied, involving plane strain tensile and anti-plane shear cracks in ideally plastic f c c and b c c crystals, analyzed within conventional small displacement gradient assumptions, the asymptotic analyses reveal striking discontinuous fields at the crack tip.For the stationary crack the stress state is found to be locally uniform in each of a family of angular sectors at the crack tip, but to jump discontinuously at sector boundaries, which are also the surfaces of shear discontinuities in the displacement field. For the quasi-statically growing crack the stress state is fully continuous from one near-tip angular sector to the next, but now some of the sectors involve elastic unloading from, and reloading to, a yielded state, and shear discontinuities of the velocity field develop at sector boundaries. In an anti-plane case studied, inclusion of inertial terms for (dynamically) growing cracks restores a discontinuous stress field at the tip which moves through the material as an elastic-plastic shock wave. For high symmetry crack orientations relative to the crystal, the discontinuity surfaces are sometimes coincident with the active crystal slip planes, but as often lie perpendicular to the family of active slip planes so that the discontinuities correspond to a kinking mode of shear.The finite element studies so far attempted, simulating the ideally plastic material model in a small displacement gradient type program, appear to be consistent with the asymptotic analyses. Small scale yielding solutions confirm the expected discontinuities, within limits of mesh resolution, of displacement for a stationary crack and of velocity for quasi-static growth. Further, the discontinuities apparently extend well into the near-tip plastic zone. A finite element formulation suitable for arbitrary deformation has been used to solve for the plane strain tension of a Taylor-hardening crystal panel containing, a center crack with an initially rounded tip. This shows effects due to lattice rotation, which distinguishes the regular versus kinking shear modes of crack tip relaxation. and holds promise for exploring the mechanics of crack opening at the tip.     

Energy partitioning for a crack under remote shear and compression

The true nature and characteristics of crack growth mechanisms in geologic materials have not been adequately described and are poorly understood. The process by which deformation energy is converted to slipping and growing cracks under compressive stresses is complex and difficult to measure. A hybrid technique employing moiré interferometry as an experimental boundary condition to a finite element method (FEM) was employed for through-cracked polycarbonate plates under remote shear and compression. Cohesive end zone and dislocation slip models are used to approximate experimentally observed displacement characteristics. Shear-driven linear elastic fracture mechanics displacement predictions are shown to be inadequate for initial displacement progression. Moiré displacement fields of relative crack face slip reveal a near tip cohesive zone. The pre-slip moiré-FEM stress fields reveal that the maximum crack tip tensile stress occurs at approximately 45 degrees and further infers cohesive zone presence. A J integral formulation uses moiré displacement data and accounts for stored energy along the crack before and after shear driven crack face slip. These energy-partitioning results track the transfer of stored energy along the crack face to the crack tip until the entire crack is actively slipping. These laboratory-scale experiments capture basic mechanical behavior and simulate thousands of years of large-scale geologic feature displacement history in just a few hours.     

线弹性断裂力学问题的扩展自然单元法

为了更加有效地求解线弹性断裂问题,提出了扩展自然单元法。该方法基于单位分解的思想,在自然单元法的位移模式中加入扩展项表征不连续位移场和裂纹尖端奇异场。通过水平集方法确定裂纹面和裂纹尖端区域,并基于虚位移原理推导了平衡方程的离散线性方程。由于自然单元法的形函数满足Kronecker delta函数性质,本质边界条件易于施加。混合模式裂纹的应力强度因子由相互作用能量积分方法计算。数值算例结果表明扩展自然单元法可以方便地求解线弹性断裂力学问题。     

Dislocation-based semi-analytical method for calculating stress intensity factors of cracks: Two-dimensional cases

Determination of the stress intensity factors of cracks is a fundamental issue for assessing the performance safety and predicting the service lifetime of engineering structures. In the present paper, a dislocation-based semi-analytical method is presented by integrating the continuous dislocation model with the finite element method together. Using the superposition principle, a two-dimensional crack problem in a finite elastic body is reduced to the solution of a set of coupled singular integral equations and the calculation of the stress fields of a body which has the same shape as the original one but has no crack. It can easily solve crack problems of structures with arbitrary shape, and the calculated stress intensity factors show almost no dependence upon the finite element mesh. Some representative examples are given to illustrate the efficacy and accuracy of this novel numerical method. Only two-dimensional cases are addressed here, but this method can be extended to three-dimensional problems.     

Cracking node method for dynamic fracture with finite elements

A new method for modeling discrete cracks based on the extended finite element method is described. In the method, the growth of the actual crack is tracked and approximated with contiguous discrete crack segments that lie on finite element nodes and span only two adjacent elements. The method can deal with complicated fracture patterns because it needs no explicit representation of the topology of the actual crack path. A set of effective rules for injection of crack segments is presented so that fracture behavior beginning from arbitrary crack nucleations to macroscopic crack propagation is seamlessly modeled. The effectiveness of the method is demonstrated with several dynamic fracture problems that involve complicated crack patterns such as fragmentation and crack branching. Copyright © 2008 John Wiley & Sons, Ltd.     

Complete Tangent Stiffness for eXtended Finite Element Method by including crack growth parameters

The eXtended Finite Element Method (XFEM) is a useful tool for modeling the growth of discrete cracks in structures made of concrete and other quasi‐brittle and brittle materials. However, in a standard application of XFEM, the tangent stiffness is not complete. This is a result of not including the crack geometry parameters, such as the crack length and the crack direction directly in the virtual work formulation. For efficiency, it is essential to obtain a complete tangent stiffness. A new method in this work is presented to include an incremental form the crack growth parameters on equal terms with the degrees of freedom in the FEM‐equations. The complete tangential stiffness matrix is based on the virtual work together with the constitutive conditions at the crack tip. Introducing the crack growth parameters as direct unknowns, both equilibrium equations and the crack tip criterion can be handled within the same standard nonlinear iterations. This new solution strategy is believed to provide the modeling capabilities to deal with simultaneous growth of several cracks. A cohesive crack modeling is used. The method is applied to a partly cracked XFEM element of linear strain triangle type with the crack length as the unknown crack growth parameter. In this paper, two examples are given. The first example verifies the theory and the implementation. The second example is the benchmark test three point bending test, where the efficiency of the complete tangential behavior is shown. Copyright © 2013 John Wiley & Sons, Ltd.     

Extended meshfree methods without branch enrichment for cohesive cracks

An extended meshless method for both static and dynamic cohesive cracks is proposed. This new method does not need any crack tip enrichment to guarantee that the crack closes at the tip. All cracked domains of influence are enriched by only the sign function. The domain of influence which includes a crack tip is modified so that the crack tip is always positioned at its edge. The modification is only applied for the discontinuous displacement field and the continuous field is kept unchanged. In addition to the new method, the use of Lagrange multiplier is explored to achieve the same goal. The crack is extended beyond the actual crack tip so that the domains of influence containing the crack tip are completely cut. It is enforced that the crack opening displacement vanishes along the extension of the crack. These methods are successfully applied to several well-known static and dynamic problems.