Equivalent stiffness calculation of composite hat stiffened laminate

Qiu, Jiabo

doi:10.21595/vp.2022.23045

Vibroengineering Procedia

Browse Procedia

Published: 12 December 2022

Check for updates

Equivalent stiffness calculation of composite hat stiffened laminate

Jiabo Qiu¹

¹Salvage and Rescue Department, Naval Submarine Academic, Qingdao, 266042, China

Cite the article Download PDF

Downloads 389

CrossRef Citations 2

Abstract

The stiffness is the main reason that restricts the large-scale of composite ships, the basis to solve this problem is to solve the stiffness of composite ships’ skeleton structure. Based on classical laminates theory, this paper calculated the equivalent in-plane stiffness and bending stiffness of composite hat stiffened laminate using equivalent stiffness principle. The hat stiffened laminate can be simplified and equivalent to orthotropic laminates using this solution idea, which solves the problem that it is difficult to obtain the analytical solutions of deformation, stress and strain of stiffened laminate. The research results can be used to calculate the stiffness of composite stiffened laminate and the local or global stiffness of the hull girder with such a skeleton structure.

1. Introduction

The single-shell skeleton structure is a typical form of composite hull, which is composed of composite panels and stiffeners. The stiffener is used to strengthen the stiffness and strength of the panels. Common stiffeners include L-shaped, T-shaped, hat shaped, etc. [1]. Although the hap stiffener is complex in technology, its strength and stiffness performance are more excellent, which is a hot research topic in composite field. For example, Yetman et al. [2-3] have studied the debonding behavior between the hat stiffener and the panel, Mo and Ge [4] have conducted experimental research and structural optimization design on the failure of the hat stiffened plate under axial load, Wan et al. [5] have studied on the buckling and post-buckling behavior of hat-stiffened composite panels under axial compressive load.

However, there are few reports on the stiffness of hat stiffened composite laminate. Hu et al. [6] have carried out research on tensile and bending stiffness of hat stiffened laminated beams, but it is only applicable to beam structures and cannot be extended to plate frame structures. This is because it is very difficult to calculate this structure’s stiffness using the classical laminated plate theory [7], and the approximate calculation can only be carried out by combining the equivalent stiffness method. For example, Zhou et al. [8] use the equivalent stiffness method to study the dynamic effects of composite sandwich plates with honeycomb structures, and Sun et al. [9] use the equivalent elastic model to study the elastic properties and lateral stiffness of multi-ribbed walls, and Wang et al. [10] used this method to study the buckling behavior of composite stiffened box structures. These studies all show that the equivalent stiffness is a good method to develop complex composite structures. Qiu and Li [11] calculated the equivalent stiffness of typical elements of composite hat stiffened laminate based on the classical laminated plate theory and the equivalent stiffness method, Wu et al. [12] verified accuracy of this method by experiment and simulation. On this basis, this paper further calculated the equivalent in-plane stiffness and bending stiffness of the entire hat stiffened laminated plate structure.

2. Solution idea of equivalent stiffness on composite hat stiffened laminate

In Fig. 1, the composite hat stiffened laminate with a length of $a_{p}$ , a width of $b_{p}$ , an $x$ -axis along the stiffener direction, and a $y$ -axis perpendicular to the stiffener direction. The cross-section dimension is shown in Fig. 2, the spacing of stiffeners is $d_{s}$ , the width of plate edge is $d_{p}$ , and the stiffeners’ number is $n$ . In the following formula, subscript $e$ represents stiffened laminated element, and $b$ represents edge panel.

The stiffness is a kind of ability to resist external loads. There is a linear relationship between deformation and load when the structure is stressed within the elastic range, as long as the deformation coordination between its constituent elements is maintained, the forces acting on it can be superimposed. The solution idea of this paper is that divide the hat stiffened laminate into $n$ stiffened plate elements and 2 edge plates, as shown in Fig. 1, calculate their stiffness respectively, and then obtain the stiffness of the stiffened laminate by superposition. The stiffness of the hat stiffened laminated element has been solved in reference [11].

Fig. 1Schematic of composite stiffened laminate

Fig. 2Section schematic of composite laminate with hat stiffeners

3. Calculation of equivalent in-plane stiffness of hat stiffened laminate

In the following formulas, the meaning and calculation formula of ${(E T)}_{x e}$ , ${(E T)}_{y e}$ , $μ_{x y e}$ and ${(G T)}_{x y e}$ can refer to reference [11].

3.1. Longitudinal in-plane stiffness

The longitudinal in-plane stiffness ${(E T)}_{x}$ of stiffened laminate is equal to the sum of the longitudinal in-plane stiffness of $n$ typical stiffened elements and 2 edge panels, namely:

1

(E T)_{x} = n (E T)_{x e} + 2 (E T)_{x b} = n (E T)_{x e} + 2 d_{p} / a_{11 p} .

3.2. Transverse in-plane stiffness

When calculating the transverse in-plane stiffness of stiffened laminate, which is decomposed into $n + 2$ segments with $n$ typical stiffened elements in the middle and 2 edge panels on both sides. When bearing the axial force $T_{y}$ per unit length, the longitudinal deformation of the stiffened laminate is the sum of the deformation of each segment, namely:

2

T_{y} b_{p} / (E T)_{y^{'}} = \sum T_{y} b_{i} / (E T)_{y i} \Rightarrow (E T)_{y^{'}} = \frac{b_{p}}{\sum b_{i} / (E T)_{y i}},

where, ${(E T)}_{y^{'}}$ is the transverse in-plane stiffness of the stiffened laminate along the x direction per unit length, $b_{i}$ is the width of each segment, and ${(E T)}_{y i}$ is the transverse in-plane stiffness of each segment. The transverse in-plane stiffness of the edge panels is calculated as follows:

3

(E T)_{y b} = E_{y p} t_{p} = \frac{1}{a_{22 p}} .

Then the transverse in-plane stiffness of the stiffened laminate:

4

(E T)_{y} = a_{p} (E T)_{y^{'}} = a_{p} b_{p} / [\frac{2 d_{p}}{(E T)_{y b}} + \frac{n d_{s}}{(E T)_{y e}}] .

3.3. In-plane Poisson's ratio

When bearing axial force $T_{x}$ , the strain of the stiffened laminate along the x direction is $ε_{x}^{0}$ , and the strain along the y direction is $ε_{y}^{0}$ , then the in-plane Poisson’s ratio $μ_{x y m}$ is defined as $- ε_{y}^{0} / ε_{x}^{0}$ . The deformation of stiffened laminated plate along the $y$ direction is the sum of deformation of each segment, namely:

5

- μ_{x y m} ε_{x}^{0} b_{p} = - \sum μ_{x y i} ε_{x}^{0} b_{i} \Rightarrow μ_{x y m} = \frac{1}{b_{p}} \sum μ_{x y i} b_{i},

where, $μ_{x y i}$ is the in-plane Poisson’s ratio of each segment, $b_{i}$ is the width of each segment. The calculation formula of edge panels is as follows:

6

μ_{x y b} = μ_{x y p} = - \frac{a_{12 p}}{a_{11 p}} .

Then the in-plane Poisson’s ratio of stiffened laminate:

7

μ_{x y m} = \frac{1}{b_{p}} (n d_{s} μ_{x y e} + 2 d_{p} μ_{x y b}) .

3.4. In-plane shear stiffness

Assume that the stiffened laminate bears the shear force along the $y$ direction, the shear strain of the stiffened laminate is $γ_{x y}$ , the shear force of each segment is $T_{x y i}$ ,and the shear strain of each segment is $γ_{x y i}$ , in order to maintain deformation coordination, the shear strain of each segment should be the same as that of the stiffened laminate, namely:

8

γ_{x y} = γ_{x y i} \Rightarrow T_{x y} / (G T)_{x y} = T_{x y i} / (G T)_{x y i},

where, ${(G T)}_{x y}$ is the in-plane shear stiffness of stiffened laminate, and ${(G T)}_{x y i}$ is the in-plane shear stiffness of each segment.

In addition, $T_{x y} = \sum T_{x y i}$ , the stiffness of stiffened laminate can be obtained:

9

(G T)_{x y} = \sum (G T)_{x y i} = n (G T)_{x y e} + 2 (G T)_{x y b} = n (G T)_{x y e} + \frac{2 d_{p}}{a_{66 p}} .

4. Calculation of equivalent bending stiffness of composite stiffened laminate

In the following formulas, the meaning and calculation formula of ${(E I)}_{x e}$ , ${(E I)}_{y e}$ , $μ_{x y e}$ and ${(G J)}_{x y e}$ can refer to reference [11].

4.1. Longitudinal bending stiffness

The longitudinal bending stiffness of stiffened laminate is equal to the sum of that of n stiffened elements and 2 the edge panels, namely:

10

(E I)_{x} = n (E I)_{x e} + 2 (E I)_{x b} = n (E I)_{x e} + 2 d_{p} [1 + 12 {(\frac{t_{p} / 2 - z_{0}}{t_{p}})}^{2}] / d_{11 p} .

4.2. Transverse bending stiffness

When calculating the transverse bending stiffness of stiffened laminate, which is decomposed into $n + 2$ segments. Under the bending moment, the corners $φ$ at both ends of the stiffened laminate should be the sum of the corners $φ_{i}$ of each segment in $y$ direction, namely:

11

φ = \sum φ_{i} \Rightarrow \frac{T_{y} b_{p}}{(E I)_{y^{'}}} = \sum \frac{T_{y} b_{i}}{(E I)_{y i}} \Rightarrow (E I)_{y^{'}} = b_{p} / \sum \frac{b_{i}}{(E I)_{y i}},

where, ${(E I)}_{y^{'}}$ is the transverse bending stiffness per unit width of the stiffened laminate along the $x$ direction, and ${(E I)}_{y i}$ is the transverse bending stiffness per unit width of each segment along the $x$ direction. The calculation formula of the transverse in-plane stiffness of the edge panels is as follows:

12

(E I)_{y b} = E_{y p} t_{p}^{3} / 12 = 1 / d_{22 p} .

Then the transverse bending stiffness of the stiffened laminate:

13

(E I)_{y} = a_{p} (E I)_{y^{'}} = a_{p} b_{p} / [\frac{2 d_{p}}{(E T)_{y b}} + \frac{n d_{s}}{(E T)_{y e}}] .

4.3. Bending Poisson’s ratio

When bearing axial force $T_{x}$ , assume that the curvature of the bending deformation of stiffened laminate in the $x$ direction is $κ_{x}$ . The curvature of bending deformation in $y$ direction is $κ_{y}$ . Then the bending Poisson’s ratio of stiffened laminate $μ_{x y b}$ is defined as $- κ_{y} / κ_{x}$ . The bending angle $φ$ at both ends of the stiffened laminate along the $y$ direction is the sum of the angles $φ_{i}$ of each segment, namely:

14

φ = \sum φ_{i} \Rightarrow - μ_{x y b} κ_{x} b_{p} = - \sum μ_{x y i} κ_{x} b_{i} \Rightarrow μ_{x y b} = \frac{1}{b_{p}} \sum μ_{x y i} b_{i},

where, $μ_{x y i}$ is the bending Poisson’s ratio of each segment and $b_{i}$ is the width of each segment. The bending Poisson’s ratio formula of the edge panels is as follows:

15

μ_{x y b} = μ_{x y p} = - d_{12 p} / d_{11 p} .

Then the bending Poisson’s ratio of stiffened laminate:

16

μ_{x y b} = \frac{1}{b_{p}} (n d_{s} μ_{x y e} + 2 d_{p} μ_{x y b}) .

4.4. Torsional stiffness

When bearing torque $T_{x y}$ , the torsion ratio of stiffened laminate is $φ^{'}$ , the torque of each segment is $T_{x y i}$ , and the torsion ratio is $φ_{i}^{'}$ . The torque of the stiffened laminate is the sum of the torque of each segment, and the torque ratio is equal to that of each segment, namely:

17

φ^{'} = φ_{i}^{'} \Rightarrow T_{x y} / (G J)_{x y} = T_{x y i} / (G J)_{x y i} \Rightarrow \sum T_{x y i} / (G J)_{x y} = T_{x y i} / (G J)_{x y i} .

Simplify the above formula to obtain the torsional stiffness of stiffened laminate:

18

(G J)_{x y} = \sum (G J)_{x y i} = n (G J)_{x y e} + 2 (G J)_{x y b} = n (G J)_{x y e} + 8 d_{p} / d_{66 p} .

5. Calculation of equivalent stiffness coefficient

After the calculation of equivalent stiffness, composite stiffened laminate can be equivalent to orthotropic composite laminate. The stiffness coefficient can be calculated according to the relationship between the mid-plane strain, curvature, internal force and internal moment of the composite laminate.

5.1. Calculation of in-plane stiffness coefficient

The relationship between in-plane strain and in-plane internal force of orthotropic composite laminate:

19

\{\begin{array}{l} ε_{x}^{0} \\ ε_{y}^{0} \\ γ_{x y}^{0} \end{array}\} = [\begin{matrix} a_{11} & a_{12} & 0 \\ a_{12} & a_{22} & 0 \\ 0 & 0 & a_{66} \end{matrix}] \{\begin{matrix} N_{x} \\ N_{y} \\ N_{x y} \end{matrix}\},

where, $ε_{x, y}^{0}$ is the in-plane strain of the laminate, $N_{x, y}$ is the in-plane internal force, $a_{i j}$ is the in-plane flexibility coefficient.

When the composite stiffened laminate bears in-plane internal force $N_{x}$ , the following formula can be obtained according to the Eq. (19):

20

ε_{x}^{0} = a_{11} N_{x} \Rightarrow P_{x} / (E T)_{x} = a_{11} N_{x} \Rightarrow N_{x} b_{p} / (E T)_{x} = a_{11} N_{x} ， ε_{y}^{0} = a_{12} N_{x} .

Simplify the above equation to obtain:

21

a_{11} = b_{p} / (E T)_{x}, μ_{x y m} = - ε_{y}^{0} / ε_{x}^{0} = - a_{12} / a_{11} .

In the same way, we can get:

22

a_{22} = a_{p} / (E T)_{y}, a_{66} = b_{p} / (G T)_{x y} .

In-plane stiffness matrix of composite stiffened laminate is the transpose of in-plane flexibility matrix.

5.2. Calculation of bending stiffness coefficient

The relationship between curvature and internal moment of orthotropic laminate:

23

\{\begin{matrix} κ_{x} \\ κ_{y} \\ κ_{x y} \end{matrix}\} = [\begin{matrix} d_{11} & d_{12} & 0 \\ d_{12} & d_{22} & 0 \\ 0 & 0 & d_{66} \end{matrix}] \{\begin{matrix} M_{x} \\ M_{y} \\ M_{x y} \end{matrix}\},

where, $κ_{x, y}$ is the curvature of the composite laminate, $M_{x, y}$ is the internal moment, $d_{i j}$ is the bending flexibility coefficient.

In the same way as section 5.1, we can get:

24

d_{11} = \frac{b_{p}}{(E I)_{x}}, \frac{d_{12}}{d_{11}} = \frac{κ_{y}}{κ_{x}} = - μ_{x y b}, d_{22} = \frac{a_{p}}{(E I)_{y}}, d_{66} = \frac{2 b_{p}}{(G J)_{x y}} .

Bending stiffness matrix of composite stiffened laminate is the transpose of bending flexibility matrix.

6. Conclusions

The research on the stiffness characteristics of composite stiffened laminate is the basis of that of composite hull. Based on the classical laminated plate theory, the equivalent in-plane stiffness and bending stiffness of composite hat stiffened laminate are obtained using the stiffness equivalence method.

Using the solution idea in this paper, the arbitrary composite stiffened laminate can be transformed into orthotropic laminate. Under the condition of the specific boundary conditions and loads, the deformation, internal force and bending moment of orthotropic laminated plates can be solved first, and then the deformation, stress and strain of each element of stiffened laminate can be further solved according to the classical laminated plate theory. The equivalent stiffness of stiffened laminate calculated by this method refers to the average distributed stiffness. Therefore, the more intensive the number of stiffeners is, the more accurate the calculation results will be.

The research results in this paper can be used to calculate the stiffness of composite stiffened laminate, and the local or global stiffness of the hull girder with such a skeleton structure.

References

M. F. Cao and L. Yang, “Application of composite materials in ship design and construction,” Naval Architecture and Ocean Engineering, Vol. 2, pp. 38–43, 2006, https://doi.org/10.3969/j.issn.1005-9962.2006.02.014

Publisher
J. E. Yetman and A. J. Sobey, “Investigation into skin stiffener debonding of top-hat stiffened composite structures,” Composite structure, Vol. 132, pp. 1168–1181, 2015, https://doi.org/10.1016/j.compstruct.2015.06

Publisher
J. E. Yetman, A. J. Sobey, J. I. R. Blake, and R. A. Shenoi, “Modelling the variability of skin stiffener debonding in post-cured top-hat stiffened panels,” Composite Structures, Vol. 211, pp. 187–195, Mar. 2019, https://doi.org/10.1016/j.compstruct.2018.11.078

Publisher
Y. Mo, D. Ge, and B. He, “Experiment and optimization of the hat-stringer-stiffened composite panels under axial compression,” Composites Part B: Engineering, Vol. 84, pp. 285–293, Jan. 2016, https://doi.org/10.1016/j.compositesb.2015.08.039

Publisher
H. L. Wan, W. M. Gong, and J. E. Qin, “Research on the buckling and post-buckling behavior of hat-stiffened composite panels under axial compressive load,” Structure and Environment Engineering, Vol. 46, No. 1, pp. 31–34, 2019, https://doi.org/10.19447/j.cnki.11-1773/v.2019.01.005

Publisher
M. Y. Hu, B. Zheng, and A. W. Wang, “Theoretical calculation of stiffness of laminated composite beam with hat-section,” Journal of Naval University of Engineering, Vol. 29, No. 6, pp. 6–12, 2017, https://doi.org/10.7495/j.issn.1009-3486.2017.06.002

Publisher
G. Q. Jiao and P. R. Jia, Mechanics of Composite Materials. Xi’an: Northwest Polytechnical University Press, 2008.

Search CrossRef
X. Q. Zhou, L. Wang, D. Y. Yu, and C. Y. Zhang, “Dynamic effective equivalent stiffness analysis on the periodical honeycomb reinforced composite laminated structure filled with viscoelastic damping material,” Composite Structures, Vol. 193, pp. 306–320, Jun. 2018, https://doi.org/10.1016/j.compstruct.2018.03.066

Publisher
J. Sun, Y. Jia, and Y. L. Mo, “Evaluation of elastic properties and lateral stiffness of multi-ribbed walls based on equivalent elastic model,” Engineering Structures, Vol. 72, pp. 92–101, Aug. 2014, https://doi.org/10.1016/j.engstruct.2014.04.014

Publisher
P. Wang, X. Huang, Z. Wang, X. Geng, and Y. Wang, “Buckling and post-buckling behaviors of a variable stiffness composite laminated wing box structure,” Applied Composite Materials, Vol. 25, No. 2, pp. 449–467, Apr. 2018, https://doi.org/10.1007/s10443-017-9643-3

Publisher
J. B. Qiu and H. D. Li, “Equivalent stiffness calculation of typical plate element of composite cap stiffened laminates,” Composites Science and Engineering, Vol. 7, pp. 33–39, 2020, https://doi.org/10.3969/j.issn.1003-0999.2020.07.005

Publisher
J. Wu et al., “Equivalent bending stiffness of composite hat-stiffened panel,” Acta Materiae Compositae Sinica, Vol. 29, No. 12, pp. 1–8, 2022, https://doi.org/10.13801/j.cnki.fhclxb.20211208.001

Publisher

Cited by

Experimental investigation on the residual bending capacity of carbon fiber corrugated sandwich structures under the repeated impact

(2024)

An improved model of fatigue life for stitched carbon fiber composites considering stiffness degradation

Xintian Liu | Xiaoying Zhang | Shuci Wang

(2023)

About this article

Received

08 November 2022

Accepted

28 November 2022

Published

12 December 2022

SUBJECTS

Materials and measurements in engineering

DOI

https://doi.org/10.21595/vp.2022.23045

Keywords

composite

laminate

hat stiffened

equivalent stiffness

Acknowledgements

The project is supported by Youth Program of National Natural Science Foundation of China: 51609252.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous article in issue Previous Next article in issue Next

Research article

2021 12 02

Hydromount elastic shell dynamic stiffness calculation using finite difference method

Andrey Daryenkov, Boris Gordeev, Sergey Okhulkov, Artem Ermolaev, Evgeniy Chernov

Research article

2018 12 31

Specific damping capacity calculation of composite plates with delamination based on higher-order Zig-Zag theory

Chaogan Gao, Chuwei Zhou

Research article

2016 12 08

Modeling and analysis of performance degradation of bolted composite lap geometries subjected to shock load

Mariyam Sattar, Cheng Wei, Awais Jalali

Research article

2016 02 15

Equivalent stiffness and dynamic response of new mechanical elastic wheel

Qiang Wang, Youqun Zhao, Xianbin Du, Mingmin Zhu, Hongxun Fu

J. Qiu, “Equivalent stiffness calculation of composite hat stiffened laminate,” Vibroengineering PROCEDIA, Vol. 47, pp. 42–48, Dec. 2022, https://doi.org/10.21595/vp.2022.23045

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/vp.2022.23045
UR  - https://doi.org/10.21595/vp.2022.23045
TI  - Equivalent stiffness calculation of composite hat stiffened laminate
T2  - Vibroengineering PROCEDIA
AU  - Qiu, Jiabo
PY  - 2022
DA  - 2022/12/12
PB  - JVE International Ltd.
SP  - 42-48
VL  - 47
SN  - 2345-0533
SN  - 2538-8479
ER  - 

Copy Ris

Copied to clipboard!

@article{Qiu_2022,
	doi = {10.21595/vp.2022.23045},
	url = {https://doi.org/10.21595/vp.2022.23045},
	year = 2022,
	month = {dec},
	publisher = {{JVE} International Ltd.},
	volume = {47},
	pages = {42--48},
	author = {Jiabo Qiu},
	title = {Equivalent stiffness calculation of composite hat stiffened laminate},
	journal = {Vibroengineering {PROCEDIA}}
}

Copy Bibtex

Copied to clipboard!

[1]J. Qiu, “Equivalent stiffness calculation of composite hat stiffened laminate,” Vibroengineering PROCEDIA, vol. 47, pp. 42–48, Dec. 2022, doi: 10.21595/vp.2022.23045.

Copy IEEE

Copied to clipboard!

Qiu, Jiabo. “Equivalent Stiffness Calculation of Composite Hat Stiffened Laminate.” Vibroengineering PROCEDIA 47 (December 12, 2022): 42–48. https://doi.org/10.21595/vp.2022.23045.

Copy Chicago

Copied to clipboard!

Equivalent stiffness calculation of composite hat stiffened laminate

Abstract

1. Introduction

2. Solution idea of equivalent stiffness on composite hat stiffened laminate

3. Calculation of equivalent in-plane stiffness of hat stiffened laminate

3.1. Longitudinal in-plane stiffness

3.2. Transverse in-plane stiffness

3.3. In-plane Poisson's ratio

3.4. In-plane shear stiffness

4. Calculation of equivalent bending stiffness of composite stiffened laminate

4.1. Longitudinal bending stiffness

4.2. Transverse bending stiffness

4.3. Bending Poisson’s ratio

4.4. Torsional stiffness

5. Calculation of equivalent stiffness coefficient

5.1. Calculation of in-plane stiffness coefficient

5.2. Calculation of bending stiffness coefficient

6. Conclusions

References

Cited by

About this article

Related Articles