Adaptive isogeometric analysis on two-dimensional trimmed domains based on a hierarchical approach

The focus of this work is on the development of an error-driven isogeometric framework, capable of automatically performing an adaptive simulation in the context of second- and fourth-order, elliptic partial differential equations defined on two-dimensional trimmed domains. The method is steered by an a posteriori error estimator, which is computed with the aid of an auxiliary residual-like problem formulated onto a space spanned by splines with single element support. The local refinement of the basis is achieved thanks to the use of truncated hierarchical B-splines. We prove numerically the applicability of the proposed estimator to various engineering-relevant problems, namely the Poisson problem, linear elasticity and Kirchhoff-Love shells, formulated on trimmed geometries. In particular, we study several benchmark problems which exhibit both smooth and singular solutions, where we recover optimal asymptotic rates of convergence for the error measured in the energy norm and we observe a substantial increase in accuracy per-degree-of-freedom compared to uniform refinement. Lastly, we show the applicability of our framework to the adaptive shell analysis of an industrial-like trimmed geometry modeled in the commercial software Rhinoceros, which represents the B-pillar of a car.

[1]  Markus Kästner,et al.  Projection and transfer operators in adaptive isogeometric analysis with hierarchical B-splines , 2018, Computer Methods in Applied Mechanics and Engineering.

[2]  K. Höllig Finite element methods with B-splines , 1987 .

[3]  Roland Wüchner,et al.  Isogeometric shell analysis with Kirchhoff–Love elements , 2009 .

[4]  Günther Greiner,et al.  Interpolating and approximating scattered 3D-data with hierarchical tensor product B-splines , 2010 .

[5]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis , 2000 .

[6]  Roland Wüchner,et al.  Analysis in computer aided design: Nonlinear isogeometric B-Rep analysis of shell structures , 2015 .

[7]  Hendrik Speleers,et al.  THB-splines: The truncated basis for hierarchical splines , 2012, Comput. Aided Geom. Des..

[8]  Annalisa Buffa,et al.  A hierarchical approach to the a posteriori error estimation of isogeometric Kirchhoff plates and Kirchhoff-Love shells , 2019, Computer Methods in Applied Mechanics and Engineering.

[9]  Giancarlo Sangalli,et al.  ANALYSIS-SUITABLE T-SPLINES OF ARBITRARY DEGREE: DEFINITION, LINEAR INDEPENDENCE AND APPROXIMATION PROPERTIES , 2013 .

[10]  Jiansong Deng,et al.  Solving higher order PDEs with isogeometric analysis on implicit domains using weighted extended THB-splines , 2019, Comput. Aided Geom. Des..

[11]  Klaus Höllig,et al.  Finite Element Analysis with B-Splines: Weighted and Isogeometric Methods , 2010, Curves and Surfaces.

[12]  B. Simeon,et al.  Adaptive isogeometric analysis by local h-refinement with T-splines , 2010 .

[13]  Rafael Vázquez,et al.  Algorithms for the implementation of adaptive isogeometric methods using hierarchical splines , 2016 .

[14]  B. Simeon,et al.  A hierarchical approach to adaptive local refinement in isogeometric analysis , 2011 .

[15]  Bert Jüttler,et al.  THB-splines: An effective mathematical technology for adaptive refinement in geometric design and isogeometric analysis , 2016 .

[16]  Andreas Keller,et al.  Finite Element Approximation with Hierarchical B-Splines , 2014, Curves and Surfaces.

[17]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[18]  Massimiliano Martinelli,et al.  Isogeometric Analysis on V-reps: First results , 2019, Computer Methods in Applied Mechanics and Engineering.

[19]  Randolph E. Bank,et al.  A posteriori error estimates based on hierarchical bases , 1993 .

[20]  Giancarlo Sangalli,et al.  Mathematical analysis of variational isogeometric methods* , 2014, Acta Numerica.

[21]  Annalisa Buffa,et al.  A Minimal Stabilization Procedure for Isogeometric Methods on Trimmed Geometries , 2019, SIAM J. Numer. Anal..

[22]  T. Hughes,et al.  ISOGEOMETRIC ANALYSIS: APPROXIMATION, STABILITY AND ERROR ESTIMATES FOR h-REFINED MESHES , 2006 .

[23]  Andrea Bressan,et al.  Some properties of LR-splines , 2013, Comput. Aided Geom. Des..

[24]  Carlotta Giannelli,et al.  Complexity of hierarchical refinement for a class of admissible mesh configurations , 2015, Comput. Aided Geom. Des..

[25]  Thomas J. R. Hughes,et al.  Hierarchically refined and coarsened splines for moving interface problems, with particular application to phase-field models of prostate tumor growth , 2017 .

[26]  David R. Forsey,et al.  Hierarchical B-spline refinement , 1988, SIGGRAPH.

[27]  Rafael Vázquez,et al.  A new design for the implementation of isogeometric analysis in Octave and Matlab: GeoPDEs 3.0 , 2016, Comput. Math. Appl..

[28]  Thomas J. R. Hughes,et al.  Variationally consistent isogeometric analysis of trimmed thin shells at finite deformations, based on the STEP exchange format , 2018, Computer Methods in Applied Mechanics and Engineering.

[29]  P. Hansbo,et al.  Fictitious domain finite element methods using cut elements , 2012 .

[30]  Alessandro Reali,et al.  Suitably graded THB-spline refinement and coarsening: Towards an adaptive isogeometric analysis of additive manufacturing processes , 2018, Computer Methods in Applied Mechanics and Engineering.

[31]  Hendrik Speleers,et al.  Strongly stable bases for adaptively refined multilevel spline spaces , 2014, Adv. Comput. Math..

[32]  Wing Kam Liu,et al.  Stress projection for membrane and shear locking in shell finite elements , 1985 .

[33]  Carlotta Giannelli,et al.  Refinement Algorithms for Adaptive Isogeometric Methods with Hierarchical Splines , 2018, Axioms.

[34]  T. Hughes,et al.  Local refinement of analysis-suitable T-splines , 2012 .

[35]  Carlotta Giannelli,et al.  Adaptive isogeometric methods with hierarchical splines: Optimality and convergence rates , 2017 .

[36]  Tom Lyche,et al.  Polynomial splines over locally refined box-partitions , 2013, Comput. Aided Geom. Des..

[37]  Ernst Rank,et al.  Efficient and accurate numerical quadrature for immersed boundary methods , 2015, Advanced Modeling and Simulation in Engineering Sciences.

[38]  Benjamin Marussig,et al.  Improved conditioning of isogeometric analysis matrices for trimmed geometries , 2018, Computer Methods in Applied Mechanics and Engineering.

[39]  Benjamin Marussig,et al.  A Review of Trimming in Isogeometric Analysis: Challenges, Data Exchange and Simulation Aspects , 2017, Archives of Computational Methods in Engineering.

[40]  John A. Evans,et al.  Isogeometric Analysis , 2010 .

[41]  F. de Prenter,et al.  Condition number analysis and preconditioning of the finite cell method , 2016, 1601.05129.

[42]  John A. Evans,et al.  Isogeometric analysis using T-splines , 2010 .

[43]  F. de Prenter,et al.  Multigrid solvers for immersed finite element methods and immersed isogeometric analysis , 2019, Computational Mechanics.

[44]  P. Gould Introduction to Linear Elasticity , 1983 .