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.
National Natural Science Foundation of China515053961. 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.
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 force FN 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.
If FNt 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:
The situation where slider separates with the sliding belt.
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 displacement yt 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 relationships z1=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.
The parameters of the friction system set in the eigenvalues analysis.
Mass (slider)
0.08 kg
k3 (coupling spring in two directions)
600 N/m
k1 (including the stiffness of Piezo-1)
890 N/m
c1 = 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/V
k2 (linear contact stiffness)
500 N/m
Cp1=Cp2 (clamped capacitance of the piezo elements)
1.79 × 10−7 F
v0 (belt velocity)
5 mm/s∼200 mm/s
F (normal load)
20 N
R1
10 kΩ
R2
10 kΩ
The real parts (a) and imaginary parts (b) of the friction system with the increase of friction coefficient.
4. Transient Dynamic Analysis4.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 velocity vr 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).
When vr>ε, 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 force Ff 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 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.
The tangential displacement (a), normal displacement (b), tangential output voltage (c), and normal output voltage (d) of the friction system with different friction coefficients.
Figure 6 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.
The tangential displacement versus tangential velocity. (a) Normal displacement versus normal velocity. (b) Phase portraits of the friction system.
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 Figure 7. 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.
The open-circuit voltage versus velocity phase trajectories.
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Ω.
The values of friction coefficients and the corresponding Eps in these three cases.
Case
μk
μs
μc
Eps (xe, ye)
1
0.4
0.6
0.85
(−0.00216, −0.0258)
2
0.7
0.9
0.85
(0.00104, −0.0246)
3
1
1.2
0.85
(0.00397, −0.0235)
Figure 8 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.
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.
Figure 9 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.
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).
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 Figure 8.
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 Figure 10. 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.
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).
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.
The contact forces of these two situations are plotted in Figure 11(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 Figures 12(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.
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.
Furthermore, the velocity versus voltage trajectory is plotted in Figures 12(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
Figure 13 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)).
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.
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 Figures 14(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).
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.
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 of R 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.
The displacement and output voltage of friction system when R is 10 KΩ (a), 50 KΩ (b), 100 KΩ (c), and 150 KΩ (d).
Additionally, the output voltage amplitudes and average generated electric power versus R in both directions are illustrated in Figure 16. The average generated electric power is evaluated by the following equation [63]:(17)Paverage=∫0TVR2/RdtT.
The output voltage amplitude and electric power of friction system with different R in the tangential direction (a) and normal direction (b).
For the curve of output voltage amplitude (red line), it increases gradually at first and then remains stable basically with the increase of R. 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:
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.
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.
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.
The external resistance R 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.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflicts of interest with respect to the research, authorship, and/or publication of this article.
Acknowledgments
This research received financial support from the National Natural Science Foundation of China (no. 51505396).
PriyaS.SongH. C.ZhouY.A Review on piezoelectric energy harvesting: materials, methods, and circuits Energy201741339SudevalayamS.KulkarniP.Energy harvesting sensor nodes: survey and implications201113344346110.1109/surv.2011.060710.000942-s2.0-79959289243PriyaS.InmanD. J.DanielJ.2008New York, NY, USASpringerStephenN. G.On energy harvesting from ambient vibration20062931-240942510.1016/j.jsv.2005.10.0032-s2.0-33644933795KimH. S.KimJ.-H.KimJ.A review of piezoelectric energy harvesting based on vibration20111261129114110.1007/s12541-011-0151-32-s2.0-83255189518WeiC.JingX.A comprehensive review on vibration energy harvesting: modelling and Realization20177411810.1016/j.rser.2017.01.0732-s2.0-85013176797AhmedR.MirF.BanerjeeS.A review on energy harvesting approaches for renewable energies from ambient vibrations and acoustic waves using piezoelectricity201726808503110.1088/1361-665x/aa7bfbWangD. W.MoJ. L.WangX. F.OuyangH.ZhouZ. R.Experimental and numerical investigations of the piezoelectric energy harvesting via friction-induced vibration20181711134114910.1016/j.enconman.2018.06.0522-s2.0-85048798245BeebyS. P.TorahR. N.TudorM. J.A micro electromagnetic generator for vibration energy harvesting20071771257126510.1088/0960-1317/17/7/0072-s2.0-34547578774KhanF. U.QadirM. U.State-of-the-art in vibration-based electrostatic energy harvesting2016261010300110.1088/0960-1317/26/10/1030012-s2.0-84989881975ErtukA.InmanD. J.2011Chichester, U.K.John Wiley & SonsWangL.YuanF. G.Vibration energy harvesting by magnetostrictive material200817411010.1088/0964-1726/17/4/0450092-s2.0-51649101733ChenJ.WangZ. L.Reviving vibration energy harvesting and self-powered sensing by a triboelectric nanogenerator20171348052110.1016/j.joule.2017.09.0042-s2.0-85034954735FuY.OuyangH.DavisR. B.Triboelectric energy harvesting from the vibro-impact of three cantilevered beams201912150953110.1016/j.ymssp.2018.11.0432-s2.0-85057354674ErturkA.InmanD. J.On mechanical modeling of cantilevered piezoelectric vibration energy harvesters200819111311132510.1177/1045389x070856392-s2.0-44349190759ErturkA.InmanD. J.An experimentally validated bimorph cantilever model for piezoelectric energy harvesting from base excitations20091802500910.1088/0964-1726/18/2/025009FakhzanM. N.MuthalifA. G. A.Harvesting vibration energy using piezoelectric material: modeling, simulation and experimental verifications2013231616610.1016/j.mechatronics.2012.10.0092-s2.0-84896332850KimM.DugundjiJ.WardleB. L.Efficiency of piezoelectric mechanical vibration energy harvesting201524505500610.1088/0964-1726/24/5/0550062-s2.0-84928636666WangJ.ZhouS.ZhangZ.YurchenkoD.High-performance piezoelectric wind energy harvester with Y-shaped attachments2019181164565210.1016/j.enconman.2018.12.0342-s2.0-85059112916LiuH. C.ZhongJ. W.LeeC. K.LeeS. W.LinL. W.A comprehensive review on piezoelectric energy harvesting technology: materials, mechanisms, and applications2018504130610.1063/1.50741842-s2.0-85059660491ChallaV. R.PrasadM. G.FisherF. T.Towards an autonomous self-tuning vibration energy harvesting device for wireless sensor network applications201120211110.1088/0964-1726/20/2/0250042-s2.0-79951588619SafaeiM.SodanoH. A.AntonS. R.A review of energy harvesting using piezoelectric materials: state-of-the-art a decade later (2008-2018)2019281111300110.1088/1361-665x/ab36e4WangX.ChenC.WangN.A frequency and bandwidth tunable piezoelectric vibration energy harvester using multiple nonlinear techniques201719036837510.1016/j.apenergy.2016.12.1682-s2.0-85008485923SongH.-C.KumarP.SriramdasR.Broadband dual phase energy harvester: vibration and magnetic field20182251132114210.1016/j.apenergy.2018.04.0542-s2.0-85047765534ZhangJ.FangZ.ShuC.ZhangJ.ZhangQ.LiC.A rotational piezoelectric energy harvester for efficient wind energy harvesting2017262112312910.1016/j.sna.2017.05.0272-s2.0-85020134659FarokhiH.GholipourA.GhayeshM. H.Efficient broadband vibration energy harvesting using multiple piezoelectric bimorph202087404100110.1115/1.4045529TangL.YangY.SohC. K.Toward broadband vibration-based energy harvesting201021181867189710.1177/1045389x103902492-s2.0-78650924737ZhouS.CaoJ.ErturkA.LinJ.Enhanced broadband piezoelectric energy harvesting using rotatable magnets20131021717390110.1063/1.48034452-s2.0-84877278091YangZ.TangL.TaoK.AwK.A broadband electret-based vibrational energy harvester using soft magneto-sensitive elastomer with asymmetrical frequency response profile2019281010LT0210.1088/1361-665x/ab3ae1ErturkA.InmanD. J.Broadband piezoelectric power generation on high-energy orbits of the bistable duffing oscillator with electromechanical couplingfing oscillator with electromechanical coupling2011330102339235310.1016/j.jsv.2010.11.0182-s2.0-79952438592YangJ.ChenJ.YangY.Broadband vibrational energy harvesting based on a triboelectric nanogenerator201446130132210.1002/aenm.2013013222-s2.0-84900475725ZhouS.CaoJ.InmanD. J.LinJ.LiuS.WangZ.Broadband tristable energy harvester: modeling and experiment verificationfication2014133333910.1016/j.apenergy.2014.07.0772-s2.0-84905638500WangX.LinL. W.Dimensionless optimization of piezoelectric vibration energy harvesters with different interface circuits201322808501110.1088/0964-1726/22/8/0850112-s2.0-84881123972ShuY. C.LienI. C.Efficiency of energy conversion for a piezoelectric power harvesting systemficiency of energy conversion for a piezoelectric power harvesting system200616112429243810.1088/0960-1317/16/11/0262-s2.0-33750602054DhoteS.YangZ.BehdinanK.ZuJ.Enhanced broadband multi-mode compliant orthoplanar spring piezoelectric vibration energy harvester using magnetic force2018135637110.1016/j.ijmecsci.2017.11.0122-s2.0-85034033127TwiefelJ.WestermannH.Survey on broadband techniques for vibration energy harvesting201324111291130210.1177/1045389x134761492-s2.0-84878365973DingW. J.2010Heidelberg, GermanySpringerIbrahimR. A.Friction-induced vibration, chatter, squeal, and chaos-Part I: mechanics of contact and friction199447720922610.1115/1.31110792-s2.0-84962088295IbrahimR. A.Friction-induced vibration, chatter, squeal, and chaos-part II: dynamics and modeling199447722725310.1115/1.31110802-s2.0-84886454789CantoniC.CesariniR.MastinuG.RoccaG.SiciglianoR.Brake comfort - a review200947890194710.1080/004231109031004322-s2.0-73349111908KangJ.KimK.Squeak noise in lead screw systems: self-excited vibration of continuous model2010329173587359510.1016/j.jsv.2010.03.0182-s2.0-77951112272Le RouzicJ.Le BotA.Perret-LiaudetJ.Friction-induced vibration by stribeck’s law: application to wiper blade squeal noise201349356357210.1007/s11249-012-0100-z2-s2.0-84892874783GhazalyN. M.El-SharkawyM.AhmedI.A review of automotive brake squeal mechanisms20141159SinouJ. J.JezequelL.On the stabilizing and destabilizing effects of damping in a non-conservative pin-disc system20081991–4435210.1007/s00707-007-0565-22-s2.0-50649104553LaiV.-V.ChielloO.BrunelJ.-F.DufrénoyP.Full finite element models and reduction strategies for the simulation of friction-induced vibrations of rolling contact systems201944419721510.1016/j.jsv.2018.12.0242-s2.0-85059807665MassiF.BailletL.CullaA.Structural modifications for squeal noise reduction: numerical and experimental validation2009511-216818910.1504/ijvd.2009.027120AriyatanapolR.XiongY. P.OuyangH.Partial pole assignment with time delays for asymmetric systems201822962619262910.1007/s00707-018-2118-22-s2.0-85045071763StenderM.TiedemannM.HoffmannN.Energy harvesting below the onset of flutter20194581721SinouJ.-J.JézéquelL.Mode coupling instability in friction-induced vibrations and its dependency on system parameters including damping200726110612210.1016/j.euromechsol.2006.03.0022-s2.0-33845303566ShinK.BrennanM. J.OhJ.-E.HarrisC. J.Analysis of disc brake noise using a two-degree-of-freedom model2002254583784810.1006/jsvi.2001.41272-s2.0-0037173820PoppK.StelterP.Stick-slip vibrations and chaos1990332162489105HelsethL. E.Excitation of energy harvesters using stick-slip motion20142308502410.1088/0964-1726/23/8/0850242-s2.0-84904437433TadokoroC.MatsumotoA.NagamineT.SasakiS.Piezoelectric power generation using friction-induced vibration20172606501210.1088/1361-665x/aa70422-s2.0-85019612282MasudaA.SawaiC.Stick-slip energy harvesting: a preliminary studyASME 2017 Conference on Smart Materials, Adaptive Structures and Intelligent SystemsSeptamber 2017Snowbird, Utah, USAHeffelE.HagedornP.Friction induce vibration for energy harvesting applications201314WangD. W.LiuM. X.QianW. J.WuX.MaQ.WuZ. Q.Parametrical investigation of piezoelectric energy harvesting via friction-induced vibration20202020132619021510.1155/2020/6190215LiZ. L.OuyangH.GuanZ. Q.Nonlinear friction-induced vibration of a slider-belt system2016138904100610.1115/1.40332562-s2.0-84974604214NiknamA.FarhangK.Friction-induced vibration due to mode-coupling and intermittent contact loss2019141202101210.1115/1.40416712-s2.0-85056718321MedioA.LinesM.2001Cambridge, UKCambridge University PressHoffmannN.FischerM.AllgaierR.GaulL.A minimal model for studying properties of the mode-coupling type instability in friction induced oscillations200229419720510.1016/s0093-6413(02)00254-92-s2.0-0036630565SoobbarayenK.SinouJ.-J.BessetS.Numerical study of friction-induced instability and acoustic radiation - effect of ramp loading on the squeal propensity for a simplified brake modelfied brake model2014333215475549310.1016/j.jsv.2014.05.0372-s2.0-84926260168BrunettiJ.MassiF.D’AmbrogioW.BerthierY.Dynamic and energy analysis of frictional contact instabilities on a lumped system201550363364710.1007/s11012-014-0020-02-s2.0-84949632780OwensB. A. M.MannB. P.Linear and nonlinear electromagnetic coupling models in vibration-based energy harvesting2012331492293710.1016/j.jsv.2011.10.0262-s2.0-80955137539