An elastoplastic model that exhibits hybrid initial and stress-induced anisotropy is developed for layered rock masses. The formulation considers the thermodynamic aspect, and is consistent and rigorous. Both initial anisotropy and stress-induced deformation anisotropy are reflected by a hybrid anisotropic stiffness matrix influenced by the deterioration development degree (DDD) in different directions. By employing a combination of the strength criterion of rock material and rock bedding plane, an anisotropic failure formulation for layered rock mass has been established. The hybrid anisotropic model has been implemented in cellular automata software for the engineering rock mass fracturing process (CASRock). The performance of the anisotropic part is demonstrated by reproducing the deformation and failure characteristics of initial or stress-induced anisotropic behaviors for layered rocks under uniaxial, conventional triaxial, and true triaxial compression and Brazilian splitting conditions. Important features, such as the strength, mechanism, deformation, DDD, and fracturing process variation, can be captured by the proposed model. In addition, a numerical simulation of tunnel excavation in a layered rock mass is performed to study the anisotropic excavation-induced damage zone (EDZ) distribution in the field. The results indicate that the model is able to reproduce the observed failure mode satisfactorily.