Research Article  Open Access
Damage Identification of Bridge Based on Chebyshev Polynomial Fitting and Fuzzy Logic without Considering Baseline Model Parameters
Abstract
The paper presents an effective approach for damage identification of bridge based on Chebyshev polynomial fitting and fuzzy logic systems without considering baseline model data. The modal curvature of damaged bridge can be obtained through central difference approximation based on displacement modal shape. Depending on the modal curvature of damaged structure, Chebyshev polynomial fitting is applied to acquire the curvature of undamaged one without considering baseline parameters. Therefore, modal curvature difference can be derived and used for damage localizing. Subsequently, the normalized modal curvature difference is treated as input variable of fuzzy logic systems for damage condition assessment. Numerical simulation on a simply supported bridge was carried out to demonstrate the feasibility of the proposed method.
1. Introduction
Bridge structures have endured progressive deterioration due to vehicle load and aging of material under the effect of external environment. In addition, the frequent occurrence of natural disasters could also cause damage and further accelerate the deterioration of bridge condition. Therefore, the damage assessment of bridge has attracted great interest of many researchers [1, 2].
Vibrationbased damage identification method has been widely applied in the past few decades [3, 4]. Theoretical background of this method is that modal parameters (frequency, mode shape, etc.) are functions of physical properties (stiffness, mass, etc.) of bridge. Therefore, the changes of modal characteristics can be treated as the damage indicators. The commonly used modal parameters for damage identification include natural frequency, mode shape, and their derivatives, such as mode shape curvature, modal strain energy, and modal flexibility [5–9].
Among these modal characteristics, natural frequency has no relation with measuring positions and can be measured most conveniently and accurately. However, the frequencybased method possesses several apparent drawbacks. For instance, it is inferior sensitive to minor damage, and it is easily affected by environment [10]. Additionally, the simply supported reinforced concrete bridges are spatial symmetric structures. Natural frequency cannot distinguish damage at symmetric locations in these symmetric structures.
Comparing with natural frequency, mode shape and its derivatives contain the spatial information with respect to location of damage, and they are found to be better indicators for damage identification. The modal curvature is recognized as a more advanced damage indicator which is calculated from the displacement mode shape and firstly proposed by Pandey et al. [7]. It is found that the absolute changes in the curvature mode shapes between undamaged and damaged structures can effectively locate the damage region for cantilever and simply supported beam model. The changes in the curvature mode shape increase with increasing size of damage. However, this method requires the baseline data from intact structures. This baseline model data can be obtained by finite element simulation. In actual simulation, the accuracy of finite element model is affected by temperature because the elastic modulus of concrete is temperature dependent [11]. However, the modulus of elasticity versus temperature curve for concrete is unclear and complicated [12, 13]. If their relationship is not reasonable to consider, it will lead to incorrect result for damage identification.
To avoid this difficulty, several researchers have conducted effective work for damage detection without baseline model data. Ratcliffe [14] proposed a modified Laplacian operatorbased method for 1D beam, which determines the damage location through a cubic curve fitting for the modal data obtained only from the damaged structure. Wu and Law [15] proposed a damage localization method based on uniform load surface (ULS) curvature for 2D plates, which can identify the damage using only the modal characteristics of damaged state if a gappedsmoothing technique is applied. Yoon et al. [16] extended the 1D gappedsmoothing method to the 2D gapped smoothing method for the damage identification of plate. And the baseline data of undamaged structure are not needed. Zhong and Oyadiji [17] presented a novel approach for crack detection based on difference between two sets of detail coefficients obtained by stationary wavelet transform without requiring the modal data of undamaged beam as a baseline. Yoon et al. [18] extended the gapped smoothing method for identifying the location of structural damage in a beam by introducing the global fitting method, which uses the modal shape data of damaged structure with an assumption that the undamaged structure is homogeneous and uniform. Cao and Qiao [10] adopted a novel Laplacian scheme for modal curvaturebased damage identification to improve the antinoise ability of standard Laplace operator. Rodríguez et al. [19] presented a Baseline Stiffness Method (BSM) to locate and evaluate the magnitude of structural damage whose baseline state is unknown.
Modal curvature can realize the damage localization and qualitatively determine the damage degree. But it is not suitable for quantitative identification. Artificial neural networks (ANNs) have been utilized by many researchers to identify damage location and severity [20–25], as they can achieve the nonlinear mapping between the inputs and outputs from certain samples training. Ko et al. [20] developed a threestage scheme for damage detection of the cablestayed Kap Shui Mun Bridge. ANNs are used in the first and third stage for damage alarming and specific damaged member(s) identification. Sahin and Shenoi [21] presented a damage detection algorithm using a combination of global (changes in natural frequencies) and local (curvature mode shapes) vibrationbased analysis data as input in ANNs for location and severity prediction of damage in beamlike structures. Lee et al. [22] proposed an ANNsbased damage detection method using the differences or the ratios of the mode shape components as the input variables, which can effectively consider the modeling errors in the baseline finite element model. Bakhary et al. [23] proposed a statistical approach which takes the effect of uncertainties into account in developing an ANNs model.
ANNs can be able to predict an output pattern when they recognize a given input pattern, but there exists an inherent inability to represent knowledge acquired by the network in an explicit form. That is the “black box” problem of ANNs [26, 27]. Fuzzy systems allow for easier understanding because they are expressed in terms of linguistic variables [28, 28]. And they are finding increasing use in structural damage identification. Ganguli et al. [29, 30] adopted the natural frequency and modal shape curvature as the input of fuzzy logic systems and genetic fuzzy logic systems for the damage detection of helicopter rotor blades. Zhao and Chen [31] proposed a method based on principal component analysis, modified mountain clustering method, descent method, and fuzzy logic systems for the damage detection of concrete bridges. Taha and Lucero [32] introduced new techniques based on Bayesian updating and fuzzy sets to improve pattern recognition and damage detection of structures.
In this study, a new acquisition method for modal curvature difference (MCD) is proposed. Firstly, the modal curvature of damaged structure can be calculated through central difference approximation based on mode shape data. Further, modal curvature of undamaged structure can be obtained by Chebyshev polynomial fitting based on the modal curvature at feature nodes of damaged structure. Therefore, MCD parameters can be calculated by use of modal curvatures before and after damage. Subsequently, a fuzzy logic system for damage condition assessment of bridge is constructed using normalized modal curvature difference (NMCD) as input variable. Numerical simulation of a simply supported beam bridge verifies the feasibility of the proposed method.
2. Theoretical Background
2.1. Modal Shape Curvature
The curvature of EulerBernoulli beam can be expressed by where is the curvature of displacement curve and and are the first and second derivatives of displacement function, respectively. Based on the assumption of small slope of beam displacement, the first derivative is close to zero; thus the denominator in (1) is equal to one. Therefore, the following expression can be used for calculating the curvature of beam:
Considering that central difference can obtain a more accurate approximation, its calculation results satisfy the project requirements [33]. Therefore, it is widely applied to estimate modal curvatures from the mode shapes in practice. It can be expressed by where represents modal curvature, represents the mode number, represents the corresponding node number, is the element length, and is the mode shape value for the th node of the th mode.
The modal curvature difference () is obtained by subtracting the undamaged modal curvature vector from the respective damaged one, given by where and are the modal curvatures for damaged and undamaged structure, respectively.
2.2. Chebyshev Polynomial Fitting
The Chebyshev polynomials are a sequence of orthogonal polynomials which are related to de Moivre’s formula and which can be defined recursively. They are important in approximation theory because they can be used for polynomial interpolation. The resulting interpolation polynomial minimizes the problem of Runge’s phenomenon and provides an approximation that is close to the polynomial of best approximation to a continuous function under the maximum norm [34].
2.2.1. Definition of Chebyshev Polynomial
The order Chebyshev polynomial can be defined by and its recursive formulas are as follows:
2.2.2. Functional Expression Using Chebyshev Polynomial
An orthogonal group can be constructed through Chebyshev polynomials, and the approximating polynomial of function at independent variable interval can be expressed by where are the orthogonal polynomials group and are the corresponding coefficients.
2.2.3. Curve Fitting Procedure Using Chebyshev Polynomial
(1) Construction of Normalized Interval. If , it can be transformed into the normalized interval by where , , and is the number of .
(2) Acquisition of Normalized Data Table. Through the construction of normalized interval, the data table for can be transformed into .
(3) Access to Zero Point of Order Chebyshev Polynomial. The zero point of order Chebyshev polynomial can be calculated by where is the degree of the polynomial. Generally, .
(4) Acquisition of Corresponding to Zero Point . The corresponding to zero point in data table can be obtained through interpolation method.
(5) Determination of Coefficients . The for orthogonal polynomials group can be calculated by
(6) Functional Relationship between and . The data table can be expressed by the following function:
(7) Fitting Error between and . The Fitting error between and can be obtained by the following mean square error function:
2.3. Fuzzy Logic Systems
A fuzzy logic system (FLS) is an information processing system that simulates the fuzzy inference ability of human brain. It is able to simultaneously handle numerical data and linguistic knowledge. It is a nonlinear mapping of an input data vector into a scholar output. Fuzzy set is characterized by a membership function which takes on values in the interval . The Mamdani system and TakagiSugeno system are the most widely used fuzzy logic systems. In the Mamdani system, fuzzy sets are used both in the antecedents and in the consequents of the rules, while the consequents are represented by functions in TakagiSugeno system. In this paper, the Mamdani system is adopted to carry out the damage identification. The fuzzy logic controller is composed of fuzzification, fuzzy inference, and defuzzification as shown in Figure 1.
3. Numerical Example
A numerical analysis is performed on a simply supported bridge as shown in Figure 2. The model consists of 21 nodes and 20 elements with equal element length ( m) and uniform bending rigidity ( kN·m^{2}).
According to the theory of modal analysis, natural frequency and modal shape are closely related with rigidity () and mass () of bridge. In general, the damage of bridge does not alter its mass but the rigidity. Therefore, the reduction of bridge rigidity caused by damage is the main factor leading to modal variation, as follows [35]: where is the number of damaged elements and , represent the rigidity of undamaged and damaged elements, respectively.
Taking the rectangular beam, for example, is the section height and is the section width. The rigidity of this section can be expressed by where is elastic modulus and is bending moment of inertia.
As can be seen from (14), the rigidity of bridge contains two aspects, that is, elasticity and section size ( and ). In the present research, there are usually two simulation methods for damage: one assumes that the section size of bridge is changed, while the other assumes that the elasticity of bridge material is altered. In this paper, the latter method is adopted and the damage parameter is defined by where is the number of damaged elements and , represent elasticity of undamaged and damaged elements, respectively.
3.1. Modal CurvatureBased Damage Localization without Baseline Model Parameters
Firstly, modal shape of damaged structure can be obtained by finite element analysis of bridge. It can be expressed by : is the node number and is the displacement mode shape, and the first vertical bending mode () is discussed in this paper. The modal curvature for damaged structure can be calculated through central difference approximation.
Secondly, the nodes related to damaged elements are excluded and the modal curvatures at remaining nodes are fitted through Chebyshev polynomial as shown in (11). And the fitting effect can be evaluated by (12). The normalized intervals are listed in Table 1. In this paper, sixorder Chebyshev orthogonal polynomials group are used for curve fitting, and the zero points are listed in Table 2.


Finally, modal curvature difference is acquired based on the above results of modal curvature for undamaged and damaged structure.
3.1.1. Damage Localization for Single Location
Taking the damage localization of element 10 with rigidity reduction (10%, 20%, 30%, and 40%), for example, modal curvatures at each node obtained by central difference approximation for damaged structure are shown in Figure 3.
The modal curvatures of undamaged structure corresponding to each damage severity are obtained by Chebyshev polynomial fitting by use of the mode shape data at nodes (1, 2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 19, 20, and 21). Taking the damage severity 30%, for example, the modal curvatures of undamaged and damaged structure are shown in Figure 4.
The fitting function of modal curvature at damage severity 30% for undamaged structure is shown as follows: where is the modal curvature at damage severity 30% for undamaged structure and is the normalized interval obtained by (8).
The fitting error is calculated according to (12), and ; it reveals that the Chebyshev polynomial fitting possesses favorable accuracy.
Modal curvature differences of element 10 with rigidity reduction (10%, 20%, 30%, and 40%) are calculated based on the above results, and they are shown in Figure 5.
3.1.2. Damage Localization for Multiple Locations
Taking the damage localization of elements 5, 10, and 16, for example, damage cases are shown in Table 3.

Modal curvature at each node calculated by central difference approximation for damaged structure is shown in Figure 6.
The modal curvatures of undamaged structure are obtained by Chebyshev polynomial fitting using the mode shape data at nodes (1, 2, 3, 4, 7, 8, 9, 12, 13, 14, 15, 18, 19, 20, and 21). Taking case 4, for example, the modal curvatures of undamaged and damaged structure are shown in Figure 7.
The fitting function of modal curvature of case 4 for undamaged structure is shown as follows:
And the fitting error .
Modal curvature differences for all cases in Table 3 are calculated based on the above results, and they are shown in Figure 8.
3.2. Fuzzy LogicBased Damage Condition Assessment
3.2.1. Calculation Process of Fuzzy Logic
(1) Inputs and Outputs. The objective of fuzzy logic is to map an input vector represented by to an output one represented by . In this paper, the vectors and can be expressed by where is the th input variable which is the normalized modal curvature difference (); is the modal curvature difference at node of mode ; is the damage severity; and and are the same with (15).
(2) Fuzzification. The fuzzification comprises the process of transforming crisp values into grades of membership for linguistic terms of fuzzy sets. The membership function is used to associate a grade to each linguistic term.
Inputs of the fuzzy logic system are and outputs are structural damage degree of elements. For input variables, we define nine membership functions to transform into linguistic terms of fuzzy sets. It is considered as a linguistic variable. It can be decomposed into a set of terms: = {negligible, very low, low, low medium, medium, medium high, high, very high, very very high} = {N, VL, L, LM, M, MH, H, VH, VVH}.
Fuzzy sets with Gaussian membership functions are used for the input variables. These fuzzy sets can be defined using the following equation and they are shown in Figure 9: where is the midpoint of the fuzzy set and is the standard deviation associated with the variable. The midpoints are selected to span the region ranging from an undamaged bridge to one with significant damage, while the standard deviations are selected in order to guarantee the enough width to capture the variations of the input variables. Therefore, the reasonable selection of and is very important for the performance of the fuzzy logic systems.
For output variables, we define five membership functions to transform the damage degree into linguistic terms of fuzzy sets. For example, “element 10” is considered as a linguistic variable and it is decomposed into a set of terms: (element 10) = {undamaged, slight damage, moderate damage, severe damage, catastrophic damage}.
The other output variables are fuzzified in a similar manner.
Fuzzy sets with Bell functions are used for the output variables. These fuzzy sets can be defined using the following equation: where are the variables that control the shape of function.
The parameters of Bell function used in this paper are shown in Table 4 and Figure 10.

(3) Rule Base. Once the input and output variables and their membership functions are defined, we can design the rule base which is composed of “if <antecedents> then <consequents>” rules. The antecedents and consequents are the input and output variables, respectively.
Rules for a fuzzy logic system are obtained by fuzzification of the numerical values obtained from finite element analysis using the following procedure.(i)The th input variable corresponding to a given structural fault is fuzzified and the corresponding membership degrees can be obtained. Therefore, each measurement has nine grades of membership based on the linguistic terms.(ii)Each input variable is then assigned to the fuzzy set with the maximum degree of membership.(iii)One rule is obtained for each fault by relating the fuzzified input variable with maximum membership degree of the fault.
3.2.2. Damage Condition Assessment for Single Location
Taking the damage condition assessment of element 10, for example, at nodes 10 and 11 are adopted as the input variables of fuzzy logic system, and the rule base is formulated and presented in Table 5.

In order to verify the inference ability of the established fuzzy system, test samples are selected and listed in Table 6 and their damage conditions corresponding to severity are also presented. The damage identification results for these test samples using fuzzy system are demonstrated in Table 7.


Damage condition of test samples can be determined through the maximum membership degree. As can be seen from Table 7, fuzzy systems can correctly identify the damage condition for all test samples. It reveals that the fuzzy systems possess favorable inference ability for damage condition assessment of single location.
3.2.3. Damage Condition Assessment for Multiple Locations
Taking the damage condition assessment of elements 5, 10, and 16, for example, at nodes 5, 6, 10, 11, 16, and 17 are adopted as the input variables of fuzzy system, and the rule base is formulated and presented in Table 8.

In order to verify the inference ability of the established fuzzy systems of damage condition assessment with multiple locations, test samples are listed in Table 9, and their damage conditions corresponding to severity are also presented.

The damage severities listed in Table 9 can be fuzzified according to the membership function shown in Figure 10. And the linguistic variable corresponding to 12% is “slight damage,” while 27% is “severe damage,” 39% is “catastrophic damage,” 28% is “severe damage,” and 32% is “severe damage.” The damage identification results for these test samples using fuzzy system are demonstrated in Table 10.

As can be seen from Table 10, fuzzy logic system can correctly identify the damage condition for all test samples; it reveals that the fuzzy system possesses favorable inference ability for damage condition assessment of multiple locations.
4. Conclusions
In practice, the baseline model parameters are difficult to obtain and their accuracy will affect the damage identification results. In order to solve this problem, Chebyshev polynomial fitting and fuzzy logicbased damage identification method are proposed for damage identification of bridge without considering baseline data. Numerical simulation of a simply supported beam bridge is used to verify the feasibility of the proposed method.
As for based damage localizing, cases are established for damage identification with single and multiple locations. Modal curvature of damaged structure can be calculated through central difference approximation. By excluding the nodes related to damaged elements, modal curvature of undamaged structure is obtained through sixorder Chebyshev polynomial fitting. Error analysis reveals that it possesses favorable accuracy. Subsequently, parameters can be calculated by use of modal curvature before and after damage. The results indicate that it is effective for damage localizing.
As for fuzzy logic systembased damage condition assessment, is used as input variable and damage severity is output variable. Gaussian and Bell functions are selected as membership functions for input and output variables, respectively. Fuzzy logic systems for damage condition assessment of single and multiple locations are constructed, and the identification results for test samples demonstrate that they possess satisfactory inference ability.
Compared with traditional finite element modeling of bridge before and after damage, the proposed method can realize damage identification depending only on the modal properties of damaged structure. It overcomes the shortcoming that the baseline model is difficult to obtain in practice. But on the other hand a higher requirement of test accuracy for the dynamic response of damaged bridge is put forward, and accurate measurement must be guaranteed in practical application.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors express their appreciation for the financial support of National Natural Science Foundation of China under Grants nos. 51408258, 51378236, and 51278222; China Postdoctoral Science Foundation funded project (2014M560237); Fundamental Research Funds for the Central Universities and Science & Technology Development Program of Jilin Province (20140203002SF).
References
 S. W. Doebling, C. R. Farrar, M. B. Prime et al., “A review of damage identification methods that examine changes in dynamic properties,” Shock and Vibration Digest, vol. 30, pp. 91–105, 1998. View at: Google Scholar
 E. P. Carden and P. Fanning, “Vibration based condition monitoring: a review,” Structural Health Monitoring, vol. 3, no. 4, pp. 355–377, 2004. View at: Publisher Site  Google Scholar
 O. S. Salawu, “Detection of structural damage through changes in frequency: a review,” Engineering Structures, vol. 19, no. 9, pp. 718–723, 1997. View at: Publisher Site  Google Scholar
 Z. R. Lu, J. K. Liu, M. Huang, and W. H. Xu, “Identification of local damages in coupled beam systems from measured dynamic responses,” Journal of Sound and Vibration, vol. 326, no. 12, pp. 177–189, 2009. View at: Publisher Site  Google Scholar
 P. Cawley and R. D. Adams, “The location of defects in structures from measurements of natural frequencies,” The Journal of Strain Analysis for Engineering Design, vol. 14, no. 2, pp. 49–57, 1979. View at: Publisher Site  Google Scholar
 M. A.B. Abdo and M. Hori, “A numerical study of structural damage detection using changes in the rotation of mode shapes,” Journal of Sound and Vibration, vol. 251, no. 2, pp. 227–239, 2002. View at: Publisher Site  Google Scholar
 A. K. Pandey, M. Biswas, and M. M. Samman, “Damage detection from changes in curvature mode shapes,” Journal of Sound and Vibration, vol. 145, no. 2, pp. 321–332, 1991. View at: Publisher Site  Google Scholar
 Z. Y. Shi, S. S. Law, and L. M. Zhang, “Structural damage localization from modal strain energy change,” Journal of Sound and Vibration, vol. 218, no. 5, pp. 825–844, 1998. View at: Publisher Site  Google Scholar
 Z. Zhang and A. E. Aktan, “Application of modal flexibility and its derivatives in structural identification,” Research in Nondestructive Evaluation, vol. 10, no. 1, pp. 43–61, 1998. View at: Publisher Site  Google Scholar
 M. S. Cao and P. Z. Qiao, “Novel Laplacian scheme and multiresolution modal curvatures for structural damage identification,” Mechanical Systems and Signal Processing, vol. 23, no. 4, pp. 1223–1242, 2009. View at: Publisher Site  Google Scholar
 A.M. Yan, G. Kerschen, P. de Boe, and J.C. Golinval, “Structural damage diagnosis under varying environmental conditions—part I: a linear analysis,” Mechanical Systems and Signal Processing, vol. 19, no. 4, pp. 847–864, 2005. View at: Publisher Site  Google Scholar
 V. Meruane and W. Heylen, “Structural damage assessment under varying temperature conditions,” Structural Health Monitoring, vol. 11, no. 3, pp. 345–357, 2012. View at: Publisher Site  Google Scholar
 S. N. Shoukry, G. W. William, B. Downie, and M. Y. Riad, “Effect of moisture and temperature on the mechanical properties of concrete,” Construction and Building Materials, vol. 25, no. 2, pp. 688–696, 2011. View at: Publisher Site  Google Scholar
 C. P. Ratcliffe, “Damage detection using a modified laplacian operator on mode shape data,” Journal of Sound and Vibration, vol. 204, no. 3, pp. 505–517, 1997. View at: Publisher Site  Google Scholar
 D. Wu and S. S. Law, “Damage localization in plate structures from uniform load surface curvature,” Journal of Sound and Vibration, vol. 276, no. 12, pp. 227–244, 2004. View at: Publisher Site  Google Scholar
 M. K. Yoon, D. Heider, J. W. Gillespie Jr., C. P. Ratcliffe, and R. M. Crane, “Local damage detection using the twodimensional gapped smoothing method,” Journal of Sound and Vibration, vol. 279, no. 12, pp. 119–139, 2005. View at: Publisher Site  Google Scholar
 S. C. Zhong and S. O. Oyadiji, “Crack detection in simply supported beams without baseline modal parameters by stationary wavelet transform,” Mechanical Systems and Signal Processing, vol. 21, no. 4, pp. 1853–1884, 2007. View at: Publisher Site  Google Scholar
 M.K. Yoon, D. Heider, J. W. Gillespie Jr., C. P. Ratcliffe, and R. M. Crane, “Local damage detection with the global fitting method using mode shape data in notched beams,” Journal of Nondestructive Evaluation, vol. 28, no. 2, pp. 63–74, 2009. View at: Publisher Site  Google Scholar
 R. Rodríguez, J. A. Escobar, and R. Gómez, “Damage detection in instrumented structures without baseline modal parameters,” Engineering Structures, vol. 32, no. 6, pp. 1715–1722, 2010. View at: Publisher Site  Google Scholar
 J. M. Ko, Z. G. Sun, and Y. Q. Ni, “Multistage identification scheme for detecting damage in cablestayed Kap Shui Mun Bridge,” Engineering Structures, vol. 24, no. 7, pp. 857–868, 2002. View at: Publisher Site  Google Scholar
 M. Sahin and R. A. Shenoi, “Quantification and localisation of damage in beamlike structures by using artificial neural networks with experimental validation,” Engineering Structures, vol. 25, no. 14, pp. 1785–1802, 2003. View at: Publisher Site  Google Scholar
 J. J. Lee, J. W. Lee, J. H. Yi, C. B. Yun, and H. Y. Jung, “Neural networksbased damage detection for bridges considering errors in baseline finite element models,” Journal of Sound and Vibration, vol. 280, no. 3–5, pp. 555–578, 2005. View at: Publisher Site  Google Scholar
 N. Bakhary, H. Hao, and A. J. Deeks, “Damage detection using artificial neural network with consideration of uncertainties,” Engineering Structures, vol. 29, no. 11, pp. 2806–2815, 2007. View at: Publisher Site  Google Scholar
 W. T. Yeung and J. W. Smith, “Damage detection in bridges using neural networks for pattern recognition of vibration signatures,” Engineering Structures, vol. 27, no. 5, pp. 685–698, 2005. View at: Publisher Site  Google Scholar
 H.B. Liu and Y.B. Jiao, “Application of genetic algorithmsupport vector machine (GASVM) for damage identification of bridge,” International Journal of Computational Intelligence and Applications, vol. 10, no. 4, pp. 383–397, 2011. View at: Publisher Site  Google Scholar
 J. M. Benítez, J. L. Castro, and I. Requena, “Are artificial neural networks black boxes?” IEEE Transactions on Neural Networks, vol. 8, no. 5, pp. 1156–1164, 1997. View at: Publisher Site  Google Scholar
 F. S. Mjalli, S. AlAsheh, and H. E. Alfadala, “Use of artificial neural network blackbox modeling for the prediction of wastewater treatment plants performance,” Journal of Environmental Management, vol. 83, no. 3, pp. 329–338, 2007. View at: Publisher Site  Google Scholar
 R. Ganguli, “A fuzzy logic system for ground based structural health monitoring of a helicopter rotor using modal data,” Journal of Intelligent Material Systems and Structures, vol. 12, no. 6, pp. 397–407, 2001. View at: Publisher Site  Google Scholar
 P. M. Pawar and R. Ganguli, “Genetic fuzzy system for online structural health monitoring of composite helicopter rotor blades,” Mechanical Systems and Signal Processing, vol. 21, no. 5, pp. 2212–2236, 2007. View at: Publisher Site  Google Scholar
 M. Chandrashekhar and R. Ganguli, “Uncertainty handling in structural damage detection using fuzzy logic and probabilistic simulation,” Mechanical Systems and Signal Processing, vol. 23, no. 2, pp. 384–404, 2009. View at: Publisher Site  Google Scholar
 Z. Zhao and C. Chen, “A fuzzy system for concrete bridge damage diagnosis,” Computers & Structures, vol. 80, no. 78, pp. 629–641, 2002. View at: Publisher Site  Google Scholar
 M. M. R. Taha and J. Lucero, “Damage identification for structural health monitoring using fuzzy pattern recognition,” Engineering Structures, vol. 27, no. 12, pp. 1774–1783, 2005. View at: Publisher Site  Google Scholar
 D. R. Lynch, Numerical Partial Differential Equations for Environmental Scientists and Engineers, 2005.
 J. H. Mathews and K. D. Fink, Numerical Methods: Using Matlab, Jay Print Pack Private Limited, New Delhi, India, 2005.
 H. B. Liu, Y. B. Jiao, and Y. F. Gong, “A fuzzy logicbased damage identification method for simplysupported bridge using modal shape ratios,” International Journal of Computational Intelligence Systems, vol. 5, no. 4, pp. 627–638, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 YuBo Jiao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.