In this paper, a unified behavioral model is used to simulate the behavior of anisotropic rocks. The proposed unified behavioral model employs the non-associated flow rule and the concept of critical state. Furthermore, the proposed model is developed based on the concept of the sub-loading surface to predict the smooth transition behavior from elastic to plastic states. To implement the model, the implicit Euler method is utilized. For modeling anisotropy, a 13-plane model is employed, with each of these planes exhibiting elastoplastic behavior. The overall behavior of the soil is calculated by assessing the separate behavior of multiple individual planes in different directions at a point. A set of unconventional constitutive equations is used separately for each plane. Using this model, the induced anisotropic effects are simulated. Subsequently, three rock samples under monotonic loading in drained conditions are simulated with the unified behavioral model and compared with laboratory data. For each of the tests conducted at different confining pressures, the deviatoric stress-axial strain and volumetric strain-axial strain graphs were plotted. Based on the results obtained from the numerical simulations, the strain softening behavior, maximum stress, and the behavior of rock samples before and after maximum stress. It was shown that the developed model can be effectively used to simulate the behavior of rocks.