TMCnet News

Kinematics Analysis and Simulation on Transfer Robot with Six Degrees of Freedom [Sensors & Transducers (Canada)]
[September 23, 2014]

Kinematics Analysis and Simulation on Transfer Robot with Six Degrees of Freedom [Sensors & Transducers (Canada)]


(Sensors & Transducers (Canada) Via Acquire Media NewsEdge) Abstract: Study focuses on transfer robot with Six Degrees of Freedom, establishing kinematic equation by D-H method, analyzing forward kinematics and obtaining inverse kinematics by using method of inverse transform. Based on vector product, it develops velocity Jacobian matrix of robot. The geometric model of robot virtual prototype is established by SolidWorks software and generates parameters such as mass and moment. Kinematic simulation for robot is performed by Mathematica software and develops curve graph of displacement, velocity and accelerated speed in x, y and z direction in end executor center of robot with measurement, analysis and assessment, which provides foundation for further kinematics analysis and structure optimization as well as motion control of robot. Copyright © 2014 IFSA Publishing, S. L.



Keywords: Transfer robot, Kinematics, Simulation, Optimization, SolidWorks, Mathematica.

(ProQuest: ... denotes formulae omitted.) 1. Preface In recent years, with deep development of technical research on robot and wider study fields, kinematic simulation analysis of robot plays more important role in spatial planning, trajectory control, motion control, off-line programming, optimization design and other aspects. Robot kinematic studies kinematic relation between each joint and rigid body, position of parts involved and relation of velocity and accelerating speed with time. Robot kinematics includes: one is direct kinematics where joint variable is given and calculate hand appearance, the other is inverse kinematics where hand appearance is given and calculate joint variable.


Transfer robot is an advanced automatic device which could increase productivity and improve working condition.

2. Robot Kinematics Model Articulated robot with 6-DOF has very complicated mechanical structure, consisting of base, waist joints, upper arm, forearm, wrist joint and hands, every motion unit of which is composed of a small mechanical system, such as shaft, bearing, bushing, key, gear, motor as well as reducer. As shown in Fig. 1.

Motion arm of 6-DOF consists of 6 arms and 6joints. Each arm is described with 4 parameters of ai_,, at_,, di, 0,., where ai_, and at_, means feature of arm i-1; di and 0; means relation of arm i-1 and arm i. For rotary joint, joint angle is joint variable and other parameters are constant; for prismatic joint, polarization is joint variable and other parameters are constant. This method to describe mechanism kinematics was brought up by Denavit and Hartenberg in 1955, called D-H method [1,2].

Establish link coordinate system by D-H method, as shown in Fig. 2. Kinematics and structure parameters of related arms are shown in Table 1.

Where ai means distance from z. to z.+I along xi ; a. means angle from z. to zi+] around xi ; di means distance from xi_, to xi along z. ; 0,. means angle from x._, to xt around z. [3].

2.1. Direct Kinematics Refer to link parameter to make sure arm transformation matrix expression In accordance with the "left to right" principle: ...(1) Put parameters in Fig. 1 into transmission matrix (1), to calculate each arm homogeneous transformation matrix °XT ,\T ,\T ,\T ,45T ,s6T . Multiply above matrixes to get robot end component homogeneous transformation matrix: ...(2) Put parameters in Fig. 1 into equation (1) and (2) to get position matrix °6T of robot end relative to basic coordinate system: ...(3) Equation (3) means position relationship of end arm coordinate system {6} relative to basic coordinate system {0}, which is the foundation of each joint kinematics analysis of robots.

Initial position: ...

Put initial value of 6t into equation (3) to get: ...

It accords with the position in Fig. 2, which proves above calculation process is correct.

By above analysis, establish kinematics equation of relation between transfer robot end actuator position and other joint variable in rectangular coordinate system, which gets forward solution of robot kinematics equation [4].

2.2. Reverse Kinematics Reverse kinematics is to get joint variable by hand position, which is used to make sure motor parameters to drive joint and plan trajectory [5].

Common reverse kinematics methods include counter transformation method, geometric method and Pieper method and so on. The paper applies for counter transformation method to solve robot reverse kinematics, that is, to multiply a reverse matrix on the left of 6 joint transformation matrixes and calculate to find the invariable element in the right, and then make these variables equal with those in the left to get a trigonometric function. According to position of vectorial combined with robot arm parameter d to get joint variables 6X, #2,..., 06.

Inverse transformation °17'"1(0I) left multiply function (2) on left and right, ... (4) to get , #3, there are two possibilities for each variable.

Inverse transformation \T~X left multiply function (2) on left and right, ...(5) Given function is equal on right and left sides to get 6^, #4, there are two possibilities for each variable.

Inverse transformation °4T~X left multiply function (2) on left and right, ... (6) Given function is equal on right and left sides to get 05, which is closed form. In the same way to get closed form of 06.

Above analysis shows inverse of kinematics function has uncertainty and has many groups of solution which are possible different space position of kinematic joint.

2.3. Calculate Jacobian Matrix Jacobian matrix of motion arm is its linear transformation of operation speed and joint speed, which can be seen as transmission ratio of movement from joint space to operation space [6], that is: ...(7) where X means operation velocity vector, q means joint velocity vector. The number of movement dimensions of robot in operation space is equal to that of rows of Jacobian matrix and joints are equal to column. Above analysis shows Jacobian matrix of robot with six degrees of freedom is a 6*6 matrix. The paper applies for vectorial method to calculate robot Jacobian matrix.

The line i of Jacobian matrix is: ...(8) where lP°n means robot end origin of coordinate in basic coordinate system {o} relative to position vectorial in coordinate system {i}, Z( means unit vectorial of z axis of coordinate system {i} in basic coordinate system {o}. °R means direction cosine matrix. From above transformation matrix, get °iR,Zi,'P5 (i = 1,2,...,6), and put them into function (8) to get line i (i = 1,2,...,6) of Jacobian matrix, so as to make robot Jacobian matrix..

Jacobian matrix of robot is shown as below: ...(9) 3. Kinematics Simulation Put model into Mathematica, which includes 7 parts and each means base, vertical shaft, base turntable, upper arm, forearm, paws and wrist. Base is fixing part and rotating joints are set between other parts around Z(. axis. Given variable to each joint, simulate, analyze and compare kinematics.

3.1. Analysis on Simulation Result Analyze point in the middle of paws, which is midpoint of paw and coincides with original point of coordinate system. Output data is shown as Fig. 3.

Fig. 3 shows that when time t=0, operation arm is in initial position, displacement in x, y, z direction is -500, 610, 1300 (unit: mm), which is the same as calculation result.

Fig. 4 shows that when time t=0, velocity in x, y, z direction is -297, 0, 245, and maximum velocity is 297, 200, 245 (unit: mm/s) which is almost the same as calculation result.

Fig. 5 shows that when time t=0, accelerated speed in x, y, z direction is 20, -113, 0 and maximum accelerated speed is 95, 113, 46 (unit: mm/s2), which is the same as theoretical analysis.

The same verification is got from Fig. 6 and Fig. 7.

Change length of each arm and movement parameter of each joint to get different displacement, velocity and accelerated speed of arm in different conditions.

4. Conclusions 1) Based on robot D-H matrix theory, establish robot kinematics mathematic model, deduce homogeneous transformation matrix of transfer robot kinematics function by combing with robot structure, by which angular of every joint and paw position can be calculated mutually and get direct and inverse kinematics.

2) Vector product method is applied to deduce speed Jacobian matrix, which could calculate displacement and speed of points in paw with each joint angular and speed mutually.

3) Simulate by software to verify given structure parameters reasonable.

4) Mathematic simulation software is applied to simulate kinematics of transfer robot with 6-degrees of freedom, measure velocity, accelerated speed, displacement curve of robot end, analyze dynamic movement features of robot. Simulation result shows: robot moves steadily in working condition which can meet working requirements.

References [1], Zixing Cai, Robotics, Tsinghua University Press, Beijing, 2000.

[2], A. K. Basu, Dayam Naveen Kumar, Nimishakavi Karthik, Vikas Padwekar, Modeling and analysis of a micro robot, Journal of Mechatronics, Vol. 1, Issue 2,2013, pp. 118-123.

[3] - You Lun Xiong et. al, Robotics, Mechanical Industry Publishing House, Beijing, 1993.

[4] . http://bbs.lcrobot.com/thread-3258-l-l.html [5] . K. G. Shin, N. D. McKay, Minimum-time control of robotic manipulators with geometricpath constraints, IEEE Transactions on Automatic Control, Vol. 30, Issue 6,1985, pp. 531-541 [6]. J. E. Bobrow, S. Dubowsky, J. S. Gibson, Timeoptimal control of robotic manipulators along specified paths, Journal of Robotic Research, Vol. 4, Issue 3,1985, pp. 554-561.

1 Yi Lu,2 Dan Liu 1 Jiangsu Jianzhu Institute, Xuzhou, Jiangsu, 221116, China 2Xuzhou Coal Mining Group, 221000, China 1 Tel.: 051683889081, fax: 051683889081 1 E-mail: [email protected] Received: 20 May 2014 /Accepted: 31 July 2014 /Published: 31 August 2014 (c) 2014 IFSA Publishing, S.L.

[ Back To TMCnet.com's Homepage ]