Next Article in Journal
An Integrated Model for Acid Fracturing without Prepad Considering Wormhole Growth
Previous Article in Journal
Performance Evaluation on Open-Graded Friction Course Reinforced by Double-Adding Fibers Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study of the Buckling Response of Stiffened FG Graphene-Reinforced Multilayer Composite Cylindrical Panels

1
School of Statistics and Big Data, Zhengzhou College of Finance and Economics, Zhengzhou 450044, China
2
Department of Innovation Engineering, University of Salento, 73100 Lecce, Italy
3
Department of Mechanical Engineering, University of Eyvanekey, Eyvanekey 99888-35918, Iran
*
Authors to whom correspondence should be addressed.
Processes 2024, 12(3), 430; https://doi.org/10.3390/pr12030430
Submission received: 9 October 2023 / Revised: 17 January 2024 / Accepted: 18 February 2024 / Published: 20 February 2024
(This article belongs to the Section Materials Processes)

Abstract

:
The present research aims at determining the axial buckling load of stiffened multilayer cylindrical shell panels made of functionally graded graphene-reinforced composites (FG-GPL RCs). Rings and stringers are applied as stiffening tools for shell panels, whose elastic properties are determined according to the Halpin–Tsai relations. The virtual work principle and finite element approach are implemented here, according to a first-order shear deformation theory (FSDT) and Lekhnitskii smeared stiffener approach, in order to determine the governing equations of the stability problem. Four different dispersions of nanofillers are assumed in the thickness direction, including the FG-X, FG-A, FG-O, and UD distributions. A large systematic investigation considers the effect of different geometric and material parameters on the buckling loads and mode shapes of the stiffened FG-GPL RC cylindrical shell panel, primarily the dispersion and weight fractions of the nanofiller, the number of rings and stringers, and the boundary conditions, with useful insights for design purposes.

1. Introduction

In the last decades, a great effort by scientists has been devoted to increase the mechanical performances and durability of many structural members, primarily resorting to polymeric or metal nanocomposite materials, where nanoparticles can be dispersed within a polymeric or metallic matrix, with different distributions, in line with design requirements. The uniform dispersion of nanoparticles as reinforcement phase is desirable, even though their high level of energy can cause some agglomeration phenomena, with low nanoparticle weight fractions. In addition, there are some practical procedures, such as the use of an ultrasonic bath instead of traditional mixing, to ensure the uniform dispersion of nanoparticles in the matrix. Another way to enhance the mechanical performance of a structure stems from the application of ribs and stringers as strengthening members. In addition, the connection type between stiffeners and original structure plays a key role in the final strength of the assembled configuration, due to the presence of stress concentrations near to the connection points. Such connections can include glue, welding, screws, or rivets that can significantly affect the overall response of the stiffened structure.
There are a lot of works in the literature that have focused on different strengthening strategies for the improvement of the mechanical performances of structures under different loading conditions, both in a static and dynamic sense. Among GPL-reinforced applications, Mollaei et al. [1] applied the 3D finite element elasticity solution and virtual work principle to assess the buckling response of FG-GPL RC cylindrical shell panels under a torsional load. Kiani [2], instead, employed a generalized differential quadrature method (GDQM) and an equilibrium criterion to determine the buckling load of FG-GPL conical shells under a thermo-mechanical loading condition, accounting for a FSDT for shell structures. Lei et al. [3] used a Ritz approach to estimate the buckling load of FG-GPL RC laminated plates in thermal conditions applying the FSDT. Shahgholian et al. [4,5] also applied the Hamiltonian principle, FSDT basics, and Rayleigh–Ritz approach to evaluate the buckling [4] and torsional buckling [5] loads of FG porous hollow cylinders reinforced by GPLs. The same FSDT theoretical basics were also applied by Dong et al. [6] to analytically evaluate the buckling load of FG cylinders with porosities and reinforced by GPLs. The same Ritz approach was also employed by Chen et al. [7] in order to control the vibration and the geometrically nonlinear buckling behavior of FG-GPL nanocomposite porous beams. In another work [8], the same authors applied the Ritz approach to determine the elastic buckling and natural frequency response of composite beams reinforced by GPLs, including FG porosities. Zhang et al. [9] presented an analytical solution for the thermo-mechanical analysis of porous FG graphene-reinforced cylindrical panels using an improved third-order shear deformation theory. Sobhani et al. [10], instead, analyzed the natural frequencies of hybrid porous nanocomposite-joined hemispherical–cylindrical–conical shells based on a FSDT and GDQM procedure. At the same time, Yang et al. [11] combined the FSDT and the Chebyshev–Ritz approach to perform a comprehensive investigation of the buckling response for porous FG-GPL nanocomposite rectangular plates, whereas Zhou et al. [12] evaluated the buckling load of FG porous cylinders reinforced with GPLs, based on a classical theory of shells and higher-order shear deformation theory (HSDT). A novel numerical model was introduced by Ansari et al. [13] to estimate the buckling load of porous FG-GPL plate-type structures, whose governing measures were determined based on the Hamiltonian principle and HSDT. Similarly, Yaghoobi and Taheri [14] proposed an analytical procedure to study the buckling behavior of sandwiched rectangular plates with a FG porous core reinforced with GPLs, whereas Tao and Dai [15] applied the isogeometric approach to accurately study the post-buckling behavior of cylindrical shell panels with a FG porous core reinforced by GPL nanofillers. This last computational strategy was also proposed by Nguyen et al. [16] to determine the buckling load and natural frequency of FG-GPL porous rectangular plates according to HSDTs.
For further applications of the HSDT and finite element approach for the buckling and post-buckling behavior of different shells and panels, the reader is referred to Refs. [17,18,19,20,21,22,23], whose limit cases can be represented by curved beams [24,25,26], as well as by a classical shell theory, or thin-walled shell theory as proposed in [27,28], respectively, for shell structures with stiffeners. For these last structural cases, we briefly recall some relevant works from the recent literature focusing on the influence of stiffeners such as ribs and stringers on the buckling response of shell members. For example, Shahani and Kiarasi [29] numerically and experimentally investigated the buckling load of stiffened cylindrical shells with cut outs under a uniform axial loading. Similarly, Wu et al. [30] focused on the buckling response of composite stiffened panels under an impact force, and Tran et al. [31] numerically computed the buckling force of stiffened curved panels under an axial force. A classical finite element approach was followed by Grondin et al. [32] to investigate the buckling response of stiffened steel rectangular plates, whereas Reddy et al. [33] proposed a layerwise theory for the buckling study of stiffened hollow cylinders. An interesting review paper about the stability response of isotropic and composite stiffened panels can be found in Ref. [34], whereas some optimization procedures are proposed in Refs. [35,36,37,38,39] to improve the mechanical response of shell structures with different shapes in the presence of stiffeners. One of the most common computational tools to treat similar problems relies on the classical finite element method (FEM), as implemented by Ifayefunmi and Ruhan [40] to evaluate the buckling response of stiffened joined conical–cylindrical shell-type structures under an axial force, but also to model the buckling of grid-stiffened composite conical shells and stiffened shell panel structures [41,42].
Despite the large use of classical FEM, there are few works theoretically defining the buckling load based on a Lekhnitskii smeared technique, where the connection between the original structure and stiffeners is considered ideal in absence of rivets, screws, and welding. For example, Duc and Thang [43] considered the nonlinear buckling response of imperfect cylindrical shells with stiffeners made of FG materials (FGMs) in thermal conditions. In another investigation [44], the same authors proposed a nonlinear analysis of stiffened FGM cylinders resting on an elastic medium and subjected to an axial force. Similarly, Duc [45] studied the nonlinear vibration response of imperfect FGM doubly curved panels with stiffeners resting on an elastic medium. Duc and Cong [46] applied the classical shell theory together with the Galerkin approach to study the nonlinear thermal response of stiffened conical shells made of FGMs on an elastic substrate.
Based on the literature overview, it has been clearly shown that nanoparticle-based reinforcements as well as the application of orthogonal stiffeners such as rings and stringers can singly provide a great effect on the buckling response of shell structures. However, the combined effect of both aspects has been not investigated yet, as proposed in the present work for the buckling study of stiffened FG graphene-reinforced composite multilayer cylindrical shell panels. In such a context, four different distributions of GPL reinforcements are assumed here throughout the thickness direction of the cylindrical panels, whereas a Lekhnitskii smeared technique is considered to model the connection between FG-GPL shell panels and stiffeners. Moreover, the Halpin–Tsai relations are established here to estimate the equivalent mechanical properties of the shell panel. The FSDT, together with the principle of virtual work and finite element procedure, are also adopted to determine the governing equations of the problem. A systematic investigation is performed to check for the effect of several parameters, including the weight fractions and patterns of GPLs, different boundary conditions, geometries, and numbers of stiffeners, on the buckling response of stiffened FG-GPL cylindrical shell panels. The work is organized as follows. After this brief introduction, we provide the main basic principles of the theoretical model in Section 2, starting with the geometrical and mechanical definition of the selected problem. The governing equations of the problem are then solved numerically according to a finite element approximation, as proposed in Section 3, from a theoretical perspective, and they are applied numerically in Section 4 within a large parametric investigation. The main concluding remarks can be found in Section 5, with useful suggestions from both a theoretical and a design standpoint for different engineering applications of novel composite structural members optimized by simply tuning the FG graphene reinforcement phase.

2. Theoretical Formulation

2.1. Geometrical and Mechanical Properties

Let us consider a stiffened FG-GPL multilayer cylindrical shell panel under a compressive axial force (Figure 1), with length L , span angle β , radius R , and thickness h . As visible in the same figure, let us assume the (x, θ, z) system of coordinates, accounting for the axial, hoop, and radial axes, respectively. We account for a shell panel stiffened with longitudinal stiffeners (stringers) and circumferential stiffeners (rings), reinforced by different GPL patterns, as reported in Figure 1.
The multilayer GPL nano-composite shell panel is characterized by perfectly bonded GPL-RC layers, where each layer combines GPLs as a nanofiller with an isotropic polymer matrix. Both rings and stringers are made of the same isotropic polymeric matrix as assumed for the shell panel in order to avoid any possible stress concentration.
GPLs can be distributed in the polymer matrix uniformly or non-uniformly. In the last case, the weight fraction of the nanofiller (GPL) is functionally graded across the radial direction of the shell panel, which is made of n L layers. More in detail, we consider four different volume fraction distributions of FG GPL-RCs, namely, a uniform pattern (UD) in addition to FG-A, O, and X patterns. In the uniform case, the GPL weight fraction remains constant for each layer, yielding to a homogeneous isotropic GPL-RC structure. For FG patterns, instead, the weight fraction of GPLs varies linearly across the thickness of the panel. More in detail, for FG-X distributions, the weight fraction of GPLs assumes the maximum value at the inner and outer layers, whereas for a FG-O distribution, the weight fraction obtains the maximum value at the mid-layers. Moreover, for an A-GPLRC-type distribution, the weight fraction is maximum at the inner layer and decreases continuously for an increased distance from it, assuming the lowest quantity of GPLs at the external surface of the shell. For different GPL distributions, the equivalent weight fraction can be evaluated as [47,48,49,50,51].
{ U D   :   V G P L ( k ) = V G P L * F G X   :   V G P L ( k ) = 2 V G P L * | 2 k n L 1 | n L F G O   :   V G P L ( k ) = 2 V G P L * ( 1 | 2 k n L 1 | n L )   F G Λ   :   V G P L ( k ) = V G P L * ( 2 ( n L k + 1 ) 1 n L )  
Here, V G P L ( k ) refers to the volume content of GPL in each layer of the panel, with k varying from 1 to n L ; V G P L * stands for the volume fraction of nanofillers, which can be evaluated as follows:
V G P L * = Δ G P L ρ m Δ G P L ρ m + ρ G P L Δ G P L ρ G P L
where Δ G P L is the weight fraction of the nanofillers, and ρ G P L and ρ m stand for the mass density of the nanofillers and polymeric matrix, respectively.
According to Halpin–Tsai micromechanics [44,45,46,47,48], the elasticity modulus of the panel can be described as
E = 3 8 ( 1 + ε L G P L η L G P L V G P L 1 η L G P L V G P L ) E m + 5 8 ( 1 + ε W G P L η W G P L V G P L 1 η W G P L V G P L ) E G P L  
ε L G P L = 2 l G P L t G P L
ε W G P L = 2 w G P L t G P L
η L G P L = E G P L E m E G P L + ε L G P L E m
η W G P L = E G P L E m E G P L + ε W G P L E m
where E m and E G P L stand for the elasticity modulus for the matrix and nanofillers, respectively. In addition, l G P L , w G P L , t G P L , and V G P L refer to the length, width, thickness, and volume fraction of the nanofillers, respectively. The rule of mixture is also employed to evaluate the Poisson’s ratio and mass density of the GPL-RC:
ρ = ρ G P L V G P L + ρ m ( 1 V G P L )
v = v G P L V G P L + v m ( 1 V G P L )
where ν G P L and ν m refer to the Poisson’s ratio of the nanofillers and matrix, respectively. Thus, the shear modulus G of the shell is determined as
G = E 2 ( 1 + ν )

2.2. Mathematical Problem

A FSDT is considered here to describe the displacement field, and a Lekhnitskii smeared stiffeners assumption is employed to define the relation between the main structure and its rings and stringers. Applying the FSDT of shells, the displacement components of the shell are defined as
u ( x , θ , z , t ) = u 0 ( x , θ , t ) + z ϕ x ( x , θ , t ) ν ( x , θ , z , t ) = ν 0 ( x , θ , t ) + z ϕ θ ( x , θ , t ) w ( x , θ , z , t ) = w 0 ( x , θ , t )
where u 0 , v 0 , w 0 are the kinematic components of the mid-plane of the panel along the axial, circumferential, and radial directions, respectively. Also, ϕ x ( x , θ , t ) ,   ϕ θ ( x , θ , t ) are the rotations of the mid-plane around the x and θ axes, respectively. The kinematic relations are also defined as
ε x = ε x 0 + z k x ε θ = ε θ 0 + z k θ γ x θ = γ x θ 0 + z k x θ γ θ z = γ θ z 0 γ x z = γ x z 0
where
ε x 0 = u 0 x ; k x = ϕ x x ; ε θ 0 = v 0 R θ + w ° R ; k θ = ϕ θ R θ ; γ x θ 0 = u ° R θ + ν ° x ; k x θ = ϕ x R θ + ϕ θ x ; γ θ z = ϕ θ + w ° R θ ν ° R ; γ x z 0 = ϕ x + w ° x ;
Based on the Lekhnitskii approach, the resultants of forces and moments are defined as [52]
N x = ( A 11 + E s A s s s ) ε x 0 + A 12 ε θ 0 + ( B 11 + E s A s z s s s ) k x + B 12 k θ N θ = A 12 ε x 0 + ( A 22 + E r A r s r ) ε θ 0 + B 12 k x + ( B 22 + E r A r z r s r ) k θ N x θ = A 66 γ x θ 0 + B 66 k x θ M x = ( B 11 + E s A s z s s s ) ε x 0 + B 12 ε θ 0 + ( D 11 + E s I s s s ) k x + D 12 k θ M θ = B 12 ε x 0 + ( B 22 + E r A r z r s r ) ε θ 0 + D 12 k x + ( D 22 + E r I r s r ) k θ M x θ = B 66 γ x θ 0 + D 66 k x θ { Q θ = k s ( A 44 + G r A r s r ) γ θ z Q x = k s ( A 55 + G s A s s s ) γ x z
with
A 11 = A 22 = h 2 h 2 E 1 υ 2 d z ; B 11 = B 22 = h 2 h 2 E 1 υ 2 z d z ; D 11 = D 22 = h 2 h 2 E 1 υ 2 z 2 d z ; A 12 = h 2 h 2 υ E 1 υ 2 d z ; B 12 = h 2 h 2 υ E 1 υ 2 z d z ; D 12 = h 2 h 2 υ E 1 υ 2 z 2 d z ; A 66 = h 2 h 2 E 2 [ 1 + υ ] d z ; B 66 = h 2 h 2 E 2 [ 1 + υ ] z d z ; D 66 = h 2 h 2 E 2 [ 1 + υ ] z 2 d z ; A 44 = A 55 = h 2 h 2 E 2 [ 1 + υ ] d z ;
I s = b s h s 3 12 + A s z s 2 ; I r = b r h r 3 12 + A r z r 2 ; z s = ± h s + h 2 ; z r = ± h r + h 2
The normal and tangential elasticity moduli of stringers and rings are defined as E s ,   G s and E r , G r , respectively, while assuming a shear correction factor k s = 5 / 6 . Both ribs and stringers as well as the main panel are considered to be perfectly joined, and each part is made of the same material. A uniaxial stress state is also assumed for both rings and stringers, which are defined geometrically by height h r , h s , and width b r , b s , respectively, for a global cross-section area A r , A s .
Moreover, the distances between two adjacent stringers and each ring are denoted by S s and S r , respectively, whereas the distances between the centroid of each stringer and ring from the mid-surface of the panel are labeled as z s and z r , respectively. The resultants of force and moment in matrix form are defined as
{ N x N θ N x θ M x M θ M x θ } = [ A ¯ 11 A ¯ 12 0 B ¯ 11 B ¯ 12 0 A ¯ 12 A ¯ 22 0 B ¯ 12 B ¯ 22 0 0 0 A ¯ 66 0 0 B ¯ 66 B ¯ 11 B ¯ 12 0 D ¯ 11 D ¯ 12 0 B ¯ 12 B ¯ 22 0 D ¯ 12 D ¯ 22 0 0 0 B ¯ 66 0 0 D ¯ 66 ] { ε x 0 ε θ 0 γ x θ 0 k x k θ k x θ } , { Q θ Q x } = k s [ A ¯ 44 0 0 A ¯ 55 ] { γ θ z γ x z }
with
A ¯ = [ A ¯ 11 A ¯ 12 0 A ¯ 12 A ¯ 22 0 0 0 A ¯ 66 ] , B ¯ = [ B ¯ 11 B ¯ 12 0 B ¯ 12 B ¯ 22 0 0 0 B ¯ 66 ] , D ¯ = [ D ¯ 11 D ¯ 12 0 D ¯ 12 D ¯ 22 0 0 0 D ¯ 66 ] , E ¯ = k s [ A ¯ 44 0 0 A ¯ 55 ]
A ¯ 11 = A 11 + E s A s s s ; A ¯ 12 = A 12 ; A ¯ 22 = A 22 + E r A r s r ; A ¯ 66 = A 66 ; A ¯ 44 = A 44 + G r A r s r ; A ¯ 55 = A 55 + G s A s s s ; B ¯ 11 = B 11 + E s A s z s s s ; B ¯ 12 = B 12 ; B ¯ 22 = B 22 + E r A r z r s r ; B ¯ 66 = B 66 D ¯ 11 = D 11 + E s I s s s ; D ¯ 12 = D 12 ; D ¯ 22 = D 22 + E r I r s r ; D ¯ 66 = D 66
The principle of virtual work is now applied to determine the governing equations of the problem:
δ U δ W = 0 U = 1 2 s [ N x ε x 0 + N θ ε θ 0 + N x θ γ x θ 0 + M x k x + M θ k θ + M x θ k x θ + Q θ γ θ z + Q x γ x z ] R d θ d x , W = 1 2 s [ P x ( w 0 x ) 2 ] R d θ d x
with U the strain energy, and W the external work associated with a uniform compressive axial load applied on the bases of the shell. A homogenous isotropic material is assumed for all stiffeners, such that E s = E r = E m .
Hence, the strain energy for the stiffened structure is defined as
U = { N x ε x 0 + N θ   ε θ 0 + N x θ γ x θ 0 + M x K x + M θ K θ + M x θ K x θ + Q x γ x z + Q θ γ θ z } R   d x   d θ
which can be arranged in the following matrix form:
U = ( [ N x , N θ , N X θ ] [ ε   x 0 ε   θ 0 γ   x θ 0 ] + [ M x , M θ , M X θ ] [ K x K θ K x θ ] + [ Q x Q ϑ ] [ γ X Z γ θ Z ] ) R d x d θ
Also, the kinematic relations (12) and (13) can be defined in the following matrix form:
[ ε x ε θ γ x θ ] = [ ε   x 0 ε   θ 0 γ   x θ 0 ] + Z [ K x K θ K x θ ] = [ x 0 0 Z x 0 0 1 R θ 1 R 0 Z 1 R θ 1 R θ x 0 Z 1 R θ Z x ] [ U 0 V 0 W 0 φ X φ θ ] = d 1 Q [ γ   x z 0 γ θ z 0 ] = [ 0 0 x 1 0 0 1 R 1 R θ 0 1 ] [ U 0 V 0 W 0 φ X φ θ ] = d 2 Q , [ ε   x 0 ε   θ 0 γ   x θ 0 ] = [ x 0 0 0 0 0 1 R θ 1 R 0 0 1 R θ x 0 0 0 ] [ U 0 V 0 W 0 φ X φ θ ] = d 3 Q [ K x K θ K x θ ] = [ 0 0 0 x 0 0 0 0 0 1 R θ 0 0 0 1 R θ x ] [ U 0 V 0 W 0 φ X φ θ ] = d 4 Q , Q = [ U 0 V 0 W 0 φ X φ θ ]
Hence, by substitution of Equation (20) into Equation (19), the strain energy can be defined in variational form as
δ U = ( ( ( d 3 Q ) T A ¯ T + ( d 4 Q ) T B ¯ T ) ( d 3 δ Q ) + ( ( d 3 Q ) T B ¯ T + ( d 4 Q ) T D ¯ T ) ( d 4 δ Q ) + ( d 2 Q ) T e ¯ T ( d 2 δ Q ) ) R   d x   d θ

3. Finite Element Modeling

A numerical procedure based on classical finite elements is now applied to solve the governing equations of a stiffened FG-GPL multilayer cylindrical panel in a discrete manner. A two-dimensional, four-noded element with 20 degrees of freedom (DOFs) is adopted to discretize the model, where the shape functions in the x and θ axes are referred to a local-coordinate system ( ξ , η ) of the following type:
ξ = 2 ( x x c ) L ( e ) η = 2 ( θ θ c ) β ( e )
where 1 ξ , η 1 in the x and θ axes. Moreover, L ( e ) and β ( e ) are the length and span angle for each element, respectively. Also, θ c and x c are the circumferential and axial coordinates of the center for each element. The approximating functions for each element are defined in terms of natural coordinates and displacement components as follows:
{ Ψ i Ψ j Ψ k Ψ m } = 1 4 { ( 1 + ξ ) ( 1 η ) ( 1 + ξ ) ( 1 + η ) ( 1 ξ ) ( 1 + η ) ( 1 ξ ) ( 1 η ) }
( ( Ψ 1 0 0 Ψ 1 ) ( Ψ 4 0 0 Ψ 4 ) ) { u 01 v 01 w 01 φ x 1 φ θ 1 u 04 v 04 w 04 φ x 4 φ θ 4 } = ψ q ( e )  
where Ψ n ,   n = 1 , 2 , 3 , 4 are the shape function components, Ψ is the shape function matrix, and u 0 i , v 0 i , w 0 i , φ x i , and φ θ i stand for the nodal DOFs defined as
u 0 = i = 1 4 Ψ i U 0 i v 0 = i = 1 4 Ψ i V 0 i w 0 = i = 1 4 Ψ i W 0 i φ x = i = 1 4 Ψ i θ x i φ θ = i = 1 4 Ψ i θ θ i
By substitution of Equation (24) into Equations (17) and (21), the virtual work principle (17) can be rewritten as
Ω 0 e [ ( ( d 3 Ψ ) T A ¯ T ( d 3 Ψ ) + ( d 4 Ψ ) T B ¯ T ( d 3 Ψ ) + ( d 3 Ψ ) T B ¯ T ( d 4 Ψ ) + ( d 4 Ψ ) T D ¯ T ( d 4 Ψ ) + ( d 2 Ψ ) T e ¯ T d 2 P x Ψ ¯ 1     T Ψ ¯ 1 Ψ ) q ( e ) ] R d x d θ = 0
where
ψ ¯ 1 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ψ 1 x 0 0 0 0 ψ 4 x 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] 5 * 20  
In addition, by replacing d 2 Ψ = B 2 ,   d 3 Ψ = B 3 ,   d 4 Ψ = B 4 into Equation (25), the following equation is derived for a cylindrical panel element:
( k 1 + k 2 + k 3 P x k g ) ( e ) q ( e ) = 0
with
k 1 e = [ B 3     T A ¯     T + B 4     T B ¯     T   B 3 ] R   d x   d θ k 2 e = [ B 3     T B ¯     T B 4 + B 4     T D ¯     T B 4 ]   R   d x   d θ k 3 e = [ B 2     T e ¯     T B 2 ]   R   d x   d θ k g e = Ψ ¯ 1     T Ψ ¯ 1 R   d x   d θ
where k g is the geometric stiffness matrix, and P x is the buckling load.
Finally, by adding the stiffness and geometric stiffness matrices, we derive the following finite element equations of the problem:
( k 1 + k 2 + k 3 P x k g ) q = 0
The buckling analysis of the shell reverts to the eigenvalue problem, such that
det ( ( k 1 + k 2 + k 3 ) P x k g ) = 0
The following displacement boundary conditions are considered
For clamped supports (CC):
v 0 , w 0 , φ x , φ θ = 0 a t ( x = 0 , θ ) , ( x = L , θ )
And for simply supports (SS):
v 0 , w 0 = 0 a t ( x = 0 , θ ) , ( x = L , θ )

4. Numerical Investigation

4.1. Validation

The numerical investigation starts by considering a validation step, where the numerical results are first derived for homogenous stiffened cylindrical panels with rings and stringers and subjected to an axial load and SS boundary conditions, while comparing them to predictions from a commercial code. To this end, in such a preliminary analysis, we assume a null value of weight fraction for GPLs, while considering the following mechanical properties for the matrix: E m = 3   GPa ,   ρ m = 1200   kg / m 3 ,   ν m = 0.34 ,   Δ = 0   wt   %   . The geometry of the cylindrical panel is defined by means of L = 2 ,   R = 0.5 ,   β = 120 ° under the same assumption for the geometrical parameters of rings and stringers: h r = h s = 0.02 ,   b r = b s = 0.04 ,   S r = 0.5 ,   S s = 0.378 ,   N r = N s = 5 .  ( N r is the number of rings, and N s is the number of stringers). Table 1 summarizes the results from our formulation in terms of the first six buckling load values, compared to predictions from literature, with a very good accordance among them. This confirms the accuracy of the proposed approach, which enables the subsequent parametric analysis for tunable geometries of stiffeners and mechanical properties of the reinforcement phase. Figure 2 also depicts the corresponding first six buckling mode shapes for the SS-supported stiffened panel, involving different symmetric and asymmetric deformations, according to the selected boundary conditions.
In addition, to ensure the correct implementation of equivalent mechanical properties of FG-GPL structures, the critical buckling load of a multilayer FG-GPL cylindrical shell ( β = 360 ° ) under an axial load without stiffeners has been determined and compared with Ref. [53] (Table 2). To this end, we set Er and Es = 0 in the Matlab Code, and the geometry and mechanical properties are assumed according to Ref. [53]. The governing equations of [53] are based on the classical Donnell shell theory, and the results of the present study are based on FSDT. The comparison between the present results and Ref. [53] for different weight fractions of GPLs and various distributions of nanofillers shows an excellent agreement (Table 2).

4.2. Numerical Results

4.2.1. Convergence Study

A convergence study is now performed for the first buckling load of stiffened FG-GPL cylindrical shell panels, as summarized in the Table 3. In this case, the critical buckling load of a stiffened multilayer cylindrical panel reinforced by GPLs with simply supported boundary conditions is obtained for a different number of elements (GPL-X, ∆ = 1 wt. %, Ss = 0.378, Sr = 0.5). The convergence study is shown in Table 3, where accurate results are obtained by applying 30 and 45 elements in the circumferential and axial directions, respectively.

4.2.2. Numerical Results of Present Investigation

The numerical analysis continues by evaluating the first six buckling loads for a stiffened multilayer cylindrical panel reinforced by GPLs, accounting for the sensitivity of the axial buckling response to different reinforcement patterns and weight fractions, but also to different boundary conditions and numbers of rings and stringers. In this case, the stiffened FG-GPL multilayer cylindrical panel assumes the following mechanical properties: E m = 3   GPa , ν m = 0.34 , ρ m = 1200   kg / m 3 (for the matrix phase), and E G P L = 1.01   TPa ,   ν G P L = 0.186 ,   ρ G P L = 1062.5   kg / m 3 ,   w G P L = 1.5   μ m ,   l G P L = 2.5   μ m ,   t G P L = 1.5   nm (for the reinforcement phase). The same geometrical properties assumed before are considered once again in this case for the cylindrical panel with stiffeners.
Table 4 summarizes the effect of different GPL patterns on the buckling load of the shell, for a SS boundary condition, under the assumption of a weight fraction ∆ = 1 wt. %. Based on this table, the GPLX pattern seems to yield the highest buckling loads of the structure compared to other GPL patterns, whereas the GPLO pattern seems to obtain the lowest buckling loads. Such differences are approximately estimated equal to 40% as acceptable for design purposes among engineers. At the same time, when the concentration of the nanofiller phase at the external sides of the panel is higher than the concentration in the core, the structure increases its mechanical stiffness. Table 4 also describes the impact of GPL weight fractions on the buckling load of the structure. More specifically, by enhancing the weight fraction of nanofillers from 0 to 0.01, for different GPL patterns, the buckling load increases significantly (approximately 400% for a GPL-X pattern). The impact of an increased weight fraction of GPLs on the buckling load of the structure is more pronounced for a GPL-X distribution than other patterns, and the impact of an increased weight fraction of nanofillers on the buckling load for GPL-O is less pronounced than the other patterns. In Table 5, we also show the sensitivity of the buckling response to varying boundary conditions, under the same assumptions for S r ,   S s ,   N r ,   N s ,   Δ .  As visible from this table, a CC boundary condition yields higher buckling loads than a SS boundary condition, due to the increased stiffness of a CC boundary condition versus a SS structure, with a relative difference of about 15%. Table 6 also shows the influence of the number of stiffeners (rings and stringers) on the buckling load of the stiffened cylindrical panel for both boundary conditions, under the assumption of uniform distribution for GPLs (GPL-UD) and a weight fraction of nanofillers ∆ = 1 wt. %. Based on results, it is worth noticing that an increased number of rings and stringers, from five to ten, provides an increased buckling load of about 15% for both boundary conditions. The first six buckling mode shapes of the shell panel are represented in Figure 2 and Figure 3 for a SS and CC boundary condition, respectively, while keeping the same number of ribs and stringers, N r = N s = 5 , whereas the last two figures (Figure 4 and Figure 5) show the corresponding mode shapes for an increased number of stiffeners, N r = N s = 10 , with a clearly different deformed configuration just related to the varied number of stiffeners, under the same boundary condition. More in detail, based on a comparative evaluation of Figure 2 with Figure 4 and Figure 3 with Figure 5, an increased number of stiffeners clearly reduces the overall deformability of the structure, as highlighted in all symmetric and asymmetric deformed configurations, for each mode shape.

5. Conclusions

This work has focused on the buckling response of stiffened cylindrical panels reinforced with GPLs subjected to an axial loading. Four different patterns of GPLs, including GPL-X, GPL-O, GPL UD, and GPL-A, have been considered here across the shell thickness. The first-order shear deformation theory and finite element approach have permitted us to define and approximately solve the governing equations of the problem. A systematic investigation has been performed to assess the sensitivity of the response to different GPL patterns and weight fractions, but also to different boundary conditions and numbers of stringers and rings. The final remarks can be summarized as follows:
(a)
The maximum and minimum values of the buckling load stem from a GPL-X and GPL-O distribution, respectively.
(b)
The maximum variation in the buckling load for different GPL patterns and weight fractions is approximately equal to 40% and 400%, respectively.
(c)
The buckling load of the structure tends to increase for increased reinforcement weight fractions, especially for a GPL-X pattern, compared to other GPL distributions.
(d)
A CC boundary condition provides higher buckling loads (of approximately 15%) compared to the results from a SS boundary condition.
(e)
The buckling load can be increased by about 15% for each fixed boundary condition, for an increased number of rings and stringers from five to ten.
Our parametric results could be a valid reference for structural design optimization purposes, due to the capability of FG GPLs to change the stiffness and strength properties of the composite panels with different boundary conditions, by simply tuning the geometry of the stiffeners, as well as the reinforcement distribution and weight fraction.

Author Contributions

Conceptualization, F.T. and R.D.; Methodology, Z.L., F.T., R.D. and M.B.; Validation, Z.L., F.T., R.D. and M.B.; Formal analysis, Z.L. and M.B.; Investigation, F.T., R.D. and M.B.; Writing—original draft, Z.L. and M.B.; Writing—review & editing, F.T. and R.D.; Visualization, Z.L.; Supervision, F.T. and R.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mollaei, S.; Babaei, M.; Asemi, K. Torsional buckling of functionally graded graphene reinforced composite laminated cylindrical panel. Arch. Appl. Mech. 2023, 93, 427–435. [Google Scholar] [CrossRef]
  2. Kiani, Y. Buckling of functionally graded graphene reinforced conical shells under external pressure in thermal environment. Compos. Part B Eng. 2019, 156, 128–137. [Google Scholar] [CrossRef]
  3. Lei, Z.; Su, Q.; Zeng, H.; Zhang, Y.; Yu, C. Parametric studies on buckling behavior of functionally graded graphene-reinforced composites laminated plates in thermal environment. Compos. Struct. 2018, 202, 695–709. [Google Scholar] [CrossRef]
  4. Shahgholian-Ghahfarokhi, D.; Safarpour, M.; Rahimi, A. Torsional buckling analyses of functionally graded porous nanocomposite cylindrical shells reinforced with graphene platelets (GPLs). Mech. Based Des. Struct. Mach. 2021, 49, 81–102. [Google Scholar] [CrossRef]
  5. Shahgholian-Ghahfarokhi, D.; Rahimi, G.; Khodadadi, A.; Salehipour, H.; Afrand, M. Buckling analyses of FG porous nanocomposite cylindrical shells with graphene platelet reinforcement subjected to uniform external lateral pressure. Mech. Based Des. Struct. Mach. 2021, 49, 1059–1079. [Google Scholar] [CrossRef]
  6. Dong, Y.H.; He, L.W.; Wang, L.; Li, Y.H.; Yang, J. Buckling of spinning functionally graded graphene reinforced porous nanocomposite cylindrical shells: An analytical study. Aerosp. Sci. Technol. 2018, 82, 466–478. [Google Scholar] [CrossRef]
  7. Chen, D.; Yang, J.; Kitipornchai, S. Nonlinear vibration and postbuckling of functionally graded graphene reinforced porous nanocomposite beams. Compos. Sci. Technol. 2017, 142, 235–245. [Google Scholar] [CrossRef]
  8. Kitipornchai, S.; Chen, D.; Yang, J. Free vibration and elastic buckling of functionally graded porous beams reinforced by graphene platelets. Mater. Des. 2017, 116, 656–665. [Google Scholar] [CrossRef]
  9. Zhang, W.; Wang, C.; Wang, Y. Thermo-mechanical analysis of porous functionally graded graphene reinforced cylindrical panels using an improved third order shear deformable model. Appl. Math. Model. 2023, 118, 453–473. [Google Scholar] [CrossRef]
  10. Sobhani, E.; Arbabian, A.; Civalek, Ö.; Avcar, M. The free vibration analysis of hybrid porous nanocomposite joined hemispherical–cylindrical–conical shells. Eng. Comput. 2021, 38, 3125–3152.
  11. Yang, J.; Chen, D.; Kitipornchai, S. Buckling and free vibration analyses of functionally graded graphene reinforced porous nanocomposite plates based on Chebyshev-Ritz method. Compos. Struct. 2018, 193, 281–294. [Google Scholar] [CrossRef]
  12. Zhou, Z.; Ni, Y.; Tong, Z.; Zhu, S.; Sun, J.; Xu, X. Accurate nonlinear buckling analysis of functionally graded porous graphene platelet reinforced composite cylindrical shells. Int. J. Mech. Sci. 2019, 151, 537–550. [Google Scholar] [CrossRef]
  13. Ansari, R.; Hassani, R.; Gholami, R.; Rouhi, H. Buckling and Postbuckling of Plates Made of FG-GPL-Reinforced Porous Nanocomposite with Various Shapes and Boundary Conditions. Int. J. Struct. Stab. Dyn. 2021, 21, 2150063. [Google Scholar] [CrossRef]
  14. Yaghoobi, H.; Taheri, F. Analytical solution and statistical analysis of buckling capacity of sandwich plates with uniform and non-uniform porous core reinforced with graphene nanoplatelets. Compos. Struct. 2020, 252, 112700. [Google Scholar] [CrossRef]
  15. Tao, C.; Dai, T. Isogeometric analysis for postbuckling of sandwich cylindrical shell panels with graphene platelet reinforced functionally graded porous core. Compos. Struct. 2021, 260, 113258. [Google Scholar] [CrossRef]
  16. Nguyen, Q.H.; Nguyen, L.B.; Nguyen, H.B.; Nguyen-Xuan, H. A three-variable high order shear deformation theory for isogeometric free vibration, buckling and instability analysis of FG porous plates reinforced by graphene platelets. Compos. Struct. 2020, 245, 112321. [Google Scholar] [CrossRef]
  17. Priyanka, R.; Twinkle, C.M.; Pitchaimani, J. Stability and dynamic behavior of porous FGM beam: Influence of graded porosity, graphene platelets, and axially varying loads. Eng. Comput. 2021, 38, 4347–4366. [Google Scholar] [CrossRef]
  18. Barati, M.R.; Zenkour, A.M. Analysis of postbuckling of graded porous GPL-reinforced beams with geometrical imperfection. Mech. Adv. Mater. Struct. 2019, 26, 503–511. [Google Scholar] [CrossRef]
  19. Anirudh, B.; Ganapathi, M.; Anant, C.; Polit, O. A comprehensive analysis of porous graphene-reinforced curved beams by finite element approach using higher-order structural theory: Bending, vibration and buckling. Compos. Struct. 2019, 222, 110899. [Google Scholar] [CrossRef]
  20. Li, Z.; Zheng, J. Nonlinear stability of the encased functionally graded porous cylinders reinforced by graphene nanofillers subjected to pressure loading under thermal effect. Compos. Struct. 2020, 233, 111584. [Google Scholar] [CrossRef]
  21. Twinkle, C.M.; Pitchaimani, J. Free vibration and stability of graphene platelet reinforced porous nano-composite cylindrical panel: Influence of grading, porosity and non-uniform edge loads. Eng. Struct. 2021, 230, 111670. [Google Scholar]
  22. Tornabene, F.; Viscoti, M.; Dimitri, R. Higher order theories for the modal analysis of anisotropic doubly-curved shells with a three-dimensional variation of the material properties. Eng. Anal. Bound. Elem. 2024, 158, 486–519. [Google Scholar] [CrossRef]
  23. Wang, Z.W.; Tang, J.; Li, S.C.; He, X.H.; Zhou, C.Y. Research on Elastic and Elastics-Plastic Buckling Load of Cylindrical Shell with an Inclined through Crack under Axial Compressive Load. Materials 2023, 16, 6123. [Google Scholar] [CrossRef] [PubMed]
  24. Yas, M.H.; Rahimi, S. Thermal buckling analysis of porous functionally graded nanocomposite beams reinforced by graphene platelets using Generalized differential quadrature method. Aerosp. Sci. Technol. 2020, 107, 106261. [Google Scholar] [CrossRef]
  25. Gao, K.; Do, D.M.; Li, R.; Kitipornchai, S.; Yang, J. Probabilistic stability analysis of functionally graded graphene reinforced porous beams. Aerosp. Sci. Technol. 2020, 98, 105738. [Google Scholar] [CrossRef]
  26. Alimoradzadeh, M.; Tornabene, F.; Dimitri, R. Nonlinear axial-lateral coupled vibration of functionally graded-fiber reinforced composite laminated (FG-FRCL) beams subjected to aero-thermal loads. Int. J. Non-Linear Mech. 2024, 159, 104612. [Google Scholar] [CrossRef]
  27. Salmani, R.; Gholami, R.; Ansari, R.; Fakhraie, M. Analytical investigation on the nonlinear postbuckling of functionally graded porous cylindrical shells reinforced with graphene nanoplatelets. Eur. Phys. J. Plus 2021, 136, 53. [CrossRef]
  28. Li, Z.; Zhang, Q.; Shen, H.; Xiao, X.; Kuai, H.; Zheng, J. Buckling performance of the encased functionally graded porous composite liner with polyhedral shapes reinforced by graphene platelets under external pressure. Thin-Walled Struct. 2023, 183, 110370. [Google Scholar] [CrossRef]
  29. Shahani, A.R.; Kiarasi, F. Numerical and Experimental Investigation on Post-buckling Behavior of Stiffened Cylindrical Shells with Cutout subject to Uniform Axial Compression. J. Appl. Comput. Mech. 2023, 9, 25–44. [Google Scholar]
  30. Wu, Q.; Hu, S.; Tang, X.; Liu, X.; Chen, Z.; Xiong, J. Compressive buckling and post-buckling behaviors of J-type composite stiffened panel before and after impact load. Compos. Struct. 2023, 304, 116339. [Google Scholar] [CrossRef]
  31. Tran, K.L.; Douthe, C.; Sab, K.; Dallot, J.; Davaine, L. Buckling of stiffened curved panels under uniform axial compression. J. Constr. Steel Res. 2014, 103, 140–147. [Google Scholar] [CrossRef]
  32. Grondin, G.Y.; Elwi, A.E.; Cheng, J.J.R. Buckling of stiffened steel plates—A parametric study. J. Constr. Steel Res. 1999, 50, 151–175. [Google Scholar] [CrossRef]
  33. Reddy, J.N.; Starnes, J.H., Jr. General buckling of stiffened circular cylindrical shells according to a layerwise theory. Comput. Struct. 1993, 49, 605–616. [Google Scholar] [CrossRef]
  34. Ni, X.Y.; Prusty, B.G.; Hellier, A.K. Buckling and post-buckling of isotropic and composite stiffened panels: A review on optimisation (2000–2015). Int. J. Marit. Eng. 2016, 158, A-251–A-267. [Google Scholar] [CrossRef]
  35. Wodesenbet, E.; Kidane, S.; Pang, S.S. Optimization for buckling loads of grid stiffened composite panels. Compos. Struct. 2003, 60, 159–169. [Google Scholar] [CrossRef]
  36. Ni, X.Y.; Prusty, B.G.; Hellier, A.K. Buckling and post-buckling of isotropic and composite stiffened panels: A review on analysis and experiment (2000–2012). Int. J. Marit. Eng. 2015, 157. [Google Scholar] [CrossRef]
  37. Kidane, S.; Li, G.; Helms, J.; Pang, S.S.; Woldesenbet, E. Buckling load analysis of grid stiffened composite cylinders. Compos. Part B Eng. 2003, 34, 1–9. [Google Scholar] [CrossRef]
  38. Wang, D.; Abdalla, M.M.; Zhang, W. Buckling optimization design of curved stiffeners for grid-stiffened composite structures. Compos. Struct. 2017, 159, 656–666. [Google Scholar] [CrossRef]
  39. Patel, S.N.; Datta, P.K.; Sheikh, A.H. Buckling and dynamic instability analysis of stiffened shell panels. Thin-Walled Struct. 2006, 44, 321–333. [Google Scholar] [CrossRef]
  40. Ifayefunmi, O.; Ruan, D. Buckling of Stiffened Cone–Cylinder Structures Under Axial Compression. Int. J. Appl. Mech. 2022, 14, 2250075. [Google Scholar] [CrossRef]
  41. Zarei, M.; Rahimi, G.H.; Hemmatnezhad, M. On the buckling resistance of grid-stiffened composite conical shells under compression. Eng. Struct. 2021, 237, 112213. [Google Scholar] [CrossRef]
  42. Yoon, J.W.; Bray, G.H.; Valente RA, F.; Childs, T.E.R. Buckling analysis for an integrally stiffened panel structure with a friction stir weld. Thin-Walled Struct. 2009, 47, 1608–1622. [Google Scholar] [CrossRef]
  43. Duc, N.D.; Thang, P.T. Nonlinear buckling of imperfect eccentrically stiffened metal–ceramic–metal S-FGM thin circular cylindrical shells with temperature-dependent properties in thermal environments. Int. J. Mech. Sci. 2014, 81, 17–25. [Google Scholar] [CrossRef]
  44. Duc, N.D.; Thang, P.T. Nonlinear response of imperfect eccentrically stiffened ceramic–metal–ceramic FGM thin circular cylindrical shells surrounded on elastic foundations and subjected to axial compression. Compos. Struct. 2014, 110, 200–206. [Google Scholar] [CrossRef]
  45. Duc, N.D. Nonlinear dynamic response of imperfect eccentrically stiffened FGM double curved shallow shells on elastic foundation. Compos. Struct. 2013, 99, 88–96. [Google Scholar] [CrossRef]
  46. Duc, N.D.; Cong, P.H. Nonlinear thermal stability of eccentrically stiffened functionally graded truncated conical shells surrounded on elastic foundations. Eur. J. Mech.-A/Solids 2015, 50, 120–131. [Google Scholar] [CrossRef]
  47. Wang, Y.; Xie, K.; Fu, T.; Shi, C. Bending and elastic vibration of a novel functionally graded polymer nanocomposite beam reinforced by graphene nanoplatelets. Nanomaterials 2019, 9, 1690. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, D. Free vibration of functionally graded graphene platelets reinforced magnetic nanocomposite beams resting on elastic foundation. Nanomaterials 2020, 10, 2193. [Google Scholar] [CrossRef] [PubMed]
  49. Zhang, J.; Lv, Y.; Li, L. Dynamic Instability of Functionally Graded Graphene Platelet-Reinforced Porous Beams on an Elastic Foundation in a Thermal Environment. Nanomaterials 2022, 12, 4098. [Google Scholar] [CrossRef]
  50. Babaei, M.; Kiarasi, F.; Asemi, K.; Dimitri, R.; Tornabene, F. Transient Thermal Stresses in FG Porous Rotating Truncated Cones Reinforced by Graphene Platelets. Appl. Sci. 2022, 12, 3932. [Google Scholar] [CrossRef]
  51. Zhou, Z.; Wang, Y.; Zhang, S.; Dimitri, R.; Tornabene, F.; Asemi, K. Numerical Study on the Buckling Behavior of FG Porous Spherical Caps Reinforced by Graphene Platelets. Nanomaterials 2023, 13, 1205. [Google Scholar] [CrossRef] [PubMed]
  52. Dong, D.T.; Van Dung, D. A third-order shear deformation theory for nonlinear vibration analysis of stiffened functionally graded material sandwich doubly curved shallow shells with four material models. J. Sandw. Struct. Mater. 2019, 21, 1316–1356. [Google Scholar] [CrossRef]
  53. Blooriyan, S.; Ansari, R.; Darvizeh, A.; Gholami, R.; Rouhi, H. Postbuckling analysis of functionally graded graphene platelet-reinforced polymer composite cylindrical shells using an analytical solution approach. Appl. Math. Mech. 2019, 40, 1001–1016. [Google Scholar] [CrossRef]
Figure 1. Geometrical scheme of a stiffened FG-GPL multilayer shell panel with four different distributions of GPLs.
Figure 1. Geometrical scheme of a stiffened FG-GPL multilayer shell panel with four different distributions of GPLs.
Processes 12 00430 g001
Figure 2. First six buckling mode shapes of a stiffened cylindrical panel (SS boundary condition, Nr = Ns = 5, GPLX Δ = 0.01 wt. %).
Figure 2. First six buckling mode shapes of a stiffened cylindrical panel (SS boundary condition, Nr = Ns = 5, GPLX Δ = 0.01 wt. %).
Processes 12 00430 g002
Figure 3. First six buckling mode shapes of a stiffened cylindrical panel (CC boundary condition, Nr = Ns = 5, GPLX ∆ = 0.01 wt. %).
Figure 3. First six buckling mode shapes of a stiffened cylindrical panel (CC boundary condition, Nr = Ns = 5, GPLX ∆ = 0.01 wt. %).
Processes 12 00430 g003aProcesses 12 00430 g003b
Figure 4. First six buckling mode shapes of a stiffened cylindrical panel (SS boundary condition, Nr = Ns = 10, GPLX Δ = 0.01 wt. %).
Figure 4. First six buckling mode shapes of a stiffened cylindrical panel (SS boundary condition, Nr = Ns = 10, GPLX Δ = 0.01 wt. %).
Processes 12 00430 g004
Figure 5. First six buckling mode shapes of a stiffened cylindrical panel (CC boundary condition, Nr = Ns = 10, GPLX Δ = 0.01 wt. %).
Figure 5. First six buckling mode shapes of a stiffened cylindrical panel (CC boundary condition, Nr = Ns = 10, GPLX Δ = 0.01 wt. %).
Processes 12 00430 g005
Table 1. Buckling loads (MN/m) as predicted by our formulation and a commercial code.
Table 1. Buckling loads (MN/m) as predicted by our formulation and a commercial code.
Buckling Load λ 1 λ 2 λ 3 λ 4 λ 5 λ 6
Commercial code10.17110.1910.4010.45110.91011.261
Present10.17510.2010.4210.45810.97711.299
Table 2. Comparison of critical buckling load of a multilayer FG-GPL cylindrical shell for present study and Ref. [53] (R = 40h, L = 2R, h = 0.045 m).
Table 2. Comparison of critical buckling load of a multilayer FG-GPL cylindrical shell for present study and Ref. [53] (R = 40h, L = 2R, h = 0.045 m).
00.10.20.30.50.81
Present (UD)21.5339.1956.7574.38109.40161.75196.60
[53] (UD)21.5839.2256.8374.41109.48161.87196.65
Present (FG-X)21.5344.8067.5290.15135.25202.50247.21
[53] (FG-X)21.5844.8267.5990.23135.30202.59247.27
Present (FG-O)21.5331.9341.4050.3567.7092.95109.53
[53] (FG-O)21.5831.9941.4250.4267.7693.01109.56
Table 3. Convergence study of the first buckling load of structure (GPL-X, ∆ = 1 wt. %).
Table 3. Convergence study of the first buckling load of structure (GPL-X, ∆ = 1 wt. %).
Number of Ements (nθ × nx)10 × 1515 × 2520 × 3025 × 3530 × 45
Buckling load70.0266.1860.1255.2254.69
Table 4. Effect of GPL weight fractions on the buckling load (MN/m) of a SS structure (Ss = 0.378, Sr = 0.5).
Table 4. Effect of GPL weight fractions on the buckling load (MN/m) of a SS structure (Ss = 0.378, Sr = 0.5).
GPL Pattern Δ G P L % λ 1 λ 2 λ 3 λ 4 λ 5 λ 6
010.17510.2010.4210.45810.97711.299
GPL-X0.531.22231.27032.08632.09635.12935.176
154.69954.77656.10356.11959.42161.046
010.17510.2010.4210.45810.97711.299
GPL-A0.525.31825.34825.87625.88227.19427.999
141.11841.16842.02542.03544.16645.474
GPL-UD010.17510.2010.4210.45810.97711.299
0.527.05027.08327.64627.65329.05529.915
143.93143.98444.89944.91047.18748.584
010.17510.2010.4210.45810.97711.299
GPL-O0.524.34524.37424.88124.88726.14926.923
139.53739.58540.40940.41942.46843.725
Table 5. Effect of boundary conditions on the buckling load (MN/m) for SS and CC structures (Ss = 0.378, Sr = 0.5, Δ = 1 wt. %).
Table 5. Effect of boundary conditions on the buckling load (MN/m) for SS and CC structures (Ss = 0.378, Sr = 0.5, Δ = 1 wt. %).
GPL PatternBoundary Condition λ 1 λ 2 λ 3 λ 4 λ 5 λ 6
CC61.62561.69864.25764.26668.72071.201
GPL-XSS54.69954.77656.10356.11959.42161.046
CC47.12047.78951.74951.79853.94253.980
GPL-UDSS43.93143.98444.89944.91047.18748.584
CC45.10245.12146.20146.21148.07049.924
GPL-OSS39.53739.58540.40940.41942.46843.725
CC46.50546.52147.78247.78849.22551.025
GPL-ASS41.11841.16842.02542.03544.16645.474
Table 6. Impact of the number of rings and stringers on the buckling loads (MN/m) of the structure (GPL-UD, Δ = 1 wt. %).
Table 6. Impact of the number of rings and stringers on the buckling loads (MN/m) of the structure (GPL-UD, Δ = 1 wt. %).
Boundary
Condition
Number of Rings and Stringers λ 1 λ 2 λ 3 λ 4 λ 5 λ
Ss = 0.378, Sr = 0.5
or Nr = Ns = 5
47.12047.78951.74951.79853.94253.980
CCSs = 0.19, Sr = 0.25
or Nr = Ns = 10
52.95353.19457.37359.29660.92560.996
Ss = 0.378, Sr = 0.5
or Nr = Ns = 5
43.93143.98444.89944.91047.18748.584
SSSs = 0.19, Sr = 0.25
or Nr = Ns = 10
49.11249.79951.47352.30753.72553.908
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Z.; Tornabene, F.; Dimitri, R.; Babaei, M. Numerical Study of the Buckling Response of Stiffened FG Graphene-Reinforced Multilayer Composite Cylindrical Panels. Processes 2024, 12, 430. https://doi.org/10.3390/pr12030430

AMA Style

Liu Z, Tornabene F, Dimitri R, Babaei M. Numerical Study of the Buckling Response of Stiffened FG Graphene-Reinforced Multilayer Composite Cylindrical Panels. Processes. 2024; 12(3):430. https://doi.org/10.3390/pr12030430

Chicago/Turabian Style

Liu, Zhihong, Francesco Tornabene, Rossana Dimitri, and Masoud Babaei. 2024. "Numerical Study of the Buckling Response of Stiffened FG Graphene-Reinforced Multilayer Composite Cylindrical Panels" Processes 12, no. 3: 430. https://doi.org/10.3390/pr12030430

APA Style

Liu, Z., Tornabene, F., Dimitri, R., & Babaei, M. (2024). Numerical Study of the Buckling Response of Stiffened FG Graphene-Reinforced Multilayer Composite Cylindrical Panels. Processes, 12(3), 430. https://doi.org/10.3390/pr12030430

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