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

TL;DR

This paper introduces continuous-assumed-strain (CAS) elements for quadratic NURBS discretizations of linear plane Kirchhoff rods, effectively removing membrane locking across various slenderness ratios with minimal computational overhead.

Contribution

The authors develop and validate CAS elements that eliminate membrane locking in NURBS-based discretizations without adding degrees of freedom or increasing computational cost.

Findings

CAS elements remove membrane locking for a broad range of slenderness ratios.

Convergence of CAS elements is independent of slenderness ratio.

CAS elements reduce spurious stress oscillations caused by membrane locking.

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.