This paper presents an embedded beam formulation for discretization independent finite element (FE) analyses of interactions between pile foundations or rock anchors and the surrounding soil in geotechnical and tunneling engineering. Piles are represented by means of finite beam elements embedded within FEs for the soil represented by 3D solid elements. The proposed formulation allows consideration of piles and pile groups with arbitrary orientation independently from the FE discretization of the surrounding soil. The interface behavior between piles and the surrounding soil is represented numerically by means of a contact formulation considering skin friction as well as pile tip resistance. The pile-soil interaction along the pile skin is considered by means of a 3D frictional point-to-point contact formulation using the integration points of the beam elements and reference points arbitrarily located within the solid elements as control points. The ability of the proposed embedded pile model to represent groups of piles objected to combined axial and shear loading and their interactions with the surrounding soil is demonstrated by selected benchmark examples. The pile model is applied to the numerical simulation of shield driven tunnel construction in the vicinity of an existing building resting upon pile foundation to demonstrate the performance of the proposed model in complex simulation environments. J. NINIĆ, J. STASCHEIT AND G. MESCHKE analytical models consider the pile-soil interaction problem by means of an elastically bedded Winkler beam, also known as the p-y approach. In this method, one-dimensional springs are used to describe the 3D soil-pile interaction effects [13][14][15]. Later, this approach was modified to account for plastic deformations of the soil by incorporating nonlinearity in the soil springs [16][17][18][19][20]. In geotechnical engineering, the p-y method has been extensively used in the analysis and the design of piles. Because of rich experimental data used to derive the shape of the p-y curves, this approach has been found effective for many practical applications. Poulos (1973) proposed another analytical, displacementbased method, to analyze the lateral pile response. In this method, a free-field soil movement profile is imposed on a pile in a simplified boundary element analysis to estimate the pile response [21]. Using such an analysis method associated with laboratory testing of piles and pile groups, a series of design charts and theoretical expressions have been developed for estimating the pile response for slope stabilization, excavation, and tunneling operations [2,9,22,23].Evidently, analytical models are based upon more or less significant simplifications of the mechanisms involved in the 3D pile-soil interaction and the interactions between soil strata.Finite element (FE) models have a comparably larger potential to realistically describe material nonlinearities, different loading cases and complex geological environments in analyses of soil-pile interaction. An early ...