**5. Numerical results**

In this section, buckling smart control of SLGS using elastically bonded PVDF nanoplate is discussed so that the effects of nonlocal parameter, mode number, Pasternak foundation, and SLGS length on the buckling of the SLGS are also considered. For this purpose, buckling load ratio is defined as follows:

$$\text{Buckling Load Ratio} = \frac{\text{Buckling load from nonlocal theory } (P\_{nl})}{\text{Buckling load from local theory } (P\_l)}$$

The orthotropic mechanical properties of SLGS with thickness *h* =0.34 *nm* are listed in Table 1 [25]. Also, mechanical and electrical characteristics of PVDF nanoplate with thickness *h* =2 *nm* are presented in Table 1 [26], which are taken from the macroscopic piezoelectric materials. In fact, the polymeric piezoelectric nanomaterials' characteristics are size-dependent and determined by experiments or molecular dynamic (MD) simulation [22]. To the best of the author's knowledge, the experimental values of piezoelectric nanomaterials' properties are not available in the literature. Therefore, in order to quantitatively analyze the present work, the characteristics of PVDF nanoplate are assumed in macroscale.

The developed nonlocal theory to date is incapable of determining the small scaling parameter *e*0*a*. However, Eringen [11] proposed *e*0*a* =0.39 *nm* by matching the dispersion curves using nonlocal theory for plane waves and Born–Karman model of lattice dynamics. For CNT, the


**Table 1.** Material properties of PVDF and SLGS.

As mentioned above, it is assumed that the SLGS plate is free from any transverse loadings. Uniform compressive edge loading along *x* and *y* axis are *Nxm* = − *P* and *N ym* = −*kP*, respectively.

> 000 1000000 0000 0000000 0000 0000000

*w*

*x y*

y

y

y

ú êú ê ú ú êú ê ú ú êú ê ú

*x y*

ê ú ë û

are defined in Appendix A.

y

f

*l P P*

*w* (32)

ê ê ê

and *L ij*

0 0 0000000 000 0000000 000 0000000 000 0000000

î þë û

*<sup>L</sup>* ) 2 <sup>+</sup> *<sup>k</sup>*( *<sup>n</sup><sup>π</sup> b* ) 2

Finally, buckling load of the system (*P*) can be calculated by solving the above equation.

In this section, buckling smart control of SLGS using elastically bonded PVDF nanoplate is discussed so that the effects of nonlocal parameter, mode number, Pasternak foundation, and SLGS length on the buckling of the SLGS are also considered. For this purpose, buckling load

Buckling load from nonlocal theory ( ) Buckling Load Ratio Buckling load from local theory ( ) <sup>=</sup> *nl*

The orthotropic mechanical properties of SLGS with thickness *h* =0.34 *nm* are listed in Table 1 [25]. Also, mechanical and electrical characteristics of PVDF nanoplate with thickness *h* =2 *nm* are presented in Table 1 [26], which are taken from the macroscopic piezoelectric materials. In fact, the polymeric piezoelectric nanomaterials' characteristics are size-dependent and determined by experiments or molecular dynamic (MD) simulation [22]. To the best of the author's knowledge, the experimental values of piezoelectric nanomaterials' properties are not available in the literature. Therefore, in order to quantitatively analyze the present work, the

The developed nonlocal theory to date is incapable of determining the small scaling parameter *e*0*a*. However, Eringen [11] proposed *e*0*a* =0.39 *nm* by matching the dispersion curves using nonlocal theory for plane waves and Born–Karman model of lattice dynamics. For CNT, the

characteristics of PVDF nanoplate are assumed in macroscale.

+

é ù é ù ê ú ê ú

<sup>ê</sup> ë û ë û

), *ηmn* =( *<sup>m</sup><sup>π</sup>*

x h

ì üé ù é ù ï ïê ú ê ú ê ú ê ú ê ú ê ú ê ú ê ú í ýê ú <sup>=</sup> ê ú ï ïú êú ê ú

*mn mn*

Substituting Eq. (31) into Eqs. (16)–(19) and (28)–(30) yields:

54 55 56 57 64 65 66 67 74 75 76 77

*LLLL LLLL LLLL*

11 12 13 14 21 22 23 31 32 33

48 Graphene - New Trends and Developments

*LLLL LLL LLL*

where *ξmn* =1 <sup>+</sup> *<sup>μ</sup>*(( *<sup>m</sup><sup>π</sup>*

**5. Numerical results**

ratio is defined as follows:

41 44 45 46 47

*<sup>L</sup>* ) 2 <sup>+</sup> ( *<sup>n</sup><sup>π</sup> b* ) 2

*L LLLL P*

*e*0*a* is found to be less than 2 *nm* [27]. Small-scale effects on the buckling behavior of SLGS and PVDF nanoplate were carried out analytically by assuming a range of values 0<*e*0*a* <2 *nm*, since its exact value is not known [27].

In the absence of similar publications in the literature covering the same scope of the problem, one cannot directly validate the results found here. However, the present work could be partially validated based on a simplified analysis suggested by Samaei et al. [17], Pradhan [28], Murmu and Pradhan [29], and Hashemi and Samaei [16] on buckling of the SLGS for which the coupled PVDF nanoplate in this paper was ignored. For this purpose, an SLGS with *E* =1 *TPa*, *ν* =0.3, *h* =0.34 *nm*, and 5< *L* <50 *nm* [17] as well as a nanoplate with *E* =1765 *GPa*, *ν* =0.25, *h* =0.34 *nm*, and 5< *L* <50 *nm* [16] are considered.

On the same basis and assuming Mindlin plate theory for buckling of SLGS embedded in a Pasternak foundation, the results obtained here are compared with those of Samaei et al. [17]. The results are shown in Fig. 2, in which buckling load ratio versus nonlocal parameter is plotted for *Kg* =0 and *Kg* =20. It is noted that in this figure the shear modulus is normalized the same as [17], *Kg* <sup>=</sup>*kgL* <sup>2</sup> / *<sup>D</sup>* in which *<sup>D</sup>* <sup>=</sup>*Eh* <sup>3</sup> / 12(1−*<sup>ν</sup>* <sup>2</sup> ). As can be seen, the two analyses agree well and show similar results. A comparison between the buckling analysis of SLGS using the theories of classical plate [29], higher-order shear deformation [28] and Mindlin plate [16] is presented in Table 2. In this table, critical buckling load for different values of nonlocal parameter and aspect ratio of length to thickness is shown. As can be seen, the present results based on Mindlin plate theory closely match with those reported by Hashemi and Samaei [16], Pradhan [28], and Murmu and Pradhan [29].

Table 1. Material properties of PVDF and SLGS.

0 *e*0*a* 2 nm, since its exact value is not known [27].

*h* 0.34 nm , and 5 *L* 50 nm [16] are considered.

[16], Pradhan [28], and Murmu and Pradhan [29].

**Fig.2** Comparison of buckling load ratio versus nonlocal parameter for 0 *Kg* and 20 *Kg* .

The developed nonlocal theory to date is incapable of determining the small scaling parameter *e a*<sup>0</sup> . However, Eringen [11] proposed *e*0*a* 0.39 nm by matching the dispersion curves using nonlocal theory for plane waves and Born– Karman model of lattice dynamics. For CNT, the *e a*<sup>0</sup> is found to be less than 2 nm [27]. Small‐scale effects on the buckling behavior of SLGS and PVDF nanoplate were carried out analytically by assuming a range of values

In the absence of similar publications in the literature covering the same scope of the problem, one cannot directly validate the results found here. However, the present work could be partially validated based on a simplified analysis suggested by Samaei et al. [17], Pradhan [28], Murmu and Pradhan [29], and Hashemi and Samaei [16] on buckling of the SLGS for which the coupled PVDF nanoplate in this paper was ignored. For this purpose, an SLGS with *E* 1TPa ,

On the same basis and assuming Mindlin plate theory for buckling of SLGS embedded in a Pasternak foundation, the results obtained here are compared with those of Samaei et al. [17]. The results are shown in Fig. 2, in which buckling load ratio versus nonlocal parameter is plotted for 0 *Kg* and 20 *Kg* . It is noted that in this figure the shear

analyses agree well and show similar results. A comparison between the buckling analysis of SLGS using the theories of classical plate [29], higher‐order shear deformation [28] and Mindlin plate [16] is presented in Table 2. In this table, critical buckling load for different values of nonlocal parameter and aspect ratio of length to thickness is shown. As can

0.25,

. As can be seen, the two

0.3, *h* 0.34 nm , and 5 *L* 50 nm [17] as well as <sup>a</sup> nanoplate with *E* 1765 GPa ,

modulus is normalized the same as [17], *Kg kgL* / *D* <sup>2</sup> in which /12(1 ) <sup>3</sup> <sup>2</sup> *D Eh*

**Figure 2.** Comparison of buckling load ratio versus nonlocal parameter for *Kg* =0 and *Kg* =20.


**Table 2.** Comparison between the buckling analysis of SLGS using the theories of classical plate, higher-order shear deformation and Mindlin plate.

The effect of the external electric voltage (*V*0) on the buckling load ratio with respect to nonlocal parameter (*μ*) is demonstrated in Fig. 3. It is shown that applying positive electric potential can increase the buckling load ratio of the SLGS and vice versa. This is because the imposed positive and negative voltages generate the axial compressive and tensile forces in the PVDF nanoplate, respectively. Meanwhile, the effect of external voltage becomes more prominent at higher *μ*. Hence, the imposed external voltage is an effective controlling parameter for buckling of the SLGS which is coupled by a smart PVDF nanoplate. It is also concluded that increasing

) is

the *μ* decreases the buckling load ratio. This is due to the fact that the increase of nonlocal parameter decreases the interaction force between graphene sheet atoms, and that leads to a softer structure. Figure 4 illustrates the effect of mode number of SLGS on the variation of buckling load ratio versus *μ*. As can be seen, buckling load ratio decreases with increasing mode numbers. Also, the small-scale effects on the buckling load ratio become more distin‐ guished at higher modes. Obviously, the difference between the buckling load ratios of the SLGS is larger at higher nonlocal parameters. Furthermore, the buckling load ratio for all mode numbers decreases by increasing the *μ*. tensile forces in the PVDF nanoplate, respectively. Meanwhile, the effect of external voltage becomes more prominent at higher . Hence, the imposed external voltage is an effective controlling parameter for buckling of the SLGS which is coupled by a smart PVDF nanoplate. It is also concluded that increasing the decreases the buckling load ratio. This is due to the fact that the increase of nonlocal parameter decreases the interaction force between graphene sheet atoms, and that leads to a softer structure. Figure 4 illustrates the effect of mode number of SLGS on the variation of buckling load ratio versus . As can be seen, buckling load ratio decreases with increasing mode numbers. Also, the small‐scale effects on the buckling load ratio become more distinguished at higher modes. Obviously, the difference between the buckling load ratios of the SLGS is larger at higher nonlocal parameters. Furthermore, the buckling load ratio for all mode numbers decreases by increasing the .

demonstrated in Fig. 3. It is shown that applying positive electric potential can increase the buckling load ratio of the SLGS and vice versa. This is because the imposed positive and negative voltages generate the axial compressive and

Table 1. Material properties of PVDF and SLGS.

50 Graphene - New Trends and Developments

0 *e*0*a* 2 nm, since its exact value is not known [27].

*h* 0.34 nm , and 5 *L* 50 nm [16] are considered.

[16], Pradhan [28], and Murmu and Pradhan [29].

0

**Dimensionless critical buckling load from higher order plate theory [28]**

*L* / *h μ*

deformation and Mindlin plate.

0.2

0.4

0.6

Buckling Load Ratio,

0.8

1

The developed nonlocal theory to date is incapable of determining the small scaling parameter *e a*<sup>0</sup> . However, Eringen [11] proposed *e*0*a* 0.39 nm by matching the dispersion curves using nonlocal theory for plane waves and Born– Karman model of lattice dynamics. For CNT, the *e a*<sup>0</sup> is found to be less than 2 nm [27]. Small‐scale effects on the buckling behavior of SLGS and PVDF nanoplate were carried out analytically by assuming a range of values

In the absence of similar publications in the literature covering the same scope of the problem, one cannot directly validate the results found here. However, the present work could be partially validated based on a simplified analysis suggested by Samaei et al. [17], Pradhan [28], Murmu and Pradhan [29], and Hashemi and Samaei [16] on buckling of the SLGS for which the coupled PVDF nanoplate in this paper was ignored. For this purpose, an SLGS with *E* 1TPa ,

On the same basis and assuming Mindlin plate theory for buckling of SLGS embedded in a Pasternak foundation, the results obtained here are compared with those of Samaei et al. [17]. The results are shown in Fig. 2, in which buckling load ratio versus nonlocal parameter is plotted for 0 *Kg* and 20 *Kg* . It is noted that in this figure the shear

analyses agree well and show similar results. A comparison between the buckling analysis of SLGS using the theories of classical plate [29], higher‐order shear deformation [28] and Mindlin plate [16] is presented in Table 2. In this table, critical buckling load for different values of nonlocal parameter and aspect ratio of length to thickness is shown. As can be seen, the present results based on Mindlin plate theory closely match with those reported by Hashemi and Samaei

**Fig.2** Comparison of buckling load ratio versus nonlocal parameter for 0 *Kg* and 20 *Kg* .

*eo a*

> **Dimensionless critical buckling load from Mindlin plate theory (present work)**

**Dimensionless critical buckling load from Mindlin plate theory [16]**

0 0.5 1 1.5 2 2.5

Nonlocal Parameter, (nm)

**Dimensionless critical buckling load from classical plate theory [29]**

0.5 9.4156 9.4031 9.4029 9.4028 1.0 8.9947 8.9807 8.9803 8.9801 1.5 8.6073 8.5947 8.5939 8.5939 2.0 8.2537 8.2405 8.2393 8.2393

0.5 9.3570 9.3455 9.3455 9.3451 1.0 8.9652 8.9528 8.9527 8.9522 1.5 8.5546 8.5421 8.5420 8.5419 2.0 8.2114 8.1900 8.1898 8.1898

100 0.0 9.8791 9.8671 9.8671 9.8671

20 0.0 9.8177 9.8067 9.8067 9.8067

**Table 2.** Comparison between the buckling analysis of SLGS using the theories of classical plate, higher-order shear

The effect of the external electric voltage (*V*0) on the buckling load ratio with respect to nonlocal parameter (*μ*) is demonstrated in Fig. 3. It is shown that applying positive electric potential can increase the buckling load ratio of the SLGS and vice versa. This is because the imposed positive and negative voltages generate the axial compressive and tensile forces in the PVDF nanoplate, respectively. Meanwhile, the effect of external voltage becomes more prominent at higher *μ*. Hence, the imposed external voltage is an effective controlling parameter for buckling of the SLGS which is coupled by a smart PVDF nanoplate. It is also concluded that increasing

 =0, Samaei et al. [8] =0, Present work =20, Samaei et al. [8] =20, Present work

**Figure 2.** Comparison of buckling load ratio versus nonlocal parameter for *Kg* =0 and *Kg* =20.

*Kg Kg Kg Kg* 0.25,

. As can be seen, the two

0.3, *h* 0.34 nm , and 5 *L* 50 nm [17] as well as <sup>a</sup> nanoplate with *E* 1765 GPa ,

modulus is normalized the same as [17], *Kg kgL* / *D* <sup>2</sup> in which /12(1 ) <sup>3</sup> <sup>2</sup> *D Eh*

**Fig.3** The effect of the external voltage on the buckling load ratio versus nonlocal parameter. **Figure 3.** The effect of the external voltage on the buckling load ratio versus nonlocal parameter. Table 3. The effect of the external voltage on the buckling load ratio versus nonlocal parameter.

higher 

ratio versus

The effects of the SLGS length and the imposed external voltage to the PVDF nanoplate on the

buckling load ratio are shown in Fig. 5. It is noted that the length of the SLGS is considered

between10 *L* 40 nm , since the maximum length of the graphene sheet taken is 45.2896 nm in

the literature by Sakhaee Pour [30], Pradhan and Murmu [13], and Samaei et al. [17]. As length of

the graphene sheet increases, the buckling load ratio increases. It is also observed that for a given

**Fig.4** The effect of mode number on the buckling load ratio versus nonlocal parameter. **Figure 4.** The effect of mode number on the buckling load ratio versus nonlocal parameter.

The effects of the SLGS length and the imposed external voltage to the PVDF nanoplate on the buckling load ratio are shown in Fig. 5. It is noted that the length of the SLGS is considered between 10< *L* <40 *nm*, since the maximum length of the graphene sheet taken is 45.2896 *nm* in the literature by Sakhaee Pour [30], Pradhan and Murmu [13], and Samaei et al. [17]. As length of the graphene sheet increases, the buckling load ratio increases. It is also observed that for a given length, the SLGS, applying negative external voltage to PVDF nanoplate, will buckle first as compared to the SLGS with positive one.

Figure 6 depicts the effects of axial half wave number (*m*) and external voltage on the buckling load ratio of the SLGS. It is obvious that the buckling load ratio decreases sharply with increasing *m*. As can be seen, with the increase of external voltage, buckling load ratio is increased. Moreover, the effect of *V*0 is not considerable for *m*<2. It means that the external voltage effect decreases with decrease of the axial half wave number.

Buckling smart control of SLGS using PVDF nanoplate versus shear modulus parameter (*Kg*) is plotted in Fig. 7. The obtained results show that at a given *Kg*, when the imposed external voltage changes from *-*1 V to 1 V, the buckling load ratio increases. It is also worth mentioning that the influences of *V*0 at higher *Kg* values are more apparent than at lower *Kg* 's. As the shear modulus parameter of the coupled system increases, generally, the buckling load ratio reduces and approaches a constant value. This is because increasing shear modulus parameter increases the structure stiffness.

Fig. 5 The effect of the external voltage on the buckling load ratio versus graphene length.

 = -1 volt = -0.5 volt = 0 volt = 0.5 volt = 1 volt

*V0 V0 V0 V0 V0*

Fig. 6 The effect of the external voltage on the buckling load ratio versus axial half wave

0 1 2 3 4 5 6 7 8 9 10

*m*

Axial half wave number,

number.

**Figure 5.** The effect of the external voltage on the buckling load ratio versus graphene length.

0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65

Buckling Load Ratio,

 = -1 volt = -0.5 volt = 0 volt = 0.5 volt = 1 volt

*V0 V0 V0 V0 V0*

10 15 20 25 30 35 40

Length, (nm)

*L*

0.25

0.3

0.35

0.4

0.45

Buckling Load Ratio,

0.5

0.55

0.6

Fig. 6 The effect of the external voltage on the buckling load ratio versus axial half wave **Figure 6.** The effect of the external voltage on the buckling load ratio versus axial half wave number.

number.

Fig. 7 The effect of the external voltage on the buckling load ratio versus shear modulus **Figure 7.** The effect of the external voltage on the buckling load ratio versus shear modulus parameter.
