Published: 01 December 2017

Dynamic displacement estimation using data fusion

Sabine Upnere1
Normunds Jekabsons2
1Technical University, Institute of Mechanics, Riga, Latvia
1Ventspils University College, Ventspils, Latvia
2University of Latvia, Faculty of Physics and Mathematics, Riga, Latvia
Corresponding Author:
Sabine Upnere
Views 326
Reads 190
Downloads 2488

Abstract

The paper describes a Kalman filtering technique for dynamic displacement estimation using accelerometer and laser sensor measurements. Data fusion of measurements from multiple sensors can give the more accurate results because of different advantages of sensors. Since the acceleration and displacement have different sampling rates, the multi-rate Kalman filter is applied. The filter is expanded with the fixed interval smoother to improve reconstruction accuracy of displacements. A modelled signal consisting of two sinus functions and Gaussian distributed noise is used to validate developed state-space model.

1. Introduction

The vibration analysis is useful for the surveillance of a structure safety and structural health monitoring. The high-velocity fluid flow interaction with structural components can induce self-excited vibrations of the system elements in many engineering applications, for example, water-cooling systems or heat exchangers. In the exploration process of flow-induced vibrations in the rod bundle, two main interesting parameters are the rod oscillation frequency depending on surrounding flow velocity and the relationship between the rod position and the flow velocity.

Most commonly measured parameter for the characterisation of the structure response on the flow excitation is acceleration due to its relatively simple implementation. Using direct numerical double-integration of acceleration measurements to obtain displacements have to encounter with an integration error, which is not linearly accumulated during integration [1]. Traditional methods, such as baseline correction techniques, two-baseline schemes or others, usually cannot give acceptable accuracy to evaluate displacements in the case when the mean value of acceleration measurements is non-zero or there are nonlinearities in displacements [2].

An alternative approach is to apply techniques of Kalman filters [3] analysing data from multiple sensors. In particular case, accelerometers and laser sensors are used to describe dynamics of a structure. Data fusion of measurements from different sensors can give the more accurate results, for example, accelerometer provides a better description of high-frequency process whereas laser sensors are more suitable to explore the low-frequency range. Due to sensors specific, the measurements of laser sensor and accelerometer are usually taken at different sampling rates. Smyth and Wu [4] used multi-rate Kalman filter data fusion to evaluate displacements. Measurement process was modelled using the state-space equation for velocity and acceleration. The technique of real-time estimation of dynamic displacement based on measured velocity and displacement from two different sensors are described by Kim and Sohn [2]. They construct a state-space model for displacement and the total error.

In this paper is described multi-rate Kalman filter with fixed interval smoother for a state-space model for displacement, velocity and the accelerometer error. A modelled signal consisting of two sinus functions and Gaussian distributed noise is used to validate developed algorithm.

2. Displacement estimation using acceleration and displacement measurements

Non-contact accelerometer and laser sensors measurements can be done to analyse the flow velocity influence on the rod oscillation frequency and amplitude. Usually, accelerometer sensors have a higher sampling rate. However double integration leads to drift of the calculated mean value. The Kalman filter can be used to avoid an integration error. We assume that Kalman filter can be applied to each axis in the plane of the motion separately.

2.1. Double integration

It is assumed that the following mathematical model Eq. (1-2) can describe the rod movement in the discrete-time domain:

1
x˙k+1=x˙k+x¨kt+εx¨(k)t,
2
xk+1=xk+x˙kt+12x¨kt2+12εx¨(k)t2,

where x˙(k) and x(k) are calculated velocity and displacement at the time step k, t is the time interval for acceleration measurements and εx¨(k) is acceleration measurement error which is a combination of offset bias and zero-mean stochastic noise process, w(k) [5].

2.2. Formulation of the Kalman filter

The dynamic displacement evaluates using displacement x(k), velocity x˙(k) and error εx¨(k) as state variables:

3
x(k)=x(k)x˙(k)εx¨(k),

Equations of the state-space model for multi-rate data fusion of laser and accelerometer measurements are as follows:

4
x(k)=Ax(k-1)+Buk-1+wk-1,
5
z(k)=Hx(k)+v(k),

where uk is measured acceleration, zk is measured displacement, vk is measurement noise of a laser sensor and:

6
A=1t0.5t201t001, B=0.5t2t0, H=100.

The system equation Eq. (4) and observation equation Eq. (5) estimates in two steps: predictor, Eq. (7-8) and corrector Eq. (10-12):

7
xfk-=Axfk-1++Buk-1,
8
Pfk-=APfk-1+AT+Q, k=1,2,, N,

where:

9
Q=qBBT,

and q is the covariance of the noise w(k):

10
xfk+=xfk-+Kfkzk-Hkxfk-,
11
Pfk+=I-KfkHkPfk-,
12
Kfk=Pfk-HkTHkPfk-HkT+r-1,

where r is the covariance of the measurement noise v(k).

Since the acceleration and displacement have different sampling rates, the optimal estimates of the state variables can process by multi-rate Kalman filter [4]. In time steps when displacement measurements are not available, only the prediction step is performed. The large sampling interval of the laser sensor can lead to drift of the displacement estimation. Smoothing can produce a better assessment [4].

Fixed-interval smoothing by Rauch-Tung-Striebel (RTS) [6] is used to correct the filtered signal. The RTS algorithm combines forward Kalman filtering and backward filtering Eqs. (15-17) with the following initialisation:

13
xNs=xfk+, PNs=Pfk+,
14
xks=xfk++Kksxk+1s-xfk+1-,
15
Pks=Pfk+-KksPfk+1--Pk+1sKksT, k=N-1,,0,

where the smoother gain is:

16
Kks=Pfk+AkTPfk+1--1.

3. Test cases

Two modelled input signals one for displacements and one for acceleration are used to validate described algorithm. Dynamic displacement due to the flow-induced vibrations of the rod is estimated using data fusion of experimentally obtained acceleration and displacement measurements.

3.1. Case with modelled signal

For displacement simulation combination of two sinusoidal signals is used. The input of displacements is approximated as a quantized signal to resemble real measurements:

17
x= -16π2sin6πt-7080π2sin80πt.

Accelerometer signals obtained using the second-order derivative of the displacement Eq. (8). Noise with a normal distribution is applied to both input signals. The modelled signals can see in Fig. 1.

Fig. 1Modelled input signals: a) the laser sensor, b) accelerometer

Modelled input signals: a) the laser sensor, b) accelerometer

a)

Modelled input signals: a) the laser sensor, b) accelerometer

b)

Fig. 2Calculated output signals: a) and corresponding error b)

Calculated output signals: a) and corresponding error b)

a)

Calculated output signals: a) and corresponding error b)

b)

Displacement sampling time is two times larger than acceleration sampling. Example of filtered signals and the difference between modelled xm and filtered xc signals are shown in Fig. 2.

Maximum amplitudes of displacement using the filter without and with smoother at two noise levels are summarised in Table 1.

In all cases, the filter gives larger maximum amplitude as a modelled signal. Comparison of maximum values shows that reconstructions of the filter with smoother provide better agreement with expected values. Smoother decreases error approximately twice.

Table 1Comparison of maximum displacements depending on noise

Maximum noise
(% of maximum amplitude)
Kalman filter
(% of maximum amplitude)
Kalman filter with smoother
(% of maximum amplitude)
13.09
+4.6
+2.0
13.66
+3.1
+1.5
13.75
+4.4
+2.0
26.04
+6.9
+3.4
24.86
+8.6
+5.1
27.83
+6.5
+4.1

3.2. Experimental case

Experimental session was started with measurements at zero flow rate. Thus we presume that only noise was presented in the time signal. Typical distributions of the signal (noise) measured by accelerometers (Acc3 and Acc4) are shown in Fig. 3. In all presented cases the signal distribution is similar to Gaussian distribution (red line in the figure).

Fig. 3Signal (noise) distribution using experimental accelerometer data

Signal (noise) distribution using experimental accelerometer data

a) Acc3

Signal (noise) distribution using experimental accelerometer data

b) Acc4

Displacements calculations using acceleration time history can also be done by so-called Omega Arithmetic (OA) method [7]. There are found no differences between displacement reconstruction using multi-rate Kalman filter or OA when flow velocity is steady. However, multi-rate Kalman filter, unlike OA, can reconstruct the rod displacements in the unsteady regime when flow velocity is slowly increasing or decreasing in time, see Fig. 4.

Fig. 4Reconstruction of displacement: a) when flow velocity is increasing using multi-rate Kalman filter, b) omega arithmetic algorithm

Reconstruction of displacement: a) when flow velocity is increasing  using multi-rate Kalman filter, b) omega arithmetic algorithm

a)

Reconstruction of displacement: a) when flow velocity is increasing  using multi-rate Kalman filter, b) omega arithmetic algorithm

b)

4. Conclusions

The fusion of an acceleration and displacement measurements using smoothing based multi-rate Kalman filter to estimate dynamic displacement is implemented and tested. Validation of proposed method is done by a modelled signal with additional normal distributed noise. Comparison of obtained results shows that smoother reduces reconstruction error approximately twice in given cases. The multi-rate Kalman filter can reconstruct time history when the flow velocity is slowly increased or decreased in unsteady flow rate measurements. The Omega Arithmetic method is less successful to reconstruct the rod coordinate oscillations in flow transients.

References

  • Thong Y. K., et al. Numerical double integration of acceleration measurements in noise. Measurement, Vol. 36, Issue 1, 2004, p. 73-92.
  • Kim K., Sohn H. Dynamic displacement estimation by fusing LDV and LiDAR measurements via smoothing based Kalman filtering. Mechanical Systems and Signal Processing, Vol. 82, 2017, p. 339-355.
  • Kalman R. E. A new approach to linear filtering and prediction problems. Journal of Basic Engineering, Vol. 82, 1960, p. 35-45.
  • Smyth A., Wu M. Multi-rate Kalman filtering for the data fusion of displacement and acceleration response measurements in dynamic system monitoring. Mechanical Systems and Signal Processing, Vol. 21, Issue 2, 2007, p. 706-723.
  • Kim J., Kim K., Sohn H. Autonomous dynamic displacement estimation from data fusion of acceleration and intermittent displacement measurements. Mechanical Systems and Signal Processing, Vol. 42, 2014, p. 194-205.
  • Rauch H. E., Tung F., Striebel C. T. Maximum likelihood estimates of linear dynamic systems. AIAA Journal, Vol. 3, Issue 8, 1965, p. 1445-1450.
  • Mercer C. Acceleration, Velocity and Displacement Spectra. Omega Arithmetic, 2006.

Cited by

Twin Rotor Damper for Human-Induced Vibrations of Footbridges
Richard Terrill | Richard Bäumer | Katrien Van Nimmen | Peter Van den Broeck | Uwe Starossek
(2020)

About this article

Received
06 November 2017
Accepted
16 November 2017
Published
01 December 2017
SUBJECTS
Mathematical models in engineering
Keywords
Kalman filter
smoother
accelerometer sensor
laser sensors
sensor data fusion
Acknowledgements

Authors would like to express gratitude for experimental data to Dr. Sergejs Dementjevs from Paul Scherrer Institute.