- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Mathematical Problems in Engineering
Volume 2013 (2013), Article ID 969021, 13 pages
Nonparametric Change Point Diagnosis Method of Concrete Dam Crack Behavior Abnormality
1College of Water Resources and Architectural Engineering, Northwest A&F University, Yangling 712100, China
2State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China
3College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China
Received 2 June 2013; Accepted 14 July 2013
Academic Editor: Xu Zhang
Copyright © 2013 Zhanchao Li 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.
The study on diagnosis method of concrete crack behavior abnormality has always been a hot spot and difficulty in the safety monitoring field of hydraulic structure. Based on the performance of concrete dam crack behavior abnormality in parametric statistical model and nonparametric statistical model, the internal relation between concrete dam crack behavior abnormality and statistical change point theory is deeply analyzed from the model structure instability of parametric statistical model and change of sequence distribution law of nonparametric statistical model. On this basis, through the reduction of change point problem, the establishment of basic nonparametric change point model, and asymptotic analysis on test method of basic change point problem, the nonparametric change point diagnosis method of concrete dam crack behavior abnormality is created in consideration of the situation that in practice concrete dam crack behavior may have more abnormality points. And the nonparametric change point diagnosis method of concrete dam crack behavior abnormality is used in the actual project, demonstrating the effectiveness and scientific reasonableness of the method established. Meanwhile, the nonparametric change point diagnosis method of concrete dam crack behavior abnormality has a complete theoretical basis and strong practicality with a broad application prospect in actual project.
China has the most water-conservancy and hydropower projects currently in the world. At the beginning of the 21st century, China entered the peak period of dam construction. There is a group of 320 m-level ultrahigh concrete arc dams completed or under construction in Southwest China. Concrete dam crack and its influence on the structure safety have always been a hot spot and difficulty in the discipline of hydraulic structure. Currently, the study is focused on whether the macrocrack of concrete structure will be expanded under load, which is an important issue for the evaluation of safety performance of concrete structure working with crack . In dam safety monitoring field, as the masses have higher and higher requirements to the safety of water conversancy projects, the problem of concrete dam crack has been more and more highlighted. It has also gradually become a hot spot and difficulty in the study on safety of hydraulic structure to study the transformation of crack behavior in concrete dam due to crack expansion and even its abnormality to cause influence on the safety of concrete structure. Such problem is defined as concrete dam crack behavior abnormality in this paper, and the emphasis is put on the study of diagnosis method of crack behavior abnormality, so as to diagnose the crack abnormality timely and rapidly.
For the hydraulic concrete structure working with crack, Li et al. , in combination with mutation theory created the grey cusp mutation model of hydraulic concrete structure crack abnormality. Bao and Yu , based on grey differential fitting modeling, created the grey state equation of concrete dam crack, analyzed the relation between the stability of grey state equation and crack subcritical expansion, and proposed the grey system method to distinguish and judge the subcritical expansion of crack. Given that concrete dam crack system is an open, dissipated, and complex nonlinear dynamical system, Gu et al.  proposed dynamical fuzzy cross-correlation factor index method to preliminarily explore the dynamic abnormality diagnosis of concrete dam crack behavior. The current study plays an important promotion role in expanding the study thought and promoting discipline crossing, and development and forecasts that applied mathematics theory will have a strong development potential and broad development space in this field.
Deep analysis shows that all abnormality diagnosis methods above are based on the sequence of concrete dam crack mouth opening displacement (CMOD), but, according to numerous researches and investigations, the author finds that whether the CMOD can be taken as the basis for diagnosis of crack behavior abnormality is not deeply studied in the existing literatures, and most of the literatures are based on the fact that CMOD generally is the only monitoring item of concrete dam crack, causing that the existing concrete dam crack abnormality diagnosis method is lack of solid theoretical foundation. Based on this, the author, according to the fracture mechanics fictitious crack model and linear asymptotic superposition assumptions, creates the double crack tip opening displacement criterion of concrete crack behavior abnormality, providing a theoretical basis for the use of CMOD for crack behavior abnormality diagnosis in actual project . Meanwhile, the basis for the existing concrete dam crack behavior abnormality diagnosis method is the parametric statistical model of safety monitoring on concrete dam crack , which is established on the basis of this statistical model and in combination with certain or some applied mathematics theories and emphasizes the study from the aspect of information analysis, and it is still required to have a deep analysis on basic issues such as scientificalness and reasonableness of the abnormality diagnosis method.
The parametric statistical model of safety monitoring on concrete dam crack takes hydraulic pressure, temperature, and timeliness as the main influential factors, in which the theoretical basis of crack timeliness component expression is the creep theory of concrete and the component of cyclic load of temperature and so forth is fitted with periodic function, and it expresses the hydraulic component with cubic or quartic polynomial continuous function of water head based on the theory of engineering mechanics . So the parametric statistical model of safety monitoring on concrete dam crack is widely applied in actual project as it can reflect the evolution law of crack behavior truthfully in a certain extent. However, when the abnormality diagnosis of crack behavior is based on the parametric statistical model of safety monitoring on concrete dam crack, as the functional relationship between the influential factors of crack and the crack effect size is set in advance and the stochastic disturbance of model needs to meet the normal consumption, the effect of crack behavior abnormality diagnosis is not reasonable and comprehensive enough. On this basis, closely combining the hot issue in the latest statistics, change point theory, based on the analysis of change point characteristics of concrete dam crack behavior and through the reduction of change point issue, establishment of basic nonparametric change point model, and progressive analysis of basic change point test method as well as the study on multichange point diagnosis method of crack abnormality, the author creates the nonparametric change point diagnosis method of concrete dam crack behavior abnormality in this paper and applies it in the actual project, so as to inspect the scientific reasonableness of nonparametric change point diagnosis method of concrete dam crack behavior abnormality.
2. Analysis on Change Point Characteristics of Concrete Dam Crack Behavior Abnormality
Crack behavior abnormality is representation that the crack evolves from stable state to instable state. The abnormity of concrete dam crack behavior is presented as the instability of model structure in crack monitoring model, that is, the dynamic relation between environment size and crack effect size changes suddenly in several time points for some reason or changes continuously and slowly in a certain time interval. For example, according to the existing crack behavior abnormality diagnosis method, the hydraulic component and temperature component in the CMOD are recoverable, while the timeliness component is unrecoverable. Thus, most of the abnormality diagnosis methods of crack behavior are against the timeliness component of CMOD. When the timeliness component meets some conditions, it is thought that the abnormality occurs to the crack behavior, for example, the grey cusp mutation model method . The general form of such methods can be described as first, the general expression of timeliness component of CMOD is defined as where is the influence factor related with timeliness and is the -dimensional regression parameter of timeliness component; for example, the expression of timeliness component of CMOD is generally .
Then the abnormality of crack behavior in the mode is presented as the regression parameter changes in some time points that is, in time point , the regression coefficient determined by samples is significantly different from the regression coefficient determined by samples statistically.
The above crack behavior abnormality is diagnosed based on the timeliness component. However, generally the hydraulic pressure and temperature components can also reflect the abnormality of crack behavior. For example, under the same water level and temperature conditions, when the concrete dam crack depth increases, the annual variation amplitude of hydraulic pressure component and temperature component in CMOD will be changed. In the statistical model of CMOD, it is reflected that the regression coefficients of hydraulic pressure component and temperature component are changed; when this change is significant statistically, the crack behavior is abnormal.
For the two given CMOD sample sequences, in order to test whether the crack is abnormal when generating the two sample sequences, the statistical two-sample analysis method can be adopted. In statistics, two-sample analysis method is mainly to analyze whether the probability distribution characteristics of two samples are significantly different statistically, for example, the Chow test . For concrete dam crack, assume that, in its monitoring model, the endogenous variable is , exogenous variables are , their sequence observation values () are from two times with different environments, in which the first () observation values are from a certain historical period and the rest ones are from another period. At this moment, if the data of the pre- period, and those of the post- period are, respectively, used to establish two linear models between and :
Obviously, if the intercepts of and are different or the coefficients of a certain exogenous variable are different, the structural relation between and is instable. On the contrary, if the intercepts and coefficients of all exogenous variable are the same, the two formulas are the same model in nature. It is called that there exists a stable structural relation between and . On this basis, it is only needed to establish one model between all the observation values between and :
It can be seen that the study on structural stability of the model can be converted into the statistical test issue of model parameter in the following form:
In order to test whether the different hypothesis on the two models is true, start from the original hypothesis that the two regression equations are the same, and then test whether the hypothesis can be refused. If the ordinary least squares method is used to estimate the models (2) and (3), respectively, as there is no limitation condition to the model coefficients, the sum of residual sum of squares of the two equations is used to represent the unconditional residual sum of squares, that is, , in which and , respectively, represent the residual sums of squares of models (2) and (3), they are respectively subject to the distribution of degrees of freedom of and , and is subject to the distribution of degree of freedom of . Then the ordinary least squares method is used to evaluate the model (4) and calculate the residual sum of squares with condition limit, which is subject to the distribution of degree of freedom . If we assume that is true, then will not exceed too much. On this basis, the hypothesis can be judged whether it is true according to whether the residual sum of squares is significant; the test statistics constructed is
Obviously, the above test statistics is subject to the distribution of first degree of freedom and second degree of freedom . If the structure between and is stable, the parameters estimated in the models (2), (3), and (4) should be approximate. Thus the difference between and is tiny; at this moment, tends to be the smaller value; otherwise, will be larger. Based on this, the structural stability of the model can be judged accordingly for the significance level given; if , is accepted; at this moment, it can be thought that the model structure is stable; that is, no abnormality occurs to the crack behavior at this moment; on the contrary, if , is refused, and thus the structure between and is instable; that is, abnormality occurs to the crack behavior at this moment.
It can be seen that when the abnormality point of concrete dam crack behavior is known, the test analysis method of two-sample issue can be adopted to diagnose whether abnormality occurs to the crack behavior. But in practice, it is difficult to know the moment when abnormality occurs to the crack behavior. When the abnormality point is unknown, the issue on crack behavior abnormality diagnosis is converted into change point issue.
The main study of change point issue is to test and analyze whether the statistical characteristics of a given time sequence are significantly changed at an unknown moment. Statistically, change point means that some statistical characteristics of the sequence or process change at an unknown moment . Assuming a sample sequence , the corresponding distribution functions are . If there exists , which makes that differs greatly from in some characteristics (such as mean value and position parameter etc.), then is called a change point of the sequence. More generally, if are divided into groups, with sample distribution in each group relatively stable, while there are mutations at , then are the change points of the sequence.
Change point is ubiquitous in nature and society, which reflects the process from quantitative change to qualitative change of things. The issue of change point, whether in theoretical exploration or actual data analysis, is very important, that is, to analyze whether the potential mechanism is changed, when it is changed, and how it is changed based on a group of data; if these problems are ignored, the accuracy of estimation and forecast will be reduced.
For the monitoring sequence of an effect size of concrete dam crack, for example, the CMOD monitoring sequence , whether abnormality occurs to the crack behavior can be expressed as the original hypothesis is true or alternative hypothesis , is true; that is, whether the distribution of CMOD sequence is changed somewhat, while the distribution change forms are diversified in theoretical exploration and practice. The distribution of observation value is changed according to a certain law and is changed to another law at an unknown moment, and this moment is the change point. Or the moment of specific change cannot be found out at all, and only the distribution law of observation is changed progressively just at the beginning. It is also possible that the sample distribution relies on a spatial parameter, while this mutation or progressive change relies on the position or interface in the space. Another possibility is that the distribution law of samples will be changed with change of system condition or an element in the system, and then the mutation or progressive change will rely on a specific nontime and space variable. Generally, the “structural change” in the parametric statistical model is often represented as the change of specific parameter . Another “structural change” is not completely decided by the change of some parameters, such issues are of diversified forms and their models are also relatively complex. The change point study on concrete dam crack behavior abnormality in this paper is mainly to detect whether the statistical behavior of the crack is significantly fluctuated, especially to detect whether the crack fluctuation exceeds its safety control range. When the crack is changed qualitatively (mainly exceeding the warning line of crack control), it is hoped that early warning can be released in time to avoid safety accident of the dam.
3. Construction of Concrete Dam Crack Behavior Change Point Diagnosis Method
It can be known from the above analysis on characteristics of concrete dam crack behavior abnormality that whether abnormality occurs to the crack behavior can be diagnosed with change point theory. According to the prior information depended, the existing change point theory diagnosis methods can be mainly classified into  parametric, semiparametric, and nonparametric diagnosis methods. In the study of change point issues, due to the uniqueness of nonparametric statistics; that is, no background of sample distribution is required, in recent years, the nonparametric methods in study on change point issue have been very active and their contents are also very rich. This paper is focused on the nonparametric change point study of concrete dam CMOD sequence, so as to establish the nonparametric diagnosis method of concrete dam crack behavior abnormality.
3.1. Reduction of Change Point Issue
Any characteristic of the random sequence is determined by its probability distribution, that is, by the probability measure in a functional space. Therefore, the change of any characteristic of the random sequence is the change of its probability distribution generally. Based on this, it can be thought that the random sequence with the change point is “bonded” by the sequence fragments with different probability distributions , while the any-probability distribution of random sequence is decided by a limited-dimension distribution set meeting the conditions of consistency . Therefore, any change of the probability characteristics is the product of some change of any limited-dimension distribution. Consider the following single change point model: assume and are two strictly stationary random sequences, command , , and define
Assume and consider the sequence , in which where is any sequence. Then it is thought that is the sequence with a change point; that is, its dimension distribution function has a change point. Here and , respectively, represent the moment when the change point begins and ends, and interval is the length of change point transient process. Assume , and , where is the sample capacity. If , it is thought that the change point is a mutation change point; if , it is thought that the change point is a progressive change point. Similarly, corresponding multichange point issues can be defined.
There are varieties of change point issues, and the different probability characteristics of stationary random sequence can be defined as the corresponding change point issues. Thus, only with the help of all finite-dimension distributions can the random sequence be fully described. Therefore, for the purpose of theoretical analysis and practical application of the change point issue, it is necessary to refine a basis change point type from so complex and various change point issues, so that other types of change point issues can be reduced to such basic type conveniently. It is thought that sequence with a change point is a basic type of change point if the mathematical expectation of sequence changes, that is, where and is any sequence.
The method to diagnose the change point of the formula (9) is a basic change point diagnosis rule. Generally, in this basic type of change point, other characteristic parameters of the random sequence also may be changed; here other parameters are called redundant parameters. Meanwhile, study shows that , as long as enough diagnosis sequence , is introduced; the test of any distribution function parameter change in sequence can be reduced to a basic type of change point test; that is, the reduction of change point issue makes it possible to adopt the same change point test method to test any information involved in the characteristics of random sequence change point as long as the basic change point test rule is analyzed.
3.2. Basic Nonparametric Change Point Model
Generally, the a posteriori change point problem can be regarded as the estimation problem of parameter vector. Assuming , is an unknown parameter vector and satisfies where and and are two prior constants. Assuming , is a determined parameter family, two types of random sequence family in the probability space should be considered:
Assume when , . The a posteriori change point problem is an estimate of change point time , of random sequence . Therefore, the change point test problem translates into the consistent estimation problem of parameter vector according to the observation of . For example, as to the change point problem with a mutation change point, when , , (), , where is an indicator function, the following statistics should be considered:
In the interval , use a sequence of points obtained to construct a continuous random process by liner interpolation, the random point in the set is the estimate of parameter , where
Under certain conditions, for random , exists. Meanwhile, when certain conditions are satisfied, . Assuming and are two strict and stationary independence random sequences, when certain conditions are satisfied, sequence , is of absolute convergence. And if , , on the continuity interval : where
where and are two standard and independent Wiener processes.
When , , and , can be simplified as
The above method is constructed in consideration of the basic type of change point. Meanwhile, the change point test method of basic change point problem mentioned above can be spread to the change point test problem of any probability characteristics.
3.3. Asymptotic Analysis of Basic Change Point Test Method
Through the asymptotic analysis of first and second types of probability errors and estimated probability error, it is able to choose the most optimal adjustment parameter in the family to meet the requirement of change point test quality in practice. Assuming is an observation of a random sequence, the following form of statistics family should be considered: where and .
When , is the statistics. When , is the maximum likelihood statistics constructed by Hinkly  for Gaussian sequence. When , is the nonparametric statistics adopted by Deshayes and Picard [13, 14]. Assume the observed sequence is from the following generating scheme: where , , , , and , is Gaussian white noise, and is the standard Wiener process. For continuous time sequence, the statistics can be expressed as where . In the interval of , , so the maximum value point of process on the interval of is the estimate of change point .
The first type of error probability can be expressed as
The second type of error probability can be expressed as
And estimated error probability can be expressed as
For the first type of error probability, the following equation exists: where . It is obvious that function is of monotone increasing on the interval of . So the asymptotic optimal method to reduce the first type of error probability is to assume and then .
For the second type of error probability, under certain conditions, the following formula exists:
It can be seen that the second type of error probability satisfies
Therefore, the asymptotic optimal method to reduce the second type of error probability is to assume , so .
For estimated error probability, any exists: where , , is a standard normal distribution function, , and . When ,
Namely, when , where satisfies ; namely, when , the estimated error probability is the minimum. Therefore, selection appropriate parameter can improve the quality of change point test.
The nonparametric diagnosis method of crack behavior abnormality of concrete dam is constructed theoretically by the reduction of any change point problems, construction of basic nonparametric change point model, and asymptotic analysis of basic change point test method. And its basic ideas are first, to construct corresponding diagnosis sequence based on the initial time sequence of CMOD according to the reduction guidelines of change point problem; then, to construct statistics sequence according to the formula (11); finally, to diagnose whether the change point exists or not according to the formula (13) under the given level of significance test.
However, in practice, for the given CMOD sequence of concrete dam, it is hard to ensure that there is only one abnormality point upon crack abnormality, so the number of crack abnormality point is unknown. Therefore, the multiple change point diagnosis method of crack behavior abnormality of concrete dam is built combing with the asymptotic property of the statistics family below.
3.4. Multiple Change Point Diagnosis Method of Crack Behavior Abnormality
For many existing change point test methods, if assuming the diagnosis sequence can be described by a mathematical model of limited parameter set, the methods belong to the scope of parametric change point test method [15–17]. When the phenomenological model of study object is known, the parametric method is effective. However, when the phenomenological model has model error during setting, the result obtained by parametric change point test method is doubtful. Therefore, this paper constructs corresponding nonparametric change point test method to diagnose the crack behavior abnormality of concrete dam which has less demand of prior information of diagnosis sequence; that is to say, the nonparametric method does not require the prior information of probability distribution of diagnosis sequence . Furthermore, for the given level of significance test, this method can almost test the abnormality characteristics of random crack behavior. Assuming a random sequence observation , satisfies where and , .
Assume and , where and . And assuming is independent identically distributed random variables sequence, and satisfies. The model (28) is the corresponding basic nonparametric change point model in case of multiple change points existed. For the actual CMOD in concrete dam, it can be turned into the form of this model by constructing appropriate diagnosis sequence. Therefore, this paper analyzes the multiple change point of this model and builds corresponding multiple change point diagnosis method of crack behavior abnormality of concrete dam.
The asymptotic analysis of basic change point test method shows that when , the first type of error probability is minimum; namely, the false alarm probability is minimum; when , the second type of error probability is minimum; namely, the false dismissal probability is minimum; when , the estimated error probability of change point is minimum; namely, the probability of change point estimated value far away from the actual change point position is minimum. Therefore, selecting different value in different change point test phase can improve the precision and credibility of the estimated result of change point . The basic steps of multiple abnormality point diagnosis method of crack behavior abnormality of concrete dam are as follows. (1) Construction of the diagnosis sequence from the original CMOD sequence: construct diagnosis sequence from original signal to make the change of various characteristics of original signal convert into the change of diagnosis sequence mathematical expectation. (2) Homogeneity test of diagnosis sequence: calculate the statistics , from diagnosis sequence. If is true, this diagnosis sequence should be deemed as homogeneity; that is to say, there is no abnormality of crack behavior. If it is not workable, this diagnosis sequence should be deemed as nonhomogeneity; that is to say, this sequence is of crack behavior abnormality point. Go to the next step.
In the formula (29), is the threshold value of crack behavior abnormality point test which relies on the specified false alarm probability, sample size, sample variance, and correlation coefficient of sample where the false alarm probability is specified in advance, the smaller this parameter value, the more abnormality point will be of missing test. In this phase, to not miss testing any possible abnormality point, the false alarm probability should be bigger. (3) Initial estimation of crack behavior abnormality point: the point that makes the statistics global maximum can be used as the estimation of first abnormality point. This point divides the diagnosis sequence into two subsample set including and where is a smaller value relying on the sample size and the maximum value in statistics. Then, test the homogeneity of each subsample and by the method in step and estimate the abnormality point by the method in step . Repeat above steps until each sample in the subsample set is homogeneous. A sequence of initial estimation of abnormality points is obtained by steps and where is the initial estimation number of abnormality points. (4) Eliminating suspicious crack behavior abnormality point: redivide the diagnosis sequence to make each subsample set only contain an abnormality point of initial estimation, that is to say
Conduct homogeneity test of each subsample set according to step . However, the false alarm probability should be smaller. If a certain subsample set is tested to be homogeneous, the corresponding abnormality points in the subsample set should be eliminated. (5) Final estimation of crack behavior abnormality point: after step , for the remaining subsample set with a size of , calculate the statistics , and the corresponding point of maximum value of the statistics is the final estimation of the abnormality point. In the above abnormality point test method, the determination of threshold is the key. To calculate the threshold , the homogeneity assumption of sequence should be workable. Then, according to the functional limit theorem of random sequence , the following formula exists under certain conditions where makes the sequence meet and is standard Brownian bridge:
When false alarm probability is given, threshold can be calculated by the following formula:
It can be seen from the above analysis that the threshold relies on and is unknown. To solve this problem, Piryatinska et al.  adopts the following method to estimate this threshold. First, when the false alarm probability is given, assuming , the initial estimation of threshold can be obtained according to the formulae (34) and (35). Then, divide the diagnosis sequence into two subsample sets including and by steps and , calculate the mean value of each subsample set, and deduct the corresponding mean value from the subsample set to form a new sample set by merging. Finally, calculate the empiric variance of the sample set as the estimation of . So the estimation of threshold is
It can be seen from the basic principle of the method that when there is no or only one abnormality point in the diagnosis sequence, the estimation is more precise and the estimation is much reliable. However, when there are multiple abnormality points in the diagnosis sequence, the error of estimation obtained by this method is larger and the estimation is suspicious. Meanwhile, it is assumed initially that can affect the threshold estimation to a certain degree.
In consideration of the problems existed, this paper improves this method and its basic ideas are under an assumed , estimate the abnormality points in the diagnosis sequence, then use empiric variance estimation according to similar method, and obtain corresponding threshold estimation according to the formula (36). Then analyze the abnormality points of diagnosis sequence again by this threshold, in this cycle, until the abnormality point estimation for two times is the same, estimation satisfies , and threshold estimation is . This improvement can eliminate the influence of multiple abnormality points and assumed on threshold to a great extent.
4. Project Case
Conduct abnormality diagnosis of the crack near 105 m elevation at the downstream face of a gravity arch dam based on the nonparametric change point diagnosis method of crack behavior abnormality of concrete dam. This dam is a concrete gravity arch dam of concentric variable radius with a crest elevation of 126.3 m and maximum dam height of 76.3 m and has 28 dam sections from left to right. It takes 12 years from dam pouring to substantial completion and is constructed by three phases. During Phase II concrete pouring, as the level rising speed is fast and interval of concrete lift is short, the contraction distortion of Phase II concrete is strongly constrained by Phase I concrete to make the top of Phase I concrete have cracks from 5# dam block to 28# dam block with a length of more than 300 m. So, the dam blocks of 14#20# were of cross-seam joint bar grouting in 1973 and of epoxy resin grouting in 1987. This paper analyzes the 1680 CMOD monitoring data of joint meter inbuilt in 18# dam block from October 10, 1974 to December 18, 2006 (as shown in Figure 1), so as to diagnose the abnormality conditions of crack behavior.
Upon the analysis of given CMOD sequence of concrete dam by adopting nonparametric change point diagnosis criteria, construct diagnosis sequence according to the original sequence first. And the CMOD means and amplitudes are the two essential characteristics where the mean value presents the base point of change of CMOD, namely, the average level of CMOD, and the amplitudes represent the amplitude of variation of crack under load in one cycle. Therefore, when there is change point in CMOD means or amplitudes, the crack of corresponding moment is deemed as abnormality. Based on this, this paper analyzes the crack behavior abnormality according to the diagnosis sequence of CMOD means, amplitudes, and correlation functions.
The diagnosis sequence based on CMOD means can be expressed as
The diagnosis sequence based on CMOD amplitudes can be expressed as where is time delay.
The diagnosis sequence based on CMOD correlation functions can be expressed as where is time delay.
This paper analyzes the diagnosis sequence according to the nonparametric change point diagnosis method of crack behavior abnormality built where Figure 5 is the diagnosis process of crack behavior abnormality based on the diagnosis sequence of CMOD means which is similar with the abnormality diagnosis process of diagnosis sequence based on CMOD amplitudes and correlation functions and not repeated here. Tables 1, 2, and 3 are the abnormality diagnosis results of diagnosis sequence based on CMOD means, amplitudes, and correlation functions.
For the diagnosis sequence of CMOD means, it can be seen from Figure 2 that the CMOD means have two obvious change points during the period of January 15, 1976 to November 25, 1976 and September 26, 1985 to August 28, 1986, that is to say, the CMOD means have abnormality. Meanwhile, the CMOD means have increased since 1987; that is to say, the CMOD means increase gradually. However, it can be seen from Table 1 that upon the abnormality analysis of the diagnosis sequence of whole means, the abnormality point of CMOD means diagnosed first occurred on August 28, 1986 which divided the diagnosis sequence of means into two subsequences. Conduct abnormality diagnosis of each subsequence and the abnormality points of two subsequences occurred on May 5, 1977 and October 20, 1994. The final estimation time of the three abnormality points was on August 7, 1986, September 15, 1977 and October 6, 1994 which was similar to the observation of abnormality point in the diagnosis sequence of means. The abnormality diagnosis results of diagnosis sequence based on CMOD amplitudes and correlation functions are also similar.
As analyzed above, the nonparametric change point diagnosis method of crack abnormality built in this paper can diagnose the abnormality condition of crack behavior. The diagnosis result not only contains obvious abnormality point of crack behavior, but also is of good diagnosis effect on the abnormality point of crack behavior that is not easy to find.
(1)Whether the concrete dam crack expands or not is an important subject of safety performance evaluation of concrete dam, and the study on the abnormality diagnosis of crack behavior of concrete dam has become the hot spot and difficult problems in hydraulic structure safety. However, as the functional relationship between influencing factor and effect size of crack is set in advance and the random noise item of model should satisfy the normality assumption, the abnormality diagnosis result of crack behavior based on the safety monitoring parametric statistics model of concrete dam crack is not reasonable enough and comprehensive. (2)The internal relationship between the crack behavior abnormality of concrete dam and statistics change point theory is deeply analyzed to make the abnormality diagnosis method of crack behavior of concrete dam based on change point theory have complete theoretical basis from the perspective of structural instability of parametric statistics model and change in the distribution law of sequence of nonparametric statistics model based on the representation of crack behavior abnormality of concrete dam in parametric statistics model and nonparametric statistics model. (3)The test of parameters change of any distribution functions is reduced as the basic type and method of change point test of sequence mathematical expectation by the reduction of change point problem and establishment of basic nonparametric change point model. The optimal adjustment parameter in the change point test statistics is selected by asymptotic analysis of basic change point test method to satisfy the requirement of change point test quality in practice. And the multiple change point diagnosis method of crack behavior abnormality of concrete dam is built considering that there may be multiple abnormality points in the crack behavior of concrete dam in practice. Finally, the nonparametric change point diagnosis method of crack behavior abnormality of concrete dam is built by integrating the above study results. (4)The nonparametric diagnosis method of crack behavior abnormality of concrete dam can be applied in practical engineering. The diagnosis result contains obvious abnormality point of crack behavior and has good diagnosis effect of crack behavior abnormality point that is not easy to find, demonstrating the effectiveness and scientific rationality of built method. Meanwhile, the nonparametric change point diagnosis method of crack behavior abnormality of concrete dam is of complete theoretic basis, strong practicability, and wide application prospect in practical engineering. (5)In accordance with the concept of dam safety monitoring, the research in this paper and many other researches on the abnormality diagnosis of crack behavior in concrete dam belong to the a posteriori category; that is, the crack monitoring data is regarded as given and the sample size is known and statistical central limit theory is used to diagnose crack behavior abnormality. However, in order to achieve the real-time monitoring and early warning of dam safety in practice, the concern about crack behavior abnormality of concrete dam is whether the current crack behavior is abnormal and whether the former monitoring model of crack safety is suitable for the current status, which corresponds to online diagnosis of crack behavior abnormality. The crack monitoring data is regarded as changing and the sample size is unknown in the online diagnosis. This will be the main research content of the authors and the research group in the future.
This work was supported by the Open Foundation of State Key Laboratory of Hydrology-Water Resources (Grant no. 2012490211), the doctoral scientific research fund of Northwest A&F University in 2012 (Grant no. 2012BSJJ007), the Fundamental Research Funds for the Central Universities (Grant no. QN2013046), and the Key Project of the National Natural Science Foundation of China (Grant no. 51139001).
- D. H. Huang, Y. P. Song, and G. F. Zhao, “Study on crack characteristics of macrocrack concrete structure,” China Civil Engineering Journal, vol. 33, no. 2, pp. 63–66, 2000.
- X. H. Li, H. Z. Xu, C. S. Gu, et al., “Crack abnormality diagnosis based on wavelet analysis and cusp catastrophic model,” Journal of Hohai University (Natural Science Edition), vol. 33, no. 3, pp. 301–305, 2005.
- T. F. Bao and H. Yu, “Detection of subcritical crack propagation for concrete dams,” Science in China Series E, vol. 52, no. 12, pp. 3654–3660, 2009.
- C. S. Gu, Z. C. Li, and B. Xu, “Abnormality diagnosis of cracks in the concrete dam based on dynamical structure mutation,” Science China Technological Sciences, vol. 54, no. 7, pp. 1930–1939, 2011.
- Z. C. Li, C. S. Gu, and Z. R. Wu, “Abnormality diagnosis of cracks in the concrete based on double crack tip opening displacement criterion,” Science China Technological Sciences, vol. 56, no. 8, pp. 1915–1928, 2013.
- C. S. Gu and Z. R. Wu, Theory and Method for Security Monitoring and Control of Dam and Dam Foundation and Their Application, Hohai University Press, Nanjing, China, 2006.
- Z. R. Wu and C. S. Gu, Theory for Security Monitoring and Control of Hydraulic Structure and Its Application, Higher Education Press, Beijing, China, 2003.
- X. H. Li, Study on Evolvement Rule and Abnormality Diagnosis Method of Major Hydraulic Concrete Structure Crack, Hohai University, Nanjing, China, 2003.
- G. C. Chow, “Test of equality between sets of coefficients in two linear regressions,” Econometrika, vol. 28, pp. 591–605, 1960.
- B. E. Brodsky and B. S. Darkhovsky, Nonparametric Methods in Change Point Problems, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1993.
- A. N. Shiryaev, Probability, Nauka, Moscow, Russia, 1980.
- D. V. Hinkly, “Time-ordered classification,” Boimetrika, vol. 59, pp. 509–523, 1972.
- J. Deshayes and D. Picard, “Off-line statistical analysis of change-point models using non parametric and likelihood methods,” in Detection of Abrupt Changes in Signals and Dynamical Systems, vol. 77 of Lecture Notes in Control and Information Sciences, pp. 103–168, 1986.
- J. Deshayes and D. Picard, “Testing for a change-point in a statistical model,” Tech. Rep., Universite de Paris-Sud, Orsay, France, 1981.
- S. Shaban, “Change-point problem and two-phase regression: an annotated bibliography,” International Statistical Review, vol. 48, pp. 83–93, 1980.
- P. Krishnaiah and B. Miao, “Review about estimation of change points,” in Handbook of Statistics, P. Krishnaiah and C. R. Rao, Eds., vol. 7, pp. 375–402, Elsevier, Amsterdam, The Netherlands, 1988.
- P. Perron, “Dealing with structural breaks,” Working Papers Series in Macroeconomics WP2005-017, Department of Economics, Boston University, 2005.
- M. Hollander and D. A. Wolfe, Nonparametric Statistical Methods, Wiley, New York, NY, USA, 1973.
- B. Brodsky, B. Darkhovsky, A. Y. Kaplan, and S. Shishkin, “A nonparametric method for the segmentation of the EEG,” Computer Methods and Programs in Biomedicine, vol. 60, no. 2, pp. 93–106, 1999.
- M. Peligrad, “Invariance principles for mixing sequences of random variables,” Annals of Probability, vol. 10, no. 4, pp. 968–981, 1982.
- A. Piryatinska, G. Terdik, W. A. Woyczynski, K. A. Loparo, M. S. Scher, and A. Zlotnik, “Automated detection of neonate EEG sleep stages,” Computer Methods and Programs in Biomedicine, vol. 95, no. 1, pp. 31–46, 2009.