The Effects of the Particle Size Ratio on the Behaviors of Binary Granular Materials

Deze Yang and Xihua Chu

School of Civil Engineering,Wuhan University,Wuhan,430072,China

ABSTRACT The particle size ratio (PSR) is an important parameter for binary granular materials, which may affect the microstructure and macro behaviors of granular materials. However, the effect of particle ratio on granular assemblies with different arrangements is still unclear.To explore and further clarify the effect of PSR in different packing structures, three types of numerical samples with regular, layered, and random packing are designed.Numerical results show that PSR has significant effects on binary granular samples with regular packing.The larger the PSR,the stronger the strength,the larger the modulus,and the smaller the angle between the shear band and the load direction.And a theoretical solution of the peak stress ratio vs.PSR is obtained for regular packing,and the results by DEM are in good agreement with the theoretical solution.Under layered packing,PSR has little effect on peak stress ratio due to similar microstructure obtained with the changing of PSR.The modulus slightly increased with the increase of PSR.Under random packing with small grain content of 50%,PSR has little effect in the range of 0.5-0.9,but in a larger range,larger PSR leads to greater modulus.

KEYWORDS Discrete element method;binary granular materials;particle size ratio;packing

1 Introduction

A granular material generally refers to a collection of a large number of discrete solid particles larger than 1 μ m in size,such as sand and soil.The unique mechanical properties of granular materials have attracted the attention of a large number of researchers.In the past two decades,a considerable number of studies have been carried out to investigate and predict the mechanical behaviors of granular materials[1-6].It has been recognized that the macroscopic behavior of granular materials is controlled by the microstructure characteristics,such as particle packing[1],particle shape[2],and particle size[3,4].

Binary granular materials are a kind of common granular materials composed of particles of two different sizes. Generally, the size ratio of small particles to large particles (PSR) and the volume fraction of small particles (VS) are key factors affecting the properties of binary granular materials[7-13]. The application of binary granular materials with different packing methods is different.For random packing, binary granular materials can be a simplified model of a well-graded mixture[7]. Many studies have reported the mechanical properties of binary particle mixture with random packing based on Discrete Element Method(DEM)[14-26].Ogarko et al.[8]found that in the case of isotropic compression,any polydisperse mixtures can be replaced by equivalent bidisperse mixtures when the size distribution moments are matched.Wiącek[23]experimentally and numerically studied the effects of the PSR and VS on the geometrical properties of binary granular mixtures,and found a strong relationship between the coordination number and the VS of binary mixtures.Kumar et al.[9],Zhu et al. [24] and Ng et al. [25] found that the modulus is related to PSR and VS. Ueda et al. [7]explored the PSR and VS of the binary mixture to the shear resistance,and identified the threshold VS,which determines whether the mechanical contribution of large or small particles can be negligible.In addition,they found that the contribution of the large particles disappears when their average spacing with respect to the small particle size is around 2.Nie et al.[26]explored the effects of the PSR and VS on the internal stability of binary mixtures,and demonstrated that the internal stability of binary mixtures generally decreases with decreasing PSR and increases with increasing VS.

For regular packing,binary granular materials are called binary granular crystals or metamaterials. At present, the properties of wave or acoustics in binary granular crystals are mainly focused[12,13], which is widely used in building shock absorption and sound insulation. However, in some engineering applications, such as for building exterior walls, a certain bearing capacity of binary granular metamaterials is required.However,there are few reports on the static properties of binary granular metamaterials. Due to the need for contact between adjacent particles and the limitation of geometric space, the PSR of binary granular metamaterials varies within a small range. In this case, whether PSR is a key factor affecting the behavior of binary granular materials remains to be investigated.

In this study, the strength and deformation characteristics of binary granular materials with different packing structures in a quasi-static state are concerned. Which is of great significance for the design of metamaterials and the study of well-graded mixtures.The PSR ranges from 0.5 to 0.9 in binary granular materials according to the actual situation in binary granular crystals.Binary granular materials with different packing structures are focused on,which are regular packing,layered packing,and random packing.Their quasi-static mechanical properties are investigated based on DEM under biaxial loading in 2-dimensional.The effect of PSR on their strength and deformation characteristics at different confining stress are discussed.

2 Discrete Element Method and Numerical Samples

2.1 Discrete Element Method

Discrete element method(DEM)was initiated by Cundall et al.in 1979[27].The cores of discrete element method are particle motion equations and contact models between particles.

The motion equation of a particle A can be written according to Newton’s second law as:

The normal force is related to the overlap of two particles:

where δcis the normal overlap,R*, ν*andG*are the effective contact radius, Poisson’s ratio, and shear modulus, respectively. The effective contact radius, Poisson’s ratio, and shear modulus are computed as:

with

whereRis the radii of the contacting particle,E,Gand υ are the Young’modulus,the shear modulus and Poisson’s ratio respectively,the superscripts A and B denote the two particles in contact.

whereµis the friction coefficient.

Figure 1:Sketch of contact model between particles

2.2 Numerical Samples

In this section,Numerical samples of binary granular materials with three packing structures were established.The size of the samples is about 1.5 m×2 m.In all samples,the radius of large particles is R=0.01 m,and the radius of small particles increases from r=0.5 R to r=0.9 R.With the decrease in the regularity of the packing structure,three different packing arrangements are generated,which are regular packing,layered packing,and random packing,as shown in Fig.2.In the regular packing,we adopted the same packing method as Wang et al.[12],filled the small particles in the square array of large particles to generate regular packing structure,as shown in Fig.2a.When PSR is 1,the packing structure is hexagonal packing with single particle size.The VS for regular binary mixtures are 18.32%,25.97%,32.28%,38.35%and 44.01%for PSR of 0.5,0.6,0.7,0.8 and 0.9,respectively.In the layered packing,the sample is generated by stacking three large particle layers and two small particle layers,and each particle layer results from a hexagonal packing of particles, as shown in Fig.2b. And the layered samples with different PSR have the same VS of 40%.In the random packing,the samples are generated by random packing of large and small particles with VS of 50%,and samples with different PSR have the same initial porosity of 0.15,as shown in Fig.2c.

Figure 2: Three arrangement methods: (a) regular packing, (b) layered packing, and (c) random packing

The biaxial test is widely used in the study of granular materials, especially in DEM numerical simulation,as it is easier to visualize particle interactions in two dimensions,and the results are valid for 3-dimensional analysis. The biaxial test is regarded as a simplification of the triaxial test, so a series of biaxial tests for granular samples is performed,as shown in Fig.3.During biaxial tests based on DEM, firstly rigid plates were manipulated by a servo-controlled mechanism to reach the initial confining pressure and complete the consolidation process.And then the loading was applied by the displacement of the rigid plates at top and bottom,and the rigid plates on the sides keep the confining pressure constant.Loading ends when the target axial strain is reached.The axial strain(ε)is defined as the ratio of the vertical displacement of the rigid plates at top and bottom(Δl)to the height of the sample(l0)after consolidation.In this study,the loading speed of the rigid plates at top and bottom was controlled as v=0.5%×l0m/s.Besides,the same microscopic parameters as Wang et al.[12]were used for simulation,as shown in Table 1.

Figure 3:The biaxial compression loading modes

Table 1: Microscopic parameters for DEM simulation

3 Numerical Results and Discussion

The numerical results for each granular sample, including the axial stress-strain curve, axialvolumetric strain curve, and the distribution of the effective strain [28,29] were presented. For the definition of the effective strain,see Appendix A.

3.1 Regular Granular Samples

A series of DEM simulations were carried out for regular granular samples with different PSR at different confining stress.The deviatoric stress-axial strain curves of regular granular samples with different PSR under different confining stress levels are shown in Figs.4a-4d, and the relationship between peak stress and PSR is shown in Fig.5a. It can be seen that with the increase of PSR, the peak deviatoric stress and corresponding axial strain increase.After reaching the peak deviatoric stress,with the increase of axial strain,the deviatoric stress decreases rapidly.The initial modulus increases with increasing PSR,and the peak deviatoric stress and corresponding axial strain increase with the increase of confining pressure stress level.Notice that regular granular samples exhibit strong brittle behavior,a similar phenomenon can be found in the study of O’Sullivan et al.[30],which is related to the rapid destruction of regular granular samples.

For the biaxial tests,there is an analytical solution for the strength of the regular packing structure.Utilizing the condition for sliding at contact points and the symmetry of the packing structure, the analytical solution of the maximum stress ratio was given by Rowe[31]:

where σ1and σ2are axial stress and lateral stress respectively,tan φuis the friction coefficient,and α is the structural slip angle,as shown in Fig.6.For regular packing structure in Fig.2a,α can be expressed by PSR as

where γ is the PSR. According to Eq.(14), the maximum stress ratio of regular packing granular samples is

Figure 4: Deviatoric stress-axial strain curves of regular granular samples under different confining stress levels:(a)100 kPa,(b)200 kPa,(c)300 kPa and(d)400 kPa

It can be seen from Fig.5b,the DEM simulation results agree well with the theoretical solution.

Figure 5:Curve of the relationship between peak deviatoric stress and PSR of regular granular samples:(a)peak deviatoric stress under different confining pressures and(b)theoretical solution by Eq.(15)

Figure 6:Stress analysis of regular packing structure

Figs.7a-7d show the relationship between the volumetric strain and the axial strain of the regular granular samples with different PSR under different confining stress.It can be seen that the regular granular samples with different PSR show the phenomenon of first volumetric contraction to volumetric dilation.The peak volumetric contraction increases with the increase of PSR and confining stress level.At the same time,the higher confining stress level and PSR,the higher the axial strain of the regular granular samples at peak volumetric contraction.

Figure 7:Volumetric strain-axial strain curves of regular granular samples under different confining stress levels:(a)100 kPa,(b)200 kPa,(c)300 kPa and(d)400 kPa

Fig.8 shows the effective strain distribution of regular granular samples at confining stress level of 100 kPa after peak deviatoric stress.Both of them show that shear bands appear after peak deviatoric stress.In the regular packing structure shown in Fig.6,the angle between the sliding particles and the loading direction is β when the particles slide,which decreases with the increase of PSR.Therefore,with the increase of PSR,the angle between the main shear band and the loading direction decreases.The formation and shape of the shear band have certain randomness, and there is a bifurcation phenomenon.

Figure 8:Effective strain distribution of regular granular samples at confining stress level of 100 kPa after peak deviatoric stress:(a)r=0.5 R,(b)r=0.6 R,(c)r=0.7 R,(d)r=0.8 R,(e)r=0.9 R and(f)r=R

3.2 Layered Granular Samples

Layered granular samples whose PSR increases from 0.5 to 0.9 are considered,as shown in Fig.2b.A series of DEM simulations were carried out at different confining stress.The deviatoric stress-axial strain curves of layered granular samples with different PSR under different confining stress levels are shown in Figs.9a-9d,and the relationship between peak deviatoric stress and PSR is shown in Fig.10.It can be seen that the layered granular samples with different PSR show similar deviatoric stress-axial strain curves and peak deviatoric stress.It is because although the granular samples are layered packing with different PSR,each particle layer is regularly hexagonal packing.Similar microstructures result in similar strengths for particle slip failure. Note that the modulus of the layered granular samples increases slightly with increasing PSR,a plausible explanation is the effect of Hertz-Mindlin contact law.According to Eq.(4),the normal contact stiffness is

Figure 9: Deviatoric stress-axial strain curves of layered granular samples under different confining stress levels:(a)100 kPa,(b)200 kPa,(c)300 kPa and(d)400 kPa

Figure 10:The influence of PSR of layered granular samples on peak deviatoric stress under different confining pressures

It can be seen that the normal contact stiffness is related to particle radius,the larger the particle radius,the greater the normal contact stiffness.With the increase of PSR,the particle radius of small particle increases,and the normal contact stiffness between small particles increase,resulting in a slight increase in modulus. To demonstrate this explanation, we use the linear contact model for layered granular samples with the same contact stiffness kn=3e7 N/m and ks=1e7 N/m.The layered granular samples are loaded in the same way at the confining stress of 100 kPa.The result is shown in Fig.11,it can be seen that the layered granular samples with different PSR have the same modulus, which demonstrated the increase in the modulus of layered granular samples is related to the Hertz-Mindlin contact model.

Figure 11:Deviatoric stress-axial strain curves of layered granular samples with linear contact model at confining stress of 100 kPa

Fig.12 shows the relationship between the volumetric strain and the axial strain of layered granular samples at different confining stress.It can be observed that similar to the regular granular samples,the volumetric contraction occurs first,and with the increase of axial strain,the volumetric dilatation occurs. Similar volumetric strain-axial strain curves are observed for layered granular samples with different PSR.

Figure 12:Volumetric strain-axial strain curves of layered granular samples under different confining stress levels:(a)100 kPa,(b)200 kPa,(c)300 kPa and(d)400 kPa

Fig.13 shows the effective strain distribution of layered granular samples at confining stress level of 100 kPa at the axial strain of 2%.Both of them show that shear bands appear.The main shear bands run through several grain layers, and several shear bands occur at the layer boundaries. Although the PSR of layered granular samples is different,the angles between the shear band and the loading direction are similar.

Figure 13:Effective strain distribution of layered granular samples at confining stress level of 100 kPa at axial strain of 2%:(a)r=0.5 R,(b)r=0.6 R,(c)r=0.7 R,(d)r=0.8 R and(e)r=0.9 R

In summary,the PSR has little influence on the strength and deformation characteristics of layered granular samples.It may be because although the PSR of the layered granular samples is different,they have very similar microstructures,resulting in similar strength and deformation properties.

3.3 Random Granular Samples

Random granular samples with the same porosity whose PSR increases from 0.5 to 0.9 are considered, as shown in Fig.2c. And VS of random granular samples is 50%. A series of DEM simulations were carried out at different confining stress. The deviatoric stress-axial strain curves of random granular samples with different PSR under different confining stress levels are shown in Figs.14a-14d.It can be seen that the particle assemblies show strain softening.The deviatoric stress tends to stabilize with the increase of axial strain.

Figure 14:Deviatoric stress-axial strain curves of random granular samples under different confining stress levels:(a)100 kPa,(b)200 kPa,(c)300 kPa and(d)400 kPa

Fig.15 shows the volumetric strain-axial strain curves of random granular samples at different confining stress. Similar to the particle assemblies with the regular and layered arrangement, the volumetric contraction occurs first, and with the increase of axial strain, the volumetric dilatation occurs. Fig.16 shows the effective strain distribution of randomly arranged particle assemblies at confining stress level of 100 kPa at the axial strain of 15%. Compared with particle assemblies with the regular and layered arrangement,particle assemblies with the random arrangement also show that shear bands appear,but the locations of shear bands are random and not very concentrated.

Figure 15:Volumetric strain-axial strain curves of random granular samples under different confining stress levels:(a)100 kPa,(b)200 kPa,(c)300 kPa and(d)400 kPa

Figure 16:Effective strain distribution of random arranged particle assemblies at confining stress level of 100 kPa at axial strain of 15%:(a)r=0.5 R,(b)r=0.6 R,(c)r=0.7 R,(d)r=0.8 R and(e)r=0.9 R

It can be seen that within the range of 0.5-0.9, PSR has little influence on the mechanical and deformation properties of random granular samples.It may be related to the similar microstructure.As indicated by Zhu et al.[24]and Gong et al.[32],Young’s modulus and the shear modulus are strongly linked to the coordination number.In this study,the coordination numberZ,mechanical coordination numberZmand volume-weighted coordination numberZvare introduced[33-34],defined as follows:

whereCis total number of contacts;Nis total number of particles;N1andN0are numbers of particles with one contact and no contacts,respectively;Vis the volume of the particle assembly;Cpis contact number of the particlep, andVpis the volume of the particlep. The relationship between the PSR and the coordination numberZ, the mechanical coordination numberZmand the volume-weighted coordination numberZvat confining stress of 100 kPa is shown in Fig.17. Similar coordination numberZand mechanical coordination numberZmcan be observed, although PSR is different.The volume-weighted coordination numberZvdecreases slightly with increasing PSR.In general,the changes in PSR do not lead to significant changes in coordination properties,which may account for similar mechanical and deformation properties of random granular samples with different PSR.

Figure 17:The relationship between PSR and coordination number

Although the mechanical and deformation properties of random granular samples are similar within such a small range of PSR, the contact behavior is changed to some extent. For binary mixture with random packing, the contacts in a particle assembly can be classified as big particlebig particle (bb), big particle-small particle (bs), and small particle-small particle (ss) contacts. The partial coordination numbers which can reflect different contacts are defined as[24]:

whereCbb,CbsandCssare numbers ofbb,bsandsscontacts,respectively,andNbandNsare numbers of big particles and small particles, respectively. The relationship between the PSR and the partial coordination at confining stress of 100 kPa is shown in Fig.17.It can be seen that the increase of PSR leads to greaterZbbandZbsand smallerZss,which is consistent with the results of Zhu et al.[24]at VS of 50%.It indicates that as the PSR decreases,bbcontacts decrease,and easier to generatebsandsscontacts.

Within a limited range of 0.5-0.9, the microstructure changes caused by PSR are very limited.We generated random granular samples with a larger range of PSR with VS of 50% in the same way,in which PSR is 1/2,1/4,1/8,and 1/16,and performed the same loadings at confining stress of 100 kPa.If the radius of large particles is fixed,the PSR will cause a large computational cost when the PSR is relatively small, so we fix the radius of small particles to 0.008 m. The deviatoric stressaxial strain curves of random granular samples with different PSR at confining stress of 100 kPa are shown in Fig.18. It can be seen that the modulus of random granular samples increases with the decrease of PSR,which is consistent with the result obtained by Zhu et al.[24].They also discussed the contribution of contact type to modulus and found that thebs,andsscontacts accounted for the main contribution in random binary granular mixtures with VS of 50%under small PSR.And it can be well reflected in the force chain distribution,as shown in Fig.19.When PSR is small,strong force chains are mainly distributed inbscontacts, and the weak chain is distributed insscontacts which have numerous quantities.Fig.20 shows the relationship between PSR and the partial coordination numbers.It can be seen that with the increase of PSR,sscontacts decrease,bsandbbcontacts increase.It is due to the number of big particles increasing with increasing PSR.

Figure 18: Deviatoric stress-axial strain curves of random granular samples at confining stress of 100 kPa

Figure 19:Force chain distribution of(a)PSR=1/8,(b)PSR=1/16

Figure 20:The relationship between PSR and the partial coordination numbers

4 Conclusions

The effect of PSR on the strength and deformation of binary granular materials with different packing methods was studied.PFC2D was used to simulate the influence of PSR on basic mechanical properties such as shear strength,volumetric strain,and so on.Besides,the effective strain distribution and the fabric tensors of the granular materials are analyzed.The main conclusions are as follows:

1) The macroscopic mechanical properties of binary granular materials with regular packing are greatly affected by the PSR. It is mainly because the effect of PSR mainly comes from the changes in microstructure induced by PSR.The theoretical formula for the peak stress ratio and PSR is obtained, which is in good agreement with the DEM simulation. The strength and modulus of binary granular materials with regular packing increase with increasing PSR,which is mainly due to the change in the structural slip angle. Clear shear bands appear in binary granular materials with regular packing after peak deviatoric stress.With the increase of PSR,the angle between the main shear band and the loading direction decreases.

2) The PSR has little influence on the macroscopic mechanical properties of binary granular materials with layered packing, which is caused by similar microstructure. The modulus of binary granular materials with layered packing increases with increasing PSR,and we suspect that it is due to the Hertz-Mindlin contact law used in the DEM simulations,which results in greater contact stiffness with an increase of PSR.

3) For random binary granular materials with VS of 50%, PSR has little effect on mechanical and deformation properties within the range of 0.5-0.9.In the larger range of PSR,with the increase of PSR,the modulus of the random binary granular materials increases.Strong force chains are mainly distributed inbscontacts, and the weak chain is distributed insscontacts which have numerous quantities.

Funding Statement:This work was supported by the National Natural Science Foundation of China(Nos.12172263,11772237).

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

Appendix A.The Definition of the Effective Strain

To describe the change of relative position between particles, Li et al. [28] and Tang et al. [29]defined the nominal strain based on the change of particle position.Fig.A1 shows the position change of particle A and one of its neighboring particles B.

Figure A1:Position change of particle A relating to a neighboring particle[28]

The coordinates of particles A and B at timetnandtn+1areandrespectively in the global coordinate system.

Refer to local coordinate systemdeformation gradient fncan be used to describe the relative change of center position of particles A and B.

where

where α1and α2correspond to the angle between the local coordinate x-axis and the global coordinatex-axis attnandtn+1.Substitute Eqs.(A3a)and(A3b)into Eq.(A4).

where

Then the displacement derivative matrix can be defined as

where I refer to the identity matrix.Then strain between particles A and B can be given as

whereDijrefers to the component of matrix D. Then the nominal strain of particle A withnAneighboring particles can be defined as