Next Article in Journal
Geometric and Colour Data Fusion for Outdoor 3D Models
Previous Article in Journal
Prediction of Gap Asymmetry in Differential Micro Accelerometers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Torque Reduction of a Three Dimensional Redundant Planar Manipulator

by
Samer Yahya
1,*,
Mahmoud Moghavvemi
1,2 and
Haider Abbas F. Almurib
3
1
Center of Research in Applied Electronics (CRAE), University of Malaya, Kuala Lumpur 50603, Malaysia
2
Faculty of Electrical and Computer Engineering, University of Tehran, P.O. Box 14399-57131, Tehran, Iran
3
Department of Electrical & Electronic Engineering, University of Nottingham Malaysia, Jalan Broga, Semenyih 43500, Malaysia
*
Author to whom correspondence should be addressed.
Sensors 2012, 12(6), 6869-6892; https://doi.org/10.3390/s120606869
Submission received: 1 April 2012 / Revised: 2 May 2012 / Accepted: 22 May 2012 / Published: 25 May 2012
(This article belongs to the Section Physical Sensors)

Abstract

: Research on joint torque reduction in robot manipulators has received considerable attention in recent years. Minimizing the computational complexity of torque optimization and the ability to calculate the magnitude of the joint torque accurately will result in a safe operation without overloading the joint actuators. This paper presents a mechanical design for a three dimensional planar redundant manipulator with the advantage of the reduction in the number of motors needed to control the joint angle, leading to a decrease in the weight of the manipulator. Many efforts have been focused on decreasing the weight of manipulators, such as using lightweight joints design or setting the actuators at the base of the manipulator and using tendons for the transmission of power to these joints. By using the design of this paper, only three motors are needed to control any n degrees of freedom in a three dimensional planar redundant manipulator instead of n motors. Therefore this design is very effective to decrease the weight of the manipulator as well as the number of motors needed to control the manipulator. In this paper, the torque of all the joints are calculated for the proposed manipulator (with three motors) and the conventional three dimensional planar manipulator (with one motor for each degree of freedom) to show the effectiveness of the proposed manipulator for decreasing the weight of the manipulator and minimizing driving joint torques.

1. Introduction

Theoretically, for a structure of the robot manipulator one actuator can be mounted on each link to drive the next link via a speed reduction unit, but actuators and speed reducers installed on the distal end become the load for actuators installed on the proximal end of a manipulator, resulting in a bulky and heavy system [1]. To reduce the weight and the inertia of a robot manipulator, many mechanisms have been proposed so far to remove the weight restriction. Some reported by [2,3] include:

  • Lightweight joint design based on a special rotary joint [46]

  • Provision of a powerful slider at the base to bear as much required driving force as possible [7]

  • The parallel mechanism is another method to reduce the mass and inertia of the manipulator [8]. A typical parallel manipulator consists of a moving platform that is connected with a fixed base by several limbs. Generally, the number of degrees of freedom of a parallel manipulator is equal to the number of its limbs. The actuators are usually mounted on or near the base, which contributes to reduce the inertia of manipulators, and

  • Concentration of the actuators at the base and transmission of the power to each joint through tendons or a special transmission mechanism [2,3,9]. This mechanism allows the actuators to be situated remotely on the manipulator base, allowing the manipulator to be made more lightweight and compact.

For a serial manipulator, direct kinematics are fairly straightforward, whereas inverse kinematics becomes very difficult. Reference [10] proposes a fused smart sensor network to estimate the forward kinematics of an industrial robot, while reference [11] measures the range data with respect to the robot base frame using the robot forward kinematics and the optical triangulation principle. The inverse kinematics problem is much more interesting and its solution is more useful, but one of the difficulties of inverse kinematics is that when a manipulator is redundant, it is anticipated that the inverse kinematics has an infinite number of solutions. This implies that, for a given location of the manipulator's end-effector, it is possible to induce a self-motion of the structure without changing the location of the end-effector. In this paper we depend on our prior works [12,13] which present a new method to solve the problem of multi-solutions of a three dimensinal planar redundant manipulator. Because this paper explains the dynamic of the manipulator and not its kinematics, the inverse kinematics methods will not be explained here. For more details about the inverse kinematics of redundant manipulators, our works [1416] can be checked.

It is mentioned earlier that the proposed manipulator could be used to reduce the weight of the manipulator which yields to a decrease in the size (power) of the motors used to control the manipulator. To show the effectiveness of the proposed manipulator in reducing the torques of its motors the inverse dynamic of the manipulator has been calculated mathematically. The inverse dynamic model provides the joint torques in terms of the joint positions, velocities and accelerations. For robot design, the inverse dynamic model is used to compute the actuator torques, which are needed to achive a desired motion [17]. Several approaches have been proposed to model the dynamics of robots. The most frequently employed in robotics are the Lagrange formulation and the Newton-Euler formulation. Because the Lagrange formulation is conceptually simple and systematic [18], it has been used in this paper. The Lagrange formulation provides a description of the relationship between the joint actuator forces and the motion of the mechanism, and fundamentally operates on the kinetic and potential energy in the system [19].

The work presented in this paper is based on our previous work [14], which presents a mechanical design for a three dimensional planar redundant manipulator, which guarantees to decrease the weight of the manipulator by decreasing the number of motors needed to control it. Because the inverse kinematics model gives an infinite number of solutions for a redundant manipulator, consequently, secondary performance criteria can be optimized [17], such as avoiding singular configurations and minimizing driving joint torques. Reference [14] studied the kinematics of the manipulator of this paper and showed in details its ability to avoid singular configurations. A comparison of the manipulability index values and the manipulability ellipsoids for the manipulator is made with the manipulability index values and the manipulability ellipsoids of the PUMA arm to show the effectiveness of using the proposed manipulator to avoid singularity. In this paper, the dynamics of this manipulator are explained in detail. The contribution of this work is to explain the ability of this manipulator for joint torque minimization. The links and motors mass distribution is studied for both the proposed (with three motors) and conventional manipulators (six motors). The driving joint torques have been studied for the proposed manipulator for each joint and the results are compared with the results of the conventional manipulators to show the effectiveness of this manipulator for minimizing driving joint torques.

2. The Mechanical Design of the Manipulator

To control the motion of the end-effector of the manipulator shown of Figure 1(a), all the motors of the manipulator should be controlled. For example, to control a five links planar redundant manipulator with the ability to rotate the entire manipulator around its vertical axis, the six motors (five motors for each joint angle and one motor to rotate the entire manipulator around its vertical axis) of the manipulator should be controlled. Using the method of our papers [12,13], the configuration of the manipulator will have three angles to be controlled instead of n angles. Figure 1(b) shows the configuration of the manipulator when there are just three angles that need to be controlled.

Because the end-effector can follow any desired path by controlling three angles (θ1, θ2 and θ3) only, therefore instead of using a motor for each joint angle, three motors can be used for controlling the manipulator. This means that for any number of degrees of freedom three dimensional planar redundant manipulators, the weight of the links will be significantly decreased using the proposed design. To make the manipulator capable of moving in a three dimensional work space, one motor will control the value of θ1—this means controlling the rotation of the entire manipulator around the vertical axis. This motor is situated in such a way as to rotate the base of the manipulator around the z-axis. The second motor controls the value of θ2, which means the rotation of the entire manipulator with its configuration. The motor is situated at the base. The third motor controls the value of θ3 and this motor is situated on the first link. This motor will rotate the second link of manipulator about the second axis, and because all the next links should rotate about their axes by the same angle θ3 therefore, there is no need to use motor for each joint angle, but the rotation of the second motor will be transferred to the next joints using gears boxes. Figure 2 shows the mechanism of the proposed manipulator.

Elaborating further, the second motor is connected to the first link using a worm gear to control the angle θ2. Figure 3 shows the position of the second motor.

The third motor is connected to the second link using a worm gear for the same reasons it was used with the first link. Controlling the third motor means controlling the angle between the first link and the second link i.e., the angle θ3. Figure 4 shows the position of the third motor.

The mechanism of the third link is shown in Figure 5. The same mechanism of the second link is used; the only one difference is that instead of using s worm as a driver and s wheel gear as a driven, two bevel gears are used. The same mechanism of the third link can be used with the next links. The last link has the mechanism shown in the Figure 6. For further details of the mechanical design of the manipulator, our reference [14] can be checked.

To ensure that all the links move at the same joint angle, the ratio between the bevel gears of each planetary gear should be equal to one. This means the bevel gears of each planetary gear should have the same diameter and number of teeth. If this arm is fixed, we get:

w 1 w 2 = N 2 N 1
where w is the angular velocity of gear and N is the number of teeth of gear. In our manipulator, it is noted that the first gear is fixed while the second gear and the arm are rotating. It is desired that both the arm and the second gear have the same angular velocity. Because the arm is not stationary, then we cannot use the previous equation. i.e., the mechanism is not an ordinary gear train but a planetary gear train. To convert this planetary gear train to an ordinary gear train, it is assumed that the arm is stationary while a first gear has an angular velocity and not fixed. This means that:
w 1 = w 1 w a
w a = w a w a = 0

And because the second gear will continue rotating with the same angular velocity, then:

w 2 = w 2

Now the Equation (1) can be rewritten as follows:

w 1 w 2 = N 2 N 1 = w 1 w a w 2

For our manipulator it is desired to move both the arm and the second gear by the same angular velocity w which means:

( 0 ) w w = N 2 N 1 N 1 = N 2

To make the manipulator to have the ability to move in a three dimensional work space, a motor is added to the base of the manipulator to make the whole manipulator capable of rotating around the z-axis. This motor controls θ1. Figure 7 shows the mechanism of the first motor.

To calculate the transformation matrix of the manipulator, the draft of the manipulator shown in Figure 8, is used. The corresponding link parameters of the manipulator are shown in Table 1. Where l1, l2, …, l5 are the length of the links, while d1 is the offset between the origin and the end-effector.

From the links parameters shown in Table 1 and using Equation (7) which defines the transformation matrix T for the links [1], we compute the individual transformations for each link:

T i 1 i = [ cos θ i sin θ i cos α i sin θ i sin α i a i cos θ i sin θ i cos θ i cos α i cos θ i cos α i a i sin θ i 0 sin α i cos α i d i 0 0 0 1 ]
where ci = cos(θi) and si = sin(θi).

T 0 1 = [ c 1 0 s 1 0 s 1 0 c 1 0 0 1 0 0 0 0 0 1 ] , T 1 2 = [ c 2 s 2 0 l 1 c 2 s 2 c 2 0 l 1 s 2 0 0 1 d 1 0 0 0 1 ] , T 2 3 = [ c 3 s 3 0 l 2 c 3 s 3 c 3 0 l 2 s 3 0 0 1 0 0 0 0 1 ] T 3 4 = [ c 4 s 4 0 l 3 c 4 s 4 c 4 0 l 3 s 4 0 0 1 0 0 0 0 1 ] , T 4 5 = [ c 5 s 5 0 l 4 c 5 s 5 c 5 0 l 4 s 5 0 0 1 0 0 0 0 1 ] , T 5 6 = [ c 6 s 6 0 l 5 c 6 s 6 c 6 0 l 5 s 6 0 0 1 0 0 0 0 1 ]

Finally we obtain the product of all six link transforms:

T 0 6 = T 0 1 T 1 2 T 2 3 T 3 4 T 4 5 T 5 6

3. Dynamics of the Manipulator

In this section, the torque of each joint is calculated. To show the effectiveness of the proposed manipulator, the joint torques are calculated using the proposed manipulator (using three motors only) and the conventional manipulators (a motor for each joint).

Let us assume for concreteness that the center of mass of each link is at its geometric center. For the manipulator used in our experiments, the mass of links without the motors are as follow: ml1 = 760 gm, ml2 = 720 gm, ml3 = 680 gm, ml4 = 640 gm, and finally ml5 = 600 gm. These masses are calculated for the manipulator with l1 = 19 cm, l2 = 18 cm, l3 = 17 cm, l4 = 16 cm, l5 = 15 cm and d2 = 21 cm.

The mass of each motor is 1,500 gm; for the manipulator of the proposed design, the first motor and the second motor are located on the base and not on the links themselves. Therefore, for our manipulator, the mass of the first link will be equal to the mass of this link (760 gm) plus the mass of the motor (1,500 gm) that controls the next links. Because there are no more motors, the mass of the links will be: m1 = 2,260 gm, m2 = 720 gm, m3 = 680 gm, m4 = 640 gm, and m5 = 600 gm. Figure 9(a) shows the mass of each link with its motor for the manipulator of the proposed design.

For the conventional three dimensional planar manipulator (one motor for each link), the mass of the first link will equal to the mass of link itself plus the mass of the motor which controls the second link position, i.e., 760 + 1,500 gm. The mass of the second link will equal to the mass of link itself plus the mass of the motor which controls the third link position, i.e., 720 + 1,500 gm. The mass of the third link will equal to the mass of third link plus the mass of the motor which controls the fourth link position, i.e., 680 + 1,500 gm. The mass of the fourth link will equal to the mass of fourth link plus the mass of the motor which controls the fifth link position, i.e., 640 + 1,500 gm, while the mass of the last link will equal to the mass of the link itself because there are no more motors, i.e., 600 gm. Figure 9(b) shows the mass of each link using the manipulator with five motors, while Table 2 shows the values of mass of the links using both the manipulator with two motors and the manipulator with five links.

It is clearly noted how the proposed method could be used to decrease the weight of manipulator. Decreasing the weight leads to a decrease of the torques of each link. The next section shows the results of the torques of each joint when the end-effector is following a desired path. Using the Lagrangian formulation, the dynamical equations of motion of the manipulator is:

j = 1 6 M i j q ¨ j + V i + G i = Q i
for i = 1,2,….,6.

The first term in this equation is the inertia forces, the second term represents the Coriolis and centrifugal forces, and the third term gives the gravitational effects [1,20,21]. Dynamics equations of the manipulator are discussed in details in the Appendix.

As shown by dynamics equations, increasing the weight of motors will increase the torques needed to control the manipulator. In order to decrease the effect of the motors weight on the inertia of manipulators, parallel manipulators are used, as we mentioned earlier. For example in reference [22], the parallel manipulator is actuated by three servo-motors located at the base which contributes to reducing the inertia of manipulators. Reference [23] shows another way to decrease the effect of the motors weight on the inertia of manipulators. This reference shows a simple configuration design, which comprises of only three joints: two at the shoulder and one at the hand. In this design, the moment of inertia of the arm is constant and independent from the joint angles. In contrast for our manipulator, we see from Equations A21A25 that the moment of inertia value is dependent on the joint angles.

4. Simulation Results

This section shows the effectiveness of using the proposed manipulator to be used when it is desired to make the end-effector follow a desired path. This section has two examples. The first example calculates the torques using both manipulators (the proposed one and the conventional three dimensional planar manipulator) and shows how effective the proposed manipulator is in decreasing the torque of each joint required to move the manipulator. To verify the estimation results and compare between them and the results measured from the manipulator itself, the second example has been shown. This example shows the results if the torque using: (1) the conventional three dimensional planar manipulator with defined desired joint angles path, (2) the proposed manipulator with the defined desired joint angles path and finally (3) the proposed manipulator with the measured joint angles path when the joint angles follow the desired joint angles path.

Case One

Torque of each joint for both the manipulators is calculated to show the effectiveness of using the proposed manipulator to decrease the torque of each joint. Using the same manipulator with l = [19,18,17,16,15]T, and d2 = 21 where all lengths are in cm, the joint angles path is defined as:

θ 1 ( t ) = 0.5 cos ( 4 t )
θ 2 ( t ) = cos ( 2 t ) + 1
θ 3 ( t ) = 4 cos ( t ) + 3

It should be remembered that when using the proposed manipulator, θ3, θ4, θ5, θ6 are equal. To show the effectiveness of the proposed manipulator in decreasing the torque, Figure 10 shows the values of the torques of the first joint using both the manipulators, the proposed manipulator (with three motors) and the manipulator of six motors.

Figure 11 shows the values of the torques of the second joint using the manipulators.

Figure 12 shows the absolute values of the torques of the third joint while Figure 13 shows the absolute values of the torques of the fourth joint using the both manipulators.

Figure 14 shows the torques of the fifth joint and finally Figure 15 shows the torques of the sixth joint angle using the both manipulators.

First of all, it is noted that the torque of the sixth joint has the same value using both the manipulators because the sixth link has the same mass for both the manipulators, in other words the mass of the sixth link is equal to the mass of the link itself only because it does not hold any motor.

Secondly, as mentioned earlier for the proposed manipulator, the third motor should balance the torque of all the third, fourth, fifth and the sixth joint. In other words, the torque of the third motor should equal to (T3 + T4 + T5 + T6) for the proposed manipulator. Figure 16 shows the power that the third motor should balance for both the manipulators. It is noted from this example that using the proposed manipulator not only decreases the number of motors used in the manipulator, but also decreases the torques of the motors used to control it.

Case Two

The trajectory applied to robot in verification experiments in this case is:

θ 1 ( t ) = e 0.7 t 25
θ 2 ( t ) = e 0.5 t 3
θ 3 ( t ) = e 0.9 t 150

Figure 17 shows the estimated (white) and measured (red) angle, angular velocity and angular acceleration of the first joint angle defined above. Figure 18 shows the estimated (white) and measured (red) angle, angular velocity and angular acceleration of the second joint angle.

Figure 19 shows the estimated (white) and measured (red) angle, angular velocity and angular acceleration of the third joint angle of the manipulator. It should be remembered again that using the proposed manipulator, θ3, θ4, θ5, θ6 are equals.

Figures 2025 show the comparison between the torque of each joint angle for: (1) the conventional three dimensional planar manipulator using the estimated joint angles path; (2) the proposed manipulator using the estimated joint angles path; and finally (3) the proposed manipulator using the measured angular position, velocity and acceleration of the manipulator joints.

The results obtained from verification experiments indicate that there is a good agreement between the torque of the joint angles for the proposed manipulator using the estimated joint angles path (green) and the measured joint angles path (red). These figures show the effectiveness of the proposed manipulator in decreasing the torque of the joint angles using the proposed manipulator.

As mentioned in the first example that for the proposed manipulator, the third motor should balance the torque of all the third, fourth, fifth and the sixth joint, i.e., the torque of the third motor should equal (T3 + T4 + T5 + T6) for the proposed manipulator, Figure 26 shows that even though that this motor (third motor) should balance the torques of four links, this motor could be smaller in size (less power) in the proposed manipulator than the third motor in the conventional three dimensional planar manipulator.

5. Conclusions

This paper presents a mechanical design for a three dimensional planar redundant manipulator. Theoretically, for each degree of freedom there should be one motor. However, in this design only three motors are needed to control any n degrees of freedom three dimensional planar redundant manipulator. Therefore, this design can be used to decrease the weight of the manipulator significantly. The design steps of this manipulator are explained in detail. The dynamical equations are calculated for both the proposed and the conventional three dimensional planar manipulators (with n motors) and it is concluded from the result, that even though the proposed manipulator has less motors, these motors could be even smaller (as regard to power) than the motors used with conventional three dimensional planar manipulators.

Appendix

This section explains the Dynamics of the manipulator used in experiments. In Equation (10):

M = i = 1 6 ( J υ i T m i J υ i + J w i T I i J w i )
where Jυi is the Jacobian submatrix assosiated with the linear velocity of the center of mass of link i and Jwi is the Jacobian submatrix assosiated with the angular velocity vector of link i.

J υ i = [ J υ i 1 , J υ i 2 , , J υ i i , 0 , 0 , , 0 ]
J w i = [ J w i 1 , J w i 2 , , J w i i , 0 , 0 , , 0 ]
where J υ i j and J w i j are the j th column vectors of Jυi and Jwi, respectively. Applying the theory of instantaneous screw motion [1] for ji, and because all the joints are revolute joints, we obtain:
J υ i j = z j 1 × j 1 p c i
J w i j = z j 1

The link Jacobian submatrices, Jυi for the proposed manipulator can be obtained as:

J υ 1 = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]

Jυ1 is a zero matrix because from the link parameters of the manipulator of Table 1, it is noted that a1 and d1 are equal to zero because the first motor is set on the base.

J υ 2 = [ d 1 c 1 l 1 2 s 1 c 2 l 1 2 c 1 s 2 0 0 0 0 d 1 s 1 l 1 2 c 1 c 2 l 1 2 s 1 s 2 0 0 0 0 0 l 1 2 c 2 0 0 0 0 ]
J υ 3 = [ d 1 c 1 s 1 ( l 1 c 2 + l 2 2 c 23 ) c 1 ( l 1 s 2 + l 2 2 s 23 ) c 1 ( l 2 2 s 23 ) 0 0 0 d 1 s 1 + c 1 ( l 1 c 2 + l 2 2 c 23 ) s 1 ( l 1 s 2 + l 2 2 s 23 ) s 1 ( l 2 2 s 23 ) 0 0 0 0 l 1 c 2 + l 2 2 c 23 l 2 2 c 23 0 0 0 ]
J υ 4 = [ d 1 c 1 s 1 ( l 1 c 2 + l 2 c 23 + l 3 2 c 234 ) c 1 ( l 1 s 2 + l 2 s 23 + l 3 2 s 234 ) c 1 ( l 2 s 23 + l 3 2 s 234 ) c 1 ( l 3 2 s 234 ) 0 0 d 1 s 1 + c 1 ( l 1 c 2 + l 2 c 23 + l 3 2 c 234 ) s 1 ( l 1 s 2 + l 2 s 23 + l 3 2 s 234 ) s 1 ( l 2 s 23 + l 3 2 s 234 ) s 1 ( l 3 2 s 234 ) 0 0 0 l 1 c 2 + l 2 c 23 + l 3 2 c 234 l 2 c 23 + l 3 2 c 234 l 3 2 c 234 0 0 ]
J υ 5 = [ d 1 c 1 s 1 ( l 1 c 2 + l 2 c 23 + l 3 c 234 + l 4 2 c 2345 ) c 1 ( l 1 s 2 + l 2 s 23 + l 3 s 234 + l 4 2 s 2345 ) c 1 ( l 2 s 23 + l 3 s 234 + l 4 2 s 2345 ) c 1 ( l 3 s 234 + l 4 2 s 2345 ) c 1 ( l 4 2 s 2345 ) 0 d 1 s 1 + c 1 ( l 1 c 2 + l 2 c 23 + l 3 c 234 + l 4 2 c 2345 ) s 1 ( l 1 s 2 + l 2 s 23 + l 3 s 234 + l 4 2 s 2345 ) s 1 ( l 2 s 23 + l 3 s 234 + l 4 2 s 2345 ) s 1 ( l 3 s 234 + l 4 2 s 2345 ) s 1 ( l 4 2 s 2345 ) 0 0 l 1 c 2 + l 2 c 23 + l 3 c 234 + l 4 2 c 2345 l 2 c 23 + l 3 c 234 + l 4 2 c 2345 l 3 c 234 + l 4 2 c 2345 l 4 2 c 2345 0 ]
J υ 6 = [ d 1 c 1 s 1 ( l 1 c 2 + l 2 c 23 + l 3 c 234 + l 4 c 2345 + l 5 2 c 23456 ) c 1 ( l 1 s 2 + l 2 s 23 + l 3 s 234 + l 4 s 2345 + l 5 2 s 23456 ) c 1 ( l 2 s 23 + l 3 s 234 + l 4 s 2345 + l 5 2 s 23456 ) c 1 ( l 3 s 234 + l 4 s 2345 + l 5 2 s 23456 ) c 1 ( l 4 s 2345 + l 5 2 s 23456 ) c 1 ( l 5 2 s 23456 ) d 1 s 1 + c 1 ( l 1 c 2 + l 2 c 23 + l 3 c 234 + l 4 c 2345 + l 5 2 c 23456 ) s 1 ( l 1 s 2 + l 2 s 23 + l 3 s 234 + l 4 s 2345 + l 5 2 s 23456 ) s 1 ( l 2 s 23 + l 3 s 234 + l 4 s 2345 + l 5 2 s 23456 ) s 1 ( l 3 s 234 + l 4 s 2345 + l 5 2 s 23456 ) s 1 ( l 4 s 2345 + l 5 2 s 23456 ) s 1 ( l 5 2 s 23456 ) 0 l 1 c 2 + l 2 c 23 + l 3 c 234 + l 4 c 2345 + l 5 2 c 23456 l 2 c 23 + l 3 c 234 + l 4 c 2345 + l 5 2 c 23456 l 3 c 234 + l 4 c 2345 + l 5 2 c 23456 l 4 c 2345 + l 5 2 c 23456 l 5 2 c 23456 ]

To formulate the submatrices of Jwi, following equations are used:

J w 1 = [ 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 ]
J w 2 = [ 0 s 1 0 0 0 0 0 c 1 0 0 0 0 1 0 0 0 0 0 ]
J w 3 = [ 0 s 1 s 1 0 0 0 0 c 1 c 1 0 0 0 1 0 0 0 0 0 ]
J w 4 = [ 0 s 1 s 1 s 1 0 0 0 c 1 c 1 c 1 0 0 1 0 0 0 0 0 ]
J w 5 = [ 0 s 1 s 1 s 1 s 1 0 0 c 1 c 1 c 1 c 1 0 1 0 0 0 0 0 ]
J w 6 = [ 0 s 1 s 1 s 1 s 1 s 1 0 c 1 c 1 c 1 c 1 c 1 1 0 0 0 0 0 ]

Now to formulate the inertia matrices of the manipulator, assuming that the moving links are homogeneous with relatively small cross section, the inertia matrix iIi of link i about its center of mass and expressed in the link link frame i th link fram is:

I i i = 1 12 m i a i 2 [ 0 0 0 0 1 0 0 0 1 ]
for I = 1,2,…,6.

Now to find Ii inertia matrix of link i about its center of mass and expressed in the base link frame, the following equation is used:

I i = R 0 i I i i ( R 0 i ) T

Using Equations (A18) and (A19), we obtain:

I 1 = [ 0 0 0 0 0 0 0 0 0 ]

I1 is a zero matrix because a1 is equal to zero.

I 2 = m 2 l 1 2 12 [ s 1 2 + c 1 2 c 2 2 s 1 c 1 ( s 2 2 1 ) c 1 ( s 2 c 2 ) s 1 c 1 ( s 2 2 1 ) c 1 2 + s 1 2 s 2 2 s 1 ( s 2 c 2 ) c 1 ( s 2 c 2 ) s 1 ( s 2 c 2 ) c 2 2 ]
I 3 = m 3 l 2 2 12 [ s 1 2 + c 1 2 c 23 2 s 1 c 1 ( s 23 2 1 ) c 1 ( s 23 c 23 ) s 1 c 1 ( s 23 2 1 ) c 1 2 + s 1 2 s 23 2 s 1 ( s 23 c 23 ) c 1 ( s 23 c 23 ) s 1 ( s 23 c 23 ) c 23 2 ]
I 4 = m 4 l 3 2 12 [ s 1 2 + c 1 2 c 234 2 s 1 c 1 ( s 234 2 1 ) c 1 ( s 234 c 234 ) s 1 c 1 ( s 234 2 1 ) c 1 2 + s 1 2 s 234 2 s 1 ( s 234 c 234 ) c 1 ( s 234 c 234 ) s 1 ( s 234 c 234 ) c 234 2 ]
I 5 = m 5 l 4 2 12 [ s 1 2 + c 1 2 c 2345 2 s 1 c 1 ( s 2345 2 1 ) c 1 ( s 2345 c 2345 ) s 1 c 1 ( s 2345 2 1 ) c 1 2 + s 1 2 s 2345 2 s 1 ( s 2345 c 2345 ) c 1 ( s 2345 c 2345 ) s 1 ( s 2345 c 2345 ) c 2345 2 ]
I 6 = m 6 l 5 2 12 [ s 1 2 + c 1 2 c 23456 2 s 1 c 1 ( s 23456 2 1 ) c 1 ( s 23456 c 23456 ) s 1 c 1 ( s 23456 2 1 ) c 1 2 + s 1 2 s 23456 2 s 1 ( s 23456 c 23456 ) c 1 ( s 23456 c 23456 ) s 1 ( s 23456 c 23456 ) c 23456 2 ]

Now to find the inertia manipulator matrix M, Equation (A1) can be rewritten as:

M = J υ 1 T m 1 J υ 1 + J w 1 T I 1 J w 1 + J υ 2 T m 2 J υ 2 + J w 2 T I 2 J w 2 + J υ 3 T m 3 J υ 3 + J w 3 T I 3 J w 3 + J υ 4 T m 4 J υ 4 + J w 4 T I 4 J w 4 + J υ 5 T m 5 J υ 5 + J w 5 T I 5 J w 5 + J υ 6 T m 6 J υ 6 + J w 6 T I 6 J w 6

Therefore, Equation (10) can be rewritten as:

[ τ 1 τ 2 τ 3 τ 4 τ 5 τ 6 ] = [ M 11 M 12 M 13 M 14 M 15 M 16 M 21 M 22 M 23 M 24 M 25 M 26 M 31 M 32 M 33 M 34 M 35 M 36 M 41 M 42 M 43 M 44 M 45 M 46 M 51 M 52 M 53 M 54 M 55 M 56 M 61 M 62 M 63 M 64 M 65 M 66 ] [ θ ¨ 1 θ ¨ 2 θ ¨ 3 θ ¨ 4 θ ¨ 5 θ ¨ 6 ] + [ V 1 ( θ , θ ¨ ) V 2 ( θ , θ ¨ ) V 3 ( θ , θ ¨ ) V 4 ( θ , θ ¨ ) V 5 ( θ , θ ¨ ) V 6 ( θ , θ ¨ ) ] + [ G 1 ( θ ) G 2 ( θ ) G 3 ( θ ) G 4 ( θ ) G 5 ( θ ) G 6 ( θ ) ]
where:
M 11 = 1 12 ( 4 m 2 ( 3 d 1 2 + l 1 2 cos 2 ( θ 2 ) ) + 4 m 3 ( 3 d 1 2 + 3 l 1 2 cos 2 ( θ 2 ) + l 2 2 cos 2 ( θ 2 + θ 3 ) + 3 l 1 l 2 cos ( θ 2 ) cos ( θ 2 + θ 3 ) ) + 12 m 4 ( d 1 2 + 1 4 ( 2 l 1 cos ( θ 2 ) + 2 l 2 cos ( θ 2 + θ 3 ) + l 3 cos ( θ 2 + θ 3 + θ 4 ) ) 2 + 12 m s ( d 1 2 + 1 4 ( 2 l 1 cos ( θ 2 ) + 2 l 2 cos ( θ 2 + θ 3 ) + 2 l 3 cos ( θ 2 + θ 3 + θ 4 ) + l 4 cos ( θ 2 + θ 3 + θ 4 + θ 5 ) ) 2 ) + 12 m 6 ( d 1 2 + 1 4 ( 2 ( l 1 cos ( θ 2 ) + l 2 cos ( θ 2 + θ 3 ) + l 3 cos ( θ 2 + θ 3 + θ 4 ) + l 4 cos ( θ 2 + θ 3 + θ 4 + θ 5 ) ) + l 5 cos ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) 2 ) + l 3 2 m 4 cos 2 ( θ 2 + θ 3 + θ 4 ) + l 4 2 m 5 cos 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) + l 5 2 m 6 cos 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) )
M 12 = M 21 = 1 2 d 1 ( l 3 m 4 sin ( θ 2 + θ 3 + θ 4 ) + 2 l 3 m 5 sin ( θ 2 + θ 3 + θ 4 ) + l 4 m 5 sin ( θ 2 + θ 3 + θ 4 + θ 5 ) + m 6 ( 2 l 3 sin ( θ 2 + θ 3 + θ 4 ) + 2 l 4 sin ( θ 2 + θ 3 + θ 4 + θ 5 ) + l 5 sin ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) + l 2 ( m 2 + 2 ( m 4 + m 5 + m 6 ) ) sin ( θ 2 + θ 3 ) + l 1 ( m 3 + 2 ( m 3 + m 4 + m 5 + m 6 ) ) sin ( θ 2 ) )
M 13 = M 31 = 1 2 d 1 ( l 4 m 5 sin ( θ 2 + θ 3 + θ 4 + θ 5 ) + m 6 ( 2 l 4 sin ( θ 2 + θ 3 + θ 4 + θ 5 ) + l 5 sin ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) sin ( θ 2 + θ 3 + θ 4 ) + l 2 ( m 3 + 2 ( m 4 + m 5 + m 6 ) ) sin ( θ 2 + θ 3 ) )
M 14 = M 41 = 1 2 d 1 ( l 5 m 6 sin ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) sin ( θ 2 + θ 3 + θ 4 + θ 5 ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) sin ( θ 2 + θ 3 + θ 4 ) )
M 15 = M 51 = 1 2 d 1 ( l 5 m 6 sin ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) sin ( θ 2 + θ 3 + θ 4 + θ 5 ) )
M 16 = 1 2 d 1 l 5 m 6 sin ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 )
M 22 = 1 192 ( 192 l 1 ( l 4 m 5 cos ( θ 3 + θ 4 + θ 5 ) + m 6 ( 2 l 4 cos ( θ 3 + θ 4 + θ 5 ) + l 5 cos ( θ 3 + θ 4 + θ 5 + θ 6 ) ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 3 + θ 4 ) + l 2 ( m 3 + 2 ( m 4 + m 5 + m 6 ) ) cos ( θ 3 ) ) + 2 l 3 2 m 4 cos ( 2 ( θ 2 + θ 3 + θ 4 ) ) l 3 2 m 4 cos ( 2 ( 2 θ 1 + θ 2 + θ 3 + θ 4 ) ) l 3 2 m 4 cos ( 4 θ 1 2 ( θ 2 + θ 3 + θ 4 ) ) + 2 l 4 2 m 5 cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) l 4 2 m 5 cos ( 2 ( 2 θ 1 + θ 2 + θ 3 + θ 4 + θ 5 ) ) l 4 2 m 5 cos ( 4 θ 1 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) + 192 l 3 l 4 m 5 cos ( θ 5 ) + 2 m 6 ( l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 96 l 5 ( l 3 cos ( θ 5 + θ 6 ) + l 4 cos ( θ 6 ) ) + 96 ( 2 l 4 l 3 cos ( θ 5 ) + l 3 2 + l 4 2 ) ) + 192 l 2 ( l 5 m 6 cos ( θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 4 + θ 5 ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 4 ) ) + 4 l 1 2 ( m 2 ( sin 2 ( 2 θ 1 ) cos ( 2 θ 2 ) + 16 ) + 48 ( m 3 + m 4 + m 5 + m 6 ) ) + 4 l 2 2 ( m 3 ( sin 2 ( 2 θ 1 ) cos ( 2 ( θ 2 + θ 3 ) ) + 16 ) + 48 ( m 4 + m 5 + m 6 ) ) + 64 l 3 2 m 4 + 192 l 3 2 m 5 + 64 l 4 2 m 5 )
M 23 = M 32 = 1 192 ( 192 l 2 ( l 5 m 6 cos ( θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 4 + θ 5 ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 4 ) ) + 2 l 3 2 m 4 cos ( 2 ( θ 2 + θ 3 + θ 4 ) ) l 3 2 m 4 cos ( 2 ( 2 θ 1 + θ 2 + θ 3 + θ 4 ) ) l 3 2 m 4 cos ( 4 θ 1 2 ( θ 2 + θ 3 + θ 4 ) ) + 2 l 4 2 m 5 cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) l 4 2 m 5 cos ( 2 ( 2 θ 1 + θ 2 + θ 3 + θ 4 + θ 5 ) ) l 4 2 m 5 cos ( 4 θ 1 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) + 192 l 3 l 4 m 5 cos ( θ 5 ) + 2 m 6 ( l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 96 l 5 ( l 3 cos ( θ 5 + θ 6 ) + l 4 cos ( θ 6 ) ) + 96 ( 2 l 4 l 3 cos ( θ 5 ) + l 3 2 + l 4 2 ) ) + 96 l 1 ( l 4 m 5 cos ( θ 3 + θ 4 + θ 5 ) + m 6 ( 2 l 4 cos ( θ 3 + θ 4 + θ 5 ) + l 5 cos ( θ 3 + θ 4 + θ 5 + θ 6 ) ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 3 + θ 4 ) + l 2 ( m 3 + 2 ( m 4 + m 5 + m 6 ) ) cos ( θ 3 ) ) + 4 l 2 2 ( m 3 ( sin 2 ( 2 θ 1 ) cos ( 2 ( θ 2 + θ 3 ) ) + 16 ) + 48 ( m 4 + m 5 + m 6 ) ) + 64 l 3 2 m 4 + 192 l 3 2 m 5 + 64 l 4 2 m 5 )
M 24 = M 42 = 1 192 ( 2 l 3 2 m 4 cos ( 2 ( θ 2 + θ 3 + θ 4 ) ) l 3 2 m 4 cos ( 2 ( 2 θ 1 + θ 2 + θ 3 + θ 4 ) ) l 3 2 m 4 cos ( 4 θ 1 2 ( θ 2 + θ 3 + θ 4 ) ) + 192 l 4 l 3 m 5 cos ( θ 5 ) + 2 l 4 2 m 5 cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) l 4 2 m 5 cos ( 2 ( 2 θ 1 + θ 2 + θ 3 + θ 4 + θ 5 ) ) l 4 2 m 5 cos ( 4 θ 1 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) + 2 m 6 ( l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 96 l 5 ( l 3 cos ( θ 5 + θ 6 ) + l 4 cos ( θ 6 ) ) + 96 ( 2 l 4 l 3 cos ( θ 5 ) + l 3 2 + l 4 2 ) ) + 96 l 2 ( l 5 m 6 cos ( θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 4 + θ 5 ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 4 ) ) + 96 l 1 ( l 5 m 6 cos ( θ 3 + θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 3 + θ 4 + θ 5 ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 3 + θ 4 ) ) + 64 l 3 2 m 4 + 192 l 3 2 m 5 + 64 l 4 2 m 5 )
M 25 = M 52 = 1 96 ( l 4 m 5 ( 48 ( l 1 cos ( θ 3 + θ 4 + θ 5 ) + l 2 cos ( θ 4 + θ 5 ) + l 3 cos ( θ 5 ) ) + l 4 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) ) ) + m 6 ( l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 48 l 5 ( l 1 cos ( θ 3 + θ 4 + θ 5 + θ 6 ) + l 2 cos ( θ 4 + θ 5 + θ 6 ) + l 3 cos ( θ 5 + θ 6 ) + 2 l 4 cos ( θ 6 ) ) + 96 l 4 ( l 1 cos ( θ 3 + θ 4 + θ 5 ) + l 2 cos ( θ 4 + θ 5 ) + l 3 cos ( θ 5 ) + l 4 ) ) )
M 26 = M 62 = 1 96 l 5 m 6 ( 48 ( l 1 cos ( θ 3 + θ 4 + θ 5 + θ 6 ) + l 2 cos ( θ 4 + θ 5 + θ 6 ) + l 3 cos ( θ 5 + θ 6 ) + l 4 cos ( θ 6 ) ) + l 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) )
M 33 = 1 96 ( 96 l 2 ( l 5 m 6 cos ( θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 4 + θ 5 ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 4 ) ) + l 4 2 m 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) ) + m 6 ( 96 l 5 l 4 cos ( θ 6 ) + l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 96 l 4 2 ) + l 3 2 ( m 4 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 ) ) ) + 96 ( m 5 + m 6 ) ) + 96 l 3 ( l 5 m 6 cos ( θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 5 ) ) + l 2 2 ( 2 m 3 ( sin 2 ( 2 θ 1 ) cos ( 2 ( θ 2 + θ 3 ) ) + 16 ) + 96 ( m 4 + m 5 + m 6 ) ) )
M 34 = M 43 = 1 96 ( l 3 2 ( m 4 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 ) ) ) + 96 ( m 5 + m 6 ) ) + 96 l 3 ( l 5 m 6 cos ( θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 5 ) ) + l 4 2 m 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( θ 2 + θ 3 + θ 4 + θ 5 ) ) ) + m 6 ( 96 l 5 l 4 cos ( θ 6 ) + l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 96 l 4 2 ) + 48 l 2 ( l 5 m 6 cos ( θ 4 + θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 4 + θ 5 ) + l 3 ( m 4 + 2 ( m 5 + m 6 ) ) cos ( θ 4 ) ) )
M 35 = M 53 = 1 96 ( l 4 m 5 ( 48 l 2 cos ( θ 4 + θ 5 ) + 48 l 3 cos ( θ 5 ) + l 4 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) ) ) + m 6 ( l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 48 l 5 ( l 2 cos ( θ 4 + θ 5 + θ 6 ) + l 3 cos ( θ 5 + θ 6 ) + 2 l 4 cos ( θ 6 ) ) + 96 l 4 ( l 2 cos ( θ 4 + θ 5 ) + l 3 cos ( θ 5 ) + l 4 ) ) )
M 36 = M 63 = 1 96 l 5 m 6 ( 48 ( l 2 cos ( θ 4 + θ 5 + θ 6 ) + l 3 cos ( θ 5 + θ 6 ) + l 4 cos ( θ 6 ) ) + l 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) )
M 44 = 1 96 ( l 3 2 ( m 4 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 ) ) ) + 96 ( m 5 + m 6 ) ) + 96 l 3 ( l 5 m 6 cos ( θ 5 + θ 6 ) + l 4 ( m 5 + 2 m 6 ) cos ( θ 5 ) ) + l 4 2 m 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) ) + m 6 ( 96 l 5 l 4 cos ( θ 6 ) + l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 96 l 4 2 ) )
M 45 = M 54 = 1 96 ( l 4 m 5 ( 48 l 3 cos ( θ 5 ) + l 4 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) ) ) + m 6 ( l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 48 l 5 ( l 3 cos ( θ 5 + θ 6 ) + 2 l 4 cos ( θ 6 ) ) + 96 l 4 ( l 3 cos ( θ 5 ) + l 4 ) ) )
M 46 = M 64 = 1 96 l 5 m 6 ( 48 l 3 cos ( θ 5 + θ 6 ) + 48 l 4 cos ( θ 6 ) + l 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) )
M 55 = 1 96 ( l 4 2 m 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 ) ) ) + m 6 ( 96 l 5 l 4 cos ( θ 6 ) + l 5 2 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) + 96 l 4 2 ) )
M 56 = M 65 = 1 96 l 5 m 6 ( 48 l 4 cos ( θ 6 ) + l 5 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) ) )
M 66 = 1 96 l 5 2 m 6 ( 32 ( cos ( 4 θ 1 ) 1 ) cos ( 2 ( θ 2 + θ 3 + θ 4 + θ 5 + θ 6 ) ) )

The second term (the Coriolis and centrifugal forces term) and third term (gravitational effects term) of Equation (A27) can be calculated as follow:

V i = j = 1 6 k = 1 6 ( M i j q k 1 2 M j k q i ) q ˙ j q ˙ k
G i = j = 1 6 m j g T J υ j i

References

  1. Tsai, L.-W. Robot Analysis: The Mechanics of Serial and Parallel Manipulators; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1999. [Google Scholar]
  2. Ma, S.; Yoshinada, H.; Hirose, S. CT ARM-I: Coupled tendon-driven manipulator model I-design and basic experiments. Proceeding of IEEE International Conference on Robotics and Automation, Nice, France, 12–14 May 1992; Volume 3. pp. 2094–2100.
  3. Ma, S.; Hirose, S.; Yoshinada, H. Design and experiments for a coupled tendon-driven manipulator. IEEE Control Syst. Mag. 1993, 13, 30–36. [Google Scholar]
  4. Feliu, V.; Ramos, F. Strain gauge based control of single-link flexible very lightweight robots robust to payload changes. Mechatronics 2005, 15, 547–571. [Google Scholar]
  5. Albu-Schaffer, A.; Haddadin, S.; Ott, Ch.; Stemmer, A.; Wimböck, T.; Hirzinger, G. The DLR lightweight robot: Design and control concepts for robots in human environments. Ind. Robot Int. J. 2007, 34, 376–385. [Google Scholar]
  6. Hagn, U.; Nickl, M.; Jorg, S.; Passig, G.; Bahls, T.; Nothhelfer, A.; Hacker, F.; Le-Tien, L.; Albu-Schaffer, A.; Konietschke, R.; et al. The DLR MIRO: A versatile lightweight robot for surgical applications. Ind. Robot: Int. J. 2008, 35, 324–336. [Google Scholar]
  7. Hirose, S.; Ma, S. Moray drive for multijoint. Proceeding of 5th International Conference of Advanced Robotics, Pisa, Italy, 19–22 June 1991; Volume 1. pp. 521–526.
  8. Merlet, J.P. Parallel Robots, 2nd ed.; Springer: Dordrecht, The Netherlands, 2006. [Google Scholar]
  9. Londi, F.; Pennestri, E.; Valentini, P.P.; Vita, L. Control and virtual reality simulation of tendon driven mechanisms. Multibody Syst. Dyn. 2004, 12, 133–145. [Google Scholar]
  10. Rodriguez-Donate, C.; Osornio-Rios, R.A.; Rivera-Guillen, J.R.; de Jesus Romero-Troncoso, R. Fused smart sensor network for multi-axis forward kinematics estimation in industrial robots. Sensors 2011, 11, 4335–4357. [Google Scholar]
  11. Lee, J.K.; Kim, K.; Lee, Y.; Jeong, T. Simultaneous intrinsic and extrinsic parameter identification of a hand-mounted laser-vision sensor. Sensors 2011, 11, 8751–8768. [Google Scholar]
  12. Yahya, S.; Moghavvemi, M.; Mohamed, H.A.F. Geometrical approach of planar hyper-redundant manipulators: Inverse kinematics, path planning and workspace. Simul. Model. Pract. Theory 2011, 19, 406–422. [Google Scholar]
  13. Mohamed, H.A.F.; Yahya, S.; Moghavvemi, M.; Yang, S.S. A new inverse kinematics method for three dimensional redundant manipulators. ICROS-SICE, Fukuoka, Japan, 18–21 August 2009; pp. 1557–1562.
  14. Yahya, S.; Moghavvemi, M.; Mohamed, H.A.F. Singularity avoidance of a six degrees of freedom three dimensional redundant planar manipulator. Comput. Math. Appl. 2012. [Google Scholar] [CrossRef]
  15. Yahya, S.; Moghavvemi, M.; Mohamed, H.A.F. A review of Singularity avoidance in the inverse kinematics of redundant robot manipulators. Int. Rev. Autom. Control 2011, 4, 807–814. [Google Scholar]
  16. Yahya, S.; Moghavvemi, M.; Mohamed, H.A.F. Redundant manipulators kinematics inversion. Sci. Res. Essays 2011, 6, 5462–5470. [Google Scholar]
  17. Khalil, W.; Dombre, E. Modeling, Identification & Control of Robots; Hermes Penton Ltd.: London, UK, 2002. [Google Scholar]
  18. Siciliano, B.; Sciavicco, L.; Villani, L.; Oriolo, G. Robotics: Modelling, Planning and Control; Springer-Verlag London Limited: London, UK, 2009. [Google Scholar]
  19. Siciliano, B.; Khatib, O. Springer Handbook of Robotics; Springer-Verlag: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  20. Angeles, J. Fundamentals of Robotic Mechanical Systems: Theory, Methods, and Algorithms, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  21. Spong, M.W.; Hutchinson, S.; Vidyasagar, M. Robot Dynamics and Control, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2004. [Google Scholar]
  22. Shang, W.W.; Cong, S. Nonlinear computed torque control for a high speed planar parallel manipulator. Mechatronics 2009, 19, 987–992. [Google Scholar]
  23. Yoshida, K.; Kurazume, R.; Umetani, Y. Torque optimization control in space robots with a redundant arm. Proceeding of IEEE/RSJ International Workshop on Intelligent Robots and Systems, Osaka, Japan, 3–5 November 1991; Volume: 3. pp. 1647–1652.
Figure 1. (a) A three dimensional planar redundant manipulator configuration; (b) A three dimensional planar redundant manipulator configuration using the method of [12,13].
Figure 1. (a) A three dimensional planar redundant manipulator configuration; (b) A three dimensional planar redundant manipulator configuration using the method of [12,13].
Sensors 12 06869f1 1024
Figure 2. The manipulator used in experiments [14]. The draft of the manipulator using the SolidWorks software (left). The mechanical design of the manipulator (right).
Figure 2. The manipulator used in experiments [14]. The draft of the manipulator using the SolidWorks software (left). The mechanical design of the manipulator (right).
Sensors 12 06869f2 1024
Figure 3. The design of the second joint angle (first link with second motor) of the manipulator [14]. The draft of the second joint angle using the SolidWorks software (left). The mechanical design of the second joint angle (right).
Figure 3. The design of the second joint angle (first link with second motor) of the manipulator [14]. The draft of the second joint angle using the SolidWorks software (left). The mechanical design of the second joint angle (right).
Sensors 12 06869f3 1024
Figure 4. The design of the third joint angle (second link with third motor) of the manipulator [14]. The draft of the third joint angle using the SolidWorks software (top left). The draft of the whole manipulator using the SolidWorks software (top right). The mechanical design of the third joint angle (bottom).
Figure 4. The design of the third joint angle (second link with third motor) of the manipulator [14]. The draft of the third joint angle using the SolidWorks software (top left). The draft of the whole manipulator using the SolidWorks software (top right). The mechanical design of the third joint angle (bottom).
Sensors 12 06869f4 1024
Figure 5. The design of the fourth joint angle (third link) of the manipulator [14]. The draft of the fourth joint angle using the SolidWorks software (top left). The draft of the whole manipulator using the SolidWorks software (top right). The mechanical design of the fourth joint angle (bottom).
Figure 5. The design of the fourth joint angle (third link) of the manipulator [14]. The draft of the fourth joint angle using the SolidWorks software (top left). The draft of the whole manipulator using the SolidWorks software (top right). The mechanical design of the fourth joint angle (bottom).
Sensors 12 06869f5 1024
Figure 6. The last joint of the manipulator.
Figure 6. The last joint of the manipulator.
Sensors 12 06869f6 1024
Figure 7. The mechanism of the first motor.
Figure 7. The mechanism of the first motor.
Sensors 12 06869f7 1024
Figure 8. The manipulator used in experiments.
Figure 8. The manipulator used in experiments.
Sensors 12 06869f8 1024
Figure 9. The position of mass for (a) the proposed manipulator; (b) the conventional manipulator.
Figure 9. The position of mass for (a) the proposed manipulator; (b) the conventional manipulator.
Sensors 12 06869f9 1024
Figure 10. The values of the torques of the first joint using the both manipulators.
Figure 10. The values of the torques of the first joint using the both manipulators.
Sensors 12 06869f10 1024
Figure 11. The values of the torques of the second joint using the both manipulators.
Figure 11. The values of the torques of the second joint using the both manipulators.
Sensors 12 06869f11 1024
Figure 12. The values of the torques of the third joint using the both manipulators.
Figure 12. The values of the torques of the third joint using the both manipulators.
Sensors 12 06869f12 1024
Figure 13. The values of the torques of the fourth joint using the both manipulators.
Figure 13. The values of the torques of the fourth joint using the both manipulators.
Sensors 12 06869f13 1024
Figure 14. The values of the torques of the fifth joint using the both manipulators.
Figure 14. The values of the torques of the fifth joint using the both manipulators.
Sensors 12 06869f14 1024
Figure 15. The values of the torques of the sixth joint using the both manipulators.
Figure 15. The values of the torques of the sixth joint using the both manipulators.
Sensors 12 06869f15 1024
Figure 16. The values of the torques of the third motor using the both manipulators.
Figure 16. The values of the torques of the third motor using the both manipulators.
Sensors 12 06869f16 1024
Figure 17. The values of the estimated and measured angular position, velocity and acceleration of the first joint angle (white: estimated, red: measured).
Figure 17. The values of the estimated and measured angular position, velocity and acceleration of the first joint angle (white: estimated, red: measured).
Sensors 12 06869f17 1024
Figure 18. The values of the estimated and measured angular position, velocity and acceleration of the second joint angle (white: estimated, red: measured).
Figure 18. The values of the estimated and measured angular position, velocity and acceleration of the second joint angle (white: estimated, red: measured).
Sensors 12 06869f18 1024
Figure 19. The values of the estimated and measured angular position, velocity and acceleration of the third joint angle (white: estimated, red: measured).
Figure 19. The values of the estimated and measured angular position, velocity and acceleration of the third joint angle (white: estimated, red: measured).
Sensors 12 06869f19 1024
Figure 20. The torque of the first joint angle.
Figure 20. The torque of the first joint angle.
Sensors 12 06869f20 1024
Figure 21. The torque of the second joint angle.
Figure 21. The torque of the second joint angle.
Sensors 12 06869f21 1024
Figure 22. The torque of the third joint angle.
Figure 22. The torque of the third joint angle.
Sensors 12 06869f22 1024
Figure 23. The torque of the fourth joint angle.
Figure 23. The torque of the fourth joint angle.
Sensors 12 06869f23 1024
Figure 24. The torque of the fifth joint angle.
Figure 24. The torque of the fifth joint angle.
Sensors 12 06869f24 1024
Figure 25. The torque of the sixth joint angle.
Figure 25. The torque of the sixth joint angle.
Sensors 12 06869f25 1024
Figure 26. The torque of the third motor.
Figure 26. The torque of the third motor.
Sensors 12 06869f26 1024
Table 1. Link parameters of the manipulator.
Table 1. Link parameters of the manipulator.
iαadθ
19000θ1
20l1d1θ2
30l20θ3
40l30θ4
50l40θ5
60l50θ6
Table 2. The mass of links for both the proposed and conventional manipulators.
Table 2. The mass of links for both the proposed and conventional manipulators.
m1(gm)1,5001,500

m2(gm)2,2602,260
m3(gm)7202,220
m4(gm)6802,180
m5(gm)6402,140
m6(gm)600600

Share and Cite

MDPI and ACS Style

Yahya, S.; Moghavvemi, M.; Almurib, H.A.F. Joint Torque Reduction of a Three Dimensional Redundant Planar Manipulator. Sensors 2012, 12, 6869-6892. https://doi.org/10.3390/s120606869

AMA Style

Yahya S, Moghavvemi M, Almurib HAF. Joint Torque Reduction of a Three Dimensional Redundant Planar Manipulator. Sensors. 2012; 12(6):6869-6892. https://doi.org/10.3390/s120606869

Chicago/Turabian Style

Yahya, Samer, Mahmoud Moghavvemi, and Haider Abbas F. Almurib. 2012. "Joint Torque Reduction of a Three Dimensional Redundant Planar Manipulator" Sensors 12, no. 6: 6869-6892. https://doi.org/10.3390/s120606869

APA Style

Yahya, S., Moghavvemi, M., & Almurib, H. A. F. (2012). Joint Torque Reduction of a Three Dimensional Redundant Planar Manipulator. Sensors, 12(6), 6869-6892. https://doi.org/10.3390/s120606869

Article Metrics

Back to TopTop