An exact solution for the free-vibration analysis of functionally graded carbon-nanotube-reinforced composite beams with arbitrary boundary conditions

We present an exact method to model the free vibration of functionally graded carbon-nanotube-reinforced composite (FG-CNTRC) beams with arbitrary boundary conditions based on first-order shear deformation elasticity theory. Five types of carbon nanotube (CNT) distributions are considered. The distributions are either uniform or functionally graded and are assumed to be continuous through the thickness of the beams. The displacements and rotational components of the beams are expressed as a linear combination of the standard Fourier series and several supplementary functions. The formulation is derived using the modified Fourier series and solved using the strong-form solution and the weak-form solution (i.e., the Rayleigh–Ritz method). Both solutions are applicable to various combinations of boundary constraints, including classical boundary conditions and elastic-supported boundary conditions. The accuracy, efficiency and validity of the two solutions presented are demonstrated via comparison with published results. A parametric study is conducted on the influence of several key parameters, namely, the L/h ratio, CNT volume fraction, CNT distribution, boundary spring stiffness and shear correction factor, on the free vibration of FG-CNTRC beams.

Carbon nanotube (CNT)-reinforced composites have shown outstanding physical, mechanical, thermal and electrical properties over traditional structural materials, drawing interest from numerous researchers. CNTs are recognized as well suited to reinforce polymer composites due to their high elastic modulus and tensile strength and low density [1][2][3][4] . As a result, enthusiasm for research activities involving CNTs has been ignited in recent years. Wagner et al. 5 performed tensile experiments on multi-walled carbon nanotubes (MWCNTs) and analysed the transformation of the elastic modulus and the break stress of nanoscale reinforced composites. Qian et al. 6 developed homologous research two years later, noting that the addition of only 1% MWCNT to polystyrene significantly improved its polymeric mechanical properties. Fiedler et al. 7 demonstrated the superiority of the CNTs as nanofillers in polymers and suggested that a distribution of CNTs should be concentrated or dispersed to realize the best possible properties. Han 8 and Wan 9 found that the introduction of low volume fractions of nanotubes in matrices can result in notable strengthening of the composite properties. Relative to micron-scale counterparts, the interfacial regions between the nanoparticles and the matrix are strongly reactive. Coleman et al. 10 compared mechanical properties and various manufacturing processes of single-walled carbon nanotube (SWCNT)-reinforced composites to MWCNT-reinforced composites. However, the collective problems of dispersion and stress transfer still lack solutions.
To overcome these issues, functionalization through chemical procedures, for example, has been adopted by researchers. Functionally graded (FG) materials, which act as primitive thermal barrier materials in the aerospace industry, are widely known for their smooth and continuous variations in material properties. In this way, the integration of single materials with different properties can be improved, and the advantages of material properties can be combined [11][12][13] . Shen 14 in 2009 first proposed a new distribution form with CNTs distributed in an FG manner in the matrix; the volume fraction of CNTs was assumed to vary along the thickness direction. The issue of nonlinear bending behaviour of functionally graded carbon-nanotube-reinforced composite (FG-CNTRC) plates was investigated in which a transverse uniform or sinusoidal load in thermal environments was taken into consideration. The results revealed that mechanical, electrical and thermal properties varied considerably with the FG distribution of CNTs.
CNT-reinforced composites, which are increasingly used, can be formed into structures such as beams, plates and shells. Beams are a fundamental and significant structure in comprehensive engineering applications in the fields of marine, aerospace, civil, and mechanical engineering, among others areas. In this regard, various studies have focused on dynamic characteristics analysis of beam structures to guide structurally reliable design of such engineering applications [15][16][17][18][19][20] . Incorporating first-order beam theory, Yas and Samadi 21 adopted the generalized differential quadrature method to analyse the issue of vibrations and buckling of carbon-nanotube-reinforced composite (CNTRC) beams on elastic foundations. Four different CNT distributions were considered, and the material properties of the nanocomposites were obtained from the rule of mixtures. In light of the von Kármán geometric nonlinearity displacement-strain relationship and Euler beam theory, the linear and nonlinear vibration behaviours of CNT-reinforced FG composite beams were presented by Rafiee et al. 22 . Numerical results showed that an increase in the CNT volume fraction led to an increase in the nonlinear-to-linear frequency ratio and the natural frequencies. Lin and Xiang 23,24 investigated the free-vibration characteristics of CNT-reinforced beams with soft-clamped and hard-clamped boundary conditions, with uniformly distributed (UD) CNTs and FG distribution being considered. The model was established based on first-order and third-order shear deformation elasticity theories and solved using the polynomial Ritz method. The research showed that ratios of nonlinear-to-linear frequency parameters and natural frequencies based on first-order and third-order shear deformation elasticity theories with soft-clamped boundary conditions showed manifest deviations. Ke et al. 25 investigated the nonlinear vibration characteristics of FG-CNTRC beams with various boundary conditions using a direct iterative method. The influences of the vibration amplitude, volume fraction of CNTs, ratio of length to thickness, boundary conditions and CNT distribution were taken into account to characterize nonlinear vibration in the beams. The response of CNT-reinforced FG composite beams under low-velocity impact was first analysed by Jam and Kiani 26 . On the basis of first-order beam theory, the behaviour of FG-CNTRC beams exposed to the impact of a small mass was solved by means of the conventional polynomial Ritz method and the Runge-Kutta method. The peak contact force was found to be proportionate to the volume fraction of CNTs and inversely proportional to the temperature, while the contact time behaved oppositely.
With the rapidly increasing industrial use of composite materials, various numerical tools and theories have been promoted to analyse the mechanical behaviour of composite structures [27][28][29][30] . The problem of nonlinear vibration of composite plates reinforced by CNTs was presented by Wang 31 . The governing equation of the CNTRC plate was derived according to higher-order shear deformation theory, and the theoretical model was solved using the improved perturbation technique. Zhang et al. 32 focused attention on the free-vibration characteristics of CNT-reinforced FG composite triangular plates and adopted the element-free IMLS-Ritz method based on first-order beam theory. In view of the first order shear deformation theory, Rafiee et al. 33 employed the Galerkin method and the harmonic balance method to investigate initially imperfect piezoelectric composite plates reinforced by SWCNTs. The vibrational characteristics and buckling of an FGM microplate with two different supports were studied by Ke et al. 34 . Wang 35,36 developed a unified semi-analytical approach and applied it to the issue of free-vibration analysis of FG-CNTRC structures of revolution, including spherical panels and doubly curved shells. The differential quadrature method and the Mindlin plate theory were applied in this research to enable scientific conclusions to be drawn. Based on the FEM and two types of shear deformation theory, Yas and Heshmati 37 established an analytical model of a FG-CNTRC beam subjected to a moving load.
Numerous studies have been conducted to illustrate the vibrational characteristics of CNTRC beams. Nevertheless, the investigations mentioned above are limited to several representative boundary conditions. A diversity of boundary restraints leads to notable changes in the free-vibration characteristics. Relatively little study has addressed the free vibration of CNT-reinforced beams with elastic supports, although various possibilities of boundary conditions appear in engineering practice. Furthermore, most actual solution procedures are customized to restricted forms of certain classes at both ends of the beam. Consequently, existent contributions are urgently sought not only to guide engineering applications but also to enhance complementary research.
Aiming at satisfying practical needs, in this investigation, a unified and satisfactorily accurate method is presented for the free-vibration analysis of FG-CNTRC beams with arbitrary boundary conditions, including various classical boundary conditions and elastic supports. The modified Fourier method was first proposed by Li to analyse the vibration of a beam 38 and was subsequently extended to plates [39][40][41][42][43][44][45][46][47] and shells [48][49][50][51][52][53] . The two functions of displacement and two functions of rotation are expressed as a linear combination of an original Fourier cosine series expansion and two complementary auxiliary polynomial functions. The supplemental items are introduced to remove the potential discontinuities of displacement components and derivatives of displacement functions at the ends of the beam and to accelerate convergence of the solution procedures. Arbitrary boundary conditions can be conveniently achieved through assigning the appropriate stiffness to four sets of boundary springs at each edge of the CNTRC beam without updating the solution procedure.
In the present work, a strong-form solution procedure of the modified Fourier method is proposed and used to solve generalized eigenvalue problems directly by submitting a modified Fourier series to the governing equations and the boundary conditions. In addition, the results obtained from the Rayleigh-Ritz technique associated with the modified Fourier method are presented here as a weak-form solution for comparison. Numerical results calculated by the present method are checked against available results published in the open literature to evaluate accuracy and validity. New results in terms of frequency parameters and mode shapes of CNTRC beams with elastic boundary conditions are presented here to provide a benchmark for future researchers in this field.
Scientific REPORTS | 7: 12909 | DOI:10.1038/s41598-017-12596-w Theory and formulation. As depicted in Fig. 1, we consider a general CNT-reinforced composite beam with a rectangular cross section. A right-handed Cartesian coordinate system is established in which the length L, width b, and thickness h are respectively defined along the x-, y-and z-directions. The main purpose of this work is to investigate the CNTRC beam with arbitrary boundary conditions, and thus, two sets of linear springs (k L u 0, and k L w 0, ) and two sets of rotational springs (K L s 0, and K L c 0, ) are artificially introduced to simulate boundary restraint forces at the two ends of the beam. By assigning a suitable stiffness to the four sets of boundary springs, an arbitrary combination of classical and elastic boundary conditions can be realized. For example, if linear and rotational restraining spring coefficients at both ends are set to infinity (or a sufficiently large number in practical numerical simulations), the perfectly clamped boundary condition can be conveniently achieved.
It is assumed that the CNTRC beam consists of a polymer matrix mixture that can be generally treated as an isotropic material and CNTs. SWCNT reinforcements are placed along the length direction and are either UD or FG in the thickness direction. Four-types of FG distribution forms are taken into account, namely, FG-X, FG-O, FG-V and FG-Λ, as illustrated in Fig. 2.
Regardless of the various distribution patterns of CNT reinforcement at the cross section, four types of FG-CNTRC beams are assumed to contain an equal CNT total weight of m tcnt and total CNT volume fraction V tcnt . The expressions of CNT volume fraction for different distributions can be written as: The effective material properties of CNTRC beams, e.g., the Young's modulus (E 11 and E 22 ), shear modulus (G 12 ) and Poisson's ratios (v 12 and v 21 ), are estimated according to matching molecular dynamics simulation results based on the rule of mixtures, which can be expressed as follows: where E cnt 11 , E cnt 22 and G cnt 12 indicate the Young's modulus along the longitudinal direction, Young's modulus in the transverse direction and shear modulus of CNTs, respectively. E m and G m indicate the Young's modulus and shear modulus of the isotropic matrix, respectively. υ cnt 12 and ρ cnt represent the Poisson's ratio and mass density of CNT, respectively, and v m and ρ m are the corresponding properties of the matrix. V m represents the matrix volume fraction, and η j (j = 1, 2, 3) denotes the CNT efficiency parameters, which are determined from the results of molecular dynamics simulations.
The displacement field for the CNTRC beam under the assumptions of first-order shear deformation elasticity theory can be expressed as follows: where u and w are the axial and transverse displacements along the x-and z-directions in the middle surface, respectively, and θ and φ represent the rotations of the normal to the section about the y-and x-axes, respectively. The strain and curvatures are defined in terms of the mid-plane displacements and rotations as: The constitutive equations are given by:  55 55 in which ε x and γ xz denote the normal and shear strain, respectively. ε x 0 , ε y 0 and γ xy represent the strain at the middle surface; k x , k y and k xy are the bending and twisting curvatures; N x , N y and N xy indicate the force resultants at the middle surface; M x , M y and M xy are the bending and twisting moment resultants; and Q xy and Q yz represent the shear force resultants.
Regarding the CNTRC beam, certain force and moment resultants, namely, N y , N xy , Q yz and M y , are equal to zero, while the corresponding strains ε y 0 , γ xy and curvature k y are assumed to be non-zero. Consequently, Eq. (5) can be expressed as:  The extensional stiffness coefficients A ij , coupling stiffness coefficients B ij , bending stiffness coefficients D ij (i, j = 1, 2, 6) and transverse shear stiffness A 55 are defined as functions of material properties, which can be written as: where κ denotes the shear correction factor. The reduced stiffness coefficients = Q i ( 1,2,6) ij are defined by the following equations:  16 26 With the aim of deriving the governing equations and boundary conditions of the CNTRC beam according to Hamilton's principle, the energy expressions are defined as follows. The total linear elastic strain energy (U s ) function of the CNTRC beam can be expressed as in which the inertia terms can be written as In addition, four sets of boundary springs are introduced at each end of the beam; the boundary springs deformation strain energy (U sp ) function is given by and the Hamilton's principle with regard to the arbitrary initial time t 1 and final time t 2 is given by Substituting Eqs (10), (11) and (13) into Eq. (14) and integrating by parts to eliminate the variational terms, the governing equations and boundary conditions can be obtained as follows: Because of the arbitrariness of the virtual displacements δu and δw, only when the values of the virtual displacements coefficients are equal to zero is Eq. (15) tractable. The governing equations can be expressed in terms of the differentials of displacement components as  and the general boundary conditions can be stated as Therefore, all classical boundary conditions and elastic supports can be directly achieved by means of the artificial spring boundary technique to assign the rigidities of the boundary springs at a certain value. The appropriate choice of admissible displacement functions plays a significant role to ensure the validity and accuracy of the proposed solution procedures. Eqs (16a-16d) indicate that the translation and rotations displacements of the CNTRC beam are required up to the second derivative. For the sake of satisfying the arbitrary boundary conditions at both ends of the beam, the displacements and rotational components are represented as 1D Fourier cosine series expansions with two supplemental auxiliary function terms. These terms are introduced to improve the convergence of the primary Fourier series representations and avoid the potential discontinuities of the displacement functions and their first-order derivatives at the boundaries. Accordingly, the functions of flexural displacements and rotation of the CNTRC beam can be universally expressed as . M represents the truncation number, A m , B m , C m and D m are the expansion coefficients of the standard Fourier series, and a i , b i and c i (i = 1, 2) denote the corresponding expansion coefficients of the auxiliary function P 1 (x) and P 2 (x), which are defined as: According to the modified Fourier series, the free-vibration characteristics of the CNTRC beam can be solved by means of strong-form solution procedures and weak-form solutions, as described below. The strong-form solution procedure is given step by step as follows.
We rewrite Eq. (18) in matrix form as Substituting Eqs (20) and (21) into Eq. (16), we obtain    Similarly, substituting Eqs (20) and (21)    Therefore, the expansion coefficients of the standard Fourier cosine series and corresponding auxiliary functions have a certain relationship according to the boundary conditions, which are expressed as: Substituting Eqs (26) and (27) into Eq. (22), then multiplying the transpose of displacement functions matrix H f on the left side and integrating both sides of the equality from 0 to L with respect to x, the partial differential equations are transformed into a standard eigenvalue problem as follows: in which K and M are the stiffness matrix and mass matrix, respectively, G is a vector that contains all undetermined coefficients of the standard Fourier series, and these matrices can be written as: The natural frequencies and modes of CNTRC beams can be obtained directly by solving the standard eigenvalue equation.
Exact solutions are often unavailable in complex vibration problems, and an approximate method is employed to complete the vibrational analysis. The Rayleigh-Ritz method associated with modified Fourier series, i.e., weak-form solution, is also presented below to compare with the strong-form solution.
In the Ritz-variational energy procedure, the accuracy of the solution will rest on how well the actual displacement can be faithfully represented by an appropriate admissible displacement field in general. Hence, auxiliary functions play a crucial role. The same displacement functions are selected for comparative purposes, and all expansion coefficients of the modified Fourier series can be regarded as generalized coordinates independently and equally.
With regard to free-vibration analysis, the Lagrange energy function of CNTRC beams consists of the strain energy, kinetic energy and boundary spring deformation strain energy as follows: Substituting Eqs (10), (11) and (13) into Eq. (30), minimizing the total expression of the Lagrange energy function via taking the derivatives of the equation with respect to the generalized coordinates and setting all expressions equal to zero to find the stationary value of the energy function, we obtain: A total of 4(M + 3) linear algebraic equations for the undetermined coefficients are achieved, which can be added and represented in a matrix form as where G * indicates the undetermined coefficients column vector, K * is treated as the total stiffness matrix of the CNTRC beams and M * indicates the corresponding mass matrix. Their expressions can be written as:  In addition, the non-dimensional frequency parameters of the natural frequency take the form of in the latter subsections unless otherwise stated. The corresponding boundary conditions at the ends of beam can be defined in terms of the spring stiffness as:

Elastically restrained case 1 (E1):
Elastically restrained case 2 (E2): The rationality of these definitions of boundary conditions in terms of assigning spring stiffness will be established through numerical examples in subsequent studies. For brevity, symbolism is applied to illustrate the boundary condition of FG-CNTRC beams, e.g., SE1 indicates a beam with S (simply supported) and E1 (elastically restrained case 1) boundary conditions at = 0 and x = L, respectively.

Convergence and validation.
As previously mentioned, modified Fourier series with infinite terms in the current solution framework are infinitely approximate in the real results. Nevertheless, the infinite terms must be numerically truncated in practical numerical simulations. Consequently, convergence studies are conducted to determine the number of series terms M used in the computation. The first four lowest frequency parameters Ω for perfectly clamped and simply supported FGV-CNT beams are considered in Table 1, in which the results obtained from strong-form solution procedures and the Rayleigh-Ritz method are given. Excellent convergence and satisfactory numerical stability of two types of current solutions can be observed. The frequency parameters Ω converge sharply as the number of series terms M increases from 4 to 15, and the results are almost invariant when the truncated number reaches a certain value (M = 10). Thus, unless otherwise illustrated, the truncated number was uniformly chosen as 10 in subsequent studies.
With satisfactory results for the convergence studies of the FG-CNTRC beam, which is assumed to feature a perfectly clamped boundary condition at both ends, the numerical validity and rationality of the mentioned definition of the boundary conditions in terms of assigning boundary spring rigidities is evaluated in this section. The first three non-dimensional frequency parameters for FG-CNTRC beams with various CNT distributions and boundary conditions are compared with those reported in the publications of Lin et al. 20 and Yas et al. 18 , as shown in Table 2 and Table 3.
The geometrical and material constants of the beams are provided as follows: GPa, E m = 2.5 GPa, υ = .  Table 2 presents the first three non-dimensional frequency parameters for UD-CNT, FGV-CNT and FGX-CNT beams with a total volume fraction V tcnt = 0.28, and two classical boundary conditions are considered, namely, S-S and C-F. The solutions from the two present numerical approaches are in outstanding agreement with the results from Lin et al. 20 and Yas et al. 18 .
The first three dimensionless natural frequencies for FG-CNTRC beams with C-S boundary conditions are presented in Table 3. The results in the present investigation are close to those in the references, and the two present numerical approaches are sufficiently accurate to enable vibrational characterization of FG-CNTRC beams subject to various boundary conditions. Moreover, a further comparison is explored to illustrate the applicability of the linear theories and assumptions in this investigation. Table 4 demonstrates the comparison of natural frequencies between the results calculated by the present method and those reported in studies 54,55 that adopted first-order beam theory along with von Karman geometric nonlinearity.
Scientific REPORTS | 7: 12909 | DOI:10.1038/s41598-017-12596-w The first three frequency parameters Ω of simply supported FG-CNTRC beams with different volume fractions are presented in Table 4. Consistency can be observed between the calculated results and the data from the literature. Thus, the mentioned numerical examples indicate that the current solutions possess rapid convergence  Table 3. Comparison of the first three dimensionless frequencies Parametric studies. Tables 2-5 indicate the accuracy and convergence of the present solution. With enhanced confidence in the present solution approach, a variety of further results for FG-CNTRC beams with different boundary conditions and material and geometrical parameters are provided in this section to serve as benchmark solutions for potential studies. In addition, the frequencies calculated by the strong-form solution procedures are close to the results obtained using the Rayleigh-Ritz method. For brevity, only the results obtained from the strong-form solution procedures are included. Figure 3 illustrates the relationship between the first three frequency parameters Ω and the length-to-thickness ratio L/h for FG-CNTRC beams with various boundary conditions and CNT distributions. The total CNT volume fraction V tcnt is fixed at 0.17, and the length-to-thickness ratio varies from 0.5 to 4. All dimensionless frequencies increase with increasing length-to-thickness ratio. The fundamental parameter rises gradually with increasing L/h ratio; in contrast, the second and third frequency parameters increase sharply. Furthermore, the graphs provide notable results regarding the influence of CNT distributions. The frequency parameters of FG-XCNT beams are always larger than the results of other distributions and beams, with FGO-CNT distributions having the smallest values regardless of the boundary conditions. The same phenomenon can also be found in the following tables. Table 5 shows the changes in the dimensionless frequencies of FG-CNTRC beams with various CNT distributions. The data in Table 5 lead us to conclude that the CNT distributions have a significant impact on the free-vibration characteristics of FG-CNTRC beams. Furthermore, symmetrical CNT distributions, i.e., FG-XCNT and FG-OCNT, play a notable role in changing the frequency parameters of the CNTRC beams relative to the uniform and asymmetric distributions through the beam thickness.
One of the primary purposes of this work is to investigate the free-vibration characteristics of FG-CNTRC beams with elastic boundary constraints. Accordingly, Fig. 4 illustrates the effects of four types of boundary spring parameters on the first three frequency parameters of FG-CNTRC beams with elastic supports. The total CNT volume fraction V tcnt and ratio L/h are 0.28 and 10, respectively. The UD-CNT, FGΛ-CNT and FGX-CNT beams are taken into consideration. The symbols Γ u , Γ w , Γ θ and Γ ϕ are defined to indicate the various kinds of boundary springs, and the boundary condition is defined as elastically restrained only at x = 0, where only one  Table 5. First three frequency parameters Ω for FG-CNTRC beams with various CNT distributions (L/h = 10, V tcnt = 0.12).
group of boundary springs is assigned to variable stiffness values ranging from 10 2 to 10 16 and the other groups are assumed to be infinite while the other boundary is clamped at x = L. Figure 4 shows that the dimensionless frequency parameters remain stable as the restraint parameters Γ u and Γ ϕ change. In contrast, the dimensionless frequencies increase sharply as the restraint parameters Γ w and Γ θ increase from 10 6 to 10 10 . Furthermore, there is little variation in the frequency parameters beyond this range. In addition, the results confirm that the influence of the restraint parameter Г w is more easily detectable than that of Г θ . Therefore, the definition of the boundary conditions mentioned above in terms of assigning the values of boundary spring stiffness is reasonable to simulate the real restraints.
The next example is focused on the influence of the total CNT volume fraction on FG-CNTRC beams with classical and elastic boundary conditions. The material properties of the beams are provided as follows:  Table 6 indicates the changes in the value of the first three dimensionless frequencies of FG-CNTRC beams as the total CNT volume fraction V tcnt increases from 0.12 to 0.28. The frequency parameters uniformly increase as the total CNT volume fraction increases.
By introducing the shear correction factor κ, the first-order shear deformation elasticity theories address the shortcomings of the Euler beam theory, which neglects the effects of transverse shear and rotary inertia. Note that all results in this study are based on the first-order beam theory and that it is necessary to study the influence of the shear correction factor on the free-vibration characteristics of FG-CNTRC beams. Table 7 presents the fundamental frequency parameters Ω for FG-CNTRC beams in the case in which the shear correction factor increases from 0.1 to 0.9 and compares the results with those calculated by Lin 21 based on Figure 3. The first three lowest frequency parameters Ω of FG-CNTRC beams with various length-to-thickness ratios. The total CNT volume fraction V tcnt is equal to 0.17, and the length-to-thickness ratio varies from 0.5 to 4. Three boundary conditions are considered: (a) F-F, (b) C-C, and (c) C-S.
Scientific REPORTS | 7: 12909 | DOI:10.1038/s41598-017-12596-w third-order shear deformation elasticity theory. The beam material properties and geometrical parameters are the same as for the FG-CNTRC beams presented in Table 5, where two classical boundary conditions and three types of CNT distribution are considered. The figures reveal that the frequency parameters monotonically increase as the shear correction factor increases from 0.1 to 0.9. To make the frequency parameters consistent with the results based on third-order beam theory, the appropriate shear correction factor κ should be selected in calculations with regard to the different boundary conditions.
Because the free-vibration results for FG-CNTRC beams with arbitrary boundary conditions are extremely limited in the literature, new results are calculated in Table 8 to provide reference data for practising engineers and to act as a benchmark for potential future studies. Finally, aiming at strengthening our understanding of vibration behaviours of FG-CNTRC beams, several selected mode shapes of the beams addressed in Table 8 are plotted in Fig. 5.

Conclusions
An accurate method is developed for the vibration analysis of FG-CNTRC beams. The distribution of CNTs through the thickness of the beam is assumed to vary continuously and smoothly, and five types of distribution, namely, UD-CNT, FGΛ-CNT, FGV-CNT, FGO-CNT and FGX-CNT, are considered. Note that this approach can be uniformly and conveniently applied in vibrational analysis of FG-CNTRC beams with arbitrary boundary conditions, including general elastic boundary conditions. The general boundary conditions can be enforced using the artificial spring technique, in which boundary springs can be assigned any value of stiffness to simulate the real boundary conditions. The energy expressions of the FG-CNTRC beams are written as functions of four displacement components based on first-order shear deformation elasticity theory. Regardless of the boundary conditions, specific geometry and material properties, the displacements and rotational components of the beam are expressed as a superposition of the standard cosine Fourier series and two auxiliary functions. The introduced auxiliary terms are intended to remove potential discontinuous displacement functions and their derivatives at each edge and to ensure the convergence of the series expansions.
By submitting modified Fourier series to governing equations and boundary conditions, the strong-form solution procedure of the modified Fourier method is proposed. For comparison, the Rayleigh-Ritz technique associated with the modified Fourier method is also presented as a weak-form solution. Numerical results obtained    Table 8. Fundamental frequency parameters Ω for a FG-CNTRC beam with various boundary conditions. by these two present methods are compared with the available results previously reported, and both accuracy and satisfactory convergence are observed. The free-vibration characteristics of FG-CNTRC beams are analysed with a variety of key parameters, for example, the L/h ratio, CNT volume fraction, CNT distribution, boundary spring stiffness and shear correction factor. New vibration results containing frequency parameters and mode shapes for the FG-CNTRC beams with classical boundary conditions and elastic supports are calculated to provide reference data for practising engineers and act as a benchmark for future studies.
Data availability statement. All data generated or analysed during this study are included in this published article.