In coal masses, joints serve as crucial components influencing mechanical responses and failure mechanisms. The joint distribution complicates the acquisition of physical parameters for engineering-scale coal masses, and traditional laboratory or in situ tests often suffer from inaccuracies and high costs. Therefore, examination of the dynamic properties of engineering-scale coal mass is generally performed by numerical simulations. This paper proposes a model construction approach using two-dimension Particle Flow Code (PFC2D) software to study the mechanical properties of engineering-scale jointed coal mass, addressing the limitations of conventional models by integrating scale effects, accuracy, and computational efficiency. Firstly, the distribution characteristics and mechanical parameters of the joints in the coal mass were obtained based on field statistics and laboratory experiments. The parameters of the laboratory-scale model were calibrated by the numerical matching method. The discrete element model for the engineering-scale coal mass was constructed by the step-by-step matching method. The confining pressure effect on the coal mass under a biaxial loading condition was studied, while the strength change, fissure evolution, and failure mechanism under different confining pressures and fissure degrees were investigated. Based on the simulation results, a quantitative relationship was established between the mechanical parameters, fissure degree, and confining pressure under compression conditions. Ultimately, the failure zone ahead of the working face and the distribution of the abutment pressure were assessed using the mechanics parameters of coal masses with diverse joint distributions.

To address severe stress concentration, excessive convergence, and instability of the roadside backfill body (RBB) in gob-side entry retaining (GER) under thick and hard roof conditions, this study investigates the control mechanism of main roof fracture position on surrounding rock stability, using the 3⁻¹01 working face of Huoluowan Coal Mine as a case study. A combined approach integrating theoretical analysis, numerical simulation, and field investigation is adopted. A statically indeterminate mechanical model based on masonry beam theory is established to characterize the lateral roof fracture behavior. The deflection and bending moment distributions are derived, and a criterion for fracture position determination is developed based on the maximum bending moment condition. The theoretical results indicate that the natural fracture position is located approximately 9.4–11.2 m inside the gob boundary. Numerical simulations using UDEC Trigon under different fracture positions (−2 m, 1 m, 5 m, and 9 m) show that fracture location significantly affects the mechanical response of GER. Fractures occurring above the roadway or RBB induce large deformation levels and more extensive plastic zones, while gob-side fracture conditions correspond to relatively lower disturbance levels and improved structural stability. The RBB exhibits shear-dominated failure characteristics, and the displacement distribution is non-uniform along height, with larger deformation in the middle-to-upper region. To improve stability, a coordinated control strategy combining anchor cable reinforcement and directional long-distance hydraulic fracturing (HF) is proposed to regulate the main roof fracture position through the formation of artificial weak planes. Field monitoring results show that the maximum displacements of the roof, floor, and ribs are 558 mm, 233.5 mm, and 71.3 mm, respectively, with a convergence ratio of 19.8%. Borehole imaging confirms the development of hydraulic fractures within the designed roof stratum, supporting the effectiveness of the proposed control approach. These results demonstrate that the fracture position of the main roof plays a key role in controlling GER stability, and its regulation provides an effective means for improving roadway performance under complex geological conditions.