Journal cover Journal topic
Journal of Sensors and Sensor Systems An open-access peer-reviewed journal
Journal topic
J. Sens. Sens. Syst., 7, 577-586, 2018
https://doi.org/10.5194/jsss-7-577-2018
J. Sens. Sens. Syst., 7, 577-586, 2018
https://doi.org/10.5194/jsss-7-577-2018

Regular research article 07 Nov 2018

Regular research article | 07 Nov 2018

# Dynamic characterization of multi-component sensors for force and moment

Dynamic characterization of multi-component force-moment sensors
Jan Nitsche1,2, Rolf Kumme1, and Rainer Tutsch2 Jan Nitsche et al.
• 1Physikalisch-Technische Bundesanstalt, Bundesallee 100, 38116 Braunschweig, Germany
• 2Technische Universität Braunschweig, Institut für Produktionsmesstechnik, Schleinitzstr. 20, 38106 Braunschweig, Germany
Abstract

An improved set-up for the characterization of multi-component sensors for force and moment is presented. It aims at calibrating such sensors under continuous sinusoidal excitation. Special focus is put on the design of load masses and adapting elements to activate uniaxial force and moment components where possible. To identify the motion and acceleration of the load mass with 6 degrees of freedom, a photogrammetric measurement system is implemented in the existing set-up. Using the set-up described, different experiments are performed to analyse a commercial multi-component sensor and perform a parameter identification for its force components.

1 Introduction

The calibration and characterization of sensors for force and moment are typically performed in a static manner . The dynamic behaviour of these sensors may differ significantly from the static behaviour. Especially the sensitivity coefficient, which is necessary to transform the sensor output signal into a force or moment signal, is subject to changes with higher frequencies of the dynamic excitation . Therefore, if such sensors are to be used in a dynamic environment, their dynamic behaviour needs to be analysed . For force and acceleration sensors, typical characterization methods are sinusoidal excitation and shock excitation . Torque sensors are currently evaluated under sinusoidal excitation (Klaus2016).

The aforementioned calibration approaches are limited to uniaxial sensors or calibrations of only one axis of multiaxial or multi-component sensors (MCSs). Such MCSs have become more popular in the last few years, which has resulted in the need for new calibration procedures . Challenges in the calibration of MCSs are, among other things, the generation of forces and moments with defined directions, the identification of the force and moment vector in a reference system and the alignment of sensor and reference coordinate systems .

For the dynamic calibration of MCSs, very few approaches exist. deployed an electrodynamic shaker set-up with a vertical excitation direction at the Physikalisch-Technische Bundesanstalt (PTB) using different adapting elements to activate sinusoidal excitations in different directions. designed a fixed sensor frame with different adapting directions for an electrodynamic plunger actuator. In both set-ups, the reference force can only be calculated in the direction of excitation. Additional components like rocking modes cannot be detected.

In the following sections, an improved set-up for the dynamic analysis of MCSs is described. The force generation is based on the periodic acceleration of a load mass connected to an MCS. The theory behind sinusoidal calibration is explained using the example of uniaxial force calibration in Sect. 2. The improved set-up includes a specially designed load mass and a three-dimensional acceleration reference based on photogrammetry. Section 3 summarizes the new set-up. A dynamic model of an MCS is described in Sect. 4. Finally, the parameters of the model are identified from experimental data as described in Sect. 5.

2 Sinusoidal calibration of force and moment sensors

The aim of the dynamic calibration of force and moment sensors is the identification of the parameters needed to describe the dynamic behaviour and the sensitivity coefficient as a function of the excitation frequency . The typical model used for the description of such a sensor is a mass-spring-damper system. The parameters needed to describe this model are the spring constant kz, the damping parameter bz and the mass m (Fig. 1). The mass m is composed of the external load mass me, the mass of adapting elements ma and the internal moving mass of the sensor m1.

Figure 1Dual-mass oscillator as a model for uniaxial force sensors.

From the given parameters, the dynamic sensitivity S can be calculated according to :

$\begin{array}{}\text{(1)}& S=\frac{c}{\sqrt{{k}_{z}^{\mathrm{2}}+\left({b}_{z}\mathit{\omega }{\right)}^{\mathrm{2}}}}\approx \frac{c}{{k}_{z}}\left(\mathrm{1}-\frac{\mathrm{1}}{\mathrm{2}}{\left(\frac{{b}_{z}}{{k}_{z}}\right)}^{\mathrm{2}}{\mathit{\omega }}^{\mathrm{2}}\right),\end{array}$

with c as a factor for unit conversion with the unit U m−1 (U as a sensor signal unit, e.g. mV V−1 for strain gauge sensors). In addition to the sensitivity S, the phase shift φs between the excitation and the sensor output can be calculated as follows:

$\begin{array}{}\text{(2)}& {\mathit{\phi }}_{\mathrm{s}}\approx -\frac{{b}_{z}}{{k}_{z}}\mathit{\omega }.\end{array}$

Both Eqs. (1) and (2) are approximations for a stiff coupling of the top mass to the sensor, which allows the influence of an additional spring-damper system describing the coupling to be neglected.

As an alternative to identifying the model parameters k and b, the dynamic sensitivity can be calculated as the ratio of the sensor signal output U and the acting dynamic force F :

$\begin{array}{}\text{(3)}& S=\frac{U}{F}.\end{array}$

In this case, only the moving mass m and its acceleration a are needed to calculate the acting force:

$\begin{array}{}\text{(4)}& F=m\cdot a.\end{array}$
3 Measurement set-up

The measurement set-up for the periodic excitation of MCSs used in this work is based on the shaker set-up located in the dynamic force calibration laboratory at the PTB . It consists of an LDS 10 kN electrodynamic shaker for frequencies up to 2 kHz and a laser scanning vibrometer as an acceleration reference. Different accelerometers can be used as additional acceleration references.

The electrodynamic shaker is only capable of generating accelerations in one direction. To activate force and moment components in all six directions, different adapting elements are needed. The design of such elements is described in Sect. 3.1.

The scanning vibrometer measures accelerations in the direction of the excitation of the shaker system. Accelerations perpendicular to the excitation direction, as well as twisting motions of the load mass, cannot be detected. To rectify this, a photogrammetric measurement system is used to detect accelerations that are not visible to the scanning vibrometer. Details of the photogrammetric system, synchronization and data analysis are described in Sect. 3.2.

In past works , different MCSs were investigated on the previously described shaker in a similar set-up. This set-up consisted of a cylindrical load mass in combination with an air-bearing guide for axial force. Transverse force components were generated using 90 angular adapters and a cubic load mass. Bending moments and torque were generated using a lever arm in combination with the aforementioned cubic load mass.

Due to its design, this set-up was only capable of generating uniaxial force components in the axial force direction Fz. Transverse force components Fx and Fy were superimposed by moment components My and Mx as a result of the distance between the sensor coordinate system and the centre of gravity of the external load mass. Torque components Mz were superimposed by bending moments Mx or My and transverse forces Fy or Fx.

To overcome these disadvantages, a special load mass was designed in order to move the centre of gravity of the load mass to the origin of the coordinate system of the sensor under test . Using this load mass, rocking modes for the axial force component Fz are reduced and transverse force components Fx and Fy can be generated without superimposed moment components. The mass of the load mass is ${m}_{\mathrm{L}}=\left(\mathrm{6.9605}±\mathrm{0.0010}\right)$kg. The angular adapter for activating transverse forces is designed in order to move the centre of gravity of the whole set-up in the axis of the shaker table to avoid moment loads on the shaker bearings. The set-up including load mass, sensor and the adapting element for transverse force components is shown in Fig. 2.

For the generation of moment components, a beam cross is used. The cross is designed in a way to move the centre of gravity of the moment generating mass components to the xy plane of the sensor coordinate system. This results in a superposition of bending moments Mx or My with only one axial force Fz and the superposition of torque Mz with only one transverse force Fx or Fy. The distance between the load mass and the sensor can be changed to generate different amplitudes of the moment component. To reduce moment loads on the shaker bearing, a counter mass is attached to the base of the angular adapter. The beam cross set-up is shown in Fig. 3.

Figure 3Beam cross for generation of moment components.

## 3.2 Three-dimensional acceleration reference

To evaluate whether signal outputs of inactive force or moment components of the sensor are a result of signal crosstalk from the active component or of bending or twisting moments, a complete three-dimensional movement of the set-up needs to be identified. The scanning vibrometer is capable of measuring the acceleration in the direction of the shaker axis at different points of the set-up. From the acceleration distribution, rocking modes around the x and y axes can be identified. Three-dimensional accelerometers can be added to specific points in the system to obtain three-dimensional acceleration information at that point. However, the mass of the accelerometer and its cable has an influence on the dynamic behaviour of the whole set-up.

To extend the acceleration measurement provided by the laser vibrometer and the accelerometers, a photogrammetric measurement set-up is installed. It consists of two stereo camera systems: one camera system observing the load mass on top of the sensor, the other one directed at the surface of the shaker. All four cameras are triggered simultaneously using a signal generator. One master camera is used to trigger two LED strobe lights to illuminate the observed area. This trigger signal is recorded by the analogue to digital converter (ADC) in the junction box of the data acquisition PC of the scanning vibrometer for synchronization purposes. The trigger signal also serves as the time stamp of the image acquisition. The six force and moment signals of the sensor under test are amplified by a 6-channel bridge amplifier (Dewetron2012), digitized by a simultaneously sampling 16-channel ADC card and saved on a separate data acquisition PC. The camera trigger signal is also recorded by this PC for synchronization. A block diagram of the extended shaker set-up is shown in Fig. 4. A photograph of the sensor with the attached load mass, the shaker and the stereo camera systems is shown in Fig. 5.

Figure 4Block diagram of the photogrammetric set-up in the dynamic shaker system.

Figure 5Photogrammetric set-up in the shaker environment.

Figure 6Vibrometer and camera trigger signals.

Random greyscale patterns are attached to the surfaces observed by the cameras. The images of these patterns are analysed using digital image correlation (DIC) algorithms . In DIC, the camera images are split into small areas and the correlation between the areas in different camera images is calculated. In this way, one area of the observed surface can be identified in the images of both cameras and at different times. Each stereo set-up is able to calculate the position of multiple areas on the observed surface in space. Those three-dimensional points can be used to calculate a three-dimensional rigid body movement as well as a deformation of the load mass in the observed area.

The rotation and translation of the rigid body transformation are calculated using the singular value decomposition (SVD) method presented by . From the resulting rotation matrix, the rotation angles are calculated in a fixed axis system. This assumption can be used as the rotations around any axis are expected to be very small and therefore do not influence each other.

Table 1Fitted sine parameters.

The cameras used for DIC are industrial CMOS cameras with frame rates of up to 166 images per second for 8-bit greyscale images. At 12-bit colour depth, the frame rate is reduced to a maximum of 82 images per second. This frame rate is too low to satisfy the Nyquist–Shannon sampling theorem for frequencies higher than 40 Hz. However, using the additional knowledge of the excitation frequency of the shaker and the precise adjustment of the sampling rate of the cameras, both frequencies can be adjusted to achieve a beat frequency.

The known frequency distribution resulting from the vibrometer measurement is used to define a fitting function for the displacement of the load mass S(t):

$\begin{array}{}\text{(5)}& S\left(t\right)=\sum _{i}\left({\stackrel{\mathrm{^}}{s}}_{i}\cdot \mathrm{sin}\left(\mathrm{2}\cdot \mathit{\pi }\cdot {f}_{\mathrm{i}}\cdot t+{\mathit{\phi }}_{i}\right)\right)+o,\end{array}$

with fi as the identified frequencies of the frequency distribution, ${\stackrel{\mathrm{^}}{s}}_{i}$ the amplitude, φi the phase shift of each frequency and o the offset of the mean. The parameters $\stackrel{\mathrm{^}}{s}$, φ and o can be estimated by a least squares fitting algorithm. The resulting function S(t) represents the axial or angular displacement of the load mass over time.

## 3.3 Experimental analysis of the set-up

An experimental evaluation of the set-up described is performed using the axial force components. The excitation frequency is set to fex=106.25Hz. The amplitude of the acceleration signal is 0.5 V. The laser vibrometer sampling frequency is 51.2 kHz and the cameras are triggered with 52.874 Hz. Figure 6 shows the acceleration signal recorded from the laser vibrometer, a fitted sine function and the camera trigger signal as a synchronization reference. The parameters of the sine fitting are listed in Table 1. The camera trigger is started manually after starting the vibrometer recording. The trigger signal at 0.97 s represents the first image acquired by the cameras.

The difference between the camera frequency and the excitation frequency results in a beat of 0.502 Hz. The first camera image is taken as a reference image against which all displacements are calculated. Each camera acquires 211 images, which results in 210 displacement measurements within 3.97 s.

For each displacement measurement, the rigid body transformation is calculated as described in Sect. 3.2. As a result, three axial translations and three rotation angles along the coordinate system of the camera set-up are obtained. For each displacement measurement, the mean of the exposure time, which is identified from the strobe trigger signal, is used as the time of acquisition. The data sets of acquisition time and axial and angular displacement are used to fit a sine function according to the excitation frequency of the shaker.

Sine fitting is performed using a nonlinear least squares fitting function, curve fit, of the Python scipy optimize package. Starting values for the fitting frequencies are taken from the frequency spectrum of the vibrometer reference. Figure 7 shows the frequencies of the fast Fourier transform (FFT) up to 1000 Hz.

Figure 7Frequency spectrum of the vibrometer signal. The y axis is scaled to make frequencies different than the excitation frequency visible. The amplitude of fex=106.25Hz is 189.3 m s−2.

Figure 8Displacement in z from DIC and the fitted sine function.

The main contributions can be seen at the excitation frequency fex=106.25Hz and its multiples 425, 212.5 and 318.75 Hz. The amplitude at the second frequency (425 Hz) reaches a value of 0.5 % of the amplitude of the excitation frequency.

Sine fitting of only the excitation frequency fex=106.25Hz results in a standard deviation of the difference between the fitting function and the measured displacements of 0.0036 mm. Using a fitting function of the two main frequencies 106.25 and 425 Hz reduces the standard deviation of the residuum to 0.0023 mm. Additional fitting frequencies show no further improvement of the results. The displacement in z, measured with the photogrammetry set-up and the fitted sine function, is shown in Fig. 8. Figure 9 shows the axial and angular displacement along the x axis. The resulting parameters from the sine fitting are listed in Table 1.

Figure 9Axial and angular displacement in x.

To calculate the resulting force and moment components on the sensor, the axial and angular acceleration of the load mass is needed. Both values can be calculated from the fitting function of the displacement as the second derivative with respect to time. For a sine function, the second derivative can easily be calculated by multiplying the function by the negative square of the angular frequency ${\mathit{\omega }}_{i}=\mathrm{2}\cdot \mathit{\pi }\cdot {f}_{\mathrm{i}}$. The acceleration of the fitting function of Eq. (5) results in

$\begin{array}{}\text{(6)}& \stackrel{\mathrm{¨}}{S}\left(t\right)=\sum _{i}\left(-{\stackrel{\mathrm{^}}{s}}_{i}\cdot {\mathit{\omega }}_{i}^{\mathrm{2}}\cdot \mathrm{sin}\left({\mathit{\omega }}_{i}\cdot t+{\mathit{\phi }}_{i}\right)\right).\end{array}$

From the values of the fitting function listed in Table 1, the amplitude of the acceleration ${\stackrel{\mathrm{^}}{a}}_{i}=-{\mathit{\omega }}^{\mathrm{2}}\cdot {\stackrel{\mathrm{^}}{s}}_{i}$ from the photogrammetric measurement is calculated to be ${\stackrel{\mathrm{^}}{a}}_{\mathrm{1}}=-\mathrm{190.6131}$m s−2 and ${\stackrel{\mathrm{^}}{a}}_{\mathrm{2}}=\mathrm{28.5161}$m s−2. While the acceleration at the excitation frequency shows a small deviation of 0.68% to the reference value of the vibrometer, the acceleration at 425 Hz differs significantly. According to the reference acceleration, the amplitude of the second frequency has a minor influence on the overall acceleration and will be neglected in the subsequent analysis. It must however be included in the pending uncertainty budget.

4 Dynamic model of multi-component sensors

For uniaxial force sensors, a dual-mass oscillator as shown in Fig. 1 can be used as an acceptable description of the dynamic behaviour . MCSs show a more complex structure which results in a more complex physical model. To reduce the complexity of the model, the observed six-component MCS is split into two independent models, one for the three force components, one for the moment components.

## 4.1 Dynamic MC force model

The dynamic model for a three-component force sensor is based on a superposition of three orthogonally aligned spring-damper systems connected to two masses m1 and m2. A simplified illustration of a two-dimensional case of this model is shown in Fig. 10. Mutual influences between the spring-damper systems are neglected. To describe this model, eight different parameters are needed: three spring constants kx, ky, kz, three damping parameters bx, by, bz and the two masses m1 and m2. The masses are assumed to be equal for excitation in the x, y and z directions.

Figure 10Two-dimensional representation of the dual-mass oscillator for multi-component force sensors.

## 4.2 Dynamic MC moment model

A dynamic model for a three-component moment sensor is more complex than the model for a three-component force sensor. Basically, the model from Fig. 10 can be adapted by changing the springs and dampers of the model from linear to rotational ones. The exchange of the masses is more difficult, as the mass moment of inertia is needed for dynamic moment measurements. This parameter depends on a rotational axis and therefore it cannot be assumed to be equal for the different excitation directions.

5 Parameter identification

## 5.1 Determination of the internal mass

To determine the internal mass m1 of the sensor, the sensor is mounted onto an indexing head with a horizontal rotation axis. The z axis of the sensor is aligned with the rotation axis of the indexing head. Sensor readings are recorded without an attached load mass. The sensor is rotated in steps of 22.5. The sensor readings are transformed into force and moment values using the calibration matrix provided by the sensor manufacturer. A sine function is fitted to the force values for Fx and Fy. The amplitude of the sine function represents the maximum force on each channel resulting from the internal mass of the sensor. The mass itself is calculated from the amplitude, divided by the local acceleration of gravity gloc:

$\begin{array}{}\text{(7)}& {m}_{\mathrm{1}}=\frac{{F}_{\mathrm{max}}}{{g}_{\mathrm{loc}}}.\end{array}$

Force readings and the sine fitting for Fx and Fy are shown in Fig. 11. The resulting amplitude of the sine function is 7.558 N for Fx and 7.532 N for Fy, resulting in internal masses of m1x=0.770kg and m1y=0.768kg.

Figure 11Force readings and sine fitting for Fx and Fy.

The described experiment is repeated with the load mass shown in Fig. 2 attached to the sensor. Including the attached load mass, the amplitude of the sine function is 75.877 N for Fx and 75.736 N for Fy. The mass of the attached load mass is calculated according to Eq. (7) to mLx=6.962 and mLy=6.951kg, which differ by 0.028% and −0.139% respectively from the weighted mass mL=6.9605kg.

## 5.2 Determination of the spring constant and damping coefficient

The spring constant is calculated according to the method presented by from the transfer function H(f) as the ratio of the acceleration on the load mass at and the acceleration of the shaker surface ab:

$\begin{array}{}\text{(8)}& H\left(f\right)=\frac{{a}_{\mathrm{t}}}{{a}_{b}}.\end{array}$

The frequency of the maximum of the transfer function represents the resonance frequency f0, from which the spring constant can be calculated:

$\begin{array}{}\text{(9)}& {f}_{\mathrm{0}}=\frac{\mathrm{1}}{\mathrm{2}\mathit{\pi }}\sqrt{\frac{k}{m}}.\end{array}$

The transfer function H(f) for Fz, resulting from a periodic chirp excitation with two different load masses m1=4.236kg and m2=8.798kg, is shown in Fig. 12. Resonance frequencies are identified at ${f}_{\mathrm{0},\mathrm{1}}=\mathrm{1477.35}$Hz and ${f}_{\mathrm{0},\mathrm{2}}=\mathrm{1092.97}$Hz. With the additional internal mass mi=0.769kg, the spring constant according to Eq. (9) is calculated to be ${k}_{\mathrm{1}}=\mathrm{4.313}×{\mathrm{10}}^{\mathrm{8}}$N m−1 and ${k}_{\mathrm{2}}=\mathrm{4.512}×{\mathrm{10}}^{\mathrm{8}}$N m−1.

Figure 12Amplitude and phase of the transfer function H(f) for excitation in Fz.

Figure 13 shows the transfer function H(f) for excitation in Fx and Fy. Because of the horizontal installation of the sensor (Fig. 2), the transfer function is calculated as the ratio of the acceleration of the load mass and the acceleration of the mounting adapter connected to the sensor. The figure shows two resonance peaks for each load mass, which can be explained by the mounting set-up. The connection between the horizontal mounting adapter connected to the sensor and the angular adapter introduces a second spring-damper system into the model. It can be seen that the resonance frequencies are identical for excitation in Fx and Fy, while the amplitudes of the resonance peaks differ for the higher resonance frequencies. The identified resonance frequencies and the resulting spring coefficients are listed in Table 2.

Figure 13Amplitude of the transfer function H(f) for excitation in Fx and Fy.

Table 2Resonance frequencies, spring constants and damping coefficients for force excitation in x and y.

In , the damping coefficient of the spring-damper model is calculated using the ratio of the resonance frequency and the width of the resonance peak of the transfer function at half the amplitude:

$\begin{array}{}\text{(10)}& \frac{{f}_{\mathrm{0}}}{\mathrm{\Delta }{f}_{\mathrm{1}/\mathrm{2}}}\approx \frac{\mathrm{1}}{b}\sqrt{k\cdot m}.\end{array}$

For the given transfer function in Fig. 12, the damping coefficients are calculated to be b1=699.18kg s−1 and b2=1566.01kg s−1. The damping coefficients for excitation in x and y are listed in Table 2.

## 5.3 Determination of the dynamic sensitivity and phase shift

To calculate the dynamic sensitivity of the sensor, sinusoidal excitations at different frequencies are used. Measurements are performed at 12 different frequencies from 53.7 to 1020.3 Hz. The sensitivity S is calculated from the ratio of the amplitude of the sensor signal $\stackrel{\mathrm{^}}{s}$, the top mass mt and the amplitude of the top mass acceleration ${\stackrel{\mathrm{^}}{a}}_{\mathrm{t}}$:

$\begin{array}{}\text{(11)}& S=\frac{\stackrel{\mathrm{^}}{s}}{{\stackrel{\mathrm{^}}{a}}_{\mathrm{t}}\cdot m}.\end{array}$

The acceleration ${\stackrel{\mathrm{^}}{a}}_{\mathrm{t}}$ is measured at eight different positions on the top mass using the laser vibrometer.

The phase shift between the acceleration of the top mass at and the sensor signal s is calculated from the phase of the sine fitting functions of the sensor and the vibrometer signal. The synchronization between the signal sources is performed using the camera trigger signal as described in Sect. 3.2.

Figure 14 shows the mean value of the dynamic sensitivity and the phase shift over the eight acceleration measurements and the resulting standard deviation. The static sensitivity from the sensor calibration data is given as a reference. A rocking mode occurs at a frequency of 375.9Hz, which leads to higher deviations for the acceleration measurements in the outer areas of the load mass. The rocking motion can also be seen in the higher standard deviation of the phase shift values. Up to the excitation frequency of 1000Hz, the sensitivity drops to 95.3% of the static value. The phase shift drops to −0.25rad at 800Hz. For higher frequencies up to 1000Hz, the phase shift does not change significantly.

Figure 14Sensitivity of Fz and phase shift between top mass acceleration and sensor signal. Error bars represent the standard deviation for the acceleration at eight different positions on the top mass.

6 Conclusions

An extended set-up for the dynamic calibration of multi-component sensors for force and moment measurement has been described. It is based on the periodic acceleration of a sensor and an attached load mass on an electrodynamic shaker. In comparison to earlier works, the design of the load mass and adapting elements was focused on activating single force and moment components where possible. Force components can be activated using a load mass with its centre of gravity in the origin of the sensor coordinate system. For the moment components, a beam cross with a movable load mass was designed.

To identify the movement of the load mass with 6 degrees of freedom, a photogrammetric set-up was installed in the existing laboratory set-up. Two stereo camera systems observe the load mass and the shaker surface. From the displacement of the observed surfaces and the time stamp of the camera images, accelerations can be calculated. With the additional information of the excitation frequency and reference measurements using a laser interferometer, the displacement of the set-up can be calculated even with camera frame rates lower than the excitation frequency. From an experimental evaluation, a deviation of 0.68% of the acceleration was achieved for an excitation frequency of fex=106.25Hz.

The dynamic parameters of the sensor are identified based on a three-dimensional mass-spring-damper system. The internal mass of the sensor was calculated in a static manner from measurements at different rotation angles around the z axis. The resulting internal masses of x=0.77 and y=0.768kg differ by 0.3%. The spring constant was calculated from a transfer function as the ratio of the top mass acceleration and the acceleration of the shaker surface. For two different load masses, the resulting spring constant differs by 4.4% for Fz. The spring constant for excitation in the x and y directions shows small deviations of less than 0.5% for the different orientations with the same load mass, while the difference at different load masses for the same orientation rises to up to 15%. It has to be pointed out that the adapting elements for excitation in the x and y directions are still subject to optimization. At frequencies from 350 to 500Hz, pitching motions occur which show an influence on the measurements.

The dynamic sensitivity of the sensor was analysed for the Fz component. In this set-up, a rocking motion is visible at 375.9Hz, which increases the standard deviation of the measured accelerations on the top mass. Up to a frequency of 1000Hz, the dynamic sensitivity drops to 95.3% of the static value. The phase shift between the sensor signal and the top mass acceleration drops to −0.25rad at frequencies of 800Hz.

As a next step, the analysis performed will be extended to the moment excitation. The previously described beam cross and moment excitation mass are currently being manufactured. Further optimization of the horizontal set-up to reduce pitching motions is suggested to reduce signal crosstalk. Alternatively, a three-dimensional shaker set-up can be used in combination with the adapting element for Fz excitation. Finally, the measurement uncertainty for the different set-ups will be analysed.

Data availability
Data availability.

Datasets used in Sects. 3.3 and 5.2 are available at the PTB's public repository https://doi.org/10.7795/720.20181105 (Nitsche et al., 2018).

Author contributions
Author contributions.

JN designed and installed the setup, performed measurements and evaluation and wrote the manuscript. RK and RT raised the research problem, supervised the research, and discussed and proofread the paper.

Competing interests
Competing interests.

Author Rainer Tutsch is a member of the editorial board of the journal.

Acknowledgements
Acknowledgements.

The authors gratefully acknowledge the funding of this work by the Deutsche Forschungsgemeinschaft (DFG) under grants Tu 135/24 and Ku 3367/1. The authors thank Thomas Bruns and Leonard Klaus of the PTB for assistance with the dynamic moment model, sine fitting and hardware.

Edited by: Ulrich Schmid
Reviewed by: two anonymous referees

References

Arun, K. S., Huang, T. S., and Blostein, S. D.: Least-squares fitting of two 3-D point sets, IEEE T. Pattern Anal., PAMI-9, 698–700, https://doi.org/10.1109/TPAMI.1987.4767965, 1987. a

Bartoli, C., Beug, M. F., Bruns, T., Elster, C., Esward, T., Klaus, L., Knott, A., Kobusch, M., Saxholm, S., and Schlegel, C.: Traceable dynamic measurement of mechanical quantities: Objectives and first results of this European project, Int. J. Metrol. Qual. Eng., 3, 127–135, https://doi.org/10.1051/ijmqe/2012020, 2012. a

Baumgarten, S., Kahmann, H., and Röske, D.: Metrological characterization of a 2 kN m torque standard machine for superposition with axial forces up to 1 MN, Metrologia, 53, 1165–1176, https://doi.org/10.1088/0026-1394/53/5/1165, 2016. a

Bruns, T., Kumme, R., Kobusch, M., and Peters, M.: From oscillation to impact: The design of a new force calibration device at PTB, Measurement, 32, 85–92, https://doi.org/10.1016/S0263-2241(01)00048-3, 2002. a

Dewetron: DAQP-BRIDGE-B Module Technical reference manual, available at: https://ccc.dewetron.com (last access: 10 October 2018), 2012. a

DIN 51309: Werkstoffprüfmaschinen – Kalibrierung von Drehmomentmessgeräten für statische Drehmomente, National standard, DIN, 2005. a

Eggert, D. W., Lorusso, A., and Fisher, R. B.: Estimating 3-D rigid body transformations: A comparison of four major algorithms, Mach. Vision Appl., 9, 272–290, https://doi.org/10.1007/s001380050048, 1997. a

ISO 376: Metallic materials: Calibration of force-proving instruments used for the verification of uniaxial testing machines, International standard, ISO, Geneva, 2011. a

Kim, G.-S.: The development of a six-component force/moment sensor testing machine and evaluation of its uncertainty, Meas. Sci. Technol., 11, 1377–1382, https://doi.org/10.1088/0957-0233/11/9/318, 2000. a

Klaus, L.: Model parameter identification from measurement data for dynamic torque calibration – Measurement results and validation, ACTA IMEKO, 5, 55–63, https://doi.org/10.21014/acta_imeko.v5i3.318, 2016. a

Kumme, R.: Untersuchung eines direkten Verfahrens zur dynamischen Kalibrierung von Kraftmessgeräten: Ein Beitrag zur Verringerung der Messunsicherheit, PTB-Bericht PTB-MA-48, Wirtschaftsverl. NW, Verl. für Neue Wiss, Bremerhaven, 1996. a

National Instruments: PXI-449x Specifications, available at: http://www.ni.com/pdf/manuals/372125f.pdf (last access: 10 October 2018), 2012. a

Nitsche, J., Baumgarten, S., Petz, M., Röske, D., Kumme, R., and Tutsch, R.: Measurement uncertainty evaluation of a hexapod-structured calibration device for multi-component force and moment sensors, Metrologia, 54, 171–183, https://doi.org/10.1088/1681-7575/aa5b66, 2017a.  a

Nitsche, J., Röske, D., and Tutsch, R.: Influence of coordinate system alignment on the calibration of multi-component force and moment sensors, in: IMEKO TC3, TC5, TC22 International Conference 2017, 2017b. a

Nitsche, J., Bruneniece, S., Kumme, R., and Tutsch, R.: Design of a calibration setup for the dynamic analysis of multi-component force and moment sensors, in: Proceedings of the XXII IMEKO World Congress, 2018a. a

Nitsche, J., Kumme, R., and Tutsch, R.: Dataset for Dynamic characterization of multi-component sensors for force and moment, Physikalisch-Technische Bundesanstalt (PTB), https://doi.org/10.7795/720.20181105, 2018b.

Pan, B., Qian, K., Xie, H., and Asundi, A.: Two-dimensional digital image correlation for in-plane displacement and strain measurement: A review, Meas. Sci. Technol., 20, 062001, https://doi.org/10.1088/0957-0233/20/6/062001, 2009. a

Park, Y.-K., Kumme, R., and Kang, D.-I.: Dynamic investigation of a three-component force-moment sensor, Meas. Sci. Technol., 13, 654–659, https://doi.org/10.1088/0957-0233/13/5/302, 2002a. a

Park, Y.-K., Kumme, R., and Kang, D.-I.: Dynamic investigation of a binocular six-component force-moment sensor, Meas. Sci. Technol., 13, 1311–1318, https://doi.org/10.1088/0957-0233/13/8/320, 2002b. a

Park, Y.-K., Kumme, R., Röske, D., and Kang, D.-I.: Column-type multi-component force transducers and their evaluation for dynamic measurement, Meas. Sci. Technol., 19, 115205, https://doi.org/10.1088/0957-0233/19/11/115205, 2008. a, b, c

Röske, D., Peschel, D., and Adolf, K.: The generation and measurement of arbitrarily directed forces and moments: The project of a multicomponent calibration device based on a hexapod structure, in: Proceedings of the 17th International Conference of IMEKO TC3, 339–349, 2001. a

Schlegel, C., Kieckenap, G., Glöckner, B., Buß, A., and Kumme, R.: Traceable periodic force calibration, Metrologia, 49, 224–235, https://doi.org/10.1088/0026-1394/49/3/224, 2012. a, b, c, d, e, f, g, h

Schleichert, J., Rahneberg, I., and Fröhlich, T.: Dynamische Kalibrierung eines Mehrkomponentensensors für Kraft und Drehmoment, Tech. Mess., 83, 131–138, https://doi.org/10.1515/teme-2015-0090, 2016. a

Vlajic, N. and Chijioke, A.: Traceable dynamic calibration of force transducers by primary means, Metrologia, 53, S136–S148, https://doi.org/10.1088/0026-1394/53/4/S136, 2016. a