Motion of fibres in sheared fibre suspensions is simulated numerically by using the lattice Boltzmann method. The orientational distributions of the fibres are presented for different Reynolds numbers, Stokes numbers, shear rate and fibre aspect ratio. Some computational results are compared with the experimental data of pipe Bow, and the qualitative agreement is achieved. The results show that the orientational distributions are greatly affected by the Reynolds numbers, while relatively insensitive to the fibre aspect ratio. The Stokes number and shear rate have obvious influence on the orientation distribution.