The settlement of a group pile foundation is the result of elastic-plastic settlement of the pile and rheological settlement of the soil. Unfortunately, most research only considers elasto-plasticity settlement of piles or rheological settlement of soils. Classical theories like the Poynting-Thompson model only provide some prediction for the elasto-plasticity behavior of pile groups, which cannot satisfactorily reflect the total settlement characteristics of pile group foundation. In the current paper, we presented a method to study settlement behaviors of pile group foundation in linear viscoelastic soil based on the pile-pile interaction factors, the calculated equation of settlement of single piles and pile groups under vertical loads that are derived, the influence law of pile end resistance, pile spacing, slenderness ratio, elastic modulus ratio, and Poisson ratio on pile end load ratio, and the interaction factors that have been presented. Finally, some typical predictive examples have been presented for a 2 × 2, 3 × 3 group pile showing settlement from a field test. It is found that the pile-pile interaction factors increase with the increase of the pile-soil elastic modulus ratio and length-to-diameter ratio but decrease with the increase of the ratio of distance to diameter, and the ratio of distance to diameter has little effect on the interaction factor when pile spacing is more than 10 d (d is the diameter of the pile). The error between the field test results and computed data is 14.9%, the settlement calculated by the present method shows good agreement with the reported results, and the correlativity coefficients R2 are greater than 0.9223, which signifies that the present method is suitable for the theoretical calculation and prediction of the composite foundation settlement of linear viscoelastic soil when piles are used to strengthen the foundation. The method would provide a unique view to quantitative prediction of the settlement characteristic of group pile foundation.