FINITE ELEMENT MODELLING OF ENERGY PILES USING DIFFERENT CONSTITUTIVE MODELS FOR THE SOIL Mouadh Rafai (mouadh.rafai@studenti.unipg.it) Università degli Studi di Perugia Dipartimento di Ingegneria Civile ed Ambientale Arianna Lupattelli (arianna.lupattelli@studenti.unipg.it) Università degli Studi di Perugia Dipartimento di Ingegneria Civile ed Ambientale Diana Salciarini (diana.salciarini@unipg.it) Università degli Studi di Perugia Dipartimento di Ingegneria Civile ed Ambientale ABSTRACT. In this paper a 2D Finite Element model is presented for the simulation of the coupled Thermo-Hydro-Mechanical (THM) processes induced by the operation of energy piles. The model, implemented in the Plaxis 2D Finite Element (FE) code, was used to simulate the behaviour of a well-known instrumented energy pile during a test performed at Lambeth College in London (Bourne-Webb et al., 2009), comparing the results obtained from different constitutive models used to describe the soil behaviour, namely: i) the simple elasto-plastic Mohr- Coulomb (MC) model; ii) the Modified Cam Clay (MCC) model; iii) the Hardening Strain (HS) model. Comparisons with the experimental results show that all the models, if correctly calibrated, are able to realistically reproduce the behaviour of the pile-soil system from a qualitative point of view, while the best prediction was found using the HS model. 1. INTRODUCTION Energy piles are an emerging technology, with several successful applications throughout the world (see, i.e., Brandl, 2006). They allow the exploiting of renewable and clean source of energy for air conditioning of buildings (Amatya et al., 2012; Murphy et al., 2014; Laloui et al., 2019). In urban areas, heating and cooling of buildings represents one of the main sources of energy consumption and one of the main causes of CO 2 production. The low enthalpy geothermal plants equipped with heat pumps provide one of the most efficient and sustainable solutions. These plants exploit the deep foundations as heat exchangers in direct contact with the soil, which is at constant temperature (10-16°C) during the year. Energy piles are ordinary reinforced concrete piles where a system of pipes is installed to allow the circulation of a heat transfer fluid to extract heat from the subsoil during winter, or release heat to the ground in summer. Many recent studies explored the thermodynamic aspects of the process and the thermal efficiency of the system (Cecinato & Loveridge, 2015; Cecinato & Salciarini, 2022), while others paid attention to the geotechnical and structural consequences associated with the cooling and heating processes of the piles, isolated or in groups (Di Donna et al., 2016; Salciarini et al., 2017; Bourne-Webb et al., 2021). Some experimental evidences obtained from instrumented real-scale energy piles subject to mechanical and thermal loading paths (Laloui et al., 2006; Bourne-Webb et al., 2009) indicate that even modest variations in temperature between the pile and the surrounding soil are able to cause significant changes in the state of stress and deformation in the pile and in the soil, particularly when the pile is embedded in a deep rigid layer. The goal of this work is to investigate the effect of using different constitutive models for the soil to reproduce the complex deformation phenomena caused by temperature variations in the vicinity of energy piles, using a FE-THM model implemented in Plaxis 2D. Experimental data collected by Bourne-Webb et al. (2009) were used to evaluate the forecasting capabilities of the models. 2. THE BENCHMARK CASE STUDY The project for a new building of Lambeth College in London (UK), founded on 143 piles, all equipped with U-shaped pipes that allow them to be used as heat exchangers, provided the benchmark case study for the numerical modelling comparison. The experimental study, led by Bourne-Webb et al. (2009), focused on the evaluation of the impact of thermal cycles on the performance of foundation piles and on the analysis of the THM response of the system. The stratigraphy of the site is shown in Figure 1a. The first layer (from the ground surface up to 1.5 m) consists of backfill, followed by terraced alluvial deposits (from 1.5 to 4.0 m), while the underlying layer consists of London clay. The groundwater table is located at a depth of about 3 m from the ground surface. The test pile XI Incontro Annuale dei Giovani Ingegneri Geotecnici. IAGIG 2022. Reggio Calabria, 12 Luglio 2022 © 2022 Associazione Geotecnica Italiana, Roma, Italia, ISBN 978-88-97517-11-5 83