Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 238714, 9 pages

http://dx.doi.org/10.1155/2015/238714

## Application of the Improved Complex Normal Form Method Based on He’s Energy Balance Method to a Noise Included Nonlinear Oscillator

School of Mechanical Engineering, Tianjin Polytechnic University, Tianjin 300192, China

Received 8 December 2014; Revised 8 April 2015; Accepted 19 April 2015

Academic Editor: P. Balasubramaniam

Copyright © 2015 Gen Ge and Zheng Chao Yun. 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

A strongly nonlinear oscillator with nonlinear damping subject to narrow band random excitation is discussed. The deterministic and the stochastic responses are studied, respectively, applying the improved complex normal form method based on He’s energy balance method. The accuracy of detecting of the response frequency has been improved significantly. The behavior of steady state responses together with the stability is analyzed by the moment method. The numerical simulation results are consistent with theoretical analysis. When the noise intensity is small enough, the analytical result given by the moment method agrees with the simulation quite well.

#### 1. Introduction

Strongly nonlinear vibration problems have drawn accelerating interest in the past two decades. People have explored various approximate methods to handle the strongly nonlinear problems. A number of them have been developed dramatically and widely applied in analysis of deterministic response problems of nonlinear systems. Actually, some of the approximate methods for deterministic response problems of nonlinear systems can be extended to random response problems. Recently, the modified Lindstedt-Poincare method (MLP) and the multiple time scale method [1–6] are the most popular ways to study noise included nonlinear problems. The basic issue of the two methods is utilizing an undetermined frequency to construct a new expanding parameter which is small enough to change strongly nonlinear systems into weak ones for convenience of using perturbing approaches. Impressive studies have been proposed by Rong and Xu et al. [7–11], and a Duffing-Rayleigh type strongly nonlinear system subject to narrow band noise is studied in detail. In [11], a parameter transformation combined with multiple time scale method is used to change a strongly nonlinear problem into a weak nonlinear problem, where is a mark denoting the strength of the nonlinear terms which are not small and is the amplitude of the first mode of the system. This type of parameter transformation is firstly used by Burton [12, 13] to solve deterministic strongly nonlinear problems. In [9, 10], another type of parameter transformation is introduced, where the undetermined frequency is expressed as and denotes the fundamental frequency of the system.

The complex normal form method proposed by Nayfeh [14] is a powerful tool to study nonlinear dynamical problems. Leung and Zhang [15] studied a Duffing oscillator’s free vibration applying this method. Then Wang and Zhang [16–21] developed this method to deal with bifurcation problems and to detect threshold parameters leading to chaotic motions. Besides that, He’s energy balance method is a powerful and convenient tool in detecting the strongly nonlinear oscillators’ frequency [22]. The classic feature of this method is constructing the system’s Hamilton function and keeping the input and output energy in balance. The striking advantage of this method is that the actual frequency does not need to be expanded into series. In other words, the steady state frequency could be solved as a single variable. Thus, the amount of calculation could be reduced significantly. The final expression of the frequency is expressed as a function containing the steady state amplitude and the nonlinear parameters. Even if the nonlinearity is driven remarkably large, the error of the approximate frequency is small enough to be neglected. In the previous study, the nonlinear terms together with the external exciting forces and parametrical exciting forces are considered as the reason of changing the frequency of system. As far as we know, studies about the complex normal form method and He’s energy method applied on stochastic nonlinear problems are few.

Although these improvements mentioned above are effective in finding frequency of strongly nonlinear systems, these studies are less than perfect. It should be noticed that in these studies the vibration phase has all been ignored. Furthermore, the accuracy of these methods essentially depends on the choice of initial values. If the initial value is chosen just near the actual amplitude the methods will offer precise enough results. But if the initial value is chosen unluckily far from the first-order amplitude, the results could be totally unacceptable. Therefore, there remains a demand for an approach to calculating the frequency, amplitude, and phase of system’s responses simultaneously.

The purpose of this paper we present is attempting to extend the complex normal form method to strongly nonlinear stochastic filed. As an example, a Duffing oscillator under narrow band noise excitation was studied in detail. Firstly, a new expression of the response was introduced to improve the performance of Wang’s method [16]. Subsequently, He’s energy balance method was applied to offer an extra equation which could be combined to the results obtained by the normal form method. Consequently, the accuracy of the frequency could be considerably improved. Thus, both frequency and phase of the response were obtained at the same time. Based on the obtained frequency, the perturbation and the moment method were used to estimate the statistical character of the response. Finally, numerical simulations were carried out to verify the accuracy of the theoretical analysis.

#### 2. The Complex Normal Form Method

We take a Duffing oscillator with nonlinear damping subject to external narrow band noise as an example which is as follows:where denotes the free frequency, denotes the nonlinear Stiffness coefficient, , are the linear and nonlinear damping coefficients, denotes the amplitude of the exciting force, denotes the center frequency of the exciting force, denotes the standard Wiener process, denotes the small strength of the , and the parameter is not small, so we call the system “strongly nonlinear.” Such type of equation is widely used in engineering fields such as vibrating beams and plates with nonlinear damping.

In this paper, we focus on the primary resonance. The actual approximate undetermined frequency is set as and the frequency of the external exciting force equals . At first, the complex variables and are introduced to rewrite the system’s variables. Thus, the system’s displacement and velocity can be expressed aswhere is unknown complex variable to be determined and overbar denotes the conjugation.

The external force can be rewritten aswhere and . Solving (2) for , yields Taking differentiation to the first function in (4) with respect to time and considering (1) and (3), one obtains

The basic issue of the complex normal form method is introducing a series of nearly identified transformation. To simplify (5), we introduce a near-identity transformation from to in the formwherewhere are coefficients to be determined. For example,

Substituting (6) and (7) into (5) yields a new equation for . Terms proportional to , , , , ,, , are resonance terms and hence cannot be eliminated while all the remaining terms should be eliminated. By eliminating the nonresonance terms of for , the simplest possible form gives

By now, similar work has been presented by Zhang and Wang [16]. The key technique of our work is using a new expression of approximate solution instead of the expression given in [16–21]. Considering the first-order guess of system’s response must include both frequency and phase variation, the approximate solution of system (1) with could be expressed asThus and could be expressed as follows: , , where + .Another expression of (10) iswhere , .

One should notice that the first-order guess of solution expressed as (11) is different from Zhang [12]. In that study, the first-order guess is given as , . The new expression ensures the phase existing.

Thus, our job is calculating the relationship of the three undetermined variables , , and .

By substituting (10) and (3) into (9), one obtainsMultiplying both sides of (12) with and separating the real and image parts yieldUntil then, the full expressions of and are proposed in (13). Theoretically, the steady state solution can be obtained by letting .

By now the method has not been finished. As we know, there are three unknown variables , , and to be determined. But we only have two equations. Chen [1, 2] offered an approach to solve this problem. With the initial displacement value given previously, the variables can be calculated by (13) subsequently. However, using this approach has limitation to some extent. As mentioned in Introduction, if the initial displacement value is chosen correctly, the solution will be precise enough. But if the initial displacement is not given properly, (13) may offer nonreal roots or totally unacceptable roots. Thus, one extra equation is needed to solve unknown variables , , and simultaneously.

#### 3. He’s Energy Balance Method

The so-called energy balance method proposed by He [22] has been used and developed by many researchers to determine frequency of nonlinear oscillations in view of the fact that it is very effective and convenient. The striking issue of He’s method is building a Hamilton and then making the input and output of it kept in balance. The validity has been verified in previous studies [23–27]. Thus, we take this approach to build an extra equation so that (13) can be solved.

The Hamilton of system (1) is given as follows:Substituting the assumed first-order approximate solution (10) into (14) and letting and , respectively, yieldThe residual can be expressed as . Then, based on Galerkin-Petrov method [28], one equation is obtained aswhere means the period.

Combing (17) with (13), three equations can be solved numerically to determine the three approximate variables , , and .

#### 4. Perturbation Procedures

By introducing a detuning parameter into , we attempt to investigate how the response of the noise included strongly nonlinear system (1) changes by the detuning parameter . Considering the parameter transformation introduced in [9, 10] and the undetermined frequency expressed as , the exciting frequency is expressed asA new form of is given aswhich is different from (10).

Substituting (3) and (19) into (9) and separating the real and image parts yield

Substituting (18) into (20), expanding (20) in Taylor series with respect to , and then ignoring the high order terms, one obtainsIntroducing a time scale transformation , (21) can be rewritten aswhere , as we know for Wiener process which is a mathematical tool to describe famous Brown’s motion; there exists , , where denotes the mathematical expectation.

Rewrite (22) as the following Ito type equation, with the parameters , , and :

The first-order steady state approximate response of noise-free system (1) is expressed as (26), where are determined by (23). When and , such thatsquaring and adding both sides of equations in (24), one obtains the frequency-amplitude response function

#### 5. The Moment Method

Next, we want to explore the probability density of the responses. Then we discuss the effect of the noise. The response of the system is determined by the steady state probability density function. As the Fokker-Planck-Kolmogorov equation for the response probability density function (PDF) is not analytically solvable, the stochastic moment method [29] is considered alternatively. The moment method is an acceptable backup to study the character of the responses approximately. We setwhere , are the small perturbation terms; substituting the above equations into (24) and ignoring the high order terms, the linearization equation is expressed aswhere denotes the standard Gauss white noise.

Taking expectation on both sides of (27) and one obtainsThus, it is obvious that , .

Based on the Floquet theory, the Jacoby matrix is used to judge the stability of the first-order moment, which is expressed as follows:The characteristic equation can be expressed as , whereThe roots of (30) areObviously, the stable condition of (27) without noise is , , which means the solution is realizable by numerical simulation.

Then the moment method is applied to study the second-order moment of the system’s responses. Similarly, we have the obvious results . Applying Ito’s rule [30] of stochastic differentiation, one obtainsThe solutions of (32) areOne obtainsForm (33), it is easy to see that, for , , there must beFunction (35) also means the solutions are stable if the inequality is satisfied (see (30)).

When the noise density is small enough, the response’s probability density function can be approximately expressed as the Gauss distribution:

#### 6. Results and Numerical Simulation

Let the coefficients be , , , , , and . Firstly, we study the noise-free system. Combing (13) and (17), the undetermined variables were obtained as , , , and . Thus the amplitude was calculated as . To compare the accuracy of this method and the improved LP method, the results were displayed as follows, respectively. With the initial displacement chosen as = −0.02, the variables and were given by the improved LP method as , , and . As seen in Figure 1, it is obvious that the accuracy of the improved complex normal form method is much better than the improved LP method. In Figure 1(a) one can see that the actual result and the approximate result are so close to each other that the difference can hardly be found. On the other hand, the result given by the improved LP method is slightly different from the actual result. Furthermore, as mentioned in Introduction, the accuracy of the improved LP method depends on the initial values. If the initial value is chosen randomly, the improved LP method may offer nonreal results.