A numerical procedure is described for predicting steady drift forces an multiple three-dimensional bodies of arbitrary shape freely floating in waves. The developed numerical approach is based on combination of a three-dimensional source distribution method, wave interaction theory art the far-field method using momentum theory. Numerical results are compared with the experimental or numerical ones, which are obtained in the literature, of steady drift forces on 33(3 by 11) floating composite vertical cylinders in waves. The results of comparison confirmed the validity of the proposed approach. Finally, the interaction effects are examined in the case of an array of 40(4 by 10) freely floating rectangular bodies in shallow water.