In order to study the movement of rigid bodies with arbitrary shapes in the fluid, a resolved CFDEM coupling algorithm is proposed. The proposed method combines the computational fluid dynamics for the fluid and the discrete element method for the rigid bodies. The behavior of the fluid phase is simulated based on the Navier-Stokes equations in the Eulcrian framework, whereas the motion of the rigid bodies is governed by the Newton's second law in the Lagrangian description. The major difficulty, namely, the representation of the moving boundaries of the solids, is solved by the immersed boundary method. The interaction forces between the fluid and the solids in different frameworks are derived according to the velocity boundary conditions on the immersed boundary points. Aiming at achieving the strong interaction effects between the fluid and the solids, several iterations are carried out to solve the coupling system. (C) 2020 Elsevier B.V. All rights reserved.