We model bedding plane slip to predict the spatial distribution of thrust faults.
We find that slip along bedding contact surfaces perturbs the local stress field.
Strain energy density and Coulomb stress are calculated as failure indicators.
Simulations predict increased propensity for failure where thrust faults are found.
This modeling approach may help predict hazardous roof conditions in mines.