Next Article in Journal
Steel Fiber-Reinforced Concrete: A Systematic Review of the Research Progress and Knowledge Mapping
Previous Article in Journal
Dynamic Compressive Mechanical Properties of UR50 Ultra-Early-Strength Cement-Based Concrete Material under High Strain Rate on SHPB Test
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Free Vibration Analysis of a Graphene-Platelet-Reinforced, Porous, Two-Cylindrical-Panel System

1
Key Laboratory of Structural Dynamics of Liaoning Province, College of Sciences, Northeastern University, Shenyang 110819, China
2
Beijing Institute of Structure and Environment Engineering, Beijing 100076, China
*
Author to whom correspondence should be addressed.
Submission received: 17 August 2022 / Revised: 27 August 2022 / Accepted: 31 August 2022 / Published: 5 September 2022

Abstract

:
In this study, a novel, dynamic model of a graphene-platelet-reinforced, porous (GPLRP) double-cylindrical-panel system is proposed. The material properties of a graphene-platelet-reinforced, porous, double-cylindrical-panel system were determined by the Halpin–Tsai micromechanics model and the typical mechanical properties of open-cell metal foams. Different types of porosity distribution and graphene platelet (GPL) distribution patterns were considered. Love’s shell theory was utilized to derive the theoretical formulation, and the Rayleigh–Ritz method was used to calculate the natural frequencies of the system. The proposed model was validated by several comparison studies with the natural frequencies in the existing literature. Finally, the effects of stiffness of Winkler springs, boundary condition, porosity coefficient, porosity distribution, GPL distribution pattern, and GPL weight fraction on the free vibration characteristics of the system were evaluated.

1. Introduction

The two-cylindrical-panel system, an important component of modern lightweight structures, is widely used in mechanical, vehicle, and marine structures [1,2]. A two-cylindrical-panel system may possess ideal structural properties, such as better vibration attenuation and lighter weight, than a single cylindrical panel [3,4]. Therefore, the in-depth study of such systems has great theoretical and practical significance.
Recently, research on the vibration characteristics of multi-beam systems and multi-plate systems has been of great interest [5,6,7,8,9,10,11,12]. Many experts have conducted a large number of studies on the vibration characteristics of multi-beam systems [13,14,15,16,17,18]. Kim et al. [19,20] established a double-beam system and analyzed the free vibration characteristics of a system. Deng et al. [21] studied the vibration characteristics of a double-functionally-graded (FG) beam system on an elastic foundation. Hao et al. [22] utilized the modified Fourier–Ritz method to analyze the free vibration characteristics of a double-beam system. Rahman and Lee [23] developed a novel harmonic balance method to analyze the nonlinear vibration characteristics of a double-beam system. Kelly and Srinivas [24] established a system of multiple beams connected by elastic layers and investigated the free vibration of the system. Han et al. [25] discussed the vibration characteristics of a double-beam system by using the improved Wittrick–Williams algorithm.
Many researchers have also conducted numerous studies on the vibration characteristics of a multi-plate system. Oniszczuk [26] investigated the free vibration of two plates connected by a Winkler elastic layer. Based on the Rayleigh–Ritz method, Jeong and Kang [27] developed a novel theoretical method to study the vibration characteristics of multiple rectangular plates coupled by a liquid. Hedrih [28,29] investigated the free vibration characteristics of double-plate systems. Stojanović et al. [30] revealed the vibration and stability characteristics of a system of multiple rectangular plates connected with elastic layers by using an analytical method. The above studies on multi-plate structures are limited to rectangular plates, though many scholars have also carried out detailed studies on circular plates. Hedrih and Simonović [31,32,33] presented a dynamic model of the double circular plate system and investigated the non-linear characteristics of the system. Noga [34] established a dynamic model of double annular and circular membranes connected with the Winkler elastic layer and studied the free vibration of the system by using analytical and numerical methods.
Metal foam materials are widely employed in energy-absorbing systems due to their energy-absorbing capability [35] and other particular characteristics [36]. Keleshteri and Jelovica [37] investigated the free vibration characteristics of the FG porous cylindrical panels. Wang and Wu [38] analyzed the influences of the porous coefficient on the free vibration of an FG porous cylindrical shell.
Due to their pores, porous metal structures are light, but they also can reduce the stiffness of such a structure, so filling materials are needed to enhance the structural stiffness. GPLs possess good mechanical properties, which is why they are employed as reinforcement materials and widely used in various composite structures [39,40]. Using the GPLRP material, Yang et al. [41] analyzed the buckling and free vibration characteristics of a plate, and Wang et al. [42] studied the nonlinear free vibration characteristics of circular cylindrical shells. Chai and Wang [43] studied the traveling wave vibration characteristics of GPLRP-joined conical-cylindrical shells in a spinning motion. Ye and Wang [44] analyzed the nonlinear forced vibration of GPLRP cylindrical shells. Xu et al. [45] investigated the free vibration of a rotating GPLRP beam by using the differential transformation method. Teng and Wang [46] studied the primary, superharmonic, and subharmonic resonances of GPLRP plates. Zhou et al. [47] established an accurate nonlinear buckling analysis of the GPLRP composite cylindrical shells under an axial compressive load. Twinkle and Jeyaraj [48] analyzed the buckling and vibration characteristics of a GPLRP, cylindrical panel. The above studies indicate that a GPLRP structure has good vibration absorption performance and high strength, so it has good prospects for engineering structures.
The above review shows that most studies mainly focus on the vibration characteristics of plane systems, and there are almost no studies on curved systems. To fill the gap in the studies of curved systems, a novel dynamic model of a graphene-platelet-reinforced, porous (GPLRP) two-cylindrical-panel system based on Love’s shell theory is established. To validate the proposed model, several comparisons of the present results with those from open literature are made. Subsequently, the effects of stiffness of the Winkler springs, boundary condition, porosity coefficient, porosity distribution, GPL distribution pattern, and GPL weight fraction on the free vibration characteristics of the system are investigated using the Rayleigh–Ritz method.

2. Theoretical Formulations

2.1. Description of the Model

The schematic diagram of the GPLRP, two-cylindrical-panel system with the radius of curvature R, thickness h, subtended angle θ0, axial length L, and arc length S (S = 0) coupled by a polymer matrix is shown in Figure 1. As shown in Figure 1, the local cylindrical coordinate system (xi, θi, zi) of the ith (i = 1, 2) GPLRP, cylindrical panel is established on the midplane of the GPLRP, cylindrical panel, the origin oi of which is located at the midpoint on the left side of the midplane. The displacements of the arbitrary point of the midplane are denoted by ui, vi, and wi along with the directions xi, θi, and zi, respectively. The polymer matrix between the GPLRP double cylindrical panels can be considered equivalent to continuously distributed Winkler springs with stiffness K. Both of the GPLRP, cylindrical panels have identical material parameters, geometric parameters, and boundary conditions.
As shown in Figure 2, three types of porosity distribution along the thickness direction, denoted by porosity-I, porosity-II, and porosity-III, are considered in this study. For different porosity distributions, the Young’s modulus, mass density, and Poisson’s ratio along the thickness direction of GPLRP, cylindrical panel can be given as follows [41]:
E ( z ) = { E * [ 1 e 1 cos ( π z h ) ] Porosity - Ι E * { 1 e 2 [ 1 cos ( π z h ) ] } Porosity - Ι Ι E * e 3 Porosity - Ι Ι Ι
ρ ( z ) = { ρ * [ 1 e m 1 cos ( π z h ) ] Porosity - Ι ρ * { 1 e m 2 [ 1 cos ( π z h ) ] } Porosity - Ι Ι ρ * e m 3 Porosity - Ι Ι Ι
μ ( z ) = μ *
in which E*, ρ*, and μ* represent Young’s modulus, mass density, and Poisson’s ratio of the GPL-reinforced cylindrical panel without pores, respectively. e1, e2, and e3 are the porosity coefficients; em1, em2, and em3 are the mass density coefficients.
According to the Halpin–Tsai micromechanics model, the effective Young’s modulus of the GPL-reinforced cylindrical panel without porosities can be given as follows [41]:
E * = 3 8 [ 1 + ξ L GPL η L GPL V GPL ( z ) 1 η L GPL V GPL ( z ) ] E m + 5 8 [ 1 + ξ w GPL η w GPL V GPL ( z ) 1 η w GPL V GPL ( z ) ] E m
in which
ξ L GPL = 2 l GPL t GPL , ξ w GPL = 2 w GPL t GPL
η L GPL = E GPL E m E GPL + ξ L GPL E m , η w GPL = E GPL E m E GPL + ξ w GPL E m
where EGPL and Em represent Young’s moduli of the GPLs and metal matrix; lGPL, wGPL, and tGPL represent the average length, width, and thickness of the GPLs, respectively; VGPL(z) is the volume fraction of the GPLs.
The mass density ρ* and Poisson’s ratio μ* of the GPL reinforced cylindrical panel without porosities can be defined by [41]
ρ * = ρ m ( 1 V GPL ( z ) ) + ρ GPL V GPL ( z )
μ * = μ m ( 1 V GPL ( z ) ) + μ GPL V GPL ( z )
in which ρGPL and μGPL represent the mass density and Poisson’s ratio of the GPLs, respectively; ρm and μm denote the mass density and Poisson’s ratio of the metal matrix, respectively.
The relationship between the modulus of elasticity and mass density of the open-cell metal foams is as follows [41]:
E ( z ) E * = ( ρ ( z ) ρ * ) 2
The relationship between the mass density coefficients and the porosity coefficients can be obtained by substituting Equations (1) and (2) into Equation (9):
{ 1 e m 1 cos ( π z / h ) = 1 e 1 cos ( π z / h ) 1 e m 2 [ 1 cos ( π z / h ) ] = 1 e 2 [ 1 cos ( π z / h ) ] e m 3 = e 3
To obtain a valuable comparison of different combinations, the total mass of the GPLRP, cylindrical panel was set to be the same for various types of porosity and GPL distribution, resulting in
0 h / 2 1 e 1 cos ( π z / h ) d z = 0 h / 2 1 e 2 [ 1 cos ( π z / h ) ] d z = 0 h / 2 e 3 d z
in which e2 and e3 can be determined by e1.
To reinforce the porous cylindrical panel, the GPLs are applied as fillers inside the porous materials. The three GPL distribution patterns in the thickness direction of the cylindrical panel are shown in Figure 3. The volume fraction VGPL(z) of the GPLs can be written as follows [41]:
V GPL ( z ) = { s j 1 [ 1 cos ( π z / h ) ] GPL A s j 2 cos ( π z / h ) GPL B s j 3 GPL C
in which sj1, sj2, and sj3 are the coefficients of the GPL volume fraction for different porosity distributions and GPL distribution patterns, and j = 1, 2, 3 correspond to the Porosity-I, -II, and -III distributions, respectively. The sj1, sj2, and sj3 can be given as follows [41]:
V GPL T h / 2 h / 2 ρ ( z ) ρ * d z = { s j 1 h / 2 h / 2 [ 1 cos ( π z / h ) ] ρ ( z ) ρ * d z s j 2 h / 2 h / 2 cos ( π z / h ) ρ ( z ) ρ * d z s j 3 h / 2 h / 2 ρ ( z ) ρ * d z
where V GPL T is the total volume fraction of the GPLs and given by [41]
V GPL T = W GPL W GPL + ( ρ GPL / ρ m ) ( 1 W GPL )
where WGPL is the weight fraction obtained by the ratio of the mass of the GPL to the mass of the cylindrical panel.

2.2. Governing Equations and Solution

According to Love’s shell theory [49], the strains of the arbitrary point at the distance z from the midplane can be written as
{ ε x ( i ) ε θ ( i ) γ x θ ( i ) } = { ε 1 ( i ) ε 2 ( i ) γ 12 ( i ) } + z { φ 1 ( i ) φ 2 ( i ) 2 φ 12 ( i ) } , ( i = 1 , 2 )
where, ε 1 ( i ) , ε 2 ( i ) , and γ 12 ( i ) ; and φ 1 ( i ) , φ 2 ( i ) , and φ 12 ( i ) are midplane strains and surface curvatures of the ith GPLRP, cylindrical panel, which can be expressed by
ε 1 ( i ) = u i x i , ε 2 ( i ) = 1 R ( v i θ i + w i ) , γ 12 ( i ) = v i x i + 1 R u i θ i
φ 1 ( i ) = 2 w i x i 2 , φ 2 ( i ) = 1 R 2 ( 2 w i θ i 2 v i θ i ) , φ 12 ( i ) = 1 R ( 2 w i x i θ i v i x i )
The stress–strain relationship of the ith GPLRP, cylindrical panel can be defined as [50]
{ σ x ( i ) σ θ ( i ) τ x θ ( i ) } = [ Q 11 ( i ) Q 12 ( i ) 0 Q 12 ( i ) Q 22 ( i ) 0 0 0 Q 66 ( i ) ] { ε x ( i ) ε θ ( i ) γ x θ ( i ) }
where
Q 11 ( i ) = Q 22 ( i ) = E ( z ) 1 μ ( z ) μ ( z ) , Q 12 ( i ) = μ ( z ) E ( z ) 1 μ ( z ) μ ( z ) , Q 66 ( i ) = E ( z ) 2 ( 1 + μ ( z ) )
The force and moment resultants of the ith GPLRP, cylindrical panel can be obtained by integrating the stresses along with the z-direction of the GPLRP, cylindrical panel, which are expressed as
[ N x ( i ) N θ ( i ) N x θ ( i ) ] = h 2 h 2 [ σ x ( i ) σ θ ( i ) τ x θ ( i ) ] d z
[ M x ( i ) M θ ( i ) M x θ ( i ) ] = h 2 h 2 [ σ x ( i ) σ θ ( i ) τ x θ ( i ) ] z d z
By substituting Equations (15) and (18) into Equations (20) and (21), we get:
N ( i ) T = ε ( i ) T S ( i )
where ε ( i ) T and N ( i ) T are given by
ε ( i ) T = { ε 1 ( i ) , ε 2 ( i ) , γ 12 ( i ) , φ 1 ( i ) , φ 2 ( i ) , 2 φ 12 ( i ) } N ( i ) T = { N x ( i ) , N θ ( i ) , N x θ ( i ) , M x ( i ) , M θ ( i ) , M x θ ( i ) }
and S ( i ) is a stiffness matrix defined by
S ( i ) = [ A 11 ( i ) A 12 ( i ) 0 B 11 ( i ) B 12 ( i ) 0 A 12 ( i ) A 22 ( i ) 0 B 12 ( i ) B 22 ( i ) 0 0 0 A 66 ( i ) 0 0 B 66 ( i ) B 11 ( i ) B 12 ( i ) 0 D 11 ( i ) D 12 ( i ) 0 B 12 ( i ) B 22 ( i ) 0 D 12 ( i ) D 22 ( i ) 0 0 0 B 66 ( i ) 0 0 D 66 ( i ) ]
in which A p q ( i ) , B p q ( i ) , and D p q ( i ) (p, q = 1, 2, 6) are stretching, coupling, and bending stiffness coefficients and can be expressed as
[ A p q ( i ) B p q ( i ) D p q ( i ) ] = h 2 h 2 Q p q ( i ) [ 1 z z 2 ]   d z
The strain energy and kinetic energy of the GPLRP, two-cylindrical-panel system can be written as
U = 1 2 i = 1 2 0 L θ 0 / 2 θ 0 / 2 ε ( i ) T S ( i ) ε ( i ) R d θ i d x i
T = 1 2 i = 1 2 0 L θ 0 / 2 θ 0 / 2 I ( i ) [ ( u i t ) 2 + ( v i t ) 2 + ( w i t ) 2 ] R d θ i d x i
where
I ( i ) = h 2 h 2 ρ ( i ) ( z ) d z
The potential energy stored in the polymer matrix between the two-cylindrical-panel system can be written as
U S = 1 2 0 L θ 0 / 2 θ 0 / 2 K ( w 1 w 2 ) 2 R d θ i d x i
where K is the stiffness of Winkler springs.
For generality and convenience, the dimensionless parameters x ¯ i and θ ¯ i are introduced
x ¯ i = 2 x i L 1 , θ ¯ i = 2 θ i θ 0
Substituting Equation (30) into Equations (16) and (17) yields
ε 1 ( i ) = 2 L u i x ¯ i , ε 2 ( i ) = 1 R ( 2 θ 0 v i θ ¯ i + w i ) , γ 12 ( i ) = 2 L v i x ¯ i + 1 R 2 θ 0 u i θ ¯ i
φ 1 ( i ) = 4 L 2 2 w i x ¯ i 2 ,     φ 2 ( i ) = 1 R 2 ( 4 θ 0 2 2 w i θ ¯ i 2 2 θ 0 v i θ ¯ i ) ,   φ 12 ( i ) = 1 R ( 4 L θ 0 2 w i x ¯ i θ ¯ i 2 L v i x ¯ i )
By separating variables, the admissible displacement functions of the ith GPLRP, cylindrical panel can be given by
u i ( x i , θ i , t ) = U i ( x ¯ i , θ ¯ i ) sin ( ω t ) v i ( x i , θ i , t ) = V i ( x ¯ i , θ ¯ i ) sin ( ω t ) w i ( x i , θ i , t ) = W i ( x ¯ i , θ ¯ i ) sin ( ω t )
where t is the time variable; the ω is the natural angular frequency; U i ( x ¯ i , θ ¯ i ) ; V i ( x ¯ i , θ ¯ i )       and W i ( x ¯ i , θ ¯ i ) are the displacement amplitude functions of the ith GPLRP, cylindrical panel.
The displacement amplitude functions of the ith GPLRP, cylindrical panel expressed by the Chebyshev polynomials can be given as [41]
U i ( x ¯ i , θ ¯ i ) = f u ( x ¯ i , θ ¯ i ) m = 1 M n = 1 N U m n ( i ) P m ( x ¯ i ) P n ( θ ¯ i ) V i ( x ¯ i , θ ¯ i )       = f v ( x ¯ i , θ ¯ i ) m = 1 M n = 1 N V m n ( i ) P m ( x ¯ i ) P n ( θ ¯ i ) W i ( x ¯ i , θ ¯ i ) = f w ( x ¯ i , θ ¯ i ) m = 1 M n = 1 N W m n ( i ) P m ( x ¯ i ) P n ( θ ¯ i )
where U m n ( i ) , V m n ( i ) , and W m n ( i ) are the unknown coefficients; M and N are the truncation numbers of Chebyshev polynomials. P m ( x ¯ i ) and P n ( θ ¯ i ) are the mth and nth Chebyshev polynomials of the first kind, which can be written as [51]
P m ( x ¯ i ) = cos [ ( m 1 ) arccos ( x ¯ i ) ] , P n ( θ ¯ i ) = cos [ ( n 1 ) arccos ( θ ¯ i ) ]
where P m ( x ¯ i ) and P n ( θ ¯ i ) are a set of orthogonal polynomials in the interval [−1,1].
In addition, in Equation (34), functions f α ( x ¯ i , θ ¯ i ) , α = u, v, w are the auxiliary functions related to boundary conditions.
f α ( x ¯ i , θ ¯ i ) = ( 1 + x ¯ i ) p ( 1 x ¯ i ) q ( 1 + θ ¯ ) r ( 1 θ ¯ ) s
in which p, q, r, and s are equal to zero, one, or two, for different boundary conditions. For instance, in the GPLRP, two-cylindrical-panel system with clamped boundaries (C) on x = 0 and x = L, free (F) at θ = −θ0/2, and simply supported (S) at θ = θ0/2, auxiliary functions can be written as
f u ( x ¯ i , θ ¯ i ) = ( 1 + x ¯ i ) 1 ( 1 x ¯ i ) 1 ( 1 + θ ¯ i ) 0 ( 1 θ ¯ i ) 1 f v ( x ¯ i , θ ¯ i ) = ( 1 + x ¯ i ) 1 ( 1 x ¯ i ) 1 ( 1 + θ ¯ i ) 0 ( 1 θ ¯ i ) 0 f w ( x ¯ i , θ ¯ i ) = ( 1 + x ¯ i ) 2 ( 1 x ¯ i ) 2 ( 1 + θ ¯ i ) 0 ( 1 θ ¯ i ) 1
The strain energy of the GPLRP, two-cylindrical-panel system can be obtained by substituting Equations (24), (31)–(33) into (26):
U = i = 1 2 1 1 1 1 [ sin ( ω   t ) ] 2 { A 22 ( i ) L θ 0 W i 2 8 R + A 66 ( i ) L 2 R θ 0 ( U i θ ¯ i ) 2 + B 22 ( i ) L W i 2 R 2 V i θ ¯ i + A 22 ( i ) L W i 2 R V i θ ¯ i + D 22 ( i ) L 2 R 3 θ 0 ( V i θ ¯ i ) 2 + B 22 ( i ) L R 2 θ 0 ( V i θ ¯ i ) 2 + A 22 ( i ) L 2 R θ 0 ( V i θ ¯ i ) 2 B 22 ( i ) L W i R 2 θ 0 2 W i θ ¯ i 2 2 D 22 ( i ) L R 3 θ 0 2 V i θ ¯ i 2 W i θ ¯ i 2 2 B 22 ( i ) L R 2 θ 0 2 V i θ ¯ i 2 W i θ ¯ i 2 + 2 D 22 ( i ) L R 3 θ 0 3 ( 2 W i θ ¯ i 2 ) 2 + A 12 ( i ) θ 0 W i 2 U i x ¯ i + A 12 ( i ) V i θ ¯ i U i x ¯ i + B 12 ( i ) R V i θ ¯ i U i x ¯ i 2 B 12 ( i ) R θ 0 2 W i θ ¯ i 2 U i x ¯ i + A 11 ( i ) R θ 0 2 L ( U i x ¯ i ) 2 + A 66 ( i ) U i θ ¯ i V i x ¯ i + 2 B 66 ( i ) R U i θ ¯ i V i x ¯ i + 2 B 66 ( i ) θ 0 L ( V i x ¯ i ) 2 + 2 D 66 ( i ) θ 0 L R ( V i x ¯ i ) 2 + A 66 ( i ) R θ 0 2 L ( V i x ¯ i ) 2 4 B 66 ( i ) R θ 0 U i θ ¯ i 2 W i x ¯ i θ ¯ i 4 B 66 ( i ) L V i x ¯ i 2 W i x ¯ i θ ¯ i 8 D 66 ( i ) L R V i x ¯ i 2 W i x ¯ i θ ¯ i + 8 D 66 ( i ) L R θ 0 ( 2 W i x ¯ i θ ¯ i ) 2 B 12 ( i ) θ 0 W i L 2 W i x ¯ i 2 2 B 12 ( i ) L V i θ ¯ i 2 W i x ¯ i 2 2 D 12 ( i ) L R V i θ ¯ i 2 W i x ¯ i 2 + 4 D 12 ( i ) L R θ 0 2 W i θ ¯ i 2 2 W i x ¯ i 2 2 B 11 ( i ) R θ 0 L 2 U i x ¯ i 2 W i x ¯ i 2 + 2 D 11 ( i ) R θ 0 L 3 ( 2 W i x ¯ i 2 ) 2 } d θ ¯ i d x ¯ i
The kinetic energy of the GPLRP, two-cylindrical-panel system can be obtained by substituting Equation (33) into Equation (27):
T = i = 1 2 1 1 1 1 1 8 L R θ 0 ω 2 [ cos ( ω t ) ] 2 [ I ( i ) ( U i 2 + V i 2 + W i 2 ) ] d θ ¯ i d x ¯ i
The potential energy of the polymer matrix can be acquired by substituting Equation (33) into Equation (29):
U S = 1 1 1 1 L R K θ 0 8 [ sin ( ω t ) ] 2 ( W 1 W 2 ) 2 d θ ¯ i d x ¯ i
To determine the solutions, the Rayleigh–Ritz method is employed. The Lagrangian energy function Π of the GPLRP, two-cylindrical-panel system can be expressed as
Π = U max + U S max T max
in which Umax and Tmax are the maximum strain energy and maximum kinetic energy of the GPLRP, two-cylindrical-panel system; USmax is the maximum potential energy of the polymer matrix.
Then, the Lagrangian energy function Π is minimized with respect to the unknown coefficient, yielding
U m n ( i ) = 0 , V m n ( i ) = 0 , W m n ( i ) = 0 , ( i = 1 , 2 ; m = 1 , , M ; n = 1 , , N )
The eigenvalue equations derived from Equation (42) can be summarized in the matrix form as
( K ω 2 M ) X = 0
in which M and K denote the mass and stiffness matrices, respectively. The explicit forms of the modal vector X can be written as follows:
X = { X ( 1 ) , X ( 2 ) } T
in which
X ( i ) = { U 11 ( i ) , , U M N ( i ) , V 11 ( i ) , , V M N ( i ) , W 11 ( i ) , , W M N ( i ) } , ( i = 1 , 2 )
For convenience, the dimensionless natural frequency of the GPLRP, two-cylindrical-panel system is defined by
ω ¯ = ω R ρ m / E m

3. Model Validation and Convergence Analysis

To prove the validity of the present formulations, the dimensionless fundamental frequency ω ¯ = ω L ρ m ( 1 μ m 2 ) / E m of a single GPLRP plate under simply supported boundary conditions is compared with the results in the existing literature [41]. The material parameters used are consistent with the literature [41]. The geometric parameters of a single GPLRP cylindrical panel are as follows: h = 0.1 m, L/S = 1, L/h = 20, WGPL = 1.0%, e1 = 0.5, i = 1, K = 0, and R = ∞. As shown in Table 1, it can be clearly seen that the results of this paper are well consistent with the literature [41].
To further clarify the correctness of the derivation of this study, Table 2 shows the comparison of the natural angular frequencies ω (rad/s) in this article with existing literature for the double-plate system under the simply supported boundary condition [26]. The geometric and material parameters are as follows: L = 1 m, S = 2 m, E = 1 × 1010 N/m2, h = 1 × 10−2 m, i = 1, 2, K = 6 × 104 N/m3, μ = 0.3, ρ = 5 × 103 kg/m3, and R = ∞. As can be observed in Table 2, very good agreement was obtained in this comparison study. The above comparisons with the open literature are sufficient to prove the accuracy of the derivation in this paper.
In this study, the Rayleigh–Ritz method was applied to solve the natural frequencies of the GPLRP, two-cylindrical-panel system. The accuracy of this method depends on the truncation numbers (M, N) of the Chebyshev polynomials. Therefore, it is necessary to discuss the influence of truncation orders M and N on the calculation results. The simulation parameters are given in Section 4. As shown in Table 3, when the truncation numbers M = N = 11, the dimensionless natural frequencies ω ¯ = ω R ρ m / E m obtained good convergence. Therefore, the truncation numbers M and N were set to 11 in the following study.

4. Results and Discussion

In this section, the material and geometric parameters of the GPLRP, two-cylindrical-panel system in this paper used are as follows: Em = 130 GPa, ρm = 8960 kg/m3, μm = 0.34, EGPL = 1.01 Tpa, ρGPL = 1062.5 kg/m3, μGPL = 0.186, tGPL = 1.5 nm, lGPL = 2.5 μm, wGPL = 1.5 μm, WGPL = 1.0%, e1 = 0.5, R = 1 m, h = 0.01R, L = R, θ0 = π/6, and K = 1 × 1014 N/m3.
Figure 4 shows the variation un dimensionless natural frequencies of the system with the spring stiffness K, considering different boundary conditions. An interesting phenomenon is that the variation of spring stiffness has no influence on the first-order natural frequency of the system. In addition, there is an elastic domain for other natural frequencies except for the first-order natural frequency, in which the natural frequency increases rapidly with the increase in the spring stiffness K. This finding is further confirmed by the results in Table 4. The reason for this phenomenon is that the natural frequency of the GPLRP, two-cylindrical-panel system can be divided into two parts, one of which is the same as the natural frequency of a single GPLRP, cylindrical panel, and the other is identical to the natural frequency of a single GPLRP, cylindrical panel on an elastic foundation with the stiffness 2K [26]. Therefore, except for the first-order natural frequency, other natural frequencies can be adjusted by adjusting the value of the spring stiffness to obtain the desired natural frequency. This means that the stiffness of the system can be changed by adjusting the spring stiffness. For the sake of brevity, only the fundamental frequency of the system is discussed below.
Figure 5 and Figure 6 demonstrate the effects of the porosity coefficient on dimensionless fundamental frequencies of the system for specific porosity distributions and boundary conditions. It can be clearly seen that the natural frequency of the system decreases significantly with the increase in the porosity coefficient, except for Porosity-I under the CCCC boundary condition. This indicates that in most cases, the greater the porosity coefficient e1, the weaker the stiffness of the system. Then, comparing Figure 5a and Figure 6a, it can also be seen that the natural frequency of the system under the SSSS boundary condition decreases rapidly with the increase in the porosity coefficient. However, the natural frequency of the system under the CCCC boundary condition changes slightly with the increase in the porosity coefficient, which indicates that the combined influence of the porosity coefficient and GPL distribution on the natural frequency is related to the boundary conditions. It is also notable that when e1 ≤ 0.5, for all porosity distributions, GPL A possesses the largest natural frequency and GPL B possesses the smallest natural frequency. However, when e1 > 0.5, GPL A does not always possess the largest natural frequency; for example, see Porosity-II under the SSSS boundary condition. It can be concluded that the effect of the porosity coefficient on the natural frequency of the system is associated with the porosity distributions, GPL distribution patterns, and boundary conditions.
To explain the effect of porosity distribution more clearly, Figure 7 and Figure 8 indicate the dimensionless fundamental frequencies versus the porosity coefficient for different porosity distributions with specific GPL distribution patterns and boundary conditions. As can be seen in Figure 7 and Figure 8, for all kinds of porosity distributions, Porosity-I possesses the largest natural frequency. This indicates that Porosity-I can significantly enhance the stiffness of the system. It is worth noting that the natural frequency of the system decreases significantly with the increase in the porosity coefficient, except for the Porosity-I under the CCCC boundary condition. This further confirms that the effect of the porosity distributions on the natural frequency of the system is associated with the boundary conditions. From the above results, it can be confirmed that the porosity coefficient, porosity distribution, GPL distribution pattern, and the boundary condition interact with each other and have coupled effect on the free vibration characteristics of the GPLRP, two-cylindrical-panel system.
Figure 9 reveals the variations in dimensionless fundamental frequencies with GPL weight fraction WGPL for different porosity distributions and GPL distribution patterns. It is clearly shown that the natural frequency increases significantly with the increase in the weight fraction WGPL for different porosity distributions and GPL distribution patterns. This means that the GPL filling material can significantly improve the effective stiffness of the system. As shown in Figure 9b, it was interesting to find that all GPL distribution patterns in Porosity-II have almost the same increase in efficiency with the increase in the GPL weight fraction WGPL. This indicates that the porosity distributions play a dominant role, and it also shows that choosing an appropriate porosity distribution is of great significance to improving the stiffness of the system. This means that the strength of the system can be improved by choosing the appropriate porosity distribution.

5. Conclusions

This study proposed a dynamic model of the GPLRP, two-cylindrical-panel system with general boundary conditions based on Love’s shell theory to study the free vibration characteristics of the system. The Rayleigh–Ritz method is used to calculate the solutions. From the results, some valuable conclusions can be drawn as follows:
(1)
The variation in the stiffness of the Winkler springs in the GPLRP, two-cylindrical-panel system has no effect on the first-order natural frequency, but other natural frequencies can be adjusted by controlling the value of the spring stiffness to obtain the desired natural frequency.
(2)
The porosity coefficient, porosity distribution, GPL distribution pattern, and boundary condition interact with each other and have a coupled influence on the vibration characteristics of the GPLRP, two-cylindrical-panel system.
(3)
Increasing the proportion of the GPL filling material can significantly improve the stiffness of the GPLRP, two-cylindrical-panel system, thereby increasing the natural frequency of the system.

Author Contributions

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

Funding

This research was supported by the National Natural Science Foundation of China (grant number 12002080).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhai, Y.; Chai, M.; Su, J.; Liang, S. Dynamics properties of composite sandwich open circular cylindrical shells. Compos. Struct. 2018, 189, 148–159. [Google Scholar] [CrossRef]
  2. Mao, Q.; Pietrzko, S. Experimental study for control of sound transmission through double glazed window using optimally tuned Helmholtz resonators. Appl. Acoust. 2010, 71, 32–38. [Google Scholar] [CrossRef]
  3. Oniszczuk, Z. Forced transverse vibrations of an elastically connected complex rectangular simply supported double-plate system. J. Sound Vib. 2004, 270, 997–1011. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Shi, D. An exact Fourier series method for vibration analysis of elastically connected laminated composite double-beam system with elastic constraints. Mech. Adv. Mater. Struc. 2021, 28, 2440–2457. [Google Scholar] [CrossRef]
  5. Stojanović, V.; Kozić, P. Vibrations and Stability of Complex Beam Systems; Springer: New York, NY, USA, 2015. [Google Scholar]
  6. Mao, Q.; Wattanasakulpong, N. Vibration and stability of a double-beam system interconnected by an elastic foundation under conservative and nonconservative axial forces. Int. J. Mech. Sci. 2015, 93, 1–7. [Google Scholar] [CrossRef]
  7. Oniszczuk, Z. Free transverse vibrations of elastically connected simply supported double-beam complex system. J. Sound Vib. 2000, 232, 387–403. [Google Scholar] [CrossRef]
  8. Oniszczuk, Z. Transverse vibrations of elastically connected double-string complex system, Part I: Free vibrations. J. Sound Vib. 2000, 232, 355–366. [Google Scholar] [CrossRef]
  9. Sobhy, M.; Zenkour, A.M. Porosity and inhomogeneity effects on the buckling and vibration of double-FGM nanoplates via a quasi-3D refined theory. Compos. Struct. 2019, 220, 289–303. [Google Scholar] [CrossRef]
  10. Jamalpoor, A.; Hosseini, M. Biaxial buckling analysis of double-orthotropic microplate-systems including in-plane magnetic field based on strain gradient theory. Compos. Part. B-Eng. 2015, 75, 53–64. [Google Scholar] [CrossRef]
  11. Xing, W.C.; Wang, Y.Q. Vibration characteristics analysis of rigid-flexible spacecraft with double-direction hinged solar arrays. Acta. Astronaut. 2022, 193, 454–468. [Google Scholar] [CrossRef]
  12. Xing, W.C.; Wang, Y.Q. Vibration characteristics of thin plate system joined by hinges in double directions. Thin-Walled Struct. 2022, 175, 109260. [Google Scholar] [CrossRef]
  13. Wu, Y.; Gao, Y. Dynamic response of a simply supported viscously damped double-beam system under the moving oscillator. J. Sound Vib. 2016, 384, 194–209. [Google Scholar] [CrossRef]
  14. Stojanović, V.; Kozić, P. Forced transverse vibration of Rayleigh and Timoshenko double-beam system with effect of compressive axial load. Int. J. Mech. Sci. 2012, 60, 59–71. [Google Scholar] [CrossRef]
  15. Han, F.; Dan, D.; Cheng, W. Exact dynamic characteristic analysis of a double-beam system interconnected by a viscoelastic layer. Compos. Part B-Eng. 2019, 163, 272–281. [Google Scholar] [CrossRef]
  16. Brito, W.; Maia, C.; Mendonca, A. Bending analysis of elastically connected Euler–Bernoulli double-beam system using the direct boundary element method. Appl. Math. Model. 2019, 74, 387–408. [Google Scholar] [CrossRef]
  17. Vu, H.; Ordonez, A.; Karnopp, B. Vibration of a double-beam system. J. Sound Vib. 2000, 229, 807–822. [Google Scholar] [CrossRef]
  18. Zhang, Y.; Lu, Y.; Ma, G. Effect of compressive axial load on forced transverse vibrations of a double-beam system. Int. J. Mech. Sci. 2008, 50, 299–305. [Google Scholar] [CrossRef]
  19. Kim, G.; Han, P.; An, K.; Choe, D.; Ri, Y.; Ri, H. Free vibration analysis of functionally graded double-beam system using Haar wavelet discretization method. Eng. Sci. Technol. 2021, 24, 414–427. [Google Scholar] [CrossRef]
  20. Kim, K.; Han, P.; Jong, K.; Jang, C.; Kim, R. Natural frequency calculation of elastically connected double-beam system with arbitrary boundary condition. AIP. Adv. 2020, 10, 055026. [Google Scholar] [CrossRef]
  21. Deng, H.; Cheng, W.; Zhao, S. Vibration and buckling analysis of double-functionally graded Timoshenko beam system on Winkler-Pasternak elastic foundation. Compos. Struct. 2017, 160, 152–168. [Google Scholar] [CrossRef]
  22. Hao, Q.; Zhai, W.; Chen, Z. Free vibration of connected double-beam system with general boundary conditions by a modified Fourier–Ritz method. Arch. Appl. Mech. 2018, 88, 741–754. [Google Scholar] [CrossRef]
  23. Rahman, M.S.; Lee, Y.-Y. New modified multi-level residue harmonic balance method for solving nonlinearly vibrating double-beam problem. J. Sound Vib. 2017, 406, 295–327. [Google Scholar] [CrossRef]
  24. Kelly, S.G.; Srinivas, S. Free vibrations of elastically connected stretched beams. J. Sound Vib. 2009, 326, 883–893. [Google Scholar] [CrossRef]
  25. Han, F.; Dan, D.; Cheng, W. An exact solution for dynamic analysis of a complex double-beam system. Compos. Struct. 2018, 193, 295–305. [Google Scholar] [CrossRef]
  26. Oniszczuk, Z. Free transverse vibrations of an elastically connected rectangular simply supported double-plate complex system. J. Sound Vib. 2000, 236, 595–608. [Google Scholar] [CrossRef]
  27. Jeong, K.-H.; Kang, H.-S. Free vibration of multiple rectangular plates coupled with a liquid. Int. J. Mech. Sci. 2013, 74, 161–172. [Google Scholar] [CrossRef]
  28. Hedrih, K. Double plate system with a discontinuity in the elastic bonding layer. Acta Mech. Sin. 2007, 23, 221–229. [Google Scholar] [CrossRef]
  29. Hedrih, K.S. Transversal vibrations of double-plate systems. Acta Mech. Sin. 2006, 22, 487–501. [Google Scholar] [CrossRef]
  30. Stojanović, V.; Kozić, P.; Ristić, M. Vibrations and stability analysis of multiple rectangular plates coupled with elastic layers based on different plate theories. Int. J. Mech. Sci. 2015, 92, 233–244. [Google Scholar] [CrossRef]
  31. Hedrih, K.R.S.; Simonović, J.D. Non-linear dynamics of the sandwich double circular plate system. Int. J. Nonlin. Mech. 2010, 45, 902–918. [Google Scholar] [CrossRef]
  32. Simonović, J.D. Multi-frequency analysis of the double circular plate system non-linear dynamics. Nonlinear Dyn. 2012, 67, 2299–2315. [Google Scholar]
  33. Hedrih, K.R.; Simonović, J.D. Energies of the dynamics in a double circular plate nonlinear system. Int. J. Bifurc. Chaos 2011, 21, 2993–3011. [Google Scholar] [CrossRef]
  34. Noga, S. Free transverse vibration analysis of an elastically connected annular and circular double-membrane compound system. J. Sound Vib. 2010, 329, 1507–1522. [Google Scholar] [CrossRef]
  35. Avalle, M.; Belingardi, G.; Montanini, R. Characterization of polymeric structural foams under compressive impact loading by means of energy-absorption diagram. Int. J. Impact. Eng. 2001, 25, 455–472. [Google Scholar] [CrossRef]
  36. Smith, B.; Szyniszewski, S.; Hajjar, J.; Schafer, B.; Arwade, S. Steel foam for structures: A review of applications, manufacturing and material properties. J. Constr. Steel Res. 2012, 71, 1–10. [Google Scholar] [CrossRef]
  37. Keleshteri, M.M.; Jelovica, J. Nonlinear vibration behavior of functionally graded porous cylindrical panels. Compos. Struct. 2020, 239, 112028. [Google Scholar] [CrossRef]
  38. Wang, Y.; Wu, D. Free vibration of functionally graded porous cylindrical shell using a sinusoidal shear deformation theory. Aerosp. Sci. Technol. 2017, 66, 83–91. [Google Scholar] [CrossRef]
  39. Sobhy, M. Magneto-electro-thermal bending of FG-graphene reinforced polymer doubly-curved shallow shells with piezoelectromagnetic faces. Compos. Struct. 2018, 203, 844–860. [Google Scholar] [CrossRef]
  40. 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]
  41. 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]
  42. Wang, Y.Q.; Ye, C.; Zu, J.W. Nonlinear vibration of metal foam cylindrical shells reinforced with graphene platelets. Aerosp. Sci. Technol. 2019, 85, 359–370. [Google Scholar] [CrossRef]
  43. Chai, Q.; Wang, Y.Q. Traveling wave vibration of graphene platelet reinforced porous joined conical-cylindrical shells in a spinning motion. Eng. Struct. 2022, 252, 113718. [Google Scholar] [CrossRef]
  44. Ye, C.; Wang, Y.Q. Nonlinear forced vibration of functionally graded graphene platelet-reinforced metal foam cylindrical shells: Internal resonances. Nonlinear Dyn. 2021, 104, 2051–2069. [Google Scholar] [CrossRef]
  45. Xu, H.; Wang, Y.Q.; Zhang, Y. Free vibration of functionally graded graphene platelet-reinforced porous beams with spinning movement via differential transformation method. Arch. Appl. Mech. 2021, 91, 4817–4834. [Google Scholar] [CrossRef]
  46. Teng, M.W.; Wang, Y.Q. Nonlinear forced vibration of simply supported functionally graded porous nanocomposite thin plates reinforced with graphene platelets. Thin-Walled Struct. 2021, 164, 107799. [Google Scholar] [CrossRef]
  47. 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]
  48. 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]
  49. Soedel, W. Vibrations of Shells and Plates; CRC Press: New York, NY, USA, 2004. [Google Scholar]
  50. Amabili, M. Nonlinear Vibrations and Stability of Shells and Plates; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  51. Zhang, W.; Niu, Y.; Behdinan, K. Vibration characteristics of rotating pretwisted composite tapered blade with graphene coating layers. Aerosp. Sci. Technol. 2020, 98, 105644. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the GPLRP, two-cylindrical-panel system. (a) two-cylindrical-panel system; (b) single-cylindrical-panel.
Figure 1. Schematic diagram of the GPLRP, two-cylindrical-panel system. (a) two-cylindrical-panel system; (b) single-cylindrical-panel.
Materials 15 06158 g001
Figure 2. Cross-section of a single GPLRP, cylindrical panel for different porosity distributions. (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Figure 2. Cross-section of a single GPLRP, cylindrical panel for different porosity distributions. (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Materials 15 06158 g002
Figure 3. Three types of GPL distribution patterns. (a) GPL A; (b) GPL B; (c) GPL C.
Figure 3. Three types of GPL distribution patterns. (a) GPL A; (b) GPL B; (c) GPL C.
Materials 15 06158 g003
Figure 4. The influences of spring stiffness on first six dimensionless natural frequencies of the GPLRP, two-cylindrical-panel system (Porosity-I, GPL A, e1 = 0.5). (a) SSSS; (b) CCCC.
Figure 4. The influences of spring stiffness on first six dimensionless natural frequencies of the GPLRP, two-cylindrical-panel system (Porosity-I, GPL A, e1 = 0.5). (a) SSSS; (b) CCCC.
Materials 15 06158 g004
Figure 5. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (SSSS). (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Figure 5. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (SSSS). (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Materials 15 06158 g005
Figure 6. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (CCCC). (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Figure 6. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (CCCC). (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Materials 15 06158 g006
Figure 7. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (SSSS). (a) GPL A; (b) GPL B; (c) GPL C.
Figure 7. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (SSSS). (a) GPL A; (b) GPL B; (c) GPL C.
Materials 15 06158 g007
Figure 8. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (CCCC). (a) GPL A; (b) GPL B; (c) GPL C.
Figure 8. The influences of porosity coefficient on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (CCCC). (a) GPL A; (b) GPL B; (c) GPL C.
Materials 15 06158 g008
Figure 9. The influences of GPL weight fraction on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (SSSS). (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Figure 9. The influences of GPL weight fraction on dimensionless fundamental frequencies of the GPLRP, two-cylindrical-panel system (SSSS). (a) Porosity-I; (b) Porosity-II; (c) Porosity-III.
Materials 15 06158 g009
Table 1. Comparison of the dimensionless fundamental frequency of a single GPLRP plate.
Table 1. Comparison of the dimensionless fundamental frequency of a single GPLRP plate.
Porosity DistributionsGPL DistributionsPresentYang et al. [41]
GPL A0.40380.3958
Porosity-IGPL B0.33330.3293
GPL C0.36330.3574
GPL A0.30890.3072
Porosity-IIGPL B0.26110.2601
GPL C0.27660.2754
GPL A0.36860.3627
Porosity-IIIGPL B0.30430.3014
GPL C0.32940.3252
Table 2. Comparison of the natural angular frequencies of the double-plate system.
Table 2. Comparison of the natural angular frequencies of the double-plate system.
ModesPresentOniszczuk [26]
152.852.8
272.072.0
384.584.5
497.797.7
5137.3137.3
6145.8145.8
Table 3. Convergence of the first six dimensionless natural frequencies of the GPLRP, two-cylindrical-panel system (Porosity-I, GPL A, e1 = 0.5).
Table 3. Convergence of the first six dimensionless natural frequencies of the GPLRP, two-cylindrical-panel system (Porosity-I, GPL A, e1 = 0.5).
Mode(M, N)
(7, 7)(8, 8)(9, 9)(10, 10)(11, 11)(12, 12)(13, 13)
10.316120.316120.316120.316120.316120.316120.31612
20.666710.666710.666670.666670.666670.666670.66667
30.693300.693300.693300.693300.693300.693300.69330
SSSS40.830130.830130.830090.830090.830090.830090.83009
50.997560.997510.997460.997460.997460.997460.99746
61.102231.101341.101261.101251.101251.101251.10125
11.041441.041431.041431.041431.041431.041431.04143
21.076921.076911.076901.076901.076901.076901.07690
CCCC31.144131.144111.144111.144111.144111.144111.14411
41.235741.235711.235701.235701.235701.235701.23570
51.286141.286091.286071.286071.286071.286071.28607
61.505741.505721.505691.505691.505681.505681.50568
Table 4. The influences of spring stiffness on dimensionless natural frequencies of the GPLRP, two-cylindrical-panel system (Porosity-I, GPL A, e1 = 0.5).
Table 4. The influences of spring stiffness on dimensionless natural frequencies of the GPLRP, two-cylindrical-panel system (Porosity-I, GPL A, e1 = 0.5).
ModeSpring Stiffness K
102104106108101010121014
SSSS10.316120.316120.316120.316120.316120.316120.31612
20.316120.316150.319200.543420.666670.666670.66667
30.666670.666670.666670.666670.693300.693300.69330
40.666670.666680.668160.693300.830090.830090.83009
50.693300.693300.693300.802720.997460.997460.99746
60.693300.693310.694720.823581.101251.101251.10125
CCCC11.041431.041431.041431.041431.041431.041431.04143
21.041431.041441.042381.076901.076901.076901.07690
31.076901.076901.076901.132671.144111.144111.14411
41.076901.076911.077841.144111.235701.235701.23570
51.144111.144111.144111.166431.286071.286071.28607
61.144111.144121.144981.228761.505681.505681.50568
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, X.; Chi, W.; Luo, J. Free Vibration Analysis of a Graphene-Platelet-Reinforced, Porous, Two-Cylindrical-Panel System. Materials 2022, 15, 6158. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15176158

AMA Style

Sun X, Chi W, Luo J. Free Vibration Analysis of a Graphene-Platelet-Reinforced, Porous, Two-Cylindrical-Panel System. Materials. 2022; 15(17):6158. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15176158

Chicago/Turabian Style

Sun, Xianguang, Weichao Chi, and Jia Luo. 2022. "Free Vibration Analysis of a Graphene-Platelet-Reinforced, Porous, Two-Cylindrical-Panel System" Materials 15, no. 17: 6158. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15176158

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