Next Article in Journal
Diamond Alpha Hilbert-Type Inequalities on Time Scales
Next Article in Special Issue
The Fractal Characteristics of Soft Soil under Cyclic Loading Based on SEM
Previous Article in Journal
Optimal H1-Norm Estimation of Nonconforming FEM for Time-Fractional Diffusion Equation on Anisotropic Meshes
Previous Article in Special Issue
Effects of Relative Density and Grading on the Particle Breakage and Fractal Dimension of Granular Materials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrete Element Modelling of Fractal Behavior of Particle Size Distribution and Breakage of Ballast under Monotonic Loading

1
School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan 430070, China
2
Faculty of Civil and Environmental Engineering, Ruhr University Bochum, 44801 Bochum, Germany
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(7), 382; https://doi.org/10.3390/fractalfract6070382
Submission received: 25 May 2022 / Revised: 1 July 2022 / Accepted: 3 July 2022 / Published: 6 July 2022
(This article belongs to the Special Issue Fractal and Fractional in Geomaterials)

Abstract

:
Breakage of particles has a great influence on the particle size distribution (PSD) and the associated mechanical behavior of ballast under train loads. A discrete element method (DEM) simulation of triaxial testing under monotonic loading was carried out using FRM (fragment replacement method) breakable particles as ballast and a flexible shell model as membrane. The coupled model was validated by comparing the load-deformation responses with those measured in previous experiments and was then used to analyze the contact orientations and the distribution of particle breakage from a micromechanical perspective. The simulation results show that higher confining pressure and larger axial strain may increase the grain breakage (Bg) and the fractal dimension (D) of ballast. It was observed that most breakage was first-generation breakage, and that the proportions of the second- to fifth-generation breakage decreased successively. Moreover, as the axial strain or confining pressure increased, the percentage of small particle fragments increased in correspondence with the PSD curves that remained concave upwards, as the fractal dimension D of PSD increased. In addition, the evolution of D exhibited a linear correlation with grain breakage Bg. Contrarily, a quadratic curve relation between D and volumetric strain was exhibited under different axial strain stages. Therefore, D has the potential to be a key indicator to evaluate the degree of ballast crushing and PSD degradation, which may contribute to better decision making concerning track bed maintenance.

1. Introduction

Ballast is a granular material used for construction of railway beds. It can be characterized as having an extremely irregular particle shape with high porosity and strength. Ballast helps to transmit and distribute axle loads from the sleepers to the sub-ballast layer [1]. The dynamic stress of a passing train degrades and fouls the ballast layer which usually experiences continuous breakage, directly resulting in track settlement, densification, and the need for frequent maintenance. Ballast breakage may lead to a sharp increase in the number of fragments, weakening interlocking and impacting mechanical performances [2]. In practice, ballast needs to be maintained and replaced when its quality deteriorates to a certain extent.
Recently, fractal theory has been widely used in describing and studying the particle size distribution (PSD) of granular material [3]. Compared with traditional methods based on Euclidean geometry, fractal theory describes PSD well because it can characterize the local shape, size, and structure of particles, and any similarities among them [4]. It is beneficial to analyze the pore structure of materials by combining more than one method, considering such factors as the relationship between fractal characteristics, or the physical or chemical attributes of materials [5]. Early studies revealed that fractal theory could be used to study the fractal characteristics of soil [4,5,6], sand [7], rockfill [8,9,10] and other material [11,12]. Compared with natural soil or quartz sand of smaller size and wider grading distribution, layers of ballast are usually composed of medium-to-coarse gravel-sized particles (10–60 mm) whose uniformity coefficients vary from 1.5 to 3.0 according to the ballast specifications [13] which can be regarded as uniformly sized granular material. As more and more small fragments were produced by ballast breakage under train loads, the PSD curve of ballast also changed, with a fractal dimension which was directly proportional to the fragment content. Therefore, the fractal theory could be used to study PSD evolution and ballast breakage. However, studies on the effect of the fractal distribution of particle size on the characteristics of ballast are still not enough in themselves.
The discrete element method (DEM) has been extensively used in railway engineering and also pavement engineering [14,15,16]. In order to analyze the evolution of fractal PSD using the fragment replacement method (FRM), the particle is represented by a group of smaller ones without bonding as implemented in [16,17,18,19,20,21]. McDowell et al. [16] adopted octahedral shear stress as the fracture criterion and replaced broken spheres with smaller ones with overlap in their successful investigation of the evolution of fractal PSD for silica sand. Ciantia et al. [17] used a multigenerational DEM approach to simulate soil breakage in a one-dimensional compression test, and obtained a limit fractal dimension of 2.47. Zhou et al. [18] found that both particle size and particle contact force were distributed in a fractal manner by one-dimensional compression test simulation, which is consistent with the findings of Vallejo and Lobo-Guerrero [19]. De Bono et al. [20] conducted one-dimensional normal compression test simulation with different fracture criteria, and pointed out that both octahedral shear stress and the largest normal contact stress can lead to the correct normal compression lines and fractal PSDs. Li et al. [21] quantified the evolution of microstructures of carbonate sands under mechanical loading and found that the fractal dimension of crack networks increased with external loading due to crack branching via cleavage.
Although the FRMs used in the above studies can accurately simulate particle breakage based on DEM, most of the parent breakable particles involved were spherical and so could hardly simulate the interlocking between particles [16,17,18,19,20,21]. However, particle shape plays a key role in the fractal behavior of ballast with irregular shapes [22,23]. Therefore, in order to reveal the fractal characteristics of ballast, a DEM–FDM (finite differential method) coupled simulation of triaxial testing under monotonic loading using an irregular-shaped breakable ballast model is carried out in this study. The effects of axial strain and confining pressure on ballast breakage and the fractal behavior of ballast PSD are investigated. The relationships between fractal dimensions, grain breakage and the stress states of ballast are then established.

2. Ballast Breakage Model Using FRM

Firstly, the DEM study using FRM to simulate particle breakage should include a reasonable fracture criterion and a fragment replacement mode. They should be defined to check (1) whether a particle should break and (2) how any individual particle does break.
It is widely accepted that the failure of spherical particles under compression is in fact a tensile failure. Jaeger [24] suggested that the ‘tensile strength’ of rock grains can be measured by diametral compression between flat platens. For a grain of diameter d under a diametral force F, a characteristic tensile stress σ induced within it can be defined as
σ = F/d2
For railway ballast with typical sizes ranging from 10 to 60 mm, Lim et al. [25] conducted indirect diametral compression and measured the tensile strengths of ballast between flat platens pointing out that Weibull strength σ0 is a function of the particle size d. This is usually expressed in the form
σ0d b
where σ0 is the value of σ such that 37% of the blocks survive the test intact; while b is a constant and a function of the statistical distribution of flaws in the material, describing the size-hardening law. Typical values of b range from −0.69 to −0.40 for different types of ballasts. Following the single-particle crushing test results [25], the characteristic tensile strength σ0 of ballast size d0 (19.7 mm) is 14.7 MPa, and the value of b is −0.56, which is used to attribute random strengths to new fragments and other initial particles.
In order to determine the correct behavior of ballast under normal compression, De Bono and McDowell [20] investigated various possibilities for the fracture criterion, including octahedral shear stress, major principal stress within a particle, mean stress, and stress calculated from the maximum contact force on a particle. The criterion based on the maximum contact force (Fmax) was shown to better match the experiments and simulations, such that
σFm = Fmax/d2
According to Raymond et al. [1], the process of ballast degradation due to loading can occur in three ways: by the grinding of small-scale asperities, by breakage of angular projections, and by the splitting of particles into approximately equal parts. In addition, an irregular clump produces more contact compared to a spherical particle, leading to a lower maximum contact force. Hence, these particles will be replaced by new fragments when the Fmax stress exceeds than half of the characteristic strength.
From the perspective of computational efficiency, in this study, each particle was allowed to split into several smaller particles, maintaining conservation of mass when the value of σFm was greater than or equal to its Weibull strength. Previous studies have used spheres or disks to simulate particle breakage within the bonded parent sphere, where new sphere fragments to be contained overlap sufficiently [16,18,20] or tangentially to each other [17,19]. As shown in Figure 1a, new sphere fragments overlap sufficiently to be contained within the bounding parent sphere, with the axes joining their centers aligned along the direction of the minor principal stress. This process produces undesirable local pressure spikes during breakage and continues until internal forces are completely released and the system reaches equilibrium [16].
It is therefore acknowledged that it is impossible to simulate a perfectly realistic ballast fracture mechanism using disks or spheres [16]; therefore, eight-pebble clumps of irregular shape were adopted for this study. The fragment replacement mode is depicted in Figure 1b, which shows that irregular particles were replaced by smaller fragments within the bounding parent particles without overlaps. During the first-generation breakage, the eight-pebble clump might split into a four-pebble clump and four single spheres. All the fragments could continue to break and generate smaller particles. Except for single spheres, most of the particles could break and generate fragments without overlap, which means internal forces could be releases and the system could reach equilibrium immediately. Hence this fragment replacement mode can greatly improve calculation efficiency, and is therefore used in this study.

3. Simulations of Monotonic Triaxial Test

Based on the experimental monotonic triaxial tests conducted by Indraratna et al. [26], DEM–FDM coupled simulations of monotonic triaxial tests were carried out under different confining pressures ranging from 60 kPa to 240 kPa. Figure 2 shows the DEM–FDM model of a triaxial test with dimensions of 300 mm diameter and 600 mm height. It consists of 971 eight-pebble clumps as breakable ballast using FRM in PFC3D (particle flow code in three dimensions) with a flexible boundary as membrane using shell model in FLAC3D (fast Lagrangian analysis of continua in three dimensions). The top and bottom plates using disk walls have the same dimensions as the triaxial test, to maintain a consistent area of contact with the ballast sample. The eight-pebble clump has an irregular shape with an aspect ratio of 1.83, as expressed by the ratio of the particle length to the particle width. In accordance with the PSD of ballast in triaxial tests, the PSD used in these simulations is shown in Table 1.
During sample generation, the particles generated at random orientation by the Ball Distribute command. Firstly, the sample was divided into five equal layers, each layer was generated in turn and allowed to fall under gravity, then compacted as performed in the lab. Once the sample attained equilibrium, a flexible membrane model was generated to replace the cylinder wall. The membrane model was linked with the top and bottom walls by each node. The shell model as a flexible membrane consisted of 1440 facet units which were assumed to be triangles of uniform thickness lying between three nodal points. Each shell structural element was defined by its geometric and material properties. The parameters of FRM ballast particles and shell elements are listed in Table 2.
In order to reduce the influence of the initial arrangement of particles on the simulation results, the same sample was used for different loading conditions. Once the assembly generation procedure was completed, artificial isotropic compression was applied to the plates and membrane to achieve the required stress state. The assigned uniform pressure, ranging from 60 to 240 kPa, was applied to all shell elements of the membrane model, so that the required confining pressure was obtained. Monotonic loading was applied to both the top and bottom disk walls with a constant speed of 0.02 m/s, and the monotonic load test continued until the axial strain reached 20%. During the loading process, data relating to porosity, the stress–strain curve, volumetric strain and particle breakage were all recorded.

4. Results and Discussion

4.1. Stress–Strain Response

Figure 3 shows the simulations and experimental results, in terms of deviator stress and of volumetric strain against axial strain under different confining pressures. Comparing the experimental results with the simulation results, a trend of stress–strain relationship is seen to be well simulated. With the development of the axial strain, the deviator stress gradually increases and attenuates after reaching a relatively stable value, showing the strain softening behavior and the maximum deviator stress increases with the increase in confining pressure. However, the shear strength of the simulated assembly is overestimated in comparison with the experimental one. This could be caused by the irregular shape of particles with an aspect ratio of 1.83, which were recognized as flaky particles. Similarly, after having conducted a set of triaxial ballast tests to investigate the effects of ballast shape on ballast performance, Roner [27] found that randomly placed flaky material exhibited higher deviator stress and higher angles of internal friction than non-flaky material at the same void ratio. This offers a possible explanation as to why the deviator stress was higher than expected. In addition, an orientation parallel to the failure plane will cause a substantial strength reduction when a significant proportion of the particles are flaky.
However, the purpose of this study is not to resolve the particle shape problem, but rather to adopt the best approach to investigate the fragment replacement mode and particle crushing under shear. Therefore, this study applies simplified and irregular shaped particles, and focus on the fractal fragmentation and gradation evolution of particles.
Figure 3b shows that the sample exhibits shear shrinkage in the initial stage, but with the increase in axial deformation, an obvious shear expansion is observed, especially under low confining pressure. As the confining pressure decreases from 240 kPa to 60 kPa, the dilatancy effect becomes more obvious. Volumetric strain agrees reasonably well with the experimental results under a high confining pressure of 240 kPa where breakage is more pronounced. This is probably attributable to the generation of smaller size fragments from the split particles. Using more and smaller fragments to form the origin clumps may provide more realistic volumetric strains at lower confining pressures. Therefore, further investigation of the effects of fragment replacement mode is needed.
Figure 4 shows the effect of different confining pressures on the peak friction angle φp of the ballast sample, and both the simulation and experimental results show that the peak friction angle decreases with the increase in confining pressure. Moreover, the simulation using uniform flat shape particles shows a higher angle of internal friction than the experimental results. This is due to the fact that the FRM particle only undergoes particle splitting at higher stress thresholds and is unable to undergo corner breakage and surface abrasion to release internal forces at the initial strain stage.
It can be seen from Figure 5a that the deformation of the sample is allowed to be uneven, and the expansion of the middle part is the most obvious after using the flexible membrane model. This verifies that the flexible membrane model is better than the servo-controlled cylinder walls for simulating the triaxial test. In order to study the development of porosity during loading, three measurement spheres with radii of 85, 120, and 85 mm were arranged in the sample, as shown in Figure 5b. Figure 5b shows the evolution of porosity with axial strain in the upper, middle and lower regions under 240 kPa confining pressure. As the axial strain increases to 12%, the porosity of each part exhibits the characteristics of first decreasing and then increasing. Subsequently, the porosities in the upper and lower parts are stable, but the porosity in the middle region continues to increase.
The evolution process of the particle displacement vector of the sample under the confining pressure of 240 kPa is shown in Figure 6. During the loading process, the particles within the top and bottom half move toward the horizontal plane in the middle. When axial strain reaches 4%, due to mutual resistance, the intermediate particles start to move outward horizontally to produce a shear sliding surface, and the angle between the sliding surface and the horizontal direction gradually increases.
Under a confining pressure of 240 kPa, the polar coordinate distribution of contact numbers and average contact forces within the sample is plotted as shown in Figure 7. With the increase in axial strain, the average contact force in the sample first increases and then decreases between 12% and 20% axial strain, which is similar to the substantial strength reduction in the deviator stress after peaking shown in Figure 3b. Most of the internal contacts of the sample are distributed in the range of ±30° from the vertical plane in the initial stage; the contact number in the vertical direction then decreases with the loading, while the contact number in the horizontal direction increases. At the end of loading, the average contact force and contact number in the vertical direction are almost 3.36 and 2.49 times larger than those in the horizontal direction.

4.2. Particle Breakage

As shown in Figure 8, the number of fragments could be ignored at the initial compaction stage (4% axial strain) before increasing linearly above 8% axial strain. This indicates that the sample becomes denser in the compaction stage due to particle rearrangement and a reduction in the number of pores. When the porosity reaches the minimum value, the particle interlocking phenomenon is further exacerbated and a large number of particles are broken off by shearing. In addition, more particle breakage occurs under higher confining pressure.
In accordance with experiment [26], the differences in the percentage by weight of each particle size fraction before and after simulation (∆Wk) are plotted against the aperture of the lower sieve corresponding to that fraction, as shown in Figure 9. The equivalent diameters and weights of each particle size are calculated by tracing all particles and extracting the volumes in simulation. For aperture sizes between 18.9 mm and 29 mm, the percentage change in particle size decreased and formed a main trough. Meanwhile, For aperture sizes between 18.9 mm and 53 mm, the percentage change of particle size increased and formed a main peak.
A measure of grain breakage (Bg) equal to the sum of positive values of ∆Wk is now plotted and expressed as a percentage, as introduced by Marsal [28] for granular materials. The Bg under 240 kPa confining pressure increases with increasing axial strain. In particular, the largest increment of Bg can be found the stage from 8% to 12% axial strain, which also corresponds to the peak in the growth rate of crushing, as shown in Figure 8. Moreover, with a large number of particles broken in this phase, the peak deviator stress is gradually deceased when the axial strain reaches 8% as shown in Figure 3a.
Figure 10b shows levels of grain breakage Bg after loading under different confining pressures. It can be observed that more particles are broken off with increasing confining pressure in both the experimental test and the simulations, which proves the acceptability of the fracture criterion and fragment replacement mode used in the simulation study. It should be noted that the breakage in simulation is underestimated compared with the experimental test. This could be explained by the fact that the simple splitting model does not take particle abrasion and small corner breakage into consideration.
For increasing levels of axial strain, the quantities and positions of fragments under 240 kPa confining pressure are plotted in Figure 11. The generations of breakage are defined from first generation to fifth generation as shown in Figure 1. At the initial stage of 8% axial strain, almost all of the fragments belong to the first generation, and the lateral bulging of the sample is not obvious. When εa = 12%, the fragments coming from the second and third generations begin to generate at locations near the top and bottom plates. As εa increases to 16%, and to 20%, new fragments coming from the fourth and fifth generations, respectively, start to appear in the sample. Apart from the area near the top and bottom plates, most of these fragments (second–third generation) are distributed along the shear surface as shown in Figure 11e, which is consistent with the distribution of particle displacement in Figure 6e. In summary, 160 of 971 eight-ball clumps were split and a total number of 1341 fragments were generated. Among these, 1237 fragments were of the first generation, accounting for 92.2% of all breakage and located all over the sample, while 104 fragments of the second to fourth generation accounted for 7.8% of all breakage and are more concentrated in the upper half of the sample. Therefore, the first-generation breakage is dominant, and the proportions of breakage decrease from the second to the fifth generations, successively

4.3. Gradation Evolution and Fractal Distribution

Turcotte et al. [3] demonstrated that any initial distribution of soil particles would tend towards a self-similar fractal distribution with increasing particle fragmentation. Railway track bed typically consists of ballast of size 10–60 mm in different ballast grading specifications, generally larger than soil or calcareous sand. Considering the ballast breakage results in the previous section, the fractal evolution of the PSD against the axial strain and confining pressure were also studied.
Figure 12a shows the evolution of PSD under monotonic loading and a confining pressure of 240 kPa. It can be seen that, as the axial strain increases, the percentage of small particle fragments ranging from 10 mm to 35 mm increases, with the grading curves remaining concave upwards. This can be attributed to the supplement of aggregates induced by particle breakage from the upper size ranges. Figure 12b shows the comparison of PSDs under different confining pressures, ranging from 60 kPa to 240 kPa. It can be seen that all the sample PSDs after test are shifted left where smaller aggregates were generated, when compared with the initial grading. With increases in confining pressure, the sample grading shifts more towards the left, especially under the highest confining pressure of 240 kPa.
According to the study by Tyler et al. [11], the fractal grading can be described as
Fu(d) = M(δ < di)/MT = (di/dmax) 3−D
where Fu(d) is the percentage finer, dmax is the maximum diameter for the given grain assembly, MT is the total mass of particles, δ is the particle size, and M(δ < di) is the cumulative mass of the particles whose grain size is smaller than di. The fractal dimension D can be determined by:
lg[M(δ < di)/MT] = (3 − D)lg(di) + (D − 3) lg(dmax)
The PSD curve presents the correlations between the percentages passing and normalized particle diameters which gradually close to a linear line with the log–log scales. Figure 13a shows that the fitting lines shift upward, and the grading curve gradually becomes linear as the axial strain increases. Meanwhile, the linear portion of the curve from which the fractal dimension can be obtained increases in length, suggesting a more reliable value. Moreover, the slope becomes steeper (i.e., the fractal dimension D increases) with increasing confining pressure, with the slope appearing to become constant as shown in Figure 13b. The size ratio of the largest and smallest particle is much smaller than that of soil or quartz sand, with values of Cu ranging from 1.5 to 3.0, which can be regarded as granular material of uniform size. This is quite different from granular materials in natural state because the grading ballast used in railway is qualified by the ballast standard, leading to a much smaller initial value of D (0.2–0.3). Predictably, a continued compression will cause a further crushing and increasing of fine particles, leading to a further fractal.
With the faster and heavier of the passing train, the ballast layer are subjetced to settle down and crushing under the long-term cyclic loading, and the fractal dimension gradually increases. Ballast will need to be maintained or replaced once the fractal dimension reaches a constant value. Further experimental research is required to predict the critical value of the fractal dimension which indicates when ballast needs to be maintained.
Figure 14 shows the evolution of fractal dimension against grain breakage under different confining pressures. It shows that the fractal dimension of the sample increases with rising confining pressure and also axial strain. The relationship between D and Bg is almost linear. Therefore, the fractal dimension can be well used to evaluate the degree of crushing and deterioration of ballast.
The relationship between D and the volumetric strain εv for these four samples under different confining pressures can be plotted in Figure 15. At the initial stage (4% axial strain), the samples are compressed with an almost constant value of D; a linear relationship between volumetric strain and D can then be observed above 8% axial strain under different confining pressures, which corresponds with the maximum sample dilation angle as shown in Figure 3b. A quadratic curve relation between the volumetric strain and fractal dimension can be obtained for different levels of axial strain, as shown in Figure 15. At 20% axial strain, the values of a, b and c are 1.48 × 10−3, 5.25 × 10−4 and 0.26, respectively.

5. Conclusions

A DEM model of irregular-shaped ballast particles using the fractal replacement method (FRM) and a FDM flexible membrane model using shell elements are presented in this paper. A series of DEM–FDM coupled simulations of triaxial testing under monotonic loading and with different confining pressures were performed to study the ballast breakage and fractal behavior of PSD. We draw the following main conclusions:
(1)
The proposed FRM model can simulate the ballast breakage and capture the essential features of the stress-dilatancy in the monotonic triaxial test. During the loading process, the particles within the top and bottom half moved toward the horizontal plane in the middle. With axial strains above 4%, due to mutual resistance, the intermediate particles started to move outward to produce a shear sliding surface, where the angle between the sliding surface and the horizontal direction gradually increased.
(2)
More particle breakage occurred under higher confining pressure or larger axial strain and most such breakage was in the second generation or earlier. The number of fragments could be ignored at the initial compaction stage and the majority of particle breakage occurred with axial strain of 8% or above. Compared with the experimental results, the simulation underestimated the degree of particle breakage. This could be explained by that the simple splitting model of FRM did not take particle abrasion and small corner breakage into consideration.
(3)
With increasing axial strain or confining pressure, the percentage of small particle fragments increased in correspondence with the PSD curves which remained concave upwards, leading to an increasing fractal dimension. However, a different situation applies to granular materials in natural state because the grading ballast used in railway is usually qualified by the ballast specification, leading to much smaller values of D.
(4)
The evolution of fractal dimension against grain breakage showed a linear relation, Conversely, a quadratic curve relation was obtained between the fractal dimension and volumetric strain under different axial strain stages. Therefore, D has the potential to be a key indicator to evaluate the degree of ballast crushing and PSD degradation, which can contribute to better decision making concerning railway track bed maintenance in practice.

Author Contributions

Conceptualization, C.C. and Y.S.; methodology, C.C. and Z.W.; software, C.C. and X.Z.; validation, C.C., X.Z. and Y.S.; formal analysis, Y.S. and R.R.; investigation, L.Z. and R.R.; resources, L.Z. and Z.W.; data curation, Y.S.; writing—original draft preparation, C.C. and X.Z.; writing—review and editing, Z.W. and Y.S.; visualization, R.R.; supervision, C.C. and Z.W.; project administration, Z.W.; funding acquisition, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 51708438) and the Self-determined and Innovative Research Funds of Wuhan University of Technology (No. 2022CG030).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Raymond, G.P.; Diyaljee, V.A. Railroad ballast load ranking classification. J. Geotech. Eng. Div. 1979, 105, 1133–1153. [Google Scholar] [CrossRef]
  2. Selig, E.T.; Waters, J.M. Track Geotechnology and Substructure Management; Thomas Telford: London, UK, 1994; pp. 16–26. [Google Scholar] [CrossRef]
  3. Turcotte, D.L. Fractals and fragmentation. J. Geophys. Res. 1986, 91, 1921–1926. [Google Scholar] [CrossRef]
  4. Zhao, W.J.; Cui, Z.; Ma, H. Fractal features of soil particle-size distributions and their relationships with soil properties in gravel mulched fields. Arab. J. Geosci. 2017, 10, 211–218. [Google Scholar] [CrossRef]
  5. Tao, G.; Wang, Q.; Chen, Q.; Nimbalkar, S.; Peng, Y.; Dong, H. Simple graphical prediction of relative permeability of unsaturated soils under deformations. Fractal Fract. 2021, 5, 153. [Google Scholar] [CrossRef]
  6. Bai, Y.; Qin, Y.; Lu, X.; Zhang, J.; Chen, G.; Li, X. Fractal dimension of particle-size distribution and their relationships with alkalinity properties of soils in the western Songnen Plain, China. Sci. Rep.-UK 2020, 10, 20603. [Google Scholar] [CrossRef]
  7. He, S.H.; Ding, Z.; Hu, H.B.; Gao, M. Effect of grain size on microscopic pore structure and fractal characteristics of carbonate-based sand and silicate-based sand. Fractal Fract. 2021, 5, 152. [Google Scholar] [CrossRef]
  8. Sui, L.; Yu, J.; Cang, D.; Miao, W.; Wang, H.; Zhang, J.; Yin, S.; Chang, K. The fractal description model of rock fracture networks characterization. Chaos Soliton Fract. 2019, 129, 71–76. [Google Scholar] [CrossRef]
  9. Xiao, Y.; Meng, M.; Daouadji, A.; Chen, Q.; Wu, Z.; Jiang, X. Effects of particle size on crushing and deformation behaviors of rockfill materials. Geosci. Front. 2020, 11, 375–388. [Google Scholar] [CrossRef]
  10. Zhang, K.; Liu, X.; Liu, W.; Zhang, S. Influence of weak inclusions on the fracturing and fractal behavior of a jointed rock mass containing an opening: Experimental and numerical studies. Comput. Geotech. 2021, 132, 104011. [Google Scholar] [CrossRef]
  11. Tyler, S.W.; Wheatcraft, S.W. Fractal scaling of soil particle-size distributions: Analysis and limitations. Soil Sci. Soc. Am. J. 1992, 56, 362–369. [Google Scholar] [CrossRef]
  12. Tang, l.; Song, Y.; Jiang, Z.; Jiang, S.; Li, Q. Pore structure and fractal characteristics of distinct thermally mature shales. Energy Fuels 2019, 33, 5116–5128. [Google Scholar] [CrossRef]
  13. Railway Ballast: TB/T2140-2018; China Railway Specifications. China Railway Administration: Beijing, China, 2018.
  14. Tavarez, F.A.; Plesha, M.E. Discrete element method for modelling solid and particulate materials. Int. J. Numer Meth. Eng. 2007, 70, 379–404. [Google Scholar] [CrossRef]
  15. Adhikari, S.; You, Z. 3D discrete element models of the hollow cylindrical asphalt concrete specimens subject to the internal pressure. Int. J. Pavement Eng. 2010, 11, 429–439. [Google Scholar] [CrossRef]
  16. McDowell, G.R.; De Bono, J.P. On the micro mechanics of one-dimensional normal compression. Géotechnique 2013, 63, 895–908. [Google Scholar] [CrossRef] [Green Version]
  17. Ciantia, M.O.; Arroyo, M.; Calvetti, F.; Gens, A. An approach to enhance efficiency of DEM modelling of soils with crushable grains. Géotechnique 2015, 65, 91–110. [Google Scholar] [CrossRef] [Green Version]
  18. Zhou, W.; Wang, D.; Ma, G. Discrete element modeling of particle breakage considering different fragment replacement modes. Powder Technol. 2020, 360, 312–323. [Google Scholar] [CrossRef]
  19. Vallejo, L.E.; Lobo-Guerrero, S. Fractal fragmentation of granular materials under compression. AIP Conf. Proc. 2009, 1145, 847–850. [Google Scholar] [CrossRef]
  20. De Bono, J.P.; McDowell, G.R. Particle breakage criteria in discrete-element modelling. Géotechnique 2016, 66, 1014–1027. [Google Scholar] [CrossRef] [Green Version]
  21. Li, H.Y.; Chai, H.W.; Xiao, X.H.; Huang, J.Y.; Luo, S.N. Fractal breakage of porous carbonate sand particles: Microstructures and mechanisms. Powder Technol. 2020, 363, 112–121. [Google Scholar] [CrossRef]
  22. Lu, M.; McDowell, G.R. The importance of modelling ballast particle shape in the discrete element method. Granul. Matter 2007, 9, 69–80. [Google Scholar] [CrossRef]
  23. Sun, Y.; Indraratna, B.; Nimbalkar, S. Three-dimensional characterization of particle size and shape for ballast. Géotechnique Lett. 2014, 4, 197–202. [Google Scholar] [CrossRef] [Green Version]
  24. Jaeger, J.C. Failure of rocks under tensile conditions. Int. J. Rock Mech. Min. Sci. 1967, 4, 219–227. [Google Scholar] [CrossRef]
  25. Lim, W.L.; McDowell, G.R. Discrete element modelling of railway ballast. Granul. Matter 2003, 7, 19–29. [Google Scholar] [CrossRef]
  26. Indraratna, B.; Ionescu, D.; Christie, H.D. Shear behaviour of railway ballast based on large scale triaxial testing. J. Geotech. Geoenviron. 1998, 5, 439–449. [Google Scholar] [CrossRef]
  27. Roner, C.J. Some Effect of Shape, Gradation and Size on the Performance of Railroad Ballast. Master’s Thesis, Report No. AAR85-324P. Department of Civil Engineering, University of Massachusetts, Amherst, MA, USA, 1985. [Google Scholar]
  28. Marsal, R.J. Mechanical properties of rockfill: In embankment-dam engineering. Int. J. Rock Mech. Min. Sci. 1975, 12, 67. [Google Scholar] [CrossRef]
Figure 1. Fragment replacement model (FRM): (a) single spheres [16]; (b) irregular-shaped particles.
Figure 1. Fragment replacement model (FRM): (a) single spheres [16]; (b) irregular-shaped particles.
Fractalfract 06 00382 g001
Figure 2. Monotonic triaxial test model: (a) test sample; (b) FLAC3D model of membrane.
Figure 2. Monotonic triaxial test model: (a) test sample; (b) FLAC3D model of membrane.
Fractalfract 06 00382 g002
Figure 3. Simulation results under monotonic loading: (a) deviator stress; (b) volumetric strain.
Figure 3. Simulation results under monotonic loading: (a) deviator stress; (b) volumetric strain.
Fractalfract 06 00382 g003
Figure 4. Influence of confining pressure on peak friction angle [26].
Figure 4. Influence of confining pressure on peak friction angle [26].
Fractalfract 06 00382 g004
Figure 5. Simulation results under confining pressure of 240 kPa: (a) sample after loading; (b) evolution of porosity.
Figure 5. Simulation results under confining pressure of 240 kPa: (a) sample after loading; (b) evolution of porosity.
Fractalfract 06 00382 g005
Figure 6. Evolution of particle displacement under 240 kPa confining pressure: (a) εa = 4%; (b) εa = 8%; (c) εa = 12%; (d) εa = 16%; (e) εa = 20%. (Dmax: the maximum displacement; Davg: the average displacement).
Figure 6. Evolution of particle displacement under 240 kPa confining pressure: (a) εa = 4%; (b) εa = 8%; (c) εa = 12%; (d) εa = 16%; (e) εa = 20%. (Dmax: the maximum displacement; Davg: the average displacement).
Fractalfract 06 00382 g006
Figure 7. The orientation of contact within the sample under 240 kPa confining pressure: (a) average contact force; (b) contact number.
Figure 7. The orientation of contact within the sample under 240 kPa confining pressure: (a) average contact force; (b) contact number.
Fractalfract 06 00382 g007
Figure 8. Number of fragments against axial strain under different confining pressures.
Figure 8. Number of fragments against axial strain under different confining pressures.
Fractalfract 06 00382 g008
Figure 9. Variation of particle distribution with grain size.
Figure 9. Variation of particle distribution with grain size.
Fractalfract 06 00382 g009
Figure 10. Effect of confining pressure on: (a) ∆Wk; (b) grain breakage Bg [26].
Figure 10. Effect of confining pressure on: (a) ∆Wk; (b) grain breakage Bg [26].
Fractalfract 06 00382 g010
Figure 11. The development of fragment distribution under confining pressure of 240 kPa: (a) εa = 4%; (b) εa = 8%; (c) εa = 12%; (d) εa = 16%; (e) εa = 20%.
Figure 11. The development of fragment distribution under confining pressure of 240 kPa: (a) εa = 4%; (b) εa = 8%; (c) εa = 12%; (d) εa = 16%; (e) εa = 20%.
Fractalfract 06 00382 g011
Figure 12. Evolution of PSDs: (a) effect of axial strain; (b) effect of confining pressure.
Figure 12. Evolution of PSDs: (a) effect of axial strain; (b) effect of confining pressure.
Fractalfract 06 00382 g012
Figure 13. Linear fitting for lg[M(δ < d)/MT] and lg(d/dM) for simulation sample: (a) effect of axial strain; (b) effect of confining pressure.
Figure 13. Linear fitting for lg[M(δ < d)/MT] and lg(d/dM) for simulation sample: (a) effect of axial strain; (b) effect of confining pressure.
Fractalfract 06 00382 g013
Figure 14. Results of fractal dimension D against grain breakage Bg.
Figure 14. Results of fractal dimension D against grain breakage Bg.
Fractalfract 06 00382 g014
Figure 15. Results of fractal dimension D against volumetric strain εv.
Figure 15. Results of fractal dimension D against volumetric strain εv.
Fractalfract 06 00382 g015
Table 1. Grain size characteristics of ballast.
Table 1. Grain size characteristics of ballast.
Parametersd10/mmd30/mmd50/mmd60/mmdmax/mmCuCcSize Ratio
Value27.132.638.941.3531.50.95.7
Table 2. Simulation parameters in PFC3D and FLAC3D.
Table 2. Simulation parameters in PFC3D and FLAC3D.
Input ParametersValue
Normal and shear stiffness of particles: N/m2 × 106
Friction coefficient of clumps0.5
Density: kg/m32650
Friction coefficient of walls and membrane0.1
Normal stiffness of walls: N/m2 × 107
Shear stiffness of walls: N/m0
Normal stiffness of membrane: N/m2 × 105
Shear stiffness of membrane: N/m0
Young’s modulus of membrane: kPa400
Sample porosity0.44
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, C.; Zhang, X.; Sun, Y.; Zhang, L.; Rui, R.; Wang, Z. Discrete Element Modelling of Fractal Behavior of Particle Size Distribution and Breakage of Ballast under Monotonic Loading. Fractal Fract. 2022, 6, 382. https://doi.org/10.3390/fractalfract6070382

AMA Style

Chen C, Zhang X, Sun Y, Zhang L, Rui R, Wang Z. Discrete Element Modelling of Fractal Behavior of Particle Size Distribution and Breakage of Ballast under Monotonic Loading. Fractal and Fractional. 2022; 6(7):382. https://doi.org/10.3390/fractalfract6070382

Chicago/Turabian Style

Chen, Cheng, Xin Zhang, Yifei Sun, Lei Zhang, Rui Rui, and Zhide Wang. 2022. "Discrete Element Modelling of Fractal Behavior of Particle Size Distribution and Breakage of Ballast under Monotonic Loading" Fractal and Fractional 6, no. 7: 382. https://doi.org/10.3390/fractalfract6070382

Article Metrics

Back to TopTop