AXY:Basics of joint parameters identification

Model based control can bring great benefit to robot control, but it requires exact model parameters. This document describes how these parameters can be identified (automatically).

Record identification motion profile

The motion profile should provoke all motion states described by the model. Ideally it should move over the whole position range of the joint, have plenty accelerations/decelerations and have different constant velocity phases.

The recorded file must contain the position feedback along with the torque feedback.

Derivative without phase shift

For parameter identification joint velocity and acceleration is needed. Standard numeric derivative of position feedback leads to phase shift, which impairs the identification. To avoid the phase shift interpolated derivative should be used instead. Hereby the numeric derivative is spline-interpolated and values at original sampling times are taken.

Interpolated derivative is implemented in aico.motionAnalyser by the function signal_deriveInterp.

Mathematical derivation

In a simple case an actuating joint can be described with the following motion equation:

$\tau =\Theta \cdot {\ddot {x}}+\rho _{v}\cdot {\dot {x}}+\rho _{c}\cdot sign({\dot {x}})\qquad \mathrm {(1)}$

Here the symbols are:

$\tau$	Torque (feedback)
${\dot {x}},{\ddot {x}}$	Velocity and acceleration (feedback)
$\Theta$	Moment of inertia
$\rho _{v}$	Viscous friction coefficient
$\rho _{c}$	Coulomb friction coefficient

It is assumed that the joints position, velocity and acceleration feedbacks together with the torque feedback are known (recorded). The unknown model parameters are $\Theta$, $\rho _{v}$ and $\rho _{c}$.

The recorded value tuples can be represented as a over-determined linear equation system:

${\begin{array}{lcl}\tau _{1}&=&\Theta \cdot {\ddot {x}}_{1}+\rho _{v}\cdot {\dot {x}}_{1}+\rho _{c}\cdot sign({\dot {x}}_{1})\\\tau _{2}&=&\Theta \cdot {\ddot {x}}_{2}+\rho _{v}\cdot {\dot {x}}_{2}+\rho _{c}\cdot sign({\dot {x}}_{2})\\\tau _{3}&=&\Theta \cdot {\ddot {x}}_{3}+\rho _{v}\cdot {\dot {x}}_{3}+\rho _{c}\cdot sign({\dot {x}}_{3})\\&&\vdots \\\tau _{n}&=&\Theta \cdot {\ddot {x}}_{n}+\rho _{v}\cdot {\dot {x}}_{n}+\rho _{c}\cdot sign({\dot {x}}_{n})\\\end{array}}$

These n equations can also be written as a single matrix equation:

${\begin{bmatrix}\tau _{1}\\\tau _{2}\\\tau _{3}\\\vdots \\\tau _{n}\end{bmatrix}}={\begin{bmatrix}{\ddot {x}}_{1}&{\dot {x}}_{1}&sign({\dot {x}}_{1})\\{\ddot {x}}_{2}&{\dot {x}}_{2}&sign({\dot {x}}_{2})\\{\ddot {x}}_{3}&{\dot {x}}_{3}&sign({\dot {x}}_{3})\\\vdots &\vdots &\vdots \\{\ddot {x}}_{n}&{\dot {x}}_{n}&sign({\dot {x}}_{n})\\\end{bmatrix}}\cdot {\begin{bmatrix}\Theta \\\rho _{v}\\\rho _{c}\\\end{bmatrix}}\qquad \Rightarrow {\boldsymbol {\tau }}=\mathbf {X} \cdot \mathbf {p} \qquad \mathrm {(2)}$

One solution of this system of equations can be obtained by using the pseudoinverse of $\mathbf {X}$.

$\mathbf {p} =\left(\mathbf {X} ^{T}\cdot \mathbf {X} \right)^{-1}\cdot \mathbf {X} ^{T}\cdot {\boldsymbol {\tau }}\qquad \mathrm {(3)}$

HowTo identify joint parameters

1. Use AF function IDM_IdentParamsJoint to move an identification profile and create a record.
2. Process the obtained record with the aico.motionAnalyser function joint_identDynamics.

Verification

See verification page.