Next Article in Journal
Improving Discrimination in Color Vision Deficiency by Image Re-Coloring
Next Article in Special Issue
Development of a Virtual Force Sensor for a Low-Cost Collaborative Robot and Applications to Safety Control
Previous Article in Journal
Electromyographic Evaluation of the Impacts of Different Insoles in the Activity Patterns of the Lower Limb Muscles during Sport Motorcycling: A Cross-Over Trial
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Parameter Identification for a Manipulator with Joint Torque Sensors Based on an Improved Experimental Design

1
School of Mechanical Engineering, Hebei University of Technology, Tianjin 300130, China
2
State Key Laboratory of Robotics and System, Harbin Institute of Technology (HIT), Harbin 150006, China
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(10), 2248; https://doi.org/10.3390/s19102248
Submission received: 10 April 2019 / Revised: 9 May 2019 / Accepted: 13 May 2019 / Published: 15 May 2019
(This article belongs to the Special Issue Sensors and Robot Control)

Abstract

:
As the foundation of model control, robot dynamics is crucial. However, a robot is a complex multi-input–multi-output system. System noise seriously affects parameter identification results, thereby inevitably requiring us to conduct signal processing to extract useful signals from chaotic noise. In this research, the dynamic parameters were identified on the basis of the proposed multi-criteria embedded optimization design method, to obtain the optimal excitation signal and then use maximum likelihood estimation for parameter identification. Considering the movement coupling characteristics of the multi-axis, experiments were based on a two degrees-of-freedom manipulator with joint torque sensors. Simulation and experimental results showed that the proposed method can reasonably resolve the problem of mutual opposition within a single criterion and improve the identification robustness in comparison with other optimization criteria. The mean relative standard deviation was 0.04 and 0.3 lower in the identified parameters than in F1 and F3, respectively, thus signifying that noise is effectively alleviated. In addition, validation experimental curves were close to the estimation model, and the average of root mean square (RMS) is 0.038, thereby confirming the accuracy of the proposed method.

1. Introduction

A manipulator is a complicated multi-input–multi-output (MIMO) system, which has strongly coupled nonlinear characteristics. Given structured and unstructured uncertainties, such as elastic deformation, assembly clearance, inertia, Coriolis, gravity, and friction torque, the repetitive positioning accuracy of a robot and the perceived accuracy of joint torque information are seriously affected. Most robot manufacturers do not provide related information or partial parameters to obtain an accurate dynamic model [1,2]. Therefore, experimental identification is the optimal choice for obtaining these types of information. However, sensor measurement and process noise will challenge experimental identification. To mitigate the effects of nonlinearity and improve the advanced model-based control performance, a complete noise reduction system must be considered to perform the accurate estimation of dynamic parameters, including statistics of noise characteristics, excitation signal optimization, and noise processing.
Some researchers have provided significant contributions to the model of robot dynamics identification. Gautier constructed a model based on energy identification in [3] and a power model in [4], depending only on the functions of joint position and velocity to avoid introducing acceleration significantly, thereby leading to including non-negligible noise in the estimation process. A direct inverse dynamic identification model, which is based on an output error, has been proposed in [5]; this model solely requires force/torque measurement to substitute a usual position output. The most common model of identification is inverse dynamic identification model (IDIM), which provides more information than other models. In addition, the model can be linearly formulated with respect to a set of minimal dynamic parameters, thereby allowing us to construct a well-conditioned over-determined regression matrix [6,7,8].
Subsequently, parameters are estimated through numerical optimization based on the abovementioned models. The least squares (LS) method [9] and maximum likelihood estimation (MLE) [10] technique are the typical robot parameter identification methods under the class of estimation algorithms. Moreover, these techniques can be expanded in several ways, such as the weighted least squares (WLS) estimation method [11] and nonlinear least squares (NLS) optimization method [5]. Another common alternative estimation approach is based on the Kalman filter and its expansions, such as the extended Kalman filter (EKF) and unscented Kalman filter [12]. Gautier compared the WLS method with EKF estimation for parameter identification of a two degrees-of-freedom (DOF) robot in [13]; the estimation results obtained through the two techniques were similar, but the EKF was sensitive to initial conditions, and its convergence velocity was slower than that of the WLS method. With the further improvement of instrument variable theory, some researchers have aimed to bridge the gap between theory and control engineering practices and complete the MLE method with profound statistical characteristics in [14,15,16]. These attempts have undeniably achieved favorable identification effects, but the choice of auxiliary variables is significantly affected by system control law. Moreover, several other identification methods for a manipulator are based on the genetic algorithm (GA) in [17] and neural networks to improve the accuracy of estimation in [18,19]. Detailed summaries of these methods can be found in [6,20] for further understanding.
Furthermore, an experimental design is a crucial link in a closed-loop identification process, because the noise on joint position and output torque measurements, friction, and other unmodeled effects seriously affect the accuracy of estimated parameters. To reduce the bias and obtain favorable results, the exact statistical properties of the uncertainty terms must be determined. An appropriate experiment is designed by reducing the sensitivity of the system to noise and minimizing the values of the variances of the estimated parameters [14,16]. This process consists of two parts, namely, excitation trajectory design and optimization criterion. In [10,21], the closed-loop identification process is combined with an optimal excitation trajectory using parameterized finite Fourier series. The optimization process implements an iterative procedure until the Cramér–Rao lower bound is reached. Atkeson [22] proposed an excitation trajectory of a fifth-order polynomial in a joint space. Moreover, the excitation signal based on modified Fourier series [2], the finite sum of harmonic sine functions [17], and the function formed by combining several ways [23] are extended well. Two main optimization criteria are adopted; one criterion is based on minimizing the condition of an observation matrix [24]; the other criterion is to minimize the variations in the Fisher information matrix [10,21]. However, noise exists in the input and output, but most optimization technologies are justified by observing a one-sided error and cannot perform full optimization. In addition, considering the time cost in the optimization process, Jingfu [25] proposed a new trajectory design method based on Hadamard’s inequality. This method is certainly a good innovation, but time is relatively loose and seems insignificant in offline optimization identification.
The present research aims to realize the dynamic parameter identification of a manipulator as accurately as possible. The primary contribution is that an improved optimization algorithm is proposed on the basis of the multi-criteria embedded nonlinear optimization method to design the exciting trajectory for minimizing the impact of noise. Considering the influence of a certain criterion on solving the parameter identification, our approach aims to reduce the condition of the observation matrix and determinant of the Fisher information matrix simultaneously in a nonlinear optimization process. In the present research, the model that will be considered is based on IDIM and is linear in the parameters to be identified, thereby implying that the input and output of the robot correspond to the output and input functions of the identification model. In addition, extracting the base set of the parameters is advantageous. Furthermore, the appropriate weight values are obtained on the basis of the statistical characteristics of noise. All joints of the robot are designed simultaneously in the identification procedure to solve the MIMO coupled problem. The simulation and experimental results show that this method can identify parameters accurately and is superior to the other nonlinear optimization methods in terms of robustness, time consumption, and other factors. In addition, the proposed optimization method is not limited to a series robot but also applicable to the optimized experimental design of parallel and exoskeleton robots. This condition is due to the fact that, to obtain a high-precision dynamic model, the input must fully stimulate the motion form of the system to ensure that the experimental data can fully reflect the physical characteristics of the system, and a good excitation design can reduce the impact of noise on the results. Examples of this include system identification of aircraft, parameter identification of sorting robots, friction identification, force-free control and collision detection.
The remainder of this paper is organized as follows. Section 2 introduces the development of a dynamic linear model, and the MLE method is used to estimate the robot parameters. Section 3 presents the different optimization criteria and describes the process of obtaining optimal excitation trajectories based on a single criterion and with the proposed multi-criteria embedded nonlinear optimization method. Section 4 discusses the identification experiment process and compares the traditional method with the proposed algorithm. Section 5 analyzes the relevant experimental results. Section 6 provides some conclusions drawn from this research.

2. Modeling and Identification of Robot Parameters

In this section, we obtained the identification procedure using the IDIM and the WLS estimation technologies to identify the base parameters of a multi-jointed robot. To reduce the system noise, we used mean and low pass filtering based on a periodic excitation trajectory to improve the signal-to-noise ratio. Moreover, we could obtain the estimated results using the linearized formulation model. The specific identification procedure is illustrated in Figure 1.

2.1. Dynamic Identification Model

Generally, we can determine the dynamic model of an n-joint rigid robot based on the Euler–Lagrange or the Newton–Euler formulation in the joint space as
M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + G ( q ) + τ f = τ
where M ( q ) n × n is the inertia matrix, C ( q , q ˙ ) n × n contains Coriolis and centripetal matrix, G ( q ) n is the vector of the gravitational matrix, τ f n is the friction force term, τ n is the input torque vector to the system, and q n is the vector of relative generalized coordinates used to model the system, q ˙ , q ¨ n represents the generalized velocities and accelerations based on differentiation of q. The friction forces are described as [26]
τ f = F c sgn ( q ˙ ) + F v q ˙
where F c , F v n × n are the diagonal matrices that describe the coefficient of Coulomb and viscous friction, correspondingly.
The mathematical model (1) in a linear parameterized form that depends on a set of parameters must be rewritten for dynamic parameter identification. Modified Denavit and Hartenberg [27] convention can help us obtain the linear model with Ns standard parameters.
τ = ϕ s ( q , q ˙ , q ¨ ) θ s
where ϕ s ( q , q ˙ , q ¨ ) n × N s is the regression matrix of a nonlinear function of joint position, velocity, and acceleration vectors, and θ s N s × 1 is the vector of standard parameters to be estimated. A total of 12 standard parameters are used by each link and joint for rigid robots, which contain the mass mj of each link j, the first mass ( m x j m y j m z j ) moments of link j, the six components of the inertia tensor of link j at the origin of frame j, ( I x x j I x y j I x z j I y y j I y z j I z z j ) and the coefficients of viscous and Coulomb friction torque of joint j ( F c j F v j ) .
The factor that truly affects robot dynamics is only part of the inertial parameters and can be estimated through identification procedures. The set of a minimum number of inertial parameters, that is, the base parameter set, can be selected on the basis of a numerical method with respect to the QR decomposition [28] or by regrouping the standard parameters through linear relations [29,30]. The dynamic equation (3) can then be regrouped into another form with Nb identifiable base parameters.
τ = ϕ b ( q , q ˙ , q ¨ ) p b
where ϕ b ( q , q , q ) n × N b is a subset of the regression matrix ϕ s , and p b N b is the base parameters by regrouping.

2.2. Data Acquisition and Signal Processing

In fact, the known system must be given continuous excitation signals to demonstrate the mechanical and physical properties fully. Although the observation matrix ϕb has full rank, and pb can be obtained through the pseudo-inverse matrix from Equation (4), the solutions of identification are local and unreliable given the complex modeling errors, measurement noise, friction, and other factors. The ideal dynamic model is inexistent. Thus, assuming that the fundamental frequency ωf, sampling frequency ωs, and the kth sampling time as tk of the excitation trajectory are known, Equation (4) can be extended in an over-determined set with M = ωs/ωf measurement points over one period T as follows:
Γ = F ( q , q ˙ , q ¨ ) p b + ε
where Γ n M × 1 is the measured force/torque vector; F ( q , q ˙ , q ¨ ) n M × N b is the observation matrix, also called the global information matrix; ε n M is the zero mean random vector of errors derived from the term of nonlinear noise.
Considering that the collected data (e.g., position and torque data) contain noise, it is a fundamental step to reducing the influence of noise on τ and F for improving the accuracy of parameter identification results. Here, velocity and acceleration are the numerical differentiation with respect to the position q. In the process of data collection, given the periodicity of excitation trajectory, we combined the mean filtering algorithm and a second-order low-pass digital Butterworth filter to process signal noise, especially acceleration and torque.

2.3. MLE for Parameter Estimation

System identification aims to determine a system model based on input and output data from a specified class of models, thereby making it equivalent to the estimated system. The principle is depicted in Figure 2. The same signal X is used to excite the system prototype Mp and the system model Me. Moreover, the output signals are τm and τp, respectively, and the error is e. An error criterion function J(a) = eTΣe, which can be taken as the function of the error, is specified to modify the parameter vector a. Recursion is repeated until the error criterion J(a) is minimal.
Robot identification handles the problem of estimating the model parameters from the measured ones based on a statistical framework during a robot excitation experiment. Assuming that the measured joint position qm(k) and torque τm(k) at kth sampling time conform to a random zero-mean Gaussian noise εq(k), ετ(k), i.e.,
q m ( k ) = q ( k ) + ε q ( k ) τ m ( k ) = τ ( k ) + ε τ ( k )
where εq(k) and ετ(k), i.e., εq(k) + ετ(k) = ε, are integrated to satisfy Equation (5). The noise-free version of these variables satisfies Equation (4).
From the perspective of statistics, the MLE method is conducted for the dynamic parameters of the robot [10]. First, a likelihood function of joint angles qm(k) and joint moment τm(k) measurement value is constructed. Second, the MLE of parameter vector p is obtained by maximizing the function. Given that the noise on different measurements are independent of each other and obey the Gaussian distribution, the following quadratic cost function can be minimized:
L ( q m , τ m | p ) = 1 2 k = 1 M j = 1 n ( ε q j 2 ( k ) σ q j 2 + ε τ j 2 ( k ) σ τ j 2 )
where L ( q m , τ m | p ) is the likelihood function, εqj(k) and ετj(k) are the noise on the joint position and torque of the jth joint, respectively. σ q i 2 and σ τ i 2 are their corresponding variances. Detailed steps for determining these variables can be obtained in [31].
Considering Equation (4), the minimization problem of Criterion (7) will be transformed into an NLS optimization problem, thereby implying the necessity to obtain the εqj(k) and ετj(k) of every estimated parameter p with respect to the given measured data qm(k) and τm(k) in practical implementation. This process is nearly infeasible using the present formulation. However, we can find that, in the actual experimental process, the noise level of the measurement value is much lower in the joint angle than in the actuator force of the joint. Therefore, no noise and error may be realized in measuring the joint angle. At this time, the maximum likelihood estimator of the parameters can be obtained by maximizing the likelihood function. Considering that the maximum likelihood function is equivalent to the logarithm of the likelihood function, the maximum likelihood parameter estimation method can be significantly simplified by taking the logarithm of the original likelihood function. Under this assumption, the MLE method reduces to the linear WLS estimation algorithm, whose weighted function is the inverse of the noise to the standard deviation of the torque value of the measured actuator [31].
p m l = ( F t Σ 1 F ) 1 F t Σ 1 Γ
where
F = [ ϕ ( q m ( 1 ) , q ˙ m ( 1 ) , q ¨ m ( 1 ) ) n × N b ϕ ( q m ( 2 ) , q ˙ m ( 2 ) , q ¨ m ( 2 ) ) n × N b ϕ ( q m ( M ) , q ˙ m ( M ) , q ¨ m ( M ) ) n × N b ]
Γ = [ τ m ( 1 ) n τ m ( 2 ) n τ m ( M ) n ]
and Σ is the diagonal covariance matrix of the measured actuator torques such that Σ = d i a g ( σ 1 2 I M σ j 2 I M σ n 2 I M ) , where IM is the (M × M) identity matrix [30]. The covariance matrix of MLE pml is equal to
( F T Σ 1 F ) 1
Furthermore, if all the noise of the measured actuator torque has the same standard deviation, then the MLE will further degrade to the standard linear LS method (e.g., p l s = ( F t F ) 1 F t Γ ). This result certainly loses the significance of the statistical framework.

3. Obtaining the Optimal Robot Excitation Trajectories

The optimal experiment design refers to the problem with respect to a nonlinear optimization with motion constraints, that is, generally, constraints on joint angles, velocities, and accelerations, and on the robot end effector positions in the Cartesian space to avoid collision with the environment and with itself. In addition to the generality, some researchers have added other conditions to the constraints, such as singular values of the observation matrix [32]. Given the actuator displacement as the control input, a criterion must be designed to represent the sensitivity of the input and output variables with respect to the solution of the estimate parameters. Considering the interaction between different optimization criteria, our main contribution is to propose a multi-criteria embedded nonlinear optimization method to obtain the optimal excitation.

3.1. Optimal Criteria for the Experiment Design

Common optimization problems are based on the deterministic framework and obtain the solution of dynamic parameters by minimizing the error term (i.e., ε in Equation (5)). To equilibrate the disturbance influence of the input and output noise on the parameter estimates, the condition number of the regression matrix is typically selected as the optimization objective function [2]. Gautier and Khalil [24] optimized a linear combination with a factor for equilibrating the observation matrix using this criterion and a fifth-order polynomial trajectory. Armstrong [33] described the optimization design of the experiment trajectory based on maximizing the minimum singular value of the square matrix FTF; some ways to expand on this condition have been used in [17,32]. Jingfu [25] proposed an approach by using the trace of the square matrix FTF based on Hadamard’s inequality. The collection of these methods can be found in [24,34].
Another common optimization problem is based on a statistical framework, and its most prominent feature is that the statistical characteristics of noise are considered in the optimization process. The representative one is the optimal objective function of the deformation of the Fisher information matrix (e.g., (FTΣ−1F)−1) presented by Swevers [10] for the MLE parameter identification of a serial robot. Afterward, it has been applied for a parallel robot [35]. Subsequently, many improvements and applications based on this criterion can be found in [15,21,23,31,36].
Notably, a new criterion based on a multi-objective optimization method is proposed in [37]. The main contribution of this approach is to find trajectories with favorable properties in accordance with more than a single criterion. This method solves the uncertainty of some solutions when the lower boundary is non-convex, but it must obtain the appropriate optimization weight value of the objective function in advance, thus undoubtedly increasing the computation and easily falling into the local solution of multiple optimizations.
From the linearized parameter model (5), we can determine that the condition number of the regression matrix represents the influence of the disturbance on the identification solution, and the physical meaning represented by the determinant value of the covariance matrix is the volume of the area of the maximum probability density function of the parameter to be estimated. Given the abovementioned characteristics, we propose a multi-criteria embedded optimization method, in which the optimization criteria (−log(det(FTΣ−1F)−1)) are embedded in the constraint conditions, and the interaction between the two criteria is compromised. In addition, the embedded optimization method makes the dual objective function optimized simultaneously in the optimization iteration process, thus avoiding the sensitivity of the multi-objective optimization approach to the goal weight value. The physical meaning is to obtain a regression matrix with good “behavior” when the solution of the estimated parameters is relatively determined, so as to reduce the sensitivity of identification results to interference.
The abovementioned optimization criteria are summarized in Table 1. Next, we will obtain the excitation trajectory based on our proposed optimization method.

3.2. Optimal Excitation Signal for the Experiment Design

Swevers proposed the excitation trajectory based on a finite Fourier series in [10]. The trajectory for each joint is a finite sum of N harmonic sine and cosine functions. The main advantages of these functions are that: First, they satisfy the numerical average in the time domain due to periodicity and conveniently improve the signal-to-noise ratio of the measurement. Second, the velocity and acceleration can be calculated using the position vector in analytical forms without using the numerical differentiation method. In addition, the effect on robot flexibility can be avoided by selecting an appropriate frequency bandwidth. The specific form for the jth of an n-link robot is defined as follows
q j ( t ) = l = 1 N a j , l ω f l sin ( ω f l t ) l = 1 N b j , l ω f l cos ( ω f l t ) + q j 0 q ˙ j ( t ) = l = 1 N a j , l cos ( ω f l t ) + l = 1 N b j , l sin ( ω f l t ) q ¨ j ( t ) = ω f l = 1 N b j , l l cos ( ω f l t ) l = 1 N a j , l l cos ( ω f l t )
where j = 1,2. t ϵ [0, T] with T is the period cycle time, ωf is the fundamental frequency, and ωf = 2π/T. N is the number of harmonics, and qj0 is the position offset of the joint of the excitation reference trajectory. In Reference [31], the fundamental frequency ωf is considered as a variable of the optimization problem to be added to the optimization design of the excitation trajectory, but the optimization results are not reflected in this process, so the selected variables are inappropriate. The base frequency must be in a good bandwidth range (e.g., the so-called bandwidth, which reflects the noise filtering characteristics of the system and also used to measure the transient response performance of the system). The trade-off of selecting the fundamental frequency is discussed in [10], it cannot be used as a variable in this optimization process. Here, 2N + 1 parameters are applied as the variables for the trajectory optimization problem of per joint. All variables include the Fourier coefficients of two joints (e.g., qi0, al = [aj,l,…,aj,N], bl = [bj,l,…,bj,N]).
Various optimization criteria are listed in Table 1. Considering the interaction between single different optimization criteria and combining with the actual physical meaning representation, the condition number of the regression matrix is adopted as the optimization objective function. In order to make the identification results insensitive to noise, the optimization problem is actually expressed as searching suitable values in constraint domain to minimize the objective function:
objFunc = F4 = minimize Cond (F)
Generally, the constraint domain mainly depends on the kinematics and geometric information of the robot system, including the position, velocity, acceleration information of the actuator and start and end points of joints. However, the optimization results are often affected by singular solutions in the optimization process, and it is easy to fall into the local optimum. In order to avoid these situations, we add F2 = (−log(det(FTΣ−1F)−1)) to the constraint as an additional constraint. Thereby, the excitation trajectory must be subjected to the following constraints:
{ ( log ( det ( F T Σ 1 F ) 1 ) ) < ω 2 ; | q j ( t ) | q max i , t | q ˙ j ( t ) | q ˙ max i , t | q ¨ j ( t ) | q ¨ max i , t q j ( t 0 ) = q j ( t f ) = 0 i , t q ˙ j ( t 0 ) = q ˙ j ( t f ) = 0 i , t q ¨ j ( t 0 ) = q ¨ j ( t f ) = 0 i , t
where ω2 represents the constraint target value of F2, which is derived from the objective value obtained when only F2 is used as the optimization objective function; q max (rad), q ˙ max (rad s−1), and q ¨ max (rad s−2) denote the bounds of joint position, velocity, and acceleration, respectively; and t0 and tf correspond to the initial and the end times; these constraints indicate that the position, velocity, and acceleration of the robot at the initial and end points of the reference trajectory are equal to 0, otherwise a strong tremor of the system will result and affect the accuracy of the parameters. Objective values are updated iteratively until the condition number of the regression matrix converges to the minimum and the optimization results satisfy the whole constraint space, the optimization terminates. At this point, the multi-objective embedded optimization method has been established.

4. Identification Implementation Process

A planar two degrees-of-freedom (2-DOF) manipulator was adopted to evaluate the effectiveness of the method, which is designed here, as demonstrated in Figure 3. Compared with a single joint, the 2-DOF manipulator considers coupling factors, and its method is applicable to the parameter identification of additional DOFs. The considered manipulator has two revolute joints equipped with brushless motors produced by Maxon, which is equipped with an encoder for position measurements. The output end of the motor is equipped with a harmonic reducer at a reduction ratio of 80, and the output end of the reducer is equipped with a spoke-type torque sensor. To protect the torque sensor from axial force crosstalk, deep-groove ball bearings and thin-wall cross roller bearings are simultaneously equipped at the connection between the link and joint for unloading the undesired load. The control unit of the robot is constituted by Platinum Maestro Network Motion Controller based on EtherCAT, which is produced by Elmo, and the control algorithms are implemented in C++. The joint coordinates are defined in accordance with Denavit–Hartenberg notation and collected in vector q = [q1 q2]T, where qi represents the angular position of joint j. The command torque vector is defined as τ = [τ1 τ2]T, where τi represents the torque applied to joint j.

4.1. Dynamic Model of the Planar Manipulator

The base parameter model of this robot contains five parameter combinations of the links with respect to the z0-axis and the first-order moments of the second link based on applying the dynamic formulation derived by Craig [38]. The friction that is modeled only in prismatic joints (independent generalized coordinates) given the friction in spherical and rotational joints (dependent generalized coordinates) can be neglected [39]. Therefore, we can obtain the vector pb of the identifiable parameters using regrouping standard parameters presented as follows:
p b = [ ( m 1 + m 2 ) d 1 2 + m 2 d 2 2 m 2 d 1 d 2 m 2 d 2 2 ( m 1 + m 2 ) g d 1 m 2 g d 2 F v 1 F c 1 F v 2 F c 2 ] T
Following the Newton–Euler approach, the planar manipulator dynamic model is expressed in Equation (4). From this equation, we can determine the regression matrix, including
F = [ q ¨ 1 c 2 ( 2 q ¨ 1 + q ¨ 2 ) s 2 ( 2 q ˙ 1 q ˙ 2 q ˙ 2 2 ) q 2 2 s 1 s 12 q ˙ 1 sgn ( q ˙ 1 ) 0 0 0 c 2 q ¨ 1 + s 2 q ˙ 1 2 q ¨ 1 + q ¨ 2 0 s 12 0 0 q ˙ 2 sgn ( q ˙ 2 ) ] ,
where si = sin(qi) and ci = cos(qi).

4.2. Simulation of the Excitation Trajectory Optimization

Several kinds of excitation trajectories are obtained by the optimization design of condition number (Cond(F)) as a criterion or (−log(det(FTΣ−1F)−1)) as the criterion or based on multi-objective optimization as the criterion. For each optimization, the trajectory is parameterized through five-term Fourier series, yielding 11 parameters for each joint (e.g., N = 5). The fundamental frequency of the trajectories is 0.05 Hz, and the sampling rate of the simulation is 100 Hz. The length of the data sequence M is 2000 data in one period. The optimization process was performed in MATLAB R2017a environment using the nonlinear constrained optimization tools “Fmincon” and “Fgoalattain”. According to the motor parameters and motion constraints, the optimization constraints are set as follows:
  • Embedded criterion F2 = (−log(det(FTΣ−1F)−1)) < −75;
  • Joint angle limits (rad): −3.14 < q1 < 3.14, −3.14 < q2 < 3.14;
  • Joint velocity limits (rad/s): −5.2 < q ˙ 1 < 5.2, −5.2 < q ˙ 2 < 5.2;
  • Joint acceleration limits (rad/s2): −4.5 < q ¨ 1 < 4.5, −4.5 < q ¨ 2 < 4.5;
  • The position, velocity, and acceleration of the two joints at the initial and end times are 0. e.g., q j ( 0 ) = q ˙ j ( 0 ) = q ¨ j ( 0 ) = q j ( 20 ) = q ˙ j ( 20 ) = q ¨ j ( 20 ) .
In addition, the covariance matrix Σ in the optimization process was derived from a random trajectory, and the parameter was estimated through the traditional linear LS method (refer to [30] (p. 294)), as plotted in Figure 4. This process aims to collect the noise characteristics of the system, that is, σ1 = 5.0832 N2m2, and σ2 = 0.0499 N2m2. The target values of the multi-objective optimization criterion were obtained from the optimization results of the former criterion. In Table 2, we selected ω1 = 8 and ω2 = −75.
To avoid special situations, the results of each criterion for four selected trajectories are summarized in Table 2. When solving the optimization problem with the F1 criterion, the values of F2 were higher than those obtained when only F2 was used. Furthermore, the opposite occurs when the optimization problems were solved by F2. Moreover, the optimization process is easily trapped in the local minimum. This phenomenon had been resolved, that is, it will not significantly increase the value of other criteria when the relevant appropriate solutions are obtained. Either by using a weighted objective function to convert two objective functions in a single function as provided in [24,33] or by using a multi-objective optimization procedure based on the goal programming in [37] (the results of multi-objective optimization are shown in the last two columns of Table 2). However, the optimization results are sensitive to the goal values given the value setting of 2-DOF, as listed in Table 3. Therefore, a multi-criteria embedded (e.g., F4) nonlinear optimization method is considered.

4.3. Multi-criteria Embedded Nonlinear Optimization

According to Equation (14), the F2 criterion is added to the constraint condition. The middle two columns of Table 3 show the results of the multi-criteria embedded optimization. In this table, this method is mainly affected by the weight values in the constraint conditions and reduces the DOF in comparison with the multi-objective optimization criteria. In addition, the optimization method is equivalent to introducing penalty function into the objective function, thereby slightly improving the robustness, and this approach has a clear physical meaning.
Similar to the optimization process of F1 and F2, the nonlinear constraint optimization function “Fmincon” was also used in the MATLAB R2017a environment. This optimization tool adopts the sequential quadratic programming (SQP) algorithm to solve the quadratic programming sub-problem in each step of the iteration process and further update the Lagrangian Hessian matrix until the value of objective function converges to the minimum and the optimization results satisfy the whole constraint space. In order to improve the optimization precision, we extend the iteration times and the error tolerance to enough to search for the solutions within the constraint space, such as “MaxFunEvals” or “TolFun,” which is set as 5000, 1e-10, respectively. Moreover, the F2 criterion in the constraint condition must satisfy the target value ω2 = −75, and the optimization must satisfy the Cramér–Rao lower bound in [10]. Figure 5 exhibits the influence of the number of points used for representing the trajectory on Cond(F). The final optimization converges to a certain value. The number of trajectory points is affected by some evaluation termination conditions in the optimization function. The evaluation conditions can be modified to obtain the optimal solution convergence. Figure 6 depicts an optimization trajectory obtained through the proposed method. The joint position, velocity, and acceleration correspond to the independent motion parameters of the actual 2-DOF manipulator.

4.4. Experimental Procedure

The parameter identification process of the 2-DOF manipulator was conducted as follows:
Experimental design. Some optimization trajectories were obtained in accordance with the optimization criteria discussed previously. To achieve an improved comparison, the optimal trajectory obtained from the same starting point was selected to satisfy the experiment requirements.
Data acquisition. The controller was equipped with the PID control law to achieve the motion. The trajectory was derived from the abovementioned optimization results (i.e., 11 harmonic functions per joint). During the actual movement, the sampling interval was 0.01 s, and the duration of one period was 20 s. To improve the signal-to-noise ratio, each trajectory continued to move for 320 s, and data of 15 consecutive periods were used for mean filtering. Torque information was collected using torque sensors, position information using a motor incremental encoder, and the velocity and acceleration information through Fourier trajectory differential analysis. All signals were processed through a second-order low-pass digital Butterworth filter to process, and the cutoff frequency was 5 Hz.
Identification. The WLS estimation method was used for the parameter identification of all experimental trajectories. The weights were derived from the noise characteristics of the random trajectory identification results, as described in Section 4.2. The position deviation was close to 0, thereby satisfying the zero-noise hypothesis.
Model validation. The accuracy of the obtained parameter estimates could be verified for a different validation trajectory by comparing the measured torque and the predicted torque based on the model and the measured position data.

5. Results

This section provides the results of the experiments described above. In Section 4.2, 12 correlation terms involved in the regression matrix are considered for the experimental design, and the 12 correlation terms are used to identify the dynamic parameters. The identification results and their relative deviations (σi) are summarized in Table 4. The relative standard deviation is used to establish a statistical analysis and characterize the quality of parameter estimation as mentioned in [40].
The results presented in Table 4 show that the mean values of the relative standard deviation are lower in each parameter than in other criteria when the F2 criterion is used because this method finds the estimation with a minimal parameter variance. By contrast, the relative standard deviation mean of each parameter estimated through the multi-criteria embedded optimization method (e.g., F4) is lower than that based on F1 and F3. However, F2 leads to a higher condition number than other methods, and the identification results are sensitive to noise, as displayed in Table 5. The experiments presented in Table 5 are deliberately designed to increase or decrease the torque residual noise. This process is allowed because obtaining fully determined noise is infeasible in the actual experiment. The results show that F4 is more robust than F2. In addition, the relative standard deviation of individual parameters obtained through each method is relatively high because the cited method fails to handle the unmodeled dynamics properly and consider the analytic derivative obtained after the trajectory reparameterization, thus resulting in the difference with the actual parameters. Evidently, all estimated parameters have physically feasible values.
Parameters obtained through the multi-criteria embedded optimization are used to generate dynamic trajectories, and then the measured torques are compared with the predicted torques, as displayed in Figure 7. Compared with the identification results of the optimized trajectory presented in Figure 4, the residual torque value is smaller than that of the latter, that is, σ1 = 0.0516 N2m2, and σ2 = 0.0161 N2m2. Despite a peak error in the predicted torque curve when the speed or acceleration is reversed, this result is reasonable due to assembly clearance or flexibility. In addition, to verify the accuracy of the identification model, an optimization trajectory that is different from F4 is selected in Section 4.2 for verification. The verification results are illustrated in Figure 8. Notably, both curves are close. Table 6 presents the root mean square (RMS) of the measured and predicted residuals.

6. Discussion and Conclusions

To achieve an ideal model control effect and improve human–robot interaction abilities, an accurate dynamic model must be designed. The experimental design, as an important part of parameter identification, must optimize the excitation trajectories. In this research, some optimization criteria for dynamic parameter identification experiments are evaluated. The results indicated that adopting only F1 or F2 criteria will significantly affect the optimization results of the opposition criterion. Although the multi-objective optimization criterion has solved the problem, this criterion was significantly affected by selecting the goal values. Combined with the advantages and disadvantages of each criterion, the multi-criteria embedded optimization method was innovatively proposed. The main contribution is to reduce the DOF of the goal values by adding the criterion to the nonlinear optimization constraints. Furthermore, the optimization results are converted to be mainly restricted by a single time. Considering the motion coupling characteristics of the joint with multiple DOFs, a 2-DOF platform was selected. Based on the optimization simulation and actual experiments, the results showed that the proposed method has achieved the rationalization of the multi-criteria optimization compared with only a single criterion optimization, and the robustness was higher than the multi-objective optimization criterion. The verification experiments presented that the identification results are considerable, and the proposed approach can be considered a suitable procedure for designing the exciting trajectories and improving the results of parameter identification. This method is not only limited to serial robots but also applicable to parallel and exoskeleton robots.
Most existing parameter identification methods focus on off-line identification. However, considering the mechanical wear, aging, temperature, and other nonlinear factors during the use of robots, real-time online identification has become a trend. Therefore, as mentioned in the literature [15,19,41,42,43,44,45], various optimization algorithms and intelligent control theories for improving identification effects must be combined to enable the robots to have the ability of autonomous identification and learning to achieve an efficient and accurate motion control, which have wide applications. In addition, future work will be expanded to include additional DOFs, such as the commonly used 6-DOF or redundant arms, rather than only a 2-DOF manipulator. Furthermore, non-planar motion planning and human–robot interaction control will be realized, and the nonlinear factors, such as flexibility and assembly clearance, will be further investigated.

Author Contributions

J.J. conceived the original ideas, designed all the experiments, and subsequently drafted the manuscript. M.Z. provided supervision to the project, formal analysis, and resources. X.Z. contributed to the construction of the experiment platform and reviewed the writing. H.Z. conducted all the experiments and provided human resources. J.Z. provided supervision and funding support for the project.

Funding

This research was supported by the Major Research Plan (Grant No. 91648201), the National Natural Science Foundation of China (Grant No. 51575157), the Self-Planned Task of the State Key Laboratory of Robotics and System (HIT) (Grant No. SKLRS201816B), and Hebei University of Technology Joint Training Program (Grant No.2017GN0008).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ayusawa, K.; Venture, G.; Nakamura, Y. Identification of humanoid robots dynamics using floating-base motion dynamics. In Proceedings of the 2008 IEEE/RSJ International Conference on Intelligent Robots and Systems, Nice, France, 22–26 September 2008; pp. 2854–2859. [Google Scholar] [CrossRef]
  2. Wu, W.; Zhu, S.; Wang, X.; Liu, H. Closed-loop dynamic parameter identification of robot manipulators using modified Fourier series. Int. J. Adv. Robot. Syst. 2012, 9, 29. [Google Scholar] [CrossRef]
  3. Gautier, M.; Khalil, W.; Restrepo, P.P. Identification of the dynamic parameters of a closed loop robot. In Proceedings of the 1995 IEEE International Conference on Robotics and Automation, Nagoya, Japan, 21–27 May 1995; pp. 3045–3050. [Google Scholar] [CrossRef]
  4. Gautier, M. Dynamic identification of robots with power model. In Proceedings of the International Conference on Robotics and Automation, Albuquerque, NM, USA, 25–25 April 1997; pp. 1922–1927. [Google Scholar] [CrossRef]
  5. Gautier, M.; Janot, A.; Vandanjon, P.O. A new closed-loop output error method for parameter identification of robot dynamics. IEEE Trans. Control. Syst. Technol. 2013, 21, 428–444. [Google Scholar] [CrossRef]
  6. Wu, J.; Wang, J.; You, Z. An overview of dynamic parameter identification of robots. Robot. Comput. Integr. Manuf. 2010, 26, 414–419. [Google Scholar] [CrossRef]
  7. Bona, B.; Curatella, A. Identification of industrial robot parameters for advanced model-based controllers design. In Proceedings of the 2005 IEEE International Conference on Robotics and Automation, Barcelona, Spain, 18–22 April 2005; pp. 1681–1686. [Google Scholar] [CrossRef]
  8. Bahloul, A.; Tliba, S.; Chitour, Y. Dynamic Parameters Identification of an Industrial Robot with and Without Payload. IFAC-PapersOnLine 2018, 51, 443–448. [Google Scholar] [CrossRef]
  9. Briot, S.; Gautier, M. Global identification of joint drive gains and dynamic parameters of parallel robots. Multibody Syst. Dyn. 2015, 33, 3–26. [Google Scholar] [CrossRef]
  10. Swevers, J.; Ganseman, C.; Tukel, D.B.; Schutter, J.D.; Brussel, H.V. Optimal robot excitation and identification. IEEE Trans. Robot. Autom. 1997, 13, 730–740. [Google Scholar] [CrossRef] [Green Version]
  11. Zak, G.; Benhabib, B.; Fenton, R.G.; Saban, I. Application of the weighted least squares parameter estimation method to the robot calibration. IEEE/ASME J. Mech. Des. 1994, 116, 890–893. [Google Scholar] [CrossRef]
  12. Poignet, P.; Gautier, M. Comparison of weighted least squares and extended kalman filtering methods for dynamic identification of robots. In Proceedings of the 2000 IEEE International Conference on Robotics and Automation, San Francisco, CA, USA, 24–28 April 2000; Volume 4, pp. 3622–3627. [Google Scholar] [CrossRef]
  13. Janot, A.; Vandanjon, P.-O.; Gautier, M. A generic instrumental variable approach for industrial robot identification. IEEE Trans. Control Syst. Technol. 2014, 22, 132–145. [Google Scholar] [CrossRef]
  14. Janot, A.; Vandanjon, P.O.; Gautier, M. An instrumental variable approach for rigid industrial robots identification. Control Eng. Pract. 2014, 25, 85–101. [Google Scholar] [CrossRef]
  15. Brunot, M.; Janot, A.; Young, P.C.; Carrillo, F. An improved instrumental variable method for industrial robot model identification. Control Eng. Pract. 2018, 74, 107–117. [Google Scholar] [CrossRef]
  16. Jiang, Z.H.; Ishida, T.; Sunawada, M. Neural network aided dynamic parameter identification of robot manipulators. In Proceedings of the 2006 IEEE International Conference on Systems, Man and Cybernetics, Taipei, Taiwan, 8–11 October 2006; pp. 3298–3303. [Google Scholar] [CrossRef]
  17. Calafiore, G.; Indri, M.; Bona, B. Robot dynamic calibration: Optimal excitation trajectories and experimental parameter estimation. J. Robot. Syst. 2001, 18, 55–68. [Google Scholar] [CrossRef]
  18. Soewandito, D.B.; Oetomo, D.; Ang, M.H., Jr. Neuro-adaptive motion control with velocity observer in operational space formulation. Robot. Comput.-Int. Manuf. 2011, 27, 829–842. [Google Scholar] [CrossRef]
  19. Hu, J.; Xiong, R. Contact force estimation for robot manipulator using semiparametric model and disturbance Kalman filter. IEEE Trans. Ind. Electron. 2018, 65, 3365–3375. [Google Scholar] [CrossRef]
  20. Urrea, C.; Pascal, J. Design, simulation, comparison and evaluation of parameter identification methods for an industrial robot. Comput. Electr. Eng. 2018, 67, 791–806. [Google Scholar] [CrossRef]
  21. Swevers, J.; Verdonck, W.; De Schutter, J. Dynamic model identification for industrial robots. IEEE Control Syst. Mag. 2007, 27, 58–71. [Google Scholar] [CrossRef]
  22. Atkeson, C.G.; An, C.H.; Hollerbach, J.M. Estimation of inertial parameters of manipulator loads and links. Int. J. Robot. Res. 1986, 5, 101–119. [Google Scholar] [CrossRef]
  23. Park, K.J. Fourier-based optimal excitation trajectories for the dynamic identification of robots. Robotica 2006, 24, 625–633. [Google Scholar] [CrossRef]
  24. Gautier, M.; Khalil, W. Exciting trajectories for the identification of base inertial parameters of robots. Int. J. Robot. Res. 1992, 11, 362–375. [Google Scholar] [CrossRef]
  25. Jin, J.; Gans, N. Parameter identification for industrial robots with a fast and robust trajectory design approach. Robot. Com.-Int. Manuf. 2015, 31, 21–29. [Google Scholar] [CrossRef]
  26. Khan, Z.A.; Chacko, V.; Nazir, H. A review of friction models in interacting joints for durability design. Friction 2017, 5, 1–22. [Google Scholar] [CrossRef] [Green Version]
  27. Khalil, W.; Dombre, E. Modeling, Identification and Control of Robots, 3rd ed.; Taylor and Francis Group: New York, NY, USA, 2002. [Google Scholar]
  28. Gautier, M. Numerical calculation of the base inertial parameters of robots. J. Robot. Syst. 1991, 8, 485–506. [Google Scholar] [CrossRef]
  29. Gautier, M.; Khalil, W. Direct calculation of minimum set of inertial parameters of serial robots. IEEE Trans. Robot. Autom. 1990, 6, 368–373. [Google Scholar] [CrossRef] [Green Version]
  30. Mayeda, H.; Yoshida, K.; Osuka, K. Base parameters of manipulator dynamic models. IEEE Trans. Robot. Autom. 1990, 6, 312–321. [Google Scholar] [CrossRef]
  31. Calanca, A.; Capisani, L.M.; Ferrara, A.; Magnani, L. MIMO closed loop identification of an industrial robot. IEEE Trans. Contr. Syst. Techol. 2011, 19, 1214–1224. [Google Scholar] [CrossRef]
  32. Do Thanh, T.; Kotlarski, J.; Heimann, B.; Ortmaier, T. Dynamics identification of kinematically redundant parallel robots using the direct search method. Mech. Mach. Theory 2012, 55, 104–121. [Google Scholar] [CrossRef]
  33. Armstrong, B. On finding exciting trajectories for identification experiments involving systems with nonlinear dynamics. Int. J. Robot. Res. 1989, 8, 28–48. [Google Scholar] [CrossRef]
  34. Presse, C.; Gautier, M. New criteria of exciting trajectories for robot identification. In Proceedings of the 1993 IEEE International Conference on Robotics and Automation, Atlanta, GA, USA, 2–6 May 1993; pp. 907–912. [Google Scholar] [CrossRef]
  35. Abdellatif, H.; Heimann, B.; Grotjahn, M. Statistical approach for bias-free identification of a parallel manipulator affected by large measurement noise. In Proceedings of the 44th IEEE Conference on Decision and Control, Seville, Spain, 15 December 2005; pp. 3357–3362. [Google Scholar] [CrossRef]
  36. Abdellatif, H.; Heimann, B. Advanced model-based control of a 6-DOF hexapod robot: A case study. IEEE/ASME Trans. Mech. 2010, 15, 269–279. [Google Scholar] [CrossRef]
  37. Díaz-Rodríguez, M.; Iriarte, X.; Mata, V.; Ros, J. On the experiment design for direct dynamic parameter identification of parallel robots. Adv. Robot. 2009, 23, 329–348. [Google Scholar] [CrossRef]
  38. Craig, J.J. Introduction to Robotics: Mechanics and Control, 3rd ed.; Pearson Education: Chennai, India, 2009. [Google Scholar]
  39. Mata, V.; Farhat, N.; Díaz-Rodríguez, M.; Manipilators, P. Dynamic parameter identification for parallel manipulators. In Parallel Manipulators, Towards New Applications; IntechOpen: London, UK, 2008. [Google Scholar]
  40. Khalil, W.; Dombre, E. Modeling, identification and control of robots. In Butterworth-Heinemann; IntechOpen: London, UK, 2004. [Google Scholar]
  41. Valdez, F.; Melin, P.; Castillo, O. An improved evolutionary method with fuzzy logic for combining particle swarm optimization and genetic algorithms. Appl. Soft Comput. 2011, 11, 2625–2632. [Google Scholar] [CrossRef]
  42. David, R.C.; Precup, R.E.; Petriu, E.M.; Rădac, M.-B.; Preitl, S. Gravitational search algorithm-based design of fuzzy control systems with a reduced parametric sensitivity. Inf. Sci. 2013, 247, 154–173. [Google Scholar] [CrossRef]
  43. Jin, H.; Liu, Z.; Zhang, H.; Liu, Y.; Zhao, J. A Dynamic Parameter Identification Method for Flexible Joints Based on Adaptive Control. IEEE-ASME Trans. Mech. 2018, 23, 2896–2908. [Google Scholar] [CrossRef]
  44. Saadat, J.; Moallem, P.; Koofigar, H. Training Echo State Neural Network Using Harmony Search Algorithm. Int. J. Artif. Intell. 2017, 15, 163–179. [Google Scholar]
  45. Vrkalovic, S.; Lunca, E.C.; Borlea, I.D. Model-free sliding mode and fuzzy controllers for reverse osmosis desalination plants. Int. J. Artif. Intell. 2018, 16, 208–222. [Google Scholar]
Figure 1. Overall procedure of parameter identification.
Figure 1. Overall procedure of parameter identification.
Sensors 19 02248 g001
Figure 2. Pml identification schematic of weighted least squares (WLS).
Figure 2. Pml identification schematic of weighted least squares (WLS).
Sensors 19 02248 g002
Figure 3. Prototype and simplified model of the two degrees-of-freedom (2-DOF) robot.
Figure 3. Prototype and simplified model of the two degrees-of-freedom (2-DOF) robot.
Sensors 19 02248 g003
Figure 4. Random trajectory identification results for obtaining the noise characteristics.
Figure 4. Random trajectory identification results for obtaining the noise characteristics.
Sensors 19 02248 g004
Figure 5. Condition number versus trajectory points for the actual 2-DOF manipulator.
Figure 5. Condition number versus trajectory points for the actual 2-DOF manipulator.
Sensors 19 02248 g005
Figure 6. Optimal trajectory found through multi-criteria embedded optimization.
Figure 6. Optimal trajectory found through multi-criteria embedded optimization.
Sensors 19 02248 g006
Figure 7. Measured, predicted, and residual torques for the two joints.
Figure 7. Measured, predicted, and residual torques for the two joints.
Sensors 19 02248 g007
Figure 8. Measured position and torque and predicted and residual torques of the validation experiment.
Figure 8. Measured position and torque and predicted and residual torques of the validation experiment.
Sensors 19 02248 g008
Table 1. Some criteria for the experiment optimization design.
Table 1. Some criteria for the experiment optimization design.
No.CriteriaReferencesFrameworks
1 1 λ min ( F T F ) Armstrong [33]deterministic
2 C o n d ( F ) Gautier and Khali [24]deterministic
3 ω 1 C o n d ( F ) + ω 2 F max F min Gautier and Khali [24]deterministic
4 ω 1 C o n d ( F ) + ω 2 1 λ min ( F T F ) Presse and Gautier [34]deterministic
5 log ( det ( F T Σ 1 F ) ) Swevers et al. [10]statistical (D-optimality)
6 1 g = 1 N b F g s Jingfu et al. [25]deterministic
(Hadamard’s inequality)
7 { c o n d ( F ) ω 1 γ F 1 * log ( det ( F T Σ 1 F ) ) ω 2 γ F 2 * Miguel et al. [37]statistical
Table 2. Values of the objective function for the exciting trajectories.
Table 2. Values of the objective function for the exciting trajectories.
TrajectoriesOptimization Criteria
F1 = Cond(F)F2 = −log(det(FTΣ−1F)−1)F3 = multi-objective
F1F2F1F2F1F2
18.9564−46.806514.0528−75.34677.9318−75.6224
24.5228−56.71279.9266−74.19397.8595−76.3168
35.4290−53.19339.236−69.62147.9149−75.7974
44.0043−58.189110.8305−77.03347.9156−75.9591
Table 3. Influence of goal value selection on the optimization results.
Table 3. Influence of goal value selection on the optimization results.
TrajectoriesGoal Weight Values
F3ω1 = 5, ω2 = −70F4ω2 = −75F4ω2 = −70
F1F2F1F2F1F2
15.0354−69.51806.9067−75.00356.5072−70.01
24.9852−70.20787.8303−75.11075.0119−70.0229
35.0671−69.05458.2943−75.00185.5365−70.0113
45.0431−69.3967.7889−75.00084.8974−70.0218
Table 4. Identified dynamic parameters and their respective deviations of the 2-DOF robot.
Table 4. Identified dynamic parameters and their respective deviations of the 2-DOF robot.
IndexF1F2F3F4
pbσipbσipbσipbσi
10.01571.990.02350.940.11610.560.10880.81
20.08940.20.02410.220.02160.950.01261.3
30.00239.50.00174.50.03630.570.03290.86
44.15140.0093.95750.0114.08470.0364.05300.049
51.21680.0181.22570.0171.23570.081.26410.11
60.12530.340.09180.380.01466.910.07421.85
70.19820.0870.06750.460.13900.910.25110.69
80.11230.370.04290.550.02733.370.02594.84
90.06920.210.10720.230.07671.650.09301.84
mean 1.41 0.81 1.67 1.37
Table 5. Comparison of sensitivity with the noise of F2 and F4.
Table 5. Comparison of sensitivity with the noise of F2 and F4.
Noise (σ1, σ2)pb
F2(5.0832,0.0499)0.02350.02410.00173.95751.22570.09180.06750.04290.1072
(0.0516,0.0161)0.04580.00650.01893.98641.21450.03400.09720.05330.1111
(8,1)0.03840.01230.01323.97761.21610.05260.08760.04860.1083
F4(5.0832,0.0499)0.10880.01260.03294.05301.26410.07420.25110.02590.0930
(0.0516,0.0161)0.10720.01000.03134.05301.24920.07470.25090.02800.0923
(8,1)0.10790.01130.03214.05281.25790.07450.25100.02690.0972
Table 6. Root mean square (RMS) of the measured and predicted residuals (Nm).
Table 6. Root mean square (RMS) of the measured and predicted residuals (Nm).
JointsF4Validation
Joint10.2270.235
Joint20.1270.195

Share and Cite

MDPI and ACS Style

Jia, J.; Zhang, M.; Zang, X.; Zhang, H.; Zhao, J. Dynamic Parameter Identification for a Manipulator with Joint Torque Sensors Based on an Improved Experimental Design. Sensors 2019, 19, 2248. https://doi.org/10.3390/s19102248

AMA Style

Jia J, Zhang M, Zang X, Zhang H, Zhao J. Dynamic Parameter Identification for a Manipulator with Joint Torque Sensors Based on an Improved Experimental Design. Sensors. 2019; 19(10):2248. https://doi.org/10.3390/s19102248

Chicago/Turabian Style

Jia, Jidong, Minglu Zhang, Xizhe Zang, He Zhang, and Jie Zhao. 2019. "Dynamic Parameter Identification for a Manipulator with Joint Torque Sensors Based on an Improved Experimental Design" Sensors 19, no. 10: 2248. https://doi.org/10.3390/s19102248

APA Style

Jia, J., Zhang, M., Zang, X., Zhang, H., & Zhao, J. (2019). Dynamic Parameter Identification for a Manipulator with Joint Torque Sensors Based on an Improved Experimental Design. Sensors, 19(10), 2248. https://doi.org/10.3390/s19102248

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop