Next Article in Journal
Selective Etching of Si versus Si1−xGex in Tetramethyl Ammonium Hydroxide Solutions with Surfactant
Next Article in Special Issue
Load Eccentricity of Compressed Composite Z-Columns in Non-Linear State
Previous Article in Journal
Relation between Ga Vacancies, Photoluminescence, and Growth Conditions of MOVPE-Prepared GaN Layers
Previous Article in Special Issue
Elastic Properties of Open Cell Metallic Foams—Modeling of Pore Size Variation Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Stability of Natural-Fiber-Reinforced Composite Cylindrical Shells with Initial Geometric Imperfection Considering Moisture Absorption and Hygrothermal Aging

1
School of Mechanical Engineering, Dalian Jiaotong University, Dalian 116028, China
2
Dalian Rail Transit Design Institute Co., Ltd., Dalian 116021, China
3
School of Civil Engineering, Dalian Jiaotong University, Dalian 116028, China
4
School of Traffic and Transportation Engineering, Dalian Jiaotong University, Dalian 116028, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(19), 6917; https://doi.org/10.3390/ma15196917
Submission received: 3 September 2022 / Revised: 28 September 2022 / Accepted: 30 September 2022 / Published: 5 October 2022
(This article belongs to the Special Issue Compressive Behavior of Materials and Structures)

Abstract

:
In this paper, the nonlinear stability of a natural-fiber-reinforced composite cylindrical shell with initial geometric imperfection is investigated. The nonlinear governing equations are established by high-order shear deformation theory. The load-edge shortening curves for different imperfection amplitudes are obtained by the Galerkin method. Several numerical examples are presented to verify the accuracy of the proposed method and to investigate the influence of initial geometric imperfection, moisture absorption, and hygrothermal aging on the post-buckling behavior of natural-fiber-reinforced composite cylindrical shells.

1. Introduction

Natural-fiber-reinforced composite (NFRC) is a kind of composite made of bio-based fibers and polymer matrix [1,2]. Compared with traditional glass/carbon fiber-reinforced composites, NFRC has many advantages, such as low cost, light weight, superior degradability, renewable, and so on [3]. NFRC is more environment friendly and has been widely used in railway coaches, automobiles, aerospace, buildings, construction, and protective armors [4,5]. However, NFRC still has the disadvantage of the significant degradation of mechanical properties caused by moisture absorption and hygrothermal aging [6,7]. This is because the natural fiber has a porous structure capable of absorbing and storing large amounts of water. Water molecules occupy the free spaces as free water and they also make connections with polymer chains as bound water. The main process of moisture absorption is the diffusion of water molecules into the micropores between polymer chains of natural fibers; the natural fibers swell dramatically and the mechanical properties of the natural fibers are significantly degraded. The second process is the capillary transport into the gaps and flaws generated at the interfaces between the fibers and the polymer matrix due to incomplete wettability and impregnation [8,9]. Therefore, it is of great significant to consider moisture absorption and hygrothermal aging in the investigation of mechanical behaviors of NFRC structures.
Some researchers have studied the mechanical properties of natural-fiber-reinforced composites. Khakpour et al. [10] investigated the effects of fiber treatment, fiber length, and fiber density on the mechanical properties of structural adhesives enhanced with natural date palm tree fibers. Akhavan-Safar et al. [11] investigated the effect of natural date palm tree fibers (DPTF) on the R-curve of E-glass/epoxy end-notch-flexure samples. Dalla Libera et al. [12] manufactured prepregs by impregnating discontinuous curaua fibers with B-stage epoxy resin and investigated the thermal, dynamic mechanical, mechanical, and morphological behavior. Ramaswamy et al. [13] provided the analytical formulation to calculate the free-edge stress values induced in natural fiber composite materials. Pawar et al. [14] studied the dynamic and mechanical characteristics of banana-epoxy, kenaf-epoxy, and hemp-epoxy natural-fiber-reinforced composite materials. Nandhakumar et al. [15] performed the experimental investigation of the mechanical properties (rigidity and hardness quality) of natural-fiber-reinforced composites. Patel and Jain [16] studied the mechanical properties in randomly oriented short natural fiber composites and discussed the different combinations of short natural fibers with reinforced plastics. Gupta et al. [17] investigated the effects of fiber orientation on the mechanical properties of natural-fiber-reinforced epoxy composite by a finite element method. The mechanical behaviors of NFRC beams, plates, and shells have also been investigated. Lim et al. [18] performed the direct measurement of the whole-field deformation of an NFRC cantilever beam by the shadow moiré method and the comparison of the experimental results with the finite element results. Saravana Bavan and Mohan Kumar [19] analyzed the deflection and stress properties of an NFRC beam by the finite element method (FEM). Jirawattanasomkul et al. [20] investigated the shear behavior of pre-damaged deep reinforced concrete beams strengthened with NFRC sheets by the non-reversed cyclic three-point bending load test. Rajeshkumar and Hariharan [21] investigated the free vibration of NFRC beams for different fiber gauge lengths under clamped-free boundary condition. Rajesh et al. [22] investigated the influence of surface pre-treatment with sodium hydroxide and hybridization effect on the free vibration of NFRC beams. Senthil Kumar et al. [23] conducted an experimental investigation on the free vibration of NFPC beams and studied the influence of fiber length and weight percentage on the mechanical properties and free vibration characteristics. Kuriakose et al. [24] investigated the free vibration of NFRC plates and analyzed the influence of carbon nanotube fillers on the vibration behavior. Selvaraj et al. [25] investigated the material, mechanical, and dynamic characteristics of natural-fiber-reinforced composite sandwich plates with multiple-core layers. Wang and Wang [26] proposed a mechanical model for NFRC under hydrothermal loading and obtained the closed-form solutions for the buckling load and natural frequency of simply supported NFRC plates. Xu et al. [27] performed the free vibration and buckling analysis of a partially or internally cracked NFRC plate with corner-point supports by using the symplectic method. Zhang et al. [28] investigated the buckling of NFRC cylindrical shells based on Reissner’s shell theory and the symplectic method.
From the above research, it can be found that there have been few investigations into NFRC structures, especially the NFRC cylindrical shell. As an important basic structure in engineering, the cylindrical shell has the advantages of light weight, strong bearing capacity, and simple processing, and the main failure form is buckling instability [29]. The existing research on the stability of NFRC cylindrical shell [28] is based on the eigenvalue buckling of a perfect shell, and experimental research [30] found that buckling instability will occur under a load lower than the theoretical critical load. Koiter [31] pointed out that this is because the actual structure inevitably has initial geometric imperfection that will significantly reduce the stability. The stability of cylindrical shells considering the initial geometric imperfection has been investigated by some researchers. Narayana et al. [32] investigated the post-buckling and imperfection sensitivity of isotropic and composite cylindrical shells by FEM. Sun et al. [33] proposed the accelerated Koiter method to analyze the post-buckling and imperfection sensitivity of isotropic, composite and stiffened cylindrical shells. Shariyat and Eslami [34] and Eslami et al. [35] investigated the dynamic post-buckling of imperfect orthotropic and laminated cylindrical shells by using the finite difference method and Runge–Kutta methods. Shen [36,37,38,39,40] proposed the boundary layer theory and analyzed the post-buckling behavior of imperfect composite cylindrical shells by using the singular perturbation technique. Using a similar method, Babaei et al. [41] investigated the post-buckling of imperfect functionally graded material (FGM) cylindrical shells. Li et al. [42] analyzed the post-buckling of imperfect anisotropic laminated cylindrical shell. Foroutan et al. [43] analyzed the static and dynamic post-buckling of imperfect stiffened FGM cylindrical shells by the Galerkin method and Runge–Kutta methods. It can be seen that the research into the stability of imperfect cylindrical shells mainly depends on numerical methods; there is little analytical research.
Motivated by these reasons, we proposed an analytical method for the post-buckling of imperfect NFRC cylindrical shells with moisture absorption and hygrothermal aging. Firstly, the nonlinear governing equations for the post-buckling of imperfect NFRC cylindrical shells are established based on higher-order shear deformation theory. Secondly, the displacement functions and imperfection function are proposed and the load-edge shortening curves are obtained by the Galerkin method. Finally, the accuracy is verified and the influences of initial geometric imperfection, moisture absorption, and hygrothermal aging are discussed.
This paper is organized as follows: Section 2 presents the model of imperfect NFRC cylindrical shells. Section 3 presents the basic equations for post-buckling problems. Section 4 presents the procedure to solve the post-buckling equilibrium path. Section 5 presents the numerical examples and discusses the effects of calculation parameters. The conclusions are drawn in Section 6.

2. Model of Imperfect NFRC Cylindrical Shells

Consider an imperfect NFRC cylindrical shell under axial compression F in Figure 1. The composite is unidirectional long sisal fiber-reinforced benzylated wood. The layer number of the NFRC cylindrical shell is 8 and the layer angle is (0/90)2S. The geometric parameters of the shell are taken as length L, radius of the mid-surface R, and thickness h. It refers to an orthogonal curvilinear coordinate system in which the x−, y−, and z− axes represent the longitudinal, circumferential, and normal of the mid-surface of the shell, respectively. The displacements of the mid-surface of the shell along the x−, y−, and z− axes are u, v, and w, respectively. Φx and Φy are the rotations of the transverse normal.
The effect of the material properties of natural-fiber-reinforced composite considering the influence of moisture absorption and hygrothermal aging can be expressed as [26,27]:
E 1 = 3 G ( 1 v f 2 β s γ s v f 3 β l γ l ) + 4 v f η G ( 1 χ s γ s ) ( c 4 + s 4 c 2 s 2 ) 1 + v f η [ ( 1 χ s γ s ) / ( 1 v f 2 β s γ s v f 3 β l γ l ) ] s 4
E 2 = 3 G ( 1 v f 2 β s γ s v f 3 β l γ l ) + 4 v f η G ( 1 χ s γ s ) ( c 4 + s 4 c 2 s 2 ) 1 + v f η [ ( 1 χ s γ s ) / ( 1 v f 2 β s γ s v f 3 β l γ l ) ] c 4
ν 12 = ( 1 v f 2 β s γ s v f 3 β l γ l ) + 2 v f η ( 1 χ s γ s ) c 2 s 2 2 ( 1 v f 2 β s γ s v f 3 β l γ l ) + 2 v f η ( 1 χ s γ s ) s 4
ν 21 = ( 1 v f 2 β s γ s v f 3 β l γ l ) + 2 v f η ( 1 χ s γ s ) c 2 s 2 2 ( 1 v f 2 β s γ s v f 3 β l γ l ) + 2 v f η ( 1 χ s γ s ) c 4
G 12 = G ( 1 v f 2 β s γ s v f 3 β l γ l ) + 4 v f η G ( 1 χ s γ s ) c 2 s 2
G 13 = G 23 = G ( 1 v f 2 β s γ s v f 3 β l γ l )
where E1, E2, G12, G13, G23, and ν12, ν21 are the Young’s moduli, shear moduli, and Poisson’s ratio, respectively, c = cosθ, s = sinθ and θ is the layer angle, γs and γl are two important coefficients describing the damage evolution induced by moisture absorption and hygrothermal aging and their expressions are γ s = 1 e η v f G t / D and γ l = 1 e v f T t v f T t e v f T t , vf is the fiber content, η is a dimensionless constant which can be determined by fitting the experiential date, D is a material parameter which governing the speed of energy dissipation, G is the initial shear modulus of the composites, t is the aging time (in h for hours), T is a material index which can be determined by fitting the empirical data, and βs, βl, and χs are three constants that can be obtained by three times multinomial fittings of the experimental results. These parameters are related to the experimental conditions (temperature, humidity, etc.). In this paper, we aim to propose an analytical method to analyze the post-buckling of imperfect NFRC cylindrical shells, so we select only the values available in the open literature [26] under specific humidity and temperature; the parameters are G = 0.104 GPa, η = 90.870, βs = 7.3459, βl = 200.0104, χs = 0.6171, D = 56.18 hN/m2, and T = 0.1/h.

3. Basic Equations for Post-Buckling of Imperfect NFRC Cylindrical Shells

According to the high-order shear deformation theory [44], the displacement components at a point of the shell can be expressed as:
u 1 ( x ,   y ,   z ) = u ( x ,   y ) 4 z 3 3 h 2 w ( x ,   y ) x + ( z 4 z 3 3 h 2 ) Φ x ( x ,   y )
u 2 ( x ,   y ,   z ) = ( 1 + z R ) v ( x ,   y ) 4 z 3 3 h 2 w ( x ,   y ) y + ( z 4 z 3 3 h 2 ) Φ y ( x ,   y )
u 3 ( x ,   y ,   z ) = w ( x ,   y ) + w 0 ( x ,   y )
Based on Novozhilov’s nonlinear shell theory [45], the nonlinear strain–displacement relations are:
ε x = u 1 x + 1 2 [ ( u 1 x ) 2 + ( u 2 x ) 2 + ( u 3 x ) 2 ]
ε y = u 2 y + u 3 R 1 + z R + ( u 1 y ) 2 + ( u 2 y + u 3 R ) 2 + ( u 2 R u 3 y ) 2 2 ( 1 + z R ) 2
γ x y = u 2 x + u 1 y 1 + z R + u 1 x u 1 y + u 2 x ( u 2 y + u 3 R ) + u 3 x ( u 2 R u 3 y )
γ x z = u 1 z + u 3 x + u 1 x u 1 z + u 2 x u 2 z + u 3 x u 3 z
γ y z = u 2 z + u 3 y u 2 R 1 + z R + u 1 y u 1 z + ( u 2 y + u 3 R ) u 2 z + ( u 3 y u 2 R ) u 3 z 1 + z R
Substituting Equations (7)–(9) into Equations (10)–(14) and retaining only the nonlinear terms with w and w0, we have:
ε = ε ( 0 ) + z χ ( 0 ) + z 3 χ ( 2 )
γ = γ ( 0 ) + z 2 χ ( 1 )
where ε = { ε x , ε y , γ x y } T , γ = { γ x z , γ y z } T , and:
ε ( 0 ) = { ε x ( 0 ) ε y ( 0 ) γ x y ( 0 ) } = { u x + 1 2 ( w x ) 2 + w x w 0 x v y + w R + 1 2 ( w y ) 2 + w y w 0 y u y + v x + w x w y + w x w 0 y + w 0 x w y }
χ ( 0 ) = { χ x ( 0 ) χ y ( 0 ) χ x y ( 0 ) } = { Φ x x Φ y y Φ x y + Φ y x }
χ ( 2 ) = { χ x ( 2 ) χ y ( 2 ) χ x y ( 2 ) } = 4 3 h 2 { Φ x x + 2 w x 2 Φ y y + 2 w y 2 Φ x y + Φ y x + 2 w x y }
γ ( 0 ) = { γ x z ( 0 ) γ y z ( 0 ) } = { Φ x + w x Φ y + w y }
χ ( 1 ) = { χ x z ( 1 ) χ y z ( 1 ) } = 4 h 2 { Φ x + w x Φ y + w y }
The constitutive equation in the k-th lamina of the imperfect NFRC cylindrical shell can be written as:
{ σ x k σ y k τ x y k τ x z k τ y z k } = [ C 11 k C 12 k 0 0 0 C 21 k C 22 k 0 0 0 0 0 G 12 k 0 0 0 0 0 G 13 k 0 0 0 0 0 G 23 k ] { ε x k ε y k γ x y k γ x z k γ y z k }
where σ x k , σ y k , τ x y k , τ x z k , and τ y z k are the in-plane stress and transverse shear stress. C 11 k = E 1 k / ( 1 ν 12 k ν 21 k ) , C 12 k = C 21 k = E 1 k ν 21 k / ( 1 ν 12 k ν 21 k ) , and C 22 k = E 2 k / ( 1 ν 12 k ν 21 k ) .
Integrating Equation (14) across the thickness, the resultant force vector N, moment vector M, high-order force vector P, shear force vector Q, and high-order shear force vector R can be obtained as:
N = { N x ,   N y ,   N x y } T = i = 1 K h i h i + 1 { σ x k ,   σ y k ,   τ x y k } T   d z
M = { M x ,   M y ,   M x y } T = i = 1 K h i h i + 1 { σ x k ,   σ y k ,   τ x y k } T z d z
P = { P x ,   P y ,   P x y } T = i = 1 K h i h i + 1 { σ x k ,   σ y k ,   τ x y k } T z 3 d z
Q = { Q x ,   Q y } T = i = 1 K h i h i + 1 { τ x z k ,   τ y z k } T d z
R = { R x ,   R y } T = i = 1 K h i h i + 1 { τ x z k ,   τ y z k } T z 2 d z
where K is the layer number.
The equilibrium equations of imperfect NFRC cylindrical shells can be expressed as:
N x x + N x y y = 0
N x y x + N y y = 0
Q x x + Q y y 4 h 2 ( R x x + R y y ) + 4 3 h 2 ( 2 P x x 2 + 2 2 P x y x y + 2 P y y 2 ) N y R + N x 2 ( w + w 0 ) x 2 + 2 N x y 2 ( w + w 0 ) x y + N y 2 ( w + w 0 ) y 2 = 0
M x x + M x y y 4 3 h 2 ( P x x + P x y y ) Q x + 4 h 2 R x = 0
M x y x + M y y 4 3 h 2 ( P x y x + P y y ) Q y + 4 h 2 R y = 0
where w0 is the initial geometric imperfection.
Airy’s stress function ψ is introduced to simplify the equilibrium equations, which satisfy:
N x = 2 ψ y 2 ,   N y = 2 ψ x 2 ,   N x y = 2 ψ x y
The first two equilibrium equations, Equations (20) and (21), are transformed into a compatible equation as follows:
2 ε x ( 0 ) y 2 + 2 ε y ( 0 ) x 2 2 γ x y ( 0 ) x y = 1 R 2 w x 2 2 w x 2 2 w y 2 + ( 2 w x y ) 2 2 w x 2 2 w 0 y 2 + 2 2 w x y 2 w 0 x y 2 w y 2 2 w 0 x 2
The equilibrium equations about stress function ψ, deflection w, and rotations of the transverse normal Φx and Φy can be rewritten as:
δ 1 4 ψ x 4 + δ 2 4 ψ x 2 y 2 + δ 3 4 ψ y 4 + δ 4 4 w x 4 + δ 5 4 w x 2 y 2 + δ 6 4 w y 4 + δ 7 3 Φ x x 3 + δ 8 3 Φ x x y 2 + δ 9 3 Φ y x 2 y + δ 10 3 Φ y y 3 1 R 2 w x 2 + 2 w x 2 2 w y 2 ( 2 w x y ) 2 + 2 w x 2 2 w 0 y 2 2 2 w x y 2 w 0 x y + 2 w y 2 2 w 0 x 2 = 0
δ 4 4 ψ x 4 + δ 5 4 ψ x 2 y 2 + δ 6 4 ψ y 4 + δ 11 4 w x 4 + δ 12 4 w x 2 y 2 + δ 13 4 w y 4 + δ 14 2 w x 2 + δ 14 2 w y 2 + δ 14 Φ x x + δ 14 Φ y y + δ 15 3 Φ x x 3 + δ 16 3 Φ x x y 2 + δ 17 3 Φ y x 2 y + δ 18 3 Φ y y 3 1 R 2 ψ x 2 + 2 ψ y 2 2 w x 2 2 2 ψ x y 2 w x y + 2 ψ x 2 2 w y 2 + 2 ψ y 2 2 w 0 x 2 2 2 ψ x y 2 w 0 x y + 2 ψ x 2 2 w 0 y 2 = 0
δ 7 3 ψ x 3 + δ 8 3 ψ x y 2 + δ 14 w x + δ 15 3 w x 3 + δ 16 3 w x y 2 + δ 14 Φ x + δ 19 2 Φ x x 2 + δ 20 2 Φ x y 2 + δ 21 2 Φ y x y = 0
δ 9 3 ψ x 2 y + δ 10 3 ψ y 3 + δ 14 w y + δ 16 3 w x 2 y + δ 17 3 w y 3 + δ 14 Φ y + δ 20 2 Φ y x 2 + δ 22 2 Φ y y 2 + δ 21 2 Φ x x y = 0
The boundary conditions at the two ends of imperfect NFRC cylindrical shells are constrained and satisfy:
w = 0 ,   w x = 0 , Φ x = 0 ,   Φ y = 0 ,   0 2 π R N x d y = F

4. Procedure for Solving the Post-Buckling Equilibrium Path

The post-buckling equilibrium path could be solved by the Galerkin method, for which the trial functions needs to be assumed. The deformation for the post-buckling is very complex and the displacement function should have more expansion terms. The displacement functions that satisfy boundary condition (39) are assumed as [46,47]:
w = n m a m n { cos [ ( m 1 ) π x L ] cos [ ( m + 1 ) π x L ] } cos ( n y N R )
Φ x = n m b m n { sin [ ( m 1 ) π x L ] sin [ ( m + 1 ) π x L ] } cos ( n y N R )
Φ y = n m c m n { cos [ ( m 1 ) π x L ] cos [ ( m + 1 ) π x L ] } sin ( n y N R )
where N is the circumferential wave number, amn, bmn, and cmn are the undetermined coefficients; m and n are selected as m = 1, 3, 5… and n = 0, 1, 2.
The initial geometric imperfection function is assumed as [46]:
w 0 = 0.5 μ h [ 1 cos ( 2 π x L ) ] cos ( y N R )
where μ is the imperfection amplitude.
Substituting Equations (40)–(43) into Equation (35), we have:
δ 1 4 ψ x 4 + δ 2 4 ψ x 2 y 2 + δ 3 4 ψ y 4 = p = 0 q = 0 ψ p q cos ( p π x L )   cos ( q y N R )  
The stress function can be expressed as a similar form to the right side of Equation (44) as follows:
ψ = p x y 2 + p = 0 q = 0 Ψ p q cos ( p π x L )   cos ( q y N R )  
where px = −(F/2π2R2) is used to satisfy the load boundary condition.
To obtain the undetermined coefficients, the Galerkin method is applied to the equilibrium Equations (36)–(38), and the following relations can be obtained:
0 2 π R 0 L L 1 ( w ,   Φ x ,   Φ y ) { cos [ ( m 1 ) π x L ] cos [ ( m + 1 ) π x L ] } cos ( n y N R ) d x d y = 0
0 2 π R 0 L L 2 ( w ,   Φ x ,   Φ y ) { sin [ ( m 1 ) π x L ] sin [ ( m + 1 ) π x L ] } cos ( n y N R ) d x d y = 0
0 2 π R 0 L L 3 ( w ,   Φ x ,   Φ y ) { cos [ ( m 1 ) π x L ] cos [ ( m + 1 ) π x L ] } sin ( n y N R ) d x d y = 0
By solving the nonlinear equations presented in Equations (46)–(48), the undetermined coefficients amn, bmn, and cmn can be determined and the end shortening can be calculated from the following equation:
= 0 L u x   d x

5. Numerical Results and Discussion

5.1. Comparison Study

There is no available research on the stability of imperfect NFRC cylindrical shells in the open literature. Firstly, the present results are compared with those of imperfect isotropic cylindrical shells to verify the accuracy of the proposed method. The computation parameters are E = 5.56 GPa, ν = 0.3, L = 71.96 mm, h = 0.247 mm, and R = 100 mm. The dimensionless axial compression is Σ = F/Fcl where Fcl = 2πEh2/[3(1 − ν2)]0.5. The dimensionless edge shortening is δ = (ΔR)/(Lh). The load-edge shortening curves for different imperfection amplitudes are shown in Figure 2. It is obvious that the present results are in good agreement with the reference results [46].
Subsequently, the present results are compared with the experimental results [48] for the buckling of imperfect carbon-fiber-reinforced composite cylindrical shells. The geometric parameters are L = 540 mm, h = 1.2 mm, R = 350 mm. The layer number is 8 and the laying angle is (0/90)2S. The load-edge shortening curves are shown in Figure 3. It is observed that the present results are in good agreement with the experimental results, and the buckling mode is very close.
Finally, the present results are compared with the finite element (FE) results by considering single-layer imperfect NFRC cylindrical shells. The geometric parameters are L = 1000 mm, h = 10 mm, and R = 1000 mm. The laying angle is 0°, the fiber content is vf = 0.1, and the aging time is t = 50 h. The load-edge shortening curves for different imperfection amplitudes are presented in Figure 4. The good agreement with the FE results can also be observed. From the above comparison, it is deduced that the proposed method can accurately analyze the stability of imperfect NFRC cylindrical shells.

5.2. Parameter Study

In this section, the effects of initial geometric imperfection, moisture absorption and hygrothermal aging on the stability of NFRC cylindrical shells are investigated. The geometrical parameters are L = 100 mm, h = 1 mm, and R = 100 mm. The dimensionless axial compression is Σ = F/Fcl, where Fcl = 2πEmh2/[3(1 − νm2)]0.5, Em = 0.312 GPa, and νm = 0.5. The dimensionless edge shortening is δ = (ΔR)/(Lh).
The load-edge shortening curves of NFRC cylindrical shells for different imperfection amplitudes are presented in Figure 5. The other parameters are t = 50 h and vf = 0.1. It is observed that the extreme point decreases with the increase in imperfection amplitude, which indicates that the NFRC cylindrical shell is sensitive to the initial geometric imperfection. When the imperfection amplitude is large (μ = 0.2 for this numerical example), the load-edge shortening curve increases monotonically without the extreme point. The variation in the critical load with the imperfection amplitude is presented in Figure 6 and the buckling modes for different imperfection amplitude are also shown. It can be observed that the critical load significantly decreases with the increasing imperfection amplitude, and decreases faster when the imperfection amplitude is small. Furthermore, the imperfection amplitude also has influence on the buckling mode.
The load-edge shortening curves of imperfect NFRC cylindrical shells for different aging time and fiber content are presented in Figure 7 and Figure 8, respectively. It can be observed that the load-end shortening curve decrease with the increase in aging time, indicating that moisture absorption and hygrothermal aging have a great influence on the stability of NFRC cylindrical shells. It can be seen from Figure 8a that when the aging time is 0 (neglecting the moisture absorption and hygrothermal aging effect), the load-end shortening curve moves up with the increase in fiber content. It is interesting, as shown in Figure 8b, that when the aging time is 100 h, the load-end shortening curve for high fiber content is lower. Here, it should be mentioned that the research aim of this paper is to qualitatively analyze the influence of moisture absorption and hygrothermal aging effect. The aging speed will be related to temperature and humidity (the coefficients in the Equations (1)–(6) will be different under different humidity and temperature), but the overall trend is similar.
In order to further discuss the effect of moisture absorption and hygrothermal aging on the stability of NFRC cylindrical shells, the variation in the critical load with aging time and fiber content is presented in Figure 9 and Figure 10, respectively, and the variation of the imperfection sensitivity with aging time and fiber content is presented in Figure 11 and Figure 12, respectively. The imperfection sensitivity is defined as λ = Σclcl0, where Σcl0 and Σcl are the critical loads for the perfect and imperfect shells, respectively. The imperfection amplitude for Figure 9, Figure 10, Figure 11 and Figure 12 is μ = 0.05. It is observed from Figure 9 that the critical load decreases with the increase in aging time, and the decrease is faster when the fiber content is large. It can be seen from Figure 10 that when the moisture absorption and hygrothermal aging is neglected, the critical load increases linearly with the increase in fiber content. When the aging time is t = 50 h, the increase in the critical load with the fiber content is very small. When the aging time is t = 100 h, the critical load first increases and then decreases with the increasing fiber content. It can be observed from Figure 11 and Figure 12 that the moisture absorption and hygrothermal aging has a very small influence on the imperfection sensitivity. The imperfection sensitivity increases slightly with increasing aging time and decreases slightly with increasing fiber content.

6. Conclusions

The nonlinear stability of natural-fiber-reinforced composite cylindrical shells with initial geometric imperfection and moisture absorption and hygrothermal aging is investigated. The nonlinear governing equations are established based on the high-order shear deformation theory. The load-edge shortening curves for different imperfection amplitudes are obtained by the Galerkin method. The effects of key influencing factors on the stability of NFRC cylindrical shells are discussed. The main conclusions are as follows: (i) NFRC cylindrical shells are sensitive to the initial geometric imperfection and small geometric imperfections will reduce the critical load; (ii) the moisture absorption and hygrothermal aging has a great influence on the stability of NFRC cylindrical shells. The critical load of NFRC cylindrical shells significantly decreases with increasing aging time; (iii) when the aging time is long, the critical load of NFRC cylindrical shells with high fiber content is even lower than those with low fiber content; (iv) the influence of moisture absorption and hygrothermal aging effect on the critical load is very obvious, while the influence on the imperfection sensitivity is very small.

Author Contributions

Conceptualization, H.Z. and Z.Z.; methodology, H.Z. and H.B.; software, H.Z.; validation, H.Z.; formal analysis, H.Z.; investigation, H.Z. and Z.Z; resources, H.B.; data curation, H.Z.; writing—original draft preparation, H.Z.; writing—review and editing, H.B. and Z.Z.; visualization, H.Z.; supervision, H.B.; project administration, H.B. and Z.Z.; funding acquisition, H.B. and Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Science and Technology Program of Liaoning Province [2021JH4/10100061].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kerni, L.; Singh, S.; Patnaik, A.; Kumar, N. A review on natural fiber reinforced composites. Mater. Today Proc. 2020, 28, 1616–1621. [Google Scholar] [CrossRef]
  2. Chichane, A.; Boujmal, R.; Barkany, A.E. Bio-composites and bio-hybrid composites reinforced with natural fibers: Review. Mater. Today Proc. 2022; in press. [Google Scholar] [CrossRef]
  3. Pan, Y.; Zhong, Z. A nonlinear constitutive model of unidirectional natural fiber reinforced composites considering moisture absorption. J. Mech. Phys. Solids 2014, 69, 132–142. [Google Scholar] [CrossRef]
  4. Rohit, K.; Dixit, S. A review—Future aspect of natural fiber reinforced composite. Polym. Renew. Resour. 2016, 7, 43–59. [Google Scholar] [CrossRef]
  5. Nayak, S.Y.; Sultan, M.T.H.; Shenoy, S.B.; Kini, C.R.; Samant, R.; Shah, A.U.M.; Amuthakkannan, P. Potential of natural fibers in composites for ballistic applications—A review. J. Nat. Fibers 2022, 19, 1648–1658. [Google Scholar] [CrossRef]
  6. Li, X.; Tabil, L.G.; Panigrahi, S. Chemical treatments of natural fiber for use in natural fiber-reinforced composites: A review. J. Polym. Environ. 2007, 15, 25–33. [Google Scholar] [CrossRef]
  7. Fuqua, M.A.; Huo, S.; Ulven, C.A. Natural fiber reinforced composites. Polym. Rev. 2012, 52, 259–320. [Google Scholar] [CrossRef]
  8. Lu, X.; Zhang, M.; Rong, M.; Shi, G.; Yang, G. All-plant fiber composites. II: Water absorption behavior and biodegradability of unidirectional sisal fiber reinforced benzylated wood. Polym. Compos. 2003, 24, 367–379. [Google Scholar] [CrossRef]
  9. Espert, A.; Vilaplana, F.; Karlsson, S. Comparison of water absorption in natural cellulosic fibres from wood and one-year crops in polypropylene composites and its influence on their mechanical properties. Compos. Part A 2004, 35, 1267–1276. [Google Scholar] [CrossRef]
  10. Khakpour, H.; Ayatollahi, M.R.; Akhavan-Safar, A.; da Silva, L.F.M. Mechanical properties of structural adhesives enhanced with natural date palm tree fibers: Effects of length, density and fiber type. Compos. Struct. 2020, 237, 111950. [Google Scholar] [CrossRef]
  11. Akhavan-Safar, A.; Salamat-Talab, M.; Delzendehrooy, F.; Zeinolabedin-Beygi, A.; da Silva, L.F.M. Effects of natural date palm tree fibers on mode II fracture energy of E-glass/epoxy plain-woven laminated composites. J. Braz. Soc. Mech. Sci. Eng. 2022, 44, 457. [Google Scholar] [CrossRef]
  12. Dalla Libera, V.; Teixeira, L.A.; Amico, S.C.; da Luz, S.M. Processing, thermal and mechanical properties of composite laminates with natural fibers prepregs. Polym. Polym. Compos. 2022, 30, 09673911221087591. [Google Scholar] [CrossRef]
  13. Ramaswamy, S.; Shariff, Z.A.; Munaf, A.A.; Leno, I.J.; Gnanaraj, S.J.P.; Jeshurun, S.B. Study on application of higher order lamination plate theory over various applications of natural fiber cross-ply composites. Mater. Today Proc. 2022, 60, 822–826. [Google Scholar] [CrossRef]
  14. Pawar, K.S.; Bagha, A.K.; Bahl, S.; Agrawal, M.K. Experimental investigation for the dynamic characteristics of short natural fiber reinforced composite materials. Indian J. Eng. Mat. Sci. 2022, 29, 366–377. [Google Scholar]
  15. Nandhakumar, S.; Mithresh Kanna, K.; Mohammed Riyas, A.; Nathin Bharath, M. Experimental investigations on natural fiber reinforced composites. Mater. Today Proc. 2021, 37, 2905–2908. [Google Scholar] [CrossRef]
  16. Patel, N.; Jain, P. An investigation on mechanical properties in randomly oriented short natural fiber reinforced composites. Mater. Today Proc. 2021, 37, 469–479. [Google Scholar] [CrossRef]
  17. Gupta, U.S.; Dharkar, A.; Dhamarikar, M.; Choudhary, A.; Wasnik, D.; Chouhan, P.; Tiwari, S.; Namdeo, R. Study on the effects of fiber orientation on the mechanical properties of natural fiber reinforced epoxy composite by finite element method. Mater. Today Proc. 2021, 45, 7885–7893. [Google Scholar] [CrossRef]
  18. Lim, J.H.; Ratnam, M.M.; Abdul Khalil, H.P.S. An experimental and finite element analysis of the static deformation of natural fiber-reinforced composite beam. Polym. Test. 2003, 22, 169–177. [Google Scholar] [CrossRef]
  19. Bavan, D.S.; Kumar, G.C.M. Finite Element analysis of a natural fiber (maize) composite beam. J. Eng. 2013, 2013, 450381. [Google Scholar] [CrossRef] [Green Version]
  20. Jirawattanasomkul, T.; Likitlersuang, S.; Wuttiwannasak, N.; Ueda, T.; Zhang, D.; Shono, M. Structural behaviour of pre-damaged reinforced concrete beams strengthened with natural fibre reinforced polymer composites. Compos. Struct. 2020, 244, 112309. [Google Scholar] [CrossRef]
  21. Rajeshkumar, G.; Hariharan, V. Free vibration characteristics of phoenix Sp fiber reinforced polymer matrix composite beams. Procedia Eng. 2014, 97, 687–693. [Google Scholar] [CrossRef]
  22. Rajesh, M.; Pitchaimani, J.; Rajini, N. Free vibration characteristics of banana/sisal natural fibers reinforced hybrid polymer composite beam. Procedia Eng. 2016, 144, 1055–1059. [Google Scholar] [CrossRef] [Green Version]
  23. Senthil Kumar, K.; Siva, I.; Jeyaraj, P.; Winowlin Jappes, J.T.; Amico, S.C.; Rajini, N. Synergy of fiber length and content on free vibration and damping behavior of natural fiber reinforced polyester composite beams. Mater. Des. 2014, 56, 379–386. [Google Scholar] [CrossRef]
  24. Kuriakose, V.M.; Rohith Sai, P.; Sanjaay Kumar, M.; Sreehari, V.M. Influence of CNT fillers in the vibration characteristics of natural fiber reinforced composite plates. Compos. Struct. 2022, 282, 115012. [Google Scholar] [CrossRef]
  25. Selvaraj, R.; Maneengam, A.; Sathiyamoorthy, M. Characterization of mechanical and dynamic properties of natural fiber reinforced laminated composite multiple-core sandwich plates. Compos. Struct. 2022, 284, 115141. [Google Scholar] [CrossRef]
  26. Wang, K.F.; Wang, B.L. A mechanical degradation model for bidirectional natural fiber reinforced composites under hydrothermal aging and applying in buckling and vibration analysis. Compos. Struct. 2018, 206, 594–600. [Google Scholar] [CrossRef]
  27. Xu, C.H.; Rong, D.L.; Zhou, Z.H.; Deng, Z.C.; Lim, C.W. Vibration and buckling characteristics of cracked natural fiber reinforced composite plates with corner point-supports. Eng. Struct. 2020, 214, 110614. [Google Scholar] [CrossRef]
  28. Zhang, J.L.; Ni, Y.W.; Li, Q.D.; Tong, Z.Z.; Zhou, Z.H.; Xu, X.S. A symplectic approach for buckling analysis of natural fiber reinforced composite shells under hygrothermal aging. Appl. Math. Mech. 2021, 42, 1238–1247. [Google Scholar]
  29. Meng, Z.; Luo, X.; Zhou, H. Lightweight design of arcuately stiffened cylindrical shells based on smeared stiffener method and active learning strategy. Thin-Walled Struct. 2022, 174, 109167. [Google Scholar] [CrossRef]
  30. Wagner, H.N.R.; Hühne, C. Robust knockdown factors for the design of cylindrical shells under axial compression: Potentials, practical application and reliability analysis. Int. J. Mech. Sci. 2018, 135, 410–430. [Google Scholar] [CrossRef]
  31. Koiter, W.T. On the Stability of Elastic Equilibrium; National Aeronautics and Space Administration: Washington, DC, USA, 1967. [Google Scholar]
  32. Narayana, Y.V.; Gunda, J.B.; Reddy, P.R.; Markandeya, R. Non-linear buckling and post-buckling analysis of cylindrical shells subjected to axial compressive loads: A study on imperfection sensitivity. Nonlinear Eng. 2013, 2, 83–95. [Google Scholar] [CrossRef]
  33. Sun, Y.; Tian, K.; Li, R.; Wang, B. Accelerated Koiter method for post-buckling analysis of thin-walled shells under axial compression. Thin-Walled Struct. 2020, 155, 106962. [Google Scholar] [CrossRef]
  34. Shairyat, M.; Eslami, M.R. Dynamic buckling and post-buckling of imperfect cylindrical shells under mechanical and thermal loads, based on the three-dimensional theory of elasticity. J. Appl. Mech. 1999, 66, 476–484. [Google Scholar] [CrossRef]
  35. Eslami, M.R.; Shariyat, M.; Shakeri, M. Layerwise Theory for dynamic buckling and postbuckling of laminated composite cylindrical shells. AIAA J. 1998, 36, 1874–1882. [Google Scholar] [CrossRef]
  36. Shen, H.S. Buckling and Postbuckling of anisotropic laminated cylindrical shells with piezoelectric fiber reinforced composite actuators. Mech. Adv. Mater. Struct. 2010, 17, 268–279. [Google Scholar] [CrossRef]
  37. Shen, H.S. Postbuckling of functionally graded fiber reinforced composite laminated cylindrical shells, Part I: Theory and solutions. Compos. Struct. 2012, 94, 1305–1321. [Google Scholar] [CrossRef]
  38. Shen, H.S. Postbuckling of functionally graded fiber reinforced composite laminated cylindrical shells, Part II: Numerical results. Compos. Struct. 2012, 94, 1322–1332. [Google Scholar] [CrossRef]
  39. Shen, H.S.; Xiang, Y. Postbuckling behavior of functionally graded graphene-reinforced composite laminated cylindrical shells under axial compression in thermal environments. Comput. Methods Appl. Mech. Eng. 2018, 330, 64–82. [Google Scholar] [CrossRef]
  40. Shen, H.S.; Xiang, Y. Effect of negative Poisson’s ratio on the postbuckling behavior of axially compressed FG-GRMMC laminated cylindrical shells surrounded by an elastic medium. Eur. J. Mech. A Solids 2021, 88, 104231. [Google Scholar] [CrossRef]
  41. Babaei, H.; Kiani, Y.; Eslami, M.R. Buckling and post-buckling analysis of geometrically imperfect FGM pin-ended tubes surrounded by nonlinear elastic medium under compressive and thermal loads. Int. J. Struct. Stab. Dyn. 2019, 19, 1950089. [Google Scholar] [CrossRef]
  42. Li, Z.M.; Liu, T.; Yang, D.Q. Postbuckling behavior of shear deformable anisotropic laminated cylindrical shell under combined external pressure and axial compression. Compos. Struct. 2018, 198, 84–108. [Google Scholar] [CrossRef]
  43. Foroutan, K.; Shaterzadeh, A.; Ahmadi, H. Static and dynamic postbuckling analysis of imperfect SSFG cylindrical shells surrounded by nonlinear elastic foundation subjected to an axial compression. Mech. Adv. Mater. Struct. 2022, 29, 1769–1781. [Google Scholar] [CrossRef]
  44. Reddy, J.N.; Liu, C.F. A higher-order shear deformation theory of laminated elastic shells. Int. J. Eng. Sci. 1985, 23, 319–330. [Google Scholar] [CrossRef]
  45. Amabili, M. Nonlinear Vibrations and Stability of Shells and Plates; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  46. Yamaki, N. Elastic Stability of Circular Cylindrical Shells; Elsevier: Amsterdam, The Netherlands, 1984. [Google Scholar]
  47. Sun, J.B.; Zhu, S.B.; Tong, Z.Z.; Zhou, Z.H.; Xu, X.S. Post-buckling analysis of functionally graded multilayer graphene platelet reinforced composite cylindrical shells under axial compression. Proc. R. Soc. A Math. Phys. Eng. Sci. 2020, 476, 20200506. [Google Scholar] [CrossRef]
  48. Bisagni, C.; Cordisco, P. An experimental investigation into the buckling and post-buckling of CFRP shells under combined axial and torsion loading. Compos. Struct. 2003, 60, 391–402. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of NFRC cylindrical shell.
Figure 1. Schematic diagram of NFRC cylindrical shell.
Materials 15 06917 g001
Figure 2. Comparison of load-edge shortening curves of imperfect isotropic cylindrical shells.
Figure 2. Comparison of load-edge shortening curves of imperfect isotropic cylindrical shells.
Materials 15 06917 g002
Figure 3. Comparison of load-edge shortening curves of imperfect carbon-fiber-reinforced composite cylindrical shells.
Figure 3. Comparison of load-edge shortening curves of imperfect carbon-fiber-reinforced composite cylindrical shells.
Materials 15 06917 g003
Figure 4. Comparison of load-edge shortening curves of single-layer imperfect NFRC cylindrical shells.
Figure 4. Comparison of load-edge shortening curves of single-layer imperfect NFRC cylindrical shells.
Materials 15 06917 g004
Figure 5. Load-edge shortening curves of imperfect NFRC cylindrical shells for different imperfection amplitudes.
Figure 5. Load-edge shortening curves of imperfect NFRC cylindrical shells for different imperfection amplitudes.
Materials 15 06917 g005
Figure 6. Variation in the critical load of imperfect NFRC cylindrical shells with imperfection amplitude.
Figure 6. Variation in the critical load of imperfect NFRC cylindrical shells with imperfection amplitude.
Materials 15 06917 g006
Figure 7. Load-edge shortening curves of imperfect NFRC cylindrical shells for different aging times.
Figure 7. Load-edge shortening curves of imperfect NFRC cylindrical shells for different aging times.
Materials 15 06917 g007
Figure 8. Load-edge shortening curves of imperfect NFRC cylindrical shells for different fiber contents: (a): t = 0; (b): t = 100 h.
Figure 8. Load-edge shortening curves of imperfect NFRC cylindrical shells for different fiber contents: (a): t = 0; (b): t = 100 h.
Materials 15 06917 g008
Figure 9. Variation in the critical load of imperfect NFRC cylindrical shells with aging time.
Figure 9. Variation in the critical load of imperfect NFRC cylindrical shells with aging time.
Materials 15 06917 g009
Figure 10. Variation in the critical load of imperfect NFRC cylindrical shells with fiber content.
Figure 10. Variation in the critical load of imperfect NFRC cylindrical shells with fiber content.
Materials 15 06917 g010
Figure 11. Variation in the imperfection sensitivity of imperfect NFRC cylindrical shells with aging time.
Figure 11. Variation in the imperfection sensitivity of imperfect NFRC cylindrical shells with aging time.
Materials 15 06917 g011
Figure 12. Variation in the imperfection sensitivity of imperfect NFRC cylindrical shells with fiber content.
Figure 12. Variation in the imperfection sensitivity of imperfect NFRC cylindrical shells with fiber content.
Materials 15 06917 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, H.; Bai, H.; Zuo, Z. Nonlinear Stability of Natural-Fiber-Reinforced Composite Cylindrical Shells with Initial Geometric Imperfection Considering Moisture Absorption and Hygrothermal Aging. Materials 2022, 15, 6917. https://doi.org/10.3390/ma15196917

AMA Style

Zhang H, Bai H, Zuo Z. Nonlinear Stability of Natural-Fiber-Reinforced Composite Cylindrical Shells with Initial Geometric Imperfection Considering Moisture Absorption and Hygrothermal Aging. Materials. 2022; 15(19):6917. https://doi.org/10.3390/ma15196917

Chicago/Turabian Style

Zhang, Hongyu, Haifeng Bai, and Zhongyi Zuo. 2022. "Nonlinear Stability of Natural-Fiber-Reinforced Composite Cylindrical Shells with Initial Geometric Imperfection Considering Moisture Absorption and Hygrothermal Aging" Materials 15, no. 19: 6917. https://doi.org/10.3390/ma15196917

APA Style

Zhang, H., Bai, H., & Zuo, Z. (2022). Nonlinear Stability of Natural-Fiber-Reinforced Composite Cylindrical Shells with Initial Geometric Imperfection Considering Moisture Absorption and Hygrothermal Aging. Materials, 15(19), 6917. https://doi.org/10.3390/ma15196917

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop