Mirza Cenanovic,Peter Hansbo,Mats G. Larson
Abstract:We construct a cut finite element method for the membrane elasticity problem on an embedded mesh using tangential differential calculus. Both free membranes and membranes coupled to 3D elasticity are considered. The discretization comes from a Galerkin method using the restriction of 3D basis funtions (linear or trilinear) to the surface representing the membrane. In the case of coupling to 3D elasticity, we view the membrane as giving additional stiffness contributions to the standard stiffness matrix resulting from the discretization of the three-dimensional continuum.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a finite - element modeling method for linear membrane problems, especially when these membranes are embedded in three - dimensional space. Specifically, the author proposes a "cut finite element method", which allows the membrane structure to pass through the pre - existing three - dimensional grid without requiring these grids to be aligned with the membrane surface. Through this method, free membranes and membrane problems coupled with three - dimensional elastomers can be dealt with more flexibly.
### Main problems and solutions:
1. **Problem description**:
- Linear elastic membrane problems usually involve two - dimensional membrane structures embedded in three - dimensional space. Traditional finite element methods require the grid to be precisely aligned with the membrane surface, which can be very complicated and time - consuming in practical applications.
- The author focuses on how to effectively simulate the behavior of these membrane structures without changing the existing three - dimensional grid.
2. **Solution**:
- **Cut finite element method**: By using high - dimensional shape functions and restricting them to the membrane surface, the author proposes a new finite - element modeling method. This method allows the membrane structure to cut through three - dimensional grid elements in an arbitrary way, thus greatly improving the flexibility of modeling.
- **Stabilization technique**: To deal with the numerical instability problems caused by cutting, the author introduces a stabilization technique, which ensures a good condition number of the discrete system and allows arbitrarily small cuts in the three - dimensional grid.
- **Coupling to three - dimensional elastomers**: This method is not only applicable to free membranes, but also can be used for membrane structures coupled with three - dimensional elastomers, such as reinforcing materials or adhesive layers. This coupling is achieved by adding the stiffness contribution of the membrane to the standard three - dimensional stiffness matrix.
### Mathematical models and formulas:
- **Tangential gradient**:
\[
\nabla_\Gamma u = P_\Gamma \nabla u
\]
where \(\nabla\) is the three - dimensional gradient operator, and \(P_\Gamma = I - n_\Gamma \otimes n_\Gamma\) is the orthogonal projection from \(\mathbb{R}^3\) to the tangential plane of the membrane surface.
- **Strain tensor**:
\[
\varepsilon_\Gamma(u) = P_\Gamma \varepsilon(u) P_\Gamma
\]
where \(\varepsilon(u)=\frac{1}{2}(\nabla u + (\nabla u)^T)\) is the three - dimensional strain tensor.
- **Stress tensor**:
\[
\sigma_\Gamma = 2\mu \varepsilon_\Gamma+\lambda \text{tr}(\varepsilon_\Gamma) P_\Gamma
\]
where \(\mu\) and \(\lambda\) are Lamé parameters respectively, defined as:
\[
\mu=\frac{E}{2(1 + \nu)}, \quad \lambda=\frac{E\nu}{1-\nu^2}
\]
\(E\) is Young's modulus and \(\nu\) is Poisson's ratio.
### Conclusion:
By introducing the cut finite element method and its stabilization technique, this paper provides an effective and flexible tool for simulating membrane structures embedded in three - dimensional elastomers. This method can be applied not only to free membranes, but also can be extended to complex structures coupled with three - dimensional elastomers, such as reinforcing materials or adhesive layers. In addition, numerical experiments show that the error level of this method is comparable to that of the traditional triangulated membrane model, but it is more efficient in dealing with complex geometries.