Static and dynamic fracture analysis in elastic solids using a multiscale extended isogeometric analysis

Published on Feb 15, 2019in Engineering Fracture Mechanics3.426
· DOI :10.1016/J.ENGFRACMECH.2018.12.024
Shuohui Yin18
Estimated H-index: 18
(XTU: Xiangtan University),
Tiantang Yu32
Estimated H-index: 32
(Hohai University)
+ 2 AuthorsShuitao Gu3
Estimated H-index: 3
(Chongqing University)
Source
Abstract
Abstract The main objective of this paper is to develop a novel multiscale approach for simulating static and stationary dynamic crack problems in two-dimensional elastic solids through an extended isogeometric analysis using the Non-Uniform Rational B-Splines (NURBS). This multiscale strategy consists of two techniques that the Nitsche’s method is applied for coupling different length-scales meshes on one hand, and on the other hand, a local enrichment technique is taken to describe the cracks independently of the mesh in the refined region by introducing the enrichment function into the isogeometric analysis displacement approximation. This method inherits the advantage of XIGA and simultaneously improves his limitation on the local mesh refinement. The proposed approach possesses several advantages in fracture modeling: (i) the singularities of the crack-tip fields can be captured precisely and conveniently by the enrichment techniques; (ii) higher-order convergence rates is achieved and (iii) stress smoothing results are achieved owing to higher continuity of NURBS basis functions. Numerical results of static and stationary dynamic crack problems validate its accuracy and efficiency.
📖 Papers frequently viewed together
157 Citations
37 Citations
29 Citations
References74
Newest
#1Yuri Bazilevs (UCSD: University of California, San Diego)H-Index: 77
#2Marco S. Pigazzini (UCSD: University of California, San Diego)H-Index: 4
Last. Hyonny Kim (UCSD: University of California, San Diego)H-Index: 13
view all 4 authors...
In this two-part paper we introduce a new formulation for modeling progressive damage in laminated composite structures. We adopt a multi-layer modeling approach, based on Isogeometric Analysis (IGA), where each ply or lamina is represented by a spline surface, and modeled as a Kirchhoff–Love thin shell. Continuum Damage Mechanics is used to model intralaminar damage, and a new zero-thickness cohesive-interface formulation is introduced to model delamination as well as permitting laminate-level ...
17 CitationsSource
#1Jiming Gu (Hohai University)H-Index: 7
#2Tiantang Yu (Hohai University)H-Index: 32
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 6 authors...
Abstract In this paper, we present an effective computational approach that combines an adaptive extended isogeometric analysis (XIGA) method with locally refined (LR) B-splines and level set methods for modeling multiple inclusions in two-dimensional (2D) elasticity problems. The advantage of XIGA is to model inclusions without considering internal inclusion interfaces by additional functions. Multiple level set functions are used to represent the location of inclusion interfaces and to define ...
37 CitationsSource
#1Yujie Guo (NUAA: Nanjing University of Aeronautics and Astronautics)H-Index: 9
#2Jason HellerH-Index: 1
Last. Dominik Schillinger (UMN: University of Minnesota)H-Index: 25
view all 5 authors...
Abstract Following a series of recent innovations, isogeometric shell analysis based on trimmed CAD surfaces is currently being developed into an accurate, efficient and mature design-through-analysis methodology. This work contributes to this emerging technology with respect to the following aspects. On the analysis side, we present a robust variationally consistent Nitsche-type formulation for thin shells at large deformations that weakly enforces coupling constraints at trimming curves. On th...
52 CitationsSource
#1Roland Kruse (Braunschweig University of Technology)H-Index: 9
#2Nhon Nguyen-Thanh (Braunschweig University of Technology)H-Index: 21
Last. L. De Lorenzis (Braunschweig University of Technology)H-Index: 22
view all 4 authors...
This paper presents an isogeometric formulation for frictionless contact between deformable bodies, based on the recently proposed concept of the third medium. This concept relies on continuum formulations not only for the contacting bodies but also for a fictitious intermediate medium in which the bodies can move and interact. Key to the formulation is a suitable definition of the constitutive behavior of the third medium. In this work, based on a number of numerical tests, the role of the mate...
8 CitationsSource
#1Cosmin Anitescu (Bauhaus University, Weimar)H-Index: 15
#2Naim HossainH-Index: 1
Last. Timon Rabczuk (Ton Duc Thang University)H-Index: 101
view all 3 authors...
Abstract An adaptive higher-order method based on a generalization of polynomial/rational splines over hierarchical T-meshes (PHT/RHT-splines) is introduced. While most problems considered in isogeometric analysis can be solved efficiently when the solution is smooth, many non-trivial simulations have rough solutions. This can be caused, for example, by the presence of re-entrant corners in the domain. For such problems, a tensor-product basis is less suitable for resolving the singularities tha...
50 CitationsSource
#1Nhon Nguyen-Thanh (NTU: Nanyang Technological University)H-Index: 21
#2Kun Zhou (NTU: Nanyang Technological University)H-Index: 63
Summary Adaptive local refinement is one of the main issues for isogeometric analysis (IGA). In this paper, an adaptive extended IGA (XIGA) approach based on polynomial splines over hierarchical T-meshes (PHT-splines) for modeling crack propagation is presented. The PHT-splines overcome certain limitations of non-uniform rational B-splines (NURBS)-based formulations; in particular, they make local refinements feasible. To drive the adaptive mesh refinement, we present a recovery-based error esti...
43 CitationsSource
#1Shuo Liu (Hohai University)H-Index: 5
#2Tiantang Yu (Hohai University)H-Index: 32
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 3 authors...
Abstract This paper presents an effective plate formulation coupling the merits of isogeometric analysis (IGA) and a new non-classical simple first-order shear deformation theory (SFSDT) for static bending, free vibration, and buckling of functionally graded (FG) moderately thick microplates. In contrast to the conventional first-order shear deformation theory (FSDT), the new SFSDT adopted here inherently owns several advantages such as free from shear-locking, capturing the shear-deformation ef...
47 CitationsSource
#1Qui X. Lieu (Sejong University)H-Index: 11
#2Jaehong Lee (Sejong University)H-Index: 43
Abstract A novel multi-resolution scheme to perform the multi-material topology optimization in the framework of isogeometric analysis (IGA) is proposed. To create high resolution optimized designs with a lower computational cost, a new variable parameter space is supplemented to define design density variables and represent the optimal material distribution utilizing the bivariate B-spline basis functions. These functions are simply obtained by k -refinement strategy in the IGA. The non-uniform...
58 CitationsSource
#1Xiaoqiong DuH-Index: 1
#2Dixiong Yang (DUT: Dalian University of Technology)H-Index: 20
Last. Shi Li (DUT: Dalian University of Technology)H-Index: 2
view all 6 authors...
This paper presents a new family of explicit time integration algorithms with controllable numerical dissipation for structural dynamic problems by utilizing the discrete control theory. Firstly, the equilibrium equation of the implicit Yu-α algorithm is adopted, and the recursive formulas of velocity and displacement for the explicit CR algorithm are used in the algorithms. Then, the transfer function and characteristic equation of the algorithms with integration coefficients are obtained by th...
13 CitationsSource
#1Zhaohui Xia (RPI: Rensselaer Polytechnic Institute)H-Index: 6
#2Yingjun Wang (SCUT: South China University of Technology)H-Index: 14
Last. Chao Mei (HUST: Huazhong University of Science and Technology)H-Index: 2
view all 4 authors...
This paper proposes a new level set-based topology optimization (TO) method using a parallel strategy of Graphics Processing Units (GPUs) and the isogeometric analysis (IGA). The strategy consists of parallel implementations for initial design domain, IGA, sensitivity analysis and design variable update, and the key issues in the parallel implementation, e.g., the parallel assembly race condition, are discussed in detail. The computational complexity and parallelization of the different steps in...
28 CitationsSource
Cited By18
Newest
#1Xu Wang (CityU: City University of Hong Kong)
#2Weiyin Ma (CityU: City University of Hong Kong)H-Index: 22
Abstract null null Analysis-suitable parameterization is a fundamental problem in IGA (IsoGeometric Analysis) implementation which significantly influences computational accuracy and efficiency. This paper proposes a variational framework to address the problem of producing a smooth parameterization of a computational domain represented in B-spline form. In order to control both angle and area distortions, a weighted and modified Liao functional is constructed. The weighting function is a modifi...
Source
#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
#1Hongting Yuan (Hohai University)H-Index: 2
#2Tiantang Yu (Hohai University)H-Index: 32
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 3 authors...
Abstract The objective of this article is to simulate crack growth of complex Mindlin–Reissner plates by developing an adaptive multi-patch extended isogeometric analysis (XIGA). Nitsche’s method is used to treat continuity between multi-patches or the coupling of non-conforming meshes, exactly describing the geometry of complex plates. The computational meshes in XIGA are independent of the cracks by introducing some enrichment functions into the displacement approximation based on the partitio...
2 CitationsSource
#1Fei Yan (CAS: Chinese Academy of Sciences)H-Index: 14
#2Wei Zhang (CAS: Chinese Academy of Sciences)H-Index: 1
Last. Shaojun Li (CAS: Chinese Academy of Sciences)H-Index: 30
view all 4 authors...
Abstract A continuous-discontinuous cellular automaton method (CDCA) is developed for dynamic crack propagation. First, a separated time integral scheme is built, in which time difference is divided into two separated time steps: one balances the change of the geometric discretization caused by crack growing, and the other balances the dynamic loading. Second, an adaptive energy conservation strategy is proposed for the CDCA that includes the change in the degree of nodal freedom for the crack t...
2 CitationsSource
#1Mansoor Shamloofard (Sharif University of Technology)H-Index: 1
#2Ahmad Assempour (Sharif University of Technology)H-Index: 16
Recently, eliminating the gap between design and formability analysis of sheet metal parts has been studied to simulate sheet metal stamping processes. In this regard, a transfer-based inverse isogeometric formulation has been proposed. This method has various advantages such as solving the governing equations in two-dimensional networks without any concern about the convergence; however, it neglects the bending effect which is a major contributor in die/punch profile radii. The present work aim...
1 CitationsSource
#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
Abstract This contribution establishes a coupled numerical method to show its effectiveness in the simulation of LEFM problems. The semi-meshless method uses the MLPG1 meshfree theory to discrete governing equations while geometry is modeled by simple CPS3 triangular cells. The cells help to propose an efficient neighboring method for meshfree interpolations. Three benchmark homogeneous cases are considered whose analytical solutions are available to carefully examine the accuracy of this method...
2 CitationsSource
#1Weihua Fang (Ministry of Water Resources)H-Index: 3
#2Xin Chen (Hohai University)H-Index: 1
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 4 authors...
Abstract This paper is devoted to numerical investigation of holes/voids effects on crack growth in solids using a locally refined (LR) B-splines extended isogeometric analysis (XIGA). We particularly focus our attention on crack-hole interaction analysis. Special enrichment functions captured discontinuity and singularity are used, by which the computational mesh are thus independent of both cracks and holes, eliminating re-meshing when modeling crack growth. Geometries of arbitrary shaped crac...
3 CitationsSource
#1Karsten Paul (RWTH Aachen University)H-Index: 3
#2Christopher Zimmermann (RWTH Aachen University)H-Index: 3
Last. Roger A. Sauer (RWTH Aachen University)H-Index: 25
view all 4 authors...
This work presents numerical techniques to enforce continuity constraints on multi-patch surfaces for three distinct problem classes. The first involves structural analysis of thin shells that are described by general Kirchhoff-Love kinematics. Their governing equation is a vector-valued, fourth-order, nonlinear, partial differential equation (PDE) that requires at least C^1continuity within a displacement-based finite element formulation. The second class are surface phase separations modele...
5 CitationsSource
#1Jiming Gu (Hohai University)H-Index: 7
#2Tiantang Yu (Hohai University)H-Index: 32
Last. Tinh Quoc Bui (TITech: Tokyo Institute of Technology)H-Index: 54
view all 6 authors...
Abstract Reliable prediction and thorough understanding of crack growth in engineering materials and structures are of challenging problems, but scientific and technical community has been continuously pursuing an efficient numerical approach to crack propagation simulation. We present in this paper a significant extension and development of the adaptive extended isogeometric analysis (XIGA) based on LR B-splines (locally refined B-splines), which facilitates the local refinement, for accurately...
11 CitationsSource