A Two-Step Method for Dynamic Parameter Identification of Indy7 Collaborative Robot Manipulator

Accurate dynamic model is critical for collaborative robots to achieve satisfactory performance in model-based control or other applications such as dynamic simulation and external torque estimation. Such dynamic models are frequently restricted to identifying important system parameters and compensating for nonlinear terms. Friction, as a primary nonlinear element in robotics, has a significant impact on model accuracy. In this paper, a reliable dynamic friction model, which incorporates the influence of temperature fluctuation on the robot joint friction, is utilized to increase the accuracy of identified dynamic parameters. First, robot joint friction is investigated. Extensive test series are performed in the full velocity operating range at temperatures ranging from 19 °C to 51 °C to investigate friction dependency on joint module temperature. Then, dynamic parameter identification is performed using an inverse dynamics identification model and weighted least squares regression constrained to the feasible space, guaranteeing the optimal solution. Using the identified friction model parameters, the friction torque is computed for measured robot joint velocity and temperature. Friction torque is subtracted from the measured torque, and a non-friction torque is used to identify dynamic parameters. Finally, the proposed notion is validated experimentally on the Indy7 collaborative robot manipulator, and the results show that the dynamic model with parameters identified using the proposed method outperforms the dynamic model with parameters identified using the conventional method in tracking measured torque, with a relative improvement of up to 70.37%.


Introduction
The control of robotic systems mainly depends on accurate dynamic models to increase performance. The topic has been covered in recent studies as it relates to control approaches (such as sliding mode control, model predictive control, passivity-based control, and state-space control) for a variety of systems [1][2][3][4][5] having common explicit dynamic models in the control loop. The accuracy of the dynamic parameters that characterize the dynamic model is the most important aspect affecting the dynamic model's accuracy. Such dynamic parameters are typically unknown or only partially known via nominal parameters approximately calculated from computer-aided design (CAD) software. Because of manufacturing limitations, such parameters' accuracy cannot be guaranteed [6]. An identification method is an alternate approach for estimating dynamic parameters (such as mass, the center of mass, and moment of inertia) employing special experiments. This method has been used and researched extensively in [6][7][8].
The most prevalent strategy for experimental identification methods is a combined Inverse Dynamics Identification Model and Least Squares (IDIM-LS) method. The accuracy of parameter estimations for such an approach is often affected by measurement noise and modeling errors. Significant research efforts have been made on identifying appropriate exciting trajectories to overcome measurement noise concerns [9][10][11]. However, the majority of modeling errors that results in considerable deterministic structural errors are caused by disregarding complex and nonlinear joint dynamics effects (such as friction), which cannot be explained by random variables. For example, strain-wave transmission methods like the Harmonic Drive (HD), which are frequently used in collaborative robot manipulators because of their benefits of high-quality characteristics of loading capability and lightweight, can cause nonlinear joint dynamics.
Friction is one of the most prominent undesired nonlinear phenomena inherent in HD, introducing modeling errors and affecting the accuracy of the dynamic model. Although assuming the friction model as Coulomb friction plus viscous friction linear to the robot joint velocity [12] is a common method for dynamic model identification, it is generally insufficient in practice. Furthermore, comparing two tests, one at robot startup and one after the thermal transient owing to friction change with robot joint temperature variation, will result in a model discrepancy. In general, since friction is influenced by temperature [13][14][15][16][17], improving the accuracy of the friction force model will significantly improve the accuracy of the robot dynamic model.
The present research aims to realize the dynamic parameter identification of a manipulator as accurately as possible. The primary contribution is using a reliable dynamic friction model to increase the accuracy of the identified dynamic parameters. Considering the influence of robot joint temperature variation on robot joint friction, the proposed method aims to increase the accuracy of dynamic parameters by minimizing non-structural uncertainty through the use of a more accurate dynamic friction model. The important contributions that this work has made can be summed up as follows: • a linear parametrization that describes the dynamics of the Indy7 collaborative robot manipulator has been developed, • the minimum set of base parameters needed to characterize the dynamic behavior of the Indy7 robot has been successfully approximated, • during robot dynamic parameter identifications, a collaborative robot's inherently nonlinear joint dynamics were taken into account, • validation of the accuracy of the torque predictions made by the identified robot dynamic model.
The remainder of the paper is organized as follows: Section 2 presents the mathematical model of the robot manipulator and the development of a linear identification model and friction model. Section 3 discusses the identification experiment process and compares the accuracy between two dynamic models with parameters identified using conventional and proposed linear identification models. After analyzing validation results, the research outcome is discussed and summarized in Sections 4 and 5, respectively.

Dynamic Modeling
The inverse dynamic model of a rigid robot composed of N moving links calculates the motor torque vector τ m (the control input) as a function of the generalized coordinates. It can be obtained from the Lagrangian or Newton Euler equation as recalled in [18,19]: where q ∈ R N represents the N joint angles, M(q) is the inertia matrix, C(q,q) captures the centrifugal and Coriolis effects, and G(q) represents the gravitational torque. The drive torque τ m applied by the joint motor controls the system.

Linear Identification Model
Using the barycentric parameters [20] or the modified Newton-Euler parameters [21] a model of the type in Equation (1) may also be parameterized and transformed into a model that is linear in a new set of unknown dynamic parameters Φ s .
where Y s (q,q,q) ∈ R N×11N is the observation matrix, which depends only on the motion data, and Φ s ∈ R 11N×1 is a vector of standard dynamic parameters: For each link i, it comprised the six components of the inertia tensor, (Ixx i , Ixy i , Ixz i , Iyy i , Iyz i , Izz i ); the three components of the first moment, (mx i , my i , mz i ); the mass, m i ; and the total inertia moment for rotor actuator and gears, Ia i .
Because friction is the primary element influencing dynamic model accuracy, a comprehensive dynamic model must contain friction in addition to the fundamental dynamic mathematical model Equation (1). In the conventional dynamic parameter identification method [12], dynamic parameters have been expanded as shown in Equation (4) to include the most familiar static friction model Equation (3).

The Proposed Linear Identification Model
Accurate dynamic models are frequently restricted to identifying important system parameters and compensating for nonlinear terms. Friction, as a primary nonlinear element in robotics, has a significant influence on model accuracy, particularly in the case of velocity reversal and low-speed disturbance, which are difficult to precisely represent. The conventional linear identification approach employs a simplified static friction model that cannot fully capture the dynamic behavior of the friction, such as the non-linearity of viscous friction seen in most HD and the temperature dependency of friction. Especially, augmented dynamic parameters in Equation (4) are highly affected by the robot joint temperature conditions. A significant variation can be exhibited in the values of these parameters depending on the temperature (hot or cold) condition of robot joints during the identification process. The experimental results presented in Table A2 in Appendix A show that for the conventional method, a significant variation in the identified dynamic parameters can be observed (highlighted part shows parameters with more than 20% relative change compared to the cold condition for conventional identification method) depending on the robot joint temperature condition during the identification process.
Therefore, to increase the robustness of conventional identification, we proposed a two-step identification procedure. In the first step, parameters of a comprehensive friction model Equation (16), which was reported in our previous work [13], are identified for each robot joint. Using the identified comprehensive friction model parameters, friction torque is computed for the measured trajectory and joint temperature during dynamic parameter identification. After subtracting friction torque from the recorded joints' torque, a modified conventional linear identification model Equation (6) is used in the second step of the dynamic parameter identification procedure.
Friction torques are considered a sum of estimates and error terms. With a reasonable assumption that the friction errorτ s τ s −τ s contain Coulomb, viscous, and offset friction contributions; that is,τ s =F c sign(q) +F vq +F o , Equation (4) is modified by replacing all friction related parameters with friction error model parameters such asF c,i ,F v,i , andF o,i for i = 1, · · · , N.
Thus, the proposed linear identification model is given by and τ m f = τ m − τ f . τ f is the friction torque vector of robot joints computed based on Eq. (16). The proposed identification method is more robust, and the variance in the identified parameters can be reduced to 20% or less in relative change as shown in Table A2 in Appendix A.

The Minimum Inertial Parameters
Not all 11N inertial parameters of a robot manipulator can be identified; some inertial parameters have no influence on the dynamic model, while others have an effect only in linear combinations. As a result, the set of standard dynamic parameters Φ s to be identified can be reduced to a minimal set of parameters known as base parameters, Φ b ∈ R b×1 , where b is the number of base parameters. These parameters can be obtained from the standard inertial parameters by eliminating those which have no effect on the dynamic model and by regrouping some others in linear relations [22][23][24]. Following numerical approach [24], a set of base parameters are computed for the Indy7 robot manipulator with coordinate systems placed as illustrated in Figure 1 and with the modified Denavit-Hartenberg (mDH) parameters in Table 1. The result is presented in Table A1 in Appendix A.
Therefore, the joint torque can be expressed as where Y B (q,q,q) ∈ R N×b is the observation matrix with full rank, and Φ B ∈ R (b+3N)×1 is base parameters,Φ b , augmented with, friction parameters for the conventional case or friction error parameters for the proposed case.

Friction Modeling
In prior work [13], we established a comprehensive friction model that included nonlinear viscous and temperature dependence. A basic overview of the formulation of the comprehensive friction model will be offered here (for further information, please see [13]).
Numerous static and dynamic friction models [9,25,26] have been investigated to characterize friction behavior. Unlike static friction models, which are discontinuous during velocity reversal, the friction phenomena have a nonlinear continuous behavior at velocity zero-crossing. As a result, for a smooth transition and a more accurate representation of the friction phenomenon, it is preferable to adopt a dynamic friction model. A dynamic friction model τ f in its generic form that takes into account internal state z (potentially multiple states), velocity, and temperature-dependent effects can be described by where δ(z i ) is the transient velocity response function, ζ(q i , T i ) is the temperature,T i , dependent velocity strengthening function including viscous friction. As all consideration, unless otherwise mentioned, is for a single joint, we shall hereafter use a plain symbol without the subscript i. Using a first-order differential equation of a generic form, the dynamics of the internal states z may be expressed as follows: with G(·) a general nonlinear function. Since the steady-state friction is a function of relative velocity only (G(z,q, q) = 0), the generic form of dynamic friction model Equation (8) converges to a static friction model Equation (10).
where the generic form of a velocity weakening function (Stribeck effect) g(q) is given by where F c is Coulomb friction, F s is stiction, and v s is Stribeck velocity. The LuGre friction model [26] is a particular instance of the generic form dynamic friction model illustrated in Equation (8), in which viscous friction is depicted as a linear proportional to relative velocity. Considering that the model is based on the average behavior of the contact between bristles, the model of the average deflection z can be given by where, σ 0 is the bristle stiffness. Therefore, the expression for the transient velocity response where σ 1 is the micro-damping coefficient.

Nonlinear Viscous Friction Model
The most typical friction contacts in robotic actuation are the lubricated friction contacts in the main HD. In contrast to the typically employed linear viscous friction, the robot actuators with HD exhibit significantly nonlinear viscous friction. In fact, changes in joint temperature cause a considerable change in viscous friction for the same joint speed. As a result, the viscous friction model is now a function of joint temperature as well as velocity.

Temperature Dependent Viscous Friction Model
At high joint temperatures, friction in HD tends to decrease considerably due to changes in contact point characteristics and lubricating layer thickness. In contrast, at low joint temperatures, the viscosity of the lubricant increases, leading to an increase in viscous friction. Considering such dependencies of viscous friction on temperature, we proposed an exponential model that maps the nonlinear viscous friction behavior over the joint velocity range, Equation (14). To account for temperature dependency, the parameters of the nonlinear viscous friction model are further parameterized as a function of joint temperature Equation (15). where

Comprehensive Friction Model
Following a thorough analysis of the nonlinear viscous and temperature dependencies, a comprehensive friction model is developed. Taking Equations (13) and (14) into account, a complete friction model can be developed as follows:

Figure 2.
Indy7 collaborative robot manipulator utilized for experimental verification, along with an exploded layout of its joint module, which includes an electric motor, harmonic drive, motor driver, encoder, and temperature sensor.
The system is actuated using a brushless DC motor with HD gear transmission. A 16-bit multi-turn absolute encoder measures each angular position. The joint torque is determined by multiplying the measured motor current by the torque constant of the motor. For each joint, a temperature sensor of type MCP9808 is placed on the housing surface of the HD in order to acquire the joint temperature. All data are logged at a sampling frequency of 4 kHz. Data is sampled at times t(k) = kT s , k = 1, 2, · · · , M, where T s = 0.25 ms is the sampling period.
The identification procedure is illustrated schematically in Figure 3. The identification procedure can generally be divided into two steps. In general, the identification procedure can be classified into two steps. First, the friction of each joint is identified. The parameters of the comprehensive friction model in Equation (16) are identified. After identifying the parameters of the comprehensive friction model, the next step is to identify the parameters of the robot's dynamic model. Estimated friction torque based on a comprehensive friction model is subtracted from measured torque during robot dynamic model parameter identification τ m . This non-friction torque τ m f will be used to identify the robot's dynamic model parameters based on the proposed identification method. These identification procedures will be described in full in the sections that follow.

Static Model Identification
An extensive experimental investigation was carried out to gain insights and distinguish relationships in the velocity and temperature during the identification of friction model parameters. To begin, the joint is rotated back and forth at various steady rates. We may presume that the centrifugal and inertial torques are negligible since only a single joint is moving at a time and data is collected for constant velocity region. As a result, the applied torque may be defined as the total of gravity-induced and friction torques, Equation (17), As all consideration, unless otherwise mentioned, is for a single joint, we shall hereafter use a plain symbol without the subscript i. An estimate of friction that is independent of direction can be obtained by taking into account forward and reverse motions for a speed level ofq. Since gravity-induced torque is directional, it may be canceled out, and the friction torque can be calculated as follows, where τ + m and τ − m are the resulting torques when the joint is moved forward and backward, respectively. The procedure can be repeated for severalq's and a friction curve can be drawn, which contains steady-state friction values. To observe temperature dependencies, the above procedure was repeated for different constant velocity points covering the considered velocity range over different constant joint temperatures. The temperature deviation during every individual experiment is maintained at ≤1 • C, which is assumed as an acceptable uncertainty.
The collected data are used to identify different temperature-independent friction curves Equation (19) where is viscous friction at constant temperature case. Then, further parameterized viscous friction model parameters, Equation (15), were identified using the estimated parameters in Equation (19) for different joint temperatures (see Figure 4a-c).
It should be noted that the temperature dependence of parameters of the velocity weakening function, g(q), is neglected due to negligible variations with temperature change.
Because the remaining joints of the robot show a similar pattern, only the result for joint one is shown in Figures 4 and 5. Figure 5 depicts the total result of the proposed velocity-temperature model Equation (14). For a more detailed identification procedure refer to [13].

Dynamic Model Identification
To identify dynamic friction model parameters, σ 0 and σ 1 , the approach used in [27] was addapted in this work. In this experiment, parameters for the dynamic friction model, σ 0 and σ 1 , were determined by adapting a method from [27]. Open-loop experiments were conducted with a sinusoidal torque input, and the resulting data were used to estimate dynamic parameters using a pattern search solver in the MATLAB optimization toolbox with a default option. The results of the identified parameters of the dynamic friction model are shown in Table 2. Because the noise in recorded data is the primary source of errors in parameter identification, signal processing is essential for ensuring the quality of the measured data. In our experiment, the measured trajectory data q,q are filtered by a 4th order Butterworth filter in both forward and reverse directions to eliminate lag of the filtered trajectories q,q. The value of angular acceleration is calculated by differential equations since it cannot be measured directly from the Indy7 collaborative robot system. To avoid phase shifts in the differential calculation, angular accelerationq, is computed through a central difference procedure.
Since the sampling frequency is much higher than the frequencies of interest in the dynamics, the data is decimated/down-sampled to reduce the required computational resources. In this experiment, a donwsampling factor of 0.8 f s /(4 f dyn ) = 80 were used, where f s = 4 kHz and f dyn = 10 Hz [28]; that is, every 80th sample is used for parameter estimation. Thus, to remove information-free samples, a down-sampling is conducted on the filtered trajectoriesq,q andq.

Torque Computation
First, the samples k = 1, · · · , M are ordered in the measurement vector Γ i and observation matrix W i for each joint i = 1, · · · , N individually; that is with Y b,i (q k ,q k ,q k ) being the ith row of the regressor evaluated in the kth sample of the filtered trajectory. The filtered data are ordered joint-wise in the measurement vector and observation matrix as Then, the measurement vector Γ is sampled down with the same downsampling factor 0.8 f s /(4 f dyn ) = 80, and the observation matrix W is ordered with Y b,i (q k ,q k ,q k ) being the ith row of the regressor evaluated in the kth sample of the filtered and downsampled trajectory instead. Finally, using the downsampled measurement vector,Γ, and the observation matrix formulated with the corresponding filtered and downsampled trajectories,Ŵ, the base parameters are estimated by solving the weighted least square (WLS) problem: where Λ is a weight matrix given by [29]: with b i being the number of base parameters related to link i.

Excitation Trajectory Design
The goal of the excitation trajectory optimization is to find a trajectory that can sufficiently excite the identified dynamic parameters. The trajectory used for parameter estimation has been designed to optimize the condition number of the regressor matrix, W. Typically, a finite Fourier series method Equation (25), which suppresses measurement noise, is used to generate periodic trajectories [30]. Each kth joint trajectory is defined as a where in this experiment w f = 0.1π is the fundamental frequency of the Fourier series, q io is the offset value of the position trajectory and L = 5 is the number of harmonics. The optimal excitation trajectory is generated, which yields in a condition number of cond(W ) = 80.25, and the result is shown in Figure 6 for a 20 seconds long trajectories.

Physical Feasibility of the Dynamic Parameters
Physical constraints on real robot dynamic parameters such as mass and inertia tensors must be considered to obtain meaningful estimations. To this end, the estimated dynamic parametersΦ B might be possibly physically infeasible (e.g., a negative link mass or not positive definite inertia tensor), and this can be caused, for instance, by modeling error or by noisy measurements. Over the last two decades, extensive study in robotics has been conducted on the idea of imposing physical feasibility as part of a parameter identification process. The physical feasibility of parameters is formulated as constraints on its mass m i , viscous-Coulomb friction F v,i , F c,i , joint module inertia I ai , and link's center-of-mass (CoM) I i parameters as follows [31][32][33]: has to be used instead of I i [18]. Therefore, the set of all physical feasibility constraints for each link i given by Equation (26) can then be written as where S(·) is the skew-symmetric matrix operator and l i is first moment-of-inertia vector, T is the center of mass relative to the link frame i (see [33]). Observing that the second inequality of Equation (28) is in fact the Schur complement of a matrix B L,i (Φ s ) allows rewriting Equation (28) as where 0 n×m , 1 3×3 respectively denote the n × m zero, and 3 × 3 identity matrices. For a robot with N links, the feasibility condition can be expressed as where B(Φ s ) is a single block-diagonal matrix given by After applying a bijective affine transformation to a convex spectrahedron set B = {Φ s ∈ R n : B(Φ s ) 0}, a new physically feasible extended base parameter set, B B , can be defined as where n = 13N, andΦ s,d is a reordered dynamic parameter according to linearly independent columns of the regressor matrix (see [33] for details). Finally, a method for feasible base parameter estimation with ordinary least square (FBPEOLS) formulated in [33] is used to compute the physically feasible base parameter vectorΦ B which minimizes the sum of squared residuals, where

Model Quality Metric
During verification experiments, the confidence of the dynamic model with identified dynamic parameters is verified using specific trajectories. A root-mean-square error (RMSE) is used to evaluate the model quality numerically: for joint i.

Results
For identification of base parameters, a persistent excitation trajectory as depicted in Figure 6 is generated. The observation matrixŴ is constructed with the measured, filtered, and decimated joint trajectories. The measurement vectorΓ is also constructed with the measured and decimated non-friction torques,τ m f . For comparison purposes, base parameters are identified using the conventional method, where the measurement vectorΓ is constructed with the measured and decimated torques,τ m . Then, a method for feasible base parameter estimation with OLS, Equation (33) is applied and the physically feasible base parameters are computed for both conventional (Equation (5)) and proposed (Equation (6)) methods.
The effectiveness of the presented identification method is demonstrated by comparing the accuracy of the dynamic model with parameters obtained using the proposed linear identification model Equation (6) to the accuracy of the dynamic model with parameters obtained using a conventional linear identification model Equation (5). Since we observed significant friction torque variations between cases of low joint temperature and cases of high joint temperature, the accuracy of the dynamic model is evaluated for two of these temperature extremes: (1) at the low joint temperature (19∼22 • C) and (2) at the high joint temperature (40∼46 • C). The scenario in which the joint temperature is low is when the robot begins operation, whereas the condition in which the joint temperature is high is when the robot experiences a thermal transient after a lengthy period of operation. The dynamic model parameters are identified for moderate robot joint temperature (29∼33 • C) condition.
To evaluate the accuracy of the estimated parameters, predicted and measured torques are compared. The experimental findings demonstrate that the estimated joint torque by the identified parameters using the proposed method can be predicted well when compared to the estimated torque using the conventional method as shown in Figures 7 and 8. For instance, from Figure 7b we can observe that the model accuracy improves 74.6% from RMSE of 13.44 Nm to that of 3.42 Nm for joint one at the low joint temperature condition. The model accuracy also improves 82.9% from RMSE of 15.95 Nm to that of 2.72 Nm for joint one at the high joint temperature condition as shown in Figure 8b. In both temperature conditions, the model accuracy is evaluated using the same excitation trajectories employed during parameter identification ( Figure 6). To confirm that the obtained dynamic model is accurate and reliable, new excitation trajectories are generated and the accuracy of the model is evaluated in the subsequent section.

Validation
To validate the results of robot dynamic parameter identification and the relevance of the proposed method, a new robot trajectory different from the trajectory used for identification is generated. Experiments of validation are conducted for three separate reference trajectories, Trajectory A, Trajectory B, and Trajectory C, with respective condition numbers of 70.3, 74.9, and 104.5. A comparison is made between the torques obtained by the dynamic model with parameters obtained using the proposed linear identification model, Equation (6), and the torques obtained by the dynamic model with parameters obtained using the conventional linear identification model, Equation (5). To further validate the accuracy of the estimated parameters, the RMSE of residuals between the measured torque and the estimated torque obtained from the identification model is introduced. The relative improvement in RMSE between the conventional and proposed methods is also examined. The dynamic model's accuracy was validated for two extreme joint temperature scenarios (low and high joint temperature). Tables 3 and 4 show the results of the root mean square error (RMSE) and relative improvement analyses for the low and high joint temperature scenarios, respectively. According to the results, the proposed identification approach has the potential to achieve a relative improvement that ranges from 38% to 78.37% when compared to the conventional one for the low joint temperature scenario. For the scenario where the joint temperature is high, it is possible to achieve a relative improvement of 15.05% to 61.64%. Table 3. RMSE of the dynamic models: with parameters obtained using conventional and proposed linear identification models. Improvements in RMSE relative to the dynamic model with parameters obtained using the conventional linear identification model. The validation experiment is performed at the lower joint temperature (19∼22 • C).  Table 4. RMSE of the dynamic models: with parameters obtained using conventional and proposed linear identification models. Improvements in RMSE relative to the dynamic model with parameters obtained using the conventional linear identification model. The validation experiment is performed at the higher joint temperature (40∼46 • C). To sum up, the dynamic model obtained based on the suggested method can reduce torque prediction error by up to 78.37% and 61.64% in the low and high joint temperature cases, respectively. This proves that the proposed parameter estimation method is effective.

Discussions
Experimental results indicate that identifying the dynamic parameters of the collaborative robot manipulator using the proposed identification approach improves the accuracy of the robot dynamic model. It is important to note that the torque prediction errors for the first and second joints are higher than those for the other joints. Due to their high self-weight, the first and second joint modules include a larger HD with a gear ratio of 121:1 and a motor power rating of 500 W. This causes these two joints to have greater joint friction than the others. For a conventional linear identification model, ignoring such a significant joint nonlinearity due to friction leads to a higher torque prediction error for the first two joints in particular. However, by integrating a more realistic friction model and adopting the proposed identification strategy during robot dynamic model parameter estimation, the proposed method can mitigate these problems and significantly reduce torque prediction errors. During robot dynamic model parameter identification, considering a realistic friction model that reflects the variance due to joint temperature changes over extended operation increases the model's accuracy. From the experimental data, it is obvious that there is a considerable variation in joint friction during robot startup and after the robot's thermal steady state is attained. Since linearity in the parameters might be lost during long time periods of robot operation, such fluctuation in friction owing to thermal factors leads to model mismatch. In contrast to friction characteristics, robot dynamic parameters should be temperature-independent.
Nevertheless, ignoring highly nonlinear and temperature-dependent friction during the identification procedure, which leads to model mismatches over the extended operation, may lead to an incorrect estimation of the robot's dynamic model parameters. This could lead to a less accurate robot dynamic model. As seen in the preceding section, experimental results also indicate that using a more realistic friction model during identification can increase the accuracy of the robot's dynamic model. The proposed method could be limited by the accessibility of robot joint temperature data. In this study, we assumed that each joint has a temperature sensor. Future work may include the development of sensorless joint temperature estimation for collaborative robot manipulators lacking a temperature sensor.

Conclusions
In this study, the problem of identifying the dynamic parameters of collaborative robots with a more realistic friction model is examined. The primary source of nonlinear joint dynamics, such as friction, in a collaborative robot manipulator, is a strain-wave transmission mechanism, such as the HD, which is commonly employed in collaborative robot manipulators due to its high-quality attributes of loading capacity and lightweight. Such nonlinear joint dynamics have a significant impact on the accuracy of the robot's dynamic parameters if not taken into account during the identification procedure. Moreover, the lubricated friction contacts in HD are greatly impacted by the variation in joint temperature, resulting in a noticeable change in friction due to the variation in joint temperature. As a result, the primary objective of this research is to incorporate nonlinear joint dynamics into the identification process of the robot's dynamic model parameters. Since the current robot dynamic model parameter identification method is based on the assumption of the dynamic model's linearity, it is hard to directly incorporate a nonlinear friction model. Therefore, we suggest a strategy consisting of two steps, the first of which is to identify a comprehensive friction model. This friction model is utilized to estimate joint friction as a function of joint velocity and temperature. Then, we proposed an identification technique in which a non-friction torque would be utilized in the robot dynamic parameter identification step. This would replace the practice of utilizing measured joint torque, which would have included friction.
Using the proposed identification approach, the parameters of the Indy7 dynamic model were obtained, and the torque prediction accuracy of the robot's dynamic model was assessed. Experiments of validation demonstrated that the torque estimation of the dynamic model obtained using identified parameters based on the proposed method can accurately represent the robot's dynamic characteristics. The results of the experiments that were done with the validation trajectories show that the proposed method is capable of achieving a reduction in torque prediction error of up to 78.37% in the low joint temperature case and up to 61.64% in the high joint temperature case This proves that the proposed parameter estimation method is effective.
Future work will be expanded to include payload estimation, and use the derived dynamic model for applications such as collision detection and model-based robot controller design. In addition, as a future study, the proposed method will be contrasted and analyzed with more advanced identification algorithms that take joint friction into account.

. Robustness of Identification Methods
The metric relative change is defined as where Φ Bi,c is ith identified base parameter in the cold (19∼22 • C) joint temperature condition and Φ Bi,h is ith identified base parameter in the hot (40∼46 • C) joint temperature condition.