Abstract
Introduction
Many load-carrying structural systems have been in use for tens or even hundreds of years. It is unavoidable for these structures to continuously accumulate damage in their service environment. Thus, it is most desirable that this damage be detected at the earliest possible stage to assure safety. Many research studies have been conducted during the past years in the area of damage detection–based structural response characteristics with different algorithms and useful databases.1–3 Generally, the existing damage detection approaches can be clarified into two major categories, that is, the dynamic detection methods and the static detection methods. Significant efforts have already been devoted to developing dynamic detection methods using various vibration characteristics, such as natural frequencies,4–6 mode shapes,7–11 modal flexibility,12–15 curvature,11,16 modal strain energy,17–19 and moving load responses.20–22 Radzieński et al. 10 defined new damage indicators based on the change in natural frequencies and any mode shape (measured or modeled) for beam-like structures. Kopsaftopoulos and Fassois 11 proposed a vibration-based statistical time series method that is capable of effective damage detection, precise localization, and magnitude estimation within a unified stochastic framework. Sung et al. 16 present a damage detection method for beam-like structure by the normalized uniform load surface (NULS) curvature obtained by modal flexibility. It has shown that changes in NULS curvature only occur at damaged elements in the structure. Xiang et al. 17 used the changes in damage-induced curvature mode shape and natural frequency to form a hybrid damage detection approach to detect damages on the surface of conical shell. Xu et al. 18 developed an energy damage detection strategy by disposing strain responses for a long-span cable-stayed bridge. Yi et al. 19 proposed a multi-stage method based on “energy-damage” theory using the wavelet packet transform (WPT) and the artificial neural network (ANN) for structural damage diagnosis. Using modal strain energy, Cha and Buyukozturk 20 developed hybrid multi-objective optimization algorithms to detect damages in various three-dimensional steel structures. Khorram et al. 21 compared the performances of two wavelet-based damage detection approaches to find the location and the size of a crack in a beam subjected to a moving load. Using Hilbert–Huang transform, Roveri and Carcaterra 22 presented a novel method for damage detection of bridge structures under a traveling load. The technique uses a single-point measurement and is able to identify the presence and the location of the damage along the beam. Cavadas et al. 23 discussed the application of data-driven methods on moving load responses to detect the occurrence and the location of damage. However, these dynamic detection algorithms are faced with some difficulties in applications. 24 One is that the precision of dynamic data are affected by many factors such as stiffness, mass, and damping. Therefore, many dynamic approaches unavoidably suppose that there is no damping and the mass is unchanged before and after damage. The next issue is that modal data for some stiff civil structures are difficult to be measured accurately, especially for high-order modes. In view of this, the static detection algorithms25–34 also attracted much attention in civil engineering area in recent years due to their simplicity for implementation. However, most of these methods need a detailed finite element model (FEM) or the baseline parameters of the intact system. It is sometimes impossible to obtain a detailed FEM or the baseline data for many existing structures. Thus, a simple and effective localization method that does not need an FEM or the baseline data of the intact structure will be welcomed.
For this reason, this article presents a static method for damage localization in beam structures by moving load. The key point of the damage location algorithm lies in the combination of the moving load technique and the pure bending theory. The first advantage of the presented algorithm is that only a few sensors are needed in the static test. The second advantage of the presented algorithm is that only the deflection parameters of the current structure are needed in damage localization process. The proposed method is verified by a numerical example to show its feasibility. It has been shown that the presented algorithm is effective on the localization of single damage or multiple damages in the beam structure.
Theory
Figure 1 presents a simple supported beam “A–B” under two symmetrical concentrated forces

Simple supported beam and corresponding bending moment diagram under two symmetrical concentrated forces.
According to the theory of material mechanics, the curvature of the deflection curve for this beam can be calculated by
where
On the other hand, the curvature can be calculated by the central difference technique of the deflection vector
where

Reciprocal theorem of displacement: (a) the load
On the other hand, according to superposition principle in linear elastic theory, the deflection of the beam in Figure 1 under two symmetrical concentrated forces will be equal to the summation of the deflections obtained by the two load cases as shown in Figure 3(a) and (b). Using the above moving load technique, the deflection vectors in Figure 3(a) and (b) can obtained as

Superposition principle of load: (a) the load
In the end, a summary of the overall process for this algorithm is given as follows:
Numerical example
Figure 4 presents a simple supported beam used to demonstrate the feasibility of the moving load method for damage localization. The beam is divided into 36 segments and the length of each segment is

The simple supported beam and sensor layout.
Damage cases in the example.
Curvature values when element 18 is damaged with 10%, 20%, and 40% stiffness reductions (no noise).
The bold values denote the curvature values corresponding to the damaged element.

Curvature values when element 18 is damaged (no noise).

Curvature values when element 18 is damaged (1% noise).
For damage cases 5 and 6, elements 9 and 18 can be determined from Figure 7 without noise to be damaged since their corresponding nodes 8, 9, 17, and 18 have larger curvatures than the other nodes. When 1% noise is considered, the damage localization result for small damage case 5 is unreliable from Figure 8. For large damage case 6, the curvature indexes in Figure 8 can still indicate the damage locations by the notable increases in the curvatures. In addition, one can conclude that the curvatures will increase as the damage extents increase from Figures 5–8. It has been shown that the proposed method is effective on the localization of single damage or multiple damages in the beam structure.

Curvature values when elements 9 and 18 are damaged (no noise).

Curvature values when elements 9 and 18 are damaged (1% noise).
It must be noted that the proposed method has some drawbacks as follows. The first is that the proposed approach can detect only damages that are between the two loading locations as shown in Figure 1. It is suggested to put the loads as close as possible to the support. This, in turn, complicates the things because deflections close to the support are extremely small and difficult to be measured with good accuracy. The second drawback is that the proposed method can detect only notable stiffness reduction in the structure. If the beam is a reinforced concrete one and rebars are cut at a certain location, the safety is dramatically reduced, while the stiffness is only slightly affected. This type of engineering disease cannot be detected by the proposed method. In addition, the support conditions of the beam structure have some impact on the proposed method. If the support conditions are non-ideal, the dynamic flexibility matrix obtained by the experimental modal parameters can be used to construct the virtual pure bending state for the beam structure. This issue will be studied in the further research.
Conclusion
A new deflection-based method has been proposed for beam damage localization by moving load. This method has two advantages: (1) only a few sensors are needed in static test and (2) only the deflection parameters of the current structure are needed in damage localization. To confirm the feasibility of the proposed method, a simple supported beam was investigated for several damage scenarios. The results showed that the presented algorithm is effective to identify structural damage locations. The presented algorithm may have good prospects in defect localization for beam structures without baseline data.
