Next Article in Journal
Investigation into the Influence of Parallel Offset Wear on Stirling Engine Piston Rod Oil-Free Lubrication Seal
Next Article in Special Issue
Robot Static Path Planning Method Based on Deterministic Annealing
Previous Article in Journal
An Energy Efficiency Tool Path Optimization Method Using a Discrete Energy Consumption Path Model
Previous Article in Special Issue
USV Application Scenario Expansion Based on Motion Control, Path Following and Velocity Planning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A New Parameter Identification Method for Industrial Robots with Friction

1
Academy for Enigineering and Technology, Fudan University, Shanghai 200433, China
2
Guanghua Lingang Engineering Application and Technology R & D (Shanghai) Co., Ltd., Shanghai 201306, China
3
School of Mechatronic Engineering and Automation, Shanghai University, Shanghai 200444, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(5), 349; https://doi.org/10.3390/machines10050349
Submission received: 11 March 2022 / Revised: 2 May 2022 / Accepted: 5 May 2022 / Published: 9 May 2022
(This article belongs to the Special Issue Advanced Control Theory with Applications in Intelligent Machines)

Abstract

:
Commonly used intelligent algorithms that are used to identify the parameters of friction of industrial robots have poor accuracy or involve complex coding, which is not conducive to their use in engineering. This paper uses the random wandering simulated annealing-based variable-step beetle antennae search (RWSAVSBAS) algorithm to identify the parameters of friction of industrial robots. The moment of friction of the third joint of the robot is experimentally obtained and used to establish a Stribeck friction model. Following this, the RWSAVSBAS algorithm is used to identify the frictional parameters of the industrial robot. These parameters can be used to accurately predict the friction-induced torque of the robot.

1. Introduction

Frictional forces in industrial robots can cause tracking errors as large as 50% to degrade the accuracy of their trajectories of motion [1]. Models of friction of the industrial robot have thus been widely developed to compensate for this [2,3,4]. Currently used empirical models of friction for the joints of robot can be divided into two categories: static and dynamic models. Dynamic models can better reflect the dynamic frictional characteristics of the joints of the robot during its motion but have complex structures and are difficult to measure. The Stribeck friction model is representative of static friction models, which can provide better predictions of the frictional force than dynamic models for the section of the robot moving at low speed. This study thus uses the Stribeck model to examine friction in industrial robots [5].
Chao Chen used the differential evolution (DE) algorithm to identify the parameters of a continuous and static friction model for industrial robots, and the experimental results showed that it can adequately respond to the frictional characteristics of such robots [6]. Zhang Jing used the genetic algorithm (GA) to identify the parameters of the Stribeck friction model for industrial robots, and the experimental results demonstrated that this method can adequately identify the parameters of friction [7]. However, coding for the DE algorithm and GA is complicated and unsuitable for use in engineering [8,9]. Wu et al. proposed an improved Stribeck model and then used the Levenberg–Marquarft (L–M) algorithm to identify the frictional parameters of the industrial robot [10]. However, the solution process of the L–M method is susceptible to the initial solution, and it requires a large memory to solve complex problems [11]. Therefore, an algorithm that is simple to code and is highly accurate is needed to identify the appropriate model of friction for industrial robots.
In our previous work, we proposed a random wandering simulated annealing variable step beetle antennae search algorithm (RWSAVSBAS) algorithm. The results of verification showed that it has a higher accuracy than commonly used affine algorithms [12,13,14]. Because this algorithm was based on the BAS algorithm, which has the advantage of simple coding, we use RWSAVSBAS in this paper to identify the frictional parameters of industrial robots.
The remainder of this paper is structured as follows: The principle of measuring frictional torque in an industrial robot is first introduced in Section 2 and is used to obtain the frictional torque in joint 3 of the robot. Following this, the principle of the RWSAVSBAS algorithm is introduced in Section 3. The performance of the algorithm in terms of identifying the parameters of friction of the robot is compared with that of commonly used improved BAS algorithms in Section 4. The conclusions of this study are summarized in Section 5.

2. Experimental Data Acquisition

2.1. Friction Torque Sampling Principle

To identify the frictional force of an industrial robot, we need to first obtain data on the speed of a certain joint of the robot and the frictional force in it. To this end, the robot joint was made to move at a uniform speed. Its inertial force, Koch force, and centrifugal force can be neglected when they are small [7]. Then, we get
τ = G ( q ) + F ( q ˙ )
When measuring the frictional force in joint i of the industrial robot, the other joints are locked. The joint of the robot was then made to move from angle θ1 to θ2 at a uniform speed. The moment of the joint consisted only of friction and gravity at this time, i.e.,
τ 1 = G ( q ) + F ( q ˙ )
Making the robot joint move from angle θ2 to θ1 with uniform speed, we have:
τ 2 = G ( q ) F ( q ˙ )
The frictional forces in the forward and reverse directions of the robot at the same speed are approximately equal. Then, we have:
F ( q ) = F ( q ˙ )
Then, after eliminating gravity from the above equation, the robot friction moment can be obtained as:
F ( q ˙ ) = 0 . 5 × ( τ 1 τ 2 )
When the robot joint movement direction and gravity are perpendicular, the direction of movement of the robot is horizontal, its gravity is 0, and the measurement of the obtained moment is the friction moment. Then, there is:
τ = F ( q ˙ )

2.2. Experimental Data Acquisition

This study used a six-joint industrial robot for experiments. This was a tandem robot consisting of a base, waist, large arm, small arm, and wrist. The third joint of the robot was “taught” to move at a uniform speed in the joint space, and information on its torque was obtained through a built-in current sensor in the servo driver. A schematic of the sampling site is shown in Figure 1.
The sampling steps are as follows
(1)
Set the trajectory of the robot in the joint space through a demonstrator.
(2)
Move the robot at a uniform speed in the space according to instructions of the trainer while obtaining information on its frictional torque through the servo drive.
(3)
The sampled data on the robot’s joint and frictional torque in it are read by the controller.
(4)
Repeat the above steps by moving the robot in the space at different speeds.
(5)
Record and process the data.
Signals of the torque in the joint and the frictional torque obtained by sampling are shown in Figure 2 and Figure 3, respectively:
Figure 3 shows that the sampled friction force increased in the low-speed section of the robot, as expected.

2.3. Friction Modeling

The Stribeck model can more accurately reflect the non-linear characteristics of friction, and can be expressed as follows:
f ( w ) = τ c + ( τ s τ c ) e ( w w s ) δ + τ v w
where w denotes the speed of rotation of the joint of the robot’s arm, f ( w ) enotes the frictional moment of the arm, τc and τs denote Coulomb friction and static friction, respectively, τc denotes the coefficient of viscous friction, δ represents the empirical parameter, which was set to one here according to Tustin’s empirical model [10,15], and ws denotes the Stribeck model.
Figure 4 shows that when the speed of the robot’s joint was zero, the intersection of the Stribeck model and the vertical coordinate was the static friction τs; when its speed was sufficiently high to form a tangent line along w, the intersection of the tangent line and the intercept of the y-axis was τc, and the slope of the tangent line was τv. The horizontal coordinate of the intersection of the tangent line at zero velocity and the line formed when τc was parallel to the x-axis was ws.

3. Related Work

The BAS is a single-particle optimization algorithm that has the advantage of a high efficiency of iteration but the disadvantage of easily falling into a local optimum for complex function optimization problems. In our previous work, we proposed the RWSAVSBAS algorithm and used it to identify geometric errors in the parameters of medical robots. This method is superior to the commonly used and improved BAS algorithm and the shuffled frog-leaping algorithm. We thus used RWSAVSBAS to identify the frictional parameters of the industrial robot [12].

3.1. RWSAVSBAS Algorithm Principle

Compared with the BAS algorithm, the RWSAVSBAS algorithm has the following main improvements:

3.1.1. “Wandering Antennae”

In the original bettle left and right Antennae base, add “Wandering Antennae”, “Wandering Antennae” iteration formula as follows:
G best new = G best + sin ( 2 π × P / h ) × s t e p
The related behavior of wolves has attracted the interest of many scholars [16,17]. Gbestnew denotes the global optimal solution of BAS after the wandering search of the wolf pack algorithm. The new solution obtained after the search is compared with the solution obtained by BAS to determine the subsequent direction of movement of the beetle. The accuracy of identification of the local optimum by BAS can thus be improved.

3.1.2. Metropolis Guidelines

The original BAS is prone to falling into local optimal solutions during the search process. Inspired by the metropolis criterion, RWSAVSBAS has an inner and outer double loop, and selectively accepts poor solutions through the metropolis criterion during iterations of the inner loop to improve the accuracy of global optimization of BAS [18].

3.1.3. Variable Step Length

The step size step is constant in the original BAS, which has better global search capability when the step size is large but can easily miss the local optimal solution. The RWSAVSBAS algorithm uses a variable step size to update the sizes of the steps of BAS. The initial step size of the RWSAVSBAS algorithm is large as it focuses on a global search, and then decreases gradually in later stages when it performs local search through dynamic adjustments to the step size of BAS. This helps to improve the search capability of the RWSAVSBAS algorithm [19].
The formula for updating the step size in RWSAVSBAS is as follows:
s t e p = s t e p × e t a
where eta denotes the step adjustment factor, and eta takes a value between 0 and 1.

3.2. RWSAVSBAS Algorithm Flow

The flow chart is as follows:
Then, the flow of RWSAVSBAS algorithm can be obtained from Figure 5 as follows [12]:
  • A: Initialize the relevant parameters.
  • B: Obtain the fitness values of the three antennae of the algorithm.
  • C: Search for the optimal solution of the algorithm in the inner loop and select it according to the metropolis criterion Q times.
  • D: Once the iteration of the inner loop is complete, use the outer loop to update the step size of the algorithm as well as the optimal solution.
  • E: If the algorithm does not satisfy the termination condition, go to step B. Otherwise, stop the iterations, and use the obtained parameters to identify frictional parameters of the industrial robot.

4. Experimental Results and Analysis

4.1. Experimental Parameter Setting

We set the number of iterations of all algorithms to 1000. The initial step size of the RWSAVSBAS algorithm was set to 0.1, the temperature for simulated annealing T was set to 10,000 [20], the scale of decay of the annealing coefficient was 0.95, the inner loop Q of the RWSAVSBAS algorithm was 100, the h-value of the “wandering antennae” of the inner loop was 100, the h-value of the “wandering antennae” of the outer loop was 1000, the rate of change in step size eta was 0.98, the number of instances of detection of the step and antennae of the beetle was set to 0.98, the antennae of the inner loop were set to 100, the ratio of the step of the beetle to the distance d0 between the antennae was set to three, and c was set to three. The parameters of the variable-step beetle antennae search (VSBAS) algorithm were the same as those of the RWSAVSBAS algorithm, 0.1 and 0.98, respectively [19]. Values of stepmax and stepmin for the improved beetle antennae search (IBAS) algorithm were set to 0.3 and 0.1, respectively, and the remaining values were the same as those of the RWSAVSBAS algorithm [21].

4.2. Adaptation Function

In this paper, the error between the friction torque predicted by the Stribeck model and the actual measured torque is defined as:
Y i = F i F i
In the above Equation (10), Fi denotes the friction moment value predicted by the Stribeck model, F i denotes the actual sampled friction moment value, i denotes the ith sampling point, and Yi denotes the friction moment error predicted by the Stribeck model for the ith time. Then, the objective function is defined as:
min ( i = 1 N Y i 2 )
In this paper, there are 17 actual sampling points, i.e., N = 17, and the objective function is the sum of the squared errors of all sampling points [8].

4.3. Experimental Results

4.3.1. Iterative Optimization Curve Comparison

Owing to the stochastic nature of intelligent algorithms, we ran each of the VSBAS, IBAS and RWSAVSBAS algorithms eight times to obtain their average curves of adaptation of the iterations eight times, as shown in Figure 6.
The RWSAVSBAS algorithm was the most efficient on average, and converged quickly to close to the optimal solution early on. The VSBAS was the slowest to converge in the early stages. Table 1 shows the average fitness, maximum fitness, and minimum fitness of each algorithm.
Table 1 shows that the average fitness of the RWSAVSBAS algorithm was 0.00013, and those of the VSBAS, IBAS were 0.60 and 0.42, respectively. Thus, the average fitness of the RWSAVSBAS algorithm was three orders of magnitude smaller than the other algorithms, and it was thus more stable. The maximum fitness of VSBAS was 1.86, which was the highest of all algorithms, followed by IBAS with 0.83, and RWSAVSBAS with only 0.00014. The RWSAVSBAS algorithm thus still delivered the best performance in terms of maximum fitness. The RWSAVSBAS algorithm had a value of 0.00013 for the minimum adaptation, which was 2–3 orders of magnitude smaller than the 0.10 recorded by VSBAS, and 0.09 by IBAS.

4.3.2. Frictional Torque Error Comparison

The frictional parameters of the four algorithms in case of minimum adaptation were selected and substituted into the above Stribeck model. Figure 7 shows the predicted errors in the minimum frictional moment for each algorithm.
It is clear from Figure 7 that the errors in the prediction of the IBAS and VSBAS algorithms were larger, with the maximum errors of VSBAS and IBAS at around 0.15 N.m and the overall errors at around 0.08 N.m. The RWSAVSBAS algorithm delivered the best performance in terms of errors in the high-speed and low-speed parts of the moving robot. Error of the frictional torque predicted by it was small and tended toward zero. Therefore, the RWSAVSBAS algorithm recorded the best performance overall.

4.3.3. Comparison of Predicted Friction Moments

A combination of the results in Table 1 and Figure 6 and Figure 7 shows that the RWSAVSBAS algorithm could accurately predict the frictional torque of the robot. We thus chose the parameters corresponding to the minimum fitness value of this algorithm over eight iterations as ones to identify friction. The values of each parameter are shown in Table 2.
The curve of frictional moment predicted by the RWSAVSBAS algorithm, which was obtained by substituting the parameters of the Stribeck friction model from Table 2 and the measured velocity data into Equation (7), is shown by the red line in Figure 8. It is also compared with the empirically acquired curve of the frictional moment.
Figure 8 shows that the difference between the frictional torque predicted by the RWSAVSBAS algorithm and the empirically measured frictional torque was not significant, and the trend of the former was consistent with that of the latter in both the low-speed and the high-speed stages. This was expected. This verifies the accuracy of the frictional parameters identified by the RWSAVSBAS algorithm.

5. Conclusions

(1)
In order to improve the accuracy of the friction parameter identification of industrial robots, this paper firstly establishes the friction moment of the robot by demonstration, and then uses the Stribeck model to establish the friction model of industrial robots. In this paper, a new Beetle Antennae search algorithm, RWSAVSBAS, is used to identify the friction model parameters of industrial robots. The experimental results show that the RWSAVSBAS algorithm maintains the advantage of BAS algorithm with fewer parameters, and at the same time has better recognition accuracy compared with other common improved Beetle Antennae search algorithm algorithms.
(2)
This study did not consider the effects of temperature, load, and other factors when measuring the frictional force. Future work in this area should focus on the effects of these external factors on the frictional force. Moreover, the RWSAVSBAS algorithm was used here only to identify the frictional parameters of industrial robots in series. Subsequent work should consider this task for industrial robots operating in parallel.
(3)
The friction model identified by RWSAVSBAS will be written into the controller, and the torque feedforward is sent down from the software algorithm layer for control. The torque feedforward can effectively improve the following performance of the robot, and the theoretically derived robot dynamics generally ignore the effect of friction, so for the actual robot control, getting an accurate friction model can also improve the robot control performance, which includes dragging and dragging demonstration and collision detection of various human–robot collaboration functions that are inseparable from friction compensation.

Author Contributions

Software, B.K., Y.H., Validation, B.K.; Writing—original draft, B.K.; Methodology, B.K., Y.H., P.W., J.Z. and S.G. funding acquisition, S.G. project administration, S.G. Writing—review & editing, B.K. and D.R. article chart editor, B.K. and D.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Shanghai Science and Technology Program (21511101701).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kermani, M.R.; Wong, M.; Patel, R.V.; Moallem, M.; Ostojic, M. Friction compensation in low and high-reversal-velocity manipulators. In Proceedings of the IEEE International Conference on Robotic and Automation, New Orleans, LA, USA, 26 April–1 May 2004; pp. 4320–4325. [Google Scholar]
  2. Iwatani, M.; Kikuuwe, R. An identification procedure for rate-dependency of friction in robotics joints with limited motion ranges. Mechatronics 2016, 36, 36–44. [Google Scholar] [CrossRef]
  3. Guo, Y.T.; Min, B.; Chen, B. Study on adaptive friction compensation of turntable servo systems based on Lucre friction model. Aerosp. Shanghai 2017, 34, 65–69. [Google Scholar]
  4. Naser, M.F.M. Hysteresis Loop of the Lucre Model; Pergamon Press, Inc.: Oxford, UK, 2015. [Google Scholar]
  5. Liu, L.L.; Liu, H.Z.; Wu, Z.Y.; Wang, Z.M. An overview of friction models in mechanical systems. Adv. Mech. 2008, 38, 201–213. [Google Scholar]
  6. Chen, C. Research on High-Performance Motion Control Method for Serial Robotic Manipulator. Ph.D. Thesis, Shandong University, Jinan, China, 2019. [Google Scholar]
  7. Zhang, J. Research on Friction Compensation Strategy of Industrial Robot Joint. Master’s Thesis, Xiamen University, Xiamen, China, 2019. [Google Scholar]
  8. Zhang, J.; Sanderson, A.C. JADE: Adaptive Differential Evolution With Optional External Archive. IEEE Trans. Evol. Comput. 2009, 13, 945–958. [Google Scholar] [CrossRef]
  9. Ding, L.; Wu, H.; Yao, Y.; Li, Y.; Xie, B.H.; Chen, B. Parameters Identification of Industrial Robots Based on WLS-ABC Algorithm. J. South China Univ. Technol. (Nat. Sci. Ed.) 2016, 44, 90–95. [Google Scholar]
  10. Wu, X.; Liu, D.; He, M.; Gao, F.; Shao, G. Research on friction modeling and compensation of robot manipulator. Chin. J. Sci. Instrum. 2018, 39, 44–50. [Google Scholar]
  11. Cao, F. Application of Improved BP-LM Algorithm in Chinese Population-aging Ratio Forecasting. J. Guangxi Univ. Financ. Econ. 2012, 25, 64–69. [Google Scholar]
  12. Kou, B.; Ren, D.; Guo, S. Geometric Parameter Identification of Medical Robot Based on Improved Beetle Antennae Search Algorithm. Bioengineering 2022, 9, 58. [Google Scholar] [CrossRef] [PubMed]
  13. Jiang, X.Y.; Li, S. BAS: Beetle antennae search algorithm for optimization problems. arXiv 2018, arXiv:1710.10724v1. [Google Scholar] [CrossRef]
  14. Jiang, X.Y.; Li, S. Beetle antennae search without parameter tun-ing (BAS WPT) for multi-objective optimization. arXiv 2017, arXiv:1711.02395v1. [Google Scholar]
  15. Marton, L.; Lantos, B. Modeling, Identification, and Compensation of Stick-Slip Friction. IEEE Trans. Ind. Electron. 2007, 54, 511–521. [Google Scholar] [CrossRef]
  16. Wu, H.S.; Zhang, F.M.; Wu, Q.S. New swarm intelligence algorithm-wolf pack algorithm. Syst. Eng. Electron. 2013, 35, 2430–2438. [Google Scholar]
  17. Haber, R.; Strzelczak, S.; Miljkovic, Z.; Castaño, F.; Fumagalli, L.; Petrović, M. Digital twin-based Optimization on the basis of Grey Wolf Method. A Case Study on Motion Control Systems. In Proceedings of the 2020 IEEE Conference on Industrial Cyberphysical Systems (ICPS), Tampere, Finland, 10–12 June 2020. [Google Scholar]
  18. Kirkpatrick, S. Optimization by simulated annealing: Quantitative studies. J. Stat. Phys. 1984, 34, 975–986. [Google Scholar] [CrossRef]
  19. Zhou, T.; Qian, Q.; Fu, Y. Fusion Simulated Annealing and Adaptive Beetle Antennae Search Algorithm. Commun. Technol. 2019, 52, 1626–1631. [Google Scholar]
  20. Gao, S. Research on annealing strategy in Simulated Annealing Algorithm. Aeronaut. Comput. Technol. 2002, 32, 20–22. [Google Scholar]
  21. Wang, H.B.; Peng, X.; Xia, X.Y.; Zhou, S.Y. Path Planning of Automated Guided Vehicle Based on Improved Beetle Antennae Search. Mach. Des. Res. 2021, 37, 27–30. [Google Scholar]
Figure 1. Schematic diagram of robot sampling site.
Figure 1. Schematic diagram of robot sampling site.
Machines 10 00349 g001
Figure 2. Robot joint torque signal.
Figure 2. Robot joint torque signal.
Machines 10 00349 g002
Figure 3. Friction torque of industrial robot.
Figure 3. Friction torque of industrial robot.
Machines 10 00349 g003
Figure 4. Stribeck model.
Figure 4. Stribeck model.
Machines 10 00349 g004
Figure 5. RWSAVSBAS algorithm flow.
Figure 5. RWSAVSBAS algorithm flow.
Machines 10 00349 g005
Figure 6. Comparison of the average adaptation curves of each algorithm optimized for 8 iterations.
Figure 6. Comparison of the average adaptation curves of each algorithm optimized for 8 iterations.
Machines 10 00349 g006
Figure 7. Comparison of minimum error of friction torque.
Figure 7. Comparison of minimum error of friction torque.
Machines 10 00349 g007
Figure 8. Friction torque prediction curves of RWSAVSBAS algorithm.
Figure 8. Friction torque prediction curves of RWSAVSBAS algorithm.
Machines 10 00349 g008
Table 1. Comparison of the average adaptation degree of each algorithm.
Table 1. Comparison of the average adaptation degree of each algorithm.
AlgorithmsAverage FitnessFitness (Max)Fitness (Min)
VSBAS0.601.860.10
IBAS0.420.830.09
RWSAVSBAS0.000130.000140.00013
Table 2. Stribeck model identification parameter values.
Table 2. Stribeck model identification parameter values.
Parametersτcτswsτv
Identifying value0.0840.6410.0660.0026
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kou, B.; Huang, Y.; Wang, P.; Ren, D.; Zhang, J.; Guo, S. A New Parameter Identification Method for Industrial Robots with Friction. Machines 2022, 10, 349. https://doi.org/10.3390/machines10050349

AMA Style

Kou B, Huang Y, Wang P, Ren D, Zhang J, Guo S. A New Parameter Identification Method for Industrial Robots with Friction. Machines. 2022; 10(5):349. https://doi.org/10.3390/machines10050349

Chicago/Turabian Style

Kou, Bin, Yao Huang, Pengpeng Wang, Dongcheng Ren, Jie Zhang, and Shijie Guo. 2022. "A New Parameter Identification Method for Industrial Robots with Friction" Machines 10, no. 5: 349. https://doi.org/10.3390/machines10050349

APA Style

Kou, B., Huang, Y., Wang, P., Ren, D., Zhang, J., & Guo, S. (2022). A New Parameter Identification Method for Industrial Robots with Friction. Machines, 10(5), 349. https://doi.org/10.3390/machines10050349

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