The construction of the new airport of Mexico City (NAICM) is challenging, as the subsoil at the site is very weak. The airport, which is under construction at the moment, is located on top of thick layers of Mexico Clay soils. These clayey layers consist of a mixture of both clay particles and pyroclastic volcanic materials. These unique soils are characterised by a high water content and compressibility. Upon loading, the soft soils will settle as a result of the consolidation process. This consolidation is defined as the dissipation of excess pore water as a result of an increase in stress. The rate of consolidation is governed by the hydraulic conductivity of the soil, which is in general very low in soft soils, and the drainage path. In order to speed up this time consuming process, vertical permeable elements can be installed in the soil. Upon installation of these elements, a reduction in drainage path is realised. This reduction in drainage path allows faster dissipation of excess pore water pressure, hence reducing the consolidation time. Due to the conditions at the site, the design of the runways is challenging. The differential- and residual settlements need to be within strict boundaries in order to keep the runways operational. Prior to the design of the runways, trial embankments have been constructed in order to measure the performance of different soil improvement techniques at this specific site. Based upon this field test results, prefabricated vertical drains (PVDs) were selected as the most efficient soil improvement technique. Besides the PVD trial, the sandpile trial also showed promising results. The PVDs and sandpiles both decrease the drainage path length inside the soft soil layers. The main difference is the relative high stiffness of the sandpiles, the stiffness of the PVDs is negligible. By installing stiff elements in the soft soil, the loads applied on the top of the soils are transferred to deeper and stiffer soils. Hence, the total settlements of the top soft soil layers are reduced. The sandpiles could therefore be feasible alternative to PVDs, applied at locations in which the total settlement needs to be limited. The investigation into the behaviour of sandpiles in Mexico Clay soils is performed in multiple stages. First, analytical models are used for gaining knowledge about the expected sandpile behaviour. In the second stage the FEM software Plaxis 2D is used for the numerical modelling of the sandpiles, using an axisymmetric model setup. The model is constructed in multiple steps, with increasing complexity. The first models simulate the material behaviour with the Mohr-Coulomb model, in the advanced models the Soft Soil Creep and Hardening Soil models are used. The final numerical model is verified by the measurement data obtained from the field trials. After fitting the numerical model to the field data, the verified model is used in the sandpile sensitivity analysis. In this analysis the sensitivity of the material properties and geometry of the sandpile on the performance are researched. The residual settlements or performance, which is defined as the difference in settlement after construction and over a period of 8 years, is affected by the length, radius and centre to centre (ctc) distance of the sandpile. Adjusting the pile stiffness has minimal effect on the performance of the sandpile. An optimum in performance is found by varying both the pile radius and ctc distance. A combination of a small pile diameter with a small ctc distance results in the best performance. When comparing the performance of both sandpile and PVD, the PVD is found to be more effective in terms of performance. The additional stiffness of the sandpile is not reducing the total settlements, on the contrary, the self-weight of the piles increases the total settlements by providing an additional load to the soft soil layers located underneath the pile tip. In conclusion the sandpiles have no additional benefit over the use of PVDs, therefore the application of sandpiles at the NAICM site is not a feasible alternative to the use of PVDs.

Due to increased traffic density, there is a tendency to construct larger immersed tunnels with larger tunnel elements both in transversal and longitudinal direction resulting in a more critical design of the joints. Especially for those tunnels with relatively soft foundation and high siltation, there is no mature experience for reference. This not only requires a more detailed analyze of the joint itself, but also of the total behavior of the tunnel and its interaction with the subsoil. Conventional deterministic designing procedure is based on average soil characteristics with sensitivity analyses including upper and lower boundaries of soil stiffness. The mentioned above trend of larger structures requires a more sophisticated probabilistic approach with deeper understanding of uncertainties and insight in probability of exceedance of requirements in serviceability limit state. In this paper, a probabilistic analysis approach has been implemented to understand the shear behavior of tunnel joints. The biggest contributor for tunnel joints shear behavior is differential settlement between adjacent parts of the tunnel structure. Two probabilistic methods, Point Estimate Method and Monte Carlo Simulation are performed on settlement determination. The former is more efficient and with a considerable accuracy while, the latter is extremely accurate with high computational costs. For simple cases, Monte Carlo Simulations are commonly used to solve soil structure interaction problems due to its straightforward process and easily understandable theory. While for large and complex soil-structure interaction problem, it is computationally intensive to complete even a single run. Such practical disadvantage can be solved by modifying the algorithm or by performing computationally efficient probabilistic methods. The soil-structure interaction analysis is researched by models set up in PLAXIS 2D. Python scripts are programmed in cooperation with PLAXIS models to perform Monte-Carlo Simulation in determination of shear forces in tunnel joints. The function of Python scripts is helping PLAXIS to select properties automatically and to storing the output after every single run. The reliability analysis is done for different conditions in 3 models. Comparison shows that reliability significantly influenced by the distribution and correlation length of soil parameters.