Research Article  Open Access
Chaotic Vibration Prediction of a FreeFloating Flexible Redundant Space Manipulator
Abstract
The dynamic model of a planar freefloating flexible redundant space manipulator with three joints is derived by the assumed modes method, Lagrange principle, and momentum conservation. According to minimal joint torque’s optimization (MJTO), the state equations of the dynamic model for the freefloating redundant space manipulator are described. The PD control using the tracking position error and velocity error in the manipulator is introduced. Then, the chaotic dynamic behavior of the manipulator is analyzed by chaotic numerical methods, in which time series, phase plane portrait, Poincaré map, and Lyapunov exponents are used to analyze the chaotic behavior of the manipulator. Under certain conditions for the joint torque optimization and initial values, chaotic vibration motion of the space manipulator can be observed. The chaotic time series prediction scheme for the space manipulator is presented based on the theory of phase space reconstruction under Takens’ embedding theorem. The trajectories of phase space can be reconstructed in embedding space, which are equivalent to the original space manipulator in dynamics. The onestep prediction model for the chaotic time series and the chaotic vibration was established by using support vector regression (SVR) prediction model with RBF kernel function. It has been proved that the SVR prediction model has a good performance of prediction. The experimental results show the effectiveness of the presented method.
1. Introduction
Redundant space robots have been receiving an increasing attention in recent years because of their ability to perform high complexity tasks. Redundant robot’s joint space dimension is greater than the task space dimension, which can overcome the shortcomings of low obstacle avoidance capacity and poor flexibility and dynamic performance. So it is of particular importance in the aerospace applications. However, the dynamic behavior of space robot systems is highly nonlinear, which in fact makes the design of controllers for such systems a challenge topic. A robust control strategy for freefloating space robots was presented in [1]. The strategy consists of a set of linear time invariant controllers coupled to the robotic manipulator through a set of high feedforward gains. The threelink space manipulator is modeled and simulated and the robustness of the controller has been demonstrated. An active damping control strategy could be used to reduce the structural vibrations of a space manipulator with flexible links during its onorbit operations. An optimized adaptive vibration control via piezoelectric devices was proposed in [2]. To implement assembling, maintenance, repair, and contingency operations in space, the notion of freeflying space manipulator systems, in which robotic manipulators are mounted on a freeflying spacecraft, was discussed in [3]. The motion control of a freefloating space manipulator with parametric uncertainties and external disturbances was investigated and performed by adaptive nonlinear H1 controllers in [4]. In [5], the theoretical formulation and the experimental validation of two novel inverse kinematics solutions for redundant space manipulators aimed at locally minimizing the torque transferred to the spacecraft were presented. The first solution is based on a weighted Jacobian pseudoinverse and is derived by using Lagrange multipliers, whereas the second one is derived according to a constrained least squares approach by exploiting the linearity of the forward kinematics and the reaction torque dynamics equations. In [6], a novel composite sliding mode control was proposed to deal with the vibration suppression and trajectory tracking of a freefloating space rigidflexible coupling manipulator with a rigid payload.
The flexible redundant space robot contains complex dynamic characteristics. The research has found that, based on the pseudoinverse control algorithm for the redundant robot, there are nonconservation and nonrepeatability in link movement. The structural flexibility inevitably causes the elastic deformation and vibration during the space manipulator’s motion. However, the selfmotion planning with redundancy to suppress vibration may result in chaotic vibration and joint angle drift phenomenon [7, 8].
Chaotic phenomena and chaotic behavior have been observed in numerous natural and model systems in physics. Chaotic dynamics appear frequently in nature and in nonlinear dynamical systems [9, 10]. In [11], the nonlinear dynamic responses of the rotating continuous flexible shaftrigid disk system with rubimpact between the disk and the stator were investigated. Time series, phase plane portrait, Poincaré map, bifurcation diagrams, and Lyapunov exponents were used to analyze the dynamic behavior of the system. In [12], an intelligent fuzzy chaotic structure controls the twolink robot arm by applying proper amplitude. This work presents the control and synchronization of chaos by designing the adaptive sliding mode controller [13]. In [14], based on MJTO, the state equations of chaos motion for the freefloating redundant space manipulator are established. The robust PD (proportional and derivative) compensatory control, the delayed feedback control, and the optimal control are designed in work space, joint space, and mode space respectively. The composite control can assure the trajectory tracking, chaotic suppression, and vibration control.
Although there exist some literatures on chaotic motions in Duffing’s oscillator, inverted pendulum, and several other systems [15–17], there are very few works on the analysis of the nonlinear chaotic dynamic behavior of space robots. If chaos cannot be predicted and controlled, it might result in disaster and collapse a space robot system on many occasions. There is important significance to predict chaos motion for space robot security. Clearly, the ability to control chaos, that is, to convert the chaotic oscillation into the desired regular motion, would be beneficial.
Nonlinear time series prediction is increasingly used in many fields of engineering [18–20]. Recently, support vector machines (SVM) and support vector regression (SVR) based on statistical theory [21, 22] have become a key machine learning method and have been applied to deal with nonlinear time series prediction [23]. They have been proved to enhance the capability of prediction. For chaotic time series prediction, a novel direct prediction method based on echostate mechanisms and SVR was proposed in [24]. The method has been tested on both the benchmark problem (MackeyGlass time series) and real life time series (Wolf monthly sunspots number) and yearly runoff time series of the Yellow River.
In this paper, we propose an analytical approach to study the chaotic vibration in space manipulator. The prediction of chaotic motion is presented to improve the space robot performance and optimize the structure parameters onorbit operations. The paper is organized as follows. In Section 2, the dynamic model for a planar redundant space manipulator with three flexible links is described. In Section 3, the phase space path graph changing over time and the time history of state variables are given in phase plane, and chaos phenomena can be determined by contrast, analysis, and synthesis. In Section 4, the prediction of chaotic motion is performed by support vector regression based on the phase space reconstruction. In Section 5, the prediction results for the chaotic vibration motion of the system are shown and discussed. Finally, conclusions are given in Section 6.
2. Dynamics Description of the FreeFloating Flexible Space Manipulator
A planar structure diagram of the freefloating flexible space manipulator with the three flexible links is shown in Figure 1. The coordinate frames have been established for system modeling based on the DenavitHartenberg’s tradition. Some parameters used in Figure 1 and in the rest of the paper are defined in Abbreviations.
The axial elongation of the links is ignored. Therefore, the links deflect only in a direction transverse to their length.
It is assumed that the longitudinal deformations of the links are negligibly small and EulerBernoulli beam theory is adequate to describe their flexible motions.
The following assumptions are made to derive a dynamical model of the space manipulator.
The system consists of a rigid body satellite and manipulator composed of three flexible links. External forces and torques are not applied to the system, and thus momentum conservation and equilibrium of forces strictly hold. The flexible motion of a link is described by a finite number of vibration modes which are constructed through component mode synthesis under the assumption of small defection. By the assumed mode technique and in the case of a pinnedpinned beam, the flexible link deflection can be expressed as where is the mode’s order () and and are the mode shape functions and generalized modalcoordinates, respectively. is the coordinate along the longitudinal axis of the link .
The position vectors of the arbitrary points of the links are, respectively, given by
The position vector of the endeffector is given bywhere , , , , , , , and .
The dynamic equations of the system are derived from the total kinetic energy and total potential energy of the system. The kinetic energy of the whole manipulator is related to the base and flexible links.
The kinetic energy of the entire system is given bywhere is the kinetic energy of the base, is the kinetic energy of the th link, is a line density of the th link, are linear velocities of the base and arbitrary point of the th link, respectively. , are the inertia matrices of the base and the th link, respectively.
Potential energy of the elastic deformations of the links is expressed by where is stiffness of flexible link.
According to D’AlembertLagrange, the dynamic equations of the manipulator are given bywhere is a Lagrange function, is the joint angle, is the joint torque, and is the generalized coordinate of the flexible links.
The dynamic model is obtained by substituting the energy expressions (6) in Lagrange’s equation. As a result, the dynamical equations can be written as where is the inertia matrix, is the damping matrix, is a stiffness matrix, , , and is the joint torque vector.
Since we assume the system is operated in freefloating mode without any external force and moment acting on the entire system, the linear momentum and angular momentum of the whole system are given bywhere are masses and inertia matrices of the base and the link, respectively, is a line density of the link , is the linear velocity of the arbitrary point of the th link, , are the inertia matrices of the base and the th link, respectively, and are linear velocities of the base and the link, respectively, and , , and are position vectors of the base, the endeffector, and the link, respectively.
At the beginning of this section we assumed that the space manipulator system is operated in freefloating mode so that total linear and angular momentum are conserved; we further assume that the total linear and angular momentum are zero. On the assumption that the total momentum is zero, the total momentum of the system can be written aswhere , , , and are the coefficient matrices.
We assumed that the planar space manipulator with an endeffector is mounted on a freefloating spacecraft. The endeffector velocity of the manipulator is expressed aswhere is the velocity of the endeffector, is the Jacobian matrix for the endeffector relative to the base, and is the Jacobian matrix for the endeffector relative to the joint.
The dynamics and vibration equations of the space manipulator can be written as According to (9)–(11), the dynamic equation of the above space manipulator system is described aswhere is inertia matrix, , and are coefficient matrices.
Utilizing (9) and (10), the relationship between joint velocities and endeffector velocities iswhere are coefficient matrices.
3. Dynamic Equations of the Space Manipulator Based on MJTO
Joint torque optimization is essential to space manipulators, in which the electric power supply is limited. According to (13), the constrained equation for the manipulator can be expressed as
Under certain conditions for the joint torque optimization and initial values, the above described system exhibits chaotic behavior. This section discusses the local optimization of the joint torque using classical Lagrange optimization technique.
The following local constrained optimization problem can be defined asUsing the Lagrange multiplier , the augmented cost function can be expressed as
In order to determine the minimum solution to the unconstrained objective function, the derivatives of with respect to and must both vanish, so that we haveThen, we getwhere , , and .
Substituting (19) into (18), we obtainwhere is a weighted MoorePenrose pseudoinverse and is an unit matrix. The second tern is possible to generate internal motions that reconfigure the manipulator structure without changing the endeffector position and orientation.
Equation (20) can be written as following the state equations: where , , , , , , and
The aim of the endeffector control law is to track the desired trajectory for the endeffector. is the command acceleration vector. The proportional and derivative (PD) control can be introduced based on the state errors in the above system. The PD control law is chosen as follows:where are position and velocity gain matrices, respectively. is the desired endeffector’s trajectory. is the actual endeffector’s trajectory. and are defined as the tracking position error and velocity error respectively.
In this section, we investigate dynamical performance of the above system. The simulation results are obtained by using MATLAB. The step size is set to be 0.001 seconds. And simulation time is 50 seconds. The parameters of the planar freefloating space redundant manipulator with 3 flexible links are given in Table 1.

In simulation, we can see that the freefloating redundant space manipulator in different initial configurations and , which can guarantee the closedloop stability may lead to chaos.
The original pose of the base is set to .
The desired trajectory for the endeffector is given by According to (22), let .
In this paper we only give a set of simulation results. We consider , as the initial conditions for the system.
According to the simulation results of nonlinear differential equations of the system, as shown in Figures 2, 3, 4, and 5, the phase space path graph changing over time and the time history of state variables are given in phase plane. Chaos phenomena can be determined by contrast, analysis, and synthesis. In phase space, the periodic motion corresponds to the closed curves, and the chaotic motion corresponds to random distribution of a certain region of the path.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
It can be seen from the phase diagrams that the trajectories in a certain random distribution neither overlap nor intersect but repeatedly fold. It can be seen that the joint’s motion performs the chaotic attractor. It can also be seen that there is a certain coverage area (domain of attraction) in which the shortterm predictable and longterm unpredictable chaotic motion exist. Preliminarily judging, there is chaotic motion in dynamics of the space robot system.
The above Poincaré sections have the infinite number of discrete points. By the nature of the Poincaré section, we can see that the system has the chaotic characteristics. By means of small data sets arithmetic, the system’s Lyapunov exponents for the time series are calculated as follows: , , and . The maximum Lyapunov exponent is positive, which demonstrates the chaotic motion in the space manipulator. Figure 5 shows the chaotic vibration curves of the generalized coordinate of the flexible links.
4. Chaotic Prediction Based on SVR
4.1. Phase Space Reconstruction
Chaotic time series prediction is based on the theory of phase space reconstruction under Takens’ embedding theorem. A reconstructed phase space is dimensional metric space in which a time series is embedded. We can reconstruct a phase space trajectory from a single observation by a time delay. Given an observation chaotic time series , selecting the embedding dimension and the delay time , the phase space can be expressed as the following form:where is a vector in the reconstruction phase space and is the dimension of the vector in the reconstructed phase space. According to Takens’ theorem, for proper embedding dimension and time delay , the trajectory of phase space can be reconstructed in embedding space, which is equivalent to original system in dynamics. Therefore, (), where is the forward prediction step and is the reconstructed prediction model.
4.2. SVR Predicting Model for Chaotic Time Series
In this section, we briefly review the support vector regression (SVR) theory [21, 22]. Given a training sample set , where is the input vector and is the desired output, the above training set can be mapped into high dimension featured space, in which optimized linear regression function can be described as follows: where is a highdimensional weight vector and the bias is a scalar and is a nonlinear mapping from the input space to a feature space. For a nonlinear regression problem, SVR model is based on Vapnik’s insensitive loss function for function approximation [21]. Given the definition of Vapnik’s insensitive loss function, we havewhere is the maximum value of tolerable error.
Therefore, by means of introduction of relaxation factors , the optimization problem can be formulated in the primal space as the following:where is a punishment parameter.
The above optimization problem can be transformed into the following Lagrange dual form:where and are Lagrange multipliers and is called the kernel function. The kernel function is defined as the following radial basis function:where is the width of the function.
After the quadratic optimization problem with inequality constraints is solved, the parameters and in (28) are obtained:Finally, the SVR function is obtained by (25) in the dual space:Setting is the support vector regression corresponding to (32) and , the onestep prediction model for chaotic time series based on SVM has the following form:where is the value of onestep forecasting; then the next point in phase space can be expressed as
5. Simulation Results and Discussions
The time series responses of the system states for the manipulator using PD control are shown in Figure 2. The 5000 data points as training and test samples from the raw data of chaotic motion in the space manipulator are collected and normalized to . In order to eliminate the influence of the initial value, the initial 2000 points were discarded, the phase space reconstruction of the chaotic time series was done by the embedding dimension , and the delay time . Then 1500 points were taken as training data; 1500 points were taken as test data. The root mean square error (RMSE) was used as integral performance index to evaluate time series prediction errors between the predicted signal and its reference signal with the SVR model obtained. The RMSE is defined aswhere is a real value, is a predicted value, and is the number of test points. The width ranges of Gaussian kernel and parameters are set to 2.8, 2.2, and 0.001, respectively. The prediction experiment was done with the SVR model obtained. The prediction and absolute error curves for the chaotic time series of the manipulator are shown in Figures 6 and 7, respectively. It can be seen from Figure 6 that the actual time series responses of the links have perfect approximation to the predicted responses. Also, it can be seen from Figure 7 that the time series responses have drastic vibration for whole motion period and prediction errors show the undesirable peaks in some sample points.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(e)
(f)
The prediction values of RMSE for the chaotic time series as shown in Tables 2 and 3 are used for evaluating the performance of chaotic vibration prediction. The smaller the value of RMSE is, the closer the predicated values are to the real values.


The results of the prediction experiment show that the proposed prediction method has good prediction accuracy and successful evaluating performance in capturing and modeling the nonlinear chaotic dynamical behavior of this space manipulator.
6. Conclusions
In this paper, chaotic behavior of a planar space flexible manipulator was investigated. Time series, phase plane portrait, Poincaré map, and Lyapunov exponents were used to analyze the chaotic dynamic behavior of the manipulator. According to the largest Lyapunov exponent, phase plane portrait, and Poincaré map figures, the simulation results reveal the existence of chaotic motion in the space manipulator with three flexible links under the condition of MJTO. Based on complex dynamics inherent characteristic of chaotic and phase space reconstruction theory for the planar freefloating flexible redundant space manipulator, the training data of the time series can be obtained. And then a model based on least squares SVR (support vector regression) is proposed, which combines features of chaotic time series’ development and change, and it obtains the training samples through phase space reconstruction of the time series by considering all the information from relative parameters’ chaotic time series. Prediction experiments show that this method has good prediction accuracy and is an effective prediction method. How to control this chaos in the space manipulator will be our future research work.
Abbreviations
:  Inertia frame 
:  The frame of the endeffector 
:  The frame of the link 
:  The pose of the base 
:  The pose of the link 
:  The vector of the joint angles 
:  The position vector of the base center with reference to the inertia frame 
:  The position vector of the th link with reference to the inertia frame 
:  The position vector of the payload center with reference to the inertia frame 
:  Joint angle of the link 
:  Attitude angle of the base 
:  Length of the flexible link 
:  Radius of the base 
:  Mass of the link 
:  Mass of the base. 
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the Foundation of Manufacturing Systems Engineering State Key Laboratory, Xi’an Jiaotong University of China (no. 201002), and the Jiangsu Province Science and Technology Support Plan (no. BE2014712).
References
 P. M. Pathak, R. P. Kumar, A. Mukherjee, and A. Dasgupta, “A scheme for robust trajectory control of space robots,” Simulation Modelling Practice and Theory, vol. 16, no. 9, pp. 1337–1349, 2008. View at: Publisher Site  Google Scholar
 M. Sabatini, P. Gasbarri, R. Monti, and G. B. Palmerini, “Vibration control of a flexible space manipulator during on orbit operations,” Acta Astronautica, vol. 73, no. 4, pp. 109–121, 2012. View at: Publisher Site  Google Scholar
 S. A. A. Moosavian and E. Papadopoulos, “Freeflying robots in space: an overview of dynamics modeling, planning and control,” Robotica, vol. 25, no. 5, pp. 537–547, 2007. View at: Publisher Site  Google Scholar
 T. F. P. A. T. Pazelli, M. H. Terra, and A. A. G. Siqueira, “Experimental investigation on adaptive robust controller designs applied to a freefloating space manipulator,” Control Engineering Practice, vol. 19, no. 4, pp. 395–408, 2011. View at: Publisher Site  Google Scholar
 S. Cocuzza, I. Pretto, and S. Debei, “Novel reaction control techniques for redundant space manipulators: theory and simulated microgravity tests,” Acta Astronautica, vol. 68, no. 1112, pp. 1712–1721, 2011. View at: Publisher Site  Google Scholar
 W. Congqing, W. Pengfei, Z. Xin, and P. Xiwu, “Composite sliding mode control for a freefloating space rigidflexible coupling manipulator system,” International Journal of Advanced Robotic Systems, vol. 10, article 124, 2013. View at: Publisher Site  Google Scholar
 S. Luo and S. Ahmad, “Predicting the drift motion for kinematically redundant robots,” IEEE Transactions on Systems, Man, and Cybernetics, vol. 22, no. 4, pp. 717–728, 1992. View at: Publisher Site  Google Scholar  MathSciNet
 F. B. M. Duarte and J. A. T. Machado, “Motion chaos in the pseudoinverse control of redundant robots,” in Proceedings of the 6th International Workshop on Advanced Motion Control, pp. 624–629, IEEE, Nagoya, Japan, 2000. View at: Publisher Site  Google Scholar
 A. S. Ravishankar and A. Ghosal, “Nonlinear dynamics and chaotic motions in feedbackcontrolled twoand threedegreeoffreedom robots,” The International Journal of Robotics Research, vol. 18, no. 1, pp. 93–108, 1999. View at: Google Scholar
 H. M. Khanlo, M. Ghayour, and S. ZiaeiRad, “Chaotic vibration analysis of rotating, flexible, continuous shaftdisk system with a rubimpact between the disk and the stator,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 1, pp. 566–582, 2011. View at: Publisher Site  Google Scholar
 R. Wang and Z. Jing, “Chaos control of chaotic pendulum system,” Chaos, Solitons and Fractals, vol. 21, no. 1, pp. 201–207, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. A. Zanjani, A. H. Jafari, and A. M. Nasrabadi, “Approach to OGY control fuzzification of 2link rigid robot arm with fixedpoints distances considering (using method similar to bellmanford algorithm),” European Journal of Scientific Research, vol. 39, no. 4, pp. 469–476, 2010. View at: Google Scholar
 S. Dadras and H. R. Momeni, “Adaptive sliding mode control of chaotic dynamical systems with application to synchronization,” Mathematics and Computers in Simulation, vol. 80, no. 12, pp. 2245–2257, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 C.Q. Wang, P.F. Wu, and X. Zhou, “Control and modeling of chaotic dynamics for a freefloating rigidflexible coupling space manipulator based on minimal joint torque's optimization,” Acta Physica Sinica, vol. 61, no. 23, Article ID 230503, 2012. View at: Publisher Site  Google Scholar
 S. M. Shahruz and L. R. Siva, “Suppression of chaos in a class of nonlinear systems by disturbance observers,” Journal of Sound and Vibration, vol. 271, no. 12, pp. 481–491, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 M. T. Arjmand, H. Sadeghian, H. Salarieh, and A. Alasty, “Chaos control in {AFM} systems using nonlinear delayed feedback via sliding mode control,” Nonlinear Analysis. Hybrid Systems. An International Multidisciplinary Journal, vol. 2, no. 3, pp. 993–1001, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Salarieh and A. Alasty, “Control of chaos in atomic force microscopes using delayed feedback based on entropy minimization,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 637–644, 2009. View at: Publisher Site  Google Scholar
 B. Samanta, “Prediction of chaotic time series using computational intelligence,” Expert Systems with Applications, vol. 38, no. 9, pp. 11406–11411, 2011. View at: Publisher Site  Google Scholar
 D. Li, M. Han, and J. Wang, “Chaotic time series prediction based on a novel robust echo state network,” IEEE Transactions on Neural Networks and Learning Systems, vol. 23, no. 5, pp. 787–797, 2012. View at: Publisher Site  Google Scholar
 M. Behmanesh, M. Mohammadi, and V. S. Naeini, “Chaotic time series prediction using improved ANFIS with imperialist competitive learning algorithm,” International Journal of Soft Computing and Engineering, vol. 4, no. 4, pp. 25–33, 2014. View at: Google Scholar
 V. N. Vapnik, The Nature of Statistical Learning Theory, Springer, New York, NY, USA, 1995. View at: Publisher Site  MathSciNet
 S. Mukherjee, E. Osuna, and F. Girosi, “Nonlinear prediction of chaotic time series using support vector machines,” in Proceedings of the 7th IEEE Workshop on Neural Networks for Signal Processing (NNSP '97), pp. 511–520, Amelia Island, Fla, USA, September 1997. View at: Publisher Site  Google Scholar
 W. Wang, C. Men, and W. Lu, “Online prediction model based on support vector machine,” Neurocomputing, vol. 71, no. 4–6, pp. 550–558, 2008. View at: Publisher Site  Google Scholar
 Z. Shi and M. Han, “Support vector echostate machine for chaotic timeseries prediction,” IEEE Transactions on Neural Networks, vol. 18, no. 2, pp. 359–372, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Congqing Wang and Linfeng Wu. 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.