Search

Quadratic NURBS-based discretizations of the Galerkin method suffer from volumetric locking when applied to nearly-incompressible linear elasticity. Volumetric locking causes not only smaller displacements than expected, but also large-amplitude spurious oscillations of normal stresses. Continuous-assumed-strain (CAS) elements have been recently introduced to remove membrane locking in quadratic NURBS-based discretizations of linear plane curved Kirchhoff rods (Casquero et al., CMAME, 2022). In this work, we propose two generalizations of CAS elements (named CAS1 and CAS2 elements) to overcome volumetric locking in quadratic NURBS-based discretizations of nearly-incompressible linear elasticity. CAS1 elements linearly interpolate the strains at the knots in each direction for the term in the variational form involving the first Lam\'e parameter while CAS2 elements linearly interpolate the dilatational strains at the knots in each direction. For both element types, a displacement vector with C1 continuity across element boundaries results in assumed strains with C0 continuity across element boundaries. In addition, the implementation of the two locking treatments proposed in this work does not require any additional global or element matrix operations such as matrix inversions or matrix multiplications. The locking treatments are applied at the element level and the nonzero pattern of the global stiffness matrix is preserved. The numerical examples solved in this work show that CAS1 and CAS2 elements, using either two or three Gauss-Legrendre quadrature points per direction, are effective locking treatments since they not only result in more accurate displacements for coarse meshes, but also remove the spurious oscillations of normal stresses., Comment: arXiv admin note: text overlap with arXiv:2311.00101

Quadratic NURBS-based discretizations of the Galerkin method suffer from membrane locking when applied to Kirchhoff-Love shell formulations. Membrane locking causes not only smaller displacements than expected, but also large-amplitude spurious oscillations of the membrane forces. Continuous-assumed-strain (CAS) elements have been recently introduced to remove membrane locking in quadratic NURBS-based discretizations of linear plane curved Kirchhoff rods (Casquero et al., CMAME, 2022). In this work, we generalize CAS elements to vanquish membrane locking in quadratic NURBS-based discretizations of linear Kirchhoff-Love shells. CAS elements bilinearly interpolate the membrane strains at the four corners of each element. Thus, the assumed strains have C0 continuity across element boundaries. To the best of the authors' knowledge, CAS elements are the first assumed-strain treatment to effectively overcome membrane locking in quadratic NURBS-based discretizations of Kirchhoff-Love shells while satisfying the following important characteristics for computational efficiency: (1) No additional degrees of freedom are added, (2) No additional systems of algebraic equations need to be solved, (3) No matrix multiplications or matrix inversions are needed to obtain the stiffness matrix, and (4) The nonzero pattern of the stiffness matrix is preserved. The benchmark problems show that CAS elements, using either 2x2 or 3x3 Gauss-Legendre quadrature points per element, are an effective locking treatment since this element type results in more accurate displacements for coarse meshes and excises the spurious oscillations of the membrane forces. The benchmark problems also show that CAS elements outperform state-of-the-art element types based on Lagrange polynomials equipped with either assumed-strain or reduced-integration locking treatments.

This paper presents the application of triangle configuration B-splines (TCB-splines) for representing and analyzing the Kirchhoff-Love shell in the context of isogeometric analysis (IGA). The Kirchhoff-Love shell formulation requires global $C^1$-continuous basis functions. The nonuniform rational B-spline (NURBS)-based IGA has been extensively used for developing Kirchhoff-Love shell elements. However, shells with complex geometries inevitably need multiple patches and trimming techniques, where stitching patches with high continuity is a challenge. On the other hand, due to their unstructured nature, TCB-splines can accommodate general polygonal domains, have local refinement, and are flexible to model complex geometries with $C^1$ continuity, which naturally fit into the Kirchhoff-Love shell formulation with complex geometries. Therefore, we propose to use TCB-splines as basis functions for geometric representation and solution approximation. We apply our method to both linear and nonlinear benchmark shell problems, where the accuracy and robustness are validated. The applicability of the proposed approach to shell analysis is further exemplified by performing geometrically nonlinear Kirchhoff-Love shell simulations of a pipe junction and a front bumper represented by a single patch of TCB-splines.

G-splines are a generalization of B-splines that deals with extraordinary points by imposing G^1 constraints across their spoke edges, thus obtaining a continuous tangent plane throughout the surface. Using the isoparametric concept and the Bubnov-Galerkin method to solve partial differential equations with G-splines results in discretizations with global C^1 continuity in physical space. Extraordinary points (EPs) are required to represent manifold surfaces with arbitrary topological genus. In this work, we allow both interior and boundary EPs and there are no limitations regarding how close EPs can be from each other. Reaching this level of flexibility is necessary so that splines with EPs can become mainstream in the design-through-analysis cycle of complex thin-walled structures. To the authors' knowledge, the two EP constructions based on imposing G^1 constraints proposed in this work are the first two EP constructions used in isogeometric analysis (IGA) that combine the following distinctive characteristics: (1) Only vertex-based control points are used and they behave as geometric shape handles, (2) any control point of the control net can potentially be an EP, (3) global C^1 continuity in physical space is obtained without introducing singularities, (4) faces around EPs are not split into multiple elements, and (5) good surface quality is attained. The studies of convergence and surface quality performed in this paper suggest that G-splines are more suitable for IGA than EP constructions based on the D-patch framework. Finally, we have represented the stiffener, the inner part, and the outer part of a B-pillar with G-spline surfaces and solved eigenvalue problems using both Kirchhoff-Love and Reissner-Mindlin shell theories. The results are compared with bilinear quadrilateral meshes and excellent agreement is found between G-splines and conventional finite elements.

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.

Analysis-suitable T-splines (AST-splines) are a promising candidate to achieve a seamless integration between the design and the analysis of thin-walled structures in industrial settings. In this work, we generalize AST-splines to allow multiple extraordinary points within the same face. This generalization drastically increases the flexibility to build geometries using AST-splines; e.g., much coarser meshes can be generated to represent a certain geometry. The AST-spline spaces detailed in this work have $C^1$ inter-element continuity near extraordinary points and $C^2$ inter-element continuity elsewhere. We mathematically show that AST-splines with multiple extraordinary points per face are linearly independent and their polynomial basis functions form a non-negative partition of unity. We numerically show that AST-splines with multiple extraordinary points per face lead to optimal convergence rates for second- and fourth-order linear elliptic problems. To illustrate a possible isogeometric framework that is already available, we design the B-pillar and the side outer panel of a car using T-splines with the commercial software Autodesk Fusion360, import the control nets into our in-house code to build AST-splines, and import the B\'ezier extraction information into the commercial software LS-DYNA to solve eigenvalue problems. The results are compared with conventional finite elements. Good agreement is found, but conventional finite elements require significantly more degrees of freedom to reach a converged solution than AST-splines.

We extend the recently introduced divergence-conforming immersed boundary (DCIB) method [1] to fluid-structure interaction (FSI) problems involving closed co-dimension one solids. We focus on capsules and vesicles, whose discretization is particularly challenging due to the higher-order derivatives that appear in their formulations. In two-dimensional settings, we employ cubic B-splines with periodic knot vectors to obtain discretizations of closed curves with C^2 inter-element continuity. In three-dimensional settings, we use analysis-suitable bi-cubic T-splines to obtain discretizations of closed surfaces with at least C^1 inter-element continuity. Large spurious changes of the fluid volume inside closed co-dimension one solids is a well-known issue for IB methods. The DCIB method results in volume changes orders of magnitude lower than conventional IB methods. This is a byproduct of discretizing the velocity-pressure pair with divergence-conforming B-splines, which lead to negligible incompressibility errors at the Eulerian level. The higher inter-element continuity of divergence-conforming B-splines is also crucial to avoid the quadrature/interpolation errors of IB methods becoming the dominant discretization error. Benchmark and application problems of vesicle and capsule dynamics are solved, including mesh-independence studies and comparisons with other numerical methods., Comment: For supplementary movies go to https://www.andrew.cmu.edu/user/hugocp/research.html