This paper introduces a mathematical model which describes the dynamics of the spread of HIV in the human body. Human immunodeficiency virus infection destroys the body immune system, increases the risk of certain pathologies, damages body organs such as the brain, kidney and heart or cause the death. Unfortunately, this infection disease currently has no cure to control the diseases. We propose a fractional order model in this paper to describe the dynamics of human immunodeficiency virus (HIV) infection. The Caputo fractional derivative operator of order ∈ (0,1] is employed to obtain the system of fractional differential equations. The basic reproductive number is derived for a general viral production rate which determines the local stability of the infection free equilibrium. The solution of the time fractional model has been procured by employing Laplace Adomian decomposition method (LADM) and the accuracy of the scheme is presented by convergence analysis Moreover, numerical simulation are performed to study the dynamical behavior of solution of the models. Simulations of different epidemiological classes at the effect of the fractional parameters revealed that most undergoing treatment join the recovered class. The results show the both viral production rate and death rate of infected cells play an important role the disease in the society.