Isogeometric dual mortar methods for computational contact mechanics

In recent years, isogeometric analysis (IGA) has received great attention in many fields of computational mechanics research. Especially for computational contact mechanics, an exact and smooth surface representation is highly desirable. As a consequence, many well-known finite element methods and algorithms for contact mechanics have been transferred to IGA. In the present contribution, the so-called dual mortar method is investigated for both contact mechanics and classical domain decomposition using NURBS basis functions. In contrast to standard mortar methods, the use of dual basis functions for the Lagrange multiplier based on the mathematical concept of biorthogonality enables an easy elimination of the additional Lagrange multiplier degrees of freedom from the global system. This condensed system is smaller in size, and no longer of saddle point type but positive definite. A very simple and commonly used element-wise construction of the dual basis functions is directly transferred to the IGA case. The resulting Lagrange multiplier interpolation satisfies discrete inf-sup stability and biorthogonality, however, the reproduction order is limited to one. In the domain decomposition case, this results in a limitation of the spatial convergence order to O(h32) in the energy norm, whereas for unilateral contact, due to the lower regularity of the solution, optimal convergence rates are still met. Numerical examples are presented that illustrate these theoretical considerations on convergence rates and compare the newly developed isogeometric dual mortar contact formulation with its standard mortar counterpart as well as classical finite elements based on first and second order Lagrange polynomials.