1. Introduction
For large-scale structures, serious structural damage may occur when vibrating deformation induced by external loads (such as traffic and winds) exceeds the designed bearing capacity. Accurate deformation information is a prerequisite for timely alerts to prevent unnecessary casualties and losses. SHM is a technology for obtaining deformation information and estimating the health condition and the structural characteristics of civil structures and infrastructures [
1,
2]. On this basis, an SHM strategy first intends to measure responses such as displacements over time, and then obtain insightful information about the current or unknown condition of a civil structure by computational techniques either in time domain or in frequency domain [
3,
4]. Various types of sensor devices such as linear variable differential transformers, optical fiber sensors, smartphones, vision cameras, and radars are used in SHM to measure the displacement responses, but dynamic deformation monitoring widely adopts GNSS and accelerometers for SHM non-stop or with high periodicity but not requiring gluing/embedding sensors into the structure [
5,
6,
7,
8]. The RTK-GNSS can achieve a subcentimeter-level measurement accuracy, and it is often selected for structural displacement estimation of large-scale structures such as long-span bridges and high-rise buildings, which usually have at least centimeter-level displacements [
9,
10,
11]. However, the accuracy of GNSS positioning is typically limited due to environmental perturbations such as occlusion, diffraction, and reflection. In addition, the GNSS-based monitoring method is insensitive to high-frequency vibration signals because the measurement noise is relatively large [
12]. The accelerometer method can obtain high-accuracy dynamical displacements over short time periods through the double integration of acceleration measurements. However, significant cumulative errors are ultimately present when the accelerometer is applied to monitor long-period quasi-static displacements [
13].
GNSS and accelerometers are complementary and therefore have been commonly integrated in monitoring applications. For example, Meng et al. [
14] successfully identified a peak vibration frequency of 10.05 Hz of Nottingham Wilford Bridge by using an GNSS/accelerometer integrated system. Kim et al. [
15] combined GPS RTK with an accelerometer in monitoring Yeongjong Grand Bridge and achieved a displacement accuracy of 2 mm in the vertical direction. Xin et al. [
16] integrated a strong motion accelerometer with GNSS PPP for seismic deformation monitoring. The root-mean-square error (RMSE) of the difference with respect to the reference was 2 mm, with a cross-correlation coefficient of 0.99. Although GNSS/accelerometer integration has been widely validated in previous studies, the tilt of the accelerometer can cause baseline errors and distort the captured peak amplitude of the displacement signal [
14,
16].
In addition to displacement, rotational deformation is another critical piece of information for structural health monitoring. Although the accelerometer can measure the rotational deformation by sensing changes in gravitational acceleration [
17], the accuracy of accelerometer-derived rotation is typically limited because accelerometers measure specific forces that involve gravitational and vibrational acceleration. To separate the different types of acceleration and finally determine the rotational deformation, it is necessary to use gyroscopes [
18]. Rossi et al. [
19] corrected the impact of rotation on GNSS/accelerometer monitoring results using rotation information recorded using an inertial measurement unit (IMU). The RMSE of the corrected displacement monitoring results was reduced by half to about 1–2 mm. de Alteriis et al. [
20] designed a low-cost and real-time monitoring device that integrated GNSS and Micro Electro Mechanical System (MEMS) IMU. Experimental results showed that the system was able to estimate position and attitude under high-frequency vibrations, and the monitoring results were in good agreement with the fiber-optic sensors. Geng et al. [
21] implemented a seismic monitoring system with six degrees of freedom by integrating an accelerometer, gyroscope, and GNSS. In a waveform simulation experiment of the 2010 MW 7.2 EI Mayor-Cucapah earthquake, the displacement obtained using the proposed method achieved an accuracy 68% higher than that obtained from the traditional accelerometer/GNSS integration.
Most previous studies of integrated GNSS/IMU deformation monitoring are based on single-antenna GNSS. Although improved accuracy and robustness have been reported, there are still some obvious limitations of single-antenna GNSS/IMU. For example, the GNSS signal environment in deformation monitoring is typically complex, which worsens the GNSS signal quality and ultimately the positioning accuracy. Moreover, the motion of the monitored structures is generally limited, resulting in a low attitude observability, especially in the heading direction [
22,
23]. The accumulation of attitude errors eventually affects the monitoring accuracy, especially when low-cost IMUs are used [
24]. Multi-antenna GNSS offers a solution to this problem. Double-difference carrier phases are formed between multiple GNSS antennas to estimate inter-antenna baselines. When carrier phase ambiguities are fixed, the baselines can be accurately inverted to attitude information [
25]. Zhang et al. [
26] evaluated the positioning performance of the traditional single-antenna and dual-antenna GNSS/IMU. They reported that the dual-antenna GNSS/IMU achieved mm level positioning accuracy, which was about 50% higher than the single-antenna GNSS/IMU. Li et al. [
27] applied the length constraints of lever arms in multi-antenna GNSS/IMU for monitoring the settlement of high-speed railway, and achieved a positioning accuracy of 1–2 mm. However, multi-antenna GNSS/IMU integration has been rarely studied in structural deformation monitoring.
In this paper, we propose a tight integration of multi-antenna GNSS and IMU for vibration monitoring. The GNSS observations are tightly coupled to the IMU by transforming their phase centers to the IMU center. Thus, the integration reduces the number of unknown parameters, enhances the geometric strength of the model, and increases the observability of attitude errors. To evaluate the performance of the proposed method, different types of vibrational signals are simulated using a shaking table and a laser rangefinder recording as a reference.
Section 2 presents the mathematical model of the tight integration of multi-antenna GNSS and IMU, including the measurement model and the system model. In
Section 3, the experimental configuration, results, and analysis are described.
Section 4 concludes this work.
3. Experiments and Analysis
To assess the performance of the tight integration of multi-antenna GNSS and IMU for vibration deformation monitoring, a simulation vibration test was conducted on 13 November 2023 on a half-open platform on the 6th floor of Hong Kong Polytechnic University.
Figure 2a illustrates the environment of the experimental site. The monitoring accuracy of integration schemes of different antennas and IMU were analyzed, and the vibration displacement information recorded using a laser rangefinder with a ranging accuracy of 0.01 mm, and a sampling rate of 20 Hz was used as the reference.
Figure 2b shows the instrument deployment for the vibration experiment. A shaking table was fixed on a tripod. A rigid triangular platform with one IMU and four GNSS antennas was mounted on the table, and the length of inter-antenna baselines was 1 m. In deformation monitoring, the observability of attitude errors is typically low because of the limited dynamics (including rotation and acceleration) of the deforming objects. The three-dimensional attitude can be determined using at least three non-collinear antennas. Therefore, this study adopts the triple-antenna triangle configuration to strengthen the GNSS/IMU positioning model so that the estimation of parameters such as attitude and IMU biases can be enhanced. A Leica GR25 GNSS geodetic receiver was attached to antenna 2, two Trimble R12 GNSS geodetic receivers were placed for rover antennas 1 and 3, respectively, and a Trimble Net R9 geodetic receiver was used for the base station. We use the configuration of the three non-collinear antennas to enhance the three-dimensional attitude determination, considering that the observability of attitude is typically low due to the limited dynamics of the deforming objects. The base station was located on the roof of the building at a distance of about 50 m. The GNSS sampling rate was 20 Hz. The IMU adopted was iXBLUE ATLANS-C, which has triple-axis accelerometers and gyroscopes, and a self-contained GNSS receiver for time synchronization (the corresponding antenna was installed at the center of the triangular platform). The error characteristics of IMU were estimated using Allan variance analysis [
29], as shown in
Table 1.
A suite of C++ libraries was developed for the data processing, where IMU records are tightly integrated with GNSS observations in double-difference relative kinematic positioning. The double-difference observation equations are formed between the primary antenna and the reference station, as well as the auxiliary antennas. GPS and BDS observations were used for positioning in the experiments. In short baselines, the frequency-dependent ionospheric delays are considerably canceled, and therefore only single-frequency (L1/B1) observations were used in the experiment to reduce the computation of high-rate data. The cut-off angle was set to 10 degrees, and observations were weighted using an elevation-based cosine function.
Figure 3 shows the sky distribution of the observed GPS and BDS satellites, where G denotes the GPS satellite and C the BDS satellite. It can be seen that during the test, about 8–9 GPS satellites and 16–17 BDS satellites were observed by Antennas 1, 3, and the base station. Satellite signals from the southeast direction with low elevation angles (below 30 degrees) were not received due to the blocking of nearby buildings (
Figure 2). Note that Antenna 2 could observe 9 GPS satellites but only 8 BDS-2 satellites because the firmware of the receiver was outdated.
Table 2 lists eleven positioning schemes that were tested in our experiment, including three single-antenna GNSS-only schemes, one triple-antenna GNSS-only scheme (using the average of the three single-antenna GNSS-only solutions), and seven GNSS/IMU integration schemes with different antennas.
Figure 4 shows the 2D horizontal displacement time series recorded with the laser rangefinder, and the horizontal displacements derived by relative positioning using schemes G1, G2, and G3. As shown in
Figure 4a, the shaking table was first kept stationary for about 10 min, which was referred to as Period 1, then generated two 1 min horizontal vibrations with frequencies of 0.50 Hz and 0.67 Hz and amplitudes of 19 mm and 17 mm, respectively. After that, a vibration with a mixed frequency of 0.11 Hz and 0.34 Hz was loaded to the shaking table for about 4 min. Then the shaking table generated a vibration with a frequency of 0.80 Hz and an amplitude of 15 mm. At last, the shaking table returned to the original position in 4 min and then repeated the vibration process. The whole vibration process was named Period 2. We can see that all eight vibration events were successfully captured by the three single-antenna GNSS-only schemes. In both Period 1 and Period 2, the G2 solution was slightly nosier than those of G1 and G3, mainly because of the relatively fewer satellites tracked by Antenna 2 (see also
Figure 3).
Figure 5 shows the errors of the horizontal displacement of different positioning schemes. We can see from
Figure 5a that during Period 1, the positioning errors of single-antenna GNSS-only schemes (i.e., G1, G2, and G3) fluctuate within about 10 mm, and the fluctuations become larger (within 20 mm) in Period 2, maybe due to the coupling of GNSS systematic errors and vibrations. The results of G2 are nosier than those of G1 and G3; for instance, large deviations are present in G2 during 123,400–123,600 s, perhaps due to relatively weaker geometry and undetected outliers. Integrating IMU with single-antenna GNSS data significantly reduces the positioning errors, which become obviously smaller and smoother. However, many low-frequency fluctuations remain visible in the error time series. This implies that the dominant influence is highly correlated with the single-antenna GNSS-only solution errors (
Figure 5a). The benefit of multiple antennas, as can be seen in GI12, GI13, GI23, and GI123 solutions, is evident when observations from more GNSS antennas are gradually integrated (
Figure 5b–d). Low-frequency systematic errors are mitigated and many error peaks are reduced (see, e.g., the period of 123,400–123,600 s). This demonstrates that the integration of multiple antennas can enhance the geometry of positioning. Typically, scheme G123, that is based on the average of three corresponding single-antenna GNSS-only solutions (
Figure 5d), can achieve similar performance with dual-antenna GNSS/IMU integration (
Figure 5c). And the triple-antenna GNSS/IMU solution GI123 further reduces the positioning errors compared with G123 by adding IMU (
Figure 5d).
Figure 6 shows the histograms of the displacement errors of different positioning schemes. For simplicity, only four typical schemes are shown (including G1, GI1, GI12, and GI123). We can see that the range of positioning errors of G1 is [−18.18, 10.59] mm, which is reduced to [−14.25, 12.79] mm when GNSS/IMU integration is adopted (i.e., GI1). When dual-antenna GNSS/IMU is used (i.e., GI12), the error range is further reduced to [−12.73, 8.38] mm. The triple-antenna GNSS/IMU solution achieves the best results, with the error range of [−8.24, 11.31] mm, and more errors are close to zero.
Table 3 presents the statistics of coordinate standard deviation errors (STDE) for all eleven positioning schemes. For the whole experiment period (including the static period and vibration period), the single-antenna GNSS-only positioning STDE is about 3.88 mm on average, with the STDEs of G1, G2, and G3 being 3.87, 4.21, and 3.57 mm, respectively. When the single-antenna GNSS is tightly integrated with IMU, the corresponding STDE only slightly declines to 3.73 mm (with the STDEs of GI1, GI2, and GI3 being 3.66, 4.12, and 3.41 mm respectively). We can find that the GNSS/IMU positioning performance is related to the positioning quality of the corresponding GNSS antenna used. This effect is less obvious when additional GNSS antennas are used in the integration. The positioning STDE is reduced to about 2.79–3.21 mm by applying dual-antenna GNSS/IMU integration, corresponding to an average accuracy improvement of 21.9% relative to the single-antenna GNSS-only solutions. The triple-antenna GNSS/IMU solution further increases the improvement rate to 33.0% by reducing the STDE to 2.60 mm. Moreover, we notice that the G123 solution is 10.6–24.2% more accurate than the single-antenna solutions (G1, G2, and G3) by means of averaging, and GI123 further reduces the STDE, which is 18.5% smaller than that of G123. This confirms that the use of both IMU and multiple antennas contributes to the improvement in accuracy. As for the statistical results of Period 1 and Period 2, they are similar with those based on the whole data period, while the STDEs of Period 2 increase by about 0.6–1.5 mm compared with those of Period 1, mainly because of the influence of vibrations.
Power spectral density (PSD) is calculated to analyze the spectrum of the displacement error sequences from the laser rangefinder and positioning schemes G1, GI1, GI12, and GI123 in the first four vibrations. The results of other positioning schemes, and other vibrations are similar, and thus, for brevity, they are not presented. As can be seen in
Figure 7, all positioning schemes can successfully identify the main frequencies of the simulated vibrations, with numerically identical peak-frequency values with respect to those based on laser records. Despite this, compared to the PSD of laser records, the GNSS-only solution G1 is obviously noisier in the high-frequency band (e.g., >1 Hz). The GNSS/IMU solutions exhibit lower noise levels compared to G1, and integrating more GNSS antennas gradually enhances such improvement. This benefits the monitoring of vibrations with high frequencies and small amplitudes, which need further study. As for vibration amplitudes, we applied Fast Fourier Transform to estimate them.
Figure 8 shows the errors of amplitudes estimated based on different positioning schemes: they are similar and mostly smaller than 3 mm; those of GI12 present the largest errors among the selected four schemes, perhaps due to the relatively poor data quality of Antenna 2.
4. Conclusions
In SHM, GNSS can be limited due to signal interference and the accelerometers being subject to rotation-induced baseline errors. Traditional single-antenna GNSS/IMU integrations generally have low observation redundancy and model geometry strength. We propose a tightly integrated multi-antenna GNSS/IMU scheme for structural health monitoring. Vibration signals were simulated with a shaking table to validate the monitoring accuracy of the proposed method, using the laser rangefinder recordings as the reference.
Experimental results show that the STDEs of single-antenna GNSS-only solutions G1, G2, and G3 were 3.87, 4.21, and 3.57 mm, respectively. Compared to the results of single-antenna GNSS-only solutions, the single-antenna GNSS/IMU solutions are affected by the GNSS antenna used; thus, the positioning STDE was only slightly reduced from about 3.88 mm to 3.73 mm on average. Specifically, the STDEs of G1, G2, and G3 were 3.66, 4.12, and 3.41 mm, respectively. Adding more GNSS antennas can effectively increase the improvements. The dual-antenna GNSS/IMU integration reduces the STDE to a level of about 3.03 mm with the STDEs of GI12, GI13, and GI23 being 2.79, 3.10, and 3.21 mm, respectively, and positioning results with an STDE of 2.60 mm were finally obtained by triple-antenna GNSS/IMU integration, which is about 32.99%, 30.29%, and 14.19% more accurate than the single-antenna GNSS-only, single-antenna GNSS/IMU, and dual-antenna GNSS/IMU solutions, respectively. The power spectral density analysis demonstrates that vibration frequencies identified based on the displacements from different positioning schemes are nearly identical. Through Fast Fourier Transform analysis, we found that the benefit of integrating multi-antenna GNSS and IMU in estimating the vibration amplitudes is not obvious. This is perhaps because of the use of poor GNSS data from some antennas.
This study validates the effectiveness of the multi-antenna GNSS/IMU integrated model, and indicates the potential of multi-antenna GNSS/IMU integration for SHM. Further study may apply the proposed method in real SHM environments where more serious signal interference, such as occlusion, multipath, and diffraction errors, is generally present.