Modeling of cohesive crack growth using an adaptive mesh refinement via the modified-SPR technique

In this paper, an adaptive finite element procedure is presented in modeling of mixed-mode cohesive crack propagation via the modified superconvergent path recovery technique. The adaptive mesh refinement is performed based on the Zienkiewicz–Zhu error estimator. The weighted-SPR recovery technique is employed to improve the accuracy of error estimation. The Espinosa–Zavattieri bilinear cohesive zone model is applied to implement the traction-separation law. It is worth mentioning that no previous information is necessary for the path of crack growth and no region of the domain is necessary to be filled by the cohesive elements. The maximum principal stress criterion is employed for predicting the direction of extension of the cohesive crack in order to implement the cohesive elements. Several numerical examples are analyzed numerically to demonstrate the capability and efficiency of proposed computational algorithm.     

A thermomechanical fracture modeling and simulation for functionally graded solids using a residual-strain formulation

This paper addresses mixed-mode crack growth in two-dimensional functionally graded solids under thermomechanical loads, and investigates the effect of mechanical and thermal loads as well as the T-stress on their crack growth behavior. A novel residual strain-based formulation in the interaction integral method is developed and used for the accurate evaluation of mixed-mode stress intensity factors and/or the T-stress. Simulation of mixed-mode crack propagation in functionally graded materials including solid oxide fuel cells under thermomechanical loads is performed by means of the finite element method and the generalized interaction integrals in conjunction with a remeshing algorithm. An iterative procedure is used for crack growth simulation including the calculation of mixed-mode stress intensity factors and/or the T-stress by means of the generalized interaction integral method, determination of crack growth direction and crack initiation condition based on selected fracture criteria, and local automatic remeshing along the crack path. The present approach employs a user-defined crack increment at the beginning of the simulation. Crack trajectories and fracture parameters obtained by the present simulation for thermomechanical loads are assessed for some numerical examples in comparison with those for mechanical loads.     

Simulation of Crack Propagation in Functionally Graded Materials Under Mixed-Mode and Non-Proportional Loading

Automatic simulation of crack propagation in homogeneous and functionally graded materials is performed by means of a remeshing algorithm in conjunction with the finite element method. The crack propagation is performed under mixed-mode and non-proportional loading. Each step of crack growth simulation consists of calculation of mixed-mode stress intensity factors by means of a novel formulation of the interaction integral method, determination of crack growth direction based on a specific fracture criterion, and local automatic remeshing along the crack path. The present approach requires a user-defined crack increment at the beginning of the simulation. Crack trajectories obtained by the present numerical simulation are compared with available experimental results.     

On damage accumulations in the cyclic cohesive zone model for XFEM analysis of mixed-mode fatigue crack growth

Predicting mixed-mode fatigue crack propagation is an important and troublesome issue in structure assessment for decades. In the present paper an extended finite element method (XFEM) combined with a new cyclic cohesive zone model (CCZM) is introduced for simulating fatigue crack propagation under mixed-mode loading conditions, which has been implemented in the commercial general purpose software ABAQUS. The algorithm allows introducing a new crack surface at arbitrary locations and directions in a finite element mesh, without re-meshing. The cyclic cohesive zone model is based on the known S–N curves and Goodman diagram for metallic materials and validated by uniaxial tension results. Furthermore, the sensitivity of the model parameter is investigated for mixed-mode fatigue. The virtual crack closure technique has been extended to the cohesive zone model and proposed to calculate the energy release rate for the generalized Paris’ law. Finally, the crack propagation rate and direction under mixed-mode fatigue loading conditions are studied.     

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.     

Adaptive mesh refinement and coarsening for cohesive zone modeling of dynamic fracture

Adaptive mesh refinement and coarsening schemes are proposed for efficient computational simulation of dynamic cohesive fracture. The adaptive mesh refinement consists of a sequence of edge‐split operators, whereas the adaptive mesh coarsening is based on a sequence of vertex‐removal (or edge‐collapse) operators. Nodal perturbation and edge‐swap operators are also employed around the crack tip region to improve crack geometry representation, and cohesive surface elements are adaptively inserted whenever and wherever they are needed by means of an extrinsic cohesive zone model approach. Such adaptive mesh modification events are maintained in conjunction with a topological data structure (TopS). The so‐called PPR potential‐based cohesive model (J. Mech. Phys. Solids 2009; 57 :891–908) is utilized for the constitutive relationship of the cohesive zone model. The examples investigated include mode I fracture, mixed‐mode fracture and crack branching problems. The computational results using mesh adaptivity (refinement and coarsening) are consistent with the results using uniform mesh refinement. The present approach significantly reduces computational cost while exhibiting a multiscale effect that captures both global macro‐crack and local micro‐cracks. Copyright © 2012 John Wiley & Sons, Ltd.     

Discrete cohesive zone model for mixed-mode fracture using finite element analysis

The discrete cohesive zone model (DCZM) is implemented using the finite element (FE) method to simulate fracture initiation and subsequent growth when material non-linear effects are significant. Different from the widely used continuum cohesive zone model (CCZM) where the cohesive zone model is implemented within continuum type elements and the cohesive law is applied at each integral point, DCZM uses rod type elements and applies the cohesive law as the rod internal force vs. nodal separation (or rod elongation). These rod elements have the provision of being represented as spring type elements and this is what is considered in the present paper. A series of 1D interface elements was placed between node pairs along the intended fracture path to simulate fracture initiation and growth. Dummy nodes were introduced within the interface element to extract information regarding the mesh size and the crack path orientation. To illustrate the DCZM, three popular fracture test configurations were examined. For pure mode I, the double cantilever beam configuration, using both uniform and biased meshes were analyzed and the results show that the DCZM is not sensitive to the mesh size. Results also show that DCZM is not sensitive to the loading increment, either. Next, the end notched flexure for pure mode II and, the mixed-mode bending were studied to further investigate the approach. No convergence difficulty was encountered during the crack growth analyses. Therefore, the proposed DCZM approach is a simple but promising tool in analyzing very general two-dimensional crack growth problems. This approach has been implemented in the commercial FEA software ABAQUS^® using a user defined subroutine and should be very useful in performing structural integrity analysis of cracked structures by engineers using ABAQUS^®.     

Modeling of crack propagation via an automatic adaptive mesh refinement based on modified superconvergent patch recovery technique

In this paper, an automated adaptive remeshing procedure is presented for simulation of arbitrary shape crack growth in a 2D finite element mesh. The Zienkiewicz-Zhu error estimator is employed in conjunction with a modified SPR technique based on the recovery of gradients using analytical crack-tip fields in order to obtain more accurate estimation of errors. The optimization of crack-tip singular finite element size is achieved through the adaptive mesh strategy. Finally, several numerical examples are illustrated to demonstrate the effectiveness, robustness and accuracy of computational algorithm in calculation of fracture parameters and prediction of crack path pattern.     

Cohesive-zone-model formulation and implementation using the symmetric Galerkin boundary element method for homogeneous solids

A new symmetric boundary integral formulation for cohesive cracks growing in the interior of homogeneous linear elastic isotropic media with a known crack path is developed and implemented in a numerical code. A crack path can be known due to some symmetry implications or the presence of a weak or bonded surface between two solids. The use of a two-dimensional exponential cohesive law and of a special technique for its inclusion in the symmetric Galerkin boundary element method allows us to develop a simple and efficient formulation and implementation of a cohesive zone model. This formulation is dependent on only one variable in the cohesive zone (relative displacement). The corresponding constitutive cohesive equations present a softening branch which induces to the problem a potential instability. The development and implementation of a suitable solution algorithm capable of following the growth of the cohesive zone and subsequent crack growth becomes an important issue. An arc-length control combined with a Newton–Raphson algorithm for iterative solution of nonlinear equations is developed. The boundary element method is very attractive for modeling cohesive crack problems as all nonlinearities are located along the boundaries (including the crack boundaries) of linear elastic domains. A Galerkin approximation scheme, applied to a suitable symmetric integral formulation, ensures an easy treatment of cracks in homogeneous media and excellent convergence behavior of the numerical solution. Numerical results for the wedge split and mixed-mode flexure tests are presented.     

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.     

Unstructured polygonal meshes with adaptive refinement for the numerical simulation of dynamic cohesive fracture

This paper presents a scheme for adaptive mesh refinement on unstructured polygonal meshes to better capture crack patterns in dynamic cohesive fracture simulations. A randomly seeded polygonal mesh leads to an isotropic discretization of the problem domain, which does not bias crack patterns, but restricts the number of paths that a crack may travel at each node. An adaptive refinement scheme is proposed and investigated through a detailed set of geometric studies. The refinement scheme is selectively chosen to optimize the number of paths that a crack may travel, while still maintaining a conforming domain discretization. The details of the refinement scheme are outlined, along with the criterion used to determine the region of refinement and the method of interpolating nodal attributes. Extrinsic cohesive elements are inserted when and where necessary, and follow the constitutive response of the Park–Paulino–Roesler cohesive model. The influence of bulk and cohesive material heterogeneity is investigated through the use of a statistical distribution of material properties. The adaptive mesh modifications are handled through a compact topological data structure. Numerical examples highlight the features of adaptive refinement in capturing physical fracture patterns while addressing computational cost. Thus, the present approach is a step towards obtaining accurate dynamic fracture patterns and fields with polygonal elements.     

A new M-integral parameter for mixed-mode crack growth in orthotropic viscoelastic material

This paper deals with a new independent path integral which provides the mixed-mode during a creep crack growth process in viscoelastic orthotropic media. The developments are based on an energetic approach using conservative laws. The mixed-mode fracture separation is introduced according to the generalization of the virtual work principle. The fracture algorithm is implemented in a finite element software and coupled with an incremental viscoelastic formulation and an automatic crack growth simulation. This M-integral provides the computation of stress intensity factors and energy release rate for each fracture mode. A numerical validation, in terms of energy release rate and stress intensity factors, is carried out on a CTS specimen under mixed-mode loading for different crack growth speeds.     

Dynamic analysis of fixed cracks in composites by the extended finite element method

This paper is dedicated to simulation of dynamic analysis of fixed cracks in orthotropic media using an extended finite element method. This work is in fact an extension to dynamic problems of the recently developed orthotropic extended finite element method for fracture analysis of composites. In this method, the Heaviside and near-tip enrichment functions are used in the framework of the partition of unity for modeling crack discontinuity and crack-tip singularities within the classical finite element method. In this procedure, elements that include a crack are not required to conform to crack edges. Therefore, mesh generation can be performed without any need to comply to crack edges and the method is capable of modeling the crack propagation without any remeshing. To determine the fracture properties, mixed-mode dynamic stress intensity factors (DSIFs) are evaluated by means of domain separation integral (J^′-integral) method. Results of the proposed method are compared with other available analytical and computational results.     

Probabilistic fracture mechanics by Galerkin meshless methods – part I: rates of stress intensity factors

 This is the first in a series of two papers generated from a study on probabilistic meshless analysis of cracks. In this paper (Part I), a Galerkin-based meshless method is presented for predicting first-order derivatives of stress-intensity factors with respect to the crack size in a linear-elastic structure containing a single crack. The method involves meshless discretization of cracked structure, domain integral representation of the fracture integral parameter, and sensitivity analysis in conjunction with a virtual crack extension technique. Unlike existing finite-element methods, the proposed method does not require any second-order variation of the stiffness matrix to predict first-order sensitivities, and is, consequently, simpler than existing methods. The method developed herein can also be extended to obtain higher-order derivatives if desired. Several numerical examples related to mode-I and mixed-mode problems are presented to illustrate the proposed method. The results show that first-order derivatives of stress-intensity factors using the proposed method agree very well with reference solutions obtained from either analytical (mode I) or finite-difference (mixed mode) methods for the structural and crack geometries considered in this study. For mixed-mode problems, the maximum difference between the results of proposed method and finite-difference method is less than 7. Since the rates of stress-intensity factors are calculated analytically, the subsequent fracture reliability analysis can be performed efficiently and accurately. Received 20 February 2001 / Accepted 19 December 2001     

Automated simulation of fatigue crack propagation for two-dimensional linear elastic fracture mechanics problems by boundary element method

In this paper, automated simulation of multiple crack fatigue propagation for two-dimensional (2D) linear elastic fracture mechanics (LEFM) problems is developed by using boundary element method (BEM). The boundary element method is the displacement discontinuity method with crack-tip elements proposed by the author. Because of an intrinsic feature of the boundary element method, a general growth problem of multiple cracks can be solved in a single-region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Local discretization on the incremental crack extension is performed easily. Further the new adding elements and the existing elements on the existing boundaries are employed to construct easily the total structural mesh representation. Here, the mixed-mode stress intensity factors are calculated by using the formulas based on the displacement fields around crack tip. The maximum circumferential stress theory is used to predict crack stability and direction of propagation at each step. The well-known Paris’ equation is extended to multiple crack case under mixed-mode loadings. Also, the user does not need to provide a desired crack length increment at the beginning of each simulation. The numerical examples are included to illustrate the validation of the numerical approach for fatigue growth simulation of multiple cracks for 2D LEFM problems.     

Fully automatic modelling of mixed-mode crack propagation using scaled boundary finite element method

The newly-developed scaled boundary finite element method (SBFEM) is able to calculate stress intensity factors directly because the singularity in stress solutions at crack tips is analytically represented. By taking this advantage, a mixed-mode crack propagation model based on linear elastic fracture mechanics (LEFM) was developed in this study. A domain is first divided into a few subdomains. Because the dimensions and shapes of subdomains can be flexibly varied and only the domain boundaries or common edges between subdomains are discretised in the SBFEM, a remeshing procedure as simple as in boundary element methods was developed with minimum mesh changes whereas the generality and flexibility of the FEM is well maintained. Fully-automatic modelling of mixed-mode crack propagation is then achieved by combining the remeshing procedure with a propagation criterion. Three mixed-mode examples were modelled. Comparisons of the numerical results with those from available publications show that the developed model is capable of predicting crack trajectories and load-displacement relations accurately and efficiently.     

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.     

用扩展有限元方法模拟混凝土的复合型开裂过程

用扩展有限元法对混凝土梁复合型开裂过程进行了数值模拟。裂纹面间的力学行为采用粘聚裂纹模型来描述,通过引入切向保留刚度考虑剪力分量的影响。开裂方向的计算采用了一种简化的最大切向应力准则。对Arrea和Ingraffea的混凝土梁复合开裂实验进行了数值模拟。计算给出了裂纹萌生、扩展的过程及破坏形态,并获得了与实验结果对比良好的荷载-裂纹开口滑移曲线。结果表明,扩展有限元法通过附加特定的位移模式,使裂纹两侧不连位移场的表达独立于网格划分,是一种能够模拟准脆性材料复合开裂问题的有效方法。     

Theory and numerics for finite deformation fracture modelling using strong discontinuities

A general finite element approach for the modelling of fracture is presented for the geometrically non‐linear case. The kinematical representation is based on a strong discontinuity formulation in line with the concept of partition of unity for finite elements. Thus, the deformation map is defined in terms of one continuous and one discontinuous portion, considered as mutually independent, giving rise to a weak formulation of the equilibrium consisting of two coupled equations. In addition, two different fracture criteria are considered. Firstly, a principle stress criterion in terms of the material Mandel stress in conjunction with a material cohesive zone law, relating the cohesive Mandel traction to a material displacement ‘jump’ associated with the direct discontinuity. Secondly, a criterion of Griffith type is formulated in terms of the material‐crack‐driving force (MCDF) with the crack propagation direction determined by the direction of the force, corresponding to the direction of maximum energy release. Apart from the material modelling, the numerical treatment and aspects of computational implementation of the proposed approach is also thoroughly discussed and the paper is concluded with a few numerical examples illustrating the capabilities of the proposed approach and the connection between the two fracture criteria. Copyright © 2005 John Wiley & Sons, Ltd.