The behavior of fiber suspensions during flow is of fundamental importance to the process simulation of discontinuous fiber reinforced plastics. However, the direct simulation of flexible fibers and fluid poses a challenging two-way coupled fluid-structure interaction problem. Smoothed Particle Hydrodynamics (SPH) offers a natural way to treat such interactions. Hence, this work utilizes SPH and a bead chain model to compute a shear flow of fiber suspensions. The introduction of a novel viscous surface traction term is key to achieve full agreement with Jeffery’s equation. Careful modelling of contact interactions between fibers is introduced to model suspensions in the non-dilute regime. Finally, parameters of the Reduced-Strain Closure (RSC) orientation model are identified using ensemble averages of multiple SPH simulations implemented in PySPH and show good agreement with literature data.