[go: up one dir, main page]

CN112179551B - Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot - Google Patents

Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot Download PDF

Info

Publication number
CN112179551B
CN112179551B CN202010843773.1A CN202010843773A CN112179551B CN 112179551 B CN112179551 B CN 112179551B CN 202010843773 A CN202010843773 A CN 202010843773A CN 112179551 B CN112179551 B CN 112179551B
Authority
CN
China
Prior art keywords
joint
tested
torque
coefficient
robot
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010843773.1A
Other languages
Chinese (zh)
Other versions
CN112179551A (en
Inventor
吴军
刘子麟
王立平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202010843773.1A priority Critical patent/CN112179551B/en
Publication of CN112179551A publication Critical patent/CN112179551A/en
Application granted granted Critical
Publication of CN112179551B publication Critical patent/CN112179551B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/22Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring the force applied to control members, e.g. control members of vehicles, triggers
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J19/00Accessories fitted to manipulators, e.g. for monitoring, for viewing; Safety devices combined with or specially adapted for use in connection with manipulators
    • B25J19/0095Means or methods for testing manipulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明提出一种机器人的关节电机转矩系数与摩擦力同步测试方法和装置,方法包括:根据机器人的动力学模型和机器人的结构参数计算机器人的驱动关节的驱动力矩解析表达式;将存储的实验数据输入驱动力矩解析表达式得到正向力矩和反向力矩、以及驱动电流进行拟合得到不同实验位置的拟合线性函数的斜率和截距;结合待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将截距取相反数得到待测试关节正向与反向运动时的摩擦力系数;计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。由此,实现关节电机转矩系数与正反向摩擦力系数的同步测试,提高测试精度。

Figure 202010843773

The invention provides a method and device for synchronizing the torque coefficient and friction force of a joint motor of a robot. The method includes: calculating an analytical expression of the driving torque of the driving joint of the robot according to the dynamic model of the robot and the structural parameters of the robot; The experimental data is input into the analytical expression of the driving torque to obtain the forward torque and reverse torque, and the driving current is fitted to obtain the slope and intercept of the fitted linear function at different experimental positions; combined with the ratio of the reduction ratio of the joint to be tested, the test is obtained. The motor torque coefficient when the joint is moving forward and reverse, take the opposite number of the intercept to get the friction coefficient when the joint to be tested is moving forward and reverse; calculate the difference between the motor torque coefficient and friction coefficient at different experimental positions The average value is obtained to obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion. Thereby, the synchronous test of the torque coefficient of the joint motor and the coefficient of forward and reverse friction force is realized, and the test accuracy is improved.

Figure 202010843773

Description

机器人的关节电机转矩系数与摩擦力同步测试方法和装置Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot

技术领域technical field

本发明公开了一种机器人的关节电机转矩系数与摩擦力同步测试方法和装置,尤其涉及工业机器人动力学与控制技术领域。The invention discloses a method and device for synchronizing the torque coefficient and friction force of a joint motor of a robot, and particularly relates to the technical field of industrial robot dynamics and control.

背景技术Background technique

工业机器人是广泛应用在各个工业制造领域的重要设备,其能在危险恶劣的工作环境中实现高效准确的运行,是实现自动化生产的基石。为了提升工业机器人的综合性能,在机器人的性能评价、优化设计、控制算法开发、参数整定等环节中,对机器人“控制-驱动-电机-机械”这一机电系统模型准确度有着较高要求,其中部分无法直接获取或不够准确的关键参数需要通过测试得到准确数值。Industrial robots are important equipment widely used in various industrial manufacturing fields. They can achieve efficient and accurate operation in dangerous and harsh working environments, and are the cornerstone of automated production. In order to improve the comprehensive performance of industrial robots, there are high requirements for the accuracy of the electromechanical system model of the robot "control-drive-motor-mechanism" in the aspects of robot performance evaluation, optimization design, control algorithm development, and parameter tuning. Some of the key parameters that cannot be directly obtained or are not accurate enough need to be tested to obtain accurate values.

从机电耦合的角度出发,工业机器人是复杂的多轴耦合机电系统,机械子系统与控制子系统之间存在相互影响的耦合作用,其中的关键过程为驱动电机的电枢电流在磁场中受到作用力产生驱动力矩,两者之间的转换系数即转矩系数是机电耦合系统的关键参数。电机转矩系数的设计值一般标注在电机铭牌上,但是当机器人应用在表面喷涂、机械加工等工业场景中时,运行环境较为恶劣,受到粉尘、温度湿度变化以及机械结构装配质量等因素的影响,关节电机的转矩系数往往与设计值有一定偏差,并且电机运转方向不同也可能带来变化。From the perspective of electromechanical coupling, industrial robots are complex multi-axis coupled electromechanical systems. There is a coupling effect between the mechanical subsystem and the control subsystem. The key process is that the armature current of the driving motor is affected in the magnetic field. The force produces the driving torque, and the conversion coefficient between the two, the torque coefficient, is the key parameter of the electromechanical coupling system. The design value of the motor torque coefficient is generally marked on the motor nameplate, but when the robot is used in industrial scenarios such as surface spraying and machining, the operating environment is relatively harsh, and it is affected by factors such as dust, temperature and humidity changes, and mechanical structure assembly quality. , the torque coefficient of the joint motor often has a certain deviation from the design value, and the different running directions of the motor may also bring changes.

驱动力矩是驱动机械结构运行的作用力,可以通过机器人动力学理论建模得到,相关方法已经较为成熟,但是由于关节实际运行时的摩擦力很难通过理论建模准确得到,机器人的实际动力学行为与理论动力学模型存在较大差距,在应用动力学模型前应当对关节的摩擦力系数也进行测试。The driving torque is the force that drives the operation of the mechanical structure, which can be obtained by modeling the robot dynamics theory. The related methods are relatively mature. However, because the friction force during the actual operation of the joint is difficult to obtain accurately through theoretical modeling, the actual dynamics of the robot is difficult to obtain. There is a big gap between the behavior and the theoretical dynamic model, and the friction coefficient of the joint should also be tested before applying the dynamic model.

通常,在电机安装前可以使用测功机对其进行转矩系数的测量,再以该系数和动力学模型为基础对安装后的机器人关节摩擦力模型进行测试,但是这一方法没有考虑到机械装配和运行环境对电机参数的影响,且步骤繁琐,不能满足机器人在安装运行一段时间后或工作环境变化后的参数快速校正的需求。Usually, a dynamometer can be used to measure the torque coefficient of the motor before it is installed, and then the friction model of the robot joint after installation is tested based on the coefficient and the dynamic model, but this method does not consider the mechanical The influence of the assembly and operating environment on the motor parameters, and the steps are cumbersome, which cannot meet the needs of the rapid correction of the parameters of the robot after installation and operation for a period of time or after the working environment changes.

目前,对于工业机器人在实际应用中驱动电机转矩系数的测试校正还缺乏相关研究,面向机器人本体的电机转矩系数与关节摩擦力的快速测试还有待完善,为满足机电系统精准建模的需要,急需提出一种面向机器人的关节电机转矩系数与摩擦力同步测试方法。At present, there is still a lack of relevant research on the test and correction of the torque coefficient of the drive motor for industrial robots in practical applications. The rapid test of the motor torque coefficient and joint friction for the robot body needs to be improved. In order to meet the needs of accurate modeling of electromechanical systems , it is urgent to propose a robot-oriented joint motor torque coefficient and friction synchronous test method.

发明内容SUMMARY OF THE INVENTION

本发明旨在至少在一定程度上解决相关技术中的技术问题之一。The present invention aims to solve one of the technical problems in the related art at least to a certain extent.

为此,本发明的目的在于提出一种机器人的关节电机转矩系数与摩擦力同步测试方法和装置,其目的是为了充分考虑机械结构与工作环境对关键机电参数的影响,实现关节电机转矩系数与正反向摩擦力系数的同步测试,提高测试精度。Therefore, the purpose of the present invention is to propose a method and device for synchronizing the torque coefficient and friction force of the joint motor of a robot. Simultaneous test of coefficient and positive and negative friction coefficient to improve test accuracy.

根据本发明的机器人的关节电机转矩系数与摩擦力同步测试方法包括:The method for synchronizing the torque coefficient and friction force of the joint motor of the robot according to the present invention includes:

基于虚功原理建立机器人的动力学模型,获取机器人的机械结构参数,根据动力学模型和机械结构参数计算机器人的各驱动关节的驱动力矩解析表达式;获取预先存储的实验数据,其中,实验数据包括:各待测试关节的角度、角速度、角加速度及驱动电流,以及按照角速度的方向将实验数据分别存储于两个数组中;将不同实验位置的两个数组中的实验数据分别输入驱动力矩解析表达式得到正向力矩和反向力矩,将正向力矩、反向力矩和驱动电流进行拟合处理,得到不同实验位置的拟合线性函数的斜率和截距;根据拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将拟合线性函数的截距取相反数得到待测试关节正向与反向运动时的摩擦力系数;计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。Establish a dynamic model of the robot based on the principle of virtual work, obtain the mechanical structure parameters of the robot, and calculate the analytical expression of the driving torque of each driving joint of the robot according to the dynamic model and the mechanical structure parameters; obtain the pre-stored experimental data, among which the experimental data Including: the angle, angular velocity, angular acceleration and driving current of each joint to be tested, and storing the experimental data in two arrays according to the direction of the angular velocity; respectively inputting the experimental data in the two arrays at different experimental positions into the driving torque analysis The forward torque and reverse torque are obtained from the expression, and the forward torque, reverse torque and driving current are fitted to obtain the slope and intercept of the fitted linear function at different experimental positions; The ratio of the reduction ratio of the joint to be tested is used to obtain the motor torque coefficient of the joint to be tested in the forward and reverse motions, and the inverse of the intercept of the fitted linear function is obtained to obtain the frictional force of the joint to be tested in the forward and reverse motions Calculate the average value of the motor torque coefficient and friction coefficient at different experimental positions, and obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion.

另外,根据本发明的机器人的关节电机转矩系数与摩擦力同步测试方法还可以具有如下附加的技术特征:In addition, the joint motor torque coefficient and friction synchronization test method of the robot according to the present invention may also have the following additional technical features:

根据本发明的一些实施例,还可以包括确定不同实验位置和待测试关节;保持其他关节不动,在不同实验位置控制待测试关节做梯形加减速往复运动,其中,待测试关节匀速运动阶段的速度与待测试关节在实际工作中的角速度大小呈正相关关系。获取在不同实验位置的实验数据。According to some embodiments of the present invention, it may also include determining different experimental positions and joints to be tested; keeping other joints still, and controlling the joints to be tested to perform trapezoidal acceleration and deceleration reciprocating motions at different experimental positions, wherein, the joints to be tested are in the uniform motion stage. There is a positive correlation between the speed and the angular velocity of the joint to be tested in actual work. Acquire experimental data at different experimental locations.

根据本发明的一些实施例,获取在不同实验位置的实验数据,包括通过绝对值编码器获取待测试关节的角度;通过光栅编码器获取待测试关节的速度;通过对瞬时速度做差分获取待测试关节的角加速度;通过霍尔传感器获取待测试关节的驱动电流。According to some embodiments of the present invention, obtaining experimental data at different experimental positions includes obtaining the angle of the joint to be tested through an absolute encoder; obtaining the speed of the joint to be tested through a grating encoder; The angular acceleration of the joint; the drive current of the joint to be tested is obtained through the Hall sensor.

根据本发明的一些实施例,驱动力矩解析表达式为:According to some embodiments of the present invention, the analytical expression of the driving torque is:

Figure GDA0003038459450000021
Figure GDA0003038459450000021

其中,τ为待测试关节实时的驱动力矩,

Figure GDA0003038459450000022
为驱动力矩中的惯性力分量,
Figure GDA0003038459450000023
为驱动力矩中的离心力与科氏力分量,G(θ)g为驱动力矩中的重力分量。另外,使用待测试关节的正向运动时的实验数据计算获得的为正向运动时的理论驱动力矩;使用待测试关节的反向运动时的实验数据计算获得的为反向运动时的理论驱动力矩。Among them, τ is the real-time driving torque of the joint to be tested,
Figure GDA0003038459450000022
is the inertial force component in the driving torque,
Figure GDA0003038459450000023
are the centrifugal force and Coriolis force components in the driving torque, and G(θ)g is the gravity component in the driving torque. In addition, using the experimental data during the forward motion of the joint to be tested is the theoretical driving torque during the forward motion; using the experimental data during the reverse motion of the joint to be tested is the theoretical drive during the reverse motion. moment.

根据本发明的一些实施例,拟合使用的公式为:According to some embodiments of the present invention, the formula used for fitting is:

Figure GDA0003038459450000031
Figure GDA0003038459450000031

其中,Kt为待测试关节电机的转矩系数,Iq为待测试关节电机的驱动电流,f为待测试关节运行所受到的摩擦力,i为待测试关节的减速比。τ为待测试关节的理论驱动力矩,包含待测试关节的正向运动时的理论驱动力矩和反向运动时的理论驱动力矩,当使用待测试关节的正向运动时的理论驱动力矩时,拟合获得的函数为正向运动拟合函数;当使用待测试关节的反向运动时的理论驱动力矩时,拟合获得的函数为反向运动拟合函数。Among them, K t is the torque coefficient of the joint motor to be tested, I q is the drive current of the joint motor to be tested, f is the friction force of the joint to be tested, and i is the reduction ratio of the joint to be tested. τ is the theoretical driving torque of the joint to be tested, including the theoretical driving torque of the joint to be tested in the forward motion and the theoretical driving torque of the reverse motion. When the theoretical driving torque of the joint to be tested is used in the forward motion, the The function obtained by fitting is the forward motion fitting function; when the theoretical driving torque of the joint to be tested is used in the reverse motion, the function obtained by fitting is the reverse motion fitting function.

根据本发明的一些实施例,在不同实验位置控制待测试关节做梯形加减速往复运动包括在机器人的预设工作空间内,以间距相等的方式均匀地选取若干个机器人末端执行器的坐标定位生成不同实验位置;在不同实验位置控制待测关节实施相同运动行程、相同速度规划的运动。According to some embodiments of the present invention, controlling the joint to be tested to perform trapezoidal acceleration and deceleration reciprocating motion at different experimental positions is included in the preset working space of the robot, and uniformly selects the coordinate positioning of several robot end effectors with equal spacing to generate Different experimental positions; control the joints to be tested in different experimental positions to implement the same motion stroke and the same speed plan.

为达到上述目的,本发明第二方面实施例提出了一种机器人的关节电机转矩系数与摩擦力同步测试装置,包括:建模模块、获取模块、第一计算模块、第二计算模块、第三计算模块,其中,建模模块用于基于虚功原理建立机器人的动力学模型,获取机器人的机械结构参数,根据动力学模型和机械结构参数计算机器人的各驱动关节的驱动力矩解析表达式;获取模块用于获取预先存储的实验数据,其中,实验数据包括:各待测试关节的角度、角速度、角加速度及驱动电流,以及按照角速度的方向将实验数据分别存储于两个数组中;第一计算模块用于将不同实验位置的两个数组中的实验数据分别输入驱动力矩解析表达式得到正向力矩和反向力矩,将正向力矩、反向力矩和驱动电流进行拟合处理,得到不同实验位置的拟合线性函数的斜率和截距;第二计算模块用于根据拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将拟合线性函数的截距取相反数得到待测试关节正向与反向运动时的摩擦力系数;第三计算模块用于计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。In order to achieve the above purpose, the embodiment of the second aspect of the present invention proposes a synchronous test device for the torque coefficient and friction force of a joint motor of a robot, including: a modeling module, an acquisition module, a first calculation module, a second calculation module, a third Three calculation modules, wherein the modeling module is used to establish a dynamic model of the robot based on the principle of virtual work, obtain the mechanical structure parameters of the robot, and calculate the analytical expression of the driving torque of each driving joint of the robot according to the dynamic model and the mechanical structure parameters; The acquisition module is used for acquiring pre-stored experimental data, wherein the experimental data includes: the angle, angular velocity, angular acceleration and driving current of each joint to be tested, and the experimental data are stored in two arrays respectively according to the direction of the angular velocity; the first The calculation module is used to input the experimental data in the two arrays of different experimental positions into the analytical expression of the driving torque respectively to obtain the forward torque and reverse torque, and fit the forward torque, reverse torque and driving current to obtain different values. The slope and intercept of the fitted linear function of the experimental position; the second calculation module is used to obtain the motor torque when the joint to be tested moves forward and backward according to the ratio of the slope of the fitted linear function to the reduction ratio of the joint to be tested The coefficient of friction is obtained by taking the inverse of the intercept of the fitted linear function to obtain the friction coefficient of the joint to be tested in forward and reverse motion; the third calculation module is used to calculate the average of the motor torque coefficient and friction coefficient at different experimental positions value to obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion.

另外,根据本发明上述实施例的机器人的关节电机转矩系数与摩擦力同步测试装置,还可以具有如下附加的技术特征:In addition, the joint motor torque coefficient and friction force synchronization test device of the robot according to the above-mentioned embodiment of the present invention may also have the following additional technical features:

进一步地,在本申请实施例的一种可能的实现方式中,还可以包括:Further, in a possible implementation manner of the embodiment of the present application, it may also include:

确定模块,用于确定不同实验位置和待测试关节;A determination module for determining different experimental positions and joints to be tested;

实验模块,用于保持其他关节不动,在不同实验位置控制待测试关节做梯形加减速往复运动,其中,待测试关节匀速运动阶段的速度与待测试关节在实际工作中的角速度大小呈正相关关系。The experimental module is used to keep other joints still, and control the joint to be tested to perform trapezoidal acceleration and deceleration reciprocating motion at different experimental positions. The speed of the joint to be tested in the uniform motion stage is positively correlated with the angular velocity of the joint to be tested in actual work. .

记录模块,用于获取在不同实验位置的实验数据。The recording module is used to obtain experimental data at different experimental locations.

进一步地,在本申请实施例的一种可能的实现方式中,记录模块,具体用于通过绝对值编码器获取待测试关节的角度、通过光栅编码器获取待测试关节的速度、通过对瞬时速度做差分获取待测试关节的角加速度、通过霍尔传感器获取待测试关节的驱动电流。Further, in a possible implementation manner of the embodiment of the present application, the recording module is specifically configured to obtain the angle of the joint to be tested through an absolute value encoder, obtain the speed of the joint to be tested through a grating encoder, and obtain the speed of the joint to be tested through an absolute value encoder. Do differential to obtain the angular acceleration of the joint to be tested, and obtain the drive current of the joint to be tested through the Hall sensor.

进一步地,在本申请实施例的一种可能的实现方式中,第一计算模块,用于计算驱动力矩解析表达式为:Further, in a possible implementation manner of the embodiment of the present application, the first calculation module is used to calculate the analytical expression of the driving torque as follows:

Figure GDA0003038459450000041
Figure GDA0003038459450000041

其中,τ为待测试关节实时的驱动力矩,

Figure GDA0003038459450000042
为驱动力矩中的惯性力分量,
Figure GDA0003038459450000043
为驱动力矩中的离心力与科氏力分量,G(θ)g为驱动力矩中的重力分量。另外,使用待测试关节的正向运动时的实验数据计算获得的为正向运动时的理论驱动力矩;使用待测试关节的反向运动时的实验数据计算获得的为反向运动时的理论驱动力矩。Among them, τ is the real-time driving torque of the joint to be tested,
Figure GDA0003038459450000042
is the inertial force component in the driving torque,
Figure GDA0003038459450000043
are the centrifugal force and Coriolis force components in the driving torque, and G(θ)g is the gravity component in the driving torque. In addition, using the experimental data during the forward motion of the joint to be tested is the theoretical driving torque during the forward motion; using the experimental data during the reverse motion of the joint to be tested is the theoretical drive during the reverse motion. moment.

本发明实施例提供的机器人的关节电机转矩系数与摩擦力同步测试方法可以包含如下有益效果:The method for synchronizing the torque coefficient and friction force of the joint motor of the robot provided by the embodiment of the present invention may include the following beneficial effects:

基于虚功原理建立机器人的动力学模型,获取机器人的机械结构参数,根据动力学模型和机械结构参数计算机器人的各驱动关节的驱动力矩解析表达式;获取预先存储的实验数据,其中,实验数据包括:各待测试关节的角度、角速度、角加速度及驱动电流,以及按照角速度的方向将实验数据分别存储于两个数组中;将不同实验位置的两个数组中的实验数据分别输入驱动力矩解析表达式得到正向力矩和反向力矩,将正向力矩、反向力矩和驱动电流进行拟合处理,得到不同实验位置的拟合线性函数的斜率和截距;根据拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将拟合线性函数的截距取相反数得到待测试关节正向与反向运动时的摩擦力系数;计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。由此,在充分考虑机械结构与工作环境对关键机电参数影响的情况下,实现关节电机转矩系数与正反向摩擦力系数的同步测试,提高测试精度。Establish a dynamic model of the robot based on the principle of virtual work, obtain the mechanical structure parameters of the robot, and calculate the analytical expression of the driving torque of each driving joint of the robot according to the dynamic model and the mechanical structure parameters; obtain the pre-stored experimental data, among which the experimental data Including: the angle, angular velocity, angular acceleration and driving current of each joint to be tested, and storing the experimental data in two arrays according to the direction of the angular velocity; respectively inputting the experimental data in the two arrays at different experimental positions into the driving torque analysis The forward torque and reverse torque are obtained from the expression, and the forward torque, reverse torque and driving current are fitted to obtain the slope and intercept of the fitted linear function at different experimental positions; The ratio of the reduction ratio of the joint to be tested is used to obtain the motor torque coefficient of the joint to be tested in the forward and reverse motions, and the inverse of the intercept of the fitted linear function is obtained to obtain the frictional force of the joint to be tested in the forward and reverse motions Calculate the average value of the motor torque coefficient and friction coefficient at different experimental positions, and obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion. Therefore, under the condition of fully considering the influence of mechanical structure and working environment on key electromechanical parameters, the synchronous test of the torque coefficient of joint motor and the coefficient of forward and reverse friction force is realized, and the test accuracy is improved.

本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the present invention will be set forth, in part, from the following description, and in part will be apparent from the following description, or may be learned by practice of the invention.

附图说明Description of drawings

本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:The above and/or additional aspects and advantages of the present invention will become apparent and readily understood from the following description of embodiments taken in conjunction with the accompanying drawings, wherein:

图1为本发明实施例的一种机器人的关节电机转矩系数与摩擦力同步测试方法的流程示意图;1 is a schematic flowchart of a method for synchronizing the torque coefficient and friction force of a joint motor of a robot according to an embodiment of the present invention;

图2为适用于本发明的一种六轴混联喷涂机器人的结构示意图;2 is a schematic structural diagram of a six-axis hybrid spraying robot suitable for the present invention;

图3是六轴混联喷涂机器人待测试关节驱动电流曲线与待测试关节角速度曲线对比图;Fig. 3 is a comparison diagram of the drive current curve of the joint to be tested and the angular velocity curve of the joint to be tested of the six-axis hybrid spraying robot;

图4是六轴混联喷涂机器人关节测试实验的计算驱动力矩曲线图;Fig. 4 is the calculated driving torque curve diagram of the joint test experiment of the six-axis hybrid spraying robot;

图5是六轴混联喷涂机器人正向运动阶段实际驱动电流散点图;Figure 5 is a scatter diagram of the actual drive current in the forward motion stage of the six-axis hybrid spraying robot;

图6是六轴混联喷涂机器人正向运动阶段计算驱动力矩散点图;Figure 6 is a scatter diagram of the calculated driving torque in the forward motion stage of the six-axis hybrid spraying robot;

图7是六轴混联喷涂机器人正向运动阶段实际驱动电流与计算驱动力矩线性拟合示意图;Fig. 7 is the linear fitting schematic diagram of the actual drive current and the calculated drive torque in the forward motion stage of the six-axis hybrid spraying robot;

图8是六轴混联喷涂机器人反向运动阶段实际驱动电流散点图;Figure 8 is a scatter diagram of the actual drive current in the reverse motion stage of the six-axis hybrid spraying robot;

图9是六轴混联喷涂机器人反向运动阶段计算驱动力矩散点图;Figure 9 is a scatter diagram of the calculated driving torque in the reverse motion stage of the six-axis hybrid spraying robot;

图10是六轴混联喷涂机器人反向运动阶段实际驱动电流与计算驱动力矩线性拟合示意图;Figure 10 is a schematic diagram of the linear fitting of the actual drive current and the calculated drive torque in the reverse motion stage of the six-axis hybrid spraying robot;

图11是六轴混联喷涂机器人关节验证实验的预测驱动电流与实际驱动电流对比图;Figure 11 is a comparison diagram of the predicted driving current and the actual driving current of the joint verification experiment of the six-axis hybrid spraying robot;

图12为本发明实施例的一种机器人的关节电机转矩系数与摩擦力同步测试装置的机构示意图。FIG. 12 is a schematic structural diagram of a device for synchronizing the torque coefficient and friction force of a joint motor of a robot according to an embodiment of the present invention.

具体实施方式Detailed ways

下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。Embodiments of the present invention are described in detail below, examples of which are illustrated in the accompanying drawings, wherein the same or similar reference numerals refer to the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary, and are intended to explain the present invention and should not be construed as limiting the present invention.

本发明的目的是为了充分考虑机械结构与工作环境对关键机电参数的影响,实现关节电机转矩系数与正反向摩擦力系数的同步测试,提高测试精度。The purpose of the invention is to fully consider the influence of the mechanical structure and the working environment on the key electromechanical parameters, realize the synchronous test of the torque coefficient of the joint motor and the positive and negative friction coefficient, and improve the test accuracy.

下面参照附图描述根据本发明实施例提出的机器人的关节电机转矩系数与摩擦力同步测试方法和装置。The following describes the method and device for synchronizing the torque coefficient and friction force of the joint motor of the robot according to the embodiments of the present invention with reference to the accompanying drawings.

图1为本发明实施例的一种机器人的关节电机转矩系数与摩擦力同步测试方法的流程示意图。如图1所示,该机器人的关节电机转矩系数与摩擦力同步测试方法,包括:FIG. 1 is a schematic flowchart of a method for synchronizing the torque coefficient and friction force of a joint motor of a robot according to an embodiment of the present invention. As shown in Figure 1, the joint motor torque coefficient and friction synchronization test method of the robot includes:

步骤101,基于虚功原理建立机器人的动力学模型,获取机器人的机械结构参数,根据动力学模型和机械结构参数计算机器人的各驱动关节的驱动力矩解析表达式。Step 101 , establishing a dynamic model of the robot based on the virtual work principle, obtaining mechanical structure parameters of the robot, and calculating an analytical expression for the driving torque of each driving joint of the robot according to the dynamic model and the mechanical structure parameters.

其中,机器人可以理解为工业机器人等,不同类型的工业机器人会应用于不同的工作环境中,由于工作环境的不同,会造成机器人运行参数的不同或改变。Among them, robots can be understood as industrial robots, etc. Different types of industrial robots will be used in different working environments. Due to different working environments, the operating parameters of the robots will be different or changed.

具体的,将机器人视作多刚体的连杆机构,基于虚功原理建立其理论动力学模型,在三维软件中建立机器人的三维模型,并通过三维软件采集机器人各构件的质量、长度、转动惯量等参数,带入基于虚功原理建立的理论动力学模型,得到机器人各个关节驱动力矩的解析表达式。其中,三维软件可以是Solidworks、Creo(Pro/E)等。Specifically, the robot is regarded as a multi-rigid linkage mechanism, its theoretical dynamic model is established based on the principle of virtual work, the 3D model of the robot is established in 3D software, and the mass, length, and moment of inertia of each component of the robot are collected through the 3D software. and other parameters, brought into the theoretical dynamic model based on the virtual work principle, and obtained the analytical expression of the driving torque of each joint of the robot. The three-dimensional software may be Solidworks, Creo (Pro/E), and the like.

步骤102,获取预先存储的实验数据,其中,实验数据包括:各待测试关节的角度、角速度、角加速度及驱动电流,以及按照角速度的方向将实验数据分别存储于两个数组中。Step 102: Acquire pre-stored experimental data, wherein the experimental data includes: the angle, angular velocity, angular acceleration and driving current of each joint to be tested, and the experimental data is stored in two arrays according to the direction of the angular velocity.

其中,预先存储的实验数据可以理解为提前确定的待测试关节,在机器人的预设工作空间内,以间距相等的方式均匀地选取若干个机器人末端执行器的坐标定位生成不同实验位置,在不同实验位置,保持其他关节不动,控制待测试关节做梯形加减速往复运动,从而获取的多组待测试关节在不同实验位置的实验数据。实验数据包括:待测试关节的角度、角速度、角加速度及驱动电流。而获得各个实验数据的方法可以是通过绝对值编码器获取待测试关节的角度;通过光栅编码器获取待测试关节的速度;通过对瞬时速度做差分获取待测试关节的角加速度;通过霍尔传感器获取待测试关节的驱动电流。另外,多组待测试关节在不同实验位置的实验数据中,还会根据角速度方向将其区分为两个数组。同时,待测试关节匀速运动阶段的速度与待测试关节在实际工作中的角速度大小呈正相关关系。Among them, the pre-stored experimental data can be understood as the joints to be tested determined in advance. In the preset working space of the robot, the coordinate positioning of several robot end effectors are evenly selected with equal spacing to generate different experimental positions. In the experimental position, keep other joints still, and control the joints to be tested to do trapezoidal acceleration and deceleration reciprocating motion, so as to obtain experimental data of multiple groups of joints to be tested in different experimental positions. The experimental data includes: the angle, angular velocity, angular acceleration and driving current of the joint to be tested. The method of obtaining each experimental data may be to obtain the angle of the joint to be tested through an absolute value encoder; to obtain the speed of the joint to be tested through a grating encoder; to obtain the angular acceleration of the joint to be tested by differentiating the instantaneous speed; Obtain the drive current of the joint to be tested. In addition, in the experimental data of multiple groups of joints to be tested at different experimental positions, they will also be divided into two arrays according to the direction of angular velocity. At the same time, the speed of the joint to be tested in the uniform motion stage has a positive correlation with the angular velocity of the joint to be tested in actual work.

具体的,从多组根据待测试关节在不同实验位置得到的实验数据中,选择其中的一组实验数据的两个数组,在两个数组中获取预先经过实验过程记录的各待测试关节的角度、角速度、角加速度及驱动电流的实验数据。Specifically, from the multiple sets of experimental data obtained according to the joints to be tested at different experimental positions, two arrays of one set of experimental data are selected, and the angles of the joints to be tested recorded in advance through the experimental process are obtained from the two arrays , angular velocity, angular acceleration and experimental data of driving current.

步骤103,将不同实验位置的两个数组中的实验数据分别输入驱动力矩解析表达式得到正向力矩和反向力矩,将正向力矩、反向力矩和驱动电流进行拟合处理,得到不同实验位置的拟合线性函数的斜率和截距。Step 103: Input the experimental data in the two arrays of different experimental positions into the driving torque analytical expression respectively to obtain the forward torque and reverse torque, and perform fitting processing on the forward torque, reverse torque and driving current to obtain different experiments. The slope and intercept of the fitted linear function of position.

其中,驱动力矩解析表达式为公式(1):Among them, the analytical expression of the driving torque is formula (1):

Figure GDA0003038459450000061
Figure GDA0003038459450000061

其中,τ为待测试关节实时的驱动力矩,

Figure GDA0003038459450000062
为驱动力矩中的惯性力分量,
Figure GDA0003038459450000063
为驱动力矩中的离心力与科氏力分量,G(θ)g为驱动力矩中的重力分量,另外,使用待测试关节的正向运动时的实验数据计算获得的为正向运动时的理论驱动力矩;使用待测试关节的反向运动时的实验数据计算获得的为反向运动时的理论驱动力矩。另外,在工业机器人中,对于没有高速运动需求的轻型机器人,在其驱动力矩中,惯性分量与重力分量将占据主要地位,因此后续计算可以忽略离心力与科氏力分量。Among them, τ is the real-time driving torque of the joint to be tested,
Figure GDA0003038459450000062
is the inertial force component in the driving torque,
Figure GDA0003038459450000063
is the centrifugal force and Coriolis force components in the driving torque, G(θ)g is the gravity component in the driving torque, in addition, the theoretical driving in the forward movement is calculated using the experimental data during the forward movement of the joint to be tested. Torque; the theoretical driving torque during reverse movement is calculated using the experimental data of the joint to be tested during reverse movement. In addition, in industrial robots, for light robots without high-speed motion requirements, inertial and gravitational components will occupy the dominant position in the driving torque, so centrifugal force and Coriolis force components can be ignored in subsequent calculations.

另外,进行拟合使用的公式为公式(2):In addition, the formula used for fitting is formula (2):

Figure GDA0003038459450000064
Figure GDA0003038459450000064

其中,Kt为待测试关节电机的转矩系数,Iq为待测试关节电机的驱动电流,f为待测试关节运行所受到的摩擦力,i为待测试关节的减速比,τ为待测试关节的理论驱动力矩,包含待测试关节的正向运动时的理论驱动力矩和反向运动时的理论驱动力矩,当使用待测试关节的正向运动时的理论驱动力矩时,拟合获得的函数为正向运动拟合函数;当使用待测试关节的反向运动时的理论驱动力矩时,拟合获得的函数为反向运动拟合函数。Among them, K t is the torque coefficient of the joint motor to be tested, I q is the drive current of the joint motor to be tested, f is the frictional force of the joint to be tested, i is the reduction ratio of the joint to be tested, and τ is the joint to be tested. The theoretical driving torque of the joint, including the theoretical driving torque in the forward motion of the joint to be tested and the theoretical driving torque in the reverse motion, when using the theoretical driving torque in the forward motion of the joint to be tested, the function obtained by fitting is the forward motion fitting function; when the theoretical driving torque of the joint to be tested is used in the reverse motion, the function obtained by fitting is the reverse motion fitting function.

具体的,将获得的多组预先存储的待测试关节在不同实验位置的角度、角速度、角加速度数据,按照角速度方向的不同,代入待测试关节的驱动力矩解析表达式,计算得到待测试关节在不同实验位置、不同运动方向的正向驱动力矩和反向驱动力矩。Specifically, the obtained multiple pre-stored angle, angular velocity, and angular acceleration data of the joint to be tested at different experimental positions are substituted into the analytical expression of the driving torque of the joint to be tested according to the direction of the angular velocity, and it is calculated that the joint to be tested is at Forward driving torque and reverse driving torque for different experimental positions and different motion directions.

然后将计算得到待测试关节在不同实验位置的正向驱动力矩和反向驱动力矩中的正向驱动力矩与采集的待测试关节电机驱动电流,通过最小二乘法做线性拟合,得到在不同实验位置的正向运动拟合线性函数的斜率与截距。Then, the forward driving torque in the forward driving torque and the reverse driving torque of the joint to be tested at different experimental positions are calculated and the collected motor driving current of the joint to be tested is linearly fitted by the least square method, and the results obtained in different experiments are obtained by linear fitting. The forward motion of the position fits the slope and intercept of a linear function.

再将计算得到待测试关节在不同实验位置的正向驱动力矩和反向驱动力矩中的反向驱动力矩与采集的待测试关节电机驱动电流,通过最小二乘法做线性拟合,得到在不同实验位置的反向运动拟合线性函数的斜率与截距。Then, calculate the forward driving torque and the reverse driving torque of the joint to be tested at different experimental positions and the collected driving current of the joint motor to be tested. The inverse motion of the position fits the slope and intercept of a linear function.

步骤104,根据拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将拟合线性函数的截距取相反数得到待测试关节正向与反向运动时的摩擦力系数。Step 104, obtain the motor torque coefficient when the joint to be tested is moving forward and backward according to the ratio of the slope of the fitted linear function and the reduction ratio of the joint to be tested, and take the inverse of the intercept of the fitted linear function to obtain the value to be tested. Coefficient of friction for joint forward and reverse motion.

具体的,将待测试关节在不同实验位置的正向拟合线性函数的斜率、反向拟合线性函数的斜率和待测试关节的减速比的比值作为待测试关节正向运动时的电机转矩系数与反向运动时的电机转矩系数,将待测试关节在不同实验位置的正向拟合线性函数的截距、反向拟合线性函数的截距的相反数作为待测试关节正向运动时的摩擦力系数与反向运动时的摩擦力系数。Specifically, the ratio of the slope of the forward fitting linear function, the slope of the backward fitting linear function and the reduction ratio of the joint to be tested at different experimental positions of the joint to be tested is taken as the motor torque when the joint to be tested moves forward Coefficient and motor torque coefficient during reverse motion, the intercept of the forward fitting linear function of the joint to be tested at different experimental positions and the opposite number of the intercept of the reverse fitting linear function as the forward motion of the joint to be tested The coefficient of friction at time and the coefficient of friction at the time of reverse motion.

步骤105,计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。Step 105: Calculate the average value of the motor torque coefficient and friction coefficient at different experimental positions, and obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion.

具体的,计算待测试关节在不同实验位置的正向与反向运动时的电机转矩系数和摩擦力系数各自的平均值,得到正向运动时与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。Specifically, calculate the respective average values of the motor torque coefficient and friction coefficient of the joint to be tested during the forward and reverse motions at different experimental positions, and obtain the target torque of the joint to be tested during the forward motion and the reverse motion. Coefficient and target friction value.

举例而言,以如图2所示的六轴混联喷涂机器人为例,六轴混联喷涂机器人共包含六个关节,其主体为由第一旋转关节1、第二旋转关节2、第三旋转关节3通过各连杆驱动喷枪7运动而组成的平面并联机构8。此外,第四旋转关节4实现平面并联机构8绕竖直轴的转动,第五旋转关节5实现第四旋转关节4及平面并联机构8绕水平轴的转动,第六旋转关节6实现第五旋转关节5、第四旋转关节4及平面并联机构8绕竖直轴的转动。For example, taking the six-axis hybrid spraying robot shown in Figure 2 as an example, the six-axis hybrid spraying robot includes six joints in total, and its main body is composed of a first rotating joint 1, a second rotating joint 2, a third rotating joint The rotating joint 3 drives the spray gun 7 to move through each connecting rod to form a plane parallel mechanism 8 . In addition, the fourth rotating joint 4 realizes the rotation of the plane parallel mechanism 8 around the vertical axis, the fifth rotating joint 5 realizes the rotation of the fourth rotating joint 4 and the plane parallel mechanism 8 around the horizontal axis, and the sixth rotating joint 6 realizes the fifth rotation The rotation of the joint 5, the fourth rotating joint 4 and the plane parallel mechanism 8 around the vertical axis.

首先基于虚功原理建立机械人的动力学模型,并将其写成如公式(3)所示的标准动力学方程形式:Firstly, the dynamic model of the robot is established based on the principle of virtual work, and it is written in the standard dynamic equation form as shown in formula (3):

Figure GDA0003038459450000081
Figure GDA0003038459450000081

其中,τ为待测试关节实时的驱动力矩,

Figure GDA0003038459450000082
为驱动力矩中的惯性力分量,
Figure GDA0003038459450000083
为驱动力矩中的离心力与科氏力分量,G(θ)g为驱动力矩中的重力分量,另外,使用待测试关节的正向运动时的实验数据计算获得的为正向运动时的理论驱动力矩;使用待测试关节的反向运动时的实验数据计算获得的为反向运动时的理论驱动力矩。由于该喷涂机器人不是有高速运动需求的轻型机器人,在其驱动力矩中,惯性分量与重力分量将占据主要地位,因此后续计算忽略离心力与科氏力分量。Among them, τ is the real-time driving torque of the joint to be tested,
Figure GDA0003038459450000082
is the inertial force component in the driving torque,
Figure GDA0003038459450000083
is the centrifugal force and Coriolis force components in the driving torque, G(θ)g is the gravity component in the driving torque, in addition, the theoretical driving in the forward movement is calculated using the experimental data during the forward movement of the joint to be tested. Torque; the theoretical driving torque during reverse movement is calculated using the experimental data of the joint to be tested during reverse movement. Since the spraying robot is not a light robot with high-speed motion requirements, inertial and gravitational components will occupy the dominant position in its driving torque, so centrifugal force and Coriolis force components are ignored in subsequent calculations.

在本实施例中,以该六轴混联喷涂机器人的关节1,即大臂驱动关节为待测试关节实施测试步骤,在三维软件中对该六轴混联喷涂机器人进行建模,将在三维软件中读取到的机器人的质量、长度、转动惯量等机械机构参数代入公式(3),得到待测试关节1的驱动力矩的解析表达式为公式(4):In this embodiment, the joint 1 of the six-axis hybrid spraying robot, that is, the big arm drive joint, is used as the joint to be tested to perform the test steps, and the six-axis hybrid spraying robot is modeled in the three-dimensional software. The mechanical parameters such as the mass, length, moment of inertia of the robot read in the software are substituted into formula (3), and the analytical expression of the driving torque of joint 1 to be tested is obtained as formula (4):

Figure GDA0003038459450000084
Figure GDA0003038459450000084

其中,θ1、θ2、θ3分别为待测试关节1、2、3的角度(单位:rad),

Figure GDA0003038459450000085
分别为待测试关节1、2、3的角加速度(单位:rad/s2),τd1为待测试关节1的驱动力矩(单位:N·m)。Among them, θ 1 , θ 2 , θ 3 are the angles of joints 1, 2, and 3 to be tested (unit: rad), respectively,
Figure GDA0003038459450000085
are the angular accelerations of the joints 1, 2 and 3 to be tested (unit: rad/s 2 ), respectively, and τ d1 is the driving torque of the joint 1 to be tested (unit: N·m).

图3为待测试关节1在实验中记录的关节驱动电流与关节角速度曲线,其中半透明虚线框中圈出的驱动电流数值原本为正值,考虑到此时驱动力矩反向而采集到的电机电流只能为正值,此处进行了去绝对值处理。由实验结果可见,在待测试关节速度由负向转为正向的时刻,其驱动电流数值存在正向的突变,在待测试关节速度由正向转为负向的时刻,其驱动电流数值存在同等幅度的反向的突变,这是由于待测试关节速度方向转变使得待测试关节受到的摩擦力方向转变导致的。Figure 3 is the joint drive current and joint angular velocity curve recorded in the experiment for the joint 1 to be tested. The drive current value circled in the semi-transparent dashed box is originally a positive value. Considering the reverse driving torque at this time, the collected motor The current can only be a positive value, and the absolute value processing is performed here. It can be seen from the experimental results that when the speed of the joint to be tested changes from negative to positive, there is a positive sudden change in the value of the driving current. When the speed of the joint to be tested changes from positive to negative, the value of the driving current exists. The reverse mutation of the same magnitude is caused by the change in the direction of the speed of the joint to be tested and the change in the direction of the friction force on the joint to be tested.

图4为将待测试关节1在实验中的角度、角速度、角加速度等数据代入式(4)得到的驱动力矩曲线,计算力矩曲线在待测试关节运动速度方向转变时保持连续。Figure 4 is the driving torque curve obtained by substituting the data of the joint 1 to be tested in the experiment, such as angle, angular velocity, angular acceleration, etc. into equation (4).

图5、图6、图8、图9,为按照待测试关节角速度方向实施数据分组存储而得到的待测试关节在正向和反向运转时的实际驱动电流和计算驱动力矩曲线。Figure 5, Figure 6, Figure 8, Figure 9 are the actual driving current and calculated driving torque curves of the joint to be tested during forward and reverse operation obtained by implementing data packet storage according to the direction of the angular velocity of the joint to be tested.

图7、图10,为对待测试关节在正向和反向运转时的实际驱动电流和计算驱动力矩以最小二乘法进行线性拟合的示意图,其中横坐标为实际驱动电流,纵坐标为计算驱动力矩,二者线性拟合效果良好,黑色实线为线性拟合结果,其表达式为公式(5):Figure 7 and Figure 10 are schematic diagrams of linear fitting of the actual drive current and calculated drive torque of the joint to be tested in forward and reverse operation by the least squares method, where the abscissa is the actual drive current, and the ordinate is the calculated drive moment, the linear fitting effect of the two is good, the black solid line is the linear fitting result, and its expression is formula (5):

Figure GDA0003038459450000086
Figure GDA0003038459450000086

其中τd+为待测试关节1在正向运行时的驱动力矩(单位:N·m),τd-为反向运行时的驱动力矩,Iq+为正向运行时的驱动电流(单位:A),Iq-为反向运行时的驱动电流。线性拟合结果表达式的斜率单位为N·m/A,截距单位为N·m。where τ d+ is the driving torque of joint 1 to be tested in forward running (unit: N m), τ d- is the driving torque in reverse running, and I q+ is the driving current in forward running (unit: A ), and I q- is the drive current during reverse operation. The unit of the slope of the linear fitting result expression is N·m/A, and the unit of the intercept is N·m.

另外,在本实施例中,该六轴混联喷涂机器人的待测试关节的减速比在该六轴混联喷涂机器人出厂的时候由厂家设置为121,且在另一些实施例中,六轴混联喷涂机器人的待测试关节的减速比也可以由厂家设置为其他指定值。所以,在本实施例中,结合该待测试关节减速比121,根据公式(2)可以得到电机转矩系数与摩擦力的实验测试结果为公式(6):In addition, in this embodiment, the reduction ratio of the joint to be tested of the six-axis hybrid spraying robot is set by the manufacturer to 121 when the six-axis hybrid spraying robot leaves the factory, and in other embodiments, the six-axis hybrid spraying robot is set to 121 by the manufacturer. The reduction ratio of the joint to be tested of the joint spraying robot can also be set to other specified values by the manufacturer. Therefore, in this embodiment, combined with the joint reduction ratio 121 to be tested, according to formula (2), the experimental test result of the motor torque coefficient and friction force can be obtained as formula (6):

Figure GDA0003038459450000091
Figure GDA0003038459450000091

其中,Kt+为待测试关节1电机在正向运行时的转矩系数,Kt-为待测试关节1电机在反向运行时的转矩系数,f+为待测试关节1电机在正向运行时的待测试关节摩擦力数值,f-为待测试关节1电机在反向运行时的待测试关节摩擦力数值。Among them, K t+ is the torque coefficient of the motor of joint 1 to be tested when it is running in the forward direction, K t- is the torque coefficient of the motor of the joint 1 to be tested when it is running in reverse, and f + is the motor of joint 1 to be tested when it is running in the forward direction. The friction force value of the joint to be tested during operation, f - is the friction force value of the joint to be tested when the motor of the joint 1 to be tested is running in the reverse direction.

在机器人不同实验位置处重复进行6组测试实验,并将实验结果取平均值,得到待测试关节1的驱动电流Iq预测表达式为公式(7):Repeat 6 groups of test experiments at different experimental positions of the robot, and take the average of the experimental results to obtain the predicted expression of the driving current I q of the joint 1 to be tested as formula (7):

Figure GDA0003038459450000092
Figure GDA0003038459450000092

其中,

Figure GDA0003038459450000093
为根据上述步骤得到的待测试关节电机正向转矩系数测试结果,
Figure GDA0003038459450000094
为根据上述步骤得到的待测试关节电机反向转矩系数测试结果,
Figure GDA0003038459450000095
为根据上述步骤得到的待测试关节正向摩擦力数值,
Figure GDA0003038459450000096
为根据上述步骤得到的待测试关节反向摩擦力数值,ω1为待测试关节角速度,其正值表示待测试关节正向运行,负值表示反向运行。in,
Figure GDA0003038459450000093
is the test result of the forward torque coefficient of the joint motor to be tested obtained according to the above steps,
Figure GDA0003038459450000094
is the test result of the reverse torque coefficient of the joint motor to be tested obtained according to the above steps,
Figure GDA0003038459450000095
is the positive friction force value of the joint to be tested obtained according to the above steps,
Figure GDA0003038459450000096
is the value of the reverse friction force of the joint to be tested obtained according to the above steps, ω 1 is the angular velocity of the joint to be tested, a positive value of which indicates that the joint to be tested is running in the forward direction, and a negative value indicates that it is running in the reverse direction.

图11,为将公式(7)中的待测试关节1电流预测模型应用于该待测试关节另一组运动实验的验证结果。其中,实际驱动电流为实验中采集得到的待测试关节的驱动电流曲线,预测驱动电流为将实验中采集得到待测试关节的位置、速度和角速度数据代入公式(7)得到的预测驱动电流曲线,从图中可以看出,两条曲线拟合良好,证明本发明提出的方法完成了对本实施例中机器人待测试关节1电机转矩系数和摩擦力的准确测试。FIG. 11 is a verification result of applying the current prediction model of the joint to be tested 1 in formula (7) to another group of motion experiments of the joint to be tested. Among them, the actual driving current is the driving current curve of the joint to be tested collected in the experiment, and the predicted driving current is the predicted driving current curve obtained by substituting the position, velocity and angular velocity data of the joint to be tested collected in the experiment into formula (7), It can be seen from the figure that the two curves are well fitted, which proves that the method proposed by the present invention has completed the accurate test of the motor torque coefficient and friction force of the joint 1 to be tested of the robot in this embodiment.

由上述可知,在本发明的实施例中,本发明的机器人的关节电机转矩系数与摩擦力同步测试方法,基于虚功原理建立机器人的动力学模型,获取机器人的机械结构参数,根据动力学模型和机械结构参数计算机器人的各驱动关节的驱动力矩解析表达式;获取预先存储的实验数据,其中,实验数据包括:各关节的角度、角速度、角加速度及驱动电流,以及按照角速度的方向将实验数据分别存储于两个数组中;将不同实验位置的两个数组中的实验数据分别输入驱动力矩解析表达式得到正向力矩和反向力矩,将正向力矩、反向力矩和驱动电流进行拟合处理,得到不同实验位置的拟合线性函数的斜率和截距;根据拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将拟合线性函数的截距取相反数得到待测试关节正向与反向运动时的摩擦力系数;计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。由此,在充分考虑机械结构与工作环境对关键机电参数影响的情况下,实现关节电机转矩系数与正反向摩擦力系数的同步测试,提高测试精度。It can be seen from the above that, in the embodiment of the present invention, the joint motor torque coefficient and friction force synchronization test method of the robot of the present invention establishes the dynamic model of the robot based on the principle of virtual work, obtains the mechanical structure parameters of the robot, and according to the dynamics The model and mechanical structure parameters are used to calculate the analytical expression of the driving torque of each driving joint of the robot; the pre-stored experimental data are obtained, wherein the experimental data includes: the angle, angular velocity, angular acceleration and driving current of each joint, and the direction of the angular velocity. The experimental data are stored in two arrays respectively; the experimental data in the two arrays at different experimental positions are respectively input into the analytical expression of the driving torque to obtain the forward torque and reverse torque, and the forward torque, reverse torque and driving current are calculated. Fitting process to obtain the slope and intercept of the fitted linear function at different experimental positions; according to the ratio of the slope of the fitted linear function to the reduction ratio of the joint to be tested, the motor torque of the joint to be tested is obtained during forward and reverse motions Take the inverse of the intercept of the fitted linear function to obtain the friction coefficient of the joint to be tested in the forward and reverse motions; calculate the average value of the motor torque coefficient and friction coefficient at different experimental positions to obtain the forward The target torque coefficient and target friction value of the joint to be tested during reverse motion. Therefore, under the condition of fully considering the influence of mechanical structure and working environment on key electromechanical parameters, the synchronous test of the torque coefficient of joint motor and the coefficient of forward and reverse friction force is realized, and the test accuracy is improved.

为了实现上述实施例,本申请还提出一种机器人的关节电机转矩系数与摩擦力同步测试装置。In order to realize the above embodiments, the present application also proposes a device for synchronizing the torque coefficient and friction force of a joint motor of a robot.

图12为本发明实施例提供的一种机器人的关节电机转矩系数与摩擦力同步测试装置的结构示意图。FIG. 12 is a schematic structural diagram of a device for synchronizing the torque coefficient and friction force of a joint motor of a robot according to an embodiment of the present invention.

如图12所示,该装置包括:建模模块1201、获取模块1202、第一计算模块1203、第二计算模块1204、第三计算模块1205。As shown in FIG. 12 , the apparatus includes: a modeling module 1201 , an acquisition module 1202 , a first calculation module 1203 , a second calculation module 1204 , and a third calculation module 1205 .

建模模块1201,用于基于虚功原理建立机器人的动力学模型,获取机器人的机械结构数,根据动力学模型和机械结构参数计算机器人的各驱动关节的驱动力矩解析表达式;The modeling module 1201 is used to establish a dynamic model of the robot based on the virtual work principle, obtain the mechanical structure number of the robot, and calculate the analytical expression of the driving torque of each driving joint of the robot according to the dynamic model and the mechanical structure parameters;

获取模块1202,用于获取预先存储的实验数据,其中,实验数据包括:各关节的角度、角速度、角加速度及驱动电流,以及按照角速度的方向将实验数据分别存储于两个数组中;The acquisition module 1202 is used to acquire pre-stored experimental data, wherein the experimental data includes: the angle, angular velocity, angular acceleration and driving current of each joint, and the experimental data is stored in two arrays respectively according to the direction of the angular velocity;

第一计算模块1203,用于将不同实验位置的两个数组中的实验数据分别输入驱动力矩解析表达式得到正向力矩和反向力矩,将正向力矩、反向力矩和驱动电流进行拟合处理,得到不同实验位置的拟合线性函数的斜率和截距;The first calculation module 1203 is used to input the experimental data in the two arrays of different experimental positions respectively into the driving torque analytical expression to obtain the forward torque and reverse torque, and fit the forward torque, reverse torque and driving current processing to obtain the slope and intercept of the fitted linear function at different experimental positions;

第二计算模块1204,用于根据拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将拟合线性函数的截距取相反数得到待测试关节正向与反向运动时的摩擦力系数;The second calculation module 1204 is configured to obtain the motor torque coefficients of the joint to be tested in forward and reverse motions according to the ratio of the slope of the fitted linear function to the reduction ratio of the joint to be tested, and to take the intercept of the fitted linear function as The opposite number obtains the friction coefficient of the joint to be tested in forward and reverse motions;

第三计算模块1205,用于计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。The third calculation module 1205 is used to calculate the average value of the motor torque coefficient and friction coefficient in different experimental positions, and obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion.

进一步地,在本申请实施例的一种可能的实现方式中,还包括:确定模块1206、实验模块1207、记录模块1208,其中,Further, in a possible implementation manner of the embodiment of the present application, it further includes: a determination module 1206, an experiment module 1207, and a recording module 1208, wherein,

确定模块1206,用于确定不同实验位置和待测试关节;a determination module 1206 for determining different experimental positions and joints to be tested;

实验模块1207,用于保持其他关节不动,在不同实验位置控制待测试关节做梯形加减速往复运动,其中,待测试关节匀速运动阶段的速度与待测试关节在实际工作中的角速度大小呈正相关关系。The experimental module 1207 is used to keep other joints still, and control the joint to be tested to perform trapezoidal acceleration and deceleration reciprocating motion at different experimental positions, wherein the speed of the joint to be tested in the uniform motion stage is positively correlated with the angular velocity of the joint to be tested in actual work. relation.

记录模块1208,用于获取在不同实验位置的实验数据。The recording module 1208 is used to acquire experimental data at different experimental locations.

进一步地,在本申请实施例的一种可能的实现方式中,记录模块1208,具体用于通过绝对值编码器获取待测试关节的角度;通过光栅编码器获取待测试关节的速度;通过对瞬时速度做差分获取待测试关节的角加速度;通过霍尔传感器获取待测试关节的驱动电流。Further, in a possible implementation of the embodiment of the present application, the recording module 1208 is specifically configured to obtain the angle of the joint to be tested through an absolute encoder; obtain the speed of the joint to be tested through a raster encoder; The angular acceleration of the joint to be tested is obtained by the speed difference; the driving current of the joint to be tested is obtained through the Hall sensor.

进一步地,在本申请实施例的一种可能的实现方式中,第一计算模块1203,用于计算驱动力矩解析表达式为:Further, in a possible implementation manner of the embodiment of the present application, the first calculation module 1203 is used to calculate the analytical expression of the driving torque as:

Figure GDA0003038459450000111
Figure GDA0003038459450000111

其中,τ为待测试关节实时的驱动力矩,

Figure GDA0003038459450000112
为驱动力矩中的惯性力分量,
Figure GDA0003038459450000113
为驱动力矩中的离心力与科氏力分量,G(θ)g为驱动力矩中的重力分量,另外,使用待测试关节的正向运动时的实验数据计算获得的为正向运动时的理论驱动力矩;使用待测试关节的反向运动时的实验数据计算获得的为反向运动时的理论驱动力矩。Among them, τ is the real-time driving torque of the joint to be tested,
Figure GDA0003038459450000112
is the inertial force component in the driving torque,
Figure GDA0003038459450000113
is the centrifugal force and Coriolis force components in the driving torque, G(θ)g is the gravity component in the driving torque, in addition, the theoretical driving in the forward movement is calculated using the experimental data during the forward movement of the joint to be tested. Torque; the theoretical driving torque during reverse movement is calculated using the experimental data of the joint to be tested during reverse movement.

需要说明的是,前述对方法实施例的解释说明也适用于该实施例的装置,此处不再赘述。It should be noted that the foregoing explanations of the method embodiment are also applicable to the apparatus of this embodiment, and details are not repeated here.

本发明实施例的机器人的关节电机转矩系数与摩擦力同步测试装置,建模模块基于虚功原理建立机器人的动力学模型,获取机器人的机械结构参数,根据动力学模型和机械结构参数计算机器人的各驱动关节的驱动力矩解析表达式;获取模块获取预先存储的实验数据,其中,实验数据包括:各待测试关节的角度、角速度、角加速度及驱动电流,以及按照角速度的方向将实验数据分别存储于两个数组中;第一计算模块将不同实验位置的两个数组中的实验数据分别输入驱动力矩解析表达式得到正向力矩和反向力矩,将正向力矩、反向力矩和驱动电流进行拟合处理,得到不同实验位置的拟合线性函数的斜率和截距;第二计算模块根据拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将拟合线性函数的截距取相反数得到待测试关节正向与反向运动时的摩擦力系数;第三计算模块计算不同实验位置的电机转矩系数和摩擦力系数的平均值,得到正向与反向运动时待测试关节的目标转矩系数与目标摩擦力数值。由此,在充分考虑机械结构与工作环境对关键机电参数影响的情况下,实现关节电机转矩系数与正反向摩擦力系数的同步测试,提高测试精度。The joint motor torque coefficient and friction force synchronization test device of the robot in the embodiment of the present invention, the modeling module establishes the dynamic model of the robot based on the virtual work principle, obtains the mechanical structure parameters of the robot, and calculates the robot according to the dynamic model and the mechanical structure parameters. The analytical expression of the driving torque of each driving joint; the acquisition module acquires the pre-stored experimental data, wherein the experimental data includes: the angle, angular velocity, angular acceleration and driving current of each joint to be tested, and the experimental data is divided according to the direction of the angular velocity. Stored in two arrays; the first calculation module inputs the experimental data in the two arrays at different experimental positions into the driving torque analytical expression respectively to obtain the forward torque and reverse torque, and the forward torque, reverse torque and drive current are calculated. Perform the fitting process to obtain the slope and intercept of the fitted linear function at different experimental positions; the second calculation module obtains the forward and reverse motions of the joint to be tested according to the ratio of the slope of the fitted linear function to the reduction ratio of the joint to be tested The friction coefficient of the joint to be tested in forward and reverse motions is obtained by taking the inverse of the intercept of the fitted linear function; the third calculation module calculates the motor torque coefficient and friction force at different experimental positions The average value of the coefficients is used to obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion. Therefore, under the condition of fully considering the influence of mechanical structure and working environment on key electromechanical parameters, the synchronous test of the torque coefficient of joint motor and the coefficient of forward and reverse friction force is realized, and the test accuracy is improved.

在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。In the description of this specification, description with reference to the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples", etc., mean specific features described in connection with the embodiment or example , structure, material or feature is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, those skilled in the art may combine and combine the different embodiments or examples described in this specification, as well as the features of the different embodiments or examples, without conflicting each other.

此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。In addition, the terms "first" and "second" are only used for descriptive purposes, and should not be construed as indicating or implying relative importance or implying the number of indicated technical features. Thus, a feature delimited with "first", "second" may expressly or implicitly include at least one of that feature. In the description of the present invention, "plurality" means at least two, such as two, three, etc., unless otherwise expressly and specifically defined.

流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。Any process or method description in the flowcharts or otherwise described herein may be understood to represent a module, segment or portion of code comprising one or more executable instructions for implementing custom logical functions or steps of the process , and the scope of the preferred embodiments of the invention includes alternative implementations in which the functions may be performed out of the order shown or discussed, including performing the functions substantially concurrently or in the reverse order depending upon the functions involved, which should It is understood by those skilled in the art to which the embodiments of the present invention belong.

在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得程序,然后将其存储在计算机存储器中。The logic and/or steps represented in flowcharts or otherwise described herein, for example, may be considered an ordered listing of executable instructions for implementing the logical functions, may be embodied in any computer-readable medium, For use with, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a system including a processor, or other system that can fetch instructions from and execute instructions from an instruction execution system, apparatus, or apparatus) or equipment. For the purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport the program for use by or in conjunction with an instruction execution system, apparatus, or apparatus. More specific examples (non-exhaustive list) of computer readable media include the following: electrical connections with one or more wiring (electronic devices), portable computer disk cartridges (magnetic devices), random access memory (RAM), Read Only Memory (ROM), Erasable Editable Read Only Memory (EPROM or Flash Memory), Fiber Optic Devices, and Portable Compact Disc Read Only Memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program may be printed, as may be done, for example, by optically scanning the paper or other medium, followed by editing, interpretation, or other suitable means as necessary process to obtain the program electronically and then store it in computer memory.

应当理解,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如,如果用硬件来实现和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。It should be understood that various parts of the present invention may be implemented in hardware, software, firmware or a combination thereof. In the above-described embodiments, various steps or methods may be implemented in software or firmware stored in memory and executed by a suitable instruction execution system. For example, if implemented in hardware as in another embodiment, it can be implemented by any one of the following techniques known in the art, or a combination thereof: discrete with logic gates for implementing logic functions on data signals Logic circuits, application specific integrated circuits with suitable combinational logic gates, Programmable Gate Arrays (PGA), Field Programmable Gate Arrays (FPGA), etc.

本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。Those of ordinary skill in the art can understand that all or part of the steps carried by the methods of the above embodiments can be completed by instructing the relevant hardware through a program, and the program can be stored in a computer-readable storage medium, and when the program is executed, Include one or a combination of the steps of a method embodiment.

此外,在本发明各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing module, or each unit may exist physically alone, or two or more units may be integrated into one module. The above-mentioned integrated modules can be implemented in the form of hardware, and can also be implemented in the form of software function modules. If the integrated modules are implemented in the form of software functional modules and sold or used as independent products, they may also be stored in a computer-readable storage medium.

上述提到的存储介质可以是只读存储器,磁盘或光盘等。尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, and the like. Although the embodiments of the present invention have been shown and described above, it should be understood that the above-mentioned embodiments are exemplary and should not be construed as limiting the present invention. Embodiments are subject to variations, modifications, substitutions and variations.

Claims (10)

1.一种机器人的关节电机转矩系数与摩擦力同步测试方法,其特征在于,包括:1. the joint motor torque coefficient of a robot and the friction force synchronization test method, it is characterized in that, comprise: 基于虚功原理建立机器人的动力学模型,获取机器人的机械结构参数,根据所述动力学模型和所述机械结构参数计算所述机器人的各驱动关节的驱动力矩解析表达式;Establish a dynamic model of the robot based on the principle of virtual work, obtain the mechanical structure parameters of the robot, and calculate the analytical expression of the driving torque of each driving joint of the robot according to the dynamic model and the mechanical structure parameters; 获取预先存储的实验数据,其中,所述实验数据包括:各待测试关节的角度、角速度、角加速度及驱动电流,以及按照所述角速度的方向将所述实验数据分别存储于两个数组中;Acquiring pre-stored experimental data, wherein the experimental data includes: the angle, angular velocity, angular acceleration and driving current of each joint to be tested, and the experimental data is stored in two arrays respectively according to the direction of the angular velocity; 将不同实验位置的所述两个数组中的实验数据分别输入所述驱动力矩解析表达式得到正向力矩和反向力矩,将所述正向力矩、所述反向力矩和所述驱动电流进行拟合处理,得到所述不同实验位置的拟合线性函数的斜率和截距;Input the experimental data in the two arrays at different experimental positions into the driving torque analytical expression to obtain the forward torque and the reverse torque, and carry out the forward torque, the reverse torque and the driving current. Fitting processing, obtaining the slope and intercept of the fitting linear function of described different experimental positions; 根据所述拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将所述拟合线性函数的截距取相反数得到所述待测试关节正向与反向运动时的摩擦力系数;According to the ratio of the slope of the fitted linear function and the reduction ratio of the joint to be tested, the motor torque coefficient of the joint to be tested is obtained when the forward and reverse motions are performed, and the intercept of the fitted linear function is taken to be the opposite number to obtain the Describe the friction coefficient when the joint to be tested moves forward and backward; 计算所述不同实验位置的所述电机转矩系数和所述摩擦力系数的平均值,得到正向与反向运动时所述待测试关节的目标转矩系数与目标摩擦力数值。Calculate the average value of the motor torque coefficient and the friction coefficient at the different experimental positions, and obtain the target torque coefficient and target friction value of the joint to be tested during forward and reverse motion. 2.如权利要求1所述的方法,其特征在于,在所述获取预先存储的实验数据之前,还包括:2. The method according to claim 1, characterized in that, before said acquiring the pre-stored experimental data, further comprising: 确定所述不同实验位置和所述待测试关节;determining the different experimental positions and the joint to be tested; 保持其他关节不动,在所述不同实验位置控制所述待测试关节做梯形加减速往复运动,其中,所述待测试关节匀速运动阶段的速度与所述待测试关节在实际工作中的角速度大小呈正相关关系;Keep other joints still, and control the joint to be tested to do trapezoidal acceleration and deceleration reciprocating motion at the different experimental positions, wherein the speed of the joint to be tested in the uniform motion stage is the same as the angular velocity of the joint to be tested in actual work. a positive correlation; 获取在所述不同实验位置的所述实验数据。Acquire the experimental data at the different experimental locations. 3.如权利要求2所述的方法,其特征在于,所述获取在所述不同实验位置的所述实验数据,包括:3. The method of claim 2, wherein the acquiring the experimental data at the different experimental positions comprises: 通过绝对值编码器获取所述待测试关节的所述角度;Obtain the angle of the joint to be tested through an absolute value encoder; 通过光栅编码器获取所述待测试关节的所述速度;Obtain the speed of the joint to be tested by using a raster encoder; 通过对瞬时速度做差分获取所述待测试关节的所述角加速度;Obtain the angular acceleration of the joint to be tested by differentiating the instantaneous velocity; 通过霍尔传感器获取所述待测试关节的所述驱动电流。The drive current of the joint to be tested is acquired through a Hall sensor. 4.如权利要求1所述的方法,其特征在于,其中,所述驱动力矩解析表达式为:4. The method according to claim 1, wherein the analytical expression of the driving torque is:
Figure FDA0003038459440000024
Figure FDA0003038459440000024
其中,τ为所述待测试关节实时的驱动力矩,
Figure FDA0003038459440000022
为所述驱动力矩中的惯性力分量,
Figure FDA0003038459440000023
为所述驱动力矩中的离心力与科氏力分量,G(θ)g为所述驱动力矩中的重力分量,另外,使用所述待测试关节的正向运动时的所述实验数据计算获得的为正向运动时的理论驱动力矩;使用所述待测试关节的反向运动时的所述实验数据计算获得的为反向运动时的理论驱动力矩。
Wherein, τ is the real-time driving torque of the joint to be tested,
Figure FDA0003038459440000022
is the inertial force component in the driving torque,
Figure FDA0003038459440000023
is the centrifugal force and the Coriolis force component in the driving torque, G(θ)g is the gravity component in the driving torque, and is calculated using the experimental data of the joint to be tested in the forward motion. is the theoretical driving torque during forward motion; the theoretical driving torque during reverse motion is calculated and obtained by using the experimental data during reverse motion of the joint to be tested.
5.如权利要求1所述的方法,其特征在于,进行拟合使用的公式为:5. The method of claim 1, wherein the formula used for fitting is:
Figure FDA0003038459440000021
Figure FDA0003038459440000021
其中,Kt为所述待测试关节电机的转矩系数,Iq为所述待测试关节电机的所述驱动电流,f为所述待测试关节运行所受到的摩擦力,i为所述待测试关节的减速比,τ为所述待测试关节的理论驱动力矩,包含所述待测试关节的所述正向运动时的理论驱动力矩和所述反向运动时的理论驱动力矩,当使用所述待测试关节的所述正向运动时的理论驱动力矩时,拟合获得的所述函数为正向运动拟合函数;当使用所述待测试关节的所述反向运动时的理论驱动力矩时,拟合获得的所述函数为反向运动拟合函数。Wherein, K t is the torque coefficient of the joint motor to be tested, I q is the drive current of the joint motor to be tested, f is the friction force experienced by the joint to be tested, and i is the friction force of the joint motor to be tested. The reduction ratio of the test joint, τ is the theoretical driving torque of the joint to be tested, including the theoretical driving torque of the joint to be tested during the forward motion and the theoretical driving torque of the reverse motion. When using all When using the theoretical driving torque of the joint to be tested in the forward motion, the function obtained by fitting is a forward motion fitting function; when using the theoretical driving torque of the joint to be tested in the reverse motion When , the function obtained by fitting is the inverse kinematic fitting function.
6.如权利要求2所述的方法,其特征在于,在所述不同实验位置控制所述待测试关节做梯形加减速往复运动,包括:6. The method of claim 2, wherein the joint to be tested is controlled to do trapezoidal acceleration and deceleration reciprocating motion at the different experimental positions, comprising: 在所述机器人的预设工作空间内,以间距相等的方式均匀地选取若干个所述机器人末端执行器的坐标定位生成所述不同实验位置;In the preset working space of the robot, evenly select the coordinate positioning of several end effectors of the robot with equal spacing to generate the different experimental positions; 在所述不同实验位置,控制所述待测试关节实施相同运动行程、相同速度规划的运动。At the different experimental positions, the joints to be tested are controlled to implement the same motion stroke and the same speed planned motion. 7.一种机器人的关节电机转矩系数与摩擦力同步测试装置,包括:7. A joint motor torque coefficient and friction force synchronization test device of a robot, comprising: 建模模块,用于基于虚功原理建立机器人的动力学模型,获取所述机器人的机械结构参数,根据所述动力学模型和所述机械结构参数计算所述机器人的各驱动关节的驱动力矩解析表达式;The modeling module is used to establish a dynamic model of the robot based on the virtual work principle, obtain the mechanical structure parameters of the robot, and calculate the driving torque analysis of each driving joint of the robot according to the dynamic model and the mechanical structure parameters. expression; 获取模块,用于获取预先存储的实验数据,其中,所述实验数据包括:各待测试关节的角度、角速度、角加速度及驱动电流,以及按照所述角速度的方向将所述实验数据分别存储于两个数组中;The acquisition module is used to acquire pre-stored experimental data, wherein the experimental data includes: the angle, angular velocity, angular acceleration and driving current of each joint to be tested, and the experimental data are respectively stored in the direction of the angular velocity according to the direction of the angular velocity. in two arrays; 第一计算模块,用于将不同实验位置的所述两个数组中的实验数据分别输入所述驱动力矩解析表达式得到正向力矩和反向力矩,将所述正向力矩、所述反向力矩和所述驱动电流进行拟合处理,得到所述不同实验位置的拟合线性函数的斜率和截距;The first calculation module is used to input the experimental data in the two arrays of different experimental positions respectively into the driving torque analytical expression to obtain a forward torque and a reverse torque, and the forward torque and the reverse torque are calculated. The torque and the drive current are fitted to obtain the slope and intercept of the fitted linear function at the different experimental positions; 第二计算模块,用于根据所述拟合线性函数的斜率和待测试关节的减速比的比值得到待测试关节正向与反向运动时的电机转矩系数,将所述拟合线性函数的截距取相反数得到所述待测试关节正向与反向运动时的摩擦力系数;The second calculation module is configured to obtain the motor torque coefficient of the joint to be tested in forward and reverse motions according to the ratio of the slope of the fitted linear function to the reduction ratio of the joint to be tested, The intercept takes the opposite number to obtain the friction coefficient when the joint to be tested moves forward and backward; 第三计算模块,用于计算所述不同实验位置的所述电机转矩系数和所述摩擦力系数的平均值,得到正向与反向运动时所述待测试关节的目标转矩系数与目标摩擦力数值。The third calculation module is used to calculate the average value of the motor torque coefficient and the friction coefficient of the different experimental positions, and obtain the target torque coefficient and target torque coefficient of the joint to be tested during forward and reverse motion. Friction value. 8.如权利要求7所述的装置,其特征在于,还包括:8. The apparatus of claim 7, further comprising: 确定模块,用于确定所述不同实验位置和所述待测试关节;a determination module for determining the different experimental positions and the joint to be tested; 实验模块,用于保持其他关节不动,在所述不同实验位置控制所述待测试关节做梯形加减速往复运动,其中,所述待测试关节匀速运动阶段的速度与所述待测试关节在实际工作中的角速度大小呈正相关关系;The experimental module is used to keep other joints still, and control the joint to be tested to perform trapezoidal acceleration and deceleration reciprocating motion at the different experimental positions, wherein the speed of the joint to be tested in the uniform motion stage is the same as that of the joint to be tested. The angular velocity during work is positively correlated; 记录模块,用于获取在所述不同实验位置的所述实验数据。a recording module for acquiring the experimental data at the different experimental positions. 9.如权利要求8所述的装置,所述记录模块,具体用于:9. The device according to claim 8, the recording module is specifically used for: 通过绝对值编码器获取所述待测试关节的所述角度;Obtain the angle of the joint to be tested through an absolute value encoder; 通过光栅编码器获取所述待测试关节的所述速度;Obtain the speed of the joint to be tested by using a raster encoder; 通过对瞬时速度做差分获取所述待测试关节的所述角加速度;Obtaining the angular acceleration of the joint to be tested by differentiating the instantaneous velocity; 通过霍尔传感器获取所述待测试关节的所述驱动电流。The drive current of the joint to be tested is acquired through a Hall sensor. 10.如权利要求9所述的装置,所述第一计算模块,用于计算所述驱动力矩解析表达式为:10. The device according to claim 9, wherein the first calculation module is used to calculate the analytical expression of the driving torque as:
Figure FDA0003038459440000031
Figure FDA0003038459440000031
其中,τ为所述待测试关节实时的驱动力矩,
Figure FDA0003038459440000032
为所述驱动力矩中的惯性力分量,
Figure FDA0003038459440000033
为所述驱动力矩中的离心力与科氏力分量,G(θ)g为所述驱动力矩中的重力分量,另外,使用所述待测试关节的正向运动时的所述实验数据计算获得的为正向运动时的理论驱动力矩;使用所述待测试关节的反向运动时的所述实验数据计算获得的为反向运动时的理论驱动力矩。
Wherein, τ is the real-time driving torque of the joint to be tested,
Figure FDA0003038459440000032
is the inertial force component in the driving torque,
Figure FDA0003038459440000033
is the centrifugal force and the Coriolis force component in the driving torque, G(θ)g is the gravity component in the driving torque, and is calculated using the experimental data of the joint to be tested in the forward motion. is the theoretical driving torque during forward motion; the theoretical driving torque during reverse motion is calculated and obtained by using the experimental data during reverse motion of the joint to be tested.
CN202010843773.1A 2020-08-20 2020-08-20 Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot Active CN112179551B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010843773.1A CN112179551B (en) 2020-08-20 2020-08-20 Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010843773.1A CN112179551B (en) 2020-08-20 2020-08-20 Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot

Publications (2)

Publication Number Publication Date
CN112179551A CN112179551A (en) 2021-01-05
CN112179551B true CN112179551B (en) 2021-06-22

Family

ID=73924148

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010843773.1A Active CN112179551B (en) 2020-08-20 2020-08-20 Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot

Country Status (1)

Country Link
CN (1) CN112179551B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114291051B (en) * 2022-01-05 2023-03-28 清华大学 Method and device for modeling internal friction force of motor servo type hydraulic line control brake system
WO2023245937A1 (en) * 2022-06-21 2023-12-28 郑州安杰莱智能科技有限公司 Conversion method and apparatus for current moment, and rehabilitation robot
CN116442220A (en) * 2023-03-30 2023-07-18 之江实验室 Parameter identification method and device for robot joint friction model and moment estimation method and device

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102269638A (en) * 2011-04-27 2011-12-07 中国科学院光电技术研究所 Integrated Measurement Method of LuGre Model Friction Parameters and Moment of Inertia of Servo Turntable
CN104011613A (en) * 2011-08-04 2014-08-27 奥尔德巴伦机器人公司 Robot with joints of variable rigidity and method for calculating said optimised rigidity
CN106994686A (en) * 2016-12-01 2017-08-01 遨博(北京)智能科技有限公司 The computational methods and device, robot of joint external force torque
CN107520258A (en) * 2016-06-22 2017-12-29 上海梅山钢铁股份有限公司 The method of testing of cold rolling mill motor losses moment coefficient
JP6272203B2 (en) * 2014-09-30 2018-01-31 オートリブ日信ブレーキシステムジャパン株式会社 Vehicle control device
CN108538381A (en) * 2018-05-10 2018-09-14 中国科学院宁波材料技术与工程研究所 The reciprocal force detection method and system of wrist device for healing and training
CN109773794A (en) * 2019-02-26 2019-05-21 浙江大学 A Neural Network-based Dynamic Parameter Identification Method for 6-axis Robot
CN110355780A (en) * 2018-04-09 2019-10-22 中国科学院宁波材料技术与工程研究所 The output torque calculation method of integrated flexible joint of robot
CN111422247A (en) * 2020-03-20 2020-07-17 合肥工业大学 Fault diagnosis and fault tolerance compensation algorithm for steer-by-wire motor and sensor
CN111438692A (en) * 2020-04-17 2020-07-24 达闼机器人有限公司 Robot control method, device, medium, equipment and robot

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4281696B2 (en) * 2005-03-14 2009-06-17 セイコーエプソン株式会社 Acceleration / deceleration control method and apparatus, and acceleration / deceleration control method program

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102269638A (en) * 2011-04-27 2011-12-07 中国科学院光电技术研究所 Integrated Measurement Method of LuGre Model Friction Parameters and Moment of Inertia of Servo Turntable
CN104011613A (en) * 2011-08-04 2014-08-27 奥尔德巴伦机器人公司 Robot with joints of variable rigidity and method for calculating said optimised rigidity
JP6272203B2 (en) * 2014-09-30 2018-01-31 オートリブ日信ブレーキシステムジャパン株式会社 Vehicle control device
CN107520258A (en) * 2016-06-22 2017-12-29 上海梅山钢铁股份有限公司 The method of testing of cold rolling mill motor losses moment coefficient
CN106994686A (en) * 2016-12-01 2017-08-01 遨博(北京)智能科技有限公司 The computational methods and device, robot of joint external force torque
CN110355780A (en) * 2018-04-09 2019-10-22 中国科学院宁波材料技术与工程研究所 The output torque calculation method of integrated flexible joint of robot
CN108538381A (en) * 2018-05-10 2018-09-14 中国科学院宁波材料技术与工程研究所 The reciprocal force detection method and system of wrist device for healing and training
CN109773794A (en) * 2019-02-26 2019-05-21 浙江大学 A Neural Network-based Dynamic Parameter Identification Method for 6-axis Robot
CN111422247A (en) * 2020-03-20 2020-07-17 合肥工业大学 Fault diagnosis and fault tolerance compensation algorithm for steer-by-wire motor and sensor
CN111438692A (en) * 2020-04-17 2020-07-24 达闼机器人有限公司 Robot control method, device, medium, equipment and robot

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
An Experimental Study on the Dynamics;Dong Wang;《IEEE》;20191231;1-11 *

Also Published As

Publication number Publication date
CN112179551A (en) 2021-01-05

Similar Documents

Publication Publication Date Title
CN112179551B (en) Method and device for synchronizing test method and device of joint motor torque coefficient and friction force of robot
KR102601888B1 (en) Vibration control of configuration-dependent dynamical systems
US10486309B2 (en) Robot controlling method, robot apparatus, program, recording medium, and method for manufacturing assembly component
JP6895242B2 (en) Robot control device, robot control method and picking device
JP5727103B2 (en) Robot apparatus and control method thereof
Duchaine et al. Computationally efficient predictive robot control
CN110076773B (en) Robot apparatus, method of manufacturing article, and motor driving apparatus
Huang et al. Modeling and elastic deformation compensation of flexural feed drive system
CN112799304B (en) A dual-motor synchronous control method and device based on time-varying friction compensation
CN109048901B (en) Traction and teaching trajectory planning method and device based on virtual friction force
JP2016083713A (en) Robot control method, robot device, program, recording medium and assembly part manufacturing method
CN109434873A (en) The measurement method of joint of robot servo motor torque constant
CN106926233B (en) A kind of planing method of planar manipulator motion path
CN112231888B (en) Ideal dynamics model establishing method based on laser tracking measurement system mechanical structure
JP2007272597A (en) Method for generating robot acceleration and deceleration pattern
Eppinger Modeling robot dynamic performance for endpoint force control
JP7158636B1 (en) Machining evaluation device, machining system, and machining evaluation method
CN113910247A (en) End-trajectory control method and system for industrial robot considering joint gap evolution
Li et al. Research on dual-driving synchronous control system modeling of gantry-type machine tools with traveling bridge
WO2020217597A1 (en) Servo control device
Saied et al. Actuator and friction dynamics formulation in control of PKMs: From design to real-time experiments
Du et al. Current-Based direct teaching for industrial manipulator
Bellakehal et al. Force/position control of parallel robots using exteroceptive pose measurements
JP4449693B2 (en) Robot control apparatus and control method thereof
CN109240269B (en) A Dynamic Performance Analysis Method for Parallel Mechanisms

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant