Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2017 / Article
!A Corrigendum for this article has been published. To view the article details, please click the ‘Corrigendum’ tab above.

Research Article | Open Access

Volume 2017 |Article ID 8975743 |

Chenming Li, Han Zhao, Shengchao Zhen, Kang Huang, Hao Sun, Ke Shao, Bin Deng, "Trajectory Tracking Control of Parallel Manipulator Based on Udwadia-Kalaba Approach", Mathematical Problems in Engineering, vol. 2017, Article ID 8975743, 12 pages, 2017.

Trajectory Tracking Control of Parallel Manipulator Based on Udwadia-Kalaba Approach

Academic Editor: Emiliano Mucchi
Received01 Aug 2017
Accepted27 Nov 2017
Published14 Dec 2017


There have been many approaches for achieving the trajectory tracking control of parallel manipulator. However, these approaches are complex for calculating Lagrangian multipliers. In this paper, unlike the former approaches, a new approach of trajectory tracking control which is based on Udwadia-Kalaba approach is presented. Using this methodology, we can obtain a concise and explicit equation of motion and consider holonomic and nonholonomic constraint whether it is ideal or nonideal simultaneously. The most important difference is that we divide constraints into structural constraints and performance constraints in this paper. Structural constraints are used to establish dynamic model without regard for trajectory control. And performance constraints are used to represent the desired trajectory. For the parallel manipulator, a nonlinear dynamics system, its constraint forces are obtained by second-order constraints. And the numerical simulation in MATLAB shows the parallel manipulator’s movement meets the requirement; tracking trajectory is exact and perfect. Through this paper, we can see that the method can simplify calculation availably.

1. Introduction

Parallel manipulator is a kind of closed-loop kinematic chain mechanism in which the end-effector is linked to the base by several independent kinematic chains [1]. In comparison with the serial manipulator which has open-loop structure, the parallel manipulator has many advantages, such as higher stiffness, higher accuracy, higher payload capacity, lower inertia, and smaller workspace. Therefore, it has attracted a great amount of attention and interest among the researchers and engineers and has been applied to many industrial applications in the past decade, like milling [2], rapid machining [3], medical resuscitation [4], and so on.

In order to achieve the great potential performances existing in the parallel manipulator system, the controller design is essential and indispensable. To deal with the trajectory tracking control problem of the parallel manipulator, it is easy to find that many advanced controllers have been proposed and successfully implemented, like proportional-integral-derivative control [5], sliding mode control [6], adaptive control [7, 8], robust control [9, 10], robust PID control [11], adaptive sliding mode control [12], adaptive robust control [13], fuzzy control [14] neural network control [15], fuzzy neural network control [16], and so on.

Parallel manipulators can be designed with various mechanical structures because the designers want them to realize different control requirements according to the actual demands. In this paper, a 3-degree-of-freedom (3-DOF) planar parallel manipulator (3-HSS, H-helix pair, S-spherical pair) has been designed and constructed. To deal with the trajectory tracking control problem of this manipulator, our methodology is based on a novel, concise approach to explicitly formulate the equations of motion for constrained systems proposed by Udwadia and Kalaba [1720].

In this paper, we will use the basic kinematics equations proposed by Udwadia and Kalaba based on Gauss’s theory for modeling and trajectory tracking of the parallel manipulators. There are several key points by using the method in this paper compared with the traditional control method. First, there is no need to obtain the complex Lagrangian multipliers during the process of dynamic modeling of the parallel manipulators. Second, we define the constraint equations by a new perspective which means the constraints are divided into structural constraints and performance constraints. Structural constraints are used to set up the dynamic model in the absence of trajectory control. And performance constraints are used to reflect the desired trajectory of the parallel manipulator. Third, we can obtain the needed motor torque for realizing the trajectory of manipulator directly by calculating the Udwadia and Kalaba equation. At last, the simulation results in this paper show that this method can control the parallel manipulator accurately.

2. Udwadia and Kalaba Equation for Constrained Mechanical Systems

In this section, we talk about the fundamental equation, called Udwadia-Kalaba equation, for constrained mechanical system. We can obtain the explicit closed-form equation of motion of the constrained mechanical system in three steps [18].

First, we consider the unconstrained mechanical system in which the coordinates are all assumed to be independent of each other. The equation of motion of this system can be expressed, by using Newtons or Lagranges equation, aswith the initial conditionswhere is the generalized coordinate -vector and is the time; is a positive definite inertia matrix; is an -vector of velocity; is an -vector of acceleration; is an -vector called the given force which could include centrifugal force, gravitational force, and control input. From (2), the generalized acceleration of the unconstrained system is given by

Second, a set of control requirement or trajectory requirements is imposed as constraints on this uncontrolled system. We shall assume that the system is subjected to holonomic constraints of the form meanwhile there are nonholonomic constraints in the form of

Then, we can differentiate these constraint equations in Lagrangian mechanics which are usually in Pfaffian form. Under the assumption of sufficient smoothness, 2-order constraint equations can be acquired by taking the time derivative once on nonholonomic constraints and twice on holonomic constraints, which can be written in the form of matrix aswhere is an matrix denoting the constraint matrix and is an vector.

Remark 1. The constraint is in either the zeroth-order or first-order form which is expressed in the form of Lagrangian equation, Maggi equation, Boltzmann and Hamel equation, and Gibbs and Appell equation. However, for further dynamic analysis and control design, second-order constraints are considered to be the most suitable form. The reason is that the acceleration will be linear by using the second-order constraint. And there is no need to worry about the loss of information during the differentiation (e.g., a constant). In fact, making the zeroth-order or holonomic constraint meet the initial condition means the initial condition retained the missing information. Thus, when we use Udwadia-Kalaba equation, the first step is to find all constraints including holonomic constraints and nonholonomic constraints in the unconstrained system. Then, we could obtain the second-order form of constraints by differentiating the constraints with respect to time once or twice. Finally, we could transform these constraints into second-order matrix equations like (6).
Eventually, the explicit equations of motion with constraints will be formed. Notice that the additional generalized forces of constraints need to be imposed on the system as a result of the presence of constraints. We therefore assume the actual explicit equation of motion of the constrained system in the form ofwhere is the additional generalized force which is imposed on the unconstrained system, arising by the holonomic and nonholonomic constraints denoting a by 1 constraint matrix. Determining the general explicit form of is our goal with the known at any time . In Lagrangian mechanics, is considered to be ideal constraints because this kind of forces of constraint does zero work under virtual displacement according to the D’ Alembert’s Principle. However nonideal constraints also exist in an actual mechanical system, and the nonideal constraint forces including friction force, electromagnetic force, will be generated. Udwadia expressed in the form ofwhere represents the ideal constraint force and represents the nonideal constraint.
Furthermore, Udwadia and Kalaba have proved the general “explicit” equations of ideal and nonideal constraint forces at any time can be expressed in the form ofwhere and the superscript “+” denotes the Moore-Penrose generalized inverse.
From (7), (8), and (9), the general “explicit” equation of motion can be expressed as

Remark 2. In particular, is a suitable -vector in (10) which is determined by the mechanical system. Now we define a given nonzero -vector virtual displacement in a constrained system at the instant time and the work done generated by under virtual displacement . And the work done which in any displacement which is subject to is the same as the work done by . Hence, we put the equation in form ofwhich means we extend Lagrange’s form of D’ Alembert’s Principle. The work done by the ideal constraint force and the nonideal constraint force under virtual displacements iswhen the constraints are ideal and the total work done under virtual displacement is zero which means equals zero; according to D’ Alembert’s Principle, (10) becomes

Therefore, at any instant of time, the constrained system is subjected to an additional constraint force , given by

If the unconstrained acceleration is zero with the constant diagonal matrix (14) becomes

Remark 3. There are two main advantages of this equation which Udwadia and Kalaba have provided. First, by using this equation, there is no need to determine the Lagrangian multipliers which are always complex. On the contrary, the parameters in the equation are all simple to obtain. Second, this equation applies to all holonomic and/or nonholonomic constrained systems including ideal or nonideal constraints.

3. Trajectory Tracking Control of the Parallel Manipulator

3.1. Dynamic Model of the Parallel Manipulator

The 3-HSS parallel manipulator is composed of a moving platform and three sets of sliding saddle chain structure. Each branched chain includes two equilong and parallel bars. One end of the bar is connected with the saddle by rolling spherical joint, and the other end is connected to the moving platform. As for the sliding saddle which is driven by the servomotor-screw-nut pair, it is moving along the ball screw which is mounted on the guide block. Hence the moving platform could achieve the 3D-motion. In addition, the platform is capable of high speed machining with the high speed motorized spindle on it.

The law of movement of each bar in branched chain is the same based on the parallelogram strut structure. Therefore, in this paper, we abstract the original mechanism into single-bar equivalent model as shown in Figure 1 during the kinematics modeling process. In this equivalent mechanism, each branched chain contains only one bar. And two ends of the bar are connected with the moving platform and the saddle by Hooke joint. Beside the saddle is driven by active pair (helix pair). Based on the above structure, the equivalent mechanism has the same degree of freedom as the original mechanism.

Next, we employ the classic universal joint, prismatic joint and spherical joint to illustrate the dynamic modeling of the delta parallel structure (as show in Figure 1). Following the procedure, we segment the delta parallel structure into four subsystems: three arm subsystems (1, 2, 3) and one platform subsystem . We shall then derive the equations of motion of the “unconstrained” subsystems.

Select reference frame . And assume the mechanism fixed on the plane . Then consider the arm-subsystem in Figure 1 separately. Denote the following.

As shown in Figure 2, the end of the branched chain which is connected with the moving platform can be described with coordinates . And the moving platform can be described with coordinates .

is height of the other end of the th branched chain. is the angle between th branched chain and -axis. is the angle between the th branched chain’s projection in the plane and -axis. is the length of the th branched chain.

Thus, when we consider the arm-subsystem separately, the motor’s coordinate of location can be expressed as and the end of the branched chain’s coordinate of location can be expressed as . Then the kinetic and potential energy of the system can be written aswhere is the system’s kinetic/potential energy, is the motor’s kinetic/potential energy, is the branched chain’s kinetic/potential energy, and is the mass of motor/branched chain.

Next, adopting the standard Lagrangian approach, we have the following equation of motion for each “unconstrained” arm subsystem :where

For the “unstrained” platform, that is, the subsystem 4 simplified as mass, by choosing the system coordinates as in Figure 2, the equation of motion is

We have

So, the dynamic equations of the system is

Now we cluster all the 4 subsystems together by lettingwhere

3.2. Structural Constraint of the Parallel Manipulator

We now introduce all the kinematic constraints among the subsystems and convert them into the second-order form.

Since the spherical joint only fixes the three translational DOF’s between the connected links, we have the kinematic constraints, for as follows:

Choosing generalized coordinate system

So, we can get

Differentiating (26) with respect to twice, we can get

Let , (27) be written as matrix in the form of :

Assume matrix, :

The total parallel manipulator system can be written in a matrix form:

4. Trajectory Tracking Simulation of the Parallel Manipulator

In this section, we aim to use the Udwadia-Kalaba approach to obtain an explicit, closed-form expression for the control force required to precisely meet the trajectory requirements of the 3-HSS parallel manipulator.

4.1. Performance Constraint of the Parallel Manipulator

Case 1. Assume that the prespecified trajectory of the mobile platform is represented by the following equation:

Differentiating the constraint equation (31) with respect to twice, we can get

Hence the system constraints can be written in the form ofwhere

Equation (31) is the performance constraints in this paper. Thus, the three servomotors which are on the saddle must make the platform satisfy (31) by controlling the three groups of branches. Finally, we consider the structural and performance constraints at the same time, which means we put , and in one matrix. So, we obtain

According to Udwadia-Kalaba approach, the closed-form control force can be expressed as

Case 2. In Case 2, we let the movement of platform be more complex in order to verify the control approach. The movement in direction is uniform motion, in direction is trigonometric function, and in direction is parabola. Thus, the desired trajectory is given as follows.

Similarly, differentiating the constraint equation (37) with respect to twice, we can get

Therefore, the matrices can be obtained. The control force can be obtained as well.

4.2. Simulation Results

In this section, the explicit analytic results given in Section 4 are verified by numerical simulations. The numerical integration through-out this paper is done in the MATLAB environment using a variable time step ode 15i integrator and the time of simulation is 1000 seconds. The values of the parameters are defined in Table 1.


Mass of saddle10 kg
Mass of branched chain.2 kg
Mass of moving platform15 kg
Moment of inertia around - and -axis
Length of branched chain m
Radius of constrained motion m
Distance between two frame4 m
Angular velocity of moving platform0.02 π

We give the initial conditions of Case 1 including initial coordinate and velocity of the generalized coordinates in Table 2.


Figures 3 and 4 show the simulated trajectory of the parallel manipulator in different views. From the two figures, we can see the moving platform can move on the given trajectory by solving the Udwadia-Kalaba equation (36). In Figures 5, 6, and 7, we show the variation of the three parameters , respectively. Figure 8 proved the distance between the platform and three saddles is constant; in other words, the length of branched chains is unchangeable which means the method in this paper is right. Figures 9 and 10 show the error of the platform in and directions, respectively. The errors are really small which means the platform can move along the desired trajectory. Figure 11 shows the forces generated by motors which verified the new approach validity.

In Case 2, we give the initial conditions in Table 3.


Figures 12 and 13 show the simulated trajectory in Case 2. From the Figures 13(a), 13(b), and 13(c), we can see the trajectory in , , and directions clearly. Figures 14, 15, and 16 show the variation of the three parameters , respectively. Figures 17, 18, and 19 show the error of the platform in , , and directions, respectively. The errors are all very small which are, respectively, in the order of , , and . Figure 20 shows the forces generated by motors in Case 2.

5. Conclusion

In this paper, Udwadia-Kalaba approach is applied for trajectory tracking control of the 3-HSS parallel manipulator for the first time. First, we establish the nonlinear dynamic model by dividing the system into four “unconstrained” subsystems which simplifies the modeling process. Next, the system structure and desired trajectory are regarded as the structure constraints and performance constraints, respectively. Finally, this approach provides an explicit, closed-form analytical expression for the control force by solving Udwadia-Kalaba equation. The simulation results of two cases show that the servoconstraint forces based on Udwadia-Kalaba equation make the platform’s movement meets the designed trajectory precisely.

In summary, there are two main advantages of this approach in this paper. First, this paper provides an easy dynamic modeling process for 3-HSS parallel manipulator. The constraints in the system which can be holonomic as well as nonholonomic are divided into structure constraints and performance constraints. What we need do is to get the second-order-form constraints. It is more convenient compared with the common approach. Second, the control performance is excellent due to the exact equation of motion by this approach. However, in another approach, there are approximations of the equation of motion or modification of the controller such as Lagrange multiplier which makes the equation complicated and the tracking performance not optimistic.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The theoretical foundation of this paper is consistent with (1) National Natural Science Foundation of China (51505116); (2) the Fundamental Research Funds for the Central Universities (JZ2016HGTB0716); (3) Natural and Science Foundation of Anhui Province (1508085SME221); (4) China Postdoctoral Science Foundation (2016M590563). In addition, the project is supported by these foundations.


  1. J. P. Merlet, Parallel Robots, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2006.
  2. J. Wu, J. Wang, T. Li, and L. Wang, “Performance analysis and application of a redundantly actuated parallel manipulator for milling,” Journal of Intelligent & Robotic Systems, vol. 50, no. 2, pp. 163–180, 2007. View at: Publisher Site | Google Scholar
  3. J. Kim, F. C. Park, S. J. Ryu et al., “Design and analysis of a redundantly actuated parallel mechanism for rapid machining,” IEEE Transactions on Robotics and Automation, vol. 17, no. 4, pp. 423–434, 2001. View at: Publisher Site | Google Scholar
  4. Y. Li and Q. Xu, “Design and development of a medical parallel robot for cardiopulmonary resuscitation,” IEEE/ASME Transactions on Mechatronics, vol. 12, no. 3, pp. 265–273, 2007. View at: Publisher Site | Google Scholar
  5. Y. Su, D. Sun, L. Ren, and J. K. Mills, “Integration of saturated PI synchronous control and PD feedback for control of parallel manipulators,” IEEE Transactions on Robotics, vol. 22, no. 1, pp. 202–207, 2006. View at: Publisher Site | Google Scholar
  6. N.-I. Kim, C.-W. Lee, and P.-H. Chang, “Sliding mode control with perturbation estimation: Application to motion control of parallel manipulator,” Control Engineering Practice, vol. 6, no. 11, pp. 1321–1330, 1998. View at: Publisher Site | Google Scholar
  7. W. Shang and S. Cong, “Nonlinear adaptive task space control for a 2-DOF redundantly actuated parallel manipulator,” Nonlinear Dynamics, vol. 59, no. 1-2, pp. 61–72, 2010. View at: Publisher Site | Google Scholar
  8. V. Vinoth, Y. Singh, and M. Santhakumar, “Indirect disturbance compensation control of a planar parallel (2-PRP and 1-PPR) robotic manipulator,” Robotics and Computer-Integrated Manufacturing, vol. 30, no. 5, pp. 556–564, 2014. View at: Publisher Site | Google Scholar
  9. K. Fu and J. K. Mills, “Robust control design for a planar parallel robot,” International Journal of Robotics & Automation, vol. 22, no. 2, pp. 139–147, 2007. View at: Google Scholar
  10. W. W. Shang and S. Cong, “Robust nonlinear control of a planar 2-DOF parallel manipulator with redundant actuation,” Robotics and Computer-Integrated Manufacturing, vol. 30, no. 6, pp. 597–604, 2014. View at: Publisher Site | Google Scholar
  11. M. A. Khosravi and H. D. Taghirad, “Robust PID control of fully-constrained cable driven parallel robots,” Mechatronics, vol. 24, no. 2, pp. 87–97, 2014. View at: Publisher Site | Google Scholar
  12. M. Zeinali and L. Notash, “Adaptive sliding mode control with uncertainty estimator for robot manipulators,” Mechanism and Machine Theory, vol. 45, no. 1, pp. 80–90, 2010. View at: Publisher Site | Google Scholar
  13. X. Zhu, G. Tao, B. Yao, and J. Cao, “Adaptive robust posture control of parallel manipulator driven by pneumatic muscles with redundancy,” IEEE/ASME Transactions on Mechatronics, vol. 13, no. 4, pp. 441–450, 2008. View at: Publisher Site | Google Scholar
  14. A. Noshadi, M. Mailah, and A. Zolfagharian, “Intelligent active force control of a 3-RRR parallel manipulator incorporating fuzzy resolved acceleration control,” Applied Mathematical Modelling, vol. 36, no. 6, pp. 2370–2383, 2012. View at: Publisher Site | Google Scholar
  15. H. Sadjadian, H. D. Taghirad, and A. Fatehi, “Neural networks approaches for computing the forward kinematics of a redundant parallel manipulator,” International Journal of Computational Intelligence, vol. 2, no. 1, pp. 40–47, 2005. View at: Google Scholar
  16. D. Zhu and Y. Fang, “Adaptive control of parallel manipulators via fuzzy-neural network algorithm,” Control Theory and Technology, vol. 5, no. 3, pp. 295–300, 2007. View at: Publisher Site | Google Scholar
  17. F. E. Udwadia and R. E. Kalaba, “A new perspective on constrained motion,” Proceedings: Mathematical and Physical Sciences, vol. 439, no. 1906, pp. 407–410, 1992. View at: Google Scholar
  18. F. E. Udwadia and R. E. Kalaba, Analytical dynamics: A new approach, Cambridge University Press, Cambridge, UK, 1996.
  19. F. E. Udwadia and R. E. Kalaba, “Explicit equations of motion for mechanical systems with nonideal constraints,” Journal of Applied Mechanics, vol. 68, no. 3, pp. 462–467, 2001. View at: Publisher Site | Google Scholar
  20. F. E. Udwadia and R. E. Kalaba, “What is the general form of the explicit equations of motion for constrained mechanical systems?” Journal of Applied Mechanics, vol. 69, no. 3, pp. 335–339, 2002. View at: Publisher Site | Google Scholar

Copyright © 2017 Han Zhao 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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.