A FDEM study on the mechanical properties and failure behavior of soft-hard interbedded rocks considering the size effect

Soft–hard interbedded rocks are widely distributed at the Earth’s surface, and their mechanical properties and failure behavior directly affect the stability of local tunnel and slope engineering projects. Previous studies have rarely considered the influence of size effects on the mechanical properties and failure behavior of such rocks. Therefore, this paper used the combined finite–discrete element numerical method (FDEM) to study the mechanical properties and failure behavior of soft–hard interbedded rock samples considering the size effect. First, appropriate input parameters were calibrated and verified by using a new parameter calibration method. Second, the effects of the element size and loading rate were studied to obtain appropriate model parameters. Finally, the effects of the layer number, sample size, and height–diameter ratio of composite rock samples on their mechanical properties and failure behavior at different layer dip angles and layer thickness ratios were investigated. The results support the following findings: (1) For the composite rock samples with a layer thickness of 10 mm, reliable simulation results can be obtained by using a 2.2 mm element size, and the loading rate should not exceed 0.2 m/s in FDEM numerical modeling. (2) The number of layers in the sample should be at least 5, and when the height–diameter ratio is a constant 2.0, the height of the sample should not be less than 110 mm. (3) As the height–diameter ratio of the composite rock samples increases, both the compressive strength and elastic modulus decrease for all layer dip angles considered, but the rock failure mode changes for layer dip angles of 15–75°; in addition, the sample size effect is most significant for layer dip angles of 60° and 75°. (4) Taking horizontally layered composite rock samples as examples, both the compressive strength and elastic modulus of samples with different layer thickness ratios decrease with increasing height–diameter ratio and their failure modes also depend on the height–diameter ratio.