Evaluación de parámetros que influyen en el transporte de cloruros en hormigón parcialmente saturado

El deterioro del hormigón debido a la presencia del ion cloruro es causa frecuente de problemas en estructuras localizadas en ambiente marino y alta montaña. Su principal efecto consiste en la despasivación del acero de refuerzo embebido en el hormigón y su consecuente inicio de la corrosión del mismo. El ingreso del ion cloruro al interior del hormigón, está condicionado por una serie de parámetros de origen medioambiental e intrínsecos del hormigón. En función de estos parámetros el ingreso de cloruros en el hormigón puede deberse principalmente a los siguientes mecanismos: difusión y succión capilar. El estudio y evaluación de la resistencia del hormigón frente a cloruros, se ha desarrollado principalmente en condiciones saturadas del hormigón. Lo que ha significado que parámetros de importancia no sean considerados. Debido a esto, distintos procesos que suceden en estructuras reales no han sido identificados y estudiados correctamente. En este trabajo, se diseñó un programa de investigación para evaluar los parámetros que influyen en el transporte cloruros en hormigones no saturados. Para esto se diseñaron tres dosificaciones diferentes de hormigón. En la primera se empleó únicamente cemento portland, para el resto se utilizaron adiciones minerales (humo de sílice y escoria de alto horno). El empleo de adiciones se debió a que tienen un papel importante en la durabilidad de hormigones frente a cloruros. Los hormigones fueron dosificados con una relación agua/material cementício de 0,40 para el hormigón elaborado únicamente con cemento portland y 0,45 para las mezclas con adiciones. Para evaluar las propiedades de los hormigones en estado fresco y endurecido se realizaron ensayos vigentes en las normativas. Con los resultados obtenidos se determinaron parámetros de resistencia mecánica, microestructurales, resistencia al transporte de cloruros e higroscópicos. Una vez caracterizados los hormigones, se diseñó una propuesta experimental para estudiar los principales parámetros presentes en estructuras reales con presencia de cloruros. Tanto la concentración de cloruro como las condiciones ambientales se han variado teniendo como referencia las situaciones reales que podrían producirse en ambientes de alta montaña en la zona centro de España. La propuesta experimental consistió en tratar de evaluar la capacidad de los hormigones al transporte de iones en ambientes de alta montaña con presencia de sales fundentes. Para esto se establecieron 5 fases experimentales donde los principales parámetros ambientales y la presencia de iones agresivos sufrieron variaciones. Al término de cada fase se obtuvieron perfiles de penetración de cloruros en los hormigones y se evaluó la influencia de los parámetros presentes en cada fase. Los resultados experimentales se implementaron en un modelo numérico basado en la teoría de elementos finitos, desarrollado por el grupo de investigación del Departamento de Materiales de Construcción. Para esto fue necesario realizar la calibración y validación del modelo numérico para cada hormigón. El calibrado del modelo precisa de datos químicos y microestructurales de cada hormigón, tales como: capacidad de combinación de cloruros y propiedades difusivas e higroscópicas. Para la validación del modelo numérico se realizaron simulaciones de la propuesta experimental. Los resultados obtenidos se compararon con los valores experimentales. Con el objeto de poder estudiar en mayor profundidad la influencia del grado de saturación del hormigón durante la difusión de cloruros, se llevó a cabo una campaña experimental que consideró distintos grados de saturación en los hormigones. Para esto se establecieron en los hormigones cuatro grados de saturación distintos (50%, 60%, 80% y 100%, aproximadamente), posteriormente se expusieron a cloruro de sodio finamente molido. Una vez transcurrido el tiempo necesario se obtuvieron experimentalmente los perfiles de penetración de cloruros para cada grado de saturación y se calcularon los coeficientes de difusión. Los datos obtenidos durante la campaña experimental han demostrado la influencia positiva que ejercen las adiciones en las mezclas de hormigón. Sus principales ventajas son el refinamiento de la red porosa y el aumento en la capacidad de combinación de cloruros, además de mejorar sus propiedades mecánicas. La porosidad total en las mezclas no presentó grandes cambios, sin embargo, el cambio en la distribución del tamaño de poros es importante en las muestras con adiciones. En especial las fabricadas con humo de sílice. Los coeficientes de difusión y migración de cloruros para las mezclas con adiciones disminuyeron significativamente, igual que los valores de resistividad eléctrica. En los ensayos de penetración del agua bajo presión, fueron las muestras con adiciones las que mostraron las menores penetraciones. Los resultados obtenidos al final de la propuesta experimental permitieron estudiar los distintos parámetros involucrados. Se observó claramente que el proceso de difusión provoca el mayor transporte de cloruros hacia el interior del hormigón. Así mismo se comprobó que el lavado superficial y el secado de las probetas, trasladan cloruros hacia las zonas externas del hormigón. El primero debido a una baja concentración de cloruros externa, mientras que el secado provoca el movimiento de la solución de poro hacia las zonas de secado depositando cloruros en ellas. Las medidas higroscópicas permitieron determinar la existencia de dos zonas distintas en el interior del hormigón. La primera se localizó en el rango de 0-10mm, aproximadamente, en ésta se puso de manifiesto una mayor sensibilidad a los cambios experimentados en el exterior de las probetas. La segunda zona se localizó a una profundidad mayor de 10mm, aproximadamente. Se observó claramente una baja influencia de los cambios externos, siendo la difusión de cloruros el principal mecanismo de transporte presente en ella. En cuanto al estudio de la influencia del grado de saturación en la difusión de cloruros, se observó claramente una marcada diferencia entre los coeficientes de difusión de cloruro obtenidos. Para grados de saturación mayores del 80% el mecanismo de penetración de cloruros por difusión existe de forma significativa. Mientras que para valores inferiores los resultados revelaron que las vías de acceso disminuyen (poros conectado con agua) considerablemente limitando en un alto grado la penetración del agresivo. Para grados de saturación inferiores del 50% los valores del coeficiente de difusión son despreciables. The deterioration of concrete due to chloride ions is a frequent problem identified in structures located in marine and high-mountain environments. After entering the outer layer of the concrete, the chlorides tend to penetrate until they reach and then depassivate the steel bars. Subsequently, this induces the deterioration process of the reinforced concrete. This chloride penetration depends on the environmental conditions and intrinsic parameters of the concrete. Several transport mechanisms, such as diffusion, capillary suction and permeability can be present into the concrete. While recent research into the study and evaluation of concretes with chloride presence has been carried out in saturated concrete, it has not considered certain parameters that can modify this condition. Consequently, at the time of writing several processes that take place in real structures have not been identified and studied. In this work a research programme is designed to evaluate the parameters that influence chloride transport into non-saturated concrete. For this, three concrete mixes were designed by using high-early-strength Portland cement and mineral admixtures (silica fume and blast-slag furnace). The water-cement ratio was 0.40 for the concrete made solely with Portland cement and 0.45 for the concretes that used mineral admixtures as a cement replacement. A set of experimental tests were performed to evaluate the concrete properties both in fresh and hardened state. In addition, an experimental simulation was carried out under laboratory conditions in which the main objective was to assess resistance of concrete to chloride penetration under high-mountain conditions with the presence of de-icing salts. The environmental conditions and surface chloride concentration of the concrete used during the experimental simulation were chosen by considering conditions found in the high-mountain environment in central Spain. For the experimental simulation five phases were designed by varying the environmental parameters and concrete surface concentration. At the end of each phase a chloride profile was obtained with the aim of assessing the influence of the parameters on chloride transport. The experimental results were then used to calibrate and validate a numerical model based on finite element theory developed by the research team from the Construction Materials Department in a previous work. In order to carry out model calibration chemical and microstructural data for the concretes was required, such as binding capacity and the diffusive and hygroscopic properties. The experimental results were compared with the numerical simulations and provided a good fit. With the objective of studying the influence of the degree of concrete saturation on chloride diffusion, an experimental programme was designed. This entailed four saturation degrees (50%, 60%, 80% and 100%) being established in several concrete samples. The samples were then exposed to ground sodium chloride. Once the time required was achieved, the chloride profiles and diffusion coefficients were obtained for each saturation degree. The results obtained from the experimental program revealed a positive influence of the mineral admixtures on the concretes. Their effects were reflected in the pore-network refinement and the increase of chloride binding capacity, together with the improvement of the mechanical properties of the concretes. Total porosity did not reveal any notable change, though the pore-size distribution showed a significant degree of change in the concretes with mineral admixtures, specifically the samples prepared through use of silica fume. The chloride diffusion and migration coefficient, as well as the electrical resistivity values, decreased significantly in the concretes with admixtures. In the water penetration under pressure test, the concretes with admixtures presented the lowest penetration depth. The results obtained in the experimental simulation allowed study of the main parameters involved during the chloride penetration processes in non-saturated concretes in the presence of chlorides. According to the results, the diffusion process was the transport mechanism that transferred the greatest amount of chlorides into the concrete samples. In addition, a chloride movement toward external zones of the concrete, caused by the washing of concrete surface and the drying processes, was identified. The washing occurred when the concrete surface came into contact with a low-chloride concentration solution which initiated the outward movement of chloride diffusion. The drying processes corresponded to a movement of pore solution launched by water evaporation from the outer layer. Furthermore, hygroscopic measurements made in the concrete allowed two areas with distinct behavioural patterns to be identified. The first one, located in the range of approximately 0-10mm, showed a greater degree of influence regarding the changes of the external conditions. The second, situated at depths greater than 10mm, displayed a low influence of external conditions. The main process in this area was diffusion. Study of the influence of the degree of concrete saturation on chloride diffusion showed a clear difference among the chloride diffusion coefficients obtained. For degrees of concrete saturation higher than 80%, chloride penetration by diffusion tends to be significant. However, in the case of extent of saturation of lower than 80%, the results revealed that the access zone through which chlorides can penetrate decreased considerably. For degrees of concrete saturation lower than 50%, the chloride diffusion coefficients were negligible.

New concepts and techniques for the development of high-efficiency concentrating photovoltaic modules

El interés por los sistemas fotovoltaicos de concentración (CPV) ha resurgido en los últimos años amparado por el desarrollo de células multiunión de muy alta eficiencia basadas en semiconductores de los grupos III-V. Estas células han permitido obtener módulos de concentración con eficiencias que prácticamente duplican las del panel plano y que llegan al 35% en los módulos récord. Esta tesis está dedicada al diseño y la implementación experimental de nuevos conceptos que permitan obtener módulos CPV que no sólo alcancen una eficiencia alta en condiciones estándar sino que, además, sean lo suficientemente tolerantes a errores de montaje, seguimiento, temperatura y variaciones espectrales para que la energía que producen a lo largo del año sea máxima. Una de las primeras cuestiones que se abordan es el diseño de elementos ópticos secundarios para sistemas cuyo primario es una lente de Fresnel y que permiten, para una concentración fija, aumentar el ángulo de aceptancia y la tolerancia del sistema. Varios secundarios reflexivos y refractivos han sido diseñados y analizados mediante trazado de rayos. En particular, utilizando óptica anidólica y basándose en el diseño de una sola etapa conocido como ‘concentrador dieléctrico que funciona por reflexión total interna‘, se ha diseñado, fabricado y caracterizado un secundario con salida cuadrada que, usado junto con una lente de Fresnel, permite alcanzar simultáneamente una elevada eficiencia, concentración y aceptancia. Además, se ha propuesto y prototipado un método alternativo de fabricación para otro de los secundarios, denominado domo, consistente en el sobremoldeo de silicona sobre células solares. Una de las características que impregna todo el trabajo realizado en esta tesis es la aproximación holística en el diseño de módulos CPV, es decir, se ha prestado especial atención al diseño conjunto de la célula y la óptica para garantizar que el sistema total alcance la mayor eficiencia posible. En este sentido muchos sistemas ópticos desarrollados en esta tesis han sido diseñados, caracterizados y optimizados teniendo en cuenta que el ajuste de corriente entre las distintas subcélulas que comprenden la célula multiunión bajo el concentrador sea muy próximo a uno. La capa antirreflectante sobre la célula funciona, en cierto modo, como interfaz entre la óptica y la célula, por lo que se ha diseñado un método de optimización de capas antirreflectantes que considera no sólo el amplio rango de longitudes de onda para el que las células multiunión son sensibles sino también la distribución angular de intensidad sobre la célula creada por la óptica de concentración. Además, la cuestión de la falta de uniformidad también se ha abordado mediante la comparación de las distribuciones espectrales y espaciales de irradiancia que crean diferentes ópticas (simuladas mediante trazado de rayos y fotografiadas) y las pérdidas de eficiencia que experimentan las células iluminadas por dichas ópticas de concentración medidas experimentalmente. El efecto de la temperatura en la óptica de concentración también ha sido objeto de estudio de esta tesis. En particular, mediante simulaciones de elementos finitos se han dado los primeros pasos para el análisis de las deformaciones que sufren los dientes de las lentes de Fresnel híbridas (vidrio-silicona), así como el cambio de índice de refracción con la temperatura y la influencia de ambos efectos sobre el funcionamiento de los sistemas. Se ha implementado un modelo que tiene por objeto considerar las variaciones ambientales, principalmente temperatura y contenido espectral de la radiación directa, así como las sensibilidades térmica y espectral de los sistemas CPV, con el fin de maximizar la energía producida por un módulo de concentración a lo largo de un año en un emplazamiento determinado. Los capítulos 5 y 6 de este libro están dedicados al diseño, fabricación y caracterización de un nuevo concepto de módulo fotovoltaico denominado FluidReflex y basado en una única etapa reflexiva con dieléctrico fluido. En este nuevo concepto la presencia del fluido aporta algunas ventajas significativas como son: un aumento del producto concentración por aceptancia (CAP, en sus siglas en inglés) alcanzable al rodear la célula con un medio cuyo índice de refracción es mayor que uno, una mejora de la eficiencia óptica al disminuir las pérdidas por reflexión de Fresnel en varias interfaces, una mejora de la disipación térmica ya que el calor que se concentra junto a la célula se trasmite por convección natural y conducción en el fluido y un aislamiento eléctrico mejorado. Mediante la construcción y medida de varios prototipos de unidad elemental se ha demostrado que no existe ninguna razón fundamental que impida la implementación práctica del concepto teórico alcanzando una elevada eficiencia. Se ha realizado un análisis de fluidos candidatos probando la existencia de al menos dos de ellos que cumplen todos los requisitos (en particular el de estabilidad bajo condiciones de luz concentrada) para formar parte del sistema de concentración FluidReflex. Por ´ultimo, se han diseñado, fabricado y caracterizado varios prototipos preindustriales de módulos FluidReflex para lo cual ha sido necesario optimizar el proceso de fabricación de la óptica multicavidad a fin de mantener el buen comportamiento óptico obtenido en la fabricación de la unidad elemental. Los distintos prototipos han sido medidos, tanto en el laboratorio como bajo el sol real, analizando el ajuste de corriente de la célula iluminada por el concentrador FluidReflex bajo diferentes distribuciones espectrales de la radiación incidente así como el excelente comportamiento térmico del módulo. ABSTRACT A renewed interest in concentrating photovoltaic (CPV) systems has emerged in recent years encouraged by the development of high-efficiency multijunction solar cells based in IIIV semiconductors that have led to CPV module efficiencies which practically double that of flat panel PV and which reach 35% for record modules. This thesis is devoted to the design and experimental implementation of new concepts for obtaining CPV modules that not only achieve high efficiency under standard conditions but also have such a wide tolerance to assembly errors, tracking, temperature and spectral variations, that the energy generated by them throughout the year is maximized. One of the first addressed issues is the design of secondary optical elements whose primary optics is a Fresnel lens and which, for a fixed concentration, allow an increased acceptance angle and tolerance of the system. Several reflective and refractive secondaries have been designed and analyzed using ray tracing. In particular, using nonimaging optics and based on the single-stage design known as ‘dielectric totally internally reflecting concentrator’, a secondary with square output has been designed, fabricated and characterized. Used together with a Fresnel lens, the secondary can simultaneously achieve high efficiency, concentration and acceptance. Furthermore, an alternative method has been proposed and prototyped for the fabrication of the secondary named dome. The optics is manufactured by direct overmolding of silicone over the solar cells. One characteristic that permeates all the work done in this thesis is the holistic approach in the design of CPV modules, meaning that special attention has been paid to the joint design of the solar cell and the optics to ensure that the total system achieves the highest attainable efficiency. In this regard, many optical systems developed in the thesis have been designed, characterized and optimized considering that the current matching among the subcells within the multijunction solar cell beneath the optics must be close to one. Antireflective coating over the cell acts, somehow, as an interface between the optics and the cell. Consequently, a method has been designed to optimize antireflective coatings that takes into account not only the broad wavelength range that multijunction solar cells are sensitive to but also the angular intensity distribution created by the concentrating optics. In addition, the issue of non-uniformity has also been addressed by comparing the spectral and spatial distributions of irradiance created by different optics (simulated by ray tracing and photographed) and the efficiency losses experienced by cells illuminated by those concentrating optics experimentally determined. The effect of temperature on the concentrating optics has also been studied in this thesis. In particular, finite element simulations have been use to analyze the deformations experienced by the facets of hybrid (silicon-glass) Fresnel lenses, the change of refractive index with temperature and the influence of both effects on the system performance. A model has been implemented which take into consideration atmospheric variations, mainly temperature and spectral content of the direct normal irradiance, as well as thermal and spectral sensitivities of systems, with the aim of maximizing the energy harvested by a CPV module throughout the year in a particular location. Chapters 5 and 6 of this book are devoted to the design, fabrication, and characterization of a new concentrator concept named FluidReflex and based on a single-stage reflective optics with fluid dielectric. In this new concept, the presence of the fluid provides some significant advantages such as: an increased concentration acceptance angle product (CAP) achievable by surrounding the cell with a medium whose refractive index is greater than one, an improvement of the optical efficiency by reducing losses due to Fresnel reflection at several interfaces, an improvement in heat dissipation as the heat concentrated near the cell is transmitted by natural convection and conduction in the fluid, and an improved electrical insulation. By fabricating and characterizing several elementary-unit prototypes it was shown that there is no fundamental reason that prevents the practical implementation of this theoretical concept reaching high efficiency. Several fluid candidates were investigated proving the existence of at least to fluids that meet all the requirements (including the stability under concentrated light) to become part of the FluidReflex concentrator. Finally, several pre-industrial FluidReflex module prototypes have been designed and fabricated. An optimization process for the manufacturing of the multicavity optics was necessary to attain such an optics quality as the one achieved by the single unit. The module prototypes have been measured, both indoors and outdoors, analyzing the current matching of the solar cells beneath the concentrator for different spectral distribution of the incident irradiance. Additionally, the module showed an excellent thermal performance.

Estructuras de hormigón sometidas a deformaciones impuestas

En este trabajo se estudia el comportamiento de elementos de hormigón estructural sometidos a deformaciones impuestas, en donde la fisuración se produce, general y fundamentalmente, en secciones en las que los esfuerzos, y por tanto las fisuras, debidos a las acciones exteriores son inexistentes o prácticamente inexistentes. Se estudia la fisuración producida, a edades tempranas, por las deformaciones debidas al calor de fraguado y la retracción autógena y a lo largo de la vida útil de la estructura, por los fenómenos reológicos del hormigón y la variación de temperatura. Se han estudiado los resultados experimentales disponibles. Se ha identificado un programa de elementos finitos, ATENA, que permite el estudio de fenómenos no lineales, especialmente la fisuración. Se ha contrastado este programa de forma detallada con los resultados experimentales disponibles. Se han estudiado las fuentes de deformaciones impuestas para su caracterización: calor de hidratación, fenómenos reológicos del hormigón y variaciones de temperaturas. Se han estudiado la evolución de las características mecánicas y de conductividad térmica del hormigón a lo largo del tiempo y especialmente a edades tempranas. Se han estudiado las evidencias experimentales disponibles y las estimaciones teóricas propuestas. Se han estudiado las recomendaciones de la CIRIA 660, que es el documento de referencia para tener en cuenta, desde el punto de vista del proyecto este tipo de fenómenos. Con el programa ATENA se han estudiado distintos ejemplos. El objetivo principal ha sido, por un lado, entender el fenómeno y, por otro, contrastar los criterios expuestos en la CIRA 660. This document studies the behavior of structural reinforced concrete structures subjected to imposed deformations, where the cracking occurs, mostly, in sections where the stresses, and thus cracks, due to the exterior actions are none-existent or practically zero. Early age cracking produced by the imposed deformations, due to heat of hydration and autogenous shrinkage, and long term effects, due to the drying shrinkage and the thermal variation, is studied. Available experimental results have been studied. A finite element program has been identified, ATENA, which allows the study of non-lineal phenomena, especially cracking. This software has been contrasted with the available experimental results. The imposed deformation sources have been studied for its characteristically sources: heat of hydration, concrete shrinkage and thermal variation. Also the evolution of the mechanical properties and the thermal conductivity of concrete in time and especially at early age have also been studied. The available experimental evidence and the proposal theoretical estimates have been studied. CIRIA 660 recommendations have been studied, which is the reference document to take into account, form the design point of view, this kind of phenomena. With the ATENA software different examples have been studied. On the one hand the main objective has been to understand the phenomena and, on the other, to contrast the CIRIA 660 proposal.

Un elemento junta de contacto entre suelo y estructura

Se desarrolla un elemento finito especial que permite la conexión entre un elemento viga (con flexión) y un elemento cuadrilátero de extensión (tensión o deformación plana). Asimismo, este elemento especial sirve para modelizar la discontinuidad existente en el contacto entre diferentes materiales. La formulación en movimientos del elemento, permite su inserción directa en un programa general de elementos finitos y, de esta forma, calcular estructuras de hormigón en contacto con el suelo (muros de contención, cimentación, estructuras enterradas como túneles y tuberías, etc.), tanto en el rango elástico como en el elastoplástico, utilizando una potencia limitada en medios computacionales.

Sobre el establecimiento de estimadores e indicadores para los métodos de contorno p-adaptables en problemas de potencial

El establecimiento de métodos numéricos capaces de garantizar una cota de exactitud previamente establecida y de autoadaptar la discretización inicial hasta conseguirlo es una de las tendencias más atractivas en la actualidad. En artículos previos se ha mostrado cómo se puede conseguir ésto con el método de los elementos de contorno y aquí se desarrollan nuevas ideas, inspiradas en tratamientos paralelos con el método de los elementos finitos, para fijar los criterios de adaptaci6n y estimación.

Thermal, stress, and traps effects in AlGaN/GaN HEMTs

GaN y AlN son materiales semiconductores piezoeléctricos del grupo III-V. La heterounión AlGaN/GaN presenta una elevada carga de polarización tanto piezoeléctrica como espontánea en la intercara, lo que genera en su cercanía un 2DEG de grandes concentración y movilidad. Este 2DEG produce una muy alta potencia de salida, que a su vez genera una elevada temperatura de red. Las tensiones de puerta y drenador provocan un stress piezoeléctrico inverso, que puede afectar a la carga de polarización piezoeléctrica y así influir la densidad 2DEG y las características de salida. Por tanto, la física del dispositivo es relevante para todos sus aspectos eléctricos, térmicos y mecánicos. En esta tesis se utiliza el software comercial COMSOL, basado en el método de elementos finitos (FEM), para simular el comportamiento integral electro-térmico, electro-mecánico y electro-térmico-mecánico de los HEMTs de GaN. Las partes de acoplamiento incluyen el modelo de deriva y difusión para el transporte electrónico, la conducción térmica y el efecto piezoeléctrico. Mediante simulaciones y algunas caracterizaciones experimentales de los dispositivos, hemos analizado los efectos térmicos, de deformación y de trampas. Se ha estudiado el impacto de la geometría del dispositivo en su auto-calentamiento mediante simulaciones electro-térmicas y algunas caracterizaciones eléctricas. Entre los resultados más sobresalientes, encontramos que para la misma potencia de salida la distancia entre los contactos de puerta y drenador influye en generación de calor en el canal, y así en su temperatura. El diamante posee une elevada conductividad térmica. Integrando el diamante en el dispositivo se puede dispersar el calor producido y así reducir el auto-calentamiento, al respecto de lo cual se han realizado diversas simulaciones electro-térmicas. Si la integración del diamante es en la parte superior del transistor, los factores determinantes para la capacidad disipadora son el espesor de la capa de diamante, su conductividad térmica y su distancia a la fuente de calor. Este procedimiento de disipación superior también puede reducir el impacto de la barrera térmica de intercara entre la capa adaptadora (buffer) y el substrato. La muy reducida conductividad eléctrica del diamante permite que pueda contactar directamente el metal de puerta (muy cercano a la fuente de calor), lo que resulta muy conveniente para reducir el auto-calentamiento del dispositivo con polarización pulsada. Por otra parte se simuló el dispositivo con diamante depositado en surcos atacados sobre el sustrato como caminos de disipación de calor (disipador posterior). Aquí aparece una competencia de factores que influyen en la capacidad de disipación, a saber, el surco atacado contribuye a aumentar la temperatura del dispositivo debido al pequeño tamaño del disipador, mientras que el diamante disminuiría esa temperatura gracias a su elevada conductividad térmica. Por tanto, se precisan capas de diamante relativamente gruesas para reducer ele efecto de auto-calentamiento. Se comparó la simulación de la deformación local en el borde de la puerta del lado cercano al drenador con estructuras de puerta estándar y con field plate, que podrían ser muy relevantes respecto a fallos mecánicos del dispositivo. Otras simulaciones se enfocaron al efecto de la deformación intrínseca de la capa de diamante en el comportamiento eléctrico del dispositivo. Se han comparado los resultados de las simulaciones de la deformación y las características eléctricas de salida con datos experimentales obtenidos por espectroscopía micro-Raman y medidas eléctricas, respectivamente. Los resultados muestran el stress intrínseco en la capa producido por la distribución no uniforme del 2DEG en el canal y la región de acceso. Además de aumentar la potencia de salida del dispositivo, la deformación intrínseca en la capa de diamante podría mejorar la fiabilidad del dispositivo modulando la deformación local en el borde de la puerta del lado del drenador. Finalmente, también se han simulado en este trabajo los efectos de trampas localizados en la superficie, el buffer y la barrera. Las medidas pulsadas muestran que tanto las puertas largas como las grandes separaciones entre los contactos de puerta y drenador aumentan el cociente entre la corriente pulsada frente a la corriente continua (lag ratio), es decir, disminuir el colapse de corriente (current collapse). Este efecto ha sido explicado mediante las simulaciones de los efectos de trampa de superficie. Por su parte, las referidas a trampas en el buffer se enfocaron en los efectos de atrapamiento dinámico, y su impacto en el auto-calentamiento del dispositivo. Se presenta también un modelo que describe el atrapamiento y liberación de trampas en la barrera: mientras que el atrapamiento se debe a un túnel directo del electrón desde el metal de puerta, el desatrapamiento consiste en la emisión del electrón en la banda de conducción mediante túnel asistido por fonones. El modelo también simula la corriente de puerta, debida a la emisión electrónica dependiente de la temperatura y el campo eléctrico. Además, también se ilustra la corriente de drenador dependiente de la temperatura y el campo eléctrico. ABSTRACT GaN and AlN are group III-V piezoelectric semiconductor materials. The AlGaN/GaN heterojunction presents large piezoelectric and spontaneous polarization charge at the interface, leading to high 2DEG density close to the interface. A high power output would be obtained due to the high 2DEG density and mobility, which leads to elevated lattice temperature. The gate and drain biases induce converse piezoelectric stress that can influence the piezoelectric polarization charge and further influence the 2DEG density and output characteristics. Therefore, the device physics is relevant to all the electrical, thermal, and mechanical aspects. In this dissertation, by using the commercial finite-element-method (FEM) software COMSOL, we achieved the GaN HEMTs simulation with electro-thermal, electro-mechanical, and electro-thermo-mechanical full coupling. The coupling parts include the drift-diffusion model for the electron transport, the thermal conduction, and the piezoelectric effect. By simulations and some experimental characterizations, we have studied the device thermal, stress, and traps effects described in the following. The device geometry impact on the self-heating was studied by electro-thermal simulations and electrical characterizations. Among the obtained interesting results, we found that, for same power output, the distance between the gate and drain contact can influence distribution of the heat generation in the channel and thus influence the channel temperature. Diamond possesses high thermal conductivity. Integrated diamond with the device can spread the generated heat and thus potentially reduce the device self-heating effect. Electro-thermal simulations on this topic were performed. For the diamond integration on top of the device (top-side heat spreading), the determinant factors for the heat spreading ability are the diamond thickness, its thermal conductivity, and its distance to the heat source. The top-side heat spreading can also reduce the impact of thermal boundary resistance between the buffer and the substrate on the device thermal behavior. The very low electrical conductivity of diamond allows that it can directly contact the gate metal (which is very close to the heat source), being quite convenient to reduce the self-heating for the device under pulsed bias. Also, the diamond coated in vias etched in the substrate as heat spreading path (back-side heat spreading) was simulated. A competing mechanism influences the heat spreading ability, i.e., the etched vias would increase the device temperature due to the reduced heat sink while the coated diamond would decrease the device temperature due to its higher thermal conductivity. Therefore, relative thick coated diamond is needed in order to reduce the self-heating effect. The simulated local stress at the gate edge of the drain side for the device with standard and field plate gate structure were compared, which would be relevant to the device mechanical failure. Other stress simulations focused on the intrinsic stress in the diamond capping layer impact on the device electrical behaviors. The simulated stress and electrical output characteristics were compared to experimental data obtained by micro-Raman spectroscopy and electrical characterization, respectively. Results showed that the intrinsic stress in the capping layer caused the non-uniform distribution of 2DEG in the channel and the access region. Besides the enhancement of the device power output, intrinsic stress in the capping layer can potentially improve the device reliability by modulating the local stress at the gate edge of the drain side. Finally, the surface, buffer, and barrier traps effects were simulated in this work. Pulsed measurements showed that long gates and distances between gate and drain contact can increase the gate lag ratio (decrease the current collapse). This was explained by simulations on the surface traps effect. The simulations on buffer traps effects focused on illustrating the dynamic trapping/detrapping in the buffer and the self-heating impact on the device transient drain current. A model was presented to describe the trapping and detrapping in the barrier. The trapping was the electron direct tunneling from the gate metal while the detrapping was the electron emission into the conduction band described by phonon-assisted tunneling. The reverse gate current was simulated based on this model, whose mechanism can be attributed to the temperature and electric field dependent electron emission in the barrier. Furthermore, the mechanism of the device bias via the self-heating and electric field impact on the electron emission and the transient drain current were also illustrated.

Los efectos de borde en problemas geotécnicos estáticos analizados por el M.E.F

En el presente artículo se plantea el problema de las condiciones de borde que deben ser consideradas, para los límites laterales de un modelo de elementos finitos, en el caso de que se pretenda analizar un problema geotécnico en condiciones estáticas. Partiendo de las soluciones existentes para el caso de problemas dinámicos de interacción suelo-estructura, los autores demuestran que puede deducirse un nuevo borde aplicable a las condiciones estáticas. Dicho borde se materializa en una matriz de rigidez que interrelaciona entre sí todos los nudos del lado que se considera. La bondad de la solución propuesta, así como la correspondiente a otras condiciones de borde, es estudiada en un problema concreto que dispone de una solución, estableciéndose unas conclusiones particulares respecto al comportamiento de las diferentes condiciones de borde analizadas.

Problemas de homogeneización en la ingeniería. Una experimentación numérica

Se expone brevemente el problema matemático de homogeneización, en sus facetas: homogeneización en el dominio y en el contorno. Con respecto a esta última, se introducen los conceptos de talla crítica y problema límite. Los resultados obtenidos se aplican a un caso particular representado por un cilindro con condiciones de borde periódicas. Se comprueba la validez del desarrollo teórico mediante una experimentación numérica utilizando elementos finitos 3-D observándose la velocidad de convergencia que se obtiene en relación con el caso límite. Finalmente se comenta la posibilidad de utilizar la técnica de homogeneización para obtener unos criterios hacia un diseño eficiente y óptimo de los resultados de las barras corrugadas en el hormigón armado.

Aplicación de algoritmos paramétricos de generación formal en el diseño estructural. La Ópera de Sídney desde una perspectiva contemporánea.

El importante desarrollo tecnológico e industrial surgido especialmente durante la segunda mitad del siglo pasado ha eliminado las históricas limitaciones técnicas en el ámbito de los pro-­? yectos arquitectónicos, desembocando en la situación actual en la que cualquier planteamiento formal puede ser analizado desde un punto de vista estructural, concluyéndose por tanto que ha desaparecido la barrera del análisis en lo que al desarrollo de un proyecto arquitectónico se refiere. En la actualidad, al igual que a finales del siglo XIX, nos encontramos en un periodo de transi-­? ción, y también, como entonces, es la tecnología la que orienta el cambio. No la tecnología de los nuevos materiales (hormigón y acero) como sucedía tras la revolución industrial sino que es la nueva tecnología digital aplicada a los sistemas de diseño, cálculo y fabricación la que están siendo el motor de la actual transformación. Hoy no es tanto el paradigma mecanicista el que prevalece en muchos casos en la concepción de los edificios sino que, nuevos elementos como la tecnología digital integrada está cambiando la forma de diseñar y concebir el entorno cons-­? truido. Ante este contexto cabría plantearse las siguientes cuestiones: ¿Puede el diseño paramétrico y la tecnología CAD-­?CAM-­?CAE en conjunción con los programas actuales de análisis estructural basados en el Método de los Elementos Finitos hacer más sencilla la construcción de estructu-­? ras ligeras y eficientes hoy en día? ¿Puede la tecnología digital ayudar a ampliar el abanico for-­? mal a la hora de diseñar edificios y a la vez permitir el uso de sistemas estructurales racionales que optimicen el consumo de materiales bajo dichas circunstancias?

Análisis del amortiguamiento en vidrio estructural bajo cargas de impacto

Los últimos avances tecnológicos y las tendencias arquitectónicas actuales han impulsado notablemente el uso de una gran variedad de productos de vidrio en la construcción que, en función de sus propias características, permiten diseñar y calcular elementos estructurales de este material bajo condiciones de seguridad. Este trabajo presenta la evaluación y el análisis de las propiedades de amortiguamiento de placas rectangulares de vidrio laminado de 1938 mm x 876 mm con distinto espesor, según el número de capas intermedias de PVB dispuestas. Mediante simulación numérica por elementos finitos y su posterior comprobación experimental, utilizando análisis modal, se determinaron las frecuencias naturales y el amortiguamiento de las placas de vidrio, tanto en condiciones de contorno libre, como en condiciones operacionales del equipo de ensayos de impacto utilizado en el programa experimental, siguiendo las especificaciones de la norma europea UNE-EN 12600:2003.

Obtención de la tensión de fractura del vidrio en ensayos de anillos concéntricos utilizando MEF

Las discrepancias entre las funciones de distribución acumulada de la tensión de fractura del vidrio es una de las dificultades encontradas para incorporar la filosofía probabilista para el diseño del vidrio. Se utilizan tres ensayos para la caracterización de la tensión de fractura del vidrio: flexión de viga en cuatro puntos, ensayo con dobles anillos concéntricos con pequeña superficie y con gran superficie. La comparación directa de las funciones de distribución presenta discrepancias claras relacionadas con el estado de tensión -uni o biaxial- y con el comportamiento de la estructura -grandes o pequeños desplazamientos. Con un procedimiento iterativo, se ajustan los resultados de los ensayos (carga - presión- probabilidad de rotura) a través de un modelo de elementos finitos que permite estimar las aéreas efectivas de cada espécimen. Los parámetros finales para la tensión de rotura mejoran la estimación de la probabilidad de la carga de rotura de cada probeta con errores inferiores al 10%.

Aplicación del M.E.F. al dibujo de curvas

Se presenta una aplicación muy concreta del Método de Elementos Finitos al trazado de carreteras, que puede ser de interés en algunos casos específicos como en el tratamiento del problema de optimización.

Un estudio sobre la exactitud del M.E.F. Aplicación a la barra recta a extensión

El objetivo de este artículo es fundamentalmente didáctico. Corresponde a un resumen de una monografía escrita sobre el tema, como report interno del Departamento de Análisis de las Estructuras de la Escuela Técnica Superior de Ingenieros de Caminos, Canales y Puertos de Santander, e intenta presentar como varía la exactitud -o velocidad de la convergencia- en los resultados obtenidos mediante el método de los elenentos finitos(MEF).Con objeto de reducir al máximo el esfuerzo computacional se ha considerado un problema de continuidad Cº, y únicamente un elemento aislado, concretamente, el caso de una barra aislada de sección variable bajo acción axil. Este ejemplo sencillo no permite deducir conclusiones definitivas y extensibles a otros casos pero sí probablenente, obtener algunas indicaciones sobre la importancia de algunos aspectos de MEF en relación con la exactitud de los resultados.

Aplicación M.E.F. al diseño de presas bóvedas

El problema de la optimización en el diseño de presas - bóveda, planteado en toda su generalidad, corresponde a la definición de este diseño, es decir, de los siguientes parámetros dela presa bóveda: 1) Forma. 2) Altura. 3) Localización. 4) Propiedades. 5) Procedimiento de construcción. 6) Tratamiento de la cimentación y juntas con el terreno. 7) Localización y tipo de obras secundarias (aliviadero, galerías, etc.). 8) Uso de armado o pretensado de modo que se minimice una cierta función mérito,que exprese el coste de la estructura en función de estos parámetros.

Dynamic analysis of payloads and structures with intermediate modal density

La necesidad de desarrollar técnicas para predecir la respuesta vibroacústica de estructuras espaciales lia ido ganando importancia en los últimos años. Las técnicas numéricas existentes en la actualidad son capaces de predecir de forma fiable el comportamiento vibroacústico de sistemas con altas o bajas densidades modales. Sin embargo, ambos rangos no siempre solapan lo que hace que sea necesario el desarrollo de métodos específicos para este rango, conocido como densidad modal media. Es en este rango, conocido también como media frecuencia, donde se centra la presente Tesis doctoral, debido a la carencia de métodos específicos para el cálculo de la respuesta vibroacústica. Para las estructuras estudiadas en este trabajo, los mencionados rangos de baja y alta densidad modal se corresponden, en general, con los rangos de baja y alta frecuencia, respectivamente. Los métodos numéricos que permiten obtener la respuesta vibroacústica para estos rangos de frecuencia están bien especificados. Para el rango de baja frecuencia se emplean técnicas deterministas, como el método de los Elementos Finitos, mientras que, para el rango de alta frecuencia las técnicas estadísticas son más utilizadas, como el Análisis Estadístico de la Energía. En el rango de medias frecuencias ninguno de estos métodos numéricos puede ser usado con suficiente precisión y, como consecuencia -a falta de propuestas más específicas- se han desarrollado métodos híbridos que combinan el uso de métodos de baja y alta frecuencia, intentando que cada uno supla las deficiencias del otro en este rango medio. Este trabajo propone dos soluciones diferentes para resolver el problema de la media frecuencia. El primero de ellos, denominado SHFL (del inglés Subsystem based High Frequency Limit procedure), propone un procedimiento multihíbrido en el cuál cada subestructura del sistema completo se modela empleando una técnica numérica diferente, dependiendo del rango de frecuencias de estudio. Con este propósito se introduce el concepto de límite de alta frecuencia de una subestructura, que marca el límite a partir del cual dicha subestructura tiene una densidad modal lo suficientemente alta como para ser modelada utilizando Análisis Estadístico de la Energía. Si la frecuencia de análisis es menor que el límite de alta frecuencia de la subestructura, ésta se modela utilizando Elementos Finitos. Mediante este método, el rango de media frecuencia se puede definir de una forma precisa, estando comprendido entre el menor y el mayor de los límites de alta frecuencia de las subestructuras que componen el sistema completo. Los resultados obtenidos mediante la aplicación de este método evidencian una mejora en la continuidad de la respuesta vibroacústica, mostrando una transición suave entre los rangos de baja y alta frecuencia. El segundo método propuesto se denomina HS-CMS (del inglés Hybrid Substructuring method based on Component Mode Synthesis). Este método se basa en la clasificación de la base modal de las subestructuras en conjuntos de modos globales (que afectan a todo o a varias partes del sistema) o locales (que afectan a una única subestructura), utilizando un método de Síntesis Modal de Componentes. De este modo es posible situar espacialmente los modos del sistema completo y estudiar el comportamiento del mismo desde el punto de vista de las subestructuras. De nuevo se emplea el concepto de límite de alta frecuencia de una subestructura para realizar la clasificación global/local de los modos en la misma. Mediante dicha clasificación se derivan las ecuaciones globales del movimiento, gobernadas por los modos globales, y en las que la influencia del conjunto de modos locales se introduce mediante modificaciones en las mismas (en su matriz dinámica de rigidez y en el vector de fuerzas). Las ecuaciones locales se resuelven empleando Análisis Estadístico de Energías. Sin embargo, este último será un modelo híbrido, en el cual se introduce la potencia adicional aportada por la presencia de los modos globales. El método ha sido probado para el cálculo de la respuesta de estructuras sometidas tanto a cargas estructurales como acústicas. Ambos métodos han sido probados inicialmente en estructuras sencillas para establecer las bases e hipótesis de aplicación. Posteriormente, se han aplicado a estructuras espaciales, como satélites y reflectores de antenas, mostrando buenos resultados, como se concluye de la comparación de las simulaciones y los datos experimentales medidos en ensayos, tanto estructurales como acústicos. Este trabajo abre un amplio campo de investigación a partir del cual es posible obtener metodologías precisas y eficientes para reproducir el comportamiento vibroacústico de sistemas en el rango de la media frecuencia. ABSTRACT Over the last years an increasing need of novel prediction techniques for vibroacoustic analysis of space structures has arisen. Current numerical techniques arc able to predict with enough accuracy the vibro-acoustic behaviour of systems with low and high modal densities. However, space structures are, in general, very complex and they present a range of frequencies in which a mixed behaviour exist. In such cases, the full system is composed of some sub-structures which has low modal density, while others present high modal density. This frequency range is known as the mid-frequency range and to develop methods for accurately describe the vibro-acoustic response in this frequency range is the scope of this dissertation. For the structures under study, the aforementioned low and high modal densities correspond with the low and high frequency ranges, respectively. For the low frequency range, deterministic techniques as the Finite Element Method (FEM) are used while, for the high frequency range statistical techniques, as the Statistical Energy Analysis (SEA), arc considered as more appropriate. In the mid-frequency range, where a mixed vibro-acoustic behaviour is expected, any of these numerical method can not be used with enough confidence level. As a consequence, it is usual to obtain an undetermined gap between low and high frequencies in the vibro-acoustic response function. This dissertation proposes two different solutions to the mid-frequency range problem. The first one, named as The Subsystem based High Frequency Limit (SHFL) procedure, proposes a multi-hybrid procedure in which each sub-structure of the full system is modelled with the appropriate modelling technique, depending on the frequency of study. With this purpose, the concept of high frequency limit of a sub-structure is introduced, marking out the limit above which a substructure has enough modal density to be modelled by SEA. For a certain analysis frequency, if it is lower than the high frequency limit of the sub-structure, the sub-structure is modelled through FEM and, if the frequency of analysis is higher than the high frequency limit, the sub-structure is modelled by SEA. The procedure leads to a number of hybrid models required to cover the medium frequency range, which is defined as the frequency range between the lowest substructure high frequency limit and the highest one. Using this procedure, the mid-frequency range can be define specifically so that, as a consequence, an improvement in the continuity of the vibro-acoustic response function is achieved, closing the undetermined gap between the low and high frequency ranges. The second proposed mid-frequency solution is the Hybrid Sub-structuring method based on Component Mode Synthesis (HS-CMS). The method adopts a partition scheme based on classifying the system modal basis into global and local sets of modes. This classification is performed by using a Component Mode Synthesis, in particular a Craig-Bampton transformation, in order to express the system modal base into the modal bases associated with each sub-structure. Then, each sub-structure modal base is classified into global and local set, fist ones associated with the long wavelength motion and second ones with the short wavelength motion. The high frequency limit of each sub-structure is used as frequency frontier between both sets of modes. From this classification, the equations of motion associated with global modes are derived, which include the interaction of local modes by means of corrections in the dynamic stiffness matrix and the force vector of the global problem. The local equations of motion are solved through SEA, where again interactions with global modes arc included through the inclusion of an additional input power into the SEA model. The method has been tested for the calculation of the response function of structures subjected to structural and acoustic loads. Both methods have been firstly tested in simple structures to establish their basis and main characteristics. Methods are also verified in space structures, as satellites and antenna reflectors, providing good results as it is concluded from the comparison with experimental results obtained in both, acoustic and structural load tests. This dissertation opens a wide field of research through which further studies could be performed to obtain efficient and accurate methodologies to appropriately reproduce the vibro-acoustic behaviour of complex systems in the mid-frequency range.