The traveling wave ultra-sonic motor (TWUSM) is operated through the frictional force between the rotor and the stator. Hence, the contact mechanism should be analyzed to estimate the motor performance. However, the nonlinearity of the contact mechanism of the TWUSM makes it difficult to propose a proper contact model and a characteristic analysis method. To address these problems, a novel contact model is proposed and be termed the cylindrical dynamic contact model (CDCM) in this research. An estimation method of the motor performance is proposed using the CDCM, an analytical method, and a numerical method. The feasibility and usefulness of the proposed characteristic analysis are verified through experimental data.
The nonlinearity of the contact mechanism of the TWUSM has been made it difficult to propose an effective characteristic analysis method
. Hence, many researchers have investigated the contact mechanisms of TWUSMs over the past few decades. Numerous types of characteristic analysis methods using diverse kinds of contact models have been proposed for the TWUSM by many researchers
. However, some of the proposed methods have to use the unreasonable values of material coefficients or inconsistent material coefficients to fit the calculated data to the simulated results. In some methods, the governing equations are complex to be solved. Some methods can analyze the partial characteristics of the motor performance. To address these conventional problems, a novel contact model for the TWUSM, which referred to CDCM, is proposed in this research. Furthermore, a characteristic analysis method using the proposed CDCM, an analytical method, and a numerical method for the TWUSM is proposed in this paper. The proposed method for the TWUSM is verified by experimental data with a prototype of the TWUSM.
2. Working Principle of the TWUSM
illustrates the working principle of the TWUSM. For the generation of the traveling wave, the piezoelectric materials are placed in the stator with a proper space difference. By applying electrical power to the piezoelectric materials with a proper time difference, the travelling wave is generated at the stator. As shown in
, by contacting the rotor to the stator and applying mechanical pressure, the rotor moves by the frictional force between the rotor and the stator.
Working principle of the TWUSM.
3. Characteristic Analysis Method using the CDCM
- 3.1 CDCM
As shown in
, we assumed the stator and the rotor of the TWUSM as cylindrical bodies with the radius of
, respectively. To reflect the dynamic contact condition, the CDCM is proposed as illustrated in
, describing the nonlinear deformation at the area of contact. It is assumed that both rolling and sliding occur between the stator and the rotor while they are in contact.
Modeling of a rotor and a stator using cylinders.
The proposed CDCM for the TWUSM.
- 3.2 Analysis of TWUSM using the FEM
The first step for the characteristic analysis of the TWUSM is the calculation of the operating frequency
(Hz) and the displacement amplitude
(m) of the TWUSM by means of 3D-FEM. We developed the 3D-FEM software and an auto-mode search program. The operating frequency can be calculated through an impedance analysis using FEM and (1), in which
) is the impedance at the excitation frequency
) (V) is the electric potential on the electrode, and
(C) is the electric charge.
- 3.3 Formula for the motion of a stator
Based on Kirchhoff’s plate theory, the neutral plane of the stator can be expressed in terms of an ideal traveling wave as (2) while assuming that the stator is not affected by the contact force
is the position in the radial direction,
is the amplitude at radius
is the number of nodal lines.
The three-dimensional problem can be replaced by a two-dimensional problem, as expressed by (3).
which is subject to
is the normal displacement amplitude of the neutral plane for a stator, λ is the wavelength,
(m) is the effective length of the circumference of the contact area,
(m) is the effective radius of the contact area, and
(m) are respectively the outer and the inner diameter of the contact area. The contact problem is described in a fixed reference (
, ẑ ) in (3). This can be transformed to a moving reference (
) using (7) and (8), where
is the wave propagation velocity (m/sec).
The moving reference follows the wave crest. Hence, (3) can be formulated into the time-independent Eq. (9) by substituting (7) and (8) into (3). The normal displacement of the neutral plane of the stator in the moving reference can be derived as (10).
- 3.4 Analysis of the contact length
When two cylinders, which are the stator and the rotor, come into contact with each other, as shown in
, the half of the contact length
can be calculated by (11), where
are the material coefficients of the stator and the rotor, respectively;
is the shape coefficient of the cylinder;
(N/m) is the normal force per unit length in the radial direction at the unit wave;
(Pa) are the equivalent elasticity modulus for the stator and the rotor, respectively;
are the Poisson’s ratios of the stator and the rotor, respectively;
(N) is the normal force applied to the TWUSM for the frictional force; and
is the contact width in the radial direction.
which is subject to
The equivalent elasticity moduli
should be calculated while taking the composite material of the TWUSM into account. The stator consists of piezoelectric ceramic material with a metal medium. The rotor is composed of the contact material and a metal medium. The equivalent elasticity modulus of the composite material can be derived by (15) when two materials are stacked in the thickness direction; their elasticity moduli are
with thicknesses of
The shape coefficient of cylinder
can be calculated by (16)
can be calculated using the displacement of the stator by (17).
is an infinite value because the rotor is a plate.
- 3.5 Analysis of the normal pressure distribution
Because the normal pressure distribution is proportional to (10), which is the normal displacement of the neutral plane of the stator, the normal pressure distribution
) on the contacting surface can be derived by (18), where
) is the amplitude of the normal pressure distribution before the deformation of the stator and where
is an unknown constant. The unknown constant
accounts for the boundary conditions, which state that
should be equal to zero when
can be found by (19). The unknown function,
, can be found using condition (20). Eq. (21) can be induced from the condition under which the normal pressure distribution is symmetric about
=0 and by substituting (18) and (19) into (20). From(21), the amplitude of the normal pressure distribution
is derived by (22).
Finally, the governing equation for the analysis of the normal pressure distribution as a function of the normal force applied to a TWUSM,
, is derived by (23). The maximum pressure
can be calculated by
- 3.6 Analysis of the stress and the strain
In the case of dynamic contact, the stress can be expressed as normal stress and tangential stress
. The normal stress,
) , which is induced by normal force can be derived by (24). The tangential stress,
) , which is generated by the tangential frictional force can be expressed by (25).
which are subject to
is the dynamic frictional coefficient.
The stress for the dynamic contact can be expressed by (31). From the Hook’s law, Eq. (32) is derived. The strain,
), can be induced by applying the infinitesimal thickness in the z direction, Δ
, as (33).
Hence, the deformation at a segmented element,
), can be determined by (34). The total deformation in the
), is derived by (35).
- 3.7 Analysis of the displacement and the velocity distribution of the stator taking the deformation into consideration
The displacement amplitude of the stator decreased nonlinearly with the deformation caused by the applied normal force and the dynamic contact condition. The deformed normal displacement equation of the stator will be an asymmetric non-linear function due to the dynamic contact condition. Hence, an approximate method is proposed in this research to describe the nonlinear traveling wave of a deformed stator. The normal displacement of the neutral plane of a deformed stator in the plus
, and in the minus
, in the moving reference can be derived by (36) and (38), respectively. The traveling wave of the neutral plane of a deformed vibrator in the plus
), and in the minus
), in the fixed reference can be induced by (37) and (39), respectively. The total number of unknown constants of (36-39) is eight, in which
can be found as the normal displacement when
equal 0. The remaining six unknown constants can be calculated by solving a 3×3 inverse matrix using the three
position and the associated normal displacement values for each plus and minus contact area.
The horizontal displacement at the contacting surface of the stator can be derived by applying Kirchhoff’s plate theory to (36-39). The horizontal displacement for the plus
area in the fixed reference,
), can be induced by (40), where
is the distance from the neutral plane to the contact surface.
The horizontal velocity of the stator at the contact surface in the fixed reference,
), can be derived by (41) for the plus
The horizontal velocity of the stator at the contact surface in the moving reference,
, can be induced by (42) for the plus
area by applying (7) and (8) to (41).
The governing equations for the minus
area can be suggested in the same way as (40-42).
- 3.8 Analysis of the normal pressure distribution considering the deformation of a stator
The normal pressure distribution induced by a nonlinear displacement function which considers the deformation of the stator at the plus
contact area of
) can be determined by (43), where
is the amplitude of the final normal pressure distribution and
is an unknown constant at the plus
contact area. For the minus
area, the normal pressure displacement can be derived in the same way as (43). The unknown constants can be solved by (44) and (45).
The unknown functions
can be induced by conditions (46) and (47).
The normal force per unit length in the radial direction at unit wave,
, is constituted with the force of the plus
and that of the minus
, as expressed by (48) and (49), respectively.
For the plus
contact area, the unknown function
can be determined using (43-49), as follows:
Eq. (50) can be modified to (51). For the minus contact area, (52) can be derived from (50).
where the integration part can be expressed by:
Finally, the normal pressure distribution taking into account the nonlinearly deformed stator is proposed in this research as (54) and (55).
- 3.9 Analysis of the motor performance
The horizontal velocity distribution of a deformed stator considering dynamic contact has an asymmetric shape, as shown in
. Hence, the slip distribution between the rotor and the stator can be described by
Distribution of the horizontal velocity at the stator and the rotor considering the nonlinear deformation induced by dynamic contact.
Slip distribution over the nonlinearly deformed contact area caused by dynamic contact.
, the frictional force acting on the surface of the stator at unit wave,
(N/m), can be induced by (56), where
(rpm) denotes the revolutions per minute of the rotor.
Free-body diagram to elucidate the forces acting on the rotor and the stator.
when (0 ≤
≤ 0) ,
which is subject to,
sign[τ] = 1, if τ > 0
sign[τ] = 0, if τ = 0
sign[τ] = 1, if τ < 0
The total frictional force acting on the surface of the stator,
(N), can be derived by (59), which is subject to (57) and (58).
The total load force acting on the rotor,
can be expressed by
. It has the same meaning as the total friction force of the rotor,
(N), as induced by (60), which also can be used for the calculation of the thrust force of the linear TWUSM.
(Nm) and the speed
(rpm) of the ring-type TWUSM can be analyzed by (61), which is subject to (57) and (58).
4. Verification of Analysis Method
As shown in
, the stator is prototyped and its impedance is measured and calculated. As displayed in
, the calculated data are in good agreement with the experimental results. Hence, it is verified that the 3D-FEM used in this research is correct.
Shape of prototype and its impedance.
The speed-torque is analyzed using the proposed characteristic analysis method derived from CDCM. These results are then compared to the experimental data, as illustrated in
. The analysis results are fitted well with the experimental data. These results prove that the proposed characteristic analysis method is correct. It takes about 1-2 seconds for the analysis using the proposed characteristic analysis method. Hence, the usefulness of the proposed analysis method is validated. As displayed in
, the stress distribution at the contact area of the stator shows an asymmetric and nonlinear shape due to the consideration of dynamic contact. The material coefficients used for the characteristic analysis are included in APPENDIX and are shown in
Speed-torque curve of TWUSM.
Calculated stress distribution at the stator using the CDCM.
Parameters for the analysis of the TWUSM
Parameters for the analysis of the TWUSM
The contact mechanism of the TWUSM is a complex nonlinear problem, which makes it difficult to propose a useful contact model and a feasible motor performance analysis method. Hence, it is noteworthy from the fact that a correct and useful characteristic analysis method using CDCM is proposed in this research.
Jong-Suk Ro He received the B.S. degree in Mechanical Engineering from Han-Yang University, Seoul, Korea, in 2001. In 2008, he earned a PhD in Electrical engineering from Seoul National University, Seoul, Korea, through the Combined Master’s and Doctorate Program. He conducted research on cellular phone modules at R&D center of Samsung Electronics as a Senior Engineer from 2008 to 2012. From 2012 to 2013, he was at Brain Korea 21 Information Technology of Seoul National University as a Post-Doctoral Fellow. Currently, he is carrying out research at Electrical Energy Conversion System Research Division of Smart Grid Team at Korea Electrical Engineering & Science Research Institute as a Researcher. His research interests are interdisciplinary research, energy conversion device, camera module, numerical analysis, and optimal design of electric machines.
Kyung-Pyo Yi He received the B.S. and M.S. degree from the School of Electrical and Electronics Engineering, Seoul National University, Seoul, Korea in 2007 and 2009, respectively. He is currently a Ph.D. student at Seoul National University, Seoul, Korea. His research focuses on analysis and design of piezoelectric devices, electric motors and superconducting magnetic energy storage(SMES).
Tae-Kyung Chung He received the B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National University in 1981, 1983, and 1987, respectively. From 1986 to 1987, he worked at Daewoo Heavy Industries as a senior researcher. Since 1988, he has been with the School of Electrical and Electronics Engineering at Chung Ang University in Korea as a professor. In 1991, he was with Harvey Mudd College in California as a visiting scholar. His main research fields are design of electric motors and actuators.
Hyun-Kyo Jung He received the B.S., M.S., and Ph.D. degree in Electrical engineering from the Seoul National University, Seoul, Korea, in 1979, 1981, and 1984, respectively. From 1985 to 1994, he was a member of the faculty with Kangwon National University. From 1987 to 1989, he was with the Polytechnic University of Brooklyn, Brooklyn, NY. From 1999 to 2000, he was a Visiting Professor with the University of California at Berkeley. He is currently a Professor at the School of Electrical Engineering and Computer Science/Electrical Engineering, Seoul National University. His research interests are the analysis and design of the electric machine.
“The effect of tangential elasticity of the contact layer between stator and rotor in traveling wave ultrasonic motors,”
International journal of NonLinear Mechanics
DOI : 10.1016/S0020-7462(01)00048-8
“Optimal design of a rod shape ultrasonic motor using sequential quadratic programming and finite element method,”
Sensors and Actuators
“Characteristics of Disk-type Linear Ultrasonic Motor for Application to x-y Stage,”
Journal of Electrical Engineering & Technology
DOI : 10.5370/JEET.2006.1.1.101
“Modal analysis and control of a bowl parts feeder activated by piezoceramic actuators,”
Journal of Sound and Vibration
DOI : 10.1016/j.jsv.2003.10.008
“Optimal design of ultrasonic motor using evolution strategy and finite element method,”
International Journal of Applied Electromagnetics and Mechanics
“Revolution speed characteristics of an ultrasonic motor estimated from the pressure distribution of the rotor,”
Proc. 12th Symp. on Ultrasonic Electronics
Bogy D. B.
“Effect of the hydrodynamic bearing on rotor/stator contact in a ring-type ultrasonic motor,”
IEEE Transactions on Ultrasononics, Ferroelectrics,. and Frequency Control
DOI : 10.1109/58.165550
“Design of a Traveling Wave Type Ultrasonic Motor,”
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
Flynn A. M.
“Piezoelectric ultrasonic micromotors,” PhD Thesis
Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology
“Estimation of the tangential stresses in the stator/rotor of traveling wave ultrasonic motors using visco-elastic foundation model,”
Proceedings of 2nd Int. Conference on Contact Mechanics
Schmidt J. P.
“A note on the contact problem in an ultrasonic travelling wave motor,”
Int. J. Nonlinear Mech.
DOI : 10.1016/S0020-7462(96)00112-6
“Performance evaluation of traveling wave ultrasonic motor based on a model with visco-elastic friction layer on stator,”
DOI : 10.1016/j.ultras.2006.05.217
Smith J. O.
Lui C. K.
“Stresses due to tangential and normal loads on an elastic solid with application to some contact stress problems,”
J. App. Mech. Trans. ASME
Norton R. L.
MACHINE DESIGN: An Integrated Approach, Chapter 6
Pearson Education. Inc.