using the code UDEC of discrete element method (DEM), with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures.