A new mixed isogeometric approach to Kirchhoff-Love shells

Katharina Rafetseder,Walter Zulehner

DOI: https://doi.org/10.1016/j.cma.2018.11.033

2018-03-07

Abstract:For Kichhoff-Love shell problems a new mixed formulation solely based on standard $H^1$ spaces is presented. This allows for flexibility in the construction of discretization spaces, e.g., standard $C^0$-coupling of multi-patch isogeometric spaces is sufficient. In terms of solution strategies, for iterative solvers efficient methods for standard second-order problems like multigrid can be used as building blocks of a preconditioner. Furthermore, a combination of the proposed mixed formulation of the bending part with a popular mixed formulation of the membrane part in order to avoid membrane locking is considered. The performance of both mixed formulations is demonstrated by numerical benchmark studies.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to construct a new hybrid isogeometric method in Kirchhoff - Love shell structure analysis to improve the efficiency and accuracy of numerical simulation. Specifically, the paper mainly focuses on the following aspects:

1. **Hybrid variational formulation**:
- A new hybrid variational formulation is proposed, in which the bending moment tensor \( \mathbf{M} \) is introduced as a new unknown quantity. This new hybrid formulation is based only on the standard \( H^1 \) space and allows the use of standard \( C^0 \) coupling between multi - patch isogeometric spaces, thus simplifying the construction of the discretization space.

2. **Avoiding membrane locking**:
- The paper considers combining the newly proposed hybrid formulation for the bending part with the popular hybrid formulation for the membrane part to avoid the membrane - locking phenomenon. Membrane locking refers to the situation where, in some cases, the shell structure will exhibit unnecessary membrane strains in a pure bending state, which will affect the accuracy of numerical results.

3. **Numerical experiments**:
- The performance of the new hybrid formulation and the combined method is verified through numerical benchmark tests. These tests include three benchmark problems in the classical shell obstacle course: the Scordelis - Lo roof, the pinched hemisphere, and the pinched cylinder. These tests not only verify the effectiveness of the new method but also show its convergence behavior at different polynomial orders.

### Main contributions
- **New hybrid variational formulation**: A new hybrid variational formulation is proposed. This formulation is based only on the standard \( H^1 \) space and allows the use of standard \( C^0 \) coupling to deal with the continuity problem between multi - patch isogeometric spaces.
- **Avoiding membrane locking**: By combining the hybrid formulation of the membrane part, the membrane - locking phenomenon is effectively avoided, and the accuracy of numerical simulation is improved.
- **Numerical verification**: The effectiveness and convergence of the new method are verified through a series of numerical experiments, showing its performance at different polynomial orders.

### Mathematical formulas
- **Bending strain tensor**:
\[
\kappa_{\alpha\beta}(u) = u_3|\alpha\beta - B^\sigma_\alpha B_{\sigma\beta} u_3 + B^\sigma_\alpha u_\sigma|\beta + B^\tau_\beta u_\tau|\alpha + B^\tau_{\beta|\alpha} u_\tau
\]
where \( u_3|\alpha\beta \) represents the second - order covariant derivative, and \( B^\tau_{\beta|\alpha} \) represents the covariant derivative of the second fundamental form.

- **Membrane strain tensor**:
\[
\epsilon_{\alpha\beta}(u) = \frac{1}{2} (u_\alpha|\beta + u_\beta|\alpha) - B_{\alpha\beta} u_3
\]

- **Hybrid variational formulation**:
\[
\int_\Omega (\hat{\mathbf{C}}^{-1}_M : \mathbf{M} : \mathbf{L}) \sqrt{a} \, d\xi - \int_\Omega \kappa(u) : \mathbf{L} \, d\xi = 0
\]
\[
-\int_\Omega \mathbf{M} : \kappa(v) \, d\xi - c(u, v) = -\langle F, v \rangle
\]
where \( \hat{\mathbf{C}}^{-1}_M=\sqrt{a}\frac{12}{t^3}\mathbf{C}^{-1} \), and \( c(u, v) \) represents the contribution of the membrane part.

### Conclusion
This paper solves several key problems in Kirchhoff - Love shell structure analysis by proposing a new hybrid isogeometric method, especially...