This work presents a modified analytical model for the bending behavior of axially functionally graded (AFG) carbon nanotubes reinforced composite (CNTRC) nanobeams. New higher order shear deformation beam theory is exploited to satisfy parabolic variation of shear through thickness direction and zero shears at the bottom and top surfaces.A Modified continuum nonlocal strain gradient theoryis employed to include the microstructure and the geometrical nano-size length scales. The extended rule of the mixture and the molecular dynamics simulations are exploited to evaluate the equivalent mechanical properties of FG-CNTRC beams. Carbon nanotubes reinforcements are distributed axially through the beam length direction with a new power graded function with two parameters. The equilibrium equations are derived with associated nonclassical boundary conditions, and Navier’s procedure are used to solve the obtained differential equation and get the response of nanobeam under uniform, linear, or sinusoidal mechanical loadings. Numerical results are carried out to investigate the impact of inhomogeneity parameters, geometrical parameters, loadings type, nonlocal and length scale parameters on deflections and stresses of the AFG CNTRC nanobeams. The proposed model can be used in the design and analysis of MEMS and NEMS systems fabricated from carbon nanotubes reinforced composite nanobeam.