Abstract

In this work, to quantitatively analyze the roughness of the surfaces of road aggregates, the contact measurement technique and contactless scanning technique were, respectively, used to capture the coordinate data of point clouds on the aggregate surface, which were then used to reconstruct the digital elevation models of aggregate particles. Then, the joint roughness coefficient (JRC) was used as an evaluation index, and the quantitative calculation methods of the two-dimensional (2D) contour line roughness and three-dimensional (3D) contour surface roughness of aggregate particles were, respectively, studied. Finally, the anisotropic characteristics and size effect of the roughness coefficients of aggregates with different lithologies were, respectively, investigated, based on which the practicability of the 3D roughness coefficient index was proven. The results demonstrate that the roughness of a road aggregate surface can be quantitatively described by the point cloud data. The 2D roughness of aggregate profile lines exhibits anisotropy, while the 3D roughness of the aggregate contour surface indicates the size effect. The subtle morphological changes of the surface textures of aggregates can be accurately described by the 3D joint roughness coefficient (JRC3D) calculated by the feature parameter method.

1. Introduction

Rock aggregates are commonly used raw materials in the field of road engineering, and their surface roughness significantly affects the performance of asphalt mixtures [1]. The rough surface is beneficial to increasing the contact area between the aggregate and asphalt, thereby creating a certain penetration depth of the asphalt, increasing the thickness of the asphalt film, and enhancing the physical adsorption of the two materials on the contact surface [26]. Additionally, there are many fine protrusions on the surfaces of aggregates, which is conducive to the formation of the aggregate skeleton structure of the asphalt mixture during the mixing process, thereby endowing asphalt concrete with good shear strength and ultimately improving the rutting resistance of the road surface [79].

From the aspect of the collection of aggregate surface morphology data, with the continuous development of measurement technology, a variety of measurement methods have been applied to obtain the characteristic information of the surfaces of rock materials. According to the acquisition method, these methods can be divided into contact contour measurement methods and noncontact optical measurement methods. Additionally, according to the dimension, the data can be divided into two-dimensional (2D) contour line coordinate data and three-dimensional (3D) contour surface coordinate data [10, 11]. To quantitatively evaluate the roughness coefficient of road aggregates, the accuracy requirements must first be met. Moreover, the measurement instruments, test methods, and calculation processes should not be too complicated and should be characterized by a fast measurement speed and convenient statistical analysis. Therefore, the contact measurement method was commonly used to quickly calculate the 2D roughness of the profile contour of the aggregate surface and the 3D topography data can also be obtained via a large amount of 2D contour data interpolation [12, 13]. In contrast, the high-precision noncontact measurement technology based on the optical principle has several advantages for obtaining data on the surface topography of rock materials [14]. A 3D scanning device can quickly obtain point cloud data describing the 3D features of the aggregate particles. These data can be used for the physical reconstruction of the aggregate and the quantification of 3D morphological parameters [15].

The quantitative evaluation methods for the characterization of aggregate surface roughness primarily include the fractal mathematics method and statistical parameter method [16]. In the fractal mathematics method, more attention is paid to the correlation between the characteristic fractal dimension of the surface topography of the rock structure and the roughness and mechanical properties [17]. In contrast, in the statistical parameter method, quantitative evaluation is conducted by establishing the regression relationship between the characteristic parameters of the surface topography of rock materials and the roughness index [15]. The joint roughness coefficient (JRC) recommended by the International Society of Rock Mechanics was commonly used as the evaluation index of the surface texture roughness of the rock aggregate [18]. In geotechnical engineering research, roughness is often used to analyze the correlations between the rock surface profile and the interface mechanics and seepage characteristics [19]. In road engineering, the roughness coefficient can also be used as an important index to analyze the surface texture roughness of aggregate particles and the correlation between asphalt and aggregate adhesion.

Therefore, in this study, the collection methods for aggregate surface morphology data and the applicability of the joint roughness coefficient (JRC) to the evaluation of the changes of the mesostructure of the aggregate surface texture are investigated. First, the measurement and quantitative evaluation method of the 2D roughness of the aggregate surface based on contact measurement is studied, and the anisotropic characteristics of aggregate surface texture roughness are analyzed. Second, the measurement and quantitative evaluation method of the 3D roughness of the aggregate surface based on 3D scanning technology is studied, and the size effect and resource characteristics of the aggregate surface texture roughness are analyzed. In addition, a polishing test of the aggregate was conducted to explore the influences of the resource characteristics and polishing degree on the surface roughness of aggregates.

2. Materials and Methods

2.1. Materials

Four types of coarse aggregates with different lithologies were obtained from the Highway Science and Technology Research Institute of Yunnan Province. The aggregate particle size range was 16–19 mm, and the aggregates were produced using jaw and counterattack crushers. The technical indices of the raw materials were obtained according to the Chinese specification “Test Methods of Aggregate for Highway Engineering” (JTG E42-2005). The technical indices of the coarse aggregates are listed in Table 1.

2.2. Evaluation Index of Surface Roughness

To quantitatively analyze the roughness of the surfaces of road aggregates, the joint roughness coefficient (JRC) was used as an evaluation index in this study. Combining engineering experience and experimental research results, Barton (1977) presented 10 typical rock surface contour lines with values ranging from 0 to 20 to evaluate the JRC value [20], as shown in Table 2.

2.3. Aggregate Surface Topography Feature Data Collection
2.3.1. Collection Method of 2D Contour Line

A Kasaka FGA51 high-precision stylus profiler, which operates on the principle of the needle-tracing method [21], was used to obtain the 2D profile curves of the aggregate surface (see Figure 1(a)). The x-axis measurement accuracy of the profiler is 1.2 μm, the grating resolution is 0.1 μm, the Z1 measurement accuracy is 1.0 μm, the resolution is 0.05 μm, and the measurement speed range is 0.02–4 mm/s; these characteristics meet the requirements of aggregate surface topography measurement (see Figure 1(b)). During the measurement process, the aggregate was fixed on a precision linear motion guideway in the x-direction so that the stylus was in contact with the surface of the specimen. Then, the stylus movement speed was set to 0.5 mm/s, the fluctuation range was set to ±0.02 mm/s, the sampling step was set to 2 μm, and the sampling line distance was set to 10 mm. By setting the parameters of the measurement system, the measurement command was issued, and the computer accepted the command to drive the control instruction to make the stylus move to the moving guideway. When the aggregate moved at the sampling interval along the x-axis, the x-coordinate of the measurement point was obtained at the position to be measured, and the stylus moved up and down with the fluctuation of the surface profile, causing the z1-axis digital sensor to output electrical signals. The signals were transmitted to the computer for processing via the photoelectric conversion processing circuit, and the z-coordinate could then be determined. This process was repeated for the observation section of the aggregate surface presented in Figure 1(c) to obtain the corresponding 2D profile curve. Finally, the data from the effective test length of 10 mm were filtered and digitized and were drawn in proportion to the coordinate grid. The roughness coefficients JRC2D of 9 contours of the aggregate surface were estimated from the visual comparison with the 10 typical rock surface contour lines shown in Table 2, as presented in Figure 1(d).

2.3.2. Collection Method of 3D Contour Surface

An EinScan-S 3D scanner was used to obtain the point cloud data that contain information about the surface geometric characteristics of the aggregate particles. The resolution of the profiler lens is 3264×2488, the scanning size is 200 × 200 mm, the measurement accuracy is 20 μm, there are 8 million effective 3D points, and omnidirectional scanning is supported; these characteristics meet the requirements of aggregate surface topography measurement (see Figure 2(a)).

The 3D scanner system operates on the principle of white-light interferometry [22, 23]. During the measurement process, the aggregate sample was placed on the stage, the white-light emitted by the light source was refracted by the lens and the beam splitter and was irradiated to the rough textured surface of the aggregate. Then, the white-light intensity of the light source was adjusted according to the light intensity and color temperature of the room, and the carrier table was simultaneously adjusted for calibration until the interference image of the aggregate surface captured by the charge-coupled device (CCD) camera in the visual interface of the software met the analysis accuracy requirements. After setting the scanning range, sampling step length, and number of frames, multiframe images of interference fringes on the 3D contour surface of the aggregate were obtained via the CCD, and the 3D scanning point of the aggregate surface in the coordinate system of the scanner was calculated by the analysis software. The principle of 3D scanning technology is presented in Figure 2(b).

After scanning the 3D contour surface of the aggregate test sample from multiple angles (see Figure 2(c)), the iterative closest point (ICP) algorithm in the Imageware software was used to register the point cloud of multiple scans. The purpose of the ICP algorithm was to determine the rotation and shift parameters of the point cloud scans from different angles to register the dataset and obtain an optimal match for the data with different incident angles due to changes in the light source [24].

The registered dataset was imported into the PCL library to eliminate outliers and noise. This library is a large cross-platform open-source C++ programming library containing many algorithms and data structures related to point clouds, including point cloud acquisition, filtering, segmentation, registration, retrieval, feature extraction, recognition, tracking, surface reconstruction, and visualization [25]. The VoxelGrid filter was used for noise reduction. A 3D mesh of the aggregate particle was generated after eliminating 96% of the point cloud data (see Figure 2(d)).

2.4. Quantitative Characterization of the Surface Roughness
2.4.1. Calculation Method of 2D Contour Line Roughness

(1) Characteristic Parameter Method. To quantitatively analyze the roughness of the two-dimensional profiles of the rock materials, a least-squares midline running through the sampling length was used as the baseline [26, 27], as shown in Figure 3.

The regression equation of the relationship between the geometric characteristic parameters and roughness proposed by Tse and Crude is generally used to calculate the 2D roughness (JRC2D) of the surface profile of rock materials [28] and is defined as follows:

The parameter Z2 is a characteristic parameter of the slope of the rock surface that is not affected by the position of the calculated reference line [29] and is expressed as follows:

When using 3D point cloud data to calculate the roughness of the 2D contour of an aggregate surface, equation (2) can be approximately expressed as follows:where L is the total length of the 2D profile contour line, n is the number of calculated samples, and xi and zi are the point cloud coordinates on the contour line.

(2) Trace Length Characterization Method. To quantitatively evaluate the roughness coefficient of road aggregates in practical engineering applications, the accuracy requirements must first be met. Moreover, the measurement instruments, test methods, and calculation processes should not be too complicated and should be characterized by a fast measurement speed and convenient statistical analysis. Therefore, the trace length characterization method proposed by Soudan was used to quickly calculate the 2D roughness of the profile contour of the aggregate surface [30] and is defined as follows:where Rp is the ratio of the actual length of the profile line to its straight-line length and xi, xi+1, yi, and yi+1 are the relative coordinates of point i and point i+1, as indicated by the subscripts.

Table 3 presents the functional relationship between Rp and JRC2D under different sampling intervals [29].

(3) Evaluation Accuracy of 2D Contour Line Roughness. According to the relative coordinate data of the measured curve exhibited in Figure 1(d), the roughness coefficients JRC2D of 9 contours of the aggregate surface were first calculated based on the characteristic parameter method. The measurement point data of each contour curve were then sampled at 0.5 mm intervals and substituted into equation (5) to determine the elongation rate Rp of the contour line. The roughness coefficient JRC2D of each 2D profile contour line was calculated by the functional relationship between Rp and JRC2D at a sampling interval of 0.50 mm. The calculation results of the two analysis methods are reported in Table 4.

From the calculation results, it can be seen that the characteristic parameter method and the trace length method can both quantitatively characterize the roughness coefficient of the 2D profile of the aggregate texture, and the calculation deviation was within 13%. Similarly, other test samples were also governed by the same law, and a good correlation was found. The comparative evaluation results demonstrate that the calculation deviation of the two roughness characterization methods was mainly caused by the difference in the sampling interval. Compared with that of the characteristic parameter method (0.002 mm), the sampling interval of the trace length method was 0.5 mm. Under the relatively large sampling interval, a portion of the mesomorphological changes of the measurement curve could not be collected, resulting in the overall low roughness value.

2.4.2. Calculation Method of 3D Contour Surface Roughness

The 2D roughness coefficient of the aggregate cannot fully reflect the geometric shape of the aggregate surface, so it is necessary to use the 3D roughness coefficient to characterize the geometric characteristics of the rough texture of the aggregate profile. In the calculation of 3D roughness, the weighted average method and characteristic parameter method are typically used.(1)Weighted Average Method. The roughness JRC2D of a large number of 2D contours on the surface of a rock material can be measured, and the 3D average roughness of this area can be calculated by the weighted average method, as follows:where is the 3D average roughness, JRCi2D is the roughness of each 2D undulating curve, and m is the number of section lines.(2)Characteristic Parameter Method. It is assumed that the coordinates of the point cloud on the rough surface of a rock are continuously differentiable. The root mean square Z2s of the 3D characteristic parameter of the relative height fluctuation can be calculated as follows [31]:When using the 3D point cloud data to calculate the 3D profile roughness coefficient JRC3D of an aggregate surface, the approximate formula is as follows:where Nx is the number of point clouds on the x-axis in the length direction of the rock profile, Nx is the number of point clouds on the y-axis in the height direction, Δx is the sampling spacing of point clouds on the x-axis, Δy is the sampling spacing of point clouds on the y-axis, and Zi and Zi+1 are the coordinates of the i-th point and the (i+1)-th point in the z-axis direction of the roughness height, respectively. The relevant parameters are presented in Figure 4.(3)Evaluation Accuracy of 3D Contour Surface Roughness. According to the roughness coefficient JRC2D of 9 contours of the aggregate surface, the 3D average roughness of the analysis area of the aggregate surface was calculated to be 14.6 via equation (6). The texture feature of the analysis area was then evaluated by the 3D contour surface roughness JRC3D. After scanning the 3D contour surface of the aggregate test sample from multiple angles, the data point cloud of the analysis area (Figure 1(c)) of the aggregate fracture surface was sampled at intervals of 10 × 10 mm, that is, Δx = Δy = 0.25 mm. Equation (8) was used to calculate the characteristic parameter Z2s, and the 3D contour surface roughness JRC3D of the analysis area of the aggregate surface was then calculated to be 13.8 via equation (1). The differences and correlations between the two calculation methods for the calculation of the 3D roughness are discussed in the subsequent section.

2.5. Change Method for the Surface Roughness of the Aggregate

An aggregate polishing test was designed to test the accuracy of the roughness coefficient JRC for the analysis of changes of the aggregate surface roughness. In the test, 1200# carborundum with a 95% alumina content was used as an abrasive, and a vibrating polishing machine was employed to change the mesostructures of the surface textures of the four aggregates with different lithologies. The steps of the aggregate polishing test were as follows:(1)After fully soaking the abrasive in water, it was drained with a basket and poured into the vibrating polishing machine.(2)The power was turned on and the wet water-drained aggregate was arranged in the abrasive in the vibrating polishing machine to be polished from different directions one-by-one.(3)After the carborundum was uniformly added to the mixture, the polishing time was set according to the polishing requirements; after the polishing test was completed, the carborundum in the mixture was washed with clean water.

At the end of each polishing test, the profile point cloud data of the particle surface were collected again, and the changes of the aggregate surface roughness under different polishing conditions were compared and analyzed.

3. Results and Discussion

3.1. Anisotropy Analysis of Aggregate Surface Roughness

To analyze the differences in the roughness of the 2D contour line of the aggregate surface in different directions and angles, coarse aggregates with a particle size between 16 and 19 mm were used as the research object, and the size of the aggregate fracture surface was 10 × 10 mm. The range was used as the analysis area, as shown in Figure 5(a). Twelve observation cross sections were, respectively, set up along the x- and y-axes at equal intervals of 0.8 mm, as shown in Figures 5(b) and 5(c). The origin of the coordinates and the direction of 0° in the analysis area was defined with the direction of 0° as the reference line, and an observation section was set up every 15° through this point for a total of 12 observation sections, as shown in Figure 5(d).

The characteristic point coordinate data of the observation section contour lines were measured from the contact stylus profiler, and the 2D contour curve was drawn. According to the functional relationship in the characteristic parameter method, the JRC2D value of each curve was programmed in MATLAB, as reported in Tables 57.

The test results demonstrate that, under different observation section sampling rules, the roughness of the 2D contour lines of the aggregate particle surface was different, and the apparent morphologies were quite different and had no correlation with the sampling direction or angle, thereby exhibiting anisotropy. When the JRC2D index is used to characterize the roughness characteristics of the texture of aggregates, the randomness of the selection of the observation cross sections will lead to deviations in the evaluation of roughness. Therefore, the 2D roughness index cannot describe the overall morphology of the aggregate surface, and it is necessary to use the 3D roughness index to comprehensively characterize the spatial geometric characteristics of the aggregate particle profile.

3.2. Influence of the Size Effect on Aggregate Surface Roughness

To analyze the roughness of the different observation areas with the same centroid and different dimensions of the aggregate fracture surface, the weighted average method and the characteristic parameter method were used to calculate the 3D average roughness of the analysis area of the aggregate fracture surface.

Via the weighted average method, equation (5) was used to calculate the 3D average roughness of the analysis area of the aggregate fracture surface. The calculation results of the sampling methods for the three observation sections along the x-axis, y-axis, and circumference were 13.2, 10.6, and 12.2, respectively.

Using the characteristic parameter method, the data point cloud of the analysis area of the aggregate surface was sampled at intervals of 0.1 mm, that is, Δx = Δy = 0.1 mm. In MATLAB, equation (7) was used to calculate the characteristic parameter Z2s, and the 3D contour surface roughness JRC3D of the analysis area of the aggregate surface was then calculated to be 12.6 by equation (1). This calculation result is closer to the average roughness coefficient of the 12 observation sections taken along the circumference calculated by the weighted average method.

The weighted average method (12 observed sections along the circumference) and the characteristic parameter method were used to study the influence of the area difference of the aggregate surface analysis on the roughness calculation results. The 3D roughness of the four analysis areas with the same centroid and different dimensions of the aggregate fracture surface was, respectively, calculated, as presented in Figure 6.

The calculation results reveal that the rough texture of the aggregate surface was characterized not only by anisotropy at the 2D level but also by the size effect at the 3D level. With the increase of the observation range of the aggregate surface, the 2D profile roughness of each observed section and the 3D profile roughness of the analysis area exhibited decreasing trends. The anisotropic characteristics displayed certain similarities in different observation ranges, but no correlation was found. Similarly, other test samples were also governed by the same law. Therefore, when evaluating the 3D roughness of an aggregate surface, the analysis area of the aggregate surface should be maintained at a constant value during the test to ensure the comparability of the evaluation results.

3.3. Influence of Lithology Characteristics on the Aggregate Surface Roughness

The point cloud data of aggregates with four different lithologies were obtained by 3D white-light interferometry, and the weighted average method (observation data were collected every 15° along the circumference) and the characteristic parameter method were used to evaluate the fracture surface of each aggregate within the range of 10 × 10 mm. The roughness results are presented in Table 8.

The test results demonstrate that when the weighted average method and the characteristic parameter method are used to evaluate the 3D roughness in the same analysis area, the difference between the calculation results can be controlled below 10%, and the results, therefore, exhibit a good correlation. Under the experimental conditions of this study, the average 3D roughness value of the limestone aggregate sample was the highest, that of granite was the lowest, and those of basalt and diabase were in the middle. When the weighted average method was used to evaluate the 3D roughness of the aggregate surface, the values of the granite samples were found to range from 8 to 16, and the roughness among the samples varied greatly. In contrast, when the characteristic parameter method was used, the JRC3D values of the granite samples ranged from 9 to 14, the variation of roughness among the samples was relatively small, and the other test samples were also governed by the same law. The comparative evaluation results demonstrate that the calculation deviation of the two roughness characterization methods was mainly caused by the differences in the test methods and data accuracy. During the process of contact measurement, factors such as the excessive hardness of the contact probe, unstable moving speeds, and easy lateral sliding will affect the measurement accuracy. In contrast, 3D noncontact measurement technology based on the optical principle is used to obtain the contour surface feature information of an object and has obvious advantages in the collection of the surface topography data of rock-like materials.

3.4. Influence of Polishing on the Aggregate Surface Roughness

To test the accuracy of the 3D roughness coefficient JRC3D for the analysis of changes of the aggregate surface roughness, an aggregate polishing test was designed. Carborundum with a 95% alumina content was used as an abrasive, and a vibrating polishing machine was employed to change the mesostructures of the surface textures of the four aggregates with different lithologies. At the end of each polishing test, the 3D profile point cloud data of the particle surface were collected again, and the changes of the aggregate surface roughness under different polishing conditions were compared and analyzed. The results are reported in Table 9. The changing trend of the roughness of each lithologic aggregate test sample under different polishing times is shown in Figure 7.

It is evident from the test data that the average roughness of the fracture surface of the limestone aggregate sample was the highest, that of granite was the smallest, and that of basalt and diabase were in the middle. With the increase in polishing time, the rough texture surface of aggregate gradually wore away, which shows that the roughness of the coarse aggregate surface of different lithologies decreases with the increase in polishing time. Among them, the average roughness of the fracture surface of limestone decreases the most, followed by basalt and diabase, while the roughness of granite decreases the least. The results show that the wear resistance of the surface texture of granite aggregate was obviously better than that of limestone aggregate, and the wear resistance of diabase aggregate was similar to that of basalt, which was consistent with the attenuation law of the mechanical index of road aggregate in practical engineering.

Under the experimental conditions of this study, when the 3D roughness of aggregate profile JRC3D ≥ 11.0, there were many pores and cracks on the aggregate surface, and the texture was uneven, which was beneficial to the formation of good physical adsorption between asphalt and aggregate on the contact surface. The test results demonstrate that the 3D profile roughness coefficient JRC3D has good sensitivity in the evaluation of the mesostructure changes of the aggregate surface texture and that this coefficient can be used as an evaluation index of the 3D roughness of road aggregate surfaces.

4. Conclusions

Based on the stylus contact measurement method and optical noncontact measurement method, the 2D and 3D quantitative evaluation methods of aggregate surface roughness were studied using the roughness coefficient JRC as the evaluation index. The following conclusions can be drawn.(1)By using contact measurement technology based on the principle of the needle-tracing method, the 2D profile curves of the aggregate surface can be obtained and can be used to quantitatively characterize the roughness coefficient of the 2D profile of the aggregate texture obtained by the characteristic parameter method. By using noncontact 3D scanning technology based on the principle of white-light interferometry, the 3D coordinate data of the analysis area of the aggregate surface can be obtained by the point cloud, which can be used to calculate the 3D roughness indexes of the aggregate surface after data processing.(2)When using the JRC2D index to evaluate the contour roughness of an aggregate surface, it is necessary to consider that the randomness of the selection of observation sections will lead to the anisotropy of the 2D roughness evaluation results. The evaluation index of 2D roughness cannot reflect the overall roughness of an aggregate surface.(3)Under the experimental conditions of this study, the average 3D roughness value of the limestone aggregate sample was the highest, that of granite was the lowest, and those of basalt and diabase were in the middle. When the weighted average method was used to evaluate the 3D roughness of the aggregate surface, the difference between the calculation results can be controlled below 10%, and the results, therefore, exhibit a good correlation. The calculation deviation of the two roughness characterization methods was mainly caused by the differences in the test methods and data accuracy.(4)The 3D profile roughness JRC3D calculated by the characteristic parameter method can be used to quantitatively evaluate the morphological characteristics of the aggregate particle surface. It is sufficiently sensitive to describe the microlevel structure of the aggregate texture and can be used as an evaluation index of the 3D roughness of road aggregate surfaces.(5)The calculation results of the 3D roughness of the aggregate surface are affected by the size effect. It is therefore suggested that when evaluating the surface roughness of the same batch of aggregates with the same particle size, the observation range of the analysis area should be fixed to avoid the deviation of the roughness evaluation results caused by the size effect.

Data Availability

The data used to support the findings of this study are available from the first author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The research was funded by the National Natural Science Fund (51708267).