Hepatitis B virus (HBV) infection is a major public health concern throughout the world. It can be treated effectively through proper medication and vaccination, although it is very difficult to cure, especially in people who have had the chronic infection. So, it is necessary to pay proper attention for its eradication. In this study, two nonlinear fractional‐order HBV infection models with recycling effects of capsids are presented via the Caputo derivative. Each model contains four compartments, namely, susceptible hepatocytes, infected hepatocytes, HBV capsids, and viruses. In the first model, the order of fractional derivatives is equal for each compartment, whereas, in the second model, the order is incommensurate. The existence and uniqueness of solutions for both models are discussed separately. The parameters are estimated in order to validate the proposed model with experimental data obtained from a chimpanzee. Stability analyses are carried out for both models theoretically. The models are solved numerically using the predictor–corrector Adams–Bashforthm–Moulton method (for commensurate order) and implicit product integration of trapezoidal type method (for incommensurate order) with various choices of fractional orders and initial conditions. All the results are presented graphically. Interestingly, the results reveal the importance of fractional‐order derivatives in capturing the dynamics of HBV transmission in the host. It is also noticed that with the decrease in the order of fractional derivative, the peak level of infection decreases, but the disease takes a long time to be cured.