Flexoelectric Effect on Vibration of Piezoelectric Microbeams Based on a Modified Couple Stress Theory.
Piezoelectricity is a well-studied electromechanical coupling effect in which the mechanical strain leads to electric polarization in piezoelectric crystals, or vice versa. Due to the excellent electromechanical characteristics at microscale , piezoelectric based microstructures have found a wide range of applications in microtechnology, like microtransducers, microgenerators, microresonantors, and so forth [2, 3]. Particularly, these above microstructures are quasi one-dimensional structures which can be efficiently characterized by simple Euler-Bernoulli beam theory. Flexoelectricity [4-7] is, however, the coupling between the mechanical strain gradient and the electric polarization, and it is a universal electromechanical mechanism in all insulators including piezoelectric materials [8-11]. Shen and Hu [12, 13] introduced an electric Gibbs function to consider both of the piezoelectricity and flexoelectricity and derived the governing equations for dielectric materials. Flexoelectricity has also been found to be useful when it is used for sensing purpose [14, 15]. Moreover, the flexoelectric structures are theoretically predicted to be more sensitive when scaled down to microdomains [16,17], yet their work did not take the size effect at microscale into account.
At microscale, size effect may be significant and even dominate the behavior of structures. Several strain gradient theories were dedicated to investigating the size effect . Among these work, Yang et al.  developed a modified couple stress theory with only one high-order material constant to account for the size effect on the microscale structures which makes the size effect more convenient to express. Afterwards, many efforts were devoted to studying the size effect on the mechanical behaviors of microbeams using the modified couple stress theory [22-25]. Among these, Park and Gao  developed a new model for the bending of a Bernoulli-Euler beam. Ma et al.  employed the Timoshenko beam model to study the size-dependence of static bending and vibration behaviors. Reddy  modified Euler-Bernoulli and Timoshenko models for functionally graded beams and studied the size-dependence of the deflection, vibration, and buckling behaviors for beam models. In these studies, the modified couple stress theory was used for elastic beams whereas the electromechanical coupling effects, for example, piezoelectric effect and flexoelectric effect, were not considered. Recently, the modified couple stress theory was utilized by Ansari et al.  for discussing the vibration characteristics of microbeams with piezoelectricity; however, flexoelectricity was not taken into account. Li et al.  reformulated a flexoelectric theory to study the bending behavior of PMBs with flexoelectricity. However, the vibration behaviors of PMBs with flexoelectricity have not been reported to date, which may be as important as the bending behavior in MEMS.
The goal of this paper is to study the natural frequency of PMBs with considering the size effect, piezoelectricity, and flexoelectricity. To achieve this goal, together with the Euler-Bernoulli beam theory, a new Gibbs function needs to be defined by employing a modified couple stress theory to derive the governing equations for the PMBs. The proposed study can be fundamental for the vibration behaviors of PMBs which may have potential application in the microelectromechanical systems (MEMS).
This paper is organized as follows. After introduction, Section 2 presents an electric Gibbs function based on the flexoelectricity theory and the modified couple stress theory. The dynamics equation and boundary conditions of piezoelectric beams were then derived by using the new Gibbs function and the Euler-Bernoulli beam model. Afterwards, the first-order natural frequencies under different boundary conditions of the current model and the traditional Euler-Bernoulli model were displayed. In particular, the natural frequency of simply supported PMBs subjected to a voltage across the thickness is solved as an example. Then the relationship between the frequency and the effective bending stiffness was discussed and employed to study the effects of couple stress, piezoelectricity, and flexoelectricity in Section 3. The new findings were summarized in Section 4.
2. Dynamics Equation of Piezoelectric Microbeams with Flexoelectricity
An Euler-Bernoulli beam in Figure 1 is considered:
[u.sub.x] = -[partial derivative]w(x, t)/[partial derivative]x z, [u.sub.y] = 0, [u.sub.z] = w(x, t). (1)
To account for the piezoelectric effect, flexoelectric effect, and nonlocal mechanical and electrical effects, an electric Gibbs energy density Gb including the energy of elasticity, piezoelectricity, flexoelectricity, and couple stress was proposed by
[mathematical expression not reproducible], (2)
where [k.sub.ij] is dielectric constant tensor, [E.sub.i] and [E.sub.j] are electric field vectors, [b.sub.ijkl] is the nonlocal electrical coupling coefficient tensor, [C.sub.ijkl] is elastic stiffness tensor, [[epsilon].sub.kl] is strain tensor, [e.sub.ijk] is piezoelectric coefficient tensor, [[mu].sub.ijkl] is a fourth- order tensor of the flexoelectric coefficient, and [[epsilon].sub.kl,j] and [E.sub.k,l] are the gradients of strain and electric field. [m.sub.ij] is deviator part of couple stress tensor and [[chi].sub.ij] is the symmetric part of curvature tensor, which, using the modified couple stress theory [20, 22], are defined as
[m.sub.ij] = 2[l.sup.2] [mu][[chi].sub.ij], [[chi].sub.ij] = 1/2 ([[theta].sub.i,j] + [[theta].sub.j,i]), [[theta].sub.i] = 1/2 [e.sub.ijk][u.sub.k,l], (3)
where [mu] is Lame's constants, I is the material length scale parameter, [[theta].sub.i] is the rotation vector, and [u.sub.i] is the gradient of displacement vector, respectively.
Thus the constitutive equations for piezoelectric Euler-Bernoulli beams can be derived as
[mathematical expression not reproducible], (4)
where [[tau].sub.xxz], [D.sub.z], and [Q.sub.zz] denote the higher-order stress tensor, electric displacement vector, and electric quadrupole, respectively. In the absence of free charges, Gauss's law requires
-[[partial derivative].sup.2][Q.sub.zz]/[partial derivative][z.sup.2] + [partial derivative][D.sub.z]/[partial derivative]z = 0. (5)
Substituting (4) into (5) leads to
[k.sub.33][E.sub.z,z] + [e.sub.31][[epsilon].sub.xx,z] + [[mu].sub.31][[epsilon].sub.xx,zz] - [b.sub.33][E.sub.z,zzz] = 0 (6)
The electrical boundary conditions are given by
[PHI] (-h/2) = 0, [PHI](h/2) = V, (7)
where [PHI] is the electric potential. Combining [E.sub.z] = -[partial derivative][PHI]/[partial derivative]z with (4)-(7) yields
[mathematical expression not reproducible] (8)
with [eta] = [square root of ([k.sub.33]/[b.sub.33])]. Neglecting the high-order nonlocal electric effect in (8) by assuming [b.sub.33] = 0 yields the following equations [27, 28]:
[mathematical expression not reproducible]. (9)
The elastic energy U , kinetic energy T, and external work are calculated by the following equations:
[mathematical expression not reproducible]. (11)
Then the governing equation for their transverse vibration can be derived by the vibrational principle
[mathematical expression not reproducible], (12)
[mathematical expression not reproducible]. (13)
For facilitating, make the substitution of [(EI).sub.e] = [C.sub.11]I, [(EI).sub.p] = [e.sup.2.sub.31] I/[k.sub.33], (EI) c = A[[mu].sup.2], and [(EI).sub.f] = [[mu].sup.2.sub.31]A/2[k.sub.33], and then we have
[(EI).sub.eff] = [(EI).sub.e] + [(EI).sub.p] + [(EI).sub.c] + [(EI).sub.f], [P.sub.eff] = [e.sub.31] Vb. (14)
If the flexoelectric effect is neglected, the effective bending rigidity can be reduced to
[(EI).sup.p.sub.eff] = ([lambda] + 2[mu])I + [e.sup.2.sub.31]I/[k.sub.33] + A[mu][l.sup.2], [P.sup.p.sub.eff] = [e.sub.31]Vb. (15)
If the working temperature of PMBs is higher than Curie's point, the piezoelectric effect will disappear and then the effective bending rigidity can be rewritten as
[(EI).sup.f.sub.eff] = ([lambda] + 2[mu])I + [[mu].sup.2.sub.31]A/2[k.sub.33] + A[mu][l.sup.2], [P.sup.f.sub.eff] = 0. (16)
The geometric parameters are defined as
[[integral].sub.A] (1, z, [z.sup.2]) dA = (A, 0, I) (17)
with A and I being the cross section area and inertia moment of beams. The natural frequencies for PMBs with different boundary conditions, such as simply supported (SS), cantilever (C-F), and clamped-clamped (C-C) types, can be calculated based on the following characteristic equations:
[mathematical expression not reproducible], (18)
[mathematical expression not reproducible] (19)
with [omega] representing the angular resonant frequency.
3. Results and Discussion
In this section the formulae derived in Section 2 will be used to study the fundamental vibration of the PMBs. The emphasis is on the effects of size, piezoelectricity, and flexoelectricity. To this end, BaTi[O.sub.3] was chosen as the material of the PMBs where the material constants are listed in Table 1. The cross section of BaTi[O.sub.3] beam is assumed to be square; that is, b = h.
Employing (13), (18), and (19) and the material constants in Table 1, we can easily get the first-order natural frequency of the PMBs Wf. Moreover, the first natural frequency of traditional Euler-Bernoulli beam model [[omega].sub.0] = [(i[pi]/L).sup.2] [square root of (EI/[rho]A)] (i = 0.5 for C-F type, i = 1 for S-S type, and i = 1.5 for C-C type) is employed for comparison. Here [[omega].sub.f] denotes the natural frequency of the present model; that is, both effects of piezoelectricity and flexoelectricity are taken into consideration. Correspondingly, [[omega].sub.p] represents the natural frequency when only piezoelectric effect is considered. The difference between the present model and the traditional model can be seen in Figure 2. Hereby PMBs have h = 10 [micro]m and I = 5 [micro]m.
From Figure 2 we can see that the differences in the frequencies between the present model and the traditional model increase significantly when the beam length decreases to less than 400 [micro]m and then turn to be nearly the same when the beam length increases to more than 800 [micro]m. In general, the trends in the frequencies of the present model resemble that of the traditional model because the frequency equation of [[omega].sub.f][|.sub.V=0] = [(in/L).sup.2] [square root of ([(EI).sub.eff]/[rho]A)] and [[omega].sub.0] are very alike with small difference between the bending stiffness [(EI).sub.eff] and (EI).
3.1. Effect of Piezoelectricity and Flexoelectricity. Throughout the following sections of this work we will calculate the first normalized natural frequencies [[bar.[omega]].sub.p] = [[omega].sub.p]/[[omega].sub.0] and [[bar.[omega]].sub.f] = [[omega].sub.f]/[[omega].sub.0], and, for facilitating the discussion, the simply supported (S-S) boundary condition was chosen as a representative. The normalized frequencies are shown in Figure 3 for PMBs where length-to-thickness ratio L/h is fixed at 20 and the length L increases from 20 [micro]m to 1000 [micro]m. Accordingly the thickness h grows from 1 [micro]m to 50 [micro]m and the material length scale parameter I is defined by I = h/2 as others do in [25, 26] which increases from 0.5 [micro]m to 25 [micro]m. In addition, the voltage V = -2 V, 0, and 2 V are applied to reveal the effects of piezoelectricity and flexoelectricity.
In Figure 3(a), where only piezoelectric effect is considered (i.e., without flexoelectricity), [[bar.[omega]].sub.p] = 1.461 is obtained at V =0 showing that the size effect increases by 46.1% for all PMBs of different length L and length scale parameter I. Equations (13) and (14) show the contribution of the size effect [(EI).sub.c] on the effective bending stiffness [(EI).sub.eff] of the PMBs. At L = 20 [micro]m, [[bar.[omega]].sub.p] increases to 1.465 when V = -2 V is applied whereas it deceases to 1.457 at V = 2 V. This result shows that the external voltage leads to the increase (or decrease) of the frequency of PMBs in the presence of a negative (or positive) voltage. A 4 V external voltage can tune the frequency of PMBs by 0.55%, and the tuning of frequency can be greater under greater external voltage. Such a piezoelectric effect can be explained by the electromechanical coupling term [(EI).sub.p] in [(EI).sub.eff] and the extra-axial force [P.sub.eff] in (13) and (14). Specifically, such a piezoelectric effect decreases monotonically with rising length L and shows significant length-dependence. As a result [[bar.[omega]].sub.p] with V = -2 V or 2 V approaches [[bar.[omega]].sub.f] = 1.461 when the length L is sufficiently large; for example, L > 900 [micro]m. In this case the piezoelectric effect turns out to be negligible.
Different from Figure 3(a), both effects of piezoelectricity and flexoelectricity are considered in Figure 3(b) for the normalized frequency [[bar.[omega]].sub.f]. Thus the difference between the results in Figures 3(a) and 3(b) [i.e., ([[bar.[omega]].sub.f] - [[bar.[omega]].sub.p])] reveals the effect of flexoelectricity on the normalized fundamental frequency. First it is noted in Figure 3(b) that, at V = 0, [[bar.[omega]].sub.f] becomes length-dependent, which is in contrast to the length-independent [[bar.[omega]].sub.p] in Figure 3(a). At L = 20 [micro]m flexoelectric effect slightly increases the normalized frequency by 0.73% relative to [[bar.[omega]].sub.p]. The relative increment then decreases with rising L and becomes negligible for long and thick PMBs. Qualitatively similar trend is observed for the flexoelectric effect at V = -2 V, which raises the normalized frequency by
0.26% at L = 20 [micro]m. The increase of the frequency decreases with rising L and almost vanishes at L = 1000 [micro]m. On the other hand, when positive voltage 2 V is applied, it is noted in Figure 3(b) that [[bar.[omega]].sub.f] first decreases with rising L and reaches the lowest value at L = 302 [micro]m. It then increases with increase of L and reaches close to [[bar.[omega]].sub.f] = 1.461 when the length is greater than 900 [micro]m. This behavior of [[bar.[omega]].sub.f] at 2 V is qualitatively different from the corresponding L-dependence of [[bar.[omega]].sub.p] in Figure 3(a) and can be attributed to the competition between the reversed effects of piezoelectricity (with positive voltage) and the flexoelectricity on the normalized frequency. Here it is noted that considering flexoelectricity always leads to the increase of the frequency independent of the voltage applied. This behavior can be attributed to the fact that flexoelectricity influences the frequency via the electromechanical coupling term [(EI).sub.f] in [(EI).sub.eff] found in (14), which always greatens the effective bending stiffness and also makes the natural frequency size-dependent. Though it can be seen from Figure 3(a) the positive voltage decreases the natural frequency with decreasing L, when L is shorter than 302 [micro]m flexoelectric effect dominates the electromechanical behaviors resulting in a significant rise in the natural frequency.
3.2. Size Effect on the Electromechanical Coupling in Vibration. After investigating the effects of piezoelectricity and flexoelectricity, we would like to further reveal the influence of the size effect of piezoelectricity and flexoelectricity. Here we first calculated the l-dependence of [[bar.[omega]].sub.p] (without flexoelectricity) in Figure 4 for PMBs subjected to a voltage of -2V, 0, and 2 V, respectively. In the calculations two groups of PMBs are considered, which have an aspect ratio L/h = 20 and 40, respectively. In addition the scale material constant I is assumed to be a constant and increases from 0 (classical theory) to 2 [micro]m. Figure 4 shows that in general [[bar.[omega]].sub.p] increases with growing I for all the PMBs under a given voltage. It is worth mentioning that at V = 0 the results remain unchanged no matter what aspect ratio L/h is considered. Specifically, the curves associated with different nonzero voltages are almost parallel with one another. These results indicate that the change in normalized frequency, that is, ([[bar.[omega]].sub.p]|[V.sub.[not equal to]0] - [[bar.[omega]].sub.p][|.sub.V=0]) due to piezoelectric effect, is not very sensitive to the variation of the scale parameter I (the subscript shows the voltage applied). On the other hand, as [[bar.[omega]].sub.p][|.sub.V=0] rises with rising I, the relative change ([[bar.[omega]].sub.p][|.sub.V[not equal to]0] - [[bar.[omega]].sub.p][|.sub.V=0])/ttp\v=0 due to piezoelectric effect decreases significantly with increasing (the subscript shows the voltage considered). For example, at I = 0 relative change due to V = -2 V (V = -2V) is 0.56% (-0.56%) for PMBs with L = 40 h. However, at I = 2 [micro]m such a relative change decreases to 0.47% atV = -2V and -0.48% at V = 2 V. In addition it is also noticed that for a constant the piezoelectric effect is less pronounced for shorter PMBs with greater L/h = 20.
Next we shall further exam the influence of the size effect on the effect of flexoelectricity by calculating the dependence of [[bar.[omega]].sub.f-p] = [[bar.[omega]].sub.f] - [[bar.[omega]].sub.p] in Figure 5. It is easy to understand that [[bar.[omega]].sub.f-p] reflects the pure effect of flexoelectricity. Figure 5 shows that, at V = 0, [[bar.[omega]].sub.f-p] decreases from around 0.98 x [10.sup.-3] to 0.90 x [10.sup.-3] when I increases from 0 to 2.0 [micro]m. These results indicate that the flexoelectricity increases the frequency of PMBs by less than 0.1%, and such an effect is found to be even smaller for greater or larger size effect. In addition, at V = -2V (or 2 V) the effect of flexoelectricity on PMB increases (decreases) slightly but the influence of voltage becomes smaller for greater l. Indeed it is seen from Figure 5 that, in terms of L/h = 40 when l = 0, [[bar.[omega]].sub.f- p][|.sub.V=-2V] is 0.32% higher than [[omega].sub.f-p][|.sub.V=0] but it declines to 0.27% when l = 2 [micro]m. In other words, by increasing I the effective bending stiffness by flexoelectricity increases; nevertheless the frequency difference of [[bar.[omega]].sub.f-p] between existing external voltage and no external voltage diminishes.
A novel electric Gibbs function was proposed to describe the size effect on the electromechanical coupling behaviors of PMBs by incorporating a modified couple stress theory into Euler-Bernoulli beam model. Based on the obtained model, the effects of piezoelectricity and flexoelectricity and the size effect were examined for the vibration behavior of PMBs. It is found that, at the microscale, the size effect can raise the effective stiffness of the structure and thus increases the natural frequency. With the assumed value of the material length scale parameter, the couple stress dominates the mechanical behavior of the PMBs and increases with rising I. Piezoelectric effect is achieved, which can increase or decrease the frequency depending on the axial force induced by external voltage. In particular, the contribution of piezoelectric effect to the natural frequency was found to be decreasing with rising or geometric size of PMBs. In the meantime, the flexoelectricity can also raise the structural stiffness via an electromechanical coupling term, which is relatively small but detectable and always increases the frequency of PMBs. With the development of material science, the higher performance of flexoelectric materials might play an more important role in MEMS.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The research reported herein was supported under NSFC (no. 11272138) and Specialized Research Fund for the Doctoral Program of Higher Education of China (no. 20123227130002).
 E. S. Hung and S. D. Senturia, "Extending the travel range of analog-tuned electrostatic actuators," Journal of Microelectromechanical Systems, vol. 8, no. 4, pp. 497-505, 1999.
 X. Li, B. Bhushan, K. Takashima, C.-W. Baek, and Y.-K. Kim, "Mechanical characterization of micro/nanoscale structures for MEMS/NEMS applications using nanoindentation techniques," Ultramicroscopy, vol. 97, no. 1-4, pp. 481-494, 2003.
 Y. Moser and M. A. M. Gijs, "Miniaturized flexible temperature sensor," Journal of Microelectromechanical Systems, vol. 16, no. 6, pp. 1349-1354, 2007.
 A. K. Tagantsev, "Piezoelectricity and flexoelectricity in crystalline dielectrics," Physical Review B, vol. 34, no. 8, pp. 5883-5889, 1986.
 L. E. Cross, "Flexoelectric effects: charge separation in insulating solids subjected to elastic strain gradients," Journal of Materials Science, vol. 41, no. 1, pp. 53-63, 2006.
 P. Zubko, G. Catalan, A. Buckley, P. R. L. Welche, and J. F. Scott, "Strain-gradient-induced polarization in SrTiO3 single crystals," Physical Review Letters, vol. 99, no. 16, Article ID 167601, 4 pages, 2007.
 N. D. Sharma, R. Maranganti, and P. Sharma, "On the possibility of piezoelectric nanocomposites without using piezoelectric materials," Journal of the Mechanics and Physics of Solids, vol. 55, no. 11, pp. 2328-2350, 2007.
 W. Ma and L. E. Cross, "Observation of the flexoelectric effect in relaxor Pb(Mg1/3Nb2/3)O3 ceramics," Applied Physics Letters, vol. 78, no. 19, pp. 2920-2921, 2001.
 W. Ma and L. E. Cross, "Flexoelectric polarization of barium strontium titanate in the paraelectric state," Applied Physics Letters, vol. 81, no. 18, pp. 3440-3442, 2002.
 W. Ma and L. E. Cross, "Flexoelectric effect in ceramic lead zirconate titanate," Applied Physics Letters, vol. 86, no. 7, Article ID 072905, 2005.
 W. Ma and L. E. Cross, "Flexoelectricity of barium titanate," Applied Physics Letters, vol. 88, no. 23, Article ID 232902, 2006.
 S. Shen and S. Hu, "A theory of flexoelectricity with surface effect for elastic dielectrics," Journal of the Mechanics and Physics of Solids, vol. 58, no. 5, pp. 665-677, 2010.
 S. Hu and S. Shen, "Variational principles and governing equations in nano-dielectrics with the flexoelectric effect," Science China Physics, Mechanics and Astronomy, vol. 53, no. 8, pp. 1497-1504, 2010.
 X. Yan, W. Huang, S. Ryung Kwon, S. Yang, X. Jiang, and F.-G. Yuan, "A sensor for the direct measurement of curvature based on flexoelectricity," Smart Materials and Structures, vol. 22, no. 8, Article ID 085016, 2013.
 W. Huang, X. Yan, S. R. Kwon, S. Zhang, F.-G. Yuan, and X. Jiang, "Flexoelectric strain gradient detection using Ba0.64Sr 0.36TiO3 for sensing," Applied Physics Letters, vol. 101, no. 25, Article ID 252903, 2012.
 S. R. Kwon, W. B. Huang, S. J. Zhang, F. G. Yuan, and X. N. Jiang, "Flexoelectric sensing using a multilayered barium strontium titanate structure," Smart Materials and Structures, vol. 22, no. 11, Article ID 115017, 2013.
 S. R. Kwon, W. Huang, L. Shu, F.-G. Yuan, J.-P. Maria, and X. Jiang, "Flexoelectricity in barium strontium titanate thin film," Applied Physics Letters, vol. 105, no. 14, Article ID 142904, 2014.
 R. D. Mindlin and N. N. Eshel, "On first strain-gradient theories in linear elasticity," International Journal of Solids and Structures, vol. 4, no. 1, pp. 109-124, 1968.
 R. A. Toupin, "Theories of elasticity with couple-stress," Archive for Rational Mechanics and Analysis, vol. 17, pp. 85-112, 1964.
 F. Yang, A. C. M. Chong, D. C. C. Lam, and P. Tong, "Couple stress based strain gradient theory for elasticity," International Journal of Solids and Structures, vol. 39, no. 10, pp. 2731-2743, 2002.
 D. C. C. Lam, F. Yang, A. C. M. Chong, J. Wang, and P. Tong, "Experiments and theory in strain gradient elasticity," Journal of the Mechanics and Physics of Solids, vol. 51, no. 8, pp. 1477-1508, 2003.
 S. K. Park and X.-L. Gao, "Bernoulli-Euler beam model based on a modified couple stress theory," Journal of Micromechanics and Microengineering, vol. 16, no. 11, pp. 2355-2359, 2006.
 H. M. Ma, X.-L. Gao, and J. N. Reddy, "A microstructure-dependent Timoshenko beam model based on a modified couple stress theory," Journal of the Mechanics and Physics of Solids, vol. 56, no. 12, pp. 3379-3391, 2008.
 J. N. Reddy, "Microstructure-dependent couple stress theories of functionally graded beams," Journal of the Mechanics and Physics of Solids, vol. 59, no. 11, pp. 2382-2399, 2011.
 R. Ansari, M. A. Ashrafi, and S. Hosseinzadeh, "Vibration characteristics of piezoelectric microbeams based on the modified couple stress theory," Shock and Vibration, vol. 2014, Article ID 598292, 12 pages, 2014.
 A. Li, S. Zhou, L. Qi, and X. Chen, "A reformulated flexoelectric theory for isotropic dielectrics," Journal of Physics D: Applied Physics, vol. 48, no. 46, Article ID 465502, 2015.
 L. Shu, F. Li, W. Huang, X. Wei, X. Yao, and X. Jiang, "Relationship between direct and converse flexoelectric coefficients," Journal of Applied Physics, vol. 116, no. 14, Article ID 144105, 2014.
 X. Liang, S. Hu, and S. Shen, "Size-dependent buckling and vibration behaviors of piezoelectric nanostructures due to flexoelectricity," Smart Materials and Structures, vol. 24, no. 10, Article ID 105012, 2015.
Xingjia Li and Ying Luo
Faculty of Civil Engineering and Mechanics, Jiangsu University, Zhenjiang
Correspondence should be addressed to Ying Luo; email@example.com
Received 12 October 2016; Accepted 12 February 2017; Published 16 March 2017
Academic Editor: Paulo B. Goncalves
Caption: FIGURE 1: The schematic of the PMB whose cross section is a rectangle with width b and height h. The origin of x, y, and z coordinate system is located at the geometric centre of the left end. Herein an electric field is applied between the top (z = h/2) and bottom (z = -h/2) surfaces. The bottom surface is grounded.
Caption: FIGURE 2: The first natural frequency of the present beam model ([[omega].sub.f]) without voltage and the traditional Euler-Bernoulli beam model ([[omega].sub.0]) versus beam length with different boundary conditions.
Caption: FIGURE 3: Length-dependence of the normalized fundamental frequency calculated for (a) [[omega].sub.p]/[[omega].sub.0] of PMBs with only piezoelectric effect and (b) [[omega].sub.f]/[[omega].sub.0] with both effects of piezoelectricity and flexoelectricity.
Caption: FIGURE 4: l-dependence of [[omega].sub.p]/[[omega].sub.0] obtained for PMBs where the effect of piezoelectricity is considered but the effect of flexoelectricity is ignored. Here PMBs have h = 10 [micro]m and L/h = 20 and 40, respectively.
Caption: FIGURE 5:1-dependence of ([[omega].sub.f] - [[omega].sub.p])/[[omega].sub.0] obtained for PMBs where h is fixed at 10 [micro]m and L/h is equal to 20 and 40, respectively.
TABLE 1: Material constants of the piezoelectric microbeam [8-11]. E v [rho] [e.sub.31] [[mu].sub.31] 150 GPa 0.33 7.5 x [10.sup.3] -4.35 C/ 10 [mu]C/m Kg/[m.sup.3] [m.sup.2] E [[epsilon].sub.r] 150 GPa 2300 E is Young's modulus, v is Poisson's ratio, [rho] is mass density, [e.sub.3l] is piezoelectric constant, is flexoelectric constant, [[epsilon].sub.r] is relative dielectric constant, and the vacuum dielectric constant [[epsilon].sub.0] = 8.85 x [10.sup./12] F/m.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Li, Xingjia; Luo, Ying|
|Publication:||Shock and Vibration|
|Date:||Jan 1, 2017|
|Previous Article:||Effects of Heating Rate on the Dynamic Tensile Mechanical Properties of Coal Sandstone during Thermal Treatment.|
|Next Article:||Period-Doubling Bifurcation of Stochastic Fractional-Order Duffing System via Chebyshev Polynomial Approximation.|