The behavior of interphase‐particle adhesion and interphase region around the nanoparticles can significantly affect the stress distribution and mechanical properties of polymeric nanocomposites. In this study, the elastic modulus of epoxy/silica nanocomposites is analyzed using the finite element method and different mathematical models. A nanoscale representative volume element including graded interphase, homogenous interphase, and no interphase model is implemented. Furthermore, the effect of interfacial adhesion is also considered. The final elastic modulus was clearly affected by the interphase modulus, especially at higher nanoparticle content. Under imperfect interfacial bonding, the existence of an interphase region leads to a slight increase in modulus, and in the absence of that area, the elastic modulus decreases to 3.28 GPa. In perfect bonding models, stress transferred from the matrix to interphase and, then, to nanoparticle, which led to a significant increase in elastic modulus. Unlike the imperfect bonding, the maximum stress was located in the elements along to the loading direction. A maximum 26% increase in elastic modulus for perfect bonding/graded interphase model with 6.54 vol% of nanosilica particles compared to bulk epoxy was achieved. Finally, on comparing the FEM analysis and theoretical results with the experimental data, good agreement between obtained results was found.