Abstract:We consider and discretize a mixed formulation for linear elasticity with weakly imposed symmetry in two and three dimensions. Whereas existing methods mainly deal with simplicial or polygonal meshes, we take advantage of isogeometric analysis (IGA) and consequently allow for shapes with curved boundaries. To introduce the discrete spaces we use isogeometric discrete differential forms defined by proper B-spline spaces. For the proposed schemes a proof of well-posedness and an error estimate are given. Further we discuss our ansatz by means of different numerical examples.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a Mixed Isogeometric Methods for linear elasticity problems with weak symmetry constraints in two - and three - dimensional spaces. Specifically, the author aims to:
1. **Introduce discrete spaces**: Use appropriate B - spline spaces to define isogeometric discrete differential forms in order to handle shapes with curved - edge boundaries.
2. **Prove the goodness of the scheme**: Provide a proof of the goodness of the proposed numerical scheme and give error estimates.
3. **Numerical verification**: Discuss and verify the proposed method through different numerical examples.
### Detailed Explanation
#### 1. Introduction of Discrete Spaces
Traditional finite element methods mainly use simplex or polygonal meshes, while this paper takes advantage of Isogeometric Analysis (IGA) to allow handling of shapes with curved - edge boundaries. To this end, the author introduces isogeometric discrete differential forms defined by appropriate B - spline spaces.
#### 2. Prove the Goodness of the Scheme
To ensure that the proposed numerical scheme is well - defined, the author provides the following two key results:
- **Well - posedness proof**: Verify an inf - sup stability result through the macro - element technique, thereby proving the well - posedness of the scheme.
- **Error estimate**: Give an error estimate to evaluate the difference between the numerical solution and the exact solution.
#### 3. Numerical Verification
The author verifies the proposed method through multiple numerical examples, demonstrating its effectiveness and accuracy in practical applications.
### Mathematical Background
This paper studies linear elasticity problems based on the mixed variational formulation and weak symmetry conditions. Specifically, the author considers the following governing equations:
\[
\begin{aligned}
A\sigma &= \varepsilon(u) & \text{in } \Omega, \\
\nabla \cdot \sigma &= f & \text{in } \Omega,
\end{aligned}
\]
where:
- \( A \) is the compliance tensor, which describes material properties;
- \( \sigma \) is the symmetric stress tensor;
- \( u \) is the displacement field;
- \( f \) is the external force distribution;
- \( \varepsilon(u)=\frac{1}{2}(\nabla u + (\nabla u)^T) \) is the symmetric gradient, corresponding to the mechanical strain in the linear theory.
To introduce the weak symmetry condition, the author proposes the following mixed weak form:
\[
\begin{aligned}
\langle A\sigma, \tau \rangle + \langle u, \nabla \cdot \tau \rangle + \langle p, \text{Skew}(\tau) \rangle &= \langle u_D, \tau \cdot \nu \rangle_{\Gamma}, & \forall \tau \in H(\Omega, \text{div}, M), \\
\langle \nabla \cdot \sigma, v \rangle &= \langle f, v \rangle, & \forall v \in L^2(\Omega, \mathbb{R}^n), \\
\langle \text{Skew}(\sigma), q \rangle &= 0, & \forall q \in L^2(\Omega, \mathbb{R}^{s(n)}).
\end{aligned}
\]
where:
- \( \text{Skew}(M) \) represents the skew - symmetric part of the matrix \( M \);
- \( s(n) = 2n - 3 \) is the dimension of the skew - symmetric matrix.
Through this method, the author avoids the limitations of the strong symmetry condition, thereby obtaining a more flexible discretization method.
### Summary
The main contribution of this paper lies in the development.