In this work, a hard sphere particle method for non-spherical particles is applied to study particle behaviors in a bubbling fluidized bed. The influence of the non-spherical particle method and the restitution coefficient are examined. The particle behaviors, including the overall fluidization state (distributions of particle translational and rotational velocities, void fraction, etc.), bubble and particle fluctuations (particle and bubble granular temperatures), are studied. Specific variation rules which are difficult to measure experimentally are presented. As the bubble boundaries in the non-spherical particle system are not as clear as in the spherical case, there is a stronger fluctuation movement of bubbles, which also dominates the particle turbulent movement. The absolute values of particle translational and rotational velocities increase when the restitution coefficient decreases. And it affects the periodic fluctuation characteristics of gas-solid interactions. The variation rules for bubbles and particles as a function of the restitution coefficient are different. (C) 2016 Elsevier B.V. All rights reserved.