Abstract

Frequency response function (FRF) data can provide considerably more information on damage in the desired frequency range as compared to modal data extracted from a very limited range around resonances. Among structural health monitoring techniques, FRF-based methods have the potential to locate structural damage. Conventional structural damage detection technology collects structural response data using contact systems, such as displacement or acceleration transducers. However, installing these contact systems can be costly in terms of labor, cost, and time. Several noncontact measurement technologies, such as optical, laser, radar, and GPS, have been developed to overcome these obstacles. Given the rapid advances in optical imaging hardware technology, the use of digital photography in structural monitoring systems has attracted considerable attention. This study develops a displacement FRF-based approach to locate damage to building structures. The proposed damage location index, CurveFRFDI, improves the sensitivity of SubFRFDI, which is a substructure FRF-based damage location index proposed by Lin et al. (2012). Moreover, the feasibility of applying the proposed approach to locate damage to building structures using displacement measured by a digital camera combined with digital image correlation techniques is also investigated in this study. A numerical example and an experimental example are presented to demonstrate the feasibility of using the proposed approach to locate damage to building structures for single and multiple nonadjacent damage locations.

1. Introduction

Vibration-based damage identification of structures refers to the in situ nondestructive sensing and analysis of the characteristics of a structure, including the structural response to external excitations, to detect changes that may indicate damage or degradation. The feasibility of applying various vibration characteristics, such as natural frequencies, mode shapes, mode shape curvatures, modal flexibility, and frequency response functions, to damage identification of structures has received considerable attention in the past few decades [1, 2].

Lee and Shin [3] identified two main advantages of using the frequency response function (FRF) data. First, modal data can be contaminated by modal extraction errors in addition to measurement errors because they are derived datasets. Second, a complete set of modal data cannot be measured in all but the simplest structures. FRF data can provide considerably more information on damage in the desired frequency range as compared to modal data that are extracted from a very limited range around resonances. Some studies have shown that FRF-based methods are highly promising tools to detect damage to building structures [48]. Ni et al. [4] identified the seismic damage of a 38-story building using measured FRFs and neural networks. Their study used principal component analysis (PCA) to reduce dimension and eliminate the noise of measured FRFs. The PCA-compressed FRF data are then used as input to neural networks for damage identification. Furukawa et al. [5] developed a method to detect damage to building structures that uses uncertain FRFs based on a statistical bootstrap method. Kanwar et al. [6] detected damage of reinforced concrete buildings using FRFs. Hsu and Loh [7] detected damage to building structures subjected to earthquake ground excitation using FRFs of intact and damaged systems as well as system matrices of the intact system to derive the damage identification equations. Lin et al. [8] proposed a substructure-based acceleration FRF approach that uses the index, SubFRFDI, to locate damage to building structures.

In addition to data analysis, data measurement is another problem that needs to be addressed in detecting damage to structures. Conventional structural damage detection technology collects structural response data by using contact systems, such as displacement or acceleration transducers. However, installing these contact systems can be costly in terms of labor, cost, and time [9]. Several noncontact measurement technologies, such as optical, laser, radar, and GPS, have been developed to overcome these obstacles [10]. The use of digital photography in structural monitoring systems has attracted considerable attention given the rapid advances in optical imaging hardware technology. Digital image correlation (DIC) is a measurement technique that extends the principles of photogrammetry to obtain full-field surface displacement measurements of an object using digital cameras. Several researchers have applied DIC to detect damage to building structures using dynamic responses [1116]. Shih et al. [11, 12] developed a low-cost digital image correlation method to measure the dynamic response of shear buildings. The accuracy of the DIC method is sufficiently high for several applications. Combe and Richefeu [13] used the DIC technique coupled with geometrical rules to develop an approach to track the nonsmooth trajectory of particles. Sieffert et al. [14] presented a digital correlation technique to capture the full-field displacement by using a high-speed camera of a full-scale structure tested on a shaking table. Lu et al. [15] presented a digital image processing approach with a unique hive triangle pattern by integrating subpixel analysis for noncontact measurement of structural dynamic response data. Hung and Lu [16] integrated this approach and GPU computing technique to save on computation time.

Many civil structures, especially high-rise buildings and long-span bridges, have low fundamental response frequencies. Measuring displacement directly, as opposed to integrating acceleration records (and potentially introducing significant error [17]), provides the opportunity to capture low-frequency response [18]. Detecting structural damage using displacement FRF should be more accurate than using acceleration FRF because the fundamental frequency of a structure becomes lower when a damage occurs. Based on these reasons, the objective of this paper is to enhance the work of Lin et al. [8]. In this paper, a damage location index (CurveFRFDI) based on SubFRFDI curvature and displacement FRF is used to locate damage to building structures to enhance the sensitivity of SubFRFDI to damage. The study also aims at investigating the feasibility of applying the proposed approach to locate damage to building structures using displacement measured by digital camera combined with the DIC technique. Moreover, a numerical example and an experimental example demonstrate the feasibility of applying the proposed approach for locating damage to building structures.

2. Damage Location Strategy

Lin et al. [8] presented a substructure-based FRF approach to locate damage to building structures. As presented in Figure 1, the ith substructure is a structure containing the ithNth floors (or degree of freedoms, DOFs) for an N-story building structure. The substructure-based FRFs of the jth DOF in the ith substructure can be simplified as follows:where and are the Fourier transforms of and , respectively; and are the absolute acceleration of the jth and (i − 1)th DOF, respectively. Theoretically, when the damage is assumed to have occurred in the columns between the ith and (i − 1)th DOFs, the substructure-based FRF is altered significantly in the ith DOF as described by . For efficiency, only one substructure-based FRF, , is determined for each substructure to reduce the computational time. Damage location is based on the FRFs, , , …, , of all substructures.

The substructure-based FRF damage location index (SubFRFDI) for the ith substructure is defined as follows:where , , , and are working parameters and the is expressed as follows:and the absolute dissimilarity is defined as follows:where and are the magnitudes of in the damaged and undamaged states, respectively. These N dissimilarities, , can be calculated correspondingly for a shear building with N floors. In equation (2), the coefficient is a control value that scales the index between 0 and 1. The damage index is not sensitive to locate damage for a small . However, the sensitivity of the damage index to locate damage will not increase for a large . The value of the coefficient is suggested to be 1∼10, with 5 used in this work. The range of selected frequencies for calculating SubFRFDI is set to be ab, where a is the starting frequency of zero and b is the end frequency, which is equal to the first modal frequency (undamaged state). The value n equals (b-a) divided by sampling time. If the properties of a structural system do not change, then is close to zero. However, if the damage to the ith floor in a shear building is severe, then the value of is high. Thus, can be regarded as the SubFRFDI corresponding to location i (the ith floor).

This study investigates the feasibility of applying displacement measured by digital camera combined with DIC technique to locate damage to building structures, and thus, the substructure-based FRF of the jth DOF in the ith substructure in equation (2) is modified to calculate :where and are the Fourier transforms of and , respectively; and are the absolute displacement of the jth and (i − 1)th DOF, respectively.

In fact, the SubFRFDI values corresponding to damaged and undamaged locations increase with the increasing of damage extent. Thus, SubFRFDI is insensitive to damage with larger damage extent. Moreover, SubFRFDI can locate single and multiple nonadjacent damages but is unable to locate multiple adjacent damages to building structures. Damage occurred at the location corresponding to a larger SubFRFDI value than adjacent locations for single and multiple nonadjacent damage cases. That is, if the ith floor is damaged, and the (i − 1)th and (i + 1)th floors are undamaged, SubFRFDIi is larger than SubFRFDIi−1 and SubFRFDIi+1, and the curve connecting SubFRFDIi−1, SubFRFDIi, and SubFRFDIi+1 is open downward. Thus, damage occurred at the location corresponding to a negative SubFRFDI curvature rather than a nonnegative SubFRFDI curvature for single and multiple nonadjacent damages. The curvature of SubFRFDIi, , can be simply defined as follows:where and are both equal to zero. In this study, a damage location index based on SubFRFDI curvature, CurveFRFDI, is used to locate damage to building structures to enhance the sensitivity of SubFRFDI to single and multiple nonadjacent damages. The CurveFRFDI for an N-story building structure is expressed as follows:where A is a control value that scales CurveFRFDIi to a very small positive value for positive SubFRFDIi curvature. The value of coefficient A is large compared with SubFRFDI curvature. It is suggested that , and A is set to be 100 in this study. The range of CurveFRFDIi is from 0 to . Like SubFRFDI, damage occurred at the location corresponding to a larger value of CurveFRFDI. Figure 2 shows the flowchart of the proposed approach for locating damage to building structures.

3. Digital Image Correlation

Digital image correlation (DIC) is an easy optical method to measure continuous deformation by tracking the position of the same physical points shown in a reference image and a deformed image. A rectangular subset of pixels is identified on the speckle pattern around a point of interest (POI) on a reference image and their corresponding location determined on the deformed image to achieve tracking.

This study employs the DIC approach presented in the work of Lu et al. [15] for displacement measurement in the experimental example. The approach used a unique hive triangle pattern by integrating subpixel analysis for noncontact measurement of structural dynamic response data. As shown in Figures 3(b) and 3(c), a fixed-size rectangular subset of pixels that make up hive triangle patterns is the region of interest (ROI), contained both within the reference (source) and within the deformed (target) images and marked with a red color box. Meanwhile, in Figure 3(a), they are designated as I and I′ in source and target images, respectively. The and are the coordinates of the points at the left-top of the ROI of the source image and the after deformed target images, respectively. Hence, u and are the relative deformations of a particular point in image space. The coordinates of and are figured out in the following steps. Two images are first selected as S and T. The image S is an undeformed (source) image, and image T is a slightly deformed image (target) relative to image S. The difference between the two images, S and T, can be estimated through a simple difference method via pixel-wise computation, as shown in the following equation:where i and j are the index of a pixel in the ith row and jth column of the source and target images based on the origin at the top-left corner. If the intensity of each pixel is from 0 to 255, the difference of each pixel between S and T is maximal when the equation value is 255, and S and T are exactly identical when the equation value is 0. If two images have the same background, the of the pixels in the background area is close to zero. Therefore, if an ROI changes the position due to deformation, the of these pixels covered in the ROI is relatively large.

Moreover, the coordinate, , of the ROI of the source image can be computed using a mean-max method represented in equation (9). The method first calculates the means of the intensities of all columns and rows pixels in the source and one of the target images having slightly deformation. The index with the maximum mean values of the columns and rows indicates the x-axis and y-axis coordinates, respectively:where and in which and h are the width and height of the source and target images in pixel. Particularly, the approach of Lu et al. employs an ideal hive triangle pattern as the ROI in calculating the correlation coefficient with the source image. The correlation coefficient used in the work is defined in the following equation:where f and are the pixel value of ROIs for the source and target images, respectively. The sign denotes the mean operator and and are the means of ROIs in the source and target images, respectively. The coordinate, , of the point at the left-top of the ROI of the after deformed target images can be figured out based on the following equation:where I and I′ denote the aforementioned ROIs for source and target images, and CC refers to the correlation coefficient function to evaluate with the two ROIs. The relative displacement in the unit pixels (i.e., the pixel displacement), , in the x-axis and y-axis coordinates can be estimated to be . The precision of the pixel displacement, , is an integer pixel value with original images. The actual length (L) of the pattern is a known quantity, and the pixel width (W) of the pattern has been evaluated by the measurement system. The actual displacement, , can be calculated by using estimated pixel displacement, , to multiply by a pixel ratio (Rp) from equation (12). At this point, the time-history displacements are estimated completely in the digital image measurement system:

The subpixel analysis can improve the precision based on the subpixel estimation of a target image. The results of Lu et al. [15] indicated that the measurement system increases the precision to a pixel value of 0.1, or even 0.01, and the precision achieves 0.01 mm if is less than 0.1 mm. Moreover, the maximum correlation coefficient is higher than 0.95 in the correct search and lower than 0.75 in the fail search. In the fail situation, another target image is select to evaluate the maximum correlation coefficient again.

By employing a noncontacted approach for measurement of structural dynamic response data, the record can be divided into 9,000 image frames if the time history of displacements is recorded as a video that contains 90 seconds with 100 fps data. The first frame will be set as a source (reference) image, and other frames are the target images and are processed sequentially to obtain the coordinates (x1, y1) in each target image; consequently, the time history of displacements can be figured out. Figure 4 graphically presents a series of frames in a video and the corresponding computed time history of displacements.

4. Examples

To confirm the feasibility of the proposed approach for locating damage to building structures, two examples, a numerical example and an experimental example, are studied.

4.1. Numerical Example

As shown in Figure 5, a six-story shear plane frame structure is studied in this example to evaluate the feasibility of the proposed approach for locating damage. Each floor consists of one beam and two columns. The cross-sectional size of the beam is 240 mm × 30 mm. The length of the beam is 360 mm. The cross-sectional size of the column is 240 mm × 15 mm. The length of the column is 180 mm. All members (beams and columns) are assumed to be made from the same material, ASTM A992 steel, with a yield stress of . In this example, SAP2000 is used for structural analysis. The natural frequencies of the first six modes are 1.653 Hz, 4.975 Hz, 8.341 Hz, 11.472 Hz, 14.130 Hz, and 15.938 Hz. The damage in this example is simulated as reduced floor stiffness. Single and multiple damage locations are studied. Table 1 presents the simulated damage cases. The 1995 Kobe earthquake record is used as the external excitation, but normalized to 100 gal as the peak ground acceleration (PGA).

4.1.1. Damage Location Using SubFRFDI

Figure 6 compares the substructure-based FRFs for the 2nd and the 5th substructures ( and ) for case Dam_2F15%. Significant changes were observed only in the substructure-based FRF of substructure 2, , which demonstrates the ability of the substructure-based FRF to characterize and locate damage.

Figures 715 present the comparison of the SubFRFDI values for damage cases listed in Table 1. Results show that although the damage locations of all cases shown in Figures 79 are predicted correctly by SubFRFDI, false detection could occur if the number of damage location is unknown. For example, the 1st and 3rd floors for case Dam_2F&4F15% (see Figure 8(a)), the 2nd floor for case Dam_3F&5F15% (see Figure 8(b)), the 1st floor for case Dam_2F&3F15% (see Figure 8(c)), and the 1st floor for case Dam_2F&3F&4F15% (see Figure 9(b)) may be detected falsely because their corresponding SubFRFDI values are large. Figures 1012 show that SubFRFDI values corresponding to damaged and undamaged locations increase with the increase in damage extent for cases of single and multiple damage locations. Thus, false detection could occur for large damage extent cases if the number of damage location is unknown.

For cases of multiple nonadjacent damage locations (cases Dam_2F&4F15% (Figure 8(a)), Dam_3F&5F15% (Figure 8(b)), and Dam_2F&4F&6F15% (Figure 9(a))), the damage locations can be predicted by the location (the lowest floor of the substructure) corresponding to the local maximum value of SubFRFDI. However, this criterion is not suitable for cases of multiple adjacent damage locations such as cases Dam_2F&3F15% (Figure 8(c)) and Dam_2F&3F&4F15% (Figure 9(b)).

The stiffness reduction of the 2nd floor (10%) is smaller than that of the 4th floor (15%), but SubFRFDI2 is larger than SubFRFDI4 for case Dam_2F10%&4F15% (see Figure 13(a)). However, the stiffness reduction of the 2nd floor (10%) is smaller than that of the 4th floor (20%) and SubFRFDI2 is smaller than SubFRFDI4 for case Dam_2F10%&4F20% (see Figure 13(b)). The results imply that the damage extent of the ith floor is not absolutely larger than that of the jth floor if SubFRFDIi is larger than SubFRFDIj. Moreover, false detection could occur for cases Dam_2F20%&4F10%, Dam_3F15%&5F10%, Dam_3F20%&5F10%, and Dam_2F15%&4F10%&6F5% because SubFRFDI1 (0.293) is a little larger than SubFRFDI4 (0.286) for case Dam_2F20%&4F10% (see Figure 13(d)), SubFRFDI2 (0.400) is just a little smaller than SubFRFDI5 (0.443) for case Dam_3F15%&5F10% (see Figure 14(c)), SubFRFDI2 (0.592) is larger than SubFRFDI5 (0.438) for case Dam_3F20%&5F10% (see Figure 14(d)), and SubFRFDI3 (0.199) is just a little smaller than SubFRFDI6 (0.214) for case Dam_2F15%&4F10%&6F5% (see Figure 15(c)). The results imply that SubFRFDI is unable to locate damage for cases of two (three) damage locations when the damage extent (stiffness reduction) is equal to and larger than 10% (5%).

4.1.2. Damage Location Using CurveFRFDI

Figures 1624 present the comparison of the CurveFRFDI values for damage cases listed in Table 1. Results show that CurveFRFDI has higher sensitivity to damage than SubFRFDI for cases of single and multiple damage locations as compared Figures 1624 with Figures 715. Nevertheless, results imply that CurveFRFDI is also unable to locate damage for cases of multiple adjacent damage locations. For example, the 3rd floor for case Dam_2F&3F15% (see Figure 17(c)) and the 3rd floor for case Dam_2F&3F&4F15% (see Figure 18(b)) are falsely detected.

For most cases in Figures 1921, the CurveFRFDI value corresponding to the damaged location increases with the increase in damage extent while that corresponding to the undamaged location is almost zero and varies slightly with the increase in damage extent (only in Figures 19(c) and 20(b), the CurveFRFDI value corresponding to the second floor wrongly increases with the damage extent). It implies that CurveFRFDI can locate damage regardless of intensity (extent) for most cases. Notably, detecting small extent damage is an important issue for early warning of structural health monitoring. The excellent ability of CurveFRFDI to locate small extent damage can be investigated from case Dam_2F5% in Figure 19(a), case Dam_3F5% in Figure 19(b), case Dam_4F5% in Figure 19(c), case Dam_2F&4F5% in Figure 20(a), case Dam_3F&5F5% in Figure 20(b), and case Dam_2F&4F&6F5% in Figure 21.

4.2. Experimental Example

In this example, the dynamic responses of an eight-story downscale steel frame (Figure 25), subjected to the Chi-Chi earthquake for different values of PGA (50 gal, 100 gal, 200 gal, 500 gal, and 1200 gal) shaking table tests, are processed to demonstrate the feasibility of the proposed approach for locating damage. Each floor of the steel frame has dimensions of 1500 mm in length, 1100 mm in width, and 1180 mm in height (Figure 26). Each floor consists of one slab and four columns. Each slab consists of two C 100 × 50 × 5 with length 1198 mm, two C 100 × 50 × 5 with length 798 mm, one Plate 1500 × 1100 × 20, and four Plate 150 × 150 × 50. Each column consists of four H 100 × 100 × 6 × 8 with length 1060 mm, four Plate 1060 × 150 × 25, and eight Plate 150 × 150 × 25. Table 2 lists the steel types used to construct the eight-story steel frame. These series of shaking table tests were undertaken by the National Center for Research on Earthquake Engineering in Taiwan. The displacements response histories of each floor are measured during the shaking table tests through linear variation differential transformation (LVDT) and a digital camera (Basler A504kc, sampling rate of 500 Hz) combined with the DIC approach presented in the work of Lu et al. [15], abbreviated as LVDT-measured data and DIC-measured data below. The damage in this example is simulated by reducing the cross section of certain columns as shown in Figure 27. Single and multiple nonadjacent damage locations are studied. Table 3 presents the simulated damage cases.

4.2.1. Damage Location Using LVDT-Measured Data

Figure 28 presents the comparison of the SubFRFDI values of LVDT-measured data for cases of single damage location, Dam50_1F and Dam100_1F. The damage locations of the two cases are predicted correctly by the SubFRFDI of LVDT-measured data. Figure 29 shows the comparison of the SubFRFDI values of LVDT-measured data for cases of two damage locations, Dam50_1F&3F, Dam100_1F&3F, Dam200_1F&3F, Dam500_1F&3F, and Dam1200_1F&3F. It shows that SubFRFDI1 and SubFRFDI3 are larger than others, and the damage locations (the 1st and the 3rd floors) of these cases are predicted correctly by the SubFRFDI if the number of damage location is known in advance. Nevertheless, SubFRFDI3 is much smaller than SubFRFDI1 and false detection of the 3rd floor may be occurred if the number of damage location is unknown.

Figure 30 presents the comparison of the CurveFRFDI values of LVDT-measured data for cases of single damage location, Dam50_1F and Dam100_1F. The damage locations of the two cases are also predicted correctly by the CurveFRFDI of LVDT-measured data. Figure 31 shows the comparison of the CurveFRFDI values of LVDT-measured data for cases of two damage locations, Dam50_1F&3F, Dam100_1F&3F, Dam200_1F&3F, Dam500_1F&3F, and Dam1200_1F&3F. Although CurveFRFDI3 is much smaller than CurveFRFDI1 and false detection of the 3rd floor may be occurred if the number of damage location is unknown, CurveFRFDI has higher sensitivity to damage than SubFRFDI as compared Figures 3031 with Figures 2829. To solve this difficulty, a threshold CurveFRFDI value is suggested to be set to locate damage. The threshold CurveFRFDI value can be determined after numerous numerical simulations of damage scenarios for a certain building structure. For example, the threshold CurveFRFDI value can set to be 0.08 for this experimental example.

4.2.2. Damage Location Using DIC-Measured Data

Figure 32 presents the comparison of the SubFRFDI values of DIC-measured data for case Dam100_1F. Figure 33 shows the comparison of the SubFRFDI values of DIC-measured data for cases of two damage locations, Dam200_1F&3F and Dam500_1F&3F. Figure 34 presents the comparison of the CurveFRFDI values of DIC-measured data for case Dam100_1F. Figure 35 shows the comparison of the CurveFRFDI values of DIC-measured data for cases of two damage locations, Dam200_1F&3F and Dam500_1F&3F. The results also show that the CurveFRFDI has higher sensitivity to damage than SubFRFDI as compared Figures 3435 with Figures 3233. The results further indicate that applying the proposed approach to locate single and multiple nonadjacent damages to building structures is feasible by using DIC-measured displacements.

The experimental example proves the feasibility of the proposed approach for location damage to building structures using DIC-measured displacement response. These data are supposed to be noise free or low noise corrupted. Future work should apply the proposed approach to measurements in the field (actual cases) to investigate its capacity to deal with DIC-measured displacement response corrupted by high noise.

5. Conclusions

The current work develops an approach for locating damage to building structures. The proposed damage location index, CurveFRFDI, improves the sensitivity of SubFRFDI, a substructure FRF-based damage location index proposed by Lin et al. [8]. The feasibility of applying the proposed approach to locate damage to building structures using DIC-measured displacement is investigated in this study. A numerical example and an experimental example are presented to demonstrate the feasibility of using the proposed approach to locate damage to building structures. The following important conclusions are drawn from the results:(1)Using SubFRFDI can predict damage location accurately for cases of single damage location. However, SubFRFDI cannot locate damage for cases of multiple damage locations with large damage extent. SubFRFDI values corresponding to damaged and undamaged locations increase with the increase in damage extent for cases of single and multiple damage locations. Thus, false detection may be occurred for large damage extent cases if the number of damage locations is unknown.(2)In most cases, the CurveFRFDI value corresponding to damaged location increases with the increase in damage extent while that corresponding to undamaged location is almost zero and varies slightly with the increase in damage extent, which implies that CurveFRFDI can locate damage regardless of intensity (extent) for most cases. Moreover, CurveFRFDI has higher sensitivity to damage than SubFRFDI. Thus, using CurveFRFDI can predict damage location more accurately than using SubFRFDI for cases of multiple damage locations with large damage extent.(3)For cases of multiple damage locations, some CurveFRFDI values corresponding to damaged location may be much smaller than others, making it difficult to detect those damage locations. To solve this difficulty, a threshold CurveFRFDI value is suggested to be set to locate damage. The threshold CurveFRFDI value can be determined after numerous numerical simulations of damage scenarios for a certain building structure.(4)Applying the proposed approach and DIC-measured displacements to locate single and multiple nonadjacent damages to building structures is feasible.(5)Both SubFRFDI and CurveFRFDI cannot locate damage for cases of multiple adjacent damage locations. How is the proposed approach applicable to the prediction of multiple adjacent damage locations should be investigated based on this research.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank the Ministry of Science and Technology of the Republic of China for financially supporting this research under Contract nos. MOST 103-2625-M-009-003 and MOST 104-2221-E-009-049-MY2.