The INTERNODES method for the treatment of non-conforming multipatch geometries in Isogeometric Analysis

In this paper we apply the INTERNODES method to solve second order elliptic problems discretized by Isogeometric Analysis methods on non-conforming multiple patches in 2D and 3D geometries. INTERNODES is an interpolation-based method that, on each interface of the configuration, exploits two independent interpolation operators to enforce the continuity of the traces and of the normal derivatives. INTERNODES supports non-conformity on NURBS spaces as well as on geometries. We specify how to set up the interpolation matrices on non-conforming interfaces, how to enforce the continuity of the normal derivatives and we give special attention to implementation aspects. The numerical results show that INTERNODES exhibits optimal convergence rate with respect to the mesh size of the NURBS spaces an that it is robust with respect to jumping coefficients.Comment: Accepted for publication in Computer Methods in Applied Mechanics and Engineerin

Bias extension test for pantographic sheets: numerical simulations based on second gradient shear energies

We consider a bi-dimensional sheet consisting of two orthogonal families of inextensible fibres. Using the representation due to Rivlin and Pipkin for admissible placements, i.e. placements preserving the lengths of the inextensible fibres, we numerically simulate a standard bias extension test on the sheet, solving a non-linear constrained optimization problem. Several first and second gradient deformation energy models are considered, depending on the shear angle between the fibres and on its gradient, and the results obtained are compared. The proposed numerical simulations will be helpful in designing a systematic experimental campaign aimed at characterizing the internal energy for physical realizations of the ideal pantographic structure presented in this paper

Efficient mass and stiffness matrix assembly via weighted Gaussian quadrature rules for B-splines

Calabr{\`o} et al. [10] changed the paradigm of the mass and stiffness computation from the traditional element-wise assembly to a row-wise concept, showing that the latter one offers integration that may be orders of magnitude faster. Considering a B-spline basis function as a non-negative measure, each mass matrix row is integrated by its own quadrature rule with respect to that measure. Each rule is easy to compute as it leads to a linear system of equations, however, the quadrature rules are of the Newton-Cotes type, that is, they require a number of quadrature points that is equal to the dimension of the spline space. In this work, we propose weighted quadrature rules of Gaussian type which require the minimum number of quadrature points while guaranteeing exactness of integration with respect to the weight function. The weighted Gaussian rules arise as solutions of non-linear systems of equations. We derive rules for the mass and stiffness matrices for uniform $C^1$ quadratic and $C^2$ cubic isogeometric discretizations. In each parameter direction, our rules require locally only $p+1$ quadrature points, $p$ being the polynomial degree. While the nodes cannot be reused for various weight functions as in [10], the computational cost of the mass and stiffness matrix assembly is comparable.RYC-2017-2264