In this manuscript, Eshelby tensor is employed to assess the strain concentration that arises in the matrix phase at the interface, offering precise values and locations of maximum strain under specific loading conditions for both spherical and cylindrical inclusions. When compared to numerical simulation results, the analytical predictions grounded in the Eshelby tensor exhibit satisfactory accuracy. Then an analytical calculation method based on Eshelby tensor for the elastic strain influence functions of reduced-order homogenization (ROH) method is developed and adopted on particle-reinforced and fibrous composites, presenting its feasibility and advantage on off-line stage calculation of ROH method. The error analyses between analytical and numerical results are conducted. The numerical results also exhibit the necessity of finer interface partitioning to obtain the response on micro-scale with higher resolution.