iBet uBet web content aggregator. Adding the entire web to your favor.
iBet uBet web content aggregator. Adding the entire web to your favor.



Link to original content: https://doi.org/10.3390/robotics10030083
Effects of Temperature and Mounting Configuration on the Dynamic Parameters Identification of Industrial Robots
Next Article in Journal
Model Predictive Control for Cooperative Transportation with Feasibility-Aware Policy
Previous Article in Journal
Semi-Autonomous Behaviour Tree-Based Framework for Sorting Electric Vehicle Batteries Components
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Temperature and Mounting Configuration on the Dynamic Parameters Identification of Industrial Robots

Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
*
Author to whom correspondence should be addressed.
Robotics 2021, 10(3), 83; https://doi.org/10.3390/robotics10030083
Submission received: 21 May 2021 / Revised: 21 June 2021 / Accepted: 25 June 2021 / Published: 29 June 2021
(This article belongs to the Topic Industrial Robotics)

Abstract

:
Dynamic parameters are crucial for the definition of high-fidelity models of industrial manipulators. However, since they are often partially unknown, a mathematical model able to identify them is discussed and validated with the UR3 and the UR5 collaborative robots from Universal Robots. According to the acquired experimental data, this procedure allows for reducing the error on the estimated joint torques of about 90% with respect to the one obtained using only the information provided by the manufacturer. The present research also highlights how changes in the robot operating conditions affect its dynamic behavior. In particular, the identification process has been applied to a data set obtained commanding the same trajectory multiple times to both robots under rising joints temperatures. Average reductions of the viscous friction coefficients of about 20% and 17% for the UR3 and the UR5 robots, respectively, have been observed. Moreover, it is shown how the manipulator mounting configuration affects the number of the base dynamic parameters necessary to properly estimate the robots’ joints torques. The ability of the proposed model to take into account different mounting configurations is then verified by performing the identification procedure on a data set generated through a digital twin of a UR5 robot mounted on the ceiling.

1. Introduction

Detailed knowledge of robot dynamic parameters can be beneficial for several applications. However, in contrast with the kinematic ones, these values are usually not fully provided by the manufacturer. As an example, Universal Robots, one of the main brands in collaborative robotics, details only the expected values of the mass and the position of the center of mass of each joint/link of its manipulators [1], for a total of 24 (4 × 6) parameters. On the other hand, the remaining 54 (9 × 6), which include links and motors inertia and Coulomb and viscous friction, remain unknown, thus preventing the definition of accurate dynamic models. According to [2,3], this is of primary importance for control algorithms currently adopted for industrial manipulators which often rely on control strategies far more complex than simple Proportional-Integral-Differential (PID) ones [4]. On the other hand, a different approach could be derived from [5], where deterministic artificial intelligence has been effectively used to learn the dynamic properties of an unmanned underwater vehicle for autonomous trajectory generation and control.
Moreover, an imprecise estimate of the torques required to execute the desired trajectory could negatively influence the effectiveness of algorithms used to provide performance indexes to evaluate the energy consumption of a manipulator [6], or to define optimal trajectories to minimize the power required by the robot without compromising its productivity [7,8]. Since these methods use the Lagrange formulation [9], more accurate knowledge of the dynamic parameters of the robot arm could lead to a better estimate of the joints torques.
High-Fidelity (HF) mathematical models could be also used for diagnostics and prognostics purposes. To do so, most of the algorithms currently employed rely on data-driven methods, which need to be trained on large data sets representative of the statistical distribution of the system behavior in both nominal and off-nominal conditions [10]. However, since the Mean Time Between Failures (MTBF) of a robot is in the order of tens of thousands of hours [11], there is a general lack of data associated with operation under degraded health conditions. To overcome this limitation, and quickly generate extensive data sets through which design and evaluate dedicated prognostic routines, simulated faults and failures can be inserted in the HF model of the manipulator. The effectiveness of this approach in industrial robots applications is discussed in [12], where prognostics routines have been applied to predict fault and failures in a roller hemming head mounted on a 5 DOF manipulator. Given these premises, reliable results can be obtained through the exact definition of both the kinematic and dynamic parameters of the robot arm, as well as through the modeling of its components.
Diagnostics and prognostics algorithms are crucial not only to minimize downtimes and economical losses in production or assembly lines, but also to guarantee the safety of the operator who is sharing his workspace with the robot. Despite the old idea of industrial automation where robots are surrounded by fences and cages to which access is usually restricted, in tasks involving Human–Robot Collaboration (HRC), the worker and the machine operate in close proximity as in the case of hand-over applications analyzed in detail in [13]. According to the norms in force, such cooperation is feasible only if the operator’s safety is always guaranteed. A possible solution consists of preventing an impact between the operator and the robot arm by implementing collision avoidance algorithms [14,15]. An example of a practical implementation of this approach has been proposed by [16], where a human–robot assembly task is disclosed. However, especially in the case of continuous cooperation of the worker and the robot, advanced algorithms for human–robot workspace sharing cannot always be developed with a satisfying safety level. As a consequence, norms require the development of the application so that the negative effects of possible undesired events would be minimized. To do so, energy approaches are often used to estimate the kinetic energy of the manipulator, which must always be under a specific threshold defined according to the possible impact zones on the worker’s body. Within this framework, a more reliable risk assessment can be reached by accurately identifying the dynamic parameters of the robot adopted. Offline simulators, effectively used to program the movements of the manipulator, often do not fully take into account its dynamic properties. Then, the kinetic energy of the robot, evaluated over its entire trajectory, does not reflect the one of the real system. This could lead to label as safe an application that is actually not. Within this framework, since the experimental campaign has been carried out using two largely adopted collaborative robots, the UR3 and UR5 manipulators from Universal Robots, the present research also aims to provide a tool for a fast and reliable classification of human–robot collaborative tasks into safe or hazardous ones.
A simple approximation of the dynamic parameters can be obtained from the CAD models of the robot arm provided by the manufacturer. However, since complex components, like the joints, are not modeled in detail and friction cannot be considered, this method is usually not effective. A better solution is proposed in [17], where a sequential identification of the dynamic properties of manipulator links is performed. Nevertheless, this is not optimal since error propagation and accumulation from one step to the next one negatively affect the algorithm performance. On the contrary, by identifying the robot dynamic parameters in a single step, the effect of the error is reduced, and the time required for the entire process is minimized. Due to the better results provided by this approach, recent studies have been using it for several industrial manipulators. As an example, in [18] the dynamic parameters of a four Degrees Of Freedom (DOF) robot have been identified, while in [19] the case of a six DOF robot arm mounted on a linear axis has been analyzed. Other studies, involving seven DOF collaborative robots, can be found in [20,21], where the inertial parameters of an ABB IRB14000 (YuMi) and a Franka Panda Emika, respectively, have been estimated.
To guarantee a complete and reliable identification of the robot dynamic parameters, specifically designed trajectories are commanded to the manipulators. To do so, different approaches, as B-splines adopted in [22] and fifth order polynomials used in [23,24] for a three DOF robot and in [25] for a six DOF manipulator, can be found in the literature. However, other authors [19,26,27] preferred Finite Fourier Series (FFS) due to their advantages in signal processing and noise reduction [28].
The present paper describes the algorithm implemented to identify the dynamic parameters of an industrial manipulator. Similar work can be found in [26], where the dynamic calibration of a UR5 collaborative robot has been carried out. However, simplifying assumptions based on the manipulator geometry have been implemented, leading to the definition of a mathematical model effective only with that specific machine. On the contrary, a more general approach, as in [29], has been used in this research, allowing the proposed methodology to be easily adapted to different industrial robots. In particular, the flexibility and the reliability of the developed mathematical model have been tested during an experimental campaign carried out with the collaborative robots UR3 and UR5.
The novelty of the present work relies on the fact that, despite other similar studies [18,19,20,21,26], the effects of temperature and mounting configuration on the robot dynamics have been taken into account and investigated with the same procedure used for the identification of the robot inertial parameters, while in studies like [30,31], friction coefficients were estimated by exciting one joint at the time through a customized trajectory. Since in the literature the dynamic parameters of a manipulator are identified only once, their variations, caused by changes in the robot working conditions, are not considered, leading to possible wrong estimates of the joints torques. As an example, joints temperatures vary during a task and they can be responsible for up to half the amount of the friction torque [32]. The effect of temperature on the dynamic parameters, and in particular on joints viscous friction coefficients, has been studied and quantified by repeating several times the identification process for both robots.
Moreover, since industrial manipulators are often positioned on the ground or on a horizontal bench, previous studies have been focused on the identification of the base dynamic parameters only in such scenarios [19,26,29]. Nevertheless, due to their high flexibility, reduced weight, and compact dimensions, collaborative robots can be easily mounted in several configurations. The mathematical model has been developed to be easily adapted to different setups to guarantee reliable estimates of the joints torques.
In addition, by mounting the robot such that the axis of its first joint is not aligned with the gravitational acceleration vector, it is possible to identify two additional base dynamic parameters, that are otherwise not identifiable. The importance of the definition of the robot mounting configuration is highlighted by estimating the joints torques of a digital twin of the UR5 mounted on the ceiling using the official offline simulator (URSim) developed by the same manufacturer.

2. Experimental Setup

The research campaign has been carried out using the collaborative robots UR3 and UR5. The commands were given to the robots using port 29999, while the feedback signals were acquired at a frequency of 125 Hz using port 30003 [33]. The data necessary for the proposed application are the joint temperatures, angular positions and velocities, and the motors currents. These last values have been used to calculate the joints torques at each trajectory point k by multiplying the motors currents, the gear ratio (G), equal to 101 for both the UR3 and the UR5, and the torque constants reported in Table 1.
Since the joints angular accelerations are not provided, they were obtained by deriving the angular velocities. To minimize the effect of noise in the data, filtering techniques, as suggested in [19,26], have been adopted for all signals.

3. Mathematical Model of the Industrial Manipulator

By considering the robot as a rigid body, its dynamic behavior can be described by a Lagrangian energy-based approach [9] which leads, for a n-DOF manipulator, to write:
M ( q ) q ¨ + C   ( q ,     q ˙ ) q ˙ + F ( q ˙ ) + G ( q ) = τ
where   M   ( q )     n × n   is the symmetric positive definite mass matrix, C   ( q ,     q ˙ )     n × n is the Coriolis and centripetal coupling matrix, F   ( q ˙ )     n . is the friction force vector, G   ( q )     n is the gravitational force vector, τ     n is the joints torques vector while q ,     q ˙ ,   q ¨     n are, respectively, the generalized joints angular positions, velocities, and accelerations.
However, for the proposed application, it is necessary to rewrite Equation (1) with respect to a set of dynamic parameters p     13 n [26] as:
Y   ( q ,   q ˙ ,   q ¨ ) ·     p = τ
where   Y     n × 13 n is called regression matrix or regressor. The robot geometric parameters used to build Y are the ones specified by the manufacturer [1]. Even though Universal Robots provides, for each robot, the corrections to the Denavit–Hartenberg (DH) parameters as a result of the kinematic calibration procedure, they cannot be implemented in a dynamic model of the manipulator. This is related to the fact that these values have been calculated to optimize the position accuracy of the manipulator Tool Center Point (TCP) via numerical optimization without taking into account their physical meaning. As an example, for the UR5 adopted in this study, some link lengths have negative values which, even if impossible from a physical point of view, allow a correct estimate of TCP poses derived using forward kinematics. Nevertheless, since possible mismatches among nominal and real DH parameters are related to production tolerances and mounting errors, their values are expected to be small. As a consequence, their influence on the calculation of the dynamic parameters and the predicted joints torques can be neglected.
Moreover, the regressor matrix is also a function of the mounting configuration. This has been implemented in Y while defining the kinematic chain of the robot arm. To do this, the first DH matrix, which describes the position and the orientation of the first robot joint with respect to its base, has been multiplied by the 3 × 3 rotation matrix WAB defined as:
A B W = R ( ψ ,     φ ,     θ ) = R o t   ( z ,     θ )   ·   R o t   ( y ,     φ )   ·   R o t   ( x ,     ψ ) A B W   =   c o s   φ     c o s   θ c o s   θ   ·   s i n   φ   ·   s i n   ψ s i n   θ   ·   c o s   ψ c o s   θ   ·   s i n   φ   ·   c o s   ψ   +   s i n   θ   ·   s i n   ψ s i n   θ   ·   c o s   φ s i n   θ   ·   s i n   φ   ·   s i n   ψ + c o s   θ   ·   c o s   ψ s i n   θ   ·   s i n   φ   ·   c o s   ψ     c o s   θ   ·   s i n   ψ s i n   φ c o s   φ   ·   s i n   ψ c o s   φ   ·   c o s   ψ
where ψ , φ , and θ are, respectively, the roll, pitch, and yaw angles describing the orientation between the robot base (B) and the fixed world (W) reference frames. So, in a common scenario where B and W are aligned, as for the robots mounted on the ground, it would be ψ = φ = θ = 0 , leading to WAB = I, where I is the 3 × 3 identity matrix. On the other hand, if the robot would be mounted on the ceiling, for example, it would be φ = θ = 0 and ψ =   π .
In more detail, for a generic trajectory point k, Y is defined as:
[ Y i d ( k ) Y c ( k ) Y v ( k ) Y I , m ( k ) ] { p i d p c p v p I , m } = τ ( k )
where Y i d ( k )     n × 10 n is the regressor block built according to [34]. Besides, for a complete model of the manipulator, Coulomb and viscous friction phenomena are described for each joint i by respectively: Y c ( k ) = d i a g   ( t a n h   ( q ˙ i 0.001 ) )     n × n and Y v ( k ) = d i a g   ( q ˙ i ) n × n . The effects of motors inertia have been implemented in Y I , m ( k ) = d i a g   ( q ¨ i )     n × n . Since each body of the robot arm can only affect the dynamics of the previous ones, the regressor matrix has an upper triangular structure, where the generic element Y a b represents the contribution of the bth body to the dynamics of the ath one. Similarly, p is composed by p i d = [ p i d , 1 ,   p i d , 2 ,   p i d , 3 ,   p i d , 4 ,   p i d , 5 ,   p i d , 6 ] T     10 n where, for the single body i,   p i d , i = [ m ,     m x ,     m y ,     m z ,     J x x ,     J y y ,   J   z z ,     J y z =   J z y ,     J x z =   J z x ,     J x y = J y x ] i     10 , which contains the information about its mass, the position of the center of mass according to the xi, yi, and zi axes, and its moments of inertia. Coulomb and viscous coefficients have been grouped inside p c = [ f c 1 ,     f c 2 ,     f c 3 ,     f c 4 ,     f c 5 ,     f c 6 ]     6   and p v   = [ f v 1 ,     f v 2 ,     f v 3 ,     f v 4 ,     f v 5 ,     f v 6 ]     6 respectively, while motor inertia are stored in p I , m = G 2 ·     [ I m 1 ,     I m 2 ,     I m 3 ,   I m 4 ,   I m 5 ,     I m 6 ]     6 .

3.1. Regressor Reduction Using SVD and QR Decomposition

Since not all the manipulator dynamic parameters are linearly independent, it is necessary to remove all the null columns from the regressor Y to obtain a reduced matrix Y B , so that Y B ·     p B = τ , where p B   is the vector of the base dynamic parameters. To do so, two different approaches, explained in detail in [35], can be adopted: QR decomposition and Singular Value Decomposition (SVD). This last method has been applied to Y ¯ 25     25 n × 13 n   obtained by stacking every regressor   Y n × 13 n   evaluated in one of the 25 random sets of angular positions, velocities, and accelerations as suggested in [19]. Y ¯ 25 is then written as Y ¯ 25 = U Σ V T   , where Σ = diag ( σ i )     13 n × 13 n is a diagonal matrix, whose non-null elements are the singular values σ i of Y ¯ 25 reported in Figure 1a.
By analyzing the trend of the singular values of Y ¯ 25 and its cumulative energy [36] in Figure 1b, it is possible to see, as also reported in [29], how 52 base dynamic parameters must be used to fully describe the robot dynamic behavior. If more are adopted, the condition number of Y ¯ 25 , defined as the ratio between the highest and the lowest singular value σ i , would drop, jeopardizing the identification process [37]. A high value of this index, in fact, makes the base dynamic parameters strongly dependent on slight variations of the signals measured from the robot that are affected by noise.
This information is then used to check the quality of the QR decomposition adopted in this work to reduce the system from Y     6 × 78 to Y B     6 × 52 and from p     78   to p B     52 . This analysis highlights that not all the 78 dynamic parameters of the manipulator (13 for each joint) can be identified. Indeed, only 30 of them are totally identifiable; 39 are identifiable with linear dependency, while the remaining 9 do not play any role in the dynamics of the manipulator.
On the other hand, when mounting the robot in different configurations (i.e., on a wall, being φ = θ = 0 and ψ =   π / 2 ), the SVD analysis applied to Y ¯ 25 assesses that the manipulator dynamics is fully described by 54 base dynamic parameters, leading to p B     54 . Since the rotation axis of the robot base is not aligned with the gravitational acceleration vector, the mass and two coordinates of the position of the center of mass of the first joint/link need to be identified. So, while the number of the totally identifiable dynamic parameters remains unchanged, 42 of them are now identifiable with linear dependency and only 6 do not affect the robot dynamics. The symbolic equations describing the 54 elements of the vectors p B of both the UR3 ( p B ,     U R 3 )   and the UR5 ( p B ,   U R 5 )   are reported in Appendix A. In the case study in which the robot is positioned on the ground or on the ceiling, the 52 base dynamic parameters are still described by Equations (A1) and (A2), but without the first two rows.

3.2. Optimized Excitation Trajectory

For proper identification of the base dynamic parameters p B , a persistent trajectory, built using 5th order Finite Fourier Series (FFS), has been commanded to both the UR3 and the UR5. Angular positions, velocities, and accelerations of each joint i are calculated according to:
q i   t   =   q i , 0   +   l   =   1 5 a i , l   s i n   ω f l t     b i , l   c o s   ω f l t q ˙ i   t   =   l   =   1 5 a i , l   ω f l   c o s   ω f l t   +   b i , l   ω f l   s i n   ω f l t q ¨ i   t   =   l   =   1 5     a i , l   ω f l 2 s i n   ω f l t   +   b i , l   ω f l 2 c o s   ω f l t
where:
  • ω f is the fundamental frequency equal for each joint to guarantee the periodicity of the robot’s movement. It is defined as ω f = 2 π / T , where T is the identification trajectory period set to 10 s as in [29], while similar or lower frequencies have been adopted in [19,26,38];
  • q i , 0 is the joint position offset equal to [0, −π/2, 0, 0, 0, 0] rad;
  • a i , l and b i , l are the coefficients of the FFS that have to be found by optimization to define a trajectory able to continuously excite all the base dynamic parameters.
Since the identification trajectory must be feasible, physical constraints, regarding the maximum joints angular positions ( q m a x ), velocities (   q ˙ m a x ), and accelerations ( q ¨ m a x ) of the manipulator, were added to the minimization problem. Such values were chosen both according to the mechanical constraints of the robots ( q ˙ m a x and q ¨ m a x   are equal to π rad/s and 5.5π rad/s2, respectively) and to avoid any collision of the robots with themselves or with the external environment by setting q m a x   = [2π, π, π, 2π, 2π, 2π] rad. By adapting the equations reported in [39], these limits have been implemented as:
q i t   =   l   =   1 5 a i , l 2 + b i , l 2   +   q i , 0     q i , m a x q ˙ i t   =   l   =   1 5 ω f l     ·   a i , l 2   +   b i , l 2     q ˙ i , m a x q ¨ i t   =   l   =   1 5   ω f l 2 ·     a i , l 2   +   b i , l 2     q ¨ i , m a x
Moreover, since non-zero values of the joints angular velocities and accelerations would lead to high vibrations at the start and the end of the commanded trajectory, their initial and final values have been set to 0 using:
q ˙ i   ( t 0 ) = q ˙ i   ( t e n d ) = l = 1 5 a i , l   ω f l = 0 q ¨ i   ( t 0 ) = q ¨ i   ( t e n d ) = l = 1 5 b i , l     ( ω f l ) 2 = 0
The FFS coefficients a i , l and b i , l can be identified using different approaches like artificial bee colony [38], particle swarm [40], or a Genetic Algorithm (GA) [18,41]. In the present work, the MATLAB Global Optimization Toolbox has been used to implement a GA with a population of 300 individuals to minimize the objective function:
min q , q ˙ , q ¨   c o n d   ( Y ¯ B   ( q ,   q ˙ ,   q ¨ ) )
where Y ¯ B = ( Y B , 1 Y B , N _ i d ) ( n x N _ i d )   x 52 is the observation matrix built by piling the single regressors Y B , k calculated for each trajectory point k of the identification trajectory. Other types of cost functions can be found in [42].
The resulting angular positions obtained for the UR5 collaborative robot are reported in Figure 2a. Similar values have been also found for the UR3.
A condition number of about half the one related to the trajectory of Figure 2a has been obtained for the one reported in Figure 2b in the case the robot would be mounted on the ceiling. This result suggests that more reliable identification of the base dynamic parameters would be obtained by mounting the manipulator on the ceiling rather than on the ground. Nevertheless, as reported in Section 5, good approximations of the joints torques of both robots have been obtained by calculating their base dynamic parameters while mounted on the ground.
To reduce the risk of an ill-conditioned Y ¯ B , the identified trajectory is repeated three times, leading to a total excitation period of 30 s. To do this, it has been necessary to use the servoj function created to directly control the joints angular positions and defined as: servoj(target joint configuration, acceleration, velocity, time, lookahead time, gain) [43]. Since they are not used in the current version of the robot software, the acceleration and the velocity parameters have been set to 0. On the contrary, the time has been set to 0.008 s being the communication frequency between the robot and the external computer equal to 125 Hz. Experimental tests showed that, for a smooth movement of the robot arm, optimal values for the lookahead time and the gain should be 0.03 and 500, respectively. This approach only allows for specifying a set of angular positions to the manipulator, while joints angular velocities and accelerations cannot be directly commanded.

4. Base Dynamic Parameters Identification

As in [19,26,44], to calculate the robots base dynamic parameters p B , the Least Squares (LS) algorithm has been adopted as:
p B = ( Y ¯ B T ·   Y ¯ B ) 1   ·   Y ¯ B T     ·   τ ¯ I D
where τ ¯ I D = ( τ I D , 1 τ I D , N _ i d ) is the joints torques vector derived from the motor currents. These values are reported, together with their standard deviations, in Figure 3a,b for the UR3 and the UR5, respectively. Additional frames have been inserted for better visualization of the small values which cannot be read at full scale.
Since both robots have been mounted on the ground, their dynamic behavior is fully described by 52 base dynamic parameters. In Figure 3a,b are reported, respectively, the numerical values associated with the elements p B , U R 3 ( 3 ) p B , U R 3 ( 54 ) and p B , U R 5 ( 3 ) p B , U R 5 ( 54 ) of the vectors reported in Equations (A1) and (A2).
Moreover, the standard deviation d i = ( [ Y ¯ B Y ¯ B ] 1 ) i , i of each element of p B has been calculated, resulting in maximum values of 0.313 for the UR3 and 0.107 for the UR5.
An alternative to Equation (9), which takes into account the different sizes of the manipulator joints, is proposed in [19]. By applying a Weighted Least Squares (WLS) method, the base dynamic parameters are calculated as:
p B w e i g h t e d = ( Y ¯ B   T ·   W T ·   W   ·   Y ¯ B ) 1   ·   Y ¯ B T ·     W T ·   W   ·   τ ¯ I D
where the coefficients W , for each trajectory point k, are W ( k ) = d i a g ( 1 G   ·   τ j , m a x ) with τ j , m a x being the maximum torques of the joint j, reported in Table 2.
In the proposed application, the LS method has been preferred to the WLS one for two main reasons:
  • While Equation (9) takes 0.05 s to be solved, 138.90 s are required to obtain the weighted base dynamic parameters for the same manipulator using Equation (10). The machine used for all the calculations is a PC with processor Intel® i7-8750H at 3.91 GHz equipped with 16 GB of RAM DDR4;
  • As reported in Table 3, better results in torques estimates have been obtained using the LS method.
A third approach is described in [46], where a Constrained Weighted Least Squares (CWLS) method has been implemented using the constraints associated with links masses and inertia tensors reported in [47]. This ensures a positive definite mass matrix for each trajectory point k. If such condition is not satisfied, the algorithm provides physically impossible results, such as the negative motors inertia ( p B   ( 49 ) p B   ( 52 ) ) of the last four joints of both manipulators. Although more accurate, this method requires a global optimization algorithm to calculate the base dynamic parameters, hence resulting in a far more time-consuming procedure than the LS and WLS. Alternative approaches can be found in [48], where an extensive review of the most used algorithms for dynamic parameters calculation is reported.
Besides, since both Coulomb ( p B   ( 37 ) p B   ( 42 ) ) and viscous friction ( p B   ( 43 ) p B   ( 48 ) ) coefficients are fully identifiable, it would be possible to use the proposed identification algorithm to estimate the health status of all the manipulator joints. An increment of the friction parameters, for example, could be related to a higher level of wear in the joint bearings and gearbox. Even though it would not be possible to locate the root cause of the fault, this information could be used, together with other health features extracted from the robot signals, to estimate the overall operating condition of the machine and its remaining useful life. The reliability of the identified friction coefficients is also highlighted by the similar values of p B   ( 39 ) and p B   ( 40 ) p B   ( 42 ) and of p B   ( 45 ) and p B   ( 46 ) p B   ( 48 ) of the UR3 and the UR5, respectively. Since, as reported in Table 2, the joints mounted at the UR3 elbow and the UR5 wrists are the same, also their Coulomb and viscous friction coefficients are similar.
The values reported in Figure 3a,b have been calculated for the robot arm only. In the case a tool is mounted, the identified dynamic parameters must be updated to allow the model to correctly estimate the joints torques. To do so, it would be only necessary to execute again the excitation trajectory and apply the LS method reported in Equation (9). A tool, in fact, can be considered as an extension of the robot last link and it does not add any additional degree of freedom to the manipulator, so the regression matrix does not change. On the other hand, in the case of a cumbersome tool, it could be necessary to find another persistent trajectory for the identification process since some changes to the robot constraints q m a x could be required to avoid possible collisions.

5. Validation of the Identified Base Dynamic Parameters

The validation of the proposed algorithm has been done by commanding the two robots with several persistent trajectories, like the ones reported in Figure 4a,b. These have been obtained by substituting pseudo-random values to the FFS coefficients ai,l and bi,l in Equation (5) so that the boundary conditions of Equations (6) and (7) are respected. This approach allows defining persistent trajectories able to excite all the identified base dynamic parameters reported in Appendix A. The corresponding joints torques of the UR3 and the UR5 are reported in Figure 5 and Figure 6, respectively.
To evaluate the accuracy of the identified base dynamic parameters, in Figure 5 and Figure 6 the measured torques (UR3 and UR5, respectively) are compared with the ones estimated with the proposed mathematical model once that the base dynamic parameters p B have been identified (Model after identification). In addition, the joints torques (Model before identification), obtained using only the dynamic parameters provided by the manufacturer [1], are also reported.
This analysis highlights how a proper identification of the robot dynamic parameters is crucial for a correct estimate of the joints torques of a manipulator. As in [26], the normalized error, defined as e N = 1 N _ v a l e T ·   e , is adopted to better quantify the mismatch among measured and modeled joint torques. In the present study, the error vector e has been calculated in three different ways:
  • e = τ ¯ U R τ ¯ B I , where τ ¯ U R = ( τ U R , 1 τ U R , N _ v a l ) and τ ¯ B I = ( τ B I , 1 τ B I , N _ v a l ) are the measured torques for each point k of the validation trajectory and the ones obtained using only the dynamic parameters found in [1];
  • e = τ ¯ U R Y ¯ B   ·   p B , with Y ¯ B = ( Y B , 1 Y B , N _ v a l )   being the observation matrix calculated for the validation trajectory;
  • e = τ ¯ U R Y ¯ B   ·   p B w e i g h t e d , where, instead of the base dynamic parameters p B , the ones obtained through the WLS method have been used.
For the validation trajectories depicted in Figure 4a,b, the normalized errors are reported in Table 3.
By both comparing the trends of the measured and the estimated torques of Figure 5 and Figure 6 and the values of the normalized errors of Table 3, the validation process has been considered to be successful.

6. Temperature Effect on the Base Dynamic Parameters

For a proper modeling and control strategy of a manipulator, the effect of temperature on the robot dynamics should be taken into account. To better show its impact on the base dynamic parameters, the identification trajectory of Figure 2a has been commanded to the UR5 just after its power-up (cold robot) and after half an hour warm-up (warm robot). Joints temperatures have been measured by built-in robot sensors and acquired at 125 Hz via TCP/IP communication using port 30003. Their mean values are hereby reported:
  • Cold robot: T = [25.2, 24.9, 22.2, 26.1, 22.4, 23.2] °C;
  • Warm robot: T = [36.5, 34.6, 32.6, 37.9, 39.7, 40.1] °C.
The identified base dynamic parameters are depicted in Figure 7, where additional frames have been inserted for better visualization of the small values which cannot be read at full scale.
The values ( p B   ( 1 ) p B   ( 36 ) ) and ( p B   ( 49 ) p B   ( 52 ) ), which should not depend on temperature, are nearly the same in both working conditions, proving the stability of the algorithm. On the other hand, there are differences in the friction parameters ( p B   ( 37 ) p B   ( 48 ) ), with a reduction of up to 31% of the viscous ones. Similar percentual variations have been also observed in [19,49].
Since viscous friction depends on joints temperature, a more detailed analysis to better quantify such a correlation has been carried out with both the UR3 and the UR5. To do so, the identification process of the base dynamic parameters has been sequentially repeated multiple times for both manipulators by commanding the same persistent trajectory of Figure 2a. The results are reported in Figure 8a,b.
Since an increment of the joint temperature reduces the viscosity of the grease used to lubricate the joints, the values of the viscous friction coefficients linearly decrease with the rise of temperature in the considered range. The overall percentage reductions of the viscous friction coefficients, between the initial and the final robot operating conditions, are reported in Table 4.
This analysis also highlights the different operating conditions of the two robots even if commanded with the same persistent trajectory. By considering the temperature ranges in Figure 8a,b, it can be seen how, because of their different heat dissipation properties, smaller joints tend to reach higher temperatures than the bigger ones. So, since grease viscosity and joint temperature are inversely proportional in the considered temperature range, this also explains why the viscous friction coefficients of the UR3 face, on average, higher percentage reductions than the ones of the UR5.
Moreover, by comparing the curves related to the UR3 elbow and the three UR5 wrist joints, which are the same, it is possible to observe how the four trends are comparable although different. The slight variations among the three wrist joints of the UR5 could be caused by their possible different lubrication conditions. On the other hand, since the UR3 has been operative for longer than the UR5, their mismatch with the UR3 elbow joint may be related to wear [50]. However, since the manufacturer does not provide any information related to the single joint usage, this assumption, although plausible, cannot be confirmed.
Moreover, to better highlight the impact of changes in grease viscosity on the joints torques necessary to execute the persistent trajectory of Figure 2a, their average and maximum differences, between the initial and the final operating conditions of both manipulators, have been reported in Table 5 and Table 6.
Even though the average variations of Table 5 could be neglected, it should be pointed out that viscous friction torques directly depend on joints angular velocities, so the values reported in this study would be different for other trajectories. This dependency is particularly relevant for industrial robots which usually reach high speeds to maximize productivity. However, by analyzing the values reported in Table 6, it becomes clear how temperature variations affect the dynamic behavior of the tested robots. Such a dependency of joints friction on temperature is also highlighted when comparing the values of Table 6 with the maximum torques, listed in Table 2, that the single joints can provide. The weighted percentage variations, reported as an average among all the joints of the same size belonging to the same robot, are reported in Table 7.
Future studies will be devoted to mapping the dependency of the robot viscous friction coefficients with temperature in the entire working range 0–50 °C defined by Universal Robots.

7. Mounting Configuration Effect on Joints Torques Estimate

To prove the adaptability and reliability of the proposed mathematical model, URSim has been adopted to estimate the torques of the digital twin of a UR5 mounted on the ceiling using the same validation trajectory depicted in Figure 4b. Despite the case in which a tool is added, since the system under analysis is only composed of the UR5, the base dynamic parameters are the same as the ones identified when the manipulator is positioned on the ground. On the contrary, the regressor Y is a function of the robot mounting configuration, so it must be recalculated. If not, there would be a mismatch among the estimated and the measured joints torques since URSim and the proposed mathematical model would describe two different physical systems. As reported in Figure 9, by updating the regression matrix, the estimated torques (Model ceiling) well approximate the ones measured from URSim (UR5 ceiling). On the contrary, if the same regressor obtained for the robot mounted on the ground would have been used, the predicted torques (Model ground) would have not correctly describe the dynamic behavior of the system under analysis, as also proved by the high value of the normalized error e N _ g r o u n d V S c e i l i n g reported in Table 7.
This can be also seen by calculating the joints torques errors, and the corresponding normalized errors reported in Table 8, of the simulated robot in three different scenarios:
  • e g r o u n d = τ ¯ g r o u n d Y ¯ B _ g r o u n d   ·   p B _ S i m , where τ ¯ g r o u n d are the torques read from URSim for each point of the validation trajectory in the case the robot is mounted on the ground, while Y ¯ B _ g r o u n d   ·   p B _ S i m are the predicted ones obtained using the base dynamic parameters identified through Equation (9) using the data from URSim;
  • e c e i l i n g = τ ¯ c e i l i n g Y ¯ B _ c e i l i n g   ·   p B _ S i m , where the acquired ( τ ¯ c e i l i n g ) and the predicted torques ( Y ¯ B _ c e i l i n g   ·   p B _ S i m ) have been calculated as in the previous point, but in the case of the robot mounted on the ceiling;
  • e g r o u n d V S c e i l i n g = τ ¯ c e i l i n g Y ¯ B _ g r o u n d   ·   p B _ S i m , where the two mounting configurations are compared to emphasize the high errors that would have been made if the robot mounting configuration would have not been taken into account.
By comparing the normalized error for the UR5 obtained after the identification of its base dynamic parameters using the LS method (Table 3) with the one in Table 8 for the simulated robot in the same operating conditions ( e N _ g r o u n d ), it can be seen that the first one is higher than the second one. This is related to the fact that there is no noise in the signals coming from the simulated robot, so the negative effect of the regressor condition number is limited. So, the main advantage of testing both the identification and the validation processes with URSim is that the quality of the manipulator base dynamic parameters and, as a consequence, the one of the predicted torques, is higher than the one obtained using the real robot. The stability of the identification process and the ability of the model to effectively predict joints torques in different mounting configurations are also confirmed by the low value of e N _ c e i l i n g , which is comparable with the one obtained when the robot is mounted on the ground.
It is also worth noting that, in the case of the UR5 mounted on the ground, the estimated torques (Model ground) and the ones reported in Figure 6 are not equal. This difference is related to the fact that URSim does not take into account link inertia and joints friction. This proves that the official offline simulator developed by the robot manufacturer cannot be effectively used for an accurate estimate of the dynamic behavior of the manipulator.

8. Conclusions

The present paper highlights how changes in the robot operating conditions affect its dynamic behavior. To do so, a mathematical model able to identify the dynamic parameters of an industrial manipulator is presented and validated with the collaborative robots UR3 and UR5. The two manipulators have been both commanded with persistent exciting trajectories based on 5th order finite Fourier series optimized using a genetic algorithm. Based on the data sampled during the experimental campaign, after the identification of the robots’ base dynamic parameters, the normalized error among the measured and the simulated joints torques is reduced by about the 90% with respect to the one obtained with only the manipulators dynamic properties provided by Universal Robots.
The identification process has been then used to investigate and quantify the effects of temperature on joints viscous friction parameters, allowing a more reliable estimate of the joints torques of the real robot in different working conditions. In the considered temperature range, experimental data showed a reduction of the viscous friction coefficients of about 20% and 17% for the UR3 and the UR5, respectively. This leads to a variation of the joints torques of up to about 8% of the maximum ones transmissible from the robots, showing how this dependency cannot be neglected in the definition of a high-fidelity model of an industrial manipulator.
Due to the possibility to fully identify both Coulomb and viscous friction coefficients and their strict relation with wear, future tests will be run to evaluate the implementation of the proposed algorithm to estimate the health status of the manipulator joints.
In addition, the proposed model has been developed to take into account different mounting configurations which affect the number of base dynamic parameters necessary to fully describe the robot dynamic behavior. The adaptability of the identification algorithm to different setups has been validated using the official simulator developed by the same manufacturer of the tested robots.
The next steps will be focused on testing the proposed algorithm on other robots mounted in different configurations. Moreover, since this work is part of a larger research project whose goal is the definition of a high-fidelity model of an industrial manipulator, more accurate and physically reliable dynamic parameters are required. To do so, the least square method used for the identification process will be replaced with a global optimization algorithm. This would allow building a digital twin of a specific machine in which diagnostics and prognostics algorithms will be implemented to study the fault to failure mechanisms in an industrial robot and their effects on its behavior.

Author Contributions

Conceptualization, A.R. and S.P.; methodology, A.R., R.G. and A.D.M.; software, A.R. and R.G.; validation, A.R. and R.G.; investigation, A.R., R.G. and A.D.M.; writing—original draft preparation, A.R. and R.G.; writing—review and editing, A.R., R.G., A.D.M., S.P., S.M. and M.S.; supervision, A.D.M., S.P., S.M. and M.S.; project administration, S.M. and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In Equations (A1) and (A2) are reported the equations associated with the base dynamic parameters of the UR3 ( p B , U R 3 )   and the UR5 ( p B , U R 5 ) for a general mounting configuration. However, if the axis of the robot base joint is aligned with the gravitational acceleration vector, it is necessary to remove the first two elements of the vectors since only 52 out of 54 base dynamic parameters are involved in the robot dynamics.
p B , U R 3 = { m 1 x 1 m 1 z 1 2.353 J 2 , x z J 1 , y y + J 2 , y y + J 3 , y y J 2 , z z J 3 , z z + J 4 , z z + J m 1 J m 2 + 0.01262 m 4 + 0.01262 m 5 + 0.01262 m 6 + 0.2247 m 4 y 4 21.99 J 3 , z z 16.84 J 2 , z z 16.84 J m 2 + m 2 m 2 x 2 4.104 J 2 , z z 4.104 J m 2 m 2 y 2 4.104 J 2 , x z 4.689 J 3 , x z + m 2 z 2 J 2 , x x J 2 , y y + J 2 , z z + J m 2 J 2 , x y J 2 , y z m 3 21.99 J 3 , z z + m 4 + m 5 + m 6 m 3 x 3 4.689 J 3 , z z m 3 y 3 4.689 J 3 , x z + 0.1123 m 4 + 0.1123 m 5 + 0.1123 m 6 + m 4 y 4 + m 3 z 3 J 3 , x x J 3 , y y + J 3 , z z J 3 , x y J 3 , y z m 4 x 4 0.08535 m 5 + 0.08535 m 6 m 5 y 5 + m 4 z 4 J 4 , x x J 4 , z z + J 5 , z z + 0.007285 m 5 + 0.007285 m 6 0.1707 m 5 y 5 J 4 , x y J 4 , x z J 4 , y y + J 5 , z z + 0.007285 m 5 + 0.007285 m 6 0.1707 m 5 y 5 J 4 , y z m 5 x 5 0.0819 m 6 + m 5 z 5 + m 6 z 6 J 5 , x x + J 6 , y y J 5 , z z + 0.006708 m 6 + 0.1638 m 6 z 6 J 5 , x y J 5 , x z J 5 , y y + J 6 , y y + 0.006708 m 6 + 0.1638 m 6 z 6 J 5 , y z m 6 x 6 m 6 y 6 J 6 , x x J 6 , y y J 6 , x y J 6 , x z J 6 , y z J 6 , z z f c , 1 f c , 6 f v , 1 f v , 6 G 2 · I m 3 G 2 · I m 4 G 2 · I m 5 G 2 · I m 6 }
p B , U R 5 = { m 1 x 1 m 1 z 1 2.353 J 2 , x z J 1 , y y + J 2 , y y + J 3 , y y J 2 , z z J 3 , z z + J 4 , z z + J m 1 J m 2 + 0.01191 m 4 + 0.01191 m 5 + 0.01191 m 6 + 0.2183 m 4 y 4 6.499 J 3 , z z 5.536 J 2 , z z 5.536 J m 2 + m 2 m 2 x 2 2.353 J 2 , z z 2.353 J m 2 m 2 y 2 2.353 J 2 , x z 2.549 J 3 , x z + m 2 z 2 J 2 , x x J 2 , y y + J 2 , z z + J m 2 J 2 , x y J 2 , y z m 3 6.499 J 3 , z z + m 4 + m 5 + m 6 m 3 x 3 2.549 J 3 , z z m 3 y 3 2.549 J 3 , x z + 0.1092 m 4 + 0.1092 m 5 + 0.1092 m 6 + m 4 y 4 + m 3 z 3 J 3 , x x J 3 , y y + J 3 , z z J 3 , x y J 3 , y z m 4 x 4 0.09465 m 5 + 0.09465 m 6 m 5 y 5 + m 4 z 4 J 4 , x x J 4 , z z + J 5 , z z + 0.008959 m 5 + 0.008959 m 6 0.1893 m 5 y 5 J 4 , x y J 4 , x z J 4 , y y + J 5 , z z + 0.008959 m 5 + 0.008959 m 6 0.1893 m 5 y 5 J 4 , y z m 5 x 5 0.0823 m 6 + m 5 z 5 + m 6 z 6 J 5 , x x + J 6 , y y J 5 , z z + 0.006773 m 6 + 0.1646 m 6 z 6 J 5 , x y J 5 , x z J 5 , y y + J 6 , y y + 0.006773 m 6 + 0.1646 m 6 z 6 J 5 , y z m 6 x 6 m 6 y 6 J 6 , x x J 6 , y y J 6 , x y J 6 , x z J 6 , y z J 6 , z z f c , 1 f c , 6 f v , 1 f v , 6 G 2 · I m 3 G 2 · I m 4 G 2 · I m 5 G 2 · I m 6 }

References

  1. Universal Robots—DH. Parameters for Calculations of Kinematics and Dynamics. Available online: https://www.universal-robots.com/articles/ur/application-installation/dh-parameters-for-calculations-of-kinematics-and-dynamics/ (accessed on 3 April 2021).
  2. Kostić, D.; de Jager, B.; Steinbuch, M.; Hensen, R. Modeling and identification for high-performance robot control: An RRR-robotic arm case study. IEEE Trans. Control Syst. Technol. 2004, 12, 904–919. [Google Scholar] [CrossRef]
  3. Bona, B.; Indri, M.; Smaldone, N. Rapid prototyping of a model-based control with friction compensation for a direct-drive robot. IEEE/ASME Trans. Mechatron. 2006, 11, 576–584. [Google Scholar] [CrossRef]
  4. Shao, X.; Wang, S.; Yang, L.; Liu, N. Research on Feedforward Control Based on Robot Dynamics Parameters Identification. In Proceedings of the 2020 IEEE International Conference on Mechatronics and Automation, ICMA 2020, Beijing, China, 13–16 October 2020; pp. 973–978. [Google Scholar]
  5. Sands, T. Development of deterministic artificial intelligence for unmanned underwater vehicles (UUV). J. Mar. Sci. Eng. 2020, 8, 578. [Google Scholar] [CrossRef]
  6. Vidussi, F.; Boscariol, P.; Scalera, L.; Gasparetto, A. Local and trajectory-based indexes for task-related energetic performance optimization of robotic manipulators. J. Mech. Robot. 2021, 13, 1–12. [Google Scholar] [CrossRef]
  7. Gadaleta, M.; Pellicciari, M.; Berselli, G. Optimization of the energy consumption of industrial robots for automatic code generation. Robot. Comput. Integr. Manuf. 2019, 57, 452–464. [Google Scholar] [CrossRef]
  8. Paryanto; Brossog, M.; Bornschlegl, M.; Franke, J. Reducing the energy consumption of industrial robots in manufacturing systems. Int. J. Adv. Manuf. Technol. 2015, 78, 1315–1328. [Google Scholar] [CrossRef]
  9. Sciavicco, L.; Siciliano, B. Robotics: Modeling, Planning and Control; Springer Science & Business Media: London, UK, 2012; ISBN 9781852332211. [Google Scholar]
  10. Vachtsevanos, G.; Lewis, F.; Roemer, M.; Hess, A.; Wu, B. Fault Diagnosis and Prognosis for Engineering Systems. Intell. Fault Diagn. Progn. Eng. Syst. 2006. [Google Scholar] [CrossRef]
  11. Majid, M.A.A.; Fudzin, F. Study on robots failures in automotive painting line. ARPN J. Eng. Appl. Sci. 2017, 12, 62–67. [Google Scholar]
  12. Grosso, L.A.; De Martin, A.; Jacazio, G.; Sorli, M. Development of data-driven PHM solutions for robot hemming in automotive production lines. Int. J. Progn. Health Manag. 2020, 11, 1–13. [Google Scholar]
  13. Melchiorre, M.; Scimmi, L.S.; Mauro, S.; Pastorelli, S.P. Vision-based control architecture for human–robot hand-over applications. Asian J. Control 2020, 23, 105–117. [Google Scholar] [CrossRef]
  14. Scimmi, L.S.; Melchiorre, M.; Mauro, S.; Pastorelli, S.P. Implementing a vision-based collision avoidance algorithm on a UR3 Robot. In Proceedings of the 23rd International Conference on Mechatronics Technology ICMT 2019, Salerno, Italy, 23–26 October 2019. [Google Scholar] [CrossRef]
  15. Melchiorre, M.; Scimmi, L.S.; Pastorelli, S.P.; Mauro, S. Collison Avoidance using Point Cloud Data Fusion from Multiple Depth Sensors: A Practical Approach. In Proceedings of the 23rd International Conference on Mechatronics Technology ICMT 2019, Salerno, Italy, 23–26 October 2019. [Google Scholar] [CrossRef]
  16. Scimmi, L.S.; Melchiorre, M.; Troise, M.; Mauro, S.; Pastorelli, S. A Practical and Effective Layout for a Safe Human-Robot Collaborative Assembly Task. Appl. Sci. 2021, 11, 1763. [Google Scholar] [CrossRef]
  17. Qin, Z.; Baron, L.; Birglen, L. A new approach to the dynamic parameter identification of robotic manipulators. Robotica 2009, 28, 539–547. [Google Scholar] [CrossRef] [Green Version]
  18. Vuong, N.D.; Ang, M.H. Dynamic model identification for industrial robots. Acta Polytech. Hung. 2009, 6, 51–68. [Google Scholar]
  19. Neubauer, M.; Gattringer, H.; Bremer, H. A persistent method for parameter identification of a seven-axes manipulator. Robotica 2015, 33, 1099–1112. [Google Scholar] [CrossRef]
  20. Taghbalout, M.; Antoine, J.F.; Abba, G. Experimental Dynamic Identification of a YuMi Collaborative Robot. IFAC PapersOnLine 2019, 52, 1168–1173. [Google Scholar] [CrossRef]
  21. Gaz, C.; Cognetti, M.; Oliva, A.; Giordano, P.R.; de Luca, A. Dynamic identification of the Franka Emika Panda Robot with retrieval of feasible parameters using penalty-based optimization. IEEE Robot. Autom. Lett. 2019, 4, 4147–4154. [Google Scholar] [CrossRef] [Green Version]
  22. Rackl, W.; Lampariello, R.; Hirzinger, G. Robot excitation trajectories for dynamic parameter estimation using optimized B-splines. In Proceedings of the IEEE International Conference on Robotics and Automation, Saint Paul, MN, USA, 14–18 May 2012; pp. 2042–2047. [Google Scholar] [CrossRef] [Green Version]
  23. An, C.H.; Atkeson, C.G.; Hollerbach, J.M. Estimation of inertial parameters of rigid body links of manipulator. In Proceedings of the 24th IEEE Conference on Decision and Control, Fort Lauderdale, FL, USA, 11–13 December 1985; pp. 990–995. [Google Scholar]
  24. Gautier, M.; Khalil, W. Exciting trajectories for the identification of base inertial parameters of robots. Int. J. Rob. Res. 1992, 11, 362–375. [Google Scholar] [CrossRef]
  25. Antonelli, G.; Caccavale, F.; Chiacchio, P. A systematic procedure for the identification of dynamic parameters of robot manipulators. Robotica 1999, 17, 427–435. [Google Scholar] [CrossRef]
  26. Kovincic, N.; Müller, A.; Gattringer, H.; Weyrer, M.; Schlotzhauer, A.; Brandstötter, M. Dynamic parameter identification of the Universal Robots UR5. In Proceedings of the Austrian Robotics Workshop 2019, Steyr, Austria, 9–10 May 2019. [Google Scholar]
  27. Swevers, J.; Ganseman, C.; Bilgin Tükel, D.; De Schutter, J.; Van Brüssel, H. Optimal robot excitation and identification. IEEE Trans. Robot. Autom. 1997, 13, 730–740. [Google Scholar] [CrossRef] [Green Version]
  28. Wu, W.; Zhu, S.; Wang, X.; Liu, H. Closed-loop dynamic parameter identification of robot manipulators using modified fourier series. Int. J. Adv. Robot. Syst. 2012, 9, 29. [Google Scholar] [CrossRef] [Green Version]
  29. Tika, A.; Ulmen, J.; Bajcinca, N. Dynamic Parameter Estimation Utilizing Optimized Trajectories. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Las Vegas, NV, USA, 25–29 October 2020. [Google Scholar]
  30. Simoni, L.; Beschi, M.; Legnani, G.; Visioli, A. Modelling the temperature in joint friction of industrial manipulators. Robotica 2019, 37, 906–927. [Google Scholar] [CrossRef] [Green Version]
  31. Hao, L.; Pagani, R.; Beschi, M.; Legnani, G. Dynamic and friction parameters of an industrial robot: Identification, comparison and repetitiveness analysis. Robotics 2021, 10, 49. [Google Scholar] [CrossRef]
  32. Nevmerzhitskiy, M.N.; Notkin, B.S.; Vara, A.V.; Zmeu, K.V. Friction Model of Industrial Robot Joint with Temperature Correction by Example of KUKA KR10. J. Robot. 2019, 2019, 6931563. [Google Scholar] [CrossRef]
  33. Remote Control Via TCP/IP 16496. Available online: https://www.universal-robots.com/articles/ur/interface-communication/remote-control-via-tcpip/ (accessed on 3 April 2021).
  34. Gabiccini, M.; Bracci, A.; De Carli, D.; Fredianelli, M.; Bicchi, A. Explicit Lagrangian formulation of the dynamic regressors for serial manipulators. In Proceedings of the XIX Aimeta Congress, Ancona, Italy, 14–17 September 2009. [Google Scholar]
  35. Gautier, M. Numerical calculation of the base inertial parameters of robots. J. Robot. Syst. 1991, 8, 485–506. [Google Scholar] [CrossRef]
  36. Brunton, S.L.; Kutz, J.N. Data-Driven Science and Engineering: Machine Learning, Dynamical Systems, and Control; Press, C.U: Cambridge, UK, 2019; Volume 60, ISBN 9781108422093. [Google Scholar]
  37. Armstrong, B. On Finding Exciting Trajectories for Identification Experiments Involving Systems with Nonlinear Dynamics. Int. J. Rob. Res. 1989, 8, 28–48. [Google Scholar] [CrossRef]
  38. Ding, L.; Wu, H.; Yao, Y.; Yang, Y. Dynamic Model Identification for 6-DOF Industrial Robots. J. Robot. 2015, 2015, 471478. [Google Scholar] [CrossRef] [Green Version]
  39. Jin, J.; Gans, N. Parameter identification for industrial robots with a fast and robust trajectory design approach. Robot. Comput. Integr. Manuf. 2015, 31, 21–29. [Google Scholar] [CrossRef]
  40. Guo, X.; Zhang, L.; Han, K. Dynamic parameter identification of robot manipulators based on the optimal excitation trajectory. In Proceedings of the 2018 IEEE International Conference on Mechatronics and Automation, Changchun, China, 5–8 August 2018; pp. 2145–2150. [Google Scholar]
  41. Calefiore, G.; Indri, M.; Bona, B. Robot dynamic calibration: Optimal excitation trajectories and experimental parameter estimation. J. Robot. Syst. 2001, 18, 55–68. [Google Scholar] [CrossRef]
  42. Jia, J.; Zhang, M.; Zang, X.; Zhang, H.; Zhao, J. Dynamic parameter identification for a manipulator with joint torque sensors based on an improved experimental design. Sensors 2019, 19, 2248. [Google Scholar] [CrossRef] [Green Version]
  43. Universal Robots. The URScript Programming Language; Universal Robots: Odense, Denmark, 2019. [Google Scholar]
  44. Benimeli, F.; Mata, V.; Valero, F. A comparison between direct and indirect dynamic parameter identification methods in industrial robots. Robotica 2006, 24, 579–590. [Google Scholar] [CrossRef]
  45. Universal Robots. Max. Joint Torques. Available online: https://www.universal-robots.com/articles/ur/robot-care-maintenance/max-joint-torques/ (accessed on 3 April 2021).
  46. Bahloul, A.; Tliba, S.; Chitour, Y. Dynamic Parameters Identification of an Industrial Robot with and Without Payload. IFAC-PapersOnLine 2018, 51, 443–448. [Google Scholar] [CrossRef]
  47. Yoshida, K.; Khalil, W. Verification of the positive definiteness of the inertial matrix of manipulators using base inertial parameters. Int. J. Rob. Res. 2000, 19, 498–510. [Google Scholar] [CrossRef]
  48. Leboutet, Q.; Roux, J.; Janot, A.; Guadarrama-Olvera, J.R.; Cheng, G. Inertial Parameter Identification in Robotics: A Survey. Appl. Sci. 2021, 11, 4303. [Google Scholar] [CrossRef]
  49. Simoni, L.; Beschi, M.; Legnani, G.; Visioli, A. Friction modeling with temperature effects for industrial robot manipulators. In Proceedings of the IEEE International Conference on Intelligent Robots and Systems, Hamburg, Germany, 28 September–2 October 2015; Volume 2015, pp. 3524–3529. [Google Scholar] [CrossRef]
  50. Bittencourt, A.C.; Axelsson, P. Modeling and experiment design for identification of wear in a robot joint under load and temperature uncertainties based on friction data. IEEE/ASME Trans. Mechatron. 2014, 19, 1694–1706. [Google Scholar] [CrossRef]
Figure 1. (a) Singular values of the Y ¯ 25 matrix as a function of the number of the dynamic parameters p ; (b) Cumulative energy of the Y ¯ 25 matrix as a function of the number of the dynamic parameters p .
Figure 1. (a) Singular values of the Y ¯ 25 matrix as a function of the number of the dynamic parameters p ; (b) Cumulative energy of the Y ¯ 25 matrix as a function of the number of the dynamic parameters p .
Robotics 10 00083 g001
Figure 2. (a) Excitation trajectory for a UR5 collaborative robot mounted on the ground; (b) Excitation trajectory for a UR5 collaborative robot mounted on the ceiling.
Figure 2. (a) Excitation trajectory for a UR5 collaborative robot mounted on the ground; (b) Excitation trajectory for a UR5 collaborative robot mounted on the ceiling.
Robotics 10 00083 g002
Figure 3. (a) Base dynamic parameters p B of the UR3 collaborative robot; (b) Base dynamic parameters p B of the UR5 collaborative robot.
Figure 3. (a) Base dynamic parameters p B of the UR3 collaborative robot; (b) Base dynamic parameters p B of the UR5 collaborative robot.
Robotics 10 00083 g003
Figure 4. (a) Validation trajectory for the UR3 collaborative robot; (b) Validation trajectory for the UR5 collaborative robot.
Figure 4. (a) Validation trajectory for the UR3 collaborative robot; (b) Validation trajectory for the UR5 collaborative robot.
Robotics 10 00083 g004
Figure 5. Measured (UR3) and estimated (Model) torques of the UR3 collaborative robot for the validation trajectory.
Figure 5. Measured (UR3) and estimated (Model) torques of the UR3 collaborative robot for the validation trajectory.
Robotics 10 00083 g005
Figure 6. Measured (UR5) and estimated (Model) torques of the UR5 collaborative robot for the validation trajectory.
Figure 6. Measured (UR5) and estimated (Model) torques of the UR5 collaborative robot for the validation trajectory.
Robotics 10 00083 g006
Figure 7. Base dynamic parameters of the UR5 collaborative robot just after power-up (cold robot) and after half an hour of warm-up (warm robot).
Figure 7. Base dynamic parameters of the UR5 collaborative robot just after power-up (cold robot) and after half an hour of warm-up (warm robot).
Robotics 10 00083 g007
Figure 8. (a) Trends of the viscous friction coefficients according to temperature for the UR3 collaborative robot; (b) Trends of the viscous friction coefficients according to temperature for the UR5 collaborative robot.
Figure 8. (a) Trends of the viscous friction coefficients according to temperature for the UR3 collaborative robot; (b) Trends of the viscous friction coefficients according to temperature for the UR5 collaborative robot.
Robotics 10 00083 g008
Figure 9. Measured (UR5 ceiling) and estimated (Model ceiling) torques of the simulated UR5 mounted on the ceiling and estimated torques of the same robot positioned on the ground (Model ground) during the validation trajectory.
Figure 9. Measured (UR5 ceiling) and estimated (Model ceiling) torques of the simulated UR5 mounted on the ceiling and estimated torques of the same robot positioned on the ground (Model ground) during the validation trajectory.
Robotics 10 00083 g009
Table 1. Torque constant for each joint of the UR3 and the UR5 used in the experimental campaign.
Table 1. Torque constant for each joint of the UR3 and the UR5 used in the experimental campaign.
RobotKt1 [Nm/A]Kt2 [Nm/A]Kt3 [Nm/A]Kt4 [Nm/A]Kt5 [Nm/A]Kt6 [Nm/A]
UR30.12490.12490.09340.04530.04530.0453
UR50.13500.13610.13550.09570.08650.0893
Table 2. Maximum joints torques of the UR3 and the UR5 collaborative robots [45].
Table 2. Maximum joints torques of the UR3 and the UR5 collaborative robots [45].
JointJoint Size Maximum   Torque   ( τ j , m a x )   [ Nm ]
UR3 base and shoulderSize 256
UR3 elbowSize 128
UR3 wrist 1, wrist 2, and wrist 3Size 012
UR5 base, shoulder, and elbowSize 3150
UR5 wrist 1, wrist 2, and wrist 3Size 128
Table 3. UR3 and UR5 normalized errors before and after the identification of their base dynamic parameters.
Table 3. UR3 and UR5 normalized errors before and after the identification of their base dynamic parameters.
Robot e N   before   Identification e N   after   Identification   ( LS ) e N   after   Identification   ( WLS )
UR30.1945 Nm0.0247 Nm0.0258 Nm
UR50.4207 Nm0.0530 Nm0.0553 Nm
Table 4. Percentage variations of viscous friction coefficients due to the increment of temperature between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
Table 4. Percentage variations of viscous friction coefficients due to the increment of temperature between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
Robot Δ f v 1 Δ f v 2 Δ f v 3 Δ f v 4 Δ f v 5 Δ f v 6
UR3−10.6%−11.6%−19.5%−26.4%−25.5%−23.9%
UR5−8.8%−21.0%−13.6%−23.5%−21.1%−15.5%
Table 5. Average variations of the joints torques between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
Table 5. Average variations of the joints torques between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
Robot Δ C 1 , a v g Δ C 2 , a v g Δ C 3 , a v g Δ C 4 , a v g Δ C 5 , a v g Δ C 6 , a v g
UR30.27 Nm0.42 Nm0.59 Nm0.40 Nm0.41 Nm0.39 Nm
UR50.70 Nm2.43 Nm1.45 Nm0.46 Nm0.72 Nm0.40 Nm
Table 6. Maximum variations of the joint torques between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
Table 6. Maximum variations of the joint torques between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
Robot Δ C 1 , m a x Δ C 2 , m a x Δ C 3 , m a x Δ C 4 , m a x Δ C 5 , m a x Δ C 6 , m a x
UR31.02 Nm1.25 Nm1.17 Nm0.87 Nm1.03 Nm0.79 Nm
UR53.09 Nm10.90 Nm4.75 Nm1.49 Nm1.57 Nm0.98 Nm
Table 7. Weighted percentage variations of the maximum joints torques between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
Table 7. Weighted percentage variations of the maximum joints torques between the initial and the final operating conditions of the UR3 and the UR5 collaborative robots.
JointJoint SizeWeighted Percentage Variation
UR3 base and shoulderSize 22.0%
UR3 elbowSize 14.2%
UR3 wrist 1, wrist 2 and wrist 3Size 07.5%
UR5 base, shoulder and elbowSize 34.2%
UR5 wrist 1, wrist 2 and wrist 3Size 14.8%
Table 8. Simulated UR5 normalized errors in the two different mounting configurations analyzed (ground and ceiling).
Table 8. Simulated UR5 normalized errors in the two different mounting configurations analyzed (ground and ceiling).
Robot e N _ g r o u n d e N _ c e i l i n g e N _ g r o u n d V S c e i l i n g
UR5 (URSim) 0.0043 Nm0.0050 Nm0.6144 Nm
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Raviola, A.; Guida, R.; De Martin, A.; Pastorelli, S.; Mauro, S.; Sorli, M. Effects of Temperature and Mounting Configuration on the Dynamic Parameters Identification of Industrial Robots. Robotics 2021, 10, 83. https://doi.org/10.3390/robotics10030083

AMA Style

Raviola A, Guida R, De Martin A, Pastorelli S, Mauro S, Sorli M. Effects of Temperature and Mounting Configuration on the Dynamic Parameters Identification of Industrial Robots. Robotics. 2021; 10(3):83. https://doi.org/10.3390/robotics10030083

Chicago/Turabian Style

Raviola, Andrea, Roberto Guida, Andrea De Martin, Stefano Pastorelli, Stefano Mauro, and Massimo Sorli. 2021. "Effects of Temperature and Mounting Configuration on the Dynamic Parameters Identification of Industrial Robots" Robotics 10, no. 3: 83. https://doi.org/10.3390/robotics10030083

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