A full Eulerian simulation method for solving fluid-structure coupling problems has been developed. It facilitates solution of dynamic interaction problems between Newtonian fluid and hyperelastic material for given initial configuration of a multi-component geometry described by voxel-based data on a fixed Cartesian mesh. The monolithic velocity field defined over the fluid and solid domains is discretized on a fixed mesh in a finite difference manner. A solid volume fraction and a left Cauchy-Green deformation tensor are temporally updated on the Eulerian frame to distinguish between two phases and to quantify the deformation level of solid, respectively. The simulation method is applied to two-dimensional motions of biconcave neo-Hookean particles in a Poiseuille flow.