Removing membrane locking in quadratic NURBS-based discretizations of linear plane Kirchhoff rods: CAS elements

Abstract

NURBS-based discretizations suffer from membrane locking when applied to primal formulations of curved thin-walled structures. We consider linear plane curved Kirchhoff rods as a model problem to study how to remove membrane locking from NURBS-based discretizations. In this work, we propose continuous-assumed-strain (CAS) elements, an assumed strain treatment that removes membrane locking from quadratic NURBS for an ample range of slenderness ratios. CAS elements take advantage of the C1 inter-element continuity of the displacement vector given by quadratic NURBS to interpolate the membrane strain using linear Lagrange polynomials while preserving the C0 inter-element continuity of the membrane strain. CAS elements are the first NURBS-based element type able to remove membrane locking for a broad range of slenderness ratios that combines the following characteristics: (1) No additional degrees of freedom are added, (2) No additional systems of algebraic equations need to be solved, and (3) The nonzero pattern of the stiffness matrix is preserved. Since the only additional computations required by the proposed element type are to evaluate the derivatives of the basis functions and the unit tangent vector at the knots, the proposed scheme barely increases the computational cost with respect to the locking-prone NURBS-based discretization of the primal formulation. The benchmark problems show that the convergence of CAS elements is independent of the slenderness ratio while the convergence of quadratic NURBS elements, local Bbar elements, and local ANS elements depends heavily on the slenderness ratio. The numerical examples also show how CAS elements remove the spurious oscillations in stress resultants caused by membrane locking while quadratic NURBS elements, local Bbar elements, and local ANS elements suffer from large-amplitude spurious oscillations in stress resultants.

References in corpus (4)

• Bézier projection: a unified approach for local projection and quadrature-free refinement and coarsening of NURBS and T-splines with particular application to isogeometric design and analysis
• Analysis-suitable unstructured T-splines: Multiple extraordinary points per face
• Isogeometric Bézier dual mortaring: Refineable higher-order spline dual bases and weakly continuous geometry
• Bézier $\bar{B}$ Projection

Cited by in corpus (7)

• Overcoming membrane locking in quadratic NURBS-based discretizations of linear Kirchhoff-Love shells: CAS elements
• A simple and efficient hybrid discretization approach to alleviate membrane locking in isogeometric thin shells
• Automated shape and thickness optimization for non-matching isogeometric shells using free-form deformation
• 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
• A new anisotropic bending model for nonlinear shells: Comparison with existing models and isogeometric finite element implementation
• A dual lumping procedure for static condensation in mixed NURBS-based isogeometric elements with optimal convergence rates for arbitrary open knot vectors