Sensitivity Analysis of Rigid Pavement Design Based on Semi-Empirical Methods: Romanian Case Study

: Due to the intensive process of road construction or rehabilitation of pavement caused by an increase in trafﬁc volume, in the ﬁeld of rigid pavement design and research in Romania, we can say that there is a need to improve the design method. In the last decade, more and more researchers have been concerned about climate change and the increase in trafﬁc volume; hence, there is a need for a renewal of the climatological, as well as trafﬁc, databases because these are part of the input data for the design process. The design method currently used in Romania for jointed plain concrete pavement design is NP081/2002. The limitation of the data and the lack of lifetime estimation of structural and functional performance are the main aspects that need to be addressed in the new design procedure. The Mechanistic–Empirical Pavement Design (MEPDG) method offers the possibility of the design of pavement structures by estimating the structural and functional performances. This paper aims to obtain a comparative study of these two methods for the analysis of the input data collected from the ﬁeld corresponding to the three failure criteria, while the symmetry of the characteristics of the material and their asymmetrical thicknesses are compared, thus contributing to the design of viable and long-lasting pavement structures using a rigid pavement with the speciﬁc characteristics of the mountainous area in northeastern Romania on the national road DN17 Suceava—Vatra Dornei. The novelty of this study consists of the implementation of the mechanistic–empirical method MEPDG instead of the old NP081/2002 method used in Romania.


Introduction
The design of road structures involves the estimation of structural and functional performances. With greater design accuracy, all maintenance and repair costs are reduced, and the corrective activities can be planned in advance. For the design of rigid pavements, over time, several methods have been developed, including predominantly theoretical (the Westergaard, Westergaard/Ioannides-USA diagram method and the Picket and Ray USA-method based on the finite element calculation scheme) and semi-empirical methods that combine theoretical design relationships with experimental results (modified Westergaard method, Highway Agency-UK method, AASHTO-American Association of State Highway and Transportation Officials method for plain concrete and for continuously reinforced concrete pavement, the MEPDG method-Mechanistic-Empirical Pavement Design USA, the Romanian method NP 081/2002, etc.) [1,2]. The semi-empirical methods have the advantage of estimating the lifetime of the pavement design on the basis of performances obtained in practice, whereas the input data allow the personalization of each case according to the climatic conditions, soil condition, pavement structure, and the traffic loads, as discussed elsewhere [3][4][5]. Gaspar et al. [6] and Plescan [7] presented comparative studies on the application of these methods, from the geographical and topographic point of view, and others discussed optimization methods of the pavement design by recalibrating the obtained models [8][9][10]. Of all the mechanistic-empirical methods, the MEPDG method is considered one of the most realistic methods, not only for the design of pavement models but also because it offers the possibility of sizing the pavement structures by modifying the parameters of the transfer function, as discussed elsewhere [11][12][13][14][15]. It was found that both the loads imposed on the ground and the time accelerate the structural and functional deterioration of rigid pavements, leading to its premature failure. Therefore, the sensitivity analysis of the most significant predictive factor of pavement performance analyzed in a known geoclimatic context and under certain traffic conditions represents the current challenge of the researchers in the field, as discussed by [16,17]. Thus, the quality of a jointed plain concrete pavement (JPCP) design is given by the degree of sensitivity of the results obtained by modifying the scenario or the input dataset, as discussed by Zhu et al. and Mu et al. [18,19]. Amin [4] presented comparatively different deterministic and stochastic approaches for calculating the pavement performance curves, highlighting the advantages and disadvantages of each approach. Chong et al. [20] presented a method of multi-objective optimization of rigid pavement design starting from the objective function of reducing energy consumption and greenhouse gas emissions, in correlation with the thickness of the pavement and the surface roughness. Although the design of JPCP is carried out under the conditions of known parameters (traffic volume, climatic conditions, material), it has been found that traffic volume increases in logarithmic progression, which affects the structural and functional performances of the pavement in a shorter time than predicted by the design methods. Bayrak and Hinislioglu [21] presented a comprehensive study that combines different analysis methods to determine the optimum thickness of the layer and the quality of the materials at a low cost-efficiency ratio.
Due to the particularities of each area, namely, environmental conditions, type of traffic, or characteristics of the base layer, the design of the road structure differs. Starting from the variability in the methods for designing the road pavement, the objective of the paper consists of a critical analysis of the efficiency of the MEPDG method for designing a jointed plain concrete pavement (JPCP). In order to achieve the proposed objective, a detailed study on the principles of the design methods for rigid pavement structures currently used in the world was undertaken, finally selecting and studying this method in detail by applying it to a significant case study in Romania located in a northeastern Romania mountain area on the National Road DN17 (E 576). This road was selected to predict the jointed plain concrete pavement models and analyze them from the performance point of view, compared to the Romanian design method according to the NP081/2002 standard. The novelty of this study consists of predicting the behavior of 72 pavement structures with different physical-mechanical characteristics, according to their response over time from the point of view of cracking, faulting, and roughness (international roughness index, IRI).  [2], the jointed plain concrete pavement is based on the following main steps: establishing the design traffic; establishing the bearing capacity of the foundation soil; conceiving the rigid pavement structure; establishing the bearing capacity at the level of the base course; determining the thickness of the surface course. Preliminary studies have been undertaken to obtain relevant data concerning the composition, intensity, and evolution of traffic, the geotechnical characteristics of the foundation soil, and the hydrologic regime on the site, provided by the Regional Roads and Bridges Directorate Iasi. The design traffic denoted Nc is expressed in millions of standard axles and is calculated with Equation (1), as a function of the average annual daily traffic (AADT) presented in Table 1, for a perspective period of 30 years, and a coefficient for the transversal distribution c rt for the lanes of 0.50.

Summary of Basic
where p p is the perspective period, c rt is the transversal distribution, MZA/AADT is the annual average daily traffic, p k is the coefficient of evolution, and f ek is the equivalence coefficient.

The Bearing Capacity of the Foundation Soil
The coefficient of subgrade reaction, K 0 , is determined in correlation with the climate and hydrologic regime and soil type, as given in Table 2. According to Romanian standard STAS 1709/2 [22], the hydrological regime is distributed as follows: hydrological regime 1, corresponding to the favorable conditions; hydrological regime 2, corresponding to the medium and unfavorable hydrological conditions, where 2a represents embankment road sectors, with a minimum height of 1.00 m, 2b is considered for sectors of road located in the mound with a height beneath 1.00 m, at ground height, at a mixed profile, and in cutting. Table 2 presents the soil type according to Romanian standard STAS 1243-1988. Because the road segment to which we applied the Romanian dimensioning method is characterized by soil type P3, climate type III, and hydrologic regime 2b, following Table 2, the coefficient of subgrade reaction is K 0 = 42 MN/m 3 . The bearing capacity at the level of the base course is expressed by the coefficient of reaction at the surface of the base course, K, as given in Figure 1a, which is determined as a function of the coefficient of subgrade reaction K 0 and the equivalent thickness of the base/subbase courses in cm, H ech . The equivalent thickness of the base/subbase courses is calculated as the sum of the equivalent thicknesses of these layers, as given by Equation (2).
where n is the number of layers, h i is the actual thickness of the layer "i", in cm, and a i is the equivalence coefficient for the layer "i". Knowing that the equivalence coefficient for ballast subbase is 0.75 (NP081/2002), it turns out that H ech = 25 × 0.75 = 18.75 cm.

The Mechanistic-Empirical Method MEPDG
The MEPDG method for JPCP design is illustrated by the flowchart presented in Figure 2 [24]. According to this flowchart, in the first stage, a JPCP structure is proposed that takes into consideration the specific site conditions of traffic, climatic conditions, and soil type. Then, the failure criteria are set on the basis of the desired performance that will be achieved at the end of the design life (for example, an acceptable level of faulting, cracking, and IRI for JPCP). Reliability levels are also selected for each of these performance indicators. Then, these input data are processed by the program using the finite element method for each axle type and the load over the design life. The key distresses, faulting, cracking, and roughness, are predicted month by month during the design period, using mechanistic-empirical performance models existing in the MEPDG methodology [25][26][27]. The roughness is then estimated on the basis of the value of the initial IRI, the distresses that occur over time, and the site conditions. The estimated performance of this structure is evaluated at the established reliability level. If the designed JPCP structure does not meet the set criteria, then it must be modified accordingly, and the design procedure applied to this modified structure is repeated. For the type of soil P3, climate type III, and hydrologic regime 2b, the dynamic elastic modulus is K = 58 MN/m 3 .

Determining the Thickness of the Concrete Slab
The design criterion is expressed as follows (Equation (3)): where σ is the tensile stress from bending in the concrete slab, determined following the design hypothesis, and σ a,t is the allowable tensile stress from bending. The allowable tensile stress from bending (σ a,t ) is determined using Equation (4).
where R k inc is the specific bending strength of the concrete at 28 days, α is the coefficient of concrete strength increase from 28-90 days, equal to 1.1, N c is the traffic for the design period, and γ is a coefficient equal to 0.05.
The design hypotheses are as follows: a for technical class I and II roads; σ = σ t + 0.8 × 0.65 × σ t∆t ≤ σ σ t,a for technical class III and IV roads; σ = σ t ≤ σ t,a , for technical class V roads. Because the analyzed road sector has a technical class III, we chose the second hypothesis. Depending on the modulus of subgrade reaction K = 58 MN/m 3 and on the allowable flexural strength σ t,a = 3.48 MPa, the slab thickness results as H slab = 23 cm, following the design diagram presented in Figure 1b.

The Mechanistic-Empirical Method MEPDG
The MEPDG method for JPCP design is illustrated by the flowchart presented in Figure 2 [24]. According to this flowchart, in the first stage, a JPCP structure is proposed that takes into consideration the specific site conditions of traffic, climatic conditions, and soil type. Then, the failure criteria are set on the basis of the desired performance that will be achieved at the end of the design life (for example, an acceptable level of faulting, cracking, and IRI for JPCP). Reliability levels are also selected for each of these performance indicators. Then, these input data are processed by the program using the finite element method for each axle type and the load over the design life. The key distresses, faulting, cracking, and roughness, are predicted month by month during the design period, using mechanistic-empirical performance models existing in the MEPDG methodology [25][26][27]. The roughness is then estimated on the basis of the value of the initial IRI, the distresses that occur over time, and the site conditions. The estimated performance of this structure is evaluated at the established reliability level. If the designed JPCP structure does not meet the set criteria, then it must be modified accordingly, and the design procedure applied to this modified structure is repeated.

Climate Data for JPCP Model Development
The Romanian National Road DN17 (E 576) placed between Suceava town and Vatra Dornei city, i.e., the road sector from km 217 to 218 (Figure 3), was selected to pre-

Climate Data for JPCP Model Development
The Romanian National Road DN17 (E 576) placed between Suceava town and Vatra Dornei city, i.e., the road sector from km 217 to 218 (Figure 3), was selected to predict the jointed plain concrete pavement models and analyzing them from the performance point of view.

Climate Data for JPCP Model Development
The Romanian National Road DN17 (E 576) placed between Suceava town and Vatra Dornei city, i.e., the road sector from km 217 to 218 (Figure 3), was selected to predict the jointed plain concrete pavement models and analyzing them from the performance point of view.  From the geological point of view, in this area, cohesive soils (such as silty and sandy clays) and the local slopes consist of the base layer, as considered from other sedimentary and metamorphic rocks. The climatic conditions in the area include severe winters (−250 • C) and hot summers (+30 • C) and a high level of seasonal precipitation between 1200 mm/m 2 and 1400 mm/m 2 , with depth of frost between 100 and 110 cm. The data were collected from the Meteorology National Agency of Romania (ANM), Călimănes , ti climate station. Knowing that the demonstrative road sector is located at 45.38 • latitude, −89.28 • longitude, 485.5 m (1593 ft) elevation, and 3 m (10 ft) depth of water table, the climate data type NARR-North American Regional Reanalysis from MEPDG software was correlated with climate data from the field.

Physical and Mechanical Characteristics of JPCP Models
According to the Romanian standard recommendation, the minimum thickness for the JPCP is 18 cm and the maximum value is 28 cm, whereas 23 cm is the thickness resulting for the case we studied with the Romanian standard NP 081/2002. For MEPDG analysis, three types of jointed plain concrete pavement models were analyzed from a thickness point of view with three asymmetric structures: 18 cm (denoted A); 23 cm (denoted B); 28 cm (denoted C) (Figure 3). In the other layers (layers 2 and 3), the geometric and structural characteristics were symmetric. Thus, the second layer was considered with the same thickness (30 mm) for all three cases and contained river-run gravel or cement-stabilized well-graded gravel. All studied variants of plain concrete pavements were analyzed assuming they have no dowels. The third layer was represented by foundation ground with infinite thickness. The length variations of rigid pavement slabs were presented in three cases: 4 m, 6 m, and 8 m. Figure 4 presents the basic matrix of JPCP.    Applying the MEPDG method, different variables such as the thickness of the concrete slab, type of the base layer, coefficient of thermal expansion, compressive strength, and the distance between the transverse joints were applied to the simulation of JPCP lifetime behavior. Thus, the AADTT type traffic value was 1776, and the simulation was conducted by varying the mechanical characteristics of jointed plain concrete pavement: two different values of the coefficient of thermal expansion (CTE), 10 × 10 −6 / • C and 13 × 10 −6 / • C; two different values of compressive strength Rc (28 days), 20 MPa and 55 MPa; unit weight, 2403 kg/m 3 ; 0.15 Poisson's ratio; cement type I; water/cement ratio of 0.45. Regarding the base foundation, the thickness was established at 30 cm, 0.27 Poisson's ratio, and two types of material were taken into account: river-run gravel and stabilized cement.
Combining these parameters into the MEPDG program, a total of 72 results were obtained, which allowed analysis of the variation of JPCP performance for 5, 10, 15, 20, 25, 30, and 35 years in correlation with the influence of each of the selected input parameters.
The models obtained using the MEPDG method were evaluated in terms of their behavior in cases of cracking, faulting, and roughness (IRI). The codification of each type of jointed plain concrete pavement resulting from the combination of parameters is presented in Table 3.

Results and Discussion
The analysis performed using the MEPDG method highlighted the response of all 72 structural variants established for the design life for the traffic loads and the predetermined climatic conditions, being useful in establishing the decisions regarding the rehabilitation of any road structure similar to that in the presented study. According to [1], the most important performance criteria for rigid pavements are cracking, faulting, and roughness of concrete slab surfaces. Thus, the obtained results are presented concerning the three performance criteria.

Estimation of Cracking
The most important degradation recorded by plain concrete pavement structures is their cracking. Thus, transverse cracks of the concrete slabs can occur at the surface of the concrete slab and propagate down (top-down cracking) or at the base of the concrete slab (bottom-up cracking) depending on the load and environmental conditions, properties of the materials, design features, and the conditions during construction. Bottom-up transverse cracking is caused by tensile bending stress at the bottom of the PCC -Pavement Concrete Compressive Strength slab, midway between two transverse joints [28,29]. Repeated heavy axle loading and a high positive temperature gradient (the top of the slab is warmer than the bottom of the slab) cause fatigue damage to occur along the bottom of the slab. This fatigue damage results in a transverse crack of the pavement that can propagate to the surface. Bottom-up transverse cracking is combined with top-down transverse cracking and calculated as a percentage of slabs cracked [30,31]. Transverse cracks of JPCP may have a small impact on overall structure performance as long as they have not evolved too much. Structural performance is influenced by transverse cracks when dynamic loads from the rigid pavements with a high roughness reduce the pavement life [18]. Transverse cracks can be caused by the uneven distribution of loads to the joints of the irregularity of the surface of the subbase, the insufficient thickness of the slab, and the shortcomings of the quality of the materials used. The influence of each parameter on the cracking of the concrete slab is presented below. Figure 5 presents the variation of cracking during the design life. With the increasing distance of the joints (length of the slab), an increase in cracking percentage of the slab of about 53.5% can be observed ( Figure 5). The increase in compressive strength of concrete from 20 MPa to 55 MPa led to a decrease in the cracking percentage of the concrete slab to 1.1% (length of 4 m). Increasing the distance between the joints (6 m, 8 m) led to a decrease in cracking resistance.
Thus, it can be considered that the structure AG1b1, with a granular base layer, showed the best behavior between the concrete strength and the distance between the joints. From the point of view of the JPCP with a stabilized base layer, a similar behavior to that of a granular base layer can be noted, with the rate of occurrence of cracks reduced. Thus, in the case of the granular base layer with a concrete slab with a length of 8 m and a compressive strength of 20 MPa (structure AG1a3), it predicts the occurrence of 40% of cracks in 2.5 years. In the case of a stabilized ballast structure, the cracking percentage of 40% occurs after approximately 15 years for the AS1a3 structure ( Figure 5). From the point of view of structural resistance to cracks, two optimal pavement types can be considered: AS1b1 and AS1b2. Figure 6 presents the sensitivity analysis of cracked concrete slabs with the smallest variation of percentage cracking in time. Thus, it can be considered that the structure AG1b1, with a granular base showed the best behavior between the concrete strength and the distance betwe joints. From the point of view of the JPCP with a stabilized base layer, a similar be to that of a granular base layer can be noted, with the rate of occurrence of cra duced. Thus, in the case of the granular base layer with a concrete slab with a leng m and a compressive strength of 20 MPa (structure AG1a3), it predicts the occurr 40% of cracks in 2.5 years. In the case of a stabilized ballast structure, the cracki centage of 40% occurs after approximately 15 years for the AS1a3 structure (Fig From the point of view of structural resistance to cracks, two optimal pavement ty be considered: AS1b1 and AS1b2. Figure 6 presents the sensitivity analysis of c concrete slabs with the smallest variation of percentage cracking in time.  4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 Cracking ( %)

AG2b1
AS1b2 AS2b1 AG1b1 AS1b1  Thus, it can be considered that the structure AG1b1, with a granular base lay showed the best behavior between the concrete strength and the distance between t joints. From the point of view of the JPCP with a stabilized base layer, a similar behav to that of a granular base layer can be noted, with the rate of occurrence of cracks duced. Thus, in the case of the granular base layer with a concrete slab with a length o m and a compressive strength of 20 MPa (structure AG1a3), it predicts the occurrence 40% of cracks in 2.5 years. In the case of a stabilized ballast structure, the cracking p centage of 40% occurs after approximately 15 years for the AS1a3 structure ( Figure  From the point of view of structural resistance to cracks, two optimal pavement types c be considered: AS1b1 and AS1b2. Figure 6 presents the sensitivity analysis of crack concrete slabs with the smallest variation of percentage cracking in time. It can be seen that the AS1b1 structure showed the lowest values of the cracki percentage, which is characterized by a stabilized base layer, a concrete slab with a leng of 4 m, CTE = 10, and a compressive strength of 55 MPa concrete. From the efficien point of view, all the variants determined using the MEPDG method in Figure 6 can  4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 Cracking ( %)

AG2b1
AS1b2 AS2b1 AG1b1 AS1b1 It can be seen that the AS1b1 structure showed the lowest values of the cracking percentage, which is characterized by a stabilized base layer, a concrete slab with a length of 4 m, CTE = 10, and a compressive strength of 55 MPa concrete. From the efficiency point of view, all the variants determined using the MEPDG method in Figure 6 can be taken into consideration in the design calculation, as the values of the cracking percentage were recorded below 5% during the 35 years. Cracks occurred due to repeated loading, leading to compressive stress. Once the cracks appeared, they propagated from the bottom up. An important factor that led to the increase in cracking rate on the Romanian traffic section was the variation in extreme temperatures during the year, which led to the expansion phenomenon and contraction of the water accumulated in the concrete pores or the interfaces between the JPCP layers. The MEPDG crack analysis resulted in a matrix of the rigid pavement with the lowest percentage of cracking, up to a maximum percentage of 50% over 35 years of the pavement life, as can be seen in Table 4. Table 4. Sensitivity matrix of pavement relative to the percentage of cracking.

Estimation of Faulting
The level difference between the edges of two adjacent slabs, at a transverse or longitudinal joint, represents another criterion for evaluating the performance of a JPCP, known as faulting. This distress can be caused by repeated heavy axle loads, poor joint load transfer efficiency, free moisture below the concrete slab, erosion of the base, subbase, subgrade, or shoulder material, and upward curling of the slab [32]. Transverse joint faulting is a measurement of the differential deflection across a joint, and the admissible value is 0.19685 inches (5 mm). The modeling results revealed that the lowest average faulting values were obtained for the 28 cm thick slab, reaching the admissible value in about 30 years. In the case of the 18 cm thick slab, faulting occurred after approximately 15 years, and, in the case of 23 cm thick slab, faulting occurred after 25 years (see Figure 7).
The stabilized base is a cemented, rigid material that distributes the load over a large area, and stresses in the subgrade are reduced, which provides excellent support for the concrete slab of a rigid pavement. The stabilized base material is stronger, uniform, and more water-resistant than a granular base. The use of cement-stabilized bases reduces the occurrence of faulting, as can be seen in Figure 8. Thus, the different behavior of the base layer is manifested by a decrease of approximately 85% in the case of the stabilized ballast layer compared to the granular layer after 5 years. Increasing the strength of concrete from 20 MPa to 50 MPa contributes to increasing the rigidity of the road structure to about 50% after 5-10 years.
x FOR PEER REVIEW 13 of 18 The stabilized base is a cemented, rigid material that distributes the load over a large area, and stresses in the subgrade are reduced, which provides excellent support for the concrete slab of a rigid pavement. The stabilized base material is stronger, uniform, and more water-resistant than a granular base. The use of cement-stabilized bases reduces the occurrence of faulting, as can be seen in Figure 8. Thus, the different behavior of the base layer is manifested by a decrease of approximately 85% in the case of the stabilized ballast layer compared to the granular layer after 5 years. Increasing the strength of concrete from 20 MPa to 50 MPa contributes to increasing the rigidity of the road structure to about 50% after 5-10 years.  The stabilized base is a cemented, rigid material that distributes the load over a large area, and stresses in the subgrade are reduced, which provides excellent support for the concrete slab of a rigid pavement. The stabilized base material is stronger, uniform, and more water-resistant than a granular base. The use of cement-stabilized bases reduces the occurrence of faulting, as can be seen in Figure 8. Thus, the different behavior of the base layer is manifested by a decrease of approximately 85% in the case of the stabilized ballast layer compared to the granular layer after 5 years. Increasing the strength of concrete from 20 MPa to 50 MPa contributes to increasing the rigidity of the road structure to about 50% after 5-10 years.  Table 5 presents the influences of each parameter taken into account related to the admissible value of faulting (0.196585 in, 0.5 cm). Thus, it can be noticed that the values lower than 100% fall within the admissible values, while the highest are variants of road structures that no longer meet this criterion at any given time. Therefore, in terms of the thickness of the concrete slabs (coded A, B, C), the tiles with the largest thickness (28 cm) recorded the lowest values of the faulting. The stabilized ballast provided a better pavement structure behavior, and the coefficient of expansion CTE = 13 × 10 −6 / • C favored the occurrence of faulting with values about 66% higher than CTE = 10 × 10 −6 / • C. Doubling the distance between the joints and the length of the concrete slab increased the compaction by approximately 50% in the first 5 years. Starting with the 15th year of pavement life, the structures with a small thickness of the slab and highest length with riverrun gravel subbase and a high coefficient expansion would require maintenance work. The results of MEPDG analysis highlighted that from 72 variants of JPCPs, 64 hypothetical cases met the condition whereby faulting was less than the nominal value (0.19685 in, 0.5 cm).

Evaluation of International Roughness Index (IRI)
Pavement roughness is used as a composite index of pavement quality and is directly related to overall ride quality, the factor of most importance to highway users. Because of these reasons, empirical roughness prediction models have been incorporated in the MEPDG, being characterized in terms of the IRI, measured in inches of roughness per miles in American standard (measured in meters/kilometers in Romanian Standard NP 081/2002). IRI depends on the initial input parameters used for pavement design and the climatic conditions that may affect roughness of the road through mechanisms such as shrinkage or swelling of subgrade soils and frost. IRI was established by  as an internationally accepted parameter, and it is determined with specialized vehicles that measured the displacement of the vehicle chassis placed on the rear axle [32,33]. The IRI threshold recommended by the FHWA (Federal Highway Administration of USA) in 1998 is 95-170 in/mi for acceptable ride quality and less than 95 in/mi for good road quality, being considered the most international index for pavement roughness measurement [3]. In 2008, DDOT (District Department of Transportation, USA) changed the IRI threshold in PCI (pavement condition index) for bituminous and concrete pavement in cases of urban road, as summarized by Arhin et al. [34] and Perez et al. [35]. Table 6 presents the IRI values relative to lifetime prediction. Analyzing the influence of the various parameters on the cumulative IRI values for 35 years, it is noted that, for an 18 cm thick concrete slab, a deviation from smoothness is recorded after 20 years, and a 38% increase in the thickness of the concrete slab leads to a 25% increase in the service life, which is within acceptable limits (see Table 6). The use of a stabilized cement base layer leads to a 75% increase in the operating life below 170 inches per mile compared to the duration of the granular subbase. The same aspect is also noticeable when a lower thermal expansion coefficient is used, being reduced by approximately 23% (e.g., CTE 10 (10 −6 / • C)). The use of a low-strength concrete (20 MPa) leads to the appearance of surface roughness exceeding acceptable limits after 23 years of service life, and doubling the distance between the joints would lead to the appearance of roughness over acceptable limits in only 20 years (Table 6). From the MEPDG simulation, it was observed that none of the road structures showed values below 95 in/mile, but there was a set of variants with an IRI index below 170 in/mi, as can be seen in Table 7.

Conclusions
The paper presented the lifetime prediction of jointed plain concrete pavement for a real traffic road segment analyzed using two methods: one according to Romanian standards and the other according to the more comprehensive MEPDG method. The estimation of the dimensions (thickness) of jointed plain concrete pavement using the Romanian Standard NP 081/2002 was very close to the MEPDG method (23 cm) but poorer and without time estimation of structural and functional performance. The MEPDG method analyzed the influence of each parameter, offering a realistic analysis as a function of the cumulative effect of these factors, such that the solutions obtained by the sensitivity method returned only the variants of the road structure that simultaneously fulfilled the three requirements of structural and functional performance-cracking, faulting, and IRI. Table 8 presents the seven solutions out of the 72 analyzed which can be chosen for the design of the pavement structures in traffic conditions and an environment specific to a temperate mountainous region, with temperature variations and annual humidity specific to the northern Romania area. In conclusion, the analyses performed in this significant study highlight the critical design phases within the MEPDG methodology, as well as possibilities for applying this method in other countries where national standards are based on classical methods. From a practical point of view, the presented study highlights how the MEPDG method can objectively lead to the elimination of possible variants of JPCP that do not simultaneously meet the requirements of structural and functional performance. From the results, the constructor or the beneficiary can choose the most convenient variant from the financial point of view.