Previsione del comportamento dinamico e sismico di un edificio a tre piani in scala reale costituito da pareti sandwich in C.A. gettato in opera da provare su tavola vibrante

INDICE INTRODUZIONE 1 1. DESCRIZIONE DEL SISTEMA COSTRUTTIVO 5 1.1 I pannelli modulari 5 1.2 Le pareti tozze in cemento armato gettate in opera realizzate con la tecnologia del pannello di supporto in polistirene 5 1.3 La connessione tra le pareti e la fondazione 6 1.4 Le connessioni tra pareti ortogonali 7 1.5 Le connessioni tra pareti e solai 7 1.6 Il sistema strutturale così ottenuto e le sue caratteristiche salienti 8 2. RICERCA BIBLIOGRAFICA 11 2.1 Pareti tozze e pareti snelle 11 2.2 Il comportamento scatolare 13 2.3 I muri sandwich 14 2.4 Il “ferro-cemento” 15 3. DATI DI PARTENZA 19 3.1 Schema geometrico - architettonico definitivo 19 3.2 Abaco delle sezioni e delle armature 21 3.3 Materiali e resistenze 22 3.4 Valutazione del momento di inerzia delle pareti estese debolmente armate 23 3.4.1 Generalità 23 3.4.2 Caratteristiche degli elementi provati 23 3.4.3 Formulazioni analitiche 23 3.4.4 Considerazioni sulla deformabilità dei pannelli debolmente armati 24 3.4.5 Confronto tra rigidezze sperimentali e rigidezze valutate analiticamente 26 3.4.6 Stima di un modulo elastico equivalente 26 4. ANALISI DEI CARICHI 29 4.1 Stima dei carichi di progetto della struttura 29 4.1.1 Stima dei pesi di piano 30 4.1.2 Tabella riassuntiva dei pesi di piano 31 4.2 Analisi dei carichi da applicare in fase di prova 32 4.2.1 Pesi di piano 34 4.2.2 Tabella riassuntiva dei pesi di piano 35 4.3 Pesi della struttura 36 4.3.1 Ripartizione del carico sulle pareti parallele e ortogonali 36 5. DESCRIZIONE DEL MODELLO AGLI ELEMENTI FINITI 37 5.1 Caratteristiche di modellazione 37 5.2 Caratteristiche geometriche del modello 38 5.3 Analisi dei carichi 41 5.4 Modello con shell costituite da un solo layer 43 5.4.1 Modellazione dei solai 43 5.4.2 Modellazione delle pareti 44 5.4.3 Descrizione delle caratteristiche dei materiali 46 5.4.3.1 Comportamento lineare dei materiali 46 6. ANALISI DEL COMPORTAMENTO STATICO DELLA STRUTTURA 49 6.1 Azioni statiche 49 6.2 Analisi statica 49 7. ANALISI DEL COMPORTAMENTO DINAMICO DELLA STRUTTURA 51 7.1 Determinazione del periodo proprio della struttura con il modello FEM 51 7.1.1 Modi di vibrare corrispondenti al modello con solai e pareti costituiti da elementi shell 51 7.1.1.1 Modi di vibrare con modulo pari a E 51 7.1.1.2 Modi di vibrare con modulo pari a 0,5E 51 7.1.1.3 Modi di vibrare con modulo pari a 0,1E 51 7.1.2 Modi di vibrare corrispondenti al modello con solai infinitamente rigidi e pareti costituite da elementi shell 52 7.1.2.1 Modi di vibrare con modulo pari a E 52 7.1.2.2 Modi di vibrare con modulo pari a 0,5E 52 7.1.2.3 Modi di vibrare con modulo pari a 0,1E: 52 7.1.3 Modi di vibrare corrispondenti al modello con solai irrigiditi con bielle e pareti costituite da elementi shell 53 7.1.3.1 Modi di vibrare con modulo pari a E 53 7.1.3.2 Modi di vibrare con modulo pari a 0,5E 53 7.1.3.3 Modi di vibrare con modulo pari a 0,1E 53 7.2 Calcolo del periodo proprio della struttura assimilandola ad un oscillatore semplice 59 7.2.1 Analisi svolta assumendo l’azione del sisma in ingresso in direzione X-X 59 7.2.1.1 Analisi svolta assumendo il modulo elastico E pari a 300000 Kg/cm2 59 7.2.1.1.1 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari ad E 59 7.2.1.1.2 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari ad E 61 7.2.1.1.3 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari ad E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 63 7.2.1.1.4 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari ad E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 66 7.2.1.2 Analisi svolta assumendo il modulo elastico E pari a 150000 Kg/cm2 69 7.2.1.2.1 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari a 0,5E 69 7.2.1.2.2 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari a 0,5E 71 7.2.1.2.3 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari a 0,5 E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 73 7.2.1.2.4 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari a 0,5 E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 76 7.2.1.3 Analisi svolta assumendo il modulo elastico E pari a 30000 Kg/cm2 79 7.2.1.3.1 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari a 0,1E 79 7.2.1.3.2 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari a 0,1E 81 7.2.1.3.3 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari a 0,1E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 83 7.2.1.3.4 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari a 0,1E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 86 7.2.2 Analisi svolta assumendo l’azione del sisma in ingresso in direzione Y-Y 89 7.2.2.1 Analisi svolta assumendo il modulo elastico E pari a 300000 Kg/cm2 89 7.2.2.1.1 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari ad E 89 7.2.2.1.2 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari ad E 91 7.2.2.1.3 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari ad E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 93 7.2.2.1.4 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari ad E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 98 7.2.2.1.5 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari ad E 103 7.2.2.1.6 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari ad E 105 7.2.2.1.7 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari ad E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 107 7.2.2.1.8 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari ad E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 112 7.2.2.2 Analisi svolta assumendo il modulo elastico E pari a 150000 Kg/cm2 117 7.2.2.2.1 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari a 0,5E 117 7.2.2.2.2 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari a 0,5E 119 7.2.2.2.3 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari a 0,5 E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 121 7.2.2.2.4 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari a 0,5 E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 126 7.2.2.2.5 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari a 0,5 E 131 7.2.2.2.6 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari ad E 133 7.2.2.2.7 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari a 0,5E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 135 7.2.2.2.8 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari a 0,5E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 140 7.2.2.3 Analisi svolta assumendo il modulo elastico E pari a 30000 Kg/cm2 145 7.2.2.3.1 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari a 0,1E 145 7.2.2.3.2 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari a 0,1E 147 7.2.2.3.3 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari a 0,1E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 149 7.2.2.3.4 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari a 0,1E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 154 7.2.2.3.5 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H e modulo elastico assunto pari a 0,1 E 159 7.2.2.3.6 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H e modulo elastico assunto pari ad E 161 7.2.2.3.7 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 2/3 H, modulo elastico assunto pari a 0,1E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 163 7.2.2.3.8 Determinazione del periodo proprio della struttura considerando la massa complessiva concentrata a 1/2 H, modulo elastico assunto pari a 0,1E, e struttura resistente costituita dai soli “maschi murari” delle pareti parallele all’azione del sisma 168 7.3 Calcolo del periodo proprio della struttura approssimato utilizzando espressioni analitiche 174 7.3.1 Approssimazione della struttura ad una mensola incastrata di peso Q=ql avente un peso P gravante all’estremo libero 174 7.3.1.1 Riferimenti teorici: sostituzione di masse distribuite con masse concentrate 174 7.3.1.2 Applicazione allo specifico caso di studio in esame con modulo elastico E=300000 kg/cm2 177 7.3.1.3 Applicazione allo specifico caso di studio in esame con modulo elastico E=30000 kg/cm2 179 7.3.2 Approssimazione della struttura ad una mensola incastrata alla base, di peso Q=ql, avente un peso P gravante all’estremo libero e struttura resistente costituita dai soli “maschi murari”delle pareti parallele all’azione del sisma 181 7.3.2.1 Applicazione allo specifico caso di studio in esame con modulo elastico E=300000 kg/cm2 181 7.3.2.2 Applicazione allo specifico caso di studio in esame con modulo elastico E=30000 kg/cm2 186 7.3.3 Approssimazione della struttura ad un portale avente peso Qp = peso di un piedritto, Qt=peso del traverso e un peso P gravante sul traverso medesimo 191 7.3.3.1 Riferimenti teorici: sostituzione di masse distribuite con masse concentrate 191 7.3.3.2 Applicazione allo specifico caso di studio in esame con modulo ellastico E=300000 kg/cm2 192 7.3.3.3 Applicazione allo specifico caso di studio in esame con modulo ellastico E=30000 kg/cm2 194 7.3.4 Approssimazione della struttura ad un portale di peso Qp = peso di un piedritto, Qt=peso del traverso e avente un peso P gravante sul traverso medesimo e struttura resistente costituita dai soli “maschi murari”delle pareti parallele all’azione del sisma 196 7.3.4.1 Applicazione allo specifico caso di studio in esame con modulo elastico E=300000 kg/cm2 196 7.3.4.2 Applicazione allo specifico caso di studio in esame con modulo elastico E=30000 kg/cm2 201 7.3.5 Approssimazione della struttura ad una mensola incastrata di peso Q=ql avente le masse m1,m2....mn concentrate nei punti 1,2….n 206 7.3.5.1 Riferimenti teorici: metodo approssimato 206 7.3.5.2 Applicazione allo specifico caso di studio in esame con modulo elastico E=300000 kg/cm2 207 7.3.5.3 Applicazione allo specifico caso di studio in esame con modulo elastico E=30000 kg/cm2 209 7.3.6 Approssimazione della struttura ad un telaio deformabile con tavi infinitamente rigide 211 7.3.6.1 Riferimenti teorici: vibrazioni dei telai 211 7.3.6.2 Applicazione allo specifico caso di studio in esame con modulo elastico E=300000 kg/cm2 212 7.3.6.3 Applicazione allo specifico caso di studio in esame con modulo elastico E=30000 kg/cm2 215 7.3.7 Approssimazione della struttura ad una mensola incastrata di peso Q=ql avente masse m1,m2....mn concentrate nei punti 1,2….n e studiata come un sistema continuo 218 7.3.7.1 Riferimenti teorici: metodo energetico; Masse ripartite e concentrate; Formula di Dunkerley 218 7.3.7.1.1 Il metodo energetico 218 7.3.7.1.2 Masse ripartite e concentrate. Formula di Dunkerley 219 7.3.7.2 Applicazione allo specifico caso di studio in esame con modulo elastico E=300000 kg/cm2 221 7.3.7.3 Applicazione allo specifico caso di studio in esame con modulo elastico E=30000 kg/cm2 226 7.4 Calcolo del periodo della struttura approssimato mediante telaio equivalente 232 7.4.1 Dati geometrici relativi al telaio equivalente e determinazione dei carichi agenti su di esso 232 7.4.1.1 Determinazione del periodo proprio della struttura assumendo diversi valori del modulo elastico E 233 7.5 Conclusioni 234 7.5.1 Comparazione dei risultati relativi alla schematizzazione dell’edificio con una struttura ad un grado di libertà 234 7.5.2 Comparazione dei risultati relativi alla schematizzazione dell’edificio con una struttura a più gradi di libertà e a sistema continuo 236 8. ANALISI DEL COMPORTAMENTO SISMICO DELLA STRUTTURA 239 8.1 Modello con shell costituite da un solo layer 239 8.1.1 Analisi dinamica modale con spettro di risposta avente un valore di PGA pari a 0,1g 239 8.1.1.1 Generalità 239 8.1.1.2 Sollecitazioni e tensioni sulla sezione di base 242 8.1.1.2.1 Combinazione di carico ”Carichi verticali più Spettro di Risposta scalato ad un valore di PGA pari a 0,1g” 242 8.1.1.2.2 Combinazione di carico ”Spettro di Risposta scalato ad un valore di 0,1g di PGA” 245 8.1.1.3 Spostamenti di piano 248 8.1.1.4 Accelerazioni di piano 248 8.1.2 Analisi Time-History lineare con accelerogramma caratterizzato da un valore di PGA pari a 0,1g 249 8.1.2.1 Generalità 249 8.1.2.2 Sollecitazioni e tensioni sulla sezione di base 251 8.1.2.2.1 Combinazione di carico ” Carichi verticali più Accelerogramma agente in direzione Ye avente una PGA pari a 0,1g” 251 8.1.2.2.2 Combinazione di carico ” Accelerogramma agente in direzione Y avente un valore di PGA pari a 0,1g ” 254 8.1.2.3 Spostamenti di piano assoluti 257 8.1.2.4 Spostamenti di piano relativi 260 8.1.2.5 Accelerazioni di piano assolute 262 8.1.3 Analisi dinamica modale con spettro di risposta avente un valore di PGA pari a 0,3g 264 8.1.3.1 Generalità 264 8.1.3.2 Sollecitazioni e tensioni sulla sezione di base 265 8.1.

Efficienza energetica dei sistemi acquedottistici: ruolo delle perdite idriche

Il lavoro di ricerca esplora il panorama dell’efficienza energetica dei sistemi acquedottistici, soffermandosi a considerare possibili indicatori che possano valutarla in maniera corretta e completa, in particolare nei confronti della presenza di perdite in rete. Si prendono in considerazione con maggiore attenzione, tra tutte le strategie per aumentare l’efficienza energetica, quelle che contemporaneamente producono anche risparmi idrici, come la riduzione della pressione e la ricerca attiva delle perdite . Dopo un inquadramento internazionale, sono stati analizzati mediante mappe tematiche di intensità energetica, i consumi energetici specifici sui sistemi acquedottistici della regione Emilia Romagna per gli anni 2006 e 2007, si è passati ad una analisi critica degli indicatori attualmente in uso. Inoltre per casi di studio sintetici e tutti i casi di studio proposti, si sono valutate curve di relazione tra percentuale di perdita idrica e aumento del consumo energetico, in grado di dare indicazioni su come ciascun sistema reagisce, in termini di aumento dell’energia consumata, all’aumentare del livello di perdita. Questa relazione appare fortemente influenzata da fattori come la modalità di pompaggio, la posizione delle rotture sulla rete e la scabrezza delle condotte. E’ emersa la necessità solo poter analizzare separatamentel’influenza sull’efficienza energeticadei sistemi di pompaggio e della rete, mostrando il ruolo importante con cui questa contribuisce all’efficienza globale del sistema. Viene proposto uno sviluppo ulteriore dell’indicatore GEE Global Energy Efficiency (Abadia, 2008), che consente di distinguere l’impatto sull’efficienza energetica dovuto alle perdite idriche e alla struttura intrinseca della rete, in termini di collocazione reciproca tra risorsa idrica e domanda e schema impiantistico.Questa metodologia di analisi dell’efficienza energetica è stata applicata ai casi di studio, sia sintetici che reali, il distretto di Marzaglia (MO) e quello di Mirabello (FE), entrambi alimentati da pompe a giri variabili.. La ricerca ha consentito di mostrare inoltre il ruolo della modellazione numerica in particolare nell’analisi dell’effetto prodotto sull’efficienza energetica dalla presenza di perdite idriche. Nell’ultimo capitolo si completa la panoramica dei benefici ottenibili attraverso la riduzione della pressione, che nei casi citati viene conseguita tramite pompe asservite ad inverter, con il caso di studio del distretto Bolognina all’interno del sistema di distribuzione di Bologna, che vede l’utilizzo di valvole riduttrici di pressione. Oltre a stimare il risparmio energetico derivante dalla riduzione delle perdite ottenuta tramite le PRV, sono stati valutati su modello i benefici energetici conseguenti all’introduzione nel distretto di turbine per la produzione di energia

Shear behavior of reinforced concrete slabs under concentrated load: an investigation through Sequentially Linear Analysis

ABSTRACT (italiano) Con crescente attenzione riguardo al problema della sicurezza di ponti e viadotti esistenti nei Paesi Bassi, lo scopo della presente tesi è quello di studiare, mediante la modellazione con Elementi Finiti ed il continuo confronto con risultati sperimentali, la risposta in esercizio di elementi che compongono infrastrutture del genere, ovvero lastre in calcestruzzo armato sollecitate da carichi concentrati. Tali elementi sono caratterizzati da un comportamento ed una crisi per taglio, la cui modellazione è, da un punto di vista computazionale, una sfida piuttosto ardua, a causa del loro comportamento fragile combinato a vari effetti tridimensionali. La tesi è incentrata sull'utilizzo della Sequentially Linear Analysis (SLA), un metodo di soluzione agli Elementi Finiti alternativo rispetto ai classici approcci incrementali e iterativi. Il vantaggio della SLA è quello di evitare i ben noti problemi di convergenza tipici delle analisi non lineari, specificando direttamente l'incremento di danno sull'elemento finito, attraverso la riduzione di rigidezze e resistenze nel particolare elemento finito, invece dell'incremento di carico o di spostamento. Il confronto tra i risultati di due prove di laboratorio su lastre in calcestruzzo armato e quelli della SLA ha dimostrato in entrambi i casi la robustezza del metodo, in termini di accuratezza dei diagrammi carico-spostamento, di distribuzione di tensioni e deformazioni e di rappresentazione del quadro fessurativo e dei meccanismi di crisi per taglio. Diverse variazioni dei più importanti parametri del modello sono state eseguite, evidenziando la forte incidenza sulle soluzioni dell'energia di frattura e del modello scelto per la riduzione del modulo elastico trasversale. Infine è stato effettuato un paragone tra la SLA ed il metodo non lineare di Newton-Raphson, il quale mostra la maggiore affidabilità della SLA nella valutazione di carichi e spostamenti ultimi insieme ad una significativa riduzione dei tempi computazionali. ABSTRACT (english) With increasing attention to the assessment of safety in existing dutch bridges and viaducts, the aim of the present thesis is to study, through the Finite Element modeling method and the continuous comparison with experimental results, the real response of elements that compose these infrastructures, i.e. reinforced concrete slabs subjected to concentrated loads. These elements are characterized by shear behavior and crisis, whose modeling is, from a computational point of view, a hard challenge, due to their brittle behavior combined with various 3D effects. The thesis is focused on the use of Sequentially Linear Analysis (SLA), an alternative solution technique to classical non linear Finite Element analyses that are based on incremental and iterative approaches. The advantage of SLA is to avoid the well-known convergence problems of non linear analyses by directly specifying a damage increment, in terms of a reduction of stiffness and strength in the particular finite element, instead of a load or displacement increment. The comparison between the results of two laboratory tests on reinforced concrete slabs and those obtained by SLA has shown in both the cases the robustness of the method, in terms of accuracy of load-displacements diagrams, of the distribution of stress and strain and of the representation of the cracking pattern and of the shear failure mechanisms. Different variations of the most important parameters have been performed, pointing out the strong incidence on the solutions of the fracture energy and of the chosen shear retention model. At last a confrontation between SLA and the non linear Newton-Raphson method has been executed, showing the better reliability of the SLA in the evaluation of the ultimate loads and displacements, together with a significant reduction of computational times.

Sviluppo di modelli motore e veicolo per l'analisi di strategie di controllo in applicazioni Software e Hardware In the Loop

This work describes the development of a simulation tool which allows the simulation of the Internal Combustion Engine (ICE), the transmission and the vehicle dynamics. It is a control oriented simulation tool, designed in order to perform both off-line (Software In the Loop) and on-line (Hardware In the Loop) simulation. In the first case the simulation tool can be used in order to optimize Engine Control Unit strategies (as far as regard, for example, the fuel consumption or the performance of the engine), while in the second case it can be used in order to test the control system. In recent years the use of HIL simulations has proved to be very useful in developing and testing of control systems. Hardware In the Loop simulation is a technology where the actual vehicles, engines or other components are replaced by a real time simulation, based on a mathematical model and running in a real time processor. The processor reads ECU (Engine Control Unit) output signals which would normally feed the actuators and, by using mathematical models, provides the signals which would be produced by the actual sensors. The simulation tool, fully designed within Simulink, includes the possibility to simulate the only engine, the transmission and vehicle dynamics and the engine along with the vehicle and transmission dynamics, allowing in this case to evaluate the performance and the operating conditions of the Internal Combustion Engine, once it is installed on a given vehicle. Furthermore the simulation tool includes different level of complexity, since it is possible to use, for example, either a zero-dimensional or a one-dimensional model of the intake system (in this case only for off-line application, because of the higher computational effort). Given these preliminary remarks, an important goal of this work is the development of a simulation environment that can be easily adapted to different engine types (single- or multi-cylinder, four-stroke or two-stroke, diesel or gasoline) and transmission architecture without reprogramming. Also, the same simulation tool can be rapidly configured both for off-line and real-time application. The Matlab-Simulink environment has been adopted to achieve such objectives, since its graphical programming interface allows building flexible and reconfigurable models, and real-time simulation is possible with standard, off-the-shelf software and hardware platforms (such as dSPACE systems).

Il linguaggio Alloy come ausilio nella specifica formale di modelli UML

UML è ampiamente considerato lo standard de facto nella fase iniziale di modellazione di sistemi software basati sul paradigma Object-Oriented; il suo diagramma delle classi è utilizzato per la rappresentazione statica strutturale di entità e relazioni che concorrono alla definizione delle specifiche del sistema; in questa fase viene utilizzato il linguaggio OCL per esprimere vincoli semantici sugli elementi del diagramma. Il linguaggio OCL però soffre della mancanza di una verifica formale sui vincoli che sono stati definiti. Il linguaggio di modellazione Alloy, inserendosi in questa fase, concettualmente può sopperire a questa mancanza perchè può descrivere con le sue entità e relazioni un diagramma delle classi UML e, tramite propri costrutti molto vicini all'espressività di OCL, può specificare vincoli semantici sul modello che verranno analizzati dal suo ambiente l'Alloy Analyzer per verificarne la consistenza. In questo lavoro di tesi dopo aver dato una panoramica generale sui costrutti principali del linguaggio Alloy, si mostrerà come è possibile creare una corrispondenza tra un diagramma delle classi UML e un modello Alloy equivalente. Si mostreranno in seguito le analogie che vi sono tra i costrutti Alloy e OCL per la definizione di vincoli formali, e le differenze, offrendo nel complesso soluzioni e tecniche che il modellatore può utilizzare per sfruttare al meglio questo nuovo approccio di verifica formale. Verranno mostrati anche i casi di incompatibilità. Infine, come complemento al lavoro svolto verrà mostrata, una tecnica per donare una dinamicità ai modelli statici Alloy.

La Ricostruzione 3D di Bologna nel XVI Secolo

Il lavoro presentato ha come oggetto la ricostruzione tridimensionale della città di Bologna nella sua fase rinascimentale. Tale lavoro vuole fornire un modello 3D delle architetture e degli spazi urbani utilizzabile sia per scopi di ricerca nell’ambito della storia delle città sia per un uso didattico-divulgativo nel settore del turismo culturale. La base del lavoro è una fonte iconografica di grande importanza: l’affresco raffigurante Bologna risalente al 1575 e situato in Vaticano; questa è una veduta a volo d’uccello di grandi dimensioni dell’intero tessuto urbano bolognese all’interno della terza cerchia di mura. In esso sono rappresentate in maniera particolareggiata le architetture civili e ecclesiastiche, gli spazi ortivi e cortilivi interni agli isolati e alcune importanti strutture urbane presenti in città alla fine del Cinquecento, come l’area portuale e i canali interni alla città, oggi non più visibili. La ricostruzione tridimensionale è stata realizzata tramite Blender, software per la modellazione 3D opensource, attraverso le fasi di modellazione, texturing e creazione materiali (mediante campionamento delle principali cromie presenti nell’affresco), illuminazione e animazione. Una parte della modellazione è stata poi testata all’interno di un GIS per verificare l’utilizzo delle geometrie 3D come elementi collegabili ad altre fonti storiche relative allo sviluppo urbano e quindi sfruttabili per la ricerca storica. Grande attenzione infine è stata data all’uso dei modelli virtuali a scopo didattico-divulgativo e per il turismo culturale. La modellazione è stata utilizzata all’interno di un motore grafico 3D per costruire un ambiente virtuale interattivo nel quale un utente anche non esperto possa muoversi per esplorare gli spazi urbani della Bologna del Cinquecento. In ultimo è stato impostato lo sviluppo di un’applicazione per sistemi mobile (Iphone e Ipad) al fine di fornire uno strumento per la conoscenza della città storica in mobilità, attraverso la comparazione dello stato attuale con quello ricostruito virtualmente.

Analisi delle tecnologie BMP per i sistemi di drenaggio urbano

I sistemi di drenaggio urbano meritano particolare attenzione per l’importante compito che devono garantire all’interno della società. L’allontanamento delle acque reflue dalle zone urbanizzate infatti è indispensabile per non creare problemi ed interferenze con attività sociali ed economiche che in esse si svolgono; inoltre un corretto funzionamento del sistema di drenaggio assicura anche un basso impatto ambientale in termini di qualità e quantità. La disponibilità di nuove tecnologie costruttive e materiali in continua evoluzione devono fare fronte alla crescente complessità idraulica, geotecnica e urbanistica delle strutture componenti i sistemi di drenaggio; in tale elaborato si vuole porre l’accento sul problema dello smaltimento delle acque meteoriche nelle zone densamente urbanizzate, la cui soluzione è costituita dall'adozione delle cosiddette BMP (Best Management Practices). Le BMP sono definite come strategie, pratiche o metodi per la rimozione, la riduzione, il ritardo o la prevenzione della quantità di costituenti inquinanti e contaminanti delle acque di pioggia, prima che giungano nei corpi idrici ricettori. Tra esse ricordiamo le cisterne o rain barrels, i pozzi asciutti, i sistemi drywell, le vasche verdi filtranti, i tetti verdi, i canali infiltranti, i fossi di infiltrazione, i pozzi perdenti, i planting container, i canali inerbiti, i bacini di infiltrazione, le gross pullutant traps (gpts), gli stagni, i sistemi di fitodepurazione estensiva (sfe), le pavimentazioni drenanti. In Italia esse risultano ancora scarsamente studiate ed applicate, mentre trovano più ampio sviluppo in realtà estere come in Inghilterra ed in Australia. Per comprendere la efficacia di tali tecniche BMP è necessario analizzarle e, soprattutto, analizzare il loro effetto sul territorio in cui esse vengono inserite e sul regime idrogeologico dell’ambiente. Questa analisi può essere svolta con prove sperimentali su aree di controllo (soluzione economicamente gravosa) o attraverso modelli matematici tramite l’utilizzo di software di calcolo che simulino il comportamento delle portate, dei tiranti idrici e degli inquinanti nei canali e nelle strutture accessorie costituenti la rete fognaria, in conseguenza ad eventi di pioggia dei quali sia nota la distribuzione spaziale e temporale. In questo elaborato si modellizza attraverso il programma Matlab un unico elemento BMP e si osservano i suoi effetti; si procede cioè alla analisi del funzionamento di un ipotetico brown roof installato nella zona di Rimini e si osservano i benefici effetti che ne derivano in termini di volumi di pioggia trattenuti dal sistema considerando diverse tipologie di pioggia e diversi parametri progettuali per il tetto (i fori dello strato inferiore del tetto, la altezza dello strato di terreno, la sua permeabilità e la sua porosità). Si procede poi con la analisi di una ipotetica zona commerciale che sorge sulle sponde di un fiume. Tali analisi vengono svolte con il software HEC-RAS per quanto riguarda la analisi dei livelli del fiume e delle zone a rischio di inondazione. I risultati forniti da questo studio preliminare vengono poi utilizzati come condizioni al contorno per una successiva analisi effettuata con il software InfoWorks in cui si valutano i benefici che derivano dalla installazione di diverse BMP nella zona commerciale oggetto di studio. Tale modello esamina un verosimile scenario inglese, la cittadina infatti si ipotizza situata in Inghilterra e anche gli eventi pluviometrici a cui ci si riferisce sono tipici eventi di pioggia inglesi. Tutti i dati di input elaborati nelle simulazioni sono stati forniti all’interno di un progetto universitario svolto presso l’università di Bradford con la supervisione dei professori Simon Tait ed Alma Schellart. Infine la parte conclusiva dell’elaborato è stata sviluppata in collaborazione con Hera S.p.a. di Rimini all’interno di un percorso di tirocinio curriculare che ha previsto la analisi delle tecnologie BMP adatte per l’ambiente cittadino di Rimini e la osservazione delle performance garantite da tecnologie esistenti sul territorio: i parcheggi permeabili. La realtà riminese infatti deve far fronte a diversi problemi che si innescano durante i periodi di pioggia ed una mitigazione dei volumi di acqua in arrivo alla fognatura grazie a questo tipo di tecnologie, senza quindi ricorrere agli usuali sistemi puntuali, porterebbe sicuramente a notevoli vantaggi. Tali tecnologie però non devono solo essere progettate e costruite in modo corretto, ma devono essere sottoposte a periodici controlli ed adeguate operazioni di manutenzione per non perdere la loro efficacia. Per studiare tale aspetto si procede quindi alla misura della permeabilità di parcheggi drenanti presenti all’interno dei comuni di Cattolica e Rimini ricercando i fattori che influenzano tale caratteristica fondamentale dei parcheggi.

Implementazione di una nuova metodologia per la simulazione di un motore sovralimentato ad alte prestazioni

In questo lavoro verrà analizzato lo sviluppo di una nuova modellazione matematica per la simulazione della dinamica del flusso all’interno del sistema di aspirazione per motori sovralimentati. Tale modellazione si basa sulla risoluzione numerica mediante la formulazione proposta da Courant, Isaacson e Rees (CIR) nel 1952 per il set delle equazioni non conservative di Eulero per il caso monodimensionale. L’applicazione attraverso il software Matlab-Simulink di tali discretizzazioni numeriche garantisce la possibilità di calcolare la dinamica del flusso all’interno del condotto. L’innovazione proposta da questo lavoro consiste nel considerare l’intero stato da iterare come un vettore, permettendo di gestire parte delle operazioni da compiere con delle matrici. Questo approccio è stato adottato sia per una maggior velocità di calcolo, sia per rendere più agevole la modifica della geometria, ad esempio in fase di progettazione. La routine di lancio del nuovo modello, infatti, gestirà autonomamente la scrittura delle matrici, a partire dai pochi parametri necessari per la definizione della geometria all’interno del codice. Si andranno quindi a presentare i passaggi più importanti che hanno portato alla scrittura del codice, con particolare attenzione poi alla fase di validazione del modello. Essa sarà basata sia su un altro codice presente in letteratura, modellato anch’esso attraverso risoluzione CIR, sia mediante dati sperimentali utilizzati per la validazione di tale implementazione. Seguirà infine un’analisi dettagliata sui fattori che influenzano, positivamente e negativamente, l’esito delle simulazioni realizzate, come la discretizzazione spaziale e quella temporale, prestando sempre particolare attenzione alla stabilità del metodo.

CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

Studio dell'intrusione salina all'interno di un acquifero freatico costiero (Ravenna; Italia)

L’acquifero freatico costiero ravennate è intensamente salinizzato fino a diversi km nell’entroterra. Il corpo dell’acquifero è formato da sabbie che poggiano su un substrato argilloso ad una profondità media di 25 m, i depositi affioranti sono sabbie e argille. Il lavoro svolto consiste in una caratterizzazione dello stato di salinizzazione con metodologie indirette (geoelettrica) e metodologie dirette (letture dei parametri fisici delle acque in pozzo). I sondaggi elettrici verticali (V.E.S.) mostrano stagionalità dovuta alle differenti quantità di pioggia e quindi di ricarica, le aree con depositi superficiali ad alta conducibilità idraulica (sabbie) hanno una lente d’acqua dolce compresa tra 0,1 e 2,25 m di spessore, al di sotto della quale troviamo una zona di mescolamento con spessori che vanno da 1,00 a 12,00 m, mentre quando in superficie abbiamo depositi a bassa conducibilità idraulica (limi sabbiosi e argille sabbiose) la lente d’acqua dolce scompare e la zona di mescolamento è sottile. Le misure dirette in pozzo mostrano una profondità della tavola d’acqua quasi ovunque sotto il livello del mare in entrambi i mesi monitorati, Giugno e Dicembre 2010, presentando una profondità leggermente maggiore nel mese di Dicembre. Dalla ricostruzione litologica risulta un acquifero composto da 4×109 m3 di sabbia, per cui ipotizzando una porosità media del 30% sono presenti 1,2×109 m3 di acqua. Dalla modellazione numerica (Modflow-SEAWAT 2000) risulta che l’origine dell’acqua salata che si trova in falda trova più facilmente spiegazione ipotizzando la sua presenza fin dalla formazione dell’acquifero, residuo delle acque marine che regredivano. Un’altra problematica analizzata è valutare l’applicazione della metodologia a minifiltri in uno studio sulla salinizzazione delle acque di falda. É stata implementata la costruzione di un transetto sperimentale, che ha permesso la mappatura dell’interfaccia acqua dolce/salmastra/salata con una precisione finora non raggiungibile.

Evoluzione Morfotettonica delle Aree Alpine "Sempione" e "Brennero" attraverso Studi Termocronologici di bassa temperatura

E’ mostrata l’analisi e la modellazione di dati termocronologici di bassa temperatura da due regioni Alpine: il Sempione ed il Brennero. Le faglie distensive presenti bordano settori crostali profondi appartenenti al dominio penninico: il duomo metamorfico Lepontino al Sempione e la finestra dei Tauri al Brennero. I dati utilizzati sono FT e (U-Th)/He su apatite. Per il Sempione i dati provengono dalla bibliografia; per il Brennero si è provveduto ad un nuovo campionamento, sia in superficie che in sotterraneo. Gli attuali lavori per la galleria di base del Brennero (BBT), hanno consentito, per la prima volta, di raccogliere dati di FT e (U-Th)/He in apatite in sottosuolo per la finestra dei Tauri occidentale. Le analisi sono state effettuate tramite un codice a elementi finiti, Pecube, risolvente l’equazione di diffusione del calore per una topografia evolvente nel tempo. Il codice è stato modificato per tener conto dei dati sotterranei. L’inversione dei dati è stata effettuata usando il Neighbourhood Algorithm (NA), per ottenere il più plausibile scenario di evoluzione morfotettonico. I risultati ottenuti per il Sempione mostrano: ipotetica evoluzione dello stile tettonico della faglia del Sempione da rolling hinge a low angle detachment a 6.5 Ma e la cessazione dell’attività a 3 Ma; costruzione del rilievo fino a 5.5 Ma, smantellamento da 5.5 Ma ad oggi, in coincidenza dei cambiamenti climatici Messiniani e relativi all’inizio delle maggiori glaciazioni; incremento dell’esumazione da 0–0.6 mm/anno a 0.6–1.2 mm/anno a 2.4 Ma nell’emisfero settentrionale. I risultati al Brennero mostrano: maggiore attività tettonica della faglia del Brennero (1.3 mm/anno), maggiore attività esumativa (1–2 mm/anno) prima dei 10 Ma; crollo dell’attività della faglia del Brennero fra 10 Ma e oggi (0.1 mm/anno) e dell’attività esumativa nello stesso periodo (0.1–0.3 mm/anno); nessun aumento del tasso esumativo o variazioni topografiche negli ultimi 5 Ma.

Studio e sperimentazione del comportamento vibratorio di una lamella azionata da un attuatore piezoelettrico

Studio di una lamella di acciaio di dimensioni relativamente ridotte, a cui è incollato un attuatore piezoelettrico che ne causa la vibrazione, sfruttando il fenomeno della risonanza. Tale struttura possiede numerose peculiarità geometriche e di funzionamento che ne determinano il comportamento dinamico. Ci si prefigge l’obiettivo di studiare, con l’ausilio di prove sperimentali, il comportamento vibratorio del sistema, al fine di realizzare un modello agli elementi finiti che ne simuli il funzionamento e che sia impiegabile nella previsione degli effetti di modifiche strutturali e dinamiche.

Oxymoron,progetto di un grattacielo direzionale a Shenzhen, Cina applicazione della Soft Kill Option (SKO) nella progettazione architettonica

Lo scopo della tesi è quello di affrontare la progettazione con un approccio,quanto più attuale e per certi versi avanguardista, chiamato Parametric design (progettazione parametrica), accoppiato efficacemente col concetto di Arte generativa (in questo caso Architettura). Già nel 1957 Luigi Moretti affrontò il tema dell’architettura parametrico-generativa fondando l’IRMOU (Istituto per la Ricerca Matematica e Operativa applicata all'Urbanistica) e oggi è una mentalità molto diffusa nei più grandi studi del mondo. Il tema non è solo tecnologico o informatico strumentale, ma è proprio un modo di pensare e immaginare il possibile, costruito o naturale che sia. E’ un modo di vivere la propria creatività. L’aggettivo “generativa” è legato al fatto che l’arte in esame è generata seguendo regole preimpostate e ben definite dal progettista, coerentemente agli obiettivi e alle finalità del progetto. L’evoluzione delle stesse, seguendo relazioni molto semplici, può dar vita a risultati sorprendenti e inaspettati, dotati di una notevole complessità che però, se letta nell’insieme, è perfettamente in armonia con l’idea progettuale di partenza. Il fascino di questa materia è il legame entusiasmante che crea tra architettura, ingegneria, poesia, filosofia, matematica, biologia, fisica, pittura ecc ecc. Questo perché i concetti di evoluzione, di relazione e di generazione appartengono a tutto ciò che ci circonda, e quindi alla concezione umana di vita. E’ possibile in questo modo permeare il costrutto progettuale con principi e regole oggettivamente riconoscibili e apprezzabili dallo spettatore perché instrisi di una forte veridicità processuale. Il titolo "Oxymoron" è la traduzione inglese della figura retorica ossimoro,la quale è strettamente connessa all’ispirazione progettuale: proviene dall’indagine approfondita di processi evolutivi (distruttivi in questo caso) caratterizzanti realtà naturali che, esplorate con sempre più accuratezza, determinano morfologie e forme aventi profonde radici strutturali. La distruzione che crea lo spazio. La genesi stessa della forma segue predominanti algoritmi matematici governati e corretti da variabili di diversa natura che definiscono l'enviroment di influenze interagenti ed agenti sul campione di studio. In questo caso la ricerca è focalizzata su processi erosivi fisici e chimici, di agenti esterni (quali vento e sali rispettivamente) ,di cui materiali inorganici, quali minerali e aggregati degli stessi (rocce), sono soggetti. In particolare, l’interesse è approfondito su fenomeni apparentemente emergenti dei tafoni e dei cosiddetti Micro canyon. A tal scopo si sfrutterà un metodo di soft kill option (SKO) di ottimizzazione topologica (optimization topology) attraverso gli strumenti informatici più idonei quali software di modellazione parametrica e di calcolo computazionale. La sperimentazione sta proprio nell'utilizzare uno strumento concepito per uno scopo, con un'ottica strettamente ingegneristica, per un'altra meta, ossia ricavare e ottenere se possibile un metodo di lavoro o anche solo un processo generativo tale da riprodurre o simulare casi e situazioni riscontrabili in natura negli eventi soggetti a erosione. Il tutto coerente con le regole che stanno alla base della genesi degli stessi. Il parallelismo tra singolarità naturale e architettura risiede nella generazione degli spazi e nella combinazione di questi. L’ambizioso obiettivo è quello di innescare un ciclo generativo, che messo in comunicazione diretta con un contesto variegato ed eterogeneo, dia vita a una soluzione progettuale dall'alto contenuto morfologico e spaziale.

Conversione per il Semantic Web di dati Turistico-Culturali: il progetto QRPlaces - Semantic Events

Il progetto QRPlaces - Semantic Events, oggetto di questo lavoro, focalizza l’attenzione sull’analisi, la progettazione e l’implementazione di un sistema che sia in grado di modellare i dati, relativi a diversi eventi facenti parte del patrimonio turistico - culturale della Regione Emilia Romagna 1, rendendo evidenti i vantaggi associati ad una rappresentazione formale incentrata sulla Semantica. I dati turistico - culturali sono intesi in questo ambito sia come una rappresentazione di “qualcosa che accade in un certo punto ad un certo momento” (come ad esempio un concerto, una sagra, una raccolta fondi, una rappresentazione teatrale e quant’altro) sia come tradizioni e costumi che costituiscono il patrimonio turistico-culturale e a cui si fa spesso riferimento con il nome di “Cultural Heritage”. Essi hanno la caratteristica intrinseca di richiedere una conoscenza completa di diverse informa- zioni correlata, come informazioni di geo localizzazione relative al luogo fisico che ospita l’evento, dati biografici riferiti all’autore o al soggetto che è presente nell’evento piuttosto che riferirsi ad informazioni che descrivono nel dettaglio tutti gli oggetti, come teatri, cinema, compagnie teatrali che caratterizzano l’evento stesso. Una corretta rappresentazione della conoscenza ad essi legata richiede, pertanto, una modellazione in cui i dati possano essere interconnessi, rivelando un valore informativo che altrimenti resterebbe nascosto. Il lavoro svolto ha avuto lo scopo di realizzare un dataset rispondente alle caratteristiche tipiche del Semantic Web grazie al quale è stato possibile potenziare il circuito di comunicazione e informazione turistica QRPlaces 2. Nello specifico, attraverso la conversione ontologica di dati di vario genere relativi ad eventi dislocati nel territorio, e sfruttando i principi e le tecnologie del Linked Data, si è cercato di ottenere un modello informativo quanto più possibile correlato e arricchito da dati esterni. L’obiettivo finale è stato quello di ottenere una sorgente informativa di dati interconnessi non solo tra loro ma anche con quelli presenti in sorgenti esterne, dando vita ad un percorso di collegamenti in grado di evidenziare una ricchezza informativa utilizzabile per la creazione di valore aggiunto che altrimenti non sarebbe possibile ottenere. Questo aspetto è stato realizzato attraverso un’in- terfaccia di MashUp che utilizza come sorgente il dataset creato e tutti i collegamenti con la rete del Linked Data, in grado di reperire informazioni aggiuntive multi dominio.