A multi-step homogenization approach is presented to predict the off-axis creep response of hybrid polymer matrix composites (HPMCs) reinforced with unidirectional carbon fibers and silica nanoparticles. The first step deals with evaluating the viscoelastic properties of silica nanoparticle-polymer nanocomposite using the Mori-Tanaka micromechanical model. Two essential features affecting the behavior, including silica nanoparticle agglomeration and interphase region generated due to the interaction between the nanoparticle and polymer are taken into account. In the second step, the off-axis viscoelastic behavior of HPMCs is extracted from the homogenized nanocomposite and carbon fiber properties using a unit cell-based micromechanical model. Some comparative studies between the predictions and available experiment are directed to verify the homogenization process. All the model predictions are in good agreement with the experimental data. The results indicate that with increasing the fiber off-axis angle from 0° to 90°, the presence of silica nanoparticles leads to a reduction in the HPMC creep compliance. Also, the proposed multi-step homogenization approach is applied to investigate the effects of volume fraction, size and agglomeration degree of silica nanoparticles; thickness and material properties of the interphase region; and off-axis angle and volume fraction of the carbon fiber on the HPMC creep response.