To effectively simulate the fracture propagation in shale, the bedding plane (BP) effect is incorporated into the augmented virtual internal bond (AVIB) constitutive relation through BP tensor. Comparing the BP-embedded AVIB with the theory of transverse isotropy, it is found the approach can represent the anisotropic properties induced by parallel BPs. Through the simulation example, it is found that this method can simulate the stiffness anisotropy of shale and can represent the effect of BPs on hydraulic fracture propagation direction. Compared with the BP-embedded virtual internal bond (VIB), this method can account for the various Poisson's ratio. It provides a feasible approach to simulate the fracture prop- agation in shale.