Abstract
This paper proposes an incipient assessment method for postconstruction settlement of highway subgrade by using statistical analysis approach with an improved static loading test method, time-extended loading test. The time-extended loading residual settlement Sr(t) and settlement rate Vs(t) are calculated from the field test data. The test procedure and the corresponding experimental validations are presented. Based on the field test data, a probabilistic model is built to bridge the final residual settlement and the actual postconstruction settlements. The proposed time-extended loading test based assessment system has been validated through three operating and well-documented road sections. According to the validation results, the construction quality can be accurately evaluated at an early stage and the corresponding remedial measures can be applied timely.
1. Introduction
Traditional incipient quality inspection methods for piled-soil composite subgrades of highways are not accurate enough. The conventional assessment system sometimes cannot reveal the true construction quality. Large postconstruction settlement can be detected and bridge bumps can be observed in lots of the road construction projects whose composite subgrades satisfied the traditional incipient quality inspection requirements [1–3]. Hence, it is necessary to develop accurate incipient assessment method for evaluating the postconstruction settlement with the consideration of the behavior of settlement process. An accurate assessment index ensures the safety and comfort of the vehicles and the passengers during the operation period of the road [4–8]. Many scholars and technicians used theoretical analysis methods [9–11], numerical simulation approaches [12–15], or field measurement data fitting methods [16–18] to study the postconstruction settlement of highways. According to the soil modulus amplification factor for creating the soil deformation, Zhou Tonghe et al. [19] proposed a theoretical method for calculating postconstruction settlement of composite subgrades. Ding Mingji [20] analyzed the relationships between the design parameters of the Cement Fly-Ash Gravel (CFG) piles and plates and the postconstruction settlement of a high-speed railway by using numerical simulation method. Kermani et al. [21] proposed two methods for estimating postconstruction settlement: using subgrade height to estimate postconstruction settlement; using stress-strain behavior of compacted layer during construction to predict the stress-strain-time change of subgrades. Chen Yuanhong et al. [10] proposed a power polynomial prediction method for postconstruction settlement of soft foundation embankment based on field monitoring data. Advanced heuristic algorithm based estimation methods [22, 23] have also been proposed for tackle problem of the multi-input impact to the settlement estimation. However, seldom of the study in the literature focuses on the investigation of the early assessment method for the postconstruction settlement of composite subgrades by using constructive quality inspection methods. At the same time, certain deficiencies exist in the traditional evaluation methods. For example, the correctness of the system parameters affects the accuracy of the theoretical analysis and numerical simulation results. The time cost of the long-term field monitoring method is expensive and it is hard to control the precision of the monitor location.
To ease the deficiencies of the traditional methods, a time-extended loading test method acted as an improvement of the static load test has been proposed in this paper. Based on the field measurement from the time-extended loading test, the incipient assessment indexes have been developed through statistical analysis approaches. The prediction method of the postconstruction settlement has been provided as well. Instead of replacing the existing evaluation method, the proposed method aims at improving the existing method as well as enriching the existing quality inspection index system. In addition to strength inspection, the proposed evaluation method can perform strength-settlement detection. Besides, early remedial actions can be implemented based on the proposed assessment indexes to reduce postconstruction settlement in the early stage.
The main contributions of this paper can be concluded as follows: a novel time-extended loading test is proposed in this paper for minimizing the construction condition impact and enhance the accuracy of the prediction results; a probabilistic model for predicting the postconstruction settlement by using time extend loading test results is developed by using 38 field data with diverse construction conditions; the proposed method and indexes have been validated through three operated road sections.
2. Time-Extended Loading Test Method
The time-extended loading (TEL) test is a supplementary test for accurately evaluating the postconstruction settlement of pile-soil composite subgrade. Apart from the approach that using the elastic cushion in the experiment mimics the functions of the cushion in practice, the additional static weight has been loaded on the cushion and the settlement changes will be recorded.
2.1. Necessity of Time-Extended Loading Test
Cushion plays a critical role in the bearing system of the pile-soil composite subgrade. Because of the mass of the subgrade, the vehicles, and the dynamic adjustment from the cushion, the stress stabilization process of the pile-soil composite subgrade generally takes a long time. Hence, the stress state of the pile-soil composite subgrade cannot reach a stable state under the conventional static loading test. Because of this deficiency of the conventional static loading test, the estimation of the postconstruction settlement may not be accurate enough. With the support from the Guangdong Province Transportation Engineering Quality Supervision Station and the Foshan Highway Bridge Engineering Monitoring Station, the authors have collected the postconstruction settlement data from 157 different high-standard roads in Guangdong area. All the monitored high-standard roads satisfied with the requirements in the static loading test and operate for 2–5 years. The pile types of sample roads include pipe piles, cement mixing piles, CFG piles, and powder jet piles. Meanwhile, the traffic flows, geological composition, and fill height of each road section are also different. Hence, the settlement data of the selected road samples have certain representativeness. The hyperbolic curve method has been used to fit the field postconstruction settlement data of each section for predicting their postconstruction settlement after 15-year operation. The corresponding prediction results are concluded in Figure 1. The predicted 15-year postconstruction settlement data are considered as the actual settlement data in this paper and denote as Sf.

According to Figure 1, the maximum and minimum postconstruction settlements are 321.15 mm and 8.78 mm, respectively. 54.78% of the total road sections have more than 100 mm postconstruction settlement. 30% of the road sections fail to satisfy the regulations. In conclusion, the traditional evaluation method, static loading test method, cannot provide an accurate prediction on the final postconstruction settlement since the stress state of the pile-soil composite subgrade is not stabilized during the static loading test and the model of the cushion is unrealistic. Therefore, it is possible to improve the detection and evaluation accuracies of the composite subgrade postconstruction settlement by extending the loading test time and obtaining settlement information under near steady-state condition. To approach the near steady-state condition of the composite subgrade, the TEL test based assessment method is proposed and discussed in the following sections.
2.2. Time-Extended Loading Test Design and Procedure
As mentioned in the previous section, the stress ratio of the pile-soil can be used to identify whether a pile-soil composite subgrade is in a stable stress state. The stress ratio of the pile-soil tends to be stable when the subgrade is loaded. In the literature, many scholars and engineers have carried out some studies [24–27] related to the behavior of the pile-soil stress ratio and they have obtained useful conclusions. Based on the observation from the field test results and corresponding theoretical analysis, the adjustment rate of the cushion can be increased when a composite subgrade is statically loaded by two times of its designed bearing capacity. Under this situation, a nearly stabilized pile-soil stress ratio under rated load can be achieved within a short time. The similar phenomenon has been observed in 90% of the author's field tests: the pile-soil final stress ratio in rated load condition of a composite subgrade can reach within 48 hours by statically loading the composite subgrade with two times of its designed bearing capacity. Therefore, the loading value in the TEL test is selected as two times the rated load and the extended test time is 48 hours.
After determining the weight and time of the TEL test, the detailed test procedures can be concluded as four parts: prepreparation, loading, data acquisition, and observation of the settlement-time relations.
Loading value determination and classification: first, the design value of the composite subgrade bearing capacity was determined. The experiment loading value was two times of the design value and it equals the maximum loading value of the single pile composite subgrade in the conventional static load test. According to the requirements in the specific regulations, the loading process can be divided into 8-12 levels, and the loading time interval of each level was controlled.
Test cushion selection: to emulate the adjustment effect of a cushion accurately, 8-10 cm sand cushions combined with an elastic rubber cushion were used. The elastic rubber cushion thickness h and the elastic modulus E are provided in Table 1. Usually, the cement mixing piles, Dry Jet Mixing Piles, and high-pressure jet grouting piles are classified as flexible piles [28, 29]. CFG piles are classified as semirigid piles [30, 31]. Pipe piles and bored piles are classified as rigid piles [32, 33].
Installation of loading devices: the test used a pressure platform with counterforce setting. The pressure platform acted as a load counterforce, and the load greater than the maximum testing load was added to the platform before the beginning of the TEL test. The hydraulic jack was used for control the test loading level. A circular plate with a diameter of 1.26 m and an area of 1.25 m2 was adopted as the load plate; see Figure 2.

(a)

(b)
Data acquisition system layout: four dial-gauges were installed at the four corners of the platen. The settlement acquisition frequency was determined based on the test requirements. The accuracy of the dial-gauges is 0.01 mm.
Observations of settlement and time relationships: during the loading phase, the data was recorded according to the requirements of the standards. The first six reading time intervals of the settlement in the TEL test are 10 minutes, 10 minutes, 10 minutes, 15 minutes, 15 minutes, and 15 minutes. Then the settlement was measured every 30 minutes.
Data analysis: from the collected test data, the corresponding settlement values at each time point during the TEL test were extracted. The detailed calculation and analysis processes are shown in Section 3.
2.3. Field Test of Time-Extended Loading Test
The proposed TEL method has potential to minimize the construction impact to the estimation of the postconstruction settlement. For validating this property of the proposed TEL method, the single pile composite subgrade field tests with 38 diverse pile types and under different test conditions were performed. The 38 test sites are selected based on the static loading test results of hundreds of projects in Guangdong Province. The pile types can be concluded as mixing cement piles representing flexible piles, CFG piles representing semirigid piles, and prefabricated pipe piles representing rigid piles. Tests nos. 1-12 are flexible piles, tests nos. 13-26 are semirigid piles, and test nos. 27-38 are rigid piles. All the settlement-time data during the TEL test is recorded for further analysis purposes.
As mentioned in the Section 2.1, the pile-soil stress ratio after applying the proposed TEL test is close to the pile-soil stress ratio under steady-state condition. To validate this argument, an example field test comparison is shown in Figure 3.

(a)

(b)
According to Figure 3(a), the stress ratio of a road section after the second year operation was 8.9 and it tended to stable around 9. According to Figure 3(b), the stress ratio after TEL test in the same road section was 9.21 which was close to the field measurement data. Hence, the proposed TEL test can effectively shrink the settlement process and minimize the impact of the construction condition to the prediction of the final postconstruction settlement.
3. Time-Extended Loading Test Based Assessment System
The proposed time-extended loading test based assessment system is a two-phase estimation approach. By using the field measurement data from TEL test, the postconstruction settlement can be estimated first by using the hyperbolic curve fitting method (final residual settlement). Then, a probabilistic model is built as a bidirectional mapping between the final residual settlement and the true postconstruction settlement. By using the proposed two-phase estimation process, the construction condition difference impact on the postconstruction settlement has been considered by the first phase, hyperbolic curve fitting process. Hence, the probabilistic approach in the second phase acts as a general model without considering the construction condition differences.
3.1. Assessment Indexes of Time-Extended Loading Test
The primary objective of this section is to develop the assessment indexes for evaluating the final postconstruction settlement implicitly. With the consideration of conventional field testing results and the cushion behaviors in practice [34, 35], the hyperbola is selected for representing the relationship between settlement S(t) and time t during the TEL test. Mathematically, we haveFormula (1) can be rewritten aswhere t is the cumulative time; t0 is the start time of the TEL test; S(t0) is the settlement after the static loading test; S(t) is the cumulative settlement at time t; a and b are the linear regression results computed from the TEL test to represent the relationship between after t0.
To obtain the final residual settlement, we havewhere Srf is the final residual settlement after the TEL test; Sr(t) is the residual settlement.
Meanwhile, the settlement rate Vs(t) can be obtained from the formula (1) by taking the derivative with respect to the time t and we have When t equals t0, the initial settlement rate Vs(t0) is .
Three field test results are analyzed in this subsection as examples for the illustration of the calculation process of the assessment indexes. The field test conditions are shown in Table 2; the selected test results with observation time are shown in Table 3.
According to (2), let , and . Based on the settlement data obtained from the TEL test, the relationship between y and x can be obtained; see Figure 4. The linear regression results for all three field test results are shown in Table 4. The parameters shown in the relationship between settlement S(t) and time t during the TEL test can be written as The settlement rate Vs(t) can be calculated from (4) and the Sr(t) can be obtained from (3). The final settlement can be obtained from (5). The detailed calculation results can be found in Table 5.

Since the final residual settlement Srf and the settlement rate Vs(t) have directly related to the actual final postconstruction settlement Sf. The Srf and Vs(t) can be used as the indicators to estimate the postconstruction settlement for the determination that whether the postconstruction settlement of a composite subgrade satisfies the specification engineering standards. The corresponding mathematic expressions are where td is the end time of the TEL test; Srt and Vrt are the assessment criteria, which usually can be obtained by using the 15-year cumulative settlement requirements of different road sections. Typically, the 15-year cumulative settlement requirements of a bridgehead transition section, a section with structures, and a general road section are 100 mm, 200 mm, and 300 mm; respectively. Different from directly using 15-year cumulative settlement requirements stated above, a probabilistic model has been built for linking the final postconstruction settlement requirements to the requirements on the final residual settlement.
3.2. Probabilistic Analysis of Time-Extended Loading Test Result
Different from the traditional evaluation system, the statistic properties of the assessment indexes Srf and Vs(td) are analyzed in this section. The corresponding assessment criteria Srt and Vrt are obtained by using the field monitoring data from 157 road sections and 38 TEL test results which are presented in the previous section. The data preprocess results for the 38 field tests are attached in Supplementary Materials.
3.2.1. Characteristics Analysis for Postconstruction Settlement Distribution
To obtain the assessment criteria, the probabilistic characteristics of the postconstruction settlement should be studied first. According to the predicted postconstruction settlement data Sf shown in Figure 1 and the state of art studies, a lognormal distribution function has been selected to represent the probabilistic distribution of the postconstruction settlement [36]. The detailed steps for obtaining the corresponding statistical parameters in the lognormal distribution function are presented as follows.
According to the assumption, if predicted postconstruction settlement Sf of the composite subgrade follows the lognormal distribution, we have . Hence, the probability density function (PDF) ff(sf) of the postconstruction settlement is An intermediate variable z is defined to convert the lognormal distribution to standard normal distribution, and we havewhere .
Since z obeys the standard normal distribution, the cumulative distribution function (CDF) Ff(sf) of the postconstruction settlement can be obtained asThe maximum likelihood estimations for and arewhere and are estimated expectation and variance, respectively; n is the number of the postconstruction data; is the postconstruction settlement value of the sample.
By plugging in the postconstruction settlement data into (9) and (10), we have and The calculation result shown above has been validated through nonparametric hypothesis test (Chi-square Fit Test). The hypothesis H0 is , where is the CDF of the hypothetical lognormal distribution. To perform Chi-square Fit Test, the sample space (settlement value) has been discretized equally into K = 17 subintervals. The length of each intervals is 20 mm. According to the postconstruction settlement data, the minimum and maximum postconstruction settlement in the sample sequence are 8.78 mm and 321.15 mm, respectively. Hence the subinterval can be represented as, ,…,,. Based on the proposed data classification, the statistical divergence degree v can be calculated aswhere nk is the number of data falling to the interval , n is the total number of samples, and is the probability that the settlement belongs to the interval , which can be approximated from the following calculations: Based on (12) and the postconstruction settlement data , we have . Meanwhile, the degree of freedom f is K – R – 1 = 14 where R is the number of unknown parameters. According to the Chi-square distribution table, we have , which means the possibility that the distribution of the postconstruction settlement which obeys the estimated lognormal distribution is larger than 90%.
Therefore, the probability density distribution of the postconstruction settlement obeys the lognormal distribution with the parameters that and .
3.2.2. Analysis of Characteristics of Final Residual Settlement Distribution
Based on the TEL test measurements and the corresponding calculation results shown in the previous section, the static load test settlement S(t0) and the final residual settlement Srf of the 38 different piles under different test conditions are plotted in Figure 5.

According to Figure 5, the maximum settlement in the static load test and the maximum final residual settlement in the TEL test are 95.32 mm and 38.76 mm, respectively (no. 24). The minimum initial settlement and final residual settlement are 3.8 mm and 0.299 mm, respectively (no. 1).
To evaluate the probabilistic property of the final residual settlement Srf, three distributions have been examined by using Chi-square Fit Test: lognormal distribution, Γ-distribution with shape parameter larger than 1, and inverse Γ-distribution with shape parameter larger than 1.
(1) Lognormal Distribution. Assume the final residual settlement Srf is a random variable with a lognormal distribution, namely, . The statistical parameters of the Srf can be computed based on the proposed method in Section 3.1, yielding we have and . The cumulative probability function of the final residual settlement Srf isSimilarly, the Chi-square Fit Test has been used to validate the results. The hypothesis H1 is , where is the cumulative lognormal probability function. Since the largest final residual settlement is 38.76 mm and the smallest Srf is 0.299 mm, the sample space has been discretized unevenly into K = 7 subintervals. The subinterval are , , , , , , . According to (12), the statistical divergence degree v is 5.547. The degree of freedom f is K - R - 1 = 4 where R is the number of unknown parameters. According to the Chi-square distribution table, we have meaning the possibility of the hypothesis H1 is larger than 90%.
(2) Gamma (Γ) Distribution. We can also assume the PDF of the final residual settlement Srf follows gamma (Γ) distribution and the shape parameter is greater than 1, namely, . According to [37], the shape and scale parameters and can be calculated and we have = 0.805 and the scale parameter = 10.05. Hence, the assumption is not satisfied since is smaller than 1. As a result, the PDF of the final residual settlement Srf does not obey the gamma (Γ) distribution.
(3) Inverse Gamma (Γ) Distribution. According to [38], the PDF of the final residual settlement Srf may follow an inverse gamma (Γ) distribution with a shape parameter greater than 1, yieldingThen the PDF and CDF of can be represented asThe maximum likelihood estimations for and arewhere and are the estimated statistical parameter of the inverse Γ distribution; is the final residual settlement of sample i. Based on field measurements, one can obtain and .
To validate the calculated parameters, the hypothesis H2 is , where is a CDF of inverse gamma distribution function with parameter and . Similarly, sample space has been discretized unevenly into K = 8 subintervals: , , , , , , , . The divergence degree v is 22.38. Based on the distribution table, we have which means the possibility of the hypothesis H2 is larger than 99.5%.
According to the Chi-square Fit Test results shown in the previous section, the PDF of the final residual settlement Srf follows the lognormal distribution with the parameters and or follows the inverse Γ distribution with the shape parameter and the scale parameter .
3.3. Assessment Criteria
3.3.1. Assessment Criterion Calculation
Based on the analysis in the previous section, the postconstruction settlement Sf is positively correlated with the final residual settlement Srf. The relationship between Srf and Sf can be generalized to the probabilistic space. Probabilistically, we havewhere Sft is the assessment criterion for the postconstruction settlement. Namely, the probability that the postconstruction settlement is smaller than its assessment criterion Sft is equal to the probability that the final residual settlement is smaller than its assessment criterion Srt. Based on this relationship, the assessment criterion value of the final residual settlement can be obtained by inverse calculation, yieldingwhere is the inverse function of the CDF . It worth pointing out that (21) can also be used to predict the postconstruction settlement.
When the allowable postconstruction settlement Sft of transitional bridge sections, structural sections, and general sections are less than 100mm, 200mm, and 300mm, the corresponding probability Ff(Sft) can be obtained from the standard Gaussian distribution table. Based on the probability value, the corresponding Srt can be obtained by using (15) and (18). Hence, the assessment criterion relationships between the final residual settlement and the postconstruction settlement are obtained; see Table 6.
Based on the similar idea, the assessment criterion for the settlement rate can be obtained. To calculate Vrt from (4), the CDF of the parameter a is necessary. The experimental cumulative distribution function of the parameter is plotted in Figure 6.

The exponential function is used to fit the experimental CDF curve of the parameter a, yielding Similar to (22), we havewhere is the inverse function of Fa(a); art is the assessment criterion for the parameter a. Based on the calculation results shown in Table 6 and (24), (3), and (4), the assessment criterion for the settlement rate Vs(td) in the TEL test can be obtained; see Table 7.
It worth pointing out that the purpose of using two assessment indexes is to enhance the reliability of the assessment results.
4. Example Analysis and Discussion
To validate the proposed assessment indexes and criteria, three representative experiment examples are selected and discussed in the following paragraphs.
4.1. Example Analysis
Three road sections with relatively complete settlement measurements were selected and presented in this section. The specific design and construction parameters are concluded in Table 8. The settlement measurements collected during the TEL test are presented in Table 9.
No. 1 and no. 2 road sections were completed and opened to traffic in July 2015; no. 3 road section was completed and opened to traffic in December 2016. The monitoring frequencies of the postconstruction settlement are shown in Table 10. The postconstruction settlement data are shown in Figure 7.

The final residual settlements in the TEL test for three target subgrade sections are collected and shown in Table 11. Meanwhile, the 15-year postconstruction settlement calculated by using hyperbolic curve fitting method and probabilistic methods presented in this paper is given in Table 11. The corresponding prediction errors are shown in Table 12
According to calculation results from Tables 6, 7, 11, and 12, the lognormal distribution evaluation result has better accuracy than inverse Γ-distribution when the final residual settlement is smaller than 4.5 mm or larger than 13.3 mm. The inverse Γ-distribution prediction has better results when the final residual settlement is smaller than 13.3 mm or larger than 4.5 mm. Hence, the construction quality can be evaluated by using the proposed novel assessment method.
In Assessment Indexes Calculation section, the first two road sections, namely, the field tests of the composite cement subgrade of the cement mixing pile and the CFG pile composite subgrade I, satisfied the proposed assessment criteria for both Srf and Vs(td). On the contrary, the field test results for the CFG pile composite subgrade II did not meet the proposed criteria. Based on this evaluation, specific remedial actions have been implemented including grouting, compaction, and prolonged preloading. Because of the probability based assessment indexes and the corresponding remedy actions, this road section has operated three years and no significant postconstruction settlement has been found.
5. Conclusion
A novel supplementary test, time-extended loading test, is proposed in this paper for minimizing the construction condition impact on the settlement estimation and increasing the accuracy of the prediction. The corresponding probabilistic model is built for establishing the TEL based assessment system. Compared with the existing estimation methods, the proposed TEL method only relies on the field test data of the first 48-hour measurements. Also, the corresponding probabilistic model further improves the accuracy of the classic hyperbolic curve fitting method by mapping the final residual settlement to the postconstruction settlement. Hence, the proposed method can identify the quality of road construction in a very early stage compared with other existing methods. Four major conclusions are listed as follows.
The settlement and settlement rate in the time-extended loading test can be obtained by using a hyperbola model. The final residual settlement and the settlement rate after 48 h can be used as the assessment indexes for predicting the final postconstruction settlement.
According to our statistical analysis results, the final residual settlement less than 4.5 mm or greater than 13.3 mm obeys the lognormal distribution. The final residual settlement greater than 4.5 mm and less than 13.3 mm follows the inverse gamma distribution.
The assessment criteria for the final residual settlement of the bridge-to-bridge transitional section, structural section, and general section corresponding are 4.5 mm, 13.3 mm, and 35.41 mm. The corresponding 48 h settlement rate criteria are 0.021 mm/h, 0.046 mm/h, and 0.055 mm/h.
For the pile-soil composite subgrade that does not satisfy the proposed assessment criteria, follow-up remedial measures, such as grouting, compaction, and prolonged preloading, should be applied in the road section based on the actual situation. According to the field experimental results, the follow-up remedial measures help to avoid the occurrences of the severe postconstruction settlements.
Data Availability
The data used to support the findings of this study are included within the supplementary information file (available here).
Disclosure
This work is an original one conducted at Hohai University.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Authors’ Contributions
All the authors have contributed to and approved this paper.
Acknowledgments
This work was supported by Postgraduate Research & Practice Innovation Program of Jiangsu Province [KYZZ15_0146], National Science Foundation for Young Scientists of China , the Fundamental Research Funds for the Central Universities , Science and Technology Project of Guangdong Provincial Transportation Department -02-013; 2015-02-01, and the China Scholarship Council (no. 201806715014).
Supplementary Materials
The 38 field test data during the TEL test are shown in Table A. (Supplementary Materials)