Non-spherical particles are widely used in the chemical and pharmaceutical industries. Often these particles are over-simplified as equivalent spherical particles for calculation of drag forces and particle-particle interactions. We have developed a three-dimensional discrete element method (DEM) for rigid cylindrical fibrous particles with a high aspect ratio. In this method, a cylindrical particle was represented as overlapped multiple spheres. The diameter of the spheres was the same as that of the fibrous particle, while the number of spheres was determined by the length of the fibrous particle. The simulations were carried out for different numbers of particles to study the dynamics of fibrous sedimentation and also to investigate the effect of particle-particle interactions on the average preferred orientation and terminal velocities of the particles. The simulation results are compared to experimental data with good agreement. It is shown that the particles interactions have a significant effect on the particles average orientation and terminal velocity in freely falling particles. However, for fibrous particles moving in a jet flow, the interactions have negligible influence on the particle orientation and terminal velocity. (C) 2018 Elsevier B.V. All rights reserved.