Abstract and Applied Analysis

Volume 2013 (2013), Article ID 813957, 6 pages

http://dx.doi.org/10.1155/2013/813957

## The Application of the Undetermined Fundamental Frequency Method on the Period-Doubling Bifurcation of the 3D Nonlinear System

^{1}School of Mechanical Engineering, Tianjin Polytechnic University, Tianjin 300072, China^{2}Department of Mechanics, Tianjin University, Tianjin 300072, China

Received 28 April 2013; Revised 10 August 2013; Accepted 10 August 2013

Academic Editor: Ren Yong

Copyright © 2013 Gen Ge and Wei 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 analytical method to predict the period-doubling bifurcation of the three-dimensional (3D) system is improved by using the undetermined fundamental frequency method. We compute the stable response of the system subject to the quadratic and cubic nonlinearity by introducing the undetermined fundamental frequency. For the occurrence of the first and second period-doubling bifurcation, the new bifurcation criterion is accomplished. It depends on the stability of the limit cycle on the central manifold. The explicit applications show that the new results coincide with the results of the numerical simulation as compared with the initial methods.

#### 1. Introduction

Period-doubling bifurcation can induce complex dynamical behavior in the nonlinear dynamic systems. It is the most classical achievement broadcasted by Feigenbaum [1]. He discovered the ratio of the difference between the values at which such successive period-doubling bifurcation occurs at a constant of around 4.6692 and then showed that the same behavior, with the same mathematical constant, would occur within a wide class of mathematical functions, prior to the onset of chaos. Since then many attempts have been made to study the period-doubling (flip) bifurcation phenomenon in the nonlinear dynamic systems. Wang and Xu [2] developed the relation between two periodic solutions analytically for a general parameter dependent dynamic system. Such relation is further confirmed by one example and shows that the 2T-periodic solution contains all the information of the T-periodic solution near the bifurcation point. From the frequency domain point of view, Floquet multipliers commonly used for the analytical bifurcations of Hopf cycles are the key to detect the appearance of a subharmonic solution. So a quasianalytical monodromy matrix approach was developed to the period-doubling bifurcation emerging near a Hopf bifurcation point [3].

As compared with the single freedom system, the dynamical behaviors surrounding the bifurcation point may become more complicated in the 3D system. Rand [4] used the center manifold theory to approximate the newly born limit cycle and then to investigate the stability of the limit cycle corresponding to the flip bifurcation. Later, Belhaq et al. [5, 6] improved the approximation of the critical value with a higher-order approximation and further solved the problem of the second period-doubling bifurcation.

In this paper, we use the center manifold theory to reduce a 3D system and then derive the critical values of the first and second period-doubling bifurcation according to the stability of the limit cycle. In terms of the undetermined fundamental frequency method, it produces more accurate results and avoids the computational complexity appending the higher-order approximation at the same time [7, 8]. Finally the whole process is precisely programmed in terms of the computer algebra Mathematica to perform the analysis more efficiently.

#### 2. Stable Response with the Undetermined Fundamental Frequency Method

In order to illustrate the analytical process, we refer to the following 3D system:

This system may be thought of as a feedback control system consisting of a damped linear oscillator in the , variables and a control variable . The origin is the equilibrium and may lose its stability at control parameter changing from to . This means that the period-doubling bifurcation appears at the value following the limit cycles.

For the value of , the center manifold theory has to be introduced to finish the reduction and obtain the equations on the center manifolds. So we set the second order polynomial of in terms of , , and

Differentiating (2) with respect to time and using (1) give Equating the same order terms on both sides of (3) produces the coefficients That leads to the following approximate flows on the center manifold:

The computational precision of the critical value depends heavily on the stable response, such as the frequency and amplitude of the 3D system. So, in order to perform the limit cycle bifurcation analysis more correctly, Belhaq et al. [6] explored the analysis through a higher-order approximation. In this paper we introduce the undetermined fundamental frequency method during the course of normal form operation.

In terms of the transformation , , (5) changes to To obtain the stable response, it demands to transform (6) into a differential equation of the first order with the complex unknown quantities . Let where is the undetermined fundamental frequency. Solving (7) obtains

Differentiating (8) with respect to and using (6) and (7), give For the simplification of (9), a third-order nonlinear transformation is considered as where , .

These transformation coefficients are suitably chosen to eliminate the nonresonance terms [9] in the final expression. So the normal form of (9) is Next, , require to be expressed in the following polar form , . Then separating the real and imaginary parts of the foregoing equation by considering the stationary condition , we have Hence, (12) produces the amplitude and the undetermined fundamental frequency

#### 3. Criterion for the Period-Doubling Bifurcation

Substituting (13) into (7), we obtain the expression of periodic solution in the trigonometric form so that the solution changes into the Cartesian form with the transformation , . Consider

Note that the limit cycle cannot show period-doubling as long as it lies in the center manifold because the latter is two-dimensional and trajectories cannot self-intersect. So it marks on the limit cycle for the expression obtained from (2). To investigate the stability, we append disturbance to in (1), that is, and linearize the variation The general solution of (16) is where . By considering the Floquet theory, the transition from stable to unstable occurs in the condition of , where is the period of the limit cycle oscillation. That produces the critical value of period-doubling bifurcation through , and the result is

In order to illustrate the accuracy of the result, the critical values obtained from different methods are presented in Table 1.

Finally, we investigate the stability to find the critical values of the second period doubling bifurcation. That is, from a bifurcation point, the asymmetric 2T-orbit born at the first period-doubling bifurcation point becomes nonstable in a flip bifurcation, where a 4T-orbit emerges. Here, we do not want to give too many details about the second period-doubling bifurcation because the computational process is very similar to the first period doubling bifurcation. We mainly follow the stability analyses of Rand [4] and Belhaq et al. [6]. The main difference exists in finding the solution of the 3D system, where we use the undetermined fundamental frequency method. It produces a better approximation of the asymptotical solution. First of all, we give the general solution in a complex form as follows:

Then, we may find the normal form of the reduced system. As for the secular terms, they can be regarded as the near resonance according to [9]. We use a near identity transformation from to , to find the normal form of the system which also includes the subharmonic components in its expression. That can be written in different order: first order , second order , , and third order: , , , , . With these secular terms we find the averaged equation of the system and then amplitude and frequency. Finally, we investigate the stability to find the critical values of the second period-doubling bifurcation. The critical values are presented in Table 2. It exhibits a better approximation than the high order analysis given by Belhaq et al. [6].

In Figure 1, the projections on , , and plane are plotted at the value of . Meanwhile the time series of trajectories , , and appear in Figure 2. A 3D phase portrait refers to the first and second period-doubling bifurcations that are portrayed at the values equal to 0.439 and 0.476 in Figure 3, respectively. Finally we programme the whole computation process in terms of the computer algebra Mathematica to accelerate the analysis more efficiently.

#### 4. Conclusion

The strategy of predicting the period-doubling bifurcation of the 3D system is presented by using the undetermined fundamental frequency method. It applies the undetermined fundamental frequency to obtain the stable response of the flows on the center manifold and then forms the criterion of period-doubling prediction by considering the stability of the limit cycle. In contrast to the result of numerical simulation, it reveals a good prediction as shown in Tables 1 and 2, compared with the analytical results of the first and second period-doubling bifurcations given by Rand and Belhaq. The whole process is constituted in terms of the computer algebra Mathematica. It enables people to research the flip bifurcation of the 3D system more accurately and efficiently.

The strategy presented in this work is sufficiently general, so it would be possible to apply the present method to consider other high-dimensional and more complicated systems, which will be the topics for further research.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 11272229) and Science Foundation of Tianjin Education Committee (Grant no. 20120902).

#### References

- M. J. Feigenbaum, “The universal metric properties of nonlinear transformations,”
*Journal of Statistical Physics*, vol. 21, no. 6, pp. 669–706, 1979. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Q. Wang and M. T. Xu, “Property of period-doubling bifurcations,”
*Chaos, Solitons and Fractals*, vol. 24, no. 2, pp. 527–532, 2005. View at Publisher · View at Google Scholar · View at Scopus - D. W. Berns, J. L. Moiola, and G. R. Chen, “Detecting period-doubling bifurcation: an approximate monodromy matrix approach,”
*Automatica*, vol. 37, no. 11, pp. 1787–1795, 2001. View at Publisher · View at Google Scholar · View at Scopus - R. H. Rand, “Analytical approximation for period-doubling following a Hopf bifurcation,”
*Mechanics Research Communications*, vol. 16, no. 2, pp. 117–123, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Belhaq and M. Houssni, “Symmetry-breaking and first period-doubling following a Hopf bifurcation in a three-dimensional system,”
*Mechanics Research Communications*, vol. 22, no. 3, pp. 221–231, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Belhaq, E. Freire, M. Houssni, and A. J. Rodríguez-Luis, “Second period-doubling in a three-dimensional system,”
*Mechanics Research Communications*, vol. 26, no. 2, pp. 123–128, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Y. T. Leung and Q. C. Zhang, “Complex normal form for strongly non-linear vibration systems exemplified by Duffing-van der Pol equation,”
*Journal of Sound and Vibration*, vol. 213, no. 5, pp. 907–914, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. C. Zhang, W. Wang, and W. Y. Li, “Heteroclinic bifurcation of strongly nonlinear oscillator,”
*Chinese Physics Letters*, vol. 25, no. 5, p. 1905, 2008. View at Publisher · View at Google Scholar - A. H. Nayfeh,
*Method of Normal Forms*, John Wiley, New York, NY, USA, 1993. View at MathSciNet