The Scientific World Journal

Volume 2012, Article ID 743138, 10 pages

http://dx.doi.org/10.1100/2012/743138

## Nonlinear Forced Response of Electromechanical Integrated Toroidal Drive to Coupled Excitation

Mechanical Engineering Institute, Yanshan University, Qinhuangdao 066004, China

Received 27 May 2012; Accepted 24 July 2012

Academic Editors: K. Prokeš and B. Yasilata

Copyright © 2012 Lizhong Xu and Fen Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The electric excitation and the parameter excitation from mesh stiffness fluctuation are analyzed. The forced response equations of the drive system to the coupled excitations are presented. For the exciting frequencies far from and near natural frequencies, the forced responses of the drive system to the coupled excitations are investigated. Results show that the nonlinear forced responses of the drive system to the coupled excitations change periodically and unsteadily; the time period of the nonlinear forced responses depends on the frequencies of the electric excitation, the mesh parameter excitation, and the nonlinear natural frequencies of the drive system; in order to improve the dynamics performance of the drive system, the frequencies of the electric excitations should not be taken as integral multiple of the mesh parameter exciting frequency.

#### 1. Introduction

Toroidal drive can transmit large torque in a very small size and is suitable for technical fields such as aviation and space flight [1–3]. As electrical and control techniques are utilized in mechanical engineering field widely, generalized composite drives become advancing edge of the mechanical science. So far, types of the generalized composite drives with integrated structure are still very limited.

The electromagnetic harmonic drive [4] and piezoelectric harmonic drive [5] are active drives in which the meshing forces between flexible gear and rigid one are controlled by electromagnetic force or piezoelectric one, and drive and power are integrated. Based on researching toroidal drive [6], the authors presented a kind of active generalized composite drive: electromechanical integrated toroidal drive. In the drive, the toroidal drive, power, and control are integrated [7].

The drive consists of four basic elements (Figure 1) (a) the central worm, (b) radially positioned planets, (c) a toroidal shaped stator, and (d) a rotor, which forms the central output shaft upon which the planets are mounted. The central worm is fixed and coils are mounted in helical grooves of its surface. The planets have permanent magnets instead of teeth. The and polar permanent magnets are mounted alternately on a planet. And the stator has helical permanent magnets instead of helical teeth. In the same manner as planet, the and polar helical permanent magnets are mounted alternately on the stator.

If a specific parameter relation is realized, pole of one element will correspond to pole of the other one all along. The attractive forces between and pole of the different elements are driving forces and the meshes without contact are realized. When the alternating voltage source is connected to the coils of the worm, a toroidal circular field is formed. It drives several planets to rotate about their own axial. And by means of magnetic forces between teeth of the planet and stator, the rotor is driven to rotate about its own axial. Thus, a power of low speed and large torque is output.

Compared with toroidal drive, the new drive is easy to produce, without wear, and does not need lubrication. It can be substituted for a servo system to simplify the structure of the existing electromechanical systems. Beside the above-mentioned fields that require compactness, the drive can be used in fields such as robots, which require accurate control.

The electromechanical integrated toroidal drive consists of a mechanical system, an electrical system and a coupled part. The mechanical vibration may occur in the mechanical system, and the electrical current oscillation may occur in the electrical system. With the coupled part, the mechanical vibration and the electrical current oscillation will influence on each other. Hence, the drive system is an electromechanical coupled dynamics system. The electromechanical coupled dynamics was first proposed for the motor [8]. Then, the electromechanical coupled dynamics of the electromechanical system consisting of the motor and mechanical system driven by the motor was developed. An electromechanical coupled dynamics model of the electromechanical system consisting of the several motors and mechanical system driven by these motors was proposed. Using the model, the natural frequency of the electromechanical system is analyzed [9]. The authors investigated nonlinear forced response of electromechanical integrated toroidal drive to voltage excitations [10]. However, not only the electric excitation occurs, but also the parameter excitations from mesh stiffness fluctuations occur. The nonlinear forced responses of the drive system to the coupled excitations consisting of the electric excitation and the parameter excitations have not been investigated yet. These nonlinear forced responses have important influence on the operating performance of the drive system. To design, evaluate, and control dynamics behavior of the drive system effectively, the nonlinear forced responses of the drive system to the coupled excitations should be developed.

In this paper, the electric excitation and the parameter excitation from mesh stiffness fluctuation are analyzed. The forced response equations of the drive system to the coupled excitations are presented. For the exciting frequencies far from and near natural frequencies, the forced responses of the drive system to the coupled excitations are investigated. The work can be used to predict the noise and dynamic load and are useful in maximizing the power density of the drive and reducing noise radiation.

#### 2. Electric Excitation of the Drive System

The magnetic energy storage of the electric system for the drive system is where is magnetic linkage of the worm coils, , is inductances of the th phase worm coils, is current of the th phase worm coils, and is phase number of the worm coils.

From (1), the electromagnetic torque on the planet is given as where is the relative rotating angle between planet and worm.

The torque consists of the static torque and the dynamic torque . The current consists of the static component and the dynamic one . The angle consists of the static angle and the dynamic angle . The torque is written in series form as If , from (3), one can give

Let denote the dynamic magnetic meshing force between a planet tooth and worm, . From (5), neglecting the high-order terms, one knows.

Here, ; it is the dynamic relative displacement between planet and worm. Thus, (6) can be changed into the following form: Here, , it is considered as linear electromagnetic mesh stiffness between a tooth of the planet and worm; , it is nonlinear component of the electromagnetic mesh stiffness.

Equation (3) shows that current fluctuation can produce electromagnetic torque fluctuation. The torque fluctuation from current fluctuation can be considered as equivalent exciting torque . Let the current periodically change as below: where is the magnitude of the fluctuation current, is the frequency of the fluctuation current, and is time.

Substituting (8) into (3), neglecting high-order terms, yields

Let denote the equivalent exciting force between a planet tooth and worm; thus

#### 3. Parametric Excitation from Stiffness Fluctuation

In operation of the drive system, the number of meshing tooth pairs between the worm and the planets is variable. It can cause fluctuation of the mesh stiffness between them. It is dependent on the number of planet teeth and the conditional face angles of the worm . At ° and , the changes of mesh stiffness along with the mesh tooth pair number are shown in Figure 2. In Figure 2, a typical stiffness variation through a mesh cycle of the drive system is given. The mesh stiffness at more contact regions is higher than that at less contact regions. If the drive rotates at appreciable speed, this time-varying stiffness will be a major excitation source of the drive system.

The mesh stiffness between a planet and worm is considered to consist of its mean value and time-varying one . The average mesh stiffness between them through one periodic time can be given by where and ; they are two and three teeth mesh stiffness, respectively.

The periodical time-varying portion of the mesh stiffness can be defined in the Fourier series form as where , is the period of the stiffness fluctuation. In Figure 2, .

For worm and planet, the periodical time-varying portion of the mesh stiffness is

Substituting (13) into (12) yields

#### 4. Forced Response Equation to Coupled Excitation

The dynamic model for the electromechanical integrated toroidal drive (see Figure 3) allows rotor and each planet to rotate about their own rotating axes and allows each planet to translate in and directions. The rotations are replaced by the corresponding translational mesh displacements as , (here, is planet number, the rotation of planet or rotor, is the rolling circle radius for planet and the radius of the circle passing through planet centers for the rotor). A displacement vector and a mass matrix are defined for each planet as and . Here, and are polar mass moment of inertia and mass for planet , respectively. is equivalent mass of rotor corresponding to its displacement . Thus, the motion equations of the drive system are where , and are relative displacements between planet- and worm, stator, or rotor, respectively, , and . and are lead angles at contact points between planet- and worm or stator, respectively, , and . Here, is center distance between worm and planet is reference circle radius of planet and are speed ratios between planet and worm or stator, respectively. and are mesh stiffness between planet- and worm or stator, respectively, and are planet support stiffness in and directions, respectively, and is torque transmitted by rotor.

Equation (15) can be written in matrix form as where and are displacement and static load vectors; respectively, , ; and are mass and stiffness matrix, respectively, ,

is nonlinear component of the force vector, is equivalent force vector caused by current fluctuation:

Consider

it is a disturbing parameter, , .

The total displacement of the each component consists of the static displacement and the dynamic one. Hence, the displacement vector is where and are static and dynamic displacement vectors; respectively,

Considering the fluctuation of the mesh stiffness caused by the changes of the mesh tooth pair number, the stiffness matrix can be expressed as where is the mean stiffness matrix, and is time-varying stiffness matrix.

Substituting (20)–(22) into (16) yields where ; it is considered as the equivalent exciting force vector caused by mesh stiffness excitation; .

Equation (24) is just the forced response equation of the drive system to coupled excitation.

#### 5. Solution of the Nonlinear Forced Response Equation

##### 5.1. Far from Natural Frequencies

First, we consider to resolve the nonlinear forced response equation of the drive system to electric excitation (let ).

For simplicity purposes, (24) should be transmitted into equations independent on each other. Then, (24) is changed into the following form: where is the diagonal mean stiffness matrix, respectively. is transmitted dynamic displacement vector. and are transmitted exciting force vectors of the forces and , respectively.

Matrix and vector are given by and . The transmitted nonlinear force vector and transmitted equivalent exciting force vector can be given as below: where is the mode matrix of (24):

At , , = , ; , + , .

Let where is natural frequency of th order mode for the nonlinear drive system, and is natural frequency of th order rotational mode for the linear drive system.

Substituting (28) and (29) into (25), let sum of the coefficients with the same-order power of the parameter equal zero, following equations can be given Here, initial conditions are , .

The solution of zero-order equation under the above initial conditions is

Substituting (31) into the second equation of (30) yields

The rotational displacement is Substituting (33) into (30) yields

In order to remove secular item, let

Substituting (35) into (34), the solutions of the first-order equations can be obtained. As the equations of the solutions are relatively complicated, it is not given here.

In a same manner, the solution of th order equation can be obtained as well. Substituting these solutions into (28) and (29), the solutions of the regular nonlinear forced response equations and natural frequencies of the drive system can be given. Then, the real solutions of the nonlinear forced responses can be calculated as below:

##### 5.2. Near Natural Frequencies

When the exciting frequency is near natural frequency, the nonlinear forced response equation can be resolved as below. Considering damping of the drive system, (25) can be changed into where is regular damping matrix, ( is damping matrix). , is relative damping coefficient.

When exciting frequency is near natural frequency, the exciting frequency can be written as

Substituting (38) and (28) into (37), let sum of the coefficients with the same-order power of the parameter equal zero, the following equations can be given: where ,

The solution of zero-order equation under the above initial conditions is

Substituting (41) into the second equation of (43) yields

Substituting rotational displacement into (42) yields

In order to remove secular item, let where .

From (44), it is known that

Thus where and .

From (46), the changes of the nonlinear vibrating magnitudes along with exciting frequencies can be given.

#### 6. Results and Discussions

When exciting frequency is far from natural frequency, from the above equations, the nonlinear forced vibrations for the drive system are analyzed. The parameters of the numerical example are shown in Table 1. Figure 4 shows changes of the forced vibrations of the transmitted variables along with nonlinear parameter . From Figure 4, the following are known.(1)Under the coupled excitations, the nonlinear forced responses of the drive system change periodically and unsteadily. The time period of the nonlinear forced responses depends on the frequencies of the electric excitation the mesh parameter excitation, and the nonlinear natural frequencies of the drive system.(2)The vibrating amplitudes of the nonlinear forced responses of the drive system to the coupled excitations are larger than those of the nonlinear forced responses to the single excitation. The vibrating amplitudes of the nonlinear forced responses of the planet are larger than that of the rotor. It is because the exciting frequencies are near to the vibrating frequency of the planet modes.(3)The vibrating amplitudes of the tangent vibration for the planet are larger than that of the axial vibration for the planet, and the frequency of the tangent vibration for the planet is smaller than that of the axial vibration for the planet.(4)The unstable periodic vibrations are harmful for the drive operation. In order to increase the dynamics performance of the drive system, the electric excitations should be avoided or the frequencies of the electric excitations should not be taken as integral multiple of the mesh parameter exciting frequency.

When exciting frequency is near natural frequency, changes of the nonlinear vibrating magnitudes along with exciting frequencies and the drive parameters are given in Figure 5. From Figure 5, the following are known.(1)As exciting frequency increase, the vibrating magnitudes of all the modes increases, and at point , the vibrating magnitudes get to the maximum, and then they decrease with increasing exciting frequency. For different modes, the previously mentioned curves bend toward the direction of the exciting frequency increase. The results are typical nonlinear character of the drive system.(2)For the first and third modes, as nonlinear parameter increases, their vibrating magnitudes increase. For the second mode, as nonlinear parameter increases, its vibrating magnitude decreases. The vibrating magnitudes are large for the first and third modes and small for the second mode. The vibrating magnitude for the first mode is larger than that for the third mode. (3)For a given exciting frequency, the vibrating amplitudes of the nonlinear forced responses of the drive system to the coupled excitations for modes 1 and 3 are larger than those of the nonlinear forced responses to the single excitation; for mode 2, the vibrating amplitude of the nonlinear forced responses of the drive system to the coupled excitations is smaller than that of the nonlinear forced responses to the single excitation. (4)Under the coupled excitations, the nonlinearity has obvious effects on the relationship between the frequency and the amplitudes; it should be considered.

#### 7. Conclusions

In this paper, the electric excitation and the parameter excitation from mesh stiffness fluctuation are analyzed. The forced response equations of the drive system to the coupled excitations are presented. For the exciting frequencies far from and near to natural frequencies, the forced responses of the drive system to the coupled excitations are investigated. Results show the following.(1)Under the coupled excitations, the nonlinear forced responses of the drive system change periodically and unsteadily. (2)The time period of the nonlinear forced responses depends on the frequencies of the electric excitation the mesh parameter excitation, and the nonlinear natural frequencies of the drive system.(3)In order to improve the dynamics performance of the drive system, the frequencies of the electric excitations should not be taken as integral multiple of the mesh parameter exciting frequency.

#### Acknowledgment

This project is supported by the National Natural Science Foundation of China (no. 51075350).

#### References

- M. R. Kuehnle, “Toroidgetriebe,” Urkunde uber die Erteilung des deutschen Patents, 1301682, 1966.
- M. R. Kuehnle, H. Peeken, C. Troeder, and S. Cerniak, “The Toroidal drive,”
*Mechanical Engineering*, vol. 103, no. 2, pp. 32–39, 1981. View at Google Scholar · View at Scopus - M. R. Kuehnle, “Toroidal transmission and method and apparatus for making and assembling same,” United States Patent, 5863273, 1999.
- Z. Delin and L. Huamin, “Side surface harmonic stepper motor,”
*Journal of Mechanical Engineering*, vol. 29, no. 5, pp. 96–98, 1993 (Japanese). View at Google Scholar - O. Barth, “Harmonic piezodrive—miniaturized servo motor,”
*Mechatronics*, vol. 10, no. 4, pp. 545–554, 2000. View at Publisher · View at Google Scholar · View at Scopus - L. Xu, Z. Huang, and Y. Yang, “Contact stress for toroidal drive,”
*Journal of Mechanical Design, Transactions of the ASME*, vol. 125, no. 1, pp. 165–168, 2003. View at Publisher · View at Google Scholar · View at Scopus - L. Xu and J. Huang, “Torques for electromechanical integrating toroidal drive,”
*Proceedings of the Institution of Mechanical Engineers, Part C*, vol. 219, no. 8, Article ID C15704, pp. 801–811, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. Shaltout and M. Al-Omoush, “Reclosing torques of large induction motors with stator trapped flux,”
*IEEE Transactions on Energy Conversion*, vol. 11, no. 1, pp. 84–90, 1996. View at Google Scholar · View at Scopus - H. Y. Kanaan, K. Al-Haddad, and G. Roy, “Analysis of the electromechanical vibrations in induction motor drives due to the imperfections of the mechanical transmission system,”
*Mathematics and Computers in Simulation*, vol. 63, no. 3–5, pp. 421–433, 2003. View at Publisher · View at Google Scholar · View at Scopus - L. Xu and X. Hao, “Nonlinear forced vibration for electromechanical integrated toroidal drive,”
*International Journal of Applied Electromagnetics and Mechanics*, vol. 28, no. 3, pp. 351–369, 2008. View at Google Scholar · View at Scopus