Patent Application
20220105628
The present invention relates to robot joint gears for robot arms comprising a plurality of robot joints connecting a robot base and a robot tool flange.
Robot arms comprising a plurality of robot joints and links where motors can rotate the joints in relation to each other are known in the field of robotics. Typically, the robot arm comprises a robot base which serves as a mounting base for the robot arm and a robot tool flange where to various tools can be attached. A robot controller is configured to control the robot joints to move the robot tool flange in relation to the base. For instance, in order to instruct the robot arm to carry out a number of working instructions.
Typically, the robot controller is configured to control the robot joints based on a dynamic model of the robot arm, where the dynamic model defines a relationship between the forces acting on the robot arm and the resulting accelerations of the robot arm. Often, the dynamic model comprises a kinematic model of the robot arm, knowledge about inertia of the robot arm and other parameters influencing the movements of the robot arm. The kinematic model defines a geometric relationship between the different parts of the robot arm and may comprise information of the robot arm such as, length, size of the joints and links and can for instance be described by Denavit-Hartenberg parameters or the like. The dynamic model makes it possible for the controller to determine which torques the joint motors shall provide in order to move the robot joints for instance at specified velocity, acceleration or in order to hold the robot arm in a static posture.
On many robot arms it is possible to attach various end effectors to the robot tool flange, such as grippers, vacuum grippers, magnetic grippers, screwing machines, welding equipment, dispensing systems, visual systems etc.
In some robots the robot joint comprises a joint motor having a motor axle configured to rotate an output axle via a robot joint gear. Typically, the output axle is connected to and configured to rotate parts of the robot arm in relation to each other. The robot joint gear forms a transmission system configured to transmit torque provided by the motor axle to the output axle for instance to provide a gear ratio between the motor axle and the output axle. The robot joint gear can for instance be provided as a spur gears, planetary gears, bevel gears, worm gears, strain wave gears or other kind of transmission systems. Commonly flexibility exist in the transmissions used for industrial robots due to the elasticity various type of transmissions. This flexibility may lead to an undesired dynamic time-varying displacement between the position of parts of the robot arm.
Taking into account the gear flexibility in the dynamic model makes the dynamic model more accurately resemble the dynamics of the real robot arm because the robot joint torque originating from the gear deformation can be known. A more accurate dynamic model can for instance allow the robot controller to control the robot arm with greater accuracy and precision. A more accurate dynamic model can also allow the robot controller to more accurately identify external disturbances, for instance human interference which is of great concern in terms of safety.
Research have been devoted to accurately identify the dynamic characteristics of the robot gear systems. However, one practical issue when taking into account the joint stiffness in the robot controller design is that the joint stiffness changes with wear due to material being worn off at the gear meshing [1], [2], [3].
Estimating the joint stiffness on industrial robots have been accomplished through off-line identification procedures by several researchers. Off-line identification procedures are not well suited for solving the problem of time-varying joint stiffness. Such procedures would have to be re-run from time to time to keep the joint stiffness information up to date. While the off-line identification procedure is running the robot is unable to conduct any other task with clear negative consequences to the user. Despite the shortcomings of the off-line calibration procedures for joint stiffness estimation a number of references will be given to such off-line identification procedures applied to solve the robot joint stiffness estimation problem. For the off-line identification of joint stiffness, one method is to apply external excitation on one robot joint at a time resulting in a reduced dynamic model hence easier identification as in [4], [5], [6], [7].
Another option is to evaluate joint stiffness values using an external laser-tracking sensor system to visually track the end effector displacements for a given applied wrench as in [8] and [9], where the applied wrench is equal to the applied forced and torques at the end effector. A third option is to incorporate two absolute position rotary encoders in the robot joint to enable direct measurement of the joint deflection and then manually impose a known force at the end-effector and using the kinematic model (the Jacobian) of the robot arm to calculate the joint torques as in [9]. Such procedure however requires external and known loading of the end-effector and the robot is not allowed to move during the procedure. A fourth option is to estimate the joint stiffness using a so-called locked-link joint procedure as in [10] where the end-effector is clamped to the environment and using either motor positions measurement and motor torques data or a force/torque sensor to measure the external wrench between the clamped end effector and the environment.
Very few works are available on the on-line estimation of time-varying nonlinear stiffness e.g. AwAs-II joint [11]. These works are focused towards Variable Stiffness Actuators (VSAs) as found in in the DLR Hand Arm System [12] by the DLR Institute of Robotics and Mechatronics, however most industrial robots do not comprise a VSA. In 2011, Flacco and De Luca [13] estimated the nonlinear stiffness of robot joints using only a motor position sensor by computing a dynamic residual based on the generalized momentum followed by a least squares algorithm to estimate the stiffness parameters.
Robustness issues were later addressed in [14] by introducing a kinematic Kalman filter to handle discretization and quantization errors and a modified recursive least squares algorithm was used to better handle poor excitation conditions. The stiffness parameters were however assumed time-invariant. Further refinements of their method were carried out in [15] modifying the stiffness estimation algorithm to deal with time-varying stiffness by using a Recursive Least Squares method based on a QR decomposition (QR-RLS) able to handle time-varying stiffness characteristics. A non-causal Savitzky-Golay (SG) filter was used to remove noise on the input/output signals. No experiments did however support their findings.
In May 2013, Menard et al. [16] developed an observer capable of on-line estimating the time-varying stiffness of a VSA. Experimental analysis on the VSA system revealed parameter uncertainties of up to 25% of the true stiffness. Friction was assumed purely viscous with a single constant coefficient. Such friction model is not capable of accurately describing real frictional characteristics of most electromechanical systems, as well-known frictional characteristics such as Coulomb friction and the Stribeck effect is then neglected completely.
In July 2013, Cirillo et al. [17] demonstrated online stiffness estimation of a VSA by measuring the elastic energy using an optoelectronic sensor built into the robot joint.
The objective of the present invention is to address the above described limitations with the prior art or other problems of the prior art. This is achieved by a method of obtaining the gear stiffness of a robot joint gear as defined by the independent claim. The dependent claims describe possible embodiments of the method according to the present invention. The advantages and benefits of the present invention are described in the detailed description of the invention. Further the objective of the present invention is addressed by a method of controlling a robot arm based on the obtained gear stiffness and a robot arm with a controller configured to control the robot arm based on the obtained gear stiffness.
A robot controller taking into account the joint flexibilities may perform well right after calibration but for a good performance over the lifetime of the robot, the joint stiffness must be estimated on-line. Further, estimating the joint stiffness on-line will allow for predictive maintenance in the case of gear unit failure. In most cases, no sensor is available for directly measuring the joint stiffness, so the stiffness information is collected by combining an accurate model of the flexibility torque with the measurements available, such as position and force/torque sensor measures, either implemented in the robot or by using external hardware.
The invention provides a simple method of obtaining the gear stiffness of a robot joint gear of a robot arm without the need to integrate expensive force/torque sensors, as the joint gear stiffness of the robot joint gear can be determined based on sensors commonly used in industrial robots, which also makes it possible to integrate such method into existing industrial robots having such sensors. Additionally, the gear stiffness can be obtained online during use of the robot and thus used as an input to the dynamic model controlling the robot whereby a more accurate controlling of the robot can be provided. Consequently, vibrations and inaccuracies of the robot's movements due to gear stiffness can be accounted for when controlling the robot. Additionally, the gear stiffness changes over time primarily due to wear of the robot joint gear and obtaining the gear stiffness online makes it possible to control the robots based on changes of gear stiffness whereby the accuracy of the robot over time is maintained/improved. Further obtaining the gear stiffness online over time makes it possible to predict failure of the robot joint gear as the gear stiffness can be used to indicate when the robot joint gear is about to fail. Consequently, the robot joint gear can be sent to service in order to exchange/repair the robot joint gear whereby unplanned downtime of the robot can be reduced.
The present invention is described in view of exemplary embodiments only intended to illustrate the principles of the present invention. The skilled person will be able to provide several embodiments within the scope of the claims. Throughout the description, the reference numbers of similar elements providing similar effects have the same last two digits. Further it is to be understood that in the case that an embodiment comprises a plurality of the same features then only some of the features may be labeled by a reference number.
The invention can be embodied into a robot arm and is described in view of the robot arm illustrated in
The robot arm comprises at least one robot controller 115 configured to control the robot joints by controlling the motor torque provided to the joint motors based on a dynamic model of the robot. The robot controller 115 can be provided as a computer comprising an interface device 117 enabling a user to control and program the robot arm. The controller can be provided as an external device as illustrated in
The joint motor is configured to rotate the motor axle by applying a motor torque to the motor axle as known in the art of motor control, for instance based on a motor control signal 233 indicating the torque, τcontrol, motor, applied by said motor axle.
The robot joint gear 229 forms a transmission system configured to transmit the torque provided by the motor axle to the output axle for instance to provide a gear ratio between the motor axle and the output axle. The robot joint gear can for instance be provided as spur gears, planetary gears, bevel gears, worm gears, strain wave gears or other kind of transmission systems. The robot joint comprises at least one joint sensor providing a sensor signal indicative of at least the angular position, q, of the output axle and an angular position, Θ, of the motor axle. For instance, the angular position of the output axle can be indicated by an output encoder 235, which provide an output encoder signal 236 indicating the angular position of the output axle in relation to the robot joint. Similarly, the angular position of the motor axle can be provided by an input encoder 237 providing an input encoder signal 238 indicating the angular position of the motor axle in relation to the robot joint. The output encoder 235 and the input encoder 237 can be any encoder capable of indicating the angular position, velocity and/or acceleration of respectively the output axle and the motor axle. The output/input encoders can for instance be configured to obtain the position of the respective axle based on the position of an encoder wheel 239 arrange on the respective axle. The encoder wheels can for instance be optical or magnetic encoder wheels as known in the art of rotary encoders. The output encoder indicating the angular position of the output axle and the input encoder indicating the angular position of the motor axle makes it possible to determine a relationship between the input side (motor axle) and the output side (output axle) of the robot joint gear.
The robot joints may optionally comprise one or more motor torque sensors 241 providing a motor torque signal 242 indicating the torque provided by the motor axle. For instance, the motor torque sensor can be provided as current sensors obtaining the current through the coils of the joint motor whereby the motor torque can be determined as known in the art of motor control. For instance, in connection with a multiphase motor, a plurality of current sensors can be provided in order to obtain the current through each of the phases of the multiphase motor and the motor torque can then be obtained based on the quadrature current obtained from the phase currents through a Park Transformation. Alternatively, the motor torque can be obtained using other kind of sensors for instance force-torque sensors, strain gauges etc.
The non-infinite gear stiffness of the robot joint results in a deflection between the input side and the output side of the robot joint gear when a torque is applied to the robot joint gear. The deflection of the robot joint gear can be indicated by a deflection variable indicating the differences between the angular position Θ of the motor axle and the angular position q of the of the output axle, thus the deflection variable is defined as:
Φjoint,i=Θi−qi eq. 1
The joint transmission torque τjoint,i defines the torque that is transferred from the motor axle to the output axle via the robot joint gear and can be modeled as a function of the deflection variable Φjoint and its time-derivative:
τjoint,i(Φjoint,i,{dot over (Φ)}joint,i)=τE,i(Φjoint,i)+τD,i({dot over (Φ)}joint,i) eq. 2
where τE,i(Φjoint,i) is a flexibility torque depending on the gear stiffness Ki and the deflection of the robot joint gear and τD,i({dot over (Φ)}joint,i) is damping torque depending on the damping coefficient Di and the first time derivative of the deflection of the robot joint gear and the time derivative of the deflection of the robot joint gear.
The τD,i({dot over (Φ)}joint,i) damping torque of the robot joint gear can for instance be obtained by following the steps:
The gear stiffness Ki of the robot joint gear can be characterized by how much the flexible torque τE causing the gear deflection changes as a function of the gear deflection. The gear stiffness can thus be expressed as:
Assuming that that no flexibility torque exists for the undeformed transmission and that the transmission has the same behavior in compression and extension, thus:
τE,i(0)=0 ∀Φjoint,i eq. 5
τE,i(−Φjoint,i)=−τE(Φjoint,i) ∀Φjoint,i eq. 6
The controller is configured to control the joint motors of the robot joints by providing motor control signals to the joint motors. The motor control signals 433i, 433i+1 . . . 433n are indicative of the motor torque τcontrol,motor,i, τcontrol,motor,i+1, and τcontrol,motor,n that each joint motor shall provide by the motor axles. The motor control signals can indicate the desired motor torque, the desired torque provided by the output axle, the currents provided by the motor coils or any other signal from which the motor torque can be obtained. The motor torque signals can be sent to a motor control driver (not shown) configured to drive the motor joint with the motor current resulting in the desired motor torque. The robot controller is configured to determine the motor torque based on a dynamic model of the robot arm as known in the prior art. The dynamic model makes it possible for the controller to calculate which torque the joint motors shall provide to each of the joint motors to make the robot arm perform a desired movement and/or be arranged in a static posture. The dynamic model of the robot arm can be stored in the memory 445.
As described in connection with
The controller is further configured to obtain the gear stiffness of at least one of the robot joint gears of the robot joints of the robot arm by:
The controller can for instance be configured to carry out the method of obtaining the gear stiffness of a robot joint gear obtained by the method illustrated in
The method comprises a step of initializing 550, a step 552 of applying a motor torque to the motor axles of the joint motors; a step 554 of obtaining the angular position of the motor axles of the joint motors; a step 556 of obtaining the angular position of the output axle of the robot joint gears and a step 560 of determining the gear stiffness of the robot joint gears based on the motor torques applied to the motor axles, the obtained angular positions of the motor axles, the obtained angular positions of the output axles and the dynamic model of the robot arm.
Step of initializing 550 comprises a step of obtaining the dynamic model Drobot of the robot arm and can be based on prior knowledge of the robot arm and robot joints, KoR [Knowledge of Robot], such as the dimensions and weight of robot joints and robot links; joint motor properties; information relating to an eventual payload attached to the robot arm, orientation of the robot arm in relation to gravity and frictional properties of the robot arm and robot joints.
The dynamic model of the robot arm can be defined and pre-stored in the memory of the controller and the user can in some embodiment be allowed to modify the dynamic model op the robot arm, for instance by providing payload information of a payload attached to the robot arm or defining the orientation of the robot arm in relation to gravity.
The dynamic model of the robot arm can be obtained by considering the robot arm as an open kinematic chain having a plurality of (n+1) rigid robot links and a plurality of n revolute robot joints, comprising a joint motor configured to rotate at least one robot link. In this example the dynamic model is provided by making the following assumptions:
A1: The rotors of the joint motors are uniform bodies having their center of mass at the axis of rotation (motor axle);
A2: The joint motors are located on the robot links preceding the driven robot links; and
A3: The angular velocity of the rotors of the joint motors are due to their own spinning.
Assumption A1 is a basic requirement for long life of a joint motor such as electric motors and also implies that the dynamics of the robot joint will be independent of the angular position of the motor axle. The kinematic arrangement of joint motors and robot links described in assumption A2 is illustrated in
The configuration of robot arm can be characterized by the generalized coordinates (q Θ)∈2N where q is a vector comprising the angular position of the output axles of the robot joint gears and Θ is a vector comprising the angular position of the motor axles as seen in “space” of the output side of the robot joint gear. Consequently:
where Θreal is the real angular position of the motor axle (e.g. as measured by an encoder) and r the gear ratio of the robot joint gear. This is the notation used throughout this application.
Similar τmotor is a vector comprising the torques of the motor axles 4 as seen in “space” of the output side of the robot joint gear. Consequently:
τmotor=r·τmotor real
where τmotor real is the real torque of the motor axles (e.g. as measured by sensors) and r the gear ratio of the robot joint gear. This is the notation used throughout this application.
Following assumptions 1-3, firstly the dynamic model as seen from the output side of the robot joint gears of the robot arm becomes [18]:
τjoint=M(q){umlaut over (q)}+C(q,{dot over (q)}){dot over (q)}+G(q)+Fq({dot over (q)})+τext eq. 7
where τjoint is a vector comprising the transmission torques τjoint,i . . . τjoint,n of each of the robot joint gears; q is a vector comprising the angular position of the output axles of the robot joint gears; {dot over (q)} is a vector comprising the first time derivative of the angular position of the output axles of the robot joint gears and thus relates to the angular velocity of the output axles; q is a vector comprising the second time derivative of the angular position of the output axles of the robot joint gears and thus relates to the angular acceleration of the output axles. M(q) is the inertia matrix of the robot arm and indicates the mass moments of inertia of the robot arm as a function of the angular position of the output axles of the robot joint gears. C(q,{dot over (q)}){dot over (q)} is the Coriolis and centripetal torques of the robot arm as a function of the angular position and angular velocity of the output axle of the robot joint gears. G(q) is the gravity torques acting on the robot arm as a function of the angular position of the output axles of the robot joint gears.
Fq({dot over (q)}) is a vector comprising the friction torques acting on the output axles of the robot joint gears. The friction torques acting on the output axle depends on angular velocity of the output axle ({dot over (q)}); however it is to be understood the friction torques acting on the output axle also can depend on other parameters such as temperatures, type of lubricants, loads to robot arm, position/orientation of the robot arm etc. Fq({dot over (q)}) can for instance be provided as linear or nonlinear functions or lookup tables (LUTs) with interpolation, and Fq({dot over (q)}) can be defined based on for instance prior knowledge of the robot, experiments, and/or be adaptively updated during robot operation. For instance, the Fq({dot over (q)}) can be obtained during calibration of the robot joints for instance by measuring the total friction torques of the robot joint gear and assuming that the friction torques act on the motor axle only thus Fq({dot over (q)})=0. τext is a vector indicating the external torques acting on the output axles of the robot joint gears. The external torques can for instance be provided by external forces and/or torques acting on parts of the robot arm. For instance, if the tool flange of the robot is subject to external forces and/or torques described by Fext, the resulting torques at the output axles of the robot joints becomes:
τext=JT(q)Fext eq. 8
where JT(q) is the transposed manipulator Jacobian of the robot arm and where Fext is a vector describing the direction and magnitude of the external forces and torques in relation to the tool flange of the robot arm.
Secondly, the dynamic model as seen from the input side of the robot joint gears becomes [18]:
τmotor=B{umlaut over (Θ)}+FΘ({dot over (Θ)})+τjoint eq. 9
where τjoint is a vector comprising transmission torque τjoint,i . . . τjoint,n of each of the robot joint gears; {umlaut over (Θ)} is a vector comprising the second time derivative of the angular position of the motor axle of the joint motor and thus relates to the angular acceleration of the motor axle. B is the positive-definite diagonal matrix indicating the mass moments of inertia of the joint motor's rotors. FΘ({dot over (Θ)}) is a vector comprising the friction torques acting on the motor axles and τmotor is a vector indicating the torque generated by the joint motors.
FΘ({dot over (Θ)}) is a vector comprising the friction torques acting on the input axles of the robot joint gears. The friction torques acting on the input axle depends on angular velocity of the input axle ({dot over (Θ)}); however, it is to be understood the friction torques acting on the input axle also can depend on other parameters such as temperatures, type of lubricants, loads to robot arm, position/orientation of the robot arm etc.
FΘ({dot over (Θ)}) can for instance be provided as linear or nonlinear functions or lookup tables (LUTs) with interpolation, and FΘ({dot over (Θ)}) can be defined based on for instance prior knowledge of the robot, experiments, and/or be adaptively updated during robot operation. For instance, FΘ({dot over (Θ)}) can be obtained for each robot joint by running the robot joint gear at a constant angular velocity for at least one whole revolution of the output axle or any positive integer multiple of a revolution without the robot arm contacting external objects, while obtaining the motor torque of the robot joint motor. Under these conditions the inertia term of the robot arm M(q){umlaut over (q)} is zero due to the constant angular velocity, the Coriolis and centripetal term C(q,{dot over (q)}){dot over (q)} cancels because only one robot joint is moved at the time, and the gravity term G(q) cancels out due to the full rotation of the output axle, and no external torques are provided to the robot arm. Consequently eq. 7 reduces to:
τjoint,average=Fq,average({dot over (q)}) eq. 10
where τjoint,average indicates the average joint transmission torque during one revolution of the output axle and Fq,average({dot over (q)}) indicates the friction force acting on the output axle of the robot joint gear during one revolution of output axle. Further, under these conditions the inertia term of the joint motors B{umlaut over (Θ)} is zero due to the constant velocity. Consequently eq. 9 reduces to:
τmotor=FΘ({dot over (Θ)})+τjoint eq. 11
combining eq. 10 and eq. 11 makes it possible to obtain the friction torques acting on the motor axles as:
F
Θ({dot over (Θ)})=τmotor−Fq,average({dot over (q)}) eq. 12
The friction torque FΘ({dot over (Θ)}) can also be mechanically modeled for instance as a Generalized Maxwell-Slip (GMS) model [19] which is based on three frictional properties
For single robot joint gear, the GMS model can be visualized as a parallel connection of M massless block-spring models (illustrated in
F
Θ=Σi=1MFi eq. 13
The dynamics of each elementary model is represented by the equations (10) and (11). If the element is sticking:
where ki is the spring stiffness of the ith element. The element remains sticking until Fi>vis(ω) where the fractional parameter v determines the maximum force Fi for each element during sticking. If the element is slipping
The attraction parameter c determines how fast the total friction force approaches s(ω) in sliding.
The Nonlinear Static Map
s(ω)=sign(ω)(FC+(FS−FC)·e−(|ω|/v
captures the Stribeck effect and nonlinear viscous friction, where FC, FS, vS and μ are, respectively, the Coulomb friction, the stiction, the Stribeck velocity, and the Stribeck shape factor, and FV
Typically the robot controller is configured to control the movements/positions of the robot arm by determining the torques generated by the motor axles τmotor of the robot joint gears based on eq. 7 and eq. 9 and desired angular positions q, angular velocity {dot over (q)} and/or desired angular accelerations {umlaut over (q)} of the output axles of the robot joint gears. The desired angular positions q, angular velocity {dot over (q)} and/or desired accelerations {umlaut over (q)} of the output axles of the robot joint gears can for instance be provided via a robot program stored in the memory or via user inputs received through the input device. The angular positions q, angular velocity {dot over (q)} and/or desired accelerations {umlaut over (q)} of the output axles can be provided directly or as other parameter from which the parameters can be obtained for instance in form of coordinates of the robot tool flange in relation to the robot base or the like. The controller is then configured to control the joint motors by regulating the current through the joint motors.
Step 552 of applying a motor torque to the motor axles of the joint motors can for instance be carried out by the controller instructing the joint motor to apply a motor torque to the motor axles of the robot joint gears. For instance, the controller can be configured to provide motor torques that maintain the robot arm in a static posture, where the motor torques are sufficient to overcome the gravity on the robot arm or by providing motor torques moving parts of the robot. Consequently, the joint motors are driven in order to generate motor torques τmotor by the motor axles.
In one embodiment the step of applying a motor torque to the motor axle using the joint motor results in movement of at least a part of the robot arm. This can be achieved by driving the joint motors with current that results in rotation of the robot joints whereby at least a part of the robot arm moves. For instance, the robot controller may be configured to apply a motor torque to the joint motors that result in movement of a part of the robot arm upward in relation to gravity and in such situation the applied motor torque is larger than the motor torque needed to compensate for gravity. Also, the robot controller may be configured to apply a motor torque to the joint motors that result in movement of a part of the robot arm downward in relation to gravity and in such situation the applied motor torque is smaller than the motor torque needed to compensate for gravity. However, it is to be understood that any motor torque resulting in movement in any direction of at least a part of the robot arm can be applied. Consequently, both the input axle of the robot joint gears and the output axle of the robot joint gears rotates. The dynamic terms of the dynamic models according to eq. 7 and eq. 9, which depends on the angular velocity {dot over (q)} of the output axle, the angular acceleration {umlaut over (q)} of the output axle, the angular velocity {dot over (Θ)} of the motor axle, and/or the angular acceleration {umlaut over (Θ)} of the motor axle can be obtained during movement of the robot arm. For instance the angular velocity {dot over (q)} of the output axle, the angular acceleration {umlaut over (q)} of the output axle, the angular velocity {dot over (Θ)} of the motor axle, the angular acceleration {umlaut over (Θ)} of the motor axle can be obtained as described in paragraphs [0053]-[0054]. This makes it possible to obtain the gear stiffness of the robot joint gear during movement and operation of the robot arm.
The method can comprise an optional step 535 of obtaining the actual motor torque τactually provided by the motor axles. The actual motor torques can for instance be obtained by obtaining the current through the joint motors whereby the actual motor torque can be obtained as known in the art of motor control. For instance, if the joint motors are provided as three-phase Permanent Magnet Synchronous Machines (PMSM) with dynamics much faster than that of the manipulator. If the joint motors are operated under their current saturation limit, the motor axle torque can be obtained by:
τactually,motor=KτImotor eq. 17
where τactually,motor is a vector comprising the actual torque provided by the motor axles of the joint motors (seen in the output space of the robot joint gear), Kτ is the positive-definite diagonal matrix of torque constants and Imotor is a vector comprising the torque-generating (quadrature) current obtained from the phase currents of the joint motors using the Park Transform.
The actual torque of the motor axle can also be obtained by using force/torque sensors such as strain gauges indicating the actual torque of the motor axle. Obtaining the actual torque τactually generated by the motor axles provides a more accurate gear stiffness as eventual deviation between the desired motor signal and the actual torque can be eliminated.
Step 554 of obtaining the angular position of the motor axles of the joint motors can be obtained by measuring the angular position of the motor axle for instance by using an encoder such as optical/magnetic encoders as known in the art of robotics. The angular position of the motor axles Θ can be stored in a memory for later usage for instance in order to store a number of angular positions of the motor axles obtained at different times. The angular position of the motor axles can for instance be used to obtain the angular velocity and/or angular acceleration of the motor axles for instance by differentiation and double differentiation of the angular position with respect to time. The angular velocity of the motor axles and angular acceleration of the motor axles can then be used to obtain the robot joint gear stiffness using the dynamic terms of the dynamic model of the robot arm. Consequently, the gear stiffness of the robot joint gear can be obtained during movement and operation of the robot arm. It is to be understood that the angular velocity of the motor axles and/or angular acceleration other motor axles alternatively can be obtained using sensors/encoders measuring the actual parameters.
Step 556 of obtaining the angular position of the output axle of the robot joint gears can be obtained by measuring the angular position of the output axle for instance by using encoders such as optical/magnetic encoders as known in the art of robotics. The angular position of the output axles q can be stored in a memory for later usage, for instance a number of angular positions of the output axles obtained at different times can be stored in a memory. Alternatively, the angular position of the output axles of the robot joint gears can be obtained as the desired angular position of the output axles upon which the controller is generating the motor torques. This makes it possible to estimate the angle of the output axles in robot joints that does not comprise encoders for measuring the angular position of the output axles. The angular position of the output axle of the robot joint gears can for instance be used to obtain the angular velocity and/or angular acceleration of the output axle of the robot joint gears for instance by differentiation and double differentiation of the angular position with respect to time. The angular velocity of the output axles and angular acceleration of the output axles can then be used to obtain the robot joint gear stiffness using the dynamic terms of the dynamic model of the robot arm. Consequently, the gear stiffness of the robot joint gear can be obtained during movement and operation of the robot arm. It is to be understood that the angular velocity of the output axles and/or angular acceleration of motor axles alternatively can be obtained using sensors/encoders measuring the actual parameters.
Step 560 of determining the gear stiffness of the robot joint gears can be embodied by considering the model of the flexibility of the robot joint illustrated in
τE(Φjoint)=M(q){umlaut over (q)}+C(q,{dot over (q)})q+G(q)+Fq({dot over (q)})+τext−τD({dot over (Φ)}joint) eq. 18
τE(Φjoint)=τmotor−B{umlaut over (Θ)}−FΘ({dot over (Θ)})−τD({dot over (Φ)}joint) eq. 19
where Φjoint is a vector comprising the deflection of each robot joint gear, {dot over (Φ)}joint is a vector comprising the first-order time derivative of the deflection of each robot joint and τD({dot over (Φ)}joint) is a vector comprising the damping torque for each robot joint gear as a function of the time derivative of the deflection of the robot joints and τE (Φjoint) is a vector comprising the flexibility torque of each of the robot joint gears as a function of the deflection of the robot joint gears.
Step 560 of determining the gear stiffness of the robot joint gears comprise thus a step 562 of obtaining the deflection Φjoint of the robot joint gears based on the angular position Θ of the motor axles and the angular position q of the output axles, where the deflection Φjoint of the robot joint gears can be indicated as the difference between the angular position Θ of the motor axle and the angular position q of the output axle as indicated by eq. 1. It is to be understood that the deflection alternatively also can be indicated as the difference between the angular position of the output axle and the angular position of the motor axle. Step 560 of determining the gear stiffness of the robot joint gears comprise a step 564 of obtaining a flexibility torque of the robot joint gear at the deflection of the robot joint gear based on at least the angular position of the motor axle, the angular position of the output axle and a dynamic model of said robot arm. The flexibility torque τE(Φjoint) of the robot joint gear at the deflection of the robot joint gear can be obtained based on eq. 18, with the angular position of the motor axles Θ and the angular position q of the output axles as input variables, as the remaining parts of eq. 18 have been obtained during the initialization and for instance stored the memory. It is noted that τext is zero in situations where no external forces/torques are acting on the robot arm, however if for instance an external force is acting on the robot tool flange the external forces/torques can then be obtained based on eq. 8. Further the contribution from external forces can also be obtained by using force/torques sensors. Additionally, or alternatively τE(Φjoint) can be obtained based on eq. 19 with the angular position of the motor axles Θ, the angular position q of the output axles and motor torques τmotor by the motor axles as input variables. It is thus to be understood that the motor torques τmotor by the motor axles can be provided as an additional/optional input to step 564. In the embodiment comprising step 553 the motor torques τmotor can be replaced by the actual motor torques τactually obtained in step 553 when obtaining the gear stiffness of the robot joint gear using eq. 19. The motor torques τmotor and the actual motor torques τactually is indicated in brackets [ ] and ( ) in
In step 566 the gear stiffness of the robot joint gears can be obtained based on the deflection of the robot joint gear obtained in step 562 and the flexibility torque τE(Φjoint) of the robot joint gear obtained in step 564. Considering the eq. 4, eq. 5 and eq. 6 makes it possible to obtain the gear stiffness as:
where τE(Φjoint) can be obtained based on eq. 18 or eq. 19 as described in connection with step 564.
The gear stiffness ki(Φjoint) obtained by eq. 20 provides can be used to obtain the average gear stiffness between a certain deflection of the robot joint gear and the situation where deflection of the robot joint gear is zero, which is typically when no motor torque is provided to the robot joint gears. In such situation eq. 20 can be simplified to
Additionally or alternatively the gear stiffness of the robot joint gears can be indicated as the flexibility torque τE(Φjoint), where a high flexibility torque at a given deflection of the robot joint gear indicates a high gear stiffness, as a higher flexibility torque is needed to achieve the deflection of the robot joint gear; reversely a low flexibility torque at a given deflection of the robot joint gear indicates a low gear stiffness as a smaller flexibility torque is needed to achieve the deflection of the robot joint gear. The gear stiffness of the robot joint gears can thus be indicated based eq. 18 and/or eq. 19.
The method can comprise an optional step 570 of comparing the obtained gear stiffness ki(Φjoint) with prior knowledge of the robot joint gear stiffness ki,limit(Φjoint) and provide a status Sjoint of the robot joint gears based on the obtained gear stiffness and said prior knowledge of the robot gear stiffness. For instance, the average gear stiffness can be obtained at different times during the life time of the robot and be compared a threshold value of the gear stiffness, where the threshold value has been determined based on knowledge about the evolution of the gear stiffness in relation to wear of the robot joint gear. The status of the robot joint gear can thus be used to indicate the wear of the robot joint gears. Typically, the gear stiffness of a robot joint gears tend to decrease before the robot joint gears breaks, as material is worn off at the gear meshing. Consequently, the controller can be adapted to provide a warning when the gear stiffness reaches a threshold value, whereby the robot joint gear can be sent for repair or replaced before the gear breaks. This makes it possible to avoided failure of the robot joint gear and thereby avoid downtime of the robot arm.
The method can comprise an optional step 580 of controlling the robot arm based on the obtained gear stiffness of the robot joint gears by providing a control signal τcontrol indicative of a desired motor torque τmotor of the motor axle. The desired motor torque τmotor can for instance be obtained based on the gear stiffness of the robot joint gears by modifying the dynamic model of the robot to include the stiffness of the robot joint gears and then update the dynamic model of the robot based on the gear stiffness obtained in step 560.
In this embodiment the methods comprise a step 665 (illustrated as an arrow) of repeating the following steps:
Step 667 is a step of fitting the obtained data points comprising the flexibility torques of the robot joint gears and the corresponding robot joint gear deflection with a mathematical function. For instance, the obtained data points for a single robot joint can be fitted to polynomial function having a linearly parametrized polynomial basis:
τE,i(Φjoint,i)=ΨT(Φjoint,i)αi eq. 22
where ΨT(Φjoint) is the linear polynomial basis of the deflection variable Φjoint, αi is a vector comprising the constants for polynomial basis functions of the ith robot joint. Based on the symmetry assumptions in eq. 5 and eq. 6 the regressor ΨT(Φjoint) can be specified to contain only odd powers of Φjoint, thus:
ΨT(Φjoint)=Φjoint2p−1, p=1,2, . . . ,P eq. 23
In summation form the polynomial function can be written as:
where αp,i is a vector comprising the constants for the polynomial basis functions for the ith robot joint. In a robot arm comprising a plurality of robot joints a polynomial function and corresponding vector of constants can be provided for each robot joint.
As described in paragraph [0027] and stated by eq. 4 the stiffness of the robot joint gear can then be determined as the change of flexibility torque of the robot joint gears with respect to the robot joint gear deflection. The method can comprise a step 668 of obtaining the gear stiffness of the robot joint robot by obtaining the slope of the flexibility torque with respect to the deflection of said robot joint gear. For instance, the slope can be determined between at least two of the data points indicated by eq. 20 and by eq. 21 for the situation where no flexibility torque exists at a deflection of zero. The slope can also be indicated by differentiation of the mathematical function fitted to the data points. In case of a polynomial function this corresponds to differentiation of eq. 24
The method of obtaining the gear stiffness of the robot joint gears of a robot joint have been demonstrated on a robot arm by obtaining the flexibility torque of a single robot joint gear as a function of robot joint gear deflection using both a method according to the present invention and an alternative method. The experiments have been performed using a UR5e robot arm provided by Universal Robots and the results are discussed in the following paragraphs.
During the experiments the robot arm 701 is arranged in a pose where the base axis 711A is parallel with the direction of gravity 723 and the shoulder axis 711b, the elbow axis 711c, the first wrist axis 711d, the second wrist axis 711e and the tool axis 711b are perpendicular to the direction of gravity 723. The base joint 703a is the only joint which is rotated during the experiment and the shoulder joint 703b, the elbow joint 703c, the first wrist joint 703d, the second wrist joint 703e and the tool joint 703b are thus fixed during the experiments.
The gear stiffness of the robot joint gear of the base joint 703a is obtained in the experiments. Firstly, the friction and the stiffness characteristics for the base joint 703a of the UR5e Robot arm is obtained using an alternative method.
The steady-state friction torque is identified by imposing different signals of constant velocity on the base joint while measuring the robot joint current. For constant velocity eq. 7, eq. 9 and eq. 17 yield Fθ+Fq=KτI. The steady-state friction torque is illustrated in
The flexibility torque of the robot joint gear of the base joint is obtained in a static setting by locking the input axle (motor axle) of the robot joint gear and imposing a set of known torques on the output axle using a Sauter FH-S 500 digital force gage while measuring the deformation of the robot joint gear using the output and input encoders of the robot.
The flexibility torque of the robot joint gear of the base joint is obtained dynamically based on the method according to the present invention by rotating the base joint according to the bang-coast-bang joint space trajectory generated with randomly generated waiting times and angular positions, velocities, and accelerations in intervals ranging from 30% to 100% of the maximum allowed values.
During the experiment the angular position of the output axle of the base joint, the angular position of motor axle of base joint, and the current trough the base joint motor are measured while the robot arm moves according the joint space trajectory illustrated in
Estimating time-varying parameters with the Recursive Least Squares (RLS) approach is obtained by incorporating a forgetting factor discounting past data. In the present example, the incoming information (the transmission deformation) is nonuniformly distributed over the parameter space. Tracking can only happen in some direction if there is an excitation in that same direction, hence the estimation algorithm tracks time-varying parameters only within the excited subspace. In particular, inspired by the results in [21], by a suitable notion of excitation subspace, the parameter vector in eq. 22 is estimated by the procedure:
where δ>0 enforces an increment of the covariance matrix P improving the algorithm alertness. More alertness is achieved by decreasing μ at the price of an increased sensitivity to disturbances. By setting μ=1 and δ=0, eq. 31 reduces to the standard RLS algorithm.
The fit of the polynomial bases (static POLY fit 989) of the flexible torque obtained statically as descried in paragraph [0063] have also been included in the
The present invention makes it thus possible to obtain the flexibility torque and the gear stiffness of the robot joint gears of a robot arm online during ordinary usage of the robot arm, where the robot arm typically performs a number of movements for instance in connection with a production setup.
| Number | Date | Country | Kind |
|---|---|---|---|
| 18194683.1 | Sep 2018 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2019/074075 | 9/10/2019 | WO | 00 |
