Kinematics of Hex-Piderix - A Six-Legged Robot

Abstract

In this paper the kinematic analysis of a six-legged robot, hereafter named Hex-piderix, is carried out. A three revolute (3R) chain has been chosen for each limb in order to mimic the leg structure of an insect. The rotation matrix, with unitary vectors, and the Denavit-Hartenberg (D-H) conventions are used to find the pose of the thorax. The problem of inverse position is solved by geometrical analysis. The direct and inverse infinitesimal kinematics are obtained by the reciprocal screw theory, considering a suction cup attached to each leg and modelling it as a UP linkage. A numerical example of the thorax pose was made by solving the equations obtained from the direct position analysis. The equations of the inverse position analysis were solved to obtain the angles of the joints. Finally, the velocity values of the thorax obtained from the infinitesimal kinematics were validated by simulating the movements of Hex-piderix using specialized software.

Keywords

Hexapod Robot Kinematics Screw Theory Suction Cup

1. Introduction

Legged robots have better performance in rough terrain than robots with wheels, thanks to the possibility of coordinated movements that provide greater flexibility and adaptability.

Legged robots can be classified as: bipeds, quadrupeds, hexapods and octopods. Most climbing robots are hexapods or quadrupeds - the former providing better static stability than the latter. They can perform static gaits by supporting the robot's body on five legs at any time, while quadrupeds can only walk steadily on a minimum of three legs. This feature makes hexapods much more stable than quadrupeds, since they can use a bigger support polygon. Notice that stability is a fundamental issue for mobile robots and speed is an important factor for robot locomotion - a quadruped robot is three times slower than a hexapod. In terms of reliability, the hexapod can continue walking even when one of its legs fails [1].

To obtain a dynamic model and to design a control algorithm of legged robots, it is important to model the kinematic behaviour of the complex multi-legged robotic mechanism.

In [2], the inverse position problem for each end-effector (EE) of the legs of the quadruped robot CLIBO is solved by homogeneous transformation matrices. The same approach is presented in [3] for a hexapod robot assuming forward motion at a constant velocity (v) in a straight path on a flat surface with alternating tripod gait. In this case, the thorax is held at a constant height and kept parallel to the ground plane during motion, and the thorax gravity centre is assumed to be located at its geometric centre. In [4], the problem of a hexapod robot is solved as in [3] but without any restriction.

A legged robot can be considered a partially parallel mechanism where the thorax is the mobile platform. In [5], the direct and inverse kinematics of the mobile platform of the quadruped robot Kamambaré, are solved by modelling it as a parallel robot. The analysis is performed using the D-H parameterization, starting at the surface and advancing towards the platform.

The work presented in [6] solves the kinematics of a four-legged robot composed by GZ-I modules, using kinematics graph theory to obtain robot posture. Two homogeneous transformation matrices describe the position of the robot, matrix R of rods and matrix J of joints. The joints are classified into horizontal and vertical joints depending on their axis of rotation.

The position problem of the tip leg has been studied for quadruped and hexapod robots [2 –6], however, the problem of the thorax pose has received little attention [5–6]. Furthermore, the infinitesimal kinematics of legged robots have been omitted.

One of the most important things in the study of legged robots is showing stability when walking. [7] conducted a study of impulse response and stability of a fractional oscillator, which determine the conditions for the fractional oscillator is strictly stable, unstable and marginally stable. Which determine the conditions to make it strictly stable, unstable and marginally stable.

In this paper, an analysis of the direct and inverse position of the hexapod robot called Hex-piderix is carried out, in addition, the infinitesimal kinematics are solved by using screw theory.

2. Description of the Hex-piderix robot

The basic configuration proposed for the robot Hex-piderix is shown in Fig. 1. The frame {0} represents the origin, while {Q} indicates the local frame attached to the robot thorax. The robot consists of six identical legs, L1, …, L6, that are distributed symmetrically around the thorax of the robot. Each limb has 3 degrees of freedom (DOF) that consist of three links, coxa, femur and tibia, connected by rotational joints. In the tip of the tibia there is a suction cup to adhere to vertical surfaces.

Figure 1.

Robot configuration.

3. Position analysis

The symmetry of the mechanism simplifies the position analysis. The inverse and direct kinematics are solved from geometrical analysis of the mechanism.

3.1 Forward position analysis

The aim of this analysis is to find the robot pose at the thorax level {Q} with respect to {0} for a given set of joint angles, q_i,j, where i represents the number of the limb (i = 1,…. 6) and j the number of the joint (j = 1, 2,3).

To determine the pose of the robot it is assumed that there are three legs touching the ground, L₁, L₃ and L₅. The corresponding coxa links, represented by the points P₁, P₃, and P₅, form an isosceles triangle composed by base a and sides b. The local frame {Q} is located at the triangle's centroid, its axis X_Q is directed to P₁ [8], as illustrated in Fig. 2.

To ensure the stability of the robot, the projection of the centre of gravity must fall within the polygon formed by the tips of the limbs.

Figure 2.

Triangle formed by coxa links.

The unitary vectors corresponding to {Q} were calculated as follows:

{\hat{X}}_{Q} = \frac{\vec{P} 1 - \vec{Q}}{| | \bar{P Q} | |},

(1)

{\hat{Y}}_{Q} = \frac{\vec{P} 2 - \vec{P} 3}{∥ \vec{P} 2 - \vec{P} 3 ∥},

(2)

{\hat{Z}}_{Q} = {\hat{X}}_{Q} \otimes {\hat{Y}}_{Q}

(3)

where, $\vec{Q}$ is the vector that includes the coordinates of the centroid.

The rotation matrix of the centroid is obtained with the unitary vectors, Eq. (4).

R = [{\hat{X}}_{Q}, {\hat{Y}}_{Q}, {\hat{Z}}_{Q}] = [\begin{matrix} u_{x} & v_{y} & w_{z} \\ u_{x} & v_{y} & w_{z} \\ u_{x} & v_{y} & w_{z} \end{matrix}]

(4)

Then, the pose of the triangle is described by the homogeneous transformation matrix 4 × 4 T^Q_C defined in Eq. (5):

T_{C}^{Q} = [\begin{matrix} R & Q \\ 0 & 1 \end{matrix}]

(5)

The rest of the analysis is performed with the D-H method for a leg. In Fig. 3 the joint variables of the mechanism leg, q_i,1 _Qi,2, >Qi,₃ are defined, as well as the lengths of the coxa, femur and tibia, d₁, d₂, d₃ respectively.

Figure 3.

D-H parameters of i-leg.

The limbs described above are similar to 3 DOF serial manipulators. From this, the homogeneous transformation matrices between frames are obtained, taking into account the D-H parameters shown in Table 1, which are: link length (a_i), link twist (α_i), joint distance (d_i) and joint angle (θ_i), required to completely describe the three joints of the legs.

Table 1.

D-H parameters for legs.

Joint	θ_i	α_i.	α_i	d_i
Coxa	q_i,1	90	d_i	0
Femur	q_i,2	0	d₂	0
Tibia	q_i,3	0	d₃,	h

The homogeneous transformation matrix that describes the translation and rotation between frames i and i-1 [9], is shown in Eq. (6).

T_{i}^{i - 1} = [\begin{matrix} \cos θ i & - \sin θ i \cos α i & \sin θ i \sin α i & a i \cos θ i \\ \sin θ i & \cos θ i \cos α i & - \cos θ i \sin α i & a i \sin θ i \\ 0 & \sin α i & \cos α i & d i \\ 0 & 0 & 0 & 1 \end{matrix}]

(6)

The homogeneous transformation matrices of the links are presented below in Eqs. (7), (8) and (9).

T_{F}^{C} = [\begin{matrix} \cos q_{i, 1} & 0 & \sin q_{i, 1} & d_{1} \cos q_{i, 1} \\ \sin q_{i, 1} & 0 & - \cos q_{i, 1} & d_{1} \sin q_{i, 1} \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}]

(7)

T_{T}^{F} = [\begin{matrix} \cos q_{i, 2} & - \sin q_{i, 2} & 0 & d_{2} \cos q_{i, 2} \\ \sin q_{i, 2} & \cos q_{i, 2} & 0 & d_{2} \sin q_{i, 2} \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}]

(8)

T_{S}^{T} = [\begin{matrix} \cos q_{i, 3} & - \sin q_{i, 3} & 0 & d_{3} \cos q_{i, 3} \\ \sin q_{i, 3} & \cos q_{i, 3} & 0 & d_{3} \sin q_{i, 3} \\ 0 & 0 & 1 & h \\ 0 & 0 & 0 & 1 \end{matrix}]

(9)

Hence, the pose of the coxa frame {C_i} with respect to {S_i}, is defined in Eq. (10).

T_{S_{i}}^{C_{i}} = T_{F_{i}}^{C_{i}} T_{T_{i}}^{F_{i}} T_{S_{i}}^{T_{i}}

(10)

Then, the frame attach to thorax {Q} is defined, with respect to {S}, in Eq. (11).

T_{c}^{s} = T_{S}^{T} T_{T}^{F} T_{F}^{C} T_{C}^{Q} .

(11)

T_{Q}^{S} = [\begin{matrix} n_{q x} & o_{q x} & a_{q x} & p_{q x} \\ n_{q y} & o_{q y} & a_{q y} & p_{q y} \\ n_{q z} & o_{q z} & a_{q z} & p_{q z} \\ 0 & 0 & 0 & 1 \end{matrix}] .

where,

\begin{array}{l} P_{q x} = r_{x} - h (c_{i, 1} v_{x} - s_{i, 1} u_{x}) + c_{i, 1} d_{1} u_{x} + s_{i, 1} d_{1} v_{x} + \\ + S_{2} d_{2} w_{x} + C_{2} d_{2} (c_{i, 1} u_{x} + s_{i, 1} v_{x}) + C_{3} d_{3} (S_{2} w_{x} + \\ + C_{2} (c_{i, 1} u_{x} + s_{i, 1} v_{x})) + S_{3} d_{3} (C_{2} w_{x} - S_{2} (c_{i, 1} u_{x} + s_{i, 1} v_{x})) \\ P_{q y} = r_{y} - h (c_{i, 1} v_{y} - s_{i, 1} u_{y}) + c_{i, 1} d_{1} u_{y} + s_{i, 1} d_{1} v_{y} + s_{i, 2} d_{2} w_{y} + \\ + c_{i, 2} d_{2} (c_{i, 1} u_{y} +_{i, 1} v_{y}) + C_{3} d_{3} (s_{i, 2} w_{y} + c_{i, 2} (c_{i, 1} u_{y} + s_{i, 1} v_{y})) + \\ + S_{3} d_{3} (c_{i, 2} w_{y} - s_{i, 2} (c_{i, 1} u_{y} + s_{i, 1} v_{y}) \\ P_{q z} = r_{z} - h (c_{i, 1} v_{z} - s_{i, 1} u_{z}) + c_{i, 1} d_{1} u_{z} + s_{i, 1} d_{1} v_{z} + s_{i, 2} d_{2} w_{z} + \\ + c_{i, 2} d_{2} (c_{i, 1} u_{z} + s_{i, 1} v_{z}) + c_{i, 3} d_{3} (s_{i, 2} w_{z} + c_{i, 2} (c_{i, 1} u_{z} + s_{i, 1} v_{z})) \\ + s_{i, 3} d_{3} (c_{i, 2} w_{z} - s_{i, 2} (c_{i, 1} u_{z} + s_{i, 1} v_{z})) \end{array}

where, [r_x,r_v,r_z] corresponds to centroid coordinates, 8_j= sin(q_i,j) and, _Cj=cos(q_i,j).

3.2 Inverse analysis position

This analysis consists of determining the joints angles, q_i,1, q_i,2 and q_i,3, starting from a given pose of the thorax.

The problem is solved for each leg separately, in a geometric way considering three restrictions: 1) all joints permit only rotation about an axis; 2) the femur and tibia always rotate around parallel axes, and 3) the physical constraints of each joint, which give a specific angular range for each active joint.

The coxa angle, Fig. 4, can be found by projecting the i-leg on the plane XY, Eq. (12).

Figure 4.

Projection of the i-leg on plane XY.

q_{i, 1} = \arctan (\frac{P_{S x}}{P_{S y}}) .

(12)

From Fig. 5 and Eq. (12), q_i,2 and q_i,3 are expressed as:

q_{i, 1} = \tan^{- 1} (c, \pm \sqrt{a^{2} + b^{2} - c^{2}}) - \tan^{- 1} (a, b) .

(13)

Figure 5.

Configuration of i-leg.

Where,

\begin{matrix} a = 2 d_{2} (\sqrt{P_{s c x}^{2} + P_{s c y}^{2}} - d_{1}), \\ b = 2 P_{s c z} d_{2}, \\ c = {(\sqrt{P_{s c x}^{2} + P_{s c y}^{2}} - d_{1})}^{2} + P_{s c z}^{2} + d_{2}^{2} - d_{3}^{2}, \\ q_{i, 3} = \cos^{- 1} [\frac{{(\sqrt{P_{s c x}^{2} + P_{s c y}^{2}} - d_{1})}^{2} + P_{s c z}^{2} - d_{2}^{2} - d_{3}^{2}}{2 d_{2} d_{3}}] . \end{matrix}

(14)

4. Infinitesimal kinematics

The velocity and acceleration analysis of the robot are solved using screw theory. As an approach to the study of the dynamics of the robot, the manipulator is assumed as a parallel mechanism. Its limbs are modelled as shown in Fig. 6.

Figure 6.

i-th kinematic chain robot.

Note that the active revolute joints are associated with the screws ³$⁴_i, ⁴$⁵_i and ⁵$⁶_i. During normal operation of the robot, the motor corresponding to twist ³$⁶_i, is used to move its limb during the locomotion. When the manipulator is in a stable position, this joint becomes a passive joint; however, since the motor has an encoder, its angular position, velocity and acceleration are known. The motors associated with screws ³$⁴_i and ⁴$⁵_i, orientate the thorax with respect to the global frame.

On the other hand, the suction cup is modelled as a UP linkage, and its corresponding screws are ¹$¹_i, ¹$²_i and ²$³_i whose axes are mutually perpendicular. In practice, the displacements of these screws rely on the elastic properties of the suction cup.

4.1 Preliminary concepts of screw theory

An infinitesimal screw $, or twist, is a vector in $ ∈ ℝ⁶ that consists of a primal part, P ($)=ŝ, and a dual part, D ($)=80, where ŝ is a unitary vector along the instantaneous axis of the screw; so is the moment produced by vector ŝ with respect to a vector τ_O/P:

s_{O} = h \hat{s} + \hat{s} \times r_{O ∕ P},

where h is the screw pitch.

Moreover, the following operations are defined [8]. Let $₁=(ŝ₁,s_O1),$₂=(ŝ₂,s_O2),$₃=(ŝ₃,s_O3) be three screws, and λ ∈ ℝ.

Addition:

$_{1} + $_{2} = ({\hat{s}}_{1} + {\hat{s}}_{2}, s_{O 1} + s_{O 2})

(15)

Multiplying by scalar:

λ $_{1} = (λ {\hat{s}}_{1}, λ s_{O 1})

(16)

Lie product

[\begin{matrix} $_{1} $_{2} \end{matrix}] = [\begin{matrix} {\hat{s}}_{1} \times {\hat{s}}_{2} \\ {\hat{s}}_{1} \times s_{O 2} - {\hat{s}}_{2} \times s_{O 1} \end{matrix}]

(17)

Killing form

K i ($_{1}, $_{2}) = {\hat{s}}_{1} \cdot {\hat{s}}_{2} .

(18)

Klein form

{$_{1}; $_{2}} = {\hat{s}}_{1} \cdot s_{O 2} + {\hat{s}}_{2} \cdot s_{O 1} .

(19)

4.2 Velocity analysis

The forward velocity analysis consists of finding the velocity state of the thorax, Vo, with respect to the global reference frame, given the joint rate velocities aω_bⁱ. It is given by the following expression:

V_{O} = [\begin{matrix} ω \\ v_{O} \end{matrix}] =_{0} ω_{1}^{i}^{0} $_{i}^{1} +_{1} ω_{2}^{i}^{1} $_{i}^{2} + \dots +_{5} ω_{6}^{i}^{5} $_{i}^{6},

(1)

where ω and v_o are, respectively, the angular and linear velocity of the thorax with respect to the global frame. The terms aω_bⁱ represent the ratios of rotational or translational velocity of body b with respect to a belonging to i-th kinematic chain. Moreover, given the actuation scheme of the robot ₅ωⁱ₆=q̇_i,1, ₄ωⁱ₅=q̇_i,2 and ₃ωⁱ₄=q̇_i,3, Eq. (20) can be written in matrix form as follows.

V_{O} = J_{i} Ω_{i},

(21)

where $J_{i} = [\begin{matrix} ^{0} $ _{i}^{1} \end{matrix}, ^{1} $ _{i}^{2}, ^{2} $ _{i}^{3}, ^{3} $ _{i}^{4}, ^{4} $ _{i}^{5}, ^{5} $ _{i}^{6}]$ is the Jacobian of the $i - t h$ kinematic chain and $Ω_{i} = {[^{0} ω _{1}^{i}, ^{1} ω _{2}^{i}, ^{2} ω _{3}^{i}, ^{3} ω _{4}^{i}, ^{5} ω _{6}^{i}]}^{T}$ is the vector of joint velocities of such a chain.

A useful concept to simplify this analysis is reciprocity. Two screws, $₁ and $₂ are reciprocal if the Klein form between them is zero, {$₁; $₂} = 0. With this definition, one can demonstrate that the screw ¹$²_i is reciprocal to all screws belonging to the i-th chain except to the ⁵$⁶_i. Then, applying the Klein form on both sides of Eq. (20), it gives:

{^{1} $_{i}^{2}; V_{O}} = {\dot{q}}_{i, 1} {^{1} $_{i}^{2};^{5} $_{i}^{6}}

(22)

On the other hand, let $ ^τ _i be a zero-pitch screw that passes by point S _i and its primal part is $P ($_{i}^{r}) = D (^{2} $ _{i}^{3}) \times P (^{1} $ _{i}^{2})$ . Then, applying the Klein form between both sides of Eq. (20) and this screw,

\begin{array}{l} {$_{i}^{r}; V_{O}} = {\dot{q}}_{i, 3} {$_{i}^{r};^{3} $_{i}^{4}} + {\dot{q}}_{i, 2} {$_{i}^{r};^{4} $_{i}^{5}} + \\ + {\dot{q}}_{i, 1} {$_{i}^{r};^{5} $_{i}^{6}} \end{array}

(23)

Eqs. (22) and (23) can be written as a lineal system as follows:

M^{T} Δ {\vec{V}}_{O} = \dot{Q},

(24)

where $M = [\begin{matrix} ^{1} $ _{1}^{2} \end{matrix}, ^{1} $ _{2}^{2}, ^{1} $ _{3}^{2}, $_{1}^{r}, $_{2}^{r}, $_{3}^{r}]$ , $Δ = [\begin{matrix} 0 & Ι \\ Ι & 0 \end{matrix}]$ is a polarity operator defined by an identity matrix $Ι_{3 \times 3}$ , and the matrix 0 of the same order. Finally,

\dot{Q} = [\begin{matrix} {\dot{q}}_{1, 1} {^{1} $_{1}^{2};^{5} $_{1}^{6}} \\ {\dot{q}}_{2, 1} {^{1} $_{2}^{2};^{5} $_{2}^{6}} \\ {\dot{q}}_{3, 1} {^{1} $_{3}^{2};^{5} $_{3}^{6}} \\ {\dot{q}}_{1, 3} {$_{1}^{r};^{3} $_{i}^{4}} + {\dot{q}}_{1, 2} {$_{1}^{r};^{4} $_{1}^{5}} + {\dot{q}}_{1, 1} {$_{1}^{r};^{5} $_{1}^{6}} \\ {\dot{q}}_{2, 3} {$_{2}^{r};^{3} $_{i}^{4}} + {\dot{q}}_{1, 2} {$_{1}^{r};^{4} $_{1}^{5}} + {\dot{q}}_{1, 1} {$_{2}^{r};^{5} $_{1}^{6}} \\ {\dot{q}}_{3, 3} {$_{3}^{r};^{3} $_{i}^{4}} + {\dot{q}}_{3, 2} {$_{3}^{r};^{4} $_{3}^{5}} + {\dot{q}}_{3, 1} {$_{3}^{r};^{5} $_{3}^{6}} \end{matrix}] .

Thus the state of velocity Vo can be calculated without knowing the ratio of passive velocities, which represents an advantage over other models.

Moreover, the inverse velocity analysis consists of finding the ratios of manipulator velocities for a given velocity state. To achieve this, it is necessary to isolate Ω_i in Eq. (21).

4.3 Acceleration analysis

Brand in 1947 introduced the concept of reduced acceleration state, A₀, of a rigid body [6]

A_{O} = [\begin{matrix} α \\ a_{O} - ω \times v_{O} \end{matrix}],

(25)

where $α = \dot{ω}$ and $a_{O}$ are the angular and linear acceleration of the thorax, respectively. Furthermore, those accelerations can be written in the form of screws as follows:

A_{O} = J_{i} {\dot{Ω}}_{i} + L_{i},

(26)

where ${\dot{Ω}}_{i}$ is the matrix whose elements are the joint acceleration ratios of the i - th. kinematic chain. Moreover L_i is the Lie screw products, which are given by:

L_{i} = \sum_{j = 0}^{4} [​^{j} w ​_{j}^{i} + ​_{1}^{j} $ ​_{i}^{j + 1} \sum_{k = j + 1}^{5} ​^{j} w ​_{k}^{i} + ​_{1}^{j} $ ​_{i}^{k + 1}] .

In the same manner as the velocity analysis, the Klein form is applied to Eq. (26) and screws ⁰$¹_i and ²$³_i. This leads to

M^{T} Δ {\vec{A}}_{O} = \overset{..}{Q}

(27)

where,

\begin{matrix} \overset{..}{Q} = {[\begin{matrix} a_{1}, a_{2}, a_{3}, b_{1}, b_{2}, b_{3} \end{matrix}]}^{T} \\ a_{i} = {\dot{q}}_{i, 1} {^{1} $_{i}^{2};^{5} $_{i}^{6}} + {^{1} $_{i}^{2}; L_{i}} \\ b_{i} = {\dot{q}}_{i, 3} {$_{i}^{r};^{3} $_{i}^{4}} + {\dot{q}}_{i, 2} {$_{i}^{r};^{4} $_{i}^{5}} + {\dot{q}}_{i, 1} {$_{i}^{r};^{5} $_{i}^{6}} + \\ {$_{i}^{r}; L_{i}} \end{matrix}

Finally, to solve the inverse analysis, it is necessary to isolated in Eq. (26).

5. Numerical example

In this section a numerical example is provided with the purpose of showing how to deal with the kinematics of this model. Hereafter the units used are millimetres, unless otherwise stated. In this case, the dimensions of the geometrical parameters of the moving triangular prism are a = 153 and b=131.203, meanwhile Q = (116.803, 216.711,160.772) is the centre of the prism; whereas the nominal positions of the revolute joints connecting the limbs to the first link are:

\begin{matrix} P_{1} = (163.838, 218.372, 192.622), \\ P_{2} = (116.497, 218.769, 70.298), \\ P_{3} = (30.387, 218.769, 196.766), \end{matrix}

The coordinates of the suction cups are:

\begin{matrix} S_{1} = (124.61, 30.166, - 39.261), \\ S_{2} = (- 68.822, 30.166, 243.952), \\ S_{3} = (253.766, 30.166, 253.893) . \end{matrix}

Links length, coxa, femur and tibia are given by: d₁ = 57.61, d₂ = 64.54 and d₃ = 11.41, respectively.

Fig. 7 shows the joints angle of the legs L₁, L₂, and L₃, for the time range 0-0.25 seconds.

Figure 7.

Position of the joints angle.

Firstly, the direct position equations were programmed to calculate the thorax pose, taking into account the initial and final positions, as well as geometric parameters. Fig. 8 shows the robot moving from the initial to the final position.

Figure 8.

Hex-piderix position.

Fig. 9 shows the thorax position, where the components x and y are kept constant during the motion of the robot, without considering the effect of the suction cup deformation.

Figure 9.

Thorax position.

With the above data, both angular and linear velocities of the numerical example are calculated. The resulting time history of the linear velocity components of the centre of thorax using the screw theory is shown in Fig. 10. Angular velocity is shown in Fig. 11.

Figure 10.

Linear velocity of the Hex-piderix robot.

Figure 11.

Angular velocity of the Hex-piderix robot.

As is shown in Figs. 8 and 9, the components of the position vector of thorax are not constants due to suction cup deformations; hence the angular velocity is different to zero.

6. Conclusion

It was developed a method to determine the rotation matrix of the thorax of a legged robot with respect to the global frame, without having to use trigonometric functions, since these could lead to ambiguous values for the orientation. Although, this method is usually applied in the study of parallel robots, it can be successfully employed in the analysis of legged robots.

The velocity and acceleration states of the thorax with respect to the fixed platform are written in screw form through each one of the three limbs of the robot. The equations obtained are linear, simple and compact through the use of the Klein form of Lie algebra. In addition, it is not necessary to know the velocities and accelerations of passive joints to solve forward velocity and acceleration analysis, respectively.

This kinematics study provides the basis for generating control strategies that allow the robot to recover from perturbations during locomotion, as well as keeping the thorax in a constant orientation.

It would be helpful to consider the suction cup attached to each tip of the tibia when modelling the kinematics of a legged robot, as this can introduce an error.

Footnotes

7. Acknowledgments

The authors acknowledge to the SIP-IPN for funding the research project 20121376

References

de Santos

P. Gonzalez

Garcia

and Estremera

Quadrupedal locomotion: an Introduction to the control of four-legged robots. Verlag, London: Springer, 2006, pp. 3–32.

Sintov

Avramovich

and Shapiro

“Design and motion planning of an autonomous climbing robot with claws.”

Robotics and Autonomous System, vol.59, pp.1008–1019, June 2011.

Roy

S. Shekhar

Singh

A. Kumar

and Pratihar

D. Kumar

. “Estimation of optimal feet forces and joint torques for on-line control of six-legged robot.” Robotics and Computer-Integrated Manufacturing, vol.27, pp.910–917, March 2011.

Garcia-Lopez

M.C.

Gorrostieta-Hurtado

Vargas-Soto

Ramos-Arreguin

J.M.

Sotomayor-Olmedo

and Moya-Morales

J.C.

“Kinematic analysis for trajectory generation in one leg of a hexapod robot.”

Procedia Technology, vol.3, pp.342–350, 2012.

Segundo-Potts

and da Cruz

J.J.

(2011). A Kinematical and Dynamical Analysis of a Quadruped Robot, Mobile Robots, Available: http://www.intechopen.com/books/mobile-robots-current-trends/a-kinematical-and-dynamical-analysis-of-a-quadruped-robot.

Fang

Zhang

, “The Mathematical Model and Control Scheme of a Four-Legged Robot Based on GZ-I and Note Module”, Lecture Notes on Artificial Intelligence, vol. 6424, pp. 300–309, 2010.

Cheng

and Chen

S.Y.

, “Exact solution of impulse response to a class of fractional oscillators and its stability”, Mathematical Problems in Engineering, Vol. 2011, AID 657839, 2011

Gallardo-Alvarado

Aguilar-Najera

C.R.

Casique-Rosas

Perez-Gonzalez

Rico-Martinez