Next Article in Journal
A Multi-Objective Reinforcement Learning Based Controller for Autonomous Navigation in Challenging Environments
Previous Article in Journal
Optimization Procedure and Toolchain for Roll Dynamic Geometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Modeling Method of Multibody System of 6-DOF Robot Based on Screw Theory

1
Robotics Institute, Beihang University, Beijing 100191, China
2
Southwest China Research Institute of Electronic Equipment, Chengdu 610036, China
3
AVIC Manufacturing Technology Institute, Beijing 100024, China
*
Author to whom correspondence should be addressed.
Submission received: 10 June 2022 / Revised: 18 June 2022 / Accepted: 21 June 2022 / Published: 22 June 2022
(This article belongs to the Section Machine Design and Theory)

Abstract

:
An accurate dynamic model is a prerequisite for realizing precise control of industrial robots. The dynamics research of multi-degree of freedom (DOF) robots is relatively unexplored and needs to be solved urgently. In this paper, a dynamic modeling method of multibody system of 6-DOF robot is proposed based on the screw theory. The established dynamic model has a more concise and unified mathematical form, and the modular matrix expression is convenient for the control of the robot. In order to ensure that the screw method is suitable for motion in a wide range of angles, quaternions are used as generalized angular coordinates, and the model established thereby eliminates singularities and improves computational efficiency. The correctness and accuracy of the screw method is verified by the simulation example, and the modeling theory and method can provide a theoretical basis for the precise control of the robot.

1. Introduction

The shortcoming of the low absolute positioning accuracy of the serial robot makes it incapable of precision machining. The fundamental way to solve this problem is to realize the precise control of the robot based on dynamics, that is to achieve precise force control and position control at the same time so as to meet the requirements of intelligence and precision. Therefore, accurate dynamic model is of great significance to the application of serial robots in high-end manufacturing.
The research of the robot dynamics based on traditional methods has been widely carried out. The traditional dynamic methods are Newton–Euler method and Lagrange method. The Newton–Euler method is a vector mechanics method and a recursive method. The motion parameters of the latter rigid body are obtained recursively from the former rigid body. The increase in the number of rigid bodies will make the derivation process cumbersome and complicated, and the parameters coupling makes the dynamic problem a nonlinear problem that is difficult to solve, so it is suitable for problems with few degrees of freedom (DOFs), but not for problems with multiple DOFs. The Lagrange method is an analytical mechanics method. The Lagrange quantity in the form of energy will become cumbersome and complicated with the increase of the number of rigid bodies, and the phenomenon of parameters coupling will intensify, which significantly increases the number of dynamic equations and the difficulty of solving them. Therefore, it is suitable for problems with few DOFs, but not for problems with multiple DOFs.
Modern mathematical method based on screw theory, Lie groups and Lie algebras, have been used in the research of robot dynamics. It is a vector mechanics method, and the dynamic equations expressed have the advantages of simplicity and unity. The existing screw method, like Newton–Euler method, is also a recursive method. Although it is possible to describe and solve dynamic problems more concisely and uniformly, when dealing with multi-DOF problems, parameters coupling also makes the calculation complicated and difficult to solve.
Dynamics of multibody systems is a branch of classical mechanics. With the development of computing technology, it has gradually become an independent discipline which studies systems composed of a large number of rigid bodies. Dynamics of multibody systems describe the interconnections between rigid bodies by graph theory. Firstly, the dynamic model of a single rigid body is established by vector mechanics, analytical mechanics, or Kane method. Then, the dynamic equations of each rigid body are combined to obtain the model of the system by graph theory, and the constraint equations are established to eliminate the redundant variables in the model and make the free variables equal to the number of DOF. The description of graph theory makes the dynamic model of the robot obtained by combining the dynamic models of each rigid body. Since the dynamic models of each rigid body are independent of each other, the coupling of the dynamic and kinematic parameters between the rigid bodies is eliminated. Decoupling is a necessary condition for solving a multi-DOF dynamic model, and the resulting dynamic model has two advantages. (1) Decoupling enables the multi-DOF dynamic model to have an analytical solution, which is convenient for the solution and control of the model. (2) Decoupling enables the dynamic equations of each rigid body to be combined and expressed as a modular matrix equation. The matrix equation, as the dynamic model of the robot, has a simple and compact form and is easy to control. It should be noted that when traditional dynamic methods and existing screw methods deal with multi-DOF dynamic problems, due to the difficulty of decoupling dynamic and kinematic parameters, the dynamic equations of each rigid body can only exist independently and cannot be unified into a modular matrix equation, which brings difficulties to the solution and control. Therefore, dynamics of multibody systems is a powerful tool for studying multi-DOF mechanisms and is suitable for analyzing complex systems composed of a large number of components in the fields of robots, vehicles, and aerospace.
Aiming at the problem of accurate dynamic modeling of multi-DOF serial robots, 6-DOF robot is selected as the research object in this paper. A dynamic modeling method of multibody system based on screw theory is proposed for the first time. This method has the advantages of simple and compact expression of screw theory and being good at solving multi-DOF dynamic problems of dynamics of multibody systems. Compared with the previous screw method, this screw method has the following characteristics. (1) The kinematic and dynamic parameters can be decoupled, and the dynamic model can be expressed in the form of a modular matrix, which is convenient for control. This method can express the dynamic model of the 6-DOF robot with only one matrix equation, so it has a more concise and unified mathematical expression. (2) By establishing and introducing constraint equations, the ideal constraining force and moment can be evaluated and analyzed. (3) Quaternions are used as generalized coordinates in this method, which eliminates singularities, improves computational efficiency, and makes the robot suitable for large-angle range motion. It should be noted that the traditional dynamics method and the existing screw method are suitable for solving problems with few DOFs. While the dynamics of multibody systems are suitable for solving problems with multiple DOFs, the applications of the two are different. The number of dynamic equations of the problem with few DOFs is few, the expression is relatively simple, and the calculation difficulty is relatively small. The calculation and solution of the problem do not have to rely on computers. The multiple-DOF problem will inevitably lead to more dynamic equations and more complex expressions. Therefore, the calculation is very difficult, and the calculation and solution must be completed with the help of computers. In fact, it is the development of computer and computing technology that promotes the development of dynamics of multibody systems. Therefore, in terms of computational efficiency, it is obvious that the problem with few DOFs is higher than the problem with multiple DOFs. However, for the multi-DOF problem, the computational efficiency can be improved by selecting the quaternion as the generalized angular coordinate. This is because when calculating the coordinate transformation matrix, the sine and cosine functions of the angle need to be calculated 29 times when the Euler angle or the Cardano angle is used as the angle coordinate, and the sine and cosine functions do not need to be calculated when the quaternion is used as the angle coordinate. The coordinate transformation matrix needs to be calculated many times in the dynamic solution process, so the quaternion can greatly reduce the calculation steps and improve the calculation efficiency. The dynamic model established by this method is correct and accurate, which is verified by a numerical simulation example. The dynamic theory and method in this paper can provide a theoretical basis for precise control and structural reliability analysis.

2. Related Work

2.1. Traditional Dynamics Method

The dynamics research based on the traditional Newton–Euler method and Lagrange method has been widely carried out and achieved corresponding results. Khalil et al. [1] studied the forward and inverse dynamics of the Gough–Stewart robot. The dynamic model of the robot is established according to the dynamic model of the legs. The dynamic model established by the Newton–Euler method has an intuitive physical explanation. In order to reduce the computational cost, the basic inertial parameters of the robot are determined through analysis. Djuric et al. [2] studied the dynamics of a reconfigurable Puma-Fanuc robot, established a dynamics model with the Newton–Euler recursion method, and each element of the inertial matrix, Coriolis moment matrix, centrifugal moment matrix, and gravitational moment vector was automatically generated using a newly developed automatic separation method. Ibrahim et al. [3] studied a hybrid robot with parallel modules connected in series and established the forward and inverse dynamics models of the robot by the Newton–Euler recursion method. The proposed algorithm can be programmed with symbolic and inertial parameters. Koopaee et al. [4] developed a two-dimensional modular snake-like robot, established the dynamic model of the robot by the Euler–Lagrange equation, and proposed an adaptive controller based on torque feedback in the gait parameter space and optimized the control gain. The effectiveness of the model and controller is verified by simulation and experiment. Kalyoncu et al. [5] studied a flexible robotic manipulator with rotating prismatic joints and established a dynamic model of the manipulator based on the Lagrange equation. Under the action of external driving force and axial force, the end of the flexible manipulator moves along a path composed of multiple straight lines. Li et al. [6] studied a 2-DOF planar parallel robot with flexible links for high-speed pick-up and placement operations. The dynamic model of the robot is established by the Euler–Lagrange method, and the dynamic characteristics of the robot are analyzed based on the model, which improves the dynamic accuracy of the end effector at high speed. Traditional dynamic methods are mostly used in the research of less-DOF robots, but they are difficult to be used in the research of multi-DOF robots. The Newton–Euler method is a recursive method. The motion parameters of the latter rigid body are obtained based on the motion parameters of the former rigid body. Therefore, the more DOFs, the more lengthy the dynamic equations and the larger the number of equations, which makes the dynamic model difficult to solve. The Lagrange method is a dynamic method based on the Lagrange quantity and derived by combining the variational method with the principle of least action. The Lagrange quantity in the form of energy becomes lengthy as the number of rigid bodies increases, which makes the Lagrange equations lengthy, and the more generalized coordinates due to the increased number of rigid bodies also increase the number of Lagrange equations. The solution of the dynamic model thus becomes very difficult. The lengthy dynamic expressions and more dynamic equations make the dynamic model difficult to solve for two reasons. (1) The coupling of dynamic parameters (mass and principal moment of inertia) and kinematic parameters make it difficult to obtain analytical solutions of the dynamic equations. (2) The ideal constraining force and moment that cannot be measured for each hinge increase the unknown quantities in the equation, which brings difficulties to the solution and control of the dynamic model. The essence of ideal constraining force and moment is the constraint condition that makes the robot move according to objective laws, and the robot without constraint conditions cannot be controlled, so the ideal constraining force and moment need to be solved.

2.2. Screw Theory

The screw theory originated in the 19th century, and with the evolution of algebra to geometric algebra, it is interrelated with the newly born Lie groups and Lie algebras, which together constitute the modern mathematical foundation of robotics. The screw theory is a vector mechanics method, so the physical meaning is clear, and it studies the translation and rotation of a straight line in space. The algebraic connotation of screw theory is embodied in screw algebra, which is a vector algebra used to describe screw, and is also a subalgebra of Lie algebra. Screw algebra, Lie group, and Lie algebra are used to describe the motion of straight lines in space and their related algebraic operations. The description has geometrical intuition and algebraic abstraction. Screw is a ray passing through the origin in Lie algebra, a geometric quantity, and an element of projective Lie algebra. A twist is a screw with velocity amplitude, which describes the motion of rigid body and is an element of Lie algebra. A wrench is a screw with force amplitude and is an element of dual Lie algebra [7,8]. The screw theory unifies velocity and angular velocity into twist, and unifies force and moment into wrench, which can succinctly and uniformly describe the motion of spatial mechanisms and solve problems. At present, the modern mathematical methods of screw theory, Lie groups and Lie algebras, have been applied by some roboticists [9,10,11,12,13,14,15]. Gallardo et al. [16] proposed a method for dynamic analysis of parallel robots, which is based on screw theory and virtual work principle, and verified the generality of the method by applying it to a 2-DOF parallel spherical mechanism and a Gough–Stewart platform. Zhao et al. [17] studied the modeling method based on the screw theory, which can simplify the kinematics and dynamics of the manipulator, and proposed a dynamic coupling index based on the Hessian matrix. By applying the method to analyze the folding parallel robot, the correctness of the modeling method and coupling index is verified. Gallardo et al. [18] conducted kinematics and dynamics research on the 2(3-RPS) manipulator based on the screw theory and virtual work principle. The screw method reduces the calculation time, and the expression is simple and compact, which can be used to solve the velocity and acceleration analysis of the space mechanism. Gallardo et al. [19] studied the 4-PRURh and 4-PRURv parallel robots, established the kinematics and dynamics models of the robot with screw theory, deduced the DOF of the parallel manipulator, and obtained the input and output equations of the velocity and acceleration of the robot in a compact form. Gallardo et al. [20] proposed a kinematics and dynamics analysis method for solving modular spatially redundant manipulators. The manipulator is made up of optional three-legged parallel mechanisms in series. The kinematics and dynamics of the basic module are studied by the screw theory, and the expression of the basic module is extended to a spatially redundant manipulator. Ding et al. [21] established a dynamic model of a robot with spatially flexible links based on screw theory and proposed a Ding–Holzer method by making some reasonable assumptions about the system. The method takes into account the spatial force and deformation of the moving flexible arm and can be used to solve the dynamic modeling problem of the flexible arm with six-dimensional force and deformation at the tip. Compared with the traditional dynamic method, the existing screw method can make the expression of the dynamic model concise and compact, however, when applied to a multi-DOF robot, the coupling of dynamic and kinematic parameters and the solution of ideal constraints cannot be solved.

3. Dynamic Modeling of Multibody System Based on Screw Theory

3.1. General Equation of Dynamics

The general equation of dynamics is the equation of dynamics expressed in generalized angular coordinates, that is, the angular coordinates are not specified as Euler angles, Cardano angles, or quaternions.
The dynamic equation derived from the screw theory has a more concise and unified expression, which is convenient for the dynamic modeling and control of the robot. The dynamics of multibody system is a powerful tool for solving multi-DOF problems. It adopts the mathematical description of graph theory to decouple the dynamic and kinematic parameters of each rigid body and solves the dynamic model with the help of a computer. In this paper, a modeling method of 6-DOF robot is proposed by combining the screw theory with the dynamics method of multibody system. This method is applicable to the dynamic modeling and solution of 6-DOF robot. The dynamic model can be expressed by only one matrix equation, and the expression is more concise and unified.
The 6-DOF robot has 6 rotational DOFs and consists of 6 rigid bodies B i   ( i = 1 , 2 , , 6 ) and 6 hinges O 1 , O 2 , , O 6 , where B 6 is the end effector, as shown in Figure 1. In order to facilitate analysis and calculation, each rigid body is simplified as a homogeneous cylinder with mass m i and central inertia tensor J i , each rigid body reference system ( O 1 , e ¯ ( 1 ) ) , ( O 2 , e ¯ ( 2 ) ) , , ( O 6 , e ¯ ( 6 ) ) is established at the position of the centroid O c 1 , O c 2 , , O c 6 , and the world reference system ( O 0 , e ¯ ( 0 ) ) is established at the centroid of the bottom surface of the base B 1 . The directions of the base vectors e 1 ( i ) , e 2 ( i ) , e 3 ( i ) are shown in Figure 1b.
The multibody system dynamic equation of the robot consists of the dynamic equations of each rigid body. Therefore, the dynamic equation of the rigid body should be established first. The dynamic equation of B i is derived according to the screw theory (the underlined ¯ represents the matrix form).
ψ ¯ i v ^ ˙ ¯ i + v ^ ˜ ¯ i ψ ¯ i v ^ ¯ i = F ^ ¯ i ( i = 1 , 2 , , 6 )
Here, ψ ¯ i , v ^ ¯ i , v ^ ˙ ¯ i , and F ^ ¯ i are the generalized inertia matrix, twist, acceleration screw, and wrench of B i   respectively, and v ^ ˜ ¯ i is a 6-order square matrix.
ψ ¯ i = m i E _ 3 × 3 J ¯ i ( i )
v ^ ¯ i = ( v ¯ i T     ω ¯ i ( i ) T ) T = ( x ˙ i y ˙ i z ˙ i ω i 1 ( i ) ω i 2 ( i ) ω i 3 ( i ) ) T ( i = 1 , 2 , , 6 )
v ^ ˙ ¯ i = ( v ˙ ¯ i T     ω ˙ ¯ i ( i ) T ) T = ( x ¨ i y ¨ i z ¨ i ω ˙ i 1 ( i ) ω ˙ i 2 ( i ) ω ˙ i 3 ( i ) ) T ( i = 1 , 2 , , 6 )
v ^ ˜ ¯ i = ω ˜ ¯ i ( i ) 0 ¯ 3 × 3 v ˜ ¯ i ω ˜ ¯ i ( i ) ( i = 1 , 2 , , 6 )
F ^ ¯ i = F ¯ i g k = 1 n e S i k e F ¯ k e k = 1 n e C ¯ i k e ( i ) × F ¯ k e ( i ) k = 1 n e S i k M ¯ k a ( i ) ( i = 1 , 2 , , 6 )
x ˙ i y ˙ i z ˙ i and x ¨ i y ¨ i z ¨ i are the components of B i   ’s velocity v i and acceleration v ˙ i in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction, respectively; ω i 1 ( i ) ω i 2 ( i ) ω i 3 ( i ) and ω ˙ i 1 ( i ) ω ˙ i 2 ( i ) ω ˙ i 3 ( i ) are the components of B i   ‘s angular velocity ω i and angular acceleration ω ˙ i in the e 1 ( i ) , e 2 ( i ) , e 3 ( i ) direction, respectively (Superscript ( i ) represents the expression relative to ( O i , e ¯ ( i ) ) , and no superscript represents the expression relative to ( O i , e ¯ ( 0 ) ) ); F i g is the gravity force of B i ; F i e , and M i a are the force element and active moment that B i 1   acts on B i   , respectively, and there are n e force elements and moments associated with B i ; the body-hinge vector C i k is the vector between the hinge point and the centroid, and the directions are from the inside to the outside, namely: C i k = 1 2 l i e 1 ( i )     ( i = 1 , 2 , , 6 ;     k = i , i + 1 ) ; the 3-order antisymmetric matrices ω ˜ ¯ i ( i ) and v ˜ ¯ i and the incidence element S i j are defined as
ω ˜ ¯ i ( i ) = 0 ω i 3 ( i ) ω i 2 ( i ) ω i 3 ( i ) 0 ω i 1 ( i ) ω i 2 ( i ) ω i 1 ( i ) 0 ,     v ˜ ¯ i = 0 z ˙ i y ˙ i z ˙ i 0 x ˙ i y ˙ i x ˙ i 0
S i j = 1 O j   has   incidence   with   B i   and   starts   with   B i 1 O j   has   incidence   with   B i   and   ends   with   B i 0 O j   has   no   incidence   with   B i             i , j = 1 , 2 , , n
In graph theory, the relationship between O j   and B i   is called incidence, O j   starts with B i   means that B i   is the starting point of O j   , and O j   ends with B i   means that B i   is the end point of O j   .

3.2. Dynamic Equation Based on Quaternion

Since the integral of angular velocity has no physical meaning, it cannot represent the angle. In practical applications, Euler angles, Cardano angles, Rodrigues parameters, or quaternions are selected as generalized angular coordinates. The first three all have singularities and are suitable for small-angle range motion. The quaternion has no singularity, is suitable for motion in a large angle range, and does not need to calculate the sine and cosine of the angle, which reduces the calculation steps and improves the calculation efficiency. For multi-DOF robots, quaternion coordinates are a better choice.
The angle in the dynamic Equation (1) of B i   is expressed by the quaternion coordinate Λ ¯ i = ( λ i 0 λ i 1 λ i 2 λ i 3 ) T , and the quaternion and the position coordinate r ¯ i = ( x i y i z i ) T form the generalized coordinate q ¯ i = ( r ¯ i T Λ ¯ i T ) T . The generalized velocity q ˙ ¯ i = ( r ˙ ¯ i T Λ ˙ ¯ i T ) T and the generalized acceleration q ¨ ¯ i = ( r ¨ ¯ i T Λ ¨ ¯ i T ) T are derived from the first derivative Λ ˙ ¯ i = ( λ ˙ i 0 λ ˙ i 1 λ ˙ i 2 λ ˙ i 3 ) T and the second derivative Λ ¨ ¯ i = ( λ ¨ i 0 λ ¨ i 1 λ ¨ i 2 λ ¨ i 3 ) T of Λ ¯ i .
Twist v ^ ¯ i and acceleration screw v ^ ˙ ¯ i have the following relationships with generalized velocity q ˙ ¯ i and generalized acceleration q ¨ ¯ i ,
v ^ ¯ i = J ¯ i v ( q ) q ˙ ¯ i ( i = 1 , 2 , , 6 )
v ^ ˙ ¯ i = J ¯ i v ( q ) q ¨ ¯ i ( i = 1 , 2 , , 6 )
Here, J i v ( q ) is the Jacobian matrix
J ¯ i v ( q ) = 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 2 λ i 1 2 λ i 0 2 λ i 3 2 λ i 2 0 0 0 2 λ i 2 2 λ i 3 2 λ i 0 2 λ i 1 0 0 0 2 λ i 3 2 λ i 2 2 λ i 1 2 λ i 0
Therefore, the dynamic Equation (1) is expressed as
ψ ¯ i J ¯ i v ( q ) q ¨ ¯ i + v ^ ˜ ¯ i ( q ) ψ ¯ i J ¯ i v ( q ) q ˙ ¯ i = F ^ ¯ i ( i = 1 , 2 , , 6 )
Here, v ^ ˜ ¯ i ( q ) is the generalized coordinate representation of v ^ ˜ ¯ i ,
v ^ ˜ ¯ i ( q ) = 2 R ¯ i R ˙ ¯ i T 0 ¯ 3 × 3 v ˜ ¯ i 2 R ¯ i R ˙ ¯ i T ( i = 1 , 2 , , 6 )
R ¯ i ,   R ˙ ¯ i T are defined as
R ¯ i = λ i 1 λ i 0 λ i 3 λ i 2 λ i 2 λ i 3 λ i 0 λ i 1 λ i 3 λ i 2 λ i 1 λ i 0 ,   R ˙ ¯ i T = λ ˙ i 1 λ ˙ i 2 λ ˙ i 3 λ ˙ i 0 λ ˙ i 3 λ ˙ i 2 λ ˙ i 3 λ ˙ i 0 λ ˙ i 1 λ ˙ i 2 λ ˙ i 1 λ ˙ i 0 ( i = 1 , 2 , , 6 )
The unconstrained rigid body has 3 rotational DOFs, while the quaternion has four parameters, so the quaternions are not independent parameters and have the following constraints,
Λ ¯ i T Λ ¯ i 1 = 0 ( i = 1 , 2 , , 6 )
In order to be combined with Equation (12), that is, to add constraint, Equation (15) needs to be derived twice with respect to time t , and rewritten into the generalized coordinate form, where the following constraint is obtained,
Λ ¯ q i T q ¨ ¯ i + Λ ˙ ¯ q i T q ˙ ¯ i = 0 ( i = 1 , 2 , , 6 )
Here,
Λ ¯ q i T = 0 ¯ 1 × 3 Λ ¯ i T ,     Λ ˙ ¯ q i T = 0 ¯ 1 × 3 Λ ˙ ¯ i T ( i = 1 , 2 , , 6 )
Combining Equation (12) with Equation (17), that is, adding the quaternion self-constraint, the dynamic equation expressed by the generalized coordinate of B i   is obtained,
ψ ¯ i J ¯ i v ( q ) Λ ¯ q i T q ¨ ¯ i + v ^ ˜ ¯   i ( q ) ψ ¯ i J ¯ i v ( q ) Λ ˙ ¯ q i T q ˙ ¯ i = F ^ ¯ i 0 ( i = 1 , 2 , , 6 )
Equation (18) can be abbreviated as
ψ ¯   i ( q )   q ¨ ¯ i = κ ¯ i ( q ) ( i = 1 , 2 , , 6 )
Here,
ψ ¯ i   ( q ) = ψ ¯ i J ¯ i v ( q ) Λ ¯ q i T ,   κ ¯ i   ( q ) = F ^ ¯ i 0 v ^ ˜ ¯   i ( q ) ψ ¯ i J ¯ i v ( q ) Λ ˙ ¯ q i T q ˙ ¯ i
The dynamic equations Equation (19) of the 6 rigid bodies are combined to obtain the dynamic equation of the robot.
ψ ¯ ( q ) q ¨ ¯ = κ ¯ ( q )
Here,
ψ ¯ ( q ) = d i a g ( ψ ¯ 1 ( q ) ψ ¯ 2 ( q )     ψ ¯ 6 ( q ) ) , q ¨ ¯ = ( q ¨ ¯ 1 T q ¨ ¯ 2 T     q ¨ ¯ 6 T ) T , κ ¯ ( q ) = ( κ ¯ 1 ( q ) T   κ ¯ 2 ( q ) T     κ ¯ 6 ( q ) T ) T
ψ ¯ ( q ) is a square matrix of order 42, and q ¨ ¯ and κ ¯ ( q ) are column matrices of order 42. The dynamic model Equation (21) of the robot has a total of 42 equations. The screw theory expresses the dynamic model of the robot as a set of equations. The form is simple and compact, which is convenient for calculation and control.

4. Quaternion-Based Constraint Equations

In the process of motion, there are certain relations between the generalized coordinates that describing the poses of each rigid body, and the analytical expressions of these relations are called constraint equations. A robot that satisfies the constraint equations can be controlled to move according to the true trajectory. The constraint equations of the robot are established by hinge constraints and quaternion self-constraints.
The robot has 6 rotation hinges, each of which has 1 rotational DOF and provides 5 constraints, namely 3 movement constraints and 2 rotation constraints. Furthermore, the angular coordinates that expressed by quaternions have their self-constraints. Specifically reflected as
(1) The hinge point of B j and its inside rigid body B j 1 are the same point O j ,
r ¯ i ( 0 ) + A ¯ ( 0 i ) C ¯ i j ( i ) r ¯ j ( 0 ) + A ¯ ( 0 j ) C ¯ j j ( j ) = 0
A ¯ ( 0 i ) is the coordinate transformation matrix that ( O i , e ¯ ( i ) ) relative to ( O 0 , e ¯ ( 0 ) ) ,
A ¯ ( 0 i ) = 2 ( λ i 0 2 + λ i 1 2 ) 1 2 ( λ i 1 λ i 2 λ i 0 λ i 3 ) 2 ( λ i 1 λ i 3 + λ i 0 λ i 2 ) 2 ( λ i 1 λ i 2 + λ i 0 λ i 3 ) 2 ( λ i 0 2 + λ i 2 2 ) 1 2 ( λ i 2 λ i 3 λ i 0 λ i 1 ) 2 ( λ i 1 λ i 3 λ i 0 λ i 2 ) 2 ( λ i 2 λ i 3 + λ i 0 λ i 1 ) 2 ( λ i 0 2 + λ i 3 2 ) 1
Equation (23) can be projected to the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction to obtain three projection expressions, corresponding to 3 movement constraints.
(2) The base vector of the rotation axis of B j and its inside rigid body B j 1 is the same base vector p j ,
A ¯ ( 0 i ) p _ j ( i ) A ¯ ( 0 j ) p _ j ( j ) = 0
Equation (25) can be projected to the direction e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) to obtain three projection expressions. Since the projections p j 1 , p j 2 , p j 3 of the base vector p j in the directions e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) have the relationship p j 1 2 + p j 2 2 + p j 3 2 = 1 , only two of the three projection expressions are independent. Only take the projection expressions in the e 1 ( 0 ) , e 2 ( 0 ) directions, corresponding to 2 rotation constraints.
(3) Quaternions are not independent parameters and are constrained by quaternion self-constraint Equation (15).
The constraint equations of the robot are obtained according to the projection expressions.
Φ 11 = x 1 + l 1 2 ( 2 ( λ 10 2 + λ 11 2 ) 1 ) = 0
Φ 12 = y 1 + l 1 ( λ 11 λ 12 + λ 10 λ 13 ) = 0
Φ 13 = z 1 + l 1 ( λ 11 λ 13 λ 10 λ 12 ) = 0
Φ 14 = 1 2 ( λ 10 2 + λ 11 2 ) = 0
Φ 15 = 2 ( λ 11 λ 12 + λ 10 λ 13 ) = 0
Φ 16 = λ 10 2 + λ 11 2 + λ 12 2 + λ 13 2 1 = 0
Φ 21 = x 1 + l 1 2 ( 2 ( λ 10 2 + λ 11 2 ) 1 ) x 2 + l 2 2 ( 2 ( λ 20 2 + λ 21 2 ) 1 ) = 0
Φ 22 = y 1 + l 1 ( λ 11 λ 12 + λ 10 λ 13 ) y 2 + l 2 ( λ 21 λ 22 + λ 20 λ 23 ) = 0
Φ 23 = z 1 + l 1 ( λ 11 λ 13 λ 10 λ 12 ) z 2 + l 2 ( λ 21 λ 23 λ 20 λ 22 ) = 0
Φ 24 = 2 ( λ 11 λ 12 λ 10 λ 13 λ 21 λ 22 + λ 20 λ 23 ) = 0
Φ 25 = 2 ( λ 10 2 + λ 12 2 λ 20 2 λ 22 2 ) = 0
Φ 26 = λ 20 2 + λ 21 2 + λ 22 2 + λ 23 2 1 = 0
Φ 31 = x 2 + l 2 2 ( 2 ( λ 20 2 + λ 21 2 ) 1 ) x 3 + l 3 2 ( 2 ( λ 30 2 + λ 31 2 ) 1 ) = 0
Φ 32 = y 2 + l 2 ( λ 21 λ 22 + λ 20 λ 23 ) y 3 + l 3 ( λ 31 λ 32 + λ 30 λ 33 ) = 0
Φ 33 = z 2 + l 2 ( λ 21 λ 23 λ 20 λ 22 ) z 3 + l 3 ( λ 31 λ 33 λ 30 λ 32 ) = 0
Φ 34 = 2 ( λ 21 λ 22 λ 20 λ 23 λ 31 λ 32 + λ 30 λ 33 ) = 0
Φ 35 = 2 ( λ 20 2 + λ 22 2 λ 30 2 λ 32 2 ) = 0
Φ 36 = λ 30 2 + λ 31 2 + λ 32 2 + λ 33 2 1 = 0
Φ 41 = x 3 + l 3 2 ( 2 ( λ 30 2 + λ 31 2 ) 1 ) x 4 + l 4 2 ( 2 ( λ 40 2 + λ 41 2 ) 1 ) = 0  
Φ 42 = y 3 + l 3 ( λ 31 λ 32 + λ 30 λ 33 ) y 4 + l 4 ( λ 41 λ 42 + λ 40 λ 43 ) = 0
Φ 43 = z 3 + l 3 ( λ 31 λ 33 λ 30 λ 32 ) z 4 + l 4 ( λ 41 λ 43 λ 40 λ 42 ) = 0
Φ 44 = 2 ( λ 30 2 + λ 31 2 λ 40 2 λ 41 2 ) = 0
Φ 45 = 2 ( λ 31 λ 32 + λ 30 λ 33 λ 41 λ 42 λ 40 λ 43 ) = 0
Φ 46 = λ 40 2 + λ 41 2 + λ 42 2 + λ 43 2 1 = 0
Φ 51 = x 4 + l 4 2 ( 2 ( λ 40 2 + λ 41 2 ) 1 ) x 5 + l 5 2 ( 2 ( λ 50 2 + λ 51 2 ) 1 ) = 0
Φ 52 = y 4 + l 4 ( λ 41 λ 42 + λ 40 λ 43 ) y 5 + l 5 ( λ 51 λ 52 + λ 50 λ 53 ) = 0
Φ 53 = z 4 + l 4 ( λ 41 λ 43 λ 40 λ 42 ) z 5 + l 5 ( λ 51 λ 53 λ 50 λ 52 ) = 0
Φ 54 = 2 ( λ 41 λ 42 λ 40 λ 43 λ 51 λ 52 + λ 50 λ 53 ) = 0
Φ 55 = 2 ( λ 40 2 + λ 42 2 λ 50 2 λ 52 2 ) = 0
Φ 56 = λ 50 2 + λ 51 2 + λ 52 2 + λ 53 2 1 = 0
Φ 61 = x 5 + l 5 2 ( 2 ( λ 50 2 + λ 51 2 ) 1 ) x 6 + l 6 2 ( 2 ( λ 60 2 + λ 61 2 ) 1 ) = 0
Φ 62 = y 5 + l 5 ( λ 51 λ 52 + λ 50 λ 53 ) y 6 + l 6 ( λ 61 λ 62 + λ 60 λ 63 ) = 0
Φ 63 = z 5 + l 5 ( λ 51 λ 53 λ 50 λ 52 ) z 6 + l 6 ( λ 61 λ 63 λ 60 λ 62 ) = 0
Φ 64 = 2 ( λ 50 2 + λ 51 2 λ 60 2 λ 61 2 ) = 0
Φ 65 = 2 ( λ 51 λ 52 + λ 50 λ 53 λ 61 λ 62 λ 60 λ 63 ) = 0
Φ 66 = λ 60 2 + λ 61 2 + λ 62 2 + λ 63 2 1 = 0
Each rigid body is subject to 5 kinematic constraints and 1 quaternion self-constraint, and a robot consisting of 6 rigid bodies has a total of 36 constraints. The dynamic Equation (21) of the robot has 42 equations and after being subjected to 36 constraints, the DOF is 6, which is consistent with the actual situation.

5. Constrained Dynamic Model of Robot

Only if the robot meets the constraints can its motion be in accord with the objective laws. The constraint equations are added to the dynamic Equation (21) by the Lagrange multiplier method, and the constrained dynamic model is obtained. The Lagrange multipliers are proportional to the ideal constraining force or moment.
The constraint equation of the robot is written as
Φ _ ( q ¯ ) = ( Φ _ 1 T     Φ _ 2 T     Φ _ 3 T     Φ _ 4 T     Φ _ 5 T     Φ _ 6 T ) T
Here, Φ _ i = ( Φ i 1     Φ i 2     Φ i 3     Φ i 4     Φ i 5     Φ i 6 ) T     ( i = 1 , 2 , , 6 ) .
In order to be combined with Equation (21), the constraint Equation (62) needs to be written in the generalized acceleration form, therefore the Equation (62) is derived twice with respect to time t , and the following equation is obtained,
Φ ¯ q q ¨ ¯ + ( Φ ¯ q q ˙ ¯ ) q q ˙ ¯ = 0 ¯
Here, Φ _ q is the Jacobian matrix of Φ _ ( q ¯ ) , and the generalized coordinates of each rigid body constitute the robot’s generalized coordinate q ¯ = ( q ¯ 1 , q ¯ 2 , , q ¯ 6 ) T = ( q 1 , q 2 , , q 42 ) T .
Φ ¯ q = Φ 11 q 1 Φ 11 q 42 Φ 66 q 1 Φ 66 q 42
( Φ ¯ q q ˙ ¯ ) q is the Jacobian matrix of Φ ¯ q q ˙ ¯ .
According to the Lagrange multiplier method, first, 36 Lagrange multipliers λ i j           ( i , j = 1 , 2 , , 6 ) are introduced into Equation (21) and combined into 6 Lagrange multiplier matrices λ ¯ i = ( λ i 1 , λ i 2 , , λ i 6 ) i = 1 , 2 , , 6 , corresponding to the rigid body B i i = 1 , 2 , , 6 , respectively. λ ¯ i constitutes the Lagrange matrix λ ¯ = ( λ ¯ 1 , λ ¯ 2 , , λ ¯ 6 ) T of the robot, from which the Lagrange equation of first kind is obtained.
ψ ¯ ( q ) q ¨ ¯ κ ¯ ( q ) + Φ ¯ q T λ ¯ = 0
Second, add constraints to Equation (65), that is, combine Equation (65) with constraint Equation (63).
ψ ¯ ( q ) Φ ¯ q T Φ ¯ q 0 ¯ q ¨ ¯ λ ¯ = κ ¯ ( q ) ζ ¯
Here, ζ ¯ = ( Φ ¯ q q ˙ ¯ ) q q ˙ ¯ .
Equation (66) is the constrained dynamic model of the robot, with a total of 78 equations, which can solve up to 78 unknowns. The expression form of this equation is compact and unified, and the modular expression is easy to control. The force control or drive control of the robot can be realized by planning the trajectory curve of q ¨ ¯ or λ ¯ , respectively. The Lagrange multiplier is proportional to the ideal constraining force or moment. By solving the Lagrange multipliers, the ideal constraining force and moment can be estimated, which can provide data and a theoretical basis for the analysis of the structural strength, reliability, and fatigue life of the robot.

6. Numerical Simulation Example

The robot is numerically simulated by Adams software to verify the correctness and accuracy of the dynamic Equation (66). As a precision machine, industrial robots are not easy to disassemble, which makes it impossible to measure the dynamic parameters of the rigid bodies that make up the robot by physical experiments, that is, the mass is obtained by weighing, and the principal moment of inertia is measured by the pendular motions. At the same time, the ideal constraint force and moment cannot be measured. Therefore, the mass m i , the principal moment of inertia J i 11 ( i )       J i 22 ( i )       J i 33 ( i ) ( J i 11 ( i )       J i 22 ( i )       J i 33 ( i ) are the components of J i in the e 1 ( i ) , e 2 ( i ) , e 3 ( i ) direction, respectively), and the Lagrange multiplier λ ¯ of each rigid body are taken as the 78 unknown parameters to be calculated, and the remaining parameters are sampled as known parameters, then the unknown parameters are solved by Equation (66).
First, the three-dimensional model of the robot is established by Adams software and simulation experiments are carried out.
Each rigid body B i that composes the robot is simplified to a homogeneous cylindrical rod, the parameters of length l i and radius R   i are provided, the material type is set to steel, and the value of gravitational acceleration is set to 9.80665. From these, the inertial parameters of the robot are generated, as shown in Table 1.
The purpose of the experiment is to solve the dynamic parameters and Lagrange multipliers of the robot. This means that the dynamic parameters and Lagrange multipliers of each rigid body B i are regarded as unknown parameters to be solved, and the other 6 parameters of each rigid body B i are regarded as known parameters, that is, the measured parameters, which are: (1) quaternion coordinates relative to ( O 0 , e ¯ ( 0 ) ) ; (2) the force element exerted on B i by O i relative to ( O 0 , e ¯ ( 0 ) ) ; (3) the active moment exerted on B i by O i relative to ( O 0 , e ¯ ( 0 ) ) ; (4) the centroid acceleration relative to ( O 0 , e ¯ ( 0 ) ) ; (5) the angular velocity relative to ( O 0 , e ¯ ( 0 ) ) ; (6) the angular acceleration relative to ( O 0 , e ¯ ( 0 ) ) .
In order to simplify the experiment and reduce the computational difficulty, the ideal experiment is to solve the dynamic parameters and Lagrange multipliers of the robot with only one movement, which means that the 6 rigid bodies that make up the robot can generate the 6 known parameters in this movement. Since the rigid bodies are connected to each other by hinges, the movement of the inner rigid body will drive the outer rigid body to move together. If the inner rigid body generates the 6 known parameters during the motion, the outer rigid body will generate the 6 known parameters at the same time. B 1 is connected to the ground and can only generate the angular acceleration in the e 3 ( 0 ) direction, but cannot generate the angular acceleration in the e 1 ( 0 ) , e 2 ( 0 ) direction, so B 2 needs to move at the same time to generate the angular acceleration in the e 1 ( 0 ) , e 2 ( 0 ) direction. B 2 can generate the 6 known parameters in this movement, and B 3 , B 4 , B 5 , B 6 also generates the 6 known parameters at the same time, so B 3 , B 4 , B 5 , B 6 does not need to move relative to the inner rigid body.
Simulation settings: The initial state is static, and the posture is shown in Figure 1. Both hinge O 1 and hinge O 2 rotate at the angular acceleration of 4.7746   r a d / s 2 (namely 15 ° / s 2 ). The direction of the axis of rotation of hinge O 1 is e 3 ( 0 ) , and the direction of the axis of rotation of hinge O 2 is e 2 ( 1 ) , and other hinges are locked so that they do not rotate. The simulation termination time is 2 s, and the change of the robot’s posture during the movement is shown in Figure 2. The trajectories of the 6 measured parameters of each rigid body over time are shown in Appendix A.
Second, the data of the motion process is sampled, and then the unknown parameters are solved.
The 1.2 s time (or any other time) is selected to sample the data of the quaternion, force element F i e , active moment M i a , angular velocity, angular acceleration, and centroid acceleration of each rigid body. The sampled data is introduced into Equation (66) to calculate the results of unknown parameters. The calculation results are shown in Table 2 and Table 3.
m i ( e 1 ( 0 ) ), m i ( e 2 ( 0 ) ), m i ( e 3 ( 0 ) ) are the calculated masses of B i in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) directions, respectively. Since B 1 is fixed to the ground, no acceleration occurs in the e 1 ( 0 ) ,     e 2 ( 0 ) direction, so the mass cannot be obtained in these directions.
Finally, the calculated results are analyzed.
Comparing the true values of inertial parameters in Table 1 with the calculated values of inertial parameters in Table 2, the maximum absolute error of the two is −0.0723 k g m 2 . The error comes from the rounding of the significant figures of the data, and the increase of the number of calculations will lead to the accumulation of errors. The error cannot be eliminated, so the very small error verifies the correctness and accuracy of the dynamic model.
The Lagrange multipliers in Table 3 have clear physical meanings, where λ i 1 , λ i 2 , λ i 3 are proportional to the ideal constraining force of B i in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction respectively; λ i 4 , λ i 5 that introduced by the constraint Equation (25) are proportional to the ideal constraining moments; while λ i 6 introduced by the quaternion self-constraint Equation (15) does not reflect the ideal constraining moment.
The posture of the robot at 1.2 s is shown in Figure 2e. Since B i i , = 1 , 2 , , 6 is affected by gravity, and the inner rigid body bears more gravity than the outer rigid body, therefore λ i 3 i = 1 , 2 , , 6 is much larger than λ i 1 , λ i 2 i = 1 , 2 , , 6 and λ i 3 i = 1 , 2 , , 6 decreases with the increase of serial number i .
The rotation axes of O 2 , O 3 , O 5 are all in the e 2 ( i ) direction and need to bear the torque generated by gravity; the rotation axes of O 1 , O 2 , O 6 are not in the e 2 ( i ) direction and do not need to bear the torque generated by gravity, so the ideal constraining moment of O 2 , O 3 , O 5 is large and the ideal constraining moment of O 1 , O 2 , O 6 is small. Therefore, the value of λ i 4 , λ i 5 i = 2 , 3 , 5 is much larger than that of λ i 4 , λ i 5 i = 1 , 4 , 6 . The number of outer rigid bodies of O 2 , O 3 , O 5 decreases, and the torque they received also decreases, so λ i 4 , λ i 5 i = 2 , 3 , 5 decreases with the increase of serial number i . In the same way, the number of outer rigid bodies of O 1 , O 4 , O 6 decreases, and the torque they received also decreases, so λ i 4 , λ i 5 i = 1 , 4 , 6 decreases with the increase of serial number i .
The data show that each rotating hinge of the robot receives great ideal constraining force in the vertical direction, while it receives very little ideal constraining force in the horizontal direction. When the direction of the hinge’s rotation axis is the same as the direction of the moment of the outer rigid body to the hinge, the hinge is subjected to a large ideal constraining moment; when the two directions are perpendicular, the ideal constraining moment of the hinge is small. The calculated values of Lagrange multipliers conform to objective physical laws, which also verifies the correctness of the dynamic model.

7. Conclusions

The dynamics of 6-DOF robot are systematically studied in this paper, and a dynamic modeling method of multibody system based on screw theory is proposed.
(1)
The dynamic and kinematic parameters of each rigid body can be decoupled by this method, and the dynamic model is concisely and compactly expressed as a modular matrix equation, which is convenient for both the solution of the dynamic model and the control.
(2)
Quaternions are used as generalized angular coordinates in this method, which eliminates singularities and makes the dynamic model suitable for motion in a large angular range. Additionally, the computational efficiency is improved.
(3)
The dynamic model established by this method is correct and accurate, which is verified by the numerical example. This method can provide a theoretical basis for the research on the dynamics of 6-DOF robot and its control method.
The future work is to carry out the research on related control methods.

Author Contributions

Conceptualization, J.C. and S.B.; methodology, J.C.; software, J.C.; validation, J.C. and S.B.; formal analysis, J.C., C.Y. and L.Z.; investigation, J.C. and C.Y.; resources, S.B.; data curation, J.C.; writing—original draft preparation, J.C.; writing—review and editing, J.C. and S.B.; visualization, C.Y.; supervision, S.B., Y.C. and Y.Y.; project administration, J.C. and C.Y.; funding acquisition, S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by “The National Key R&D Program of China (No. 2017YFB1301700)”.

Acknowledgments

We thank Weilong Fan for technical support.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The trajectories of the 6 measured parameters of each rigid body over time.
Figure A1. The quaternion coordinates of each rigid body B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the quaternion coordinates of B 1 ; (b) the quaternion coordinates of B 2 ; (c) the quaternion coordinates of B 3 ; (d) the quaternion coordinates of B 4 ; (e) the quaternion coordinates of B 5 ; (f) the quaternion coordinates of B 6 .
Figure A1. The quaternion coordinates of each rigid body B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the quaternion coordinates of B 1 ; (b) the quaternion coordinates of B 2 ; (c) the quaternion coordinates of B 3 ; (d) the quaternion coordinates of B 4 ; (e) the quaternion coordinates of B 5 ; (f) the quaternion coordinates of B 6 .
Machines 10 00499 g0a1
Here, quaternion0, quaternion1, quaternion2, quaternion3 are the four components of the quaternions, respectively.
Figure A2. The force element exerted on B i by O i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the force element exerted on B 1 by O 1 ; (b) the force element exerted on B 2 by O 2 ; (c) the force element exerted on B 3 by O 3 ; (d) the force element exerted on B 4 by O 4 ; (e) the force element exerted on B 5 by O 5 ; (f) the force element exerted on B 6 by O 6 .
Figure A2. The force element exerted on B i by O i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the force element exerted on B 1 by O 1 ; (b) the force element exerted on B 2 by O 2 ; (c) the force element exerted on B 3 by O 3 ; (d) the force element exerted on B 4 by O 4 ; (e) the force element exerted on B 5 by O 5 ; (f) the force element exerted on B 6 by O 6 .
Machines 10 00499 g0a2
Here, FX, FY, FZ are the three components of the force element in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction, respectively.
Figure A3. The active moment exerted on B i by O i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the active moment exerted on B 1 by O 1 ; (b) the active moment exerted on B 2 by O 2 ; (c) the active moment exerted on B 3 by O 3 ; (d) the active moment exerted on B 4 by O 4 ; (e) the active moment exerted on B 5 by O 5 ; (f) the active moment exerted on B 6 by O 6 .
Figure A3. The active moment exerted on B i by O i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the active moment exerted on B 1 by O 1 ; (b) the active moment exerted on B 2 by O 2 ; (c) the active moment exerted on B 3 by O 3 ; (d) the active moment exerted on B 4 by O 4 ; (e) the active moment exerted on B 5 by O 5 ; (f) the active moment exerted on B 6 by O 6 .
Machines 10 00499 g0a3
Here, TX, TY, TZ are the three components of the active moment in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction, respectively.
Figure A4. The centroid acceleration of B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the centroid acceleration of B 1 ; (b) the centroid acceleration of B 2 ; (c) the centroid acceleration of B 3 ; (d) the centroid acceleration of B 4 ; (e) the centroid acceleration of B 5 ; (f) the centroid acceleration of B 6 .
Figure A4. The centroid acceleration of B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the centroid acceleration of B 1 ; (b) the centroid acceleration of B 2 ; (c) the centroid acceleration of B 3 ; (d) the centroid acceleration of B 4 ; (e) the centroid acceleration of B 5 ; (f) the centroid acceleration of B 6 .
Machines 10 00499 g0a4
Here, cm_ACCX, cm_ACCY, cm_ACCZ are the three components of the centroid acceleration in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction, respectively.
Figure A5. The angular velocity of B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the angular velocity of B 1 ; (b) the angular velocity of B 2 ; (c) the angular velocity of B 3 ; (d) the angular velocity of B 4 ; (e) the angular velocity of B 5 ; (f) the angular velocity of B 6 .
Figure A5. The angular velocity of B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the angular velocity of B 1 ; (b) the angular velocity of B 2 ; (c) the angular velocity of B 3 ; (d) the angular velocity of B 4 ; (e) the angular velocity of B 5 ; (f) the angular velocity of B 6 .
Machines 10 00499 g0a5
Here, WX, WY, WZ are the three components of the angular velocity in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction, respectively.
Figure A6. The angular acceleration of B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the angular acceleration of B 1 ; (b) the angular acceleration of B 2 ; (c) the angular acceleration of B 3 ; (d) the angular acceleration of B 4 ; (e) the angular acceleration of B 5 ; (f) the angular acceleration of B 6 .
Figure A6. The angular acceleration of B i relative to ( O 0 , e ¯ ( 0 ) ) . (a) the angular acceleration of B 1 ; (b) the angular acceleration of B 2 ; (c) the angular acceleration of B 3 ; (d) the angular acceleration of B 4 ; (e) the angular acceleration of B 5 ; (f) the angular acceleration of B 6 .
Machines 10 00499 g0a6
Here, WDX, WDY, WDZ are the three components of the angular acceleration in the e 1 ( 0 ) , e 2 ( 0 ) , e 3 ( 0 ) direction, respectively.

References

  1. Khalil, W.; Guegan, S. Inverse and direct dynamic modeling of Gough-Stewart robots. IEEE Trans. Robot. 2004, 20, 754–761. [Google Scholar] [CrossRef] [Green Version]
  2. Djuric, A.M.; El Maraghy, W.H. Automatic separation method for generation of reconfigurable 6R robot dynamics equations. Int. J. Adv. Manuf. Technol. 2010, 46, 831–842. [Google Scholar] [CrossRef]
  3. Ibrahim, O.; Khalil, W. Inverse and direct dynamic models of hybrid robots. Mech. Mach. Theory 2010, 45, 627–640. [Google Scholar] [CrossRef]
  4. Koopaee, M.J.; Pretty, C.; Classens, K.; Chen, X. Dynamical modeling and control of modular snake robots with series elastic actuators for pedal wave locomotion on uneven terrain. J. Mech. Des. 2020, 142, 031120. [Google Scholar] [CrossRef]
  5. Kalyoncu, M. Mathematical modelling and dynamic response of a multi-straight-line path tracing flexible robot manipulator with rotating-prismatic joint. Appl. Math. Model. 2008, 32, 1087–1098. [Google Scholar] [CrossRef]
  6. Li, H.; Yang, Z.; Huang, T. Dynamics and elasto-dynamics optimization of a 2-DOF planar parallel pick-and-place robot with flexible links. Struct. Multidiscip. Optim. 2009, 38, 195–204. [Google Scholar] [CrossRef]
  7. Dai, J.S. Screw Algebra and Kinematic Approaches for Mechanisms and Robotics; Springer: London, UK, 2014. [Google Scholar]
  8. Dai, J.S. Geometrical Foundations and Screw Algebra for Mechanisms and Robotics; Higher Education Press: Beijing, China, 2014. [Google Scholar]
  9. Sun, Y.; Wang, S.; Mills, J.K.; Zhi, C. Kinematics and dynamics of deployable structures with scissor-like-elements based on screw theory. Chin. J. Mech. Eng. 2014, 27, 655–662. [Google Scholar] [CrossRef]
  10. Müller, A. Screw and lie group theory in multibody dynamics. Multibody Syst. Dyn. 2018, 42, 219–248. [Google Scholar] [CrossRef] [Green Version]
  11. Cibicik, A.; Egeland, O. Dynamic modelling and force analysis of a knuckle boom crane using screw theory. Mech. Mach. Theory 2019, 133, 179–194. [Google Scholar] [CrossRef]
  12. Liu, W.; Gong, Z.; Wang, Q. Investigation on Kane dynamic equations based on screw theory for open-chain manipulators. Appl. Math. Mech. 2005, 26, 627–635. [Google Scholar]
  13. Pardos-Gotor, J.M. Screw Theory in Robotics: An Illustrated and Practicable Introduction to Modern Mechanics; CRC Press: Boca Raton, FL, USA, 2021. [Google Scholar]
  14. Lynch, K.M.; Park, F.C. Modern Robotics; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  15. Ajwad, S.A.; Ullah, M.I.; Islam, R.U.; Iqbal, J. Modeling robotic arms–A review and derivation of screw theory based kinematics. In Proceedings of the International Conference on Engineering and Emerging Technologies, Lahore, Pakistan, 1–2 December 2014; p. 98. [Google Scholar]
  16. Gallardo, J.; Rico, J.M.; Frisoli, A.; Checcacci, D.; Bergamasco, M. Dynamics of parallel manipulators by means of screw theory. Mech. Mach. Theory 2003, 38, 1113–1131. [Google Scholar] [CrossRef]
  17. Zhao, T.; Geng, M.; Chen, Y.; Li, E.; Yang, J. Kinematics and dynamics Hessian matrices of manipulators based on screw theory. Chin. J. Mech. Eng. 2015, 28, 226–235. [Google Scholar] [CrossRef]
  18. Gallardo-Alvarado, J.; Aguilar-Nájera, C.R.; Casique-Rosas, L.; Rico-Martínez, J.M.; Islam, M.N. Kinematics and dynamics of 2 (3-RPS) manipulators by means of screw theory and the principle of virtual work. Mech. Mach. Theory 2008, 43, 1281–1294. [Google Scholar] [CrossRef]
  19. Gallardo-Alvarado, J.; Rodríguez-Castro, R.; Delossantos-Lara, P.J. Kinematics and dynamics of a 4-PRUR Schönflies parallel manipulator by means of screw theory and the principle of virtual work. Mech. Mach. Theory 2018, 122, 347–360. [Google Scholar] [CrossRef]
  20. Gallardo-Alvarado, J.; Aguilar-Nájera, C.R.; Casique-Rosas, L.; Pérez-González, L.; Rico-Martínez, J.M. Solving the kinematics and dynamics of a modular spatial hyper-redundant manipulator by means of screw theory. Multibody Syst. Dyn. 2008, 20, 307–325. [Google Scholar] [CrossRef]
  21. Ding, X.; Selig, J.M. Dynamic modelling of a compliant arm with 6-dimensional tip forces using screw theory. Robotica 2003, 21, 193–197. [Google Scholar] [CrossRef]
Figure 1. 6-DOF robot and its schematic of the reference systems (a) 6-DOF robot; (b) schematic of the robot’s reference systems.
Figure 1. 6-DOF robot and its schematic of the reference systems (a) 6-DOF robot; (b) schematic of the robot’s reference systems.
Machines 10 00499 g001
Figure 2. The postures at each time of the robot’s simulation model. (a) 0 s; (b) 0.3 s; (c) 0.6 s; (d) 0.9 s; (e) 1.2 s; (f) 1.5 s; (g) 1.8 s; (h) 2 s.
Figure 2. The postures at each time of the robot’s simulation model. (a) 0 s; (b) 0.3 s; (c) 0.6 s; (d) 0.9 s; (e) 1.2 s; (f) 1.5 s; (g) 1.8 s; (h) 2 s.
Machines 10 00499 g002aMachines 10 00499 g002b
Table 1. Physical parameters of robot.
Table 1. Physical parameters of robot.
Rigid   Body   B i B 1 B 2 B 3 B 4 B 5 B 6
Length   l i ( m )1.0451.2451.210.3150.290.3
Radius   R   i ( m )0.50.30.20.20.20.3
m i   ( k g )6402.60122746.07261186.1661308.7953284.2877661.7042
J i 11 ( i )   ( k g m 2 )800.3251123.573323.72336.17595.685829.7767
J i 22 ( i )   ( k g m 2 )982.8126416.4934156.58385.64134.835319.8511
J i 33 ( i )   ( k g m 2 )982.8126416.4934156.58385.64134.835319.8511
Table 2. The solved inertial parameters.
Table 2. The solved inertial parameters.
Rigid   Body   B i B 1 B 2 B 3 B 4 B 5 B 6
m i   ( e 1 ( 0 ) )   ( k g )2746.07261186.1658308.7950284.2880661.7042
m i   ( e 2 ( 0 ) )   ( k g )2746.06851186.1663308.7954284.2878661.7042
m i   ( e 3 ( 0 ) )   ( k g )6402.59822746.07271186.1656308.7955284.2878661.7042
J i 11 ( i )   ( k g m 2 )800.3261123.568123.71406.17935.684729.7768
J i 22 ( i )   ( k g m 2 )982.8122416.4211156.59205.64174.836419.8510
J i 33 ( i )   ( k g m 2 )982.8122416.5364156.58405.64064.835319.8511
Table 3. The solved Lagrange multipliers.
Table 3. The solved Lagrange multipliers.
Lagrange   Multipliers   λ ¯ i λ ¯ 1 λ ¯ 2 λ ¯ 3 λ ¯ 4 λ ¯ 5 λ ¯ 6
λ i 1 1.5166 × 10 6 1.5166 × 10 6 4.1854 × 10 7 1.6008 × 10 5 3.8787 × 10 5 2.6473 × 10 5
λ i 2 5.7217 × 10 4 5.7217 × 10 4 1.3463 × 10 4 5.8943 × 10 5 1.8120 × 10 5 1.3748 × 10 6
λ i 3 1.1366 × 10 5 5.0867 × 10 4 2.3938 × 10 4 1.2305 × 10 4 9.2770 × 10 3 6.4891 × 10 3
λ i 4 4.6041 × 10 5 7.7488 × 10 13 7.5058 × 10 13 3.4288 × 10 4 8.5585 × 10 12 5.0122 × 10 3
λ i 5 4.8692 × 10 6 4.0621 × 10 14 3.9347 × 10 14 6.5407 × 10 3 4.4865 × 10 13 9.5612 × 10 2
λ i 6 8.1974 × 10 14 2.5700 × 10 13 7.9403 × 10 14 9.0539 × 10 13 9.0540 × 10 13 1,0207 × 10 4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cheng, J.; Bi, S.; Yuan, C.; Cai, Y.; Yao, Y.; Zhang, L. Dynamic Modeling Method of Multibody System of 6-DOF Robot Based on Screw Theory. Machines 2022, 10, 499. https://0-doi-org.brum.beds.ac.uk/10.3390/machines10070499

AMA Style

Cheng J, Bi S, Yuan C, Cai Y, Yao Y, Zhang L. Dynamic Modeling Method of Multibody System of 6-DOF Robot Based on Screw Theory. Machines. 2022; 10(7):499. https://0-doi-org.brum.beds.ac.uk/10.3390/machines10070499

Chicago/Turabian Style

Cheng, Jun, Shusheng Bi, Chang Yuan, Yueri Cai, Yanbin Yao, and Ling Zhang. 2022. "Dynamic Modeling Method of Multibody System of 6-DOF Robot Based on Screw Theory" Machines 10, no. 7: 499. https://0-doi-org.brum.beds.ac.uk/10.3390/machines10070499

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop