A new simulation method for solving fluid-structure coupling problems has been developed. All the basic equations are numerically solved on a fixed Cartesian grid system in a finite difference scheme. A volume-of-fluid approach, which has been developed for computing multiphase flows, is applied to describing the multi-component geometry. The temporal change in the solid deformation is described on the Eulerian frame by updating a left Cauchy-Green deformation tensor, which represents constitutive equations for the hyperelastic Cauchy stress. The present simulation method is applied to problems of hydrodynamic interactions among several neo-Hookean materials.