In this research, distribution of residual stresses in three-phase representative volume element (RVE) consisting of carbon fiber, carbon nanotube and matrix, as three dimensional has been determined. Finite element analysis and cylinder analytical model have been used to determination of residual stresses in every phase. Cylinder analytical model has been developed to three-phase RVE. Coefficient of thermal expansion has been considered differently at various directions in carbon fiber and carbon nanotube phases. Placement of carbon nanotubes has been considered parallel to carbon fiber and as ring around carbon fiber using electrophoresis method. Different volume fraction of carbon nanotube has been considered as 0%, 1%, 2%, 3%, 4% and 5% and distribution of residual stresses of different volume fraction was compared together, but volume fraction of carbon fiber is considered as constant equal 60%. Radial residual stress, tangential residual stress, axial residual stress and residual stress invariant have been determined in different phases. Results of residual stress invariant of two different analysis as finite element and cylinder analytical model have been compared in different phases that are in good agreement together.