This thesis presents a numerical analysis of pile groups deflections, under axial load in layered soils. The method uses a simplified version of the boundary element technique, as proposed by Poulos and Davis (1968) for a homogeneous, isotropic and linear elastic half-space, based on Mindlin s equation. In this case, the major step consisted in replacing the multi-layered system for an equivalent homogeneous half-space using the well known procedures suggested by Poulos (1979)and Yamashita et al (1987). Results of the deflection of pile groups obtained with the proposed method were compared to results obtained using finite element analysis and to results published in the technical literature obtained by several authors. It has been shown that the proposed method provides results which compare favorably with existing solution and basically presents the advantage of the simplification in handling the soil and geometrical datas necessary for the problem description. Interaction factors for groups of pile were produced for some specific cases consisting of a double layer system.