The present article provides the description of the solvation forces between large spheres in a fluid. The molecular dynamics (MD) method was applied to the relatively simple systems in which a pair of structureless macroparticles, either solvophobic or solvophilic, is immersed in a simple fluid of two types, either a soft-sphere or a Lennard-Jones fluid. When a pair of solvophobic macroparticles was in the attractive Lennard-Jones fluid, no dense layer of the solvent particles formed near the surface of the macroparticles and the strong attractive forces were induced between them. In the other combinations of macroparticles and fluids, the dense layers formed and the solvation forces oscillated, exhibiting the attraction and repulsion, whose periodic distance was about the diameter of solvent particles. Our results agreed well with those of the other simulation and theoretical studies with respect to the solvent density profile near a macroparticle and the force-distance profile between macroparticles. The benefit of our approach would be the simplicity in specifying or finding the bulk condition that is in equilibrium with the thin film of molecules between large surfaces. The present method can be applied straightforward to macroparticles immersed in mixtures and complex fluids described by the bead-spring model, to which the conventional grand canonical ensemble Monte Carlo (GCEMC) method is hardly accessible.