Abstract
Dynamic modulus is a key parameter in the pavement structure design and analysis. A haversine loading waveform is the most widely recommended waveform in the laboratory to obtain the dynamic modulus of asphalt mixtures by test protocols, which may not completely represent all field loading conditions. The aim of this study is to investigate the suitable vertical compressive stress pulse waveforms at different depths under different pavement structures and to obtain the best fitting waveforms by using the least squares method. Specifically, the vertical compressive stress of different pavement structures was calculated by utilizing a three-dimensional finite element program incorporating with the measured wheel loads. The results revealed that the vertical stress pulse waveforms at different depths in asphalt layer were different within different pavement structure combinations. Generally, the square waveform fitted the vertical stress pulse waveform better for the shallow depth of the surface layer. The haversine waveform became more suitable while the depth increased. The bell waveform was better when the depth went deeper (i.e., 10 cm). In addition, the method of choosing waveforms at different depths of different pavement structures was provided, and the equation of calculating wavelength was recommended. Moreover, dynamic modulus under different loading waveforms were analysed through laboratory test, and the process for obtaining dynamic modulus of asphalt layer in any depth was presented.
1. Introduction
The mechanistic-empirical (M-E) is one of the fundamental methods that is routinely used for the asphalt pavement structural design and analysis. The most recent and successful achievement in the application and further enhancement of the M-E method is the development of the Mechanistic-Empirical Pavement Design Guide (M-EPDG) [1]. The dynamic modulus is one of the basic parameters as the input material property [1]. In laboratory, the dynamic modulus value of the hot mixture asphalt (HMA) is typically determined as a function of the test temperature, loading waveform, and frequency [2].
A haversine waveform is utilized as the loading waveform to obtain dynamic modulus corresponding to laboratory test methods [3–7]. In reality, the vertical compressive stress varies along with the pavement depth, which leads to the waveform of the vertical compressive stress pulse also vary at different depths within the asphalt concrete (AC) layer. The haversine loading waveform may not fully be simulative for the field loading conditions [8]. Based on the research result of Barksadale, the vertical compressive stress pulse could be fitted by the haversine waveform at pavement surface position [9]. Furthermore, according to the literature [10], the vertical compressive stress pulse waveform tended to approach a square waveform while the point of interest was close to the surface of the pavement and changed to a haversine as the depth increases. The shape of the vertical compressive stress pulse waveform approached a triangular one at deeper position in the AC layer. Moreover, haversine and bell waveforms could fit the vertical compressive stress pulse waveforms well for a moving vehicle [11].
Loading frequency can be calculated by loading time. Previous research has shown that loading time is closely related to the wavelength of vertical compressive stress, depth of interest, and vehicle speed. Hu et al. introduced an iterated method to calculate the loading time [10]. Brown derived an equation to calculate the loading time, and the relationship between loading time, t (s), depth, d (m), and vehicle speed, (km/h), was defined as follows [12]:
Huang proposed a haversine waveform to fit moving load in the VESYS and Kenlayer program, and loading time (di) was calculated as follows [13]:where ai is the tire contact radius for each axle type and s is the vehicle speed.
Several study results have shown that the tire contact area is not a simple circular shape, and the vertical tire-pavement contact pressure (TPCP) is not uniformly distributed [14–17]. The measured TPCP is nonuniformly distributed and changes with the discrepancy in axle load, tire inflation pressure, tire type, and tire tread patterns [15]. As a result, the mechanical responses are inaccurate if a circular uniformly distributed load is considered, which significantly affects the vertical compressive stress in the pavement structure. Therefore, the measured TPCP can be used to calculate more accurate vertical compressive stress and obtain an accurate vertical compressive stress pulse waveform. Furthermore, it provides an accurate loading waveform on obtaining dynamic modulus in the laboratory test and precise parameters for pavement structure design.
This study focuses on calculating the vertical compressive stress by using the measured wheel loads and then obtaining waveforms of the vertical compressive stress pulse. A three-dimensional finite element (3D-FE) program incorporated with the measured vertical TPCP and horizontal forces (i.e., lateral pressure) is utilized to calculate the vertical compressive stress. The objectives are to investigate the accurate loading waveforms, wavelength, and loading time at different depths in the AC layer under different combinations of layer thickness and AC modulus. Besides, this study demonstrates the method to acquire precisely dynamic modulus at different depths for the field AC layer through the laboratory test.
2. Methodology
2.1. 3D-FE Model
For mechanistic response analyses, it is more reasonable to adopt nonlinear viscoplastic models to fit the field pavement conditions [18–20]. However, it is proved that using multilayer linear elastic theory for modelling the AC layer, base layer, and subgrade is satisfactorily verified with measured field data [10, 21]. In addition, it is less time-consuming. Therefore, the 3D-FE multilayer linear elastic model was chosen in this study.
Figure 1 illustrates the 3D-FE model for analysis. The dimension of this model was 5 m × 10 m × 6 m. The boundary conditions were assumed to be as follows: (1) zero Z-directional displacement at the bottom, (2) zero X-directional displacement on both sides, and (3) zero Y-directional displacement on both sides. As previously mentioned, the aim of this study was to investigate the vertical compressive stress under the measured vertical TPCP and horizontal forces; thus, only the maximum vertical compressive stress within the pavement structure was modelled and analysed in the following part.

(a)

(b)
2.2. Loading Parameters
According to the literature [15], the longitudinal tread pattern tire is one of the most common types of tires in China, which was chosen for TPCP modelling. Moreover, 3D loads (vertical and horizontal loads) were used to represent the tire-pavement pressure of the moving load under a constant speed in this study. The obtained moving load was uniformly distributed on the node of each meshed unit in the 3D-FE model. The simplified 3D-FE model for tire-pavement contact area and the loading pattern on the 3D-FE model were illustrated, as shown in Figure 2. As shown in Figure 2(c). It should be noted that the vertical and horizontal TPCP directions acting on the nodes of each rib is completely symmetrical. Total weight of the rear axle of 100 kN, tire pressure of 600 kPa, and a single wheel load of 25 kN were implemented to the model. The measured vertical TPCP was consistent with the tire tread pattern, which was used to compute the vertical and horizontal stress [15]. Based on the available literature, the lateral tire contact stress was generally in the perpendicular direction and approximately 15% to 50% of the maximum vertical compressive stress [22–25]. Besides, lateral tire contact stress was taken as 30% of the maximum vertical compressive stress. The longitudinal stress was considered to be 12% of the maximum vertical compressive stress [25].

(a)

(b)

(c)
2.3. Pavement Structure
Typical AC layer thickness, namely, 2.5 cm, 3.8 cm, 5 cm, 7.5 cm, 10 cm, 15 cm, 20 cm, and 30 cm, was investigated. The thickness of two base layers was 20 cm. The literature [10] showed that the modulus ratio (R) between the layer of interest and the immediate succeeding layer below has a significant influence on the vertical compressive stress waveform. R values were considered as 0.02, 0.05, 0.5, 1, 5, 20, and 50 in this study. It should be noted that when calculating the vertical stress, the thickness equivalent conversion of the multilayer asphalt pavement structure was still valid [9, 10]. On this basis, in order to simplify the model, only one layer of AC was selected. The minimum and maximum moduli of the AC layer are 1,000 MPa and 10,000 MPa, respectively. The moduli of two base layers were set the same, and the value was between 10 MPa and 200,000 MPa. Similarly, the moduli of subgrade were between 10 MPa and 1,000 MPa. The modulus and thickness of various pavements were considered which covered most of the commonly used pavement structures. More than 1,000 cases were analysed in the study.
3. Analysis of Vertical Compressive Stress Pulse Waveform
The vertical compressive stress pulse is different at various depths, and three shapes of waveforms were utilized to fit the curves. Moreover, the least squares method was used to choose the optimal fitting waveform.
3.1. Equations for Fitting the Vertical Compressive Stress
Researchers utilized square, haversine, triangle, and bell waveforms to fit the vertical compressive stress pulse waveform of the AC layer. The results demonstrated that the vertical compressive stress pulse waveform changed at different depths in the AC layer [8–11]. In addition, according to the analysis results, the triangle waveform could not fit the vertical compressive stress pulse waveform appropriately. However, square, haversine, and bell waveforms could fit the vertical compressive stress pulse waveform at different depths in different pavement structures. As a result, square, haversine, and bell waveforms were selected. Fitting equations for the square, haversine, and bell waveforms are listed in the following equations, respectively:where d (m) is the length of the vertical compressive stress pulse and x (m) is the horizontal distance between the point of interest and the centre of the loading.
3.2. Process of Choosing the Best Fitting Waveform
As previously stated, three different waveforms were utilized to fit the vertical compressive stress pulse waveforms. The method of normalizing the vertical compressive stress was beneficial for obtaining suitable fitting curves. The steps for selecting the best fitting waveform curves were as follows:(1)The vertical compressive stress was calculated by the 3D-FE model at each depth. The maximum vertical compressive stress distribution along the traffic direction was selected as the representative for analysis.(2)The normalized vertical compressive stress pulse waveform was calculated by the ratio of vertical compressive stress value to its maximum value along the traffic direction at the same depth.(3)The normalized vertical compressive stress pulse waveform were fitted by using equations (3)–(5), respectively.(4)The optimal curve was chosen by the least squares method.
3.3. Vertical Compressive Stress Pulse Waveform in the Same Pavement Structure
Due to the length of the article, in this part and the following part, only parts of the fitting results were demonstrated.
S-Z indicates the normalized vertical compressive stress pulse waveform. S and D indicate the wavelength of the square waveform, and the best fitting waveform at the given depth was the square waveform. H and D represent the wavelength of the haversine waveform, and the best fitting waveform at the given depth was the haversine waveform. B and D indicate the wavelength of the bell waveform, and the best fitting waveform at the given depth was the bell waveform. Z (cm) is the vertical depth between the point of interest and the pavement surface.
The vertical compressive stress at different depths in the AC layer was calculated, and the pavement structure was as follows: the thickness of the AC layer was 20 cm, and the moduli of AC layer, base, and subgrade were 1,000 MPa, 20,000 MPa, and 50 MPa, respectively. The fitting results were at the depths of 0 cm, 5 cm, 10 cm, and 17.8 cm, as shown in Figure 3.(1)The vertical compressive stress pulse varied with the depth of interest changed. Meanwhile, the length of the vertical compressive stress pulse increased as the depth increased. Optimal fitting waveform for the vertical compressive stress was not always haversine. As the depth of interest increased, the optimal fitting waveform was square at shallow position, in Figure 3(a). Then, it changed to haversine, in Figure 3(b). On a comparison basis, the bell waveform was the best fitting shape curve at deeper position (Z = 17.8 cm), in Figure 3(d). In conclusion, the depth was a main factor influencing the shape and wavelength of vertical compressive stress pulse.(2)Two waveforms could be used to fit and define the vertical compressive stress at a particular depth, in Figure 3(c). Both the haversine and bell waveforms could be utilized to fit the vertical compressive stress waveforms at Z = 10 cm; additionally, the wavelengths of the two kinds of waveforms were close.

(a)

(b)

(c)

(d)
3.4. Vertical Compressive Stress Pulse Waveform in Different Pavement Structures
Four groups of modulus combinations of the surface AC layer, base, and subgrade were selected for this analysis at two different depths. R values of four pavement structures were 0.5, 0.5, 1, and 20, respectively. The fitting results of vertical compressive stress pulse are shown in Figure 4. In the figure, the following can be observed: Figure 4(a): the thickness of the AC layer was 15 cm, moduli of the pavement structure were 1,000 MPa, 2,000 MPa, and 40 MPa for the AC layer, base layer, and subgrade, respectively, and the position of interest was at a depth of 10 cm. Figure 4(b): the thickness of the AC layer was 15 cm, moduli of the pavement structure were 10,000 MPa, 20,000 MPa, and 100 MPa for the AC layer, base layer, and subgrade, respectively, and the position of interest was at a depth of 10 cm. Figure 4(c): the thickness of the AC layer was 30 cm, moduli of the pavement structure were 1,000 MPa, 1,000 MPa, and 30 MPa for the AC layer, base layer, and subgrade, respectively, and the position of interest was at a depth of 20 cm. Figure 4(d): the thickness of the AC layer was 30 cm, moduli of the pavement structure were 10,000 MPa, 500 MPa, and 20 MPa for the AC layer, base layer, and subgrade, respectively, and the position of interest was at a depth of 20 cm.

(a)

(b)

(c)

(d)
In Figure 4, the following conclusions could be drawn:(1)Comparing Figures 4(a), 4(c), and 4(d) and Figures 4(b), 4(c), and 4(d), the vertical compressive stress pulse waveform changed with different R values, thickness of pavement structures, and different depths; additionally, wavelength was not the same.(2)In Figures 4(a) and 4(b), among different pavement structures, when the R value and the thickness of AC layer were the same, regardless of the modulus of AC value, the vertical compressive stress pulse waveform was the same. Moreover, haversine and bell waveforms could fit the vertical compressive stress pulse waveform well at a certain depth, and the wavelength of two waveforms was close. The R value was one of the key factors that affect the waveform and wavelength, and the modulus of subgrade had a little effect on waveform and wavelength.(3)In Figures 4(c) and 4(d), the bell waveform fitted the vertical compressive stress pulse waveform well in deeper position of pavement structure, i.e., at a depth of 20 cm. Wavelength increased with the increase in the R value within the same thickness of the pavement structure.
4. Recommended Waveforms and Wavelength for Each Depth
Numerous computations and analyses were conducted in this study. Based on these analyses, the appropriate loading waveform and wavelength were proposed for different depths in the AC layer. The depth of interest and the R value were the two main factors that affected the waveforms and wavelength. The following part stated the method to obtain appropriate loading waveform and wavelength at any depth in the AC layer.
HAC (cm) represented the thickness of the AC layer, h (cm) represented the depth of the position of interest, the R value was the modulus ratio between the layer of interest and the immediate succeeding layer below, and D (cm) was the wavelength of the vertical compressive stress pulse.
5. Recommended Proper Waveform for Each Depth
According to more than 1,000 cases of computation and analyse results, appropriate loading waveform was suggested. The following steps were the methods to get a proper waveform at different depths in different pavement structures:(1)A square waveform was suggested for a pavement depth within the range 0 < h < 2.5 cm, regardless of the AC layer thickness.(2)A haversine waveform was recommended for a pavement depth range of 2.5 cm ≤ h < 5 cm when 2.5 cm ≤ HAC < 5 cm.(3)A haversine waveform was proposed when 5 cm ≤ HAC ≤ 10 cm and for a pavement depth range of 2.5 cm ≤ h < 0.5 HAC. Furthermore, for 0.5HAC ≤ h ≤ HAC and 0.02 ≤ R ≤ 5, a haversine waveform was proposed, and bell waveform was applicable for the other cases.(4)A haversine waveform was found to be appropriate when 10 cm < HAC, at a position of 2.5 cm ≤ h ≤ 5 cm. For the other case scenarios, the following steps were used to select the appropriate loading waveform:(4.1)A haversine waveform was recommended when 10 cm < HAC < 20 cm. 5 cm < h ≤ 7.5 cm. Moreover, when h > 7.5 cm, both the haversine and bell waveforms could satisfactorily fit the vertical stress pulse when R ≤ 1. However, when R > 1, the bell waveform exhibited better curve fitting tendency than the haversine waveform.(4.2)A haversine waveform was recommended when 20 cm ≤ HAC and 5 cm < h ≤ 10 cm. For other depths, the bell waveform was more suitable.
5.1. Recommended Wavelength for Each Depth
Depending on the optimal fitting waveforms that were recommended at different depths, simultaneously, wavelength could be obtained. The wavelength could be divided into four conditions for analysis, namely, (1) condition 1: 0 < h < 2.5 cm; (2) condition 2: 2.5 cm ≤ HAC ≤ 3.8 cm; (3) condition 3: 3.8 cm < HAC and h ≤ 0.5757HAC + 0.0424; and (4) condition 4: else. The fitting equations of wavelength for each condition were proposed, as presented in Table 1. The wavelength was mainly affected by the depth of interest in the shallow position (i.e., 3.8 cm) based on the equations of wavelength.
The wavelength was 21.7 cm for condition 1, equal to the maximum length of the tire tread. Figure 5(a) illustrates the wavelength under conditions 2 and 3, and the wavelengths could be fitted accurately by equations (∗) and (∗∗), respectively. As presented in Figures 5(b) and 5(c), the wavelength corresponds to different R and h. Besides, equation (∗∗∗) fitted the wavelength of condition 4 properly.

(a)

(b)

(c)
6. Dynamic Modulus of Different AC Layers
The dynamic modulus of AC at various temperatures and frequencies could easily be obtained under different loading waveforms by the laboratory test [2]. In this section, the method and example of obtaining dynamic modulus for random depth of the AC layer was presented.
6.1. Method for Obtaining Dynamic Modulus
Figure 6 illustrates a typical pavement structure. The AC layer was divided into several sub-AC layers, the thickness was hACi, the modulus was ESub-ACi, and the depth of interest was hi.

The modulus was calculated from the lower position to the upper position. The iteration method was used to obtain the loading waveforms, wavelength, frequency, and dynamic modulus. The difference between two calculated modulus values was less than 1%, which meant the modulus was appropriate. The method of obtaining dynamic modulus depends on the waveform and wavelength, as shown in Figure 7.

Loading time (T) could easily be calculated by using equation (6). Loading frequency (f) could be expressed by using equation (7):where D is the wavelength and is the vehicle speed.
In addition, dynamic modulus curves under different loading waveforms at a certain temperature could be fitted by using the following equation:
6.2. Example of Obtaining Dynamic Modulus
The common pavement structure in China and the probable loading waveform for each layer are presented in Table 2.
All the test samples of AC-13/20/25 were molded a target air void of 4 ± 1%, as specified by “Technical specification for construction of highway asphalt pavements” [26]. As shown in Figure 8, the dynamic moduli of AC-13, AC-20, and AC-25 under different loading waveforms at the temperature of 21°C were obtained. The dynamic modulus curves under different loading waveforms were different. Within the same test condition, dynamic modulus values under the square loading waveform were much smaller than those under haversine and bell loading waveforms. Each E∗ curve was fitted as presented in equations (9)∼(14) ( is the dynamic modulus, and f is the loading frequency):

The pavement structure has, as shown in Table 2, HAC = 18 cm. According to “Specifications for design of highway asphalt pavements” [27], the modulus of the base layer: EB = 18,000 MPa and the modulus of the subgrade: Es = 50 MPa. In general, considering the minimum speed of China’s high-grade highway is 60 km/h, the speed was set as 72 km/h = 20 m/s.
Based on the methodology of obtaining dynamic modulus for the depth of interest as presented in Figure 7, the dynamic modulus, wavelength, frequency, and loading waveforms at the depth of h1 = 2 cm, h2 = 7 cm, and h3 = 14 cm were calculated as follows:(1)h3 = 14 cm. The modulus at a depth of 14 cm was assumed as E3 = 20,000 MPa. 1 < R3 = E3/EB. Loading waveform: bell waveform. Via equation (∗∗∗), wavelength: D3 = 61.870 cm, loading time: T3 = 0.0309 s, and loading frequency: f3 = 32.326 Hz. Via equation (12), E3-1 = 13,945 MPa. R3-1 = E3-1/EB < 1. Loading waveform: haversine and bell waveforms. Via equation (∗∗∗), wavelength: D3-1 = 60.138 cm, loading time: T3-1 = 0.0301 s, and loading frequency: f3-1 = 33.257 Hz.(1.1)Loading waveform: bell waveform. Via equation (12), E3-1-1 = 14,001 MPa. Then, |E3-1-1-E3-1|/E3-1 ≤ 1%; therefore, E3 = E3-1-1 = 14,001 MPa.(1.2)Loading waveform: haversine waveform. Via equation (9), E3-2-1 = 16,872 MPa. R3-2-1 = E3-2-1/EB < 1. Via equation (∗∗∗), wavelength: D3-2-1 = 60.980 cm, loading time: T3-2-1 = 0.0305 s, and loading frequency: f3-2-1 = 32.798 Hz. Via equation (∗∗∗), E3-2-1-1 = 16,837 MPa. Then, |E3-2-1-1-E3-2-1|/E3-2-1 ≤ 1%; therefore, E3 = E3-2-1-1 = 16,837 MPa. Therefore, when the loading waveform was bell, E3 = 14,001 MPa, wavelength: D3 = 60.138 cm, and loading frequency: f3 = 33.257 Hz. When the loading waveform was haversine, E3 = 16,837 MPa, wavelength: D3 = 60.980 cm, and loading frequency: f3 = 32.798 Hz.(2)h2 = 7 cm. The modulus at the a depth of 7 cm was assumed as E2 = 15,000 MPa. Haversine was the loading waveform at the depth of 7 cm.(2.1)When E3 = 14,001 MPa, 1 < R2-1 = E2/E3. Via equation (∗∗∗), wavelength: D2-1 = 45.650 cm, loading time: T2-1 = 0.0228 s, and loading frequency: f2-1 = 43.812 Hz, and via equation (10), E2-1 = 15,864 MPa. 1 < R2-1-1 = E2-1/E3. Via equation (∗∗∗), wavelength: D2-1-1 = 45.830 cm, loading time: T2-1-1 = 0.0229 s, and loading frequency: f2-1-1 = 43.639 Hz, and via equation (10), E2-1-1 = 15,855 MPa. Then, |E2-1-1-E2-1|/E2-b-1 ≤ 1%; therefore, E2 = E2-1-1 = 15,855 MPa.(2.2)When E3 = 16,837 MPa, R2-2 < 1. Via equation (∗∗∗), wavelength: D2-2 = 45.121 cm, loading time: T2-2 = 0.0226 s, and loading frequency: f2-2 = 44.325 Hz. Via equation (10), E2-2 = 15,892 MPa. R2-2-1 = E2-2/E3 < 1. Via equation (∗∗∗), wavelength: D2-2-1 = 45.278 cm, loading time: T2-2-1 = 0.0226 s, and loading frequency: f2-2-1 = 44.172 Hz, and via equation (10), E2-2-1 = 15,884 MPa. Then, |E2-2-1-E2-2|/E2-2 ≤ 1%; therefore, E2 = E2-2-1 = 15,884 MPa. Therefore, when E3 = 14,001 MPa, E2 = 15,855 MPa, wavelength: D2 = 45.8306 cm, and loading frequency: f2 = 43.639 Hz; When E3 = 16,837 MPa, E2 = 15,884 MPa, wavelength: D2 = 45.278 cm, and loading frequency: f2 = 44.172 Hz;(3)h1 = 2 cm. h1 < 2.5 cm, the square waveform was suitable, wavelength: D1 = 21.7 cm, loading time: T1 = 0.00868 s, and loading frequency: f1 = 92.166 Hz, and via equation (14), E1 = 2,728 MPa.
7. Conclusions
The 3D-FE model incorporating with the measured TPCP was utilized for calculating the vertical compressive stresses in AC layers. Based on extensive calculation and numerical fitting, proper loading waveforms and new wavelength calculation equations were proposed for the laboratory dynamic modulus test. The conclusions could be drawn as follows:(1)Vertical compressive stress pulse is not always the same. The vertical compressive stress pulse tends to approach a square waveform when the position of interest is near the pavement surface. Then, it approaches to haversine with the increasing of depth. At greater pavement depths (i.e., 10 cm), the shape of the vertical compressive stress pulse changes to a bell waveform.(2)The square waveform could be suitable merely with depths less than 2.5 cm; however, haversine and bell waveform could be more widely applicable. Additionally, haversine and bell waveforms could be suitable in some depth simultaneously.(3)The waveform and wavelength vary at different positions beneath the pavement surface. The method to select a suitable waveform in the AC layer has been given. The equations to obtain loading wavelength within different loading waveforms at different depths are proposed.(4)According to the method of recommendation of suitable waveforms and wavelength, R and depth are the two main factors that influence the waveform and wavelength in the AC layer.(5)The process to calculate the accurate dynamic modulus at different depths for the field AC layer under different waveforms and loading time through the laboratory test is presented.
Data Availability
Previously reported (pictures or tables) data used to support this study are included within the article. These prior studies (and datasets) are cited at relevant places within the text as references [15].
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors wish to thank Dr. Yalong Pi from A&M University and Dr. Yao Zhang from Southeast University for their suggestions for this study. This study was funded by the National Nature Science Foundation of China (no. 51878168) and Open Research Fund of National Engineering Laboratory for Advanced Road Materials (no. NLARM-ORF-2018-01).