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.     

A consistent partly cracked XFEM element for cohesive crack growth

Present extended finite element method (XFEM) elements for cohesive crack growth may often not be able to model equal stresses on both sides of the discontinuity when acting as a crack‐tip element. The authors have developed a new partly cracked XFEM element for cohesive crack growth with extra enrichments to the cracked elements. The extra enrichments are element side local and were developed by superposition of the standard nodal shape functions for the element and standard nodal shape functions for a sub‐triangle of the cracked element. With the extra enrichments, the crack‐tip element becomes capable of modelling variations in the discontinuous displacement field on both sides of the crack and hence also capable of modelling the case where equal stresses are present on each side of the crack. The enrichment was implemented for the 3‐node constant strain triangle (CST) and a standard algorithm was used to solve the non‐linear equations. The performance of the element is illustrated by modelling fracture mechanical benchmark tests. Investigations were carried out on the performance of the element for different crack lengths within one element. The results are compared with previously obtained XFEM results applying fully cracked XFEM elements, with computational results achieved using standard cohesive interface elements in a commercial code, and with experimental results. The suggested element performed well in the tests. Copyright © 2007 John Wiley & Sons, Ltd.     

A multicorner variable order singularity triangle to model neighbouring singularities

A six noded triangular finite element has been proposed to model variable order singularities in the first derivatives of field variables occurring at multiple locations (up to three) lying in close proximity. The element satisfies the convergence requirements fully. The element will be very useful in the analysis of neighbouring cracks, zig-zag cracks, kinked cracks, etc., involving mechanical and/or thermal loadings. The first degenerate form, which is a two point variable order singularity element, can be used, for example, to cover the distance between two adjacent crack tips. The second degenerate form, which is a one point variable order singularity element, can be used around a crack tip like any other crack tip elements. Six case studies have been reported, three involving mechanical loadings and the other three thermal loadings. These have helped to check the performance of the element and its two degenerate forms. The accuracy of the results in all the examples is good.     

A new crack tip element for the phantom‐node method with arbitrary cohesive cracks

We have developed a new crack tip element for the phantom‐node method. In this method, a crack tip can be placed inside an element. Therefore, cracks can propagate almost independent of the finite element mesh. We developed two different formulations for the three‐node triangular element and four‐node quadrilateral element, respectively. Although this method is well suited for the one‐point quadrature scheme, it can be used with other general quadrature schemes. We provide some numerical examples for some static and dynamic problems. Copyright © 2007 John Wiley & Sons, Ltd.     

An algorithm to generate quadrilateral or triangular element surface meshes in arbitrary domains with applications to crack propagation

A new hybrid algorithm for automatically generating either an all-quadrilateral or an all-triangular element mesh within an arbitrarily shaped domain is described. The input consists of one or more closed loops of straight-line segments that bound the domain. Internal mesh density is inferred from the boundary density using a recursive spatial decomposition (quadtree) procedure. All-triangular element meshes are generated using a boundary contraction procedure. All-quadrilateral element meshes are generated by modifying the boundary contraction procedure to produce a mixed element mesh at half the density of the final mesh and then applying a polygon-splitting procedure. The final meshes exhibit good transitioning properties and are compatible with the given boundary segments which are not altered. The algorithm can support discrete crack growth simulation wherein each step of crack growth results in an arbitrarily shaped region of elements deleted about each crack tip. The algorithm is described and examples of the generated meshes are provided for a representative selection of cracked and uncracked structures.     

Finite element modelling of cracked inelastic shells with large deflections: two‐dimensional and three‐dimensional approaches

Higher utilization of structural materials leads to a need for accurate numerical tools for reliable predictions of structural response. In some instances, both material and geometrical non‐linearities are allowed for, typically in assessments of structural collapse or residual strength in damaged conditions. The present study addresses the performance of surface‐cracked inelastic shells with out‐of‐plane displacements not negligible compared to shell thickness. This situation leads to non‐linear membrane force effects in the shell. Hence, a cracked part of the shell will be subjected to a non‐proportional history of bending moment and membrane force. An important point in the discretization of the problem is whether a two‐dimensional model describes the structural performance sufficiently, or a three‐dimensional model is required. Herein, the two‐dimensional modelling is performed by means of a Mindlin shell finite element. The cracked parts are accounted for by means of inelastic line spring elements. The three‐dimensional models employ eight‐noded solid elements. These models also account for ductile crack growth due to void coalescence by means of a modified Gurson–Tvergaard constitutive model, hence providing detailed solutions that the two‐dimensional simulations can be tested against. Using this, the accuracy of the two‐dimensional approach is checked thoroughly. The analyses show that the two‐dimensional modelling is sufficient as long as the cracks do not grow. Hence, using fracture initiation as a capacity criterion, shell elements and line springs provide acceptable predictions. If significant ductile tearing occurs before final failure, the line spring ligaments have to be updated due to crack growth.     

Modeling arbitrary crack propagation in coupled shell/solid structures with X‐FEM

In this paper, the extended finite element method (X‐FEM) formulation for the modeling of arbitrary crack propagation in coupled shell/solid structures is developed based on the large deformation continuum‐based (CB) shell theory. The main features of the new method are as follows: (1) different kinematic equations are derived for different fibers in CB shell elements, including the fibers enriched by shifted jump function or crack tip functions and the fibers cut into two segments by the crack surface or connecting with solid elements. So the crack tip can locate inside the element, and the crack surface is not necessarily perpendicular to the middle surface. (2) The enhanced CB shell element is developed to realize the seamless transition of crack propagation between shell and solid structures. (3) A revised interaction integral is used to calculate the stress intensity factor (SIF) for shells, which avoids that the auxiliary fields for cracks in Mindlin–Reissner plates cannot satisfy exactly the equilibrium equations. Several numerical examples, including the calculation of SIF for the cracked plate under uniform bending and crack propagation between solid and shell structures are presented to demonstrate the performance of the developed method. Copyright © 2015 John Wiley & Sons, Ltd.     

New crack‐tip elements for XFEM and applications to cohesive cracks

An extended finite element method scheme for a static cohesive crack is developed with a new formulation for elements containing crack tips. This method can treat arbitrary cracks independent of the mesh and crack growth without remeshing. All cracked elements are enriched by the sign function so that no blending of the local partition of unity is required. This method is able to treat the entire crack with only one type of enrichment function, including the elements containing the crack tip. This scheme is applied to linear 3‐node triangular elements and quadratic 6‐node triangular elements. To ensure smooth crack closing of the cohesive crack, the stress projection normal to the crack tip is imposed to be equal to the material strength. The equilibrium equation and the traction condition are solved by the Newton–Raphson method to obtain the nodal displacements and the external load simultaneously. The results obtained by the new extended finite element method are compared to reference solutions and show excellent agreement. Copyright © 2003 John Wiley & Sons, Ltd.     

An over‐deterministic method for calculation of coefficients of crack tip asymptotic field from finite element analysis

An over‐deterministic method has been employed for calculating the stress intensity factors (SIFs) as well as the coefficients of the higher‐order terms in the Williams series expansions in cracked bodies, using the conventional finite element analysis. For a large number of nodes around the crack tip, an over‐determined set of simultaneous linear equations is obtained, and using the fundamental concepts of the least‐squares method, the coefficients of the Williams expansion can be calculated for pure mode I, pure mode II and mixed mode I/II conditions. A convergence study has been conducted to examine the effects of the number of nodes used, the number of terms in Williams expansion and the distance of the selected nodes from the crack tip, on the accuracy of the results. It is shown that the simple method presented in this paper, yields accurate results even for coarse finite element meshes or in the absence of singular elements. The accuracy of SIFs and the coefficients of higher‐order terms are validated by using the available results in the literature.     

Effects of CFRP bond locations on the Mode I stress intensity factor of centre‐cracked tensile steel plates

The fatigue life of cracked steel members can be greatly extended by externally attached carbon fibre reinforced plastics (CFRP), which reduces the stress intensity factors (SIFs) at the crack tip. Access to cracks is sometimes limited and the CFRP has to be attached away from the cracks. There is a lack of knowledge on SIFs for such strengthening scheme. This paper presents the effects of CFRP bond locations on the Mode I SIF of centre‐cracked tensile (CCT) steel plate. The Mode I SIF at the crack tip is calculated using the finite element (FE) models. A correction factor is introduced as a function of CFRP bond location and crack length. The FE results are compared and agree well with experimental tests conducted by the authors. By combining with another two factors (one considering CFRP mechanical properties and the other considering CFRP bond width) derived previously by the authors, SIF formulae are proposed for CFRP reinforced CCT steel plates.     

Mechanisms and modeling of low cycle fatigue crack propagation in a pressure vessel steel Q345

The fatigue process near crack is governed by highly concentrated strain and stress in the crack tip region. Based on the theory of elastic–plastic fracture mechanics, we explore the cyclic J-integral as breakthrough point, an analytical model is presented in this paper to determine the CTOD for cracked component subjected to cyclic axial in-plane loading. A simple fracture mechanism based model for fatigue crack growth assumes a linear correlation between the cyclic crack tip opening displacement (ΔCTOD) and the crack growth rate (da/dN). In order to validate the model and to calibrate the model parameters, the low cycle fatigue crack propagation experiment was carried out for CT specimen made of Q345 steel. The effects of stress ratio and crack closure on fatigue crack growth were investigated by elastic–plastic finite element stress–strain analysis of a cracked component. A good comparison has been found between predictions and experimental results, which shows that the crack opening displacement is able to characterize the crack tip state at large scale yielding constant amplitude fatigue crack growth.     

Oscillatory crack tip triangular elements for finite element analysis of interface cracks

In this paper a special crack tip element has been developed in which displacements and stresses have the same behaviour as those of bi‐material interface cracks with open tips. The element degenerates into a traditional triangular quarter point element in cases of homogeneous cracks. An isoparametric co‐ordinate system (ρ, t) is defined in this study, and numerical techniques using these co‐ordinates to evaluate Jacobian matrices, shape function derivatives, and element stiffness matrices are developed. Also, equations calculating the complex stress intensity factor using displacements are obtained in this study. Numerical results are in good agreement with known analytical solutions in two examples. Copyright © 2003 John Wiley & Sons, Ltd.     

Analysis of three‐dimensional crack initiation and propagation using the extended finite element method

An Erratum has been published for this article in International Journal for Numerical Methods in Engineering 2005, 63(8): 1228. We present a new formulation and a numerical procedure for the quasi‐static analysis of three‐dimensional crack propagation in brittle and quasi‐brittle solids. The extended finite element method (XFEM) is combined with linear tetrahedral elements. A viscosity‐regularized continuum damage constitutive model is used and coupled with the XFEM formulation resulting in a regularized ‘crack‐band’ version of XFEM. The evolving discontinuity surface is discretized through a C⁰ surface formed by the union of the triangles and quadrilaterals that separate each cracked element in two. The element's properties allow a closed form integration and a particularly efficient implementation allowing large‐scale 3D problems to be studied. Several examples of crack propagation are shown, illustrating the good results that can be achieved. Copyright © 2005 John Wiley & Sons, Ltd.     

The extended finite element method for fracture in composite materials

Methods for treating fracture in composite material by the extended finite element method with meshes that are independent of matrix/fiber interfaces and crack morphology are described. All discontinuities and near‐tip enrichments are modeled using the framework of local partition of unity. Level sets are used to describe the geometry of the interfaces and cracks so that no explicit representation of either the cracks or the material interfaces are needed. Both full 12 function enrichments and approximate enrichments for bimaterial crack tips are employed. A technique to correct the approximation in blending elements is used to improve the accuracy. Several numerical results for both two‐dimensional and three‐dimensional examples illustrate the versatility of the technique. The results clearly demonstrate that interface enrichment is sufficient to model the correct mechanics of an interface crack. Copyright © 2008 John Wiley & Sons, Ltd.     

Element‐local level set method for three‐dimensional dynamic crack growth

An approximate level set method for three‐dimensional crack propagation is presented. In this method, the discontinuity surface in each cracked element is defined by element‐local level sets (ELLSs). The local level sets are generated by a fitting procedure that meets the fracture directionality and its continuity with the adjacent element crack surfaces in a least‐square sense. A simple iterative procedure is introduced to improve the consistency of the generated element crack surface with those of the adjacent cracked elements. The discrete discontinuity is treated by the phantom node method which is a simplified version of the extended finite element method (XFEM). The ELLS method and the phantom node technology are combined for the solution of dynamic fracture problems. Numerical examples for three‐dimensional dynamic crack propagation are provided to demonstrate the effectiveness and robustness of the proposed method. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluation of the two‐parameter fracture criterion for various crack configurations made of 7075‐T6 aluminium alloy using finite‐element fracture simulations

For many years, a two‐parameter fracture criterion (TPFC) has been used to correlate and predict failure loads on cracked metallic fracture specimens. The current study was conducted to evaluate the use of the TPFC on a high‐strength aluminium alloy, using elastic‐plastic finite‐element (FE) analyses with the critical crack‐tip‐opening angle (CTOA) fracture criterion. In 1966, Forman generated fracture data on middle‐crack tension, M(T), specimens made of thin‐sheet 7075‐T6 aluminium alloy, which is a quasi‐brittle material. The fracture data included a wide range of specimen half‐widths (w) ranging from 38 to 305 mm. A two‐dimensional FE analysis code (ZIP2D) with a “plane‐strain core” option was used to model the fracture process with a critical CTOA chosen to fit the M(T) test data. Fracture simulations were then conducted on other M(T), single‐edge‐crack tension, SE(T), and bend, SE(B), specimens over a wide range in widths (w = 19‐610 mm). No test data were available on the SE‐type specimens. The results supported the TPFC equation for net‐section stresses less than the material proportional limit. However, some discrepancies in the FE fracture simulations results were observed among the numerical analyses made on the three specimen types. Thus, more research is needed to improve the transferability of the TPFC from the M(T) specimen to both the SE(T) and SE(B) specimens for quasi‐brittle materials.     

Stress intensity factors for three‐dimensional surface cracks using enriched finite elements

The analysis of three‐dimensional crack problems using enriched crack tip elements is examined in this paper. It is demonstrated that the enriched finite element approach is a very effective technique for obtaining stress intensity factors for general three‐dimensional crack problems. The influence of compatibility, integration, element shape function order, and mesh refinement on solution convergence is investigated to ascertain the accuracy of the numerical results. It is shown that integration order has the greatest impact on solution accuracy. Sample results are presented for semi‐circular surface cracks and compared with previously obtained solutions available in the literature. Good agreement is obtained between the different numerical solutions, except in the small zone near the free surface where previously published results have often neglected the change in the stress singularity at the free surface. The enriched crack tip element appears to be particularly effective in this region, since boundary conditions can be easily imposed on the stress intensity factors to accurately represent the correct free surface condition. Copyright © 2002 John Wiley & Sons, Ltd.     

基于响应的梁损伤识别

采用扭转弹簧模拟悬臂梁的损伤，导出了损伤梁位移模态和转角模态的近似公式，获得了损伤梁在单点激励下，零初始条件的位移和转角响应。发现损伤梁的转角响应在损伤处发生阶跃变化，而损伤梁的转角响应对位置的一阶偏导数在损伤处有脉冲两数型突变的特点，从而提出了基于损伤梁转角响应的损伤判据函数。通过对损伤梁和损伤桁架结构的数值模拟，表明提出的判据函数不仅可以利用简谐激励下的响应识别梁的损伤，电可以利用冲击激励下的响应识别桁架的损伤。实验结果表明利用所提出的判据函数，可以同时判别损伤的位置和损伤的程度。     

Combined effect of shear and large deformation on stress intensity factors

An attempt has been made to study the influence of large deformation on the stress intensity factor in a cracked plate subjected to bending including shear deformation. It is assumed that singular terms for stress resultants and strains in the case of large deformation have the same angular distribution and order of singularity as in the case of a linear problem. With this in view the small deformation singular element has been used at the crack tip region surrounded by large deformation plate bending elements. The finite element analysis, based on total Lagrangian formulation combined with the modified Newton–Raphson technique, has been used to get numerical results. Several examples connected with large deformation of cracked plates subjected to bending are studied. Using the above technique stress intensity factors for linear and non-linear cases have been compared.     

Analytical and computational investigation into the influence of the compressive stress cycles on crack growth under variable amplitude loading using CTOD

Fatigue crack growth rate of 6061 T651 aluminum alloy centre‐cracked plates, subject to a variable amplitude loading (VAL), is established using a coupled analytical/computational approach. The method utilises the generalised Willenborg retardation model in conjunction with the Walker model, with several of the models parameters obtained through nonlinear finite element analyses. A two‐parameter approach is also used in this study to explore the influence of compressive stress cycles in a VAL scenario on the crack tip opening displacement, and in turn, on the resulting fatigue crack growth rate of the material. The work includes the evaluation of the crack tip opening displacement and residual stresses in the vicinity of a crack by nonlinear finite element analysis, and application of the Generalised Willenborg model to evaluate the fatigue crack growth; particularly, under the influence of the compressive stress cycles, tensile overloads and underloads have been investigated. The finite element analysis results are compared with the experimental results reported in earlier studies by the authors. The results further demonstrate two important phenomena, that is, (i) the influence of FCG retardation effects due to the tensile overloads and (ii) the acceleration effect due to the applied compressive underloads of a VAL stress–time history.