`Journal of Applied MathematicsVolume 2012, Article ID 465023, 10 pageshttp://dx.doi.org/10.1155/2012/465023`
Research Article

## An Optimal Approach to Study the Nonlinear Behaviour of a Rotating Electrical Machine

1Faculty of Mechanical Engineering, “Politehnica” University of Timisoara, M. Viteazu Boulevard, No. 1, 300222 Timisoara, Romania
2Center of Advanced Research in Engineering Sciences, Romanian Academy, Timisoara Branch, M. Viteazu Boulevard, No. 24, 300223 Timisoara, Romania

Received 1 August 2012; Accepted 13 September 2012

Academic Editor: Livija Cveticanin

Copyright © 2012 Nicolae Herişanu and Vasile Marinca. 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

We introduce a novel analytical approximate technique, called Optimal Variational Method (OVM), to investigate the nonlinear behaviour of a rotating electrical machine modelled as an oscillator with cubic elastic restoring force and time variable coefficients. The proposed procedure involves the presence of some initially unknown convergence-control parameters whose values are later optimally determined. Comparisons between the obtained results and exact ones reveal that OVM is very effective and convenient providing highly accurate results.

#### 1. Introduction

Mathematical modelling of many physical and engineering systems often leads to nonlinear ordinary or partial differential equations. In order to analyze such mathematical models, an effective method is required to provide solutions conforming to physical reality. In some cases, inherent difficulties are overcome by replacing a nonlinear differential equation with a corresponding linear differential equation that approximates the original one close enough to give useful results, but this is only a harsh approximation.

There are some known approaches intended to obtain approximate solutions to nonlinear dynamical systems. The most common and widely studied methods for solving nonlinear differential equations are the perturbation methods [13], but almost all perturbation methods are based on such an assumption that a small parameter must exist in the equation under investigation. The approximate solutions obtained through perturbation methods are valid, in most cases, only for small values of parameters, and moreover, there is no criterion on how small the parameters should be.

In order to overcome these limitations related to small parameters, several effective methods were developed, such as the method of Lie group [4], the Adomian decomposition method [5], the weighted-linearization method [6], and the method of harmonic balance [7].

The aim of this paper is to propose a novel variational approach to study the nonlinear behaviour of a rotating electrical machine modelled as an oscillator with cubic elastic force and time variable coefficients.

Rotating electrical machines are complex dynamical systems exhibiting typical problems from rotor dynamics, especially for high speeds or heavy duty rotor shafts with high inertial loads [8, 9]. An important engineering challenge in this field is to develop such mathematical models and methods, which allows quite accurate predictions of the dynamical phenomena at the stage of design, when compared against experimental measurements obtained after the stage of machine building. These models and methods should be made available in order to provide an easy-to-use tool to operate design modifications when these may be required to suppress undesired dynamic phenomena.

The most common sources of vibration in rotating machines are related to the unbalanced forces of the rotor, shaft misalignment and nonlinearity of the bearing stiffness, variable elasticity, bad bearings and mechanical looseness, and other electrical and mechanical faults which generate nonlinear vibration in the system.

These phenomena should be controlled in order to make the machine run smoothly and reliably. It is well known that rotating parts cannot be perfectly balanced. From engineering point of view, it is impossible to make any rotor perfectly mass balanced. Hence, residual unbalance is always present to some extent, even though the rotating structure is well constructed, but if it deteriorates, then damaging vibration occurs. Supplementary problems could arise in case of some horizontal rotating machines, when the gravity effect is not negligible for certain stiffness conditions. In case of gravity deflection, the shaft center leaves the bearing centerline, which leads to vibration occurrence. Misalignment could also occur in the electrical rotating machine after some amount of running. All these could be highly detrimental, affecting the integrity of the system [8, 9].

In this paper, a new analytical procedure, namely, the Optimal Variational Method (OVM) is employed to study the problem of nonlinear vibrations of an electric machine supported by nonlinear bearings characterized by nonlinear stiffness of Duffing type while the entire system is subject to a parametric excitation caused by an axial thrust and a forcing excitation caused by an unbalanced force of the rotor. In these conditions, the dynamical behaviour of the investigated electrical machine will be governed by the following second-order strongly nonlinear differential equation [8]: with the initial conditions which can be written in the more convenient form where and the dot denotes derivative with respect to time and is the initial amplitude of the oscillations. Note that it is unnecessary to assume the existence of any small or large parameters in (1.3).

Equation (1.3) describes a system oscillating with an unknown period . We switch to a scalar time . Under the transformation the original Equation (1.3) becomes with the initial conditions where the prime denotes derivative with respect to .

#### 2. A Novel Variational Method and Solutions

In order to show the basics of OVM, we consider the following differential equation: The variational principle for (2.1) can be easily established if there exists a functional which admits as extremals the solutions of (2.1), where is the Lagrangian of the system (2.1) and [] is the domain of interest.

This problem is based on the study of the conditions under which there exists a function   such that Euler’s equation of the functional (2.2) coincide with the system (2.1), that is,

On physical grounds, the primary significance of this problem, called “the inverse problem in Newtonian mechanics” [10] rests on the fact that the acting forces in Newtonian system (2.1) need not necessarily be derivable from a potential. Equation (2.2) is called action functional or action for short.

In our procedure, we assume that the approximate solution of (2.1) depends on several parameters : such that the action functional (2.2) becomes

The parameters from (2.5), called convergence-control parameters, can be determined optimally applying the Ritz method [8]:

From (2.6) and from the initial condition (1.7)1 which becomes we can obtain optimally the parameters and the frequency of the system (2.1).

We remark that the condition (1.7)2 is identically verified by the solution (2.4). On the other hand, the expression of the solution (2.4) is not unique.

The validity of the proposed approach is illustrated on (1.6). In this case, the Lagrangian of (1.6) can be written as

If we consider in (2.4), then the approximate solution of (1.6) can be written as

Also, we can choose this approximate solution in the form or and so on.

Substituting (2.9) into (2.8) and this into (2.10) we have the following results for and : where

The values of the parameters and the frequency Ω are obtained from (2.6) which become: and from the initial condition (2.7) which can be written as

#### 3. Test Examples

The validity of the proposed procedure for solving the investigated problem is illustrated on three examples, considering different parameters and initial amplitudes.

Case a. As a first example, we consider the following set of parameters: or for the corresponding coefficients of (1.6):
The parameters and the frequency Ω are obtained from (2.14) and (2.15):
The approximate solution of (1.1) in this case becomes
Figure 1 shows the comparison between the present solution (3.4) and numerical integration results obtained using a fourth-order Runge-Kutta scheme.

Figure 1: Comparison between the numerical solution of (1.1) and the approximate solution (3.4) for , , , , , , , .

Case b. As a second example, we consider the following parameters: Keeping the rest of the coefficients unchanged and following the same procedure, we obtain the optimal values of the convergence-control parameters ,  , and the frequency Ω: and consequently, the approximate solution of (1.1) becomes in this case:
The comparison between the present solution (3.7) and numerical integration results in the second case is presented in Figure 2.

Figure 2: Comparison between the numerical solution of (1.1) and the approximate solution (3.7) for , , , , , , , .

Case c. Finally, for the last example we consider the following set of parameters: In this case, following the procedure described above, we obtain and the approximate solution of (1.1) becomes
The comparison between the obtained approximate analytical solution (3.10) and numerical integration results in the third case is presented in Figure 3.

Figure 3: Comparison between the numerical solution of (1.1) and the approximate solution (3.10) for , , , , , , , .

It can be observed from Figures 13 that the approximate analytical results obtained using the proposed procedures are in good agreement with the numerical ones.

#### 4. Conclusions

In this work, an optimal variational approach is employed to propose a new analytic approximate solution for nonlinear conservative oscillations of an electrical machine. The construction of our variational approach is different from traditional approach especially concerning the involvement of some initially unknown parameters called convergence-control parameters, whose optimal values ensure a fast convergence of the approximate analytical solutions. This is in fact the main advantage of the OVM, which provides us with a simple and rigorous way to control and adjust the convergence of the solutions. The proposed procedure is valid even if the nonlinear equation does not contain any small or large parameters. Three test examples illustrate that the proposed analytical approach is very effective and yields accurate results comparing to those obtained via numerical integration using a fourth-order Runge-Kutta method.

#### References

1. A. H. Nayfeh, Problems in Perturbation, John Wiley & Sons, New York, NY, USA, 1985.
2. J. J. Stoker, Nonlinear Vibrations in Mechanical and Electrical Systems, John Wiley & Sons, New York, NY, USA, 1986.
3. M. P. Cartmell, S. W. Ziegler, R. Khanin, and D. I. M. Forehand, “Multiple scales analysis of the dynamics of weakly nonlinear mechanical systems,” Applied Mechanics Reviews, vol. 56, no. 5, pp. 455–492, 2003.
4. P. J. Olver, Application of Lie Group to Differential Equation, Springer, New York, NY, USA, 1986.
5. G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.
6. V. P. Agrwal and N. N. Denman, “Weighted linearization technique for period approximation in large amplitude nonlinear oscillators,” Journal of Sound and Vibration, vol. 99, pp. 463–473, 1985.
7. R. E. Mickens, Oscillations in Planar Dynamics Systems, vol. 37, World Scientific, Singapore, 1996.
8. V. Marinca and N. Herişanu, Nonlinear Dynamical Systems in Engineering. Some Approximate Approaches, Springer, Berlin, Germany, 2011.
9. V. Marinca and N. Herişanu, “Nonlinear dynamic analysis of an electrical machine rotor-bearing system by the optimal homotopy perturbation method,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 2019–2024, 2011.
10. R. M. Santilli, Foundations of Theoretical Mechanics I, The Inverse Problem in Newtonian Mechanics, Springer, Berlin, Germany, 1978.