Next Article in Journal
High Precision Thermoforming 3D-Conformable Electronics with a Phase-Changing Adhesion Interlayer
Previous Article in Journal
Design of a Micromachined Z-axis Tunneling Magnetoresistive Accelerometer with Electrostatic Force Feedback
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Vibration of a Micro Piezoelectric Precision Drive System

1
School of Mechanical Engineering, Jiangsu University of Science and Technology, Zhenjiang 212003, China
2
Robotics and Microsystems Center, Scoochow University, Suzhou 215021, China
*
Author to whom correspondence should be addressed.
Micromachines 2019, 10(3), 159; https://doi.org/10.3390/mi10030159
Submission received: 17 January 2019 / Revised: 20 February 2019 / Accepted: 22 February 2019 / Published: 26 February 2019

Abstract

:
A micro piezoelectric precision drive system is proposed, which is advantageous due its small size, large transmission ratio, and large output torque. The working principle of the proposed piezoelectric precision drive system is presented, and the nonlinear dynamic model and equations of the system are established. Using the Linz Ted-Poincaré and perturbation methods, the nonlinear approximate solutions of the dynamic equations are calculated. The results indicate that the nonlinear intensity of the drive system is inversely proportional to the number of meshing movable teeth. It was also noted that the rotor is most affected by the nonlinear phenomenon. These results can be utilized both to optimize the dimensions of the piezoelectric precision drive system and to reduce the intensity of vibrations during operation.

1. Introduction

With the development of science and technology, the precision transmission has been widely utilized in society. At the same time, numerous novel transmissions such as the contactless air conveyor system, magnetic coupling transmission, and piezoelectric drive system have been developed [1,2,3,4]. Among them, the piezoelectric drive system is attracting increasing attention from scholars.
The piezoelectric actuator has superior performances, such as a fast response, a large output force, a small size, and no electromagnetic interference. As a result, it has been successfully used in the fields of precision positioning, robot actuating, energy harvesting, and vibration reduction [5,6,7,8].
Numerous developments in piezoelectric driving have been achieved. Ozaki et al. [9] proposed a bioinspired flapping-wing micro aerial vehicle with piezoelectric direct-driven actuation, and a two-wing prototype with a wingspan of 114 mm was designed and fabricated. The total mass of the aerial vehicle was 598 mg, and the maximum measured lift force was 6.52 mN at a driving voltage of 100 V. Utilizing a piezoelectric bimorph, David et al. [10] developed an improved tactile sensor. As the sensor can distinguish soft materials, it was used for brain tumor resection. Besides, Zhang et al. [11] designed a finger joint based on a hybrid multi-degree of freedom (DOF) piezoelectric ultrasonic motor. The stator of the piezoelectric motor consisted of a multi-layered piezoelectric longitudinal vibrator and a sandwich bending vibrator, which could generate a high longitudinal vibration velocity at a low input voltage. Through a testing experiment, a maximum torque of 23.5 mNm was obtained. Tajitsu [12] designed a piezoelectric poly-L-lactic acid fabric. The fabric could feel people’s movement. Thus, it was used for controlling a humanoid robot. Moreover, a novel piezoelectric-driven robot was proposed by Dharmawan et al. [13]. The robot consisted of a piezoelectric unimorph actuator and a four-bar linkage. Driven by the piezoelectric force, two four-bar linkages generated forward and backward movements. Fang et al. [14] proposed a precision position stage with a piezoelectric actuator. In addition, Zou et al. [15] developed an insect-scale flapping-wing robot with a wingspan of 35 mm. The robot had piezoelectric actuators. Li [16] proposed a novel piezoelectric-based unmanned underwater vehicle with the functions of piezoelectric vibration reduction and piezoelectric energy harvesting.
Regarding the nonlinear dynamic characteristic of the piezoelectric system, numerous works have also been conducted. Considering the inertial, geometric, and piezoelectric nonlinearities, Abdelkefi et al. [17] presented a nonlinear calculating method for an energy harvester. Using the Galerkin technique, Gauss law, and extended Hamilton principle, a reduced-order model of the harvester was built. The effect of nonlinear piezoelectric coefficients on the system response was also investigated. Moreover, Meesala et al. [18] developed a method of multiple scales as an approximate solution, as well as amplitude and phase modulation equations, to capture the nonlinear behavior of an energy harvester. With this method, the parameters of the nonlinear piezoelectric constitutive relations were estimated. Takahashi et al. [19] studied the nonlinear behavior of the ferroelectric ceramic, which excited at a resonant mode. It was found that an induced vibration stress had a greater effect on the nonlinearity than an electric field. In addition, Ozaki et al. [20] studied nonlinear piezoelectric vibration in resonant piezoelectric devices. The purpose of the study was to introduce a nonlinear parameter to the finite elements method (FEM) and to establish a method for measuring the nonlinear parameter by evaluating a nonlinear model for the piezoelectric vibration. The abovementioned studies have focused on the nonlinear vibration of the piezoelectric material itself. However, the effects of the piezoelectric nonlinear characteristic on the driving system have had fewer studies.
In this study, a micro piezoelectric precision drive system is proposed. The system integrates piezoelectric driving and a movable tooth transmission, and it has the advantages of being small in size, with a large transmission ratio, and a large output torque. Compared with the electromagnetic drive, it has merits such as a compact structure, a fast response, electromagnetic interference resistance. The nonlinear dynamic performance results in a drive system with a deteriorated load-carrying capability. For the proposed piezoelectric precision drive system, the nonlinear free vibration of the transmission part [21] and the nonlinear dynamics of the driving part under piezoelectric excitation [22] are investigated. Meanwhile, using the numerical method, the chaotic vibrations of the driving part, transmission part, and nonlinear axis force are analyzed separately [23,24,25]. In [24], the chaotic vibration of a movable tooth system was a special nonlinear vibration, but emphasis was placed on the chaotic characteristic and the influence of parameters on chaotic vibrations. However, the nonlinear response characteristic of a movable tooth transmission under piezoelectric excitation was not investigated. Unlike previous works, this new research mainly focuses on the nonlinear solution process, the influence law of parameters on the frequency characteristic, and force vibration rules, which are the different from the previous chaotic characteristics.
The aim of this study is to reveal the nonlinear vibration response law under the nonlinear piezoelectric effect. Using the dynamic theory of a discrete system, the nonlinear dynamic model of a robot-used piezoelectric precision drive system is built. With the Linz Ted-Poincaré and perturbation methods, the approximate nonlinear dynamic responses are solved. Furthermore, the effects of the system parameters on the nonlinear natural characteristic and dynamic responses are analyzed. These results can be used to reduce the vibration of the proposed drive system.

2. Operating Principle of the Drive System

The operating principle of the proposed micro piezoelectric precision drive system is shown in Figure 1. It consists of the following parts: Piezoelectric actuator (1), displacement magnifying mechanism (2), swaying rod (3), rotor (4), central gear (5), movable tooth (6), harmonic generator (7) and preload spring (8). The system is driven by two piezoelectric actuators, with dimensions of 5 mm × 5 mm × 30 mm, which are installed at the lower shell, with a phase difference of 90°.
In the initial state, the two piezoelectric actuators maintain the original length, the harmonic rod contacts tightly with the movable teeth under the force of the preloaded spring, and the swaying rod sways toward the displacement magnifying mechanism - the random movable tooth located at the A-site in Figure 1a.
When the proposed piezoelectric precision drive system works, sinusoidal signals with a 90° phase difference are applied to the two piezoelectric actuators. Then, the actuators are lengthened in their axis direction. Driven by the deformation force, a harmonic wave is generated under the displacement amplification of the swaying rod. The harmonic force of the harmonic generator pushes the movable teeth slide along the central gear profile.
When the swaying rod sways to the 90° position, the movable tooth moves into the middle position of the adjacent addendum and root (see Figure 1b) and then moves into the adjacent roots at the 180° position (see Figure 1c). Similarly, the swaying rod sways to the 360° position and returns to the initial position, the movable tooth moves into B-site of the adjacent addendum in Figure 1a, and it accomplishes the φ angle movement of a period. Thus, driven by the continuous input signal, a continuous rotation of the rotor can be generated.

3. Dynamic Models and Equations

The driving source of the proposed micro piezoelectric precision drive system is derived from the two piezoelectric actuators. Thus, its exciting force is piezoelectric excitation. From the piezoelectric equations [26], the output force of the piezoelectric actuator can be obtained. The forces of the displacement amplification mechanism and the swaying rod are shown in Figure 2, where Fc is the exciting force of the movable tooth system. From the force balance theorem, the force Fc can be written as:
F C = ( F p l 1 M O ) l 5 k δ k l 3 ( l 5 + l 6 ) l 3 ( l 5 + l 6 + l 7 ) = [ ( 2 l p c 33 s 33 P + c 33 d 33 A p U p ) l 1 2 l p M O ] l 5 2 l p k δ k l 3 ( l 5 + l 6 ) l 3 ( l 5 + l 6 + l 7 ) ,
where c33, d33, and s33 are elasticity modulus, elastic flexibility factor, and piezoelectric strain constant of the piezoelectric actuator, respectively; P is preload of piezoelectric actuator; lp and Ap are the per piezoelectric pieces thickness and cross-sectional area of piezoelectric actuator, respectively; Up is the exciting voltage, U p = U p p [ 1 + sin ( ω t ) ] , U p p is the peak-peak value of voltage, ω is the exciting frequency; Mo is the torque of O point; li is the length of each part, here, l1 = 6 mm, l2 = 11 mm, l3 = 35 mm, l5 = 6 mm, l6 = 16 mm, l7 = 44 mm; k and δk are the stiffness and deformation of preloaded spring, respectively.
Figure 3 shows the dynamic model of the drive system, where subscripts s, c, r and p represent the harmonic rod, central gear, rotor and movable tooth, respectively. Variables xj, yj, uj represent the x direction, y direction, and circumferential direction linear displacement, respectively. The nonlinear dynamic equations were built in [24], and the modeling process is simplified in this report. The nonlinear equations are used to analyze the nonlinear response characteristic. Substituting the nonlinear meshing stiffness into the dynamic model, the nonlinear dynamic equation of the drive system can be given as follows:
M q ¨ + K q = F + Δ F + Δ F e ,
where M, K, and q are the mass matrix, stiffness matrix and generalized coordinate array, respectively. F, ΔF, and ΔFe are the outer force array, nonlinear meshing force increment array and exciting meshing force array, respectively (here, the force comes from the piezoelectric actuator).
The regularization dynamic equation can be written as:
Δ q ¨ N + K N Δ q N = Δ F N + Δ F e N ,
where KN is regular stiffness matrix, ΔFN and ΔFeN are the regular meshing force array and regular outer force array, respectively, and Δ F N = A N T Δ F = B k u i ε P N T , ΔqN is the regular coordinate array.
First, the nonlinear free vibration was investigated through the Linz Ted-Poincaré method. In this condition, the outer force array was zero, that is ΔFeN = 0.
Assuming a regular coordinate array ΔqN, natural frequency ωi can be expressed in the following forms:
{ Δ q N = q N 0 + ε q N 1 + ε 2 q N 2 + ω i 2 = ω 0 i 2 ( 1 + ε σ 1 + ε 2 σ 2 + ) ,
Substituting Equation (9) into Equation (8), with ΔFeN = 0, the nonlinear dynamic equation can be changed into a linear equation set as follows:
{ q ¨ N 0 + ω 0 i 2 q N 0 = 0 q ¨ N 1 + ω 0 i 2 q N 1 = σ 1 q ¨ N 0 + P N B k u 0 i q ¨ N 2 + ω 0 i 2 q N 2 = σ 1 q ¨ N 1 σ 2 q ¨ N 0 + σ 1 P N B k u 1 i ,
Solving the above equation set with a zero initial condition, the solutions of the zero-order equation, first-order equation and second-order equation can be expressed by:
q N 0 i = A N 0 i cos ( ω 0 i t ) ,
{ q N 1 1 = P N 1 B k ( A N 1 2 A N 0 2 cos ω 02 t cos ω 01 t ω 01 2 ω 02 2 + A N 0 3 A N 1 3 cos ω 03 t cos ω 01 t ω 01 2 ω 03 2 ) q N 1 2 = P N 2 B k ( A N 0 1 A N 1 1 cos ω 02 t cos ω 01 t ω 01 2 ω 02 2 + A N 0 3 A N 1 3 cos ω 03 t cos ω 02 t ω 03 2 ω 02 2 ) q N 1 3 = P N 3 B k ( A N 0 1 A N 1 1 cos ω 03 t cos ω 01 t ω 01 2 ω 03 2 + A N 0 2 A N 1 2 cos ω 03 t cos ω 02 t ω 02 2 ω 03 2 ) ,
{ q N 2 1 = P N 1 2 B k 3 A N 1 1 ( cos ω 02 t cos ω 01 t ) ω 02 2 ω 01 2 ( P N 2 A N 0 1 A N 1 1 A N 1 2 ω 02 2 ω 01 2 + P N 2 A N 0 3 A N 1 2 A N 1 3 ω 02 2 ω 03 2 + P N 3 A N 0 2 A N 1 2 A N 1 3 ω 02 2 ω 03 2 ) + P N 1 2 B k 3 A N 1 1 ( cos ω 03 t cos ω 01 t ) ω 03 2 ω 01 2 ( P N 2 A N 0 3 A N 1 2 A N 1 3 ω 03 2 ω 02 2 + P N 3 A N 0 1 A N 1 1 A N 1 3 ω 03 2 ω 01 2 + P N 3 A N 0 2 A N 1 2 A N 1 3 ω 03 2 ω 02 2 ) q N 2 2 = P N 2 2 B k 3 A N 1 2 ( cos ω 02 t cos ω 01 t ) ω 02 2 ω 01 2 ( P N 1 A N 0 2 A N 1 2 A N 1 1 ω 01 2 ω 02 2 + P N 1 A N 0 3 A N 1 1 A N 1 3 ω 01 2 ω 03 2 + P N 3 A N 0 1 A N 1 1 A N 1 3 ω 02 2 ω 03 2 ) + P N 2 2 B k 3 A N 1 2 ( cos ω 03 t cos ω 02 t ) ω 03 2 ω 02 2 ( P N 1 A N 0 3 A N 1 1 A N 1 3 ω 03 2 ω 01 2 + P N 3 A N 0 1 A N 1 1 A N 1 3 ω 03 2 ω 02 2 + P N 3 A N 0 2 A N 1 2 A N 1 3 ω 03 2 ω 02 2 ) q N 2 3 = P N 3 2 B k 3 A N 1 3 ( cos ω 02 t cos ω 03 t ) ω 02 2 ω 03 2 ( P N 1 A N 0 2 A N 1 2 A N 1 1 ω 02 2 ω 01 2 + P N 3 A N 0 1 A N 1 1 A N 1 2 ω 02 2 ω 01 2 + P N 3 A N 0 3 A N 1 2 A N 1 3 ω 02 2 ω 03 2 ) + P N 1 2 B k 3 A N 1 3 ( cos ω 01 t cos ω 03 t ) ω 01 2 ω 03 2 ( P N 1 A N 0 3 A N 1 2 A N 1 1 ω 01 2 ω 02 2 + P N 1 A N 0 3 A N 1 1 A N 1 3 ω 01 2 ω 03 2 + P N 3 A N 0 1 A N 1 2 A N 1 1 ω 01 2 ω 02 2 ) ,
From q N 0 i , q N 1 i and q N 2 i , the response displacement can be written as:
q = A N ( q N 0 i + ε q N 1 i + ε 2 q N 2 i ) ,
From the first-order equation and second-order equation, to eliminate the secular term, σ 1 i can be achieved:
σ 1 i = A N 1 i P N i B k ω 0 i 2 ,
{ σ 2 1 = P N 1 B k A N 1 1 ω 01 2 A N 0 1 [ ( P N 1 2 B k 2 A N 1 1 ω 01 2 P N 1 B k ) ( A N 0 2 A N 1 2 ω 01 2 ω 02 2 + A N 0 3 A N 1 3 ω 01 2 ω 03 2 ) + P N 1 B k 2 A N 0 1 A N 1 1 ω 01 2 ( P N 2 ω 01 2 ω 02 2 + P N 3 ω 01 2 ω 03 2 ) ] σ 2 2 = P N 2 B k A N 1 2 ω 2 2 A N 0 2 [ ( P N 2 2 B k 2 A N 1 2 ω 02 2 P N 2 B k ) ( A N 0 1 A N 1 1 ω 02 2 ω 01 2 + A N 0 3 A N 1 3 ω 02 2 ω 03 2 ) + P N 2 B k 2 A N 0 2 A N 1 2 ω 02 2 ( P N 1 ω 02 2 ω 01 2 + P N 3 ω 02 2 ω 03 2 ) ] σ 2 3 = P N 3 B k A N 1 3 ω 03 2 A N 0 3 [ ( P N 3 2 B k 2 A N 1 3 ω 03 2 P N 3 B k ) ( A N 0 1 A N 1 1 ω 03 2 ω 01 2 + A N 0 2 A N 1 2 ω 03 2 ω 02 2 ) + P N 3 B k 2 A N 0 3 A N 1 3 ω 03 2 ( P N 1 ω 03 2 ω 01 2 + P N 2 ω 03 2 ω 02 2 ) ] ,
Thus, the nonlinear natural frequency can be written in the following form:
ω i 2 = ω 0 i 2 ( 1 + ε σ 1 + ε 2 σ 2 ) ,
In Equation (3), when ΔFeN is not a zero value, the equation changes into a nonlinear forced vibration equation. Using the perturbation method, the nonlinear forced vibration equation can be solved. Firstly, transfer Equation (3) into a linear equation set as follows:
{ q ¨ N 0 + ω 0 i 2 q N 0 = 0 q ¨ N 1 + ω 0 i 2 q N 1 = σ 1 q ¨ N 0 + P N B k u 0 i + C k cos ( ω e t ) P e N q ¨ N 2 + ω 0 i 2 q N 2 = σ 1 q ¨ N 1 σ 2 q ¨ N 0 + σ 1 P N B u 1 i + C k σ 1 cos ( ω e t ) P e N ,
where Ck is constant parameter.
Solving the above linear equation set, the solution can be given as:
q N 0 i = A N 0 i cos ( ω 0 i t ) ,
{ q N 1 1 = P N 1 B k ( A N 1 2 A N 0 2 cos ω 02 t cos ω 01 t ω 01 2 ω 02 2 + A N 0 3 A N 1 3 cos ω 03 t cos ω 01 t ω 01 2 ω 03 2 ) + C k P e N 1 ω 01 2 ( cos ω e t cos ω 01 t ) ω 01 2 ω e 2 q N 1 2 = P N 2 B k ( A N 0 1 A N 1 1 cos ω 02 t cos ω 01 t ω 01 2 ω 02 2 + A N 0 3 A N 1 3 cos ω 03 t cos ω 02 t ω 03 2 ω 02 2 ) + C k P e N 2 ω 02 2 ( cos ω e t cos ω 02 t ) ω 02 2 ω e 2 q N 1 3 = P N 3 B k ( A N 0 1 A N 1 1 cos ω 03 t cos ω 01 t ω 01 2 ω 03 2 + A N 0 2 A N 1 2 cos ω 03 t cos ω 02 t ω 02 2 ω 03 2 ) + C k P e N 3 ω 03 2 ( cos ω e t cos ω 03 t ) ω 03 2 ω e 2 ,
{ q N 2 1 = P N 1 B k 2 A N 1 2 ( cos ω 02 t cos ω 01 t ) ω 02 2 ω 01 2 [ P N 2 A N 0 1 A N 1 1 ω 02 2 ω 01 2 + ( P N 2 A N 0 3 + P N 3 A N 0 2 ) A N 1 3 ω 02 2 ω 03 2 ] + [ P N 3 A N 0 1 A N 1 1 ω 03 2 ω 01 2 + ( P N 2 A N 0 3 + P N 3 A N 0 2 ) A N 1 2 ω 03 2 ω 02 2 ] P N 1 B k 2 A N 1 3 ( cos ω 03 t cos ω 01 t ) ω 03 2 ω 01 2 + P N 1 B k P e N 2 C k ( cos ω 02 t cos ω 01 t ) ( ω 01 2 ω 02 2 ) ( ω 02 2 ω e 2 ) + P N 1 B k P e N 3 C k ( cos ω 03 t cos ω 01 t ) ( ω 03 2 ω 01 2 ) ( ω e 2 ω 03 2 ) + ( P e N 1 + P N 1 B k P e N 1 ω 02 2 ω e 2 + P N 1 B k P e N 3 ω 03 2 ω e 2 ) C k ( cos ω e t cos ω 01 t ) ω 01 2 ω e 2 q N 2 2 = P N 2 B k 2 A N 1 1 ( cos ω 02 t cos ω 01 t ) ω 02 2 ω 01 2 [ P N 1 A N 0 2 A N 1 2 ω 01 2 ω 02 2 + ( P N 1 A N 0 3 + P N 3 A N 0 1 ) A N 1 3 ω 01 2 ω 03 2 ] + [ P N 3 A N 0 2 A N 1 2 ω 03 2 ω 02 2 + ( P N 1 A N 0 3 + P N 3 A N 0 1 ) A N 1 1 ω 03 2 ω 01 2 ] P N 2 B k 2 A N 1 3 ( cos ω 03 t cos ω 02 t ) ω 03 2 ω 02 2 + P N 2 B k P e N 1 C k ( cos ω 02 t cos ω 01 t ) ( ω 01 2 ω 02 2 ) ( ω 01 2 ω e 2 ) + P N 2 B k P e N 3 C k ( cos ω 03 t cos ω 02 t ) ( ω 03 2 ω 02 2 ) ( ω e 2 ω 03 2 ) + ( P e N 2 + P N 2 B k P e N 1 ω 01 2 ω e 2 + P N 2 B k P e N 3 ω 03 2 ω e 2 ) C k ( cos ω e t cos ω 02 t ) ω 02 2 ω e 2 q N 2 3 = P N 3 B k 2 A N 1 2 ( cos ω 02 t cos ω 03 t ) ω 02 2 ω 03 2 [ P N 3 A N 0 3 A N 1 3 ω 02 2 ω 03 2 + ( P N 1 A N 0 2 + P N 3 A N 0 1 ) A N 1 1 ω 02 2 ω 01 2 ] + [ P N 1 A N 0 3 A N 1 3 ω 01 2 ω 03 2 + ( P N 1 A N 0 2 + P N 3 A N 0 1 ) A N 1 2 ω 01 2 ω 03 2 ] P N 3 B k 2 A N 1 1 ( cos ω 03 t cos ω 01 t ) ω 03 2 ω 01 2 + P N 3 B k P e N 2 C k ( cos ω 03 t cos ω 02 t ) ( ω 02 2 ω 03 2 ) ( ω 02 2 ω e 2 ) + P N 3 B k P e N 1 C k ( cos ω 03 t cos ω 01 t ) ( ω 01 2 ω 03 2 ) ( ω e 2 ω 01 2 ) + ( P e N 3 + P N 3 B k P e N 1 ω 01 2 ω e 2 + P N 3 B k P e N 2 ω 02 2 ω e 2 ) C k ( cos ω e t cos ω 03 t ) ω 03 2 ω e 2 ,
Substituting Equations (14)–(16) into Equation (9), the displacement of nonlinear forced vibration can be acquired.

4. Results and Discuss

4.1. Nonlinear Free Vibration Analysis

Table 1 presents the system parameters of the drive system. By substituting the parameters into the frequency equation, Equation (12), the nonlinear natural frequencies of the proposed piezoelectric precision drive system were found (see Table 2). Meanwhile, by selecting the harmonic offset a, movable tooth radius rp, and harmonic plate radius rs as research parameters, the effect of the parameters on nonlinear natural frequencies was investigated, as presented in Figure 4. In Table 2 and Figure 4, ωi and Δωi represent nonlinear natural frequencies and the frequency difference between nonlinear and linear frequencies. From Table 2 and Figure 4, it can be observed that:
(1) Under the influence of the nonlinear effect, the first-order and second-order frequencies are larger than their linear natural frequencies, whereas the third-order frequency is smaller than its linear natural frequency. As the frequency order increases, the frequency difference between the nonlinear and linear frequencies grows. When the number of meshing movable teeth is 16, the frequency difference is less than that with 15 meshing teeth. Hence, the less meshing teeth there are, the stronger the nonlinear phenomenon.
(2) As the harmonic offset a grows, the nonlinear frequency of each order increases. This is because the meshing forces between the movable tooth and each element become larger with an increase in a, resulting in growth of their meshing stiffness. For the first-order and second-order conditions, the frequency difference Δωi increases when a increases, while it decreases for the third-order condition. Thus, the nonlinear phenomenon of the first-order and second-order frequencies strengthens with the increase of harmonic offset a.
(3) The variation rules of nonlinear frequencies ωi and frequency difference Δωi are similar to the change in movable tooth radius rp and harmonic plate radius rs. As rp and rs grow, the nonlinear frequencies decrease, and the reason for this lies in the decreasing meshing stiffness. The frequency differences of the first-order and second-order conditions are reduced with the increase in rp and rs, while the value is increased for the third-order condition. Hence, the nonlinear intensities of the first-order and second-order conditions are directly proportional to changes in rp and rs, while that of the third-order condition is inversely proportional.
The influence of parameters on nonlinear free vibration responses was also investigated. For this purpose, the parameters of the number of meshing movable teeth f, harmonic offset a, movable tooth radius rp, and harmonic plate radius rs were selected to be investigated. The acquired results are shown in Figure 5, Figure 6, Figure 7 and Figure 8. In the figures, Δus, Δuc, and Δur represent the torsional responses of the harmonic rod, central gear, and rotor. In Figure 5, Figure 6, Figure 7 and Figure 8, it can be observed that:
(1) The number of meshing movable teeth f has an important influence on the nonlinear displacement responses of the novel drive system. When the number of meshing movable teeth f changes from 16 to 15, the amplitude of Δus reduces by a third, whereas the values of Δuc and Δur increase 5- and 6-fold. Hence, with less meshing movable teeth, the nonlinear torsional vibration of the harmonic rod weakens, and the nonlinear torsional vibration of the central gear and rotor increases.
(2) As the parameter a increases, the amplitude of Δus increases, while that of Δuc and Δur decreases. Thus, the nonlinear vibration of the harmonic rod enhances as a increases, and those of the central gear and rotor weaken. From the large variation of Δuc and Δur, the influence of the nonlinear effect on the central gear and rotor is larger than that of the harmonic rod.
(3) The effect laws of movable tooth radius rp and harmonic plate radius rs on the displacement response are similar. As rp and rs increase, the amplitude of Δus decreases, the amplitude of Δuc firstly increases and then decreases, and that of Δuc firstly decreases and then increases. In addition, the range of Δus is smaller to that of Δuc and Δur.
In summary, under all conditions, the effect of the number of meshing movable teeth f on the displacement response is the most obvious. Further, for the nonlinear torsional vibration of each element, the rotor is most affected by the nonlinear phenomenon.

4.2. Nonlinear Forced Vibration Analysis

From Equations (14)–(16) and Equation (9), the nonlinear forced vibration responses of the proposed drive system can be determined. When the outer exciting frequency (ωe = 1000 rad/s) was far from the resonance frequency, this effected the displacement responses of the drive system, as shown in Figure 9. When the outer exciting frequency was close to the first-order resonance frequency, this effected the close resonance responses, which are presented in Figure 10. From Figure 9 and Figure 10, it is known that:
(1) When the outer exciting frequency is far from the resonance frequency, the amplitude of each element is quite small.
(2) The translational vibration displacements of the harmonic rod and movable teeth are larger than those of other elements. In addition, these characteristics are consistent with the motion feature of the proposed piezoelectric precision drive system. When the drive system operates, the swaying movement of harmonic rod occurs. Then, the harmonic forces push the movable teeth to generate two-dimensional motion.
(3) When the outer exciting frequency is close to the resonance frequency, the response amplitudes in the xs and ys directions are obviously larger than that far from the resonance frequency. At this point, resonances in the xs and ys directions occur, and the resonance in the xs direction is much more severe. Hence, the vibration mode is translational vibration of the harmonic rod for the first-order resonance.

4.3. Simulation Verification Analysis

To verify the correctness of the theoretical model, a finite element method (FEM) software, ANSYS, was used to simulate the nonlinear frequencies of the system. When building the 3D model, the number of meshing movable teeth was set to 16. The FEM model and first-order vibration mode are shown in Figure 11. The theoretical and simulated nonlinear natural frequencies are presented in Table 3.
The results show that the first-order and second-order simulated frequencies were larger than their theoretical values, while the third-order simulation frequency was smaller than its theoretical value. As the frequency order increased, the errors between the theoretical and simulated frequencies were reduced. The largest error was less than 4%.

5. Conclusions

In this study, a micro piezoelectric precision drive system was proposed, and its operating principle was presented. Considering the nonlinear effect of the change in the number of meshing movable teeth, the nonlinear dynamic model and equations of the proposed drive system were established. Using the Linz Ted-Poincaré and perturbation methods, the solutions of the nonlinear responses were calculated. The results indicate that:
(1) The nonlinear intensity of the drive system is inversely proportional to the number of meshing movable teeth.
(2) Regarding the nonlinear torsional vibration of each element, the rotor is most affected by the nonlinear phenomenon.
(3) For the first-order resonance, translational vibration is obvious in the system.
These results can be used for optimizing the structure and improving the dynamic properties of the proposed micro piezoelectric precision drive system.

Author Contributions

Conceptualization, C.L.; Investigation, W.Z.; Validation, J.F.; Writing—original draft, C.L.; Writing—review and editing, L.S.

Funding

This research was funded by the China Postdoctoral Science Foundation (grant number 2018M640515), the Research Project of State Key Laboratory of Mechanical System and Vibration (grant number MSV201808), the Natural Science Foundation of the Jiangsu Higher Education Institutions (grant number 18KJB460008) and the Science and Technology Innovation Project of Jiangsu University of Science and Technology for Youth.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, X.R.; Zhong, W.; Li, C.; Fang, J.; Liu, F. Development of a contactless air conveyor system for transporting and positioning planar objects. Micromachines 2018, 9, 487. [Google Scholar] [CrossRef] [PubMed]
  2. Zhong, W.; Ji, X.; Li, C.; Fang, J.; Liu, F. Determination of permeability and inertial coefficients of sintered metal porous media using an isothermal chamber. Appl. Sci. 2018, 8, 1670. [Google Scholar] [CrossRef]
  3. Li, C.; Xing, J.C. Dynamic optimization of an electromechanical integrated harmonic piezoelectric motor. J. Mech. Sci. Technol. 2018, 32, 2517–2526. [Google Scholar] [CrossRef]
  4. Li, C. Design and mechanical analysis of a mesh-type rotary harmonic piezoelectric motor. Int. J. Appl. Electromagnet. Mech. 2018, 56, 461–478. [Google Scholar] [CrossRef]
  5. Zhu, X.B.; Wen, Z.J.; Chen, G.Z.; Liang, J.; Liu, P. A decoupled flexure-based rotary micropositioning stage with compact size. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2018, 232, 4167–4179. [Google Scholar] [CrossRef]
  6. Yang, Y.L.; Lou, J.Q.; Wu, G.H.; Wei, Y.D.; Fu, L. Design and position/force control of an S-shaped MFC microgripper. Sens. Actuators A Phys. 2018, 282, 63–78. [Google Scholar] [CrossRef]
  7. Yoon, S.; Kim, J.; Cho, K.; Ko, Y.H.; Lee, S.K.; Koh, J.H. Piezoelectric energy generators based on spring and inertial mass. Materials 2018, 11, 2163. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, Y.X.; Niu, S.K.; Hu, Y.T. A periodic piezoelectric smart structure with the integrated passive/active vibration-reduction performances. Smart Mater. Struct. 2017, 26, 065015. [Google Scholar] [CrossRef]
  9. Ozaki, T.; Hamaguchi, K. Bioinspired flapping-wing robot with direct-driven piezoelectric actuation and its takeoff demonstration. IEEE Robot. Autom. 2018, 3, 4217–4224. [Google Scholar] [CrossRef]
  10. David, O.U.; Johan, S.; Ralf, S. Improved tactile resonance sensor for robotic assisted surgery. Mech. Syst. Signal Process 2018, 99, 600–610. [Google Scholar] [Green Version]
  11. Zhang, X.F.; Zhang, G.B.; Nakamura, K.; Ueha, S. A robot finger joint driven by hybrid multi-DOF piezoelectric ultrasonic motor. Sens. Actuators A Phys. 2011, 169, 206–210. [Google Scholar] [CrossRef]
  12. Tajitsu, Y. Piezoelectric Poly-L-lactic acid fabric and its application to control of humanoid robot. Ferroelectrics 2017, 515, 44–58. [Google Scholar] [CrossRef]
  13. Dharmawan, A.G.; Hariri, H.H.; Soh, G.S.; Foong, S.; Wood, K.L. Design, analysis, and characterization of a two-legged miniature robot with piezoelectric-driven four-bar linkage. J. Mech. Robot. 2018, 10, 021003. [Google Scholar] [CrossRef]
  14. Fang, J.W.; Zhang, L.F.; Long, Z.L.; Wang, M.Y. Fuzzy adaptive sliding mode control for the precision position of piezo-actuated nano positioning stage. Int. J. Precis. Eng. Manuf. 2018, 19, 1447–1456. [Google Scholar] [CrossRef]
  15. Zou, Y.; Zhang, W.; Ke, X.; Lou, X.; Zhou, S. The design and micro fabrication of a sub 100 mg insect-scale flapping-wing robot. Micro Nano Lett. 2017, 12, 297–300. [Google Scholar] [CrossRef]
  16. Li, C. Frequency sensitivity and mode characteristics of body structure for a piezoelectric-based UUV. Appl. Sci. 2018, 8, 1774. [Google Scholar] [CrossRef]
  17. Abdelkefi, A.; Nayfeh, A.H.; Hajj, M.R. Effects of nonlinear piezoelectric coupling on energy harvesters under direct excitation. Nonlinear Dyn. 2012, 67, 1221–1232. [Google Scholar] [CrossRef]
  18. Meesala, V.C.; Hajj, M.R. Identification of nonlinear piezoelectric coefficients. J. Appl. Phys. 2018, 124, 065112. [Google Scholar] [CrossRef]
  19. Takahashi, S.; Yamamoto, M.; Sasaki, Y. Nonlinear piezoelectric effect in ferroelectric ceramics. Jpn. J. Appl. Phys. Part 1 Regul. Pap. Short Note Rev. Pap. 1996, 37, 5292–5296. [Google Scholar] [CrossRef]
  20. Ozaki, R.; Liu, Y.; Hosaka, H.; Morita, T. Piezoelectric nonlinear vibration focusing on the second-harmonic vibration mode. Ultrasonics 2018, 82, 233–238. [Google Scholar] [CrossRef] [PubMed]
  21. Li, C.; Xu, L.Z.; Gao, L.C. Nonlinear free vibration of driving system of an electromechanical integrated harmonic piezoelectric motor. J. Vib. Shock 2016, 35, 7–12. [Google Scholar]
  22. Li, C.; Xing, J.C.; Xu, L.Z.; Sun, L. Nonlinear dynamics of a piezoelectric precision drive system used for robot joint. Int. J. Appl. Electromagnet. Mech. 2018, 58, 57–77. [Google Scholar] [CrossRef]
  23. Li, C.; Xing, J.C. Effects of nonlinear axial forces on chaotic vibration for a mesh-type harmonic piezoelectric motor. Int. J. Appl. Electromagnet. Mech. 2018, 56, 505–622. [Google Scholar] [CrossRef]
  24. Li, C.; Xing, J.C.; Fang, J.W.; Zhao, Z. Numerical analysis on chaotic vibration of drive system for a movable tooth piezoelectric motor. Shock Vib. 2017, 2017, 3216010. [Google Scholar] [CrossRef]
  25. Li, C.; Xu, L.Z.; Zhang, J.Q. Bifurcation and chaotic vibration for an electromechanical integrated harmonic piezodrive system. J. Mech. Sci. Technol. 2016, 30, 2961–2970. [Google Scholar] [CrossRef]
  26. Zhao, C.S. Ultrasonic Motors Technologies and Applications; Science Press: Beijing, China, 2010. [Google Scholar]
Figure 1. Operating principle of the micro piezoelectric precision drive system. (a) Initial position and 360° position; (b) 90° position; (c) 180° position; (d) 270° position.
Figure 1. Operating principle of the micro piezoelectric precision drive system. (a) Initial position and 360° position; (b) 90° position; (c) 180° position; (d) 270° position.
Micromachines 10 00159 g001
Figure 2. Forces of displacement amplification mechanism and swaying rod.
Figure 2. Forces of displacement amplification mechanism and swaying rod.
Micromachines 10 00159 g002
Figure 3. Dynamic model of the drive system. (a) Movable tooth structure; (b) model simplification.
Figure 3. Dynamic model of the drive system. (a) Movable tooth structure; (b) model simplification.
Micromachines 10 00159 g003
Figure 4. Effect of the parameters on nonlinear natural frequencies. (a) Harmonic offset a change; (b) movable tooth radius rp changes; (c) harmonic plate radius rs.
Figure 4. Effect of the parameters on nonlinear natural frequencies. (a) Harmonic offset a change; (b) movable tooth radius rp changes; (c) harmonic plate radius rs.
Micromachines 10 00159 g004
Figure 5. Effect of the meshing movable teeth f on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Figure 5. Effect of the meshing movable teeth f on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Micromachines 10 00159 g005
Figure 6. Effect of the harmonic offset a on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Figure 6. Effect of the harmonic offset a on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Micromachines 10 00159 g006
Figure 7. Effect of the movable tooth radius rp on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Figure 7. Effect of the movable tooth radius rp on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Micromachines 10 00159 g007
Figure 8. Effect of the harmonic plate radius rs on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Figure 8. Effect of the harmonic plate radius rs on nonlinear free vibration responses. (a) Δus response; (b) Δuc response; (c) Δur response.
Micromachines 10 00159 g008
Figure 9. Displacement responses of the drive system far from resonance. (a) Δxs response; (b) Δys response; (c) Δus response; (d) Δxc response; (e) Δyc response; (f) Δuc response; (g) Δxr response; (h) Δyr response; (i) Δur response; (j) Δxp1 response; (k) Δyp1 response.
Figure 9. Displacement responses of the drive system far from resonance. (a) Δxs response; (b) Δys response; (c) Δus response; (d) Δxc response; (e) Δyc response; (f) Δuc response; (g) Δxr response; (h) Δyr response; (i) Δur response; (j) Δxp1 response; (k) Δyp1 response.
Micromachines 10 00159 g009aMicromachines 10 00159 g009b
Figure 10. Displacement responses of the drive system close to resonance. (a) Δxs response; (b) Δys response; (c) Δus response; (d) Δxc response; (e) Δyc response; (f) Δuc response; (g) Δxr response; (h) Δyr response; (i) Δur response; (j) Δxp1 response; (k) Δyp1 response.
Figure 10. Displacement responses of the drive system close to resonance. (a) Δxs response; (b) Δys response; (c) Δus response; (d) Δxc response; (e) Δyc response; (f) Δuc response; (g) Δxr response; (h) Δyr response; (i) Δur response; (j) Δxp1 response; (k) Δyp1 response.
Micromachines 10 00159 g010aMicromachines 10 00159 g010b
Figure 11. Finite elements method (FEM) simulation of nonlinear frequencies. (a) Mesh generation; (b) first-order vibration mode.
Figure 11. Finite elements method (FEM) simulation of nonlinear frequencies. (a) Mesh generation; (b) first-order vibration mode.
Micromachines 10 00159 g011
Table 1. Parameters of the drive system.
Table 1. Parameters of the drive system.
ParametersRotorCentral GearHarmonic RodMovable Tooth
mj (Kg)1.31 × 1025.64 × 10−22.59 × 10−23.30 × 10−5
Ij (Kg)9.34 × 10−38.41 × 10−21.30 × 10−21.32 × 10−5
rj (mm)31.633.2292
Table 2. Nonlinear natural frequencies of the drive system.
Table 2. Nonlinear natural frequencies of the drive system.
Number of Meshing TeethFrequencies (rad/s)ω01ω02ω03
16ω0i347,066233,801230,295
ωi349,023249,219191,179
Δωi195715,41839,116
Δωi/ω0i × 1000.566.5916.99
15ω0i346,874235,535229,473
ωi336,116269,738168,614
Δωi10,75834,20260,858
Δωi/ω0i × 1003.1014.5226.52
Table 3. Comparison of nonlinear frequency simulation results.
Table 3. Comparison of nonlinear frequency simulation results.
Orderω1ω2ω3
Theoretical value (rad/s)349,023249,219191,179
Simulation value (rad/s)362,690254,236190,217
Error (%)3.922.010.50

Share and Cite

MDPI and ACS Style

Li, C.; Zhong, W.; Fang, J.; Sun, L. Nonlinear Vibration of a Micro Piezoelectric Precision Drive System. Micromachines 2019, 10, 159. https://doi.org/10.3390/mi10030159

AMA Style

Li C, Zhong W, Fang J, Sun L. Nonlinear Vibration of a Micro Piezoelectric Precision Drive System. Micromachines. 2019; 10(3):159. https://doi.org/10.3390/mi10030159

Chicago/Turabian Style

Li, Chong, Wei Zhong, Jiwen Fang, and Lining Sun. 2019. "Nonlinear Vibration of a Micro Piezoelectric Precision Drive System" Micromachines 10, no. 3: 159. https://doi.org/10.3390/mi10030159

APA Style

Li, C., Zhong, W., Fang, J., & Sun, L. (2019). Nonlinear Vibration of a Micro Piezoelectric Precision Drive System. Micromachines, 10(3), 159. https://doi.org/10.3390/mi10030159

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop