A dual lumping procedure for static condensation in mixed NURBS-based isogeometric elements with optimal convergence rates for arbitrary open knot vectors

Abstract

Locking is a common effect in finite element and isogeometric analysis. In the case of plates, transverse shear locking is most prominent, for shells several other types of locking exist. A common cure are mixed methods that introduce additional fields of unknowns into the variational formulation. These fields reduce constraints and thus alleviate locking significantly. As a drawback, the discretized additional fields increase computational costs significantly. These fields are often eliminated by static condensation, which requires the inverse of a part of the stiffness matrix. In Lagrange-based finite elements, this inverse is computed on element level, due to a discontinuous interpolation of additional fields. Since isogeometric analysis features higher continuity, static condensation must be performed on patch level, which requires a costly matrix inversion on that level. In this contribution, the virtual shear parameters of a mixed isogeometric plate formulation are interpolated by enhanced approximate dual basis functions. This allows to conduct row-sum lumping of the relevant matrix part at a minimal loss of accuracy, since this part becomes diagonal dominant. For a properly chosen integration space, this lumped matrix becomes the identity matrix. Thus, the proposed condensation procedure does not require an inversion anymore. The crucial and novel point is the proposed treatment of knot vectors with limited internal continuity. With the help of several single- and multi-patch examples, both with full and with limited internal continuity, we show that the proposed procedure obtains optimal error convergence rates in all cases, while without these alterations, convergence rates are significantly deteriorated.

arXiv comment52 pages, 18 figures, 2 tables, Preprint submitted to Computer Methods in Applied Mechanics and Engineering

References in corpus (10)

• An isogeometric method for the Reissner-Mindlin plate bending problem
• Isogeometric Bézier dual mortaring: Refineable higher-order spline dual bases and weakly continuous geometry
• Towards higher-order accurate mass lumping in explicit isogeometric analysis for structural dynamics
• Biorthogonal splines for optimal weak patch-coupling in isogeometric analysis with applications to finite deformation elasticity
• Overcoming membrane locking in quadratic NURBS-based discretizations of linear Kirchhoff-Love shells: CAS elements
• Removing membrane locking in quadratic NURBS-based discretizations of linear plane Kirchhoff rods: CAS elements
• Bézier $\bar{B}$ Projection
• Higher order accurate mass lumping for explicit isogeometric methods based on approximate dual basis functions
• Vanquishing volumetric locking in quadratic NURBS-based discretizations of nearly-incompressible linear elasticity: CAS elements
• An objective isogeometric mixed finite element formulation for nonlinear elastodynamic beams with incompatible warping strains