Research Article  Open Access
Prediction Methods for Routine Maintenance Costs of a Reinforced Concrete Beam Bridge Based on Panel Data
Abstract
In order to make informed decisions on routine maintenance of bridges of expressways, the hierarchical regression analysis method was used to quantify factors influencing routine maintenance cost. Two calculation models for routine maintenance cost based on linear regression and timeseries analysis were proposed. The results indicate that the logarithm of the historical routine maintenance cost is the dependent variable and the bridge age is the independent variable. The linear regression analysis was used to obtain a cost prediction model for routine maintenance of a beam bridge, which was combined with the quantity and price, and verified by a physical engineering example. In order to cope with the cost changes and future demands brought about by the emergence of new maintenance technologies, the timeseries analysis method was used to obtain a model to predict the engineering quantities for the routine maintenance of a bridge based on standardized minor repair engineering quantities. Taking into account the actual cost of the minor repair project as well as the timeseries analysis’ sample size demands, the annual engineering quantity was randomly decomposed into four quarterly data quantities, and the timeseries analysis result was verified by physical engineering. These results can improve the calculation accuracy of the routine maintenance costs of reinforced concrete beam bridges. Furthermore, it can have a certain application value for improving the cost measurement module of bridge maintenance management systems.
1. Introduction
Maintenance cost planning is an important component of highway infrastructure asset management. In the classification system of highway maintenance engineering in China, routine maintenance has the characteristics of regularity, repeatability, and timeliness, as compared with medium maintenance, heavy maintenance, preventive maintenance, emergency maintenance, or special maintenance. It is the basic work of highway maintenance management and plays an important role in prolonging the periodic overhaul life and reducing the cost of road users [1, 2]. The measurement and calculation of annual maintenance expenditures (AMEX) is an important module of an asset management system [3–8]. Figure 1 shows the valuable uses of AMEX prediction models [9–13].
The objects of the asset management of highway infrastructure involve subgrade, pavement, bridge and culvert, and tunnel and traffic facilities. The structural characteristics and working principles of each asset are quite different, so each should be treated differently in the process of management and maintenance. Many scholars have devoted themselves to the study and development of the maintenance cost budget of the asset management of highway infrastructure. A large number of reports have focused on the calculation of routine maintenance. Most of the research methods are based on regression analysis, and the dependent variable is the routine maintenance cost [9, 14–24] or routine maintenance quantity [25]. According to the choice of independent variables, the approaches can be summarized into the following four categories:(1)Using performance as an independent variable, for instance, AlMansour and Sinha [17] determined the regression equation of road annual routine maintenance cost and the pavement service performance index. Lu and Tolliver [18] determined the regression equation of annual routine maintenance cost and the international roughness index.(2)Taking the natural environment, traffic environment, and road age as independent variables, the representative results are as follows: Hagood [19] concluded that the independent variables of five models of pavement routine maintenance cost in Nevada included the length of road section, area, road age, altitude, climate, temperature, cumulative equivalent axle load, and the cost of the previous year. The California pavement maintenance cost model studied by Gibby et al. [20] includes the annual average daily traffic (AADT) of heavy vehicles, AADT of passenger cars, road age, temperature, presence of shoulders, presence of bridges, and road area. Volovski et al. [21] considered the road age, AADT, average annual rainfall, length of the road, number of lanes, nature of the road (new or reconstructed pavement), and other factors when analyzing the cost of road maintenance in Indiana, USA. Qiu et al. [22] et al. constructed an indicator system of highway routine maintenance rates and considered six calculation variables such as AADT, annual average daily weight of vehicles, routine maintenance engineering amount, and special maintenance engineering amount.(3)Considering both the road performance and the road technical conditions, Woldemariam et al. [9] used the ANN method to predict the cost of pavement maintenance. The input variables were annual average temperature, annual average rainfall, total network size, IRI, and truck vehiclekilometers traveled. Wu [23] reported on the quantitative relationship between the annual maintenance cost of highway asphalt pavement in Zhejiang Province, China, as well as the pavement quality index (PQI), traffic volume, truck ratio, and operating life. Yue [25] summarized the influencing factors of damage maintenance amounts in the routine maintenance of asphalt pavement using pavementengineering quality grade, width, thickness, performance index, service life after conversion, annual average daily heavy traffic (AADT), annual average rainfall, and annual average temperature difference.(4)Directly providing the routine maintenance cost intervals without considering any independent variables: for instance, the US Highway Economic Requirements System State (HERSST) calculated that the average daily maintenance cost of a road unit length varies from $285 to $7830 per lanekilometer. Ola and Ugiomoh [24] obtained the criteria of the annual maintenance cost of a road, which states that the cost of a paved twolane rural main road is from $870/km to $1730/km and that the cost of a paved fourlane main road is from $1682/km to $6743/km, and finally that that of a twolane main road without paving is $703/km to $1407/km. It is obvious that the notion of a cost interval is affected by various factors.
Research on the maintenance cost of bridge assets involves the preventive maintenance opportunity of a bridge based on a costbenefit analysis [26], the maintenance cost model based on lifecycle cost [3–8, 27–32], the maintenance cost model [28] included in the development of bridge maintenance management systems [5, 6, 33, 34], and the research specifically on the prediction of bridge maintenance cost [10, 31, 35–39]. The factors influencing the maintenance cost of bridges are the same as those influencing the technical status of bridges [40–43]. Studies in the literature [44] have reported on analyses of the influence of bridge superstructure form, completion time, load grade, bridge length, and other factors. Other studies in the literature [45] have considered the influence of environment, bridge age, maintenance measures, and other factors when evaluating the technical status of bridges. The AADT, structure location, and other factors were taken into account when verifying the adaptability of a bridge deck in Maryland, USA. In the relevant literature [46–48], the influence of environmental and traffic factors on bridges was considered, respectively. In the research literature of Hegazy et al. [34], factors such as bridge deck type (steel bridge or concrete bridge), year, highway type (interstate highway or other), average daily traffic (ADT), width, length, and technical status were considered. In the report of Ali et al. [10], climatic factors (including annual average temperature and annual average rainfall), average daily truck traffic (ADTT), technical status, total bridge deck area of each state, average bridge deck area of each state, and other factors were considered.
The existing research can be summarized as follows:(1)Transportation infrastructure maintenance has obvious regional characteristics. Researchers analyzed the problem with a specific area as their research object [19–23].(2)Most researchers concentrate on pavement maintenance costs because this constitutes the largest portion of the road maintenance cost, but they concentrate less on the study of bridge maintenance cost.(3)Factors affecting the asset condition should be taken into account in the calculation of the maintenance cost of pavement and bridges. Different scholars consider different influencing factors in different periods. In particular, traffic influencing factors, such as AADT, ESAL, heavytruck AADT, and heavytruck ratio, are considered. They all choose a certain factor directly and do not provide a basis for selecting the corresponding factor.(4)There are more studies to be found on repair expenditure [10], preservation expenditure [14, 35, 36], or maintenance expenditure specifically [9], but fewer studies on routine maintenance [15, 16, 23, 25] or daily maintenance costs.
We believe that there are significant differences between different types of maintenanceimplementation occasion and implementation content, while routine maintenance engineering has the characteristics of regularity, and the corresponding cost must be paid every year. Therefore, there is both theoretical and practical significance to focus study on maintenance engineering. This paper draws lessons from an idea that is based on the research on the routine cost of pavement maintenance. The data in this paper are derived from the historical routine maintenance cost and details of reinforced concrete (RC) beam bridges since the opening of the 11 expressways in Shaanxi Province. With the aim of setting research boundary conditions and adopting quantitative analysis methods, the influencing factors of routine maintenance cost of this type of bridge asset are selected and determined. Based on these, this paper puts forward two calculation methods: the regression model with routine maintenance cost as the dependent variable, and the prediction model with routine maintenance details as the dependent variable. Their respective applicable conditions are delineated, which are verified by physical engineering.
In this paper, the routine maintenance tasks for reinforced concrete beam bridges include the maintenance of drainage facilities, expansion joints, abutment, the maintenance of guardrails, antithrowing nets, expansion devices, and drainage facilities. Among these items, the maintenance of the guardrail involves painting for concrete guardrails and rustremoving paint for the handrails. The maintenance of the antithrowing mesh refers to replacing the pillars and bases of the antithrowing mesh. The maintenance of expansion devices involves repairing the concrete, replacing rubber strips, and welding the steel of expansion joints. Maintenance of drainage facilities involves replacing the centralized drainage pipe of a bridge, replacing the fourway base of PVC and repairing the cover of the drainage hole.
2. Materials and Methods
2.1. Setting Boundary Conditions
Expressway routine maintenance costs are affected by a number of factors. Quantitative studies on their cost investment are required to limit the boundary conditions, which include the following.
2.1.1. Delimiting “Special Service Supply Level”
In the current maintenance management system of China, different maintenance types influence each other along the same highway. The cost of small repairs and their quality of implementation affects the timing of preventive maintenance. The implementation effect of preventive maintenance affects the time required for large and medium repairs and whether timely preventive maintenance and overhaul will affect routine maintenance expenditures. To maintain an expressway’s operating service level, the expressway operating company sets a minimum acceptable level for the performance indicators of the expressway assets. This paper defines this level as the level of maintenanceservice supply. The maintenance costs reflected by different levels of maintenanceservice supply levels vary. Therefore, it is necessary to set the boundary conditions of the maintenanceservice supply level to ensure that bridge assets have the same maintenanceservice level. In this case, the research is valuable in practical applications.
2.1.2. Climate and Environmental Factors Normalized to Regional Factors
In the course of evaluating the assets of expressway bridges, it is important to note that they are influenced by natural environmental factors. This paper normalizes these natural environmental factors into regional factor variables. In other words, bridge assets located in the same area are situated in the same range of temperature variation, precipitation, and altitude. In this research, bridge assets located in the same climate region are treated as a group of analysis objects.
2.1.3. Traffic Factor Indicators
Traffic volume factors can be reflected by the AADT, heavytruck ratio, ESAL, and cumulative ESAL. In the literature, different scholars have selected various different indicators of traffic volume.
To determine the traffic factors affecting the cost of routine maintenance of bridge assets scientifically and reasonably, this paper divides the traffic volume indicators into the following categories: AADT (hereinafter referred to as traffic volume indicator 1), AADT + heavytruck volume (hereinafter referred to as traffic volume indicator 2), AADT + heavytruck ratio (hereinafter referred to as traffic volume indicator 3), ESAL (hereinafter referred to as traffic volume indicator 4), and (hereinafter referred to as traffic volume indicator 5),
Using the above, the appropriate traffic indicator is determined by quantitative analysis.
2.1.4. Years of Usage
Due to the different opening years of each highway project, the analysis converts the different opening years into bridge age and uses the same number of years of traffic in the same area as a set of sample data.
2.1.5. Road Segmentation
This paper takes the expressway in Shaanxi Province, China, as the analysis object. According to the management system of Shaanxi Province, a specific branch is responsible for the operational management of each expressway. The maintenance budget is also formulated on the basis of each individual expressway. This paper uses the historic cost data of each operating branch of the expressway for the underlying analysis. As a result, this work considers a highway managed by each branch as a specific analysis object. The maintenance management level of each branch is influenced by the maintenance policy of the central office, and each branch is under the unified jurisdiction of the group company. Therefore, the maintenanceservice supply level of each branch is the same.
2.2. Study Methodology
2.2.1. Applicability of Multivariate Regression Analysis Methods
Multiple regression analysis can be used for predictive regression as well as interpreted regression [49]. Predictive regression can be used to predict the trend of the dependent variables, where the interpreted regression mainly focuses on the relationship of the variables. This interpreted regression meets the needs of this paper. The influencing factors of this study have continuous variables, such as years of usage and traffic volume, and also include categorical variables, such as regional factors. When there are categorical variables in a multiple regression analysis, the method requires hierarchical regression analysis [50], which satisfies the requirements of analyzing multiple influencing factors. At the same time, the available mathematical tools can be used to judge the degree of influence of each factor on the routine maintenance of each asset.
2.2.2. Preprocessing of Variables
(1) Preprocessing of Dependent Variables. In the analysis, the total routine maintenance investment for each bridge is taken as the dependent variable (interpreted variable). The total cost includes the sum of the daily cleaning and routine maintenance under the boundary conditions of an identical level of maintenanceservice supply. Because the routine maintenance cost of each asset is affected by the size of the bridge [9, 10, 14, 34] and the year of implementation, this research adopts the following methods:(1)In order to eliminate the impact of the asset scales on the routine maintenance costs, the annual actual cost of each asset is standardized by equation (1).(2)In order to eliminate the impact of price factors of manpower, materials, and machinery in different opening years, the routine maintenance cost data for different years are converted into the cost value of the base year (2015) with price indices. In other words, in equation (1) is the routine maintenance cost that is converted into the base year:where is the routine maintenance cost per lanekilometer in the t year of road bridge i (yuan/lane·km), is the routine maintenance cost in the tyear of road bridge i (yuan), is the number of lanes of road i (usually 4, 6, or 8 lanes), and is length of the road i in the lane k (km).
In introductory econometrics, when the dependent variable is positive and the unit is yuan, the logarithm is given to the dependent variable first, using as the dependent variable. The logarithm can eliminate or reduce the skewness of the mathematical distribution of the dependent variable and can narrow its range. This brings it closer to the classical, linear model assumption instead of the direct use of as the dependent variable. In this study, the calculation results of equation (1) are taken as the dependent variable after taking the logarithm.
(2) Preprocessing of Independent Variables. The independent variables used in this paper include traffic volume, regional factors, and age. They also include factors of scale involving the length, width, and area of bridges [10, 34] because the historical data are accumulated for each road, and it is impossible to identify the cost of the bridges of different scales on the road; hence, such factors are not included in the independent variables. The regional factor in the independent variable is the category variable, and the other two factors are continuous variables. The most important hypothesis in multiple regression analysis is the linear relationship hypothesis. This requires that the independent variables in the model must be continuous, and further indicates that the categorical variables are not suitable for multiple regression analysis [50]. In reality, the categorical variables are extremely important. In order to enable the use of the regional factors in the regression analysis with other continuous variables, it is necessary to virtualize the categorical variables. This requires that the categorical variables be converted into dummy variables first.
The key step to virtualize the categorical variables is the recoding of them. Usually, the recoding method employs virtual coding, which converts the categorical variables into 0 and 1 dichotomous variables. After conversion, the dummy variables can be taken as source data into the multiple regression equation for regression analysis, similar to that for the other continuous variables. When the type of the categorical variable is k, the virtualized categorical variables need to be split into k1 virtual variables. These can then be used in the hierarchical regression analysis function module in a statistical software package to input the entire set of k1 dummy variables into the multiple regression analysis.
There are three regional variables analyzed in this paper as categorical variables: Shanbei, Guanzhong, and Shannan. When virtualizing these three categorical variables, one first must choose a reference group. This study regards Shannan as the reference group, which enables the division of the three categorical variables into two dummy variables.
(3) Independent Variable Classification and Order. When there are categorical variables in a multiple regression analysis, the best choice is the hierarchical regression analysis [50]. In the hierarchical regression analysis, one first groups the independent variables and then brings the independent variables into the multiple regression equations in the order of grouping. In this paper, the hierarchical group is classified as follows: the regional variable is used as group 1; the age is taken as group 2; and the traffic volume indicator is used as group 3. It should be noted that, for traffic volume indicators 1 through 4 and for traffic volume indicator 5, only a single traffic volume indicator can be input at a given time for hierarchical regression analysis.
2.3. Parametric Analysis Using Multivariate Regression
2.3.1. Multivariate Regression Analysis Models and Indicators
The multivariate hierarchical regression model is shown in the following equation:where is the routine maintenance cost per lanekilometer in the tyear of road bridge i (yuan/lane·km), is the regional factor variable (Shanbei, Guanzhong, and Shannan), is the age, and is the traffic volume indicator.
Regression coefficient indicates the interpretation degree of the corresponding independent variable for the dependent variable . The regression coefficient is a nonstandardized statistical parameter with units. Although it can reflect the number of influences of the corresponding variables of the independent variable, it cannot convincingly indicate the importance of each independent variable according to their sizes. Therefore, it cannot be used for comparison purposes. If one wants to compare the regression coefficients and get a more intuitive interpretation of each independent variable against the dependent variables, all variables in equation (2) can be standardized by calculating Z scores; the intercept term a disappears. Equation (2) then becomes equation (3), and it is then very straightforward to judge the interpretation degree of independent variables for various dependent variables. in equation (3) is defined by equation (4):where is the standard deviation of the independent variable and is the standard deviation of dependent variables after taking the logarithm.
2.3.2. Analysis Content and Criteria for Assessing the Influencing Factors
The purpose of the regression analysis in this paper is to explore the interpretation degree of independent variables against the dependent variables. Therefore, the overall explanatory power, R^{2}, of the model is reported first. The statistical significance of R^{2} is explained by the test data of the F test. In the process of hierarchical regression analysis, when a variable enters the group to participate in the regression analysis before other variables, the first variable plays the role of the control variable in the group, while the subsequent variables are added. represents the additional explanatory power that the group is dedicated under the influence of the previous set of variables. When is as large as possible and reaches statistical significance, it is considered that the addition of the group can effectively enhance the explanatory power of the model.
The research focuses next on the standardized coefficient. In Tables 1 and 2, each independent variable has a corresponding standardized regression coefficient, the β coefficient, and a significance test result. Standardizing the regression coefficient can unify the units of each independent and dependent variable and reduce the error caused by conflicting units, making the results more accurate.


For the β coefficient, the larger the absolute value of the ttest is, the smaller the p value becomes. The p value indicates the significance of the ttest. In statistics, if , then the coefficient test is generally considered to be significant, indicating that the independent variable can effectively predict the variation of the dependent variable. In the case, the possibility of an incorrect conclusion is 5% and that of a correct one is 95%.
In the coefficient table of regression analysis results, the β coefficient can be positive or negative. This reflects that the standard deviation—led by the doubling of an independent variable—causes the standard deviation of the dependent variable to increase by the multiple (or decrease by the fraction) of the absolute value of β, respectively. A larger absolute value of β of an independent variable means that the independent variable affects the dependent variable to a greater degree. Therefore, in this paper, we only focus on the absolute value of β.
Based on the regression analysis results of all the models and the actual scenario, this paper addresses whether an independent variable can be the main influencing factor of routine maintenance cost. The result is as follows:(1)When β > 0.30 and the p value is statistically significant, the explanatory variable is the main influencing factor.(2)When β > 0.30 and the p value are not statistically significant, the result might be caused by insufficient sample data. However, the practical significance of the variable cannot be ignored, so the explanatory variable is also considered to be the main influencing factor.(3)When β is approximately 0.20, it is necessary to assess the importance of the explanatory variable based on its practical significance. If it is of great practical significance, then it can be considered as the main influencing factor; if not, then it cannot be considered.
3. Results and Discussion
3.1. Analysis and Discussion of the Model Results
This paper analyzes 11 expressways managed by Shaanxi Provincial Communication Construction Group. The longest operation duration among them is 15 years, and the shortest is 3 years. According to the regional characteristics of Shaanxi Province, 11 highways are divided into those of three regions, namely, Northern, central, and Southern Shaanxi, which are taken as three independent variables of region. During research, we collected the routine beam bridge maintenance cost of all 11 expressways for each operation duration. Preprocessing of Dependent Variables shows the price index adjustment and normalization of the data. We take its natural logarithm as the dependent variable. Independent Variable Classification and Order classifies the independent variables into different groups and defines their time series for hierarchical regression analysis. The analysis results are shown in Tables 1 and 2, and Figure 2 shows the corresponding standardized regression coefficients.
From Tables 1 and 2, the following could be obtained:(1)The explanatory capability of the overall model of all independent variables is visibly high, indicating that they can explain the change in the dependent variable.(2)From the perspective of the individual explanatory capability of each hierarchy, the explanatory capability of Group 1 is relatively low, while that of Group 2 is significantly high ( = 0.375, ). This demonstrates that taking into consideration the roadoperation duration can effectively improve the explanatory capability of the model.(3)The increment regression result of the traffic factor in Group 3 shows that the increment of any traffic factor has no statistical significance.
As shown in Figure 2, the values of the β coefficient of operation duration are high in each group, remaining larger than 0.6, and the ttest results are of high statistical significance. However, for other independent invariables, the β coefficients are small and the ttest results are of no statistical significance.
Therefore, we conclude the analysis on the influencing factors of routine beam bridge maintenance costs as follows:(1)Regarding the regional factor, there is no obvious difference in the routine beam bridge maintenance cost between Northern, central, and Southern Shaanxi.(2)The roadoperation duration has a strong explanatory capability and is statistically significant. Therefore, it can be an important influencing factor.(3)For the five traffic factors, although the model reflects that they have great overall explanatory capability, the values of the β coefficient are small and the statistical significances are low. As a result, the traffic volume indicators are considered to be unimportant influencing factors.
3.2. Calculation and Verification of Routine Maintenance Cost of the RC Beam Bridge
According to the analysis results of the influencing factors of routine RC beam bridge maintenance costs, the regression model of routine maintenance cost prediction can be established in two manners: first, routine maintenance cost is set as the dependent variable, and roadoperation duration is set as the independent variable. Through regression analysis, their functional relation can be determined and then used to predict the annual maintenance cost in the future. Second, based on the analysis results of influencing factors and the timeseries forecasting method, the amount of routine maintenance is taken as the dependent variable to predict its annual value in the future. In this way, according to the current maintenance quota, the cost can be calculated.
3.2.1. Prediction and Verification of Routine Maintenance Cost Based on Linear Regression
Here, we continue to use the data in the quantitative analysis on influencing factors, namely, the historical maintenance cost data (for 2015) after price index adjustment, standardization, and logarithm processing. It is set as the independent variable and expressed by y. The road operation duration (beam bridge age) is taken as the independent variable. The final regression formula is
To verify the formula accuracy, we compare the predicted value of the regression formula with the actual routine beam bridge maintenance costs in 2016 and 2017, as shown in Table 3. Data in the table are normalized to 2017 prices.

Nonparametric testing does not require samples to obey an assumption of normal distribution, so it is usually applied to the data test of a smallscale sample, classified data, and an ordinal number. As part of the doublecorrelation nonparametric test, the Wilcoxon signedrank test meets the requirement of the consistency check between the predicted value and actual value in this study.
The Wilcoxon signedrank test is a nonparametric hypothesis testing method. It checks the data in two data sets in pairs (a paired difference test) to determine whether the two data sets come from the same distributed collection.
The Wilcoxon signedrank test, as a type of doublecorrelation sample test, is based on correlation detection. It analyzes the correlation between two groups of data first. When the correlation reaches a certain degree, the Wilcoxon signedrank test is carried out. The test does not require the samples to obey the assumption of normal distribution and is usually applied to test the data of smallscale samples.
In analysis of the result of the correlation coefficient, apart from the coefficient, the SIG value (namely, the p value), which indicates the significance, also exists.
The original hypothesis of the Wilcoxon signedrank test was that there is no significant difference between the predicted value and the actual value, and the alternative hypothesis was that they are different.
The test result shows that the predicted value is highly correlated with the actual value and is statistically significant, indicating the feasibility of the prediction formula.
3.2.2. Prediction and Verification of the Engineering Quantity of Routine Maintenance Based on Time Series
The cost prediction model built on nonlinear regression and logarithm processing makes a forecast on the basis of historical data. However, the reference value of historical data is severely affected owing to the impact of new materials, new technologies, and new equipment on cost. This is especially true for the period during which policy changed from business tax to VAT; the items of cost changed significantly. Therefore, considering the separation of quantity and price, the prediction of the engineering quantity of routine maintenance can not only eliminate the influence of new technology on cost forecasts but can also exempt the impact of new taxes and feepricing policies.
According to the analysis’ results of influencing factors, the engineering quantity of RC routine beam bridge maintenance is affected by various natural and accidental factors which will gradually show a certain time characteristic in the long term. Timeseries analysis is usually used to study and analyze a series of data arranged in an order. It can reveal the inherent law of the data and make shortterm predictions. Therefore, this paper uses timeseries prediction with the assistance of the SAS analysis software to study and make shortterm predictions on the engineering quantity of RC routine beam bridge maintenance.
The modeling steps for timeseries analysis are shown in Figure 3.
(1) Data Preprocessing. Data preprocessing in timeseries analysis includes standardization of engineering quantities, elimination of outliers, and supplementation of missing values.
In order to eliminate the influence of the scale of the bridge structure on the routine engineering quantity data of the RC beam bridge and according to the relationship between details and the length and width of the bridge, unitization of routine maintenance in each detail is as follows. The routine maintenance related to the length and width of bridges is divided by the product of the bridge length and the number of lanes, for example, the replacement of expansion joints. The amount of routine maintenance related only to the length of the bridge is divided by the length of the bridge, for example, the maintenance of steel handrails.
In timeseries analysis, due to the limitation of actual sample size, the representative value of one year’s routine maintenance is divided into four seasons’ data, which not only meets the sample size requirements predicted by the time series but also conforms to the actual situation of quarterly payment and settlement in the project.
Aggregating the actual data, this paper describes the performance of timeseries analysis on three routine maintenance items: replacement of expansion joints, concrete repair of expansion joints, and maintenance of steel handrails.
(2) Data Verification. ① Stationarity Test. The run length method is used to test the stationarity of the sample data, and the test results are used as the basis for selecting the appropriate model. Results are shown in Table 4.
 
is timeseries length; is the number of values greater than the average of the time series; is the number of values less than the average of time series; is the total number of run lengths of the sequence; ; ; is the test statistics; ; . 
Given the significance level α = 0.05 and that the values of in Table 4 are all less than 1.96, the selected sample sequences all meet the stationarity requirements.
② Pure Randomness Test. The purpose of the pure randomness test is to judge the interdependence of sample sequence data according to the hypothesis significance test, so as to determine whether there is research value in modeling analysis. If the sample sequence is a nonpure random sequence, it shows that the sequence has the value of further analysis. The test results are shown in Table 5.
 
is the significance of the chisquared test in a pure random test. 
After the pure randomness test, the significance of the chisquared test for each sample sequence in Table 5 is less than 0.05, so all sequences are nonpure randomness sequences.
(3) Model Identification and Rating. For the tested data, the sequence diagram is drawn to determine whether there are singularities and, if so, make corrections.
The BIC criterion is a method for determining the order of the best criterion function in a timeseries analysis. According to the BIC minimum criterion and combining with the autocorrelation function graph and partial autocorrelation function graph of the sample sequence, the model is selected and the order of the model is determined. Through optimization and adjustment, it was finally determined that the timeseries models of the three routine maintenance cases were all firstorder models.
(4) Parameter Estimation and Testing. The maximum likelihood method was used to estimate the model parameters, and the significance test was performed on the parameters. The results are shown in Table 6.
 
is means of a sequence of samples; is the coefficient of the autoregressive parameter term; is significance of the ttest for parameter estimation. 
Through the test, the ttest significance of all model parameters given in Table 6 is less than 0.05, and the autocorrelation test results of the residuals all meet the requirements, so the model passes the test.
(5) TimeSeries Prediction. According to the timeseries model, the various routine maintenance quantities of 8 seasons in total in 2016 and 2017 were predicted, and the reliability of the model was verified with the actual values. The results are shown in Table 7. The difference in the table is the difference between the forecasted value and the actual values.

It can be seen from Table 7 that the differences between the predicted values of the three engineering quantities and the actual values are small. In particular, the predicted value of the concrete repair of expansion joints in 2017 is in good agreement with the actual value, and the absolute error is only 0.02. Thus, it can be seen that the prediction results of the timeseries model is in accordance with engineering practice, and the deviation from the actual data is within a reasonable range, which indicates that the model is feasible.
4. Conclusions
Reasonable determination of expressway routine maintenance cost is an important issue in highway asset management. Previous studies have mainly focused on the routine maintenance of pavement assets. In order to explore the cost model of the routine maintenance of bridge assets, in this study, we applied hierarchical regression analysis to interpretative regression analysis. This paper relies upon the historic data of 518.254 km of RC beam bridges from 11 expressways managed by the Shaanxi Province Transportation Group in China. Two kinds of routine maintenance cost prediction models of RC beam bridge assets based on the combination of quantity and price, and on the separation of quantity and price, were obtained.
The research basis of this paper is the preprocessing of the dependent and independent variables. Preprocessing mainly involves taking the logarithm, virtualizing the categorytype independent variables, classifying the independent variable groups, and ordering the entrance hierarchical regression. Analysis results show that the maintenance cost of RC beam bridge assets is only affected by bridge age.
The influencing factors analysis results were used to calculate the routine maintenance expenditures of a RC beam bridge. A forecasting model of a combination of quantity and price was put forward, based on the principle of nonlinear regression and a forecasting model of routine maintenance amount based on the principle of a time series and the idea of the separation of quantity and price. The engineering examples were verified. The results showed that the forecasted values of the two forecasting methods are highly correlated with the actual values, and the differences between the forecasted values and the actual ones were of no statistical significance.
The research conclusion has important practical application value for improving the accuracy of the cost calculation of routine maintenance of RC beam bridges in Shaanxi Province and also for improving the cost calculation module and data module processing method of the bridge maintenance management system.
As the research conclusions in this paper were obtained based on the maintenance data of RC beam bridges of expressways in Shanxi Province, whether the results can be applied directly to other provinces should be verified. Other provinces can adopt the method in this paper to determine the appropriate influencing factors to calculate the expenditures.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work is part of a research project on “Research on the Annual Consumption Index of Routine Maintenance of Highway,” funded by the Shaanxi Provincial Department of Transportation under Grant no. 1634k and the Central University Basic Research Business Funding Project under Grant no. 300102218112. Data support was provided by the Shaanxi Provincial Communications Group Corporation. Their support is gratefully acknowledged.
References
 T. F. Fwa, K. C. Sinha, and J. D. N. Riverson, “Highway routine maintenance programming at network level,” Journal of Transportation Engineering, vol. 114, no. 5, pp. 539–554, 1988. View at: Publisher Site  Google Scholar
 T. F. Fwa, C. Y. Tan, and W. T. Chan, “Roadmaintenance planning using genetic algorithms. II: analysis,” Journal of transportation engineering, vol. 120, no. 5, pp. 710–722, 1994. View at: Publisher Site  Google Scholar
 Y. M. Wu and J. Chang, “Study on life cycle cost model of existing bridges based on maintenance management system,” Journal of Suzhou University of Science and Technology: Engineering and Technology, vol. 23, no. 1, pp. 44–49, 2010. View at: Google Scholar
 J. L. Smith, “Lifecycle cost analysis of reinforced concrete bridges rehabilitated with CFRP,” Theses and DissertationsCivil Engineering, 2015, https://uknowledge.uky.edu/ce_etds/33. View at: Google Scholar
 Z. J. Song, Research on Information Management of Integrated Construction and Maintenance of Highway Bridge, Southeast University, Nanjing, China, 2015.
 Y. H. Dong, The Research of Bridge Maintenance Management System, Chang’an University, Xi’an, China, 2014.
 N. M. Okasha and D. M. Frangopol, “Computational platform for the integrated lifecycle management of highway bridges,” Engineering structures, vol. 33, no. 7, pp. 2145–2153, 2011. View at: Publisher Site  Google Scholar
 G. Barone and D. M. Frangopol, “Lifecycle maintenance of deteriorating structures by multiobjective optimization involving reliability, risk, availability, hazard and cost,” Structural Safety, vol. 48, pp. 40–50, 2014. View at: Publisher Site  Google Scholar
 W. Woldemariam, J. MurilloHoyos, and A. M. Samuel Labi, “Estimating annual maintenance expenditures for infrastructure: artificial neural network approach,” Journal of Infrastructure Systems, vol. 22, no. 2, Article ID 04015025, 2016. View at: Google Scholar
 G. S. Ali, M. Volovski, S. Alqadhi, and M. Alinizzi, “Estimation of annual repair expenditure for interstate highway bridges,” Infrastructure Asset Management, vol. 6, no. 1, pp. 40–47, 2019. View at: Publisher Site  Google Scholar
 M. Moloney, T. McKenna, K. Fitzgibbon, and E. McKeogh, “Quality data for strategic infrastructure decisions in Ireland,” Infrastructure Asset Management, vol. 4, no. 2, pp. 40–49, 2017. View at: Publisher Site  Google Scholar
 A. Adetola and J. Goulding, “Collaborative framework for road infrastructure management,” Infrastructure Asset Management, vol. 3, no. 2, pp. 71–80, 2016. View at: Publisher Site  Google Scholar
 A. Taggart, L. Tachtsi, M. Lugg, and H. Davies, “UKRLG framework for highway infrastructure asset management,” Infrastructure Asset Management, vol. 1, no. 1, pp. 10–19, 2014. View at: Publisher Site  Google Scholar
 M. Irfan, M. B. Khurshid, A. Ahmed, and S. Labi, “Scale and condition economies in asset preservation cost functions: case study involving flexible pavement treatments,” Journal of Transportation Engineering, vol. 138, no. 2, pp. 218–228, 2012. View at: Publisher Site  Google Scholar
 ““Estimating cost per lane mile for routine highway operations, and maintenance,” Midwest Regional University Transportation Center, University of Wisconsin, Madison, WI, USA, 2011, Project 0712. View at: Google Scholar
 B. T. Adey, B. García de Soto, and C. Senn, “Evaluating the operation and routine maintenance of municipal roads in Switzerland,” Infrastructure Asset Management, vol. 4, no. 3, pp. 96–108, 2017. View at: Publisher Site  Google Scholar
 A. AlMansour and K. C. Sinha, “Economic analysis of effectiveness of pavement preventive maintenance,” Transportation Research Record, vol. 1442, 1994. View at: Google Scholar
 P. Lu and D. Tolliver, “Multiobjective pavementpreservation decision making with simulated constraint boundary programming,” Journal of Transportation Engineering, vol. 139, no. 9, pp. 880–888, 2013. View at: Publisher Site  Google Scholar
 M. Hagood, Highway Routine Maintenance Cost Estimation for Nevada, University of Nevada, Las Vegas, NV, USA, 2014.
 R. Gibby, R. Kitamura, and H. Zhao, “Evaluation of truck impacts on pavement maintenance cost,” Transportation Research Record, vol. 1262, pp. 48–56, 1990. View at: Google Scholar
 M. Volovski, J. MurilloHoyos, T. U. Saeed, and S. Labi, “Estimation of routine maintenance expenditures for highway pavement segments: accounting for heterogeneity using randomeffects models,” Journal of Transportation Engineering, Part A: Systems, vol. 143, no. 5, p. 04017006, 2017. View at: Publisher Site  Google Scholar
 Z. W. Qiu, M. J. Feng, and L. J. Li, “Research on calculation methods of highway maintenance funds,” Road Construction Machinery and Construction Mechanization, vol. 27, no. 4, pp. 5556, 2010. View at: Google Scholar
 C. L. Wu, Research on Estimation Model of Highway Asphalt Pavement Maintenance Cost Based on Target Maintenance Cycle, Chongqing Jiao tong University, Chongqing, China, 2016.
 S. A. Ola and I. A. Ugiomoh, “Cost of road maintenance in developing countries,” Journal of Transportation Engineering, vol. 112, no. 4, pp. 440–445, 1986. View at: Publisher Site  Google Scholar
 R. Yue, Research of Maintenance Amount Prediction on the Highway with Intelligent Algorithm, Department of Computer Science and Technology, Wuhan Univ. of Technology, Wuhan, China, 2015.
 C. X. Zhang and W. Wang, “Study on the method of determining preventive maintenance time for bridges based on optimal cost effectiveness,” Highway, vol. 7, pp. 159–162, 2013. View at: Google Scholar
 J. G. Qiao and H. X. Hou, “Study on optimization of bridge transformation plan based on life cycle cost,” Journal of China and Foreign Highway, vol. 34, no. 5, pp. 93–95, 2014. View at: Google Scholar
 S. Kubota and I. Mikami, “Development of product data model for maintenance in concrete highway bridges,” Applied Computational Intelligence and Soft Computing, vol. 2013, Article ID 148785, 12 pages, 2013. View at: Publisher Site  Google Scholar
 I.C. Wu and S.H. Hsieh, “A framework for facilitating multidimensional information integration, management and visualization in engineering projects,” Automation in Construction, vol. 23, pp. 71–86, 2012. View at: Publisher Site  Google Scholar
 D. M. Frangopol, Y. Dong, and S. Sabatino, “Bridge lifecycle performance and cost: analysis, prediction, optimisation and decisionmaking,” Structure and Infrastructure Engineering, vol. 13, no. 10, pp. 1239–1257, 2017. View at: Publisher Site  Google Scholar
 Y. Zhang, D. A. Novick, A. Hadavi, and R. J. Krizek, “Lifecycle cost analysis of bridges and tunnels,” Construction Research Congress, vol. 183, pp. 1–9, 2005. View at: Publisher Site  Google Scholar
 B. S. Yanev, “Bridge maintenance life cycle cost assessment,” in Proceedings of the First USJapan Workshop on LifeCycle Cost Analysis and Design of Civil Infrastructure Systems, Honolulu, HI, USA, August 2000. View at: Google Scholar
 K. C. Sinha, S. A. Labi, B. G. McCullouch, A. Bhargava, and Q. Bai, “Updating and enhancing the Indiana bridge management system (IBMS),” Indiana Department of Transportation and Purdue University, West Lafayette, IN, USA, 2009, Publication FHWA/IN/JTRP2008/30, Joint Transportation Research Program. View at: Google Scholar
 T. Hegazy, E. Elbeltagi, and H. ElBehairy, “Bridge deck management system with integrated lifecycle cost optimization,” Transportation Research Record: Journal of the Transportation Research Board, vol. 1866, no. 1, pp. 44–50, 2004. View at: Publisher Site  Google Scholar
 S. Labi, M. M. Rodriguez, and K. C. Sinha, “Assessing preservation needs for a bridge network: a comparison of alternative approaches,” Structure and Infrastructure Engineering, vol. 4, no. 3, pp. 221–235, 2008. View at: Publisher Site  Google Scholar
 K. C. Sinha, S. Labi, and M. Rodriguez, “Procedures for the estimation of pavement and bridge preservation costs for fiscal planning and programming,” Indiana Department of Transportation and Purdue University, West Lafayette, IN, USA, 2005, Publication FHWA/IN/JTRP2005/17, Joint Transportation Research Program. View at: Publisher Site  Google Scholar
 J. Jebreen and D. W. Johnston, “Bridge maintenance level of service optimization based on an economic analysis approach,” 1995. View at: Google Scholar
 D. Lee, “Comparative maintenance costs of different bridge types,” in Bridge Management, pp. 145–153, Springer, Boston, MA, USA, 1990. View at: Google Scholar
 L. C. Neves and D. M. Frangopol, “Condition, safety and cost profiles for deteriorating structures with emphasis on bridges,” Reliability Engineering and System Safety, vol. 89, no. 2, pp. 185–198, 2005. View at: Publisher Site  Google Scholar
 C. Zhou, Research on Condition Assessment of Existing Bridges and the Maintenance Managing Method, Zhejiang University, Hangzhou, China, 2008.
 S.Y. Ok, S.Y. Lee, and W. Park, “Robust multiobjective maintenance planning of deteriorating bridges against uncertainty in performance model,” Advances in Engineering Software, vol. 65, pp. 32–42, 2013. View at: Publisher Site  Google Scholar
 F. Biondini and D. M. Frangopol, “Lifecycle performance of deteriorating structural systems under uncertainty: review,” Journal of Structural Engineering, vol. 142, no. 9, Article ID F4016001, 2016. View at: Publisher Site  Google Scholar
 D. Kifokeris, J. A. Campos e Matos, Y. Xenidis, and L. Bragança, “Bridge quality appraisal methodology: application in a reinforced concrete overpass roadway bridge,” Journal of Infrastructure Systems, vol. 24, no. 4, Article ID 04018034, 2018. View at: Publisher Site  Google Scholar
 J. L. Qin, Research on Key Technology to Bridge Maintenance Planning, Beijing University of Civil Engineering and Architecture, Beijing, China, 2014.
 Z. J. Hu, Study of RC Beam Bridge Condition State Assessment System Based on Al, Tongji University, Shanghai, China, 2016.
 G. Du, M. Safi, L. Pettersson, and R. Karoumi, “Life cycle assessment as a decision support tool for bridge procurement: environmental impact comparison among five bridge designs,” The International Journal of Life Cycle Assessment, vol. 19, no. 12, pp. 1948–1964, 2014. View at: Publisher Site  Google Scholar
 Z. H. Zong, Z. G. Yang, Y. F. Xia, and W. W. Yuan, “Vehicle load model of Xinyihe bridge under congestion operation,” China Journal of Highway and Transport, vol. 29, no. 2, pp. 44–51, 2016. View at: Google Scholar
 W. S. Han, Y. Zhao, and H. J. Liu, “Research status and development trend of windvehiclebridge coupled vibration,” China Journal of Highway and Transport, vol. 31, no. 7, pp. 1–23, 2018. View at: Google Scholar
 J. M. Wooldridge, Introductory Econometrics, China Renmin University Press, Beijing, China, 2014.
 H. Z. Qiu, Quantitative Research and Statistical Analysis, Chongqing University Press, Chongqing, China, 2015.
Copyright
Copyright © 2019 Xiaoli Shi 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.