Modern high speed printing machines are able to print up to 700 m/min. At this rate, little excitations lead to vibrations, which may lead to loss of contact between the rollers (bouncing). This bouncing results in white stripes, being visible on the printed image. To enable the simulation of the whole printing process, including effects like bouncing, a discrete multibody model is developed. The rollers are modeled by several rigid bodies. These bodies are connected to each other by rotational springs, which allow simulation of the first bending eigenmodes of each roller. The contact area between the rollers is modeled by several nonlinear translational springs and damping elements. These elements change their stiffness and damping values depending on the distance between the rollers. If a defined distance is exceeded, the values become zero, which represents the loss of contact (bouncing). The unknown spring and damping elements of this model are parametrized with help of an experimental modal analysis. This paper presents the development of a flexible multibody model to simulate nonlinear effects in printing process.