Timoshenko Beam Model Predicts When Thick Hard Mine Roofs Will Crack

Deep beneath the Kouzidong Coal Mine in China, a slab of fine sandstone nearly nine metres thick hangs over the coal face like a loaded spring. When it finally breaks, the release is violent: rib spalling more than a metre deep, roof falls between hydraulic supports, and roadways deforming so severely that steel straps fracture and bolts pull free of the shattered coal. A new study published in Discover Geoscience offers a way to anticipate exactly where and at what angle such a roof will fail, combining a century-old structural mechanics theory with modern digital imaging to turn roof fracture from a surprise into a calculation.

The research team, led by Hao Yue and Xinghai Lei with colleagues Qiang Li and Fengqi Liu, tackled a long-standing blind spot in mining engineering. For decades, engineers have modelled the rock layers above a longwall working face using the Euler–Bernoulli beam theory, which assumes a beam’s cross-section stays planar and perpendicular to its neutral axis as it bends. That assumption works well for thin, slender roof strata. But it quietly ignores shear deformation, and as roof thickness grows, the overburden stops behaving like a slender beam and starts behaving like a short, stubby one. The researchers show that with a periodic weighting interval of 30 metres, once roof thickness exceeds roughly 3.75 metres, the slenderness ratio falls outside the range where Euler–Bernoulli theory remains valid.

Their solution was to rebuild the roof model on Timoshenko beam theory, which allows the cross-section to rotate independently of the beam axis and explicitly accounts for shear deformation and rotational inertia. Crucially, the team did not stop at a single two-dimensional slice. Previous studies typically analysed fracture either along the direction of face advance or along the face length, but rarely both at once. By establishing fracture models in both the advancing and lengthwise directions, the researchers could reconstruct the full three-dimensional fracture morphology of the overlying strata, characterised by two parameters for each layer: the advanced fracture distance, where the break occurs ahead of the face, and the fracture angle, which governs how the broken blocks interlock.

The mechanical setup couples the Timoshenko beam with a Winkler elastic foundation, treating the interaction between the main roof and the softer immediate roof beneath it as a bed of linear springs whose reaction is proportional to deflection. This is a deliberate improvement over the rigid-foundation assumptions of earlier models, because it better reflects the real stress state of the rock mass, even though rock is ultimately elastoplastic rather than purely elastic. During the first weighting, when the roof spans an unsupported opening for the first time, the roof above the coal is modelled as a beam on an elastic foundation while the span above the goaf is treated separately, with boundary and continuity conditions stitching the two solutions together. During periodic weighting, the goaf side loses its lateral constraint entirely, so the roof is treated as a plate with three elastically supported edges and one free edge.

Solving the resulting differential equations yields closed-form expressions for the bending moment along the roof, and the fracture position follows from the point of maximum moment. The sensitivity analysis reveals trends with direct operational consequences. The advanced fracture distance decreases as the elastic foundation coefficient and the roof’s elastic modulus increase, but grows with roof thickness, meaning a thicker roof breaks farther ahead of the working face and demands a longer advance support distance to protect the coal wall. The fracture angle tells the opposite story in places: it decreases with increasing weighting interval and foundation coefficient, but increases with roof thickness and elastic modulus, while Poisson’s ratio barely matters. Because the fracture angle controls whether the broken roof forms a stable voussoir beam or fails by sliding or rotation, these trends translate directly into hazard prediction. A larger fracture angle favours sliding instability; a smaller one favours rotational collapse.

Stage differences matter too. Both the advanced fracture distance and the fracture angle during periodic weighting come out smaller than during first weighting, a consequence of the goaf side’s lost constraint giving the roof more rotational freedom. The practical implication is that advance support ranges should be larger for the first weighting than for subsequent ones. The model also echoes field observations from China’s Shendong Mining Area, where weighting intervals under 20 metres correspond to stepped subsidence and sliding instability with severe dynamic loading, exactly what the theory predicts for short intervals and large fracture angles.

To test the theory against physical reality, the team built a two-dimensional similarity model on a 180 by 140 by 16 centimetre platform, scaled at 1:150 geometrically, with materials of sand, cement, and gypsum proportioned to reproduce the strength hierarchy of the Kouzidong strata. A digital image correlation system with two five-megapixel CCD cameras filming at 75 frames per second tracked the principal strain field as the model face advanced. The results were strikingly quantised: the immediate roof failed at a principal strain threshold of about 2 percent, the main roof at about 6 percent, and the key stratum at about 8 percent. Higher thresholds correspond to stronger layers, which means the key stratum accumulates the most strain energy before failing and releases the most energy when it does, a finding that flags key-stratum fracture as the moment demanding pre-emptive reinforcement.

The experiment reproduced the full failure sequence in miniature. The immediate roof collapsed first at 30 centimetres of excavation, the main roof broke for the first time at 45 centimetres and then fractured periodically at roughly 15-centimetre intervals, and at 125 centimetres the key stratum failed with an audible crack, dragging extensive overlying failure with it. When the researchers imported images of the final fracture morphology into MATLAB and binarised them using Otsu’s method, the overlying strata formed an approximately trapezoidal fracture pattern, with the open-off cut side fracturing at about 75 degrees and the advancing side at about 70 degrees, a five-degree difference attributed to the hydraulic supports altering the effective foundation stiffness rather than to image-processing error, which a threshold sensitivity test bounded at plus or minus 2 degrees.

The agreement between theory and experiment was close enough for engineering use. For the initial weighting, the model predicted an advanced fracture distance of 2.1 metres and a fracture angle of 76 degrees, against experimental values of 2.25 metres and 72 degrees. For the overall fracture angle, the theoretical 76 degrees versus the measured 70 degrees represents a relative error of about 8.6 percent, comfortably within the 10 to 20 percent natural fluctuation of field weighting intervals, and small enough that both values point to the same instability mode and the same control strategy.

The study closes the loop from prediction to control. At Kouzidong, the team designed a pressure-relief blasting scheme with fan-patterned boreholes in the roadways, inclined at 45, 60, and 75 degrees with depths exceeding 18 metres, and numerical simulation showed the peak abutment stress coefficient dropping from 2.1 to 1.92 after blasting. They further proposed dividing the 266-metre-long face into three zones: reinforced end control zones within 10 metres of each boundary using special supports with 40 to 70 percent more support area, a high-resistance support zone where intact roof demands high initial support forces to prevent bed separation, and a central pressure-stabilised zone covering the middle half of the face where fragmented roof is protected by moving supports in coordinated groups of two or three, orchestrated through a remote intelligent control platform. Together, the Timoshenko-based fracture model and the zonal control scheme offer mines facing thick, hard roofs a quantitative route to safer, more predictable extraction.

Subject of Research: Fracture mechanics and stability control of thick hard roof strata over longwall coal mining faces using Timoshenko beam theory and physical simulation

Article Title: Study on fracture mechanism and stability control of longwall working face with thick hard roof in coal mine

Article References: Yue, H., Lei, X., Li, Q., & Liu, F. (2026). Study on fracture mechanism and stability control of longwall working face with thick hard roof in coal mine. Discover Geoscience, 4(1), Article 393. https://doi.org/10.1007/s44288-026-00763-4

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00763-4

Keywords: longwall mining, thick hard roof, Timoshenko beam theory, Winkler elastic foundation, roof fracture angle, first weighting, periodic weighting, digital image correlation, voussoir beam, pressure-relief blasting, hydraulic supports, coal mine safety