Ensuring the stability of slopes is critical to the safe operation of geotechnical engineering. Evaluating slope stability to minimize geologic risks induced by destabilization is significant in reducing casualties and property damage. A conventional, single-coefficient strength reduction method is widely applied in slope stability analyses, but this method ignores the attenuation degree of different parameters in the slope destabilization. A new double-strength reduction method considering different contributions of the mechanics’ parameters is proposed in this study for evaluating the stability of nonhomogeneous slope. First, the role of each mechanic’s parameters in the slope destabilization was investigated theoretically and numerically using ABAQUS software 2022. The results indicate that the effect of elasticity (E), Poisson’s ratio (v), and soil gravity (γ) on the evolution of factor of safety (FOS) are insignificant and can be neglected compared with cohesive force (c), and angle of internal friction (φ). Next, an improved method was constructed to correlate the FOS with cohesive force (c) and the angle of internal friction (φ). Then, a numerical method was constructed based on the computation of the mathematical–mechanical relationship between FOS and the mechanical parameters, and the stability of slope is estimation based on the Mohr–Coulomb yield criterion. Finally, the double-strength reduction coefficient method proposed in this study, the limit equilibrium method, and the traditional finite element strength reduction coefficient method were applied to nonhomogeneous slopes and slopes containing a soft underlying layer for comparison, and the difference between them was within the range of ±5%. The results indicate that both the limit equilibrium method and the traditional finite element strength reduction method tend to overestimate the FOS of intricate slopes compared with the evaluated method proposed in this study. This parallel comparison serves to validate the accuracy of the double-strength reduction method proposed in the present study. Further, based on the proposed method, the relationship between slope stability and slope displacement is established, which provides a theoretical basis for the safety assessment of slope engineering.

The mechanical evolution of coal is evident when the pore pressure and the surrounding stress alone influence it. However, the evolution of the mechanical response of saturated coal under the coupling effect of pore pressure and confining pressure needs further investigation. This study identifies the mechanical behaviors of burst tendency dry and saturated coal under the stress condition where confining and pore pressure simultaneously increase but keep the constant difference by conducting a series of triaxial compressions on high burst tendency dry and saturated coal samples. The results show that the elastic modulus (E) and strength (σ_peak) of dry coal increase from 3.4 to 4.8 GPa and 78.5 to 92.6 MPa, respectively, and the macro shear failure angle decreases from 64.2° to 56.5° when the confining pressure increases from 9 to 15 MPa. However, these parameters show the opposite evolution law when the pore pressure increases. Furthermore, the E and σ_peak of saturated coal decrease from 3.84 to 2.75 GPa and 73.4 to 60.3 MPa, respectively, and the macro shear failure angle of saturated coal increases from 64.7° to 72.4° when the confining pressure and pore pressure increase simultaneously. The coefficient μ is proposed to reveal the evolution of strength at the effective confining pressure. Furthermore, the Mohr–Coulomb failure criterion, including μ, is ameliorated for application in coal under pore pressure conditions. In addition, a model was developed to reveal the effect of a pore-rich layer on the angle of macrocracks, which was confirmed by acoustic emission. The research reveals the mechanical behavior of coal under high pore pressure. Improved Mohr–Coulomb criterion criteria provide new guidance and vision for the analysis of coal instability in high pore pressure coal seams.

Understanding the behavior of fluid flow and solute transport in fractured rock is of great significance to geoscience and engineering. The discrete fracture network is the predominate channel for fluid flow through fractured rock as the permeability of fracture is several magnitudes higher than that of the rock matrix. As the basic components of the fracture network, investigating the fluid flow in crossed fractures is the prerequisite of understanding the fluid flow in fractured rock. First, a program based on the successive random addition algorithm was developed to generate rough fracture surfaces. Next, a series of fracture models considering shear effects and different surface roughness were constructed. Finally, fluid dynamic analyses were performed to understand the role of flowrate and surface roughness in the evolution of flow field, concentration field, solute breakthrough, and solute mixing inside the crossed fractures. Results indicated that the channeling flow at the fracture intersection became more pronounced with the increasing Péclet number (Pe) and Joint Roughness Coefficient (JRC), the evolution of the concentration field was influenced by Pe and the distribution of the concentration field was influenced by JRC. For Pe < 10, the solute transport process was dominated by molecular diffusion. For 100 > Pe > 10, the solute transport process was in the complete mixing mode. In addition, for Pe > 100, the solute transport process was in the streamline routing mode. The concentration distribution was affected by the local aperture at the fracture intersection corresponding to different surface roughness. Meanwhile, the solute mixing equation was improved based on this result. The research results are beneficial for further revealing the mechanism of fluid flow and solute transport phenomenon in fractured rock.

Equivalent hydraulic aperture and fracture surface roughness are two significant factors affecting the fluid flow behaviors in rock fractures. To understand the role of fracture surface roughness and aperture in the fluid flow through 3D self-affine rough fractures, roughness fracture surfaces with joint roughness coefficients equal to 2.5, 7.5, 12.5, and 17.5 were established, and the Navier–Stokes equation was used to compute the fluid flow in these 3D self-affine rough fractures with a mechanical aperture increase from 0.2 mm to 0.8 mm with a gradient of 0.2 mm. The results show that when the fracture mechanical aperture is 0.2 mm, the impact of fracture surface roughness on fluid flow is considerable, while this effect decreases obviously with the increase of fracture mechanical aperture. Comparing the permeability obtained by the Navier–Stokes equation with the cube law under different hydraulic gradients, we found that their deviation increased with the increase of hydraulic gradients. This allows for the definition of a critical hydraulic gradient (), below which the permeability can be properly predicted using the cubic law for its simplicity, and above which Forchheimer’s equation should be applied, and Forchheimer’s coefficients , , and can be calculated by the prediction equations established in this study.

Freeze–thaw action is a key factor in the deterioration of the dynamic mechanical behavior of rocks in cold regions. This study used yellow sandstone, which is prevalent in the seasonally cold region of Xinjiang, China. The yellow sandstone samples were subjected to various temperatures and a range of freeze–thaw cycles. Impact mechanical tests were performed using a Split Hopkinson Pressure Bar (SHPB) system on the treated samples. The effects of freezing temperature and changes in impact load on the mechanical properties of frozen–thawed sandstone were examined. Additionally, the damage fractal characteristics of the sandstone were analyzed using fractal theory. The results indicate that as the freezing temperature decreases, the stress–strain curves of frozen–thawed specimens exhibit a clear initial compaction stage. The dynamic strength of the specimens decreases with lower freezing temperatures and shows a logarithmic relationship with the loading strain rate; however, the dynamic deformation modulus exhibits no significant correlation with the strain rate. The fractal dimension is positively correlated with the strain rate, indicating that lower freezing temperatures correspond to a higher rate of increase in the fractal dimension. These findings offer valuable insights into the damage deterioration characteristics of frozen–thawed rocks under varying temperature conditions.