Purpose. Creating a generalized algorithm to account for factors (coal seam thickness, enclosed rock mechanical properties, the dimension and bearing capacity of artificial support patterns) causing a gateroad state under the effect of longwall face and goaf. Methods. The assessment of the gateroad stability is based on numerical simulation of the rock stress-strain state (SSS). The finite element method is used to find out the changes in the SSS of surrounding rocks at various stages of longwall mining. The elastic-plastic constitutive model and Hoek-Brown failure criterion implemented in codes RS2 and RS3 (Rocscience) are applied to determine rock displacements dependently on the coal seam thickness, enclosed rock strength, width and strength of artificial support (a packwall comprised of hardening mixture “BI-lining”). To specify the mechanical properties of the packwall material a series of experimental tests were conducted. A computational experiment dealing with 81 combinations of affecting factors was carried out to estimate the roof slag and floor heaving in the gateroad behind the longwall face. A group method of data handling (GMDH ) is employed to generalize the relationships between rock displacements and affecting factors. Findings. The roof-to-floor closure in the gateroad has been determined at the intersection with the longwall face and goaf dependently on the coal seam thickness, enclosed rock strength, width of the packwall, and strength of hardening material. It is revealed that the support material gains the strength value of 30 MPa on the 3rd day from its beginning to use which is fully corresponding to the requirements of protective element bearing capacity. The possibility of using untreated mine water to liquefy the mixture is proved, that allows simplifying and optimizing the solute mixing and pumping technology. Originality. This study contributes to improving the understanding of the factors that influence the stability of underground mining operations and highlights the importance of utilizing numerical simulations in optimizing mining designs. The impact of each factor on the resulting variable (decrease in cross-section of gate road by height) based on the combinatorial algorithm of structural identification of the model is estimated as follows: the packwall width is 48%, the thickness of coal seam is 25%, the strength of enclosing rocks is 23%, and the strength of the packwall material is 4%. Practical implications. The findings provide stakeholders with a technique to determine reasonable parameters for support and protective systems, and the predictive model developed can be used to mitigate potential instability issues in longwall mining excavations. The results have implications under similar geological settings and can be valuable for mine design and optimization in other regions.