Reduced-order homogenization (ROH) and related methods are important computational tools for simulating the material behavior of composites. These methods generally sacrifice accuracy in exchange for superior computational efficiency, relative to methods such as classical computational homogenization (CCH). In this study, building on the recently developed solver-free CCH, we propose a fine-scale solver-free reduced order homogenization approach that avoids solving the fine-scale equilibrium equations and approximates the phase-average eigenstrains by sampling the fine-scale eigenstrain at a small number of points. The proposed method, which we call solver-free ROH, works by pre-computing history-dependent eigenstrain influence function tensors and sampling point contribution factors based on training data from a small set of CCH simulations. Then, during the online stage of the computation, phase-average eigenstrains are computed from sampling point eigenstrains and used to compute homogenized stresses and strains. Focusing on small-deformation problems, this paper formulates and verifies the solver-free ROH approach for nonlinear periodic heterogeneous media. First, in the formulation, we delineate a sampling point approximation of the phase eigenstrains and describe the use of this approximation within the coarse-scale stress update. Next, we verify the solver-free ROH using loading cases outside the training data set. Finally, we use the proposed method to simulate a multilayer composite plate in three point bending (3pt-bend) and open hole tension (OHT), demonstrating the method’s efficiency and accuracy relative to the CCH.