In this paper, a humoral and cellular immunity virus dynamics model with the Beddington-DeAngelis incidence rate is set up. We derive the basic reproductive number R0, the cytotoxic T lymphocytes immune response reproductive number R1, the humoral immune response reproductive number R2, humoral immune response competitive reproductive number R3, and cytotoxic T lymphocytes immune response competitive reproductive number R4, and a full description of the relation between the existence of the equilibria and reproductive numbers is given. The global properties of the five equilibria are obtained by constructing Lyapunov functions. Copyright © 2014 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]