Next Article in Journal
Exploration of Zero-Valent Iron Stabilized Calcium–Silicate–Alginate Beads’ Catalytic Activity and Stability for Perchlorate Degradation
Next Article in Special Issue
Effect of Magnesium Salt (MgCl2 and MgSO4) on the Microstructures and Properties of Ground Granulated Blast Furnace Slag (GGBFS)-Based Geopolymer
Previous Article in Journal
On Growth and Morphology of TiO2 Nanotubes on CP-Ti by Anodic Oxidation in Ethylene Glycol Electrolyte: Influence of Electrolyte Aging and Anodization Parameters
Previous Article in Special Issue
Experimental Investigation on the Influencing Factors of Compressive Strength of Foamed Lightweight Material Utilizing Completely Decomposed Granite
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 15
Issue 9
10.3390/ma15093337
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation of Triaxial Shear Behaviors of Cemented Sands with Different Sampling Conditions Using Discrete Element Method

by
Xuqun Zhang
1,
Zhaofeng Li
2,*,
Pei Tai
2,
Qing Zeng
2 and
Qishan Bai
2
1
Guangzhou Metro Design & Research Institute Co., Ltd., Guangzhou 510080, China
2
School of Civil and Environmental Engineering, Harbin Institute of Technology, Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(9), 3337; https://0-doi-org.brum.beds.ac.uk/10.3390/ma15093337
Submission received: 15 April 2022 / Revised: 2 May 2022 / Accepted: 3 May 2022 / Published: 6 May 2022
(This article belongs to the Special Issue Microstructures and Mechanical Properties of Cement-Based Composites)
Browse Figures
Review Reports Versions Notes

Abstract

:
In cemented sand, the influences of the sampling factors (i.e., the curing time, cement–sand ratio, and initial void ratio) on the triaxial shear behavior were investigated using discrete element method. Cemented sand samples with different initial conditions were prepared and subjected to the consolidated drained triaxial shearing test. In the simulations, the peak strength, residual strength, and pre-peak stiffness of cemented sand were enhanced by increasing the curing time and cement–sand ratio, and the enhancements could be explained by the increases in bond strength and bond number. Resulting from the increases of these two sampling factors, bond breakage emerged at a greater axial strain but lower intensity. However, some uncommon phenomena were generated; that is, the contractive but strain-softening response occurred in the sample with a curing time of 3 days, and the shear band and the strain-hardening behavior coexisted in the sample with a cement–sand ratio of 1%. The peak strength and pre-peak stiffness were also enhanced by decreasing the initial void ratio, more distinctly than by increasing the curing time and cement–sand ratio. However, the residual strength, bond breakage, and failure pattern with the persistence of shear band were insensitive to this change.

1. Introduction

In engineering practice, the strength of loose sand can be efficiently enhanced by adding a cementation agent, such as Portland cement, gypsum, and calcite precipitated microbially [1,2,3,4]. The loose sand is then turned into the cemented sand, with bonds usually formed at the contacts between sand particles. Therefore, to provide reliable instructions for sand stabilization by cementation, the behavior of cemented sand should be well understood.
Laboratory tests on cemented sand have shown that the behavior of cemented sand is sensitive to the sampling condition, thus presenting difficulties when developing the associated constitutive model [5,6,7] and evaluating the performance of cement-stabilized soil in ground improvement [8,9]. Consoli et al. [10] suggested the curing time as one of the influential factors for the strength of cemented sand. Li et al. [11] revealed the importance of the cement–sand ratio (i.e., the mass of cement divided by the mass of dry sand) or cement content to the stress–strain response of cemented sand, and, in their tests, the strain-softening behavior was produced when only 1% cement was used in sample preparation. The significance of the initial void ratio was also recognized, and a variable named the cement–void ratio (i.e., the porosity over the cement–sand ratio) was further proposed and found to be positively correlated with strength and stiffness [12,13]. According to these laboratory tests, the sampling factors of the curing time, cement–sand ratio, and initial void ratio appear to play a major role in the mechanical responses of cemented sand, while other factors, such as the curing stress and particle size, may be correlated with the three major ones and only play a minor role [14,15,16].
The above experimental studies have helped to identify the sampling factors which are most significant for cemented sand. However, these studies mainly focused only on the peak strength of cement sand. In fact, to establish the constitutive model, other properties such as the stiffness, residual strength, and volumetric strain are also needed [17,18]. Furthermore, merely by experimentation, these works are insufficient to unravel the regulation mechanisms of the sampling factors toward the soil responses, which are related to the bond breakage at the microscale and are also important for constitutive modeling [19,20,21,22]. Therefore, in parallel with laboratory tests, simulations using the discrete element method (DEM) have been frequently carried out in the study of cemented sand [23,24,25,26]. This is because DEM simulations, on one hand, can provide microscopic insight into the soil behavior [27,28,29] and, on the other hand, can control the sampling and shearing process in an accurate and idealized manner, eliminating the bias from manual operation [30,31,32,33,34]. Using this technique, Li et al. [26] preliminary discussed the influence of curing time (essentially the bond strength) on the stress–strain response of cemented sand; that is, the strain-softening response and the shear band emerged in the long-term cured sample, while the strain-hardening response and the bulging-type failure pattern emerged in the short-term cured sample. Furthermore, the evolution of bond breakage was captured to explain the influence of curing time, which could not be achieved in the laboratory test. However, for other important sampling factors, i.e., the cement–sand ratio and initial void ratio, such a numerical investigation has still not been carried out.
Therefore, in this study, the influences of the sampling factors (i.e., the curing time/bond strength, cement–sand ratio, and initial void ratio) on the triaxial shear behaviors of cemented sand (i.e., the stress–strain relationship, failure pattern, and evolution of bond breakage) were systematically investigated using the numerical technique of the discrete element method, owing to the advantages of this method over the laboratory test. In the remainder of the paper, details of the DEM simulation, including the elements used in the simulation and the simulation process of the triaxial shearing test, are described first. Then, the shear behaviors in different samples, including the stress–strain response, failure pattern, and bond breakage rate, are compared and discussed to understand the influences of the sampling factors.

2. Details of DEM Simulation

Simulations by means of the discrete element method were carried out to produce the shear behaviors of cemented sands with different sampling conditions, which were subjected to the consolidated drained shearing test. Details of the DEM simulations, including the elements used in the simulation and the simulation process, are described below.

2.1. Elements Used in the Simulation

The laboratory tests on cemented sand by Wang and Leung [35] served as the reference to the DEM simulations here. In the laboratory tests, samples were prepared by mixing the Portland cement and Ottawa 20–30 sand, and they were enclosed laterally by a flexible membrane and vertically by rigid walls. In order to reproduce these configurations, each DEM simulation used the elements of rigid walls and three kinds of particles, i.e., sand particles, cement particles, and membrane particles.
Sand particles in the samples followed the size distribution in Figure 1, which was scaled up from that of Ottawa 20–30 sand by 3.0 times in order to make the computational cost affordable [35]. As shown in Table 1, the particle density was set as 2650 kg/m3, the same as that of Ottawa 20–30 sand. A linear contact model was assigned between the sand particles for its simplicity, following the practice in Guo and Zhao [36]. Therefore, the interparticle contact stiffness was linearized using Equation (1), which was based on the Hertz–Mindlin theory and derived by Li et al. [26].
k n = R ⋅ 2 2 π 3 p c G 1 − υ 2 / 3   and   k s k n = 2 1 − υ 2 − υ ,
where kn and ks are the contact normal and tangential stiffness, respectively, pc is the confining pressure and was set as 50 kPa for all the samples, R is the radius of the soil particle and was evaluated as 0.8 mm, i.e., half of D50 in Figure 1, and G and υ are the shear modulus and Poisson’s ratio of soil particles, respectively; the measured values of quartz crystal with G = 29 GPa and υ = 0.2 were used here [37]. Accordingly, the contact normal and tangential stiffnesses were estimated in a linear fashion as 5.0 × 105 and 4.0 × 105 N/m, respectively. The coefficient of friction was given as 0.5, which is also the measured value of quartz crystal [38].
Cement particles were generated around the sand particles, and the parallel bond model was assigned between them to form the cementation in the sample. The radii of the cement particle and cementation bond were set as 1/4 of D50 of the sand particles, in accordance to the X-ray tomography image of cemented sand [39,40]. As can be seen in Table 1, the particle density was given as 3150 kg/m3, i.e., the density of Portland cement. Using these configurations, the real sample could be better reproduced in DEM; that is, different amounts of cementation could be controlled by changing the number of cemented particles, and different curing times could be represented by changing the bond strength. For instance, as Portland cement was simulated here, the bond strengths of 5.0 MPa, 2.5 MPa, and 1.25 MPa could be used to represent the curing times of 28 days, 7 days, and 3 days, respectively [41]. Then, stiffnesses for these three kinds of bonds were determined as 20.5, 31.1, and 82.1 GPa/m, respectively; hence, all bonds broke with a tensile displacement of ~60 μm in the simulation, which is in line with the stretching test results on the cemented bonds with various thicknesses by Jiang et al. [42].
Membrane particles were assembled and linked together by the contact bond model to intimate a thin wall-like membrane boundary in DEM. This particle-based membrane boundary was then used to laterally wrap the cemented sample, i.e., the mixture of sand particles and cemented particles. Particle density was given as 1800 kg/m3, according to the experimentally measured value by Li et al. [43]. Particle size was set as moderate (~1.0 mm) to eliminate the membrane penetration [44]. Contact bond stiffness was assigned as 2.5 × 103 N/m, i.e., 1/200 of the contact stiffness of sand particles, which could eliminate the lateral constraint on the sample during triaxial test [45].
In addition to the membrane particles, rigid walls were used to vertically wrap the sample, which was frictionless to provide the free end boundary condition. Accordingly, the constraint from the rigid wall was also eliminated, and the strain localization or shear band could develop more freely in the sample [46,47,48,49]. The wall stiffness was identical to that of the sand particle, plausibly turning the particle–wall contact into an interparticle contact.

2.2. Simulation Process

The DEM simulation started with the preparation of the cemented sand samples, as shown in Figure 2a. Similar to the practice in the laboratory test, the sand particles and cemented particles were first generated in a rigid mold with a height of 14 cm and a radius of 3.5 cm. The total number of these two kinds of particles was more than 20,000, enabling the sample to serve as a representative element volume and produce the stable shear behavior, according to the simulation results in Kuhn and Bagi [50]. Note that the mold for sample preparation was essentially the combination of a cylindrical rigid wall as the lateral edge and two flat rigid walls as the top and bottom edges. The pressure acting from these rigid walls onto the mixture of sand particles and cemented particles was then adjusted to around 0 kPa. The parallel bond model with the parameters summarized in Table 1 was assigned to the contacts between the two kinds of particles, thereby forming the cemented sand sample.
Since this study was targeted at the influences of three sampling factors (i.e., curing time, cement–sand ratio, and initial void ratio) on the shear behavior of cemented sand, seven different samples were prepared, as summarized in Table 2. Specifically, Samples C28_R5L, C7_R5L, and C3_R5L, which had the same cement–sand ratio of 5% and the same initial void ratio of 0.659, but different bond strength varying from 1.25 to 5.0 MPa, were prepared to investigate the influence of curing time. Note that, as aforementioned, the value of 5.0 MPa could represent a curing time of 28 days for Portland cement, while 2.5 MPa and 1.25 MPa could represent curing times of 7 days and 3 days, respectively. Then, in terms of the cement–sand ratio, Sample C28_R5L served as the reference, and Samples C28_RS3L and C28_R1L were further generated using different numbers of cement particles. Lastly, C28_R5M and C28_R5D were also prepared in order to unravel the influence of initial void ratio on the shear behavior of cemented sand.
After sample preparation in the mold, the rigid lateral wall was removed and replaced by the particle-based membrane, as shown in Figure 2b. To eliminate the leakage of soil particles from the pores of the membrane, the membrane particles were arranged with a hexagon-like pattern, i.e., having each one connected with six neighbors. Details of the establishment process were previously described by Li et al. [45]. Note that the top and bottom ends of the membrane were allowed to deform, which is different to the laboratory setting where these two ends are fastened by O-rings. The purpose was again to reduce the constraint from the boundary condition. Thereafter, with the establishment of the membrane boundary, all samples were consolidated under the confining pressure of 50 kPa, as exemplified in Figure 2c. The confining pressure acting from the rigid wall boundary was realized using the servo-control algorithm proposed by Thornton [51], while that from the membrane boundary was applied using the algorithm proposed by Li et al. [26]. Note that, at the end of consolidation, the void ratio of the sample was measured, taken as the initial value in Table 1.
Finally, the triaxial compression tests were carried out. Samples were sheared under the drained condition by moving the top and bottom walls at a constant rate of 10−3 m/s, as indicated by Figure 2d. With this rate, the inertia number of each sample was maintained at a level far lower than 10−3, demonstrating that the quasi-static condition was well preserved throughout the test. The shearing process was stopped when the axial strain reached 20%. The stress–strain response, failure pattern, and bond breakage rate in the seven prepared samples were measured. Accordingly, in the next section, the influences of the sampling factors (i.e., the curing time, cement–sand ratio, and initial void ratio) on these shearing responses are discussed.

3. Simulation Results and Discussion

3.1. Influence of Curing Time

In order to investigate the influence of the curing time (essentially the bond strength), the shearing responses of Samples C3_R5L, C7_R5L, and C28_R5L are compared in this section. Figure 3 shows the stress–strain curves of these three samples. Obviously, both the peak strength and the pre-peak stiffness could increase with the curing time. The strain-softening phenomenon could be observed in Sample C3_R5L (see the orange solid line in Figure 3), which was cured for only 3 days, and the peak stress of 117 kPa was located at the axial strain of 0.96%. Then, for Sample C7_R5L, the peak stress was 178 kPa and the corresponding axial strain was 1.20%, while these two variables for Sample C28_R5L were 270 kPa and 1.50%, respectively (see the green and black solid lines). A similar effect of curing time on the stress–strain relationship was observed in the laboratory results by Li et al. [26]. Comparing the values of these two variables here and those in Table 2, it can be confirmed that the regulation of bond strength to peak strength exists and is linear to some extent. Similar to the peak strength, the residual strength at the axial strain of 20% could also be enhanced by the increase in curing time, and the associated values of this variable are given in Table 3. All the above improvements by the curing time indicate that stronger bonds can store more energy and can better resist the force-chain buckling, thereby developing a higher level of stress in the sample [52,53,54].
In addition to the soil strength, the volumetric response is affected significantly by the curing time. In Sample C3_R5L, although strain softening occurred, the sample volume still contracted (see the orange dashed line). Note that, after the peak, the sample intended to dilate but failed eventually. The reason for this evolution pattern is that the bonds in Sample C3_R5L were weak and were consumed during shearing. According to Kuhn and Bagi [55], the rigid-body rotation of particle pairs dominates the dilation. Then, owing to the weak bonds being gradually consumed, the particle pairs in the sample might be loose and cannot contribute to dilation effectively and persistently. From this viewpoint, stronger bonds may result in a greater tendency of volumetric dilation. As expected, Sample C7_R5L began to dilate slightly by an ultimate value of −0.43%, while Sample C28_R5L dilated more distinctly by −1.65% (see the green and black dashed lines). Clearly, these observations demonstrate that the influence of the curing time/bond strength on the stress–strain response is complicated, which may pose difficulties in the constitutive modeling of cemented sand.
In addition to the stress–strain response, the curing time plays an important role in the failure pattern. Figure 4 displays the appearances of the three samples at the end of the shearing test. As shown in Figure 4a of Sample C3_R5L with a strain-softening but contractive behavior, the deformation was homogenous. Note that, as mentioned in Section 2.2, the top and bottom ends of the membrane were not fixed to eliminate the constraint from the boundary. It is suggested that, if these two ends were fastened, a bulging failure pattern, which was observed in the laboratory test with the early-age sample by [26], would have been reproduced in Sample C3_R5L. Instead, strain localization was observed of Sample C7_R5L, attributed to the stronger bonds (see Figure 4b). In Sample C28_R5L with a strain-softening and dilative behavior, a shear band near the bottom was further generated, as shown in Figure 4c. This is because the breakage of a strong bond is comparable to the force-chain buckling to some extent, which can initialize the shear banding, according to Tordesillas [52].

3.2. Influence of Cement–Sand Ratio

Responses of Samples C28_R1L, C28_R3L, and C28_R5L are discussed in this section to understand the influence of the cement–sand ratio on the shear behavior of the cemented sand. Figure 5 depicts the stress–strain responses of these three samples. Similar to the effect of increasing the curing time, the peak/residual strength and the pre-peak stiffness could also be enhanced by increasing the cement–sand ratio. In Sample C28_R1L with a ratio of 1%, although the cement was cured for 28 days, a strain-hardening and contractive behavior seemingly appeared, with a peak stress of only 60 kPa; that is, little improvement was produced by the cementation (see the orange lines in Figure 5). The same observation was also found in the laboratory tests of Li et al. [11]. Then, upon increasing the cement–sand ratio from 1% to 3% in Sample C28_R3L, the bond number was increased from 5246 to 17,106, as given in Table 2. Due to this fundamental change at the microscale, the soil response was improved notably, and strain hardening and volumetric dilation were generated (see the green lines). As the cement–sand ratio increased to 5% in Sample C28_R5L, the bond number was increased to 16,141 and the soil response was further improved, manifested by a higher peak/residual strength and volumetric dilation (see the black lines). Note that the specific values of the soil strengths and volumetric strain in these three samples are summarized in Table 3.
Figure 6 presents the failure patterns of the three samples at the axial strain of 20%. Undoubtedly, the pattern was impacted by the cement–sand ratio. Although Sample C28_R1L exhibited a strain-hardening and contractive behavior, strain localization emerged and a shear band seemingly developed, as shown in Figure 6a. Note that the expected result is similar to that of Sample C3_R5L in Figure 4a, where the deformation was homogeneous. Essentially, even though 1% cement was used in the sand, the breakage of strong bonds could still occur through force-chain buckling, resulting in strain localization. As the cement–sand ratio was increased to 3% in Sample C28_R3L, the shear band became discernable in Figure 6b, in line with its strain-hardening and dilative behavior. As the ratio was further increased to 5% in Sample C28_R5L, the shear band became much more distinct in Figure 6c. Therefore, unlike the curing time/bond strength, the influence of cement–sand ratio is complex, related to the aspect of failure pattern, instead of the stress–strain response. This is due to the effect of bond breakage on localization at the microscale, which cannot be revealed by most constitutive models.

3.3. Influence of Initial Void Ratio

To unravel the influence of the initial void ratio, the shear behaviors of Samples C28_R5L, C28_R5M, and C28_R5D are compared in this section. Figure 7 shows the stress–strain responses of these three samples, whereby all samples generated the strain–softening behavior (see all the solid lines). Obviously, with the decrease in the initial void ratio, the peak strength and the pre-peak stiffness could be enhanced. As mentioned in Section 3.1, Sample C28_R5L with an initial void ratio of 0.659 exhibited a peak strength of 178 kPa. In Sample C28_R5M with the initial void ratio decreased by 20% to 0.618, the peak strength was significantly increased to 331 kPa, double that in Sample C28_R5L, while that in Sample C28_R5D further increased to 575 kPa. All these values are summarized in Table 3. Comparing the results in Table 3, it can be found that the peak strength was more sensitive to the initial void ratio than the curing time and cement–sand ratio. However, unlike the peak strength, the residual strength at the end of the shearing test was almost unchanged, i.e., independent of the initial void ratio. This is because the bonds were destroyed in the shearing interface, converting the cemented sand back to be a pure mixture of cement particles and sand particles, which then behaved following the critical state theory.
In line with the strain-softening behavior, all samples exhibited dilative behavior (see the dashed lines). The consistency among these samples implies that once sufficient cement is used and well cured, e.g., using the curing time of 28 days and the cement–sand ratio of 5%, an improvement in the soil, whether loose or dense, can be guaranteed using Portland cement. In contrast, when using an insufficient curing time, e.g., Sample C3_R5L, or insufficient among of cement, e.g., Sample C28_R1L, contractive and strain-hardening behavior may occur. Nevertheless, as the initial void ratio decreased, the volumetric dilation could be further enhanced, and the ultimate values for Samples C28_R5L, C28_R5M, and C28_R5D were −1.65%, −2.70%, and −5.96%, respectively.
Figure 8 shows the failure mode of the three samples. Again, response consistency could be observed among the samples; that is, the shear band persisted near the bottom of each sample. Nevertheless, more soil was involved in the rigid bottom part in Samples C28_R5M and C28_R5L than in Sample C28_R5L, manifested by the shearing interface inclined at a higher angle. To some extent, such a difference in the failure mode is comparable to that in the case of the rigid footing penetrating into the soil foundation; that is, punching failure occurs in the loose sand, which is similar to Figure 8a, and general shear failure occurs in dense sand, which is similar to Figure 8b,c [56].

3.4. Accumulative Bond Breakage

In order to provide more insight into the constitutive modeling of cemented sand, the evolutions of bond breakage for the seven samples in Table 2 are also discussed. Here, the rate of accumulative bond breakage Br is defined in Equation (2) and used in the discussion.
B r = N b N t × 100 % ,
where Nb is the number of accumulative broken bonds, and Nt is the total number of bonds before shearing. Note that Nt for each sample is given in Table 2, and Nb was measured during the shearing process.
Figure 9 plots the evolutions of Br for the seven samples, allowing the influences of different sampling factors to be recognized. Firstly, with the increase in curing time/bond strength, bond breakage may emerge at a greater axial strain and Br can be lowered throughout the shearing process. To be specific, at the axial strain of 5%, Br in Sample C3_R5L reached 19.6% (see the cyan line), while the values in Sample C7_R5L and C28_R5L were reduced to 13.9% and 10.5%, respectively (see the earth-yellow and black lines). Secondly, similar to the curing time, by increasing the cement-sand ratio, Br could also be lowered. Here, at the beginning of the test, the evolution patterns among Samples C28_R1L, C28_R3L, and C28_R5L were similar. However, after the axial strain of 3%, the difference grew. At the end of shearing, Br developed to 44.3%, 26.2%, and 21.5% for these three samples, respectively (see the blue, red, and black lines). Lastly, the evolution of Br appeared to be insensitive to the initial void ratio, since the responses of Samples C28_R5L, C28_R5M, and C28_R5D were similar to each other.

4. Conclusions

The influences of the sampling factors (i.e., the curing time, cement–sand ratio, and initial void ratio) on the triaxial shear behavior of cemented sand were investigated using DEM simulations. The salient findings are as follows:
  • Peak strength, residual strength, and pre-peak stiffness were enhanced by either increasing the curing time or increasing the cement–sand ratio. The enhancements were fundamentally attributed to the increases in bond strength and bond number.
  • Curing time complicated the stress–strain relationship of cemented sand, since strain-softening but contractive behavior was generated in the sample with a curing time of 3 days. Cement–sand ratio disrupted the correlation between the failure pattern and stress–strain evolution pattern, since the shear band occurred in the sample with strain-softening and contractive behavior, which had a cement–sand ratio of only 1%.
  • By decreasing the initial void ratio, the peak strength and pre-peak stiffness can be significantly enhanced, and the shear band may incline at a higher angle. However, the residual strength and failure pattern are insensitive to this change.
  • Bond breakage may emerge later and be less intensive when increasing the curing time. It can also be intensified due to the medium shearing strain by increasing the cement–sand ratio. However, the whole evolution pattern is insensitive to the change in the initial void ratio.
  • Overall, the mechanical behaviors of cemented sand, in terms of the strength, stiffness, and volumetric dilation, were found to be significantly enhanced by increasing the curing time, cement–sand ratio, and packing density. The failure pattern was also changed, attributed to the regulation of the bond breakage at the microscale. These results provide important insight into other cementation methods, such as using gypsum, biopolymer, or MICP.

Author Contributions

Conceptualization, Z.L.; data curation, Z.L.; formal analysis, Z.L.; funding acquisition, X.Z.; investigation, Z.L. and Q.Z.; methodology, Z.L.; resources, P.T.; software, P.T.; supervision, X.Z.; validation, P.T. and Q.Z.; visualization, Z.L.; writing—original draft, Z.L. and Q.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China (2019YFC1511101), Shenzhen Science and Technology Program (Grant No. RCBS20200714114855326, JCYJ20200109112816582, GXWD20201230155427003-20200824104811001, and GXWD20201230155427003-20200824153508001), and Guangdong Basic and Applied Basic Research Foundation (Grant No. 2019A1515110512).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Jun Kang Chow and Yu-Hsing Wang from Hong Kong University of Science and Technology for their support and valuable comments on this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, D.; Park, K.; Kim, D. Effects of Ground Conditions on Microbial Cementation in Soils. Materials 2013, 7, 143–156. [Google Scholar] [CrossRef]
  2. Park, S.-S.; Woo, S.-W.; Jeong, S.-W.; Lee, D.-E. Durability and Strength Characteristics of Casein-Cemented Sand with Slag. Materials 2020, 13, 3182. [Google Scholar] [CrossRef]
  3. Konstantinou, C.; Biscontin, G.; Logothetis, F. Tensile Strength of Artificially Cemented Sandstone Generated via Microbially Induced Carbonate Precipitation. Materials 2021, 14, 4735. [Google Scholar] [CrossRef]
  4. Qin, Y.; Wang, Q.; Xu, D.; Chen, W. Mechanical Behavior and Healing Efficiency of Microcapsule-Based Cemented Coral Sand under Various Water Environments. Materials 2021, 14, 5571. [Google Scholar] [CrossRef]
  5. Gao, Z.; Zhao, J. Constitutive modeling of artificially cemented sand by considering fabric anisotropy. Comput. Geotech. 2012, 41, 57–69. [Google Scholar] [CrossRef]
  6. Diambra, A.; Ibraim, E.; Peccin, A.; Consoli, N.C.; Festugato, L. Theoretical Derivation of Artificially Cemented Granular Soil Strength. J. Geotech. Geoenviron. Eng. 2017, 143, 04017003. [Google Scholar] [CrossRef] [Green Version]
  7. Tan, Y.; Gu, Q.; Ning, J.; Liu, X.; Jia, Z.; Huang, D. Uniaxial Compression Behavior of Cement Mortar and Its Damage-Constitutive Model Based on Energy Theory. Materials 2019, 12, 1309. [Google Scholar] [CrossRef] [Green Version]
  8. Kasama, K.; Whittle, A.J.; Zen, K. Effect of spatial variability on the bearing capacity of cement-treated ground. Soils Found. 2012, 52, 600–619. [Google Scholar] [CrossRef] [Green Version]
  9. Wei, X.; Ku, T. New design chart for geotechnical ground improvement: Characterizing cement-stabilized sand. Acta Geotech. 2019, 15, 999–1011. [Google Scholar] [CrossRef]
  10. Consoli, N.C.; Cruz, R.C.; Floss, M.F. Variables Controlling Strength of Artificially Cemented Sand: Influence of Curing Time. J. Mater. Civ. Eng. 2011, 23, 692–696. [Google Scholar] [CrossRef]
  11. Li, D.; Liu, X.; Liu, X. Experimental Study on Artificial Cemented Sand Prepared with Ordinary Portland Cement with Different Contents. Materials 2015, 8, 3960–3974. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Consoli, N.; Fonseca, A.; Silva, S.; Cruz, R.; Fonini, A. Parameters controlling stiffness and strength of artificially cemented soils. Géotechnique 2012, 62, 177–193. [Google Scholar] [CrossRef] [Green Version]
  13. Consoli, N.C.; Cruz, R.C.; Da Fonseca, A.V.; Coop, M.R. Influence of Cement-Voids Ratio on Stress-Dilatancy Behavior of Artificially Cemented Sand. J. Geotech. Geoenviron. Eng. 2012, 138, 100–109. [Google Scholar] [CrossRef]
  14. Rotta, G.V.; Consoli, N.C.; Prietto, P.D.M.; Coop, M.R.; Graham, J. Isotropic yielding in an artificially cemented soil cured under stress. Geotechnique 2003, 53, 493–501. [Google Scholar] [CrossRef]
  15. Liu, T.; Wang, Z.; Zou, D.; Zhou, A.; Du, J. Strength enhancement of recycled aggregate pervious concrete using a cement paste redistribution method. Cem. Concr. Res. 2019, 122, 72–82. [Google Scholar] [CrossRef]
  16. Guan, W.; Qi, Q.; Zhang, Z.; Nan, S. Effect of Sand Particle Size on Microstructure and Mechanical Properties of Gypsum-Cemented Similar Materials. Materials 2020, 13, 765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Jiang, M.J.; Yu, H.S.; Harris, D. Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses. Int. J. Numer. Anal. Methods Géoméch. 2006, 30, 723–761. [Google Scholar] [CrossRef]
  18. Estrada, N.; Taboada, A. Yield surfaces and plastic potentials of cemented granular materials from discrete element simulations. Comput. Geotech. 2013, 49, 62–69. [Google Scholar] [CrossRef]
  19. Rahimi, M.; Chan, D.; Nouri, A. Bounding Surface Constitutive Model for Cemented Sand under Monotonic Loading. Int. J. Géoméch. 2016, 16, 04015049. [Google Scholar] [CrossRef]
  20. Chen, R.; Huang, J.; Leung, A.K.; Chen, Z.; Chen, Z. Experimental investigation on water release and gas emission of evapotranspirative capillary barrier landfill covers. Soil Sci. Soc. Am. J. 2021, 86, 311–323. [Google Scholar] [CrossRef]
  21. Zhou, C.; Chen, R. Modelling the water retention behaviour of anisotropic soils. J. Hydrol. 2021, 599, 126361. [Google Scholar] [CrossRef]
  22. Gai, X.; Sánchez, M. An elastoplastic mechanical constitutive model for microbially mediated cemented soils. Acta Geotech. 2018, 14, 709–726. [Google Scholar] [CrossRef]
  23. Brugada, J.; Cheng, Y.; Soga, K.; Santamarina, J.C. Discrete element modelling of geomechanical behaviour of methane hydrate soils with pore-filling hydrate distribution. Granul. Matter 2010, 12, 517–525. [Google Scholar] [CrossRef]
  24. Jiang, M.; Yan, H.; Zhu, H.; Utili, S. Modeling shear behavior and strain localization in cemented sands by two-dimensional distinct element method analyses. Comput. Geotech. 2011, 38, 14–29. [Google Scholar] [CrossRef]
  25. Obermayr, M.; Dressler, K.; Vrettos, C.; Eberhard, P. A bonded-particle model for cemented sand. Comput. Geotech. 2012, 49, 299–313. [Google Scholar] [CrossRef]
  26. Li, Z.; Wang, Y.H.; Ma, C.H.; Mok, C.M.B. Experimental characterization and 3D DEM simulation of bond breakages in artificially cemented sands with different bond strengths when subjected to triaxial shearing. Acta Geotech. 2017, 12, 987–1002. [Google Scholar] [CrossRef]
  27. Li, Z.; Wang, Y.; Chow, J.; Su, Z.; Li, X. 3D pore network extraction in granular media by unifying the Delaunay tessellation and maximal ball methods. J. Pet. Sci. Eng. 2018, 167, 692–701. [Google Scholar] [CrossRef]
  28. Chow, J.K.; Tai, P.; Li, J.; Li, Z.; Wang, W. Over-stiff and over-damped problem of multi-sphere approach for ellipse-wall collision using discrete element method. Powder Technol. 2021, 394, 735–747. [Google Scholar] [CrossRef]
  29. Tai, P.; Chen, Z.; Li, Z.; Chen, R.; Lu, H.; Li, Y. Experimental Investigation on the Influencing Factors of Compressive Strength of Foamed Lightweight Material Utilizing Completely Decomposed Granite. Materials 2022, 15, 1060. [Google Scholar] [CrossRef]
  30. Ketterhagen, W.R.; Am Ende, M.T.; Hancock, B.C. Process Modeling in the Pharmaceutical Industry using the Discrete Element Method. J. Pharm. Sci. 2009, 98, 442–470. [Google Scholar] [CrossRef]
  31. O’Sullivan, C. Particle-Based Discrete Element Modeling: Geomechanics Perspective. Int. J. Géoméch. 2011, 11, 449–464. [Google Scholar] [CrossRef]
  32. Niu, Q.; Revil, A.; Li, Z.; Wang, Y.-H. Relationship between electrical conductivity anisotropy and fabric anisotropy in granular materials during drained triaxial compressive tests: A numerical approach. Geophys. J. Int. 2017, 210, 1–17. [Google Scholar] [CrossRef]
  33. Liu, T.; Wei, H.; Zou, D.; Zhou, A.; Jian, H. Utilization of waste cathode ray tube funnel glass for ultra-high performance concrete. J. Clean. Prod. 2019, 249, 119333. [Google Scholar] [CrossRef]
  34. Zhou, A.; Zhang, W.; Wei, H.; Liu, T.; Zou, D.; Guo, H. A novel approach for recycling engineering sediment waste as sustainable supplementary cementitious materials. Resour. Conserv. Recycl. 2021, 167, 105435. [Google Scholar] [CrossRef]
  35. Wang, Y.-H.; Leung, S.-C. A particulate-scale investigation of cemented sand behavior. Can. Geotech. J. 2008, 45, 29–44. [Google Scholar] [CrossRef]
  36. Guo, N.; Zhao, J. Local fluctuations and spatial correlations in granular flows under constant-volume quasistatic shear. Phys. Rev. E 2014, 89, 042208. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, N.; Ciantia, M.O.; Arroyo, M.; Gens, A. A contact model for rough crushable sand. Soils Found. 2021. [Google Scholar] [CrossRef]
  38. Sandeep, C.S.; Senetakis, K. Effect of Young’s Modulus and Surface Roughness on the Inter-Particle Friction of Granular Materials. Materials 2018, 11, 217. [Google Scholar] [CrossRef] [Green Version]
  39. DeJong, J.T.; Fritzges, M.B.; Nüsslein, K. Microbially Induced Cementation to Control Sand Response to Undrained Shear. J. Geotech. Geoenviron. Eng. 2006, 132, 1381–1392. [Google Scholar] [CrossRef]
  40. Mujah, D.; Shahin, M.A.; Cheng, L. State-of-the-Art Review of Biocementation by Microbially Induced Calcite Precipitation (MICP) for Soil Stabilization. Geomicrobiol. J. 2016, 34, 524–537. [Google Scholar] [CrossRef]
  41. Zimbelmann, R. A contribution to the problem of cement-aggregate bond. Cem. Concr. Res. 1985, 15, 801–808. [Google Scholar] [CrossRef]
  42. Jiang, M.; Liu, F.; Zhou, Y. A bond failure criterion for DEM simulations of cemented geomaterials considering variable bond thickness. Int. J. Numer. Anal. Methods Géoméch. 2014, 38, 1871–1897. [Google Scholar] [CrossRef]
  43. Li, Z.; Wang, Y.-H.; Li, X.; Yuan, Q. Validation of discrete element method by simulating a 2D assembly of randomly packed elliptical rods. Acta Geotech. 2017, 12, 541–557. [Google Scholar] [CrossRef]
  44. Membrane penetration effects in triaxial testing: Baldi, G; Nova, R J Geotech Engng Div ASCEV110, NGT3, March 1984, P403–420. Int. J. Rock Mech. Min. Sci. Géoméch. Abstr. 1984, 21, 125. [CrossRef]
  45. Li, Z.; Chow, J.K.; Li, J.; Tai, P.; Zhou, Z. Modeling of flexible membrane boundary using discrete element method for drained/undrained triaxial test. Comput. Geotech. 2022, 145. [Google Scholar] [CrossRef]
  46. Rowe, P.W.; Barden, L. Importance of Free Ends in Triaxial Testing. J. Soil Mech. Found. Div. 1964, 90, 1–27. [Google Scholar] [CrossRef]
  47. Chow, J.K.; Li, Z.; Wang, Y.-H. Comprehensive microstructural characterizations of 1-D consolidated kaolinite samples with fabric tensors and pore elongation factors. Eng. Geol. 2018, 248, 22–33. [Google Scholar] [CrossRef]
  48. Chen, R.; Liu, L.; Li, Z.; Deng, G.; Zhang, Y.; Zhang, Y. A novel vertical stress-controlled apparatus for studying suffusion along horizontal seepage through soils. Acta Geotech. 2021, 16, 2217–2230. [Google Scholar] [CrossRef]
  49. Yuan, Q.; Li, Z.; Gao, Y.; Wang, Y.H.; Li, X. Local responses in 2D assemblies of elliptical rods when subjected to biaxial shearing. Acta Geotech. 2019, 14, 1685–1697. [Google Scholar] [CrossRef]
  50. Kuhn, M.R.; Bagi, K. Specimen Size Effect in Discrete Element Simulations of Granular Assemblies. J. Eng. Mech. 2009, 135, 485–492. [Google Scholar] [CrossRef]
  51. Thornton, C. Numerical simulations of deviatoric shear deformation of granular media. Géotechnique 2000, 50, 43–53. [Google Scholar] [CrossRef]
  52. Tordesillas, A. Force chain buckling, unjamming transitions and shear banding in dense granular assemblies. Philos. Mag. 2007, 87, 4987–5016. [Google Scholar] [CrossRef]
  53. Li, Z.; Wang, Y.-H.; Chow, J.K. Density Effect and Associated Unjamming Events on the Aging-Induced Stiffness Increase in Sand. Int. J. Geomech. 2018, 18, 04018173. [Google Scholar] [CrossRef]
  54. Liu, C.; Tai, P.; Li, Z.; Hu, W. Mechanism of Packing Rigidity Gain in Sand Aging: From the Perspective of Structural Order Evolution. KSCE J. Civ. Eng. 2022, 1–12. [Google Scholar] [CrossRef]
  55. Kuhn, M.R.; Bagi, K. Contact rolling and deformation in granular media. Int. J. Solids Struct. 2004, 41, 5793–5820. [Google Scholar] [CrossRef] [Green Version]
  56. Liang, W.; Zhao, J. Multiscale modeling of large deformation in geomechanics. Int. J. Numer. Anal. Methods Géoméch. 2019, 43, 1080–1114. [Google Scholar] [CrossRef]
Materials 15 03337 g001 550
Figure 1. Grain size distribution of sand particles.
Figure 1. Grain size distribution of sand particles.
Materials 15 03337 g001
Materials 15 03337 g002 550
Figure 2. Simulation process: (a) sample preparation; (b) membrane boundary establishment; (c) consolidation; (d) triaxial compression test.
Figure 2. Simulation process: (a) sample preparation; (b) membrane boundary establishment; (c) consolidation; (d) triaxial compression test.
Materials 15 03337 g002
Materials 15 03337 g003 550
Figure 3. Influence of the curing time on the stress–strain response.
Figure 3. Influence of the curing time on the stress–strain response.
Materials 15 03337 g003
Materials 15 03337 g004 550
Figure 4. Influence of the curing time on the failure mode: (a) Sample C3_R5L; (b) Sample C7_R5L; (c) Sample C28_R5L.
Figure 4. Influence of the curing time on the failure mode: (a) Sample C3_R5L; (b) Sample C7_R5L; (c) Sample C28_R5L.
Materials 15 03337 g004
Materials 15 03337 g005 550
Figure 5. Influence of the cement–sand ratio on the stress–strain response.
Figure 5. Influence of the cement–sand ratio on the stress–strain response.
Materials 15 03337 g005
Materials 15 03337 g006 550
Figure 6. Influence of the cement–sand ratio on the failure mode: (a) Sample C28_R1L; (b) Sample C28_R3L; (c) Sample C28_R5L.
Figure 6. Influence of the cement–sand ratio on the failure mode: (a) Sample C28_R1L; (b) Sample C28_R3L; (c) Sample C28_R5L.
Materials 15 03337 g006
Materials 15 03337 g007 550
Figure 7. Influence of the initial void ratio on the stress–strain response.
Figure 7. Influence of the initial void ratio on the stress–strain response.
Materials 15 03337 g007
Materials 15 03337 g008 550
Figure 8. Influence of the initial void ratio on the failure mode: (a) Sample C28_R5L; (b) Sample C28_R5M; (c) Sample C28_R5D.
Figure 8. Influence of the initial void ratio on the failure mode: (a) Sample C28_R5L; (b) Sample C28_R5M; (c) Sample C28_R5D.
Materials 15 03337 g008
Materials 15 03337 g009 550
Figure 9. The rate of accumulative bond breakage for the seven samples.
Figure 9. The rate of accumulative bond breakage for the seven samples.
Materials 15 03337 g009
Table
Table 1. Parameters of three types of particles and rigid wall for all samples.
Table 1. Parameters of three types of particles and rigid wall for all samples.
ElementsParametersValues
Sand particlesDensity2650kg/m3
Particle radius0.9–3.54mm
Contact normal stiffness5 × 105N/m
Contact tangential stiffness4 × 105 N/m
Coefficient of friction0.5
Cement particleDensity3150kg/m3
Particle radius0.62mm
Coefficient of friction0.5
Bond radius0.62mm
Parallel bond strength1.25–5.0MPa
Parallel bond stiffness20.5–82.1GPa/m
Membrane particlesDensity1800kg/m3
Particle radius1mm
Contact bond stiffness2.5 × 103N/m
Coefficient of friction0.0
Rigid wallsNormal stiffness5 × 105N/m
Coefficient of friction0.0
Table
Table 2. Parameters of the samples with different initial conditions.
Table 2. Parameters of the samples with different initial conditions.
Sample LabelCuring TimeBond StrengthCement-Sand RatioSoil Particle NumberCement Particle NumberBond NumberInitial Void Ratio
(Days)(MPa)(%)(-)(-)(-)(-)
C28_R5L285.05.0561116,14130,4636.59 × 10−1
C7_R5L72.55.0561116,14130,4636.59 × 10−1
C3_R5L31.255.0561116,14130,4636.59 × 10−1
C28_R3L285.03.05705984317,1066.80 × 10−1
C28_R1L285.01.05798333552467.00 × 10−1
C28_R5M285.05.0561116,14131,0636.19 × 10−1
C28_R5D285.05.0561116,14132,2044.87 × 10−1
C: curing time; R: cement–sand ratio; L/M/D: loose/medium dense/dense sample.
Table
Table 3. Responses of the samples with different initial conditions.
Table 3. Responses of the samples with different initial conditions.
Sample LabelPeak Deviator StressAxial Strain of the PeakResidual StressUltimate Volumetric Strain
(kPa)(%)(kPa)(%)
C28_R5L2701.50143−1.65
C7_R5L1781.2095−0.43
C3_R5L1170.96770.86
C28_R3L1591.33103−1.15
C28_R1L601.04541.20
C28_R5M3311.59120−2.70
C28_R5D5751.67169−5.96
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, X.; Li, Z.; Tai, P.; Zeng, Q.; Bai, Q. Numerical Investigation of Triaxial Shear Behaviors of Cemented Sands with Different Sampling Conditions Using Discrete Element Method. Materials 2022, 15, 3337. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15093337

AMA Style

Zhang X, Li Z, Tai P, Zeng Q, Bai Q. Numerical Investigation of Triaxial Shear Behaviors of Cemented Sands with Different Sampling Conditions Using Discrete Element Method. Materials. 2022; 15(9):3337. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15093337

Chicago/Turabian Style

Zhang, Xuqun, Zhaofeng Li, Pei Tai, Qing Zeng, and Qishan Bai. 2022. "Numerical Investigation of Triaxial Shear Behaviors of Cemented Sands with Different Sampling Conditions Using Discrete Element Method" Materials 15, no. 9: 3337. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15093337

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop