IBM (Immersed boundary method) and DEM (Discrete element method) coupling method were used to simulate the flow of cylinder-shaped particles in a fluidized bed. The greatest advantage of IBM-DEM is that it can reveal the microscopic characteristics of dense-phase gas-particle flow in Cartesian grids. Large cylinder-shaped particles are very difficult to fluidize, and slugging flow can be observed even if the static bed height is low. The gas flow field around the particle in fluidized bed is analyzed, and the formation and development of vortex behind the particle is affected obviously by the neighboring particles. Particle trajectory is obtained, and the effect of gas phase on particle rotation becomes active as particle size increases. Compared with the experimental results, the gas-solid force in simulation results is higher. This calculation error may be reduced by decreasing the grid size.