# Influence of the velocity on quasi-static deflections of industrial articulated robots

M.K. Gonzalez, N. Theissen, N. Agirre, Jon Larrañaga, Patxi Hacala

► To cite this version:

M.K. Gonzalez, N. Theissen, N. Agirre, Jon Larrañaga, Patxi Hacala. Influence of the velocity on quasi-static deflections of industrial articulated robots. The International Journal of Advanced Manufacturing Technology, 2023, 125, pp.1429-1438. 10.1007/s00170-022-10661-x. hal-04145406

HAL Id: hal-04145406

https://hal.science/hal-04145406

Submitted on 29 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Influence of the velocity on quasi-static deflections of industrial articulated robots

Monica Katherine Gonzalez¹ ID • Nikolas Theissen¹ • Nora Agirre² • Jon Larrañaga³ • Patxi Hacala⁴,⁵ • Andreas Archenti¹

Received: 25 August 2022 / Accepted: 3 December 2022 / Published online: 7 January 2023
© The Author(s) 2023

## Abstract

This article presents the measurement and analysis of the influence of velocity on the quasi-static deflections of industrial manipulators of three different manufacturers. Quasi-static deflection refers to the deflection of the end effector position of articulated robots during movement at low velocity along a predefined trajectory. Based on earlier reported observations by the authors, there exists a difference in the static and quasi-static deflections considering the same points along a trajectory. This work investigates this difference to assess the applicability of robotic compliance calibration at low velocities. For this assessment, the deflections of three industrial articulated robots were measured at different speeds and loads. Considering the similarity among the robot models used in this investigation, this work also elaborates on the potential influence of the measurement procedure on the measured deflections and its implications for the compliance calibration of articulated robots. For all industrial articulated robots in this investigation, the quasi-static deflections are significantly larger than the static ones but similar in trend. Additionally, the magnitude of the quasi-static deflections presents a proportional relationship to the Cartesian velocity.

**Keywords** Industrial robot · Contact applications · Quasi-statics · Position error · Compliance

## 1 Introduction

Industrial articulated robots are one of the key enablers of modern manufacturing environments. They are characterised by flexibility and lower cost per unit of working

space compared to other specialised machinery such as machining centres. These capabilities have allowed the integration of industrial robots across important industries in a variety of applications, and commodities such as consumer electronics, industrial machinery, or vehicles, since 1961 [1].

Some of the most significant disadvantages of articulated industrial manipulators in contact applications, such as robotic milling [2], are their comparably lower accuracy [3] and higher vibrations [4] in contrast to the benchmark machinery used for these applications. The improvement of positioning accuracy is the subject of the study of robot calibration [5]. Positioning accuracy is the distance between the commanded and the attained position of an articulated robot's end effector (EE) [6]. As described by Mooring [7], robot calibration intends to improve the positioning accuracy by:

1. Modelling the kinematics (imperfect geometries and configurations [8] as well as gear backlash [9]) and non-kinematic parameters (joint and link compliance [10]) and thermo-elasticity errors [11]

2. Measuring the positioning accuracy of the actual robot [12]

Monica Katherine Gonzalez and Nikolas Theissen contributed equally to this work.

✉ Monica Katherine Gonzalez
mgon@kth.se

¹ Department of Production Engineering, KTH Royal Institute of Technology, Brinellvägen 68,114 28, Stockholm, Sweden

² Research and Development, MMAEN S.L., Pol. Plazaola, Manzana E, Nave 2, 31195, Berrioplano, Spain

³ Faculty of Engineering, Mondragon University, Loramendi, 4,20500, Gipuzkoa, Spain

⁴ Dynamics and Control Department, IDEKO Research Centre, Arriaga Kalea, 2, 20870, Elgoibar, Spain

⁵ Robotics Department, ESTIA Institute of Technology, Technopole Izarbel, 64210, Bidart, France
3. Identifying its model parameters from the measure-
   ments [13]

4. Implementing the model parameters to optimise the accuracy of the manipulator [14].

Robot original equipment manufacturers (OEMs) offer working space calibration services to improve the mean and maximum positioning accuracy to typically less than 0.5 mm and 1.0 mm [12] for all the robot working space.

Contact applications [15] such as machining [16] and especially high material removal rate operations [17], require tight workpiece tolerances. They commonly demand an accuracy better than ±0.25 mm on the location of features. This may require calibration for operating spaces rather than the OEM calibration for the working space [18]. This operating space calibration, done for a limited set of manipulator configurations in close vicinity, should yield a local optimum in positioning accuracy.

Garnier and Subrin [19] provided recently a summary of
measurement instruments and procedures used for compli-
ance calibration. All the described measurement procedures
evaluate the industrial manipulator in static configurations
and follow ISO 9283 [6] recommendations. These mea-
surement procedures are less complex and minimise the
uncertainty associated with the manipulator's configuration.
However, when the actual contact application is performed,
geometric and non-geometric (motion dependent) errors
contribute to the robot's positioning accuracy [20]. There-
fore, Theissen et al. [21] proposed a measurement scheme
in which the industrial manipulator moved at velocities
from 50 to 250 mm s⁻¹. In preceding research work, the
authors reported an apparent difference between static and
quasistatic (QS) deflections of industrial manipulators. To
investigate this phenomenon further, this article analyses
and quantifies the effect of Cartesian velocity on the mea-
sured deflections of industrial manipulators to evaluate the
suitability of quasi-static stiffness measurements for robotic
compliance calibration.

This paper is organised as follows: Section 2 discusses the concept and the rationale for measurements under motion. Section 3 summarises the experimental measurement procedures carried out to analyse and quantify the effect of velocity. Section 4 focuses on the measurement results as well as an analysis of the validity of the measurement results from a metrological point of view.

2 Quasi-static compliance

The mechanical stiffness of a system can be defined as its
capacity to sustain loads, which results in a change in its
geometry [22]. The magnitude of a structure’s mechanical
stiffness depends on the frequency of the displacement,

vibration, or excitation force. This magnitude is minimal,
close to or at an eigenfrequency of the mechanical structure.
The inverse of stiffness is named compliance.

This article differentiates between static and quasi-
static deflections. The static deflections are measured at
discretised static configurations along trajectories, while the
quasi-static deflections are measured under motion along
the same trajectories. Thus, in the static state, the wrenches
and the deflections have a frequency component at 0 Hz.
In the quasi-static state, the wrenches and the deflections
have frequency components higher than 0 Hz and lower
than approximately 10% of the first eigenfrequency of
the industrial manipulator [23, 24]. In other words, this
article uses the terminology quasi-static to indicate that
the machine under investigation is moving. At the same
time, the frequency components of the displacement and the
excitation force are at least five to ten times smaller than
the first eigenfrequency of the system. In this frequency
range, it should be feasible to approximate the behaviour
of the mechanical system accurately through static models
instead of dynamic models. This approximation can be
visualised in Fig. 1. The dynamic compliance of the
industrial manipulator remains almost constant close to 0
Hz and up to 5 Hz, while it starts changing and then
increases significantly close to the first eigenfrequency at
approximately 14 Hz.

Other authors have used the term quasi-static to indicate
that the robot [2] or machine tool [25] are moving during
the measurement phase. In that sense, the term quasi-static
compliance refers to the parameterisation of a compliant
manipulator model by using a tuple, which consists of the
configuration in terms of joint angles Θ, the quasi-static

<figure><img src="image_2.png" /><figcaption>Fig. 1 Dynamic compliance of an ABB IRB 6700 measured at [625 –950 1300]mm in RBCS</figcaption></figure>
wrench $W_{QS}$, and the deflection $\Delta X_{QS}$ vector, as follows:

$$ \| K_{\theta}^{-1} J(\Theta)^{T} \cdot W_{QS} - J(\Theta) \Delta X_{QS} = 0 \|^2 \quad (1) $$

Here $K_{\theta}$ denotes the diagonal joint stiffness matrix, and the term $J(\Theta)$ expresses the configuration of the articulated robot as a function of the joint angles. In short, the quasi-static compliance modelling phase may adopt any model, from the basic six degree of freedom (DOF) torsional spring [26] to the elaborate 258 DOF link and joint compliance model [10], but uses measurement data from a quasi-static measurement instead of a static one to identify the model parameters.

# 3 Methodology

## 3.1 Measurement procedure — description

The proposed measurement procedure quantified the QS deflections. Those are the EE deflections while moving along a trajectory. The trajectories under investigation were circular due to the design constraints of the utilised measurement instrument, the loaded double ball bar (LDBB) [27]. A more detailed description of the measurement procedure can be found in [21] and [28].

This work aimed at quantifying the influence of varying velocity of a QS measurement on a single industrial articulated robot (ABB IRB 6700). Furthermore, it applied the same QS measurement for a single Cartesian velocity on two more manipulators (KUKA KR 270 and a Stäubli TX 200) to evaluate the transferability of the findings.

The experiments comprised measurements with an ABB IRB 6700 at 0, 1, 5, 50, 150, and 500 mm s⁻¹. The label 0 mm s⁻¹ was meant to express a static measurement in accordance with the ISO standard. There were two additional sets of measurements performed on a KUKA KR 270 and a Stäubli TX 200 at 0 and 50 mm s⁻¹. The measurement of the deflections on the ABB IRB 6700 was recorded using both a laser tracker (LT) and the linear variable displacement transducer (LVDT) located inside the LDBB. The LT’s maximum sampling rate is 1000 Hz. Thus, the potential incorrect attribution of the movement as deflection may lead to significant errors at velocities from 150 onwards, e.g. close to 500 µm for LT measurements at the velocity of 500 mm s⁻¹. Therefore, no data is presented from the LT for the velocities of 150 and 500 mm s⁻¹.

In short, the measurement procedure provided data to identify the wrench $W_{QS}$ and the quasi-static deflections $\Delta X_{QS}$ at a manipulator measurement configuration $\Theta$. These data were synchronised to avoid a biased attribution of the former two to the latter. This is discussed in detail starting in the next Section 3.2.

The manipulator measurement configurations (MMC) were generated using offline programming. The circular trajectory for these experiments and the poses along the trajectory can be visualised through coordinate systems in Fig. 2. For the static measurement, these MMCs correspond to the discretised Cartesian poses along the circular trajectory. Configuration or target T01 and T37 were at the same position, meaning that the static circle path had the same start and end point.

For the quasi-static measurement, there were only four nominal MMCs. These were the vertices on two arcs that described a circle using the circular movement functions, namely, ABB’s RAPID MoveC, KUKA’s KRL CIRC, and Stäubli’s VAL3 movec. The whole trajectory was created to obtain several MMC in between the vertices, which equaled the number of measurement points along the trajectory. First, the manipulator executed a trigger movement of 3 mm in the Z-axis direction at the start point of the circle. Then, the manipulator moved two times clock-wise (CW) about the same circular trajectory (see Fig. 3).

The quasi-static wrench exerted multiple force components simultaneously. The Z-axis component’s magnitude equalled approx. 1/3 of the total load, while the remaining 2/3 were split between the X- and Y-axis components. Their contribution depended on the position of the LDBB along the circular trajectory [28]. For the static measurement, once

<figure><img src="image_3.png" /><figcaption>Fig. 2 The targets along the circle for the static (above) and QS measurement (below). The targets along the circle’s first quadrant are exemplified for the static measurement</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 3** LDBB quasi-static manipulator measurement; CW signal was the sum of the measurement data and the angular overshoot (AOS) movement</figcaption></figure>

the manipulator reached a stable position, the static loads were induced at the tool centre point (TCP).

The measurement application point (MAP) for the LVDT is located at the TCP of the EE. The spherically mounted retroreflector (SMR), target for the MAP of the LT, is located close to the TCP. The MAPs for the LVDT and the LT can be visualised in Fig. 4. The experiments on the KUKA KR 270 and Stäubli TX 200 were measured using only the LVDT. Additional parameters such as the position of the centre of the circular trajectory in the robot base coordinate system (RBCS), AOS, and apparent load (AL) for each of the three industrial manipulators can be seen in Table 1. Each industrial manipulator has different boundary constraints in terms of time, operating space availability, and control system. This led to varying implementations of the measurement procedure on each manipulator. Each measurement was repeated five times.

## 3.2 Measurement procedure — error separation analysis

This subsection evaluates all factors associated with the measurement procedure that could have introduced bias to the measurement data.

### 3.2.1 Transient measurement data

In the context of this work, transient measurement data can be defined as data that was acquired from the industrial manipulator during the change from equilibrium or steady state to a moving state. There were mainly two different kinds of transient phases in the measurement stages associated with loaded circular testing, which corresponded to acceleration and deceleration as well as trigger movements. For instance, for the measurement of

<figure><img src="image_5.png" /><figcaption>**Fig. 4** Visualisations of the three industrial manipulators with their EEs and MAPs. Robot (a) has a dummy tool which has one MAP for the LT close to the TCP, and the MAP for the LVDT is in the TCP. Robots (b) and (c) are equipped with motor spindles for machining</figcaption></figure>
<table><caption>**Table 1** Experimental conditions for the comparison of static and quasi-static deflections of three different industrial manipulators</caption><thead><tr><th>Manipulator</th><th>ABB IRB 6700</th><th>KUKA KR 270</th><th>Stäubli TX200</th></tr></thead><tbody><tr><td>Centre point in RBCS in mm</td><td>[625 - 950 1300]</td><td>[1250 450 550]</td><td>[750 1150 - 200]</td></tr><tr><td>Circle diameter in mm</td><td>250 mm</td><td>250 mm</td><td>250 mm</td></tr><tr><td>AOS in °</td><td>180°</td><td>180°</td><td>180°</td></tr><tr><td>AL</td><td>125:125:500 N</td><td>125:125:500 N</td><td>125:125:500 N</td></tr><tr><td>MAP</td><td>TCP LVDT [100 0 180]<br>SMR LT [145 0 165]</td><td>TCP LVDT<br>[155 0 245]</td><td>TCP LVDT<br>[160 0 275]</td></tr><tr><td>No. static configurations</td><td>37</td><td>21</td><td>37</td></tr><tr><td>Cartesian velocity in mm s<sup>-1</sup></td><td>1, 5, 50, 150, 500</td><td>50</td><td>50</td></tr><tr><td>No. QS configurations</td><td>784250, 157900,<br>15 650, 5 350, 1 550</td><td>19276</td><td>15874</td></tr></tbody></table>

quasi-static stiffness, the industrial manipulator is moving along a trajectory at a constant velocity [21]. At the beginning of the test, the manipulator had to accelerate until it reached constant velocity and, at its end, decelerate back to zero velocity. The measurement data captured during these stages contained contributions from the manipulator's dynamics. This would have introduced bias in the identification of the static mechanical stiffness, as the model does not look to quantify the effect of the inertia and velocity terms. Hence, an overshoot travel, also known as AOS, is employed. The magnitude of the AOS selected was 180° as the transient behaviour was assumed to be contained within this domain.

Furthermore, the controller of the industrial manipulators under investigation and the measurement instrument did not communicate bilaterally. Even if they could communicate, there could have existed a significant delay [29]. Thus, trigger signals were employed to facilitate data synchronization. Trigger signals are usually mechanical movements that involve the fast movement of the TCP by a significant distance compared to the measurand from the referenced starting position and orthogonal to the commanded trajectory. The trigger movement is usually employed at the measurement's beginning and end to highlight the measurement data. Both the AOS and the trigger movement can be visualised in Fig. 3.

### 3.2.2 Mechanical base load reference

All measurements presented in this analysis featured a mechanical base load reference (MBLR). The deflections of the industrial manipulators were not derived as the difference between the unloaded and loaded configuration but as the difference of a loaded configuration with respect to another loaded configuration at the MBLR. This approach can also be considered as a pre-loading of the components

in the force loop [25]. This procedure was meant to reduce errors in the measurement resulting from play or backlash (out of the scope of this work). Moreover, some industrial manipulators may exhibit hysteresis that could be reflected in the measurements if the comparison included both CW and counter clock-wise (CCW) trajectories. Therefore, this comparison was excluded from this investigation.

In general, the measurements were intended to capture the compliance of the industrial manipulator without contributions from other error sources. Thus, a suitable MBLR needed to be selected. By default, the MBLR in all tests equalled 125 N. The MBLR was adjusted depending on the payload and reach of the industrial manipulator if additional information was available. The wording AL is used to highlight the idea that a MBLR has been used. An AL of 250 N indicated that the system is loaded with 375 N, but the deflection was estimated with respect to the MBLR of 125 N. This is visualised in Fig. 5. The figure is a schematic reproduction of the measurement data, which means that the data is qualitatively accurately represented, but not quantitatively. The magnitudes were changed to highlight the observations. The figure shows the two reference trajectories at 50 and 500 mm s⁻¹ as well as the trajectories for the AL of 250 and 500 N. It can be assumed that at lower velocities, the kinematic errors predominantly affect the position of the reference trajectory. In comparison, dynamic and controller errors should be more dominant at higher velocities. The concept of the MBLR should ensure that these error sources do not bias the compliance evaluation, as one would expect the kinematic and control errors for the trajectory at AL of 125 and 50 mm s⁻¹ and AL of 250 and 50 mm s⁻¹ to have the same magnitude. Therefore, their difference can be considered as an accurate reflection of the compliance of the industrial manipulator. In addition, it could be seen in Fig. 5 that the reference trajectories were offset from one another based on
<figure><img src="image_6.png" /><figcaption>**Fig. 5** Schematic representation of the measurement data at a velocity of 50 and 500 mm s⁻¹ (solid and dashed line, respectively) and ALs of 250 and 500 N. The circle centers of the 50 and 500 mm s⁻¹ circular trajectories appear as offset due to the change in velocity (a phenomenon which was observed from the measurements)</figcaption></figure>

their velocity. This shift of the reference trajectory did not introduce an error because the deflections of the industrial manipulators are calculated as the relative difference to each corresponding (velocity-wise) reference trajectory.

### 3.2.3 Measurement instruments

The measurement instruments require the ability to exert a controllable mechanical load while measuring the deflection in terms of distance and position. This subsection evaluates the influences of the measurement procedure on the deflection measurement and the mechanical load exertion.

This article presented measurement data, which were obtained from a LVDT as well as a LT. The measurements were performed using a Micro-Epsilon DTA-3G8-3-CA® [30] LVDT for all velocities and with a Leica AT901® [31] LT up to a velocity of 50 mm s⁻¹. The sampling rate of the measurement instruments was selected to ensure that at least one data point was measured at a nominal travel distance of 50 µm, i.e. at a Cartesian velocity of 50 mm s⁻¹ the sampling rate equals 1000 Hz. For both measurement instruments, there was no technical specification about their accuracy for dynamic measurements such as the QS measurements used in this work. The uncertainty associated with the position measurement of the LT and the distance measurement of the LVDT are stated in Table 2.

The LDBB induces load by controlling the pressure in the chamber using a proportional pressure control valve (PPCV) [32]. The PPCV and its measurement systems are not affected by this static load, as the load is not exerted

on them. However, the load fluctuated for the quasi-static measurement at higher velocities. For a trajectory at 50 mm s⁻¹, there existed a steady-state control error but no fluctuations. The load fluctuated within approximately 2% at 500 mm s⁻¹. The force along the static and quasi-static trajectory was also measured on the side of the link (TL), by mounting the TL to a Kistler 9255C Dynamometer® [33]. The test measured the loads on the TL for both the static and quasi-static measurements. The results differed by less than 1 N on average for all ALs. This experiment was repeated five times. The uncertainty associated with the control of the magnitude of the static force to push between the two objects is stated in Table 2.

## 4 Results and discussion

The effect of the velocity of the movement on the difference between static and QS compliance of industrial manipulators is shown in Figs. 6 and 7 and Table 3. The figures illustrate the average Cartesian deflections and uncertainty for the static and QS procedures for the ABB IRB 6700. Figure 7 shows that the static and QS deflections are significantly different in magnitude but similar in trend. The same phenomenon was observed in two other manipulators included in the study. Table 3 summarises the average deflections and their uncertainty for the static and the QS measurement data for the three industrial manipulators. The static measurement data were considered the reference, see the column *Static* 0 mm s⁻¹. It can be

<table><caption>**Table 2** Investigated type B uncertainty contributors associated with the proposed measurement procedure. It is assumed that these uncertainty contributors are normally distributed</caption><thead><tr><td>Contributor</td><td>Magnitude</td><td>Assmt. type</td></tr></thead><tbody><tr><td>LDBB load variation</td><td>Max. 2.8 N</td><td>B</td></tr><tr><td>LVDT dist. accuracy</td><td>5 µm full-scale</td><td>B</td></tr><tr><td>LT posi. accuracy</td><td>10 µm + 5 µm m<sup>-1</sup></td><td>B</td></tr></tbody></table>

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Comparison of the average static deflection and QS deflections measured at different speeds for the IRB 6700. The static measurement data were considered as the reference (100%)</figcaption></figure>

seen that for each QS measurement, the deflections are greater than for the static deflections. At 50 mm s⁻¹, the QS deflections are 8 to 17% bigger than the static deflections.

For the ABB IRB 6700, see Fig. 6, the results suggest that the Cartesian velocity has a proportional relationship to the measured deflections. The LVDT data imply that beyond a velocity of 50 mm s⁻¹, the higher the velocity, the higher the QS deflections. Additionally, both the LVDT and LT data suggest that for velocities between 1 and 50 mm s⁻¹, a constant change in magnitude can be expected.

Additionally, as shown in Fig. 5, the positions of the circular trajectory varied at different velocities. The smaller the levers created by the configuration, the higher the robot's Cartesian stiffness [10]. Other authors have attributed the difference in the deflection of the robot to link weights and to the effect of the gravity compensator [34]. Therefore, a change is expected in the Cartesian stiffness and measured deflections as the manipulator is not always measured at the same position. The same notion of position shift also applies to the repetition of the circular trajectory at each velocity as defined by the path positioning repeatability of the industrial manipulator [6]. Under the assumption that it is only the change in configuration θ that influences the Cartesian stiffness Kₓ, the change in Cartesian stiffness dependent on the change of the configuration can be quantified according to Eq. 1. For the investigation of the ABB IRB 6700, the change in configuration from the QS trajectories of 1 to 500 mm s⁻¹ equalled on average a positioning change of 1.5 mm. Equation 1 was used to simulate the effect of the change in position on the expected Cartesian deflection with the following set of parameters: normally distributed translational offsets with an average of ±2 mm and normally distributed orientation offsets with an average of ± 20 m deg in the Cartesian components at

100 000 samples around the circular trajectory. According to the simulation results, this systematic change in the configuration should have contributed with less than 1% to the difference in deflection between the static and QS trajectories.

Theoretically, as the excitation frequencies in the QS wrench or the manipulator's movement are well below the first eigenfrequency, one would expect no significant differences between static and QS deflections. The frequency components of the quasi-static measurements at the different velocities have been identified using a fast Fourier transform (FFT) of the LVDT and LT measurement data. The maxima of the spectral magnitude occur at frequencies less than 1 Hz. Nevertheless, the performed measurements have shown the opposite. Thus, it could be argued that the combined effect of manipulators mechanics in terms of the inertial forces, and backlash as well as the control in terms of the control loop gains, and the criterion of passage of the position control might explain the difference between the static and QS deflections. This difference can be perceived as a relative loss of stiffness between static and moving industrial articulated robots. An error separation analysis of these effects can support the understanding of the systematic difference between static and QS deflections. This work presented how the QS measurement procedure might have influenced the measurand, and it was not possible to find an associated systematic bias.

## 4.1 Uncertainty

The investigated type B uncertainty contributors, presented in Table 2, can be considered insignificant in comparison to the measurands. This is expressed through the combined expanded (k = 2) uncertainty presented in Table 2. The
<figure><img src="image_9.png" /></figure>

(a) Mean Cartesian deflection and uncertainty in the direction of the X-axis of the RBCS.

<figure><img src="image_10.png" /></figure>

(b) Mean Cartesian deflection and uncertainty in the direction of the Y-axis of the RBCS.

<figure><img src="image_11.png" /></figure>

(c) Mean Cartesian deflection and uncertainty in the direction of the Z-axis of the RBCS.

<caption>**Fig. 7** Mean Cartesian deflections of the ABB IRB 6700 and their uncertainty in the RBCS</caption>

<table><caption>**Table 3** Comparison of average static and QS deflections of three different industrial manipulators</caption><thead><tr><th></th><th></th><th>Static 0 mm s<sup>-1</sup></th><th>QS 50 mm s<sup>-1</sup></th></tr></thead><tbody><tr><td rowspan="2">ABB IRB 6700</td><td>LT def. in %</td><td>100%</td><td>114 ± 5%</td></tr><tr><td>LVDT def. in %</td><td>100%</td><td>108 ± 5%</td></tr><tr><td>KUKA KR 270</td><td>LVDT def. in %</td><td>100%</td><td>108 ± 1%</td></tr><tr><td>Stäubli TX200</td><td>LVDT def. in %</td><td>100%</td><td>117 ± 2</td></tr></tbody></table>

measurand is assumed to be normally distributed. This information has been used to create the visualisation of the error bars in Fig. 7 and to calculate the uncertainty presented in Table 2.

## 5 Conclusion

The work presented in this article quantifies the influence of Cartesian velocity on quasi-static deflections of industrial manipulators and analyses how the QS measurement procedure could have influenced the presented measurands. For all industrial articulated robots in this experimental investigation, the static and quasi-static deflections were significantly different in magnitude but similar in trend, being the quasi-static deflections larger. No systematic error in the measurement procedure could be found to explain this difference.

From the authors' points of view, this implies that the manipulator's calibrated compliant model obtained according to ISO 9283 may have systematic inaccuracies associated with the robot motion. Therefore, a calibrated compliant manipulator model derived from a QS measurement may yield more accurate estimates for trajectory optimisation compared with conventional static stiffness identification procedures. A successful implementation of quasi-static compliance calibration requires a measurement procedure that produces accurate results with fewer resource requirements. For this scheme to be economically viable, it is preferential to perform the calibration faster, i.e. reducing the required downtime of the industrial manipulator. In a quasi-static compliance calibration, the downtime consists of the setup and measurement time. The measurement time is proportional to the velocity of the motion. Therefore, the effect of the Cartesian velocity on the measured deflections was investigated. According to the results obtained, the magnitude of the quasi-static deflections had a proportional relationship to the Cartesian velocity above above 50 mm s¹. Based on the observed influence of the velocity on the load-induced deflections, the suggestion is to perform a quasi-static compliance calibration that considers not only the combination of force and motion, but pays special
attention to the Cartesian velocity of the trajectory for the
intended application.

**Author contribution** Monica Katherine Gonzalez: methodology, investigation, formal analysis, writing — original draft. Nikolas Alexander Theissen: methodology, investigation, formal analysis, writing — original draft. Nora Agirre: resources, investigation. Patxi Hacala: resources, investigation. Jon Larrañaga: resources, investigation, writing — review and editing. Andreas Archenti: supervision, writing — review and editing.

**Funding** Open access funding provided by Royal Institute of Technology. The authors would like to thank VINNOVA (Sweden's Innovation Agency), CDTI (Spanish Centre for the Development of Industrial Technology), and the SMART Advanced Manufacturing Cluster for funding this study as a part of the COMACH project (Grant Agreement ID: S0120-COMACH). This study has also been conducted partially under the framework of the project MIRAGED (CER-20191001), supported by CDTI-Acreditación y Concesión de Ayudas Destinadas a Centros Tecnológicos de Excelencia Cervera, with partial funding from the European Union's Horizon 2020 research and innovation programme (Grant agreement number: 971442 FIBREMACH project), ARISTARCO project (Government of Navarra- R&D 2020 call), and supported by the European Regional Development Fund (Fondos Feder). Finally, the authors would also like to express their gratitude to the Center for Design and Management of Manufacturing Systems and Excellence in Production Research for their financial support.

Declarations

**Conflict of interest** The authors declare no competing interests.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

References

1. International Federation of Robotics. History

2. Tunc LT, Gonul B (2021) Effect of quasi-static motion on the dynamics and stability of robotic milling. CIRP Annals

3. Bo L, Tian W, Zhang C, Hua F, Cui G, Li Y (2021) Positioning error compensation of an industrial robot using neural networks and experimental study. Chinese Journal of Aeronautics

4. Cvitanic T, Nguyen V, Melkote SN (2020) Pose optimization in robotic machining using static and dynamic stiffness models. Robot Comput Integr Manuf 66:101992

5. Nubiola A, Bonev IA (2013) Absolute calibration of an ABB IRB 1600 robot using a laser tracker. Robot Comput Integr Manuf 29(1):236–245

6. International Organization for Standardization. ISO 9283:1998 Manipulating industrial robots - performance criteria and related test methods

7. Mooring B, Roth ZS, Driels MR (1991) Fundamentals of manipulator calibration. Wiley, New York

8. Aoyagi S, Suzuki M, Takahashi T, Fujioka J, Kamiya Y (2012) Calibration of kinematic parameters of robot arm using laser tracking system: compensation for non-geometric errors by neural networks and selection of optimal measuring points by genetic algorithm. Int J Autom Technol 6(1):29–37

9. Jawale HP, Thorat HT (2013) Positional error estimation in serial link manipulator under joint clearances and backlash. J Mech Robot 5:2

10. Klimchik Alexandr, Furet Benoit, Caro S, Pashkevich A (2015) Identification of the manipulator stiffness model parameters in industrial environment. Mech Mach Theory 90:1–22

11. Theissen NA, Mohammed A, Archenti A (2019) Articulated industrial robots: an approach to thermal compensation based on joint power consumption. In: Blunt L, Knapp W (eds) Laser metrology and machine performance XIII, pages 81–90. European Society for Precision Engineering and Nanotechnology, Bedfordshire, UK

12. Ibaraki S, Theissen NA, Archenti A, Alam M (2021) Moktadir evaluation of kinematic and compliance calibration of serial articulated industrial manipulators. Int J Autom Technol 15(5):567–580

13. Dumas C, Caro S, Garnier S, Furet B (2011) Joint stiffness identification of six-revolute industrial serial robots. Robot Comput Integr Manuf 27(4):881–888

14. Zaeh M. F., Roesch O. (2014) Improvement of the machining accuracy of milling robots. Prod Eng 8(6):737–744

15. Siciliano B, Khatib O (2007) Handbook of robotics. Springer-Verlag, Berlin

16. Verl A, Valente A, Melkote S, Brecher C, Ozturk E, Tunc LT (2019) Robots in machining. CIRP Ann Manuf Technol 68(2):799–822

17. Ji W, Wang L (2019) Industrial robotic machining: a review. Int J Adv Manuf Technol 103(1):1239–1255

18. Theissen NA, Gonzalez MK, Archenti A (2022) Working vs. operating space kinematic calibration of serial articulated industrial manipulators. In: EUSPEN GA 22t International Conference and Exhibition 30th May – 3rd June 2022

19. Garnier S, Subrin K (2022) A metrological device for robot identification. Robot Comput Integr Manuf 73:102249

20. Schneider U, Drust M, Ansaloni M, Lehmann C, Pelliccari M, Leali F, Gunnink JW, Verl A (2016) Improving robotic machining accuracy through experimental error investigation and modular compensation. Int J Adv Manuf Technol 85(1):3–15

21. Theissen NA, Gonzalez MK, Barrios A, Archenti A (2021) Quasi-static compliance calibration of serial articulated industrial manipulators. Int J Autom Technol 15(5):590–598

22. Rivin EI (2010) Handbook on stiffness & damping in mechanical design. ASME Press, New York

23. Weck M, Brecher C (2006) Statisches Verhalten von Werkzeugmaschinen. Springer, Berlin, pp 163–177

24. International Organization for Standardization. ISO 230-1:2012 Test code for machine tools — Part 1: Geometric accuracy of machines operating under no-load or quasi-static conditions

25. Laspas T, Theissen N, Archenti A (2020) Novel methodology for the measurement and identification for quasi-static stiffness of five-axis machine tools. Precision Eng 65:164–170

26. Kenneth Salisbury J (1980) Active stiffness control of a manipulator in cartesian coordinates. 95–100
27. Archenti A, Nicolescu M (2013) Accuracy analysis of machine tools using elastically linked systems. CIRP Ann Manuf Technol 62(1):503–506

28. Gonzalez M, Hosseini A, Theissen NA, Archenti A (2020) Quasi-static loaded circular testing of serial articulated industrial manipulators. In: Verl A, Parisel N (eds) 52nd International Symposium on Robotics // ISR 2020, Berlin. VDE VERLAG

29. Zhou R, Kauschinger B, Ihlenfeldt S (2021) Data synchronization by continuous spatial measurement with double ballbar. Measurement 174:108909

30. Micro-Epsilon (2017) induSENSOR-LVDT

31. Leica. User Manual AbsoluteTracker AT901

32. Theissen NA (2021) Precision measurement instruments for machinery’s mechanical compliance: design and operation : Mea-

surement instruments for physics-based calibration of advanced manufacturing machinery. PhD thesis, KTH, Production Engineering

33. Kistler Instrumente AG (2021) Multi-component dynamometer up to 60 kn: quartz 3-component dynamometer type 9255c

34. Peng X, Yao X, Liu S, Wang H, Liu K, Senthil Kumar A, Wen Feng L, Bi G (2021) Stiffness modeling of an industrial robot with a gravity compensator considering link weights. Mech Mach Theory 161:104331

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>