**1. Introduction**

Piezoelectric polymers have been widely used in sensors, actuators, electronic information and intelligent structures because of its grea<sup>t</sup> force–electric coupling characteristics [1–6]. The piezoelectric polymers usually participate in the work of piezoelectric instruments in the form of piezoelectric sheets which usually are simplified to a piezoelectric cantilever beam [7–9].The problems of piezoelectric cantilever beams are usually difficult to be solved analytically due to the existence of the force–electric coupling constitutive relation. It is known that the design of piezoelectric instruments often requires the analytical expression of the problem of piezoelectric cantilever beams as a theoretical reference. Therefore, it is necessary and meaningful to find an efficient analytical method for solving the problem of piezoelectric cantilever beams and giving their analytical solutions.

In the past twenty years, many researchers have studied the problem of piezoelectric cantilever beams and obtained some corresponding solutions. Wang and Chen [10] obtained a general solution of the control equation for the three-dimensional problem of transverse isotropic piezoelectric material by means of a set of new potential functions representing displacement component and potential function, and solved the problem of spatial piezoelectric material under the action of concentrated transverse shear force. Lin et al. [11] derived the analytical expressions of displacement, potential, and stress

distribution of piezoelectric beams which were simply supported at both ends under a uniform load. According to the plane stress problem, Mei and Zeng [12] directly derived the equation of state of piezoelectric beams from the piezoelectric physical equation, and on this basis, the exact state equation solution of electromechanical coupling e ffect of simply supported piezoelectric beams at both ends under a uniform load was given. On the basis of three-dimensional constitutive equations and their simplified equations of elastic piezoelectric materials, Zhu [13] derived the analytic solution to a piezoelectric cantilever beam with concentrated force at the free end in terms of displacements and voltage. For the orthotropic piezoelectric plane problem, Ding et al. [14–16] solved a series of piezoelectric beam problems and obtained the corresponding exact solutions with the trial and error method on the basis of the general solution in the case of three distinct eigenvalues, and expressed all displacements, electrical potential, stresses, and electrical displacements by three displacement functions in terms of harmonic polynomials. Yang and Liu [17] investigated the bending of transversely isotropic cantilever beams under an end load, and derived the simplified linear elastic equations of piezoelectric cantilever beams according to the characters of the problem. Pang et al. [18] manufactured a typical Li- and Ta/Sb-modified, alkaline niobate-based, lead-free piezoelectric ceramics by two-step sintering and investigated the sintering condition dependence of dielectric constants and piezoelectric properties. Zhu et al. [19] studied the active vibration control of piezoelectric cantilever beams, where an adaptive feed forward controller (AFC) was utilized to reject the vibration with unknown multiple frequencies. Peng et al. [20] presented time-delayed feedback control to reduce the non-linear resonant vibration of a piezoelectric elastic beam and examined three single-input linear time-delayed feedback control methodologies: displacement, velocity, and acceleration time-delayed feedback. Liu and Yang [21] studied the bending problem of a cantilever beam made of a transversely isotropic piezoelectricity medium under uniformly distributed loads. Shi et al. [22,23] studied the analytical solution of a density functionally gradient piezoelectric cantilever under axial and transverse uniform loads and applied DC voltages and then, solved the force–electric coupling plane strain problem of simply supported beams under a uniform load by the inverse method. Wang et al. [24] dealt with the vibration analysis of a circular plate surface bonded by two piezoelectric layers, based on the Kirchho ff plate model. Recently, Lian et al. [25] studied the problem of a functionally graded piezoelectric cantilever beam under combined loads, but non-dimensionalization was not considered in solving the problem. There is still a lot of research performed in this field, which will not be elaborated here. The summation of results of existing research shows that there are still some unsolved problems. First, non-dimensionalization was not considered in the existing research. We know that piezoelectric materials have not only mechanical properties, but also electrical properties. So, there are both mechanical units and electrical units to be solved, which may lead to computational errors. Second, the existing research basically provides theoretical solutions, but there are a few related experimental verifications. Therefore, the reliability of theoretical solutions cannot be guaranteed. Besides, there has been no unified and e ffective method for solving the problems of piezoelectric structure.

Parameter perturbation method is a general analysis method for solving approximate solutions of non-linear mechanical problems. It has been successfully applied to various fields of non-linear structural analysis, such as non-linear bending and post-buckling, and has become a powerful tool for solving non-linear problems of structures. Generally speaking, the perturbation method is based on a selected small parameter. In order to solve the problem of parameter selection, Chen and Li [26] put forward the concept of free parameter perturbation method, that is, there is no need to point out the physical meaning of perturbation parameters during perturbation, which provides a new idea for solving the parameter selection problem of parameter perturbation method. Lian et al. [27] solved the Hencky membrane problem without a small-rotation-angle assumption by the single-parameter perturbation method. The successful application of perturbation method depends, to a large extent, on the reasonable choice of small parameters, but the selection of perturbation parameters does not have a set of step-by-step procedures, which can only rely on deep understanding and multiple attempts. To avoid the di fficulty in the selection of perturbation parameters, researchers can select multiple parameters, that is, the so-called "multi-parameter perturbation method". For multi-parameter perturbation method, Nowinski and Ismail [28] solved the cylindrical orthotropic circular plate problem under a uniform load by using the two-parameter perturbation method. The application of the multi-parameter perturbation method in beam problem was proposed by Chien [29] in 2002, the classical Euler–Bernoulli equation of bending beams was solved by using load and beam height differences as perturbation parameters. Later, He and Chen [30] simplified the bending moment by using the quasi-linear analysis method, so that the parameter perturbation process was directly aimed at the algebra equation rather than the integral equation, and the two-parameter perturbation solution of the large deflection bending problem of a cantilever beam was obtained, and the integrity of the two-parameter perturbation solution was analyzed. Recently, He et al. [31,32] comprehensively analyzed the large deflection problem of beams with height difference under various boundary conditions, put forward the so-called "two-parameter perturbation method", and successfully applied this method to the solution of bimodular von-Kármán thin plate equation. But so far, the perturbation method of three or more parameters has only a few reports.

In this study, we will derive the theoretical solution of the bending problem of piezoelectric cantilever beams by the multi-parameter perturbation method. The whole paper is organized as follows. In Section 2, the mechanical model of the problem solved here will be established, and the governing equations will be given and dimensionless. In Section 3, the three piezoelectric parameters will be selected as perturbation parameters, and the dimensionless governing equations will be solved by the multi-parameter perturbation method. The solution presented in this paper will be compared with the existing analytical solution from Yang and Liu [17] in Section 4. Next, in Section 5, we will show the related experiments of the piezoelectric cantilever beam, compare the experimental results with the solution presented here, and also discuss the effect of the piezoelectric properties on the deformation of piezoelectric cantilever beams. According to the results mentioned above, some main conclusions will be drawn in Section 6.

#### **2. Mechanical Model and Basic Equations**

In this study, the mechanical model of the transversely isotropic piezoelectric cantilever beam is established by using two-dimensional elastic beam theory and neglecting shear deformation. As shown in Figure 1, an transversely isotropic piezoelectric cantilever beam is fixed at its right end and subjected to a uniformly distributed load *q* on its upper surface, a concentrated force *P* and a bending moment *M* at its left end, in which *l*, *b*, and *h* denote the length, width, and height of the beam, respectively, and *O* denotes the origin of the coordinates. A rectangular coordinate system is introduced with the upper and lower surfaces of the beam lying in *z* = −*h*/2 and *z* = *h*/2.

**Figure 1.** Scheme of a piezoelectric cantilever beam.

*Polymers* **2019**, *11*, 1934

Supposing that the polarization direction is the forward direction of the *z*-axis, let us take a microelement in the piezoelectric cantilever beam, and from the balance of the force, we may obtain, by neglecting the body force

$$\begin{cases} \frac{\partial \sigma\_x}{\partial x} + \frac{\partial \tau\_{zx}}{\partial z} = 0\\ \frac{\partial \tau\_{zx}}{\partial x} + \frac{\partial \sigma\_z}{\partial z} = 0 \end{cases} \tag{1}$$

where σ*x*, σ*z* and τ*zx* are the stress components. The equation of Maxwell electric displacement conservation is

$$
\frac{\partial D\_x}{\partial \mathbf{x}} + \frac{\partial D\_z}{\partial z} = 0,\tag{2}
$$

where *Dx* and *Dz* are the electric displacement components. The constitutive equations of piezoelectric polymeric materials considered are

$$\begin{cases} \varepsilon\_{x} = s\_{11}\sigma\_{x} + s\_{13}\sigma\_{z} + d\_{31}E\_{z} \\ \varepsilon\_{z} = s\_{13}\sigma\_{x} + s\_{33}\sigma\_{z} + d\_{33}E\_{z} \\ \gamma\_{zx} = s\_{44}\tau\_{zx} + d\_{15}E\_{x} \\ D\_{x} = d\_{15}\tau\_{zx} + \lambda\_{11}E\_{x} \\ D\_{z} = d\_{31}\sigma\_{x} + d\_{33}\sigma\_{z} + \lambda\_{33}E\_{z} \end{cases} \tag{3}$$

where ε*x*, ε*z*, and γ*zx* are the strain components; and *Ex* and *Ez* are the electric field intensity components. The geometric equations of the piezoelectric cantilever beam are

$$
\varepsilon\_x = \frac{\partial u}{\partial \mathbf{x}} , \varepsilon\_z = \frac{\partial w}{\partial z} , \gamma\_{z\mathbf{x}} = \frac{\partial u}{\partial z} + \frac{\partial w}{\partial \mathbf{x}} , \tag{4}
$$

where *u* and *w* are the displacement components. From Equation (4), the strain consistency equation is obtained as follows:

$$
\frac{
\partial^2 \varepsilon\_x
}{
\partial z^2
} + \frac{
\partial^2 \varepsilon\_z
}{
\partial x^2
} - \frac{
\partial^2 \gamma\_{zx}
}{
\partial z \partial x
} = 0.
\tag{5}
$$

The relationship between electric field intensity and electric potential are

$$E\_{\mathbf{x}} = -\frac{\partial \phi}{\partial \mathbf{x}} , E\_z = -\frac{\partial \phi}{\partial z} , \tag{6}$$

where φ is the electric potential function. By introducing the Airy stress function *<sup>U</sup>*(*<sup>x</sup>*, *<sup>z</sup>*), the stress components can be expressed as

$$
\sigma\_{\mathbf{x}} = \frac{\partial^2 \mathcal{U}}{\partial \mathbf{z}^2}, \sigma\_{\mathbf{z}} = \frac{\partial^2 \mathcal{U}}{\partial \mathbf{x}^2}, \tau\_{\mathbf{z}\mathbf{x}} = -\frac{\partial^2 \mathcal{U}}{\partial \mathbf{z}\partial \mathbf{x}}.\tag{7}
$$

The boundary conditions of the problem of the piezoelectric cantilever beam are

$$\begin{cases} \int\_{-h/2}^{h/2} \tau\_{x\overline{x}} dz = \int\_{-h/2}^{h/2} \frac{\partial^2 \mathcal{U}}{\partial z \partial x} dz = -\frac{\mathcal{P}}{\mathcal{b}},\\ \int\_{-h/2}^{h/2} \sigma\_x dz = \int\_{-h/2}^{h/2} \frac{\partial^2 \mathcal{U}}{\partial z^2} dz = 0, \quad \text{, at } x = 0, \\ \int\_{-h/2}^{h/2} z \sigma\_x dz = \int\_{-h/2}^{h/2} z \frac{\partial^2 \mathcal{U}}{\partial z^2} dz = \frac{\mathcal{M}}{\mathcal{b}} \end{cases} \tag{8}$$

$$
\sigma\_z = \frac{\partial^2 \mathcal{U}}{\partial \mathbf{x}^2} = 0,\\
\tau\_{z\mathbf{x}} = -\frac{\partial^2 \mathcal{U}}{\partial \mathbf{x} \partial z} = 0,\text{ at } z = \hbar/2,\tag{9}
$$

$$
\sigma\_z = \frac{\partial^2 \mathcal{U}}{\partial \mathbf{x}^2} = -q\_\prime \tau\_{z\mathbf{x}} = -\frac{\partial^2 \mathcal{U}}{\partial \mathbf{x} \partial \mathbf{z}} = 0,\text{ at }\mathbf{z} = -\hbar/2,\tag{10}
$$

$$\int\_{-h/2}^{h/2} D\_x dz = 0, \text{ at } \mathbf{x} = 0 \text{ and } \mathbf{x} = l\_\prime \tag{11}$$

*Polymers* **2019**, *11*, 1934

$$E\_z = \frac{\partial \phi}{\partial z} = 0, \text{ at } z = \pm h/2,\tag{12}$$

and

$$u = 0, w = 0, \frac{\partial w}{\partial \mathbf{x}} = 0, \text{ at } z = 0 \text{ and } \mathbf{x} = l. \tag{13}$$

Substituting Equations (3), (6), and (7) into Equations (2) and (4), we may obtain two equations of the stress function *<sup>U</sup>*(*<sup>x</sup>*, *z*) and the potential function φ

$$d\_{31}\frac{\partial^3 \mathcal{U}}{\partial z^3} + d\_{33}\frac{\partial^3 \mathcal{U}}{\partial x^2 \partial z} - d\_{15}\frac{\partial^3 \mathcal{U}}{\partial x^2 \partial z} = \lambda\_{33}\frac{\partial^2 \phi}{\partial z^2} + \lambda\_{11}\frac{\partial^2 \phi}{\partial x^2} \tag{14}$$

and

$$s\_{11}\frac{\partial^4 \mathcal{U}}{\partial z^4} + (2s\_{13} + s\_{44})\frac{\partial^4 \mathcal{U}}{\partial \mathbf{x}^2 \partial z^2} + s\_{33}\frac{\partial^4 \mathcal{U}}{\partial \mathbf{x}^4} = d\_{31}\frac{\partial^3 \phi}{\partial z^3} + d\_{33}\frac{\partial^3 \phi}{\partial \mathbf{x}^2 \partial z} - d\_{15}\frac{\partial^3 \phi}{\partial \mathbf{x}^2 \partial z}.\tag{15}$$

Equations (14) and (15) are usually called governing equations. Let us introduce the following dimensionless quantities:

$$\begin{aligned} X &= \frac{x}{h}, Z = \frac{z}{h}, S\_{13} = \frac{s\_{13}}{s\_{11}}, S\_{33} = \frac{s\_{33}}{s\_{11}}, S\_{44} = \frac{s\_{44}}{s\_{11}}, \overline{d}\_{31} = \frac{d\_{31}}{\sqrt{s\_{11}\lambda\_{11}}}, \overline{d}\_{33} = \frac{d\_{33}}{\sqrt{s\_{11}\lambda\_{11}}}, \overline{d}\_{15} = \frac{d\_{15}}{\sqrt{s\_{11}\lambda\_{11}}}, \\\ \Phi &= \frac{\phi \cdot \sqrt{s\_{11}\lambda\_{11}}}{h}, \overline{\lambda}\_{33} = \frac{\lambda\_{33}}{\lambda\_{11}}, \overline{P} = \frac{p}{h^2} s\_{11}, \overline{b} = \frac{h}{h}, \overline{M} = \frac{M}{h^3} s\_{11}, \overline{\eta} = qs\_{11}, \overline{u} = \frac{\mu}{h}, \overline{w} = \frac{w}{h^2}, \overline{U} = \frac{\|L\|\_{1}}{h^2} \end{aligned} \tag{16}$$

From Equation (16), Equations (14) and (15) can be transformed into

$$
\overline{d}\_{31} \frac{\partial^3 \overline{\mathcal{U}}}{\partial Z^3} + \overline{d}\_{33} \frac{\partial^3 \overline{\mathcal{U}}}{\partial X^2 \partial Z} - \overline{d}\_{15} \frac{\partial^3 \overline{\mathcal{U}}}{\partial X^2 \partial Z} = \overline{\lambda}\_{33} \frac{\partial^2 \Phi}{\partial Z^2} + \frac{\partial^2 \Phi}{\partial X^2} \tag{17}
$$

and

$$\frac{\partial^4 \overline{\mathcal{U}}}{\partial \mathcal{Z}^4} + (2\mathcal{S}\_{13} + \mathcal{S}\_{44}) \frac{\partial^4 \overline{\mathcal{U}}}{\partial X^2 \partial \mathcal{Z}^2} + \mathcal{S}\_{33} \frac{\partial^4 \overline{\mathcal{U}}}{\partial X^4} = \overline{d}\_{31} \frac{\partial^3 \Phi}{\partial \mathcal{Z}^3} + \overline{d}\_{33} \frac{\partial^3 \Phi}{\partial X^2 \partial \mathcal{Z}} - \overline{d}\_{15} \frac{\partial^3 \Phi}{\partial X^2 \partial \mathcal{Z}}.\tag{18}$$

The boundary conditions can be transformed into

$$\int\_{-1/2}^{1/2} \frac{\partial^2 \overline{\mathcal{U}}}{\partial Z \partial X} dZ = -\frac{\overline{P}}{\overline{b}}, \int\_{-1/2}^{1/2} \frac{\partial^2 \overline{\mathcal{U}}}{\partial Z^2} dZ = 0,\\ \int\_{-1/2}^{1/2} Z \frac{\partial^2 \overline{\mathcal{U}}}{\partial Z^2} dZ = \frac{\overline{M}}{\overline{b}}, \text{ at } X = 0,\tag{19}$$

$$\frac{\partial^2 \overline{U}}{\partial X^2} = -\frac{\partial^2 \overline{U}}{\partial X \partial Z} = 0, \text{ at } Z = 1/2,\tag{20}$$

$$\frac{\partial^2 \overline{U}}{\partial X^2} = -\overline{q}\_\prime - \frac{\partial^2 \overline{U}}{\partial X \partial \overline{Z}} = 0, \text{ at } \overline{Z} = -1/2,\tag{21}$$

$$\int\_{-1/2}^{1/2} (-\overline{d}\_{15}\frac{\partial^2 \overline{U}}{\partial Z \partial X} - \frac{\partial \Phi}{\partial X}) dZ = 0, \text{ at } X = 0 \text{ and } X = l/h,\tag{22}$$

$$\frac{\partial \Phi}{\partial Z} = 0, \text{ at } Z = \pm 1/2,\tag{23}$$

and

$$
\overline{u} = 0, \overline{w} = 0, \frac{\partial \overline{w}}{\partial X} = 0, \text{ at } Z = 0 \text{ and } X = l/h. \tag{24}
$$

#### **3. Multi-parameter Perturbation Solution**

Equations (17) and (18) are two partial differential equations which are usually difficult to solve analytically. Here, we use the multi-parameter perturbation method to solve them. The piezoelectric coefficients are usually very small [33], thus, they can be selected as perturbation parameters to meet the requirement of convergence in perturbation expansions. From the point of view of the perturbation idea, if the cantilever beam without piezoelectric properties is regarded as an unperturbed system, the

piezoelectric cantilever beam can be looked upon as a perturbed system. Selecting *d*31, *d*33 and *d*15 as the perturbation parameters, the Φ and *U* can be expanded as

$$\begin{aligned} \Phi &= \Phi\_0^0 + \Phi\_1^1 \overline{d}\_{31} + \Phi\_2^1 \overline{d}\_{33} + \Phi\_3^1 \overline{d}\_{15} + \Phi\_1^{\text{II}} (\overline{d}\_{31})^2 + \Phi\_2^{\text{II}} (\overline{d}\_{33})^2 \\ &+ \Phi\_3^{\text{II}} (\overline{d}\_{15})^2 + \Phi\_4^{\text{II}} \overline{d}\_{31} \overline{d}\_{33} + \Phi\_5^{\text{II}} \overline{d}\_{31} \overline{d}\_{15} + \Phi\_6^{\text{II}} \overline{d}\_{33} \overline{d}\_{15} \end{aligned} \tag{25}$$

and

$$\begin{aligned} \overline{\boldsymbol{U}} &= \overline{\boldsymbol{U}}\_0^0 + \overline{\boldsymbol{U}}\_1^1 \overline{\boldsymbol{d}}\_{31} + \overline{\boldsymbol{U}}\_2^1 \overline{\boldsymbol{d}}\_{33} + \overline{\boldsymbol{U}}\_3^1 \overline{\boldsymbol{d}}\_{15} + \overline{\boldsymbol{U}}\_1^{\Pi} (\overline{\boldsymbol{d}}\_{31})^2 + \overline{\boldsymbol{U}}\_2^{\Pi} (\overline{\boldsymbol{d}}\_{33})^2 \\ &+ \overline{\boldsymbol{U}}\_3^{\Pi} (\overline{\boldsymbol{d}}\_{15})^2 + \overline{\boldsymbol{U}}\_4^{\Pi} \overline{\boldsymbol{d}}\_{31} \overline{\boldsymbol{d}}\_{33} + \overline{\boldsymbol{U}}\_5^{\Pi} \overline{\boldsymbol{d}}\_{31} \overline{\boldsymbol{d}}\_{15} + \overline{\boldsymbol{U}}\_6^{\Pi} \overline{\boldsymbol{d}}\_{33} \overline{\boldsymbol{d}}\_{15} \end{aligned} \tag{26}$$

where Φ00 and *U*00, ΦI*i* and *U*I*i* (*i* = 1, 2, 3), and ΦII*i* and *U*II*i* (*i* = 1, 2, ... , 5, 6) are unknown functions of *X* and *Z*.

First, we solve the zero-order perturbation equations. Substituting Equations (25) and (26) into Equations (17) and (18) and comparing the coefficients of (*d*31)0, (*d*33)<sup>0</sup> and (*d*15)0, we may obtain the zero-order perturbation equations

$$\begin{cases} \begin{aligned} \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_0^0}{\partial Z^2} + \frac{\partial^2 \Phi\_0^0}{\partial X^2} &= 0 \\ \frac{\partial^4 \overline{\Pi}\_0^0}{\partial Z^4} + (2S\_{13} + S\_{44}) \frac{\partial^4 \overline{\Pi}\_0^0}{\partial X^2 \partial Z^2} + S\_{33} \frac{\partial^4 \overline{\Pi}\_0^0}{\partial X^4} &= 0 \end{aligned} \end{cases} \tag{27}$$

The corresponding boundary conditions are

$$\int\_{-1/2}^{1/2} (-\frac{\partial^2 \overline{\mathcal{U}\_0}}{\partial Z \partial X}) dZ = -\frac{\overline{\mathcal{P}}}{\overline{\mathfrak{b}}}, \int\_{-1/2}^{1/2} \frac{\partial^2 \overline{\mathcal{U}\_0}}{\partial Z^2} dZ = 0,\\ \int\_{-1/2}^{1/2} Z \frac{\partial^2 \overline{\mathcal{U}\_0}}{\partial Z^2} dZ = \frac{\overline{\mathcal{M}}}{\overline{\mathfrak{b}}}, \text{ at } X = 0,\tag{28}$$

$$\frac{\partial^2 \overline{U}\_0^0}{\partial X^2} = -\frac{\partial^2 \overline{U}\_0^0}{\partial X \partial Z} = 0, \text{ at } Z = 1/2,\tag{29}$$

$$\frac{\partial^2 \overline{\mathcal{U}}\_0^0}{\partial X^2} = -\overline{\eta}\_\prime - \frac{\partial^2 \overline{\mathcal{U}}\_0^0}{\partial X \partial Z} = 0,\text{ at}\, Z = -1/2,\tag{30}$$

$$\int\_{-1/2}^{1/2} (-\frac{\partial \Phi\_0^0}{\partial X}) dZ = 0, \text{ at } X = 0 \text{ and } X = l/h \tag{31}$$

and

$$\frac{\partial \Phi\_0^0}{\partial Z} = 0, \text{ at } Z = \pm 1/2. \tag{32}$$

Suppose,

$$\begin{cases} \Phi\_0^0 = X^2 g\_1^0(Z) + Xg\_2^0(Z) + g\_3^0(Z) \\ \overline{L}\_0^0 = \frac{X^2}{2} f\_1^0(Z) + Xf\_2^0(Z) + f\_3^0(Z) \end{cases} \tag{33}$$

where *g*0*i* (*Z*) and *f* 0*i* (*Z*) (*i* = 1, 2, 3) are unknown functions of *Z* which can be determined by Equations (27) and (33), please see Appendix A.

Next, let us solve the first-order perturbation equations. Comparing the coefficients of (*d*31)1,(*d*33)<sup>1</sup> and (*d*15)1, we may obtain the first-order perturbation equations as follows.

For term (*d*31)1:

$$\begin{cases} \frac{\partial^3 \overline{\mathcal{U}}\_0^0}{\partial \mathcal{Z}^3} - \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_1^1}{\partial \mathcal{Z}^2} - \frac{\partial^2 \Phi\_1^1}{\partial \mathcal{X}^2} = 0\\ \frac{\partial^4 \overline{\mathcal{U}}\_1^l}{\partial \mathcal{Z}^4} + (2S\_{13} + S\_{44}) \frac{\partial^4 \overline{\mathcal{U}}\_1^l}{\partial \mathcal{X}^2 \partial \mathcal{Z}^2} + S\_{33} \frac{\partial^4 \overline{\mathcal{U}}\_1^l}{\partial \mathcal{X}^4} - \frac{\partial^3 \Phi\_0^0}{\partial \mathcal{Z}^3} = 0 \end{cases} \tag{34}$$

*Polymers* **2019**, *11*, 1934

for term (*d*33)1:

$$\begin{cases} \frac{\partial^3 \overline{\mathcal{U}}\_0^0}{\partial X^2 \partial \overline{\mathcal{U}}^2} = \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_2^1}{\partial \overline{\mathcal{U}}^2} + \frac{\partial^2 \Phi\_2^1}{\partial X^2} \\\ \frac{\partial^4 \overline{\mathcal{U}}\_2^1}{\partial \overline{\mathcal{U}}^4} + (2S\_{13} + S\_{44}) \frac{\partial^4 \overline{\mathcal{U}}\_2^1}{\partial X^2 \partial \overline{\mathcal{U}}^2} + S\_{33} \frac{\partial^4 \overline{\mathcal{U}}\_2^1}{\partial X^4} - \frac{\partial^3 \Phi\_0^0}{\partial X^2 \partial \overline{\mathcal{U}}^2} = 0 \end{cases} \tag{35}$$

and for term (*d*15)1:

$$\begin{cases} -\frac{\partial^3 \overline{U}\_0^0}{\partial X^2 \partial \overline{Z}} = \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_3^1}{\partial \overline{Z}^2} + \frac{\partial^2 \Phi\_3^1}{\partial X^2} \\\ \frac{\partial^4 \overline{U}\_3^1}{\partial \overline{Z}^4} + (2S\_{13} + S\_{44}) \frac{\partial^4 \overline{U}\_3^1}{\partial \overline{X^2 \partial \overline{Z}^2}} + S\_{33} \frac{\partial^4 \overline{U}\_3^1}{\partial X^4} + \frac{\partial^3 \Phi\_0^0}{\partial \overline{X^2 \partial \overline{Z}}} = 0 \end{cases} \tag{36}$$

The corresponding boundary conditions are

⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪⎩

$$\begin{cases} \int\_{-1/2}^{1/2} -\frac{\partial^2 \overline{U}\_1}{\partial Z \partial X} dZ = 0, \\ \int\_{-1/2}^{1/2} \frac{\partial^2 \overline{U}\_2}{\partial Z^2} dZ = 0, \quad \text{, at } X = 0, \\ \int\_{-1/2}^{1/2} Z \frac{\partial^2 \overline{U}\_3}{\partial Z^2} dZ = 0 \end{cases} \tag{37}$$

$$\frac{\partial^2 \overline{U\_i}^l}{\partial X^2} = -\frac{\partial^2 \overline{U\_i}^l}{\partial X \partial Z} = 0 (i = 1, 2, 3), \text{ at } Z = 1/2,\tag{38}$$

$$\frac{\partial^2 \overline{\mathcal{U}\_i^1}}{\partial X^2} = 0, -\frac{\partial^2 \overline{\mathcal{U}\_i^1}}{\partial X \partial Z} = 0 (i = 1, 2, 3), \text{ at } Z = -1/2,\tag{39}$$

$$\begin{cases} \int\_{-1/2}^{1/2} (-\frac{\partial \Phi^1\_1}{\partial X}) dZ = 0\\ \int\_{-1/2}^{1/2} (-\frac{\partial \Phi^1\_2}{\partial X}) dZ = 0\\ \int\_{-1/2}^{1/2} (-\frac{\partial^2 \overline{\Pi}^0\_0}{\partial Z dX} - \frac{\partial \Phi^1\_3}{\partial X}) dZ = 0 \end{cases}, \text{ at } X = 0 \text{ and } X = l/h,\tag{40}$$

and

$$\frac{\partial \Phi\_i^1}{\partial Z} = 0(i = 1, 2, 3), \text{ at } Z = \pm 1/2. \tag{41}$$

Similarly, suppose

$$\begin{cases} \Phi\_i^1 = X^2 g\_{3i-2}^1(Z) + X g\_{3i-1}^1(Z) + g\_{3i}^1(Z) \\ \overline{\mathcal{U}}\_i^1 = \frac{X^2}{2} f\_{3i-2}^1(Z) + X f\_{3i-1}^1(Z) + f\_{3i}^1(Z) \end{cases} \quad (i = 1, 2, 3), \tag{42}$$

where *<sup>g</sup>*<sup>I</sup>*i*(*Z*) and *f* I*i* (*Z*) (*i* = 1, 2, 3, ... , 9) are unknown functions of *Z* which can be determined by Equations (34)–(36) and (42), please see Appendix A.

Then, we solve the second-order perturbation equations. Comparing the coefficients of (*d*31)2, (*d*33)2, (*d*15)2, *d*31*d*33, *d*31*d*15 and *d*33*d*15, we may obtain the two-order perturbation equations as follows. For term (*d*31)2:

$$\begin{cases} \frac{\partial^3 \overline{\mathcal{U}}\_1^1}{\partial \mathcal{Z}^3} = \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_1^{\text{II}}}{\partial \mathcal{Z}^2} + \frac{\partial^2 \Phi\_1^{\text{II}}}{\partial \mathcal{X}^2} \\\ \frac{\partial^4 \overline{\mathcal{U}}\_1^{\text{II}}}{\partial \mathcal{Z}^4} + (2S\_{13} + S\_{44}) \frac{\partial^4 \overline{\mathcal{U}}\_1^{\text{II}}}{\partial \mathcal{X}^2 \partial \mathcal{Z}^2} + S\_{33} \frac{\partial^4 \overline{\mathcal{U}}\_1^{\text{II}}}{\partial \mathcal{X}^4} = \frac{\partial^3 \Phi\_1^{\text{I}}}{\partial \mathcal{Z}^3} \end{cases} \tag{43}$$

for term (*d*33)2:

$$\begin{array}{ll}\frac{\partial^{3}\overline{U}\_{2}^{\mathrm{I}}}{\partial X^{2}\partial Z} = \overline{\lambda}\_{33}\frac{\partial^{2}\Phi\_{2}^{\mathrm{II}}}{\partial Z^{2}} + \frac{\partial^{2}\Phi\_{2}^{\mathrm{II}}}{\partial X^{2}}\\\frac{\partial^{4}\overline{U}\_{2}^{\mathrm{II}}}{\partial Z^{4}} + \left(2S\_{13} + S\_{44}\right)\frac{\partial^{4}\overline{U}\_{2}^{\mathrm{II}}}{\partial X^{2}\partial Z^{2}} + S\_{33}\frac{\partial^{4}\overline{U}\_{2}^{\mathrm{II}}}{\partial X^{4}} = \frac{\partial^{3}\Phi\_{2}^{\mathrm{I}}}{\partial X^{2}\partial Z}\end{array} \tag{44}$$

*Polymers* **2019**, *11*, 1934

for term (*d*15)2:

$$\begin{cases} -\frac{\partial^3 \Pi\_3^\mathrm{fl}}{\partial X^2 \partial Z} = \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_3^\mathrm{fl}}{\partial Z^2} + \frac{\partial^2 \Phi\_3^\mathrm{fl}}{\partial X^2} \\ \frac{\partial^4 \Pi\_3^\mathrm{fl}}{\partial Z^4} + (2S\_{13} + S\_{44}) \frac{\partial^4 \Pi\_3^\mathrm{fl}}{\partial X^2 \partial Z^2} + S\_{33} \frac{\partial^4 \Pi\_3^\mathrm{fl}}{\partial X^4} = -\frac{\partial^3 \Phi\_3^\mathrm{fl}}{\partial X^2 \partial Z} \end{cases} \tag{45}$$

for term *d*31*d*33:

$$\begin{aligned} \frac{\partial^3 \Pi\_2^\mathrm{I}}{\partial Z^3} + \frac{\partial^3 \Pi\_1^\mathrm{I}}{\partial X^2 \partial Z} &= \overline{\Lambda}\_{33} \frac{\partial^2 \Phi\_4^\mathrm{II}}{\partial Z^2} + \frac{\partial^2 \Phi\_4^\mathrm{II}}{\partial X^2} \\ \frac{\partial^4 \Pi\_4^\mathrm{II}}{\partial Z^4} + \left(2S\_{13} + S\_{44}\right) \frac{\partial^4 \Pi\_4^\mathrm{I}}{\partial X^2 \partial Z^2} + S\_{33} \frac{\partial^4 \Pi\_4^\mathrm{I}}{\partial X^4} &= \frac{\partial^3 \Phi\_2^\mathrm{I}}{\partial Z^3} + \frac{\partial^3 \Phi\_1^\mathrm{I}}{\partial X^2 \partial Z} \end{aligned} \tag{46}$$

for term *d*31*d*15:

$$\begin{aligned} \frac{\partial^3 \overline{\Pi}\_3^\ddagger}{\partial Z^3} - \frac{\partial^3 \overline{\Pi}\_1^\ddagger}{\partial X^2 \partial Z} &= \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_5^\ddagger}{\partial Z^2} + \frac{\partial^2 \Phi\_5^\ddagger}{\partial X^2} \\ \frac{\partial^4 \overline{\Pi}\_5^\ddagger}{\partial Z^4} + \left( 2S\_{13} + S\_{44} \right) \frac{\partial^4 \Pi\_5^\ddagger}{\partial X^2 \partial Z^2} + S\_{33} \frac{\partial^4 \Pi\_5^\ddagger}{\partial X^4} &= \frac{\partial^3 \Phi\_3^\ddagger}{\partial Z^3} - \frac{\partial^3 \Phi\_1^\ddagger}{\partial X^2 \partial Z} \end{aligned} \tag{47}$$

and for term *d*33*d*15:

$$\begin{cases} \frac{\partial^3 \overline{\Pi}\_3^\mathbf{I}}{\partial X^2 \partial \overline{Z}} - \frac{\partial^3 \overline{\Pi}\_2^\mathbf{I}}{\partial X^2 \partial \overline{Z}} = \overline{\lambda}\_{33} \frac{\partial^2 \Phi\_6^\mathbf{I}}{\partial \overline{Z}^2} + \frac{\partial^2 \Phi\_6^\mathbf{II}}{\partial X^2} \\ \frac{\partial^4 \Pi\_6^\mathbf{I}}{\partial X^4} + \left(2S\_{13} + S\_{44}\right) \frac{\partial^4 \Pi\_6^\mathbf{I}}{\partial X^2 \partial \overline{Z}^2} + S\_{33} \frac{\partial^4 \Pi\_6^\mathbf{I}}{\partial X^4} = \frac{\partial^3 \Phi\_3^\mathbf{I}}{\partial X^2 \partial \overline{Z}} - \frac{\partial^3 \Phi\_2^\mathbf{I}}{\partial X^2 \partial \overline{Z}} \end{cases} \tag{48}$$

 The corresponding boundary conditions are

⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪⎩

⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪⎩

> ⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪⎩

$$\begin{cases} \int\_{-1/2}^{1/2} \frac{\partial^2 \overline{U}\_1^{\parallel}}{\partial Z \partial X} dZ = 0, \\ \int\_{-1/2}^{1/2} \frac{\partial^2 \overline{U}\_2^{\parallel}}{\partial Z^2} dZ = 0, \quad \text{, at } X = 0, \\ \int\_{-1/2}^{1/2} Z \frac{\partial^2 \overline{U}\_3^{\parallel}}{\partial Z^2} dZ = 0 \end{cases} \tag{49}$$

$$\frac{\partial^2 \overline{U}\_i^{\Pi}}{\partial X^2} = -\frac{\partial^2 \overline{U}\_i^{\Pi}}{\partial X \partial Z} = 0\\(i = 1, 2, 3, 4, 5, 6), \text{ at } Z = 1/2,\tag{50}$$

$$\frac{\partial^2 \overline{\mathcal{U}}\_i^{\Pi}}{\partial X^2} = 0,\\ -\frac{\partial^2 \overline{\mathcal{U}}\_i^{\Pi}}{\partial X \partial \mathcal{Z}} = 0\\ (i = 1, 2, 3, 4, 5, 6), \text{ at } \mathcal{Z} = -1/2,\tag{51}$$

$$\begin{cases} \int\_{-1/2}^{1/2} \frac{\partial \Phi^{\mathrm{I}}}{\partial \mathcal{X}} dZ = 0, \int\_{-1/2}^{1/2} \frac{\partial \Phi^{\mathrm{I}}\_{2}}{\partial \mathcal{X}} dZ = 0\\ \int\_{-1/2}^{1/2} \left(\frac{\partial^{2} \mathcal{U}\_{3}}{\partial \mathcal{X} \partial \mathcal{Z}} + \frac{\partial \Phi^{\mathrm{I}}\_{3}}{\partial \mathcal{X}}\right) dZ = 0, \int\_{-1/2}^{1/2} \frac{\partial \Phi^{\mathrm{I}}\_{4}}{\partial \mathcal{X}} dZ = 0\\ \int\_{-1/2}^{1/2} \left(\frac{\partial^{2} \mathcal{U}\_{1}}{\partial \mathcal{X} \partial \mathcal{Z}} + \frac{\partial \Phi^{\mathrm{I}}\_{5}}{\partial \mathcal{X}}\right) dZ = 0, \int\_{-1/2}^{1/2} \left(\frac{\partial^{2} \mathcal{U}\_{2}}{\partial \mathcal{X} \partial \mathcal{Z}} + \frac{\partial \Phi^{\mathrm{II}}\_{6}}{\partial \mathcal{X}}\right) dZ = 0 \end{cases}, \text{ at } X = 0 \text{ and } X = l/h, \tag{52}$$

and

$$\frac{\partial \Phi\_i^{\rm II}}{\partial Z} = 0, \text{ at } Z = \pm 1/2. \tag{53}$$

Suppose, 
$$\begin{cases} \Phi\_i^{\text{II}} = \mathbf{X}^2 g\_{3i-2}^{\text{II}}(Z) + \mathbf{X} g\_{3i-1}^{\text{II}}(Z) + g\_{3i}^{\text{II}}(Z) \\ \overline{\Pi}\_i^{\text{II}} = \frac{\mathbf{X}^2}{2} f\_{3i-2}^{\text{II}}(Z) + \mathbf{X} f\_{3i-1}^{\text{II}}(Z) + f\_{3i}^{\text{II}}(Z) \end{cases} (i = 1, 2, 3, 4, 5, 6), \tag{54}$$

where *g*II*i* (*Z*) and *f* II *i* (*Z*) (*i* = 1, 2, 3, ... , 18) are unknown functions of *Z* which can be determined by Equations (43)–(48) and (54), please see Appendix A.

Thus, we can obtain

Φ = *B*06 + *<sup>d</sup>*31[*X*<sup>2</sup>(− 3λ33 *qZ*<sup>2</sup> + 1 4λ33 *q*) + *<sup>X</sup>*(− 6 *b*λ33 *PZ*<sup>2</sup> + 1 2*b*λ33 *P*) + 1 <sup>2</sup>(<sup>λ</sup>33)<sup>2</sup> *qZ*<sup>4</sup> − 1 <sup>4</sup>(<sup>λ</sup>33)<sup>2</sup> *qZ*<sup>2</sup> +(2*S*13+*S*44) 2λ33 *qZ*<sup>4</sup> + 6 *b*λ33 *MZ*<sup>2</sup> − <sup>3</sup>(2*S*13+*S*44) 20λ33 *qZ*<sup>2</sup> + *B*I6] + *d*33[− 1 2λ33 *qZ*<sup>4</sup> + 3 4λ33 *qZ*<sup>2</sup> − 1 2λ33 *qZ* + *B*I12] + *<sup>d</sup>*15[−12 *qX*<sup>2</sup> − *Pb X* + 1 2λ33 *qZ*<sup>2</sup> + 1 2λ33 *qZ*<sup>4</sup> − 3 4λ33 *qZ*<sup>2</sup> + *B*I18] +(*d*31)<sup>2</sup>*B*II6 + (*d*33)<sup>2</sup>*B*II12 + (*d*15)<sup>2</sup>*B*II18 + *<sup>d</sup>*31*d*33*<sup>B</sup>*II24 + *<sup>d</sup>*31*d*15*<sup>B</sup>*II30 + *<sup>d</sup>*33*d*15*<sup>B</sup>*II36 (55)

*U* = *X*22 (−2*qZ*<sup>3</sup> + 32 *qZ* − *q*2 ) + *<sup>X</sup>*(−2*bPZ*<sup>3</sup> + 32*bPZ* + *C*08) + (2*S*13+*S*44) 10 *qZ*<sup>5</sup> + 2*bMZ*<sup>3</sup> − (2*S*13+*S*44) 20 *qZ*<sup>3</sup> + *C*011*Z* + *C*012 + *<sup>d</sup>*31(*XC*I8 + *C*I11*Z* + *C*I12) + *<sup>d</sup>*33(*XC*I20 + *C*I23*Z* + *C*I24) <sup>+</sup>*d*15(*XC*I32 + *C*I35*Z* + *<sup>C</sup>*I36)+(*d*31)<sup>2</sup>[*XC*II8 + 1 <sup>10</sup>(<sup>λ</sup>33)<sup>2</sup> *qZ*<sup>5</sup> + (2*S*13+*S*44) 10λ33 *qZ*<sup>5</sup> − 1 <sup>20</sup>(<sup>λ</sup>33)<sup>2</sup> *qZ*<sup>3</sup> − (2*S*13+*S*44) 20λ33 *qZ*<sup>3</sup> + *C*II11*Z* + *<sup>C</sup>*II12]+(*d*33)<sup>2</sup>(*XC*II20 + *C*II23*Z* + *<sup>C</sup>*II24)+(*d*15)<sup>2</sup>(*XC*II32 + *C*II35*Z* + *C*II36) <sup>+</sup>*d*31*d*33(*XC*II44 − 1 5λ33 *qZ*<sup>5</sup> + 1 10λ33 *qZ*<sup>3</sup> + *C*II47*Z* + *C*II48) + *<sup>d</sup>*31*d*15(*XC*II56 + 1 5λ33 *qZ*<sup>5</sup> − 1 10λ33 *qZ*<sup>3</sup> +*<sup>C</sup>*II59*<sup>Z</sup>* + *C*II60) + *<sup>d</sup>*33*d*15(*XC*II68 + *C*II71*Z* + *C*II72) . (56)

Finally, from Equations (55) and (56), we can obtain the expression of displacement components, stress components, and electric displacement components. The detailed derivation is shown in Appendix B. Thus, the bending problem of a piezoelectric cantilever beam under combined loads is solved. It can be seen from the derivation above that the piezoelectric effect is not shown in the zero-order perturbation solution, that is, the zero-order perturbation solution is the solution of the cantilever beam without piezoelectric properties which is regarded as the unperturbed system. The piezoelectric properties are only shown in the first-order and second-order perturbation solutions. In other words, the mechanical meaning of the first-order and second-order solutions is the influence of piezoelectric properties on the deformation of piezoelectric cantilever beams. This phenomenon is consistent with the basic idea of perturbation method.

#### **4. Comparison of the Solution Presented Here and the Existing Solution**

The theoretical solution for a piezoelectric cantilever beam under combined loads is given in this paper by a new method which is usually called the multi-parameter perturbation method. The validity of the theoretical solution should further be verified. For this purpose, we compare the solution presented here with the solution given in reference [17].

Before the comparison, we need to make a degradation of the solution presented here. In reference [17], only the concentrated force is considered. In this paper, however, the concentrated force, bending moment, and uniformly distributed load are all considered. Thus, for the convenience of comparison, we let the bending moment and uniformly distributed load equal to zero, that is, let

$$q = 0,\\
\mathcal{M} = 0.\tag{57}$$

Substituting Equation (57) into Equations (A44) and (A45), the displacement components can be transformed into

$$\begin{split} w &= \left(\frac{d\_{33}d\_{31}}{\lambda\_{33}} - s\_{13}\right) \frac{\delta P}{\delta \mathbf{h}^{3}} \mathbf{x} \mathbf{z}^{2} - \left(\frac{d\_{31}d\_{31}}{\lambda\_{33}} - s\_{11}\right) \frac{\Delta P}{\delta \mathbf{h}^{3}} \mathbf{x}^{3} \\ &+ \left(\frac{d\_{31}d\_{31}}{\lambda\_{33}} - s\_{11}\right) \frac{\delta P}{\delta \mathbf{h}^{3}} \mathbf{l}^{2} \mathbf{x} - \left(\frac{d\_{31}d\_{31}}{\lambda\_{33}} - s\_{11}\right) \frac{4P}{\delta \mathbf{h}^{3}} \mathbf{l}^{3} \end{split} \tag{58}$$

and

$$\begin{array}{l} \mu = \left(\frac{d\_{31}d\_{31}}{\lambda\_{33}} - s\_{11}\right)\frac{\delta P}{\delta h^{3}}x^{2}z - \left(\frac{d\_{33}d\_{31}}{\lambda\_{33}} - s\_{13} - s\_{44} - \frac{d\_{15}d\_{31}}{\lambda\_{33}}\right)\frac{2P}{\delta h^{3}}z^{3} \\\ -\left(\frac{d\_{31}d\_{31}}{\lambda\_{33}} - s\_{11}\right)\frac{\delta P}{\delta h^{3}}l^{2}z - \left(\frac{d\_{15}d\_{31}}{\lambda\_{33}} + s\_{44}\right)\frac{\delta P}{2\delta h}z \end{array} . \tag{59}$$

Similarly, substituting Equation (57) into Equations (A46), (A47), and (A48), the stress components can be written as

$$
\sigma\_x = -\frac{12P}{bh^3}xz\_\prime \tag{60}
$$

$$
\sigma\_z = 0,\tag{61}
$$

and

$$
\pi\_{zx} = \frac{\epsilon P}{bh^3} z^2 - \frac{3P}{2bh}.\tag{62}
$$

Substituting Equation (57) into Equations (A49) and (A50), the expressions of electric displacement components are

$$D\_x = (d\_{15} + \frac{\lambda\_{11} d\_{31}}{\lambda\_{33}}) (\frac{6P}{bh^3} z^2 - \frac{3P}{2bh}) \tag{63}$$

and

$$D\_z = 0.\tag{64}$$

By comparing Equations (58)–(64) with the expressions of displacement components, stress components, and electric displacement components in reference [17], it can be found that they are exactly the same, which indicates that the solution obtained here is correct. It should be mentioned that the structures studied in this paper and in reference [17] are both piezoelectric cantilever beams, but the structure in this paper is subjected to combined loads and the structure in reference [17] is subjected only to a concentrated force. In addition, non-dimensionalization is considered, the relevant experiments are carried out, and a new method called the multi-parameter perturbation method is given in this paper. These di fferences mentioned above constitute the advancements of this paper, compared with reference [17].

## **5. Experimental Verification**

To further verify the validity of the theoretical solution presented here, we carry out the relevant experiments of piezoelectric cantilever beams. The mechanical model of the theoretical part is shown in Figure 1, it can be seen that it is a piezoelectric cantilever beam subjected to three kinds of loads. In the experiment, it is very di fficult to apply these three kinds of loads at the same time. Therefore, we apply only the concentrated force at the cantilever end to carry out the experiments, that is, this experiment corresponds only to the case where the bending moment and the uniformly distributed load in the theoretical solution are zero. The details of the experiments are as follows. The main experimental equipments include a non-contact laser displacement sensor (ZSY Group Ltd, London, UK), a bench clamp (a cantilever beam clamping device), weights, and the ZLDS10X measuring software (ZSY Group Ltd, London, UK). The measuring range of the non-contact laser displacement sensor is 1 m, the accuracy is 0.01%, and the sampling frequency is 2 kHz. The experimental specimens consist of two groups of PbZrTiO3-5 (Generally abbreviated as PZT-5) piezoelectric ceramic sheets in which one group has piezoelectric properties and the other group has no piezoelectric properties. The size of the experimental specimens is 60 mm × 10 mm × 1 mm. The experimental specimen and non-contact laser displacement sensor are shown in Figure 2, the experimental device is shown in Figure 3, and the material constants are shown in Table 1.

The clamping length of the experimental specimens is 10 mm, therefore, the length of the piezoelectric cantilever beam is 50 mm. The deformations of the free end of piezoelectric cantilever beam are measured at the applied load 0.49 N, 0.98 N, and 1.96 N. The measured experimental data and theoretical calculation results are shown in Tables 2 and 3, respectively. It should be noted that the self-weight of the piezoelectric cantilever beam is 0.0367 N, and the ratio of the self-weight to the

minimum applied load is 0.075, which indicates that the self-weight of the piezoelectric cantilever beam is very small and thus may be ignored.

**Figure 2.** Scheme of experimental specimens and measuring instruments: (**a**) PbZrTiO3-5 (Generally abbreviated as PZT-5) piezoelectric ceramic specimens. (**b**) The non-contact laser displacement sensor.

(**a**) (**b**) 

**Figure 3.** Scheme of experimental device: (**a**) The cantilever beam device. (**b**) Theintegralmeasuring device.

**Table 1.** Physical properties of PZT-5 materials [33].




**Table 3.** Comparison of deformation test results between piezoelectric cantilever beam and cantilever beam without piezoelectric properties.

From Table 2, it can be seen that the theoretical results are in good agreemen<sup>t</sup> with the experimental results, and the relative errors under every level load are less than 15% allowed in engineering. This indicates that the analytical solution presented in this paper is reliable.

Table 3 shows that the deformation of the piezoelectric cantilever beam is smaller than the cantilever beam without piezoelectric properties. This means that the piezoelectric properties have a certain e ffect on the deformation of the piezoelectric cantilever beam, and its e ffect is, to a certain extent, hindering the deformation of the cantilever beam. This phenomenon can be explained by energy conservation. For piezoelectric cantilever beams, part of the work done by external forces is transformed into the elastic strain energy of piezoelectric cantilever beams, while the other part is transformed into the electric energy due to the existence of piezoelectric properties. For cantilever beams without piezoelectric properties, the work done by external forces is basically transformed into the elastic strain energy of cantilever beams. Therefore, the deformation of cantilever beams without piezoelectric properties is larger than that of cantilever beams with piezoelectric properties. The phenomenon mentioned above is commonly known as the piezoelectric sti ffening e ffect peculiar to piezoelectric materials and structures.
