The mechanical properties and failure modes of shale that contains bedding were analyzed under coupled seepage pressure and stress. In the numerical model, seven sets of shale samples with different dip‐bedding values were considered, and a seepage‐stress coupling numerical simulation was carried out. The results showed that the compressive strength of shale exhibited obvious anisotropy with increase in the bedding angle, and the compressive strength decreased with an increase in the osmotic pressure. The elastic modulus of shale showed an increasing trend with increase in the bedding angle, while the osmotic pressure had little effect on the elastic modulus of shale. The results also showed that weak bedding had a significant impact on the shale failure mode. When the osmotic pressure was 4 MPa, the shale fracture mode showed five failure modes. When the osmotic pressure was 8 MPa, the shale fracture mode exhibited four failure modes. When the osmotic pressure was 12 MPa, the shale fracture mode showed four failure modes. As the osmotic pressure increased, the osmotic pressure played a leading role in the shale rupture process. The spatial distribution of the acoustic emission demonstrated some self‐similarity patterns, and its fractal characteristics were analyzed. It was found that the fractal dimension better reflects the rupture damage to the shale. The larger the fractal dimension is, the more severe the rupture of shale.