About this Journal Submit a Manuscript Table of Contents
Journal of Applied Mathematics
Volume 2013 (2013), Article ID 252808, 5 pages
Research Article

Determination of Pavement Rehabilitation Activities through a Permutation Algorithm

1Road Management Division, Seoul Metropolitan Government, Republic of Korea
2Department of Civil Engineering, Seoul National University of Science & Technology, Republic of Korea
3ROADTECH Corporation, Ansan-si, Gyunggi-do, Republic of Korea

Received 15 April 2013; Accepted 29 May 2013

Academic Editor: Zong Woo Geem

Copyright © 2013 Sangyum Lee 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.


This paper presents a mathematical programming model for optimal pavement rehabilitation planning. The model maximized the rehabilitation area through a newly developed permutation algorithm, based on the procedures outlined in the harmony search (HS) algorithm. Additionally, the proposed algorithm was based on an optimal solution method for the problem of multilocation rehabilitation activities on pavement structure, using empirical deterioration and rehabilitation effectiveness models, according to a limited maintenance budget. Thus, nonlinear pavement performance and rehabilitation activity decision models were used to maximize the objective functions of the rehabilitation area within a limited budget, through the permutation algorithm. Our results showed that the heuristic permutation algorithm provided a good optimum in terms of maximizing the rehabilitation area, compared with a method of the worst-first maintenance currently used in Seoul.

1. Introduction

This study involves the objective function of the rehabilitation area to be maximized, using the evaluation model of pavement performance, which is the structural and functional behavior of a road section or length of pavement; furthermore, using budget limitations for pavement rehabilitation, the more deteriorated areas are fixed, based on solving a maximizing problem in this paper. Regression analysis was used to characterize the structural performance under rutting and fatigue deterioration conditions (Figure 1).

Figure 1: Pavement distress: (a) rutting and (b) fatigue cracking.

The functional behavior of the pavement is based on the passenger’s perception of the level of service of the pavement and is related to the comfort quality (i.e., smoothness and safety) of the ride. Pavement riding quality was determined in terms of the international roughness index (IRI), as shown in Figure 2. The measurement of riding serviceability of the pavement used the quarter-car simulation model, which accounts for the sprung mass of the vehicle body, the suspension spring and damper constants, the unsprung mass of the suspension, tire, and wheel, and the spring constant of the tire [1]. For the calculation of IRI, the road profile was filtered using the quarter-car simulation, with specific parameter values that corresponded to a simulation speed of 80  (Figure 2) [2].

Figure 2: International roughness index (IRI) determination through measurement of pavement roughness and quarter-car simulation.

In this study, we created an efficient optimization method, based on the heuristic permutation algorithm, to determine the optimum pavement rehabilitation area for a specified budget. The proposed algorithm, based on the harmony search algorithm [36], was developed from field measurement data obtained from a road survey vehicle and quarter-car simulations. Section 2 presents a review of the relevant literature of structural and functional deterioration models. Section 3 describes the determination of the optimum rehabilitation area, using the permutation algorithm. Finally, Section 4 presents the application of a case study.

2. Structural and Functional Deterioration Models

Much attention has been paid to the deterioration model of pavement performance in terms of fatigue cracking and rut depth as permanent deformation [5]. In this approach, the initial value of fatigue cracking is assigned a value of 10 (i.e., undamaged status), which decreases as described in (1): where is the Seoul pavement index (SPI) of the crack deterioration model for Seoul city, and is the percentage of crack area of the pavement surface. The decreases as the crack area increases.

The deterioration model related to rut depth is based on the averaged measurement of permanent deformation in specific pavement sections and is given as where is the SPI of the rutting deterioration model in Seoul city, and is the rut depth (mm). The intact status of 10 (the initial value) decreases as the rut depth increases.

The functional pavement condition of the IRI value has a similar form as the structural deterioration models specified in (1) and (2): where is the functional pavement deterioration index of the IRI () value, related to ride quality, and is the international roughness index. The unit for the IRI of meters per kilometer () is based on the slope measurement of the pavement roughness.

Finally, the combination model for determining pavement integrity considers both structural and functional deterioration: where is the integral SPI consisting of a combination of the structural and functional SPIs, including fatigue cracking, rutting, and IRI conditions, for surface deterioration of the pavement in Seoul. The threshold SPI value for pavement rehabilitation is 5.

Figure 3 shows three cases; each case is dependent on one of the indices, , , or , while the other two indices remain constant. For example, the distress curve of the integral SPI can be calculated for a fixed value of 3 mm for the rut depth and 2  for the IRI value, as shown in Figure 3(a). The distress curves of the integral SPIs shown in Figures 3(b) and 3(c) are obtained from fixed values of 1% fatigue cracking and 2  IRI and 1% fatigue cracking and 3 mm rutting, respectively.

Figure 3: Distress models: (a) fatigue cracking, (b) rutting, and (c) IRI.

If the SPI integral is <5, then the rut depth, cracking area, and IRI are checked to determine the appropriate rehabilitation technique. The procedure is given as follows: check whether the rut depth is ≥18 mm; if so, a modified asphalt overlay of 15 cm thickness is applied; check whether the crack area on the pavement surface is ≥30%; if so, a modified asphalt overlay of 10 cm thickness is applied; if 17% ≤ crack area <30%, then a modified asphalt overlay of 5-cm thickness is applied; check whether the IRI is ≥7.3 ; if so, a modified asphalt overlay of 5 cm thickness is applied; for all other cases, a traditional asphalt overlay of 5 cm thickness is applied. The costs of the above rehabilitation overlays are shown in Table 1.

Table 1: Costs of the rehabilitation overlays.

3. Permutation Algorithm

The permutation algorithm conceptualizes a behavioral phenomenon of randomly generated permutations, in which a randomly selected permutation set continues to improve the optimization solution, similar to the search for a better state of harmony [6, 7]. This section describes the proposed permutation algorithm, based on the heuristic algorithm, which searches for an optimized solution. First, we present an application overview of the permutation algorithm, used to formulate the solution vectors for optimization. The object function is then evaluated. Finally, the application setup of the permutation algorithm is explained in terms of pavement rehabilitation.

3.1. Algorithm Procedure

Budget limitations for pavement rehabilitation in Seoul city require that the worst deteriorated areas are fixed first (worst-first maintenance) with an asphalt overlay, followed by subsequent areas with slightly less wear. The worst-first maintenance is currently used in Seoul city. However, this rehabilitation procedure has not been evaluated to determine whether or not it provides an increase in the overall averaged value, an indicator of the overall pavement surface condition for Seoul city. The more deteriorated pavement surface areas should be rehabilitated through the permutation algorithm to optimize resources within a limited budget, compared with the current rehabilitation method, such as worst-first maintenance.

The procedure for the permutation algorithm consists of Steps 1 through 4. Algorithm preparation is set up in Step 1, as follows: the sectional values of (4) are calculated, regarding the number of deteriorated sections; these values are then checked to determine whether or not they are below the threshold value of 5, to decide the appropriate rehabilitation method; the deteriorated sections below an SPI value of 5 are then grouped; the deteriorated sections are sorted in ascending order, in terms of their values; and finally, the worst-first maintenance sequence is set up.

In Step 2, the optimization problem is specified as follows: where is the th deteriorated area, which is computed by multiplying the length of the area by its width. For example, the worst deteriorated area, which shows the lowest value of , can be found in DeArea1 (e.g., ). is the specified budget limitation.

In Step 3, a sorted column vector of the worst-first sequence in Step 1 is filled in the first row of a permutation memory (PM) matrix. The second row of the PM matrix is filled with a randomly generated permutation vector, consisting of the columns of the rehabilitation sections covered by the limited budget, as well as the corresponding total area, , as shown in where the first row represents the sorted column vector of the worst-first maintenance (e.g., is the first row (superscript) and the second worst deterioration area (subscript)); the second row represents the randomly generated permutation vector; is the total number of deteriorated sections which show an value <5, based on satisfying the limited city budget; and is the total deteriorated area. Some of the deteriorated sections will not be rehabilitated due to budget constraints.

In Step 4, a comparison process is carried out between the first and second rows in (6). Firstly, the second row in the PM matrix is randomly permutated until is larger than in (6). Secondly, a new permutation vector is randomly generated until the new permutation set provides a larger rehabilitation area in comparison with the second row, as given in If the new permutation set satisfies the above criterion of the larger rehabilitation area then, it can be added to the third row in the PM matrix, as shown in (7) and

Additionally, based on an evaluation of the objective function value of in (5), another new permutation set is randomly searched to update the fourth row of the PM matrix, consisting of a larger rehabilitation area, compared with the third row. Finally, if the stopping criterion (i.e., the maximum number of iterations) is satisfied, the procedure of adding better permutation sets into the PM matrix is terminated. Otherwise, Step 4 is repeated.

3.2. Application Setup for the Maximized Rehabilitation Area of Deteriorated Pavement Sections

In terms of applying the permutation algorithm to the area maximization problem of pavement rehabilitation, the costs of rehabilitation overlays shown in Table 1 are considered, as well as the specified budget limitations of this study. The optimization problem can be solved using the permutation algorithm, based on the setup shown in Table 2. For the simplicity of applying the permutation algorithm, the section numbers shown in Table 2 are sequentially assigned in ascending order, in terms of the integral SPI values.

Table 2: Setup of the measured field data used for the heuristic permutation algorithm.

4. Results

Our permutation algorithm was used to evaluate 150 pavement sections being considered for rehabilitation, as shown in Table 2. According to Step 3, the worst-first sequence was stored in the PM matrix, and the second row was obtained by randomly permutated selection, as shown in (9).

Permutation memory (PM) matrix generation for worst-first maintenance, along with the first permutation, is as follows:

There were 50 sections that had an integral SPI value of <5. The 30 sections in the first row can be rehabilitated, using a limited budget of seven billion won (Korean currency), as shown in (9). The second row was updated until the second row chosen by the permutation algorithm provided more rehabilitated areas, compared with the first row of the worst-first sequence. The rehabilitated sections were limited to 25 sections, which provided more rehabilitated area than the approach using worst-first maintenance, as shown in (10).

PM matrix consisting of a better solution in the second row, compared with the worst-first maintenance, is as follows:

The columns above 25th in the second row of the PM matrix were assigned a zero value, as shown in (10). Additionally, a third row was added to the PM matrix when the randomly generated third row resulted in better selection than the second row group, in terms of rehabilitation area covered by the limited budget, as shown in (11). Equation (12) and Figure 4 show the best selection resulting in the largest rehabilitation area, after ~1000 iterations.

Figure 4: Rehabilitated area as a function of the number of iterations.

PM matrix updating procedure

Final PM matrix obtained using the permutation algorithm

5. Conclusion

A permutation algorithm has been described to determine the optimal pavement rehabilitation planning which maximizes the rehabilitation area, using the empirical models of fatigue and rutting deterioration, as well as the functional model of IRI. To evaluate the suggested permutation algorithm, the worst-first maintenance method used in the current rehabilitation technique was compared. Our results indicated that the developed permutation algorithm, described in this study, provided a better solution in terms of covering a larger deteriorated area within the limited city budget, when compared with the worst-first maintenance, which is currently used in Seoul city.


This work was supported by the National Research Foundation of Korea (NRF) Grant funded by the Korea government (MEST) (no. 2011-0030848).


  1. R. Hass, W. R. Hudson, and J. Zaniewski, Modern Pavement Engineering, Krieger Publishing Company, Malabar, Fla, USA, 1994. View at Zentralblatt MATH
  2. M. W. Sayers and S. M. Karamihas, “Interpretation of road roughness profile data,” Final Report Contract DTFH 61-92-C00143, Federal Highway Administration, 1996.
  3. S. Mun and Z. W. Geem, “Determination of viscoelastic and damage properties of hot mix asphalt concrete using a harmony search algorithm,” Mechanics of Materials, vol. 41, no. 3, pp. 339–353, 2009. View at Publisher · View at Google Scholar · View at Scopus
  4. S. Mun and Z. W. Geem, “Determination of individual sound power levels of noise sources using a harmony search algorithm,” International Journal of Industrial Ergonomics, vol. 39, no. 2, pp. 366–370, 2009. View at Publisher · View at Google Scholar · View at Scopus
  5. Y. Suh, S. Mun, and I. Yeo, “Fatigue life prediction of asphalt concrete pavement using a harmony search algorithm,” KSCE Journal of Civil Engineering, vol. 14, no. 5, pp. 725–730, 2010. View at Publisher · View at Google Scholar · View at Scopus
  6. Z. W. Geem, “Particle-swarm harmony search for water network design,” Engineering Optimization, vol. 41, no. 4, pp. 297–311, 2009. View at Publisher · View at Google Scholar · View at Scopus
  7. J. Fourie, R. Green, and Z. W. Geem, “Generalized adaptive harmony search: a comparative analysis of modern harmony search,” Journal of Applied Mathematics, vol. 2013, Article ID 380985, 13 pages, 2013. View at Publisher · View at Google Scholar