Abstract: Based on Voxel finite element mesh, the stress analysis of spherical inclusion composites was conducted. The obvious stress concentration phenomenon of elements neighboring interface was induced by the stepped interface within Voxel mesh. A local stress averaging method was proposed to treat the stress concentration induced by Voxel finite element mesh. In addition, the effects of stress averaging domain, stress averaging weight function and mesh density were taken into account. It is found that the local stress averaging method can effectively reduce the stress concentration of elements neighboring interface. But the averaging domain size should be not too big or too small. The two depth layer Voxel elements are the best choice as the averaging domain through calculation, which also has mesh independence. The method can be further used in the progressive damage evolution analysis of the spherical inclusion composites.