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

Authors: Qingkai Han; Xueyan Zhao; Xiaoguang Yang; Bangchun Wen
Journal: Mathematical Problems in Engineering
 (2010)

---

## 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.


---

## Body

## 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 matricesA 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:(2.1)θ̈=F(θ,θ̇,t),
where θ={θ1,θ2,…,θn}T is the state variable, and n is the number of DOF of the system.The given motion goal of the system is(2.2)g={g1,g2,…,gn}T.Define a tracking error of the two-bar linkage as(2.3)e=θ-g={θ1-g1,θ2-g2,…,θn-gn}T={e1,e2,…,en}T.Equation (2.1) is linearized and expanded in the neighborhood of the goal via Taylor series and it becomes(2.4)θ̈=F(θ,θ̇,t)=F(g+e,ġ+ė,t)=F(g,ġ,t)+(∂F(g,ġ,t)∂g)e+(∂F(g,ġ,t)∂ġ)ė+o2(g,ġ)=F(g,ġ,t)+Jge+Jġė+o2(g,ġ),
where Jg and Jġ are Jacobian matrices of F(g,ġ,t) with respect to g and ġ, respectively.An OPCL controller for the system is designed as [8](2.5)U=g̈-F(g,ġ,t)-Jge-Jġė+Aė+Be,
where the term of g̈-F(g,ġ,t) is the open-loop part, and the term of -Jge-Jġė+Aė+Be 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:(2.6)θ̈=F(θ,θ̇,t)+U=F(g,ġ,t)+Jge+Jġė+o2(g,ġ)+g̈-F(g,ġ,t)-Jge-Jġė+Aė+Be=g̈+Aė+Be+o2(g,ġ).
The tracking error equation defined as (2.3) can be deduced as follows with omitting higher-order term:(2.7)ë=θ̈-g̈=g̈+Aė+Be+o2(g,ġ)-g̈=Aė+Be.It is noticed that, if the controlled system of (2.6) is asymptotically stable when coefficient matrices Aand B are constant and with negative real parts of eigenvalues, that is, their elements aii and bii should be negative, (2.7) can be rewritten in an expanded form as(2.8){ë1ë2⋮ën}=diag(a11,a22,…,ann){ė1ė2⋮ėn}+diag(b11,b22,…,bnn){e1e2⋮en}={a11ė1+b11e1a22ė2+b22e2⋮annėn+bnnen}.Thus, ifëi=aiiėi+biiei (i=1,2,…,n) is asymptotically stable, the system is also stable. Therefore, the aforementioned assumption is true.A Lyapunov functionVis defined as(2.9)V(ėi,ei)=bii2aiiei2-12aiiėi2>0.
Then(2.10)V̇(ėi,ei)=biieiėi-ėiëiaii=biieiėi-ėi(aiiėi+biiei)aii=-ė2i<0.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 Figure1. Its two joints can be driven, respectively, at o1 and o2, so that the upper link (Link 1) and the lower link (Link 2) can rotate, respectively, around their own joints of o1 and o2 within the range of [-π,π]. The upper joint o1 is fixed on the ground.Figure 1
The schematic diagram of a two-bar linkage.In coordinate systemO1xy, as shown in Figure 1, θ1 is the rotating angle of Link 1 with respect to y-axis, and θ2 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:(2.11)[M11M12M21M22]{θ̈1θ̈2}+{C1(θ1,θ2,θ̇1,θ̇2)C2(θ1,θ2,θ̇1,θ̇2)}+{K1(θ1,θ2)K2(θ1,θ2)}={τ1τ2},
where(2.12)M11=m1d12+m2(l12+d22+2l1d2cosθ2)+I1+I2,M21=m2(d22+l1d2cosθ2)+I2,M12=m2(d22+l1d2cosθ2)+I2,M22=m2d22+I2,C1(θ1,θ2,θ̇1,θ̇2)=-m2l1d2θ̇22sinθ2-2m2l1d2θ̇1θ̇2sinθ2,C2(θ1,θ2,θ̇1,θ̇2)=m2l1d2θ̇12sinθ2,K1(θ1,θ2)=(m1d1+m2l1)gsinθ1+m2d2gsin(θ1+θ2),K2(θ1,θ2)=m2gd2sin(θ1+θ2).In (2.11), mi is the mass of Link i, i=1,2; Ii is the moment of inertia of Link iwith respect to its mass center, i=1,2; di is the distance between Link i and jointi, i=1,2; gis the acceleration of gravity; τi is driving moment in joint i,i=1,2.Equation (2.11) can be rewritten in the following form, which is also in the form of   (2.1):(2.13){θ̈1θ̈2}=[H11H12H21H22]{R1R2},
where
(2.14)[H11H12H21H22]=[M11M12M21M22]-1,{R1R2}={τ1-C1(θ1,θ2,θ̇1,θ̇2)-K1(θ1,θ2)τ2-C2(θ1,θ2,θ̇1,θ̇2)-K2(θ1,θ2)}.


## 2.1. A Dynamical System with an OPCL Controller

Let a typical dynamical system be as follows:(2.1)θ̈=F(θ,θ̇,t),
where θ={θ1,θ2,…,θn}T is the state variable, and n is the number of DOF of the system.The given motion goal of the system is(2.2)g={g1,g2,…,gn}T.Define a tracking error of the two-bar linkage as(2.3)e=θ-g={θ1-g1,θ2-g2,…,θn-gn}T={e1,e2,…,en}T.Equation (2.1) is linearized and expanded in the neighborhood of the goal via Taylor series and it becomes(2.4)θ̈=F(θ,θ̇,t)=F(g+e,ġ+ė,t)=F(g,ġ,t)+(∂F(g,ġ,t)∂g)e+(∂F(g,ġ,t)∂ġ)ė+o2(g,ġ)=F(g,ġ,t)+Jge+Jġė+o2(g,ġ),
where Jg and Jġ are Jacobian matrices of F(g,ġ,t) with respect to g and ġ, respectively.An OPCL controller for the system is designed as [8](2.5)U=g̈-F(g,ġ,t)-Jge-Jġė+Aė+Be,
where the term of g̈-F(g,ġ,t) is the open-loop part, and the term of -Jge-Jġė+Aė+Be 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:(2.6)θ̈=F(θ,θ̇,t)+U=F(g,ġ,t)+Jge+Jġė+o2(g,ġ)+g̈-F(g,ġ,t)-Jge-Jġė+Aė+Be=g̈+Aė+Be+o2(g,ġ).
The tracking error equation defined as (2.3) can be deduced as follows with omitting higher-order term:(2.7)ë=θ̈-g̈=g̈+Aė+Be+o2(g,ġ)-g̈=Aė+Be.It is noticed that, if the controlled system of (2.6) is asymptotically stable when coefficient matrices Aand B are constant and with negative real parts of eigenvalues, that is, their elements aii and bii should be negative, (2.7) can be rewritten in an expanded form as(2.8){ë1ë2⋮ën}=diag(a11,a22,…,ann){ė1ė2⋮ėn}+diag(b11,b22,…,bnn){e1e2⋮en}={a11ė1+b11e1a22ė2+b22e2⋮annėn+bnnen}.Thus, ifëi=aiiėi+biiei (i=1,2,…,n) is asymptotically stable, the system is also stable. Therefore, the aforementioned assumption is true.A Lyapunov functionVis defined as(2.9)V(ėi,ei)=bii2aiiei2-12aiiėi2>0.
Then(2.10)V̇(ėi,ei)=biieiėi-ėiëiaii=biieiėi-ėi(aiiėi+biiei)aii=-ė2i<0.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 Figure1. Its two joints can be driven, respectively, at o1 and o2, so that the upper link (Link 1) and the lower link (Link 2) can rotate, respectively, around their own joints of o1 and o2 within the range of [-π,π]. The upper joint o1 is fixed on the ground.Figure 1
The schematic diagram of a two-bar linkage.In coordinate systemO1xy, as shown in Figure 1, θ1 is the rotating angle of Link 1 with respect to y-axis, and θ2 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:(2.11)[M11M12M21M22]{θ̈1θ̈2}+{C1(θ1,θ2,θ̇1,θ̇2)C2(θ1,θ2,θ̇1,θ̇2)}+{K1(θ1,θ2)K2(θ1,θ2)}={τ1τ2},
where(2.12)M11=m1d12+m2(l12+d22+2l1d2cosθ2)+I1+I2,M21=m2(d22+l1d2cosθ2)+I2,M12=m2(d22+l1d2cosθ2)+I2,M22=m2d22+I2,C1(θ1,θ2,θ̇1,θ̇2)=-m2l1d2θ̇22sinθ2-2m2l1d2θ̇1θ̇2sinθ2,C2(θ1,θ2,θ̇1,θ̇2)=m2l1d2θ̇12sinθ2,K1(θ1,θ2)=(m1d1+m2l1)gsinθ1+m2d2gsin(θ1+θ2),K2(θ1,θ2)=m2gd2sin(θ1+θ2).In (2.11), mi is the mass of Link i, i=1,2; Ii is the moment of inertia of Link iwith respect to its mass center, i=1,2; di is the distance between Link i and jointi, i=1,2; gis the acceleration of gravity; τi is driving moment in joint i,i=1,2.Equation (2.11) can be rewritten in the following form, which is also in the form of   (2.1):(2.13){θ̈1θ̈2}=[H11H12H21H22]{R1R2},
where
(2.14)[H11H12H21H22]=[M11M12M21M22]-1,{R1R2}={τ1-C1(θ1,θ2,θ̇1,θ̇2)-K1(θ1,θ2)τ2-C2(θ1,θ2,θ̇1,θ̇2)-K2(θ1,θ2)}.


## 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 Table1.Table 1
The values of structural parameters of the two-bar linkage.
Parameterm1m2I1I2lc1lc2l1l2gValue110.0830.330.51129.8The goal trajectories of the two rotating angles of the two links,θ1 and θ2, are designed as    (3.1)θ1(t)=sin(t),θ2(t)=sin(t).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 simulationn.(2)
Set the control parameter matricesA and B 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 ofAandB 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.Table 2
The values of control parameters for different motions.
Control parameterAControl parameterBMotion patternsSimulation resultsA=diag(-10,-10)B=diag(-20,-20)Single-periodic motionsSection3.2A=diag(-16,-16)B=diag(-3,-9)Multiple-periodic motionsSection3.3A=diag(-8,-8)B=diag(-4,-8)Quasiperiodic motionsSection3.4A=diag(-2.5,-2.5)B=diag(-7,-7)Chaotic motionsSection3.5The 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

GivenA=diag(-10,-10) and B=diag(-20,-20), the two-bar linkage can achieve single-periodic motions. The simulated motions with initial conditions of θ1=1, θ̇1=0, θ2=1, and θ̇2=0  are shown in Figure 2.Single-periodic motions.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2In Figures2(a) and 2(b), the motions of θ1 and θ2 are obviously periodic, that is, harmonic. The two phase plane portraits of θ1 and θ2, 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

GivenA=diag(-16,-16) and B=diag(-3,-9), 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 θ1=1, θ̇1=0, θ2=1, and θ̇2=0.The simulated multiple-periodic motions.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2In Figures3(a) and 3(b), the motions of θ1 and θ2 are also periodical but with other harmonic components. The two phase plane portraits of θ1 and θ2, 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

GivenA=diag(-16,-16) and B=diag(-3,-9), 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 θ1=1, θ̇1=0, θ2=1, and θ̇2=0.The simulated quasiperiodic motion.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2It can be seen from Figure4 that the motions of θ1 and θ2 are complex. Both the time histories and the phase plane portraits of θ1 and θ2 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 θ1 and θ2 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

GivenA=diag(-2.5,-2.5) and B=diag(-7,-7), 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 θ1=1, θ̇1=0, θ2=1, and θ̇2=0.The simulated chaotic motions.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2From Figures5(a) and 5(b), the simulated responses of θ1 and θ2 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 θ1 and θ2 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.


## 3.1. Conditions and Simulation Steps

In numerical simulations, the dimensionless structural parameters of the two-bar linkage are used and listed in Table1.Table 1
The values of structural parameters of the two-bar linkage.
Parameterm1m2I1I2lc1lc2l1l2gValue110.0830.330.51129.8The goal trajectories of the two rotating angles of the two links,θ1 and θ2, are designed as    (3.1)θ1(t)=sin(t),θ2(t)=sin(t).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 simulationn.(2)
Set the control parameter matricesA and B 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 ofAandB 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.Table 2
The values of control parameters for different motions.
Control parameterAControl parameterBMotion patternsSimulation resultsA=diag(-10,-10)B=diag(-20,-20)Single-periodic motionsSection3.2A=diag(-16,-16)B=diag(-3,-9)Multiple-periodic motionsSection3.3A=diag(-8,-8)B=diag(-4,-8)Quasiperiodic motionsSection3.4A=diag(-2.5,-2.5)B=diag(-7,-7)Chaotic motionsSection3.5The 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

GivenA=diag(-10,-10) and B=diag(-20,-20), the two-bar linkage can achieve single-periodic motions. The simulated motions with initial conditions of θ1=1, θ̇1=0, θ2=1, and θ̇2=0  are shown in Figure 2.Single-periodic motions.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2In Figures2(a) and 2(b), the motions of θ1 and θ2 are obviously periodic, that is, harmonic. The two phase plane portraits of θ1 and θ2, 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

GivenA=diag(-16,-16) and B=diag(-3,-9), 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 θ1=1, θ̇1=0, θ2=1, and θ̇2=0.The simulated multiple-periodic motions.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2In Figures3(a) and 3(b), the motions of θ1 and θ2 are also periodical but with other harmonic components. The two phase plane portraits of θ1 and θ2, 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

GivenA=diag(-16,-16) and B=diag(-3,-9), 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 θ1=1, θ̇1=0, θ2=1, and θ̇2=0.The simulated quasiperiodic motion.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2It can be seen from Figure4 that the motions of θ1 and θ2 are complex. Both the time histories and the phase plane portraits of θ1 and θ2 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 θ1 and θ2 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

GivenA=diag(-2.5,-2.5) and B=diag(-7,-7), 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 θ1=1, θ̇1=0, θ2=1, and θ̇2=0.The simulated chaotic motions.
(a)
Time history ofθ1(b)
Time history ofθ2(c)
Phase-plane portrait ofθ1(d)
Phase-plane portrait ofθ2(e)
Poincare mapping ofθ1(f)
Poincare mapping ofθ2(g)
Amplitude spectrum ofθ1(h)
Amplitude spectrum ofθ2From Figures5(a) and 5(b), the simulated responses of θ1 and θ2 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 θ1 and θ2 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 θ1 and θ2 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 θ1 and θ2 are shown in Table 4.Table 3
The Lyapunov exponents of the simulated chaotic motions of Figure5.
FirstSecondThirdFourthFifthθ10.54960.0960-0.0027-0.1008-2.0600θ20.1431-0.0371-0.0787-0.1567-2.4240Table 4
The calculated values with surrogate data method for chaotic motions.
QDμsχPθ15.52411.58514.51526.3252×10-6θ28.97151.55618.94083.8644×10-19As shown in Table3, the maximum Lyapunov exponents of the motions of angle θ1 and θ2 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θ1 and θ2, that is, 6.3252×10-6 and 3.8644×10-19, 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 ofA and Bare diagonal and with negative real parts of eigenvalues. It is reliable to force a two-bar linkage to achieve different stable motions of θ1 and θ2.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 ofA and B of the OPCL controller for θ1 and θ2 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 are6.3252×10-6 and 3.8644×10-19, smaller than  .05.



---

*Source: 986319-2010-07-27.xml*