Aging of buildings during their service life has attracted the attention of researchers on structural health monitoring (SHM). This paper is related with detecting damage in building structures at the earliest possible stage during seismic activity to facilitate decision-making on evacuation before physical inspection is possible. For this, a simple method for damage assessment is introduced to identify the damage story of multistory buildings from acceleration measurements under a wave propagation approach. In this work, damage is assumed as reduction in shear wave velocities and changes in damping ratios that are directly related with stiffness loss. Most damage detection methods are off-line processes; this is not the case with this method. First, a real-time identification system is introduced to estimate the current parameters to be compared with nominal values to detect any changes in the characteristics that may indicate damage in the building. In addition, this identification system is robust to constant disturbances and measurement noise. The time needed to complete parameter identification is shorter compared to the typically wave method, and the damage assessment can keep up with the data flow in real time. Finally, using a robust threshold, postprocess of the compared signal is performed to find the location of the possible damage. The performance of the proposed method is demonstrated through experiments on a reduced-scale five-story building, showing the ability of the proposed method to improve early stage structural health monitoring.

1. Introduction

Structures like buildings are frequently exposed to various kinds of damages during their service life mainly due to the continuous deterioration caused by material degradation, corrosion, fatigue, and unexpected heavy loads induced by earthquakes and strong winds. In this sense, structural health monitoring (SHM) allows to track changes in the dynamic characteristics of the structural system for inspecting structure’s integrity and detecting any probable structural damage. Most SHM methods in civil engineering are based on modal analysis, which studies global changes in natural frequencies and modal shapes [1, 2]. Under this idea, extensive reviews of vibration-based damage identification methods are presented in [35] and subsequently every 5 years since 2003. Recent reviews on vibration-based damage detection can be consulted in [6, 7]. A limitation of this group of methods is that they require exciting a building at high frequencies, condition that is not easy to achieve, and therefore, damage may go unnoticed. In the same research direction, an improved method to detect the location of the damage is the curvature mode shape proposed in [8, 9] that is more sensitive to stiffness loss and strain change based on a flexibility index [10]. Damage can be localized by comparing changes in displacement mode shapes. However, if the damage is distributed throughout the structure and a baseline measurement data set of the undamaged structure is not available, the algorithm performance may decrease as a damage indicator.

Other studies for structural health monitoring are developed around a pattern classification framework; a review of these methods is presented in [11]. For instance, Sohn and Farrar and Nair and Kiremidjian [12, 13] provide a procedure for damage detection and localization in a mechanical system analyzing the measured structural-vibration response with autoregressive models (ARX) and autoregressive moving average (ARMA), respectively. The assumption is that when the structure is damaged, the error between a measured signal and that obtained from these two models increases; thus, changes in the coefficients of the autoregressive part of the model can indicate damage. Following this line, convolutional neural network is an algorithm that allows to realize damage diagnosis using image classification [14, 15]. Other applications can be found in [16] that are based on a machine learning perspective to detect the damage of a bridge pier. Like in [17], damage detection is carried out by using a support vector machine (SVM), and Grubbs’ test method is adopted to delete the outliers for the SVM. An artificial bee colony (ABC) algorithm via modal and frequency analyses is presented in [18] to detect structural damages under different identify conditions such as grouping, noise effect, mode shape order, and sensor location and is also evaluated by several explicit test functions. Damage detection is treated as an optimization problem. Moreover, note that all previous pattern recognition methods demand large time histories from the undamaged structure and intensive data processing.

Some new techniques employ model-based methods assuming that a detailed well-correlated numerical model of the structure is available for damage identification. Among them, finite element (FE) model is typically employed to minimize the differences between analytical and experimental results and is usually proposed as an optimization problem that minimizes the discrepancy by adjusting the unknowns of the FE mode, employing finite element model updating (FEMU). An extensive bulk of research studies of structural dynamic model updating techniques based on direct and iterative techniques of FEMU can be found in [19, 20]. The authors highlight the current issues, applications, and observations for further advancements in the field. Following this line, Sehgal and Kumar [21] presented the theory and numerical application of a damped updating technique, which is based on response models and Derringer’s function approach. Updating process is formulated as an optimization problem, wherein desirability functions are formulated, based on natural frequencies and modal assurance criterion (MAC) values. Similarly, in [22], the application of Derringer’s function-based weighted model updating method for damage detection in a simulated cantilever beam structure is investigated. In this method, FEMU is treated as a multiobjective optimization problem, where the number of objectives needs to be defined in such a way as to reduce errors in responses predicted by the FE model. Maximization of an individual desirability function leads to minimization of error in the corresponding finite element response to which the individual desirability function is attached. Experimental validation of this technique can be found in [23], where structural dynamic model updating has been utilized in experimentally evaluating the extent of damage at six locations of a damaged cantilever beam structure. Accurate damage identification was successfully performed by evaluating percent reduction in flexural rigidity. Similarly, a damped updating method to identify damage as well as damping parameters in a structure under real conditions is presented in [24]. Novelty of the work corresponds to experimental validation of the desirability function-based multiobjective optimization approach to identify the structural damage as well as damping parameters accurately in the real-life structure. It is important to note that the success of finite element model updating depends heavily on the selection of updating parameters since an ill-conditioned numerical problem may produce an updated model that becomes unsatisfactory or unrealistic [25].

Recently, the wave propagation approach has become relevant as an alternative method to identify building local characteristics. The underlying hypothesis is that, after damage occurred, shear wave velocity of propagation within the damaged zone is reduced. In the same research direction, [2628] presented different wave methods for structural health monitoring and discussed advantages and limitations of this approach. An extension of studies on pulse propagation in uniform and layered shear beams is presented in [29, 30], where system identification of a layered shear building model is developed by using a direct algorithm based on ray theory. Similarly, Rahmani and Todorovska [31, 32] presented system identification with a waveform inversion algorithm for characterizing the shear velocities in layered shear beam and torsional shaft models. Seismic records from a 54-story steel-frame building are analyzed, and from them, wave velocity is derived as a damage indicator. In [33, 34], a method is proposed to monitor changes in the velocity of wave propagation through the structure identified by a least squares fit of a Timoshenko beam model. The author reported that the wave method could be successfully applied to a damaged structure with very dispersed wave propagation, resulting from bending deformation, in the absence of significant soil-structure interaction effects. From report results, wave propagation methods are promising in practice due to their robustness when applied to real damaged structures and insensitivity to the effects of soil-structure interaction [35, 36]. One common characteristic of all the above methods under this approach is that the analysis is performed off-line. This allows data processing or filtering that is normally computationally intensive; however, it prevents real-time response to changes in structures’ health condition.

This paper extends the results presented in [3032], where the off-line estimation of wave velocities in a layered shear beam model is studied. Here, two important features are included: (i) a new parameterization of the identification system problem to estimate velocities and damping coefficient and (ii) a methodology to locate possible damage, which is based on comparing changes in the wave parameters, processing these changes, and assessing damage using appropriate thresholds. Although the wave propagation approach can be applied at large and small scales, in this paper, an intermediate scale is proposed that focuses on the damage location at a story level. The method is applied to structures built with frames under the action of shear forces in one direction. The developed identification algorithm is intended for use to SHM systems to provide early alerts during seismic activity. Therefore, if a structure is not safe, this system identification facilitates decision-making on evacuation before physical inspection is possible. Such timely damage assessment will help to avoid endangering human lives and to prevent injuries caused by potential collapse of a weakened structure during shaking from aftershocks, as well as unnecessary evacuation and loss of function of important structures such as hospitals when damage is nonexistent. Moreover, real-time damage detection allows fast response, avoids data storage, and discards the use of off-line data processing for eliminating acceleration measurement bias. This last one is achieved by using linear integral filters (LIF), as in [37], which eliminate constant disturbances, attenuate measurement noise, and avoid wave dispersion. This last feature is convenient as there is evidence of dispersion in structural systems [32, 38]. The system identification uses the normalized least square algorithm (NLSA) with forgetting factor, which provides real-time estimation of current parameters in the discretized shear beam model in contrast with [33, 39]. Therefore, the contribution of this paper can be summarized as follows:(i)Proposing an online identification system that allows to recover in real time the best possible set of parameters in the least square sense that minimizes the error between the real and predicted output. This approach avoids dealing directly with wave reflection, transmission, and dispersion, as in previous wave methods. Moreover, the parametric projection algorithm guarantees estimates with a physical sense.(ii)Introducing a new parameterization for the identification system that improves estimations of wave parameters considering measurements from a specific story and its adjacent floors. This improves the precision of the estimated parameters compared with previous works. Moreover, the discrete shear beam model only requires two free parameters per floor in contrast with the discrete shear building model.(iii)Reducing the time needed to complete the parameter identification in comparison with a typical wave propagation method, and allowing the damage assessment to keep up with the data flow in real time.(iv)Providing a real-time damage detection algorithm that is easier to implement in comparison with other alternative methodologies already available in the literature. When there is damage, ratios between parameters of shear wave velocity and damping in two adjoining floors always change at the damaged story. When there is no damage, these ratios would not change with time. Moreover, the selected damage detection threshold is robust at certain close values.

The specimen used in this study is a five-story small-scale building, where first, shear wave velocities and damping coefficients are obtained when the structure is damage-free. Later, two different scenarios are presented. In the first one of them, damage is produced at the second story and in the second one, in the second and fifth floors. In both cases, damage is performed by reducing the cross section of the columns, which results in stiffness loss. The experimental setup allows to easily calculate both analytically and experimentally the new expected stiffness values, as this is a useful parameter to evaluate the estimated shear wave velocities and damping coefficient, because these last parameters cannot be experimentally measured. However, from the results published in [30], it is known that wave propagation methods are sensitive to a wide type of damages. Moreover, considering nominal and currently estimated parameters, a methodology to locate the damage is introduced. The most important challenge is in defining the threshold range for damage detection, which was determined using the data collected in the tests to find the maximum change level. This allows us to calculate the standard deviation that each parameter may suffer.

This paper is organized as follows: the employed mathematical model is briefly presented in Section 2, and the real-time system identification design is presented in Section 3. The damage location methodology is described in Section 4, while experimental results are presented in Section 5. Finally, concluding remarks are presented in Section 6.

2. Mathematical Model of the Building

The dynamics of a multidegree-of-freedom (MDOF) building structure seismically excited is modeled employing a one-dimensional wave equation with Kelvin damping to study the wave propagation characteristics. The model and proper boundary conditions are [37]where β is the shear wave velocity, η is a constant proportional to damping in the structure, and is the shear module, whereas and are signals of displacement and velocity, and is the ground motion. Furthermore, according to equations (2a2d), the following is assumed.

2.1. Remarks

(1)The building is initially at rest according to equations 2a and 2b(2)Only the base of the structure is excited during seismic activity, see equation (2c)(3)The building is stress-free at the roof in accordance with equation (2d)(4)From equation (1), the shear building model can be treated as an elastic continuous shear beam, and thus, wave propagation effects can be studied and used as indicators of structural health

In order to propose a real-time identification system, the continuous shear beam model is discretized in points that coincide with instrumented stories. Therefore, different shear wave velocities and damping coefficients are produced, as shown in Figure 1. Discretization is carried out employing three-point central-difference approximation for stories 1 to (n − 1) and a first-order backward approximation for the remaining ones. Moreover, in order to improve discretization accuracy, the shear wave velocity and damping coefficient at each point are obtained with the average values of the nearest points, i.e.,which will induce a different parameterization than the one presented in [38].

Furthermore, the resulting state matrix is expanded to use the acceleration measurement as the excitation signal, which is applied only at the Dirichlet boundary condition, in contrast to shear building models, where excitation is applied to all stories at the same time [40]. Eventually, the new ordinary differential and algebraic system equations have the possibility to be solved on real time [41]. Therefore, wave equation (1) is approximated bywherewith ; is the total number of floors. Moreover, definingthen and are rewritten aswithwith

However, equation (12) can be rewritten as , i.e.,

Similarly, .

On the contrary, , where is the building height. Matrices and are the result of discretization, and they are related to the number of floors of the building. Considering that , the acceleration measurement at the boundary condition is available and equivalent to seismic action. Similarly, the displacement and velocity of the basement are denoted by , and , respectively, while the remaining , and with are the absolute displacements, velocities, and accelerations of each floor. Therefore, equation (5) does not need a coordinate transformation to be solved compared with other works [42]. Additionally, defining allows expressing equation (5) in the space form aswithwhose output equation is given bywhere and denote, respectively, the null and identity matrices of size .

3. System Identification

Defining and , equation (18) is rearranged in a suitable form for identification:where elements of diagonal matrices and are used as elements in vectors and ; matrices and vectors , are defined as

Hence, equation (19) is now of the form

However, equation (21) still depends on the nonmeasurable position and velocity signal contained in vectors and because displacement and velocity cannot be directly measured during seismic excitation. To overcome this problem, the following procedure is introduced. First, the model (21) is expressed in the frequency domain using the Laplace transform:where . Now, deriving three times equation (22), constant disturbances are removed, producingwhere is given by

Note that while equation (22) depends on nonmeasurable signals, equation (24) depends only on measurable acceleration signals. Furthermore, (24) can be written in the time domain aswhere the superscript represents the j-th time derivative of the corresponding variable. Moreover, considering that acceleration signals are contaminated with offset and measurement noise, linear filter integrals (LIF) are introduced to eliminate this constant disturbance and to attenuate measurement noise [42]. This technique overcomes limitations that affect the performance of numerical integrators, as it is the case with bias and amplification of low-frequency noise. LIF is defined by the operatorwhere is the number of integrations over finite time intervals of signal . The time interval, , is defined by the size of the moving window aswhere is the sampling period and determines the length of the time integration window, which must be properly selected by making the system bandwidth and the filter cutoff frequencies to coincide. Applying the LIF operator to both sides of equation (26), which corresponds to the fifth-order low-pass filter, yields terms depending on the measured variables and , i.e.,


Thus, , and , is the binomial coefficient, and finally, is a degree of freedom to normalize the amplitude of the filtered signals. From equation (29), the system output is parameterized aswhere is a vector containing the real parameters of the system that will be estimated with the normalized least square algorithm (NLSA) and is the regressor matrix formed by filtered acceleration signals.

4. Damage Location Methodology

In this section, a two-stage damage location methodology is described. First, the parametric estimation algorithm developed is presented; then, once the current parameters that characterize the damage behavior are recovered, they are compared with nominal parameters, and a damage location algorithm (DLA) is developed.

4.1. Parameter Estimation

Let be the vector containing estimated parameters in equation (33) which are identified by using the normalized least square algorithm (NLSA) with forgetting factor given by equations (34 and 35) [43]:where is the forgetting factor, is the covariance matrix, with , , , and , which satisfies and assures that the normalized estimated error tends to zero, i.e.,

Moreover, NLSA is improved using a parameter projection scheme that ensures only positive estimates and is defined bywhere and are the projected estimates of and , respectively, with , and are the lower bounds for estimated parameters, whereas and allow to improve the estimation, and finally, . Note that produce an approximation of the sign function, e.g.,which allows to assume in equation (37) that when ; moreover, when . Similar conditions apply to , ensuring that always lie within , as shown in Figure 2. Therefore,

Moreover, projected estimates in (37 and 38) have the advantage to be continuously differentiable functions, a useful feature for assuring that the estimation is stable (see Appendix A).

4.2. Damage Location Algorithm

Assuming that the reference parameters that characterize a nondamage building have been previously estimated or calculated from experiments or material properties, the system identification proposed in the above sections is used to identify real-time parameters that may indicate damage. These parameters are then compared to the nominal parameters in order to locate structural damage. However, an important question in this procedure is how to define an appropriate threshold for damage detection? For this, next, some considerations that must be considered in its design are enumerated:(a)Taking into account that, for several experiments under the same conditions, results may vary, it is necessary to obtain the average of all involved parameters to obtain a nominal model.(b)Because structural damage causes changes in all structural parameters, it is necessary to find the maximum change level. This allows to calculate the standard deviation that each parameter may suffer.(c)Since all the parameters change, the ratio of parameters corresponding to 2 contiguous floors cancels this effect and only occurs when the rate of change is greater than the previously calculated standard deviation.(d)Finally, due to the damping, which exhibits a random behavior and generally inverses to the wave propagation velocity, the weighted sum of both parameters defines the damage detection threshold.

Therefore, the developed procedure is as follows:

For ,

The relationships between pairs of nominal parameters are obtained.

The relationships between pairs of identified parameters that could indicate damage are also calculated:

For ,

The relative errors between the reference parameters and those possibly indicating damage are computed, i.e.,

The standard deviation of the relative errors is determined bywhere and are the mean values of and , respectively. Moreover, if the relative error , ; otherwise, the previous difference is considered null. Similar conditions are applied to values of .

Finally, damage is localized using the following expression:

Note that when , a value of is assigned to prevent indetermination in equation (47). Figure 3 shows the damage detection and localization scheme developed here. This algorithm has the advantage to operate in real time, and therefore, the current condition of structures can be estimated during a seismic action, and an early warning can be issued. Moreover, the damage assessments can be stored and compared, and this information can be useful to decide whether the building should be evacuated, repaired, demolished, or still kept in use.

5. Experimental Results

A five-story building prototype, as shown in Figure 4(a), isused to evaluate the performance of the proposed damage location algorithm. The building is made from aluminum, with the exception of three columns that are made from brass. The dimensions of the columns of the first and the remaining fourth floors are 0.635 × 0.635 × 31.5 cm and 0.635 × 0.635 × 36 cm, respectively. Considering that the first story has an attachment to the base with an height of 4.5 cm, the structure has dimensions of 60 × 50 × 180 cm, and it is mounted on a shaking table actuated by servomotors from Parker, model 406T03LXR, which can be moved in the x-axis and y-axis of the horizontal plane, respectively. However, experiments are based only on the y-axis motion. In addition, the structure is equipped with accelerometers based on MEMS Technologies, model ADXL203E, with a range of ±1.7 g and a bandwidth of [0–50] Hz, which allow to measure absolute acceleration at every story and ground level. Data acquisition was taken using two PCI-6221-M series, electronic boards from National Instruments. Communication between these boards and Simulink were carried out using the Matlab Real-Time Windows Target toolbox. Experiments were executed using a sampling time of 0.001 s.

The North-South component of the 1985 México City earthquake, measured at Secretaria de Comunicaciones y Transportes (SCT), is employed as a source of seismic excitation for the experimental structure. Taking into account the reduced-scale of the building, the amplitude of excitation is adjusted to match with the structure maximum motion, as depicted in Figure 5(a). It is important to note that the used signal has enough richness to identify the structural parameters.

To have preliminary information about the building characteristic, the structure is excited using a chirp signal with a swept-frequency from 0.1 to 15 Hz. The resulting building frequencies are  Hz,  Hz,  Hz,  Hz, and Hz, as depicted in Figure 5(b). Considering that the largest frequency is 11.51 Hz, the linear integral filter is designed with a cutoff frequency of 25 Hz to avoid modifying the structural response. Hence, the integration period is  s, and the integration window length results in . The NLSA uses the following parameters; and . For parameter estimation, the building model is discretized at 5 points that match the stories, where acceleration measurement is recorded.

The seismic excitation is applied only at the building-base level, and shear wave velocities and Kelvin damping coefficients are identified for every floor that characterize the nondamage building, as presented in Table 1. From this table, it is clear that the algorithm estimates similar structural properties for floors 2–5, while the first floor differs because its height is slightly smaller than that of other floors, increasing its wave velocity. Evolution of identified parameters is shown in Figures 6(a) and 6(b). However, as real wave velocities and damping coefficients are unknown, they cannot be compared with their estimates. The parameters are validated by comparing the stiffness from the estimated parameters with that obtained from material properties. To carry out this procedure, it is useful to define the following stiffness and damping equations:

Equations (48) and (49) allow to compute lateral stiffness and damping in the structure. Note that is a diagonal matrix, where each one of its elements is the mass of the respective story. For the experiments, it was possible to measure each one of these masses that resulted in  kg and  kg. Table 2 presents the identified lateral stiffness compared with that calculated directly using nominal values of mechanical properties [44], showing that the estimated stiffness is sufficiently close to analytical values on floors 2–5. The first floor is affected by the clamping mechanism increasing its stiffness. In consequence, an error of approximately 5.5% is generated between analytical and estimated stiffness in the worst case. For a real-time experiment, results are considered satisfactory. Damping ratio is more difficult to compare because reference values are not precise [40].

The parameters obtained in this section will be used as reference parameters in the next sections for the purpose of damage detection. For convenience, parameters and are renamed as and , while and are renamed as and , respectively, where subscript r refers to nominal or reference values.

5.1. Damage Location Task

In this section, the problem of damage detection and location in the five-story reduced-scale building described in the previous section is studied. Here, the building is affected by the presence of damage with two different damage settings: (1) damage at the second story and (2) damage at both second and fifth stories, as in Figures 4(b) and 4(c), respectively. Experiments consist of 2 total tests. In both cases, induced damage was through stiffness loss, and the square cross section of these columns was reduced from of 6.35 to 5 mm, i.e., approximately 23% less with respect to nominal conditions. The intention of this setup controlled experiments is only to ensure a good and appropriate evaluation of the proposed algorithm, with respect to analytical values. Nevertheless, methods under the wave propagation approach are robust to the type and amount of damage, as discussed in [30]. It is important to note that, during experiments, brass and aluminum columns were numbered and referenced at specific positions within the building in order to assure that, after machining, they were replaced in their original position. Moreover, to rebuild the experimental prototype, the torque exerted on the union screws was calibrated to 3 N/m. This helped to guarantee similar conditions for all experiments. Moreover, the 1985 earthquake, as shown in Figure 5(a), was also used as the seismic excitation for experiments.

Note that the damage produces a decrease in the shear wave velocity and Kelvin damping parameters and with respect to reference values, as shown in Figures 7(a), 8(a), 7(b), and 8(b), respectively, where subscript d is used to describe parameters obtained with damage presence. Details of these can be found in Tables 3 and 4.

Moreover, vibration frequencies are calculated from estimated parameter matrix . A comparison between nominal and estimated bandwidths after damage can be found in Table 5. From the results in Table 5, it is clear that all vibration frequencies changed and that the structural response bandwidth is reduced, which is to be expected because damage produces a different structural system in accordance with [4, 5]. This confirms that shear wave velocity and Kelvin damping are successful damage indicators since they are sensitive to structural changes and intrinsically contain information about the damage. On the contrary, considering that the mass-loss in columns is small with respect to the mass of floors, it is valid to assume that the masses here are the same as those of the experiment one without damage, i.e.,  kg and  kg. Therefore, the stiffness and damping values are recalculated using equations (48 and 49). Table 6 presents the results in more detail and their comparison with the nominal values.

Despite the information provided by the algorithm, shear wave velocities, Kelvin damping, vibration frequencies, and stiffness values, damage location is still a hard task to realize. For instance, observing results in Table 3, it is evident that the parameters corresponding to the first floor are the most affected because they present a difference of 30.8% between nominal wave propagation velocities and those that identify a damaged behavior. However, considering that the rest of the parameters also changed, it is difficult to assure that damage is located on the first story. For example, the second floor presents a difference of 29% between parameters corresponding to the same level, while on the third floor, there is a difference of 23% for the 2 different damage settings. A similar situation happens analyzing the vibration frequencies and stiffness values presented in Tables 5 and 6, respectively. In contrast to these parameters, the Kelvin damping coefficients presented in Table 4 vary in a random way, and this even more complicates damage location. To overcome these problems, the algorithm proposed in Section 4.2 uses a threshold for the damage location. The simplicity of the method makes it quite promising for real-time applications since no database is needed. Moreover, in order to evaluate the performance of the damage location algorithm, the method developed in Section 4.2 is applied; in addition, a sensitivity analysis is carried out by changing the detection threshold, in particular varying the standard deviation from 70% to 110% of the value calculated in equation (46). Figures 9(a) and 9(b) confirm that the structural damage occurs on floor 2 and then on both floors 2 and 5, for configurations 1 and 2, respectively. The results are still satisfactory varying the standard deviation from 70% to 110%, indicating that the selected damage detection threshold is robust to the threshold selection. All the previous results were obtained using commercial computers.

6. Conclusion

A real-time damage location algorithm based on the wave propagation approach has been presented. Damage location is achieved by using an adaptive detection algorithm and a threshold design scheme that detects changes in the parameters. The threshold has been adjusted to obtain the best performance. The aim is to complement methods based on the wave propagation approach, which are traditionally developed in the frequency domain and are not suitable for real-time applications. Here, damage is assumed as stiffness loss and is identified by comparing the nominal shear wave velocities and the Kelvin damping coefficients with parameters obtained after a damage occurs. An advantage of the proposed method is the low computing time required for damage detection and location. Since structural parameters are identified in real time, the assessment of structural health can be carried out precisely at the same time. In addition to this, the method can be modified if a more precise location of the damage is desired by placing an additional sensor at the damaged story. Therefore, elements with structural damage can be identified in a specific floor. Furthermore, the proposed method has analytical properties that help to obtain robust results. In fact, the inclusion of linear integral filters in the parameterization removes constant disturbances and attenuates measurement noise from acceleration signals. Moreover, the parametric projection produces stronger convergence properties that improve the identification algorithm. For real-life engineering structures, damage identification, utilizing output-only acceleration measurements, is attractive and promising for practical applications. In this sense, the use of wave methods provides an interesting alternative to modal analysis methods.

We also recognized that it is necessary to carry out more and extensive research to assess the potential of this approach. For instance, damage produced by cracks, cuts, or drilled holes in columns needs to be explored, as well as experiments in x and y coordinates. Nevertheless, the experimental results presented here are a first step in the proper direction. It is convenient to emphasize that the use of controlled experiments ensures a fair and appropriate evaluation of the proposed algorithm, with respect to analytical values. Nevertheless, the proposed wave propagation approach is robust to several types and amounts of damage, as discussed in [30]. At this time, damage quantification is out of the scope of this work, and it is expected to be studied in the near future as structural damage may exhibit nonlinear behavior. This would imply that the amount of damage could not be directly proportional to the estimated wave propagation velocities.


A. Bounded and Differentiable Parameters

Led and , with , , be positive and bounded estimated parameters. They are differentiable, and their time derivatives are given bywhere and are provided by the NRLS algorithm. Additionally, variables and because the right-hand side term inside the brackets in equations (A.1) and (A.2) is also bounded. Consequently, is differentiable because this depends on estimated parameters and , which are also positive, bounded, and differentiable, such that .

Additionally, since and are matrices positive definite matrices and assuming is a diagonally dominant matrix, this ensures that all eigenvalues of matrix have strictly negative real part such that for each time, where .

B. Notation

Symbols used in this paper are shown in Table 7

Data Availability

Some of the data generated or used during the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


The first author was grateful to the Cátedras CONACYT. This research was sponsored by grant UNAM-PAPIIT IT101420.