In the current development of flexible multibody dynamics the efficient and accurate consideration of distributed and non-linear forces are an active area of research. Examples are forces due to body-body contact or due to elastohydrodynamics. This leads to many additional modes for representing the local deformations, which describe the potential area on which the forces act. Recent publications show that these can be several hundred to several thousand additional modes. A conventional, monolithic numerical time integration scheme would lead to unacceptable computing times. This paper presents a method for an efficient time integration of such systems. The core idea is to treat the equations associated with modes representing local deformations separately. Using the Newmark formulas, a fixed point iteration is proposed for these separated equations, which can always be stabilized with decreasing step size. The concluding examples underline this property, as well as the fact that the proposed method massively outperforms the conventional, monolithic time integration with increasing number of modes.