A robust patch coupling method for NURBS-based isogeometric analysis of non-conforming multipatch surfaces

Abstract This paper presents a simple yet flexible method for coupling non-conforming NURBS patches in isogeometric frameworks. It consists of virtually refining bordering patches to derive coupling constraints in a master–slave formulation. These constraints are then enforced via a substitution method for condensation of the slave variables, thereby reducing the model size. In the special case of hierarchical meshes, the method results in an exact connection. For an arbitrary non-conforming configuration, the master–slave formulation takes the interface constraints into account in a least-squares sense. This yields an accurate weak coupling without overconstraining the interface. The coupling relationships only depend on the mesh itself and not on any problem-dependent parameters, allowing them to be generated in a pre-processing step with very limited numerical efforts. Besides this low computational cost, the main merit of the proposed method lies in its simplicity and robustness, yielding good results for arbitrarily strongly non-conforming patch configurations. Several numerical examples are studied for different problem types requiring C 0 - or C 1 -continuity, in particular time-harmonic acoustics and (dynamic) thin plate bending. Benchmarking of the proposed approach against existing similar techniques illustrates its superior accuracy and robustness.

[1]  John A. Evans,et al.  Isogeometric finite element data structures based on Bézier extraction of NURBS , 2011 .

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

[3]  Jiansong Deng,et al.  Polynomial splines over hierarchical T-meshes , 2008, Graph. Model..

[4]  Per Christian Hansen,et al.  Some Applications of the Rank Revealing QR Factorization , 1992, SIAM J. Sci. Comput..

[5]  Anthony T. Patera,et al.  Domain Decomposition by the Mortar Element Method , 1993 .

[6]  I. Babuska,et al.  Finite Element Solution of the Helmholtz Equation with High Wave Number Part II: The h - p Version of the FEM , 1997 .

[7]  Wei Wang,et al.  Nitsche method for isogeometric analysis of Reissner-Mindlin plate with non-conforming multi-patches , 2015, Comput. Aided Geom. Des..

[8]  Barbara Wohlmuth,et al.  Isogeometric mortar methods , 2014, 1407.8313.

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

[10]  Faker Ben Belgacem,et al.  The Mortar finite element method with Lagrange multipliers , 1999, Numerische Mathematik.

[11]  Alessandro Reali,et al.  GeoPDEs: A research tool for Isogeometric Analysis of PDEs , 2011, Adv. Eng. Softw..

[12]  J. W. Brown,et al.  Fourier series and boundary value problems , 1941 .

[13]  Andrew J. Kurdila,et al.  『Fundamentals of Structural Dynamics』(私の一冊) , 2019, Journal of the Society of Mechanical Engineers.

[14]  Les A. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communication.

[15]  Bert Jüttler,et al.  IETI – Isogeometric Tearing and Interconnecting , 2012, Computer methods in applied mechanics and engineering.

[16]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[17]  M. Scott,et al.  Acoustic isogeometric boundary element analysis , 2014 .

[18]  Saeed Shojaee,et al.  Free vibration analysis of thin plates by using a NURBS-based isogeometric approach , 2012 .

[19]  Ernst Rank,et al.  Weak coupling for isogeometric analysis of non-matching and trimmed multi-patch geometries , 2014 .

[20]  I. Babuska,et al.  Finite element solution of the Helmholtz equation with high wave number Part I: The h-version of the FEM☆ , 1995 .

[21]  Peter Betsch,et al.  Isogeometric analysis and domain decomposition methods , 2012 .

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

[23]  C. D. Boor,et al.  On Calculating B-splines , 1972 .

[24]  Onur Atak,et al.  An isogeometric indirect boundary element method for solving acoustic problems in open-boundary domains , 2017 .

[25]  Louis Jezequel,et al.  A C0/G1 multiple patches connection method in isogeometric analysis , 2015 .

[26]  Roland Wüchner,et al.  A Nitsche‐type formulation and comparison of the most common domain decomposition methods in isogeometric analysis , 2014 .

[27]  Thomas J. R. Hughes,et al.  Truncated hierarchical Catmull–Clark subdivision with local refinement , 2015 .

[28]  Wim Desmet,et al.  A performance study of NURBS-based isogeometric analysis for interior two-dimensional time-harmonic acoustics , 2016 .

[29]  Alessandro Reali,et al.  Studies of Refinement and Continuity in Isogeometric Structural Analysis (Preprint) , 2007 .

[30]  J. Reddy Theory and Analysis of Elastic Plates and Shells , 2006 .

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

[32]  Eugen J. Skudrzyk,et al.  The Foundations of Acoustics: Basic Mathematics and Basic Acoustics , 1972 .

[33]  Yuri Bazilevs,et al.  The bending strip method for isogeometric analysis of Kirchhoff–Love shell structures comprised of multiple patches , 2010 .

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

[35]  Bernd Hamann,et al.  Iso‐geometric Finite Element Analysis Based on Catmull‐Clark : ubdivision Solids , 2010, Comput. Graph. Forum.

[36]  M. Cox The Numerical Evaluation of B-Splines , 1972 .

[37]  Vinh Phu Nguyen,et al.  Isogeometric analysis: An overview and computer implementation aspects , 2012, Math. Comput. Simul..

[38]  P. Bar-Yoseph,et al.  Mechanically based models: Adaptive refinement for B‐spline finite element , 2003 .

[39]  Vinh Phu Nguyen,et al.  Nitsche’s method for two and three dimensional NURBS patch coupling , 2013, 1308.0802.

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

[41]  P. Bouillard,et al.  Error estimation and adaptivity for the finite element method in acoustics: 2D and 3D applications , 1999 .

[42]  Sven Klinkel,et al.  The weak substitution method – an application of the mortar method for patch coupling in NURBS‐based isogeometric analysis , 2015 .

[43]  Roland Wüchner,et al.  Isogeometric analysis of trimmed NURBS geometries , 2012 .