Next Article in Journal
A High-Strain-Rate Superplasticity of the Al-Mg-Si-Zr-Sc Alloy with Ni Addition
Next Article in Special Issue
Correction: Wang et al. Parametric Formula for Stress Concentration Factor of Fillet Weld Joints with Spline Bead Profile. Materials 2020, 13, 4639
Previous Article in Journal
One-Step Etching Characteristics of ITO/Ag/ITO Multilayered Electrode in High-Density and High-Electron-Temperature Plasma
Previous Article in Special Issue
Dynamics of Space-Fractional Euler–Bernoulli and Timoshenko Beams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Prediction for Nonlinear Buckling of Elastically Supported FG-GPLRC Arches under a Central Point Load

1
College of Urban and Rural Construction, Zhongkai University of Agriculture and Engineering, Guangzhou 510225, China
2
Wind and Vibration Engineering Research Center, Guangzhou University, Guangzhou 510006, China
3
School of Engineering, RMIT University, P.O. Box 71, Bundoora, VIC 3083, Australia
4
Department of Civil and Environmental Engineering, The Hong Kong Polytechnic University, Kowloon, Hong Kong, China
5
Hong Kong Branch of National Rail Transit Electrification and Automation Engineering Technology Research Center, The Hong Kong Polytechnic University, Kowloon, Hong Kong, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(8), 2026; https://doi.org/10.3390/ma14082026
Submission received: 24 March 2021 / Revised: 15 April 2021 / Accepted: 16 April 2021 / Published: 17 April 2021

Abstract

:
In this paper, we present an analytical prediction for nonlinear buckling of elastically supported functionally graded graphene platelet reinforced composite (FG-GPLRC) arches with asymmetrically distributed graphene platelets (GPLs). The effective material properties of the FG-GPLRC arch are formulated by the modified Halpin–Tsai micromechanical model. By using the principle of virtual work, analytical solutions are derived for the limit point buckling and bifurcation buckling of the FG-GPLRC arch subjected to a central point load (CPL). Subsequently, the buckling mode switching phenomenon of the FG-GPLRC arch is presented and discussed. We found that the buckling modes of the FG-GPLRC arch are governed by the GPL distribution pattern, rotational restraint stiffness, and arch geometry. In addition, the number of limit points in the nonlinear equilibrium path of the FG-GPLRC arch under a CPL can be determined according to the bounds of successive inflexion points. The effects of GPL distribution patterns, weight fractions, and geometric configurations on the nonlinear buckling behavior of elastically supported FG-GPLRC arches are also comprehensively discussed.

1. Introduction

Functionally graded material (FGM) structures, characterized by a continuous change in the material compositions along one or multiple directions, have attracted extensive attention from both research and industrial communities owing to their excellent stiffness and strength-to-weight properties as compared with homogeneous composite structures [1,2,3]. To better understand the performance of FGM structures in practical engineering applications, researchers have conducted a series of investigations on the structural behavior of FGM structures. For instance, Ke et al. [4] presented analytical solutions for the nonlinear vibration responses of FGM beams with different end supports and discussed the influence of bending-stretching coupling on the nonlinear vibration of FGM beams. Librescu et al. [5] analyzed the vibration and stability of FGM thin-walled beams under a high-temperature environment. Yan et al. reported analytical solutions for the dynamic instability [6] and dynamic responses [7] of FGM beams with open-edge cracks. They investigated the crack depth and location effects on the mechanical behavior of such beams. Nguyen et al. [8] studied the mechanical buckling of stiffened FGM plates by using the finite element method. They found that the addition of stiffeners to the FGM plate could significantly reduce the weight of the FGM plate. Chen et al. [9] investigated the imperfection sensitivity in the nonlinear vibration of initially stressed FGM plates. According to their results, the effects of the initial stress, geometric imperfection, and volume fraction index were quite significant on the nonlinear vibration behavior of the FGM plate. Hao et al. [10] employed an asymptotic perturbation method to analyze the nonlinear oscillations, bifurcations, and chaotic motions of FGM plates. More relevant studies to investigate the significant performance of FGM structures could be found from the open literature [11,12,13,14,15,16,17].
Graphene, as an emerging high-performance nanofillers, has attracted considerable attention in aerospace, mechanical, thermal, and electrical engineering fields. It has been demonstrated by researchers that the reinforced performance of graphene nanoplatelets (GPLs) is significantly superior to other reinforcement materials. Adding a low concentration of GPLs can improve the stiffness and strength of reinforced composites significantly [18,19,20,21,22,23]. By introducing GPLs to FGM materials, novel FG GPLs-reinforced composite (FG-GPLRC) structures have been developed recently and have since attracted extensive attention in both research and engineering communities [24]. Yang and his co-workers conducted pioneering studies on the mechanical behaviors of FG-GPLRC structures, such as beams [25,26,27,28], shells [29,30], and plates [31,32]. Focusing on the stability analysis of FG-GPLRC arches, the authors devoted extensive efforts to the investigation of such structures, including the characteristics of nonlinear static buckling, dynamic buckling, and free vibration for FG-GPLRC arches with different boundary conditions and external loads [33,34,35,36,37,38,39]. In addition, Liu et al. [40] analyzed the nonlinear behavior and stability of functionally graded porous (FGP) arches reinforced by GPLs and obtained the critical buckling load under a uniform load. Zhao et al. [41] discussed the linear buckling, fundamental frequency, and dynamic instability of porous arches using an analytical method. Li et al. [42] further investigated the mechanics of the confined porous arches by including the thermal effect. Hitherto, to the best of our knowledge, far too little attention has been paid to the nonlinear buckling behavior of elastically supported FG-GPLRC arches with asymmetric distributed GPLs under a central point load (CPL).
Therefore, to fill this research gap, the nonlinear buckling behavior of elastically supported FG-GPLRC arches under a CPL is investigated in this study. The effective material properties of the FG-GPLRC arch are formulated by the modified Halpin–Tsai micromechanical model, because the Young’s modulus of the FG-GPLRC arch predicted by this model agreed well with the experimental results [19]. Subsequently, the principle of virtual work is employed to derive the nonlinear buckling load for the limit point buckling and bifurcation buckling modes, from which the critical geometric parameters to identify the buckling load and the number of limit points can also be determined. In addition, the influences of GPL distribution patterns, GPL weight fractions, geometric configurations, as well as the flexibility of the rotational constraints on the nonlinear buckling behavior of the arch are discussed in detail. The bifurcation stability criteria presented in this work can provide essential information for the structural/material design of FG-GPLRC arches that have great potential in various engineering applications, for example, arch-shaped micro-electromechanical devices as sensors and transducers [43], and dielectric elastomer actuators as lightweight speakers [44,45]. In addition, the presence of analytical solutions is useful to engineers and researchers for benchmarking the convergence and validity of numerical methods for arch buckling analysis.

2. Effective Material Properties

In Figure 1, we consider an elastically restrained FG-GPLRC arch with multiple layers, NL, under a CPL Q. The rectangular cross section of the arch is b × h (width × thickness).
It is assumed that the reinforcement of GPLs is uniformly distributed in the isotropic polymer matrix, so each individual GPLRC layer can be regarded as an isotropic and homogeneous material. The variation of GPLs is continuous along the arch thickness in accordance with the power law distribution. The GPL volume fraction of the kth GPLRC layer V GPL k is formulated as follows:
V GPL k = V GPL * [ 1 2 + ( k 1 N L 1 ) n ]
where NL is the total number of the GPLRC layer, and n is the power law index that characterizes the distribution of GPLs. When n = 0, it corresponds to a uniform distribution of GPLs reinforcements in the thickness direction. V GPL * is the total GPL volume fraction which can be determined using:
V GPL * = W GPL W GPL + ( ρ GPL / ρ m ) ( 1 W GPL )
where WGPL is the GPLs weight fraction and ρ GPL and ρ m are the mass densities of GPLs and matrix, respectively.
According to the modified Halpin–Tsai micromechanical model [43], the effective Young’s modulus Ek of the kth GPLRC layer is given by
E k = 3 ( 1 + ξ L η L V GPL k ) 8 ( 1 η L V GPL k ) × E m + 5 ( 1 + ξ T η T V GPL k ) 8 ( 1 η T V GPL k ) × E m
with
η L = ( E GPL / E m ) 1 ( E GPL / E m ) + ξ L ,   η T = ( E GPL / E m ) 1 ( E GPL / E m ) + ξ T ξ L = 2 ( a GPL / b GPL ) × ξ T / 2 ,   ξ T = 2 ( b GPL / t GPL )
where EGPL and Em are the Young’s moduli of GPLs and matrix, respectively. aGPL, bGPL, and tGPL are the length, width, and thickness of GPLs, respectively.
The Poisson’s ratio υ k of each GPLRC layer is determined by the rule of mixture as follows:
υ k = υ GPL V GPL k + υ m ( 1 V GPL k )
where υ GPL and υ m are the Poisson’s ratios of GPLs and matrix, respectively.
Consider GPLs with the material properties EGPL = 1010 GPa, ρ GPL = 1062.5 kg / m 3 , υ GPL = 0.186 , aGPL = 2.5 µm, bGPL = 1.5 µm, and tGPL = 1.5 nm, and epoxy (polymer matrix) with Em = 3 GPa, υ m = 0.34 , and ρ m = 1200 kg / m 3 [25], the effective Young’s moduli of the FG-GPLRC arch with 10 GPLRC layers under various GPLs distributions (different power law indices) are shown in Figure 2.

3. Mathematical Modeling

As shown in Figure 1, an elastically supported FG-GPLRC arch with an angle 2Θ, a radius R, and an arc length S under the effect of a CPL is studied in this section. The stiffness of the elastic rotational constraints at both ends is k, the radial and axial displacement of the arch are v and w, respectively, and θ is the angular coordinate. By using the Donnell’s shallow shell theory [39], the nonlinear strain-displacement relations for the FG-GPLRC arch are adopted as follows:
ε = w ˜ v ˜ + v ˜ 2 2 z R v ˜   with   v ˜ = v / R   and   w ˜ = w / R
where (   ) = d (   ) / d θ .
According to the principle of virtual work, the governing equation is established as:
δ W = Θ Θ R b h / 2 h / 2 E k ε δ ε   d z d θ Θ Θ δ D ( θ ) Q R δ v ˜   d θ + i = ± Θ k v ˜ i δ v ˜ i = 0
where δ D ( θ ) is the Dirac delta function [27].
Substituting Equation (6) into Equation (7), one obtains
δ W = Θ Θ [ N R ( δ w ˜ δ v ˜ + v ˜ δ v ˜ ) M δ v ˜ δ D ( θ ) Q R δ v ˜ ]   d θ + i = ± Θ k v ˜ i δ v ˜ i = 0
with the axial force
N = A 11 ( w ˜ v ˜ + v ˜ 2 2 ) + B 11 R v ˜
and the bending moment
M = B 11 ( w ˜ v ˜ + v ˜ 2 2 ) D 11 R v ˜ = ( D 11 B 11 2 A 11 ) v ˜ R B 11 A 11 N
where
{ A 11   B 11   D 11 } = b h / 2 h / 2 E k { 1   z   z 2 }   dz = b k = 1 N L z k z k + 1 E k { 1   z   z 2 }   dz .
As the present model is derived in accordance with the classical Euler–Bernoulli theory, so that A 11 , B 11 and D 11 are Ek-based constants.
From Equation (8), the boundary conditions for the elastically supported FG-GPLRC arch are given by
v ˜ ( Θ ) = v ˜ ( Θ ) = w ˜ ( Θ ) = w ˜ ( Θ ) = 0 M ( Θ ) + k v ˜ Θ = M ( Θ ) k v ˜ Θ = 0
Substituting Equations (9) and (10) into Equation (8), obtains the equilibrium equations as follows:
( N R ) = 0
N R + ( N R v ˜ ) M δ D ( θ ) Q R = 0
Inserting Equations (9), (10) and (13) into Equation (14) yields the radial equilibrium equation as follows:
v ˜ i v μ 2 + v ˜ = δ D ( θ ) Q R 2 κ μ 2 1   with   κ = D 11 B 11 2 A 11   and   μ 2 = N R 2 κ
where κ is effective bending stiffness of the FG-GPLRC arch.
Solving Equation (11) and considering the boundary conditions of the elastically support of the arch, the analytical solution of the dimensionless radial displacement can be determined as follows:
v ˜ = β 2 μ 2 θ 2 2 μ 2 + β ( cos μ θ cos β ) K 1 μ 2 + B ( cos μ θ cos β ) K 12 R + P μ 2 β [ K 2 cos μ θ β + K 3 + H ( θ ) ( μ θ sin μ θ ) ]
where β = μ Θ , B = B 11 / A 11 , and α is the flexibility of the rotational constraints defined as:
α = κ k S
The dimensionless CPL is defined as:
P = Q R 2 Θ 2 κ
and the parameters K1, K2, K3, and K12 are given as:
K 1 = 2 α + 1 2 α β cos β + sin β ,   K 2 = 2 α β sin β cos β + 1 2 α β cos β + sin β K 3 = 1 cos β 2 α β cos β + sin β ,   K 12 = 2 α β 2 α β cos β + sin β
In Equation (16), H ( θ ) is a step function defined by [33] as:
H ( θ ) = { 1       θ < 0 1           θ 0
Substituting Equation (16) into Equation (9), obtains the nonlinear equilibrium relationship as:
A 1 P 2 + B 1 P + C 1 = 0
where A1, B1, and C1 are determined as:
A 1 = Ξ 1 K 2 2 Ξ 2 K 2 + Ξ 3
B 1 = 2 β 2 Ξ 1 K 1 K 2 β 2 Ξ 2 K 1 K 2 cos β β 3 + sin β K 3 β 3 + B 0 λ [ ( 2 K 2 Ξ 1 Ξ 2 ) β 3 K 12 + K 2 sin β + cos β 1 β 2 ]
C 1 = β 4 Ξ 1 K 1 2 1 6 + β 2 r 11 2 λ 2 + B 0 2 λ 2 ( β 6 Ξ 1 K 12 2 + β K 12 sin β ) + B 0 λ ( 2 β 5 Ξ 1 K 12 K 1 + K 1 sin β + 1 )
with B0 = B11/A11h, r 11 = r x / h , r x = κ / A 11 , and λ is the arch geometrical parameter as
λ = R Θ 2 h = S h Θ 2
The parameters Ξ 1 , Ξ 2 , and Ξ 3 are given as
Ξ 1 = β sin β cos β 4 β 5 ,   Ξ 2 = ( cos β 1 ) 2 2 β 5 ,   Ξ 3 = sin β cos β + 3 β 4 sin β 4 β 5

4. Nonlinear Buckling Analysis

4.1. Limit Point Buckling

When the FG-GPLRC arch buckles in a limit point instability mode, the limit point loads may be either local maxima or local minima on the nonlinear equilibrium path, which are derived from Equation (21) as:
A 2 P 2 + B 2 P + C 2 = 0
with
{ A 2 B 2 C 2 } = β { A 1 B 1 C 1 }
Accordingly, the solutions of the limit point load and nonlinear equilibrium path of FG-GPLRC arches under a CPL can be determined by solving Equations (21) and (27).

4.2. Bifurcation Buckling

When the FG-GPLRC arch buckles in a bifurcation mode, the equilibrium equation for the arch can be derived by substituting the critical states { v ˜ + v ˜ b ,   w ˜ + w ˜ b ,   N + N b ,   M + M b } for the bifurcation buckling into Equation (14), we obtain:
v ˜ b i v μ 2 + v ˜ b = 0
The general solution of Equation (29) is solved as:
v ˜ b = E 1 sin μ θ + E 2 cos μ θ + E 3 θ + E 4
where E1, E2, E3, and E4 are unknown coefficients.
Substituting the boundary conditions of the arch into Equation (30) yields the following characteristic equation for the coefficients E1, E2, E3, and E4 as:
[ sin β 2 α + β 2 ( sin β cos β 2 α β ) ] ( cos β + sin β 2 α β ) = 0
From Equation (31), for the first case, it is true that the first term becomes zero, from which the critical axial force parameter β b for the bifurcation buckling is solved as:
β b = η b π
and the corresponding axial force is obtained as:
N b = μ b 2 κ R 2 = ( η b π ) 2 κ ( S / 2 ) 2 = N E 2
Substituting Equation (32) into Equation (21), the bifurcation buckling equation is determined as:
A 3 P 2 + B 3 P + C 3 = 0
where
A 3 = A 1 ,   β = η b π ,   B 3 = B 1 ,   β = η b π ,   C 3 = C 1 ,   β = η b π
from which the bifurcation load is solved as:
P = B 3 ± B 3 2 4 A 3 C 3 2 A 3
The existence of real solutions in Equation (36) requires B 3 2 4 A 3 C 3 0 , this condition gives a critical geometric parameter λ b 1 that can trigger the bifurcation buckling of the arch.
Consider the second case of Equation (31) (i.e., the second term is zero), we may obtain the critical axial force parameter β s n as follows:
β s n = η s n π   ( s n = s 1 ,   s 2 ,   s 3 ,   , )
Following similar procedures, a critical geometric parameter λ s n that governs the number of inflection points in the nonlinear equilibrium path of the FG-GPLRC arch can also be determined. It is worth noting that the CPL corresponding to the first (as well as the lowest) axial force parameter β s 1 is also the lowest buckling load of the arch whose geometric parameter is λ s 1 . When the FG-GPLRC arch has λ λ s 1 , the arch becomes a slightly curved beam and does not perform typical nonlinear buckling behavior.

5. Numerical Studies and Discussion

In this section, the nonlinear buckling behavior of the FG-GPLRC arch with a cross section of b × h = 0.03 m × 0.025 m and 10 GPLRC layers is investigated. The material properties of GPLs and matrix are adopted as the same in Section 2. To verify the present solutions, a finite element (FE) analysis conducted by ANSYS 17.0 is used to obtain the nonlinear equilibrium path of the FG-GPLRC arch. In numerical modeling, the SHELL181 element is used to establish a multi-layer FG-GPLRC arch model. The COMBIN14 element is adopted to model elastic end restraints whose rotational restraint stiffness k is defined by Equation (17). The section commands are used to define the layered structures which can provide the input options for specifying the thickness, material, and orientation of each layer. The FG-GPLRC arch model is meshed by 100 elements along the arch length direction, and 10 layers are considered to be in the thickness direction, as depicted in Figure 3.
As a large deflection and rotation may occur during the whole deformation of the arch, the NLGEOM command is used for nonlinear analysis when considering the effect of geometric nonlinearities. In addition, the arc-length method activated by the ARCLEN command is, then, used to trace the nonlinear equilibrium path of the FG-GPLRC arch with a load step of 200, as specified by the NSUBST command. By setting the minimum and maximum multipliers for the arc-length radius, the subsequent displacement and load proportional factors and the increment size can be adjusted and computed automatically using the arc-length radius.
For comparison, the FE and present results for the limit point buckling path and bifurcation buckling path are plotted in Figure 4, in which NE0 = the critical bifurcation axial force NE2 of the pure epoxy arch, and vc/f = the dimensionless displacement of the arch crown with f being the rise of the arch. It should be mention that, for triggering the bifurcation buckling behavior of the FG-GPLRC arch, an antisymmetric geometric imperfection of 0.1% arch length is introduced to the arch in the bifurcation buckling analysis. It is observed that the present analytical solutions agree very well with the FE results.
Next, we discuss the buckling load of the elastically supported arch, i.e., the limit point buckling load (the first upper limit point load) for brevity, unless stated otherwise. Figure 5 shows the influence of the power law index on the buckling load and the effective bending stiffness of the arch for α = 2 . In Figure 5b, “D110” is the bending stiffness of the pure epoxy arch. It is found that the buckling load and the effective bending stiffness significantly decrease as the power law index increases, but the effect of the power law index tends to be less pronounced when the power law index is higher than three. This is because a higher content of GPLs would be concentrated on the bottom and less on other layers when the FG-GPLRC arch has a higher power law index, thereby leading to a lower bending stiffness of the arch.
Table 1 lists the buckling load of the arch having different GPL weight fractions. It is found that, by comparing with a pure epoxy arch (i.e., WGPL = 0.0%), GPLs have a remarkably reinforced effect on improving the buckling bearing capacity of the FG-GPLRC arch and the buckling load of the arch increases with an increase in WGPL. In addition, the effects of the geometry and size of GPLs on the buckling load are depicted in Figure 6. For bGPL/tGPL = 103, the reinforced effect of GPLs on improving the buckling bearing capacity of the FG-GPLRC arch becomes stable under different ratios of aGPL/bGPL. Figure 7 presents the flexibility of the rotational constraints on the buckling load of the FG-GPLRC arch. It is noted that, in all cases, the buckling load decreases as the flexibility of the rotational constraints increases, and the fixed FG-GPLRC arch has the highest one, as expected.
Figure 8 displays the nonlinear equilibrium path of the FG-GPLRC arch in a load-displacement curve form to discuss the phenomenon of buckling mode switching. It is observed from Figure 8a that, when the arch has a smaller geometric parameter (e.g., λ = 2.6 ), the bifurcation point is located behind the limit point, so the buckling mode of the arch is the limit point instability mode. However, when the FG-GPLRC arch has a higher geometric parameter (e.g., λ = 5.5 ), its bifurcation point is located before the limit point, as shown in Figure 8b, leading to a bifurcation buckling mode of the arch instead of a limit point buckling mode. Accordingly, a critical geometric parameter λ b 2 is defined to identify the buckling mode of the arch, whose value is determined by solving Equations (27) and (34) at β = β b simultaneously. When the FG-GPLRC arch has λ > λ b 2 , the buckling mode of the arch is bifurcation mode.
The effects of the rotational constraints and the power law index on the critical geometric parameters to distinguish the buckling mode are presented in Figure 9. As shown in this figure, in Region 1 ( λ > λ b 2 ), the buckling mode of the arch is bifurcation mode. In Region 2 ( λ b 1 λ λ b 2 ), the buckling mode of the arch is either limit point instability mode or bifurcation mode, depending on which mode would occur first. In Region 3 ( λ s 1 λ < λ b 1 ), the buckling mode of the arch is limit point instability mode. In Region 4 ( λ < λ s 1 ), no buckling occurs in the arch. It should be mentioned that there is no solution to λ b 2 for α = 0 which refer to the fixed support, it implies that the FG-GPLRC arch with fixed ends cannot buckle in a bifurcation mode.
The nonlinear equilibrium paths of the FG-GPLRC arch, having the geometric parameters λ = λ s 2 = 4.0861 and λ = λ s 3 = 5.7552 , are given in Figure 10. The number of limit points of the arch, in which the flexibility of the rotational constraints is set as two, is indicated in the figure. Note that we can determine the number of limit points within the bounds of successive inflexion points. When λ < λ s 1 , there is no limit point in the nonlinear equilibrium path. When λ s 1 < λ < λ s 2 , there are two limit points in the nonlinear equilibrium path. When λ s 2 < λ < λ s 3 , there are four limit points in the nonlinear equilibrium path, and so on. Hence, the first 18 limit points and the corresponding critical geometric parameters λ s n of the arch with the power law index n = 1.5 and n = 3 are plotted in Figure 11 for further clarification.

6. Conclusions

In this paper, we studied the nonlinear buckling behavior of elastically supported FG-GPLRC arches with asymmetric distributed GPLs under a CPL. Analytical solutions for the limit point buckling, bifurcation buckling, and mode switching were derived. A finite element analysis was also employed to verify the present results, showing good accuracy of the present method for predicting the nonlinear buckling of the FG-GPLRC arch. According to the numerical studies, the influences of the distribution of GPLs, weight fractions, geometric configurations, as well as boundary restraints are discussed. It is found that the nonlinear buckling load of the FG-GPLRC arch decreases as the power law index or the flexibility of the rotational constraints increases. Switching of the buckling mode of the FG-GPLRC arch is quite sensitive to the power law index, the flexibility of rotational constraints, and the arch geometry. It is also found that, according to the present analytical method, the number of limit points in the nonlinear equilibrium path of the FG-GPLRC arch under a CPL can be easily determined.

Author Contributions

Conceptualization, J.F. and A.L.; Methodology, Z.Y., J.Y., and J.L.; Software, Z.Y. and J.L.; Validation, Z.Y., A.L., and J.Y.; Formal analysis, Z.Y. and S.-K.L.; Investigation, Z.Y., J.L., and S.-K.L.; Resources, A.L., J.F., and J.Y.; Data curation, Z.Y. and J.L.; Writing—original draft preparation, Z.Y.; Writing—review and editing, J.F., S.-K.L., and J.Y.; Visualization, Z.Y. and S.-K.L.; Supervision, J.F., A.L., and J.Y.; Project Administration, J.F. and A.L., J.Y., and S.-K.L.; Funding acquisition, J.F., A.L., Z.Y., J.Y., and S.-K.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Natural Science Foundation of China (no. 51925802, 11972123, and 51878188), the Technology Planning Project of Guangdong Province (no. 2020A1414010319), the Technology Planning Project of Guangzhou City (no. 201807010021), the General Research Fund (project no. PolyU 152008/19E) from the Research Grants Council of Hong Kong, and Discovery Project from the Australian Research Council (DP210103656). The authors are grateful for these supports.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data can be obtained from the first author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qiu, J.; Sahmani, S.; Safaei, B. On the NURBS-based isogeometric analysis for couple stress-based nonlinear instability of PFGM microplates. Mech. Based Des. Struc. 2020, 1–25. [Google Scholar] [CrossRef]
  2. Wu, H.; Yang, J.; Kitipornchai, S. Mechanical Analysis of Functionally Graded Porous Structures: A Review. Int. J. Struct. Stab. Dy. 2020, 20, 2041015. [Google Scholar] [CrossRef]
  3. Tho Hung, V.; Thuy Dong, D.; Thi Phuong, N.; Le Ly, N.; Quang Minh, T.; Trung, N.T.; Hoai Nam, V. Nonlinear buckling behavior of spiral corrugated sandwich FGM cylindrical shells surrounded by an elastic medium. Materials 2020, 13, 1984. [Google Scholar] [CrossRef]
  4. Ke, L.L.; Yang, J.; Kitipornchai, S. An analytical study on the nonlinear vibration of functionally graded beams. Meccanica 2010, 45, 743–752. [Google Scholar] [CrossRef]
  5. Librescu, L.; Oh, S.Y.; Song, O. Thin-walled beams made of functionally graded materials and operating in a high temperature environment: Vibration and stability. J. Therm. Stresses 2005, 28, 649–712. [Google Scholar] [CrossRef]
  6. Yan, T.; Kitipornchai, S.; Yang, J. Parametric instability of functionally graded beams with an open edge crack under axial pulsating excitation. Compos. Struct. 2011, 93, 1801–1808. [Google Scholar] [CrossRef]
  7. Yan, T.; Kitipornchai, S.; Yang, J.; He, X.Q. Dynamic behaviour of edge-cracked shear deformable functionally graded beams on an elastic foundation under a moving load. Compos. Struct. 2011, 93, 2992–3001. [Google Scholar] [CrossRef]
  8. Nguyen, H.N.; Tan, T.C.; Luat, D.T.; Phan, V.D.; Thom, D.V.; Minh, P.V. Research on the buckling behavior of functionally graded plates with stiffeners based on the third-order shear deformation theory. Materials 2019, 12, 1262. [Google Scholar] [CrossRef] [Green Version]
  9. Chen, C.S.; Tan, A.H. Imperfection sensitivity in the nonlinear vibration of initially stresses functionally graded plates. Compos. Struct. 2007, 78, 529–536. [Google Scholar] [CrossRef]
  10. Hao, Y.X.; Chen, L.H.; Zhang, W.; Lei, J.G. Nonlinear oscillations, bifurcations and chaos of functionally graded materials plate. J. Sound Vib. 2008, 312, 862–892. [Google Scholar] [CrossRef]
  11. Yang, J.; Chen, Y. Free vibration and buckling analyses of functionally graded beams with edge cracks. Compos. Struct. 2008, 83, 48–60. [Google Scholar] [CrossRef]
  12. Kitipornchai, S.; Ke, L.L.; Yang, J.; Xiang, Y. Nonlinear vibration of edge cracked functionally graded Timoshenko beams. J. Sound Vib. 2009, 324, 962–982. [Google Scholar] [CrossRef]
  13. Huang, C.S.; Huang, S.H. Analytical solutions based on fourier cosine series for the free vibrations of functionally graded material rectangular mindlin plates. Materials 2020, 13, 3820. [Google Scholar] [CrossRef]
  14. Bateni, M.; Eslami, M.R. Non-linear in-plane stability analysis of FGM circular shallow arches under central concentrated force. Int. J. Nonlin. Mech. 2014, 60, 58–69. [Google Scholar] [CrossRef]
  15. Bateni, M.; Eslami, M.R. Non-linear in-plane stability analysis of FG circular shallow arches under uniform radial pressure. Thin Wall Struct. 2015, 94, 302–313. [Google Scholar] [CrossRef]
  16. Czechowski, L.; Kołakowski, Z. Analysis of the functionally step-variable graded plate under in-plane compression. Materials 2019, 12, 4090. [Google Scholar] [CrossRef] [Green Version]
  17. Zaczynska, M.; Kazmierczyk, F. Multi-mode buckling analysis of FGM channel section beams. Materials 2020, 13, 2567. [Google Scholar] [CrossRef] [PubMed]
  18. Kim, H.; Abdala, A.A.; Macosko, C.W. Graphene/polymer nanocomposites. Macromolecules 2010, 43, 6515–6530. [Google Scholar] [CrossRef]
  19. Rafiee, M.A.; 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] [PubMed]
  20. Rafiee, M.A.; Rafiee, J.; Srivastava, I.; Wang, Z.; Song, H.; Yu, Z.Z.; Koratkar, N. Fracture and fatigue in graphene nanocomposites. Small 2010, 6, 179–183. [Google Scholar] [CrossRef] [PubMed]
  21. Qi, X.; Zhang, S.; Wang, T.; Guo, S.; Ren, R. Effect of High-Dispersible Graphene on the Strength and Durability of Cement Mortars. Materials 2021, 14, 915. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, J.; Yan, H.; Jiang, K. Mechanical properties of graphene platelet-reinforced alumina ceramic composites. Ceram. Int. 2013, 39, 6215–6221. [Google Scholar] [CrossRef]
  23. Tang, L.C.; Wan, Y.J.; Yan, D.; Pei, Y.B.; Zhao, L.; Li, Y.B.; Lai, G.Q. The effect of graphene dispersion on the mechanical properties of graphene/epoxy composites. Carbon 2013, 60, 16–27. [Google Scholar] [CrossRef]
  24. Zhao, S.; Zhao, Z.; Yang, Z.; Ke, L.; Kitipornchai, S.; Yang, J. Functionally graded graphene reinforced composite structures: A review. Eng. Struct. 2020, 10, 110339. [Google Scholar] [CrossRef]
  25. 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]
  26. Yang, Z.; Tam, M.; Zhang, Y.; Kitipornchai, S.; Lv, J.; Yang, J. Nonlinear dynamic response of FG graphene platelets reinforced composite beam with edge cracks in thermal environment. Int. J. Struct. Stab. Dy. 2020, 2043005. [Google Scholar] [CrossRef]
  27. Zhang, L.H.; Lai, S.K.; Wang, C.; Yang, J. DSC regularized Dirac-delta method for dynamic analysis of FG graphene platelet-reinforced porous beams on elastic foundation under a moving load. Compos. Struct. 2021, 255, 112865. [Google Scholar] [CrossRef]
  28. 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] [Green Version]
  29. Dong, Y.; Li, Y.; Li, X.; Yang, J. Active control of dynamic behaviors of graded graphene reinforced cylindrical shells with piezoelectric actuator/sensor layers. Appl. Math. Model. 2020, 82, 252–270. [Google Scholar] [CrossRef]
  30. Dong, Y.; Li, X.; Gao, K.; Li, Y.; Yang, J. Harmonic resonances of graphene-reinforced nonlinear cylindrical shells: Effects of spinning motion and thermal environment. Nonlinear Dynam. 2020, 99, 981–1000. [Google Scholar] [CrossRef]
  31. Song, M.; Kitipornchai, S.; Yang, J. Free and forced vibrations of functionally graded polymer composite plates reinforced with graphene nanoplatelets. Compos. Struct. 2017, 159, 579–588. [Google Scholar] [CrossRef]
  32. Wu, H.; Yang, J.; Kitipornchai, S. Parametric instability of thermo-mechanically loaded functionally graded graphene reinforced nanocomposite plates. Int. J. Mech. Sci. 2018, 135, 431–440. [Google Scholar] [CrossRef] [Green Version]
  33. Yang, Z.; Liu, A.; Yang, J.; Fu, J.; Yang, B. Dynamic buckling of functionally graded graphene nanoplatelets reinforced composite shallow arches under a step central point load. J. Sound Vib. 2020, 465, 115019. [Google Scholar] [CrossRef]
  34. Yang, Z.; Feng, C.; Yang, J.; Wang, Y.; Lv, J.; Liu, A.; Fu, J. Geometrically nonlinear buckling of graphene platelets reinforced dielectric composite (GPLRDC) arches with rotational end restraints. Aerosp. Sci. Technol. 2020, 107, 106326. [Google Scholar] [CrossRef]
  35. Yang, Z.; Zhao, S.; Yang, J.; Lv, J.; Liu, A.; Fu, J. In-plane and out-of-plane free vibrations of functionally graded composite arches with graphene reinforcements. Mech. Adv. Mater. Struc. 2020, 1–11. [Google Scholar] [CrossRef]
  36. Yang, Z.; Xu, J.; Lu, H.; Lv, J.; Liu, A.; Fu, J. Multiple equilibria and buckling of functionally graded graphene nanoplatelet-reinforced composite arches with pinned-fixed end. Crystals 2020, 10, 1003. [Google Scholar] [CrossRef]
  37. Huang, Y.; Yang, Z.; Liu, A.; Fu, J. Nonlinear buckling analysis of functionally graded graphene reinforced composite shallow arches with elastic rotational constraints under uniform radial load. Materials 2018, 11, 910. [Google Scholar] [CrossRef] [Green Version]
  38. Yang, Z.; Huang, Y.; Liu, A.; Fu, J.; Wu, D. Nonlinear in-plane buckling of fixed shallow functionally graded graphene reinforced composite arches subjected to mechanical and thermal loading. Appl. Math. Model. 2019, 70, 315–327. [Google Scholar] [CrossRef]
  39. Yang, Z.; Yang, J.; Liu, A.; Fu, J. Nonlinear in-plane instability of functionally graded multilayer graphene reinforced composite shallow arches. Compos. Struct. 2018, 204, 301–312. [Google Scholar] [CrossRef]
  40. Liu, Z.; Yang, C.; Gao, W.; Wu, D.; Li, G. Nonlinear behaviour and stability of functionally graded porous arches with graphene platelets reinforcements. Int. J. Eng. Sci. 2019, 137, 37–56. [Google Scholar] [CrossRef]
  41. Zhao, S.; Yang, Z.; Kitipornchai, S.; Yang, J. Dynamic instability of functionally graded porous arches reinforced by graphene platelets. Thin. Wall Struct. 2020, 147, 106491. [Google Scholar] [CrossRef]
  42. Li, Z.; Zheng, J.; Zhang, Z. Mechanics of the confined functionally graded porous arch reinforced by graphene platelets. Eng. Struct. 2019, 201, 109817. [Google Scholar] [CrossRef]
  43. Younis, M.I.; Ouakad, H.M.; Alsaleem, F.M.; Miles, R.; Cui, W. Nonlinear dynamics of MEMS arches under harmonic electrostatic actuation. J. Microelectromech. Syst. 2010, 19, 646–656. [Google Scholar] [CrossRef]
  44. Youn, J.-H.; Jeong, S.M.; Hwang, G.; Kim, H.; Hyeon, K.; Park, J.; Kyung, K.-U. Dielectric elastomer actuator for soft robotics applications and challenges. Appl. Sci. 2020, 10, 640. [Google Scholar] [CrossRef] [Green Version]
  45. Hosoya, N.; Baba, S.; Maeda, S. Hemispherical breathing mode speaker using a dielectric elastomer actuator. J. Acoust. Soc. Am. 2015, 138, EL424–El428. [Google Scholar] [CrossRef]
Figure 1. Configuration and coordinate system of an elastically supported functionally graded graphene platelet reinforced composite (FG-GPLRC) arch.
Figure 1. Configuration and coordinate system of an elastically supported functionally graded graphene platelet reinforced composite (FG-GPLRC) arch.
Materials 14 02026 g001
Figure 2. Effective Young’s modulus of an elastically restrained FG-GPLRC arch.
Figure 2. Effective Young’s modulus of an elastically restrained FG-GPLRC arch.
Materials 14 02026 g002
Figure 3. Mesh and deformation of an elastically restrained FG-GPLRC arch.
Figure 3. Mesh and deformation of an elastically restrained FG-GPLRC arch.
Materials 14 02026 g003
Figure 4. Comparison of FE results and analytical solutions.
Figure 4. Comparison of FE results and analytical solutions.
Materials 14 02026 g004
Figure 5. Effect of the power law index on (a) the buckling load and (b) the effective bending stiffness.
Figure 5. Effect of the power law index on (a) the buckling load and (b) the effective bending stiffness.
Materials 14 02026 g005
Figure 6. Effect of geometry and size of GPLs on the buckling load.
Figure 6. Effect of geometry and size of GPLs on the buckling load.
Materials 14 02026 g006
Figure 7. Effect of the flexibility of the rotational constraints on the buckling load.
Figure 7. Effect of the flexibility of the rotational constraints on the buckling load.
Materials 14 02026 g007
Figure 8. Nonlinear load-displacement curve of an elastically restrained FG-GPLRC arch.
Figure 8. Nonlinear load-displacement curve of an elastically restrained FG-GPLRC arch.
Materials 14 02026 g008
Figure 9. Critical geometric parameters for the buckling mode of an elastically restrained FG-GPLRC arch.
Figure 9. Critical geometric parameters for the buckling mode of an elastically restrained FG-GPLRC arch.
Materials 14 02026 g009
Figure 10. Nonlinear load-axial force curve of an elastically restrained FG-GPLRC arch.
Figure 10. Nonlinear load-axial force curve of an elastically restrained FG-GPLRC arch.
Materials 14 02026 g010
Figure 11. Number of limit points of an elastically restrained FG-GPLRC arch under different power law indices.
Figure 11. Number of limit points of an elastically restrained FG-GPLRC arch under different power law indices.
Materials 14 02026 g011
Table 1. Buckling loads of the FG-GPLRC arch having different weight fractions of distributed graphene platelets (GPLs).
Table 1. Buckling loads of the FG-GPLRC arch having different weight fractions of distributed graphene platelets (GPLs).
WGPL (%)n = 0n = 0.5n = 1n = 3
0.00.5088
0.21.01570.86790.83130.7846
0.41.52271.22131.14541.0503
0.62.03001.57301.45671.3123
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, Z.; Liu, A.; Yang, J.; Lai, S.-K.; Lv, J.; Fu, J. Analytical Prediction for Nonlinear Buckling of Elastically Supported FG-GPLRC Arches under a Central Point Load. Materials 2021, 14, 2026. https://doi.org/10.3390/ma14082026

AMA Style

Yang Z, Liu A, Yang J, Lai S-K, Lv J, Fu J. Analytical Prediction for Nonlinear Buckling of Elastically Supported FG-GPLRC Arches under a Central Point Load. Materials. 2021; 14(8):2026. https://doi.org/10.3390/ma14082026

Chicago/Turabian Style

Yang, Zhicheng, Airong Liu, Jie Yang, Siu-Kai Lai, Jiangen Lv, and Jiyang Fu. 2021. "Analytical Prediction for Nonlinear Buckling of Elastically Supported FG-GPLRC Arches under a Central Point Load" Materials 14, no. 8: 2026. https://doi.org/10.3390/ma14082026

APA Style

Yang, Z., Liu, A., Yang, J., Lai, S. -K., Lv, J., & Fu, J. (2021). Analytical Prediction for Nonlinear Buckling of Elastically Supported FG-GPLRC Arches under a Central Point Load. Materials, 14(8), 2026. https://doi.org/10.3390/ma14082026

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