Buckling Analysis of Functionally Graded Tapered Microbeams via Rayleigh–Ritz Method

Abstract

In the present study, the buckling problem of nonhomogeneous microbeams with a variable cross-section is analyzed. The microcolumn considered in this study is made of functionally graded materials in the longitudinal direction and the cross-section of the microcolumn varies continuously throughout the axial direction. The Bernoulli–Euler beam theory in conjunction with modified strain gradient theory are employed to model the structure by considering the size effect. The Rayleigh–Ritz numerical solution method is used to solve the eigenvalue problem for various conditions. The influences of changes in the cross-section and Young’s modulus, size dependency, and non-classical boundary conditions are examined in detail. It is observed that the size effect becomes more pronounced for smaller sizes and differences between the classical and non-classical buckling loads increase by increasing the taper ratios.

1. Introduction

Engineering is a profession in which the principles of nature are applied to create useful objects. A mechanical engineer develops a new engine. A mechatronic engineer designs a robot. Civil engineers design a bridge, a dam, or a building. An electronic engineer designs a computer or an integrated circuit. For many reasons, an engineer is not only interested in a design that works at a nominal level but also one which is somehow the best design. Simply put, the process of determining the best design is called optimization. Consequently, we may wish to design the lowest-cost bridge for the site.

On the other hand, it is not always possible or economical to design and build structures/prismatic structures with one-type (homogeneous) of, or a constant, cross-section. For this reason, it is useful to use two or more materials with variable cross-sections in many engineering applications.

Functionally graded materials (FGMs) can be described as relatively new types of composite materials that are used in various engineering applications, such as turbine blades, cutting tools, rocket bodies, engine cylinders, artificial bones, and dental implants. Due to the gradual and smooth variation of their material properties, FGMs have some desired properties, for example, in avoiding cracking and delamination problems and diminishing undesired stress concentrations between laminates. These improved materials were introduced by a group of Japanese scientists in 1984 to propose a thermal barrier material that can withstand a surface temperature of 2000 °K and a temperature range of 1000 °K for a section thinner than 10 mm in aircrafts and space shuttles [1].

With rapid developments in technology, FGMs have recently been applied in micro-/nano-electromechanical systems (MEMS/NEMS), for example, as a structural element in shape memory alloy films [2], electrically actuated MEMS devices [3], and atomic force microscopes [4]. The characteristic dimensions of these microstructures are specified in micrometers and/or nanometers. As has been demonstrated through experimentation, the size effect phenomenon plays a key role in the mechanical characteristics of microstructures [5,6,7,8,9,10,11,12]. Experimental studies to determine the static and dynamic behaviors of structural elements are very important and valid. However, testing these structures under different conditions is very laborious, expensive, and not feasible, especially at the micro and nanoscales.

Solid mechanics is a branch of continuum mechanics that deals with the static and dynamic responses of structures by modelling them mathematically. The mathematical modelling of structures can be simpler, cheaper, and more useful than experimental works. Consequently, continuum mechanics may be considered as a potential approach to determine the mechanical behaviors of small-sized structures. Unfortunately, classical continuum theory fails to estimate the responses of microstructures due to inability to consider the size effect phenomena. Thus, the utilization of size-dependent theories containing at least one material length scale parameter is required. Several higher-order elasticity theories have been proposed such as couple stress, nonlocal elasticity, strain gradient, nonlocal strain gradient, and doublet mechanics theories [13,14]. These theories have been frequently employed to analyze microstructures in the past [15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30].

Free vibration and buckling behaviors of axially functionally graded (AFG) tapered beams were numerically investigated by utilizing differential quadrature [31], differential transform [32], and direct integral methods [33]. The transverse vibration of functionally graded thick beams was examined using the Timoshenko beam theory [34] and a refined beam theory considering thickness-stretching effect [35]. The interpolation matrix method was applied to obtain the critical buckling loads of AFG-tapered thick beams based on the first-order shear deformation beam theory [36]. A general solution for the free vibration of AFG-tapered cantilevers with a tip mass was presented based using the Bernoulli–Euler beam theory and the Myklestad method [37]. The thermo-mechanical buckling of bi-directional functionally graded porous microbeams was studied [38]. Material properties were changed throughout the height and length of the microbeam according to the power law. The free vibrational response of non-prismatic AFG Timoshenko beams were subjected to thermal variation in humid environments and examined via the harmonic differential quadrature method [39].

Size-dependent stability analysis of homogeneous microbeams with varying cross-sections was carried out based on the Bernoulli–Euler beam theory and modified strain gradient theory [40]. Moreover, the buckling response of Bernoulli–Euler AFG-tapered microstructures was examined using a modified couple stress theory for different boundary conditions [41]. Nonlocal strain gradient theory was used to model the buckling response of AFG nanobeams lying on a variable elastic foundation [42]. An analytical solution for the static and dynamic responses of functionally graded tapered microbridges was introduced on the basis of modified couple stresses and three-dimensional beam theories [43]. Stability and free vibration responses of tapered functionally graded composite microbeams reinforced by carbon nanotube were subjected to a magnetic field and investigated using the finite element method [44]. The free vibration of viscoelastic tapered microbars embedded in a visco-Pasternak foundation was studied based on strain gradient theory and the Kelvin–Voigt model [45]. Size-dependent nonlinear vibrations of functionally graded microbeams with porosities was carried out based on modified couple stress and Bernoulli–Euler beam theories. The numerical results were obtained by utilizing generalized differential quadrature and direct iterative methods [46]. A buckling analysis of bi-directional functionally graded tapered Bernoulli–Euler porous nanobeams was performed via nonlocal elasticity theory [47]. Moreover, the thermo-mechanical vibrational response of rotary functionally graded tapered Bernoulli–Euler nanobeams was perused for cantilever and propped cantilever beams on the basis of nonlocal elasticity theory [48].

From this literature survey, it can be seen that there are many studies on the linear and nonlinear mechanical behaviors of functionally graded tapered micro and nanobeams based on nonlocal elasticity, modified couple stress, and nonlocal strain gradient theories in particular. To the best of the authors’ knowledge, there is no study on the size-dependent buckling analysis of axially functionally graded tapered microstructures based on the modified strain gradient theory by considering higher-order boundary conditions and by possessing variable length scale parameters. In this work, microstructure-dependent buckling of AFG-tapered microbeams is examined. Euler-Bernoulli beam and modified strain gradient theories are utilized to formulate the AFG-tapered microbeams. The Rayleigh–Ritz method is applied to obtain the critical buckling loads of cantilever and propped cantilever microbeams at various taper ratios and gradient indices. Effects of taper ratios, material property gradient indices, length scale parameters, and nonclassical boundary conditions are investigated in detail.

2. Theory and Formulation

Modified strain gradient theory (SGT) is introduced by Lam et al. [5]. This nonclassical elasticity theory includes three additional material length scale parameters to take into consideration when investigating size effect. According to this theory, the strain energy $U$ in a linear elastic isotropic material can be expressed as

$$U=\frac{1}{2}{{\displaystyle \int }}_0^L{{\displaystyle \int }}_A\left({\sigma }_{ij}{\epsilon }_{ij}+p_i{\gamma }_i+{\tau }_{ijk}^{\left(1\right)}{\eta }_{ijk}^{\left(1\right)}+m_{ij}^s{\chi }_{ij}^s\right)dAdx$$

(1)

$${\epsilon }_{ij}=\frac{1}{2}\left(u_{i,j}+u_{j,i}\right)$$

(2)

$${\epsilon }_{mm,i}={\gamma }_i$$

(3)

$$\begin{array}{c}{\eta }_{ijk}^{\left(1\right)}=\frac{1}{3}\left({\epsilon }_{jk,i}+{\epsilon }_{ki,j}+{\epsilon }_{ij,k}\right)-\frac{1}{15}[({\delta }_{ij}\left({\epsilon }_{mm,k}+2{\epsilon }_{mk,m}\right)\\ +{\delta }_{jk}\left({\epsilon }_{mm,i}+2{\epsilon }_{mi,m}\right)+{\delta }_{ki}\left({\epsilon }_{mm,j}+2{\epsilon }_{mj,m}\right))]\end{array}$$

(4)

$${\chi }_{ij}^{\left(s\right)}=\frac{1}{2}\left({\theta }_{i,j}+{\theta }_{j,i}\right)$$

(5)

$${\theta }_i=\frac{1}{2}{e}_{ijk}{u}_{k,j}$$

(6)

where $u_i$, ${\theta }_i$, ${\epsilon }_{ij}$, ${\gamma }_i$, ${\eta }_{ijk}^{\left(1\right)}$, and ${\chi }_{ij}^s$ represent the components of the displacement and rotation vectors, the classical strain tensor, the dilatation gradient vector, the deviatoric stretch gradient, and the symmetric rotation gradient tensors, respectively. $\delta$ is Kronecker delta and $e_{ijk}$ is an alternate symbol. Furthermore, the components of the classical stress tensor and the higher-order stress tensors can be defined as [5]

$${\sigma }_{ij}=\lambda {\epsilon }_{mm}{\delta }_{ij}+2G{\epsilon }_{ij}$$

(7)

$$p_i=2G{l}_0^2{\gamma }_i$$

(8)

$${\tau }_{ijk}^{\left(1\right)}=2G{l}_1^2{\eta }_{ijk}^{\left(1\right)}$$

(9)

$$m_{ij}^s=2G{l}_2^2{\chi }_{ij}^s$$

(10)

where $l_0,l_1,l_2$ are the additional material length scale parameters. $\lambda$ and $G$ are the Lamé constants and can be defined in terms of Young’s modulus and Poisson’s ratio as

$$\lambda =\frac{Ev}{\left(1+v\right)\left(1-2v\right)}$$

(11)

$$G=\frac{E}{2\left(1+v\right)}$$

(12)

The displacement components of an initially straight Bernoulli–Euler beam can be defined as

$$u\left(x,z\right)=-z\frac{\mathrm{d}w\left(x\right)}{\mathrm{d}x}$$

(13)

$$v\left(x,z\right)=0$$

(14)

$$w\left(x,z\right)=w\left(x\right)$$

(15)

where $u, v$, and $w$ are the components of the displacement vector in $x-, y-,$ and $z-$ axes, respectively. Using Equations (13)–(15) in Equations (2)–(5) yields the following nonzero components of the classical and higher-order strains as

$${\epsilon }_{xx}=-z\frac{{\mathrm{d}}^{2}w}{\mathrm{d}{x}^2}$$

(16)

$${\gamma }_x=-z\frac{{\mathrm{d}}^{3}w}{\mathrm{d}{x}^3}$$

(17)

$${\gamma }_z=-\frac{{\mathrm{d}}^{2}w}{\mathrm{d}{x}^2}$$

(18)

$${\eta }_{xxx}^{\left(1\right)}=-\frac{2}{5}\left(z\frac{{\mathrm{d}}^{3}w}{\mathrm{d}{x}^3}\right)$$

(19)

$${\eta }_{xxz}^{\left(1\right)}={\eta }_{xzx}^{\left(1\right)}={\eta }_{zxx}^{\left(1\right)}=-\frac{4}{15}\left(\frac{d^{2}w}{d{x}^2}\right)$$

(20)

$${\eta }_{xyy}^{\left(1\right)}={\eta }_{xzz}^{\left(1\right)}={\eta }_{yxy}^{\left(1\right)}={\eta }_{yyx}^{\left(1\right)}={\eta }_{zxz}^{\left(1\right)}={\eta }_{zzx}^{\left(1\right)}=\frac{1}{5}\left(z\frac{{\mathrm{d}}^{3}w}{\mathrm{d}{x}^3}\right)$$

(21)

$${\eta }_{yyz}^{\left(1\right)}={\eta }_{yzy}^{\left(1\right)}={\eta }_{zyy}^{\left(1\right)}=\frac{1}{15}\left(\frac{d^{2}w}{d{x}^2}\right)$$

(22)

$${\eta }_{zzz}^{\left(1\right)}=\frac{1}{5}\left(\frac{{\mathrm{d}}^{2}w}{\mathrm{d}{x}^2}\right)$$

(23)

$${\chi }_{xy}^s={\chi }_{yx}^s=-\frac{1}{2}\left(\frac{d^{2}w}{d{x}^2}\right)$$

(24)

By substituting Equations (16)–(24) into Equations (7)–(10), the nonzero components of macro and micro stresses can be obtained as (by neglecting Poisson effect [29])

$${\sigma }_{xx}=-Ez\frac{{\mathrm{d}}^{2}w}{\mathrm{d}{x}^2}$$

(25)

$$p_x=-2G{l}_0^{2}z\frac{{\mathrm{d}}^{3}w}{\mathrm{d}{x}^3}$$

(26)

$$p_z=-2G{l}_0^2\frac{d^{2}w}{\mathrm{d}{x}^2}$$

(27)

$${\tau }_{xxx}^{\left(1\right)}=-\frac{4}{5}G{l}_1^{2}z\frac{{\mathrm{d}}^{3}w}{d{x}^3}$$

(28)

$${\tau }_{xxz}^{\left(1\right)}={\tau }_{xzx}^{\left(1\right)}={\tau }_{zxx}^{\left(1\right)}=-\frac{8}{15}G{l}_1^2\frac{d^{2}w}{\mathrm{d}{x}^2}$$

(29)

$${\tau }_{xyy}^{\left(1\right)}={\tau }_{xzz}^{\left(1\right)}={\tau }_{yxy}^{\left(1\right)}={\tau }_{yyx}^{\left(1\right)}={\tau }_{zxz}^{\left(1\right)}={\tau }_{zzx}^{\left(1\right)}=\frac{2}{5}G{l}_1^{2}z\frac{{\mathrm{d}}^{3}w}{d{x}^3}$$

(30)

$${\tau }_{yyz}^{\left(1\right)}={\tau }_{yzy}^{\left(1\right)}={\tau }_{zyy}^{\left(1\right)}=\frac{2}{15}G{l}_1^2\frac{d^{2}w}{\mathrm{d}{x}^2}$$

(31)

$${\tau }_{zzz}^{\left(1\right)}=\frac{2}{5}G{l}_1^2\frac{d^{2}w}{\mathrm{d}{x}^2}$$

(32)

$$m_{xy}^s=m_{yx}^s=-G{l}_2^2\frac{d^{2}w}{\mathrm{d}{x}^2}$$

(33)

Using above equations into Equation (1) yields an expression for the strain energy $U$ as

$$U=\frac{1}{2}{{\displaystyle \int }}_0^L\left[\left(EI+GA\left(2l_0^2+\frac{8}{15}{l}_1^2+l_2^2\right)\right){\left(\frac{d^{2}w}{d{x}^2}\right)}^2+\left(GI\left(2l_0^2+\frac{4}{5}{l}_1^2\right)\right){\left(\frac{d^{3}w}{d{x}^3}\right)}^2\right]dx$$

(34)

where $A$ and $I$ are the cross-section and second moment of area, respectively.

3. Buckling of an AFG-Tapered Microcolumn

In the present study, it is considered that the cross-section $A_{\left(x\right)}$ and second moment of area $I_{\left(x\right)}$ are linearly varied along the longitudinal direction. Moreover, the material properties of Young’s modulus $E_{\left(x\right)}$, the shear modulus $G_{\left(x\right)}$, and the length scale parameters $l_{i_{\left(x\right)}}$ change gradually in this direction. They can be described as

$$A\left(x\right)=A_L\left(1-\alpha \frac{x}{L}\right)$$

(35)

$$I_{\left(x\right)}=I_L\left(1-\alpha \frac{x}{L}\right)$$

(36)

$$E_{\left(x\right)}=E_L+\left(E_R-E_L\right)V_R$$

(37)

$$G_{\left(x\right)}=G_L+\left(G_R-G_L\right)V_R$$

(38)

$$l_{i_{\left(x\right)}}=l_{i_L}+\left(l_{i_R}-l_{i_L}\right)V_R(i=0,1,2)$$

(39)

where $\alpha$ represents the taper ratio, and the subscripts $L$ and $R$ denote left and right sides of the microbeam, respectively. $V_R$ defines the volume fraction of the right side’s material equals ${\left(\frac{x}{L}\right)}^n$ in which $n$ is the non-negative material property gradient index. It is notable that if $\alpha =0$, the microbeam will have a uniform section. On the other hand, if $n=0 or \infty$, the microbeam will be homogeneous. Variation of the volume fraction of right side’s material throughout the longitudinal direction for various gradient indices is depicted in Figure 1. It is seen from this figure that $V_R$ becomes bigger for lower values of the material property gradient index and approaches to zero by increasing the material property gradient index.

The strain energy $U$ in Equation (34) can be rewritten for AFG-tapered microbeams as following

$$\begin{array}{c}{U}^{*}=\frac{1}{2}{{\displaystyle \int }}_0^L[\left(E_{\left(x\right)}{I}_{\left(x\right)}+G_{\left(x\right)}{A}_{\left(x\right)}\left(2l_{0\left(x\right)}^2+\frac{8}{15}{l}_{1\left(x\right)}^2+l_{2\left(x\right)}^2\right)\right){\left(\frac{d^{2}w}{d{x}^2}\right)}^2\\ +\left(G_{\left(x\right)}{I}_{\left(x\right)}\left(2l_{0\left(x\right)}^2+\frac{4}{5}{l}_{1\left(x\right)}^2\right)\right){\left(\frac{d^{3}w}{d{x}^3}\right)}^2]dx\end{array}$$

(40)

4. Application of Rayleigh–Ritz Method

The total potential energy of the AFG-tapered microbeam can be expressed as

$$\Pi =U^{*}-W$$

(41)

where $U^{*}$ is the strain energy of the AFG-tapered microbeam and $W$ is the work done by axial compressive load as

$$W=\frac{P}{2}{{\displaystyle \int }}_0^L{\left(\frac{dw}{dx}\right)}^{2}dx$$

(42)

where $P$ is the axial compressive load. The total potential energy of the AFG-tapered microbeam is written by using Equations (40) and (42) in Equation (41) as

$$\begin{array}{c}\Pi =\frac{1}{2}{{\displaystyle \int }}_0^L[\left(E_{\left(x\right)}{I}_{\left(x\right)}+G_{\left(x\right)}{A}_{\left(x\right)}\left(2l_{0\left(x\right)}^2+\frac{8}{15}{l}_{1\left(x\right)}^2+l_{2\left(x\right)}^2\right)\right){\left(\frac{d^{2}w}{d{x}^2}\right)}^2\\ +\left(G_{\left(x\right)}{I}_{\left(x\right)}\left(2l_{0\left(x\right)}^2+\frac{4}{5}{l}_{1\left(x\right)}^2\right)\right){\left(\frac{d^{3}w}{d{x}^3}\right)}^2]dx-\frac{P}{2}{{\displaystyle \int }}_0^L{\left(\frac{dw}{dx}\right)}^{2}dx\end{array}$$

(43)

Utilizing Rayleigh–Ritz method with assumed trial function, $w\left(x\right)$ can be defined as

$$w\left(x\right)={\displaystyle \sum }_{i=1}^N{a}_i{\phi }_i\left(x\right)$$

(44)

where $a_i$ are the constants and ${\phi }_i\left(x\right)$ is the admissible function which is necessary to satisfy only the geometric (essential) boundary conditions. The geometric boundary conditions of clamped and pinned ends for classical case are

$$\mathrm{Clamped}\left(\mathrm{C}\right):w=0,\frac{dw}{dx}=0$$

(45)

$$\mathrm{Pinned}\left(\mathrm{P}\right):w=0$$

(46)

For nonclassical boundary conditions, two enable boundary conditions are possible at fixed ends as [5,15]

$$M_{nc}=\frac{d^{3}w}{d{x}^3}=0$$

(47)

$$\frac{d^{2}w}{d{x}^2}=0$$

(48)

In the present study, two admissible functions are chosen to be dependent on cantilever (clamped-free) and propped cantilever (clamped-pinned) microbeams as following

$$\mathrm{Cantilever}\left(\mathrm{BC}1\right):\hspace{1em}{\phi }_i^1\left(x\right)={\left(x\right)}^{i+1},i=1,2,\dots ,10$$

(49)

$$\mathrm{Cantilever}\left(\mathrm{BC}2\right):\hspace{1em}{\phi }_i^2\left(x\right)={\left(x\right)}^{i+2}, i=1,2,\dots ,10$$

(50)

$$\mathrm{Propped}\mathrm{cantilever}\left(\mathrm{BC}1\right):{\phi }_i^1\left(x\right)={\left(x\right)}^{i+1}\left(L-x\right),i=1,2,\dots ,10$$

(51)

$$\mathrm{Propped}\mathrm{cantilever}\left(\mathrm{BC}2\right):{\phi }_i^2\left(x\right)={\left(x\right)}^{i+2}\left(L-x\right), i=1,2,\dots ,10$$

(52)

The stationary points of the total potential energy are the solutions which satisfy equilibrium as

$$\frac{\partial \Pi }{\partial a_i}=0,i=1,2,\dots ,10$$

(53)

For a non-trivial solution, the determinant of the coefficient matrix must be equal to zero. Consequently, ten roots of this eigenvalue problem can be assessed and the smallest of them indicates the critical buckling load.

For illustration purposes, the solution of Equation (53) is presented for homogeneous uniform cantilever beams with $i=2$. In this case, $w\left(x\right)$ can be written according to Equations (44) and (49) as

$$w\left(x\right)=a_1{x}^2+a_2{x}^3$$

(54)

The total potential energy of a homogeneous uniform beam can be expressed as

$$\Pi =\frac{1}{2}{{\displaystyle \int }}_0^L\left[EI{\left(\frac{d^{2}w}{d{x}^2}\right)}^2-P{\left(\frac{dw}{dx}\right)}^2\right]dx$$

(55)

Substituting Equation (54) in Equation (55) yields

$$\begin{array}{c}\Pi =\frac{1}{2}{{\displaystyle \int }}_0^L\left[EI{\left(2a_1+6a_{2}x\right)}^2-P{\left(2a_{1}x+3a_2{x}^2\right)}^2\right]dx\\ =2EIL\left(a_1^2+3a_1{a}_{2}L+3a_2^2{L}^2\right)-\frac{P{L}^3}{30}\left(20a_1^2+45a_1{a}_{2}L+27a_2^2{L}^2\right)\end{array}$$

(56)

Use of Equation (56) in Equation (53) gives the following relations

$$\frac{\partial \Pi }{\partial a_1}=EIL\left(4a_1+6a_{2}L\right)-P{L}^3\left(\frac{4}{3}{a}_1+\frac{3}{2}{a}_{2}L\right)=0$$

(57)

$$\frac{\partial \Pi }{\partial a_2}=EIL\left(6a_{1}L+12a_2{L}^2\right)-P{L}^3\left(\frac{3}{2}{a}_{1}L+\frac{9}{5}{a}_2{L}^2\right)=0$$

(58)

Equations (57) and (58) can be written in matrix form as

$$\left[\begin{array}{cc}{K}_{11}& K_{12}\\ K_{21}& K_{22}\end{array}\right]\left\{\begin{array}{c}{a}_1\\ a_2\end{array}\right\}=0$$

(59)

where

$$K_{11}=4L\left(EI-\frac{P{L}^2}{3}\right),K_{12}=K_{21}=3L^2\left(2EI-\frac{P{L}^2}{2}\right),K_{22}=3L^3\left(4EI-\frac{3P{L}^2}{5}\right)$$

(60)

Two dimensionless roots are obtained by equaling to zero the determinant of the coefficient matrix in Equation (59) and multiplying $\frac{L^2}{EI}$ as

$${\lambda }_1=2.4860\mathrm{and}{\lambda }_2=32.1807$$

(61)

The smallest root (${\lambda }_1$) indicates the critical buckling load. It is notable that Equation (56) is written for $i=10$ as

$$\left[\begin{array}{ccc}{K}_{11}& \cdots & K\_110\\ ⋮& \ddots & ⋮\\ K_{101}& \cdots & K_{1010}\end{array}\right]\left\{\begin{array}{c}{a}_1\\ ⋮\\ a_{10}\end{array}\right\}=0$$

(62)

Similarly, ten roots can be assessed and the smallest of them indicates the critical buckling load.

5. Results and Discussion

First, in order to prove the validity and accuracy of the current results, comparisons and convergences for the dimensionless classical critical buckling loads of tapered homogeneous columns are presented in

Table 1. Comparison and convergence of dimensionless critical buckling loads for homogeneous tapered cantilever microbeams ($l_0=l_1=l_2=0$).

Taper Ratio, $\mathit{\alpha }$	Present					Exact
	N = 2	N = 4	N = 6	N = 8	N = 10	[49]
0	2.4860	2.4674	2.4674	2.4674	2.4674	2.467
0.1	2.4165	2.3928	2.3928	2.3928	2.3928	2.393
0.3	2.2702	2.2351	2.2351	2.2351	2.2351	2.235
0.5	2.1078	2.0621	2.0621	2.0621	2.0621	2.062
0.7	1.9140	1.8655	1.8653	1.8653	1.8653	1.865
0.9	1.6462	1.6229	1.6212	1.6211	1.6211	1.621

and

Table 2. Comparison and convergence of dimensionless critical buckling loads for homogeneous tapered propped cantilever microbeams ($l_0=l_1=l_2=0$).

Taper Ratio, $\mathit{\alpha }$	Present					Exact
	N = 2	N = 4	N = 6	N = 8	N = 10	[49]
0	20.9187	20.1943	20.1907	20.1907	20.1907	20.1907
0.1	19.9908	19.1736	19.1686	19.1686	19.1686	19.17
0.3	18.0621	17.0426	17.0354	17.0353	17.0353	17.03
0.5	15.9584	14.7452	14.7394	14.7394	14.7394	14.74
0.7	13.4459	12.1775	12.1773	12.1772	12.1772	12.18
0.9	9.7113	9.0608	9.0307	9.0295	9.0294	9.029

for various taper ratios and cantilever and propped cantilever microbeams, respectively. It is revealed from the table that the exact and present results agree very well for two boundary conditions and all taper ratios.

Because the comparative results above are based on classical theory, another comparison of size-dependent dimensionless critical buckling loads for homogeneous tapered cantilever microbeams made of epoxy is given in

Table 3. Comparison of size-dependent dimensionless critical buckling loads for tapered cantilever epoxy microbeams ($l_0=l_1=l_2=h_1$).

Taper Ratio, $\mathit{\alpha }$	BC1		BC2
	Present	[40]	Present	[40]
0	40.38752	40.3875	41.52736	41.5274
0.2	35.20555	35.2055	36.04467	36.0447
0.4	29.81004	29.8100	30.38676	30.3868
0.6	24.04835	24.0483	24.40082	24.4008
0.8	17.51253	17.5125	17.67871	17.6787

. For the purpose of comparison, the material and geometric properties are used as indicated in Ref. [40] for this table only. It is observed that there is an excellent agreement between the compared results.

Table 4. Dimensionless buckling loads for nonhomogeneous tapered microbeams for various gradient indices (α = 0.25, $l_L=h_1$).

n	Cantilever			Propped Cantilever
	CT	SGT-BC1	SGT-BC2	CT	SGT-BC1	SGT-BC2
0	2.2757	11.5570	11.8369	17.5814	89.5126	91.6124
0.5	3.2955	29.7361	30.0206	22.8377	175.9789	178.2837
1	3.7659	43.0086	43.5267	25.4402	232.7055	236.3154
2	4.1851	59.4975	60.5222	28.5113	314.3887	319.9357
10	4.5350	77.8720	79.8187	34.3698	564.3991	577.0945

and

Table 5. Dimensionless buckling loads for nonhomogeneous tapered microbeams for various gradient indices (α = 0.75, l_L = h₁).

n	Cantilever			Propped Cantilever
	CT	SGT-BC1	SGT-BC2	CT	SGT-BC1	SGT-BC2
0	1.8104	9.1937	9.3599	11.4694	58.3994	59.4905
0.5	2.5358	21.4730	21.6085	14.5021	107.0497	108.1194
1	2.8923	30.5769	30.8118	16.0267	138.5213	140.1774
2	3.2404	43.0156	43.4926	17.8791	183.6129	186.1580
10	3.5961	62.2555	62.3944	22.0821	348.4897	354.3562

show the variation of dimensionless critical buckling loads for AFG-tapered microbeams for various gradient indices. It is assumed in the analysis that $E_L=2E_R$, $l_L=2l_R$ and $v_L=v_R=0.3$. The tables demonstrate that an increase in the material property gradient index gives rise to an increment in the critical buckling loads. Furthermore, it can be said that the critical buckling loads obtained based on SGT are always bigger than those evaluated by CT. In addition, it is seen from these two tables that the critical buckling loads for $\alpha =0.25$ are greater than those for $\alpha =0.75$. Moreover, it is notable that the critical buckling loads obtained for BC2 are larger than those for BC1 and that this situation is more prominent for propped cantilever microbeams.

Variations of dimensionless buckling loads for cantilever and propped cantilever AFG-tapered microbeams with respect to h₁/l for different taper ratios are illustrated in Figure 2 and Figure 3, respectively. It can be seen from these figures that size-dependent buckling loads decrease by increasing h₁/l_L, approaching classical results. Moreover, it is observed that the critical buckling loads for $\alpha =0.3$ are greater than those for $\alpha =0.6$.

Influences of nonclassical boundary conditions on the size-dependent dimensionless buckling loads with respect to various taper ratios for cantilever and propped cantilever microbeams are depicted in Figure 4 and Figure 5, respectively. It can be found from the figures that a decrement occurs in the buckling loads when increasing the taper ratio. Furthermore, the difference between the critical buckling loads evaluated for BC1 and BC2 is more evident for smaller taper ratios, especially for a uniform cross-section when $\alpha =0$. This difference diminishes by increasing the taper ratio.

6. Conclusions

In this study, buckling behavior of axially functionally graded non-prismatic microbeams is examined via modified strain gradient theory. The Rayleigh–Ritz method is implemented to solve the problem with cantilever and propped cantilever microbeams for various taper ratios and gradient indices. A detailed parametric study is performed to investigate the influences of taper ratio, material gradient index, length scale parameter, and boundary conditions on the critical buckling loads. It can be concluded that size effect is more important for smaller h₁/l_L. Futhermore, non-classical boundary conditions are more effective for propped cantilever microbeams than cantilever ones.