The effective thermal conductivity of composite materials with thermal contact resistance at interfaces is studied by lattice Boltzmann modeling in this work. We modified the non-dimensional partial bounce-back scheme, proposed by Han et al. [Int. J. Thermal Sci., 2008. 47: 1276-1283], to introduce a real thermal contact resistance at interfaces into the thermal lattice Boltzmann framework by re-deriving the redistribution function of heat at the phase interfaces for a corrected dimensional formulation. The modified scheme was validated in several cases with good agreement between the simulation results and the corresponding theoretical solutions. Furthermore, we predicted the effective thermal conductivities of composite materials using this method where the contact thermal resistance was not negligible, and revealed the effects of particle volume fraction, thermal contact resistance and particle size. The results in this study may provide a useful support for materials design and structure optimization.