Mathematical Problems in Engineering

Volume 2010 (2010), Article ID 986319, 15 pages

http://dx.doi.org/10.1155/2010/986319

## Periodic and Chaotic Motions of a Two-Bar Linkage with OPCL Controller

^{1}School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110004, China^{2}School of Mechanical Engineering, University of Birmingham, Birmingham B15 2TT, UK

Received 16 December 2009; Accepted 25 June 2010

Academic Editor: Irina N. Trendafilova

Copyright © 2010 Qingkai Han et al. 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 two-bar linkage, which is described in differential dynamical equations, can perform nonlinear behaviors due to system parameters or external input. As a basic component of robot system, the investigation of its behavior can improve robot performance, control strategy, and system parameters. An open-plus-close-loop (OPCL) control method therefore is developed and applied to reveal and classify the complicated behaviors of a two-bar linkage. In this paper, the conception and stability of OPCL are addressed firstly. Then it is applied to the dynamical equations of two-bar linkage. Different motions including single-periodic, multiple-periodic, quasiperiodic, and chaotic motions are unfolded by numerical simulations when changing the controller parameters. Furthermore, the obtained chaotic motions are sorted out for qualitative and quantificational study using Lyapunov exponents and hypothetic possibilities of surrogate data method.

#### 1. Introduction

A two-bar linkage, as a basic component of mechanical system, can perform nonlinear motions, among which chaotic motion is the most typically complicated one. It is known that the performance of such a mechanism is influenced by system parameters such as mass and friction coefficient, the initial states, and external input such as driving torque controlled by specially designed controller. Conversely, study on these motions can provide a novel way to improve system’s performance, optimize structure design, and develop new control strategy.

For a two-bar linkage, the motions of its two rotating links can be single periodic, multiple periodic, quasiperiodic, and chaotic. In the past decades, many control strategies were explored to obtain certain motions of a two-bar linkage mechanism. A neural controller was employed to achieve two typical synchronous motions, that is, giant rotating motion and small swing motion as stated in [1, 2]. A PD feedback controller to obtain chaotic motions for a two-bar linkage was reported in [3, 4]. The bifurcation characteristics of motions of a two-bar linkage were presented in [5]. The transferring process of a two-bar linkage from period-doubling bifurcation to chaotic motions was studied by changing control variables in [6]. In recent research, it is known that the open-plus-close-loop (OPCL) control strategy is powerful for complicated dynamic systems, as stated in [7, 8], which has been used for chaotic control of chaotic systems [9] and synchronous systems [10, 11]. Besides, the parametric OPCL method [12] and the nonlinear OPCL method were also applied for motion control of some dynamic systems [13].

In this paper, an OPCL controller is proposed for a two-bar linkage to achieve different motions as stated above. The dynamics of the two-bar linkage is modeled with nonlinear differential equations, followed by the Lyapunov stability analysis of the controlled system. By changing the two coefficient matrices **A** and **B** of the OPCL controller, we obtain typical motions of the two-bar linkage in numerical simulations including single-periodic, multiple-periodic, quasiperiodic, as well as chaotic motions. These different motions are described in both qualitative and quantificational ways, such as phase-space portraits, frequency spectra, Lyapunov exponents, and the hypothetic possibilities of surrogate data.

#### 2. Dynamical Equations of a Two-Bar Linkage with OPCL Controller

##### 2.1. A Dynamical System with an OPCL Controller

Let a typical dynamical system be as follows:
where is the state variable, and *n* is the number of DOF of the system.

The given motion goal of the system is

Define a tracking error of the two-bar linkage as

Equation (2.1) is linearized and expanded in the neighborhood of the goal via Taylor series and it becomes where and are Jacobian matrices of with respect to and , respectively.

An OPCL controller for the system is designed as [8]
where the term of is the open-loop part, and the term of is the closed-loop part. In this controller, the coefficient matrices **A** and **B** are assumed to be diagonal.

The controlled system, that is, the dynamical system with its OPCL controller, is re-written as follows: The tracking error equation defined as (2.3) can be deduced as follows with omitting higher-order term:

It is noticed that, if the controlled system of (2.6) is asymptotically stable when coefficient matrices **A **and **B** are constant and with negative real parts of eigenvalues, that is, their elements and should be negative, (2.7) can be rewritten in an expanded form as

Thus, if () is asymptotically stable, the system is also stable. Therefore, the aforementioned assumption is true.

A Lyapunov function *V *is defined as
Then

According to the Lyapunov stability theory, it is proved that the error equation of (2.7) is asymptotically stable if the real parts of the eigenvalues of the coefficient matrices **A** and **B** are negative. It can be found that such an OPCL controller is reliable enough to be applied on a two-bar linkage in order to get different motions as expected and these motions are asymptotically stable.

##### 2.2. Dynamic Equations of Two-Bar Linkage

A two-bar linkage is shown in Figure 1. Its two joints can be driven, respectively, at and , so that the upper link (Link 1) and the lower link (Link 2) can rotate, respectively, around their own joints of and within the range of [-,]. The upper joint is fixed on the ground.

In coordinate system *O _{1}xy*, as shown in Figure 1, is the rotating angle of Link 1 with respect to

*y-*axis, and is the rotating angle of Link 2 with respect to the centerline of Link 1. The dynamical equations of the two-bar linkage are given as follows: where

In (2.11), is the mass of Link *i*, ; is the moment of inertia of Link *i *with respect to its mass center, ; is the distance between Link *i* and joint* i*, ; *g *is the acceleration of gravity; is driving moment in joint *i, *.

Equation (2.11) can be rewritten in the following form, which is also in the form of (2.1): where

#### 3. Simulations on Different Motions of Two-Bar Linkage

##### 3.1. Conditions and Simulation Steps

In numerical simulations, the dimensionless structural parameters of the two-bar linkage are used and listed in Table 1.

The goal trajectories of the two rotating angles of the two links, and , are designed as

The simulation steps, mainly referring to the system of (2.6), are listed herein.(1)Set the initial states and the structural parameter values of the two-bar linkage and the total number of simulation *n*.(2)Set the control parameter matrices and of the OPCL controller of (2.5).(3)Set the goal trajectory of (3.1).(4)Calculate the joint angles based on (2.6) using the fourth-order Runge-Kutta method.

With different values of **A **and** B** of the OPCL controller, the two-bar linkage can achieve different motions including single-periodic, multiple-periodic, quasiperiodic, and chaotic motions. The typical parameter values of **A** and **B** and their corresponding motions are listed in Table 2. They are illustrated in the following sections, which are also listed in the fourth column of Table 2.

The obtained motions of link joints can be investigated qualitatively and quantificationally in different ways, including observation of frequency spectra and estimation of nonlinear parameters. Single- or multiamplitude lines in frequency spectra indicate a periodic motion or a multiple-periodic motion. Wide-range frequency distributions are often generated from quasiperiodic motion and chaotic motion. In Poincare mapping plots, one single point or some scattered points are often mapped by periodic motions whereas attractors with concentrated area and so-called strange attractors are mostly from either quasiperiodic or chaotic motions. In addition, for chaotic motions, there are also some critical methods to judge, for example, the well-known nonlinear parameter estimation of positive maximum Lyapunov exponent [14, 15] and the checking possibility of surrogate data method [16] which is powerful to distinguish a chaotic motion from a random one.

##### 3.2. Single-Periodic Motions

Given and , the two-bar linkage can achieve single-periodic motions. The simulated motions with initial conditions of , , , and are shown in Figure 2.

In Figures 2(a) and 2(b), the motions of and are obviously periodic, that is, harmonic. The two phase plane portraits of and , shown in Figures 2(c) and 2(d), are closed curves, which prove that the motions of the two-bar linkage are stable. Furthermore, the Poincare mapping portrait of each has only one isolated point, as shown in Figures 2(e) or 2(f). The frequency spectra of the two rotating angles show that there exists only single dominant frequency of about 0.16 Hz, as shown in Figures 2(g) and 2(h). The simulations can also show that periodic motions of the two-bar linkage are stable in this case.

##### 3.3. Multiple-Periodic Motions

Given and , the two-bar linkage can realize multiple-periodic motions. The typical simulated motions of the two rotating angles are shown in Figure 3 with initial conditions of , , , and .

In Figures 3(a) and 3(b), the motions of and are also periodical but with other harmonic components. The two phase plane portraits of and , shown in Figures 3(c) and 3(d), are closed curves, which can indicate that the motions of the two-bar linkage are stable. Accordingly, the Poincare mapping portrait of each, shown in Figure 3(e) or 3(f), has some isolated points. There also exist more than two and four obvious frequency lines, as shown in Figures 3(g) and 3(h), which show that the multiperiodic motions of the two-bar linkage are stable in this case.

##### 3.4. Quasiperiodic Motions

Given and , the two-bar linkage can realize quasiperiodic motions. The simulated motions of the two joints are shown in Figure 4 with the initial conditions of , , , and .

It can be seen from Figure 4 that the motions of and are complex. Both the time histories and the phase plane portraits of and are difficult to distinguish the motion type rather than traditional harmonics. The Poincare mapping portrait of each, shown in Figures 4(e) and 4(f), is the concentrated stick-like area. More than six and eight frequency lines appear in the corresponding frequency spectra of and shown in Figures 4(g) and 4(h). The simulations can also show that the motions of the two-bar linkage are quasiperiodic in this case.

##### 3.5. Chaotic Motions

Given and , the two-bar linkage can realize chaos motions. The simulated chaotic motions are shown in Figure 5, in the case of the initial conditions of , , , and .

From Figures 5(a) and 5(b), the simulated responses of and are irregular without obvious periods. The two phase plane portraits (shown in Figures 5(c) and 5(d)) and the Poincare mapping (shown in Figures 5(e) and 5(f)) of and illustrate irregular shape or strange attractors. Their corresponding amplitude spectra also unfold multifrequency lines and explicit broadband ranges. According to the qualitative theory of chaos, these motions are chaotic in this case.

According to the above simulation results, even for the same initial conditions, controlled behaviors vary with the parameters in the controllers. This is because the basins of entrainment, whose counterparts in maps are addressed in [9, 10], depend on the parameters of controller.

#### 4. Discussion on the Simulated Chaotic Motions

In order to quantitatively describe the simulated chaotic motions of the two-bar linkage, Lyapunov exponent and the hypothesis possibility with surrogate data method are used based on nonlinear theory [14–17]. The first five order Lyapunov exponents of the time series of and are calculated according to algorithm in [14] and shown in Table 3. The calculated hypothesis possibilities with the surrogate data method [17] for the time series of and are shown in Table 4.

As shown in Table 3, the maximum Lyapunov exponents of the motions of angle and are 0.5496 and 0.1431, respectively. The positive values indicate that the motions of the two angles are chaotic.

From the calculated hypothesis possibilities with the surrogate data method for the motions of and , that is, and , it is demonstrated that both of them are chaotic due to the checking possibility values which are smaller than .05, referring to an empirical value in [17].

#### 5. Conclusions

In this paper, the dynamical model of a two-bar linkage with OPCL controller is proposed in order to obtain different motions. It is verified that the OPCL controlled system is asymptotically stable based on the Lyapunov theory, when the control coefficient matrices of **A** and **B **are diagonal and with negative real parts of eigenvalues. It is reliable to force a two-bar linkage to achieve different stable motions of and .

Numeral simulations are conducted to demonstrate that the two-bar linkage can achieve single-periodic, multiple-periodic, quasiperiodic, and chaotic motions by changing the control parameters of **A** and **B** of the OPCL controller for and successfully. The proposed OPCL approach works on the given conditions.

Furthermore, for the simulated chaotic motions, the maximum Lyapunov exponents are positive, that is, 0.5496 and 0.1431, respectively. The calculated hypothesis possibilities of them with the surrogate data method for the same chaotic motions are and , smaller than .05.

#### Acknowledgment

The authors gratefully acknowledge that the work was supported by the Key Project of Science and Technology Research Funds of Chinese Ministry of Education (Grant no. 108037).

#### References

- K. Matsuoka, N. Ohyama, A. Watanabe, and M. Ooshima, “Control of a giant swing robot using a neural oscillator,” in
*Proceedings of the 1st International Conference on Natural Computation (ICNC '05)*, vol. 3611 of*Lecture Notes in Computer Science*, pp. 274–282, August 2005. View at Scopus - Q. Han, Z. Qin, X. Yang, and B. Wen, “Rhythmic swing motions of a two-link robot with a neural controller,”
*International Journal of Innovative Computing, Information and Control*, vol. 3, no. 2, pp. 335–342, 2007. View at Google Scholar · View at Scopus - S. Lankalapalli and A. Ghosal, “Possible chaotic motions in a feedback controlled 2R robot,” in
*Proceedings of the 13th IEEE International Conference on Robotics and Automation*, pp. 1241–1246, IEEE, April 1996. View at Scopus - A. S. Ravishankar and A. Ghosal, “Nonlinear dynamics and chaotic motions in feedback-controlled two- and three-degree-of-freedom robots,”
*International Journal of Robotics Research*, vol. 18, no. 1, pp. 93–108, 1999. View at Google Scholar · View at Scopus - F. Verduzco and J. Alvarez, “Bifurcation analysis of A 2-DOF robot manipulator driven by constant torques,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 9, no. 4, pp. 617–627, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - K. Li, L. Li , and Y. Chen, “Chaotic motion of a planar 2-dof robot,”
*Machine*, vol. 29, no. 1, pp. 6–8, 2002 (Chinese). View at Google Scholar - E. A. Jackson and I. Grosu, “An open-plus-closed-loop (OPCL) control of complex dynamic systems,”
*Physica D*, vol. 85, no. 1-2, pp. 1–9, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L.-Q. Chen and Y.-Z. Liu, “A modified open-plus-closed-loop approach to control chaos in nonlinear oscillations,”
*Physics Letters A*, vol. 245, no. 1-2, pp. 87–90, 1998. View at Publisher · View at Google Scholar · View at Scopus - L.-Q. Chen, “An open-plus-closed-loop control for discrete chaos and hyperchaos,”
*Physics Letters A*, vol. 281, no. 5-6, pp. 327–333, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - L.-Q. Chen and Y.-Z. Liu, “An open-plus-closed-loop approach to synchronization of chaotic and hyperchaotic maps,”
*International Journal of Bifurcation and Chaos*, vol. 12, no. 5, pp. 1219–1225, 2002. View at Publisher · View at Google Scholar · View at Scopus - Q.-K. Han, X.-Y. Zhao, and B.-C. Wen, “Synchronization motions of a two-link mechanism with an improved OPCL method,”
*Applied Mathematics and Mechanics*, vol. 29, no. 12, pp. 1561–1568, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - L.-Q. Chen, “The parametric open-plus-closed-loop control of chaotic maps and its robustness,”
*Chaos, Solitons and Fractals*, vol. 21, no. 1, pp. 113–118, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - Y.-C. Tian, M. O. Tadé, and J. Tang, “Nonlinear open-plus-closed-loop (NOPCL) control of dynamic systems,”
*Chaos, solitons and fractals*, vol. 11, no. 7, pp. 1029–1035, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, “Determining Lyapunov exponents from a time series,”
*Physica D*, vol. 16, no. 3, pp. 285–317, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. A. Tsonis and J. B. Elsner, “Nonlinear prediction as a way of distinguishing chaos from random fractal sequences,”
*Nature*, vol. 358, no. 6383, pp. 217–220, 1992. View at Publisher · View at Google Scholar · View at Scopus - M. T. Rosenstein, J. J. Collins, and C. J. De Luca, “A practical method for calculating largest Lyapunov exponents from small data sets,”
*Physica D*, vol. 65, no. 1-2, pp. 117–134, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. Theiler, S. Eubank, A. Longtin, B. Galdrikian, and J. Doyne Farmer, “Testing for nonlinearity in time series: the method of surrogate data,”
*Physica D*, vol. 58, no. 1–4, pp. 77–94, 1992. View at Publisher · View at Google Scholar · View at Scopus