Abstract: In order to reveal the influence of bedding properties on the stability of coal mine roadway, the study established a 3D numerical model combined finite differential method (FDM) coupling discrete element method (DEM). In the model, the two sides of the roadway, the floor and the far-field surrounding rock are simulated by the continuous medium domain created by FLAC^3D, while the layered roof rock mass of the roadway is created by the bonded particle model (BPM) in PFC^3D, and the bedding planes are realized by adding parallel weak planes with specific stiffness and strength. In the study, the properties of the bedding were simulated by varying the strength of the weak plane to investigate the fracture characteristics of the layered roof under the action of excavation and unloading. The results showed that the fracture distribution in the intact roof was dome-shaped, while the fractures in the layered roof were distributed along the weak plane. The failure of the layered roof started from the layer adjacent to the goaf and expanded inward. The increase in bedding strength will change the stress field distribution within the roof, and to alter the development process of fractures in the roof. The coupling strategy of PFC^3D/FLAC^3D can effectively reproduce the fracture development process of the layered roof in coal mines under the action of roadway excavation and unloading.