Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Computational Science, Engineering & Technology Series
INNOVATION IN COMPUTATIONAL STRUCTURES TECHNOLOGY
Edited by: B.H.V. Topping, G. Montero, R. Montenegro
Enrichment Schemes and Accuracy of the Extended-Generalised Finite Element Method for Modelling Traction-Free and Cohesive Cracks
Q.Z. Xiao and B.L. Karihaloo
School of Engineering, Cardiff University, United Kingdom
Q.Z. Xiao, B.L. Karihaloo, "Enrichment Schemes and Accuracy of the Extended-Generalised Finite Element Method for Modelling Traction-Free and Cohesive Cracks", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Innovation in Computational Structures Technology", Saxe-Coburg Publications, Stirlingshire, UK, Chapter 13, pp 265-286, 2006. doi:10.4203/csets.14.13
Keywords: cohesive crack, crack, enrichment scheme, extended-generalised finite element method (XFEM), stress recovery.
In the modelling of crack problems using the extended-generalised finite element method (XFEM) [1,2,3,4], the crack faces behind the crack tip are modelled by enrichment of discontinuous Heaviside functions and the crack tip region can be modelled by enrichment of proper branch functions such as the true crack tip asymptotic displacement fields [1,2,5,6], or discontinuous Heaviside functions with proper partition of the element including the crack tip . This partition divides the partially cracked element into fully cracked and uncracked parts, and makes the enriched discontinuity conform to the crack.
This study investigates the advantages and disadvantages of various crack tip enrichment schemes, and the accuracy of the resulting crack tip displacement and stress fields. Three methods for evaluating the stresses are compared: direct differentiation of the displacements, simple interpolation of the averaged nodal stress values evaluated from adjacent elements by bi-linear extrapolation from the Gauss points using shape functions (AVG), and the statically admissible stress recovery (SAR) . Both cracks with traction-free faces and cohesive cracks are considered. For cohesive cracks, the accuracy of load-deformation curves, evolution of the fracture process zone (FPZ), crack opening profile, and distribution of the traction in the FPZ are also studied.
If the crack tip region is enriched with branch functions such as the true crack tip asymptotic displacement field, partially cracked elements can be handled directly without further partition. However, sophisticated quadrature rules are required to handle the possible singularity at the tip and/or the nature of the angular oscillations of the enrichment functions . On the other hand, if only the jump function is used in the crack tip region, the partially cracked elements need to be partitioned into a fully cracked part and an uncracked part. However, this enrichment requires only the standard Gauss-Legendre quadrature.
If the true crack tip asymptotic displacement field is used as the enrichment function at the crack tip, but the coefficients appearing in it are assumed to be independent at each enriched node, the accuracy of the method is no different from that obtained by the use of the jump function only for enrichment. However, when the enriched fields at different nodes are enforced to be the same, the enrichment approximation reduces to the real crack tip asymptotic field. The accuracy is then improved and the SIFs can be obtained directly for linear problems [5,6].
Although the above discussion is mainly based on homogeneous crack problems, it is obviously not limited to these problems.
purchase the full-text of this chapter (price £20)