Numerical investigation of the plastic deformation behaviour of graded crushed stone

A numerical model of the dynamic triaxial test of graded crushed stone was established based on DEM (Discrete Element Method) to study its dynamic characteristics. The influence of test conditions on simulation results was analysed through numerical simulation. A method for determining the test conditions was proposed, and the reliability of the simulation was verified. We studied the accumulation rule and failure standard of permanent deformation. We then determined the critical failure stress level. The prediction model of permanent deformation cumulative failure was then established. When the calculation time step is greater than 1E−4 s per step, the stability of the dynamic triaxial numerical test of graded crushed stone is good. The plastic deformation of the simulated specimen tends to be stable under 10,000 dynamic loading cycles. When the specimen height and diameter are greater than 20 and 10 cm, respectively, the specimen size has little influence on the simulated axial strain. The recommended specimen size is 10 cm (Ф) × 20 cm (h). The action time curve results are consistent with indoor measurement results, which proves the simulation’s reliability. The critical failure stress is approximately linearly correlated with the confining pressure, and the cumulative failure equation of plastic deformation is established.


Introduction
Graded crushed stone is often used in road pavement bases. It is immensely significant to study the dynamic and deformation characteristics of graded crushed stone to improve pavement durability [1]. Janoo and Bayer experimentally analysed the properties of graded crushed stone raw materials, divided the evaluation index and grade of raw materials and analysed the influence of different factors on the dynamic modulus and shear strength of graded crushed stone through triaxial tests [2]. Brown and Chan studied the nonlinear characteristics of graded crushed stone using many laboratory tests and fitted the recommended regression equation and parameter range through a large amount of data [3]. Rada and Witczak studied the behavior law of plastic deformation and rebound deformation of graded gravel through an indoor dynamic triaxial test [4]. Azam  accumulate to the limit value when the number of loading actions tends to be infinite [5]. Werkmeister et al. established a calculation model of graded crushed stone pavement using finite element software and combined it with an indoor dynamic triaxial test to calculate the permanent deformation. The orthogonal method was used to analyse the change law of water content, gradation, raw material characteristics, load and stress level [6]. Wang et al. studied the development law and distribution state of plastic deformation of graded crushed stone base under long-term vehicle load [7,8]. Wang et al. analysed the dynamic deformation characteristics of graded gravel mixtures [9]. Yang et al. explored the influence of fine particle content on the dynamic deformation behaviour and cumulative plastic strain of graded gravel [10]. Because of the obvious nonlinear mechanical properties and discontinuous medium structure of graded crushed stone, the relationship between mechanical parameters, experimental conditions and deformation performance of graded crushed stone cannot be theoretically analysed in detail owing to the limitation of current mechanical theory. In addition, the traditional indoor test method has low efficiency and high cost, so the research depth and breadth of the traditional indoor test method for studying the dynamic characteristics of graded gravel are limited. To solve the abovementioned problems, the DEM is widely used for studying mechanical properties of road materials, which provides a theoretical basis and method support for the study of dynamic characteristics of graded crushed stone [11][12][13][14][15]. Therefore, based on DEM and particle flow theory, the dynamic triaxial numerical test method of graded crushed stone is constructed in this work to study the plastic deformation law and critical failure stress of graded gravel, and the cumulative failure equation of plastic deformation of graded gravel is proposed.

Materials
The raw material used in this work was Yinghu limestone from Ankang, Shaanxi Province. The apparent density of the raw material is given in Table 1; the mineral aggregate gradation is given in Table 2. A, B and C are the median values of skeleton dense gradation, continuous gradation 1 and continuous gradation 2, respectively, in JTG D50-2015 code for the design of highway asphalt pavement (hereinafter 'the specification').

Test conditions.
(1) Test instrument. In this paper, UTM-25 (universal testing machine) was used to conduct the indoor dynamic triaxial test, as shown in Fig 1. (2) Loading time and loading waveform. Loading time and loading waveform are important parameters for dynamic testing of permanent deformation of loose aggregates. The load waveform and loading time of the dynamic test should be determined according to the response law of the pavement under the action of the traffic load. When the load is quite far away from a certain point on the road, the point's stress is 0, and when the load directly acts on the point, the stress reaches the maximum. It is assumed that with change in time, the load strength exhibits the law of half-sine function; the expression for load strength is given by formula (1): where d is the load action time (s) and q is the load strength (MPa). When the load is far away or L t ¼ � d 2 , then L t = 0; when the load directly acts on t = 0, the load strength reaches the maximum value.
The time of load action depends on vehicle speed v and tyre contact area radius a. Previous studies have shown that when the load distance is six times that of the tyre grounding area radius, it can be considered that the load does not affect the road surface. The specific expression is given by formula (2): where v is the vehicle speed (m s−1) and a is the tyre contact area radius (m).
The smaller the vehicle speed v is, the longer the load acts on the road surface, and the greater the impact is. The known v is 64 km�h −1 , and the load action time is 0.1 s. Therefore, considering the most adverse effect, this paper selects the half-sine wave with dynamic load action time t z = 0.1 s and dynamic load action time interval t j = 0.9 s and cycle 1.0 s. The load waveform curve is shown in Fig 2. (3) Confining pressure. Before the dynamic triaxial test, ANSYS was used to calculate the stress distribution of each structural layer of an asphalt pavement when the normal tyre ground pressure is 0.7-1.1 MPa within the elastic-plastic theory range. The horizontal stress distribution in the granular layer is generally between 20 and 150 kPa. To simulate the actual confining pressure of graded crushed stone, the maximum σ3 is slightly greater than the maximum effective stress actually borne by the structural layer. The minimum σ3 should not be less than the pressure experienced by the structural layer. The confining pressure ranges from 50 to 150 kPa (50, 75, 100, 125 and 150 kPa).
(4) Deviator stress. The ultimate strength in static triaxial tests under different confining pressures was divided into five stress levels to compare the plastic deformation characteristics of crushed stones under long-term dynamic loading. For the elastic-plastic loose aggregate of graded crushed stone, on the one hand, its stress-strain curve exhibits nonlinear change: the greater the load, the more prominent the plastic performance of mixture. On the other hand, the loading stress level selection should consider the stress state of a graded gravel layer asphalt pavement structure under general vehicle loads. Therefore, for the repeated loading test of permanent deformation of graded crushed stone, the stress level was selected between 0.1 and 0.9, and the stress level S = 0.1, 0.3, 0.5, 0.7 and 0.9 was selected in this work. S is the ratio of maximum cyclic stress to shear strength, as given by formula (3): where σ max is the maximum cyclic stress (kPa) and σ s is the ultimate strength of the static triaxial test (kPa).
(5) Number of loading actions. The number of loading actions is an important factor for studying the permanent deformation behaviour of graded crushed stone. Based on existing research, the number of loadings in the repeated loading test in this paper is 10,000 [16][17][18][19][20][21][22][23][24].
2.2.2 Test procedure. This test refers to the AASHTO DESIGNATION T307-99 test method, and the specific steps are as follows: 1. After the specimen is prepared, place it on the test bench with both ends of the rubber membrane tightly tied to the base; install the pressure chamber cover, connect the loading shaft and the actuator and inject silicone oil until the entire pressure chamber is filled.
2. Install the axial differential displacement metre and connect the piezoelectric pressure sensor.
3. Set the loading scheme and set the confining pressure until the required value is stable for 30 min.
4. When the deformation of the specimen is stable under static axial stress, the dynamic loading module is used to apply the dynamic load, and the stress control method is used; that is, the deviator stress remains unchanged during the test. When the specimen is loaded 10,000 times, the deformation stops. Remove the confining pressure, remove the test piece and prepare for the next round of tests.

5.
Change the next specimen, change the deviator stress and confining pressure, repeat steps (1)-(4) and sort the data.

Construction of the numerical test method
3.1.1 Generation of graded crushed stone. The DEM built-in command 'wall' was used to generate four 'walls' to form a closed rectangular area for simulating the test mould. Assuming there are n types of crushed stones in the range of particle size to form the mixture according to a certain gradation, the surface mass m of the specimen in a two-dimensional plane state equivalent to M is as follows: where M is the total mass of the test piece (g), , D is the diameter of the test piece (cm) and h is the height of the test piece (cm), m is the surface quality of a specimen (g) and S is the two-dimensional mapping area of the specimen (cm 2 ).
We calculate the two-dimensional mapping area S i of the type i gravel as where S i is the two-dimensional mapping area of type i gravel (cm 2 ) and ρ i is the density of It can be seen from Fig 3 that because the centre coordinates of aggregate particles are generated randomly, particle overlapping will inevitably occur in the finite space of the simulated test mould, and significant additional stress will be generated between particles, as shown in Fig 3(A). To eliminate this adverse effect, the generated 'wall' should keep a certain distance from the particles so that they have enough space to dissipate their own additional stress at the beginning of the calculation cycle. See Fig 3(B) for the moulded specimen after additional stress dissipation.

Formation of the indenter and the diaphragm.
A row of spheres is generated on both sides of the 'wall' on the top and bottom surfaces. The sphere and sphere are connected as a whole by bond model. The dynamic triaxial test indenter of graded gravel is simulated. The original two 'wall' are given acceleration to the horizontal two sides, and the two 'wall' are pulled out to both sides along the horizontal direction. We kept the left and right 'walls' still and simulated the dynamic triaxial test rubber membrane of graded gravel to simulate the application of confining pressure, as shown in

Sensitivity and determination of test conditions
To determine the standard test conditions for the dynamic triaxial numerical test of graded crushed stone, the sensitivity of the results of the dynamic triaxial numerical test of graded crushed stone to the calculation time step, loading times and specimen size was investigated in this work. See Fig 5 for the loading diagram of the dynamic triaxial numerical test of graded crushed stone.
3.2.1 Influence of calculation time step. The calculation time step reflects the relationship between computer simulation time and the actual time. In principle, the calculation time step should not exceed the default value of the computer. However, the larger the calculation time step, the more information contained in the same time step and the greater the amount of calculation. It is proven that the selection of the calculation time step has a great influence on the stability of the numerical test [25]. Different calculation time steps dt were selected to stabilise the virtual test piece, and the influence of the calculation time step on the stability of the numerical test was analysed to select a reasonable calculation time step. See Fig 6 for the simulation diagram. It can be seen from Fig 6 that when the calculation time step is greater than 1E−4 s per step, the numerical test effect is stable and the simulated indenter has no major disturbance. When the calculation time step is 1E−3 s per step, the upper simulated indenter flies away from the virtual specimen range. In the specimen range, the lower simulated indenter tilts, and a small number of aggregates fly out to both sides of the specimen. When the calculation time step is 1E−2 s per step, the upper and lower indenters fly away from the virtual specimen range, and many aggregate flies out to both sides of the specimen. To ensure the stability of the numerical experiment while considering the computational efficiency of the numerical experiment, the numerical experiment calculation step was taken as 1E−4 s per step.

Influence of load times.
The dynamic triaxial simulation test method was used to study the influence of the number of loading actions on the permanent deformation of the graded crushed stone. The microcontact force state of the virtual sample of the graded crushed stone upon application of a dynamic load is shown in Fig 7. The number of loading actions is 100,000. The relationship between the permanent strain ε of the graded crushed stone and the number of loading actions N is shown in Fig 8.  Fig 7 shows that when the virtual sample of graded crushed stone is not subjected to dynamic loading, the direction of load transmission is not obvious. The distribution of the contact force network is disorderly in the entire specimen. When the virtual sample of graded gravel is subjected to dynamic load, the contact force concentration phenomenon appears at the top and bottom of the specimen. By contrast, the distribution of contact force networks in the virtual specimen does not change. With increase in the number of loadings, we found that the contact force is transmitted along the main skeleton in a strip form. When the virtual specimen of graded crushed stone is in the dynamic load interval stage, the contact stress concentration inside the specimen increases and the stress concentration at the top and bottom decreases. It can be seen that the dynamic load not only acts on the surface of the virtual specimen but also has a significant effect on the interior of the specimen.
It can be seen from Fig 8 that the axial permanent strain of graded crushed stone materials increases significantly with the increase in loading times (especially the first 1000 times). When the loading times reach 10,000, the axial permanent strain growth is no longer obvious and gradually tends to be stable. The continuous increase in loading times has little effect on the axial permanent strain. Therefore, to simplify the test and comprehensively monitor the influence of loading times on axial permanent strain, the dynamic load action times of 10,000 was adopted in the indoor permanent deformation test of graded crushed stone.

Influence of specimen size.
(1) Height of the test piece. Fig 9 shows the change rule of ε 10000 of graded crushed stone with respect to the height of the test piece (h). The stress condition of the test is that the confining pressure is 50 kPa, the diameter of the specimen is Ф = 10 cm, and the nominal maximum particle size of the aggregate is abbreviated as D max.

PLOS ONE
Plastic deformation behaviour of graded crushed stone It can be seen from Fig 9 that the influence of the height (h) of the graded crushed stone specimens with different maximum nominal sizes D max on ε 10000 is similar. With increase in specimen height, ε 10000 decreases first and then tends to be stable. When D max � 31.5 mm, when ε 10000 tends to be stable, H � 20 cm; when D max = 37.5 mm, when ε 10000 tends to be stable, h = 15 cm. Therefore, when h � 20 cm, the gravel particle size tends to be stable, and the size of the specimen has little effect on the numerical test results.
For further analysis, this paper presents the appropriate displacement field of the virtual specimen at different specimen heights, as shown in Fig 10.  It can be seen from Fig 10 that when h � 15 cm, the aggregate in the virtual specimen has obvious displacement, and the aggregate in the middle of the test mould moves violently and accumulates towards the middle of the specimen. This is because the upper and lower parts of the aggregate are disturbed when the indenter acts on the aggregate, resulting in opposite vertical displacement. When 20 cm < h < 25 cm, because of the increase in specimen height, the

PLOS ONE
Plastic deformation behaviour of graded crushed stone displacement vector has been partially consumed, which results in the aggregate in the middle of the specimen moving laterally but not to the middle. When h � 25 cm, a 'weak influence area' appears in the middle area of the test mould, and the aggregate displacement is minimal in this area. This is due to the continuous increase in the specimen height and because the displacement vector has been exhausted before the opposite moving aggregates meet, which prevents the indenter from disturbing the aggregate at this time. The axial strain only comes from the compaction of the upper and lower parts of the aggregate, so the axial strain is small.
To sum up, when the height of the specimen is small, the transverse displacement of the aggregate and the axial strain are too large, which obviously cannot give full play to the performance of the graded crushed stone material. When the height of the test piece is too large, there is a 'weak influence area' in the middle of the test mould. Therefore, with respect to the height of the test piece, it is necessary to ensure that the axial strain is within the reasonable range and avoid the occurrence of a 'weak influence area'. It is suggested that the specimen height h should be 20-25 cm.
(2) Diameter of the test piece. Fig 11 shows the change rule of ε10000 of graded crushed stone with respect to the diameter of the test piece (Ф). The stress condition of the test is that the confining pressure is 50 kPa and the height of the specimen is h = 20 cm.
It can be seen from Fig 11 that the influence of the diameter of graded crushed stone specimens with different nominal maximum particle sizes D max on ε 10000 is similar. With the increase in specimen diameter, ε 10000 first decreases and then tends to be stable. When the maximum particle size D max � 9.5 mm, when ε 10000 tends to be stable, Ф � 15 cm; when the maximum particle size D max � 9.5 mm, when ε 10000 tends to be stable, Ф = 10 cm. At this time, when the particle size is larger than or equal to 15 cm, the test results are slightly affected by the particle size.
For further analysis, the microdisplacement vector field of the virtual specimen with different specimen diameters is shown in Fig 12.  It can be seen from Fig 12 that the aggregate particles in specimens with different diameters have different degrees of lateral movement. When the diameter is less than or equal to 15 cm, the upper part of the aggregate particles in the specimen clearly exhibits a 45˚oblique downward shear plane. Because of the constraints of the test moulds on both sides of the specimen, the vertical displacement of the particles on both sides of the specimen is large. When the specimen diameter is greater than or equal to 20 cm, the transverse movement of the aggregate particles mainly occurs in the middle of the specimen, and the displacement of the aggregate particles on both sides of the specimen is minimal, which is due to the action of the indenter. However, because of the large diameter of the specimen, the displacement vector is largely dissipated before the displacement vector is transferred to the two sides of the test mould. The displacement of the aggregate on both sides of the specimen is minimal, and the 'weak influence area' appears on both sides of the specimen, so the axial strain is small.
To sum up, when Ф � 20 cm, there is a 'weak influence area' on both sides of the specimen. Although the axial strain is stable, it cannot give full play to the role of mould confining pressure. Therefore, with respect to the diameter of the test piece, we should not only ensure the role of the test mould confining pressure but also avoid the occurrence of a 'weak influence area'. It is suggested that the diameter of the specimen should be 10-20 cm.
When the height of the test piece is h � 20 cm and the diameter of the test piece is Ф � 10 cm, the influence of the gravel particle size on the test result tends to be stable, and the size effect of the test piece has little influence on the numerical test results. At the same time, considering the simulation degree and calculation speed of the numerical test, it is suggested that the size of the specimen should be 10 cm (Ф) × 20 cm (h).

Reliability verification
The reliability of the numerical simulation method was verified by taking a class A configuration in Table 2 as an example.

Simulation of the trial mould.
According to the research on the dynamic triaxial numerical test of graded crushed stone, the size of the test model should be 10 cm (Ф) × 20 cm (h). The DEM command 'wall' was used to generate two vertical 'walls' with a length of 20 cm and two horizontal 'walls' with a length of 10 cm to form a closed rectangle to simulate the test mould.

Generation of graded crushed stone.
The equivalent conversion method was used to convert the aggregate volume into a two-dimensional mapping area in numerical experiments. Here, the two-dimensional mapping area of 19-31.5 mm gravel is calculated as an example. The DEM command 'ball' was used to generate particles with diameters ranging from 19 to 31.5 mm. When the total area reached 44.9 cm 2 , the particles were stopped. By analogy, the generation of graded crushed stone is completed.

Reliability verification.
The indoor dynamic triaxial test results and simulation results of graded crushed stone are shown in Fig 13. It can be seen from Fig 13 that the maximum error of the numerical and laboratory tests with a confining pressure of 50 kPa is 7%, and the average error is 3%. The maximum error in the numerical and indoor tests with a confining pressure of 100 kPa is 5%, and the average error is 2%. The maximum error in the numerical and indoor tests with a confining pressure of 150 kPa is 8%, and the average error is 1%. The results of the numerical and laboratory tests are in good agreement, which proves that the dynamic triaxial numerical test method is reliable.  Fig 14; the confining pressure is taken as 50 kPa.
It can be seen from Fig 14 that the relationship between the plastic deformation curve and the number of loading actions varies with the stress. The curve of the relationship between the axial strain of graded crushed stone and the number of loading actions n is mainly distributed in three types of deformation regions: stable section, critical section and failure zone. The strain curves in each zone are called stable, critical and failure curves, respectively.

Cumulative equation of plastic deformation of graded crushed stone.
Through the linear analysis of plastic deformation, the plastic deformation of graded gravel increases with the increase in number of loadings. The equation should meet the boundary conditions of N = 0, ε = 0 and N = 1, ε = ε 1 assuming a cumulative plastic deformation equation for graded gravel. According to the boundary conditions, the cumulative equation of the plastic deformation of graded crushed stone is analysed and established.
where ε is the axial strain of graded crushed stone (%), ε 1 is the ultimate permanent strain of graded crushed stone (%), ξ is the cumulative coefficient of plastic deformation and N is the dynamic load action times. According to formula (6), the plastic deformation data of graded crushed stone are fitted, and the ultimate permanent strain and cumulative plastic deformation coefficient of graded crushed stone under three confining pressures (50, 100 and 150 kPa) are obtained. The results are given in Table 3. R 2 is the correlation coefficient.

Critical failure stress
Based on the above analysis, the upper limit of dynamic stress corresponding to the critical interval is the critical failure stress σ p . If the dynamic stress of the graded crushed stone base is  less than the critical failure stress, the plastic deformation of the graded crushed stone base will be effectively controlled. In this work, the plastic strain accumulation equation of graded crushed stone was fitted through the formation mechanism of plastic deformation and the dynamic triaxial test results of plastic deformation. The ultimate strain of graded gravel was regressed by SPSS software, and the critical failure stress was determined.

Plastic deformation failure criteria.
To determine the value of plastic deformation failure standard of graded crushed stone, the stress-strain relationship of graded crushed stone was studied, as shown in Fig 15.  Fig 15 shows that the stress of graded crushed stone increases with the increase in strain initially, then decreases with the increase in strain after reaching the peak value and finally tends to be stable. The failure strain corresponding to the peak value is mainly in the range of 4%-5%, which indicates that the antideformation performance of graded crushed stone decreases rapidly when the axial strain is 4%-5%. According to a previous study, an axial plastic strain of more than 4% in the pavement is unacceptable [26]. Therefore, 4% axial plastic deformation was used as the failure criterion ε p to determine the critical failure stress of graded crushed stone.

Determination of critical failure stress.
The schematic of the critical failure stress of graded crushed stone according to the limited plastic strain data of the regression graded crushed stone in Table 3 is shown in Fig 16. The failure strain line of graded broken stone is drawn in the figure, as shown by the dotted line in the figure.
It can be seen from Fig 16 that under the same confining pressure, the ultimate plastic strain increases with the increase in deviator stress, and the strain rate decreases gradually; under the same deviator stress, the ultimate plastic strain decreases with the increase in confining pressure. Confining pressure has a great influence on critical failure stress. The greater the confining pressure is, the greater the deviator stress required to reach the same strain. The critical failure stress and confining pressure are approximately linear. See Table 4 for the critical failure stress corresponding to different confining pressures. To sum up, to effectively control the plastic deformation of graded crushed stone under repeated dynamic loading, the deviator stress under different confining pressures σ d should be less than the corresponding critical failure stress σ p . The critical failure stress level corresponding to the critical failure stress S p is approximately 0.310-0.424.

Cumulative failure equation of plastic deformation
4.3.1 Cumulative failure stress condition of plastic deformation. The horizontal stress range of the graded crushed stone transition layer should be considered in the selection of confining pressure. According to ANSYS mechanical analysis, the horizontal stress range of the graded crushed stone transition layer is 48.9-70.8 kPa, and the horizontal stress range of the graded crushed stone base is 26.3-56.7 kPa. Therefore, in the dynamic triaxial simulation test of graded gravel, the range of confining pressure was taken as 20-100 kPa to cover the horizontal stress range of the graded gravel transition layer. We did not take the four levels 20, 50, 80 and 100 kPa.
It can be seen from Table 4 that to ensure failure, the cumulative failure stress level of plastic deformation used for calculation should be greater than the upper limit of critical failure stress level, i.e. S d > 0.424. Test efficiency and accuracy should also be considered in the determination of S d . If S is too small, the number of dynamic load actions required to reach the failure strain of graded crushed stone specimens is very large, and the simulation test time is too long. If S is too large, it is inconsistent with the actual stress condition of the graded gravel transition layer. In view of this, S d is proposed to be 0.5, 0.55, 0.6, 0.65 and 0.7.

Determination of the cumulative failure equation of plastic deformation.
Assuming that there is no cumulative failure limit of plastic deformation for graded crushed stone materials, the cumulative failure equation of plastic deformation should satisfy the following boundary conditions: where S d is the cumulative failure stress level of different plastic deformations and N d is the number of dynamic load actions required to reach the failure strain of graded crushed stone simulation specimens.
No plastic deformation cumulative failure equation strictly satisfies these two boundary conditions and can well fit all the test results. Therefore, these two boundary conditions are often relaxed to obtain a more reasonable cumulative failure equation of plastic deformation: where α and β are the undetermined regression coefficients. Eq (9) satisfies the boundary condition given by Eq (8)

Conclusions
1. In this work, the dynamic triaxial numerical test method of graded crushed stone based on particle flow theory and its application were studied. On the basis of DEM, a dynamic triaxial numerical test method for graded crushed stone was constructed, and the sensitivity of test conditions to numerical test results was studied. Then, the standard test conditions of the numerical test were proposed. Results show that the stability of the dynamic triaxial numerical test of graded crushed stone is better when the calculation time step is greater than 1E−4 s per step. The dynamic load has a disturbing effect on all parts of the aggregates of the virtual specimen. With increase in the number of dynamic loadings, the plastic deformation of graded crushed stone gradually accumulates, and the plastic deformation of the simulated specimen tends to be stable under 10,000 cyclic dynamic loadings. When the specimen's size is 10 cm (F) × 20 cm (h), the effect of its size can be ignored.
2. The dynamic triaxial numerical test method was used to study the plastic deformation law of graded crushed stone. The horizontal stress of the graded crushed stone transition layer was calculated using ANSYS. The peak deviator and critical failure stresses of graded crushed stone were determined using the numerical test method. Results show that the horizontal stress range of the graded crushed stone transition layer is 48.9-70.8 kPa, the horizontal stress range of the graded crushed stone base is 26.3-56.7 kPa, and the cumulative failure test stress level of plastic deformation should be greater than 0.424.
3. On the basis of the dynamic triaxial test method, the cumulative failure equation of the plastic deformation of graded crushed stone was established, and the coefficient of the cumulative failure equation was regressed. Results show that confining pressure has a great influence on critical failure stress. The greater the confining pressure is, the greater the deviator stress required to reach the same strain. The critical failure stress and confining pressure are approximately linear, and the plastic deformation of graded crushed stone is established. The cumulative failure equation is lgS d = lgα-βlgN d .