Complex jointed rock masses are common in nature. When a roadway is located in complex jointed rock mass, numerous joints can weaken the strength of surrounding rock, resulting in large deformation of the roadway and impact on mine production. In this paper, a DFN–DEM jointed roadway model was constructed based on discrete element method, with the failure characteristics and layout optimization of roadway analyzed. Also, the mechanical properties of jointed rock masses with different sizes and rotation angles were numerically investigated. And the failure and deformation characteristics of jointed roadway with different layout angles were analyzed. The results demonstrated that the UCS of jointed rock masses under different loading directions are significantly different, with an RVE size of 6 m. UCS of the jointed rock masses reaches maximum value at the rotation angles of 30° and −60°, and the average UCS of rock masses has a positively relationship with the joint density P₃₀. The shear instability of joints results in the macro-failure of rock masses, and brings about the asymmetric failure of jointed roadway. The layout angle of jointed roadway can significantly affect its failure and deformation characteristics. With the roadway axial direction adjusted to −60° relative to the direction of principal stress, the deformation degree of the surrounding rock reached minimum values. In engineering practice, an optimal layout angle and supporting of jointed roadway are effective for avoiding the large deformation of surrounding rock.

The prevention and control of stress relief gas have crucial influence on safety of coal mine, it is not only the inevitable requirement of safe production, but also an effective way to realize reasonable utilization of gas. Taking the 1103 working face of the Weijiadi coal mine as the background, the mining stress field was analysed by means of numerical simulation, theoretical analysis and field practice, and the control technology of stress relief gas was studied. The results showed that scope of stress relief zone drops gradually associated with an increase of distance from the roof (or floor) to the working face. Additionally, the shape of the stress relief body exhibited a ring-shaped distribution, while four corners of the goaf roof and floor underwent high permeability zones due to a deep stress relief body, where the permeability of two corners near the transportation roadway of the floor was higher. The results provided good information for W-shaped ventilation mode with two inlets and one return which was adopted in the working face. More importantly, the optimized layout of boreholes was put forward, which eventually were useful for solving the gas overrun of the working face. The technology used in 1103 working face has an attractive and practical background with other extensive applications for the prevention and control of relief gas.

Rock masses with a distinct structure may present a transversely isotropic character; thus, the stress state in a transversely isotropic elastic half-plane surface is an important way to assess the behavior of the interaction between the distributed loading and the surroundings. Most previous theoretical analyses have considered a loading direction that is either vertical or horizontal, and the stress distribution that results from the effect of different loading directions remains unclear. In this paper, based on the transversely isotropic elastic half-plane surface theory, a stress solution that is applicable to distributed loading in any direction is proposed to further examine the loading effect. The consistency between the analytical solution and numerical simulations showed the effectiveness of the proposal that was introduced. Then, it was utilized to analyze the stress distribution rule by changing the Poisson’s ratio and Young’s modulus of the model. The effects of the formation dip angle on the stress state are also examined. The stress distribution, depending on the physical property parameters and relative angle, is predicted using an analytical solution, and the mechanisms associated with the transversely isotropic elastic half-plane surface subjected to the loading in any direction are clarified. Additionally, extensive analyses regarding this case study, with respect to the mechanical behavior associated with changes in the stress boundary, is available. Hence, the proposed analytical solution can more realistically account for the loading problem in many engineering practices.

Jointed rocks under local load are ubiquitous in civil engineering. The instability and failure of jointed rocks are fatal to engineering safety. This paper numerically investigated the effects of loading area and joint angle on the strength dividing points, energy evolution, and crack distribution characteristics of non-persistent jointed rocks. The results demonstrated that the closer the absolute value of joint angle to 45° and the smaller the loading area, the lower the strength dividing points of rocks. The curves of rock joint angle versus total energy at peak and of elastic energy versus amplitude of post-peak abrupt energy change render a W-shape distribution. Meanwhile, compared with joint angle, loading area has more influence on rock energy input. The larger the loading area, the higher the crack fractal dimension, the crack entropy, and the penetration rate. Tensile cracks outnumber shear cracks when jointed rocks are damaged, and shear cracks increases significantly at the post-peak stage.