Soil-rock mixtures with large particle size variations are often used as fill materials for expressway construction in mountainous areas. Conventional testing methods do not enable fast and nondestructive monitoring of real-time changes in the compaction quality of soil-rock filled subgrades. Selecting an appropriate evaluation method is the key to controlling the compaction quality of a soil-rock filled subgrade. In this study, three-dimensional DEM models of subgrade materials were reconstructed by a spherical harmonic series whose harmonization degree was fixed at 15. The macroscopic and mesoscopic behaviours and characteristics of the subgrade under vibratory rolling were analysed. The results showed that the porosity, contact force and coordination number of the subgrades tended to be stable in the last two passes. The subgrades with 4 filler combinations presented the similar mechanical anisotropy and meso-mechanical states. On-site monitoring of subgrades under vibratory rolling and settlement after construction was performed, and the results were considered. An evaluation method and criterion to control the compaction quality of the SRM subgrade was proposed, i.e., whether the average value of the vibration compaction value from the second to last pass differed by more than 2% from the average value in the last pass.

Spatial distribution orientations of blocks can cause significant errors in the discrete element model (DEM) calculation of soil-rock mixture (SRM). To avoid this error, spherical harmonic (SH) series whose harmonization degrees fixed at 15 were proposed for block reconstruction. This research refers to the case-history of a deep excavation rift valley spanning from the Mabian to Zhaojue of the Leshan-Xichang Expressway, mainly containing moderately-weathered silty mudstone, in the Leshan City, Sichuan Province, China. The appropriate degree of finite-term SH series is selected by the volume, surface area. 100 blocks were scanned on site, and sphericity and angularity of the blocks were calculated. The sphericity and angularity of 50 reconstructed blocks were considered for the error analysis of SH method. Moreover, stochastic polyhedron method was considered for comparing different block reconstructions. The maximum block placement angle was defined to control the spatial distribution orientations of the blocks. Large scale direct tests were carried. Numerical simulations of large-scale direct shear tests were conducted to study the influence of the spatial distribution orientations of the blocks on the mechanical properties of the SRMs. The results revealed that the finite-term SH series fixed at 15 accurately reflected the shape characteristics and mechanical behaviors of actual blocks. The spatial distribution orientations of the blocks had a minimal impact on the friction angle and cohesion of SRM constructed through the SH method. The SRMs developed via the SH method exhibited marginal variations in contact force and anisotropy index of contact across diverse block placement strategies. The evolution of coordination number was closer when employing the SH method under varied block placement methods. Blocks reconstructed by the SH method, could mitigate errors in DEM calculation caused by the spatial distribution orientations of the blocks.