This study employs similar simulation testing and discrete element simulation coupling to analyze the failure and deformation processes of a model coal seam's roof. The caving area of the overburden rock is divided into three zones: the delamination fracture zone, broken fracture zone, and compaction zone. The caving and fracture zones' heights are approximately 110 m above the coal seam, with a maximum subsidence of 11 m. The delamination fracture zone's porosity range is between 0.2 and 0.3, while the remainder of the roof predominantly exhibits a porosity of less than 0.1. In addition, the numerical model's stress analysis revealed that the overburden rock's displacement zone forms an 'arch-beam' structure starting from 160 m, with the maximum and minimum stress values decreasing as the distance of advancement increases. In the stress beam interval of the overburden rock, the maximum value changes periodically as the advancement distance increases. Based on a comparative analysis between observable data from on-site work and numerical simulation results, the stress data from the numerical simulation are essentially consistent with the actual results detected on-site, indicating the validity of the numerical simulation results.