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.