Simulación numérica de elementos estructurales de hormigón armado de alta resistencia sometidos a explosiones

En los últimos años ha crecido el interés por el uso del hormigón de alta resistencia en la construcción de estructuras resistentes a explosiones y de elementos de protección frente a impactos y explosiones. En el presente trabajo se analiza la capacidad de los modelos numéricos existentes a la hora de reproducir resultados experimentales obtenidos sobre hormigones de alta resistencia sometidos a cargas explosivas y cargas estáticas. Para ello se emplea una campaña experimental desarrollada por la Agencia de Defensa de Suecia sobre vigas de este material. Los resultados obtenidos demuestran la capacidad de las herramientas numéricas para reproducir adecuadamente el comportamiento del hormigón si bien, paradójicamente, proporcionan mejores resultados en régimen dinámico que en régimen estático. Las simulaciones numéricas, junto con el análisis de los resultados experimentales, permiten obtener conclusiones acerca del diferente comportamiento de los elementos estructurales sometidos a altas velocidades de deformación frente a su comportamiento en régimen estático.

Experimental and numerical analysis of reinforced concrete elements subjected to blast loading

El hormigón es uno de los materiales de construcción más empleados en la actualidad debido a sus buenas prestaciones mecánicas, moldeabilidad y economía de obtención, entre otras ventajas. Es bien sabido que tiene una buena resistencia a compresión y una baja resistencia a tracción, por lo que se arma con barras de acero para formar el hormigón armado, material que se ha convertido por méritos propios en la solución constructiva más importante de nuestra época. A pesar de ser un material profusamente utilizado, hay aspectos del comportamiento del hormigón que todavía no son completamente conocidos, como es el caso de su respuesta ante los efectos de una explosión. Este es un campo de especial relevancia, debido a que los eventos, tanto intencionados como accidentales, en los que una estructura se ve sometida a una explosión son, por desgracia, relativamente frecuentes. La solicitación de una estructura ante una explosión se produce por el impacto sobre la misma de la onda de presión generada en la detonación. La aplicación de esta carga sobre la estructura es muy rápida y de muy corta duración. Este tipo de acciones se denominan cargas impulsivas, y pueden ser hasta cuatro órdenes de magnitud más rápidas que las cargas dinámicas impuestas por un terremoto. En consecuencia, no es de extrañar que sus efectos sobre las estructuras y sus materiales sean muy distintos que las que producen las cargas habitualmente consideradas en ingeniería. En la presente tesis doctoral se profundiza en el conocimiento del comportamiento material del hormigón sometido a explosiones. Para ello, es crucial contar con resultados experimentales de estructuras de hormigón sometidas a explosiones. Este tipo de resultados es difícil de encontrar en la literatura científica, ya que estos ensayos han sido tradicionalmente llevados a cabo en el ámbito militar y los resultados obtenidos no son de dominio público. Por otra parte, en las campañas experimentales con explosiones llevadas a cabo por instituciones civiles el elevado coste de acceso a explosivos y a campos de prueba adecuados no permite la realización de ensayos con un elevado número de muestras. Por este motivo, la dispersión experimental no es habitualmente controlada. Sin embargo, en elementos de hormigón armado sometidos a explosiones, la dispersión experimental es muy acusada, en primer lugar, por la propia heterogeneidad del hormigón, y en segundo, por la dificultad inherente a la realización de ensayos con explosiones, por motivos tales como dificultades en las condiciones de contorno, variabilidad del explosivo, o incluso cambios en las condiciones atmosféricas. Para paliar estos inconvenientes, en esta tesis doctoral se ha diseñado un novedoso dispositivo que permite ensayar hasta cuatro losas de hormigón bajo la misma detonación, lo que además de proporcionar un número de muestras estadísticamente representativo, supone un importante ahorro de costes. Con este dispositivo se han ensayado 28 losas de hormigón, tanto armadas como en masa, de dos dosificaciones distintas. Pero además de contar con datos experimentales, también es importante disponer de herramientas de cálculo para el análisis y diseño de estructuras sometidas a explosiones. Aunque existen diversos métodos analíticos, hoy por hoy las técnicas de simulación numérica suponen la alternativa más avanzada y versátil para el cálculo de elementos estructurales sometidos a cargas impulsivas. Sin embargo, para obtener resultados fiables es crucial contar con modelos constitutivos de material que tengan en cuenta los parámetros que gobiernan el comportamiento para el caso de carga en estudio. En este sentido, cabe destacar que la mayoría de los modelos constitutivos desarrollados para el hormigón a altas velocidades de deformación proceden del ámbito balístico, donde dominan las grandes tensiones de compresión en el entorno local de la zona afectada por el impacto. En el caso de los elementos de hormigón sometidos a explosiones, las tensiones de compresión son mucho más moderadas, siendo las tensiones de tracción generalmente las causantes de la rotura del material. En esta tesis doctoral se analiza la validez de algunos de los modelos disponibles, confirmando que los parámetros que gobiernan el fallo de las losas de hormigón armado ante explosiones son la resistencia a tracción y su ablandamiento tras rotura. En base a los resultados anteriores se ha desarrollado un modelo constitutivo para el hormigón ante altas velocidades de deformación, que sólo tiene en cuenta la rotura por tracción. Este modelo parte del de fisura cohesiva embebida con discontinuidad fuerte, desarrollado por Planas y Sancho, que ha demostrado su capacidad en la predicción de la rotura a tracción de elementos de hormigón en masa. El modelo ha sido modificado para su implementación en el programa comercial de integración explícita LS-DYNA, utilizando elementos finitos hexaédricos e incorporando la dependencia de la velocidad de deformación para permitir su utilización en el ámbito dinámico. El modelo es estrictamente local y no requiere de remallado ni conocer previamente la trayectoria de la fisura. Este modelo constitutivo ha sido utilizado para simular dos campañas experimentales, probando la hipótesis de que el fallo de elementos de hormigón ante explosiones está gobernado por el comportamiento a tracción, siendo de especial relevancia el ablandamiento del hormigón. Concrete is nowadays one of the most widely used building materials because of its good mechanical properties, moldability and production economy, among other advantages. As it is known, it has high compressive and low tensile strengths and for this reason it is reinforced with steel bars to form reinforced concrete, a material that has become the most important constructive solution of our time. Despite being such a widely used material, there are some aspects of concrete performance that are not yet fully understood, as it is the case of its response to the effects of an explosion. This is a topic of particular relevance because the events, both intentional and accidental, in which a structure is subjected to an explosion are, unfortunately, relatively common. The loading of a structure due to an explosive event occurs due to the impact of the pressure shock wave generated in the detonation. The application of this load on the structure is very fast and of very short duration. Such actions are called impulsive loads, and can be up to four orders of magnitude faster than the dynamic loads imposed by an earthquake. Consequently, it is not surprising that their effects on structures and materials are very different than those that cause the loads usually considered in engineering. This thesis broadens the knowledge about the material behavior of concrete subjected to explosions. To that end, it is crucial to have experimental results of concrete structures subjected to explosions. These types of results are difficult to find in the scientific literature, as these tests have traditionally been carried out by armies of different countries and the results obtained are classified. Moreover, in experimental campaigns with explosives conducted by civil institutions the high cost of accessing explosives and the lack of proper test fields does not allow for the testing of a large number of samples. For this reason, the experimental scatter is usually not controlled. However, in reinforced concrete elements subjected to explosions the experimental dispersion is very pronounced. First, due to the heterogeneity of concrete, and secondly, because of the difficulty inherent to testing with explosions, for reasons such as difficulties in the boundary conditions, variability of the explosive, or even atmospheric changes. To overcome these drawbacks, in this thesis we have designed a novel device that allows for testing up to four concrete slabs under the same detonation, which apart from providing a statistically representative number of samples, represents a significant saving in costs. A number of 28 slabs were tested using this device. The slabs were both reinforced and plain concrete, and two different concrete mixes were used. Besides having experimental data, it is also important to have computational tools for the analysis and design of structures subjected to explosions. Despite the existence of several analytical methods, numerical simulation techniques nowadays represent the most advanced and versatile alternative for the assessment of structural elements subjected to impulsive loading. However, to obtain reliable results it is crucial to have material constitutive models that take into account the parameters that govern the behavior for the load case under study. In this regard it is noteworthy that most of the developed constitutive models for concrete at high strain rates arise from the ballistic field, dominated by large compressive stresses in the local environment of the area affected by the impact. In the case of concrete elements subjected to an explosion, the compressive stresses are much more moderate, while tensile stresses usually cause material failure. This thesis discusses the validity of some of the available models, confirming that the parameters governing the failure of reinforced concrete slabs subjected to blast are the tensile strength and softening behaviour after failure. Based on these results we have developed a constitutive model for concrete at high strain rates, which only takes into account the ultimate tensile strength. This model is based on the embedded Cohesive Crack Model with Strong Discontinuity Approach developed by Planas and Sancho, which has proved its ability in predicting the tensile fracture of plain concrete elements. The model has been modified for its implementation in the commercial explicit integration program LS-DYNA, using hexahedral finite elements and incorporating the dependence of the strain rate, to allow for its use in dynamic domain. The model is strictly local and does not require remeshing nor prior knowledge of the crack path. This constitutive model has been used to simulate two experimental campaigns, confirming the hypothesis that the failure of concrete elements subjected to explosions is governed by their tensile response, being of particular relevance the softening behavior of concrete.

Validación numérica del comportamiento cíclico de conexiones interiores en estructuras con forjados reticulares

En este artículo se analizan distintas alternativas para modelar conexiones interiores forjado reticularpilar de hormigón armado y se propone una solución que con un coste computacional razonable proporcione resultados muy próximos a los experimentales. El modelo propuesto se valida numéricamente con los resultados obtenidos previamente de un ensayo de laboratorio realizado sobre un espécimen a escala 3/5 de la misma tipología estructural, sometido a ensayo cuasi-estático hasta el colapso. Para tal fin, se aplica la teoría avanzada de mecánica no lineal de hormigón armado sobre un modelo tridimensional de elementos finitos en análisis no-lineal: elementos tipo fibra para componentes unidimensionales y tipo capa para los bidimensionales. Los resultados principales de la simulación fueron: a) la obtención de la curva de capacidad mediante un análisis de empuje incremental (push-over) con control de desplazamientos; b) las curvas histeréticas en la losa, y c) los patrones de fisuras, logrando una concordancia revelada con los resultados experimentales.

Análisis de metodologías para la medida de la resistencia a tracción por spalling

En el presente trabajo se presenta un estudio teórico y experimental de la determinación de la resistencia a tracción de materiales frágiles a partir de la técnica experimental de spalling. Se utilizan diferentes metodologías propuestas por varios autores para determinar la resistencia a tracción por spalling y se lleva a cabo un estudio mediante simulaciones numéricas de las diversas variables que influyen en este tipo de ensayos. Además, se realiza una campaña experimental de ensayos a una alúmina del 99,5% de pureza cuyos resultados son utilizados para la determinación de la resistencia a tracción de este material a partir de tres métodos diferentes propuestos por varios autores. Se expone el estudio y comparación de los resultados experimentales obtenidos de resistencia a tracción de la alúmina empleando técnicas de fotografía a alta velocidad y un sistema de correlación digital de imágenes. Los resultados muestran que la resistencia a tracción obtenida difiere en función de las diferentes metodologías propuestas.A theoretical and experimental study of the tensile strength of brittle materials using the experimental procedure of spalling of long bars is presented in this article. Different methodologies proposed by several authors are used to obtain the tensile strength of Al2O3 monolithic ceramic. The hypotheses needed for the experimental set-up are also checked, and the requirements of the set-up and the variables are also studied by means of numerical simulations. An experimental campaign has been carried out including high speed video and a digital image correlation system to obtain the tensile strength of alumina. Finally, a comparison of the test results provided by three different methods proposed by different authors are presented. The tensile strength obtained from three different methods on the same specimens provides different results.

La barrera del análisis estructural y la representación gráfica en el desarrollo de los proyectos arquitectónicos: el caso de la Ópera de Sidney

En numerosas ocasiones a lo largo de la historia la imaginación de los creadores ha ido por delante de las posibilidades técnicas de cada momento. Así, muchas de estas nuevas ideas han requerido largos periodos de tiempo para materializarse como realidad construida, hasta que el desarrollo tecnológico e industrial hubo alcanzado un grado de madurez suficiente. En el campo de la arquitectura, estas limitaciones técnicas se han ido acotando paulatinamente hasta desembocar en la situación actual en la que cualquier planteamiento formal puede ser representado gráficamente y analizado desde un punto de vista estructural, superádose de este modo la barrera existente históricamente en el tratamiento de las formas. A lo largo del presente tesis doctoral se analiza cómo la formulación del Método de los Elementos Finitos en la década de los cincuenta y las curvas de Bézier en la década de los sesenta del siglo pasado y la posterior generalización de los ordenadores personales y de los programas informáticos asociados (C.A.D. y F.E.M. principalmente) en los estudios de arquitectura e ingeniería a partir de la década de los noventa, posibilitó el desarrollo de cualquier propuesta arquitectónica, por compleja que ésta fuese, provocando una verdadera revolución a nivel formal en el mundo de la arquitectura, especialmente en el campo de la edificación singular o icónica. Se estudia este proceso a través de ocho edificios; cuatro anteriores y otros tantos posteriores a la desaparición de la barrera anteriormente referida, establecida de forma simbólica en la década de los años ochenta del siglo XX: Frontón de Recoletos en Madrid, Edificio Seagram en Nueva York, Habitat ’67 en Montreal, Ópera de Sídney, museo Guggenheim de Bilbao, ampliación del Victoria & Albert Museum en Londres, tanatorio “Meiso no Mori” en Gifu y nueva sede de la CCTV en Pekín. De entre ellos, la Ópera de Sídney, obra del arquitecto danés Jørn Utzon, condensa gran parte de los aspectos relevantes investigados en relación a la influencia que los métodos de representación y análisis estructural ejercen en la concepción y construcción de las obras de arquitectura. Por este motivo y por considerarse un hito de la arquitectura a nivel global se toma como caso de estudio. La idea general del edificio, que data de 1956, se enmarca en una época inmediatamente anterior a la del desarrollo científico y tecnológico anteriormente referido. Esta ausencia de herramientas de diseño disponibles acordes a la complejidad formal de la propuesta planteada condicionó enormente la marcha del proyecto, dilatándose dramáticamente en el tiempo y disparándose su coste hasta el punto de que el propio arquitecto danés fue separado de las obras antes de su conclusión. Además, la solución estructural finalmente construida de las cubiertas dista mucho de la prevista por Utzon inicialmente. Donde él había imaginado unas finas láminas de hormigón flotando sobre el paisaje se materializó una estructura más pesada, formada por costillas pretensadas de hormigón con unas secciones notablemente mayores. La forma también debió ser modificada de modo ostensible respecto a la propuesta inicial. Si este edificio se pretendiese construir en la actualidad, con toda seguridad el curso de los acontecimientos se desarrollaría por senderos muy diferentes. Ante este supuesto, se plantean las siguientes cuestiones: ¿sería posible realizar un análisis estructural de la cubierta laminar planteada por Utzon inicialmente en el concurso con las herramientas disponibles en la actualidad?; ¿sería dicha propuesta viable estructuralmente?. A lo largo de las siguientes páginas se pretende dar respuesta a estas cuestiones, poniendo de relieve el impacto que los ordenadores personales y los programas informáticos asociados han tenido en la manera de concebir y construir edificios. También se han analizado variantes a la solución laminar planteada en la fase de concurso, a través de las cuales, tratando en la medida de lo posible de ajustarse a las sugerencias que Ove Arup y su equipo realizaron a Jørn Utzon a lo largo del dilatado proceso de proyecto, mejorar el comportamiento general de la estructura del edificio. Por último, se ha pretendido partir de cero y plantear, desde una perspectiva contemporánea, posibles enfoques metodológicos aplicables a la búsqueda de soluciones estructurales compatibles con la forma propuesta originalmente por Utzon para las cubiertas de la Ópera de Sídney y que nunca llegó a ser construida (ni analizada), considerando para ello los medios tecnológicos, científicos e industriales disponibles en la actualidad. Abstract On numerous occasions throughout history the imagination of creators has gone well beyond of the technical possibilities of their time. Many new ideas have required a long period to materialize, until the technological and industrial development had time to catch up. In the architecture field, these technical limitations have gradually tightened leading to the current situation in which any formal approach can be represented and analyzed from a structural point of view, thus concluding that the structural analysis and the graphical representation’s barrier in the development of architectural projects has dissappeared. Throughout the following pages it is examined how the development of the Finite Element Method in the fifties and the Bezier curves in the sixties of the last century and the subsequent spread of personal computers and specialized software in the architectural and engineering offices from the nineties, enabled the development of any architectural proposal independently of its complexity. This has caused a revolution at a formal level in architecture, especially in the field of iconic building. This process is analyzed through eight buildings, four of them before and another four after the disappearance of the above mentioned barrier, roughly established in the eighties of the last century: Fronton Recoletos in Madrid, Seagram Building in New York Habitat '67 in Montreal, Sydney Opera House, Guggenheim Museum Bilbao, Victoria & Albert Museum extension in London, Crematorium “Meiso no Mori” in Gifu and the new CCTV headquarters in Beijing. Among them, the Sydney Opera House, designed by Danish architect Jørn Utzon, condenses many of the main aspects previously investigated regarding the impact of representation methods and structural analysis on the design and construction of architectural projects. For this reason and also because it is considered a global architecture milestone, it is selected as a case study. The building’s general idea, which dates from 1956, is framed at a time immediately preceding the above mentioned scientific and technological development. This lack of available design tools in accordance with the proposal’s formal complexity conditioned enormously the project’s progress, leading to a dramatic delay and multiplying the final budget disproportionately to the point that the Danish architect himself was separated from the works before completion. Furthermore, the built structure differs dramatically from the architect’s initial vision. Where Utzon saw a thin concrete shell floating over the landscape a heavier structure was built, consisting of prestressed concrete ribs with a significantly greater size. The geometry also had to be modified. If this building were to built today, the course of events surely would walk very different paths. Given this assumption, a number of questions could then be formulated: Would it be possible to perform a structural analysis of Utzon’s initially proposed competition-free-ways roof’s geometry with the tools available nowadays?; Would this proposal be structurally feasable?. Throughout the following pages it is intended to clarify this issues, highlighting personal computers and associated software’s impact in building design and construction procedures, especially in the field of iconic building. Variants have also been analyzed for the laminar solution proposed in the competition phase, through which, trying as far as possible to comply with the suggestions that Ove Arup and his team did to Jørn Utzon along the lengthy process project, improving the overall performance of the building structure. Finally, we have started from scratch and analyzed, from a contemporary perspective, possible structural solutions compatible with Utzon’s Opera House’s original geometry and vision –proposal that was never built (nor even analyzed)-, taking into consideration the technological, scientific and industrial means currently available.

Depreciación diferencial y estructura urbana : análisis tipológico de la edificción residencial en la almendra central de Madrid

En la práctica de la valoración inmobiliaria se han desarrollado diferentes métodos para la obtención del valor, ajustados a las características del bien a tasar. El procedimiento para la determinación del valor del inmueble difiere si se trata de un suelo vacante, o una vivienda de segunda mano. La aplicación de los distintos métodos técnicos de valoración necesita de información concreta sobre las condiciones del mercado. Sin embargo, esta información no siempre está disponible. Por ello, se hace precisa la búsqueda de procedimientos de valoración alternativos, que permitan aplicar métodos de valoración inusuales en los diferentes submercados. Por otra parte, la aplicación del método del coste y, más concretamente, la cuantificación de la depreciación, desvincula completamente la pérdida de valor de un inmueble de la localización del mismo. Esta tesis doctoral se plantea la búsqueda de un modelo que relacione la pérdida de valor con variables de localización, partiendo de la hipótesis de que tal relación existe, y además es determinante en la pérdida de valor. Para ello, se ha seleccionado un ámbito de análisis muy concreto, tipológica y geográficamente: la edificación residencial colectiva en la Almendra Central de Madrid. El espectro es suficientemente acotado como para ser abarcable, pero suficientemente complejo como para permitir la existencia de una variabilidad que permita extraer conclusiones de interés para el campo de la valoración inmobiliaria. Con carácter previo a la aplicación de un análisis econométrico, se ha estudiado en profundidad el proceso de formación del espacio urbano y la evolución de su tejido residencial. Se ha prestado especial atención a aquellos aspectos que pudieran ser más relevantes para la caracterización de la construcción, ya sea por factores de localización o cronológicos. En este sentido, se demostrará que es especialmente relevante la identificación de un edificio con una determinada época, y no tanto con su fecha de construcción. Esta relación no es lineal a lo largo del tiempo, identificándose periodos más penalizados, desde el punto de vista de la valoración. La utilización de modelos de precios hedónicos permite confirmar la estrecha relación entre el mercado secundario y el de obra nueva: tiene más peso en la cuantificación de la depreciación el valor de la vivienda nueva, que el estado de conservación del inmueble a valorar.

Estudio Hidrogeológico de la Llanura Manchega mediante simulación numérica

El objetivo del presente estudio es el análisis de un sistema complejo, el Acuífero de la Mancha Occidental, mediante un modelo numérico de simulación que represente de la manera más rigurosa posible la evolución del Sistema Acuífero 23. Este modelo se realiza en régimen permanente con una rigurosa configuración del sistema desde el punto de vista geológico y geométrico. De esta forma se deja iniciado y planeado un modelo y su estructura que será una base real de futuras formulaciones transitorias, base de los sucesivos análisis de explotación y predicción. . Las distintas situaciones que en las últimas décadas ha experimentado el Sistema 23, soluciones que se han dado para las mismas, y los posibles planes de actuación que se podrían llevar a cabo en un futuro frente a condiciones cambiantes de clima y explotación se podrán estudiar a partir del presente modelo de simulación numérica. ABSTRACT The main purpose of this study is to analyses a complex system such as Western La Mancha Aquifer by the use of a numeric model that simulates as accurately as possible the evolution of Aquifer 23. This model is made in steady state with a thorough configuration of the system from both the geological and geometric point of view. Therefore, it is left initiated with a planned model and its structure which will be used as a real base for following transient flow simulations that also will be the foundation of subsequent prediction and exploitation analysis. The different situations that Aquifer 23 has experienced during the last decades; the solutions that have been given for them; and the possible plans that would be implemented in the future in order to deal with the changing environmental and exploitation conditions will be able to be studied using this model.

Análisis del comportamiento en servicio de los puentes de fábrica ante las nuevas condiciones de explotación

Los puentes arco de fábrica representan una parte importante del inventario de puentes en España y en Europa. De aquí, la importancia cuantitativa de estas estructuras y el interés que representa tanto para las Administraciones, organismos públicos y privados como para la sociedad en general, mantener estos puentes en servicio. Para poder alargar su vida útil, aún más si cabe, se plantea imprescindible llevar a cabo una labor de conservación y mantenimiento adecuada. La importancia de la evaluación estructural de los puentes de fábrica se basa en la conveniencia de conocer el comportamiento estructural tanto en condiciones de servicio como en agotamiento. El comportamiento en Estado Límite Último ha sido estudiado en profundidad y como resultado, la seguridad frente a agotamiento ha quedado definida con cierta confianza. Sin embargo, el comportamiento en Estado Límite de Servicio no es tan conocido. Se considera necesaria una revisión del concepto de ELS en puentes arco de fábrica puesto que muchos de ellos, cumpliendo las comprobaciones establecidas para ELU, presentan daños asociados a las actuales condiciones de explotación existentes actuales. El presente trabajo se enmarca dentro de un proyecto que está llevando a cabo la UIC (Union Internationale de Chemins de Fer) desde 2004. El objetivo general del mismo es conocer el comportamiento de los puentes arco de fábrica y mejorar los métodos de evaluación, mantenimiento y reparación existentes. Con este estudio, se pretende contribuir a mejorar la caracterización del comportamiento en servicio de estas estructuras. Para ello se ha realizado un análisis del estado tensional de la bóveda y el relleno estudiando la influencia de la configuración geométrica de los puentes y las propiedades de los materiales que los componen. Entender el funcionamiento de estas estructuras y conocer sus peculiares características ha sido el paso previo al desarrollo del estudio. Para ello, se ha acudido al minucioso trabajo desarrollado por diversos autores, desde los primeros análisis mediante la línea de presiones hasta los actuales métodos basados en elementos finitos. En primer lugar, se ha realizado un estudio paramétrico de diferentes configuraciones geométricas. El objetivo es comprobar la influencia de las dimensiones absolutas y las relaciones existentes entre las mismas en el comportamiento en servicio de los puentes. A continuación, se ha estudiado la influencia de las propiedades mecánicas de los rellenos rígido y granular (modificando su módulo de deformabilidad) en el estado tensional de la bóveda y el propio relleno. Para las dos etapas anteriores se emplea un modelo simplificado de puente arco de fábrica. A continuación, se crea un modelo más completo del mismo, para estudiar la influencia de las propiedades mecánicas del substrato de cimentación en el comportamiento estructural en servicio de la estructura. El objetivo final de este trabajo fin de máster, es conocer el comportamiento en servicio de los puentes arco de fábrica, para poder establecer las variables que condicionan el mismo. El siguiente paso en este campo de estudio, consiste en fijar unos valores límite para las mismas. Como síntesis de las conclusiones obtenidas, se menciona que el comportamiento estructural de los puente arco de fábrica, en mayor o menor medida, está influenciado por su geometría y las propiedades de los rellenos y substrato de cimentación que forman parte del mismo. Masonry arch bridges represent an important part of the total bridges, both in Spain and Europe. For this reason, these structures are so important in terms of quantity. They play an essential role for public administration and people in general. In order to remain those bridges serviceable, suitable inspections and repairs are required. It is necessary to have knowledge of the bridge condition so that an assessment can be made with confidence. It is therefore necessary to check that bridges behave properly under Ultimate Limit State criterion and Serviceability Limit State criterion. ULS has been studied thoroughly and as result, requirements under collapse are well described. However, structural behaviour under conditions of serviceability is not well evaluated. Serviceability Limit State criterion should be redefined for masonry arch bridges assessment. This is because many bridges that fulfil SLS criterion have damages. The purpose of this work is to take part in the improvement of assessment under serviceability. A vault and fill structural behaviour analysis has been performed from two points of view: geometry of the bridge and materials’ properties. The starting point was learning about masonry arch bridges behaviour and their special features. For this task, a revision of several authors’ thorough study has been made, from the line of thrust analysis to the current finite element analysis. Firstly, a parametric study of typologies of bridge has been made. The aim is to know how vault and fill behaviour changes modifying both absolute dimensions and relation between them. In the next step, a study of both surface fill and backfill properties has been done. The aim is to know how vault and fill behaviour changes modifying Young’s modulus. The principal conclusion achieved along this work is that the structural behaviour of masonry bridges depends of the bridge geometry and its material properties.

Nuevos desarrollos en la caracterización geométrica y simulación mecánica de paredes arteriales

El presente trabajo resume los desarrollos realizados por el Grupo de Mecánica Computacional de la ETSICCP de la UPM en el campo de la biomecánica. Se presentan estrategias de caracterización geométrica de paredes arteriales a partir de tomografía computarizada y de reconstrucciones tridimensionales del cayado aórtico, se describe una metodología de implementación de las tensiones iniciales dentro del contexto del análisis por elementos finitos y se expone un modelo constitutivo de daño continuo que busca reproducir la degradación material del tejido aórtico.

Desarrollo de métodos de procesado de imágenes laríngeas en color para la detección del espacio glotal y de las cuerdas vocales

La medicina ha evolucionado de forma que las imágenes digitales tienen un papel de gran relevancia para llevar a cabo el diagnóstico de enfermedades. Son muchos y de diversa naturaleza los problemas que pueden presentar el aparato fonador. Un paso previo para la caracterización de imágenes digitales de la laringe es la segmentación de las cuerdas vocales. Hasta el momento se han desarrollado algoritmos que permiten la segmentación de la glotis. El presente proyecto pretende avanzar un paso más en el estudio, procurando asimismo la segmentación de las cuerdas vocales. Para ello, es necesario aprovechar la información de color que ofrecen las imágenes, pues es lo que va a determinar la diferencia entre una región y otra de la imagen. En este proyecto se ha desarrollado un novedoso método de segmentación de imágenes en color estroboscópicas de la laringe basado en el crecimiento de regiones a partir de píxeles-semilla. Debido a los problemas que presentan las imágenes obtenidas por la técnica de la estroboscopia, para conseguir óptimos resultados de la segmentación es necesario someter a las imágenes a un preprocesado, que consiste en la eliminación de altos brillos y aplicación de un filtro de difusión anisotrópica. Tras el preprocesado, comienza el crecimiento de la región a partir de unas semillas que se obtienen previamente. La condición de inclusión de un píxel en la región se basa en un parámetro de tolerancia que se determina de forma adaptativa. Este parámetro comienza teniendo un valor muy bajo y va aumentando de forma recursiva hasta alcanzar una condición de parada. Esta condición se basa en el análisis de la distribución estadística de los píxeles dentro de la región que va creciendo. La última fase del proyecto consiste en la realización de las pruebas necesarias para verificar el funcionamiento del sistema diseñado, obteniéndose buenos resultados en la segmentación de la glotis y resultados esperanzadores para seguir mejorando el sistema para la segmentación de las cuerdas vocales. ABSTRACT Medicine has evolved so that digital images have a very important role to perform disease diagnosis. There are wide variety of problems that can present the vocal apparatus. A preliminary step for characterization of digital images of the larynx is the segmentation of the vocal folds. To date, some algorithms that allow the segmentation of the glottis have been developed. This project aims to go one step further in the study, also seeking the segmentation of the vocal folds. To do this, we must use the color information offered by images, since this is what will determine the difference between different regions in a picture. In this project a novel method of larynx color images segmentation based on region growing from a pixel seed is developed. Due to the problems of the images obtained by the technique of stroboscopy, to achieve optimal results of the segmentation is necessary a preprocessing of the images, which involves the removal of high brightness and applying an anisotropic diffusion filter. After this preprocessing, the growth of the region from previously obtained seeds starts. The condition for inclusion of a pixel in the region is based on a tolerance parameter, which is adaptively determined. It initially has a low value and this is recursively increased until a stop condition is reached. This condition is based on the analysis of the statistical distribution of the pixels within the grown region. The last phase of the project involves the necessary tests to verify the proper working of the designed system, obtaining very good results in the segmentation of the glottis and encouraging results to keep improving the system for the segmentation of the vocal folds.

Estudio y aplicación de redes neuronales artificiales al análisis y diseño de antenas reflectarrays

En los últimos años, estamos siendo testigos de la alta implantación en la sociedad de dispositivos de comunicación. Lo que hace años estaba reservado a un público reducido, con claras necesidades en comunicación, se ha trasladado al público general, dado la amplia variedad de servicios que sobre los nuevos medios de comunicación se han desarrollado. De hecho, el mayor tráfico de datos en la actualidad no se produce al hilo de necesidades de máxima importancia, sino como producto de nuevos hábitos cotidianos. En este contexto de renovación tecnológica constante en busca de la eficiencia, las antenas reflectoras reflectarray (o, simplemente, los reflectarrays, RAs, [1]) se presentan como una opción competitiva contra los reflectores parabólicos metálicos. En su versión más simple, una antena reflectarray se trata de una estructura compuesta de un elemento alimentador radiante, como puede ser una bocina, y de una superficie plana, consistente en multitud de elementos individuales dispuestos en una rejilla periódica. Sobre esta superficie plana, los frentes de onda provenientes del alimentador son reflejados formando frentes de ondas planas, de una manera análoga a como lo hace un reflector parabólico. A partir de la configuración inicial, y centrándose en el principio de funcionamiento, se ha ido modificando el tipo de elemento RA empleado, consiguiendo RA cada vez más operativos. Es, sobre todo, con el nacimiento de la tecnología impresa cuando las antenas RAs vuelven a cobrar interés. Aunque el uso de tecnología impresa supuso un gran impulso en los RAs, también abrió otros desafíos en lo que al diseño de ellos se refiere. Desde el punto de vista del análisis, es común suponer que el elemento RA se encuentra en un ambiente infinitamente periódico, de forma que se puedan aplicar las condiciones de contorno de Floquet (suposición de periodicidad local). Desde un punto de vista funcional, en general, los elementos RA de tecnología impresa presentan un ancho de banda reducido, que condiciona el ancho de banda del RA completo. Entre las soluciones aportadas, es comúnmente aceptado que las estructuras multicapa, con resonadores a distintas frecuencias cercanas, pueden mitigar en parte el problema del ancho de banda. Por ello, en la actualidad, los elementos RA más comunes están compuestos por varios elementos resonadores, cuyas dimensiones constituyen los parámetros de diseño libres. Es decir, en función de dichas dimensiones, el elemento RA tendrá un valor del coeficiente de reflexión u otro. Esto supone un aumento en la complejidad a la hora de analizar dicho elemento por los métodos numéricos conocidos, como el Método de los Momentos (MoM) o el Método de Elementos Finitos (FEM, por las siglas de su traducción inglesa Finite Element Method), que redundará en un mayor tiempo de cómputo en el análisis. Por otra parte, como se muestra en la Figura R.1, el diseño de un RA conlleva analizar multitud de veces el elemento RA considerado. En efecto, se trata de un método de diseño indirecto, en donde las dimensiones de los parámetros geométricos libres de cada elemento RA se obtienen de manera iterativa usando un optimizador. Se ve claro, entonces, que el aumento en tiempo de análisis del elemento RA repercute en gran medida en el tiempo de diseño total, por lo que una reducción en el tiempo de análisis del elemento RA podría ser muy beneficioso. Uno de los métodos para conseguir reducir el tiempo de diseño de un RA, que se puede encontrar en la literatura, es emplear un modelo de la respuesta del elemento RA en función de los parámetros libres. La cuestión que aflora es cuál es la herramienta idónea para modelar la respuesta del elemento RA. En los últimos años se han propuestos varias formas. La primera de ellas consistía en encontrar un equivalente circuital. Esta aproximación está bien extendida para otras estructuras EM, donde los equivalentes circuitales con componentes LC ofrecen respuestas muy precisas con respecto a las que ofrecen las estructuras EM en sí. A raíz del carácter no lineal de la respuesta, hay autores que han propuesto para el diseño de RAs la creación de tablas de datos (look up tables) que, para cada parámetro de diseño de interés (suele ser el desfase introducido por el elemento) guardan las dimensiones de los parámetros geométricos libres asociados. De esta forma, consiguen un diseño rápido, pero poco versátil, ya que la tabla ofrece un único valor para cada entrada, por lo que es difícil jugar con más de una restricción de diseño. Más recientemente, se está comenzando a utilizar, para la caracterización de estructuras EM, unos sistemas llamados Redes Neuronales Artificiales (ANN, por sus siglas en inglés Artificial Neural Network). El uso fundamental de los mismos en EM es el de servir como interpoladores no lineales. Se trata de sistemas que admiten múltiples parámetros de entradas y múltiples parámetros de salida. Antes de poder ser usados como interpoladores, deben ser entrenados. Para ello, necesitan de un conjunto de pares de los parámetros de entrada a la red, con los valores de las salidas asociados. Algunos usos en electromagnetismo de las ANNs que se pueden encontrar en la literatura son: el modelado de filtros; la caracterización de dispositivos activos; la obtención de modelos que aceleran los algoritmos que calculan la dirección de llegada en antenas de radar; o el diseño de arrays de antenas. Volviendo al modelado de elementos RA, en este trabajo haremos uso de las ANNs para caracterizar distintos tipos de elementos RA. A lo largo de estos últimos años, se ha considerado esta posibilidad como una de las más prometedoras. De hecho, podemos encontrar algunas pocas referencias al respecto, varias de las cuales han sido publicadas por distintos autores durante la elaboración del trabajo recogido en esta Tesis. Como veremos, los resultados que vamos a presentar aportan novedades con respecto a la citada literatura. Particularmente, en este trabajo se ha realizado la caracterización de un elemento RA de tres capas, considerando hasta 9 parámetros de entrada (seis parámetros geométricos, las dos coordenadas del ángulo de incidencia, y la frecuencia) y 4 parámetros de salida complejos (los coeficientes de reflexión para dos polarizaciones ortogonales lineales). Haciendo uso de esta caracterización en el flujo de diseño de RAs, se ha realizado el análisis y el diseño de varias antenas RA con restricciones de diseño de comunicaciones espaciales. Los resultados fueron exitosos comparados con los resultados obtenidos por los métodos tradicionales. De manera puntualizada, podríamos resumir las aportaciones que se verán en esta Tesis como: Caracterización de distintos elementos RA mediante ANNs basadas en el Perceptrón Multicapa (MLP). En concreto, se ha realizado con éxito la caracterización de un elemento RA de parche acoplado a línea de retardo a través de apertura; la caracterización de un elemento RA basado en dipolos sobre substratos de distintas características eléctricas en el rango de centenas de GHz; y la caracterización de un elemento RA basado en 3 parches apilados, con 9 parámetros libres en su caracterización. Uso del FEM, de la técnica de segmentación en subdominios y de la generación y manipulación de accesos MAM para el análisis y la caracterización de elementos RA mediante ANNs. Desarrollo de una nueva técnica de obtención de muestras, para el caso de estructura multicapa cuyo estudio EM se pueda dividir en dos pasos: estudio de cada capa y conexión de capas. De esta forma, se ha podido reducir en varios órdenes de magnitud el tiempo necesario para obtener el set de entrenamiento de las ANNs. Valoración del uso de distintos métodos de entrenamiento de segundo orden para el entrenamiento de redes ANN MLP, en la caracterización de elementos RA. Desarrollo de una nueva técnica para realizar el entrenamiento de redes ANNs basadas en el MLP, denominada como Entrenamiento en Cascada. Dado el alto número de parámetros a caracterizar, era difícil conseguir una red que, partiendo del número de entradas deseado, proporcionara convergencia con precisión suficiente. Con el algoritmo propuesto y probado en esta Tesis, se consiguió entrenar redes de 8 parámetros de entradas (el noveno parámetro, la frecuencia, correspondía a redes diferentes para cada valor) con gran precisión. Desarrollo de un método adaptativo para mejorar la precisión de las ANNs en el análisis de antenas RA. Este método, basado en re-entrenar las ANNs para sub rangos de los parámetros de entrada en donde el error es mayor, aporta una precisión mayor, al mejorar el entrenamiento global de las ANNs, en un tiempo aceptable, ya que solo se incluyen nuevas muestras en torno a los valores donde el error es mayor. Análisis de antena RA completa, con cobertura según especificaciones de la misión AMAZONAS (haz conformado, banda Ku), usando las caracterización el elemento RA obtenida mediante ANNs. La mejora en tiempo de análisis conseguida con respecto al uso del MoM está en un factor 102, con precisiones comparables. Diseño de antenas RA completas, con especificaciones de haz pincel y EuTELSAT (banda Ku). De nuevo, la mejora en tiempo de diseño conseguida están en torno a 102. De todos los puntos anteriores, son de destacar los dos últimos, que forman el objetivo principal de esta Tesis. Esto es, el uso de modelos rápidos de elementos RA mediante ANNs para el análisis y el diseño de antenas para comunicaciones por satélite.

Estudio Acústico y Electroacústico de la sala de conciertos Ritmo y Compás (Madrid)

En el presente proyecto se ha realizado un estudio de las condiciones acústicas y del refuerzo sonoro de la sala de conciertos Ritmo y Compás (Madrid), repasando las principales magnitudes de calidad acústica. El proyecto combina medidas y simulaciones. Con las medidas se obtienen los parámetros acústicos de la sala, que permiten cuantificar la sensación que producen en la audiencia las cualidades acústicas de un espacio. Mediante software de simulación electroacústica (EASE) se ha generado el modelo geométrico y acústico de la sala, ajustado en base a los resultados experimentales. El objetivo fundamental del trabajo es generar un modelo virtual que refleje fielmente las condiciones reales de la sala de conciertos, de forma que se puedan estudiar en profundidad las características acústicas y electroacústicas del recinto. En el proceso de análisis se ha querido resaltar la importancia de la acústica en salas de conciertos de música amplificada, muchas veces relegada a un segundo plano por el diseño electroacústico. La calidad de una sala de este tipo depende estrechamente de la sonorización, pero su inevitable relación con las cualidades acústicas del recinto obliga a establecer una serie de criterios acústicos mínimos que aseguren las mejores condiciones para el sistema sonoro. Dado el peso de las simulaciones, en el proyecto se presentan los resultados obtenidos con distintos métodos y recursos de análisis software para apoyar el estudio, con una información completa que muestre, con la máxima claridad, el potencial de la sala Ritmo y Compás como sala de conciertos. SUMMARY. This project is about the acoustical and electroacustical studio of “Ritmo y Compás”, an important concert hall in Madrid, analyzing the main magnitudes of acoustic quality. The project combines measurements and simulations. With the measurements it is obtained the acoustic parameters of the hall, allowing the quantification of the sensation produced in the audience by the acoustic attributes of the hall. With the acoustic simulation software (Ease), based on the experimental results, the geometric and acoustic model has been created. The main purpose of this study is to generate a virtual model that accurately reflects the real conditions of the concert hall, allowing the deep study of the acoustic and electroacoustic features of the hall. In the process of the analysis, the importance of the acoustic characteristic in the amplified music concert halls was emphasized, which is often underestimate because of the electroacoustic design. The quality of this kind of hall strictly depends on the sound electrical system, but its inevitable relation with the acoustic characteristics of the hall, forces to establish a series of minimum acoustic rules that assures the best conditions for the sound system. Due to the importance of the simulations in the project, the results are presented with different methods and analysis software resources to back up the study with complete information that shows the maximum quality and clarity of the potential of the enclosure Ritmo y Compás as a concert hall.

Análisis del acoplamiento vibro-acústico en resonadores de membrana

El presente Trabajo Fin de Máster pretende llevar a cabo el análisis del comportamiento vibratorio de resonadores de membrana, consistentes en un panel delgado y ligero montado a cierta distancia de un elemento constructivo rígido y pesado. Este tipo de sistemas resonantes son empleados habitualmente como absorbentes de media-baja frecuencia en aplicaciones de acondicionamiento acústico de salas. El análisis hará especial hincapié en la influencia del acoplamiento mecánico-acústico entre la placa vibrante (estructura) y el colchón de aire (fluido) encerrado entre la misma y la pared rígida. En primer lugar, realizaremos el análisis modal experimental del resonador objeto de ensayo a partir de las mediciones de su respuesta vibratoria, con el fin de caracterizar su comportamiento en base a sus primeros modos propios acoplados de flexión. El análisis de las señales vibratorias en el dominio de la frecuencia para la identificación de dicho modos se realizará en el entorno de programación MATLAB, haciendo uso de una herramienta propia que implementa los métodos de cálculo y los algoritmos necesarios para tal fin. Asimismo, simularemos el comportamiento del resonador mediante el método de elementos finitos (FEM), utilizando las aplicaciones ANSYS y SYSNOISE, considerando diferentes condiciones frontera en el modelo generado. Los resultados aquí obtenidos serán de utilidad para complementar aquellos obtenidos de forma experimental a la hora de extraer conclusiones prácticas del análisis realizado. SUMMARY. This Master's Thesis intends to carry out the analysis of the vibratory behaviour of resonance absorbers, consisting of a thin and lightweight panel mounted at a distance from a rigid wall. Such systems are commonly used as sound absorption systems for mid-low frequency in room acoustics applications. The analysis will emphasize the influence of mechanical-acoustic coupling between the vibrating plate (structure) and the air cushion (acoustic element) enclosed behind it. First of all, we are performing the experimental modal analysis of the resonance absorber under test from the vibrational response measurements, in order to characterize its behaviour based on its first bending coupled-modes. The analysis of vibration signals in the frequency domain for the identification of such modes will be made in MATLAB programming environment, using a proprietary tool that implements the calculation methods and algorithms needed for this purpose. Furthermore, we are simulating the behaviour of the resonance absorber applying the Finite Element Method (FEM) – using ANSYS and SYSNOISE applications - considering different boundary conditions in the model created. The results from the simulation will be useful to complement those obtained experimentally when drawing practical conclusions from this analysis.

Análisis acústico del estudio de grabación Ritmo & Compás

En la actualidad, el éxito de una grabación musical depende, en gran medida, de la acústica que caracteriza la sala de control y la sala de grabación de un estudio donde se haya realizado ésta, así como del equipo utilizado. La acústica de un estudio de grabación viene determinada por la absorción, reverberación y aislamiento acústico de ese recinto. Estos son los factores que le dan la personalidad las salas. En este proyecto, se analiza el comportamiento acústico de las dos salas que componen estudio de grabación de „Ritmo & Compas‟, una sala de control y una sala de grabación, mediante la realización de medidas in situ. Se comparan los resultados con los valores aconsejados y se observa si ambas salas cumplen las condiciones óptimas que corresponden a cada una de éstas. En el caso de haber deficiencias acústicas, se proporcionan soluciones de acondicionamiento acústico en las salas, que se validan mediante simulación. Además, se realiza el estudio sobre el aislamiento a ruido aéreo entre los dos recintos colindantes que componen el estudio de grabación. ABSTRACT. Nowadays, the success of a musical recording depends, to a large extent, on the good acoustics that characterizes the recording studio where it has been made, as well as the equipment used. The acoustics of a recording studio is determined by absorption, reverberation and sound insulation of the enclosure. These are the factors that give personality to the site. This project aims to study the acoustic behavior of the two rooms that make the recording studio 'Ritmo & Compas', a control room and a recording room, by performing in situ measurements. Results will be compared with the recommended values and examined if both rooms meet the optimal conditions that correspond to each of these. In the case of having acoustic deficiencies will be provided acoustic solutions in the rooms, which are validated by simulation. In addition, there will be the study of airborne sound insulation between the two neighboring enclosures that make up the recording studio.

Modelización y simulación numérica de sistemas geotérmicos planos

En el presente proyecto se estudiará la conveniencia de utilizar simulación numérica para evaluar, diseñar y mejorar el diseño de instalaciones geotérmicas planas. Para ello, se modelizará el conjunto de terreno y sistema geotérmico mediante el acoplamiento de un modelo tridimensional (terreno) y un modelo unidimensional, formado por la red de tubería geotérmica, que estará dispuesta en un plano cuya orientación dependerá de las aplicaciones. Para la simulación numérica se utilizará un software de análisis y resolución por elementos finitos tipo comercial. Finalmente, se definirán unos parámetros geométricos de diseño y, teniendo en cuenta las propiedades térmicas de cada material, se realizarán estudios paramétricos para determinar la influencia de cada grupo de parámetros. Ello requerirá la realización de cientos de simulaciones numéricas que permitirán optimizar el rendimiento de sistemas geotérmicos planos tanto horizontales (soleras), como verticales (muros pantalla). ABSTRACT In this project it is going to be studied the convenience of numerical simulation as to evaluate, design and improve the design of geothermal planar installations. With this purpose, the terrain and the geothermal system will be modeled by coupling a three-dimensional model (terrain) and a one-dimensional model formed by the geothermal piping network disposed in a plane with an orientation which will depend on the applications. For the numerical simulation commercial software of analyzing and resolving by finite elements will be used. Finally, geometric design parameters will be defined, and taking into account the thermal properties of each material, parametric studies will be performed to determine the influence of each group of parameters. This will require the completion of hundreds of numerical simulations which will enhance the performance of both horizontal (slabs) and vertical (diaphragm walls) geothermal systems.