Development of Deflection Prediction Model for Concrete Block Pavement Considering the Block Shapes and Construction Patterns

Concrete block pavement (CBP) is distinct from typical concrete or asphalt pavements. It is built by using individual blocks with unique construction patterns forming a discrete surface layer to bear traffic loadings.The surface structure of CBP varies depending on the block shapes and construction patterns, so it is hard to apply a general equivalent elastic modulus estimation method to define the surface structural strength. In this study, FEM analysis and dynamic loading test were carried out to develop a deflection prediction model for CBP considering the block shapes and construction patterns. Based on the analysis results, it was found that block shapes did not havemuch effect on load distribution, whereas construction patterns did. By applying the deflection prediction model to the rutting model for CBP proposed by Sun, the herringbone bond pattern showed the best performance comparing with stretcher bond or basket weave bond pattern. As the load repetition increased to 1.2 million, the rutting depth of CBP constructed by herringbone bond pattern was 2mm smaller than those constructed by the other two patterns.


Introduction
Concrete block pavement (CBP) has been widely used in sidewalks, motorways, and ports because of its aesthetic features and easy maintenance.Unlike continuous concrete pavement or asphalt pavement, CBP can bear traffic loadings with individual blocks constructed in a uniform pattern.The blocks are connected in a continuous structure with joint sand to minimize the rotation and movement of blocks for structural stability.Therefore, CBP should be viewed as a single system in which all components, the surface layer, bedding sand under the joint sand, and underlying layers are taken into account, rather than considering only the blocks as the layer carrying the applied load or providing a wearing layer or roughness [1].
The load bearing capacity of the CBP increases mainly because of shear force at the joints and compressive force caused by dilatancy of the joint sand.In general, CBP bearing capacity also changes over time compared to the initial service life; with traffic loading increase the pavement strain decreases and load bearing capacity increases due to compaction of the sublayers [2,3].According to Kuipers [4], the repeated traffic loading causes the rotation of blocks; the construction patterns would affect the performance of CBP when the rotation of blocks at narrow joint intervals reaches its limit.CBP lockup occurs only when a sufficient load is applied over a flexible base.As interlocking efficiency increases, stresses between blocks decrease and are transferred to the substructure.Kuipers [4] proposed that once interlock between the blocks is achieved, the elastic modulus of the joint sand rises around 10 times compared with the initial value.Also the interaction between paving blocks with heavy loading was determined that it could provide increased pavement stiffness and thus increased load dissipation resulting in lower transmitted stress on the subgrade [5].

Advances in Materials Science and Engineering
It is difficult to estimate the elastic modulus of the surface layer because of its discrete property.On the basis of findings from previous studies, Shackel [6] recorded the complex elastic modulus of concrete blocks and bedding sand layers measured in FWD test and laboratory test.The elastic modulus varied widely from 500 MPa to 4,000 MPa depending on conditions such as block shape and testing method.CMAA [7] suggested an initial complex elastic modulus of the surface layer (blocks and bedding sand layer) of 350 MPa (50,750 psi); the complex elastic modulus would rise to 3,100 MPa (450,000 psi) after 10,000 ESALs were applied.Although numerous researchers have presented a wide range of elastic modulus values for the layer of blocks, it seems undesirable to apply only one of them to design CBP structure [8][9][10].
One of the major structural distresses of a CBP is rutting caused by traffic loading.Rutting of CBPs is the product of horizontal and vertical movements of the blocks.Yasuhisa et al. [11] assessed the conditions of 48 CBPs laid on motorways in service and found that rutting accounted for 40% of all distresses.Panda and Ghosh [12] proved that the vertical load distribution in pavement structures was not significantly affected by construction patterns.It was found that shape, size, and thickness of the blocks have a significant influence on the behavior of CBP [13].However, Miura et al. [14] found that pavement performance was influenced more by the block shape and laying pattern than the block thickness by using rutting depth as an indicator.Based on the analysis results of finite-element model, Mampearachchi and Gunarathna [15] pointed that the performance of CBP was found to be affected more by the construction pattern.Then the viewpoint was demonstrated by the deflection basins data measured in the field test [16].Variation of deflection among different laying patterns was significantly high compared with deflection variations of block shapes.The herringbone bond had the lowest and the stack bond showed the largest deflection.Sun [17] developed a rutting prediction model based on a deflection by considering factors such as traffic load, accumulated traffic volume, resilient deflection of pavement, and block thickness.Several prediction models presented in other studies also did not consider the block shapes and construction patterns [18,19].
In this study, to develop a deflection prediction model for CBP considering the block shapes and construction patterns, a FEM analysis and a dynamic loading test were performed to analyze load-deflection behavior.By using multiregression analysis, a load distribution model based on the deflection was developed and applied to the rutting model proposed by Sun [17].Then the calculated rut depth was compared with the rutting results measured from Accelerated Pavement Testing (APT).

Rutting Prediction Model
Rutting of CBP mainly happens when using a granular aggregate base, and it is divided into permanent deformation at the base course and at the subgrade.Rutting is caused by repeated loading and is related to permanent deformation which accumulated in additional compaction of bedding sand at an early stage of service and continued traffic loads at the aggregate base, subbase, and subgrade.If a pavement's substructure (aggregate base, subbase, and subgrade) is relatively weak, shear stress at the joints from a concentrated load exceeds the shear strength at the interlock of the filling sand and blocks, which leads to rut.Rutting at the sublayers deteriorates the performance of the pavement structure and causes poor drainage and roughness.In general, rut depth is a major design standard for CBPs using an aggregate base [17,20,21].Houben et al. [18] presented the rut depth prediction model by taking traffic volume into consideration in the following equation: where RD is rut depth (mm),  is the number of load repetitions, and   ,   are model constants.
Huurman and Boomsma [19] defined rut depth occurring at the CBP through strain at the sand layers and the aggregate base.This empirical model mainly reflects the characteristics of materials in the environment condition of Netherlands, as shown in  perm at sand layer =  (  1000 ) perm at aggregate layer =  (  1000 ) where  is the number of load repetitions and , , , and  are model constants.
It is difficult to accurately estimate the stresses at sublayers of CBP due to the discrete property of the surface layer.Therefore Sun [17] presented a rutting prediction model of CBP by considering factors such as the magnitude of loads, the bearing capacity of structure (resilient deflection), the thickness of block layer, and the number of load repetitions, while the rut depth is relative to both the accumulation of permanent deformation and the lateral distribution of vehicle loads.Based on the regression analysis of test results, the following equation is used for the estimation of the rut depth: where RD is rut depth (mm),  0 is resilient deflection (mm),  0 is load pressure (MPa),  0 is radius of loaded areas (cm), ℎ is thickness of block (cm),  is number of load repetitions, and  is wheel configuration coefficient;  = 1 for single wheel and  = 1.17, 1.25, and 1.34 corresponding to ℎ = 8, 10, and 12 cm for dual wheels.
The rutting model of Huurman and Boomsma [19] is an empirical model reflecting the properties of materials suited for Netherlands.However, Sun's model can predict rutting by calculating the deflection of the pavement structure under diverse conditions for incorporation at the design stage.The allowable design rut depth for CBP with a flexible base varies depending on the location of construction and its function.Allowable rut depth for ports in China is 30 mm [17] and that for motorways is 15 mm [20] in Netherlands and 35 mm [21] in Japan.

Methodology to Develop Deflection Prediction Model
3.1.Laboratory Loading Test.As shown in Figure 1, the loading test was conducted in a 1,500 mm × 1,500 mm × 700 mm test pit.In this test, a section with 80 mm blocks for the surface layer, 30 mm for bedding sand, 200 mm aggregate layer, and 300 mm subgrade was used.In order to avoid the loss of bedding sand, a geotextile was placed between the base and bedding sand.Three construction patterns considered in most of the CBP design guide (ASCE 2015) [7,21] were assessed: stretcher, basket weave, and herringbone bond, which were the same test variables in Mampearachchi and Senadeera's study [16].The blocks used in this test were 200 mm × 100 mm × 80 mm rectangular blocks.Flexural strength and dynamic elastic modulus tests were in accordance with ASTM C 78 and ASTM C 215 [22,23] and were found to be 4.77 MPa and 12.09 GPa, respectively.The bulk density of the blocks was 2.24 kg/mm 3 .The particle size analysis was performed for the joint sand, bedding sand, and permeable base material [24].Figure 2 indicates the results of the passing ratio analysis of the base material, bedding sand, and joint sand; the passing ratios of the bedding sand and joint sand met ASTM C 33 and ASTM C 144, respectively, and the base material was uniformly distributed in aggregate sizes ranging from 40 mm to 2 mm [25,26].
A light falling weight deflectometer (LFWD) model KFD-100A produced by Tokyo Sokki Kenkyujo Co., Ltd., was used to investigate the load-deflection behavior.A 9000 N load was applied on a 150 mm radius plate in the center of the test pit.In order to identify the scope of the deflection caused by loading, deflection sensors were installed which can measure deflection of up to 1 mm at the loading point and at spots 200 mm and 500 mm away to right and left from the loading point.

Structural Analysis Model.
A three-dimensional (3D) finite-element analysis of the CBP using ABAQUS 6.10 [27] was performed with C3D8R as the element used in the analysis.Figure 3 shows the element configuration and system for the structural analysis considering three layers, which were the block surface layer, bedding sand layer, and base layer in an area of 1.5 m × 1.5 m.The blocks were connected with a joint sand element.The joint width was fixed to 3 mm and the sides of the test pit were in roller boundary condition to create the same condition as in the loading test.All of the elements were modeled as solid-deformable bodies having homogeneous and elastic property.The model was subjected to a linear static loading analysis.The subgrade is defined as a set of spring elements beneath the base layer to simulate the subgrade reaction by using the spring coefficient .The equidistant meshes were used for different analytical factors.The blocks, joint sand element, and layers of the pavement were fully bonded.The loading on the concrete blocks was applied uniformly on the circular steel plates with a diameter of 300 mm.
Input values from prior researches were used to select the block size, soil support condition, block thickness, construction pattern, loading scale, and contact area as factors affecting the vertical behavior of CBPs [17].Yaginuma et al. [28] generalized block shapes as a shape factor (ratio of total side area to top area, SF) to determine the proper sizes for a block that can be used for roads and argued that a larger SF leads to higher load transfer efficiency between blocks.The SF concept was designed to simulate block shapes.Five levels of SF used in the analysis were 2.3 to 2.7 at intervals of 0.1, calculated according to blocks in wide use in Korea.The ratio of length to width of the blocks was fixed at a 2 : 1.
Three construction patterns were simulated as the same with lab test.The elastic modulus of the base was assessed for four levels which were 100 MPa, 300 MPa, 500 MPa, and 740 MPa.Further, subgrade reaction was considered at three levels which were high (300 pci), medium (200 pci), and low (100 pci).Three different block thicknesses were used in the model, which were 60 mm, 80 mm, and 100 mm blocks.The base was built from 150 mm to 300 mm at intervals of 50 mm.Since previous studies had reported that the elastic modulus of the blocks had an insignificant effect on pavement performance, this was not considered in the analysis [29].
Input data for the 3D FEM analysis include load, the material properties and thickness of the pavement layer, and the soil support value.Table 1 summarized the properties of the materials employed for the structural analysis.The elastic modulus of the block was evaluated by the laboratory test; therefore the elastic modulus of the joint and bedding sands was used as presented in previous studies [15].

Load-Deflection Behavior
To evaluate the load-deflection behavior, 10 times of the test loads were dropped for each experimental variable.and the stretcher bond showed the largest in maximum deflection.Figures 4-6 show the comparison results of deflection curves from FEM analysis and lab test for each construction pattern.Both the test results and analysis results in a 9,000 N load condition appeared similar in all construction patterns, but deflection basin differed among construction patterns.By considering maximum deflection, the stretcher bond had an estimated 1.2 mm of deflection, and the value fell by 25% to 0.9 mm in the herringbone bond with a noticeably different deflection curvature.The curvature is proportional to stress; it was presumed that the stress levels are from high to low, followed by stretcher, basket weave, and herringbone.
In the FEM analysis, tensile stress on the base ranged from 60 to 100 kPa depending on the construction patterns in the condition of a 9,000 N load with a 300 mm radius.The heavier load led to a stronger effect on the deflection of the patterns.Because the lab test results and analytical results were matched well, the analytical model was set up to estimate the deflection for the remaining variables such as construction patterns and block shapes.

Multiregression Analysis of Deflection
Although the deflection obtained from a 3D FEM analysis can be reflected into the design as it is, determining the value   under various conditions and conducting a regression analysis reduces design time, which is a recent trend in pavement design [30].In the load distribution model for CBP, the maximum deflection obtained from the analysis was considered as resilient deflection and a dependent variable in a regression analysis using PASW Statistics 18.Among the independent variables were shape factor (SF), contact area ( 0 ), applied load ( 0 ), subgrade reaction coefficient , pattern ( ℎ ), block thickness (ℎ 1 ), base thickness (ℎ 2 ), and elastic modulus of the base (  ).Since construction pattern is on a nominal scale that is difficult to use in a regression analysis, a binary variable was used in values "0" without applying a herringbone pattern and "1" with applying a herringbone pattern.A total of 200 FEM analyses were performed with the variables range of SF of 2.3 to 2.7, traffic loading of 0.125 MPa to 0.6 MPa, base thickness of 20 to 30 cm, and block thickness of 6 to 10 cm.The estimation method for the multiregression equation was using a simultaneous input to examine the influence of specific independent variables with the other independent variables controlled.The results of the multiregression analysis, excluding variables rejected by significant probability, are summarized in Table 3.This multiregression analysis model can explain 91.5% of the deflection caused by the independent variables.The Durbin-Watson value verifies the satisfaction of the independence of error.In Table 4, which shows the coefficients from the multiregression analysis on deflection prediction, the significance probability (sig.) for all of the dependent variables except pattern 3 is less than 0.05, which statistically fits the regression model.Although the value is bigger than 0.05, the significance probability (sig.) of pattern 3, which is 0.07, also fits the regression model.
The equation to calculate resilient deflection considering the various block shapes, thicknesses, and patterns is given as where  0 is maximum deflection (mm), SF is shape factor,  0 is diameter of grounding area (mm),  0 is applied load (N),    ℎ 1 is block thickness (mm), ℎ 2 is base thickness (mm),  ℎ :  ℎ = 1 for herringbone pattern and  ℎ = 0 for nonherringbone pattern, and   is elastic modulus of base (MPa).
Deflection as a dependent variable increases or declines with  value as the values of the independent variables change.Block shapes and construction patterns influence changes in the value of the deflection.If load  0 increases, deflection increases by 5.356 − 5, and as the value of additional independent variables grows, deflection drops.Figure 7 compares the deflection from the analysis and the value calculated in the model.As deflection grows, the model appears unstable, which can be utilized in designing.

Application of Rutting Prediction Model
The developed deflection prediction model was implemented at the CBP rutting model suggested by Sun [17] in (3).APT results were used to verify the validity of the rutting model of Sun [17].The specific methodology for the APT and findings are presented in Lin et al. [31].
The test sections were consistent with an 80 mm surface layer laying herringbone pattern, a 30 mm sand layer, and 200 mm/300 mm thickness base layers.Resilient deflection was measured with FWD, which was applied to the rut depth prediction model to verify the validity of the developed equation.With a load condition of 40 kN, 0.871 mm and 0.736 mm of deflection occurred in cross sections with a base thickness of 200 mm and 300 mm, respectively.The rut depth was evaluated in five phases depending on the number of load repetitions for each pavement in the transverse direction.APT was operated at 100,000 cycles, representing 1.08 million total ESALs.The rut depth is shown in Figures 8 and 9.Under the total load of 1,081,000 ESALs, the rutting depths of block pavement with a granular aggregate base of 30 cm and 20 cm were 21 mm and 25 mm, respectively.
The accumulated rutting measured in the APT and the rut depths calculated in Sun's model are shown in Figure 10, with a similar trend of rutting caused by growing traffic load; the rut depth calculated using the developed deflection estimation model was the largest, and the result obtained from applying the FWD deflection data to Sun's model was the smallest.The results have shown that the prediction model was reasonable from a more conservative perspective, and the error in terms of methodology was small as 5 mm of rutting in 1,000,000 ESALs.Therefore, the proposed model is adequate to be applied in the design.To assess rutting by the thickness of the block and base, block shapes, and construction patterns, (4) was applied to Sun's [17] model for sensitivity analysis.Block thicknesses used were 60 mm, 80 mm, and 100 mm; base thicknesses were 200 mm, 250 mm, and 300 mm; and block shapes of SF were 2.2, 2.4, and 2.6.The contact radius was 150 mm.A 40 kN load was applied to a dual wheel.Figure 11 depicts differences in the rut depth depending on block thickness.During the early age of service, rut depth rose sharply with all variables.Rutting varied largely as the traffic load increased.As block thickness increases from 60 mm to 100 mm, rut depth dropped by 25%.Laying in a herringbone pattern was set to be relatively lower rut depth up to 2 mm.These findings are similar to those presented by Hiroshi et al. [32].
Figure 12 shows differences in rut depth depending on base thickness with block thickness fixed to 80 mm.As base thickness rose from 200 mm to 300 mm, rut depth decreased by 22%.Laying blocks in a herringbone pattern decreased rut depth by up to 2 mm, as it did when the block thickness was used as a variable.Increasing the block thickness by 20 mm or the base thickness from 50 mm to 100 mm also resulted in similar performances.Figure 13 depicts differences in rut depth depending on block shapes (SF) with the block and base thickness fixed to 80 mm and 200 mm, respectively.As shown in Figure 13, the difference in rut depth depending on block shape was Advances in Materials Science and Engineering   insignificant; rut depth was found to be affected more by the construction pattern than by block shape.These findings are similar to those from a field test [16].It is inconclusive, however, that construction pattern is the only determinant in load distribution effect.Because CBP performs as a complex interlocking mechanism of horizontal, vertical, and rotational behaviors, it seems unreasonable to evaluate its performance by separating construction patterns and block shapes too strictly.

Conclusions
In this study, a resilient deflection prediction model was developed to estimate rutting of CBP.Diverse sets of block shapes and construction patterns were assessed using FEM analysis and dynamic loading test.LFWD was used for the loading test to evaluate correlations between loading and deflection depending on block shape and construction pattern.In order to verify the validity of the test results and develop a deflection prediction model for CBP, 3D FEM modeling simulations were performed using ABAQUS 6.10 [27].
In addition, a resilient deflection estimation model was developed for CBP that considers construction patterns and block shapes via a multiregression analysis.Among the independent variables of the model were contact radius ( 0 ), load ( 0 ), block shape and construction pattern, block thickness (ℎ 1 ), base thickness (ℎ 2 ), and elastic modulus (  ).The deflection prediction model for CBP was applied to the existing rutting model and results show that as block thickness increases from 60 mm to 100 mm, rut depth dropped by 25% with the load repetition increasing to 1.2 million.By increasing the thickness of base from 200 mm to 300 mm, rut depth could be reduced by 22%.Increasing the thickness of the surface layer by 20 mm has similar rutting resistance performances by increasing the base thickness from 50 mm to 100 mm.It was found out that the block shape did not have significant effects on load distribution, but construction pattern was a relevant parameter.Laying in a herringbone pattern was set to be relatively lower rut depth than the other up to 2 mm.However, it is inconclusive that construction pattern is the only determinant in load distribution effect.Because CBP performs as a complex interlocking mechanism of horizontal, vertical, and rotational behaviors, it seems unreasonable to evaluate its performance by separating construction patterns and block shapes too strictly.

Figure 1 :
Figure 1: Specification of dynamic laboratory loading test.

Figure 2 :
Figure 2: Passing ratios of sand and base material.

Figure 4 :
Figure 4: Comparison of deflection from the analysis and test (stretcher).

Figure 5 :
Figure 5: Comparison of deflection from the analysis and test (basket weave).

Figure 6 :
Figure 6: Comparison of deflection from the analysis and test (herringbone).

Figure 7 :
Figure 7: Comparison of analysis and calculated deflection results.

Figure 10 :Figure 11 :
Figure 10: Difference in rutting depth between the APT data and the model.

Figure 13 :
Figure 13: Difference in rut depth depending on block shapes.

Table 2
summarized the average, standard deviation, and normalized value calculated by load-deflection data obtained in dynamic loading test.The herringbone bond pattern had the lowest

Table 1 :
Properties of materials used for the structural analysis.

Table 2 :
Summary of average, standard deviation, and normalized deflection for each pattern.

Table 4 :
Coefficients of multiregression analysis for deflection.