Next Article in Journal
Real-Time Weld Quality Prediction Using a Laser Vision Sensor in a Lap Fillet Joint during Gas Metal Arc Welding
Next Article in Special Issue
Ultra-Sensitive Minute Mass Sensing Using a Microcantilever Virtually Coupled with a Virtual Cantilever
Previous Article in Journal
Challenge Accepted? Individual Performance Gains for Motor Imagery Practice with Humanoid Robotic EEG Neurofeedback
Previous Article in Special Issue
Advances in Sonothrombolysis Techniques Using Piezoelectric Transducers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design and Analysis of Electromagnetic-Piezoelectric Hybrid Driven Three-Degree-of-Freedom Motor

1
School of Electrical Engineering, Hebei University of Science and Technology, Shijiazhuang 050018, China
2
National Engineering Laboratory of Energy-Saving Motor & Control Technique, Anhui University, Hefei 230601, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(6), 1621; https://doi.org/10.3390/s20061621
Submission received: 16 February 2020 / Revised: 12 March 2020 / Accepted: 13 March 2020 / Published: 14 March 2020
(This article belongs to the Special Issue Development of Piezoelectric Sensors and Actuators)

Abstract

:
Multi-DOF movement actuators are widely used in industry, mainly in the fields of bionics and precision machining. With the non-stop improvement of modern-day industry, the requirements for the precision, integration and flexibility of multi-degree-of-freedom motion actuators in the industrial field have progressively increased. This paper presents a novel electromagnetic–piezoelectric hybrid driven three-degree-of-freedom motor. The driving method of the hybrid drive motor can be divided into electromagnetic driving and piezoelectric driving. The motor structure and working principle are analyzed. The structural parameters are obtained by modal analysis of the stators and rotor. The rationality of the stator structure is proved by using the transient analysis of the piezoelectric stators. The magnetic field characteristics of the motor are analyzed by both analytical method and the finite element method. The contact pressure and displacement between the piezoelectric stator and the rotor are analyzed by the analytical method. A motor drive model is established, which provides the basis for motor optimization design and control. Finally, a motor prototype and its test platform were built, and the experimental results are presented to verify the rationality of the motor design.

1. Introduction

A multi-degree-of-freedom (DOF) motor is a device whose rotor can achieve two degrees of freedom or more degrees of freedom of motion. Multi-DOF motors are widely used in bionic fields and multi-degree-of-freedom motion fields. With the further development of modern intelligent industry, the demand for high-precision, high-integration and high-flexibility multi-DOF motors is rapidly growing. In general, a multi-DOF actuator consists of a few traditional motors [1]. However, the mechanism has a complicated structure, a large volume and low positioning control precision, which cannot meet the requirements of the micro-intelligent drive system [2]. Therefore, the multi-DOF integrated actuator drive technology has attracted the attention of experts and scholars.
In the past decade, many multi-DOF motion devices were developed, which have been widely used in bionic robots and multi-DOF motion platforms [3,4,5,6,7]. Mashimo et al. proposed a spherical ultrasonic motor consisting of three annular stators and a spherical rotor. This design enables three-DOF motion, and simplifies the structure of multi-DOF motors [8]. Guo et al. proposed a small three-DOF ultrasonic motor. It is driven by the first longitudinal and second bending vibration modes. The motor performance was improved through a simple piezoelectric lead zirconatetitanate (PZT)-tube stator [9].
Wang, et al., proposed a novel two-DOF spherical ultrasonic motor with three traveling-wave type annular stators to improve the tolerance and controllability of multi-free ultrasonic motors [10]. Chen et al. proposed a novel two-DOF planar linear ultrasonic motor that can achieve the 2-DOF drive function through a variable-mode excitation method [11]. Shi et al. designed a new multi-DOF ultrasonic motor with a ring-shaped composite stator. Large output torque and compact structure are obtained through two orthogonal axial bending modes and a radial bending mode [12]. Liu et al. proposed a sandwich multi-DOF ultrasonic motor with hybrid excitation. The motor is driven through a combination of three resonant modes [13]. The piezoelectric multi-DOF motor motors have the advantages of compact structure, flexible design, easy realization of miniaturization, self-locking, fast response speed, and no external electromagnetic interference [14], so it is suitable for low-speed and high-torque applications. However, it has a short operation life, and it is not suitable for continuous operation. Therefore, Bai et al. proposed a motion platform based on a spherical motor for high-performance nozzle-substrate negotiation in conformal printing of curved electrons. The ball joint type ball motor provides continuous three-DOF rotation, which greatly improves printing efficiency [15]. Li et al. proposed a Halbach array permanent magnet (PM) spherical motor consisting of a spherical rotor and a spherical shell stator. The Halbach array PM structure is used to improve the motor magnetic field distribution and the performance [16]. Lee et al. introduced a spherical multi-DOF motor with double air gaps and proposed a method to improve the motor output characteristics and stability [17]. Li et al. proposed a synthesis strategy for the stator magnetic field of a PM spherical motor to support the further control of the motor. The electromagnetically-driven multi-DOF motor has the characteristics of high mechanical integration and high speed, which is suitable for continuous operation [18]. However, the structure of the electromagnetically-driven multi-DOF motor is relatively complicated, and its magnetic field is difficult to control [19].
The development of hybrid drive multi-DOF motors has gradually become a research trend [20]. The hybrid driven multi-DOF motor mainly combines various drive modes, and hence it does not have the problems of the conventional one drive mode. For example, electromagnetic hybrid-driven multi-DOF motors have problems, such as their difficulty to control the magnetic field, low positioning accuracy and difficulty to transmit high torque through output shaft [21,22], and the multi-DOF motor with piezoelectric hybrid driven has the problems of low operating life and unstable performance [23]. In order to compensate the weaknesses of each type, a novel electromagnetic-piezoelectric hybrid-driven three-DOF is proposed. This hybrid-driven three-DOF motor combines two different driving types. The coarse movement of the motor over a large range is driven by electromagnetic drive. The small movement of the motor is driven by the piezoelectric drive. The electromagnetic-piezo-hybrid driving method not only improves the accuracy of the multi-degree-of-freedom motor, but also realizes the continuous and stable operation of the motor.
In this paper, the basic structure, working principle and hybrid drive strategy of the motor are described. Then, the motor structural model is established by the finite element analysis, and the modal analysis, harmonic response analysis and transient analysis of the motor are carried out. The air gap magnetic field characteristics of the motor are analyzed by analytical method, and the rationality of the results is verified by the finite element method. The contact pressure and contact displacement between the piezoelectric stator and the rotor are analyzed through the analytical method, and the motor drive model is established. Finally, a motor prototype and its motor test platform are built. The design rationality of the motor is proven by the experimental results on the prototype.

2. Structure and Principle of the Motor

As shown in Figure 1, the electromagnetic-piezoelectric hybrid drive three-DOF motor is mainly composed of the motor base, pre-pressure regulating system, PMs, electromagnetic stator core, electromagnetic stator coils, piezoelectric stators (No.1 stator, No.2 stator and No.3 stator) and rotor. The number of permanent magnets is four, two of which are N poles and two are S poles, which are alternately arranged.
The permanent magnet is a spherical structure, which is attached to the rotor. The angle (α) between the central axis of the three piezoelectric stators and the horizontal plane is 17.596°. The included angle of the central axis of each stator is 120°. There are four spherical permanent magnets attached to the rotor. The pre-pressure between the stator and the rotor is adjusted by a pre-pressure regulating system.
As shown in Figure 2a, the pre-pressure regulating system is mainly composed of a pre-pressure base, a pre-pressure bracket, a pre-pressure rod, a pre-pressure sleeve and a pre-pressure slider. The pre-pressure sleeve and the motor base are connected by a screw, and the bottom of the pre-pressure rod and the pre-pressure sleeve are connected by another screw. In Figure 2b, the force generated by rotating the pre-pressure lever is converted into the pressure of the piezoelectric stator on the rotor spherical shell. When the pre-pressure bar moves upward, the pre-pressure slider moves in a direction close to the rotor case. In this way, the pre-pressure between the piezoelectric stators and rotor is applied.
As shown in Figure 3, the piezoelectric stator is composed of the piezoelectric base body and the piezoelectric ceramic sheets. Alternating voltages with the same frequency and amplitude, but a phase difference of 90° is applied to the two phases of the piezoelectric ceramic sheets.
The piezoelectric stator responds to a pulsating wave with a phase difference of 90° in both time and space, which is in accordance with the principle of a piezoelectric inverse effect. A circularly traveling wave is formed by the superposition of the two-phase standing waves. Then the rotor rotates along the tangential direction of the particle surface of the piezoelectric stator surface, and the direction of rotor rotation is opposite to the direction of the traveling wave.
As shown in Figure 4a,b, the line passing through the No. 2 stator (S2) is defined as the Y-axis, the axis perpendicular to the Y-axis is the X-axis, and the line perpendicular to the X-Y plane and passing through the origin is defined as the Z-axis. The angle between the piezoelectric stators and the X-Y plane is α°.
As shown in Figure 4c, the electromagnetically-driven stator is mounted outside the spherical rotor. The rotor torque is generated by the interaction between the energized coil and the PM magnetic field. The 3-DOF motion is achieved by the synthesis of four stators motions. ωx, ωy, and ωz are defined as the angular velocity vectors of the rotor rotating around the X, Y, and Z axes, respectively. According to the motor structure, the relation among ωx, ωy, ωz] and the angular velocity vectors ω1, ω2, ω3, ω4 of the four stators can be expressed as:
ω = [ ω x ω y ω z ] = [ 3 2 ω 1 cos α 3 2 ω 3 cos α 1 2 ω 1 cos α + ω 2 cos α 1 2 ω 3 cos α ω 1 sin α + ω 2 sin α + ω 3 sin α + ω 4 ]

3. Motor Design and Analysis

3.1. Modal Analysis

Modal analysis is an important issue for the motor dynamic properties. The modal analysis can predict and optimize the structural dynamics in product design.
Through modal analysis, the pattern shape and the motor structure in the working frequency domain are identified. The basic structure of the motor is then determined.
The COMSOL finite element model of the motor is established in this paper. The structural parameters of the motor are as follows:
The piezoelectric stator is formed by bonding the elastic metal body and the piezoelectric ceramic. The number of teeth in the piezoelectric stator is 72. The stator materials are listed in Table 1. The elastic metal body (piezoelectric base body) material is phosphor bronze and the piezoelectric material is PZT-8.
The outer radius of the piezoelectric stators is 30 mm, the tooth height of the piezoelectric stators is 2 mm, the thickness of the piezoelectric ceramics is 0.5 mm and the thickness of the intermediate elastic body is 2.5 mm. The model of a piezoelectric stator is shown in Figure 3. Each tooth occupies a tooth height of 2 mm, a groove angle of 5° between the teeth, a circumferential angle of 3°, a tooth gap width of 5.5 mm, an elastic base thickness of 2.5 mm, and a piezoelectric ceramic thickness of 0.5 mm. The outer radius of the electromagnetic stator is 93.6 mm. The air gap between the electromagnetic stator and rotor is 1.5 mm. The inner radius of the rotor is 49 mm. The PM thickness is 10 mm. In order to predict the dynamic characteristics of the motor structure, the modal analysis of the stators and rotor is carried out.
Based on the finite element model of the motor, the vibration frequency of the piezoelectric stator is obtained by modal analysis of a single stator. The specific steps are as follows: The solid mechanics module is used in the COMSOL finite element analysis software, and the material of the selected piezoelectric stator is shown in Table 1. The inner ring of the piezoelectric stator is set fixed, and the frequency analysis range is set to 20–60 kHz. The mode shape and frequency of the 9th order bending vibration mode of the piezoelectric stator are obtained through modal analysis. Figure 5 is a 9th-order modal shape diagram of a piezoelectric stator.
As shown in Figure 5, the modal frequency of the 9th-order mode of the motor is 40–41 kHz. When the motor rotor and the electromagnetic stator have a modal at 40–41 kHz, the motor rotor and the electromagnetic stator resonate, which causes the structural deformation of the rotor and reduces the performance of the motor. Therefore, in order to avoid the above problems, the solid mechanics module in the COMSOL finite element analysis software is used to perform modal analysis on the motor rotor and the electromagnetic stator. The specific steps are: the definition of the material parameters is shown in the table, and the boundary conditions are set to free boundary conditions. The frequency range is defined as 40–41 kHz. The outer diameter variable of the rotor is defined as wqR. The scanning parameter values of wqR are 51.5 mm, 52.5 mm, 53.5 mm, 54.5 mm, 55.5 mm.
As shown in Figure 6, when the outer radius is 53.5 mm, the rotor structure has two resonant modal modes with frequencies of 39.987 kHz and 41.72 kHz. The motor does not have any natural (Eigen) frequency between 40–41 kHz when the wqR is 53.5 mm. Then it can be concluded that the working frequency of the motor is not close to the resonant frequency.
As shown in Figure 7, the electromagnetic stator has no vibration mode in the frequency range of 40–41 kHz. Therefore, the electromagnetic stator does not resonate, and the structure of the rotor will not be deformed.

3.2. Frequency Domain Analysis

The frequency domain analysis of the piezoelectric stator determines the steady-state response when it is excited by a sinusoidal voltage. The impedance and frequency variation and the electrical resonance frequency of the piezoelectric stator can also be obtained.
By using the solid mechanical analysis module in the COMSOL finite element analysis software, a voltage load is applied to the piezoelectric ceramics of the stator, and the form of the load is a plural form. The magnitude of the voltage on one side of the piezoelectric substrate is 0 V, and the magnitude of the voltage on the other side is 200 V. The layout of the applied load is shown in Figure 8a. The excitation voltage applied by phase A is 200 × exp(π/2 × i), and the excitation voltage applied by phase B is 200 × exp(0 × i). From the modal analysis results in Figure 5, the piezoelectric stator frequency is about at 40 kHz, so the frequency domain analysis is performed with 39.5–41.5 kHz. The frequency domain analysis result of the piezoelectric stator obtained by post-processing is shown in Figure 8.
As shown in Figure 8b, the piezoelectric stator reaches the resonance frequency at 40.365 kHz, and the frequency corresponding to the maximum operating speed of the motor is obtained.

3.3. Transient Analysis

The transient analysis of the piezoelectric stator utilizes a state equation that describes the continuous state change of the motor. The trajectory and vibration displacement of the piezoelectric stator are investigated by finite element analysis. Based on the frequency domain analysis, the phase A voltage excitation is defined as 200 × sin(2π × f), the phase B voltage excitation is defined as 200 × cos(2π × f), and f is 40.365 kHz. The material damping [24] is design by the damping ratio and frequency for modal analysis. The material parameters [25] of the piezoelectric ceramic sheet are:
The ceramic piezoelectric constant matrix is
[ e ] = [ 0 0 3.875 0 0 3.875 0 0 13.91 0 0 0 0 10.344 0 10.344 0 0 ]   C / m 2 ,
the dielectric constant matrix is
[ ε s   ] = [ 8.002 0 0 0 8.002 0 0 0 4.968 ] × 10 9 F / m ,
and the piezoelectric ceramic short-circuit stiffness matrix is
[ c E ] = [ 14.7 8.11 8.11 0 0 0 8.11 13.2 7.3 0 0 0 8.11 8.11 13.2 0 0 0 0 0 0 3.29 0 0 0 0 0 0 3.13 0 0 0 0 0 0 3.13 ] × 10 10 Pa .
The movement of the motor stator at 0–5 ms is shown in the Figure 9. Four points in the piezoelectric stator drive boundary are selected, which is used to analyze the stator vibration characteristics, as shown in Figure 9a. It can be seen from Figure 9b–d that the trajectory of piezoelectric stator particle displacement is not elliptical at the beginning of the applied excitation. The displacement of the piezoelectric stator particle on the X-Z and Y-Z planes gradually increases, eventually reaching the maximum. Currently, the motion track of the particle becomes an elliptical.
It can be seen from Figure 10 that the displacement in the X, Y or Z direction is 0 µm at the time of t = 0. The displacement of the piezoelectric stator in the X, Y or Z direction reaches the maximum value within 0.1 ms, and the maximum value does not change. In other words, the time for the motor to start is 0.1 ms. The maximum displacement of the piezoelectric stator surface mass in the Z direction is about 0.5263 µm, and the results show that the displacement of the piezoelectric stator in the X and Y directions is about 0.1831 µm and 0.2443 µm, respectively. In order to generate a traveling wave that drives the rotor to rotate, the piezoelectric ceramic sheet is polarized along the thickness direction, that is, the Z direction, and the X and Y directions are not polarized. Therefore, when the piezoelectric ceramic sheet is excited, the Z-direction displacement is greater than the X- and Y-directions due to the inverse piezoelectric effect principle. The piezoelectric stator particle displacement in the X or Y direction is very small, so side shifting does not occur during the piezoelectric stator operation. It can be concluded that the piezoelectric stator structure is reasonable.

3.4. Piezoelectric Stator and Rotor Contact Analysis

In order to simplify calculation, the three-dimensional contact between the piezoelectric stator and the rotor is reduced to a one-dimensional problem (as shown in Figure 11). A wave contact between the piezoelectric stator and the rotor is analyzed.
The contact between the piezoelectric stator and the rotor is equivalent to the contact of the parabolic indenter with the friction contact layer, as shown in Figure 11. The piezoelectric stator friction contact layer is set to be displaced in the Z-axis direction to the maximum displacement when the piezoelectric stator is in contact with the rotor. The piezoelectric stator is a rigid body and does not displace, so the friction between the piezoelectric stator and rotor is ignored.
The waveform of the piezoelectric stator at a certain moment is:
u 1 = Acos ( k x ω t )
The curvature radius Rc of this waveform is [26]:
R c = | [ 1 + ( d u 1 d x ) 2 ] 3 2 / d 2 u 1 d x 2 |
The pre-pressure applied by each stator is F, and then the one-dimensional pre-stress experienced by each peak is the fixed-rotor contact zone displacement function, which can be defined as:
w = x 2 2 R c + d 0   , z = 0 , | x | a
where, d0 is an unknown constant and, a is half the contact distance of the fixed rotor within one wavelength of the piezoelectric stator.
The friction between the piezoelectric stator and the rotor is neglected, and the shear force at the piezoelectric stator and rotor boundary is zero. That is x = 0, τxy = 0.
The piezoelectric stator and rotor elastic contact normal stress distribution is p(x), z = 0, |x| ≤ a.
According to the balance of the piezoelectric stator and rotor:
a a p ( x ) d x = F n b
where, n is the number of traveling waves, and b is the equivalent width of the piezoelectric stator contact.
When the fixed rotor contact force p(ξ) acts, the vertical displacement expression [26] is:
1 π E a a ln ( x ξ ) 2 p ( ξ ) d ξ = x 2 2 R c + d 0
where, E is the Young’s modulus of the friction layer.
By solving the first type of Fredholm equation and solving the boundary conditions [27,28,29], the rotor contact force distribution is as follows:
p ( x ) = ( P E π R c a 2 ) 1 2 ( a 2 x 2 ) 1 2
The piezoelectric stator-rotor contact surface displacement is:
u = x 2 2 R c + P π E ( 2 ln 2 + 1 )
The parameters of the piezoelectric stator are as shown in Table 2.
According to the data in Table 2, from Equations (7) and (8), we use MATLAB software to numerically simulate the contact pressure and contact surface distribution of the stator and rotor of the motor. The piezoelectric stator-rotor contact displacements calculated with the parameters are shown in Figure 12a–c. It can be seen from Figure 12a–c that the greater the pre-pressure, the greater the length of the rotor contact. The greater the pre-pressure, the greater the displacement of the pressure contact point. The greater the pre-pressure, the greater the pressure at the contact point.

3.5. Modeling and Analysis of Rotor Air Gap Magnetic Field

The rotor PM structure is shown in Figure 13. According to the magnetic medium boundary conditions, the magnetic field generated by the rotor permanent magnet can be divided into two parts:
{ B 1 = μ 0 μ r H 1 r R o B 2 = μ 0 μ m H 2 + μ 0 M r R i r < R o
where, μ0 is the magnetic permeability in vacuum, its value is 4π × 10−7 H/m, μm is the relative magnetic permeability of the permanent magnet, and its value is close to 1, and μr is the relative permeability of air, its value is also close to 1.
The PM material is NdFeB35, Mr is the residual magnetization. |Mr| = Br/μ0 (A/m), Br is the residual magnetic flux density with a value of 1.17 T.
The parallel permanent magnet magnetization vector is:
M r = { M r r M r θ M r ϕ } = ( 1 ) p 1 | M r | { cos ( φ α p ) sin θ cos ( φ α p ) cos θ sin ( φ α p ) }
α p = α r / 2 + 2 π ( p 1 ) / L
π / 2 β r / 2 θ π / 2 + β r / 2
π ( p 1 ) / L φ α r + π ( p 1 ) / L
where p is the permanent magnet number (p = 1, 2, 3, 4), L is the total number of permanent magnets, αr is the permanent magnet with a longitude angle of 41° and βr is the permanent magnet with a latitude angle of 85°.
The PM mounting positions are symmetrically distributed, hence
M r = 0
There is no current excitation in the rotor area, and we have:
× H i = 0
H i = Φ i
where, i = 1, 2.
The magnetic scalar potentials in the two regions can be expressed as:
2 Φ 1 = 0
2 Φ 2 = 0
The Laplace equation of scalar magnetic potentials Φi can be expressed as:
1 r 2 r ( r 2 Φ i r ) + 1 r 2 sin θ θ ( sin θ Φ i θ ) + 1 r 2 sin 2 θ 2 Φ i φ 2 = 0
The solution of the magnetic potential equation by using the variable separation method is obtained as:
Φ i ( r , θ , φ ) = R ( r ) Y ( θ , φ ) = n = 0 + m = 0 n ( C n r n + D n r n + 1 ) P n m ( cos θ ) ( A m cos m φ + B m sin m φ ) = n = 0 + m = n n ( C n r n + D n r n + 1 ) Y n m ( θ , φ )
where, i = 1,2. Y n m ( θ , φ ) = S n m P n m ( cos θ ) e i m φ (m = 0, 1, … n), S n m = 2 n + 1 4 π ( n m ) ! ( n m ) ! , Y n m ( θ , φ ) is a spherical harmonic function, P n m ( cos θ ) is the associated Legendre function.
The boundary conditions for the solution are:
Φ 1 | r = 0
B 1 r = B 2 r | r = r o
H 1 θ = H 2 θ | r = r o
H 1 ϕ = H 2 ϕ | r = r o
The magnetic field strength in spherical coordinates is:
H i = [ H r e r , H θ e θ , H φ e φ ] T = [ Φ i r , 1 r Φ i θ , 1 r sin θ Φ i φ ] T
where, i = 1,2.
The fundamental component of flux density can be obtained by solving the boundary conditions as
B 11 = [ B 11 r B 11 θ B 11 φ ] = μ 0 μ r r 4 D n , 1 m S n m [ 3 ( sin θ ) 2 sin 2 φ 2 sin θ cos θ sin 2 φ 2 sin θ cos 2 φ ]
In order to verify the accuracy of the magnetic field analysis, the air gap magnetic field at the distance of 0.5 mm from the permanent magnet is taken, and the comparative analysis is carried out by the finite element method (FEM). The results are shown in Figure 14. The results of the analytical method and the FEM are in relative agreement. The feasibility of the analytical method is verified, which provides a basis for the fast analysis and optimization of the motor prototype.

4. Driver Model of Electromagnetic-Piezoelectric-Hybrid Drive Motor

The mathematical driving model of the electromagnetic-piezoelectric-hybrid drive three-degree-of-freedom motor is composed of two parts: the electromagnetic drive and the piezoelectric drive.

4.1. Piezoelectric Drive Modeling

According to the Mindlin theory [30,31,32,33], the tangential force distribution of the stator and rotor is:
When Tmax μ ( PE π R ) , the tangential force distribution is provided by the following equation:
{ f ( x ) = μ ( PE π R ) 1 2 ( 1 x 2 a 2 ) 1 2 c | x | a f ( x ) = μ ( PE π R ) 1 2 [ ( 1 x 2 a 2 ) 1 2 2 c a ( 1 x 2 c 2 ) ] c max | x | a f ( x ) = μ ( PE π R ) 1 2 [ ( 1 x 2 a 2 ) 1 2 2 c a ( 1 x 2 c 2 ) + c max a ( 1 x 2 c max 2 ) ] | x | c max
where cmax is c max a = ( 1 T max μ PE π R ) 1 2 , c a = ( 1 T max T μ PE π R ) 1 2 , μ is the dynamic friction, T is the resistance of the motor movement, Tmax is the maximum resistance of the motor operation and c is a function about T.
When Tmax = μ ( PE π R ) , the tangential force distribution is provided by the following equation:
{ f ( x ) = μ ( PE π R ) 1 2 ( 1 x 2 a 2 ) 1 2 c | x | a f ( x ) = μ ( PE π R ) 1 2 [ ( 1 x 2 a 2 ) 1 2 2 c a ( 1 x 2 c 2 ) ] c max | x | a
where, c a = ( 1 T max T μ PE π R ) 1 2 .
Under the normal contact pressure in the contact area, the tangential thrust Ft on the resultant rotor is:
F t = a a f ( x ) d x
The output torque is
T j = n R F t
where n is the number of traveling wave peaks, R is the outer radius of the piezoelectric stators, and j (j = 1, 2, 3) is the number of the piezoelectric stators.

4.2. Electromagnetic Drive Modeling

The magnetic field generated by the energized coil interacts with the magnetic field generated by the PMs, producing an interaction Lorentz force. The electromagnetic force generated by a coil dl micro-element can be obtained by the Lorentz force law. Then, the electromagnetic torque generated by the dl micro-element is calculated, and the torque micro-element is integrated to obtain the torque generated by the single coil. Finally, the rotation torque of the electromagnetic stator drive is obtained [29]:
T 4 = i = 1 N J k r 1 r 2 δ 0 δ 1 0 2 π B 11 r ( r , θ i , φ i , ϕ , δ ) sin δ sin ϕ d ϕ d r d θ
where k is the coil number, and N is the total number of coils.

4.3. Electromagnetic-Piezoelectric Hybrid Motor Modeling

As shown in Figure 4, the line passing through the No. 2 stator (S2) is defined as the Y-axis, the axis perpendicular to the Y-axis is the X-axis, and the line perpendicular to the X-Y plane and passing through the origin is defined as the Z axis, and the angle between the piezoelectric stators and the X-Y plane is α°, the torque of the three-degree-of-freedom motor is T, which is synthesized by the torque vectors T1, T2, T3 and T4 of the four stators:
T = [ T x T y T z ] = [ 3 2 T 1 cos α 3 2 T 3 cos α 1 2 T 1 cos α + T 2 cos α 1 2 T 3 cos α T 1 sin α + T 2 sin α + T 3 sin α + T 4 ]

5. Experimental Verification and Analysis

In order to verify the rationality of the motor structure and performance analysis, a prototype experimental test system is built as shown in Figure 15. The system includes a 3-DOF prototype of an electromagnetic-piezoelectric-hybrid drive motor, a driver and a host computer. The driver used in this experiment is UMD-3II. The power line of the DC12V power supply is connected to the UMD-3II driver power port, the A phase of the ceramic chip is connected to the sine port of the driver, and the cosine port of the B phase driver of the ceramic chip. The UMD-3II driver uses RS232 serial communication. The serial communication line of the driver is connected to the host computer. Adjust the frequency of the excitation voltage to measure the speed of the motor on the host computer. The rotation speed measured by this platform is the rotation speed (n1, n2, n3) around the central axis of the stator. nx, ny, can be calculated by Equation (33).
n = [ n x n y n z ] = [ 3 2 n 1 cos α 3 2 n 3 cos α 1 2 n 1 cos α + n 2 cos α 1 2 n 3 cos α n 1 sin α + n 2 sin α + n 3 sin α + n 4 ]
Figure 16 is a graph of the rotation speed n1 with frequency. The FEM rotation speed n1 is obtained according to a previous method [34]. The experimental and FEM methods are in good agreement. It can be seen from Figure 16 that the rotation speed of the motor decreases as the frequency increases.
According to Equation (33), the electromagnetic waves of No. 1 and No. 3 stators (S1 and S3) are driven in the forward direction, and that of No. 2 stator (S2) is reversely driven to obtain the rotation speed of the motor at different voltage excitation frequencies, as shown in Figure 17. As can be seen from Figure 17, the maximum speed of motor is about 40.2 kHz, which is the approximate resonant frequency of a piezoelectric stator. It can be known from the experimental results that the speed of the motor decreases as the frequency increases, so the speed can be adjusted by changing the voltage excitation frequency.

6. Conclusions

A novel electromagnetic-piezoelectric-hybrid drive three-degree-of-freedom motor was proposed. The structural model was established by finite element analysis software, and the modal analysis, harmonic response analysis and transient analysis of the motor were conducted. The magnetic field characteristics are analyzed and compared by analytical and finite element methods. The contact pressure and contact displacement between the piezoelectric stator and the rotor were analyzed by the analytical method, and the motor drive model was established. The rationality and feasibility of the design are verified by the experimental results on a prototype. The technical scheme of the electromagnetic-piezoelectric-hybrid three-degree-of-freedom motor is feasible, and it has research value. In the next step, it is planned to investigate further the electromagnetic-piezoelectric-hybrid drive mechanism for more effective modeling and system-level multi-discipline optimization for high-performance three-DOF motor drives.

Author Contributions

Writing—review and editing, supervision, project administration, funding acquisition, Z.L.; methodology, visualization, MATLAB software, writing—original draft preparation, P.G.; COMSOL software, data curation, Z.W.; ANSOFT software, L.Z. formal analysis, Q.W. All authors have read and agree to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (51877070, 51577048, 51637001); Natural Science Foundation of Hebei Province of China (E2018208155); Hebei Province Higher Education Science and Technology Research Key Project (ZD2018228); the Talent Engineering Training Support Project of Hebei Province, grant No. A201905008; High-energy-saving Motor and Control Technology National and Local Joint Engineering Laboratory Open Project Funded Project (KFKT201901), Hebei Province Graduate Innovation Funding Project (CXZZSS2018085, CXZZSS2019084).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, X.; Li, S.; Wang, Q.; Wen, Y.; Gong, N. Dynamic analysis and current calculation of a permanent magnet spherical motor for point-to-point motion. IET Electr. Power Appl. 2019, 13, 426–434. [Google Scholar] [CrossRef]
  2. Xie, L.; Si, J.; Hu, Y.; Wang, Z. Overview of 2-degree-of-freedom rotary-linear motors focusing on coupling effect. IEEE Trans. Magn. 2019, 55, 1–11. [Google Scholar] [CrossRef]
  3. Nakamura, K.; Zhang, X.; Ueha, S. A multi-degrees-of-freedom ultrasonic motor design for robotics applications. In Proceedings of the 2006 IEEE Ultrasonics Symposium, Vancouver, BC, Canada, 2–6 October 2006; pp. 2281–2284. [Google Scholar]
  4. Bansevicius, R.; Blechertas, V. Multi-degree-of-freedom ultrasonic motors for mass-consumer devices. J. Electroceramics 2008, 20, 221–224. [Google Scholar] [CrossRef]
  5. Hoshina, M.; Toyama, S. Development of spherical ultrasonic motor as a camera actuator for pipe inspection robot. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, St. Louis, MO, USA, 10–15 October 2009; pp. 2379–2384. [Google Scholar]
  6. Shen, S.C.; Tsai, P.C.; Wang, Y.J.; Huang, H.J. A new type of multi-degree-of-freedom miniaturization actuator using symmetric piezoelectric pusher element for a pocket sun-tracking system. Sens. Actuators A Phys. 2012, 182, 114–121. [Google Scholar] [CrossRef]
  7. Hoshina, M.; Mashimo, T.; Fukaya, N.; Matsubara, O.; Toyama, S. Spherical ultrasonic motor drive system for camera orientation in pipe inspection. Adv. Robot. 2013, 27, 199–209. [Google Scholar] [CrossRef]
  8. Mashimo, T.; Toyama, S.; Ishida, H. Design and implementation of spherical ultrasonic motor. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2009, 56, 2514–2521. [Google Scholar] [CrossRef] [PubMed]
  9. Guo, M.; Hu, J.; Zhu, H.; Zhao, C.; Dong, S. Three-degree-of-freedom ultrasonic motor using a 5-mm-diameter piezoelectric ceramic tube. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2013, 60, 1446–1452. [Google Scholar]
  10. Wang, J.; Hu, X.; Wang, B.; Guo, J. A novel two-degree-of-freedom spherical ultrasonic motor using three travelling-wave type annular stators. J. Cent. South Univ. 2015, 22, 1298–1306. [Google Scholar] [CrossRef]
  11. Zhou, X.; Chen, W.; Liu, J. Novel 2-dof planar ultrasonic motor with characteristic of variable mode excitation. IEEE Trans. Ind. Electron. 2016, 63, 6941–6948. [Google Scholar] [CrossRef]
  12. Shi, S.; Xiong, H.; Liu, Y.; Chen, W.; Liu, J. A ring-type multi-DoF ultrasonic motor with four feet driving consistently. Ultrasonics 2017, 76, 234–244. [Google Scholar] [CrossRef]
  13. Yang, X.; Liu, Y.; Chen, W.; Liu, J. Sandwich-type multi-degree-of-freedom ultrasonic motor with hybrid excitation. IEEE Access 2017, 4, 905–913. [Google Scholar] [CrossRef]
  14. Liu, Y.; Deng, J.; Su, Q. Review on multi-degree-of-freedom piezoelectric motion stage. IEEE Access 2018, 6, 59986–60004. [Google Scholar] [CrossRef]
  15. Bai, K.; Xu, R.; Lee, K.-M.; Dai, W.; Huang, Y. Design and development of a spherical motor for conformal printing of curved electronics. IEEE Trans. Ind. Electron. 2018, 65, 9190–9200. [Google Scholar] [CrossRef]
  16. Li, H.; Li, T. End-effect magnetic field analysis of the halbach array permanent magnet spherical motor. IEEE Trans. Magn. 2018, 54, 1–9. [Google Scholar] [CrossRef]
  17. Cho, S.; Lim, J.-S.; Oh, Y.J.; Jeong, G.; Kang, D.-W.; Lee, J. A study on output characteristics of the spherical multi-DoF motor according to the number of phases and pole pitch angles. IEEE Trans. Magn. 2018, 54, 1–5. [Google Scholar] [CrossRef]
  18. Li, B.; Zhang, S.; Li, G.; Li, H. Synthesis strategy for stator magnetic field of permanent magnet spherical motor. IEEE Trans. Magn. 2018, 54, 1–5. [Google Scholar] [CrossRef]
  19. Li, Z.; Wang, Y.; Ge, R.; Zhao, G. The summary and latest research of PM spherical m-DoF motor. Micromotors 2011, 44, 66–70. [Google Scholar]
  20. Zhang, J.; Yuan, L.; Chen, S.-L.; Liang, Y.; Huang, X.; Zhang, C.; Yang, G. A survey on design of reaction spheres and associated speed and orientation measurement technologies. ISA Trans. 2019, 1–14. [Google Scholar] [CrossRef]
  21. Li, Z.; He, W.; Liu, L.; Wang, Q. Calculation and optimization of a permanent magnetic hybrid driven 3-dof motor. Int. J. Appl. Electromagn. Mech. 2018, 58, 1–33. [Google Scholar]
  22. Li, Z.; He, W.; Liu, L.; Wang, Q. Magnetic field calculation and dynamics simulation of a permanent magnetic hybrid driven 3-dof motor. J. Electr. Eng. Technol. 2019, 14, 711–718. [Google Scholar] [CrossRef]
  23. Liu, Y.; Yan, J.; Wang, L.; Chen, W. A two- DoF ultrasonic motor using a longitudinal–bending hybrid sandwich transducer. IEEE Trans. Ind. Electron. 2019, 66, 3041–3050. [Google Scholar] [CrossRef]
  24. Pons, J.L.; Rodríguez, H.; Rocon, E.; Fernández, J.F.; Villegas, M. Practical consideration of shear strain correction factor and rayleigh damping in models of piezoelectric transducers. Sens. Actuators A Phys. 2004, 115, 202–208. [Google Scholar] [CrossRef]
  25. Sheldon, I. Conversion of Piezoelectric Material Data; Collaborative Solutions Inc.: Redondo Beach, CA, USA, 1999. [Google Scholar]
  26. Kong, X. Solid Contact Mechanics; China Railway Publishing House: Beijing, China, 1999. [Google Scholar]
  27. Aizikovich, S.M.; Vasil’ev, A.S.; Volkov, S.S. The axisymmetric contact problem of the indentation of a conical punch into a half-space with a coating inhomogeneous in depth. J. Appl. Math. Mech. 2015, 79, 500–505. [Google Scholar] [CrossRef]
  28. Vasiliev, A.S.; Volkov, S.S.; Aizikovich, S.M.; Mitrin, B.I. Plane contact problem on indentation of a flat punch into a transversely-isotropic half-plane with functionally graded transversely-isotropic coating. Z. Angew. Math. Phys. 2016, 68, 4. [Google Scholar] [CrossRef]
  29. Li, Z.; Guo, P.; Han, R.; Wang, Q. Torque modeling and characteristic analysis of electromagnetic piezoelectric hybrid-driven 3-degree-of-freedom motor. Adv. Mech. Eng. 2018, 10, 1–15. [Google Scholar] [CrossRef] [Green Version]
  30. Lucas, B.N.; Hay, J.C.; Oliver, W.C. Using multidimensional contact mechanics experiments to measure poisson’s ratio. MRS Proc. 2004, 19, 58–65. [Google Scholar] [CrossRef]
  31. Lei, G.; Zhu, J.; Guo, Y.; Liu, C.; Ma, B. A review of design optimization methods for electrical machines. Energies 2017, 10, 1962. [Google Scholar] [CrossRef] [Green Version]
  32. Lei, G.; Wang, T.; Zhu, J.; Guo, Y.; Wang, S. System-level design optimization method for electrical drive systems—Robust approach. IEEE Trans. Ind. Electron. 2015, 62, 4702–4713. [Google Scholar] [CrossRef]
  33. Zhao, D.; Li, Y.; Jernej, B. 6-DoF haptic rendering of static coulomb friction using linear programming. IEEE Trans. Haptics 2018, 11, 325–337. [Google Scholar] [CrossRef]
  34. Hatsuzawa, T.; Toyoda, K.; Tanimura, Y. Speed control characteristics and digital servosystem of a circular traveling wave motor. Rev. Sci. Instrum. 1986, 57, 2886–2890. [Google Scholar] [CrossRef]
Figure 1. Structure of electromagnetic-piezoelectric hybrid drive motor.
Figure 1. Structure of electromagnetic-piezoelectric hybrid drive motor.
Sensors 20 01621 g001
Figure 2. Pre-pressure regulating system: (a) Schematic diagram of pre-pressure system (b) Pre-pressure regulating diagram.
Figure 2. Pre-pressure regulating system: (a) Schematic diagram of pre-pressure system (b) Pre-pressure regulating diagram.
Sensors 20 01621 g002
Figure 3. Piezoelectric stator.
Figure 3. Piezoelectric stator.
Sensors 20 01621 g003
Figure 4. Motor structure: (a) Top view of piezoelectric stators/rotor structure, (b) Front view of piezoelectric stators/rotor structure, and (c) Structure of electromagnetic stator/rotor.
Figure 4. Motor structure: (a) Top view of piezoelectric stators/rotor structure, (b) Front view of piezoelectric stators/rotor structure, and (c) Structure of electromagnetic stator/rotor.
Sensors 20 01621 g004
Figure 5. 9th-order mode shape diagram of piezoelectric stator.
Figure 5. 9th-order mode shape diagram of piezoelectric stator.
Sensors 20 01621 g005
Figure 6. Rotor mode shape diagram.
Figure 6. Rotor mode shape diagram.
Sensors 20 01621 g006
Figure 7. Electromagnetic stator mode shape diagram.
Figure 7. Electromagnetic stator mode shape diagram.
Sensors 20 01621 g007
Figure 8. Motor frequency domain analysis chart: (a) layout of the applied load, (b) Motor frequency impedance diagram.
Figure 8. Motor frequency domain analysis chart: (a) layout of the applied load, (b) Motor frequency impedance diagram.
Sensors 20 01621 g008
Figure 9. Trajectory diagram of the motor: (a) Schematic diagram of the point on the piezoelectric stator, (b) Trajectory on X-Y plane, (c) Trajectory on X-Z plane, and (d) Trajectory on Y-Z plane.
Figure 9. Trajectory diagram of the motor: (a) Schematic diagram of the point on the piezoelectric stator, (b) Trajectory on X-Y plane, (c) Trajectory on X-Z plane, and (d) Trajectory on Y-Z plane.
Sensors 20 01621 g009
Figure 10. Piezoelectric stator time-displacement diagram.
Figure 10. Piezoelectric stator time-displacement diagram.
Sensors 20 01621 g010
Figure 11. Piezoelectric Stator and rotor contact diagram.
Figure 11. Piezoelectric Stator and rotor contact diagram.
Sensors 20 01621 g011
Figure 12. Piezoelectric stator and rotor contact: (a) Pressure vs. displacement of the contact under different pre-pressures, (b) Displacement of contact surface vs. contact length under different pre-pressures, and (c) Displacement of contact vs. pre-pressure.
Figure 12. Piezoelectric stator and rotor contact: (a) Pressure vs. displacement of the contact under different pre-pressures, (b) Displacement of contact surface vs. contact length under different pre-pressures, and (c) Displacement of contact vs. pre-pressure.
Sensors 20 01621 g012
Figure 13. Permanent magnet structure diagram of the rotor.
Figure 13. Permanent magnet structure diagram of the rotor.
Sensors 20 01621 g013
Figure 14. Finite element and analytical comparison analysis of air gap magnetic field.
Figure 14. Finite element and analytical comparison analysis of air gap magnetic field.
Sensors 20 01621 g014
Figure 15. Prototype experimental test system.
Figure 15. Prototype experimental test system.
Sensors 20 01621 g015
Figure 16. Speed vs. frequency diagram of piezoelectric stator.
Figure 16. Speed vs. frequency diagram of piezoelectric stator.
Sensors 20 01621 g016
Figure 17. Synthetic speed vs. frequency diagram.
Figure 17. Synthetic speed vs. frequency diagram.
Sensors 20 01621 g017
Table 1. Main Parameters of Stators and Rotor Materials.
Table 1. Main Parameters of Stators and Rotor Materials.
Motor ComponentMaterialYoung’s Modulus (N/m2)Poisson’s RatioDensity (kg/m3)
Piezoelectric ceramicsPZT-8112 × 1090.337600
Piezoelectric base bodyPhosphor bronze8 × 1090.38800
Rotor materialHard aluminum72 × 1090.332700
Rotor PMNdFeB3516 × 10100.37500
Electromagnetic stator coilCopper11 × 10100.358960
Electromagnetic stator coreSilicon steel2 × 10110.257600
Table 2. Structural Parameters of piezoelectric stator.
Table 2. Structural Parameters of piezoelectric stator.
ParametersValue (units)
b2 mm
λ6.67 π
f40.365 kHz
V200 V
n9
E72 × 109 N/m2

Share and Cite

MDPI and ACS Style

Li, Z.; Guo, P.; Wang, Z.; Zhao, L.; Wang, Q. Design and Analysis of Electromagnetic-Piezoelectric Hybrid Driven Three-Degree-of-Freedom Motor. Sensors 2020, 20, 1621. https://doi.org/10.3390/s20061621

AMA Style

Li Z, Guo P, Wang Z, Zhao L, Wang Q. Design and Analysis of Electromagnetic-Piezoelectric Hybrid Driven Three-Degree-of-Freedom Motor. Sensors. 2020; 20(6):1621. https://doi.org/10.3390/s20061621

Chicago/Turabian Style

Li, Zheng, Peng Guo, Zhe Wang, Liang Zhao, and Qunjing Wang. 2020. "Design and Analysis of Electromagnetic-Piezoelectric Hybrid Driven Three-Degree-of-Freedom Motor" Sensors 20, no. 6: 1621. https://doi.org/10.3390/s20061621

APA Style

Li, Z., Guo, P., Wang, Z., Zhao, L., & Wang, Q. (2020). Design and Analysis of Electromagnetic-Piezoelectric Hybrid Driven Three-Degree-of-Freedom Motor. Sensors, 20(6), 1621. https://doi.org/10.3390/s20061621

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