THREE-DIMENSIONAL COUPLED FEM MODELLING AND PROGRAMMING OF PARTIALLY SATURATED POROUS MEDIA

The purpose of the work presented in this thesis is to investigate the fully and partially saturated behaviour of soils, behaviour that can be extended also to geomaterials like concrete. The physical - mathematical approach proposed within this manuscript is a coupled thermo-hydro-mechanical model, suitable for consolidation / subsidence analyses of unsaturated soils. This coupled formulation, can therefore be qualified as u – pw – pg – (T), by the introduction of basic state variables involved in the processes, that here are: the displacements field u, the liquid (water) pressure field pw, the gas (dry air and water vapour) pressure field pg, and eventually the temperature T that is involved on the modelling of non – isothermal process. Due to the coexistence of two different fluid phases, liquid and gaseous one, this model can be regarded as a multiphase approach to a deforming porous medium as proposed by Lewis and Schrefler in the framework of the hybrid mixture theory for porous media firstly presented by Hassanizadeh and Gray and Zienkiewicz et al. The evolution at macroscopic scale of the state variables above mentioned, in particular of pressures of both liquid and gas, is basically influenced by the microstructure of the material that characterizes the behaviour of a soil with relation on capillary effects and deformability. The physical approach proposed here is based on averaging techniques applied to the physical quantities that can be estimated in a representative elementary volume (REV) . With the addition of water retention functions that provide a description of the relation that exists among capillary pressure and the degree of water saturation, a complete set of fluid balance equations and mechanical and thermodynamic equilibrium equations can be obtained for the medium in a macroscopic scale. A coupled (thermo)-hydro-mechanical formulation u – p – (T) that deals with a fully saturated porous medium has been implemented with success in the past in the F.E. two-dimensional program PLASCON and its further extensions to three- dimensionality with PLASCON3D. The present work focused on the extension and upgrading of the relative simple single phase theory along with its numerical implementations, towards a more realistic multiphase description of the porous material, where voids may be filled up with both liquid and gas that interacts each other by mean of the concept of capillary pressure. An improved code PLASCON3D_PS based on the fully coupled u – pw – pg – (T) formulation and developed from previous versions has been realized. Due to the lack in literature of three-dimensional coupled numerical and experimental tests, some numerical results of benchmark tests and real case problems, that derive from two-dimensional domains, will be presented.

Lo scopo del lavoro presentato con questa tesi è di studiare il comportamento dei terreni in stato di totale e parziale saturazione, comportamento che può essere anche esteso ad altri geomateriali come il calcestruzzo. Il modello fisico - matematico proposto in questo lavoro è denominato termo- idro- meccanico in formulazione accoppiata, adatto ad analisi di consolidazione / subsidenza per suoli non saturi. La formulazione accoppiata può essere quindi qualificata come u – pw – pg – (T), attraverso l’introduzione delle variabili di stato, che qui sono: il campo degli spostamenti u, il campo di pressione liquida (acqua) pw e quello di pressione del gas (aria secca e vapore acqueo) pg ed eventualmente la temperatura T coinvolta nella modellazione di processi non isotermi. Data la coesistenza di due differenti fasi fluide, liquida e gassosa, il modello può essere visto come un approccio multifase al mezzo poroso deformabile, come già proposto da Lewis and Schrefler nel contesto della teoria delle miscele ibride per mezzo poroso, inizialmente presentata da Hassanizadeh and Gray e Zienkiewicz et al. L’evoluzione a scala macroscopica delle variabili di stato menzionate sopra, in particolare delle pressioni del liquido e del gas, sono fondamentalmente influenzate dalla microstruttura del materiale che caratterizza il comportamento del suolo in relazione agli effetti capillari e alla deformabilità. L’approccio fisico proposto qui è basato sull’uso di tecniche di media applicate alle grandezze fisiche stimate su un volume rappresentativo elementare (REV). Con l’aggiunta di funzioni di ritenzione dell’acqua che forniscono una descrizione della relazione che esiste tra pressione capillare e grado di saturazione dell’acqua, si può ottenere un set completo di equazioni di bilancio di massa per i fluidi e di equilibrio meccanico e termodinamico per il mezzo a scala macroscopica. Una formulazione accoppiata termo-idro-meccanica u – p – (T) che tratta mezzi porosi completamente saturi è stata implementata con successo in passato nel codice bidimensionale a elementi finiti PLASCON e le successive estensioni alla tridimensionalità con PLASCON3D. Il presente lavoro si focalizza sull’estensione e aggiornamento della relativamente semplice teoria a singola fase, con relative implementazioni numeriche, verso una più realistica descrizione multifase del materiale poroso, dove i vuoti possono essere riempiti da liquido e gas che interagiscono tra loro attraverso il concetto di pressione capillare. Sulla base delle versioni precedenti è stato realizzato un codice chiamato PLASCON3D_PS che implementa la formulazione accoppiata u – pw – pg – (T). Data la mancanza in letteratura di test sperimentali e numerici su domini tridimensionali, verranno presentati alcuni risultati numerici di benchmark test e casi reali derivanti da domini bidimensionali.