Estudio de efectos de interacción estructura-suelo-estructura sobre la respuesta sísmica de edificios pilotados en función del tipo de onda sísmica y su ángulo de incidencia

Máster Universitario en Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI)

Reduced-order modelling, circuit-level design and SOI fabrication of microelectromechanical resonators

This thesis deals with two important research aspects concerning radio frequency (RF) microresonators and switches. First, a new approach for compact modeling and simulation of these devices is presented. Then, a combined process flow for their simultaneous fabrication on a SOI substrate is proposed. Compact models for microresonators and switches are extracted by applying mathematical model order reduction (MOR) to the devices finite element (FE) description in ANSYS c° . The behaviour of these devices includes forms of nonlinearities. However, an approximation in the creation of the FE model is introduced, which enables the use of linear model order reduction. Microresonators are modeled with the introduction of transducer elements, which allow for direct coupling of the electrical and mechanical domain. The coupled system element matrices are linearized around an operating point and reduced. The resulting macromodel is valid for small signal analysis around the bias point, such as harmonic pre-stressed analysis. This is extremely useful for characterizing the frequency response of resonators. Compact modelling of switches preserves the nonlinearity of the device behaviour. Nonlinear reduced order models are obtained by reducing the number of nonlinearities in the system and handling them as input to the system. In this way, the system can be reduced using linear MOR techniques and nonlinearities are introduced directly in the reduced order model. The reduction of the number of system nonlinearities implies the approximation of all distributed forces in the model with lumped forces. Both for microresonators and switches, a procedure for matrices extraction has been developed so that reduced order models include the effects of electrical and mechanical pre-stress. The extraction process is fast and can be done automatically from ANSYS binary files. The method has been applied for the simulation of several devices both at devices and circuit level. Simulation results have been compared with full model simulations, and, when available, experimental data. Reduced order models have proven to conserve the accuracy of finite element method and to give a good description of the overall device behaviour, despite the introduced approximations. In addition, simulation is very fast, both at device and circuit level. A combined process-flow for the integrated fabrication of microresonators and switches has been defined. For this purpose, two processes that are optimized for the independent fabrication of these devices are merged. The major advantage of this process is the possibility to create on-chip circuit blocks that include both microresonators and switches. An application is, for example, aswitched filter bank for wireless transceiver. The process for microresonators fabrication is characterized by the use of silicon on insulator (SOI) wafers and on a deep reactive ion etching (DRIE) step for the creation of the vibrating structures in single-crystal silicon and the use of a sacrificial oxide layer for the definition of resonator to electrode distance. The fabrication of switches is characterized by the use of two different conductive layers for the definition of the actuation electrodes and by the use of a photoresist as a sacrificial layer for the creation of the suspended structure. Both processes have a gold electroplating step, for the creation of the resonators electrodes, transmission lines and suspended structures. The combined process flow is designed such that it conserves the basic properties of the original processes. Neither the performance of the resonators nor the performance of the switches results affected by the simultaneous fabrication. Moreover, common fabrication steps are shared, which allows for cheaper and faster fabrication.

Comparison between 2.5-D and 3-D realistic models for wind field adjustment

[EN]In previous works, many authors have widely used mass consistent models for wind field simulation by the finite element method. On one hand, we have developed a 3-D mass consistent model by using tetrahedral meshes which are simultaneously adapted to complex orography and to terrain roughness length. In addition, we have included a local refinement strategy around several measurement or control points, significant contours, as for example shorelines, or numerical solution singularities. On the other hand, we have developed a 2.5-D model for simulating the wind velocity in a 3-D domain in terms of the terrain elevation, the surface temperature and the meteorological wind, which is consider as an averaged wind on vertical boundaries...

Two-dimensional numerical approach for the vibration isolation analysis on thin walled wave barriers in poroelastic soils.

[EN]This paper is concerned with the vibration isolation efficiency analysis of total or partially buried thin walled wave barriers in poroelastic soils. A two-dimensional time harmonic model that treats soils and structures in a direct way by combining appropriately the conventional Boundary Element Method (BEM), the Dual BEM (DBEM) and the Finite Element Method es developed to this aim.

Local scale air quality model with several pollutant sources

[EN]This paper shows a finite element method for pollutant transport with several pollutant sources. An Eulerian convection–diffusion–reaction model to simulate the pollutant dispersion is used. The discretization of the different sources allows to impose the emissions as boundary conditions. The Eulerian description can deal with the coupling of several plumes. An adaptive stabilized finite element formulation, specifically Least-Squares, with a Crank-Nicolson temporal integration is proposed to solve the problem. An splitting scheme has been used to treat separately the transport and the reaction. A mass-consistent model has been used to compute the wind field of the problem…

Some advances in open problems of isogeometric analysis

[EN]Isogeometric analysis (IGA) has arisen as an attempt to unify the fields of CAD and classical finite element methods. The main idea of IGA consists in using for analysis the same functions (splines) that are used in CAD representation of the geometry. The main advantage with respect to the traditional finite element method is a higher smoothness of the numerical solution and more accurate representation of the geometry. IGA seems to be a promising tool with wide range of applications in engineering. However, this relatively new technique have some open problems that require a solution. In this work we present our results and contributions to this issue…

Crustal heterogeneities and complexities of fault processes

In this work we study the relation between crustal heterogeneities and complexities in fault processes. The first kind of heterogeneity considered involves the concept of asperity. The presence of an asperity in the hypocentral region of the M = 6.5 earthquake of June 17-th, 2000 in the South Iceland Seismic Zone was invoked to explain the change of seismicity pattern before and after the mainshock: in particular, the spatial distribution of foreshock epicentres trends NW while the strike of the main fault is N 7◦ E and aftershocks trend accordingly; the foreshock depths were typically deeper than average aftershock depths. A model is devised which simulates the presence of an asperity in terms of a spherical inclusion, within a softer elastic medium in a transform domain with a deviatoric stress field imposed at remote distances (compressive NE − SW, tensile NW − SE). An isotropic compressive stress component is induced outside the asperity, in the direction of the compressive stress axis, and a tensile component in the direction of the tensile axis; as a consequence, fluid flow is inhibited in the compressive quadrants while it is favoured in tensile quadrants. Within the asperity the isotropic stress vanishes but the deviatoric stress increases substantially, without any significant change in the principal stress directions. Hydrofracture processes in the tensile quadrants and viscoelastic relaxation at depth may contribute to lower the effective rigidity of the medium surrounding the asperity. According to the present model, foreshocks may be interpreted as induced, close to the brittle-ductile transition, by high pressure fluids migrating upwards within the tensile quadrants; this process increases the deviatoric stress within the asperity which eventually fails, becoming the hypocenter of the mainshock, on the optimally oriented fault plane. In the second part of our work we study the complexities induced in fault processes by the layered structure of the crust. In the first model proposed we study the case in which fault bending takes place in a shallow layer. The problem can be addressed in terms of a deep vertical planar crack, interacting with a shallower inclined planar crack. An asymptotic study of the singular behaviour of the dislocation density at the interface reveals that the density distribution has an algebraic singularity at the interface of degree ω between -1 and 0, depending on the dip angle of the upper crack section and on the rigidity contrast between the two media. From the welded boundary condition at the interface between medium 1 and 2, a stress drop discontinuity condition is obtained which can be fulfilled if the stress drop in the upper medium is lower than required for a planar trough-going surface: as a corollary, a vertically dipping strike-slip fault at depth may cross the interface with a sedimentary layer, provided that the shallower section is suitably inclined (fault "refraction"); this results has important implications for our understanding of the complexity of the fault system in the SISZ; in particular, we may understand the observed offset of secondary surface fractures with respect to the strike direction of the seismic fault. The results of this model also suggest that further fractures can develop in the opposite quadrant and so a second model describing fault branching in the upper layer is proposed. As the previous model, this model can be applied only when the stress drop in the shallow layer is lower than the value prescribed for a vertical planar crack surface. Alternative solutions must be considered if the stress drop in the upper layer is higher than in the other layer, which may be the case when anelastic processes relax deviatoric stress in layer 2. In such a case one through-going crack cannot fulfil the welded boundary conditions and unwelding of the interface may take place. We have solved this problem within the theory of fracture mechanics, employing the boundary element method. The fault terminates against the interface in a T-shaped configuration, whose segments interact among each other: the lateral extent of the unwelded surface can be computed in terms of the main fault parameters and the stress field resulting in the shallower layer can be modelled. A wide stripe of high and nearly uniform shear stress develops above the unwelded surface, whose width is controlled by the lateral extension of unwelding. Secondary shear fractures may then open within this stripe, according to the Coulomb failure criterion, and the depth of open fractures opening in mixed mode may be computed and compared with the well studied fault complexities observed in the field. In absence of the T-shaped decollement structure, stress concentration above the seismic fault would be difficult to reconcile with observations, being much higher and narrower.

A model for diffusive and displacive phase transitions in steels

Heat treatment of steels is a process of fundamental importance in tailoring the properties of a material to the desired application; developing a model able to describe such process would allow to predict the microstructure obtained from the treatment and the consequent mechanical properties of the material. A steel, during a heat treatment, can undergo two different kinds of phase transitions [p.t.]: diffusive (second order p.t.) and displacive (first order p.t.); in this thesis, an attempt to describe both in a thermodynamically consistent framework is made; a phase field, diffuse interface model accounting for the coupling between thermal, chemical and mechanical effects is developed, and a way to overcome the difficulties arising from the treatment of the non-local effects (gradient terms) is proposed. The governing equations are the balance of linear momentum equation, the Cahn-Hilliard equation and the balance of internal energy equation. The model is completed with a suitable description of the free energy, from which constitutive relations are drawn. The equations are then cast in a variational form and different numerical techniques are used to deal with the principal features of the model: time-dependency, non-linearity and presence of high order spatial derivatives. Simulations are performed using DOLFIN, a C++ library for the automated solution of partial differential equations by means of the finite element method; results are shown for different test-cases. The analysis is reduced to a two dimensional setting, which is simpler than a three dimensional one, but still meaningful.

Plasmonenresonanzen von sichelförmigen metallischen Nanoobjekten

Metallische Objekte in der Größenordnung der optischen Wellenlänge zeigen Resonanzen im optischen Spektralbereich. Mit einer Kombination aus Kolloidlithographie, Metallfilmbedampfung und reaktivem Ionenstrahl¨atzen wurden Nanosicheln aus Gold bzw. Silber mit identischer Form und Orientierung in Sichelform mit einer Größe von 60nm bis 400nm hergestellt. Der Öffnungswinkel der Nanosicheln lässt sich kontinuierlich einstellen. Durch die einheitliche Orientierung lassen sich Messungen am Ensemble direkt auf das Verhalten des Einzelobjektes übertragen, wie ein Vergleich der Extinktionsspektren einer Ensemblemessung am UV/Vis/NIR-Spektrometer mit einer Einzelpartikelmessung in einem konfokalen Mikroskop zeigt. Die optische Antwort der Nanosicheln wurde als zwei-dimensionales Modell mit einer Finite Elemente Methode berechnet. Das Ergebnis sind mehrere polarisationsabhängige Resonanzen im optischen Spektrum. Diese lassen sich durch Variation des Öffnungswinkels und der Gr¨oße der Nanosichel verschieben. Durch Beleuchten lassen sich plasmonische Schwingungen anregen, die ein stark lokalisiertes Nahfeld an den Spitzen und in der Öffnung der Nanosicheln erzeugen. Das Nahfeld der Partikelresonanz wurde mit einer Fotolackmethode nachgewiesen. Die Untersuchungen am UV/Vis/NIR-Spektrometer zeigen mehrere polarisationsabhängige Resonanzen im Spektralbereich von 300 nm bis 3200 nm. Die Resonanzen der Nanosicheln lassen sich durch den Öffnungswinkel und den Durchmesser in der Größenordnung der Halbwertbreite im optischen Spektrum verschieben. In der Anwendung als Chemo- bzw. Biosensor zeigen Gold-Nanosicheln eine ähnliche Empfindlichkeit wie vergleichbare Sensoren auf der Basis von dünnen Metallstrukturen. Das Nahfeld zeichnet sich durch eine starke Lokalisierung aus und dringt, je nach Multipolordnung, zwischen 14 nm und 70 nm in die Umgebung ein. Quantenpunkte wurden an das Nahfeld der Nanosicheln gekoppelt. Die Emission der Quantenpunkte bei einer Wellenlänge von 860nm wird durch die Resonanz der Nanosicheln verstärkt. Die Nanosicheln wurden als optische Pinzette eingesetzt. Bei einer Anregung mit einem Laser bei einer Wellenlänge von 1064 nm wurden Polystyrolkolloide mit einem Durchmesser von 40 nm von den resonanten Nanosicheln eingefangen. Die Nanosicheln zeigen außergewöhnliche optische Eigenschaften, die mithilfe der Geometrieparameter über einen großen Bereich verändert werden können. Die ersten Anwendungen haben Anknüpfungspunkte zur Verwendung in der Sensorik, Fluoreszenzspektroskopie und als optische Pinzette aufgezeigt.

Magma injections and destabilization of basaltic volcanoes: a numerical study. Application to La Reunion (Indian ocean, France) and Stromboli (Tyrrhenian sea, Italy)

Most basaltic volcanoes are affected by recurrent lateral instabilities during their evolution. Numerous factors have been shown to be involved in the process of flank destabilization occurring over long periods of time or by instantaneous failures. However, the role of these factors on the mechanical behaviour and stability of volcanic edifices is poorly-constrained as lateral failure usually results from the combined effects of several parameters. Our study focuses on the morphological and structural comparison of two end-member basaltic systems, La Reunion (Indian ocean, France) and Stromboli (southern Tyrrhenian sea, Italy). We showed that despite major differences on their volumes and geodynamic settings, both systems present some similarities as they are characterized by an intense intrusive activity along well-developed rift zones and recurrent phenomena of flank collapse during their evolution. Among the factors of instability, the examples of la Reunion and Stromboli evidence the major contribution of intrusive complexes to volcano growth and destruction as attested by field observations and the monitoring of these active volcanoes. Classical models consider the relationship between vertical intrusions of magma and flank movements along a preexisting sliding surface. A set of published and new field data from Piton des Neiges volcano (La Reunion) allowed us to recognize the role of subhorizontal intrusions in the process of flank instability and to characterize the geometry of both subvertical and subhorizontal intrusions within basaltic edifices. This study compares the results of numerical modelling of the displacements associated with high-angle and low-angle intrusions within basaltic volcanoes. We use a Mixed Boundary Element Method to investigate the mechanical response of an edifice to the injection of magmatic intrusions in different stress fields. Our results indicate that the anisotropy of the stress field favours the slip along the intrusions due to cointrusive shear stress, generating flank-scale displacements of the edifice, especially in the case of subhorizontal intrusions, capable of triggering large-scale flank collapses on basaltic volcanoes. Applications of our theoretical results to real cases of flank displacements on basaltic volcanoes (such as the 2007 eruptive crisis at La Reunion and Stromboli) revealed that the previous model of subvertical intrusions-related collapse is a likely mechanism affecting small-scale steeply-sloping basaltic volcanoes like Stromboli. Furthermore, our field study combined to modelling results confirms the importance of shallow-dipping intrusions in the morpho-structural evolution of large gently-sloping basaltic volcanoes like Piton de la Fournaise, Etna and Kilauea, with particular regards to flank instability, which can cause catastrophic tsunamis.

Comparison between analytical and numerical methods for the assessment of masonry arch bridges: Case study of Clemente Bridge on Savio river

This research has focused on the study of the behavior and of the collapse of masonry arch bridges. The latest decades have seen an increasing interest in this structural type, that is still present and in use, despite the passage of time and the variation of the transport means. Several strategies have been developed during the time to simulate the response of this type of structures, although even today there is no generally accepted standard one for assessment of masonry arch bridges. The aim of this thesis is to compare the principal analytical and numerical methods existing in literature on case studies, trying to highlight values and weaknesses. The methods taken in exam are mainly three: i) the Thrust Line Analysis Method; ii) the Mechanism Method; iii) the Finite Element Methods. The Thrust Line Analysis Method and the Mechanism Method are analytical methods and derived from two of the fundamental theorems of the Plastic Analysis, while the Finite Element Method is a numerical method, that uses different strategies of discretization to analyze the structure. Every method is applied to the case study through computer-based representations, that allow a friendly-use application of the principles explained. A particular closed-form approach based on an elasto-plastic material model and developed by some Belgian researchers is also studied. To compare the three methods, two different case study have been analyzed: i) a generic masonry arch bridge with a single span; ii) a real masonry arch bridge, the Clemente Bridge, built on Savio River in Cesena. In the analyses performed, all the models are two-dimensional in order to have results comparable between the different methods taken in exam. The different methods have been compared with each other in terms of collapse load and of hinge positions.

Micro-FEM models based on micro-CT reconstructions for the in vitro characterization of the elastic properties of trabecular bone tissue.

This master’s thesis describes the research done at the Medical Technology Laboratory (LTM) of the Rizzoli Orthopedic Institute (IOR, Bologna, Italy), which focused on the characterization of the elastic properties of the trabecular bone tissue, starting from october 2012 to present. The approach uses computed microtomography to characterize the architecture of trabecular bone specimens. With the information obtained from the scanner, specimen-specific models of trabecular bone are generated for the solution with the Finite Element Method (FEM). Along with the FEM modelling, mechanical tests are performed over the same reconstructed bone portions. From the linear-elastic stage of mechanical tests presented by experimental results, it is possible to estimate the mechanical properties of the trabecular bone tissue. After a brief introduction on the biomechanics of the trabecular bone (chapter 1) and on the characterization of the mechanics of its tissue using FEM models (chapter 2), the reliability analysis of an experimental procedure is explained (chapter 3), based on the high-scalable numerical solver ParFE. In chapter 4, the sensitivity analyses on two different parameters for micro-FEM model’s reconstruction are presented. Once the reliability of the modeling strategy has been shown, a recent layout for experimental test, developed in LTM, is presented (chapter 5). Moreover, the results of the application of the new layout are discussed, with a stress on the difficulties connected to it and observed during the tests. Finally, a prototype experimental layout for the measure of deformations in trabecular bone specimens is presented (chapter 6). This procedure is based on the Digital Image Correlation method and is currently under development in LTM.

Locating the source of volcanic tremor at stromboli volcano, italy

We have developed a method for locating sources of volcanic tremor and applied it to a dataset recorded on Stromboli volcano before and after the onset of the February 27th 2007 effusive eruption. Volcanic tremor has attracted considerable attention by seismologists because of its potential value as a tool for forecasting eruptions and for better understanding the physical processes that occur inside active volcanoes. Commonly used methods to locate volcanic tremor sources are: 1) array techniques, 2) semblance based methods, 3) calculation of wave field amplitude. We have choosen the third approach, using a quantitative modeling of the seismic wavefield. For this purpose, we have calculated the Green Functions (GF) in the frequency domain with the Finite Element Method (FEM). We have used this method because it is well suited to solve elliptic problems, as the elastodynamics in the Fourier domain. The volcanic tremor source is located by determining the source function over a regular grid of points. The best fit point is choosen as the tremor source location. The source inversion is performed in the frequency domain, using only the wavefield amplitudes. We illustrate the method and its validation over a synthetic dataset. We show some preliminary results on the Stromboli dataset, evidencing temporal variations of the volcanic tremor sources.

Detecting interfaces in a parabolic-elliptic problem from surface measurements

Assuming that the heat capacity of a body is negligible outside certain inclusions the heat equation degenerates to a parabolic-elliptic interface problem. In this work we aim to detect these interfaces from thermal measurements on the surface of the body. We deduce an equivalent variational formulation for the parabolic-elliptic problem and give a new proof of the unique solvability based on Lions’s projection lemma. For the case that the heat conductivity is higher inside the inclusions, we develop an adaptation of the factorization method to this time-dependent problem. In particular this shows that the locations of the interfaces are uniquely determined by boundary measurements. The method also yields to a numerical algorithm to recover the inclusions and thus the interfaces. We demonstrate how measurement data can be simulated numerically by a coupling of a finite element method with a boundary element method, and finally we present some numerical results for the inverse problem.