An efficient variable-node XFEM for modeling multiple crack growth: A Matlab object-oriented implementation

Published on Feb 1, 2020in Advances in Engineering Software4.141
路 DOI :10.1016/J.ADVENGSOFT.2019.102750
Tinh Quoc Bui5
Estimated H-index: 5
(Duy Tan University)
Sources
Abstract
Abstract We present a Matlab object-oriented implementation of an efficient variable-node extended finite element method (XFEM) for multiple crack growth simulation in brittle materials. To enhance the accuracy of the present method, not only the local mesh refinement strategy is used, but also the variable-node elements are adopted to directly link different scale elements, saving the computational cost. In order to efficiently conduct coalescence detection, a modification is also presented. A stability analysis is performed to determine active crack tips, and the crack length control scheme is employed to model multiple crack growth. Mesh is locally refined simultaneously with the extension of cracked area via a straightforward manner. A general computer program structure for the Matlab object-oriented programming code for the variable-node XFEM is described. All those desirable features of the present work are demonstrated through five numerical examples with branching in linear-elastic fracture mechanics. In addition, this work provides an efficient computer code written by Matlab programming language for computational fracture simulation, which is for educational and academic uses.
馃摉 Papers frequently viewed together
References68
Newest
#1Manik Bansal (IITR: Indian Institute of Technology Roorkee)H-Index: 4
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 32
Last. St茅phane Bordas (University of Luxembourg)H-Index: 67
view all 4 authors...
We propose a parallel and computationally efficient multi-split XFEM approach for 3-D analysis of heterogeneous materials. In this approach, multiple discontinuities (pores and reinforcement particles) may intersect any given element (we call those elements multi-split elements). These discontinuities are modeled by imposing additional degrees of freedom at the nodes. The main advantage of the proposed scheme is that the mesh size remains independent of the relative distance among the heterogene...
Source
#1Konstantinos Agathos (ETH Zurich)H-Index: 8
#2St茅phane Bordas (China Medical University (Taiwan))H-Index: 67
Last. Eleni Chatzi (ETH Zurich)H-Index: 29
view all 3 authors...
Abstract Partition of unity enrichment is known to significantly enhance the accuracy of the finite element method by allowing the incorporation of known characteristics of the solution in the approximation space. However, in several cases it can further cause conditioning problems for which a number of remedies have been proposed in the framework of the extended/generalized finite element method (XFEM/GFEM). Those solutions often involve significant modifications to the initial method and resul...
Source
#1Subhendu Chakraborty (IISc: Indian Institute of Science)H-Index: 5
#2Sundararajan Natarajan (IISc: Indian Institute of Science)H-Index: 33
Last. St茅phane Bordas (Cardiff University)H-Index: 67
view all 5 authors...
AbstractThis paper presents optimal numerical integration schemes for a family of Polygonal Finite Elements with Schwarz鈥揅hristoffel (SC) conformal mapping with improved convergence properties. With Wachspress interpolation, the arbitrary polygons of a finite element mesh are first mapped to regular polygons using isoparametric mapping and then to the unit-disk using SC mapping. Numerical integration is done over this unit-disk by determining integration points which are optimal with respect to ...
Source
#2Eleni Chatzi (ETH Zurich)H-Index: 29
Last. St茅phane Bordas (Cardiff University)H-Index: 67
view all 3 authors...
A numerical scheme is proposed for the detection of multiple cracks in three dimensional (3D) structures. The scheme is based on a variant of the extended finite element method (XFEM) and a hybrid optimizer solution. The proposed XFEM variant is particularly well-suited for the simulation of 3D fracture problems, and as such serves as an efficient solution to the so-called forward problem. A set of heuristic optimization algorithms are recombined into a multiscale optimization scheme. The introd...
Source
#1Tiantang Yu (Hohai University)H-Index: 33
#2Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 57
Abstract The present paper is concerned with numerical simulation of two-dimensional (2-D) cracks and material interfaces by an effective computational approach. A local mesh refinement in terms of extended finite element method is thus described. The new approach combines a posteriori error estimation algorithm, a local non-conformal mesh connection strategy, and local enrichment. An error estimator based on recovery strain for adaptivity is used; allowing the mesh where it is needed is subsequ...
Source
#1Steffen Beese (Leibniz University of Hanover)H-Index: 2
#2Stefan Loehnert (Leibniz University of Hanover)H-Index: 13
Last. Peter Wriggers (Leibniz University of Hanover)H-Index: 78
view all 3 authors...
In this contribution we present a gradient enhanced damage based method to simulate discrete crack propagation in 3D polycrystalline microstructures. Discrete cracks are represented using the eXtended finite element method. The crack propagation criterion and the crack propagation direction for each point along the crack front line is based on the gradient enhanced damage variable. This approach requires the solution of a coupled problem for the balance of momentum and the additional global equa...
Source
#1Gao Yi (Hohai University)H-Index: 3
#2Tiantang Yu (Hohai University)H-Index: 33
Last. Satoyuki Tanaka (Hiroshima University)H-Index: 23
view all 4 authors...
Abstract In this paper, a new effective approach based on conservative integral approach associated with extended finite element method (XFEM) is developed for evaluating stress intensity factors (SIFs) of bi-material V-notched structures. The XFEM model for bi-material V-notches is established, which owns various features: (a) jump enrichment functions are taken for describing the intersection of notch-faces; (b) eight (real eigenvalue) or sixteen (complex eigenvalue) branch functions are emplo...
Source
#2Giulio Ventura (Polytechnic University of Turin)H-Index: 18
Last. St茅phane Bordas (Duy Tan University)H-Index: 67
view all 4 authors...
Summary phane Bordas thanks funding provided by the European Research Council Starting Independent Research Grant (ERC Stg grant agreement No. 279578)鈥淩ealTCut towards real time multiscale simulation of cutting in non-linear materials with applications to surgical simulation and computer guided surgery鈥. Stephane Bordas also thanks the financial support of the Fonds National de la Recherche Luxembourg FWO-FNR grant INTER/FWO/15/10318764. Eleni Chatzi would like to acknowledge the support of the ...
Source
#2Leandro Lopes da Silva (UFMG: Universidade Federal de Minas Gerais)H-Index: 2
Abstract This work presents an object-oriented implementation of the G/XFEM to model the crack nucleation and propagation in structures made of either linear or nonlinear materials. A discontinuous function along with the asymptotic crack-tip displacement fields are used to represent the crack without explicitly meshing its surfaces. Different approach are explained in detail that are used to capture the crack nucleation within the model and also determine the crack propagation path for various ...
Source
#1Danas Sutula (Cardiff University)H-Index: 5
#2Pierre Kerfriden (Cardiff University)H-Index: 22
Last. St茅phane Bordas (Duy Tan University)H-Index: 67
view all 4 authors...
Abstract The three-part paper deals with energy-minimal multiple crack propagation in a linear elastic solid under quasi-static conditions. The principle of minimum total energy, i.e. the sum of the potential and fracture energies, which stems directly from the Griffith鈥檚 theory of cracks, is applied to the problem of arbitrary crack growth in 2D. The proposed formulation enables minimisation of the total energy of the mechanical system with respect to the crack extension directions and crack ex...
Source
Cited By4
Newest
#1Shuohui Yin (XTU: Xiangtan University)H-Index: 18
#2Ning Zhang (XTU: Xiangtan University)
Last. Yu Cong (Universit茅 Paris-Saclay)
view all 0 authors...
Abstract null null This paper aims to develop a variable-node multiscale extended finite element method (V-XFEM) for dynamic fracture analysis of the linearly uncoupled and coupled physical phenomena in a compact formula. The general governing equations for the linearly uncoupled and coupled physical phenomena are presented in a compact form. The local mesh refinement technique for modeling cracks is used to improve the accuracy and efficiency, in which variable-node elements without modifying t...
Source
#1Raphael Suda (TU Wien: Vienna University of Technology)H-Index: 2
#2Thomas Kiefer (TU Wien: Vienna University of Technology)H-Index: 6
Last. Josef F眉ssl (TU Wien: Vienna University of Technology)H-Index: 18
view all 4 authors...
Abstract Since lateral tensile stresses trigger failure of vertically perforated clay block masonry under vertical loading, reinforcement of the bed joints introduces a new way to improve the vertical resistance of masonry. The aim of this work is to estimate the feasible increase of the vertical compressive strength by means of the eXtended Finite Element Method (XFEM). Using a unit cell approach and a stochastic strength distribution, the increase of the masonry strength鈥檚 5%-quantile could be...
Source
#2Ran Tian (HBUT: Hubei University of Technology)
Last. Fenghua Nie (WUT: Wuhan University of Technology)
view all 6 authors...
One of the most important tasks in the numerical analysis of the fracture problem is to achieve a high-precision approximation of the singular stress fields at crack tips. Partially or fully enrich...
Source
#1Kelson Pothin Wolff (PUC-MG: Pontif铆cia Universidade Cat贸lica de Minas Gerais)H-Index: 1
#2Roque Luiz da Silva Pitangueira (UFMG: Universidade Federal de Minas Gerais)H-Index: 5
Last. Rodrigo Guerra Peixoto (UFMG: Universidade Federal de Minas Gerais)H-Index: 5
view all 3 authors...
Abstract The displacement-based methods for the extraction of the Stress Intensity Factors (SIFs) were created decades ago and, since them, they have been used in a wide range of numerical analysis, most of them based on standard 2D Finite Element Method (FEM). Although the accuracy of the SIFs computed by the displacement-based methods generally be inferior compared with the energy release methods, the use of the former methods have been motivated by its simplicity, ease of implementation and l...
Source
The purpose of this paper is to evaluate some numerical integration strategies used in generalized (G)/extended finite element method (XFEM) to solve linear elastic fracture mechanics problems. A range of parameters are here analyzed, evidencing how the numerical integration error and the computational efficiency are improved when particularities from these examples are properly considered.,Numerical integration strategies were implemented in an existing computational environment that provides a...
Source
This website uses cookies.
We use cookies to improve your online experience. By continuing to use our website we assume you agree to the placement of these cookies.
To learn more, you can find in our Privacy Policy.