© 2022, The Author(s), under exclusive licence to Società Italiana di Fisica and Springer-Verlag GmbH Germany, part of Springer Nature.Antiferromagnetic spin-1 XYZ model is examined by using a mean-field approach with the introduction of spin operators on the simple cubic lattice. The model includes the crystal field interaction (Dz) along the z-axis, the Dzyaloshinskii–Moriya interaction (Δ m) and an external magnetic field with the components of Hx= Hy= Hz= H. The bilinear exchange interaction parameter (Jz) is taken as a scaling parameter chosen to be negative to simulate the antiferromagnetic interactions between the nearest-neighbor spins. Thermal variations of the total magnetization and its components are investigated in detail to obtain the phase diagrams on the (H/ | Jz| , T/ | Jz|) , (Dz/ | Jz| , T/ | Jz|) and (Δ m/ | Jz| , T/ | Jz|) planes. The model exhibits antiferromagnetic, paramagnetic and random phase regions. Very interesting various phase lines and critical points are observed including the tricritical points, bicritical points, critical end points and two more. The reentrant behavior is also observed for appropriate values of the system parameters.