# Noninvasive Remote Sensing Techniques for Infrastructures Diagnostics

Authors: Angelo Palombo; Stefano Pignatti; Angela Perrone; Francesco Soldovieri; Tony Alfredo Stabile; Simone Pascucci
Journal: International Journal of Geophysics
 (2011)

---

## Abstract

The present paper aims at analyzing the potentialities of noninvasive remote sensing techniques used for detecting the conservation status of infrastructures. The applied remote sensing techniques are ground-based microwave radar interferometer and InfraRed Thermography (IRT) to study a particular structure planned and made in the framework of the ISTIMES project (funded by the European Commission in the frame of a joint Call “ICT and Security” of the Seventh Framework Programme). To exploit the effectiveness of the high-resolution remote sensing techniques applied we will use the high-frequency thermal camera to measure the structures oscillations by high-frequency analysis and ground-based microwave radar interferometer to measure the dynamic displacement of several points belonging to a large structure. The paper describes the preliminary research results and discusses on the future applicability and techniques developments for integrating high-frequency time series data of the thermal imagery and ground-based microwave radar interferometer data.


---

## Body

## 1. Introduction

During the last decades, the advances in nondestructive electromagnetic sensing techniques (optical, thermal, and microwave) operating at different spectral, spatial and time scales make a multidepth, multiresolution and multiscale diagnostics and monitoring activities thus to be used for the status assessment of existing structures [1].These nondestructive techniques are primarily suited for the detection and characterisation of alterations and defects in the near surface of structures [2] and for the displacement monitoring and verification of the structural integrity of buildings and structures [3, 4].As regards the sensing techniques used for this study, the thermal camera technology is a non-destructive investigation technique commonly used in different applications and in particular in the infrastructures and architectural heritage diagnostics [5]. Thermal imagery is used to identify subsurface ununiformities since the heat flow from the surface to the inner structure is affected by the defects and inhomogeneities present in the structure. The InfraRed Termography (IRT) analysis allows diagnostic results ranging from qualitative characterization (e.g., detection) to quantitative analysis (e.g., inhomogeneities characterization). It is worth noting that the application of IRT for detecting and characterizing defects in construction materials can be quite complicated as different physical properties affecting the thermal and optical properties of the structure [6] have to be considered. Among them, the most important are the conductivity, the diffusivity, and the specific heat of the investigated material, whereas the spectral properties such as emissivity, absorption, reflection, and transmission play a key role on the heat distribution in the material. Other properties/characteristics that influence the thermal behaviour of a material are related to its porosity, volumetric mass, and water content.Furthermore, the use of high-frequency imagery thermal camera gives the opportunity of measuring full field deformation and vibration characteristics (i.e., frequency and attenuation). From these measurements different modal parameters of the structure (e.g., natural frequencies, mode shapes, and damping ratios) can be extracted and displayed and used for the definition of a baseline set of dynamic characteristics of the structure. This kind of analysis can be subsequently used for the structural health monitoring and the application of vibration-based damage detection techniques (see e.g., [7, 8]).For the present study, we used high-frequency imagery (HFI) acquisitions for the retrieving of fast displacements and vibrations in the structure.On the other hand, concerning the microwave interferometry technology applied in this study, it has proven to be a powerful remote sensing tool for vibration measurement of structures [9]. It was recently applied for tracking the vibration of bridges excited by vehicular traffic [8, 10] to monitor the displacements of heritage architectural structures such as the two bell towers,Giotto’s Tower and Arnolfo’s Tower, in Florence [11] or the Tower of Pisa [12], for deflection measurements on vibrating stay cables [13], and for in-field dynamic monitoring of engineering structures [4, 14]. All these studies have been carried out by the use of an innovative noncontact microwave interferometer, named IBIS-S (Image By Interferometric Survey of Structures), manufactured by the Italian company IDS S.p.A. This innovative radar system has been also applied for this study and it is based on the Stepped-Frequency Continuous Wave (SF-CW) technique [15] and on the Differential Interferometric technique [16]. The former providing the system with a range resolution capability and the latter allowing the system to evaluate the displacement response of each target detected in the illuminated scenario. IBIS-S sensor is fast and easy to be installed, allowing to operate in all weather conditions, both day and night, even at a significant distance and without the need of installing and wiring sensors. A further advantage obtained by the IBIS-S instrument is the direct measurement of the displacement of a large number of targets simultaneously, in real-time, and with high accuracy.In this paper, the application of different sensors and methodologies to several noninvasive infrastructures near surface diagnostics is presented. The analysis of acquired data on different structures and materials included also time series recorded by the thermal camera and radar on several points of a test bed cement beam.


## 2. Data and Methods

### 2.1. Test Beds

The data for this study were acquired in Montagnole (French Alpes) site in an innovative facility owned by LCPC (Laboratoire Central des Ponts et Chaussées), that is one of the ISTIMES partners.The Montagnole test site is mainly used to certify metallic protection nets that are used in mountains to prevent catastrophic rockslides. For ISTIMES project, it was exploited as a purely research-oriented facility to verify the behaviour of noninvasive sensing techniques during the progressive damage of an on purpose built concrete beam structure. In particular, the experiment has regarded the progressive damage, in different stages, of a concrete beam by means of falling blocks thus allowing the different techniques to be tested in the presence of hazards.In the Montagnole experiment an iron ball of 2.5 Tons was used (Figure1(a)) and for this study we analyze only the four falling actions where the ball has impacted directly on the beam. In particular, drops 1 to 3 have been performed from an altitude of 1 m with respect to the beam and drop 4 from a 5 m of height. The last drop (i.e., number 4) was used to study and diagnose the significant structural damage as induced by a heavy direct impact of the ball onto the armed concrete beam in order to understand the progressive energy release. Therefore, for this experiment, we assumed that such drop rate was convenient to have a refined monitoring of the beam at the beginning and at the end of each daily experiment session.(a) Picture of the experimental site with the heavy load (iron ball) suspended at the height of 1 m over the beam. (b) Example of a frame acquired by HFI FLIR camera. The frame depicts the heavy load just before the impact on the concrete beam. The black arrow indicates the mixed pixels (i.e., background and beam) used for calculating the oscillation frequency after the impact.
(a)(b)


### 2.2. Data

For this study we used the following instrumentation for a rapid and noninvasive diagnostic of the infrastructure status. Specifically, as regards the IRT instrumentation, we used a FLIR SC7900-VL thermal camera (LWIR; 7.7–11.5μm), whereas an IBIS-S microwave radar interferometer operating in the Ku frequency band with a central frequency of 16.75 GHz was used for the microwave interferometry.The FLIR camera features has a “snap shot” integration range from 10μs to 10 ms, which incorporates a high quantic efficiency MCT focal plane array thus ensuring a very high-spectral resolution with a Noise Equivalent Temperature Difference (NEDT) <25 mK. The frame rates are up to 200 FPS (frames per second) in full frame mode (320×256) and 6000 FPS in subwindowing mode. In this study, we used the following optimal configuration of the LWIR camera: 160×128 frames, thus allowing 800 FPS.The IBIS-S microwave radar consists of a sensor module installed on a tripod with a 3D rotating head. An USB interface allows the connection with a control PC, which is used to configure the acquisition parameters, to view in real time the first results, and to perform the storage of the measurements. The two horn antennas transmit the electromagnetic signals in the frequency band of16.75±0.30 GHz and receive the echoes from the target. Sensor module and PC are connected to a battery unit.IBIS-S system can provide a range resolution of 0.5 m, a displacement sensitivity of 0.01 mm, a maximum operational distance of 500 m, and a maximum sampling frequency of 200 Hz [17].


### 2.3. Methods

#### 2.3.1. HFI Thermal Imagery

In the survey of measurements on the Montagnole test site, we took advantage of the thermal camera (LWIR) to obtain very high-frequency thermal imagery on the concrete beam during the direct impact of the falling blocks. The FLIR camera was mounted on a tripod at about 40 m (i.e., in box 3) and 70 m (i.e., near box 5) of distance from the beam. Box 3 was used for the first three drops (height of the falling ball of 1 m); for the drop 4 (height of the falling ball of 5 m), due to safety reasons, the camera was installed near the box 5 (70 m). For all the drops, the acquisition time was of 46 s for a total of 37000 frames.The thermal imagery acquisitions were used to calculate the beam oscillation frequency as derived from the difference between the beam temperature and the background temperature (see Figure1(b)).The method is based on the presence of mixed pixels composed by both the beam and the background as shown by the black arrow in Figure1(b), labelled as i-pixel. The i-pixel temperature measured by the FLIR camera accounts for the time-varying temperature relative to the beam and the background. The following equations were used to retrieve the beam oscillation:(1)c1Tb+c2Tt=Tm,c1+c2=1,
where, c1 and Tb are the fraction of the i-pixel occupied by the background and its temperature, respectively; c2 and Tt are the fraction of the i-pixel occupied by the beam and its temperature, respectively. Tm is the temperature measured by the camera for the i-pixel.The first equation of (1) represents the linear mixing of the temperatures of the beam and the background, the second is the constrain used for the total pixel fraction. The Tb and Tt temperatures, which are unknowns in (1), can be extracted from the temperature values of the neighbouring pixels of the i-pixel. The neighbouring pixels used for calculating Tb and Tt have to be pure pixels (i.e., made up of only beam or target).In (1), c2 represents the oscillation quantity expressed in terms of pixel fraction that varies from 0 to 1. Moreover, when c2 oscillation reaches values near to the maximum and minimum, this represents situations in which the amplitude of oscillation is close to or greater than the maximum measurable for a pixel, respectively. In the case that a pixel is not sufficient to encompass the oscillation, it is possible to solve the problem by considering, using the sum of the neighbouring pixels, “equivalent pixels” having a larger dimension compared to the original one.


#### 2.3.2. Ground-Based Microwave Radar Interferometer

IBIS-S is based on two well-known radar techniques: the Stepped-Frequency Continuous Wave (SF-CW) and the Differential Interferometric [12, 18].The SF-CW technique allows the system to resolve the scenario in the range direction [15, 19], that is, to detect the position of target surfaces placed at different distances from the sensor. It is based on the synthesis and transmission of a burst of N monochromatic pulses equally and incrementally spaced in frequency (with fixed frequency step of Δf) leading to a work bandwidth B:(2)B=(N-1)Δf,
The N monochromatic pulses sample the scenario in the frequency domain similarly to a short pulse with a large bandwidth B. In a SF-CW radar, the signal source dwells at each frequency fx=f0+kΔf(k=0,1,2,…,N-1), the single tone signal is sufficiently long to allow the echoes backscattered by the target to reach the receiver. Hence, the duration of each monochromatic pulse (Tpulse) depends on the maximum distance (Rmax) to be observed in the scenario and is given by,(3)Tpulse≥2Rmax⁡c,
where c represents the velocity of the light.At each sampled time instant, bothI (In-phase) and Q (Quadrature) components of the received signals are acquired so that the resulting data consist of a vector of N complex samples, representing the frequency response measured at N discrete frequencies.By taking the Inverse Discrete Fourier Transform (IDFT) the radar response is reconstructed in the time domain: each complex sample in this domain represents the signal (echo) from a range (distance) interval of lengthcTpulse/2.The amplitude range profile of the radar echoes is then obtained by calculating the magnitude of each bin of the IDFT of acquired vector samples. This range profile gives a one-dimensional mapping of scattering objects in the “view space” as function of their relative distance from the sensor.The differential interferometric technique allows the system to measure the displacements of the structure illuminated by the antenna main beam by comparing the phase information of the backscattered electromagnetic waves collected at different times (see, e.g., [16]).Generally speaking, when a target surface moves with respect to the sensor module (emitting and back-receiving the electromagnetic wave) at least a phase shift arises between the signals reflected by the target surface at different times. Hence, the displacement of the investigated object is determined from the phase shift measured by the radar sensor at the discrete acquisition times. The radial displacementdp (i.e., the displacement along the direction of wave propagation) and the phase shift Δφ are linked by the following:(4)dp∝λ4πΔφ,
where λ is the wavelength of the probing electromagnetic signal.The radar emits a series of electromagnetic waves for the entire measurement period and processes phase information at regular time intervals to find any displacement occurring between one emission and the next one. It is worth underlining that the interferometric technique in (4) provides a measurement of the displacement along the line of sight for all the range bins of the structure illuminated by the antenna beam; once dp has been evaluated, the vertical displacement d can be easily found by making some geometric projections.IBIS-S system was applied in both Montagnole and Musmeci test beds. In the first case it was installed on a hill at a height of about 15 m from the ground level and at a distance of about 70 m from the target. In order to improve the accuracy and the signal-to-noise ratio of the measurements, a corner reflector was fixed on the beam. To acquire the measurements a range sampling frequency of 200 Hz and an acquisition time of about 60 s for each drop were used.In the Musmeci test bed IBIS-S system was placed just under the bridge deck with a perpendicular direction at a distance of about 30 m. A sampling frequency of about 100 Hz and an acquisition time of about 5 minutes were applied.In both cases the use of IBIS-S system allowed the dynamic monitoring of the investigated target without any contact between the sensor and the investigated target. This characteristic was very useful in Montagnole test bed where very restricted security conditions were required in order to carry out the experiment.


## 2.1. Test Beds

The data for this study were acquired in Montagnole (French Alpes) site in an innovative facility owned by LCPC (Laboratoire Central des Ponts et Chaussées), that is one of the ISTIMES partners.The Montagnole test site is mainly used to certify metallic protection nets that are used in mountains to prevent catastrophic rockslides. For ISTIMES project, it was exploited as a purely research-oriented facility to verify the behaviour of noninvasive sensing techniques during the progressive damage of an on purpose built concrete beam structure. In particular, the experiment has regarded the progressive damage, in different stages, of a concrete beam by means of falling blocks thus allowing the different techniques to be tested in the presence of hazards.In the Montagnole experiment an iron ball of 2.5 Tons was used (Figure1(a)) and for this study we analyze only the four falling actions where the ball has impacted directly on the beam. In particular, drops 1 to 3 have been performed from an altitude of 1 m with respect to the beam and drop 4 from a 5 m of height. The last drop (i.e., number 4) was used to study and diagnose the significant structural damage as induced by a heavy direct impact of the ball onto the armed concrete beam in order to understand the progressive energy release. Therefore, for this experiment, we assumed that such drop rate was convenient to have a refined monitoring of the beam at the beginning and at the end of each daily experiment session.(a) Picture of the experimental site with the heavy load (iron ball) suspended at the height of 1 m over the beam. (b) Example of a frame acquired by HFI FLIR camera. The frame depicts the heavy load just before the impact on the concrete beam. The black arrow indicates the mixed pixels (i.e., background and beam) used for calculating the oscillation frequency after the impact.
(a)(b)


## 2.2. Data

For this study we used the following instrumentation for a rapid and noninvasive diagnostic of the infrastructure status. Specifically, as regards the IRT instrumentation, we used a FLIR SC7900-VL thermal camera (LWIR; 7.7–11.5μm), whereas an IBIS-S microwave radar interferometer operating in the Ku frequency band with a central frequency of 16.75 GHz was used for the microwave interferometry.The FLIR camera features has a “snap shot” integration range from 10μs to 10 ms, which incorporates a high quantic efficiency MCT focal plane array thus ensuring a very high-spectral resolution with a Noise Equivalent Temperature Difference (NEDT) <25 mK. The frame rates are up to 200 FPS (frames per second) in full frame mode (320×256) and 6000 FPS in subwindowing mode. In this study, we used the following optimal configuration of the LWIR camera: 160×128 frames, thus allowing 800 FPS.The IBIS-S microwave radar consists of a sensor module installed on a tripod with a 3D rotating head. An USB interface allows the connection with a control PC, which is used to configure the acquisition parameters, to view in real time the first results, and to perform the storage of the measurements. The two horn antennas transmit the electromagnetic signals in the frequency band of16.75±0.30 GHz and receive the echoes from the target. Sensor module and PC are connected to a battery unit.IBIS-S system can provide a range resolution of 0.5 m, a displacement sensitivity of 0.01 mm, a maximum operational distance of 500 m, and a maximum sampling frequency of 200 Hz [17].


## 2.3. Methods

### 2.3.1. HFI Thermal Imagery

In the survey of measurements on the Montagnole test site, we took advantage of the thermal camera (LWIR) to obtain very high-frequency thermal imagery on the concrete beam during the direct impact of the falling blocks. The FLIR camera was mounted on a tripod at about 40 m (i.e., in box 3) and 70 m (i.e., near box 5) of distance from the beam. Box 3 was used for the first three drops (height of the falling ball of 1 m); for the drop 4 (height of the falling ball of 5 m), due to safety reasons, the camera was installed near the box 5 (70 m). For all the drops, the acquisition time was of 46 s for a total of 37000 frames.The thermal imagery acquisitions were used to calculate the beam oscillation frequency as derived from the difference between the beam temperature and the background temperature (see Figure1(b)).The method is based on the presence of mixed pixels composed by both the beam and the background as shown by the black arrow in Figure1(b), labelled as i-pixel. The i-pixel temperature measured by the FLIR camera accounts for the time-varying temperature relative to the beam and the background. The following equations were used to retrieve the beam oscillation:(1)c1Tb+c2Tt=Tm,c1+c2=1,
where, c1 and Tb are the fraction of the i-pixel occupied by the background and its temperature, respectively; c2 and Tt are the fraction of the i-pixel occupied by the beam and its temperature, respectively. Tm is the temperature measured by the camera for the i-pixel.The first equation of (1) represents the linear mixing of the temperatures of the beam and the background, the second is the constrain used for the total pixel fraction. The Tb and Tt temperatures, which are unknowns in (1), can be extracted from the temperature values of the neighbouring pixels of the i-pixel. The neighbouring pixels used for calculating Tb and Tt have to be pure pixels (i.e., made up of only beam or target).In (1), c2 represents the oscillation quantity expressed in terms of pixel fraction that varies from 0 to 1. Moreover, when c2 oscillation reaches values near to the maximum and minimum, this represents situations in which the amplitude of oscillation is close to or greater than the maximum measurable for a pixel, respectively. In the case that a pixel is not sufficient to encompass the oscillation, it is possible to solve the problem by considering, using the sum of the neighbouring pixels, “equivalent pixels” having a larger dimension compared to the original one.


### 2.3.2. Ground-Based Microwave Radar Interferometer

IBIS-S is based on two well-known radar techniques: the Stepped-Frequency Continuous Wave (SF-CW) and the Differential Interferometric [12, 18].The SF-CW technique allows the system to resolve the scenario in the range direction [15, 19], that is, to detect the position of target surfaces placed at different distances from the sensor. It is based on the synthesis and transmission of a burst of N monochromatic pulses equally and incrementally spaced in frequency (with fixed frequency step of Δf) leading to a work bandwidth B:(2)B=(N-1)Δf,
The N monochromatic pulses sample the scenario in the frequency domain similarly to a short pulse with a large bandwidth B. In a SF-CW radar, the signal source dwells at each frequency fx=f0+kΔf(k=0,1,2,…,N-1), the single tone signal is sufficiently long to allow the echoes backscattered by the target to reach the receiver. Hence, the duration of each monochromatic pulse (Tpulse) depends on the maximum distance (Rmax) to be observed in the scenario and is given by,(3)Tpulse≥2Rmax⁡c,
where c represents the velocity of the light.At each sampled time instant, bothI (In-phase) and Q (Quadrature) components of the received signals are acquired so that the resulting data consist of a vector of N complex samples, representing the frequency response measured at N discrete frequencies.By taking the Inverse Discrete Fourier Transform (IDFT) the radar response is reconstructed in the time domain: each complex sample in this domain represents the signal (echo) from a range (distance) interval of lengthcTpulse/2.The amplitude range profile of the radar echoes is then obtained by calculating the magnitude of each bin of the IDFT of acquired vector samples. This range profile gives a one-dimensional mapping of scattering objects in the “view space” as function of their relative distance from the sensor.The differential interferometric technique allows the system to measure the displacements of the structure illuminated by the antenna main beam by comparing the phase information of the backscattered electromagnetic waves collected at different times (see, e.g., [16]).Generally speaking, when a target surface moves with respect to the sensor module (emitting and back-receiving the electromagnetic wave) at least a phase shift arises between the signals reflected by the target surface at different times. Hence, the displacement of the investigated object is determined from the phase shift measured by the radar sensor at the discrete acquisition times. The radial displacementdp (i.e., the displacement along the direction of wave propagation) and the phase shift Δφ are linked by the following:(4)dp∝λ4πΔφ,
where λ is the wavelength of the probing electromagnetic signal.The radar emits a series of electromagnetic waves for the entire measurement period and processes phase information at regular time intervals to find any displacement occurring between one emission and the next one. It is worth underlining that the interferometric technique in (4) provides a measurement of the displacement along the line of sight for all the range bins of the structure illuminated by the antenna beam; once dp has been evaluated, the vertical displacement d can be easily found by making some geometric projections.IBIS-S system was applied in both Montagnole and Musmeci test beds. In the first case it was installed on a hill at a height of about 15 m from the ground level and at a distance of about 70 m from the target. In order to improve the accuracy and the signal-to-noise ratio of the measurements, a corner reflector was fixed on the beam. To acquire the measurements a range sampling frequency of 200 Hz and an acquisition time of about 60 s for each drop were used.In the Musmeci test bed IBIS-S system was placed just under the bridge deck with a perpendicular direction at a distance of about 30 m. A sampling frequency of about 100 Hz and an acquisition time of about 5 minutes were applied.In both cases the use of IBIS-S system allowed the dynamic monitoring of the investigated target without any contact between the sensor and the investigated target. This characteristic was very useful in Montagnole test bed where very restricted security conditions were required in order to carry out the experiment.


## 2.3.1. HFI Thermal Imagery

In the survey of measurements on the Montagnole test site, we took advantage of the thermal camera (LWIR) to obtain very high-frequency thermal imagery on the concrete beam during the direct impact of the falling blocks. The FLIR camera was mounted on a tripod at about 40 m (i.e., in box 3) and 70 m (i.e., near box 5) of distance from the beam. Box 3 was used for the first three drops (height of the falling ball of 1 m); for the drop 4 (height of the falling ball of 5 m), due to safety reasons, the camera was installed near the box 5 (70 m). For all the drops, the acquisition time was of 46 s for a total of 37000 frames.The thermal imagery acquisitions were used to calculate the beam oscillation frequency as derived from the difference between the beam temperature and the background temperature (see Figure1(b)).The method is based on the presence of mixed pixels composed by both the beam and the background as shown by the black arrow in Figure1(b), labelled as i-pixel. The i-pixel temperature measured by the FLIR camera accounts for the time-varying temperature relative to the beam and the background. The following equations were used to retrieve the beam oscillation:(1)c1Tb+c2Tt=Tm,c1+c2=1,
where, c1 and Tb are the fraction of the i-pixel occupied by the background and its temperature, respectively; c2 and Tt are the fraction of the i-pixel occupied by the beam and its temperature, respectively. Tm is the temperature measured by the camera for the i-pixel.The first equation of (1) represents the linear mixing of the temperatures of the beam and the background, the second is the constrain used for the total pixel fraction. The Tb and Tt temperatures, which are unknowns in (1), can be extracted from the temperature values of the neighbouring pixels of the i-pixel. The neighbouring pixels used for calculating Tb and Tt have to be pure pixels (i.e., made up of only beam or target).In (1), c2 represents the oscillation quantity expressed in terms of pixel fraction that varies from 0 to 1. Moreover, when c2 oscillation reaches values near to the maximum and minimum, this represents situations in which the amplitude of oscillation is close to or greater than the maximum measurable for a pixel, respectively. In the case that a pixel is not sufficient to encompass the oscillation, it is possible to solve the problem by considering, using the sum of the neighbouring pixels, “equivalent pixels” having a larger dimension compared to the original one.


## 2.3.2. Ground-Based Microwave Radar Interferometer

IBIS-S is based on two well-known radar techniques: the Stepped-Frequency Continuous Wave (SF-CW) and the Differential Interferometric [12, 18].The SF-CW technique allows the system to resolve the scenario in the range direction [15, 19], that is, to detect the position of target surfaces placed at different distances from the sensor. It is based on the synthesis and transmission of a burst of N monochromatic pulses equally and incrementally spaced in frequency (with fixed frequency step of Δf) leading to a work bandwidth B:(2)B=(N-1)Δf,
The N monochromatic pulses sample the scenario in the frequency domain similarly to a short pulse with a large bandwidth B. In a SF-CW radar, the signal source dwells at each frequency fx=f0+kΔf(k=0,1,2,…,N-1), the single tone signal is sufficiently long to allow the echoes backscattered by the target to reach the receiver. Hence, the duration of each monochromatic pulse (Tpulse) depends on the maximum distance (Rmax) to be observed in the scenario and is given by,(3)Tpulse≥2Rmax⁡c,
where c represents the velocity of the light.At each sampled time instant, bothI (In-phase) and Q (Quadrature) components of the received signals are acquired so that the resulting data consist of a vector of N complex samples, representing the frequency response measured at N discrete frequencies.By taking the Inverse Discrete Fourier Transform (IDFT) the radar response is reconstructed in the time domain: each complex sample in this domain represents the signal (echo) from a range (distance) interval of lengthcTpulse/2.The amplitude range profile of the radar echoes is then obtained by calculating the magnitude of each bin of the IDFT of acquired vector samples. This range profile gives a one-dimensional mapping of scattering objects in the “view space” as function of their relative distance from the sensor.The differential interferometric technique allows the system to measure the displacements of the structure illuminated by the antenna main beam by comparing the phase information of the backscattered electromagnetic waves collected at different times (see, e.g., [16]).Generally speaking, when a target surface moves with respect to the sensor module (emitting and back-receiving the electromagnetic wave) at least a phase shift arises between the signals reflected by the target surface at different times. Hence, the displacement of the investigated object is determined from the phase shift measured by the radar sensor at the discrete acquisition times. The radial displacementdp (i.e., the displacement along the direction of wave propagation) and the phase shift Δφ are linked by the following:(4)dp∝λ4πΔφ,
where λ is the wavelength of the probing electromagnetic signal.The radar emits a series of electromagnetic waves for the entire measurement period and processes phase information at regular time intervals to find any displacement occurring between one emission and the next one. It is worth underlining that the interferometric technique in (4) provides a measurement of the displacement along the line of sight for all the range bins of the structure illuminated by the antenna beam; once dp has been evaluated, the vertical displacement d can be easily found by making some geometric projections.IBIS-S system was applied in both Montagnole and Musmeci test beds. In the first case it was installed on a hill at a height of about 15 m from the ground level and at a distance of about 70 m from the target. In order to improve the accuracy and the signal-to-noise ratio of the measurements, a corner reflector was fixed on the beam. To acquire the measurements a range sampling frequency of 200 Hz and an acquisition time of about 60 s for each drop were used.In the Musmeci test bed IBIS-S system was placed just under the bridge deck with a perpendicular direction at a distance of about 30 m. A sampling frequency of about 100 Hz and an acquisition time of about 5 minutes were applied.In both cases the use of IBIS-S system allowed the dynamic monitoring of the investigated target without any contact between the sensor and the investigated target. This characteristic was very useful in Montagnole test bed where very restricted security conditions were required in order to carry out the experiment.


## 3. Results and Discussion

### 3.1. HFI Results

For Montagnole experiment, the drops of the heavy iron ball on the beam were used as mechanical solicitation to measure the dynamic response of the beam. For this study we used a simple model by considering the behaviour of the beam as the one of a damped harmonic oscillator. The equation describing the oscillation is(5)z(t)=A0e-ςω0tsin(1-ςω0t+φ),
where A0 is the amplitude and φ the phase that dictates the initial conditions, ω0 is the undamped angular frequency of the oscillator and ζ is damping ratio that are given, respectively, by:(6)ω0=km,ς=c2mω0,
where k is the elastic constant, m is the mass, and c is the viscous damping coefficient.The value ofζ determines the behaviour of the system in fact, the damped harmonic oscillator can be (a) for ζ>1, the system is represented by exponentially decays; (b) for ζ=1, the system returns to equilibrium as quickly as possible without oscillating; (c) for ζ<1, the system is represented by an underdamped harmonic oscillator (Figure 2).Figure 2
The graph shows the under-dumping decay attained for the fraction of thei-pixel occupied by the beam (i.e., blue line C2) used for the Montagnole experiment. Red line depicts the theoretical dumper oscillator calculated from the parameters of Table 1.The quantity that describes the beam oscillation is thec2 value represented in (1). This value was calculated for all the frames of each acquisition performed on the beam. Figure 2 depicts the comparison between the quantity c2 (blue line) calculated for the third drop and the theoretical underdamped harmonic oscillator (red line) that approximates the c2 oscillations. The parameters used for evaluating the underdamped harmonic oscillation, that is, those that better approximates the c2 oscillation, are reported in Table 1.Table 1
Parameters used for calculating the underdamped harmonic oscillator from (5).
A00.700φ1.1ω028.049ζ0.021Freq.4.5 HzMoreover, in order to obtain the resonance frequencies of the beam we applied the Fast Fourier Transform (FFT) on the quantityc2 evaluated for the four drops. Figure 3 depicts the time behaviour of the c2 for the four drops, whereas Figure 4 reports the amplitude of the FFT of c2. Table 2 shows the oscillation frequencies of the two most significant FFT coefficients.Table 2
Significant oscillation frequencies of the beam under the four direct impacts.
DropFirst frequency (Hz)Secondary frequency (Hz)14.79.324.513.734.513.244.68.4The graphs show the underdamped decay attained for the fraction of thei-pixel occupied by the beam used for the four Montagnole drops.
(a)(b)(c)(d)FFT absolute coefficients relative to thec2 oscillations.
(a)(b)(c)(d)


### 3.2. Ground-Based Microwave Radar Interferometer Results

IBIS-S data acquired in Montagnole test site were processed by using the commercial IBIS DATA VIEWER (IBISDV) software. The software works by using MATLAB libraries and is based on (i) Inverse Discrete Fourier Transform (IDFT), and (ii) Differential Interferometric technique.The IBIS-S results allowed the measure of the dynamic response in time and frequency domain of the beam affected by the direct impact of the heavy iron ball.Figures5 and 6 show the response in time and frequency domain of the beam oscillation, due to the direct impact of the heavy iron ball, for the four drops, respectively.The graphs show the time domain response of the beam after each drop. A time window of about 20 s was considered.
(a)(b)(c)(d)The graphs report the frequency domain response of the beam during the four drops, a frequency range up to 10 Hz is considered.
(a)(b)(c)(d)In the time domain a window of about 20 sec after each direct impact was considered in order to study the dynamic behaviour of the beam. The results show a negative line of sight (LOS) permanent displacement for each drop. In particular, the displacement increases from the first to the last drop as reported in Table3.Table 3
LOS permanent displacements of the beam after each drop.
Drop numberPermanent displacement (mm)1−0.52−2.53−3.04−10.0The different LOS permanent displacement could be due to a combination between horizontal and vertical movement of the beam towards the IBIS-S sensor after each drop. This displacement is less than 3 mm and grows slightly for the first three drops where the falling block conditions (drop height equal to 1 m) are the same. The permanent displacement is amplified after the fourth drop where the iron ball falls from 5 m height.Before applying the Fourier Transform to the signals of Figure5, the mean value and the linear trend were removed from the signal in the selected window. The amplitude of the spectra in a frequency range up to 10 Hz is reported in Figure 6. For all the drops, two main resonant peak frequencies arise, the first one at about 1.5 Hz and the second one at about 4.5 Hz. As can be argued from Table 4, both the main resonant peak frequencies decrease after any direct impact of the iron ball on the beam.Table 4
Oscillation frequencies of the beam under the four direct impacts measured by Ground-based microwave radar interferometer.
DropFirst frequency (Hz)Secondary frequency (Hz)11.64.821.64.631.44.541.43.8


## 3.1. HFI Results

For Montagnole experiment, the drops of the heavy iron ball on the beam were used as mechanical solicitation to measure the dynamic response of the beam. For this study we used a simple model by considering the behaviour of the beam as the one of a damped harmonic oscillator. The equation describing the oscillation is(5)z(t)=A0e-ςω0tsin(1-ςω0t+φ),
where A0 is the amplitude and φ the phase that dictates the initial conditions, ω0 is the undamped angular frequency of the oscillator and ζ is damping ratio that are given, respectively, by:(6)ω0=km,ς=c2mω0,
where k is the elastic constant, m is the mass, and c is the viscous damping coefficient.The value ofζ determines the behaviour of the system in fact, the damped harmonic oscillator can be (a) for ζ>1, the system is represented by exponentially decays; (b) for ζ=1, the system returns to equilibrium as quickly as possible without oscillating; (c) for ζ<1, the system is represented by an underdamped harmonic oscillator (Figure 2).Figure 2
The graph shows the under-dumping decay attained for the fraction of thei-pixel occupied by the beam (i.e., blue line C2) used for the Montagnole experiment. Red line depicts the theoretical dumper oscillator calculated from the parameters of Table 1.The quantity that describes the beam oscillation is thec2 value represented in (1). This value was calculated for all the frames of each acquisition performed on the beam. Figure 2 depicts the comparison between the quantity c2 (blue line) calculated for the third drop and the theoretical underdamped harmonic oscillator (red line) that approximates the c2 oscillations. The parameters used for evaluating the underdamped harmonic oscillation, that is, those that better approximates the c2 oscillation, are reported in Table 1.Table 1
Parameters used for calculating the underdamped harmonic oscillator from (5).
A00.700φ1.1ω028.049ζ0.021Freq.4.5 HzMoreover, in order to obtain the resonance frequencies of the beam we applied the Fast Fourier Transform (FFT) on the quantityc2 evaluated for the four drops. Figure 3 depicts the time behaviour of the c2 for the four drops, whereas Figure 4 reports the amplitude of the FFT of c2. Table 2 shows the oscillation frequencies of the two most significant FFT coefficients.Table 2
Significant oscillation frequencies of the beam under the four direct impacts.
DropFirst frequency (Hz)Secondary frequency (Hz)14.79.324.513.734.513.244.68.4The graphs show the underdamped decay attained for the fraction of thei-pixel occupied by the beam used for the four Montagnole drops.
(a)(b)(c)(d)FFT absolute coefficients relative to thec2 oscillations.
(a)(b)(c)(d)


## 3.2. Ground-Based Microwave Radar Interferometer Results

IBIS-S data acquired in Montagnole test site were processed by using the commercial IBIS DATA VIEWER (IBISDV) software. The software works by using MATLAB libraries and is based on (i) Inverse Discrete Fourier Transform (IDFT), and (ii) Differential Interferometric technique.The IBIS-S results allowed the measure of the dynamic response in time and frequency domain of the beam affected by the direct impact of the heavy iron ball.Figures5 and 6 show the response in time and frequency domain of the beam oscillation, due to the direct impact of the heavy iron ball, for the four drops, respectively.The graphs show the time domain response of the beam after each drop. A time window of about 20 s was considered.
(a)(b)(c)(d)The graphs report the frequency domain response of the beam during the four drops, a frequency range up to 10 Hz is considered.
(a)(b)(c)(d)In the time domain a window of about 20 sec after each direct impact was considered in order to study the dynamic behaviour of the beam. The results show a negative line of sight (LOS) permanent displacement for each drop. In particular, the displacement increases from the first to the last drop as reported in Table3.Table 3
LOS permanent displacements of the beam after each drop.
Drop numberPermanent displacement (mm)1−0.52−2.53−3.04−10.0The different LOS permanent displacement could be due to a combination between horizontal and vertical movement of the beam towards the IBIS-S sensor after each drop. This displacement is less than 3 mm and grows slightly for the first three drops where the falling block conditions (drop height equal to 1 m) are the same. The permanent displacement is amplified after the fourth drop where the iron ball falls from 5 m height.Before applying the Fourier Transform to the signals of Figure5, the mean value and the linear trend were removed from the signal in the selected window. The amplitude of the spectra in a frequency range up to 10 Hz is reported in Figure 6. For all the drops, two main resonant peak frequencies arise, the first one at about 1.5 Hz and the second one at about 4.5 Hz. As can be argued from Table 4, both the main resonant peak frequencies decrease after any direct impact of the iron ball on the beam.Table 4
Oscillation frequencies of the beam under the four direct impacts measured by Ground-based microwave radar interferometer.
DropFirst frequency (Hz)Secondary frequency (Hz)11.64.821.64.631.44.541.43.8


## 4. Conclusions

The proposed multitechniques approach, which is the main aim of the ISTIMES Project, can be useful to individuate suitable parameters which will allow a continuous and/or rapid remote control of the condition of infrastructures to plan appropriate interventions for their conservation. This is even more important when analyzing complex structures where the use of a single diagnostic tool can provide results leading to ambiguous interpretations.However, it is difficult to have at disposal a complete setup of instruments measuring the same parameter, for example, as in this case where HFI thermal imagery and ground-based (G-B) microwave radar were applied to study the dynamic behaviour of a beam affected by changeable strains.G-B microwave radar was applied with the aim of measuring the dynamic displacement of the beam. This technique provided good results, giving information about the behaviour of the target investigated in time and frequency domains. In particular, it allowed the estimation of the displacement affecting the beam after each drops and the measurement of the main frequencies characterizing the structure. Two main peaks at different frequencies were highlighted; both the frequencies decrease after the two last drops due to the increase of the number of severe cracks observed on the beam. The cracks damaging the beam have induced a loss of stiffness with a consequent raise of the period of the eigen modes and a decrease of the eigen-frequencies.HFI thermal imagery camera was applied with the aim of measuring the structure oscillations by high frequency analysis. Also in this case, the results allowed to discriminate two main peaks at different frequencies after each drop. The first peak measured by HFI at around 4.5 Hz is comparable with the second peak measured by G-B radar. However, HFI frequencies do not highlight a decrease after the last drop as well as the G-B results (see Table5). This could be explained considering the different line of view of the two techniques; indeed, G-B microwave radar was placed at a different height respect to the HFI camera and it observed the beam by a specific angle. HFI thermal camera was placed at the same height of the beam and it observed the target perpendicularly. Moreover, the information provided by G-B microwave radar is related to displacement of the beam along the line of sight, while the information from HFI camera is related to the displacement of the beam in the plane perpendicular to the sensor. The comparison between HFI and G-B results also suggests that G-B radar better points out the low peak frequencies, whereas the HFI thermal camera highlightes the high frequencies.Table 5
Comparison between the oscillation frequencies measured by HFI and GB microwave Radar with the relative percentage difference.
DropHFI (HZ)GB microwave radar (Hz)Relative percentage difference14.74.82.124.54.62.234.54.50.044.63.817.4About the lowest frequencies retrieved by G-B radar (see Table4), two possible explanation could be given as: a vibration mode of the beam or the main oscillation frequency of dropping machinery after each drop.Future work will be focused on (a) laboratory and field measurements to better understand the different results attained by the HFI and GB microwave Radar (b) a numerical modelling of the beam to estimate its main frequencies, and (c) the development and analysis of the integration of the data acquired by different sensor technologies on infrastructures in order to evaluate their effectiveness and suitability for this kind of application.



---

*Source: 204976-2011-08-11.xml*