Accurate equilibrium-based interlaminar stress recovery for isogeometric laminated composite Kirchhoff plates

Alessia Patton,Pablo Antolin,John-Eric Dufour,Josef Kiendl,Alessandro Reali

DOI: https://doi.org/10.48550/arXiv.2007.10029

2020-07-20

Abstract:In this paper, we use isogeometric Kirchhoff plates to approximate composite laminates adopting the classical laminate plate theory. Both isogeometric Galerkin and collocation formulations are considered. Within this framework, interlaminar stresses are recovered through an effective post-processing technique based on the direct imposition of equilibrium in strong form, relying on the accuracy and the higher continuity typically granted by isogeometric discretizations. The effectiveness of the proposed approach is proven by extensive numerical tests.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in the isogeometric analysis framework based on the Classical Laminated Plate Theory (CLPT) in composite laminates, accurately recover the interlaminar stresses through an effective post - processing technique. Specifically, the paper aims to:

1. **Accurately evaluate the three - dimensional stress state**: Due to the mismatch of material properties between different layers, laminated composites may exhibit complex mechanical behaviors under external loads, and even lead to delamination failure. Therefore, accurately evaluating the three - dimensional stress state inside the laminate is crucial for designing and evaluating its structural response.

2. **Overcome the limitations of existing theories**:
- **Equivalent Single - Layer Theory (ESL)**: Although it is simple to calculate and has low cost, it cannot provide sufficient accuracy to describe the inter - layer stresses, especially near the edges and notch boundaries.
- **Layer - Wise Theory (LW)**: It can approximate the inter - layer stresses relatively well, but as the number of layers increases, the calculation cost rises significantly.

3. **Combine the advantages of Isogeometric Analysis (IGA)**: Utilize the high - order continuity and high precision provided by isogeometric analysis, and accurately recover the out - of - plane stress distribution from the rough two - dimensional displacement solution by directly applying the equilibrium conditions in the strong form.

### Specific problems
- **How to improve the prediction accuracy of the interface stress in laminated composite plates while maintaining computational efficiency?**
- **How to ensure the accurate recovery of the three - dimensional stress field from the rough two - dimensional displacement solution by using the high - order continuous shape functions in the isogeometric analysis framework?**

### Solution
The paper proposes a post - processing technique based on equilibrium conditions. This technique can accurately recover the three - dimensional stress distribution between layers from the two - dimensional displacement solution obtained by the Classical Kirchhoff Laminated Plate Theory (CLPT) by directly applying the equilibrium equations in the strong form. This method not only improves the computational efficiency but also ensures the accuracy of the results.

### Key formulas
- **Displacement field**:
\[
u_1(x_1, x_2, x_3)=-x_3 w,_{1}, \quad u_2(x_1, x_2, x_3)=-x_3 w,_{2}, \quad u_3(x_1, x_2)=w
\]
where \(w\) is the transverse deflection.

- **Bending moment**:
\[
M_{\alpha\beta}=-\int_{-\frac{t}{2}}^{\frac{t}{2}} x_3 C_{\alpha\beta\gamma\delta}(x_3) \kappa_{\gamma\delta} \, dx_3
\]
where \(C_{\alpha\beta\gamma\delta}\) is the elastic tensor and \(\kappa_{\gamma\delta}\) is the curvature.

- **Equilibrium equation**:
\[
\nabla \cdot \sigma + b = 0
\]
where \(\sigma\) is the stress tensor and \(b\) is the body force.

Through these formulas and methods, the paper shows how to accurately recover the interface stress distribution in laminated composite plates while maintaining computational efficiency.