Abstract

The quality of industrial robots essentially depends on the properties of their kinematic couples. This research has involved conducting an experimental study of the friction torque in a joint of the KUKA KR10 industrial robot and building its model. It has been established that the largest impact on friction in the joint is caused by its axial load and velocity, as well as the temperature of the mechanism, which is generally not homogeneous. It is not possible to measure temperature fields in the joints of a serial industrial robot directly. This study has set forth a method to estimate friction torque taking into account the temperature factor indirectly. For this, we have used the motor temperature available for measuring in combination with special periodical motions, performed by the robot, during which we estimated the actual friction torque in the joint and calculated a temperature correction based on our findings.

1. Introduction

The mechanical rigidity of an industrial robot-manipulator is largely determined by the rigidity of its joints. To increase rigidity and minimize backlash, the joints of industrial robots are designed to have high kinematic precision and high load capacity but increased friction on the downside. The authors of the studies [1, 2] point out a number of factors affecting the friction force in robot joints: joint design, manufacturing and assembly precision, material of the joint’s components, lubricants used, geometry and permanency of contact between the mating surfaces, wear and tear of the joint’s parts, backlash, vibration, etc.

Taking friction forces into account may considerably improve the efficiency of controlling the robot [3, 4] and allows assessing the state of lubricant in the joint [5], the mechanical wear and tear, predicting its service life [1], and implementing force-torque control algorithms using the robot’s axial torque data [6]. The literature contains references to many identification techniques [69] and friction models. Both experimental [10, 11] and theoretical friction force determination methods [2, 12, 13] are available.

The most comprehensive friction models include dependencies based on such variables as angular velocity, load, and temperature [5, 10, 12]. It has been experimentally demonstrated [10, 14] that temperature has a substantial effect on friction torque. However, it cannot always be measured since the joint’s interior is not accessible due to its design [2, 5, 12, 15]. In this case, the joint’s friction area temperature can be obtained either by indirect measurements, or through estimation. In the study [10], the temperature of the lubricant was used as the joint’s own temperature. When estimating friction area temperature with the aid of state observers, motor winding temperature [12] or joint housing temperature [5] can be used. The observers are based on various heat transmission models based in turn on the law of conservation of thermal energy and taking into account the motor’s thermal output [12, 15]. Moreover, in some cases lubricant viscosity can be taken into account, as well as thermal effect of friction on it [5]. Temperature estimation methods based on models may have limited accuracy which deteriorates with time due to various unmeasured factors [4].

This study explores friction in one of the joints of the KUKA KR10 R900 sixx robot-manipulator [16]. Its control system allows measuring the motor temperature only. To account for the thermal factor influence on friction in the robot’s joint, a technique is proposed based on introducing a motor temperature correction, which is regularly calculated using the friction torque registered when special axial motions are performed by the robot.

2. Overview

2.1. Friction Function and Its Components

Axial control torque of the industrial robot can be described bywhere is inertia matrix; is vector of centrifugal and Coriolis forces; , , and are angular positions, velocities, and accelerations in the robot’s axles, respectively; is axial load torque; is friction torque counter-directed to the axial velocity.

Friction is a very complex process and to some degree unique to every mechanism [14]. The studies dedicated to tribology use a wide range of friction function formalizations [2, 17]. One of the fullest descriptions of the friction torque in a robot-manipulator joint is presented in the study [10]. It establishes a dependence of friction torque on load torque , angular velocity , and the joint’s temperature and can be presented as follows:where is friction component dependent on load; is Stribeck friction; is viscous friction; is Coulomb friction. Let us consider these components in greater detail.

Coulomb friction is related to the force of intermolecular interaction of the mating surfaces and is a constant [14].

The way, load has an effect on the friction component , is linked with the pressing force of the mating surfaces, their roughness, the presence of lubricant between them, and the material properties of the parts [14]. Axial load of the robot consists of holding torque (caused by the gravitational force, and dependent on the robot’s kinematics) and load torque (caused by the external forces applied to the robot at the current time ):To describe the dependence of friction on load, both linear [5, 7, 10] and nonlinear functions [2] can be used. Most often, this dependence is described as follows:where is a friction coefficient.

The viscous friction component arises when the joint’s mechanical elements move in viscous fluid. This dependence is often expressed linearly in terms of velocity [6, 9, 15, 17]; in some cases a nonlinear function approximation is used: power [12], trigonometrical [18], and exponential [11] functions. The value and characteristics of viscous friction are substantially dependent on the viscosity of lubricant used [14], which in turn depends on temperature. For this reason, a more general description of viscous friction involves the joint’s both velocity and temperature . In the studies [2, 5, 10], this dependence is formalized as follows:where are coefficients to characterize viscosity; are coefficients to take temperature into account.

The nature of the Stribeck effect is related to the catching of peaks on rough surfaces at the microlevel through a layer of liquid lubricant. The effect depends on the finish of the surfaces, lubricant viscosity, and quantity, and it is registered at low velocities and can manifest itself to a varying degree. The Stribeck effect is identified in the mechanisms explored in the studies [3, 18, 19], where it is described as follows:where – Stribeck’s velocity; and – coefficients to characterize the amplitude and slope of the nonlinearity of Stribeck friction, respectively. The expression (6) reflects the dependence of Stribeck friction on velocity only; however, there are studies which, to describe this component further, have taken into account the influence of the joint’s temperature and load [10]. There are also a number of joint designs in which the Stribeck effect has been either not detected, or neglected by the researchers [2, 4, 5, 12, 13, 20].

2.2. Friction Measurement Methods

The most common solution for estimating friction torque in robot-manipulator joints is separating it, using measurable control torque , where the contribution of the remaining components of the equation (1) is calculated with due regard to the mass-inertial characteristics of the robot’s links, and then eliminated [12, 1922]. If these characteristics are unknown, various additional computational approaches and methodologies can be used to identify them [24, 69, 13].

The study [1] proposes an alternative method that does not require identifying mass-inertial characteristics of the robot’s links. In a nutshell, it consists in paired counter-directed motions, performed by one robot axle at constant velocity , without external load . At constant speed there is no influence of inertia forces. Centrifugal and Coriolis forces are directed radially from considered joint and create torque acting only on subsequent joints and links. Then, according to (1), allow the components of the inertia forces and those of centrifugal and Coriolis forces , while the axial control torque counteracts only the friction forces and the load torque , caused by the gravitational forces. Both components depend on the angle . Notably, the gravitational component of the load torque does not depend on the motion direction , whilst the friction torque component changes sign as the direction reverses:where and are control torques registered at the robot axis, moving in the forward and reverse direction, respectively. Then the friction torque and load torque areThis method and its variations [11] do not require a priori knowledge of the robot’s mass-inertial characteristics, and the identified dependences and have high repeatability [18, 19, 23].

2.3. Experiment Preparation

Figure 1 presents an overall view of the industrial robot-manipulator KUKA KR 10 R900 sixx and its kinematic structure with six degrees of freedom. The robot’s load capacity is 10 kg, with maximum reach of 900 mm.

The friction model identification experiments are performed based on the robot’s second axle. To identify friction torque and holding torque the method (8) is used, involving paired motions in counter side directions, as shown in Figure 2. For load torque , the study uses gravitational torque . To take into account the thermal factor, the study uses motor temperature available from the robot’s control system; additionally, the robot’s joint housing temperature is measured. In the study the friction model is considered as the expression (2). Then the friction torque components dependent on load , velocity , and temperature are identified consecutively.

3. Measurement Results

3.1. Dependence of Friction Torque upon Load Torque

To estimate the friction torque component , which is dependent on load torque in the joint, a series of experiments are carried out where, according to the method (8) and scheme in Figure 2, robot motions can be reviewed at various velocities (Algorithm 1).

Inputs:
velocity = [] // vector of velocity values, deg/s
// control torque,
// joint angular position, deg/s
Outputs:
joint angular velocity, deg/s
// control torque until mowing clockwise,
// control torque until mowing counterclockwise,
// friction torque,
// load torque,
// motor temperature, °C
// CYCLE OF CHANGING VELOCITY VALUES
for i = 1 to length(velocity) step 1
= velocity[i] // assigning the velocity value to the robot joint velocity
Moving to = 180.0 // joint clockwise motion
// assigning until mowing clockwise
Moving to = 0.0 // joint counterclockwise motion
// assigning until mowing counterclockwise
and calculating by eq. (8)
data measuring
end for

The experiment results are presented in Figure 3. Taking into account that each experimental curve is obtained for and assuming that during a single experiment (pair of counter-directed motions) , the expression (2) can be reduced towhere   if  . Accordingly, the effects of the component on friction are linked with the shape of the curves presented in Figure 3. Further, their different, but invariable within the same experiment, vertical shifts are determined by the other components in the expression (2). Assuming that for each curve, it is possible to identify its vertical shift and define the law of the total quadratic dependence of friction torque upon the load in the joint:In Figure 3, the experimental dependences are juxtaposed with their approximation obtained in accordance with (10) and plotted with various shifts . The revealed dependence (10) is expressly nonlinear which is in contrast with its traditionally linear shape (4).

3.2. Viscous Friction

To identify the viscous friction component , let us consider experiments at different velocities at varying temperature . Algorithm 2 presents the pseudocode for experiments motion.

Inputs:
velocity = [;;;] // vector of velocity values, deg/s
// motor temperature, °C
// control torque, Nm
Temp_max = // temperature parameter, °C
// joint angular position, deg/s
Outputs:
   joint angular velocity, deg/s
// control torque until mowing clockwise,
// control torque until mowing counterclockwise,
// friction torque,
// load torque,
// motor temperature, °C
// CYCLE OF CHANGING VELOCITY VALUES IN JOINT HEATING-UP MODE
while < 65°C
// Cycle of changing velocity values
for i = 1 to 4 step 1
= velocity[i] // assigning the velocity value for the robot joint velocity
Moving to = 180.0 // joint clockwise motion
// assigning until mowing clockwise
Moving to = 0.0 // joint counterclockwise motion
// assigning until mowing counterclockwise
and calculating by eq. (8)
data measuring
end for
Temp_max = Temp_max + 1 // changing temperature parameter to °C up
// Heating-up cycle by °C making intensive motions at maximum velocity
while < Temp_max
= // assigning the velocity maximum value for the robot joint velocity
Moving to = 180.0 // joint clockwise motion
Moving to = 0.0 // joint counterclockwise motion
end while
end while
// CYCLE OF CHANGING VELOCITY VALUES IN JOINT COOLING-DOWN MODE
while > 25°C
// Cycle of changing velocity values
for i = 1 to 4 step 1
= velocity[i] // assigning the velocity value for the robot joint velocity
Moving to = 180.0 // joint clockwise motion
// assigning until mowing clockwise
Moving to = 0.0 // joint counterclockwise motion
// assigning until mowing counterclockwise
and calculating by eq. (8)
data measuring
end for
Temp_max = Temp_max – 1 // changing temperature parameter to °C down
wait for ( <= Temp_max) // waiting for natural cooling to °C down
end while

Figure 4 presents time plots of motor temperature and joint’s housing temperature , including two phases: heat-up through intensive motions at maximum velocity and natural cool-down. As temperature changed, in order to accumulate experimental data, motions were regularly performed by the robot at four velocities. The obtained experimental dependences , plotted for (points, corresponding to curves minimum, shown on Figure 3) in relation to motor temperature and joint housing temperature , are presented in Figures 5 and 6, respectively.

In all of the presented experiments, a considerable reduction in friction torque with the rise of temperature could be observed. However, it should be noted that dependences of friction torque on temperatures, registered during heat-up and cool-down, are differ. That difference is true for both motor and housing temperatures, being the hottest and coolest points in the joint’s design, respectively, and any of their linear combination that reflects the temperature in some middle point of the joint’s mechanism (not shown in the plot). The observed ambiguity can be explained by the joint’s complex design where friction concurrently occurs on many mating surfaces that are within a heterogeneous temperature field. Thus, the experiments presented in Figures 5 and 6, on the one hand, demonstrate a significant effect of joint temperature on friction and, on the other hand, the complexity of the task of its objective monitoring.

For further experiments, let us perform identification of the viscous friction component in terms of mean temperature and mean friction torque :where – time, and – friction torques registered in heat-up and cool-down modes, respectively.

Figure 7 presents mean experimental data and their approximation, obtained in accordance with the simplified expression (2):where the components excluded, whose effects according to the terms of the experiment are either nonexistent ( for ) or small, are eliminated (the Stribeck friction component manifests itself only at low velocities, which are not used in the experiment in Figure 7).

The viscous friction component in the study is identified according to the formula:which is obtained based on the expression (5) but differs from the latter by a nonlinear dependence on velocity , identified for the mechanism being studied.

3.3. Stribeck Effect

Consider the component of friction torque (2) caused by the Stribeck effect. This effect occurs at low velocities at the pass between friction at rest and friction of motion. The effect manifests itself in increased friction at the beginning of motion and is characterized by the fact that as velocity rises, total friction does not increase, as it follows, for example, from (12), but decreases.

To identify the Stribeck component of friction torque, let’s turn again to the experiment presented in Figure 3, where successive motions are considered in a wide range of velocities starting from 2.3 deg/s. In the process of conducting the experiment the joint’s mechanism heated up inevitably. The temperature factor has a substantial influence on viscous friction torque; however, the magnitude of this influence was identified and formalized as (13) and hence can be accounted for by this correction:which at a given velocity reflects the change in viscous friction torque , caused by the joint’s temperature changing from to . This allows identifying and eliminating the influence of uncontrolled temperature changes from the experimental data, obtaining an estimate of friction torque dependence upon velocity for any temperature :where is friction torque at actual temperature .

Figure 8 presents two dependences of friction torque on velocity at constant temperature °С and load : experimental and model ones. The experimental dependence is built using the expression (15), where for temperature the motor temperature is used. The model dependence is presented by the expression (13), which was obtained earlier by approximation of experimental data at velocities, ranging from 19.1 to 101.2 deg/s. As can be seen from Figure 8, the model dependence provides an acceptable accuracy of describing the experimental data over the whole range of working velocities, including in the proximity of . Therefore, we can conclude that the experimental dependence at low velocities has no manifestations such as the Stribeck effect. So when describing friction in the mechanism under study, this effect can be neglected, accepting that .

3.4. Mean Friction Torque Function

Combining the components of friction torque included in (2), dependent upon load torque L (10), velocity , and temperature (13), as well as taking into account Coulomb friction component , we get an expression for friction torque in the second joint of the KUKA KR 10 robot:where .

It is noteworthy that the expression (16) with presented parameters is obtained, using mean values (12) under the conditions of established temperature ambiguity, and is a qualitative description of friction torque, in the first place. The following section of this paper discusses the use of (16) for quantitative estimates of friction torque.

3.5. Friction Function Temperature Parameter Adjustment

The effects of temperature on friction in the robot joint are substantial. In the case (see Figures 5 and 6) temperature may be responsible for up to half the amount of friction torque, so its effects must not be neglected. However, the study showed that it is impossible to establish a clear relationship between friction torque and temperature in a single point of the joint’s mechanism (the study considered motor and joint housing temperatures, as well as a set of their linear combinations). This can be explained by the joint’s complex design where friction occurs simultaneously between many paired parts, which operate in contact with each other under different conditions, including temperature. Building a friction model, based on distributed temperature, is difficult due to the complex task of registering thermal fields inside the closed joint’s mechanism.

Let us review a technique that allows a practical estimation of friction torque, using the expression (16). It is implemented by making a regularly updated correction to the registered motor temperature :where is temperature parameter in the expression (16), is time, and is correction calculation time.

To calculate the correction , let us turn to the expression (16) and express the temperature parameter from it:where , , and are friction torque, load torque, and velocity in the joint at the time . Thenwhere is motor temperature, registered at the correction calculation time.

Friction torque and load torque included in (18) can be obtained using the method (8), which involves performing two counter-directed motions at constant velocity. The need for performing special motions may necessitate interrupting the robot’s control program, which is generally undesirable, so the interruption procedure should have reasonable periodicity.

Our study has shown that the more intensive thermal processes in the joint currently are (see Figure 4), the larger a change of the temperature correction turns out to be and the more often its adjustment is required. This observation has been used to develop a heuristic periodicity criterion:where is scale coefficient regulating the periodicity of corrections and is parameter setting the time, upon expiration of which correction is made forcibly.

The temperature correction adjustment mechanism is launched when the condition (20) is true. As a measure of temperature change intensity, we used the time which it takes the motor temperature to change by °C. This choice was made in view of the peculiar features of temperature diagnostics in the control systems of KUKA robots, where registration step is 1°C. Parameter provides extra operational security in established temperature conditions. Algorithm 3 presenting pseudocode for realising described technique.

Inputs:
// joint angular position, deg/s
joint angular velocity, deg/s
// joint angular velocity value for temperature adjustment, deg/s
// motor temperature, °C
// time of robot’s control program executing, s
// time parameter, upon expiration of which temperature correction is
made forcibly, s
// identified parameters of friction function eq. (16)
// periodicity criterion parameter
Outputs:
// time which it takes the motor temperature to change by °C, s
// adaptation time parameter, s
// motor temperature at time , °C
// friction torque at time ,
// load torque at time ,
// temperature parameter at time , °C
// temperature correction at time , °C
// friction torque,
// ROBOT’S CONTROL PROGRAM CYCLE
while robot’s control program executing
// Temperature change monitoring cycle
if >= 1 // if motor temperature changed by °C
// time for temperature change
end if
// Temperature parameter adjustment cycle
if 0 or // compliance of heuristic periodicity criterion eq. (20)
Pause of robot’s control program executing
// rewriting temporary variable
// ascending time value to adaptation time parameter
// assigning the velocity value for the robot joint velocity
Moving to = 180.0 // joint clockwise motion
// assigning until mowing clockwise
Moving to = 0.0 // joint counterclockwise motion
// assigning until mowing counterclockwise
and calculating by eq. (8)
by eq. (18) // temperature parameter adjustment
// temperature correction calculation by eq. (19)
end if
// temperature parameter calculation online by eq. (17)
by eq. (16) // friction torque estimation online using with temperature
correction calculated after adjustment cycle
Motions of robot’s control program executing
end while

Let’s demonstrate the proposed friction torque estimation technique, using the experimental data shown earlier in Figure 5. These data present a pairs sequence of axle counter-directed motions, performed at four velocities at different temperatures in heat-up and cool-down modes. Each pair of motions allows, using the method (8), getting a friction torque, which can be used both for calculating the temperature correction and for comparison purposes in analyzing the accuracy of friction torque estimates.

The experiment in Figure 5 showed a substantial difference in the dependences of friction torque upon motor temperature in heat-up and cool-down modes. In Figures 9 and 10 these modes are presented separately and labelled with friction torque estimates. The estimates are obtained using the expression (16) with temperature correction in accordance with (17) at times, determined by the condition (20) for , = 60 minutes. In the presented experiments to calculate temperature corrections a velocity = 67.5 deg/s is used. The labels show moments of time when adjustments were made and corresponding values of the temperature correction .

The statistical study showed that when choosing the velocity for calculating the temperature correction , it is advisable to use the observations below for guidance. Turning to (16) it can be noted that the higher the velocity , the more pronounced the contribution of temperature influence on the friction torque compared to its other components. However, in practice, as velocity grows it is becoming more difficult to stabilize it. And, as a result, the registered torque (1) has a growing component related to acceleration , which according to the method (8) should be absent. Consequently, it is recommended that, in the general case, mid-range velocities should be used for calculating the temperature correction.

The accuracy of friction torque estimates in the proposed technique is linked with how frequently the temperature correction is adjusted . In the experiments in Figures 9 and 10, six adjustments were made in each of them for the values of the parameters and = 60 minutes. Using the same experimental data and changing the values of the parameters and , let us consider the effect of temperature adjustment frequency on the accuracy of friction torque estimates. Let us use standard deviation, taken for the whole range of velocities under study as an accuracy test:where is number of experiments for a given velocity, is number of velocities under study, is friction torque obtained with the use of the method (8) and smoothed out using polynomial approximation (as shown in Figures 9 and 10), and is friction torque estimate obtained using the proposed technique.

The findings of a comparative experiment are presented in Figure 11 in the form of a deviation histogram depending on the number of corrections for heat-up and cool-down modes. It is noteworthy that, in both modes, the joint’s temperature varies quite widely—around °C (see Figure 4). Also, the cooling process, lasting for 300 minutes, runs three times as slow as the heating process, with the adjustment procedure carried out three times as rare, with comparable friction torque estimation accuracy. This confirms the universality of the proposed criterion (20). The presented histogram illustrates that even in the modes of intensive natural heat-up and cool-down of the robot, about seven corrections are required for a temperature variation range of °C. The corresponding standard deviation does not exceed 2 Nm.

On the whole, it should be noted that the correction can be updated in relatively long time periods owing to high inertia of thermal processes and in view of the temperature estimation in (17), which indirectly takes into account its change rate. This periodicity does not affect the robot’s performance to any significant degree and ensures high accuracy of friction torque estimation in its joints in real time without using any auxiliary equipment.

4. Conclusions

The study has explored the friction function in the second joint of the KUKA KR10 R900 sixx robot depending on load, velocity, and temperature. To separate friction torque from an axial control torque signal, an experimental technique, not requiring knowledge of the robot’s mass-inertial characteristics, has been used. Nonlinear nature of the friction torque dependence upon load torque has been confirmed experimentally and its description proposed. The combined influence of velocity, motor winding temperature, and joint housing temperature on friction torque during heat-up and cool-down has been explored. A function that approximates mean experimental data describing the dependence has been proposed. The study has demonstrated the absence of the Stribeck effect in the joint.

The experimental data do not allow building a definite friction torque dependence on the friction area temperature, since the robot’s design does not allow measuring temperature that determines friction in the joint. To solve this problem, a technique to estimate temperature in friction area has been proposed. In a nutshell, the technique consists in indirect estimation of temperature using the axial torque of the robot in the process of its performing special motions and adjusting the friction function temperature correction. An adjustment frequency criterion, independent of temperature mode, has been proposed. The effect of adjustment frequency on the accuracy of torque estimation has been studied. The conducted experiments have shown that the proposed solution allows estimating friction torque with a standard deviation not exceeding 2 Nm. Recommendations for choosing special motion velocities and adjustment periodicity have been given.

The friction function estimation technique is universal and allows developing a friction model depending on the joint’s angular velocity, load, and temperature in any serial industrial robot with a serial kinematic structure. The friction-area temperature estimation technique can be tested in robot operation modes, when substantial changes in the joint temperature occur.

Precise knowledge of axial friction forces will allow estimating the external forces acting on the robot using data on its axial torques. The availability of mass-inertial characteristics and the robot’s kinematic model along with axial friction models will help estimate the external forces, acting on the process tools, as well as implement force-torque control algorithms. Within the latter task some identification methods described in articles [6, 8] potentially can be combined with described nonlinear friction model. Nevertheless, the temperature factor compensation is important in the modes of significant temperature changes.

Data Availability

The plots, graphics, and estimation results data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was financially supported by the Ministry of Education and Science of the Russian Federation [state contract 02.G25.31.0173].