In this paper, we construct a numerical method for a delayed viral infection model with general incidence rate. We prove that the obtained discrete model has the same dynamics as the corresponding continuous model, such as positivity, boundedness and global behaviors of solutions with no restriction on the time step size. Furthermore, numerical simulations are given to illustrate and confirm our main analytical results.