# Investigation of Piezoelectric Energy Harvesting via Nonlinear Friction-Induced Vibration

Authors: D. W. Wang; M. X. Liu; X. Wu; W. J. Qian; Q. Ma; Z. Q. Wu
Journal: Shock and Vibration
 (2020)

---

## Abstract

In this work, piezoelectric energy harvesting (PEH) performance via friction-induced vibration (FIV) is studied numerically. A nonlinear two-degree-of-freedom friction system (mass-on-belt) with piezoelectric elements, which simultaneously considers the stick-slip motion, model coupling instability, separation, and reattachment between the mass and belt, is proposed. Both complex eigenvalue analyses and transient dynamic analysis of this nonlinear system are carried out. Results show that it is feasible to convert FIV energy to electrical energy when the friction system is operating in the unstable vibration region. There exists a critical friction coefficient (μc) for the system to generate FIV and output visible voltage. The friction coefficient plays a significant role in affecting the dynamics and PEH performance of the friction system. The friction system is able to generate stronger vibration and higher voltage in the case that both the kinetic friction coefficient and static friction coefficient are larger than μc. Moreover, it is seen that the separation behavior between contact pair can result in overestimating or underestimating the vibration magnitude and output voltage amplitude, and the overestimate or underestimate phenomenon is determined by the located range of friction coefficient. Furthermore, it is confirmed that an appropriate value of external resistance is beneficial for the friction system to achieve the highest output voltage. The obtained results will be beneficial for the design of PEH device by means of FIV.


---

## Body

## 1. Introduction

Recently, with the problem of energy shortage becoming more and more serious, energy harvesting from ambient sources has become an important research field and received lots of attention from academia and industry [1–3]. Among the diverse new and renewable ambient energy sources, vibration, as a typical and common mechanical motion, has become an attractive source by virtue of its ubiquity and ease of access through MEMS technology. Over last decades, energy harvesting from ambient vibration has been extensively studied [4–7].Generally, five main techniques can realize the conversion of vibration energy to electric energy, i.e., electromagnetic, electrostatic, piezoelectric, magnetostrictive, and triboelectric [8–14]. All these techniques have been successfully demonstrated experimentally and numerically in the past. An overview of the basics in energy harvesting has been provided by Priya et al. [3]. Among these techniques, piezoelectric energy harvesting (PEH) has attracted much more attention as it is easy to be implemented and does not rely on external voltage or magnetic field [15–20]. However, widespread usage of PEH technique is likely to be shadowed by a certain limitation; this limitation is mainly due to the fact that the natural frequency of PEH device is often not consistent with the vibration frequency from the ambient. As a matter of fact, an energy harvester can effectively implement its own functions only if ambient vibration frequency well matches with the resonant frequency of the harvester, which indicates that, with a slight deviation between these two frequencies, the ability of harvesters in energy harvesting would be reduced rapidly [6, 21–23]. Considering this, many new approaches, such as structure modification, circuit improvement, and nonlinear techniques, have been proposed to broaden the frequency response spectrum of harvester to match the ambient vibration frequency as much as possible [24–35]. However, most ambient vibrations are random and their frequency spectrums drift over time, which consequently brings about great challenges for the design of PEH devices [36].To solve this problem, a new approach is proposed: seeking a kind of vibration energy which is independent of the ambient frequency and can be used for harvesting and conversion, such as friction-induced vibration (FIV). FIV is a typical self-excited vibration, which does not require any exterior force excitation, but it can be triggered under specific working conditions [37]. FIV phenomenon is commonly seen in our daily lives [38, 39]. Examples include chatter or squeak from brake systems, tool cutting systems, mechanical gear systems, window wiper blades, and lead screw drives [40–42]. Considering that the FIV is uncontrollable and unpredictable due to its inherently irregular nature and it is undesirable and harmful in most cases; thus a large amount of research on FIV mainly focuses on understanding its mechanism and searching for potential effective approaches to improve FIV of friction systems [43–47]. Meanwhile, from another perspective, sometimes the FIV can be useful if this kind of vibration energy is collected in a certain way and converted to electrical energy, which can be thought of as “makes the best out of an undesirable situation” [48]. Furthermore, FIV occurs simultaneously in both the normal and tangential directions of the friction system, indicating that the vibration energy in both directions can be collected [8].Up to now, many works have pointed out that the most probable reasons for the occurrence of FIV are stick-slip motion and model coupling phenomenon [49–51]. To study the energy harvesting performance related to the stick-slip motion, Helseth [52] established a hook and loop system and found that the stick-slip motion between hook and loop would cause the fluctuations of contact force. These force fluctuations could accordingly excite different energy harvesting systems and generate fluctuating electric signals. Tadokoro et al. [53] and Masuda and Sawai [54] experimentally verified the feasibility of PEH via stick-slip motion; moreover, Tadokoro et al. [53] also established a one-degree-of-freedom (DOF) mathematical model to simulate the PEH by means of stick-slip motion. On the other hand, to study the energy harvesting related to the model coupling instability, Wang et al. [8] experimentally and numerically studied the mode-coupling instability between a pad and disc and verified the feasibility of energy harvesting from high-amplitude FIV. By using a wobbling disc energy harvester, Heffel and Hagedorn [55] found that converting FIV energy to electrical energy could lead to a broadband applicability. In addition, in a recent work performed by Stender et al. [48], external forcing was verified to be able to pull part of the friction energy out of the system and thus make energy harvesting possible in certain parameter ranges. These studies are definitely necessary for realizing energy harvesting via FIV. However, the knowledge on this aspect is still limited, and the PEH behavior via FIV remains unknown and somehow mysterious. In our previous work to detect the PEH via FIV, a single-degree-of-freedom friction system is proposed to discuss the possible PEH via stick-slip vibration. Secondly, a two-degree-of-freedom friction system is used to discuss the possible PEH via model-coupling vibration. Visibly, these two parts are almost independent [56]. However, in a real world of FIV, the stick-slip motion and model coupling instability usually occur simultaneously, and sometimes separation and subsequent reattachment behaviors between the contact pair may appear as well during the friction process [57, 58], which indicates that a more specific and complicated mathematical model should be established, to reveal the relationship between energy harvesting and FIV numerically. This paper makes such a contribution.In this work, the PEH performance by means of FIV is analyzed numerically. A nonlinear two-degree-of-freedom (DOFs) friction system (mass-on-belt) with piezoelectric elements, which simultaneously considers the stick-slip motion, model coupling, separation, and reattachment between the mass and belt, is proposed. Through introducing piezoelectric elements into the friction system, the energy harvesting behaviors along both the tangential and normal directions by means of FIV are verified. Moreover, the effects of friction coefficient and the separation behavior between contact pair on the dynamics and energy harvesting are studied. Finally, the dynamics behaviors of the friction system with different external electric resistances are analyzed, and the role of electric resistances in affecting PEH performance is revealed. Some conclusions on the piezoelectric energy harvesting related to FIV are presented.


## 2. Theoretical Modelling

A friction system with two DOFs is employed to study the PEH via FIV; this system is extended from the work of Tadokoro et al. [53]. In their work, a ball sample rubs against a sliding flat sample and stick-slip motion occurs; then a piezoelectric film is attached with the fixture of the test setup and successfully transfers the FIV energy to electric energy. Based on the test configuration, a one-degree-of-freedom friction system is proposed to investigate the PEH via stick-slip friction motion. Inspired by this system, a more realistic two-DOF friction system which simultaneously accounts for the stick, slip, loss of contact, and reattachment between the contact pair due to mode-coupling is proposed in this work, as shown in Figure 1.Figure 1
A 2-DOF friction system.The slider (m) rubs against a moving belt, which has a constant siding velocity v0. A number of linear springs and dampers link with the slider, which are used to store and dissipate energy of the friction system. The spring k1 and the damper c1 connect with the slider along the tangential direction (x direction), while a damper c2 connects with the slider along the normal direction (y direction). Additionally, piezoelectric elements named Piezo-1 and Piezo-2 (loaded with external electric resistance R) are, respectively, introduced into the friction system along the tangential and normal directions. These two piezo elements can respond to the vibration of the system through strain, so as to convert the vibration energy in both directions. A spring k3 is linked to the slider at a 135° angle relative to the tangential direction, which provides coupling in the tangential and normal directions. Moreover, spring k2 and knl, which are characterized by linear and cubic spring, respectively, are used to simulate the contact stiffness between contact pair. The effects of impact at reattachment between the slider and belt are neglected [57]. A constant normal compressive force, F, directly acts on the slider along the y direction.If the slider is still in contact with the belt, the contact forceFN between the contact pair is expressed as(1)FN=−k2y−knly3.The equations of friction motion for the slider yield(2)m0000m0000000000x¨y¨V¨xV¨y+c10000c200α10−Cp100α20−Cp2x˙y˙V˙xV˙y+k1+k32−k32+μk2α10−k32k2+k320α200−1R10000−1R2xyVxVy+0μknl000knl0000000000x3y3Vx3Vy3=0−F00.In the above equation, the parametersα1 and α2are defined as the force factors of Piezo-1 and Piezo-2, respectively; Vxand Vy represent the generated voltages from Piezo-1 and Piezo-2, respectively. It should be emphasized that equation (2) is valid only if the contact force still exists between contact pairs; that is,(3)FNt>0.IfFNt is found to be smaller than 0, separation is identified to happen between the contact pair. The normal displacement of slider is larger than zero in this moment, and the springs k2 and knl will not work on the contact stiffness. As a consequence, the situation where the slider separates with the belt can be depicted as shown in Figure 2, and the dynamical equation of the system can be expressed as the following equation:Figure 2
The situation where slider separates with the sliding belt.(4)m0000m0000000000x¨y¨V¨xV¨y+c10000c200α10−Cp100α20−Cp2x˙y˙V˙xV˙y+k1+k32−k32α10−k32k320α200−1R10000−1R2xyVxVy=0−F00.Assuming that the belt surface is rigid, the condition of separation between the contact pair is determined by the normal displacement of slider, which is written as(5)yt>0.When the normal displacementyt becomes less than zero, it is identified that the slider reattaches with the belt again. Then the dynamical equation of the friction is determined by equation (2) until the next round of separation appears again. Therefore, the velocity, displacement, and other state parameter values determined by equation (2) at the last step in contact are used for the initial conditions for equation (4), and the states of the mass before reattachment are the last moment in the separation motion (equation (2)).


## 3. Eigenvalue Analysis

It is meaningful to study the stability of the slider in the vicinity of the origin, which will be beneficial for detecting the critical friction coefficientμc for the generation of unstable vibration. In the eigenvalue analysis procedure, the slider is assumed to stay in the steady sliding stage; thus the separation and reattachment behaviors between the contact pair will not be accounted, and only (2) is employed to calculate the eigenvalue of the friction system.The equilibrium points (EPs) (xe, ye, Vxe, and Vye) of the friction system are found by solving (2) in the case where the first and second derivatives of these four variables are assumed to be zero, which become(6)k1+k32−k32+μk2α10−k32k2+k320α200−1R10000−1R2xyVxVy+0μknl000knl0000000000x3y3Vx3Vy3=0−F00.Define the transformation relationshipsz1=x,z2=x˙,z3=y,z4=y˙,z5=Vx,z6=Vy. Thus equation (6) can be expressed in the following form:(7)z˙1z˙2z˙3z˙4z˙5z˙6=010000−k1+k3/2M−c1M−μk2−k3/2M0−α1M0000100k3/2M0−k2+k3/2M−c2M0−α2M0α1Cp100−1R1Cp10000α2Cp20−1R2Cp2z1z2z3z4z5z6+00000000−μknlM00000000000−knlM000000000000000z13z23z33z43z53z63+000−FM00.The Jacobian matrix can be derived via expanding the nonlinear equation (7) in a first-order truncated Taylor series at the equilibrium point, as written as follows:(8)J=010000−k1+k3/2M−c1M−μk2−k3/2M−3μknlMye20−α1M0000100k3/2M0−k2+k3/2M−3knlMye2−c2M0−α2M0α1Cp100−1R1Cp10000α2Cp20−1R2Cp2.The stability of the system can be evaluated by calculating the eigenvaluesλ of the Jacobian matrix [59]. The eigenvalues λ of the matrix J can be obtained by solving the following equation:(9)detλE−J=0.The parameters set in this analysis are from the work performed by Tadokoro et al. [53], which are listed in Table 1. Figure 3 depicts the variation of real and imaginary parts of the eigenvalues with the increase of friction coefficients. Visibly, there exists a specified friction coefficient value (0.85) at which system will become unstable. This value of friction coefficient is called critical friction coefficient μc. The eigenvalues are a pair of complex conjugates with identical negative real parts for friction coefficient less than the critical value. At and beyond this value, the real parts of two eigenvalues branch out to one positive and one negative, while the corresponding imaginary parts start to merge into one value. Such mode merging phenomena indicate that the friction system has a strong tendency to generate self-excited instability, and accordingly the vibration energy may be transferred to electric energy through the action of Piezo elements. Based on the critical friction coefficient μc, the nonlinear transient dynamic analysis of the friction system will be performed in the next section.Table 1
The parameters of the friction system set in the eigenvalues analysis.
Mass (slider)0.08 kgk3 (coupling spring in two directions)600 N/mk1 (including the stiffness of Piezo-1)890 N/mc1 = c2 (including the damping values of piezo elements)0.077 N/(m/s)knl (nonlinear contact stiffness)1200 N/m3α1=α2 (force factor)0.00036 N/Vk2 (linear contact stiffness)500 N/mCp1=Cp2 (clamped capacitance of the piezo elements)1.79 × 10−7 Fv0 (belt velocity)5 mm/s∼200 mm/sF (normal load)20 NR110 kΩR210 kΩFigure 3
The real parts (a) and imaginary parts (b) of the friction system with the increase of friction coefficient.
(a)(b)


## 4. Transient Dynamic Analysis

### 4.1. Procedure of the Transient Dynamic Analysis

In the transient dynamic analysis, the stick, slip, separation, and reattachment behaviors between the slider and belt are taken into consideration. Assume that the relative sliding velocityvr is expressed as(10)vr=v0−x˙t.Therefore, the friction system has two motion phases when the slider is sliding on the belt (y < 0), that is, slip phase vr≠0 and stick phase vr=0. For convenience of numerical computation, the criterion conditions for the relative velocity vr≠0 and vr=0 are changed to vr>ε and vr≤ε, respectively, where ε represents a threshold value for the sliding velocity that is sufficiently close to zero (10−6).Whenvr>ε, that is, the friction system stays in the slip stage, the friction force Ffis expressed as(11)Ff=−μkk2yt+knlyt3sgnvr.Here the parameterμk represents the kinetic friction coefficient of the friction system. When vr≤ε, it indicates that the friction system is in the stick stage; the friction force Ff is determined by the maximum static friction force Fstatic and the applied horizontal force Fapplied on the slider. Here Fstatic is expressed as(12)Fstatic=−μsk2yt+knlyt3sgnvr.The applied horizontal force is expressed as(13)Fapplied=k1xt+c1x˙t+αV1t+12k3xt−12k3yt,and therefore, in the stick phase, the friction forceFf can be summarized in the following form:(14)Ff=k1xt+c1x˙t+α1V1t+12k3xt−12k3yt,whenFapplied<Fstatic,−μsk2yt+knly3tsgnvr,whenFapplied>Fstatic.Then equation (2) can be transformed to the following equation:(15)z˙1z˙2z˙3z˙4z˙5z˙6=010000−k1+k3/2M−c1Mk3/2M0−α1M0000100k3/2M0−k2+k3/2M−c2M0−α2M0α1Cp100−1R1Cp10000α2Cp20−1R2Cp2z1z2z3z4z5z6+00000000000000000000−knlM000000000000000z13z23z33z43z53z63+0FfM0−FM00.During sliding stage, the dynamical equation of the system is determined by equation (15), and the contact force between contact pair is estimated at the last time step; once the contact force FN<0, it is suggested that the slider is separating with belt. At this moment, the equation of motion is changed to equation (4). The critical point is determined, which corresponds to the moment when the contact force FN=0. At this critical point, the vibration state is switched from sliding stage to separation stage. During the separation process, the normal displacement of the slider yt is detected, as long as yt<0, indicating that the slider is touching with the belt. The critical point yt=0 suggests the generation of reattachment between contact pair; at this moment, the equation of motion is changed to equation (2) until the stage of separation is satisfied again.Equation (4) is used to calculate the dynamic behavior of the friction system in the separation process, which can be expressed as follows:(16)z˙1z˙2z˙3z˙4z˙5z˙6=010000−k1+k3/2M−c1Mk3/2M0−α1M0000100k3/2M0−k2+k3/2M−c2M0−α2M0α1Cp100−1R1Cp10000α2Cp20−1R2Cp2z1z2z3z4z5z6+000−FM00.Visibly, it is a typical nonlinear, nonsmooth, and discontinuous friction system, which brings some difficulties in determining the moment that the motion of the slider transits from one stage to another stage. The flow chart of the numerical computation is shown in Figure4.Figure 4
The flow chart of the transient dynamic analysis.


### 4.2. Transient Dynamic Analysis Results

The critical friction coefficient value for the generation of self-excited vibration is determined in the above section, which is equal to 0.85. In this section, friction coefficient is set as 0.6, 0.85, and 1.2, respectively, and used in the transient dynamic analysis. These friction coefficients are lower, equal to, and larger than the critical friction coefficient, respectively. The belt velocity (v0) is set as 0.1 m/s. The corresponding dynamical behaviors and output voltage are shown in Figure 5. No visible sustained vibration signals can be seen for the friction system at a lower friction coefficient; hence no visible output voltages can be detected. With the friction coefficient increasing to the critical value, strong self-excited vibration appears and consequently visible voltage output is found to be generated in both directions. This phenomenon indicates that the FIV energy converts to electric energy when the FIV appears. In the case of friction coefficient equal to 1.2, the system exhibits larger amplitude of vibration. As a consequence, higher voltages are outputted in both directions. In addition, it can be seen that the vibration in the normal direction is larger than that in the tangential direction, which accordingly results in a higher output voltage in the normal direction. This is because the normal contact stiffness will not work on the system during the separation stage, which results in the larger displacement in the normal direction, while the friction force disappears and will not act on the slider. Thus, the transient dynamic analysis results further prove the results exhibited in the stability analysis and confirm the flexibility of the PEH by means of FIV.Figure 5
The tangential displacement (a), normal displacement (b), tangential output voltage (c), and normal output voltage (d) of the friction system with different friction coefficients.
(a)(b)(c)(d)Figure6 shows the displacement versus velocity trajectories of friction system with different coefficients in tangential and normal directions, respectively. As can be seen from the periodic orbits appearing in this figure, the steady-state vibration amplitude of the friction system with larger friction coefficient is much larger than that of the friction system with lower friction coefficient, which further indicates that larger friction coefficient leads to a stronger vibration level. Additionally, due to the consideration of separation and reattachment behavior between the slider and belt, the slider’s tangential stick-slip motion becomes unconventional: sometimes the speed of the slider can be even larger than the belt speed during the friction process and forms an interesting stick-slip limit cycle.Figure 6
The tangential displacement versus tangential velocity. (a) Normal displacement versus normal velocity. (b) Phase portraits of the friction system.
(a)(b)Expectedly, the large-amplitude periodic (limit-cycle) response on the high-energy orbit of the friction system with larger coefficient is also observed in the velocity versus voltage trajectory shown in Figure7. For a larger friction coefficient, the orbit exhibits a large-amplitude limit cycle, which indicates that a sustained and stable voltage can be outputted during the friction process, while for the friction system with lower friction coefficient, the orbit gradually evolves to a point with the increase of time, suggesting that no visible FIV energy is transferred to electric energy. This phenomenon further verifies the role of friction coefficient in affecting the stability and energy harvesting behaviors of friction system and also supports the opinion that the FIV can be transferred to electric energy by piezo elements.Figure 7
The open-circuit voltage versus velocity phase trajectories.
(a)(b)


## 4.1. Procedure of the Transient Dynamic Analysis

In the transient dynamic analysis, the stick, slip, separation, and reattachment behaviors between the slider and belt are taken into consideration. Assume that the relative sliding velocityvr is expressed as(10)vr=v0−x˙t.Therefore, the friction system has two motion phases when the slider is sliding on the belt (y < 0), that is, slip phase vr≠0 and stick phase vr=0. For convenience of numerical computation, the criterion conditions for the relative velocity vr≠0 and vr=0 are changed to vr>ε and vr≤ε, respectively, where ε represents a threshold value for the sliding velocity that is sufficiently close to zero (10−6).Whenvr>ε, that is, the friction system stays in the slip stage, the friction force Ffis expressed as(11)Ff=−μkk2yt+knlyt3sgnvr.Here the parameterμk represents the kinetic friction coefficient of the friction system. When vr≤ε, it indicates that the friction system is in the stick stage; the friction force Ff is determined by the maximum static friction force Fstatic and the applied horizontal force Fapplied on the slider. Here Fstatic is expressed as(12)Fstatic=−μsk2yt+knlyt3sgnvr.The applied horizontal force is expressed as(13)Fapplied=k1xt+c1x˙t+αV1t+12k3xt−12k3yt,and therefore, in the stick phase, the friction forceFf can be summarized in the following form:(14)Ff=k1xt+c1x˙t+α1V1t+12k3xt−12k3yt,whenFapplied<Fstatic,−μsk2yt+knly3tsgnvr,whenFapplied>Fstatic.Then equation (2) can be transformed to the following equation:(15)z˙1z˙2z˙3z˙4z˙5z˙6=010000−k1+k3/2M−c1Mk3/2M0−α1M0000100k3/2M0−k2+k3/2M−c2M0−α2M0α1Cp100−1R1Cp10000α2Cp20−1R2Cp2z1z2z3z4z5z6+00000000000000000000−knlM000000000000000z13z23z33z43z53z63+0FfM0−FM00.During sliding stage, the dynamical equation of the system is determined by equation (15), and the contact force between contact pair is estimated at the last time step; once the contact force FN<0, it is suggested that the slider is separating with belt. At this moment, the equation of motion is changed to equation (4). The critical point is determined, which corresponds to the moment when the contact force FN=0. At this critical point, the vibration state is switched from sliding stage to separation stage. During the separation process, the normal displacement of the slider yt is detected, as long as yt<0, indicating that the slider is touching with the belt. The critical point yt=0 suggests the generation of reattachment between contact pair; at this moment, the equation of motion is changed to equation (2) until the stage of separation is satisfied again.Equation (4) is used to calculate the dynamic behavior of the friction system in the separation process, which can be expressed as follows:(16)z˙1z˙2z˙3z˙4z˙5z˙6=010000−k1+k3/2M−c1Mk3/2M0−α1M0000100k3/2M0−k2+k3/2M−c2M0−α2M0α1Cp100−1R1Cp10000α2Cp20−1R2Cp2z1z2z3z4z5z6+000−FM00.Visibly, it is a typical nonlinear, nonsmooth, and discontinuous friction system, which brings some difficulties in determining the moment that the motion of the slider transits from one stage to another stage. The flow chart of the numerical computation is shown in Figure4.Figure 4
The flow chart of the transient dynamic analysis.


## 4.2. Transient Dynamic Analysis Results

The critical friction coefficient value for the generation of self-excited vibration is determined in the above section, which is equal to 0.85. In this section, friction coefficient is set as 0.6, 0.85, and 1.2, respectively, and used in the transient dynamic analysis. These friction coefficients are lower, equal to, and larger than the critical friction coefficient, respectively. The belt velocity (v0) is set as 0.1 m/s. The corresponding dynamical behaviors and output voltage are shown in Figure 5. No visible sustained vibration signals can be seen for the friction system at a lower friction coefficient; hence no visible output voltages can be detected. With the friction coefficient increasing to the critical value, strong self-excited vibration appears and consequently visible voltage output is found to be generated in both directions. This phenomenon indicates that the FIV energy converts to electric energy when the FIV appears. In the case of friction coefficient equal to 1.2, the system exhibits larger amplitude of vibration. As a consequence, higher voltages are outputted in both directions. In addition, it can be seen that the vibration in the normal direction is larger than that in the tangential direction, which accordingly results in a higher output voltage in the normal direction. This is because the normal contact stiffness will not work on the system during the separation stage, which results in the larger displacement in the normal direction, while the friction force disappears and will not act on the slider. Thus, the transient dynamic analysis results further prove the results exhibited in the stability analysis and confirm the flexibility of the PEH by means of FIV.Figure 5
The tangential displacement (a), normal displacement (b), tangential output voltage (c), and normal output voltage (d) of the friction system with different friction coefficients.
(a)(b)(c)(d)Figure6 shows the displacement versus velocity trajectories of friction system with different coefficients in tangential and normal directions, respectively. As can be seen from the periodic orbits appearing in this figure, the steady-state vibration amplitude of the friction system with larger friction coefficient is much larger than that of the friction system with lower friction coefficient, which further indicates that larger friction coefficient leads to a stronger vibration level. Additionally, due to the consideration of separation and reattachment behavior between the slider and belt, the slider’s tangential stick-slip motion becomes unconventional: sometimes the speed of the slider can be even larger than the belt speed during the friction process and forms an interesting stick-slip limit cycle.Figure 6
The tangential displacement versus tangential velocity. (a) Normal displacement versus normal velocity. (b) Phase portraits of the friction system.
(a)(b)Expectedly, the large-amplitude periodic (limit-cycle) response on the high-energy orbit of the friction system with larger coefficient is also observed in the velocity versus voltage trajectory shown in Figure7. For a larger friction coefficient, the orbit exhibits a large-amplitude limit cycle, which indicates that a sustained and stable voltage can be outputted during the friction process, while for the friction system with lower friction coefficient, the orbit gradually evolves to a point with the increase of time, suggesting that no visible FIV energy is transferred to electric energy. This phenomenon further verifies the role of friction coefficient in affecting the stability and energy harvesting behaviors of friction system and also supports the opinion that the FIV can be transferred to electric energy by piezo elements.Figure 7
The open-circuit voltage versus velocity phase trajectories.
(a)(b)


## 5. Impact Factors Analysis

It has been reported that the dynamics of friction system are very sensitive to the contact parameters [38–40], which indicates that the variation of working parameters can significantly change the energy harvesting behaviors of friction system as well. In this section, the effect of friction coefficient on the dynamics and energy harvesting is studied. Then, the dynamics and energy harvesting behaviors of the friction system with and without consideration of the separation between slider and belt are analyzed, respectively. Finally, the dynamics behaviors of the friction system with different external electric resistances are analyzed, and the role of electric resistances in affecting PEH performance is revealed.


### 5.1. Friction Coefficient

According to the value of critical friction coefficient (0.85), the dynamic behaviors of the friction system with kinetic and static friction coefficients, denoted byμk and μs, located within three ranges (μk<μs<μc,μk<μc<μs, and μc<μk<μs), are examined and compared. The values of friction coefficients in these three cases are listed in Table 2, and concurrently listed in Table 2 are the Eps of the friction system under different friction coefficients. The initial condition is set at (xe, 0.1, ye, 0, 0, 0), with the tangential velocity of the slider equaling the velocity of the belt. The external electric resistance in this analysis is set at 10 KΩ.Table 2
The values of friction coefficients and the corresponding Eps in these three cases.
CaseμkμsμcEps (xe, ye)10.40.60.85(−0.00216, −0.0258)20.70.90.85(0.00104, −0.0246)311.20.85(0.00397, −0.0235)Figure8 depicts the dynamical behaviors and output voltages of the friction system with three different cases of friction coefficients. Visibly, the vibration level and output voltage are strongly related to the friction coefficient. In the tangential direction, the friction system exhibits the strongest self-excited vibration level when the friction coefficient is located in Case 3. Accordingly, the voltage generated in this situation exhibits the highest value. Meanwhile, for the lowest friction coefficient of Case 1, the tangential vibration amplitude and output voltage are visibly lower than those generated from the other two cases (Figures 8(a) and 8(c)). Similarly, for the normal direction, a larger coefficient is found to be able to generate stronger vibration and higher output voltage as well. Therefore, the vibration level and voltage values generated in Case 3 are larger than those from Case 2, and vibration and voltage signals generated from Case 1 show the lowest values.Figure 8
The tangential displacement (a), normal displacement (b), tangential output voltage (c), and normal output voltage (d) of the friction system with three different cases of friction coefficients.
(a)(b)(c)(d)Figure9 shows the phase plots of friction system with different friction coefficient cases. For Case 1, the tangential stick-slip limit cycle is very complicated and exhibits lots of different cycles of stick-slip motion. Similarly, the normal vibration also shows a multiperiod limit cycle. Meanwhile, for Case 2, the phase trajectory in tangential direction is very unconventional, in which several limit cycles are characterized by appearing and repeating themselves; additionally, the tangential velocity of slider exceeds the belt velocity (v > 0.1 m/s) and causes overshoot phenomenon during the vibration. Additionally, the normal vibration also forms a multiperiod limit cycle. Case 3 represents the situation in which both the static and kinetic friction coefficients are larger than the critical friction coefficient; it is seen that the tangential stick-slip and overshoot always happen. In this case, the friction system has a large-amplitude (limit-cycle) response in both the normal direction and the tangential direction, compared with the other two cases.Figure 9
The displacement verus velocity trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d).
(a)(b)(c)(d)The velocity versus voltage trajectory of the friction system in three cases of friction coefficient is plotted. In the tangential direction, it is visible that the orbit of friction system with Case 3 possesses biggest limit cycle, which suggests that the friction system can generate higher amplitude voltage than the other two cases. Meanwhile the output voltage amplitude in Case 1 shows the smallest limit cycle. Indeed, both trajectories from Case 1 and Case 2 exhibit very complicated orbits. In the normal direction, as expected, the orbit of friction system with Case 3 exhibits the largest amplitude of limit cycle, which indicates that a sustained and high voltage is outputted in this direction during the sliding process. Additionally, the orbit of friction system with Case 1 forms a smallest limit cycle, which reflects that relative lower voltage is generated in this situation. This phenomenon can well reflect the time domain signal of voltage shown in Figure8.To sum up, the dynamic and output voltages are strongly related to the relationship among the kinetic friction coefficientμk, static friction coefficient μs, and critical friction coefficient μc. The friction system can generate stronger vibration and output higher voltage when the friction coefficient is located in the range of μc<μk<μs. In contrast, lower vibration and less voltage are generated when the friction coefficient locates in the range of μk<μs<μc.


### 5.2. Comparison and Analysis of Ignoring/considering Separation

In most previous studies of FIV of small models, it is assumed that the slider is still in contact with belt during vibration process [50, 60]. Although the possibility of loss of contact was mentioned in [61, 62], the effect of separation was not studied. As the unstable FIV caused by the friction force at the contact interface grows, it is important to consider separation during the vibration, when friction temporarily disappears. The loss of contact will affect not only the contact behaviors of friction system but also the energy harvesting behaviors. Therefore, in this section, the dynamics and energy harvesting behaviors of the friction system in both cases, that is, ignoring separation and considering separation between slider and the belt, are analyzed and compared.The flow chart of transient dynamic analysis when the friction system ignores the separation between slider and belt is shown in Figure10. Due to ignoring the separation, the equations of motion are given by equation (2). The iteration process starts from the stick mode, and slider velocity is checked at the end of every single time step during numerical computation. Once the slider velocity is not equal to the belt velocity (0.1 m/s), the iteration process switches to the slide mode, and initial value of slide mode is the last value of the friction system in the stick mode. Then, the slider rubs on the moving belt until the next stick happens.Figure 10
The flow chart of transient dynamic analysis when the friction system ignores the separation between belt and slider.In the previous study, it is found that the separation can occur in Case 2 and Case 3; thus, it is worth noting that comparison between ignoring separation and considering separation will be performed in these two cases.


#### 5.2.1.μk<μc<μs

Numerical results withμk = 0.7, μs = 0.9, and initial condition (xe, 0.1, ye, 0, 0, 0) are shown in Figure 11. Visibly, the separation and reattachment behaviors make both the tangential and normal dynamics of friction system different from their counterparts when separation is ignored. For the tangential direction, the displacement mainly reflects in the positive direction when considering separation, while when separation is ignored, the displacement fluctuates between the positive and negative directions, suggesting that the slider moves back and forth with the moving belt. In addition, the amplitude of displacement fluctuation in the case of considering separation is larger than that in the case of ignoring separation, as shown in Figure 11(a); thus the output voltage in this direction is larger in the case of considering separation (Figure 11(c)). Meanwhile, for the normal direction, the situation changes significantly. It can be seen that the normal displacement amplitude when considering separation is smaller compared with the situation of ignoring separation; accordingly the output voltage is smaller when separation is taken into consideration, as shown in Figures 11(b) and 11(d).Figure 11
The tangential displacement (a), normal displacement (b), tangential voltage (c), normal voltage (d), and contact force (e) of the friction system in the time domain whenμk = 0.7 and μs = 0.9.
(a)(b)(c)(d)(e)The contact forces of these two situations are plotted in Figure11(e). For the case when separation is considered, multiple separation events can be observed as the values of contact force vary from zero to positive and then drop to zero during the vibration. Meanwhile, when ignoring separation, the value of contact force gradually grows and becomes negative, and the amplitude is larger than that in the case of considering separation. This phenomenon well reflects the results shown in the displacement and output voltage signals.The displacement verus velocity phase trajectories of the friction system in both directions are shown in Figures12(a) and 12(b). For the case of ignoring separation, the tangential phase trajectory exhibits a typical stick-slip limit cycle, while the phase trajectory in the case of considering separation is very unconventional, in which several limit cycles appear and repeat themselves; in addition, the slider is found to be able to slide on the belt with its speed passing through the speed of the belt without stick motion. For the normal direction, when separation is ignored, the phase trajectory exhibits a larger-amplitude limit cycle compared to the cases of considering separation.Figure 12
The displacement versus velocity phase trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d) whenμk = 0.7 and μs = 0.9.
(a)(b)(c)(d)Furthermore, the velocity versus voltage trajectory is plotted in Figures12(c) and 12(d). Visibly, the trajectory in the tangential direction when considering separation exhibits a larger-amplitude and more complicated limit cycle, compared to the case of ignoring separation. This is attributed to the fact that the velocity of slider can sometimes exceed the belt velocity. As a consequence, the tangential output voltage when considering separation is visibly larger compared to the situations of ignoring separation (Figure 12(c)). Meanwhile, for the trajectory in the normal direction, for the case of ignoring separation, the trajectory exhibits a larger-amplitude limit cycle, which indicates that a sustained and stable voltage can be outputted during the friction process. In contrast, when separation is taken into consideration, the trajectory shows a smaller limit cycle, suggesting that less voltage can be outputted in this situation (Figure 12(d)).


#### 5.2.2.μc<μk<μs

Figure13 shows the simulation results in the case where μk = 1 and μs = 1.2, and the initial condition of system is set at (xe, 0.1, ye, 0, 0, 0). It can be seen that the dynamics and output voltage behaviors show significant difference between these two cases (ignoring and considering separation). For the tangential direction, the displacement signal shows larger amplitude compared with the situation of considering separation. However, the displacement amplitude in the normal direction shows totally opposite results (Figure 13(a) and 13(b)). Consequently, the output voltage in tangential direction when ignoring separation is visibly larger than that in the case of considering separation, while the normal output voltage is smaller compared with the situation of considering separation (Figure 13(c) and 13(d)). This phenomenon is also reflected in the signals of contact force in these two cases (ignoring and considering separation); it can be seen that the contact force in the case of considering separation is larger than that in the case of ignoring separation (Figure 13(e)).Figure 13
The tangential displacement (a), normal displacement (b), tangential voltage (c), normal voltage (d), and contact force (e) of the friction system in the time domain whenμk = 1 and μs = 1.2.
(a)(b)(c)(d)(e)Expectedly, when separation is ignored, the displacement versus velocity trajectory shows larger limit cycle compared to the case of considering separation in the tangential direction, and thus the corresponding velocity versus voltage trajectory when ignoring separation possesses larger limit cycle in this direction, as shown in Figures14(a) and 14(b). This phenomenon suggests that ignoring separation will overestimate vibration level and energy harvesting ability in the tangential direction of friction system. Conversely, when ignoring separation, the displacement versus velocity limit cycle in the normal direction exhibits much smaller amplitude than the case of considering separation, and the corresponding velocity versus voltage trajectory possesses smaller limit cycle in this direction, which indicates that ignoring separation will underestimate the vibration level and energy harvesting ability in the normal direction of friction system, as shown in Figures 14(c) and 14(d).Figure 14
The displacement versus velocity phase trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d) whenμk = 1 and μs = 1.2.
(a)(b)(c)(d)Combining the puzzling results shown in these two situations, a possible physical explanation is provided. When ignoring the separation behavior of the friction interface, the normal spring is stretched and provides continuous resistance to the slider, which accordingly prevents the slider from continuously moving upward, while the friction force still contributes to the system vibration in the tangential direction. In contrast, when separation is considered in this system, the normal contact stiffness will lose its effect when separation occurs; at this time, the friction force disappears, and the normal contact stiffness will not provide resistance force to the normal vibration. Since the tangential and normal vibration is coupled, the vibration of this system when separation is considered is actually more complex. The friction coefficient also plays a significant role in affecting the dynamics and energy harvesting performance of the friction system; visibly, the different friction coefficient will cause complete opposite results. These are very interesting findings and they are useful for designing PEH devices via FIV.


### 5.3. Comparison and Analysis of the Friction System with Different Electric Resistances (R)

In this section, the effect ofR on the dynamic and output voltage is studied. Figure 15 exhibits the dynamical behaviors and output voltages of the friction system under different values of R. It is found that the variation of R will not cause vibration level of friction system generating significant variation. However, R is found to be able to affect the output voltages values in both directions. Visibly, larger R is beneficial for outputting higher voltages within limit. Meanwhile, with R further increasing, the output voltage will not show noticeable variation.Figure 15
The displacement and output voltage of friction system whenR is 10 KΩ (a), 50 KΩ (b), 100 KΩ (c), and 150 KΩ (d).
(a)(b)(c)(d)Additionally, the output voltage amplitudes and average generated electric power versusR in both directions are illustrated in Figure 16. The average generated electric power is evaluated by the following equation [63]:(17)Paverage=∫0TVR2/RdtT.Figure 16
The output voltage amplitude and electric power of friction system with differentR in the tangential direction (a) and normal direction (b).
(a)(b)For the curve of output voltage amplitude (red line), it increases gradually at first and then remains stable basically with the increase ofR. Meanwhile, for the electric power generated from the piezo elements, the curves (black line) show a trend of rising firstly and then falling with the further increase of R, indicating that there is an optimal external resistance that gives the maximum output power, which is in agreement with the experimental results reported by Wang et al. [8]. Moreover, it is worth noting that the external resistance of 40 KΩ in the tangential direction is beneficial for the system to obtain the maximum average output power of 0.2112 mW, while the external resistance of 60 KΩ in the normal direction gives the maximum average output power of 38 mW, which indicates that the optimal external resistance to generate highest power in the tangential direction is relatively less than that in the normal direction. The reason is that the output voltage in the tangential direction can arrive at a stable value at a lower optimal external resistance, compared to that of the output voltage in the normal direction. Therefore, it can be concluded that R has little effect on the system’s dynamical behaviors; however, it can significantly affect the energy harvesting behavior in a certain degree, and appropriate external electric resistances values in both directions are beneficial for obtaining the highest harvested energy.


## 5.1. Friction Coefficient

According to the value of critical friction coefficient (0.85), the dynamic behaviors of the friction system with kinetic and static friction coefficients, denoted byμk and μs, located within three ranges (μk<μs<μc,μk<μc<μs, and μc<μk<μs), are examined and compared. The values of friction coefficients in these three cases are listed in Table 2, and concurrently listed in Table 2 are the Eps of the friction system under different friction coefficients. The initial condition is set at (xe, 0.1, ye, 0, 0, 0), with the tangential velocity of the slider equaling the velocity of the belt. The external electric resistance in this analysis is set at 10 KΩ.Table 2
The values of friction coefficients and the corresponding Eps in these three cases.
CaseμkμsμcEps (xe, ye)10.40.60.85(−0.00216, −0.0258)20.70.90.85(0.00104, −0.0246)311.20.85(0.00397, −0.0235)Figure8 depicts the dynamical behaviors and output voltages of the friction system with three different cases of friction coefficients. Visibly, the vibration level and output voltage are strongly related to the friction coefficient. In the tangential direction, the friction system exhibits the strongest self-excited vibration level when the friction coefficient is located in Case 3. Accordingly, the voltage generated in this situation exhibits the highest value. Meanwhile, for the lowest friction coefficient of Case 1, the tangential vibration amplitude and output voltage are visibly lower than those generated from the other two cases (Figures 8(a) and 8(c)). Similarly, for the normal direction, a larger coefficient is found to be able to generate stronger vibration and higher output voltage as well. Therefore, the vibration level and voltage values generated in Case 3 are larger than those from Case 2, and vibration and voltage signals generated from Case 1 show the lowest values.Figure 8
The tangential displacement (a), normal displacement (b), tangential output voltage (c), and normal output voltage (d) of the friction system with three different cases of friction coefficients.
(a)(b)(c)(d)Figure9 shows the phase plots of friction system with different friction coefficient cases. For Case 1, the tangential stick-slip limit cycle is very complicated and exhibits lots of different cycles of stick-slip motion. Similarly, the normal vibration also shows a multiperiod limit cycle. Meanwhile, for Case 2, the phase trajectory in tangential direction is very unconventional, in which several limit cycles are characterized by appearing and repeating themselves; additionally, the tangential velocity of slider exceeds the belt velocity (v > 0.1 m/s) and causes overshoot phenomenon during the vibration. Additionally, the normal vibration also forms a multiperiod limit cycle. Case 3 represents the situation in which both the static and kinetic friction coefficients are larger than the critical friction coefficient; it is seen that the tangential stick-slip and overshoot always happen. In this case, the friction system has a large-amplitude (limit-cycle) response in both the normal direction and the tangential direction, compared with the other two cases.Figure 9
The displacement verus velocity trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d).
(a)(b)(c)(d)The velocity versus voltage trajectory of the friction system in three cases of friction coefficient is plotted. In the tangential direction, it is visible that the orbit of friction system with Case 3 possesses biggest limit cycle, which suggests that the friction system can generate higher amplitude voltage than the other two cases. Meanwhile the output voltage amplitude in Case 1 shows the smallest limit cycle. Indeed, both trajectories from Case 1 and Case 2 exhibit very complicated orbits. In the normal direction, as expected, the orbit of friction system with Case 3 exhibits the largest amplitude of limit cycle, which indicates that a sustained and high voltage is outputted in this direction during the sliding process. Additionally, the orbit of friction system with Case 1 forms a smallest limit cycle, which reflects that relative lower voltage is generated in this situation. This phenomenon can well reflect the time domain signal of voltage shown in Figure8.To sum up, the dynamic and output voltages are strongly related to the relationship among the kinetic friction coefficientμk, static friction coefficient μs, and critical friction coefficient μc. The friction system can generate stronger vibration and output higher voltage when the friction coefficient is located in the range of μc<μk<μs. In contrast, lower vibration and less voltage are generated when the friction coefficient locates in the range of μk<μs<μc.


## 5.2. Comparison and Analysis of Ignoring/considering Separation

In most previous studies of FIV of small models, it is assumed that the slider is still in contact with belt during vibration process [50, 60]. Although the possibility of loss of contact was mentioned in [61, 62], the effect of separation was not studied. As the unstable FIV caused by the friction force at the contact interface grows, it is important to consider separation during the vibration, when friction temporarily disappears. The loss of contact will affect not only the contact behaviors of friction system but also the energy harvesting behaviors. Therefore, in this section, the dynamics and energy harvesting behaviors of the friction system in both cases, that is, ignoring separation and considering separation between slider and the belt, are analyzed and compared.The flow chart of transient dynamic analysis when the friction system ignores the separation between slider and belt is shown in Figure10. Due to ignoring the separation, the equations of motion are given by equation (2). The iteration process starts from the stick mode, and slider velocity is checked at the end of every single time step during numerical computation. Once the slider velocity is not equal to the belt velocity (0.1 m/s), the iteration process switches to the slide mode, and initial value of slide mode is the last value of the friction system in the stick mode. Then, the slider rubs on the moving belt until the next stick happens.Figure 10
The flow chart of transient dynamic analysis when the friction system ignores the separation between belt and slider.In the previous study, it is found that the separation can occur in Case 2 and Case 3; thus, it is worth noting that comparison between ignoring separation and considering separation will be performed in these two cases.


### 5.2.1.μk<μc<μs

Numerical results withμk = 0.7, μs = 0.9, and initial condition (xe, 0.1, ye, 0, 0, 0) are shown in Figure 11. Visibly, the separation and reattachment behaviors make both the tangential and normal dynamics of friction system different from their counterparts when separation is ignored. For the tangential direction, the displacement mainly reflects in the positive direction when considering separation, while when separation is ignored, the displacement fluctuates between the positive and negative directions, suggesting that the slider moves back and forth with the moving belt. In addition, the amplitude of displacement fluctuation in the case of considering separation is larger than that in the case of ignoring separation, as shown in Figure 11(a); thus the output voltage in this direction is larger in the case of considering separation (Figure 11(c)). Meanwhile, for the normal direction, the situation changes significantly. It can be seen that the normal displacement amplitude when considering separation is smaller compared with the situation of ignoring separation; accordingly the output voltage is smaller when separation is taken into consideration, as shown in Figures 11(b) and 11(d).Figure 11
The tangential displacement (a), normal displacement (b), tangential voltage (c), normal voltage (d), and contact force (e) of the friction system in the time domain whenμk = 0.7 and μs = 0.9.
(a)(b)(c)(d)(e)The contact forces of these two situations are plotted in Figure11(e). For the case when separation is considered, multiple separation events can be observed as the values of contact force vary from zero to positive and then drop to zero during the vibration. Meanwhile, when ignoring separation, the value of contact force gradually grows and becomes negative, and the amplitude is larger than that in the case of considering separation. This phenomenon well reflects the results shown in the displacement and output voltage signals.The displacement verus velocity phase trajectories of the friction system in both directions are shown in Figures12(a) and 12(b). For the case of ignoring separation, the tangential phase trajectory exhibits a typical stick-slip limit cycle, while the phase trajectory in the case of considering separation is very unconventional, in which several limit cycles appear and repeat themselves; in addition, the slider is found to be able to slide on the belt with its speed passing through the speed of the belt without stick motion. For the normal direction, when separation is ignored, the phase trajectory exhibits a larger-amplitude limit cycle compared to the cases of considering separation.Figure 12
The displacement versus velocity phase trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d) whenμk = 0.7 and μs = 0.9.
(a)(b)(c)(d)Furthermore, the velocity versus voltage trajectory is plotted in Figures12(c) and 12(d). Visibly, the trajectory in the tangential direction when considering separation exhibits a larger-amplitude and more complicated limit cycle, compared to the case of ignoring separation. This is attributed to the fact that the velocity of slider can sometimes exceed the belt velocity. As a consequence, the tangential output voltage when considering separation is visibly larger compared to the situations of ignoring separation (Figure 12(c)). Meanwhile, for the trajectory in the normal direction, for the case of ignoring separation, the trajectory exhibits a larger-amplitude limit cycle, which indicates that a sustained and stable voltage can be outputted during the friction process. In contrast, when separation is taken into consideration, the trajectory shows a smaller limit cycle, suggesting that less voltage can be outputted in this situation (Figure 12(d)).


### 5.2.2.μc<μk<μs

Figure13 shows the simulation results in the case where μk = 1 and μs = 1.2, and the initial condition of system is set at (xe, 0.1, ye, 0, 0, 0). It can be seen that the dynamics and output voltage behaviors show significant difference between these two cases (ignoring and considering separation). For the tangential direction, the displacement signal shows larger amplitude compared with the situation of considering separation. However, the displacement amplitude in the normal direction shows totally opposite results (Figure 13(a) and 13(b)). Consequently, the output voltage in tangential direction when ignoring separation is visibly larger than that in the case of considering separation, while the normal output voltage is smaller compared with the situation of considering separation (Figure 13(c) and 13(d)). This phenomenon is also reflected in the signals of contact force in these two cases (ignoring and considering separation); it can be seen that the contact force in the case of considering separation is larger than that in the case of ignoring separation (Figure 13(e)).Figure 13
The tangential displacement (a), normal displacement (b), tangential voltage (c), normal voltage (d), and contact force (e) of the friction system in the time domain whenμk = 1 and μs = 1.2.
(a)(b)(c)(d)(e)Expectedly, when separation is ignored, the displacement versus velocity trajectory shows larger limit cycle compared to the case of considering separation in the tangential direction, and thus the corresponding velocity versus voltage trajectory when ignoring separation possesses larger limit cycle in this direction, as shown in Figures14(a) and 14(b). This phenomenon suggests that ignoring separation will overestimate vibration level and energy harvesting ability in the tangential direction of friction system. Conversely, when ignoring separation, the displacement versus velocity limit cycle in the normal direction exhibits much smaller amplitude than the case of considering separation, and the corresponding velocity versus voltage trajectory possesses smaller limit cycle in this direction, which indicates that ignoring separation will underestimate the vibration level and energy harvesting ability in the normal direction of friction system, as shown in Figures 14(c) and 14(d).Figure 14
The displacement versus velocity phase trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d) whenμk = 1 and μs = 1.2.
(a)(b)(c)(d)Combining the puzzling results shown in these two situations, a possible physical explanation is provided. When ignoring the separation behavior of the friction interface, the normal spring is stretched and provides continuous resistance to the slider, which accordingly prevents the slider from continuously moving upward, while the friction force still contributes to the system vibration in the tangential direction. In contrast, when separation is considered in this system, the normal contact stiffness will lose its effect when separation occurs; at this time, the friction force disappears, and the normal contact stiffness will not provide resistance force to the normal vibration. Since the tangential and normal vibration is coupled, the vibration of this system when separation is considered is actually more complex. The friction coefficient also plays a significant role in affecting the dynamics and energy harvesting performance of the friction system; visibly, the different friction coefficient will cause complete opposite results. These are very interesting findings and they are useful for designing PEH devices via FIV.


## 5.2.1.μk<μc<μs

Numerical results withμk = 0.7, μs = 0.9, and initial condition (xe, 0.1, ye, 0, 0, 0) are shown in Figure 11. Visibly, the separation and reattachment behaviors make both the tangential and normal dynamics of friction system different from their counterparts when separation is ignored. For the tangential direction, the displacement mainly reflects in the positive direction when considering separation, while when separation is ignored, the displacement fluctuates between the positive and negative directions, suggesting that the slider moves back and forth with the moving belt. In addition, the amplitude of displacement fluctuation in the case of considering separation is larger than that in the case of ignoring separation, as shown in Figure 11(a); thus the output voltage in this direction is larger in the case of considering separation (Figure 11(c)). Meanwhile, for the normal direction, the situation changes significantly. It can be seen that the normal displacement amplitude when considering separation is smaller compared with the situation of ignoring separation; accordingly the output voltage is smaller when separation is taken into consideration, as shown in Figures 11(b) and 11(d).Figure 11
The tangential displacement (a), normal displacement (b), tangential voltage (c), normal voltage (d), and contact force (e) of the friction system in the time domain whenμk = 0.7 and μs = 0.9.
(a)(b)(c)(d)(e)The contact forces of these two situations are plotted in Figure11(e). For the case when separation is considered, multiple separation events can be observed as the values of contact force vary from zero to positive and then drop to zero during the vibration. Meanwhile, when ignoring separation, the value of contact force gradually grows and becomes negative, and the amplitude is larger than that in the case of considering separation. This phenomenon well reflects the results shown in the displacement and output voltage signals.The displacement verus velocity phase trajectories of the friction system in both directions are shown in Figures12(a) and 12(b). For the case of ignoring separation, the tangential phase trajectory exhibits a typical stick-slip limit cycle, while the phase trajectory in the case of considering separation is very unconventional, in which several limit cycles appear and repeat themselves; in addition, the slider is found to be able to slide on the belt with its speed passing through the speed of the belt without stick motion. For the normal direction, when separation is ignored, the phase trajectory exhibits a larger-amplitude limit cycle compared to the cases of considering separation.Figure 12
The displacement versus velocity phase trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d) whenμk = 0.7 and μs = 0.9.
(a)(b)(c)(d)Furthermore, the velocity versus voltage trajectory is plotted in Figures12(c) and 12(d). Visibly, the trajectory in the tangential direction when considering separation exhibits a larger-amplitude and more complicated limit cycle, compared to the case of ignoring separation. This is attributed to the fact that the velocity of slider can sometimes exceed the belt velocity. As a consequence, the tangential output voltage when considering separation is visibly larger compared to the situations of ignoring separation (Figure 12(c)). Meanwhile, for the trajectory in the normal direction, for the case of ignoring separation, the trajectory exhibits a larger-amplitude limit cycle, which indicates that a sustained and stable voltage can be outputted during the friction process. In contrast, when separation is taken into consideration, the trajectory shows a smaller limit cycle, suggesting that less voltage can be outputted in this situation (Figure 12(d)).


## 5.2.2.μc<μk<μs

Figure13 shows the simulation results in the case where μk = 1 and μs = 1.2, and the initial condition of system is set at (xe, 0.1, ye, 0, 0, 0). It can be seen that the dynamics and output voltage behaviors show significant difference between these two cases (ignoring and considering separation). For the tangential direction, the displacement signal shows larger amplitude compared with the situation of considering separation. However, the displacement amplitude in the normal direction shows totally opposite results (Figure 13(a) and 13(b)). Consequently, the output voltage in tangential direction when ignoring separation is visibly larger than that in the case of considering separation, while the normal output voltage is smaller compared with the situation of considering separation (Figure 13(c) and 13(d)). This phenomenon is also reflected in the signals of contact force in these two cases (ignoring and considering separation); it can be seen that the contact force in the case of considering separation is larger than that in the case of ignoring separation (Figure 13(e)).Figure 13
The tangential displacement (a), normal displacement (b), tangential voltage (c), normal voltage (d), and contact force (e) of the friction system in the time domain whenμk = 1 and μs = 1.2.
(a)(b)(c)(d)(e)Expectedly, when separation is ignored, the displacement versus velocity trajectory shows larger limit cycle compared to the case of considering separation in the tangential direction, and thus the corresponding velocity versus voltage trajectory when ignoring separation possesses larger limit cycle in this direction, as shown in Figures14(a) and 14(b). This phenomenon suggests that ignoring separation will overestimate vibration level and energy harvesting ability in the tangential direction of friction system. Conversely, when ignoring separation, the displacement versus velocity limit cycle in the normal direction exhibits much smaller amplitude than the case of considering separation, and the corresponding velocity versus voltage trajectory possesses smaller limit cycle in this direction, which indicates that ignoring separation will underestimate the vibration level and energy harvesting ability in the normal direction of friction system, as shown in Figures 14(c) and 14(d).Figure 14
The displacement versus velocity phase trajectory of the friction system in tangential direction (a) and normal direction (b). The velocity versus voltage trajectory of the friction system in tangential direction (c) and normal direction (d) whenμk = 1 and μs = 1.2.
(a)(b)(c)(d)Combining the puzzling results shown in these two situations, a possible physical explanation is provided. When ignoring the separation behavior of the friction interface, the normal spring is stretched and provides continuous resistance to the slider, which accordingly prevents the slider from continuously moving upward, while the friction force still contributes to the system vibration in the tangential direction. In contrast, when separation is considered in this system, the normal contact stiffness will lose its effect when separation occurs; at this time, the friction force disappears, and the normal contact stiffness will not provide resistance force to the normal vibration. Since the tangential and normal vibration is coupled, the vibration of this system when separation is considered is actually more complex. The friction coefficient also plays a significant role in affecting the dynamics and energy harvesting performance of the friction system; visibly, the different friction coefficient will cause complete opposite results. These are very interesting findings and they are useful for designing PEH devices via FIV.


## 5.3. Comparison and Analysis of the Friction System with Different Electric Resistances (R)

In this section, the effect ofR on the dynamic and output voltage is studied. Figure 15 exhibits the dynamical behaviors and output voltages of the friction system under different values of R. It is found that the variation of R will not cause vibration level of friction system generating significant variation. However, R is found to be able to affect the output voltages values in both directions. Visibly, larger R is beneficial for outputting higher voltages within limit. Meanwhile, with R further increasing, the output voltage will not show noticeable variation.Figure 15
The displacement and output voltage of friction system whenR is 10 KΩ (a), 50 KΩ (b), 100 KΩ (c), and 150 KΩ (d).
(a)(b)(c)(d)Additionally, the output voltage amplitudes and average generated electric power versusR in both directions are illustrated in Figure 16. The average generated electric power is evaluated by the following equation [63]:(17)Paverage=∫0TVR2/RdtT.Figure 16
The output voltage amplitude and electric power of friction system with differentR in the tangential direction (a) and normal direction (b).
(a)(b)For the curve of output voltage amplitude (red line), it increases gradually at first and then remains stable basically with the increase ofR. Meanwhile, for the electric power generated from the piezo elements, the curves (black line) show a trend of rising firstly and then falling with the further increase of R, indicating that there is an optimal external resistance that gives the maximum output power, which is in agreement with the experimental results reported by Wang et al. [8]. Moreover, it is worth noting that the external resistance of 40 KΩ in the tangential direction is beneficial for the system to obtain the maximum average output power of 0.2112 mW, while the external resistance of 60 KΩ in the normal direction gives the maximum average output power of 38 mW, which indicates that the optimal external resistance to generate highest power in the tangential direction is relatively less than that in the normal direction. The reason is that the output voltage in the tangential direction can arrive at a stable value at a lower optimal external resistance, compared to that of the output voltage in the normal direction. Therefore, it can be concluded that R has little effect on the system’s dynamical behaviors; however, it can significantly affect the energy harvesting behavior in a certain degree, and appropriate external electric resistances values in both directions are beneficial for obtaining the highest harvested energy.


## 6. Conclusions

A 2-DOF friction system (slider-on-belt model) with piezoelectric elements, which accounts for the stick, slip, separation, and reattachment between contact pairs, is employed to study the energy harvesting performance via friction-induced vibration (FIV). Both the eigenvalues analysis and transient dynamic analysis are carried out to study the dynamics and energy harvesting behavior of this friction system, and the main conclusions are summarized as follows:(1)
Eigenvalue analysis results show that friction coefficient value (μc=0.85) is the critical friction coefficient of the system which will generate self-excited vibration. Transient dynamic analysis results verify the feasibility of piezoelectric energy harvesting (PEH) by means of FIV. A larger friction coefficient is proved to be able to generate more intense vibration and output higher voltage.(2)
The dynamic and output voltages are crucially dependent on the relationship among the kinetic friction coefficientμk, static friction coefficient μs, and critical friction coefficient μk. When the friction coefficient locates in the range of μc<μk<μs, the system will generate stronger vibration and output higher voltage. In contrast, lower vibration and less voltage will be generated when the friction coefficient locates in the range of μk<μs<μc.(3)
If the separation between contact pair is ignored, it is likely to overestimate or underestimate the vibration magnitude and output voltage amplitude. The overestimation or underestimation phenomenon is mainly determined by the located range of friction coefficient. These are very interesting findings and can be exploited in design of energy harvester by means of FIV.(4)
The external resistanceR has little effect on the dynamical behaviors of the friction system, but this will affect the output voltages of friction system within limits. In addition, there exists an appropriate value of R for achieving the largest output voltage and electric power.



---

*Source: 8891756-2020-07-27.xml*