Adaptive extended isogeometric analysis based on PHT-splines for thin cracked plates and shells with Kirchhoff–Love theory

Published on Jul 5, 2019in Applied Mathematical Modelling3.633
· DOI :10.1016/J.APM.2019.07.002
H.S. Yang2
Estimated H-index: 2
(BIT: Beijing Institute of Technology),
H.S. Yang3
Estimated H-index: 3
(BIT: Beijing Institute of Technology),
Chunying Dong19
Estimated H-index: 19
(BIT: Beijing Institute of Technology)
Source
Abstract
Abstract In this paper, a posteriori error estimation and mesh adaptation approach for thin plate and shell structures of through-the-thickness crack is presented. This method uses the extended isogeometric analysis (XIGA) based on PHT-splines (Polynomial splines over Hierarchical T-meshes), which is abbreviated as XIGA-PHT. In XIGA-PHT, the isogeometric displacement approximation is locally enriched with enrichment functions, which efficiently capture the displacement discontinuity across the crack face as well as the stress singularity in the vicinity of the crack tip. On the one hand, the rotational degrees of freedom (RDOFs) are not required in Kirchhoff-Love theory, which drastically reduces the complexity of enrichment mode and computational scale for crack analysis. On the other hand, the PHT-splines basis functions can automatically satisfy the requirement of C 1 -continuity for the Kirchhoff-Love theory. Moreover, the PHT-splines facilitate the local refinement, which is the deficiency of NURBS-based isogeometric formulations. The local refinement is highly suitable for adaptive analysis. The stress recovery-based posteriori error estimator combined with the superconvergent patch recovery (SPR) technique is used to evaluate the approximate local discretization error. A new strategy for selecting enriched recovered functions in the enriched areas was proposed. Special functions extracted from the asymptotic stress solutions are applied to obtain the recovered stress field in the enriched area. The results of stress intensity factors or J-integral values obtained by the adaptive XIGA-PHT are compared with reference solutions. Several thin plate and shell illustrative examples demonstrate the effectiveness and accuracy of the proposed adaptive XIGA-PHT.
📖 Papers frequently viewed together
15 Citations
301 Citations
1,008 Citations
References81
Newest
#1Javier Videla (University of Chile)H-Index: 5
#2Cosmin Anitescu (Bauhaus University, Weimar)H-Index: 15
Last. Elena Atroshchenko (University of Chile)H-Index: 11
view all 5 authors...
Abstract In this work, we demonstrate the application of PHT-splines for time-harmonic acoustic problems, modeled by the Helmholtz equation. Solutions of the Helmholtz equation have two features: global oscillations associated with the wave number and local gradients caused by geometrical irregularities. We show that after a sufficient number of degrees of freedom is used to approximate global oscillations, adaptive refinement can capture local features of the solution. We compare residual-based...
24 CitationsSource
#1Konstantinos Agathos (ETH Zurich)H-Index: 8
#2Stéphane Bordas (China Medical University (Taiwan))H-Index: 65
Last. Eleni Chatzi (ETH Zurich)H-Index: 26
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...
29 CitationsSource
#1Peng Yu (Cardiff University)H-Index: 19
#2Cosmin Anitescu (Bauhaus University, Weimar)H-Index: 15
Last. Pierre Kerfriden (Cardiff University)H-Index: 21
view all 5 authors...
Abstract This paper presents a novel methodology of local adaptivity for the frequency-domain analysis of the vibrations of Reissner–Mindlin plates. The adaptive discretization is based on the recently developed Geometry Independent Field approximaTion (GIFT) framework, which may be seen as a generalization of the Iso-Geometric Analysis (IGA). Within the GIFT framework, we describe the geometry of the structure exactly with NURBS (Non-Uniform Rational B-Splines), whilst independently employing P...
22 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 5 authors...
Abstract We present an effective adaptive analysis procedure in terms of isogeometric analysis (IGA) based on locally refined (LR) B-splines for two-dimensional elasticity problems. One major advantage of the adopted LR B-splines over B-splines and non-uniform rational B-splines (NURBS) is the local refinement, which makes it highly suitable for adaptive analysis. The Zienkiewicz–Zhu estimation is used on error estimator based local refinement. The local refinement is implemented by the structur...
41 CitationsSource
#1Annalisa Buffa (EPFL: École Polytechnique Fédérale de Lausanne)H-Index: 49
#2Eduardo M. GarauH-Index: 7
In this paper, we develop a function-based a posteriori error estimators for the solution of linear second-order elliptic problems considering hierarchical spline spaces for the Galerkin discretiza...
10 CitationsSource
#1Elena Atroshchenko (University of Chile)H-Index: 11
#2Satyendra Tomar (University of Luxembourg)H-Index: 15
Last. Stéphane Bordas (Cardiff University)H-Index: 65
view all 4 authors...
This paper presents an approach to generalize the concept of isogeometric analysis by allowing different spaces for the parameterization of the computational domain and for the approximation of the solution field. The method inherits the main advantage of isogeometric analysis, ie, preserves the original exact computer‐aided design geometry (for example, given by nonuniform rational B‐splines), but allows pairing it with an approximation space, which is more suitable/flexible for analysis, for e...
63 CitationsSource
#1Nhon Nguyen-Thanh (NTU: Nanyang Technological University)H-Index: 21
#2Weidong Li (NTU: Nanyang Technological University)H-Index: 6
Last. Kun Zhou (NTU: Nanyang Technological University)H-Index: 63
view all 3 authors...
This paper develops a coupling approach which integrates the meshfree method and isogeometric analysis (IGA) for static and free-vibration analyses of cracks in thin-shell structures. In this approach, the domain surrounding the cracks is represented by the meshfree method while the rest domain is meshed by IGA. The present approach is capable of preserving geometry exactness and high continuity of IGA. The local refinement is achieved by adding the nodes along the background cells in the meshfr...
22 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...
47 CitationsSource
#1Lin Chen (University of Sheffield)H-Index: 11
#2Erik Jan LingenH-Index: 1
Last. René de Borst (University of Sheffield)H-Index: 44
view all 3 authors...
Summary Adaptive hierarchical refinement in isogeometric analysis is developed to model cohesive crack propagation along a prescribed interface. In the analysis, the crack is introduced by knot insertion in the NURBS basis, which yields C−1 continuous basis functions. To capture the stress state smoothly ahead of the crack tip, the hierarchical refinement of the spline basis functions is used starting from a coarse initial mesh. A multi-level mesh is constructed, with a fine mesh used for quanti...
22 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...
38 CitationsSource
Cited By12
Newest
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...
Source
#1Chuang Xu (BIT: Beijing Institute of Technology)H-Index: 1
#2Chunying Dong (BIT: Beijing Institute of Technology)H-Index: 19
Abstract The isogeometric analysis boundary element method (IGABEM) has great potential in the simulation of heat conduction problems due to its exact geometric representation and good approximation properties. In this paper, the radial integration IGABEM (RI-IGABEM) is proposed to solve isotropic heat conduction problems in inhomogeneous media with heat source. Similar to traditional BEM, the domain integrals cannot be avoided since the foundational solution for the Laplace equation is used to ...
1 CitationsSource
#1Amir Reza Isazadeh (Sharif University of Technology)H-Index: 1
#2Mansoor Shamloofard (Sharif University of Technology)H-Index: 1
Last. Ahmad Assempour (Sharif University of Technology)H-Index: 16
view all 3 authors...
Abstract Recently, isogeometric methodology has been successfully implemented in one-step inverse analysis of sheet metal stamping processes. However, these models are not capable of analyzing forming processes which require severe deformation and/or several forming stages. This paper presents a multi-step inverse isogeometric methodology to enhance the precision of one-step models in predictions of the initial blank, strain distributions, and drawability of the formed parts. This methodology de...
1 CitationsSource
#1Hao Shi (THU: Tsinghua University)
#2Haiyu Qiao (HUST: Huazhong University of Science and Technology)
Last. Houfa Shen (THU: Tsinghua University)H-Index: 12
view all 4 authors...
Source
#1Zhenyu Liu (ZJU: Zhejiang University)H-Index: 13
#2Minglong Yang (ZJU: Zhejiang University)H-Index: 3
Last. Jianrong Tan (ZJU: Zhejiang University)H-Index: 22
view all 4 authors...
ABSTRACT A new stochastic isogeometric analysis method based on reduced basis vectors (SRBIGA) is proposed for engineering structures with random field material properties and external loads. Based on the Galerkin isogeometric functions, the proposed SRBIGA applies the Karhunen-Loeve expansion to discretize the random field uncertainties. Inspired by the stochastic Krylov subspace theory, the structural responses of linear elasticity structures with random field uncertainties are represented bas...
2 CitationsSource
#1H.S. Yang (BIT: Beijing Institute of Technology)H-Index: 3
#2Chunying Dong (BIT: Beijing Institute of Technology)H-Index: 19
Last. Yonghong Wu (BIT: Beijing Institute of Technology)H-Index: 63
view all 3 authors...
Abstract In this paper, a three-dimensional (3D) isogeometric finite element (FE) and boundary element (BE) coupling approach based on non-uniform rational B-splines (NURBS) is developed to simulate the deformed body. Based on the geometric exactness and higher smoothness of isogeometric analysis (IGA), NURBS discretization is applied to the approximation of geometric representation, displacement field and traction field simultaneously. An improved power series expansion method (M-PSEM) is used ...
1 CitationsSource
#1Yonghong Wu (BIT: Beijing Institute of Technology)H-Index: 63
#2Chunying Dong (BIT: Beijing Institute of Technology)H-Index: 19
Last. H.S. Yang (BIT: Beijing Institute of Technology)H-Index: 3
view all 3 authors...
Abstract In this paper, Finite Element-Indirect Boundary Element (FE-IBE) coupling approach based on the Isogeometric Analysis (IGA) framework is proposed to deal with Acoustic-Structural Interaction (ASI) problems. Numerical simulation of thin-walled structures is carried out by using isogeometric Reissner-Mindlin shell theory. The isogeometric indirect boundary element method (IBEM) is used to simulate the fluid domain, and the coupling formula is established by connecting the primary variable...
1 CitationsSource
#1Yonghong Wu (BIT: Beijing Institute of Technology)H-Index: 63
#2Chunying Dong (BIT: Beijing Institute of Technology)H-Index: 19
Last. H.S. Yang (BIT: Beijing Institute of Technology)H-Index: 3
view all 3 authors...
Abstract This paper introduces two NURBS (Non-Uniform Rational B-Splines)-based Isogeometric FE (Finite Element)-BE (Boundary Element) coupling approaches, i.e. DCM (Directed Coupling Method) and SICM (Symmetric-Iterative Coupling Method), to solve the structural-acoustic interaction problems. IGA (IsoGeometric Analysis) can eliminate the substantial geometrical errors and time-consuming meshing steps. Meanwhile, CAD (Computer Aided Design) only makes use of surface representations to model geom...
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...
9 CitationsSource