Fast and accurate elastic analysis of laminated composite plates via isogeometric collocation and an equilibrium-based stress recovery approach

Alessia Patton,John-Eric Dufour,Pablo Antolin,Alessandro Reali
DOI: https://doi.org/10.48550/arXiv.1901.08797
2019-01-25
Abstract:A novel approach which combines isogeometric collocation and an equilibrium-based stress recovery technique is applied to analyze laminated composite plates. Isogeometric collocation is an appealing strong form alternative to standard Galerkin approaches, able to achieve high order convergence rates coupled with a significantly reduced computational cost. Laminated composite plates are herein conveniently modeled considering only one element through the thickness with homogenized material properties. This guarantees accurate results in terms of displacements and in-plane stress components. To recover an accurate out-of-plane stress state, equilibrium is imposed in strong form as a post-processing correction step, which requires the shape functions to be highly continuous. This continuity demand is fully granted by isogeometric analysis properties, and excellent results are obtained using a minimal number of collocation points per direction, particularly for increasing values of length-to-thickness plate ratio and number of layers.
Numerical Analysis
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to solve the problem of rapid and accurate elastic analysis of laminated composite plates. Specifically, the author proposes a new method that combines the isogeometric collocation method and the stress - recovery technique based on equilibrium to efficiently analyze the mechanical behavior of laminated composite plates. #### Main problems 1. **Limitations of traditional methods**: - **Two - dimensional theory**: It is usually not sufficient to accurately describe delamination damage and interface failure. - **Layer - by - layer theory**: Although it can provide more accurate results, the computational cost is high, especially when the number of layers increases. 2. **Improving computational efficiency and accuracy**: - **Reducing degrees of freedom**: By using a single element throughout the thickness and adopting homogenized material properties, the number of degrees of freedom is significantly reduced. - **Stress recovery**: In order to obtain an accurate out - of - plane stress state, the equilibrium equation is imposed in a strong form in the post - processing step, ensuring the accuracy of the stress results. #### Method overview - **Isogeometric Collocation (IGA - C)**: This is a strong - form method that replaces the standard Galerkin method and can achieve a high - order convergence rate and significantly reduce the computational cost. - **Homogenized single - element method**: The model is simplified by modeling the entire plate thickness with one element and using homogenized material properties. - **Equilibrium - based stress recovery**: In the post - processing step, the out - of - plane stress components are recovered by directly integrating the equilibrium equation to ensure the accuracy of the stress distribution. #### Key formulas - **Homogenized expression of the elastic tensor**: \[ C_{ij}=\sum_{k = 1}^{N}t_kC_{ij}^{(k)}+\text{correction term} \] where \(C_{ij}\) is the homogenized elastic tensor component, \(t_k\) is the volume fraction of the \(k\)-th layer, and \(C_{ij}^{(k)}\) is the elastic tensor component of the \(k\)-th layer. - **Out - of - plane stress recovery formulas**: \[ \sigma_{13}(X_3)=-\int_{\bar{X}_3}^{X_3}(\sigma_{11,1}(\zeta)+\sigma_{12,2}(\zeta)+b_1(\zeta))d\zeta+\sigma_{13}(\bar{X}_3) \] \[ \sigma_{23}(X_3)=-\int_{\bar{X}_3}^{X_3}(\sigma_{12,1}(\zeta)+\sigma_{22,2}(\zeta)+b_2(\zeta))d\zeta+\sigma_{23}(\bar{X}_3) \] \[ \sigma_{33}(X_3)=-\int_{\bar{X}_3}^{X_3}(\sigma_{13,1}(\zeta)+\sigma_{23,2}(\zeta)+b_3(\zeta))d\zeta+\sigma_{33}(\bar{X}_3) \] Through this method, the paper shows that while maintaining high accuracy, the computational complexity is significantly reduced, which is suitable for the mechanical analysis of various composite plates.