- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
International Journal of Distributed Sensor Networks
Volume 2013 (2013), Article ID 818967, 12 pages
Health Condition Evaluation of Cable-Stayed Bridge Driven by Dissimilarity Measures of Grouped Cable Forces
1Department of Bridge Engineering, Tongji University, Shanghai 200092, China
2Technology Center of Shanghai Construction Group, Co., Ltd., Shanghai 200431, China
3Shanghai Urban Construction Design Research Institute, Shanghai 200125, China
Received 28 May 2013; Revised 25 August 2013; Accepted 26 August 2013
Academic Editor: Ting-Hua Yi
Copyright © 2013 Dan-Hui Dan 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.
An index system for the health-status evaluation of cable-stayed bridges is presented based on a set of dissimilarity measures (DMs) of the cable force of grouped cables. The DMs and their labels corresponding to health status derived by the serviceability limit states of a bridge are determined by the combination of the finite element method and the influence matrix method. The Monte-Carlo method is used to determine the rational values of thresholds of proposed DMs. By comparing the indices with the thresholds, the holistic health status assessment of a cable-stayed bridge can be easily and reasonably determined. Based on eight years of historical data on the cable forces investigated from a real cable-stayed bridge, the proposed methods are applied to evaluate the health status of this cable-stayed bridge. The results show the validity of the proposed methods.
The relationship between the current internal force state (action effects) and the structure’s limit state (resistance forces) bridges the gap between measureable quantities and structural health status. Thus, the availability of these quantities can facilitate the evaluation of the health status of the long-span cable-stayed bridges . Stayed cables comprise the main structural components of a cable-stayed bridge, and the changes in inclined cable force significantly affect the state of the whole structure . After years of bridge service, deviation attributed to tower-top displacement, foundation settlement, shrinking and creeping, and so on are likely to occur between real and designed cable forces . These deviations are apparently linked to specific and holistic working statuses when the structure is in service. To detect and explore the deviation of grouped cable forces rationally, the changes in the corresponding internal force state of the whole structure can be determined [4, 5]. By associating these changes in the internal force state with the known relationships between such changes and a certain structural limit states, the health status of a cable-stayed bridge structure could be evaluated. Therefore, monitoring the grouped cable forces of long-span cable-stayed bridges and developing a safety evaluation method based on such monitoring data have become important tasks in bridge structural safety evaluation.
A simple comparison method continues to be widely used for the evaluation of cable forces. By computing the relative changes in single-cable force over time to reflect the state of cable forces, such state can be considered normal when the relative change value is less than a certain value (e.g., 5%) . This estimation method has several limitations. On one hand, a rational threshold system has never been built, which makes quantitative evaluation impossible; on the other hand, considering the stress redistribution in the cable plane, the significant changes in single-cable force may not directly affect the safety of the whole structure. However, even when the cable force deviations occurred in each cable is inconsiderable and hereby jumped to a positive conclusion, the holistic cable forces would sometimes fail to meet the safety requirements of the whole structure. Therefore, to evaluate bridge safety based on cable force monitoring data, the whole cable group should be considered simultaneously.
The present synthetic evaluation methods like variable weight method  and gray related degree method  continue to be the most commonly used in evaluation of grouped cable forces. The two methods rely on certain parameters with poor physical significance, thus yielding unsatisfactory results. The former method depends on a balance coefficient representing the degree of correlation of structural components, whereas the latter is constrained by the correlation of the cable forces, where all such methods have no significant relationship with the structures’ working state and thus have the disadvantage of poor interpretability. The statistical methods and pattern recognition methods recently are used in structural status evaluation  and health monitoring [9, 10], with a capacity of dealing with multiple quantities measured by integrated monitoring system or comprehensive bridge detection project synchronously, and so the ability to conduct synthetic evaluation of bridge condition can therefore be expected. For lack of effective channel of communication between patterns and the definition of health status of target structures, the interpretability of these kinds of methodologies still are unsatisfactory. Thus, developing a rational and interpretable technology for the evaluation of the working state of whole structures based on grouped cable forces is necessary.
Based on a comprehensive analysis of grouped cable forces, this paper proposes a rational evaluation method that considers the dissimilarity measure (DM) of grouped cable forces. The threshold system of DMs can then be obtained by using the Monte-Carlo method. Finally, based on the historical measurement data of cable forces of a long-span, cable-stayed bridge located in Guangdong, China, this paper conducts a health evaluation for same bridge and confirms the validity of the evaluation method.
2. Measure of Grouped Cable Forces
A holistic assessment of structural health status can be achieved by analyzing the data samples of grouped cable forces. However, two issues must be addressed: first, the assumption that grouped cable forces are directly related to structural health status should be proven. Second, valid measures of grouped cable forces should be determined to achieve an accurate assessment.
2.1. Correlation of Grouped Cable Forces and Structural Inner Force
As a space-stressed system, a cable-stayed bridge is characterized by a condition wherein loads are borne by the combination of the main girder and the cable tower as well as the stayed cables. Among all loads, vertical loads are primarily balanced by the vertical component of the cable forces. Considering that all cables are anchored to the same tower and main girder, the functionality of a single cable is restricted by that of other cables. All cables are closely related and thus serve a joint function of balancing the loads. Therefore, the inner force of structural components aside from the cables correlates with the holistic cable forces but not with the force of a single cable. Changes in single-cable force, whether present or not, have a negligible effect on the inner force of the whole structure.
The relationship between the distribution of grouped cable forces and the behavior of the structure is significantly affected by factors such as traffic loads, uneven settlement, and temperature. The effects of wind load and nonlinear factors on the internal forces of grouped cables are negligible . Thus, the influence matrix method based on linear theory can be used to describe the effect of grouped cable forces on the behavior of the bridge structure .
Figure 1 shows that if the effect of wind load is disregarded, the relationship between the grouped cable forces and the target mechanical variables of the bridge structure (e.g., external load, inner force, stress, or displacement) could be described as where is a column vector of the grouped cable forces of a finished bridge. , , , and , respectively, denote the cable forces produced by the dead loads, traffic loads, uneven settlement, and temperature; stands for the additional cable forces of the group cables that emerge after reaching the destination through cable force adjustment. is used to denote the passive controlled vectors.
stands for the actions applied to the structure. Index represents the dead loads, represents the traffic loads, represents the uneven settlement, and represents the temperature effect. represents the initial cable force obtained after cable force is adjusted. All elements in are row vectors, and is column vector, which is used to denote the comprehensive active controlled vectors.
stands for the influence matrix and is a constant under the linear hypothesis. , and the submatrices in relate some factors to corresponding cable forces. The referred factors contain dead loads, initial cable force after cable force adjustment, traffic loads, uneven settlement, temperature effect, and corresponding cable forces. The size of first dimension of is equal to the number of grouped cables, and the size of second dimensions is equal to that of .
In the same way, the relationship between the target mechanical variables and the comprehensive active controlled vectors under the ultimate limit state (ULS) and serviceability limit state (SLS) is given by where is column vector of the target mechanical variables, generally refers to typical mechanical variables in typical sections of a cable-stayed bridge, including internal forces, stress and displacement; is the influence matrix between the comprehensive active controlled vectors and the target mechanical variables. The target mechanical variables are typically used in function equations which describe a certain conditions under the ULS and SLS in the design stage. Let represent the function equation or function margin, which is given by where denotes the thresholds or resistance column vectors of the corresponding target mechanical variables.
According to the influence matrix features, (1) could be rewritten as follows: where refers to the grouped cable forces under the combination of the permanent loads; denotes the grouped cable forces under the combination of quasipermanent loads; and stands for the grouped cable forces under the combination of live loads, which has a random value. “” and “”, respectively, represent the average annual temperature effects and the momentary temperature effects of the structure. The former is recognized as a quasipermanent load and has high periodicity and certainty, whereas the latter is characterized by randomness.
The target mechanical variables , which reflect the structure state, can also be rewrite as
Thus, the relationship between the target mechanical variables and the grouped cable forces could be written as the following linear combination:
Equations (1) to (6) show that the state of the inner forces of the structure is closely related to the grouped cable forces. Moreover, the target mechanical variables that reflect the state of the structure can be approximately expressed as a linear combination of grouped cable forces. Thus, through the grouped cable forces, the condition of the target mechanical variables can be obtained; moreover, the function margin can be expressed by the following equation:
The resistance and the effects in (3) are difficult to monitor, but the grouped cable forces in (7) can be measured easily and thus can serve as an agent in function equation. Thus, the original function equation based on resistance and effects is transformed into one that is based on the agent of grouped cable forces, which facilitates to form a structural reliability evaluation frame based on aforementioned monitored measurement agent.
2.2. DMs of Grouped Cable Forces
To determine structure’s states accurately, rational measures should be given for the vectors of grouped cable forces. Grouped cable forces are high dimensional and are thus difficult to deal with. Therefore, low-dimensional measures of the vectors should be introduced to simplify the analysis.
In practical applications, the most common method used for vector measurement is DM. DM determines the dissimilarity among multidimension data samples, which are essentially a group of special dissimilar functions. The differences between two sample vectors can possibly be transformed into scalars. Common DM functions include Euclidian distance, street block distance, Chebyshev distance, Canberra distance, angle separation degree, and correlation coefficient [13, 14]. The two data sample vectors can be set as
Thus, could be written as shown in Table 1.
All DMs could constitute the index system used to evaluate the health status of a structure. When using the given DMs to measure grouped cable forces, observations under specific conditions could be selected as the reference, and together with the rest observations , DMs can be derived via operation on DM operators, that is, . Designed cable forces and finished-bridge cable forces could both be used as .
DMs are used to classify the data points and to identify the boundaries between data points under different categories. Their capability to distinguish such categories could be measured based on the value of the criteria , , and . The given criteria are defined in (9) . A larger value of these criteria denotes a smaller within-class spread and a lager between-class spread of DMs. And a larger between-class spread of DMs results in better classification performance: where is the within-class scatter matrix. is the sample between-class covariance matrix, where is the population covariance matrix, and is the covariance matrix of class . The maximum likelihood estimates of and are, respectively, and : where is the sample mean of the class , which is given by is the sample population mean
3. Threshold Systems of DMs
To achieve an accurate and rational evaluation of the health status of bridge structures, the threshold system of the DMs of grouped cable forces should be determined. From (7), it is known that the thresholds of DMs can be obtained by comparing the real-time target mechanical variables of the significant sections with its corresponding thresholds under two limit states: the ULS and SLS . The solving method is given in the following subsections.
3.1. Definition of Threshold Systems of DMs
The threshold systems of the DMs of grouped cable forces can be obtained by using (3), (6), and (7). Cable-stayed bridges typically have two limit states, and each limit state has upper and lower limits in real conditions. Equation (3) could thus be rewritten as where , denote the lower and upper limit thresholds of the target mechanical variables under ULS, respectively, whereas , denote the same thresholds under SLS, respectively. The four limit thresholds can be obtained based on structure design specifications or through finite element analysis. , stand for the target mechanical variables calculated with grouped cable forces under conditions of ULS and SLS, which can be obtained as follows:
For a cable-stayed bridge, the target mechanical variables that have clear design limits include the displacement of tower top, deflection in the middle of the main span, rotation angle of the girder beam end, and bending moment of the significant girder sections, as well as the bending moment of the tower and pylon-beam junction. These target mechanical variables could be set as an index for evaluation, whereas their design limits can be set as thresholds.
It should be mentioned that if the direct measurement is good enough, it is a more accurate way to monitor and evaluate the status of the bridge. It is reported by Yi et al. that the GPS techniques can help to get the accurate deformation information of bridge; so if there have GPS devices in the monitoring system, it is valuable to use these direct deformation data to evaluate the bridge [15, 16]. However, by now, the price of GPS monitoring system is so expensive that some bridge health monitoring systems do not adopt this technique [17, 18]. In our study, we consider this situation and build our condition evaluation framework on the basis of these two quantities measured indirectly.
Equations (6), (7), and (15) show that the threshold systems of grouped cable force could be indirectly obtained based on the target mechanical variables that reflect the structure state, that is, where subscript “” represents the indexes , , , , , , and , . The thresholds of the DMs of grouped cable forces can then be derived as where is a specific vector of grouped cable forces used as the reference, typically the design cable force or finished-bridge cable force.
Therefore, the evaluation of grouped cable forces could be converted to the comparison of the DMs with its corresponding thresholds,
, which is the threshold of the structural target mechanical variable , is usually a scalar or low-dimensional vector with dimensions that are far less than the number of group cables. Given that (17) does not have a definite solution, it cannot be used to calculate the thresholds. In practical applications, the Monte-Carlo methods can be used to generate mass of samples of grouped cable forces, and then all samples can be labeled into two classes: the safe and unsafe label classes. Subsequently, by distinguishing these labeled DMs samples, the boundary of safe and unsafe samples can be derived, and the threshold system of DMs can be further obtained.
3.2. Flowchart for Threshold Solving Using the Monte-Carlo Method
The initial comprehensive active controlled vectors and the comprehensive influence matrix are obtained through cable force adjustment process firstly, which is operated under the combinations of action effects under SLS, including dead loads, traffic loads, temperature, uneven settlement, and so on. Then, the values of the comprehensive active controlled subvectors corresponding to dead loads, initial cable force, uneven settlement, and global temperature effect are fixed as constants; the values of comprehensive active controlled subvectors, and , corresponding to traffic loads and transient temperature effect, are varied within a rational range randomly, which are then determined using the Monte-Carlo method . The random comprehensive active controlled vectors are then established and used to calculate the random grouped cable force vectors and target mechanical variables via influence matrix method. The DMs of all cable forces vectors are calculated against its references, and (15) is used to determine the sample labels of all cable force vectors by comparing the target mechanical variables with their limits or thresholds. Finally, the threshold system of the DMs of grouped cable forces is determined according to the boundaries between DMs with different labels.
The steps employed by the Monte-Carlo method to determine the threshold system of DMs are shown in Figure 2.
4. Engineering Application
4.1. Engineering Background
Taking a long-span cable-stayed bridge with two towers in China’s Guangdong province as a case study, seen in Figure 3, the feasibility of the method discussed above will be verified in these sections. The bridge has four cable planes, with 40 cables in each plane for a total of 160 cables. Considering that cable force detection has been conducted eight times and that the bridge is not entirely longitudinally symmetrical, the cable forces of the 80 cables of two cable planes at the upstream side are selected as data samples. Totally 16 data samples can be gathered.
All cable forces are regularly detected by the same bridge detection organization and thus the consistency and credibility of the data are guaranteed. Through an analysis of the variation of single cable force, the results show significant similarities. Except for the variations of a few cable forces of nearly ±10% in the previous year, the range of changes of rest cable forces was within −5% to 5%. At present, there are not apparent indications of abnormality occurred during its service history, which means the bridge is in a good condition.
As previously mentioned, the key to apply the Monte-Carlo method to obtain the threshold values of grouped cable force is to bond the class label to a certain limit state. This paper uses Midas Civil Software to build the finite element model (FEM) of the bridge, as shown in Figure 4. In this model, the living loads are arranged at the most unfavorable position. The bending moment of the unilateral bottom section of the tower is chosen as the mechanical variable to calculate the ultimate limit state and the mid-span deflection of girders is taken as the mechanical variable to calculate the serviceability limit state. The bending moment arrives at its ceiling with counterclockwise direction and floor value with clockwise direction when the two limit live load distributions are applied to the girder. While, for serviceability limit state, the two extremes live loads distributions are arranged in this way that the deflection of middle span has exceeded the design specifications or codes both upward and downward. The threshold values of the two mechanical variables are shown in Table 2. From the FEM model of the bridge, the initial active controlled vectors of cable force and the influence matrix can also be obtained. The two limit states are shown in Figure 5.
4.2. Monte-Carlo Statistical Test
To obtain the mass of samples of grouped cable forces and to ensure that the samples can fully reflect the structural behavior of the target bridge under real operation conditions, a traffic load model must be properly built. In order to simplify the calculation, such modeling should differ from lane load modeling in terms of design specification in that it should be randomly distributed on the deck, but it should not stand on any certain vehicle loading manners or any random distribution types of traffic loads.
Therefore, a model of traffic with 80 independent concentrated loads located on the anchorages of the cables and girder are built and applied to the bridge FEM models. Considering the short-term temperature effect, the lower limit of loads should be negative. Thus, the preliminary design is that each concentrated load is evenly distributed within an interval of −70 and 210 tons. To avoid mass repetition of the FEM calculations, the influence matrix method is used to derive the cable forces in the Monte-Carlo simulations operation, under the postulate that the influence matrix can properly represent the complex functional relations among random live loads, grouped cable forces, and target mechanical variables. A total of 180,000 rounds of simulation are conducted, and the same number of cable force samples was obtained. To simplify the illustration, 5000 samples are randomly selected from the population. For every sample, the corresponding values of the target mechanical variables are also derived from the simulation simultaneously.
Figures 6 and 7 show that the grouped cable forces provided by the simulation cover the design cable forces, the region of cable forces of determined by the SLS, as well as the historical cable forces and the scope of these samples approached to the limit determined by ULS. It can be concluded that the parameters used in the simulation have good representativeness and can reflect the possible conditions of real cable forces. Similarly, the simulation values of target mechanical variables overspread the feasible engineering region, and the sample set contains the full information on structural behavior.
4.3. Comparison of DMs and Determination of Threshold System
It is reasonable to designate the finished-bridge state as the reference, which implies that the target bridge is in a health status at the very beginning of its service history; choosing the grouped cable forces on phase of finished bridge as standard, the DM of the samples of cable forces are calculated. The corresponding class labels can then be determined based on the comparison between the simulation values and the SLS threshold values of the target mechanical variables. Therefore, when using DMs to describe grouped cable force samples, each DM is related to a class label.
According to (9), we can assess the classification capacity of DMs based on their corresponding class labels, and as shown in Table 3, it can be seen from that the best DM is Canberra distance, followed by street block distance, and Euclidean distance. However, the values of Chebyshev distance, angular separation, and correlation coefficient are extremely small and are thus inappropriate for classifications.
Figure 8 shows that two clear distinctions exist between the points of DMs with two different class labels in Figure 8(a) Canberra distance, Figure 8(b) street block distance, and Figure 8(c) Euclidean distance, which divides the whole measurement space into three regions: completely safe region, completely unsafe region, and mixed region. For Figure 8(d) Chebyshev distance, the completely unsafe line is clear, but the completely safe line is difficult to identify. For Figure 8(e) angular separation and Figure 8(f) correlation coefficient, the points with safe and unsafe labels are fully mixed, and almost the entire measurement space is a mixed area. Thus, the two DMs lose the capability of classification.
Based on the formulas of DMs, a lesser similarity of the samples denotes larger values of Euclidean distance, street block distance, Chebyshev distance, and Canberra distance but smaller values of angular separation and correlation coefficient. Combined with the criterion, the maximum of DMs with safe labels and the minimum of DMs with unsafe labels are chosen as the thresholds, as shown in Table 4.
4.4. Health Status Evaluation
To validate the rationality and validity of the given threshold system, the DMs in practical application are derived by samples of historical measured cable forces and finished-bridge cable forces. With eight years’ measurement results and considering 160 cables divided into two samples (80 upstream cables as the first sample and 80 downstream cables as the second), the number of samples of historical measured cable forces comes up to 16. Thus, 16 results can be obtained for each DM, which are drawn in scatter plots in Figure 9. The green points are DMs of the measures at the finished-bridge state (given that they are compared with themselves, DMs equal zero), whereas the red points are the maximum DMs, which are the closest points to the safe boundaries. The results show that all data are located at the safe region of the corresponding measurement space. This demonstrates that both the safety requirements and serviceability requirements of the bridge are met. This deduction from these three figures is consistent with the current service status of the bridge.
Figures 9(a)–9(c) show that whether Canberra distance, street block distance, or Euclidean distance, the maximum values of DMs occur in the most recent monitoring year (2010). With the increasing lifespan of the bridge, each DM value gradually increases and approximates the upper limit of the safe region. The tendency discovered by DM history curves demonstrates that the overall performance of the bridge is gradually degrading and the disquieting departure from baseline status is increasing continuously, which should call the attention of the bridge owner.
This paper proposed an evaluation method for the health status of a cable-stayed bridge based on a set of DMs of grouped cables forces. By mapping the relationship between grouped cable forces and health status and building the measure space of grouped cable forces, the evaluation of structure health status is realized. Likewise, the structure’s health status can be grasped in a macro sense. The main research conclusions are as follows.
(1) The status of grouped cable forces corresponds to the specific inner force state of a bridge. Thus, a holistic state evaluation of bridge health can be accomplished by evaluating the state of grouped cable forces as the agent, and the measure of grouped cable forces, which are DMs, can further service for this evaluation.
(2) Based on the Monte-Carlo method and combined with the threshold of target mechanical variables under the bridge limit state, a threshold system of the DMs of grouped cable forces can be established.
(3) Among all the DMs, Canberra distance, street block distance, and Euclidean distance exhibit better classification performance and can thus be used for similar cases.
(4) With increasing service time, the DMs values of the grouped cable forces continuously increase, which indicates that the performance of the bridge structure degrades annually.
Based on the finding of this paper, historical data on cable forces can be used to establish a safety evaluation system for cable-stayed bridges. Likewise, the data can also be used to facilitate the functions of online early warning and real-time evaluation of bridges based on the fact that the health monitoring system of bridges can acquire sufficient real-time data on cable forces.
However, this paper does not consider the influence of nonlinear factors in calculating the influence matrix. Relevant research considering more other factors will be conducted in the future.
This paper was sponsored by the project of the National Key Technology R&D Program in the 12th Five Year Plan of China (Grant no. 2012baj11b01), the National Natural Science Foundation of China (Grant no. 50978196), the Key State Laboratories Freedom Research Project (Grant no. sldrce09-d-01), the Fundamental Research Funds for the Central Universities, the State Meteorological Administration Special Funds of Meteorological Industry Research (Grant no. 201306102), and Research of the Provincial Department of the Ministry of Education of Guangdong Province (2008b090500222).
- D. F. Wisniewski, J. R. Casas, and M. Ghosn, “Codes for safety assessment of existing bridges—current state and further development,” Structural Engineering International, vol. 22, no. 4, pp. 552–561, 2012.
- H. M. Irvine, Cable Structures, Dover Publications, New York, NY, USA, 1991.
- F. N. Catbas and A. E. Aktan, “Condition and damage assessment: issues and some promising indices,” Journal of Structural Engineering, vol. 128, no. 8, pp. 1026–1036, 2002.
- N. Wang, C. O'Malley, B. R. Ellingwood, and A.-H. Zureick, “Bridge rating using system reliability assessment. I: assessment and verification by load testing,” Journal of Bridge Engineering, vol. 16, no. 6, pp. 854–862, 2011.
- N. Wang, B. R. Ellingwood, and A.-H. Zureick, “Bridge rating using system reliability assessment. II: improvements to bridge rating practices,” Journal of Bridge Engineering, vol. 16, no. 6, pp. 863–871, 2011.
- H. Wu, Y. Xiang, and Y. Li, Research on Safety Evaluation of Large Span Cable Stayed Bridge, China Municipal Engineering, 2006, Chinese.
- X. Shang and Y. Xu, “Safety-based cables' condition evaluation of cable stayed bridge with grey theory,” Journal of Chang'an University, vol. 24, no. 1, p. 52, 2004.
- M. Gul and F. N. Catbas, “Identification of structural changes by using statistical pattern recognition,” in Proceedings of the 6th International Workshop on Structural Health Monitoring, Stanford, 2007.
- H. Sohn, C. R. Farrar, N. F. Hunter, and K. Worden, “Structural health monitoring using statistical pattern recognition techniques,” Journal of Dynamic Systems, Measurement and Control, Transactions of the ASME, vol. 123, no. 4, pp. 706–711, 2001.
- M. Gul and F. Necati Catbas, “Statistical pattern recognition for Structural Health Monitoring using time series modeling: theory and experimental verifications,” Mechanical Systems and Signal Processing, vol. 23, no. 7, pp. 2192–2204, 2009.
- X. Hua and G. Xu, Nonlinear Analysis of Bridge Structure, vol. 6, People's Communication Press, 1997.
- R. Xiao and H. Xiang, “Influence matrix optimization method on cable force of cable stayed bridge,” Journal of Tongji University, no. 3, pp. 235–240, 1998.
- P. C. Mahalanobis, “On the generalized distance in statistics,” in Proceedings of the National Institute of Sciences of India, pp. 49–55, 1936.
- R. W. Andrew, Statistical Pattern Recognition, 2nd edition, 2004.
- T. Yi, H. Li, and M. Gu, “Full-scale measurements of dynamic response of suspension bridge subjected to environmental loads using GPS technology,” Science China Technological Sciences, vol. 53, no. 2, pp. 469–479, 2010.
- T. H. Yi, H. N. Li, and M. Gu, “Recent research and applications of GPS based technology for bridge health monitoring,” Science China Technological Sciences, vol. 53, no. 10, pp. 2597–2610, 2010.
- T.-H. Yi, H.-N. Li, and M. Gu, “Experimental assessment of high-rate GPS receivers for deformation monitoring of bridge,” Measurement, vol. 46, no. 1, pp. 420–432, 2013.
- T.-H. Yi, H.-N. Li, and M. Gu, “Wavelet based multi-step filtering method for bridge health monitoring using GPS and accelerometer,” Smart Structures And Systems, An International Journal, vol. 11, no. 4, pp. 331–348, 2013.