Dual and approximate dual basis functions for B-splines and NURBS – Comparison and application for an efficient coupling of patches with the isogeometric mortar method

Abstract

This contribution defines and compares different methods for the computation of dual basis functions for B-splines and Non-Uniform Rational B-splines (NURBS). They are intended to be used as test functions for the isogeometric mortar method, but other fields of application are possible, too. Three different concepts are presented and compared. The first concept is the explicit formula for the computation of dual basis functions for NURBS proposed in the work of Carl de Boor. These dual basis functions entail minimal support, i.e., the support of the dual basis functions is equal to the support of the corresponding B-spline basis functions. In the second concept dual basis functions are derived from the inversion of the Gram matrix. These dual basis functions have global support along the interface. The third concept is the use of approximate dual basis functions, which were initially proposed for the use in harmonic analysis. The support of these functions is local but larger than the support of the associated B-spline basis functions. We propose an extension of the approximate dual basis functions for NURBS basis functions. After providing the general formulas, we elaborate explicit expressions for several degrees of spline basis functions. All three approaches are applied in the frame of the mortar method for the coupling of non-conforming NURBS patches. A method which allows complex discretizations with multiple intersecting interfaces is presented. Numerical examples show that the explicitly defined dual basis functions with minimal support severely deteriorate the global stress convergence behavior of the mechanical analysis. This fact is in accordance with mathematical findings in literature, which state that the optimal reproduction degree of arbitrary functions is not possible without extending the support of the dual basis functions. The dual basis functions computed from the inverse of the Gram matrix yield accurate numerical results but the global support yields significantly higher computational costs in comparison to computations of conforming meshes. Only the approximate dual basis functions yield accurate and efficient computations, where neither accuracy nor efficiency is significantly deteriorated in comparison to computations of conforming meshes. All basic cases of T-intersections and star-intersections are studied. Furthermore, an example which combines all basic cases in a complex discretization is given. The applicability of the presented method for the nonlinear case and for shell formulations is shown with the help of one numerical example.