Grain Size Effect on the Mechanical Behavior of Cohesionless Coarse-Grained Soils with the Discrete Element Method

Biaxial compression tests with the same specimen size and different maximum grain sizes were simulated for coarse-grained soils using the discrete elementmethod to study the influence of grain size on the mechanical properties and force chain.*emaximum grain sizes were 40, 20, 10, and 5mm, respectively. *e grading with self-similar fractal structure in mass is designed to ensure the same pore structure for soils. *e shear strength increased with the increase in maximum grain size. Evident increase in shear strength and significant size effect were observed when the ratio of the specimen diameter to maximum grain size was less than five. *e shear dilation of coarse-grained soils increases with the increase in maximum grain size. *e contact force distribution was uniform when maximum grain size was small but tends to be uneven with the increase in maximum grain size, thereby causing the increase in shear strength by stable strong force chains. *is finding demonstrates size effect on the mechanical properties and force chain of cohesionless coarse-grained soils under the biaxial compression condition.


Introduction
e size effect on cohesionless coarse-grained soils is an important topic in geotechnical engineering.Size effect, including sample and particle size effect, can be studied through tests and numerical simulations.A large number of test results for coarse-grained soils have been obtained.When the maximum grain size is large or the specimen size is small, the shear strength is high [1][2][3][4][5][6][7].However, a few factors in laboratory tests, such as particle shape, gradation, uniformity, and particle breakage of soil, affect the test results [8][9][10][11][12].
e grading with the self-similar fractal structure in mass was used in triaxial compression tests to ensure the same pore characteristics for specimens with different maximum grain sizes and avoid the influence of grading [7].However, the laboratory tests can only obtain a few macroscopic behaviors.With the increase in numerical simulation methods, the size effect of coarse-grained soils was simulated based on numerical tests, thereby controlling the error from the test material.Furthermore, the load transfer and fabric change in the deformation and failure process can be obtained, and the mechanism of size effect can be analyzed.erefore, analyzing the size effect on coarse-grained soils using numerical tests is significant.
e discrete element method (DEM), which is an important method in numerical simulation, has been widely used to investigate coarse-grained soils since it was proposed by Cundall and Strack [13].e size effect of coarse-grained soil is one of the main contents of DEM.Potyondy and Cundall [14] studied the influence of the ratio of specimen size to particle size on the elastic modulus and internal friction angle of the sample using particle flow code (PFC).e results showed that the mechanical parameters are later the same in the two-dimensional tests, whereas the parameters in the threedimensional tests are related to the size.Matthew and Katalin [15] used the DEM to carry out the biaxial simulation tests for four cohesionless particle materials with different specimen sizes.eir results showed that sample shear strength reduces with the increase in specimen size.A series of triaxial compression tests were carried out on coarse-grained soil through PFC. e results demonstrated that the change in the ratio of specimen size to particle size has no evident influence on the shear strength of coarse-grained soils when the ratio is greater than 40 [16].Rao and Chang [17] studied the influence of grain size on stress-strain-strength behaviors by using PFC 2D for soils with the maximum grain sizes of 60, 80, and 120 mm and the same ratios of specimen size to maximum particle size.
Scholars have attempted to explain the phenomenon of size effect in terms of its formation mechanism due to the complexity and uncertainty of the size effect of coarse-grained soils.Zhu et al. [12] believed that boundary constraints of the specimen can be attributed to the size effect.However, the specific action of the constraint on the formation of the size effect is not explained.Zhang and Hou [18] explained the size effect of granular materials due to the boundary constraint based on force chain theory.e contact force between the particles is evidently affected by the boundary when the specimen size is smaller than the characteristic length of the force chain.
e biaxial compression tests for cohesionless granular soils with different maximum grain sizes and self-similar fractal structures were performed based on the DEM to analyze the effect of maximum grain size on the mechanical properties and force chain of coarse-grained soils with the same specimen size.
e mechanism of size effect of cohesionless coarse-grained soils was preliminarily discussed by the concept of the force chain.

Introduction of Contact Stiffness Model.
e linear contact stiffness model of PFC 2D was adopted in this paper [19].In the contact plane, the contact force vector between two contact entities can be decomposed into normal and tangential components.e normal contact force component is the normal contact stiffness product and the overlap of particles and the normal unit vector of the contact surface, which can be written as follows [19]: where F n i is the normal contact force component, k n is the contact normal secant stiffness, U n is the relative contact displacement in the normal direction, and n i is the normal vector unit that defines the contact plane.
e shear force is computed in an incremental fashion.When the contact is formed, total shear force is initialized to zero.Each subsequent relative shear displacement increment leads to an increment of elastic shear force that is added to the current value, which can be written as follows [19]: where ΔF s i is the shear contact force component, k s is the shear stiffness, and ΔU s i is the relative contact displacement in the shear direction.
Based on the assumption that the stiffness of the two contacting entities acts in succession, the contact stiffness for the linear contact model is computed.e contact normal secant stiffness is given by the following equation: e contact shear tangent stiffness is calculated as follows: where superscripts [A] and [B] denote the two particles in contact.

Specimen.
Figure 1 shows the specimen of cohesionless coarse-grained soils with a diameter of 101 mm and a height of 200 mm.Among which, #1 and #3 walls were the lower and upper loading plates, respectively, and #2 and #4 were flexible boundaries, which can simulate the rubber membrane in triaxial tests.e test soils with self-similar fractal structure in mass as used in the biaxial compression test were utilized to ensure the same pore characteristics for specimens with different maximum grain sizes.e grading was achieved as follows [20]: where M(d < d i ) is the mass of grains with the grain size d less than a certain grain size d i , M t is the gross mass of grains, d max is the maximum grain size, and D is the fractal dimension.Four different scaling methods were adopted to scaled rockfill materials.e results showed that the range of the  Advances in Civil Engineering fractal dimension value is between 1.463 and 1.783 [21].erefore, D was set to 1.5 in the numerical tests.A total of four specimens with the similar grading were used in the numerical tests.e maximum grain size d max was 5, 10, 20, and 40 mm, and the minimum grain size d min was 0.25, 0.5, 1, and 2 mm, respectively.e particles of size less than d min were dealt with the equivalent replacement method.e average grain size d 50 was 3.1, 6.3, 12.6, and 25.2 mm, respectively.
e coefficient of uniformity C u and the coefficient of curvature of the specimens are 3.3 and 1.30, respectively.e graining curves are shown in Figure 2.

Parameters.
us far, the parameters used in DEM simulation cannot be directly obtained by laboratory tests because no definite relation exists between the mesostructure   Advances in Civil Engineering and macroparameters.e trial and error method is used to ensure that the numerical test results are in agreement with the laboratory test results.After repeated trial and error, the mesostructure parameters of biaxial numerical tests were obtained.e normal and tangential contact stiffness of the particles was 4 × 10 7 N/m, the normal and tangential contact stiffness of the wall was 2 × 10 8 and 0 N/m, respectively, the friction coefficient between particles was 2.4, the particle density was 2600 kg/m 3 , the porosity of specimens was 0.15, and the loading rate was 0.3 m/s.

Size Effects on Strength.
e deviator stress q-axial strain ε a curves with different maximum grain sizes under the confining pressure of 300 kPa are shown in Figure 3. e curves under other confining pressures possess similar behaviors.e specimen was found to exhibit a strain softening behavior.e peak deviator stress is considered to be the failure deviator stress q f .e q f − D/d max curves are shown in Figure 4. e deviator stress q f increases with the decrease in D/d max , and the increasing tendencies of q f are evident when D/d max ≤ 5. is result is consistent with the laboratory test under the same conditions [7].

Size Effects on Deformation.
e volumetric strain ε v -axial strain ε a curves with different maximum grain sizes under a confining pressure of 300 kPa are shown in Figure 5. e curves under other confining pressures have similar behaviors.
e volume of the specimen contracts initially and then dilates and eventually reaches a steady state when the axial strain is up to 15%.e behaviors between the volumetric strain corresponding to the axial strain of 15% and D/d max are shown in Figure 6. e shear dilatancy of cohesionless coarse-grained soils is evident with the decrease in D/d max .

Size Effects on Force Chain.
e mechanical properties of granular materials are determined by the force chain network formed by the mutual contact of particles [22][23][24].
e size effect mechanism of coarse-grained soils can be explained by the distribution and transfer law of contact force.Figure 7 is the force chain distribution of coarsegrained soils with different maximum particle sizes under the peak stress.As shown in Figure 7, "black" is denoted as the contact force which is larger than the average value, and its width is proportional to the contact force.
Figure 7 shows that the contact force is transferred mainly through the strong force chain.
e strong force chain is mainly distributed along the direction of the major principal stress, and the weak force chain is distributed around the strong chain and acts as support.Both chains form a force chain network structure.With the increase in confining pressure, the lateral strong force chains gradually develop to support the axial strong force chains, and both chains form a stable force chain network.e macroscopic behavior indicates that when the confining pressure is high, the shear strength is also high, thus indirectly proving the rationality of the numerical simulation.
e breakage of the contact points causes the entire force chain to disintegrate and recombine.Maintaining stability is difficult when more particles exist on the force chain.In general, maintaining the stability of the force chain is difficult when its particle number exceeds five [25,26].When the maximum particle size is small, more particles exist in the sample, and the contact force is uniformly distributed; thus, forming the strong force chain is difficult.erefore, the sliding and dislocation between particles easily occur, and the macroscopic behavior of particles shows that the soils have low shear strength.e increase in the maximum particle size results in evident sample boundary restraint.e contact force becomes uniform, and the stable strong chain becomes long due to the decrease in the particle number.Extending the strong force chains from one side of the sample to the other is easy, thereby forming a stable structure that is failure resistant.Macroscopically, the sample shows high shear strength and evident size effect.

Conclusion
e DEM simulation of biaxial compression tests for coarsegrained soils with different maximum grain sizes and selfsimilar fractal structures in mass was carried out.e effect of grain size on the mechanical properties and force chain of coarse-grained soils was also discussed.e main conclusions are as follows.
(1) e shear strength of coarse-grained soils increases with the increase in the maximum particle size.e shear strength of soils evidently increases when the ratio of the sample diameter to the maximum particle size is less than five.(2) When the maximum particle size increases, the shear dilatancy of cohesionless coarse-grained soil also increases.(3) When the maximum particle size is small, many particles exist in the sample, and the contact force distribution is uniform; thus, forming a strong force chain is difficult.When the maximum particle size increases, the particle number decreases, the contact force distribution of particles tends to be uneven, and the stable strong force chain increases the shear strength, thereby resulting in size effect.