Research Article  Open Access
Qingkai Han, Xueyan Zhao, Xiaoguang Yang, Bangchun Wen, "Periodic and Chaotic Motions of a TwoBar Linkage with OPCL Controller", Mathematical Problems in Engineering, vol. 2010, Article ID 986319, 15 pages, 2010. https://doi.org/10.1155/2010/986319
Periodic and Chaotic Motions of a TwoBar Linkage with OPCL Controller
Abstract
A twobar 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 openpluscloseloop (OPCL) control method therefore is developed and applied to reveal and classify the complicated behaviors of a twobar linkage. In this paper, the conception and stability of OPCL are addressed firstly. Then it is applied to the dynamical equations of twobar linkage. Different motions including singleperiodic, multipleperiodic, 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 twobar 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 twobar 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 twobar 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 twobar linkage was reported in [3, 4]. The bifurcation characteristics of motions of a twobar linkage were presented in [5]. The transferring process of a twobar linkage from perioddoubling bifurcation to chaotic motions was studied by changing control variables in [6]. In recent research, it is known that the openpluscloseloop (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 twobar linkage to achieve different motions as stated above. The dynamics of the twobar 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 twobar linkage in numerical simulations including singleperiodic, multipleperiodic, quasiperiodic, as well as chaotic motions. These different motions are described in both qualitative and quantificational ways, such as phasespace portraits, frequency spectra, Lyapunov exponents, and the hypothetic possibilities of surrogate data.
2. Dynamical Equations of a TwoBar 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 twobar 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 openloop part, and the term of is the closedloop 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 rewritten as follows: The tracking error equation defined as (2.3) can be deduced as follows with omitting higherorder 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 twobar linkage in order to get different motions as expected and these motions are asymptotically stable.
2.2. Dynamic Equations of TwoBar Linkage
A twobar 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 yaxis, and is the rotating angle of Link 2 with respect to the centerline of Link 1. The dynamical equations of the twobar 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 TwoBar Linkage
3.1. Conditions and Simulation Steps
In numerical simulations, the dimensionless structural parameters of the twobar 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 twobar 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 fourthorder RungeKutta method.
With different values of A and B of the OPCL controller, the twobar linkage can achieve different motions including singleperiodic, multipleperiodic, 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 multipleperiodic motion. Widerange 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 socalled 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 wellknown 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. SinglePeriodic Motions
Given and , the twobar linkage can achieve singleperiodic motions. The simulated motions with initial conditions of , , , and are shown in Figure 2.
(a) Time history of
(b) Time history of
(c) Phaseplane portrait of
(d) Phaseplane portrait of
(e) Poincare mapping of
(f) Poincare mapping of
(g) Amplitude spectrum of
(h) Amplitude spectrum of
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 twobar 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 twobar linkage are stable in this case.
3.3. MultiplePeriodic Motions
Given and , the twobar linkage can realize multipleperiodic motions. The typical simulated motions of the two rotating angles are shown in Figure 3 with initial conditions of , , , and .
(a) Time history of
(b) Time history of
(c) Phaseplane portrait of
(d) Phaseplane portrait of
(e) Poincare mapping of
(f) Poincare mapping of
(g) Amplitude spectrum of
(h) Amplitude spectrum of
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 twobar 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 twobar linkage are stable in this case.
3.4. Quasiperiodic Motions
Given and , the twobar linkage can realize quasiperiodic motions. The simulated motions of the two joints are shown in Figure 4 with the initial conditions of , , , and .
(a) Time history of
(b) Time history of
(c) Phaseplane portrait of
(d) Phaseplane portrait of
(e) Poincare mapping of
(f) Poincare mapping of
(g) Amplitude spectrum of
(h) Amplitude spectrum of
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 sticklike 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 twobar linkage are quasiperiodic in this case.
3.5. Chaotic Motions
Given and , the twobar linkage can realize chaos motions. The simulated chaotic motions are shown in Figure 5, in the case of the initial conditions of , , , and .
(a) Time history of
(b) Time history of
(c) Phaseplane portrait of
(d) Phaseplane portrait of
(e) Poincare mapping of
(f) Poincare mapping of
(g) Amplitude spectrum of
(h) Amplitude spectrum of
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 twobar 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 twobar 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 twobar linkage to achieve different stable motions of and .
Numeral simulations are conducted to demonstrate that the twobar linkage can achieve singleperiodic, multipleperiodic, 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: Google Scholar
 Q. Han, Z. Qin, X. Yang, and B. Wen, “Rhythmic swing motions of a twolink 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
 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: Google Scholar
 A. S. Ravishankar and A. Ghosal, “Nonlinear dynamics and chaotic motions in feedbackcontrolled two and threedegreeoffreedom robots,” International Journal of Robotics Research, vol. 18, no. 1, pp. 93–108, 1999. View at: Google Scholar
 F. Verduzco and J. Alvarez, “Bifurcation analysis of A 2DOF 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 Site  Google Scholar  Zentralblatt MATH
 K. Li, L. Li , and Y. Chen, “Chaotic motion of a planar 2dof robot,” Machine, vol. 29, no. 1, pp. 6–8, 2002 (Chinese). View at: Google Scholar
 E. A. Jackson and I. Grosu, “An openplusclosedloop (OPCL) control of complex dynamic systems,” Physica D, vol. 85, no. 12, pp. 1–9, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L.Q. Chen and Y.Z. Liu, “A modified openplusclosedloop approach to control chaos in nonlinear oscillations,” Physics Letters A, vol. 245, no. 12, pp. 87–90, 1998. View at: Publisher Site  Google Scholar
 L.Q. Chen, “An openplusclosedloop control for discrete chaos and hyperchaos,” Physics Letters A, vol. 281, no. 56, pp. 327–333, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L.Q. Chen and Y.Z. Liu, “An openplusclosedloop 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 Site  Google Scholar
 Q.K. Han, X.Y. Zhao, and B.C. Wen, “Synchronization motions of a twolink mechanism with an improved OPCL method,” Applied Mathematics and Mechanics, vol. 29, no. 12, pp. 1561–1568, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L.Q. Chen, “The parametric openplusclosedloop control of chaotic maps and its robustness,” Chaos, Solitons and Fractals, vol. 21, no. 1, pp. 113–118, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y.C. Tian, M. O. Tadé, and J. Tang, “Nonlinear openplusclosedloop (NOPCL) control of dynamic systems,” Chaos, solitons and fractals, vol. 11, no. 7, pp. 1029–1035, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar
 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. 12, pp. 117–134, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar
Copyright
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.