Abstract
In this paper, we propose an (n + 3)-dimensional nonlinear virus dynamics model that describes the interactions of the virus, uninfected target cells, n-stages of infected cells and B cells. We assume that the incidence rate of infection, the generation and removal rates of all compartments are given by general nonlinear functions. We derive two threshold parameters, the basic reproduction number, R-0 and the humoral immunity number, R1 and establish a set of conditions on the general functions which are sufficient to determine the global dynamics of the model. Utilizing Lyapunov functions and LaSalle's invariance principle, the global asymptotic stability of all steady states of the model is proved. Numerical simulations are conducted for specific forms of the general functions in order to illustrate the dynamical behavior.