Abstract

A model of a segmented electrode multilayer cantilever piezoelectric actuator was established to predict its actuation performance, and then, theoretical and numerical analyses of the strain nodes were performed based on normalized deflection and strain distributions. The segmented electrodes instead of the continuous electrodes are applied in a multilayer cantilever piezoelectric actuator which can avoid the modal displacement offsets at the high vibration modes, thereby enhancing the tip deflection. The theoretical analysis and simulation results show that the tip deflection of the segmented electrode at the second mode was almost 100% larger than that of the continuous electrode. At the second mode, the maximum error between the theoretical calculation value of the tip deflection and the simulation result is 6.8%. It is because the segmented electrode is optimally designed at the strain node, which avoids the modal displacement offsets of a multilayer cantilever piezoelectric actuator at the high vibration modes; meanwhile, the theoretical results are closer to the FEM simulation results. It reveals that the tip deflection of a multilayer cantilever piezoelectric actuator can be precisely estimated by the proposed model. This research can provide some useful guidance improving the actuation performance and optimizing the design of a multilayer cantilever piezoelectric actuator.

1. Introduction

Piezoactuators are widely used in microelectromechanical systems (MEMS) because of the characteristics of small size, thinness, and high displacement, such as atomic force microscopes [1], biosensors [2], microelectromechanical switches [3], and micropositioning platforms [4], etc. The actuation performance improvement and design optimization of such devices have always been the main focus of many researchers. In particular, there are many reports in the theoretical research of piezoelectric actuators. Based on the Euler model, Wang et al. [5] and Zhang et al. [6] presented the governing equations for the piezoelectric actuators with a sandwich layer. In a study conducted by Zhang et al. [7], a simple MCPAs distributed parameter model is developed to simulate the fundamental wave of piezoelectricity in thickness-extension mode. In order to reduce the poor piezoelectric effect caused by the damage of piezoelectric materials, some researchers have designed multilayer piezoelectric actuators to improve the flexibility and compactness of the structure. Afonin [8] constructed a generalized structural parameter model of nanomechatronics multilayer electromagnetic elastic actuators. Shivashankar and Gopalakrishnan [9] reported a d33 mode surface-bondable multilayer actuator that can provide large braking force and stroke for driving large, thick, and stiffer structures. Peng et al. [10] proposed a piezoelectric multilayer actuator considering buffer layers and analyzed the dependence of the resonance frequency at the first mode and tip deflection on different layer thicknesses (buffer layer, electrode layer, and substrate layer). The contributions of the above research mainly focused on the first mode while ignoring the other higher-frequency modes, because the tip deflection of the cantilever beam at higher modes is smaller than the tip deflection at the first mode. However, piezoelectric cantilever beams can also provide superior performance at higher modes [11], which is more common in energy harvesters [12]. Ly et al. developed a 31-effect piezoelectric bending cantilever based on the Euler-Bernoulli beam theory. Their results indicated that the voltage and bandwidth at the second mode of resonant frequency were much larger than those at the first mode [13]. Except for the first mode, the cantilevered piezoelectric energy harvester (PEH) has fixed strain nodes in other vibration modes, and there are dynamic strain distributions with opposite strain signs on both sides of the node. Theoretical and experimental results demonstrate that covering the strained nodes with continuous electrodes may lead to a strong cancelation of the electrical output [14]. In order to improve the piezoelectric performance, Zizys et al. investigated the segmentation of a vibration-shock cantilevered PEH working in higher transverse vibration modes [15]. Rafique et al. used segmented electrodes to enhance the output power of the PEH [16]. Liu et al. employed the first and second bending vibration modes to design a novel bonded-type piezoelectric actuator, which obtained higher power density than previous designs [17]. Although segmented electrodes could be applied to MCPAs at the high vibration modes, there is currently no complete electromechanical equation that can describe MCPA with segmented electrodes (MCPA-S).

In this paper, by optimizing the structure of segmented electrodes, a MCPA with segmented electrodes based on the strain nodes is designed to improve the actuation performance [18]. The strain nodes are determined by the normalized deflection and strain distribution [14]. Based on Euler-Bernoulli beam theory and piezoelectric constitutive equation, a complete electromechanical coupling model is developed for the MCPA-S. Here, the electrodes are connected in series at the first mode and in parallel at the second mode, which is different from those reported in [11]. It can prevent the displacement offset in the electric potential, increase the modal electromechanical coupling term, and improve the applied capacity. In order to understand whether segmented electrodes can eliminate the influence of strain node at higher modes, we have studied the relationship of tip deflection of segmented and continuous electrodes MCPAs with excitation frequency, excitation voltage, and beam length under different modes. In addition, the structural parameters of the MCPAs were optimized by simulating different thicknesses of the substrate, piezoelectric, and buffer layers, as well as the different Young’s modulus ratios. The proposed model and prediction results can provide useful guidance for optimizing the construction and efficiency of MCPAs.

2. Design and Modeling

2.1. Design

Except for the fundamental modes, the dynamic strain distribution of the cantilevered beam changes direction at fixed strain nodes. The modal actuation capability of a cantilevered beam is closely related to the position of the piezoelectric actuators. To increase the driving force, the use of segmented actuators to control adaptive structures is proposed by avoiding the position of dynamic strain phase changes. At high modes, when the top surface of the entire piezoelectric layer was covered by continuous electrodes, the actuation capability had been significantly reduced. Therefore, we apply electrode segmentation at the nodes to a multilayer actuator (considering buffer and electrode layers), which is different from the traditional sandwich structure. This paper takes the MCPA in the second mode as an example to analyze the actuation capability. There is one strain node in the second-order mode, and the electrode is cut at the node and divided into two sections of electrodes.

Figure 1 depicts a two-dimensional schematic diagram of the MCPA-S. One end of the multilayer cantilever piezoelectric actuator is attached to the base composed of five different layers from bottom to top: the substrate, the buffer, the second electrode, the piezoelectric layer, and the first electrode as shown in Figure 1(a). The first electrode and the piezoelectric layer are cut at the strained node to form segmented electrodes for the MCPA. The contact between the second electrode and the lower surface of the piezoelectric film is continuous, but the contact between the first electrode and the upper surface of the piezoelectric layer is discontinuous. The polarization direction is reversed after passing through the strain node, and the split position is the strain node position L1. The different electrode connections are adopted under the different modes of the cantilever beam. At the first mode, the strain distribution is in the same phase for the MCPA, because there is no strain node [10]. Here, the electrode directions in the L1 and L1-L regions are opposite, and the electrode wires are connected in series with the applied voltage, as shown in Figure 1(b). At the second mode, there is a strain node and the strain distributions in L1 and L1-L regions are 180 degrees out of phase. The connection of the electrodes is described in Figure 1(c), which is arranged in parallel to prevent modal displacement in the electric potential. And then the applied voltage should be applied on the L1 and L1-L regions of the cantilever beam to generate the deflection for the MCPA.

Z1, Z2, Z3, and Z4 depict the vertical coordinates of the bottom-surface of substrate, buffer, piezoelectric layer, and the first electrode, respectively. Zi and Z5 indicate the vertical coordinates of the top surface of buffer and the second electrode. The length and width of the cantilever beam are denoted by L and b. In Cartesian coordinate system, the x- and z-axes are consistent with the directions 1 and 3, respectively, and the z-axis represents the polarization direction of the piezoelectric layer. The coordinate origin of the x-z plane corresponds to the leftmost point of the MCPA. The mid-plane of the substrate is denoted by the dotted line. The neutral plane is located at z0 from the mid-plane of the substrate. Moreover, h is used to describe the thickness of each layer, and its subscripts , and indicate the piezoelectric, substrate, buffer, and electrodes layers, respectively. The transverse deflection of the MCPA occurs along the z-axis and is a function with the x value and time t.

2.2. Modeling
2.2.1. Electromechanical Coupling Equation

The constitutive equation for piezoelectric materials can be described in the form [10]where stress, strain, and electric field are denoted by T, S, and E, respectively, elastic constant of the piezoelectric material is described by c, and piezoelectric coupling coefficient under steady electric field is depicted by ep. Here, the axial strain and polarization direction are marked as subscripts 1 and 3, respectively. The elastic stiffness component can be calculated by c11,p = 1/s11,p according to the plane-stress presumption of the MCPA. Under a constant electric field, s11,p denotes the elastic compliance. In addition, e31 can be expressed as ep31 = d31/s11,p with the commonly used piezoelectric constant d31. The constitutive equations applied to the substrate layer and the buffer layer can be given as [19]

The constitutive equation for the electrode layer is [19]and the axial strain at value x and time t can be obtained by [20]

According to our previous research [7], the position of the neutral plane z0 can be given by

According to the moment balance equation (the beam’s cross section), the bending moment is given by

At the first mode, the polarization and electric field are consistent along the entire beam length; thus, the uniform electric field E3(t) can be given in terms of voltage across the piezoelectric layer and the thickness hp as . At the second mode, the polarization direction and electric field of the two segmented piezoelectric layers are different. Therefore, for the 0-L1 area (same direction) the electric field can be obtained by ; and for the L-L1 area (opposite direction), .

By substituting equations (1)–(4) into equation (5) and integrating equation (5), it can be further simplified to

EI represents the bending stiffness of the composite structure’s cross section, which is expressed by the following equation:

represents the spatial distribution of the electric potential and is related to the modes and the structure of the piezoelectric cantilever beam. For the MCPA with continuous electrodes (MCPA-C), is expressed as [21]

For the MCPA with segmented electrodes (Figure 1), the uniform electric fields are applied at different modes by conducting the different electrode lines. The potential spatial distribution can be regarded as the sum of the two electrode regions. At the first mode, the electric field is exerted to the segmented electrodes in series, and the can be derived by

At the second mode, the electric field is exerted to the segmented electrodes in parallel, and the can be derived bywhere H(x) denotes the Heaviside function and the coupling term is given as

The configurations of the microcantilever conform to the Euler-Bernoulli beam hypothesis, which have been presented in our previous research [10]. Considering viscous air (medium) damping and Kelvin-Voigt (or strain rate) damping, the governing equation of the cantilever beam can be written as [21]

Here, the viscous damping coefficient is denoted by ca, the inertia moment of the cross section area is described by I, the strain rate damping term is expressed as csI, and the mass per unit length of the MCPA is represented by m, which is obtained bywhere the uniform densities of the different layers are ps (Si substrate), pi (buffer), pe (electrodes), and pp (piezoelectric), respectively. The mass per unit length of the different layers is bps (Si substrate), bpi (buffer), bpe1 (electrodes), bpp (piezoelectric), and bpe2 (electrodes), respectively. Finally, by inserting equation (6) into equation (11), the electromechanical coupling equation of the MCPA at the first two modes can be obtained:

2.2.2. Modal Analysis

Based on the standard modal expansion approach, a series of the absolutely uniformly convergent eigenfunctions are used to describe the transverse deflection of the cantilever beam [11]:

The mass normalized eigenfunction and modal coordinate of the clamped-free beam at the rth vibration mode are, respectively, given by and . The deflection of the cantilever beam can be exactly obtained from this equation. is written as [11]and it satisfies the orthogonality conditions [9].

Here, (dimensionless frequency number) of the rth vibration mode can be expressed as follows:and is given by

Bending strain distribution can be measured directly by the curvature eigenfunction that is the second derivative of the displacement eigenfunction (equation (16)). For a positive definite system (), the positions of the strain nodes can be determined by calculating the roots of equation (20) at where denotes the length position (dimensionless) for the MCPA. By combining equation (20) with equations (18) and (19), the strain nodes positions (dimensionless) of the first three modes can be obtained in Table 1.

In addition, is the undamped natural frequency of the rth mode, which is written as

Equation (16) is simplified by using the orthogonal condition of equation (17) and then substituted into equation (15). The mechanical motion equation in modal coordinates can be derived as follows:where denotes the modal mechanical damping ratio. The coupling term (modal electromechanical) is defined as

(23) can be further rewritten as [11]where is the Dirac function.

By substituting equation (8) into equation (23), for the continuous electrodes can be rewritten as

For the segmented electrodes, is related to the spatial distribution of the electric potential at the vibration modes. By substituting equation (9) into equation (25), at the first mode can be expressed as

At the second mode, it can be obtained by

Employing the separating variables method, we record and , where and represent the amplitudes. By substituting them into equation (22), can be calculated by

By substituting equation (28) into equation (15), the transverse deflection can be redescribed as the following formula:

Finally, the tip deflection that occurred at the free end of the cantilever is expressed as

3. Verification by FEM Simulation

3.1. Material Properties and Structural Parameters

The MCPA finite element model consists of Si substrate, SiO2 buffer, Pt first electrode, piezoelectric, and Pt second electrode. The geometrical dimensions of the segmented electrode MCPA include length L, width b, thicknesses hs, hi, he1, hp, he2, and segmented length L1. All dimensions are listed in Table 2. At the first and second mode, the mechanical damping ratios were and , respectively. Poisson’s ratio was set in this paper, and other material property parameters are shown in Tables 3 and 4.

3.2. Strain Distribution of the Cantilevers with MCPA

In the simulations, the piezoelectric material is modeled by “solid5” composed of 3D 8-node hexahedral coupled-field elements, and the nonpiezoelectric materials are modeled by “solid45” including 8-node linear structural elements. At the beam’s fixed end, the freedom degree of displacement is limited to be zero. The electrode connection of the first and second electrodes is implemented using coupling commands. For the upper surface of the second electrode layer and the lower surface of the piezoelectric layer, the voltage is coupled and constrained to be zero. The applied voltage is coupled to the upper surface of the first electrode layer. Figures 2(a), 2(b), 3(a) and 3(b) show the raw and fine mesh of the finite element (FEM); the circled area marks the position of the strain node. The strain distributions of the continuous and segmented electrode MCPAs at the first two modes are, respectively, presented in Figures 2(c), 2(d), 3(c), and 3(d). The black line of MCPA is the cutting position of the segmented electrode (L1), as shown in Figure 3. It can be observed that, at the first mode, the strain of the continuous and segmented electrodes decreases monotonously from the fixed end to the free end in the strain contours of the MCPAs. However, at the second mode, there is a minimum strain magnitude at a certain region of the beam, which is much lower than at the fixed end, indicating the presence of a strain node.

4. Results and Discussion

4.1. Determination of Strain Nodes

Figure 4 shows the dependence of normalized deflection and normalized strain distribution of the MCPAs with continuous and segmented electrodes on the dimensionless position along the beam axis . Here, the normalized deflection curve is monotonically decreasing, and there is no zero point at the first mode, as shown in Figure 4(a). At the second mode there is a zero point on the normalized deflection curve, indicating the existence of strain nodes. In Figure 4(b), at the first mode, the strain distribution curve is monotonically increasing without strain node. At the second mode, there is a strain node at , which is similar to the reported result [11]. Therefore, the existence of strain nodes at the second mode is analyzed by both the theoretical and simulation models.

4.2. Effect of Segmented Electrode Length on Tip Deflection

To verify the influence of the segmented electrode length L1 on the tip deflection of the MCPA, the L1-tip deflection curves under different resonance frequency were simulated when V = 1–5 V, L1 = 0–1.0 mm, and L = 1 mm, as shown in Figure 5. Figure 5(a) shows that tip deflection at the first mode does not change with the increase of L1 under constant applied voltages. Figure 5(b) indicates that, at the second mode, as L1 increases from 0 to 0.216 mm, the tip deflection increases monotonically, but as L1 increases from 0.216 to 1 mm, the tip deflection decreases nonmonotonically. The tip deflection increases and reaches a peak at L1 = 0.216 mm (strain node); then the tip deflection decreases and reaches the lowest value at L1 = 0.620 mm. Meanwhile, the maximum value of tip deflection increases as the applied voltage increases. At different modes, the electrode is segmented at the strain node and the electrodes are connected in different ways, in which the modal electromechanical coupling coefficient can reach a large value. For the modal displacement in the spatial potential, the cancelation is prevented to increase the tip deflection at the second mode. It indicates that the mechanism of the segmented electrode to avoid the modal displacement offsets at the high modes has been verified for the MCPAs. By adjusting the length of the segmented electrode and positioning segmentation location at the strain node, the larger tip deflection value can be obtained to improve the actuation performance.

4.3. Tip Deflections with the Segmented Electrode under Different Excitation Frequencies

Under different excitation frequencies, the tip deflection at the first and second modes of the MCPA with continuous and segmented electrodes is described in Figures 6 and 7. The tip deflection reaches a peak at the resonance frequency, and the peak value increases with the growth of the applied voltages, as shown in Figure 6. For continuous electrodes and segmented electrodes MCPAs, the tip deflections at the first mode are almost the same. At the second mode, the tip deflection of the MCPA-S at the second mode is almost 100% larger than that of the MCPA-C, as depicted in Figure 7. The maximum error of the tip deflection at the second mode is 6.8% between the theoretical and simulation results, and the theoretical calculation results of tip deflection are close to its simulation results. It reveals that, under different excitation frequencies, there is zero/one strain node at the first/second mode for the MCPAs, which has different degrees of influence on the tip deflection. Furthermore, the series/parallel connection is valid for MCPA-S at the first/second mode. At the second mode, the dynamic strain distribution of the beam can change the strain direction on both sides of the strain node [11]. When the strain nodes are covered by the continuous electrodes, the tip deflection of the beam is canceled. Therefore, the optimized segmented electrode is cut at the strain node, and the wires of the segmented electrodes are connected in series at the first mode and in parallel at the second mode. It can prevent the modal displacement in the electric potential from canceling out [21] and improve the actuation performance of MCPA-S.

4.4. Dependence of Tip Deflection on Applied Voltage

To perceive the dependence of tip deflection on applied voltage, theoretical calculations and simulations are carried out on the tip deflection of MCPA-C and MCPA-S under different applied voltages. Here, the tip deflection increases linearly with the increase of the applied voltage from 1 to 5 V, as indicated in Figure 8. The theoretical slope of the applied voltage-tip deflection curve at the first mode is 12.85 μm/V at 24 kHz in Figure 8(a). In Figure 8(b), at 150 kHz, the theoretical slopes of the applied voltage-tip deflection curves at the second mode are 0.88 μm/V and 1.73 μm/V, respectively, for MCPA-C and MCPA-S. The theoretical and simulation tip deflections of the MCPA-S are 97% larger than that of MCPA-C, and the maximum relative error is 7.8%. Obviously, the theoretical values are in good agreement with the simulation values. It shows that the presented theoretical model can correctly estimate the dynamic tip deflection of the MCPA under the applied voltage.

4.5. Dependence of the Tip Deflection on the Beam Length

In order to investigate the relationship between the tip deflection and the beam length [18], the values of tip deflection of the MCPA were obtained when L = 0–2.0 mm and V = 1–5 V. The theoretical/simulation results of the L-tip deflection curves of MCPA-C and MCPA-S at first and second mode are given in Figures 9 and 10. Under different constant voltages, as the beam length L increases from 0 to 2 mm, the tip deflections at the first and second modes increase nonlinearly, which is consistent with [22]. At the second mode, the tip deflection of the segmented electrode is obviously greater than that of the continuous electrode, as shown in Figures 10(a) and 10(b). By analyzing the relative deviations of the tip deflection when V = 1 V (low applied voltage), the theoretical and simulation results are expressed by lines and dots, respectively, as shown in Figure 10(c). The maximum deviations of the tip deflections between the theoretical and simulation results are 4.9% for the MCPA-S and 4.1% for the MCPA-C. It is similar to the reported results [21]. The validity of the proposed model has been verified for segmented electrode MCPA. Therefore, it is a useful strategy to improve the actuation performance by adjusting the beam length L, which can be used to drive microelectromechanical switches [3].

4.6. Dependence of Tip Deflection on the Substrate/Piezoelectric Layer Thickness Ratio

In this section, r = Es/Ep and h = hs/hp are defined to describe Young’s modulus ratio and thickness ratio between the substrate and the piezoelectric layer [23]. In order to investigate dependence of tip deflection on the thickness ratio between the substrate and piezoelectric layer, when V = 1 V, he = 1 μm, L = 1 mm, b = 200 μm, and hi = 1 μm, the h-tip deflection curves of MCPAs at the first two modes under different Young’s modulus ratios were simulated and analyzed. They are shown in Figures 11 and 12. Two monotonously changing regions are formed on two sides of the maximum point of each curve. At the first and second modes, as h and r increase, the tip deflections increase/decrease in the upward/downward region. In particular, when the thickness ratio is constant (h = 0–2.3), the tip deflection increases with the increase of Young’s modulus (r = 0.8–2.4). In the downward region, when h is constant (h = 2.3–5), the tip deflection decreases as r increases (r = 0.8–2.4). For MCPAs at the first and second mode, the change tendencies of tip deflections with h and r are similar through theoretical calculation and simulation. The tip deflection of MCPA-S is good with that of MCPA-C at the first mode. Obviously, with a smaller thickness ratio and a larger Young’s modulus ratio, the curve in the upward region has a larger slope; i.e., in this region the tip deflection is more sensitive to the thickness change, which is beneficial for position sensing [24]. At the second mode, the slopes of the tip deflection curves for MCPA-S are significantly larger than that of the curves for MCPA-C, as shown in Figure 12. It indicates that the tip deflection of MCPA-S is more sensitive to change of thickness ratio and Young’s modulus ratio. A smaller thickness ratio and larger Young’s modulus ratio can achieve the maximum tip deflection, thereby improving the actuation performance of the MCPAs [25].

5. Conclusions

In summary, by considering the influence of the strain node on the beam at the high modes for the MCPAs, the optimal position of segmented electrode should be at the strain node. At the high modes of MCPAs, the tip deflections of the MCPA-S are larger than that of the MCPA-C, so the actuation performance can be improved. The strain node of the MCPA-S is 0.216 mm at the second mode, and the developed complete electromechanical coupling model can predict the tip deflection at the first two modes. The results of theoretical calculation and simulation indicate that the tip deflections of the MCPA-S are consistent with the MCPA-C at the first mode. The tip deflection of the MCPA-S is almost 100% larger than that of the MCPA-C at the second mode, and the tip deflection at the second mode has a maximum error of 6.8% between the theoretical and simulation results. Obviously, the theoretical values and the simulation results are relatively close, indicating that the proposed model can precisely estimate the tip deflection of MCPAs. The reliance of tip deflection on segmented electrode length, actuation frequency, voltage, beam length, substrate/piezoelectric thickness, and Young’s modulus ratio are also discussed. The results indicate that, for the MCPA-S, under a certain beam length and high voltage pressure, a smaller thickness ratio and a larger Young’s modulus ratio of the substrate/piezoelectric layer are beneficial to gain a larger tip deflection. The proposed model verifies the mechanism of MCPA-S to avoid the modal displacement offsets at the high modes, and predicted results can provide valuable guidance for optimizing the construction and efficiency of MCPAs.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest with respect to the research, authorship, and/or publication of this article.

Acknowledgments

This work was supported by NNSF of China (11832016 and 51775471), Hunan Innovative Province Construction Special Major Landmark Innovation Demonstration Project, Changsha Zhuzhou Xiangtan Landmark Engineering Technology Project (2019XK2303 and 2020GK2014), Hefei General Machinery Research Institute Co., Ltd Project (2021ZKKF043), and Postgraduate Scientific Research Innovation Project of Hunan Province (CX20200642).