Summary
This paper presents a simple three‐dimensional (3D) Distinct Element Method (DEM) for numerical simulation of the mechanical behavior of bonded sands. First, a series of micro‐mechanical tests on a pair of aluminum rods glued together by cement with different bond sizes were performed to obtain the contact mechanical responses of ideally bonded granular material. Second, a 3D bond contact model, which takes into account the influences of bond sizes, was established by extending the obtained 2D experimental results to 3D case. Then, a DEM incorporating the new contact model was employed to perform a set of drained triaxial compression tests on the DEM bonded specimens with different cement contents under different confining pressures. Finally, the mechanical behavior of the bonded specimens was compared with the available experimental results. The results show that the DEM incorporating the simple 3D bond contact model is able to capture the main mechanical behavior of bonded sands. The bonded specimen with higher cement content under lower confining pressure exhibits more pronounced strain softening and shear dilatancy. The peak and residual strengths, the apparent cohesion and peak/residual friction angles, and the position and slope of the critical state line increase with increase in cement content. Microscopically, bond breakage starts when the system starts to dilate and the maximum rate of bond breakage coincides with the maximum rate of dilation. Bond breakage is primarily due to tension‐shear failure and the percentage of such failures is independent of both confining pressure and cement content. Copyright © 2015 John Wiley & Sons, Ltd.