A Bézier extraction based XIGA approach for three-dimensional crack simulations

Published on Nov 1, 2018in Advances in Engineering Software3.884
· DOI :10.1016/J.ADVENGSOFT.2018.08.014
S.K. Singh6
Estimated H-index: 6
(IITR: Indian Institute of Technology Roorkee),
Indra Vir Singh31
Estimated H-index: 31
(IITR: Indian Institute of Technology Roorkee)
+ 1 AuthorsB.K. Mishra27
Estimated H-index: 27
(IITR: Indian Institute of Technology Roorkee)
Sources
Abstract
Abstract In this work, Bezier extraction based extended isogeometric analysis (XIGA) has been performed for three-dimensional (3-D) crack simulations in linear elastic materials. The Bezier extraction represent the original NURBS basis functions in terms of Bernstein shape functions defined over C 0 continuous isogeometric Bezier elements, and provides an element wise data structure of NURBS which makes the implementation similar to extended finite element method (XFEM). A crack in the domain is modeled by Heaviside and crack tip enrichment functions. The values of mixed-mode stress intensity factors (SIFs) are obtained by interaction integral approach using virtual cuboidal domain. Various 3-D crack problems are solved by Bezier extraction based XIGA, and the results are compared with XFEM and analytical solutions.
📖 Papers frequently viewed together
3,633 Citations
46 Citations
200 Citations
References95
Newest
#1Tiantang Yu (Hohai University)H-Index: 32
#2Wenjiang Lai (Hohai University)H-Index: 2
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 3 authors...
The isogeometric analysis (IGA), which employs non-uniform rational B-splines (NURBS) basis as shape functions for both representation of geometry and approximation of the field variables, owns several inherent advantages to become an effective numerical method such as an exact geometry description with fewer control points, high-order continuity, and high accuracy. Unlike the C0-continuity shape functions in the conventional finite element method (FEM), the high-order basis functions in IGA are...
8 CitationsSource
#1V.B. Pandey (IITR: Indian Institute of Technology Roorkee)H-Index: 4
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 31
Last. Vikas Kumar (DMRL: Defence Metallurgical Research Laboratory)H-Index: 14
view all 6 authors...
In the present work, elasto-plastic creep crack growth simulations are performed using continuum damage mechanics and extended finite element method. Liu–Murakami creep damage model and explicit time integration scheme are used to evaluate the creep strain and damage variable for various materials at different temperatures. Compact tension and C-shaped tension specimens are selected for the simulation of crack growth analysis. For damage evaluation, both local and nonlocal approaches are employe...
35 CitationsSource
#1Manish Kumar (IITR: Indian Institute of Technology Roorkee)H-Index: 5
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 31
Last. Vikas Kumar (DMRL: Defence Metallurgical Research Laboratory)H-Index: 14
view all 6 authors...
Abstract In this work, an extended finite element based fully implicit algorithm is proposed to predict the creep deformation and creep crack growth in elasto-plastic-creeping solids. The effect of stress relaxation and redistribution due to creep strain is captured through elasto-plastic-creep analysis. In creep crack growth analysis, mostly C ∗ is used to compute the creep crack growth which includes only extensive creep whereas C ( t ) -integral includes small-scale creep, transition creep an...
28 CitationsSource
#1S.K. Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 6
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 31
Last. G. Bhardwaj (Thapar University)H-Index: 12
view all 4 authors...
Abstract In this work, an extended isogeometric analysis (XIGA) is used for the analysis of through-thickness crack in a homogeneous and isotropic plate. In isogeometric analysis (IGA), non-uniform rational B-splines (NURBS) are used as a basis function. The plate kinematics is modelled by Reddy’s higher-order shear deformation theory (HSDT). The C 1 continuity requirement of HSDT can be easily fulfilled by the NURBS basis functions. In order to obtain the plate fracture parameters (moment inten...
29 CitationsSource
#2Giulio Ventura (Polytechnic University of Turin)H-Index: 18
Last. Stéphane Bordas (Duy Tan University)H-Index: 65
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)“RealTCut 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 ...
37 CitationsSource
#1Himanshu Pathak (IIT Mandi: Indian Institute of Technology Mandi)H-Index: 63
Abstract In the present technological scenario, functionally graded materials (FGMs) are being introduced by scientific community for a wide range of effective engineering applications especially in aerospace, nuclear reactor, bio-medical and military uses. Safety and long service life of engineering components made from FGMs are one of the major objectives of engineering design process. The prediction of fracture behaviour is quite essential for safe designing of such critical components. In th...
20 CitationsSource
#1Manik Bansal (IITR: Indian Institute of Technology Roorkee)H-Index: 4
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 31
Last. I.A. Khan (BARC: Bhabha Atomic Research Centre)H-Index: 5
view all 5 authors...
Abstract In the present study, a two-scale stochastic framework has been proposed for predicting the failure strength probability of heterogeneous materials. The analysis at both scales (meso and macro) is performed under plane stress condition. The meso-scale analysis is performed by XFEM whereas the macro-scale analysis is performed by FEM. The heterogeneities (pores and reinforced particles) are considered at meso-scale. The effect of shape, size, clustering and volume fraction of heterogenei...
13 CitationsSource
#1Danas Sutula (Cardiff University)H-Index: 5
#2Pierre Kerfriden (Cardiff University)H-Index: 21
Last. Stéphane Bordas (Duy Tan University)H-Index: 65
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’s 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...
38 CitationsSource
#1Pedro M. A. Areias (University of Évora)H-Index: 37
#2Timon Rabczuk (IST: Instituto Superior Técnico)H-Index: 101
Abstract Realistic 3D finite strain analysis and crack propagation with tetrahedral meshes require mesh refinement/division. In this work, we use edges to drive the division process. Mesh refinement and mesh cutting are edge-based. This approach circumvents the variable mapping procedure adopted with classical mesh adaptation algorithms. The present algorithm makes use of specific problem data (either level sets, damage variables or edge deformation) to perform the division. It is shown that glo...
100 CitationsSource
#1Danas Sutula (Cardiff University)H-Index: 5
#2Pierre Kerfriden (Cardiff University)H-Index: 21
Last. Stéphane Bordas (Duy Tan University)H-Index: 65
view all 4 authors...
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’s 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 extension l...
50 CitationsSource
Cited By23
Newest
#1Keke Li (Hohai University)H-Index: 1
#2Tiantang Yu (Hohai University)H-Index: 32
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 3 authors...
Abstract null null Limit analysis can be used to directly calculate the ultimate load of structures without considering cumbersome elastic–plastic analysis. We present an adaptive mesh refinement strategy by mean of extended isogeometric analysis (XIGA) in association with second-order cone programming (SOCP) for kinematic limit analysis of hole and inclusion problems. The location of hole and inclusion can be captured by level set function without regard for its interfaces. The local refined no...
Source
#1Huilong Ren (Bauhaus University, Weimar)H-Index: 11
#2Xiaoying Zhuang (Tongji University)H-Index: 55
Last. Timon Rabczuk (Ton Duc Thang University)H-Index: 101
view all 5 authors...
The derivation of nonlocal strong forms for many physical problems remains cumbersome in traditional methods. In this paper, we apply the variational principle/weighted residual method based on nonlocal operator method for the derivation of nonlocal forms for elasticity, thin plate, gradient elasticity, electro-magneto-elasticity and phase-field fracture method. The nonlocal governing equations are expressed as an integral form on support and dual-support. The first example shows that the nonloc...
1 CitationsSource
#1Wenbin Hou (DUT: Dalian University of Technology)H-Index: 10
#2Kai Jiang (DUT: Dalian University of Technology)H-Index: 1
Last. Ping Hu (DUT: Dalian University of Technology)H-Index: 21
view all 5 authors...
Abstract Strong imposing Dirichlet boundary conditions remain a challenge for eXtended Finite Element Methods (XFEMs) or eXtended Isogeometric Analysis (XIGA) of discontinuous problems. Moreover, the physical meaning of the additional Degrees of Freedom (DOFs) in the displacement expression of XFEMs or XIGA is not clarified. To address these issues, we proposed a new method that combines XIGA and B++ splines for modeling the fracture behaviors of single and multiple cracks in 2D elasticity solid...
Source
#1Xiaofei Hu (DUT: Dalian University of Technology)H-Index: 15
#2Jiaxi Li (DUT: Dalian University of Technology)
Last. Weian Yao (DUT: Dalian University of Technology)H-Index: 12
view all 4 authors...
Abstract null null This study further extends a developed symplectic analytical singular element (SASE) library to include three-dimensional (3D) cracks in dissimilar materials. We construct a trial symplectic eigen solution for 3D bimaterial interfacial cracks and use it as the crack tip enrichment to define the interior displacement field of a new singular element. The symplectic approach is applied for the numerical modelling of 3D bimaterial cracks for the first time. In the derivation, a no...
Source
#1Hongjun Yu (HIT: Harbin Institute of Technology)H-Index: 15
#2Meinhard Kuna (Freiberg University of Mining and Technology)H-Index: 31
Abstract A critical overview of past research on the interaction integral (I-integral) method is presented. The I-integral is the two-state mutual energy release rate decided by a designable auxiliary field and the actual field. Due to the designability of the auxiliary field, the I-integral has been used extensively in extracting the individual stress intensity factors (SIFs) and the T-stress of a crack in single materials or at a bi-material interface. The practical implementation of the I-int...
1 CitationsSource
#1S.K. Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 6
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 31
Abstract In this work, we present an extended isogeometric analysis (XIGA) for cracked functionally graded magneto-electro-elastic (FGMEE) material. The material properties inside the FGMEE domain vary exponentially. An electrically and magnetically impermeable crack in the FGMEE domain is modeled by adding appropriate enrichment functions into the isogeometric analysis (IGA) approximation. Heaviside function is utilized to model the displacement, electric, and magnetic potential along the crack...
3 CitationsSource
#1Huilong Ren (Bauhaus University, Weimar)H-Index: 11
#2Xiaoying ZhuangH-Index: 55
Last. Timon RabczukH-Index: 101
view all 5 authors...
The derivation of nonlocal strong forms for many physical problems remains cumbersome in traditional methods. In this paper, we apply the variational principle/weighted residual method based on nonlocal operator method for the derivation of nonlocal forms for elasticity, thin plate, gradient elasticity, electro-magneto-elasticity and phase field fracture method. The nonlocal governing equations are expressed as integral form on support and dual-support. The first example shows that the nonlocal ...
#1V.B. Pandey (IITR: Indian Institute of Technology Roorkee)H-Index: 4
#2Indra Vir Singh (IITR: Indian Institute of Technology Roorkee)H-Index: 31
Last. B.K. Mishra (IITR: Indian Institute of Technology Roorkee)H-Index: 27
view all 3 authors...
Abstract In this paper, a strain-based damage model is developed to account the effect of different strain ratios on the fatigue failure. Thereafter, continuum damage mechanics and extended finite element method based framework is proposed to simulate the low cycle fatigue crack growth simulations. The proposed damage model is imposed in the developed framework to evaluate the fatigue life of cracked specimen. Chaboche mixed hardening model is implemented to determine the stress–strain accuratel...
Source
#1S. S. Durga Rao (Indian Institute of Technology Madras)H-Index: 1
#2Sethuraman Raju (Indian Institute of Technology Madras)H-Index: 2
Abstract Ill conditioning of system matrices is one of the drawbacks of Generalized/eXtended IsoGeometric Analysis (GIGA/XIGA), which results in increase in computational cost in the case of iterative solvers or erroneous results in the case of direct solvers. Linear dependency between the basis functions which can be measured in terms of smallest angle among all the bases is one of the attributes of ill conditioning in XIGA. In the present study, a simplified method to find the smallest angle a...
Source
#1Aanchal Yadav (Thapar University)H-Index: 2
#2R.U. Patil (IITs: Indian Institutes of Technology)H-Index: 7
Last. G. Bhardwaj (Thapar University)H-Index: 12
view all 5 authors...
In this work, the extended isogeometric analysis (XIGA) is presented for the analysis of cracks in elastic material subjected to thermal, mechanical, and thermo-mechanical loading. The comparative ...
4 CitationsSource