Next Article in Journal
Fabrication and Characterization of a Metallic–Dielectric Nanorod Array by Nanosphere Lithography for Plasmonic Sensing Application
Next Article in Special Issue
Effects of Graphene Oxide-Gold Nanoparticles Nanocomposite on Highly Sensitive Foot-and-Mouth Disease Virus Detection
Previous Article in Journal
Facile Synthesis of Antimony Tungstate Nanosheets as Anodes for Lithium-Ion Batteries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bending and Elastic Vibration of a Novel Functionally Graded Polymer Nanocomposite Beam Reinforced by Graphene Nanoplatelets

1
Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Beijing Key Laboratory of CO2 Utilization and Reduction Technology, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, China
2
Institute of Systems Engineering, China Academy of Engineering Physics, Mianyang 621900, China
3
Beijing Key Lab of MFPTS, China Academy of Safety Science & Technology, Beijing 100029, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2019, 9(12), 1690; https://doi.org/10.3390/nano9121690
Submission received: 22 October 2019 / Revised: 19 November 2019 / Accepted: 22 November 2019 / Published: 26 November 2019
(This article belongs to the Special Issue Multifunctional Polymer-Based Nanocomposite Materials)

Abstract

:
A novel functionally graded (FG) polymer-based nanocomposite reinforced by graphene nanoplatelets is proposed based on a new distribution law, which is constructed by the error function and contains a gradient index. The variation of the gradient index can result in a continuous variation of the weight fraction of graphene nanoplatelets (GPLs), which forms a sandwich structure with graded mechanical properties. The modified Halpin–Tsai micromechanics model is used to evaluate the effective Young’s modulus of the novel functionally graded graphene nanoplatelets reinforced composites (FG-GPLRCs). The bending and elastic vibration behaviors of the novel nanocomposite beams are investigated. An improved third order shear deformation theory (TSDT), which is proven to have a higher accuracy, is implemented to derive the governing equations related to the bending and vibrations. The Chebyshev–Ritz method is applied to describe various boundary conditions of the beams. The bending displacement, stress state, and vibration frequency of the proposed FG polymer-based nanocomposite beams under uniformly distributed loads are provided in detail. The numerical results show that the proposed distributions of GPL nanofillers can lead to a more effective pattern of improving the mechanical properties of GPL-reinforced composites than the common ones.

1. Introduction

High-strength and multifunctional polymer-based nanocomposites with low-content carbon-based nanofillers have attracted widespread attentions due to their remarkable mechanical properties such as large deformation, stretchability, good chemical compatibility, and energy storage capacity [1,2,3,4,5,6,7,8,9,10]. That is because the carbon-based nanofillers, such as carbon nanotubes (CNTs), graphene nanoplatelets (GPLs), and graphene oxides (GOs), can provide unique nanostructures and superior properties in the nanocomposites. The functionally graded material (FGM) is a novel type of composite material whose mechanical properties smoothly and continuously vary in a preferred direction [11,12,13,14]. To better use the superior mechanical properties of carbon-based nanofillers and inspire from the concept of FGM, the functionally graded carbon nanotube-reinforced composite (FG-CNTRC) [15,16,17] and functionally graded graphene nanoplatelet reinforced composite (FG-GPLRC) [18,19,20] have been introduced, where the weight fractions of the CNTs and GPLs vary in the thickness direction.
The dispersion of nanofillers in the polymer matrix based on a reasonable graded distribution can effectively make use of reinforcements. To address the effects of nanofiller distributions on the mechanical behaviors of FG polymer-based nanocomposites, different types of distributions, such as the uniform distribution (UD), FG-V shape, FG-O shape and FG-X shape, were introduced and employed in many reports, e.g., references [21,22,23,24,25,26]. In addition, the distribution laws in forms of general polynomials have also been implemented [27,28,29]. All existing functions to describe the nanofiller distribution law are not adjustable since no adjustable parameter is included. Until now, the distribution law having adjustable parameters, which can lead to continuously graded mechanical properties for the nanofiller-reinforced polymer nanocomposites, has not been reported. It is believed that an adjustable distribution has great potential for introducing a novel type of nanocomposite, and also can be used to optimize the mechanical performances of the nanofiller-reinforced structures. The first contribution of the current work is to propose a new distribution law with an adjustable parameter as the gradient index.
In recent years, the graphene nanoplatelet (GPL) has been considered as an ideal reinforcement because of its excellent material properties, including high stiffness and strength. A polymer with a low concentration of GPLs exhibits distinct improvement in mechanical properties. Yang and his coauthors devoted great efforts to the bending, buckling, and vibration behaviors of FG-GPLRC structures [30,31,32]. Shen et al. [33,34] discussed various results of the nonlinear bending, vibration, and buckling behaviors of FG-GPLRC plates, panels, and shells. Mao and Zhang [35] examined the linear and nonlinear vibration responses of graphene-reinforced piezoelectric composite plates under external voltage excitation. Wang et al. [36] studied the vibration and bending behaviors of functionally graded nanocomposite doubly curved shallow shells reinforced by GPLs. Gholami et al. [37] investigated the nonlinear, harmonically excited vibration of a third-order shear deformable FG-GPLRC rectangular plate. Although the studies concerning on the mechanical behaviors of FG-GPLRC structures are massive, few studies focus on improving the distribution of the GPLs. Furthermore, the bending and vibration performances are the most basic and critical properties for structural components. Lezgy-Nazargah and Salahshuran [38] proposed a novel mixed-field theory with relatively low number of unknown variables for investigating bending and vibration analysis of multi-layered composite plates. Wattanasakulpong and Ungbhakorn [39] studied the bending, buckling and vibration behaviors of carbon-nanotube-reinforced composite (CNTRC) beams. Ghannadpour et al. [40] investigated the bending, buckling and vibration analyses of nonlocal Euler beams. Zhao et al. [41] analyzed the bending and vibrations of functionally graded trapezoidal nanocomposite plates reinforced with GPLs. The second contribution of our work is to introduce a novel FG-GPLRC beam using the proposed distribution law and investigate its bending and vibration behaviors.
In present work, a novel functionally graded graphene nanoplatelet reinforced composite (FG-GPLRC) is proposed by developing a new distribution law. The variation of the gradient index results in a continuous variation of the weight fraction of graphene nanoplatelets (GPLs) and forms graded mechanical properties within composite according to the thickness direction. The modified Halpin–Tsai micromechanics model is employed to evaluate the effective Young’s modulus. The bending and elastic vibration behaviors of the novel FG-GPLRC beams are investigated. An improved third-order shear deformation theory (TSDT), which is proved to have a higher accuracy, is implemented to derive the governing equations. Chebyshev–Ritz method is applied to describe various boundary conditions of the beam. Parametric studies are performed and present the bending displacement, stress state, and vibration frequency of the new FG-GPLRC beam under the uniformly distributed loads. The effects of gradient index of the proposed GPL distributions on the bending and vibration performance are addressed.

2. A Novel GPL-Reinforced Nanocomposite

2.1. Evaluation of Effective Mechanical Properties

It is well known that the Eshelby-Mori-Tanaka theory and the Halpin–Tsai micromechanics model are the most common methods to determine the effective elastic moduli of nanocomposite materials. Much literature suggest that the Eshelby-Mori-Tanaka theory is able to yield a more accurate prediction of the nanocomposite mechanical properties [42,43,44,45,46], and this conclusion was supported by experiments for CNT-reinforced nanocomposite.
The Halpin–Tsai model is based on the generalized self-consistent method. Although this model was initially developed for traditional fiber composites, successful efforts were undertaken to apply this model to nanofillers-reinforced nanocomposites by introducing orientation and shape factors that account for the geometry of the filler phase. The Halpin–Tsai model was used to calculate the effective Young’s modulus of a polymer nanocomposite reinforced by graphene nanoplatelets (GPLs), and the comparisons between the theoretical predictions and experimental results were also performed [47,48]. Due to the simple form in mathematics, the Halpin–Tsai micromechanics model has been widely employed to estimate the effective Young’s modulus of functionally graded graphene nanoplatelets reinforced composites [27,28,29,30,31,32,33,34,35,36,37]. The main objective of the current work is to propose an adjustable distribution law to find a more effective way to use the GPL reinforcements. It is believed that the effective material properties resulted from the Halpin–Tsai model can provide the results with some reference values. Thus, the Halpin–Tsai micromechanics model is selected.
A beam composed of FG-GPLRCs is considered herein. The length, width, and thickness of the beam are L, b, and h, respectively. The origin of the coordinate system is fixed at the center at the left end of the beam. The matrix of the nanocomposite is a polymer, and the GPL nanofillers are uniformly or non-uniformly dispersed across the thickness direction of the beam. A modified Halpin–Tsai model is used to calculate the effective Young’s modulus of the GPL/polymer composite. We assume that GPLs are effective rectangular solid fillers dispersed in a polymer matrix, and the effective Young’s modulus EC of the GPL/polymer composite is approximated by Voigt–Reuss model as follows [18]:
E C = 3 8 E L + 5 8 E T
where longitudinal modulus EL and transverse modulus ET are determined by the Halpin–Tsai model as follows:
E L = 1 + ξ L η L V G P L 1 η L V G P L E M
E T = 1 + ξ W η W V G P L 1 η W V G P L E M
Substituting Equations (2) and (3) into Equation (1), we obtain the following expression:
E C = 3 8 1 + ξ L η L V G P L 1 η L V G P L E M + 5 8 1 + ξ W η W V G P L 1 η W V G P L E M
where
η L = ( E G P L / E M ) 1 ( E G P L / E M ) + ξ L
η W = ( E G P L / E M ) 1 ( E G P L / E M ) + ξ W
where EM and EGPL are Young’s moduli of the polymer matrix and GPLs, respectively; VGPL is the volume fraction of the graphene nanoplatelets; ξL and ξW characterize the geometry and size of the GPL nanofillers, respectively, and are defined as follows:
ξ L = 2 ( l G P L h G P L )
ξ W = 2 ( w G P L h G P L )
where lGPL, wGPL, and hGPL are the average length, width, and thickness of GPLs, respectively. Mass density ρC and Poisson’s ratio vC of the GPL/polymer nanocomposite are calculated by the rule of mixtures as follows:
ρ C = V G P L ρ G P L + V M ρ M
v C = V G P L v G P L + V M v M
where VM is the volume fraction of the polymer matrix; subscripts “GPL”, “M”, and “C” denote the GPLs, polymer matrix, and GPL/polymer nanocomposite, respectively. The volume fraction of GPLs is:
V G P L = g G P L g G P L + ( ρ G P L / ρ M ) ( 1 g G P L )
where gGPL is the total weight fraction of the GPLs in the nanocomposite.

2.2. A New GPL Distribution

Previous works proved that different reinforcement distributions significantly affect the mechanical properties of the composites. FG-X, FG-O and FG-V are the most common types of distributions for the GPLs or CNTs, where the total weight fractions or volume fractions of the GPLs or CNTs are constant regardless how the reinforcements dispersed in the polymer matrix.
The motivation of the present work is to find a GPL distribution with an adjustable parameter that controls the weight fraction of GPLs in a preferred direction and maintain a constant total weight fraction. This will be helpful to address the impacts of the GPL distributions on the structural responses and provide a new distribution for designing polymer nanocomposites and optimizing the mechanical properties of FG-GPLRC composites in different applications.
The new distribution function is constructed by using an error function:
{ g G P L ( z ) = g G P L 0 2 r π Erf ( r ) e 4 r 2 ( z h 1 2 ) 2 0 z h 2 g G P L ( z ) = g G P L 0 2 r π Erf ( r ) e 4 r 2 ( z h + 1 2 ) 2 h 2 z 0
where Erf ( r ) = 2 π 0 r e t 2 d t is the error function; r is the gradient index and takes only positive values (r > 0); g G P L 0 is the total weight fraction of the GPLs. Although the value of r cannot be zero, when the selected r approaches 0, the value of Equation (8) approaches 1 as its limit, which represents the uniform GPL distribution pattern. From Equation (8), the following features of the proposed GPL distribution are obtained:
  • The total weight fraction of GPLs g G P L 0 remains constant with the r variations.
  • The GPLs symmetrically disperse in the matrix about the mid-plane of the beam.
  • When r increases, many GPLs are increasingly dispersed to the upper and lower surfaces of the beam.
  • This GPL graded distribution is continuous along the thickness and makes the mechanical properties of the FG-GPLRCs continuously vary.
To clearly illustrate the mentioned characteristics for the proposed distribution law, Figure 1 presents the variation of GPL weight fraction with respect to r. The total weight fraction of the GPLs in Figure 1 is assumed to be 0.5%.
For further illustrating, Figure 2 presents the schematic of the GPL distributions of the proposed model with different values of the gradient index. The features of the new FG-GPLRCs are clearly observed.
In fact, a functionally graded GPL-reinforced composite (FG-GPLRC) structure is ideal in combining the advantages of both FGMs and GPLs. However, the fabrication of such functionally graded structures with a continuous and smooth variation of GPLs across the thickness is extremely difficult due to the constraint of manufacture technology. For overcoming this problem, the functionally graded GPL reinforced multilayer structures are introduced, because a functionally graded GPL reinforced multilayer nanocomposite structure in which each individual layer is made from a mixture of uniformly distributed GPL reinforcements and polymer matrix with GPL concentration incrementally varying layer by layer is much easier to fabricate. The published results have provided authentic evidence that such a multilayer structure is an excellent approximation to the ideal functionally graded structure with a continuous and smooth variation of GPLs across the thickness direction when the total number of layers is sufficiently large [18,31,32,49,50,51].
In the following, the dimensions of the graphene nanoplatelets are lGPL = 2.5 µm, wGPL = 1.5 µm, and hGPL = 1.5 nm. The material properties of GPLs and epoxy are ρGPL = 1.06 g/cm3, EGPL = 1.01 TPa, ρM = 1.2 g/cm3, EM = 3.0 GPa, vM = 0.34, and vGPL = 0.186.

3. Theory and Formulations

In the present work, an accurate and efficient modeling is established based on an improved third-order shear deformation theory, and the governing equations associated with the static and vibration problems of the beams with various boundary conditions are derived using the minimum energy methodology and Chebyshev–Ritz method.

3.1. An Improved Third Order Shear Deformation Theory

The improved third-order shear deformation theory (TSDT) was originally proposed by Shi [52] based on rigorous kinematic of displacements. The published literature has proven that the results obtained by the improved TSDT are more reliable and accurate than many other theories because the kinematics of displacements is derived from elasticity theory instead of the hypothesis of displacement like other existing approaches.
The displacement fields of the improved third-order shear deformation theory (TSDT) are expressed as follows [52,53]:
{ u = u 0 ( x , t ) + 5 4 ( z 4 3 h 2 z 3 ) ϕ x ( x , t ) + ( 1 4 z 5 3 h 2 z 3 ) w 0 ( x , t ) x w = w 0 ( x , t )
where u0 and w0 define the generalized displacements at the mid-plane of the beam in the x and z directions, and φx is the rotation of the beam. From the above displacement fields, the small normal strain εxx and transverse shear strain γxz can be written as
{ ε x x = u x = u 0 x + z 4 ( 5 ϕ x x + 2 w 0 x 2 ) 5 z 3 3 h 2 ( ϕ x x + 2 w 0 x 2 ) γ x z = u z + w x = 5 4 ( ϕ x + w 0 x ) 5 z 2 h 2 ( ϕ x + w 0 x )
Then, the stresses are obtained under the assumption of Hook’s law:
{ σ x x = Q 11 ( z ) ε x x σ x z = Q 55 ( z ) ε x z
where Q11(z) and Q55(z) are the elastic constants that continuously vary along the beam thickness and expressed as
{ Q 11 ( z ) = E ( z ) 1 ν 2 Q 55 ( z ) = E ( z ) 2 ( 1 + ν )
The strain energy U of the beam is
U = 1 2 0 L h / 2 h / 2 b ( σ x x ε x x + σ x z ε x z ) d z d x
Substituting Equations (9)–(12) into Equation (13), we obtain the strain energy expression as a function of the material stiffness and strain components as
U = 1 2 0 L { A 11 ( u 0 x ) 2 + B 11 2 ( 5 u 0 x ϕ x x + u 0 x 2 w 0 x 2 ) + D 11 16 ( 5 ϕ x x + 2 w 0 x 2 ) 2 10 E 11 3 h 2 ( u 0 x ϕ x x + u 0 x 2 w 0 x 2 ) 5 F 11 6 h 2 [ 5 ( ϕ x x ) 2 + 6 ϕ x x 2 w 0 x 2 + ( 2 w 0 x 2 ) 2 ] + 25 9 h 4 H 11 ( ϕ x x + 2 w 0 x 2 ) 2 + 25 16 A 55 ( ϕ x + w 0 x ) 2 D 55 25 2 h 2 ( ϕ x + w 0 x ) 2 + 25 h 4 F 55 ( ϕ x + w 0 x ) 2 } d x
where A11, B11, D11, F11, H11, A55, D55 and F55 are the stiffness constants, and they are defined as
{ ( A 11 , B 11 , D 11 , E 11 , F 11 , H 11 ) = b h / 2 h / 2 Q 11 ( z ) ( 1 , z , z 2 , z 3 , z 4 , z 6 ) d z ( A 55 , D 55 , F 55 ) = b h / 2 h / 2 Q 55 ( z ) ( 1 , z 2 , z 4 ) d z
The kinetic energy of the composite beam is
T = 1 2 0 L h / 2 h / 2 b { ρ ( z ) [ ( u t ) 2 + ( w t ) 2 ] } d z d x
where ρ(z) is the mass density of the beam, which varies in the thickness direction.
Substituting Equation (9) into Equation (16), we can rewrite the kinetic energy as
T = 1 2 0 L { I 0 [ ( u 0 t ) 2 + ( w 0 t ) 2 ] + I 1 2 ( 5 u 0 t ϕ x t + u 0 t 2 w 0 x t ) + I 2 16 [ 25 ( ϕ x t ) 2 + 10 ϕ x t 2 w 0 x t + ( 2 w 0 x t ) 2 ] 10 I 3 3 h 2 ( u 0 t ϕ x t + u 0 t 2 w 0 x t ) 5 I 4 6 h 2 [ 5 ( ϕ x t ) 2 + 6 ϕ x t 2 w 0 x t + ( 2 w 0 x t ) 2 ] + 25 I 6 9 h 4 [ ( ϕ x t ) 2 + 2 ϕ x t 2 w 0 x t + ( 2 w 0 x t ) 2 ] } d x
where I i = b h / 2 h / 2 ρ ( z ) z i , i = 0 , 1 , 2 , 3 , 4 , 6 is the inertia terms.
The work done by a uniform distributed load q is:
V = 0 L q w ( x ) d x
The total energy function (Π) of the FG-GPLRC beam for the bending problem can be written as:
Π = U + V
The total energy functional (Π) of FG-GPLRC beams for the free-vibration analysis is expressed as:
Π = U T

3.2. Chebyshev–Ritz Method

Ritz method is known as an effective tool to analyze the structural behavior. Since the adoption of trial functions only depends on the essential type of boundary condition [54], various functions may be selected as the admissible functions. In the present work, each displacement function can be written in the form of triplicate series of Chebyshev polynomials multiplied by a boundary function, which ensures that the displacement component satisfies the essential geometric boundary condition of the beam, i.e.,
{ u 0 ( x , t ) = B u ( x ) i = 1 N U ¯ i ( t ) P i ( x ) w 0 ( x , t ) = B w ( x ) i = 1 N W ¯ i ( t ) P i ( x ) ϕ x ( x , t ) = B φ ( x ) i = 1 N V ¯ i ( t ) P i ( x )
where N is the truncated number of Chebyshev polynomials. BΞ (x) (Ξ = u, w and ϕ) are the boundary functions; ( U ¯ n , W ¯ n , V ¯ n ) denote the unknown coefficients corresponding to time and are expressed as follows: ( U ¯ n , W ¯ n , V ¯ n ) = ( U n , W n , V n ) e i ω n t , i = 1 , where ωn is the vibration frequency.
Pi(x) in Equation (21) is the i-th Chebyshev polynomial of the first kind, which is commonly known as “the most optimal expansion” [55,56], and defined in the interval of [−1, 1] as:
P i ( x ) = cos ( ( i 1 ) arccos ( 2 x L 1 ) ) , i = 1 , 2 , 3
The recursive relationship is
{ P 0 ( χ ) = 1 P 1 ( χ ) = χ P i + 1 ( χ ) = 2 χ P i ( χ ) P i 1 ( χ )
There are two distinct advantages of selecting Chebyshev polynomial series as the admissible functions for each displacement component [57,58]: (1) Pi(x) is a set of complete and orthogonal series in the interval of [−1, 1] and has more rapid convergence and better numerical stability in computation than other polynomials; (2) Pi(x) can be expressed in a simple and unified form of cosine function as shown in Equation (22), which reduces the coding efforts.
The boundary functions BΞ(x) (Ξ = u, w and ϕ) that correspond to u, w and ϕ are provided in the following uniform formula:
B Θ ( x ) = ( x L ) L Θ ( 1 x L ) R Θ
where LΘ and RΘ are indices from the following essential geometric boundary conditions:
(1)
Hinged-Hinged (H-H)
x = 0: u0 = 0; w = 0; φx ≠ 0
x = L: u0 = 0; w = 0; φx ≠ 0
(2)
Clamped-Clamped (C-C)
x = 0: u0 = 0; w = 0; d w d x = 0 ; φx ≠ 0
x = L: u0 = 0; w = 0; d w d x = 0 ; φx ≠ 0
(3)
Clamped-Hinged (C-H)
x = 0: u0 = 0; w = 0; d w d x = 0 ; φx = 0
x = L: u0 = 0; w = 0; φx = 0
Table 1 lists the values of the indices for various boundary conditions.
Substituting Equation (21) into the total energy functional (Π) for the bending in Equation (19) and free vibration in Equation (20) and taking the derivative with respect to the unknown coefficients in the procedure of finding minimization requires
Π U n = 0 , Π W n = 0 , Π V n = 0
The aforementioned procedure produces a system of simultaneous equations with an equal number of unknown coefficients (Un, Wn, Vn).
The equation system of solving the static bending problem of the FG-GPLRC beam under a distributed load q can be given as:
[ K ] Δ ¯ = F
where [K] is the structural stiffness matrix; F is a column vector associated with the external load from Equation (17); and Δ ¯ is the solution column vector for the static problem.
The generalized eigenvalue problem for free vibration is expressed as follows:
( [ K ] ω 2 [ M ] ) Δ = 0
where [K] is the structural stiffness matrix, [M] is the mass matrix, and ω is the natural frequency. Vector Δ is the eigenvector from the displacement functions, which represents the modal shapes of the structures. The dimensions of the aforementioned matrices are 3 N × 3 N. The dimension of Δ is 3 N × 1.

4. Convergence and Validation Studies

The validation and accuracy of the current model developed based on the improved third-order shear deformation theory and Chebyshev–Ritz method are verified by comparing the transverse deflections, stresses and vibration frequencies of the FG and FG sandwich beams reported in the existing literature.
A hinged-hinged (H-H) sandwich beam with power-law type FG face sheets and a homogeneous hardcore, which is subjected to a uniform distributed load, is selected for the first example. The power-law index of the facesheets is 0.5. The top and bottom surfaces of the beam are metal-rich, while its core is made of ceramic. The material properties of the metal (Al) are E = 70 GPa and ρ = 2702 kg/m3, and those of ceramic (Al2O3) are E = 380 GPa and ρ = 3960 kg/m3. Both Poisson’s ratios are 0.3. The length-to-height ratio of the sandwich beam is set to be 5 (L/h = 5), and its layer thickness ratio is 2:1:2.
The transverse deflection (w) at mid-span, normal stresses (σxx) at the top surface of the mid-span, transverse shear stress (σxz) at mid-plane at the left end, and natural frequencies of the sandwich beam are calculated and compared with the results from [59] in the following dimensionless form:
w ¯ = w h 3 12 384 E m 5 q 0 L 4 , σ ¯ x x = σ x x ( L 2 , z ) h q 0 L , σ ¯ x z = σ x z ( 0 , z ) h q 0 L , ω i = Ω i L 2 h ρ m E m
Table 2 tabulates the dimensionless results, which show excellent consistency between the present results and the published results based on a high-order shear deformation theory. Meanwhile, 6 is a reasonable number of Chebyshev polynomial terms, which is sufficiently large to obtain accurate results.
In addition, Table 3 and Table 4 list the dimensionless mid-span displacements and fundamental frequencies of the sandwich beams with the layer thickness ratio of 1:2:1. The H-H and C-C boundary conditions are used, and various power-law indices (p) are considered. The results are compared with the available ones obtained by the first-order shear deformation theory (FSDT), third-order shear deformation theory (TSDT) and quasi-3D beam theory from the literature [60,61]. The present results well match the published results. It is important to note that the numerical results from the present model are closer to the results from the quasi-3D theory than the other beam theories, which indicates that the present model is accurate and efficient.

5. Results and Discussion

To facilitate the presentation, the following dimensionless parameters are introduced as follows:
Vertical   displacement :   w ¯ ( x ) = w ( x ) 100 E M h 3 q L 4
Dimensionless   normal   stress :   σ ¯ x x = σ x x ( L 2 , z ) h q L
Dimensionless   shear   stress :   σ ¯ x z = σ x z ( 0 , z ) h q L
Dimensionless   frequencies :   Ω i = ω i L 2 h ρ M E M
where ρM and EM are the density and elastic modulus of the polymer matrix, respectively.

5.1. Bending Analysis

5.1.1. Bending Deflection

In this section, the effects of the gradient indices and GPL weight fractions on the bending behaviors of the new FG-GPLRC beams are investigated. Table 5, Table 6 and Table 7 present the numerical results of the maximum dimensionless displacement w ¯ max of the FG-GPLRC beams with various boundary conditions.
Figure 3 shows the relationship between maximum dimensionless displacement w ¯ max and GPL weight fractions. There is a nonlinear variation between the values of w ¯ max with respect to the GPL weight fractions. Regardless of the gradient index, w ¯ max decreases when the total weight fraction of the GPLs g G P L 0 increases. The curves in Figure 3 clearly indicate that the decrease in w ¯ max increases when the GPL content is relatively low, while the decrease tends to be slight when the GPL content increases to a relatively high level, which suggests that the bending response of the beam is more sensitive to the variations of the GPL contents when g G P L 0 is relatively low.
Another concerning point is the effect of the gradient index on the bending displacements of the beam. Figure 4 presents the relationship between the maximum dimensionless displacements w ¯ max and the values of gradient index r. L/h = 5 and L/h = 20 are considered. The curves in Figure 4 indicate that parameter r has remarkable effects on the bending performances of the FG-GPLRC beam because the variation of r changes the variation of GPLs in the thickness direction of the beam when the weight fraction is fixed and affects the beam stiffness.
For the thick beam (L/h = 5), as shown in Figure 4a, the regulation between r and w ¯ max becomes complex and interesting. When the total weight fraction g G P L 0 is very small, e.g., 0.25%, the increase in r will result in a decrease of w ¯ max implies that dispersing more reinforcements near the surfaces can improve the bending resistance of the beam. However, when the total weight fraction g G P L 0 is relatively high, e.g., 1.25%, when gradient index r increases, w ¯ max first reduces and subsequently increases, which suggests that there is a gradient index value that can minimize w ¯ max . This r value can produce the best GPL distribution that provides the greatest improvement to the FG-GPLRC beam.
For the thin beam (L/h = 20), as shown in Figure 4b, the maximum dimensionless displacements w ¯ max decrease when gradient index r increases regardless of the geometric parameters and boundary conditions. As mentioned, when the total GPL weight fractions remain constant, the increase in r implies that more GPLs are dispersed in the top and bottom surfaces of the beam. Therefore, Figure 4b indicates that dispersing more GPL nanofillers near the top and bottom surfaces of the beam is most effective in enhancing the beam stiffness. The numerical results in Table 5, Table 6 and Table 7 also support this conclusion.
In addition, for the thin beam, Figure 4b demonstrates that there is a nonlinear descent in the maximum dimensionless displacements w ¯ max versus the increase in gradient index r. For comparison, the maximum dimensionless displacement w ¯ max of the polymer composite beam of dispersing the GPLs in a FG-X type and UD type, are also displayed in Figure 4b. Figure 4b shows an intersection point of the two curves, which implies that the proposed FG-GPLRC beam and conventional FG-X composite beams have identical bending deflections. The corresponding value of r for this point will be of significance because when r is beyond this point, the new FG-GPLRC beams can provide better deformed-resistant capability than the conventional FG-X ones. In addition, we find that the value in Figure 4b is equal to 1.6 (r = 1.6), and this r at the intersection point is constant regardless of the boundary conditions, geometric parameters and total GPL weight fractions. This result can be mathematically explained that when r = 1.6, the values of the function from Equation (8) are closer to those from the FG-X function ( g ( z ) = 4 | z h | ) at a given z-position, i.e., the GPL distribution along the thickness direction is similar to the GPL distribution from the FG-X one when r = 1.6.

5.1.2. Stress State

A hinged-hinged FG-GPLRC beam with a length-to-height ratio of 5 is used to evaluate the normal and shear stresses in the thickness direction, which reveals the effects of the gradient index and GPL weight fraction on the stress distributions of the beam.
The effect of the gradient index on the normal stress distributions is displayed in Figure 5. Figure 5 shows that the normal stress linearly varies with the thickness coordinates when r approaches 0 because in this situation, the GPLs are uniformly dispersed in the polymer matrix, and the mechanical properties of the FG-GPLRCs have no change in the thickness direction. When r increases, nonlinear change occurs and become more remarkable because a greater r indicates that more GPLs are concentrated on the top and bottom surfaces, which makes a significantly uneven GPL distribution and produces the graded mechanical properties of the FG-GPLRCs. In addition, the increasingly concentrated dispersion of GPLs to the surfaces increases the normal stresses near the top and bottom of the beam when the gradient index increases but reduces those in the middle part of the beam. As mentioned, when r increases, the new FG-GPLRC beam can become a sandwich structure, and the stress distributions in Figure 5 prove the presence of the smooth and continuous variation of stresses throughout the entire cross-section.
Meanwhile, the effects of the total weight fractions of the GPLs on the normal stress distributions are presented in Figure 5. The cases of r = 1 and r = 4 are considered in Figure 6a,b, respectively. The variations of normal stresses are linear when g G P L 0 = 0 for a pure polymer beam, but they nonlinearly change when g G P L 0 > 0 because the uneven distribution of GPLs produces the gradient variation of material properties. Further addition of GPLs increases the normal stresses near the top and bottom surfaces of the beam and reduces those in the internal part far away from the two surfaces.
Figure 7 shows the variation of the shear stress distributions in the thickness direction of the FG-GPLRC beams with various gradient index values. To satisfy the traction-free condition on the surfaces of the improved TSDT, the shear stresses on the top and bottom surfaces are zero, and the maximum shear stress appears near the mid-surface of the beam. The curves in Figure 7 show that the change in gradient index significantly affects the shear stress distribution. When the gradient index increases, the maximum shear stress moves to the beam surfaces.
Figure 8 shows the effects of the total weight fraction on the shear stress distributions. As observed in Figure 8, when the gradient index is fixed, the increase in GPL weight fraction increases the shear stresses near the top and bottom surfaces and reduces those in the center portion of the beam.

5.2. Elastic Vibration Behavior

Table 8, Table 9 and Table 10 show the dimensionless fundamental frequencies of FG-GPLRC beams with different gradient indices and weight fractions. The total weight fraction of the GPL is 0–1.5% with an increment of 0.25%. Three immovable boundary conditions are considered.
The numerical results in Table 8, Table 9 and Table 10 indicate that the change in gradient index can significantly affect the fundamental frequency and consequently affect the elastic vibration characteristics of the composite beam. The effect becomes increasingly obvious when the total weight fraction increases. For a fixed concentration of GPLs in the FG-GPLRC beam, when the gradient index increases, the dimensionless frequency increases. As mentioned, the increase in gradient index implies that more GPLs are dispersed near the top and bottom surfaces of the beam. Therefore, the above phenomenon proves that dispersing more GPL nanofillers near the top and bottom surfaces of the beam is most effective in enhancing the beam stiffness.
The curves in Figure 9 reflect the relationship between the fundamental frequencies and the weight fractions. The results in Figure 9 indicate that the fundamental frequencies increase when the GPL concentration increases regardless of the gradient index. The aforementioned phenomenon clearly indicates that low-content GPL additions significantly improve the stiffness of the beam.
Figure 10 shows the change in dimensionless fundamental frequencies with respect to the gradient index. The total weight fraction of GPLs is 1.25%.
Figure 10a shows the results of the thick beam (L/h = 5), and Figure 10b is related to the thin one (L/h = 20). We can reach a contrary conclusion with that of the bending behaviors: for the thick beam, when the gradient index increases, the fundamental frequencies first increase and subsequently decrease. There is a suitable r value that can result in a peak fundamental frequency, which indicates that the FG-GPLRC beam with this r value is the stiffest. However, this value is not constant and related to the geometric, weight fraction and boundary conditions.
For the thin beam (L/h = 20), the fundamental frequency increases with the increase in gradient index. A remarkable variation is observed when r is relatively small, and the variation becomes slight when r increases. This result demonstrates that the frequency is more sensitive to the gradient index when the gradient index is relatively small.
As mentioned before, the multilayer structures are introduced to solve the manufacturing problems caused by the distribution laws in practice. It is expected that as a total layer number increase, the functionally graded graphene reinforced multilayer beam will have the some static and vibration response to the monolithic FG-GPLRC beams. Therefore, we calculated the static displacements and fundamental frequencies of the multilayered nanocomposite beams with various layer numbers and made a comparison with those from the monolithic ones. It is worth to note that the weight fractions of each layer GPL-reinforced nanocomposites are determined by substituting the corresponding z-coordinate of mid-plane of each layer into Equation (8). In this example, the length-to-height ratio is set to be 5, and the total weight fraction of GPLs is 1.0%. The clamped-clamped boundary condition is taken into account. Table 11 lists the maximum displacements and fundamental frequencies of GPL-reinforced multilayer beams for different gradient index r. One can observe that when the total layer number is beyond 12, the results of multilayer beams are very closed to the results of monolithic ones, implying that a multilayer GPLRC beam with 12 or more layers is an excellent approximation for an ideal functionally graded beam structure with a continuous and smooth variation in both material composition and properties. In addition, this result also demonstrates that the new FG-GPLRC structures proposed in present work can be realized easily using the multilayer alternatives in the manufacturing processing.

6. Conclusions

A new distribution law based on the error function has been proposed and is further employed to introduce a novel FG-GPLRC. The GPL distributions in the thickness direction can be adjusted by a gradient index to maintain a constant total GPL weight fraction. The modified Halpin–Tsai micromechanics model is used to evaluate the effective Young’s modulus of the FG-GPLRCs.
A computational modeling based on an improved third-order shear deformation theory and Chebyshev–Ritz methodology is developed. The comparisons prove that the developed modeling is more accurate and efficient. The bending deflections, stresses, and natural frequencies of a novel FG-GPLRC beam are investigated using the proposed modeling. The parameter studies are performed, and some conclusions are summarized based the numerical results:
(1)
With respect to the GPL distribution patterns, low-content GPL additions can significantly improve the stiffness and bending resistance of the beam.
(2)
For the thick FG-GPLRC beam, the new distribution law can produce the stiffest FG-GPLRC beams with the lowest bending displacements and highest fundamental frequencies. In other words, by adjusting the gradient index, the most optimized distribution law for the new FG-GPLRCs can be found, which results in an FG-GPLRC beam with the greatest bending resistance and vibration stiffness.
(3)
For the thin FG-GPLRC beam, the increase in gradient index reduces the bending displacements and increases the fundamental frequencies. In addition, if the gradient index is beyond 1.6, the new FG-GPLRCs exhibit better capabilities in the static and vibration analysis than the common FG-X ones.
(4)
A multilayer GPLRC beam with 12 or more layers is an ideal alternative structure for fabricating the new FG-GPLRC structures.

Author Contributions

Conceptualization, Y.W. and K.X.; Methodology, Y.W. and K.X.; Formal analysis, Y.W. C.S. and K.X.; Validation, Y.W. and K.X.; Writing—original draft, Y.W.; Writing—review & editing, T.F. and C.S., Supervision, T.F. and C.S.

Funding

The research was funded by the National Natural Science Foundation of China grant number 51976097, and the Science Fund for Creative Research Groups of NSFC grant number 51621062.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. García-Macías, E.; Rodríguez-Tembleque, L.; Sáez, A. Bending and free vibration analysis of functionally graded graphene vs. carbon nanotube reinforced composite plates. Compos. Struct. 2018, 186, 123–138. [Google Scholar] [CrossRef]
  2. Coleman, J.N.; Khan, U.; Blau, W.J.; Yurii, K.G. Small but strong: A review of the mechanical properties of carbon nanotube–polymer composites. Carbon 2006, 44, 1624–1652. [Google Scholar] [CrossRef]
  3. Viculis, L.M. Intercalation and exfoliation routes to graphite nanoplatelets. J. Mater. Chem. 2005, 15, 974–978. [Google Scholar] [CrossRef]
  4. Lu, W.; Lin, H.; Wu, D.; Chen, G. Unsaturated polyester resin/graphite nanosheet conducting composites with a low percolation threshold. Polymer 2006, 47, 4440–4444. [Google Scholar] [CrossRef]
  5. Si, Y.; Samulski, E.T. Synthesis of water soluble graphene. Nano Lett. 2008, 8, 1679–1682. [Google Scholar] [CrossRef] [PubMed]
  6. Spitalsky, Z.; Tasis, D.; Papagelis, K.; Galiotis, C. Carbon nanotube—Polymer composites: Chemistry, processing, mechanical and electrical properties. Prog. Polym. Sci. 2010, 35, 357–401. [Google Scholar] [CrossRef]
  7. Dang, Z.M.; Yuan, J.K.; Zha, J.W.; Zhou, T.; Li, S.T.; Hu, G.H. Fundamentals, processes and applications of high-permittivity polymer–matrix composites. Prog. Mater. Sci. 2012, 57, 660–723. [Google Scholar] [CrossRef]
  8. Ramanathan, T.; Abdala, A.A.; Stankovich, S.; Dikin, D.A.; Herrera-Alonso, M.; Piner, R.D.; Adamson, D.H.; Schniepp, H.C.; Chen, X.R.R.S.; Ruoff, R.S.; et al. Functionalized graphene sheets for polymer nanocomposites. Nat. Nanotechnol. 2008, 3, 327–331. [Google Scholar] [CrossRef]
  9. Formica, G.; Milicchio, F.; Lacarbonara, W. Hysteretic damping optimization in carbon nanotube nanocomposites. Compos. Struct. 2018, 194, 633–642. [Google Scholar] [CrossRef]
  10. Sengupta, R.; Bhattacharya, M.; Bandyopadhyay, S.; Bhowmick, A.K. A review on the mechanical and electrical properties of graphite and modified graphite reinforced polymer composites. Prog. Polym. Sci. 2011, 36, 638–670. [Google Scholar] [CrossRef]
  11. Koizumi, M. FGM Activities in Japan. Compos. Part B Eng. 1997, 28, 1–4. [Google Scholar] [CrossRef]
  12. Wang, Y.; Wu, D. Thermal effect on the dynamic response of axially functionally graded beam subjected to a moving harmonic load. Acta Astronaut. 2016, 127, 171–181. [Google Scholar] [CrossRef]
  13. Gupta, A.; Talha, M. Recent development in modeling and analysis of functionally graded materials and structures. Prog. Aerosp. Sci. 2015, 79, 1–14. [Google Scholar] [CrossRef]
  14. Sahmania, S.; Safaei, B. Nonlinear free vibrations of bi-directional functionally graded micro/nano-beams including nonlocal stress and microstructural strain gradient size effects. Thin Walled Struct. 2019, 140, 342–356. [Google Scholar] [CrossRef]
  15. Mohammadi, M.; Arefi, M.; Dimitri, R.; Tornabene, F. Higher-Order Thermo-Elastic Analysis of FG-CNTRC Cylindrical Vessels Surrounded by a Pasternak Foundation. Nanomaterials 2019, 9, 79. [Google Scholar] [CrossRef]
  16. Liew, K.M.; Lei, Z.X.; Yu, J.L.; Zhang, L.W. Postbuckling of carbon nanotube-reinforced functionally graded cylindrical panels under axial compression using a meshless approach. Comput. Methods Appl. Mech. Eng. 2014, 268, 1–17. [Google Scholar] [CrossRef]
  17. Frikha, A.; Zghal, S.; Dammak, F. Dynamic analysis of functionally graded carbon nanotubes-reinforced plate and shell structures using a double directors finite shell element. Aerosp. Sci. Technol. 2018, 78, 438–451. [Google Scholar] [CrossRef]
  18. Song, M.; Kitipornchai, S.; Yang, J. Free and forced vibrations of functionally graded polymer composite plates reinforced with graphene nanoplatelets. Compos. Struct. 2016, 159, 579–588. [Google Scholar] [CrossRef]
  19. Wang, Y.; Xie, K.; Shi, C.; Fu, T. Nonlinear bending of axially functionally graded microbeams reinforced by graphene nanoplatelets in thermal environments. Mater. Res. Express 2019, 6, 085615. [Google Scholar] [CrossRef]
  20. Mao, J.J.; Zhang, W. Buckling and post-buckling analyses of functionally graded graphene reinforced piezoelectric plate subjected to electric potential and axial forces. Compos. Struct. 2019, 216, 392–405. [Google Scholar] [CrossRef]
  21. Arefi, M.; Bidgoli, E.M.R.; Dimitri, R.; Tornabene, F. Free vibrations of functionally graded polymer composite nanoplates reinforced with graphene nanoplatelets. Aerosp. Sci. Technol. 2018, 81, 108–117. [Google Scholar] [CrossRef]
  22. Ke, L.L.; Yang, J.; Kitipornchai, S. Nonlinear free vibration of functionally graded carbon nanotube-reinforced composite beams. Compos. Struct. 2010, 92, 676–683. [Google Scholar] [CrossRef]
  23. Tam, M.; Yang, Z.; Zhao, S.; Yang, J. Vibration and Buckling Characteristics of Functionally Graded Graphene Nanoplatelets Reinforced Composite Beams with Open Edge Cracks. Materials 2019, 12, 1412. [Google Scholar] [CrossRef] [PubMed]
  24. Wu, Z.; Zhang, Y.; Yao, G.; Zhou, Y. Nonlinear primary and super-harmonic resonances of functionally graded carbon nanotube reinforced composite beams. Int. J. Mech. Sci. 2019, 153–154, 321–340. [Google Scholar] [CrossRef]
  25. Bahaadini, R.; Saidi, A.R. Aeroelastic analysis of functionally graded rotating blades reinforced with graphene nanoplatelets in supersonic flow. Aerosp. Sci. Technol. 2018, 80, 381–391. [Google Scholar] [CrossRef]
  26. Banić, D.; Bacciocchi, M.; Tornabene, F.; Ferreira, A.J.M. Influence of Winkler-Pasternak Foundation on the Vibrational Behavior of Plates and Shells Reinforced by Agglomerated Carbon Nanotubes. Appl. Sci. 2017, 7, 1228. [Google Scholar] [CrossRef]
  27. Yang, B.; Mei, J.; Chen, D.; Yu, F.; Yang, J. 3D thermo-mechanical solution of transversely isotropic and functionally graded graphene reinforced elliptical plates. Compos. Struct. 2018, 184, 1040–1048. [Google Scholar] [CrossRef]
  28. Yang, B.; Kitipornchai, S.; Yang, Y.F.; Yang, J. 3D thermo-mechanical bending solution of functionally graded graphene reinforced circular and annular plates. Appl. Math. Model. 2017, 49, 69–86. [Google Scholar] [CrossRef]
  29. Wang, Y.; Xie, K.; Fu, T.; Shi, C. Vibration response of a functionally graded graphene nanoplatelet reinforced composite beam under two successive moving masses. Compos. Struct. 2019, 209, 928–939. [Google Scholar] [CrossRef]
  30. Feng, C.; Kitipornchai, S.; Yang, J. Nonlinear free vibration of functionally graded polymer composite beams reinforced with graphene nanoplatelets (GPLs). Eng. Struct. 2017, 140, 110–119. [Google Scholar] [CrossRef]
  31. Wu, H.; Yang, J.; Kitipornchai, S. Dynamic instability of functionally graded multilayer graphene nanocomposite beams in thermal environment. Compos. Struct. 2017, 162, 244–254. [Google Scholar] [CrossRef] [Green Version]
  32. Song, M.; Yang, J.; Kitipornchai, S.; Zhu, W. Buckling and postbuckling of biaxially compressed functionally graded multilayer graphene nanoplatelet-reinforced polymer composite plates. Int. J. Mech. Sci. 2017, 131, 345–355. [Google Scholar] [CrossRef] [Green Version]
  33. Shen, H.S.; Xiang, Y.; Lin, F. Thermal buckling and postbuckling of functionally graded graphene-reinforced composite laminated plates resting on elastic foundations. Thin Walled Struct. 2017, 118, 229–237. [Google Scholar] [CrossRef]
  34. Shen, H.S.; Xiang, Y.; Fan, Y. Nonlinear vibration of functionally graded graphenereinforced composite laminated cylindrical shells in thermal environments. Compos. Struct. 2017, 182, 447–456. [Google Scholar] [CrossRef]
  35. Mao, J.J.; Zhang, W. Linear and nonlinear free and forced vibrations of graphene reinforced piezoelectric composite plate under external voltage excitation. Compos. Struct. 2018, 201, 551–565. [Google Scholar] [CrossRef]
  36. Wang, A.; Chen, H.; Hao, Y.; Zhang, W. Vibration and bending behavior of functionally graded nanocomposite doubly-curved shallow shells reinforced by graphene nanoplatelets. Results Phys. 2018, 9, 550–559. [Google Scholar] [CrossRef]
  37. Gholami, R.; Ansari, R. Nonlinear harmonically excited vibration of third-order shear deformable functionally graded graphene platelet-reinforced composite rectangular plates. Eng. Struct. 2017, 156, 197–209. [Google Scholar] [CrossRef]
  38. Lezgy-Nazargah, M.; Salahshuran, S. A new mixed-field theory for bending and vibration analysis of multi-layered composite plate. Arch. Civ. Mech. Eng. 2018, 18, 818–832. [Google Scholar] [CrossRef]
  39. Wattanasakulpong, N.; Ungbhakorn, V. Analytical solutions for bending, buckling and vibration responses of carbon nanotube-reinforced composite beams resting on elastic foundation. Comput. Mater. Sci. 2013, 71, 201–208. [Google Scholar] [CrossRef]
  40. Ghannadpour, S.A.M.; Mohammadi, B.; Fazilati, J. Bending, buckling and vibration problems of nonlocal Euler beams using Ritz method. Compos. Struct. 2013, 96, 584–589. [Google Scholar] [CrossRef]
  41. Zhao, Z.; Feng, C.; Wang, Y.; Yang, J. Bending and vibration analysis of functionally graded trapezoidal nanocomposite plates reinforced with graphene nanoplatelets (GPLs). Compos. Struct. 2017, 180, 799–808. [Google Scholar] [CrossRef]
  42. Formica, G.; Lacarbonara, W.; Alessi, R. Vibrations of carbon nanotube-reinforced composite. J. Sound Vib. 2010, 329, 1875–1889. [Google Scholar] [CrossRef]
  43. Formica, G.; Talò, M.; Lacarbonara, W. Nonlinear modeling of carbon nanotube composites dissipation due to interfacial stick–slip. Int. J. Plast. 2014, 53, 148–163. [Google Scholar] [CrossRef]
  44. Odegard, G.M.; Gates, T.S.; Wise, K.E.; Park, C.; Siochi, E.J. Constitutive modeling of nanotube–reinforced polymer composites. Compos. Sci. Technol. 2003, 63, 1671–1687. [Google Scholar] [CrossRef]
  45. Talò, M.; Krause, B.; Pionteck, J.; Lanzara, G.; Lacarbonara, W. An updated micromechanical model based on morphological characterization of carbon nanotube nanocomposites. Compos. Part B Eng. 2017, 115, 70–78. [Google Scholar] [CrossRef]
  46. Formica, G.; Talò, M.; Lanzara, G.; Lacarbonara, W. Parametric identification of carbon nanotube nanocomposites constitutive response. J. Appl. Mech. 2019, 86, 041007. [Google Scholar] [CrossRef]
  47. Rafiee, A.M.; Rafiee, J.; Wang, Z.; Song, H.; Yu, Z.Z.; Koratkar, N. Enhanced mechanical properties of nanocomposites at low graphene content. ACS Nano 2009, 3, 3884–3890. [Google Scholar] [CrossRef]
  48. Shokrieh, M.; Esmkhani, M.; Shokrieh, Z.; Zhao, Z. Stiffness prediction of graphene nanoplatelet/epoxy nanocomposites by a combined molecular dynamics—Micromechanics method. Comput. Mater. Sci. 2014, 92, 444–450. [Google Scholar] [CrossRef]
  49. Tha, H.C.; Ferreira, A.J.M.; Tran, T.D.; Phung-Van, P. Free vibration, buckling and bending analyses of multilayer functionally graded graphene nanoplatelets reinforced composite plates using the NURBS formulation. Compos. Struct. 2019, 220, 749–759. [Google Scholar] [CrossRef]
  50. Gholami, R.; Ansari, R. Large deflection geometrically nonlinear analysis of functionally graded multilayer graphene platelet-reinforced polymer composite rectangular plates. Compos. Struct. 2017, 180, 760–771. [Google Scholar] [CrossRef]
  51. Zhang, Z.; Li, Y.; Wu, H.; Zhang, H.; Wu, H.; Jiang, S.; Chai, G. Mechanical analysis of functionally graded graphene oxide-reinforced composite beams based on the first-order shear deformation theory. Mech. Adv. Mater. Struct. 2018. [Google Scholar] [CrossRef]
  52. Shi, G. A new simple third-order shear deformation theory of plates. Int. J. Solids Struct. 2007, 44, 4399–4417. [Google Scholar] [CrossRef] [Green Version]
  53. Nuttawit, W.; Gangadhara, B.P.; Donald, W.K. Thermal buckling and elastic vibration of third-order shear deformable functionally graded beams. Int. J. Mech. Sci. 2011, 53, 734–743. [Google Scholar]
  54. Reddy, J.N. Mechanics of Laminated Composite Plates and Shells: Theory and Application, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  55. Wang, L.; Liu, Y.; Liu, Y. An inverse method for distributed dynamic load identification of structures with interval uncertainties. Adv. Eng. Softw. 2019, 131, 77–89. [Google Scholar] [CrossRef]
  56. Liang, K.; Sun, Q.; Liu, X.R. Investigation on imperfection sensitivity of composite cylindrical shells using the nonlinearity reduction technique and the polynomial chaos method. Acta Astronaut. 2018, 146, 349–358. [Google Scholar] [CrossRef]
  57. Fox, L.; Parker, I.B. Chebyshev Polynomials in Numerical Analysis; Oxford University Press: London, UK, 1968. [Google Scholar]
  58. Zhou, D.; Lo, S.H.; Au, F.T.K.; Cheung, Y.K.; Liu, W.Q. 3-D vibration analysis of skew thick plates using Chebyshev–Ritz method. Int. J. Mech. Sci. 2006, 48, 1481–1493. [Google Scholar] [CrossRef]
  59. Nguyen, T.K.; Nguyen, B.D. A new higher-order shear deformation theory for static, buckling and free vibration analysis of functionally graded sandwich beams. J. Sandw. Struct. Mater. 2015, 17, 613–631. [Google Scholar] [CrossRef]
  60. Vo, T.P.; Thai, H.T.; Nguyen, T.K.; Inam, F.; Lee, J. Static behaviour of functionally graded sandwich beams using a quasi-3D theory. Compos. Part B Eng. 2015, 68, 59–74. [Google Scholar] [CrossRef] [Green Version]
  61. Nguyen, T.K.; Vo, T.P.; Nguyen, B.D.; Lee, J. An analytical solution for buckling and vibration analysis of functionally graded sandwich beams using a quasi-3D shear deformation theory. Compos. Struct. 2016, 156, 238–252. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Variation of the weight fraction of the graphene nanoplatelets (GPLs) versus r.
Figure 1. Variation of the weight fraction of the graphene nanoplatelets (GPLs) versus r.
Nanomaterials 09 01690 g001
Figure 2. Schematic of GPL distributions of the proposed model: (a) without GPLs; (b) uniform distribution; (c) FG r = 2; (d) FG r = 5.
Figure 2. Schematic of GPL distributions of the proposed model: (a) without GPLs; (b) uniform distribution; (c) FG r = 2; (d) FG r = 5.
Nanomaterials 09 01690 g002
Figure 3. Variation in the maximum dimensionless displacement w ¯ max with respect to the GPL weight fraction g G P L 0 .
Figure 3. Variation in the maximum dimensionless displacement w ¯ max with respect to the GPL weight fraction g G P L 0 .
Nanomaterials 09 01690 g003
Figure 4. Relationship between maximum dimensionless displacement w ¯ max and gradient index r: (a) L/h = 5; (b) L/h = 20.
Figure 4. Relationship between maximum dimensionless displacement w ¯ max and gradient index r: (a) L/h = 5; (b) L/h = 20.
Nanomaterials 09 01690 g004
Figure 5. Normal stress distributions of the FG-GPLRC beam with different r values.
Figure 5. Normal stress distributions of the FG-GPLRC beam with different r values.
Nanomaterials 09 01690 g005
Figure 6. Effects of the total weight fractions of GPLs on the normal stress distributions: (a) r = 1; (b) r = 4.
Figure 6. Effects of the total weight fractions of GPLs on the normal stress distributions: (a) r = 1; (b) r = 4.
Nanomaterials 09 01690 g006
Figure 7. Effects of the gradient index on the shear stress distributions of the FG-GPLRC beam.
Figure 7. Effects of the gradient index on the shear stress distributions of the FG-GPLRC beam.
Nanomaterials 09 01690 g007
Figure 8. Effects of the total weight fractions of GPLs on the shear stress distribution: (a) r = 1; (b) r = 4.
Figure 8. Effects of the total weight fractions of GPLs on the shear stress distribution: (a) r = 1; (b) r = 4.
Nanomaterials 09 01690 g008
Figure 9. Variation of the dimensionless fundamental frequency with respect to the gradient index.
Figure 9. Variation of the dimensionless fundamental frequency with respect to the gradient index.
Nanomaterials 09 01690 g009
Figure 10. Variations of the maximum dimensionless displacements with respect to gradient index r: (a) L/h = 5; (b) L/h = 20.
Figure 10. Variations of the maximum dimensionless displacements with respect to gradient index r: (a) L/h = 5; (b) L/h = 20.
Nanomaterials 09 01690 g010
Table 1. Values of the indices for various boundary conditions.
Table 1. Values of the indices for various boundary conditions.
Boundary ConditionsLuLϕLwRuRϕRw
H-H101101
C-C112112
C-H112101
Table 2. Deflection w ¯ , normal stress σ ¯ x x , shear stress σ ¯ x z and natural frequency of FG sandwich beam.
Table 2. Deflection w ¯ , normal stress σ ¯ x x , shear stress σ ¯ x z and natural frequency of FG sandwich beam.
N w ¯ σ ¯ x x ( 0.5 L , 0.5 h ) σ ¯ x z ( 0 , 0 ) ω1ω2ω3
20.26650.77930.84954.667728.135835.2954
30.26640.77830.91044.667218.262227.9520
40.32811.16250.84764.236318.184627.9520
50.32811.16240.84484.236315.150727.9518
60.32811.16240.84164.235115.099027.9518
70.32811.16240.84164.235115.099027.9518
80.32811.16240.84164.235115.099029.5854
Ref. [59]0.32821.16320.83714.3052
Table 3. Comparisons of the vertical dimensionless mid-span displacements of functionally graded (FG) sandwich beams.
Table 3. Comparisons of the vertical dimensionless mid-span displacements of functionally graded (FG) sandwich beams.
SourceH-HC-C
p = 1p = 5p = 10p = 1p = 5p = 10
FSDT [60]5.44088.14099.02321.37702.06352.2614
TSDT [60]5.41228.57629.48001.33721.98962.1747
Qusi-3D [60]5.36128.51379.40501.30771.94162.1211
Present5.38228.51409.40491.31261.96372.1487
Table 4. Comparisons of the dimensionless fundamental frequencies of the FG sandwich beams.
Table 4. Comparisons of the dimensionless fundamental frequencies of the FG sandwich beams.
SourceH-HC-C
p = 1p = 2p = 5p = 1p = 2p = 5
HSDT [61]4.11053.73343.37718.37477.71497.0723
TSDT [61]4.11053.73343.37718.37057.71147.0691
Qusi-3D [61]4.11853.74103.38408.46537.80087.1550
Present4.11463.73743.37918.41637.77317.1265
Table 5. Maximum dimensionless displacements of the functionally graded graphene nanoplatelet reinforced composite (FG-GPLRC) beams (H-H).
Table 5. Maximum dimensionless displacements of the functionally graded graphene nanoplatelet reinforced composite (FG-GPLRC) beams (H-H).
L/hGradient Index0%0.25%0.5%0.75%1.0%1.25%1.5%
5uniform15.42478.01965.41914.09273.28822.74832.3609
r = 215.42476.45314.15063.07422.44602.03291.7402
r = 415.42475.73563.75192.85762.33501.98651.7351
r = 615.42475.49023.64682.83942.37272.06181.8364
r = 815.42475.36023.58762.82742.39402.10811.9019
r = 1015.42475.27863.54692.81462.40192.13231.9394
20uniform13.91877.23664.89013.69312.96722.48002.1304
r = 213.91875.46333.40462.47441.94421.60151.3618
r = 413.91874.57282.75281.97521.54291.26731.0760
r = 613.91874.28642.55601.83031.43021.17641.0007
r = 813.91874.14792.46281.76271.37851.13550.9676
r = 1013.91874.06662.40871.72381.34901.11240.9492
Table 6. Maximum dimensionless displacements of the FG-GPLRC beams (C-C).
Table 6. Maximum dimensionless displacements of the FG-GPLRC beams (C-C).
L/hGradient Index0%0.25%0.5%0.75%1.0%1.25%1.5%
5uniform4.27182.22101.50081.13350.91070.76110.6538
r = 24.27182.06301.40991.07980.87760.74020.6404
r = 44.27182.04891.52061.24911.07340.94660.8494
r = 64.27182.02971.56951.34271.19631.08911.0048
r = 84.27182.00891.58361.38341.25751.16631.0947
r = 104.27181.99121.58531.40111.28841.20831.1461
20uniform2.85951.48671.00460.75870.60960.50950.4377
r = 22.85951.14240.71840.52500.41400.34190.2913
r = 42.85950.97280.60060.43920.34820.28940.2482
r = 62.85950.91750.56570.41650.33320.27950.2419
r = 82.85950.89020.54880.40570.32640.27570.2402
r = 102.85950.87390.53860.39920.32240.27340.2393
Table 7. Maximum dimensionless displacements of the FG-GPLRC beams (C-H).
Table 7. Maximum dimensionless displacements of the FG-GPLRC beams (C-H).
L/hGradient Index0%0.25%0.5%0.75%1.0%1.25%1.5%
5uniform7.51353.90642.63971.99361.60171.33871.1500
r = 27.51353.40082.26511.71091.37841.15570.9956
r = 47.51353.22362.27611.81711.53281.33441.1859
r = 67.51353.14922.29601.89721.65171.47891.3474
r = 87.51353.10022.29461.93031.71051.55721.4409
r = 107.51353.06512.28631.94271.73901.59901.4938
20uniform5.85863.04602.05831.55451.24891.04390.8967
r = 25.85862.31781.45031.05670.83170.68590.5837
r = 45.85861.95521.19080.86210.67810.56020.4779
r = 65.85861.83791.11310.80760.63820.53000.4546
r = 85.85861.78071.07590.78200.62000.51700.4454
r = 105.85861.74681.05410.76700.60950.50950.4403
Table 8. Dimensionless fundamental frequencies of the FG-GPLRC beams (H-H).
Table 8. Dimensionless fundamental frequencies of the FG-GPLRC beams (H-H).
L/hGradient Index0%0.25%0.5%0.75%1.0%1.25%1.5%
5uniform2.75103.81584.64275.34325.96216.52257.0386
r = 22.75104.26515.32566.19356.94787.62478.2444
r = 42.75104.53505.62416.45727.15387.76438.3151
r = 62.75104.63925.71356.49137.11417.64218.1065
r = 82.75104.69695.76446.51127.09027.56727.9764
r = 102.75104.73405.79956.52927.08277.52917.9046
20uniform2.92684.05974.93945.68466.34306.93937.4884
r = 22.92684.67255.92026.94577.83728.63669.3677
r = 42.92685.10776.58507.77598.80029.712310.5423
r = 62.92685.27586.83438.07899.141810.082810.9349
r = 82.92685.36326.96278.23289.312510.264011.1221
r = 102.92685.41667.04058.32569.414310.370811.2306
Table 9. Dimensionless fundamental frequencies of the FG-GPLRC beams (C-C).
Table 9. Dimensionless fundamental frequencies of the FG-GPLRC beams (C-C).
L/hGradient Index0%0.25%0.5%0.75%1.0%1.25%1.5%
5uniform5.31977.37888.977810.332411.529112.612913.6109
r = 25.31977.68109.302510.636711.804312.858213.8276
r = 45.31977.72268.98049.917810.705811.405512.0454
r = 65.31977.76448.84809.576910.153610.647811.0904
r = 85.31977.80768.81359.44129.910610.297210.6339
r = 105.31977.84418.81189.38549.795610.121510.3976
20uniform6.53869.069611.034912.699914.170915.503016.7297
r = 26.538610.347213.051115.269717.198018.927220.5089
r = 46.538611.214014.276916.699918.760320.582022.2321
r = 66.538611.547414.711217.151419.183220.948422.5245
r = 86.538611.723314.937517.378719.381221.096422.6079
r = 106.538611.832215.078017.520119.503621.185422.6525
Table 10. Dimensionless fundamental frequencies of the FG-GPLRC beams (C-H).
Table 10. Dimensionless fundamental frequencies of the FG-GPLRC beams (C-H).
L/hGradient Index0%0.25%0.5%0.75%1.0%1.25%1.5%
5uniform3.97675.51596.71127.72388.61849.428610.1746
r = 23.97675.93417.28358.38869.351710.218211.0134
r = 43.97676.11267.29688.18048.91689.564810.1527
r = 63.97676.19017.27568.02008.60649.10389.5442
r = 83.97676.24177.28267.95718.46468.87989.2376
r = 103.97676.27887.29837.93528.39868.76699.0768
20uniform4.54066.29827.66308.81929.840710.765811.6176
r = 24.54067.22039.129610.697612.060613.282914.4007
r = 44.54067.862210.078111.848813.363514.707615.9278
r = 64.54068.109610.425212.244513.779515.126016.3366
r = 84.54068.239110.604512.444413.981915.318116.5093
r = 104.54068.318510.714412.566214.103715.431416.6070
Table 11. Maximum displacements w ¯ max and fundamental frequencies Ω1 of multilayer beams.
Table 11. Maximum displacements w ¯ max and fundamental frequencies Ω1 of multilayer beams.
Total Layer Number w ¯ max Ω1
r = 1r = 3r = 5r = 1r = 3r = 5
20.88101.84314.106711.96508.26865.5381
40.86330.94091.322112.085011.55819.7511
80.86390.96741.119012.080911.391010.5805
120.86430.97201.132312.077811.363710.5198
160.86450.97431.136812.076411.350410.4995
200.86460.97751.142512.074811.332410.4871
500.86480.97771.144212.074611.330410.4656
Monolithic0.86480.97771.145212.074611.330210.4610

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/nano9121690

AMA Style

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(12):1690. https://doi.org/10.3390/nano9121690

Chicago/Turabian Style

Wang, Yuewu, Ke Xie, Tairan Fu, and Congling Shi. 2019. "Bending and Elastic Vibration of a Novel Functionally Graded Polymer Nanocomposite Beam Reinforced by Graphene Nanoplatelets" Nanomaterials 9, no. 12: 1690. https://doi.org/10.3390/nano9121690

APA Style

Wang, Y., Xie, K., Fu, T., & Shi, C. (2019). Bending and Elastic Vibration of a Novel Functionally Graded Polymer Nanocomposite Beam Reinforced by Graphene Nanoplatelets. Nanomaterials, 9(12), 1690. https://doi.org/10.3390/nano9121690

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