Seismic Response Assessment of RC Wall Structures. Applications of Non-Linear Finite Element Modelling Techniques.

Abstract (in Italian) Il problema del consumo di risorse energetiche è divenuto oggi di grandissima importanza in relazione anche al fenomeno del riscaldamento globale legato all’inquinamento provocato da centrali a combustibile fossile. Nel mondo, oltre alle centrali a combustibile tradizionale, gli impianti nucleari sono diffusi in trenta nazioni. In questo caso il problema energetico è correlato anche al più importante problema della sicurezza delle popolazioni nei confronti del rischio di fuoriuscita di radiazioni e quindi al problema dello sviluppo di metodi razionali e sicuri per il progetto di nuove strutture e per la valutazione della sicurezza di quelle esistenti. Infatti una considerevole quantità di energia viene prodotta ogni anno da centrali nucleari che necessitano di opere di rinnovamento a breve termine (Hsu et al. 2014) ed in particolare il 27% dell’energia elettrica europea viene prodotto in strutture nucleari. I componenti in calcestruzzo rinforzato delle centrali nucleari (come ad esempio vessel di contenimento e sistemi strutturali a pareti) sono di estrema importanza per la sicurezza e per l’operatività delle centrali stesse. Tali componenti in calcestruzzo armato devono essere verificati nei confronti di elevate sollecitazioni dovute all’azione di forti eventi sismici, impatti ed esplosioni. Le strutture, se sottoposte a tali azioni, mostrano un forte comportamento non lineare con conseguente apertura di fessure. È dunque di fondamentale importanza poter prevedere correttamente il comportamento tridimensionale di queste strutture in modo tale da poter valutarne sia la risposta globale (ad esempio in termini di spostamenti) sia locale (fessurazione multi - direzionale). Appare dunque chiara la necessità di sviluppare strumenti numerici adeguati ed affidabili modelli per i materiali al fine di descrivere con buona approssimazione tutti questi aspetti così complessi e difficili da cogliere con una modellazione semplificata equivalente. La validazione di tali procedure con prove sperimentali costituisce un aspetto chiave del problema. La tesi affronta dunque il problema della modellazione in campo non lineare di strutture a pareti attraverso lo sviluppo di modelli calibrati su recenti ed innovative campagne sperimentali internazionali. Attualmente sono a disposizione degli ingegneri civili numerosi strumenti numerici per la modellazione di strutture a pareti di impianti nucleari, (Richard et al. 2013). La tesi si concentra dunque sulla modellazione del comportamento dinamico non lineare di un provino in scala ¼ di un edificio di servizio per l’alloggiamento di apparecchiature elettriche per reattori nucleari CPY testato nell’ambito della blind prediction organizzata all’interno del benchmark internazionale SMART2013, (www.smart2013.eu). Inoltre, è stato analizzato e modellato il comportamento ciclico di un vessel di contenimento in C.A. per l’alloggiamento di reattori nucleari testato in scala 1/13 presso il laboratorio NCREE di Taipei (Taiwan). La modellazione effettuate con elementi shell multi-layered ed il modello fessurativo tipo “smeared” PARC_CL. Nelle analisi il modello PARC_CL è stato utilizzato per la valutazione delle matrice di rigidezza non lineare nei punti di integrazione degli elementi shell. Il modello è implementato nella user subroutine UMAT.for del software commerciale ABAQUS. Tale approccio è stato applicato in quanto è stato osservato che la modellazione con elementi shell e fessurazione diffusa può predire adeguatamente gli indicatori di danno. È possibile in particolare stimare gli EDP (Engineering Demand Parameter) sia locali (aperture di fessura, deformazione nelle barre) che globali (spostamenti di interpiano), parametri fondamentali per la progettazione sismica definita “performance-based”. La modellazione con elementi shell e modello fessurativo PARC_CL si è rivelata adatta per la valutazione del comportamento dinamico non lineare e di entrambe le strutture presentate nei casi studio. Nel primo capitolo a carattere introduttivo sono descritte le principali problematiche che sono alla base dello studio e la metodologia utilizzate nel corso della ricerca. Obiettivo della ricerca è la valutazione delle performance sismiche di strutture a pareti (in particolare edifici nucleari) attraverso modellazione agli elementi finiti di veri test di laboratorio. I capitoli da 2 a 4 sono raggruppati in una parte prima relativa agli strumenti di indagine di maggiore interesse utilizzati nel corso della ricerca. Il primo di tali aspetti è la modellazione del materiale calcestruzzo armato attraverso il modello non lineare a fessurazione diffusa sviluppato per carichi ciclici PARC_CL. Il modello PARC (Phyisical Approach on Reinforced Concrete) è un modello a fessurazione ad inclinazione fissa ed armatura “smeared” (diffusa). Il modello è particolarmente efficiente nella valutazione nella risposta di strutture a pareti in quanto calibrato su problemi dominati da comportamento a taglio. Tali problemi sono infatti tipici di strutture strategiche quali strutture off shore, nuclei di contenimento, sistemi di pareti sismo resistenti. Il modello inoltre è funzionale alla valutazione delle performance locali della struttura in quanto in grado di valutare importanti parametri quali l’apertura di fessure nel calcestruzzo (parametro particolarmente importante per strutture di contenimento, particolarmente impiegate nell’ambito della progettazione nucleare, in quanto in questo la struttura deve anche prevenire la fuoriuscita di radiazioni o altri agenti contaminanti in caso di sollecitazioni elevate, dovute non solo ad eventi sismici ma anche ad esempio impatti ed esplosioni). Inoltre è possibile valutare tensioni e deformazioni nel calcestruzzo e nelle armature. Il modello è applicato alla modellazione con elementi shell multistrato. Dal punto di vista globale è dunque in grado di cogliere anche il comportamento tridimensionale di strutture complesse, come dimostrato nel capitolo 6 della tesi. Grazie all’implementazione del comportamento ciclico dei legami costitutivi di calcestruzzo ed acciaio il modello è inoltre in grado di tenere conto in maniera realistica della dissipazione dell’energia per isteresi, aspetto particolarmente delicato nella valutazione della risposta dinamica della struttura. Il capitolo 3 costituisce una sintesi dell’approccio probabilistico all’analisi strutturale. In particolare, è presentato e descritto attraverso esempi applicativi un metodo probabilistico utilizzato in ambito di progettazione di strutture nucleari per la valutazione di curve di capacità in funzione di diverse misure di intensità (intensity measures) e parametri di valutazione della domanda (engineer demand parameters), che costituiscono un importante aspetto dei moderni approcci di verifica e progettazione impiegati non solo in ambito accademico. Tale ricerca è particolarmente importante ai fini di valutare la robustezza e l’affidabilità degli strumenti numerici impiegati nell’ingegneria strutturale. Tale procedura è stata applicata al fine di valutare le curve di fragilità nel capitolo 6. Il capitolo 4 costituisce una breve review di importanti aspetti relativi all’utilizzo di accelerogrammi e alla generazione di spettri di risposta affrontati durante la ricerca. In particolare, è affrontato nel dettaglio il tema dell’utilizzo del cosiddetto approccio energetico per la realizzazione di spettri di risposta ai fini ad esempio per la progettazione di organi dissipativi ad attrito. Inoltre viene discussa la problematica dell’impiego di accelero grammi reali in analisi dinamiche al passo e sono forniti esempi pratici relativi a campionamento, filtro di frequenze, correzione della base-line. Infine un paragrafo è dedicato alla generazione di spettri di Fourier di segnali sismici, aspetto fondamentale dell’identificazione strutturale. In quella che costituisce la seconda parte della tesi è presentata una review relativa all’utilizzo del modello PARC per la modellazione di strutture a pareti e pareti singole, prove sperimentali su tali strutture disponibili in letteratura ed è posta attenzione anche sulle procedure di calcolo analitiche disponibili per la determinazione della resistenza degli elementi a parete nei confronti delle diverse modalità di rottura. La terza parte della tesi raggruppa in sé i capitoli 6, 7 ed 8 e costituisce il nucleo di maggiore interesse nel corso della trattazione. Nel capitolo 6 è descritta nel dettaglio la modellazione di una struttura a pareti con forte comportamento torsionale. Tale struttura è stata testata in scala ¼ presso i laboratori di Saclay (Francia) su tavola vibrante con l’applicazione di una sequenza di accelerogrammi di intensità crescente fino a 1.6g. La risposta non lineare della struttura è stata prima analizzata “in cieco” ovvero senza conoscere i risultati della prova sperimentale ma avendo noti solo gli input in termini di spostamento e accelerazione e le caratteristiche dei materiali. Il modello sviluppato con elementi shell e PARC_CL è stato in grado di cogliere correttamente il campo di spostamenti e i fenomeni locali di fessurazione, taglio, snervamento delle barre d’armatura. Un ulteriore modello è stato realizzato in fase di post-analisi per valutare l’influenza di alcune variabili (introduzione della tavola vibrante nel modello globale, modellazione delle masse addizionali, definizione dell’input sismico, calibrazione del damping alla Rayleigh). Tuttavia il modello non è stato in grado di cogliere la rottura per crushing del calcestruzzo osservata in corrispondenza dell’interfaccia tra fondazione ed una delle pareti che costituivano il sistema sismo resistente. Nel capitolo 7 è invece descritta la modellazione di un vessel di contenimento per reattori nucleari testato per carico ciclico quasi - statico in scala 1/13 presso il laboratorio NCREE di Taipei (Taiwan). Il modello implementato è stato in grado di cogliere il comportamento si in fase lineare che fessurata del vessel. Inoltre localmente sono state stimate tensioni e deformazioni nelle barre d’armatura. Fenomeni di fessurazione sono stati osservati sia numericamente che sperimentalmente e globalmente è stato possibile stimare la forza di picco raggiunta in fase sperimentale. Durante il test è stata identificata una modalità di rottura per taglio scorrimento all’interfaccia tra vessel e fondazione. Tale osservazione ha portato alla conclusione della necessità di sviluppare un modello locale da implementare nel codice di calcolo per tenere in considerazione appunto dal comportamento tra interfacce in calcestruzzo armato. Nell’ottavo capitolo è dunque riportata come estensione dello studio lo sviluppo e la calibrazione preliminare di una tecnica di modellazione locale in grado di cogliere le rotture di tipo taglio scorrimento. La filosofia alla base di tale modellazione è la stessa alla base delle formulazioni analitiche utilizzate in fase di verifica che valutano la resistenza a scorrimento come somma di due contributi principali per la trasmissione del taglio. La trasmissione del taglio per attrito viene considerata attraverso opportuna modellazione e definizione delle proprietà delle superfici di contatto tra parete e fondazione (dunque considerando incidenze diverse per gli elementi di fondazione e gli elementi che definiscono la parete). Il contributo dato dalla dowel action viene introdotto invece attraverso lo sviluppo di un elemento multi-spring spaziale in grado di tenere conto anche della risposta ciclica delle barre d’armatura attraverso l’implementazione di un legame sperimentale taglio-scorrimento. Per concludere, nella moderna pratica progettuale una delle principali problematiche è la valutazione degli “engineering demand parmaters” in modo da stimare correttamente le performance della struttura allo stato limite di danno (soprattutto per garantire l’operatività degli impianti e degli elementi non strutturali) ed allo stato limite di collasso. Questo aspetto risulta molto importante se si tiene in considerazione il fatto che le pareti in C.A. costituiscono spesso il sistema resistente di edifici strategici. Recenti programmi sperimentali su strutture a parete, in particolare il programma SMART2013 (Saclay, Francia), www.smart2013.eu, e la campagna di prove su una struttura cilindrica di contenimento presso il laboratorio NCREE di Taipei (Taiwan), hanno evidenziato che attualmente esistono ancora incertezze che dovrebbero essere considerate nella valutazione della sicurezza sismica delle strutture. In particolare entrambe le prove sono state condizionate da una modalità di rottura per taglio scorrimento. I modelli proposti nella presente tesi si sono dimostrati efficienti nella valutazione delle performance locali e globali di strutture a parete in C.A. di edifici a pareti in calcestruzzo rinforzato ed in particolare di centrali nucleari. Tuttavia, ulteriori studi saranno portati avanti a partire dai risultati ottenuti per l’implementazione e la calibrazione di un modello numerico in grado ci cogliere le rotture per taglio scorrimento negli elementi a pareti.