Research Article  Open Access
Xiaoping Ji, Jia Li, Zhifei Cui, Shouwei Li, Yue Xiong, Jianming Hu, Yingjun Jiang, "A Prediction Method for the California Bearing Ratio of SoilRock Mixture Based on the Discrete Element Method and CT Scanning", Advances in Civil Engineering, vol. 2020, Article ID 9794756, 12 pages, 2020. https://doi.org/10.1155/2020/9794756
A Prediction Method for the California Bearing Ratio of SoilRock Mixture Based on the Discrete Element Method and CT Scanning
Abstract
Because of the large amount of gravel with particle sizes over 40 mm in the soilrock mixture (SRM), it is impossible to determine its California Bearing Ratio (CBR) via the indoor test method, which is a key parameter for designing the backfill in underground mined cavities or the road subgrade constructed with SRM. In this paper, Xray computed tomography (CT) scanning and 3D reconstruction technology were used to construct the 3D structure of SRM particles with a particle size greater than 5 mm. Based on the vertical vibration test method (VVTM) and PFC^{3D}, the numerical simulation method (NSMCBR) of SRM was established. The CBR of the SRM with a maximum particle size over 40 mm (SRMG) was studied by NSMCBR, and the effects of factors such as maximum particle size, soil content, and largesize particle content (d ≥ 40 mm) on the CBR were investigated via NSMCBR. Based on the laboratory tests and NSMCBR, the prediction model and the determining method of CBR of SRMG were established and verified. The results show that the maximum deviation between the CBR obtained from NSMCBR and laboratory tests was 7.4%. The CBR of SRMG decreases linearly with the increase in soil content and increases with the increase in maximum particle size and largesize particle content. The practical project shows that the maximum deviation between the predictive and measured values of the CBR of SRMG was less than 1.5%, indicating that the prediction model and the method established in this paper have high reliability.
1. Introduction
The soilrock mixture (SRM) is a heterogeneous mixture of soil particles, rock particles, water, and pores, which is widely used in highway subgrade construction in mountainous areas [1]. Due to the complexity of material composition and the irregularity of structural distribution, the SRM presents the complex mechanical properties. The California Bearing Ratio (CBR) is an important parameter of the subgrade structure design, which can well express the ability of the subgrade to resist permanent deformation; the mechanical properties of SRM are closely related to the compaction method. Moreover, the laboratory testing of SRM can only be achieved by mimicking the practical conditions. It is not economical to drill cores from a road construction site for mechanical behavior testing. Therefore, SRM samples are compacted and formed in the laboratory with the same compaction level as in the field. The correlation of mechanical properties between the laboratoryproduced samples and field drilled cores is the key indicator for evaluating the quality of the laboratory compaction method. Throughout the world, forming CBR specimens is done mainly by three methods, namely, the static pressure compaction method, the shaking table method, and the surface compaction method. However, these methods of compaction do not conform to modern heavy rolling equipment, which results in low indoor standard density. The compactness can easily satisfy the requirements of the code, which causes the unsatisfactory compactness of the soilrock mixed subgrade and greatly affects the strength of the subgrade.
To overcome the shortcomings of the traditional methods, Chinese researchers fully simulated the effect of modern heavy rolling equipment and developed the vertical vibration test method (VVTM) based on the principle of the directional vibration roller. Compared with the traditional methods, the VVTM can reduce the optimal moisture content and increase the maximum dry density and strength of the mixtures. The correlation between the mechanical properties of an indoor test piece and an onsite core sample can reach 90%. Therefore, the VVTM has been widely used, especially in semirigid base material design [2–4]. However, the VVTM is only applicable to cohesive finegrained soil with particle sizes below 40 mm or noncohesive freedraining giantgrained soil and coarsegrained soil. Because SRM contains a large number of particles with particle sizes greater than 40 mm (SRMG), it is impossible to accurately determine its maximum dry density and forming specimens by the laboratory methods including the VVTM.
To solve the problem that the mechanical properties of SRM with maximum particle sizes above 40 mm (SRMG) cannot currently be tested in laboratory tests, the virtual experiments such as the discrete element modeling (DEM) are one of the most effective solutions [5–7]. The DEM can not only complete various experimental projects as in the real environment, but also the experimental results obtained can be equivalent to or even better than those achieved in the real environment. In addition, due to its high resolution and nondestructive nature [8, 9], the numerical simulation of SRM from the microperspective based on Xray computed tomography (CT) has gradually become the focus. Wu et al. [10] obtained the threedimensional structure of sand via an Xray synchrotron radiation scanner and simulated the triaxial test of sand by DEM, which verified the effectiveness of CT technology in the study of SRM. Yan et al. [11] used sphere unit and sphere unit combination to simulate soil particles and rock particles, respectively, and studied the influence of spatial distribution of soilrock particles on SRM macromechanics. Graziani et al. [12] determined the relationship between the micro and macromechanical properties via DEM.
It is can be summarized from the above literature reviews that the shortcomings remain in the numerical simulation of SRM. (1) The existing numerical simulation methods were mostly aimed to study the mesostructure of SRM and its impact on the macroperformance, and almost no related virtual experiments have been established for the traditional laboratory tests. (2) There are few studies on the SRM with largesize particles in engineering applications. (3) Most of the existing molding methods of SRM in DEM are still static pressure molding, and there is no research on the SRM model based on the VVTM. To addresses these gaps, the authors have performed the following work: the numerical simulation method of the VVTM for SRM (NSMSRM) based on the DEM and CT scanning was developed; based on NSMSRM, the CBR of SRMG was studied comprehensively, and the CBR predictive model of SRMG was established; the reliability of the model was verified by practical engineering.
2. Testing Methods and Materials
2.1. Test Methods
2.1.1. Vertical Vibration Test Method
To fully simulate the rolling effect of a modern heavyduty vibratory roller, Chinese scholars [4, 13] have developed the indoor vertical vibration test method (VVTM) based on the principle of the directional vibratory roller, as shown in Figure 1. The VVTM method includes the vibration compaction method (VCM) and the vibration shaping method (VSM). The VCM was used to determine the maximum dry density and optimum water content, and the VSM was used to shape the mechanical specimens. By studying the effect of the vibration frequency, exciting force, and static pressure on the dry density of SRM, the vibration compaction parameters of SRM are recommended as shown in Table 1 [14].

2.1.2. Test of the CBR
According to the T01341993 bearing ratio test in “Test Methods of Soils for Highway Engineering” (JTG E402007) [15], the vibration forming test specimen was placed on the lifting platform of the pavement material strength tester, and the load with 45 N was applied on the penetration rod; then the dial indicator pointer of force and deformation measurement was adjusted to integer; then the initial reading was recorded. The speed of pressing in the specimen was 1 mm/min. The forces of 0.3 mm, 0.6 mm, and 1.0 mm pressure sensors on the data display instrument were recorded, and then the pz relation curve with the unit pressure (p) as the abscissa and the penetration (z) as the ordinate was drawn. The ratio of unit pressure and standard unit pressure was obtained when the penetration is 2.5 mm as the bearing ratio of the materials, which is recorded as the CBR.
2.2. Materials
The SRM is usually divided into soil and rock. Particles smaller than 5 mm in diameter are considered soil, and those larger than 5 mm are considered rock [16]. The soilrock ratio is the mass ratio of soil to rock. The corresponding soil content P_{5} is the mass ratio of soil to SRM by mass. The used SRM was taken from the construction site of subgrade excavation in Wuyi County, Jinhua City, Zhejiang Province. The characteristics of raw materials are shown in Tables 2 and 3.


3. Establishment of the Numerical Simulation for the CBR
The indoor CBR test method is only suitable for the mechanical property test of SRM_{d≤40 mm} and cannot analyze the mechanical behavior of SRMG. Therefore, this paper presents a numerical simulation method of the CBR (NSMCBR) for SRMG.
3.1. Methodology to Establish the NSMCBR for SRMG
The method to establish NSMCBR for SRMG is shown in Figure 2.
3.2. Procedures to Establish the NSMCBR for SRMG
3.2.1. CT Scanning of the Shape Characteristic of Rock Particles
In this paper, CT scanning technology was used to randomly collect 50 representative shape features of rock particles, obtain the digital information of the shape of the particles, and reconstruct the 3D contour of the rock particles. The experimental equipment is an industrial CT scanning device, the image reconstruction matrix is 512 × 512 pixels, the maximum spatial resolution is 0.02 mm, the density resolution is 0.3%, the image display matrix is 1024 × 1024 pixels, and the time of image reconstruction is less than 3 s, as shown in Figure 3. First, according to the coordinate matrix of aggregate surface contour points, the cloud point data of particle contour (Figure 4) were reconstructed by using the MATLAB program and converted into an STL file. Then, the STL file was used to generate the aggregate surface geometry model by using the Geometry import command in PFC3D, and then some spherical particles were used to fill the geometric model of aggregate to create the clump template via the command of Clump Template in PFC3D, so as to generate aggregate models of different shapes and sizes, as shown in Figure 5.
3.2.2. Constitutive Model of SRM
Considering the engineering properties of the materials, a linear contact model with damping was used to simulate the properties of SRM [17, 18]. The normal and tangential stiffness F_{n} of two contact entities in the linear contact model were realized. The contact between two entities occurred at a very small point, and only the transfer force cannot resist the bending moment, so the contact moment M_{c} was equal to zero. Contact force F_{c} was divided into elastic force F_{l} and inelastic force F_{d}. The elastic and inelastic properties of the contact elements were usually represented by the spring constant K and damper β, and the friction between particles was represented by a slider with a friction coefficient. The mechanical relationship between particle elements is shown in Figure 6.
In Figure 6, K_{n} and K_{s} represent the normal and tangential elastic constants between two particle elements, β_{n} and β_{s} represent the normal and tangential damping coefficients, respectively, and u represents the friction coefficient of the slider. The linear spring is parallel to buffer , and there is a gap on the surface, which is called the distance between contacts. Theoretically, the particle surface contact is shown in Figure 7. When the distance between them is less than or equal to zero, the contact surface is active.
3.2.3. Numerical Simulation for the CBR Test
According to the bearing ratio test in “Test Methods of Soils for Highway Engineering” (JTG E402007) [15], the CBR particle discrete element model was established. Step 1: the model of the test tube, bearing plate, and base was generated: The modulus of the test tube, bearing plate, and base was 30 GPa, and Poisson’s ratio was 0.2. The inner wall diameter of the test tube was 150 mm, the outer wall was 160 mm, and the height was 150 mm; the diameter of the bearing plate was 149 mm and the height was 25 mm; the diameter of the base was 180 mm and the height was 25 mm. The DEM model should also have the same size as the physical experiment model so that a meaningful comparison can be made between them, as shown in Figure 8. Step 2: formation of soilrock particles: According to the particle composition, geometry, and physical parameters, the corresponding model was generated. The soil particle model was simplified as a sphere, and the rock particle model was input according to the CT scanning results. In order to avoid generating too many particles leading to inefficient simulation calculation, the minimum particle size of the model was controlled as 1 mm [19]. The physical parameters of the SRM are set in the PFC^{3D} simulation as shown in Table 4. Step 3: formation of SRM specimens: First, the gravity was applied so that all particles can reach a static equilibrium state, as shown in Figure 9. The vibration frequency was 25 Hz, the force peak was 5.3 kN, and the vibration compaction time was selected according to the compaction degree K, as shown in Table 5. Figure 10 shows the SRM virtual specimen after vibration compaction. Step 4: determination of the CBR: In the loading process, the load plate with a total mass of 10 kg was set first. The load of 45 N was applied on the penetration rod, and the load was started after the test piece is stable, and the force on the penetration rod was recorded. The uniaxial compression loading method was adopted in the numerical simulation. The bearing plate was loaded axially at the loading rate of 1 mm/min for 5 minutes until the penetration reached 5.0 mm; the historical records of the unit pressure (p) on the penetration rod were recorded, and the loading time was output every 30 seconds (i.e., the penetration is 0.5 mm). In addition, the penetration unit pressure p at 2.5 mm penetration was obtained, as shown in Figure 11. The CBR value of the test piece was calculated by the formula: CBR = p_{2.5}/7 × 100%.

(a)
(b)

(a)
(b)
(a)
(b)
3.3. Reliability Verification
The numerical and laboratory tests on the CBR of SRM were performed for nine kinds of soilrock ratios: 100 : 0, 80 : 20, 70 : 30, 60 : 40, 50 : 50, 40 : 60, 30 : 70, 20 : 80, and 0 : 100. The gradation of the test is shown in Table 6. The value of laboratory CBRp and simulation CBRp,s was obtained; then the cumulative relative deviation δ was calculated as the following equation:

During the process of numerical testing, the friction coefficient between soil particles and rock particles was repeatedly adjusted. The results show that δ = 0.35 reached the minimum when the friction coefficient of the rock particle is 0.41 and the damping coefficient is 0.65 and the friction coefficient of the soil particle is 0.31 and the damping coefficient is 0.34. The forcedeformation curves of SRM with different soilrock ratios obtained by NSMCBR are shown in Figure 12, and the curve between the numerical and experimental CBR of the SRM is shown in Figure 13. The graph shows that the numerical values correspond well with the experimental values. The maximum deviation is 7.4% and the average deviation is 1.1%, which show that the NSMCBR has high reliability. This result verifies the feasibility and accuracy of the numerical simulation and lays a foundation for the next NSMCBR for SRMG.
4. Application of NSMCBR to Investigate CBR of SRMG
4.1. Research Programme
NSMCBR was carried out on SRMG with three kinds of maximum particle size D_{max} (60 mm, 80 mm, and 100 mm), four kinds of largesize particle content P (0%, 10%, 20%, and 30%), and five kinds of soilrock ratios (30 : 70, 40 : 60, 50 : 50, 60 : 40, and 70 : 30). P refers to the ratio of rock particles within 40–100 mm to SRM by mass, and the soilrock ratio refers to the ratio of soil content to rock content by mass. The influence rules of maximum particle size, soilrock ratio, and largesize particle content on the CBR of SRMG were investigated.
4.2. Results and Discussion
4.2.1. Effect of the SoilRock Ratio on the CBR of SRMG
Figure 14 shows the variation rule of the CBR of SRMG obtained by NSMCBR with the soilrock ratio. The results show that the CBR of SRMG increases with the decrease in the soilrock ratio in the range of 70 : 30–30 : 70. This phenomenon can be explained from two aspects: on the one hand, with the gradual increase in rock particles, rock particles are embedded and contacted with each other to form a stable skeleton structure. The stress will be transmitted by the force chain through the contact point of rock particles, which can resist larger deformation [20]. On the other hand, the overall deformation resistance of the mixture depends on the stiffness of its constituent materials, so the CBR value of SRMG will increase with the increase in rock particles with larger stiffness and smaller deformation.
(a)
(b)
(c)
4.2.2. Effect of the LargeSize Particle Content on the CBR of SRMG
Figure 15 shows the variation of the CBR of SRMG obtained by NSMRM with largesize particle content (P). The results show that the CBR of SRMG increases with the increase in P. The reason is that in the simulation scheme, the size of the test cylinder used is certain. In a certain range of model volume, with the increase in the content of largesize particles, the largesize particles are more likely to be close to each other to form a skeletonlike dense structure. Therefore, the contact force between particles and the CBR value can be increased.
(a)
(b)
(c)
4.2.3. Effect of the Maximum Particle Size on the CBR of SRMG
Figure 16 shows that the CBR value of SRMG increases with the increase in maximum particle size D_{max.} The reason is that the SRMG model is uniformly generated in a certain volume range, so the gradation of largesize particles is identical. With the increase in the maximum particle size, the total number of particles decreased, and the total surface area and incomplete contact among the particle decreased. As a result, the porosity decreases, the model becomes denser, and the CBR increases. It can be seen that the skeleton structure formed by largesize particles can improve the CBR value of the soilrock mixture. Therefore, while using the soilrock mixture to fill the subgrade, a largesize rock can be considered to enhance the permanent deformation resistance of the subgrade.
(a)
(b)
(c)
5. Prediction Model of the CBR and Its Application
5.1. Establishment of CBR Prediction Model
Given that the CBR of SRMG has as a positive correlation and a nonlinear relationship with the maximum particle size, content of largesize particles, and CBR, a predictive model was simplified as shown in the following equation:where A, B, C, and D are constants; CBR_{40} is the indoormeasured CBR for SRM_{d≤40 mm}, %; P is the content of the largesize particle, %; D_{max} is the maximum particle size of SRMG, mm; P_{5} is the soil content of SRM, %; CBR is the predicted CBR of SRMG, %.
According to the simulated results, the model parameters were calculated as follows: A = 0.18346, B = 0.17317, C = 0.49944, and D = 0.04678. The simplified prediction model of the CBR of SRMG was obtained as follows:
To verify the validity of the prediction model for the CBR, the predictive results and NSMRM results are compared in Figure 17. The ratio of the predicted and simulated CBR basically fluctuated around the value of 1, which showed that the prediction model has high accuracy and good predictive effect.
5.2. Determination Method of the Predictive CBR Value
The CBR of SRMG is estimated by taking materials from the construction test section. The steps are as follows:(i)SRMG samples are divided into A1 (d ≤ 40 mm) and A2 (40 mm < d < 100 mm) using a 40 mm sieve; then the contents of largesize particle P and soil content P5 were calculated.(ii)The maximum size (D_{max}) in A2 was determined.(iii)The VVTM was carried out on part A1, and the maximum dry density ρ_{dmax} and optimum water content ω_{opt} of part A1 were determined. Then, the indoor cylindrical specimens of 150 mm × 100 mm were prepared via the VVTM according to the compactness K and determined ρ_{dmax} and ω_{opt}.(iv)The CBR value of SRMG was calculated by substituting the parameters CBR_{40}, P, P_{5}, and D_{max} into equation (3).
5.3. Field Application and Verification
The Wuyi section of WuchengWuyi Highway, 235 National Road, Jinhua City, Zhejiang Province, has a total length of 4.472 km and a design standard of a twoway sixlane firstclass highway. The design speed is 80 km/h. The subgrade was filled with SRM. The maximum particle size of used SRM in the test was smaller than 100 mm. The soil particles were quartzite with 42.6 MPa of the average saturated uniaxial compressive strength. There were two sections in the test. The maximum particle sizes of Section I and Section IIsec2 were 60 mm and 80 mm, respectively.
The field CBR test was carried out on two sections, respectively, according to the field CBR value test method of “Field Test Methods of Subgrade and Pavement for Highway Engineering” (JTG E602008) [21]. The values of the field and prediction are shown in Table 7.

The results show that the predicted value basically coincided with the field value in the field with a maximum deviation of 1.5%. Thus, the established prediction model was reliable and can provide more precise design parameters of the CBR for the design subgrade filled with SRM.
6. Conclusions
(1)The numerical simulation method of the CBR of soilrock mixtures was established based on CT scanning and DEM. The maximum deviation in the CBR between numerical simulation and the laboratory test was 7.4%, which showed that the established NSMCBR has high reliability.(2)Based on the NSMCBR, the effects of the maximum particle size, soilrock ratio, and largesize particle content on the CBR of SRMG were investigated. The results showed that the CBR of SRMG linearly increased with the decrease in soil content and the increase in maximum particle size and content of the largesize particle.(3)The prediction model of the CBR of SRMG was established. The comparison of the predicted value and the field value showed that the maximum deviation between the two values was 1.5%. It showed that the prediction model established in this study has high reliability and can be widely applied.
Data Availability
The data used to support the findings of this study are obtained directly from the tests and included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China under project no. 51908460 and the Scientific Research of Central Colleges of China for Chang'an University under project nos. 300102219520 and 310821172008.
References
 W.J Xu and H.Y Zhang, “Research status and development trend of soilrock mixture,” Advances in Science and Technology of Water Resources, vol. 33, no. 01, pp. 80–88, 2013. View at: Google Scholar
 X. P. Ji, “Study on vibration test method for soil of earthrock mixed fill subgrade,” Highway, vol. 64, no. 02, pp. 27–32, 2019. View at: Google Scholar
 X. Ji, Y. Jiang, and Y. Liu, “Evaluation of the mechanical behaviors of cementstabilized cold recycled mixtures produced by vertical vibration compaction method,” Materials and Structures, vol. 49, no. 6, pp. 2257–2270, 2016. View at: Publisher Site  Google Scholar
 Y. Jiang and J. Xue, “Investigation into physical and mechanical properties of SRXstabilised crushed rock using different compaction methods,” International Journal of Pavement Engineering, vol. 20, no. 7, pp. 866–873, 2019. View at: Publisher Site  Google Scholar
 D. Bao and X. Zhang, “Granular matter and granular flow,” Journal of Zhejiang University (Science Edition), vol. 21, no. 05, pp. 514–517, 2003. View at: Google Scholar
 L. Jin, “Study of effect of mesomechanical parameters on macromechanical character of soilrock mixture,” Engineering Journal of Wuhan University, vol. 51, no. 05, pp. 409–417, 2018. View at: Google Scholar
 Y. X. Zhao, Z. X. Liu, and C. Mazzotti, “Numerical experiments on triaxial compression strength of soilrock mixture,” Advances in Civil Engineering, vol. 2019, Article ID 8763569, 15 pages, 2019. View at: Publisher Site  Google Scholar
 Z. Cheng and J. Wang, “Experimental investigation of interparticle contact evolution of sheared granular materials using Xray microtomography,” Soils And Foundations, vol. 58, no. 6, pp. 1492–1510, 2018. View at: Publisher Site  Google Scholar
 Z. Cheng and J. Wang, “A particletracking method for experimental investigation of kinematics of sand particles under triaxial compression,” Powder Technology, vol. 328, pp. 436–451, 2018. View at: Publisher Site  Google Scholar
 M. Wu, J. Wang, and A. Russell, “DEM modelling of minitriaxial test based on onetoone mapping of sand particles,” Géotechnique, vol. 1–14, 2020. View at: Publisher Site  Google Scholar
 Y. Yan, J.F Zhao, and S.Y Ji, “Discrete element analysis of the influence of rock content and rock spatial distribution on shear strength of rocksoil mixtures,” Engineering Mechanics, vol. 34, no. 06, pp. 146–156, 2017. View at: Google Scholar
 A. Graziani, C. Rossini, and T. Rotonda, “Characterization and DEM modeling of shear zones at a large dam foundation,” International Journal Of Geomechanics, vol. 12, no. 6, pp. 648–664, 2012. View at: Publisher Site  Google Scholar
 M. J. Li, “Study on vibration compaction characteristics of road materials,” Chang’an University, Xi’an, China, 2002, Ph.D. thesis. View at: Google Scholar
 X.P Ji, S.W Li, and Y. Xiong, “Research on vibration test method for soilrock mixture,” Highway, vol. 64, no. 02, pp. 27–32, 2019. View at: Google Scholar
 Transport, “Research institute of highway ministry of test methods of soil for highway engineeringJTG E402007,” 2007. View at: Google Scholar
 Y. Dong, “Experimental study on intensity character of rocksoil aggregate mixture,” Rock and Soil Mechanics, vol. 06, pp. 1269–1274, 2007. View at: Google Scholar
 Y. L. Tu, “The macro and meso mechanical characteristics and nonlinear elastoplastic constitutive model of soilrock aggregate,” Chongqing University, Chongqing, China, 2017, Ph.D. thesis. View at: Google Scholar
 P. Wang, “Soilrockmixture largescale triaxial test and mesoscopicmechanics characteristics of pfc simulation,” Chongqing University, Chongqing, China, 2017, Master thesis. View at: Google Scholar
 Yu Li, L. Li, and K. Chen, “Study on the relationship between particle size distribution and stressstrain of soil,” Journal of Railway Science and Engineering, vol. 15, no. 07, pp. 1722–1729, 2018. View at: Google Scholar
 L. Fan, “Study on the rapid detection method of compaction quality with earthrock subgrade based on PFWD,” Chang’an University, Xi’an, China, 2011, Master thesis. View at: Google Scholar
 Transport, “Research institute of highway ministry of, filed test method of subgrade and pavement for high engineeringJTG E602008,” 2008. View at: Google Scholar
Copyright
Copyright © 2020 Xiaoping Ji 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.