Modelización de reacciones químicas en presencia de un campo electromagnético variable

En esta tesis se aborda el estudio del proceso de isomerización del sistema molecular LiNC/LiCN tanto aislado como en presencia de un pulso láser aplicando la teoría del estado de transición (TST). Esta teoría tiene como pilar fundamental el hecho de que el conocimiento de la dinámica en las proximidades de un punto de silla de la superficie de energía potencial permite determinar los parámetros cinéticos de la reacción objeto de estudio. Históricamente, existen dos formulaciones de la teoría del estado de transición, la versión termodinámica de Eyring (Eyr38) y la visión dinámica de Wigner (Wig38). Ésta última ha sufrido recientemente un amplio desarrollo, paralelo a los avances en sistemas dinámicos que ha dado lugar a una formulación geométrica en el espacio de fases que sirve como base al trabajo desarrollado en esta tesis. Nos hemos centrado en abordar el problema desde una visión fundamentalmente práctica, ya que la teoría del estado de transición presenta una desventaja: su elevado coste computacional y de tiempo de cálculo. Dos han sido los principales objetivos de este trabajo. El primero de ellos ha sido sentar las bases teóricas y computacionales de un algoritmo eficiente que permita obtener las magnitudes fundamentales de la TST. Así, hemos adaptado con éxito un algoritmo computacional desarrollado en el ámbito de la mecánica celeste (Jor99), obteniendo un método rápido y eficiente para la obtención de los objetos geométricos que rigen la dinámica en el espacio de fases y que ha permitido calcular magnitudes cinéticas tales como el flujo reactivo, la densidad de estados de reactivos y productos y en última instancia la constante de velocidad. Dichos cálculos han sido comparados con resultados estadísticos (presentados en (Mül07)) lo cual nos ha permitido demostrar la eficacia del método empleado. El segundo objetivo de esta tesis, ha sido la evaluación de la influencia de los parámetros de un pulso electromagnético sobre la dinámica de reacción. Para ello se ha generalizado la metodología de obtención de la forma normal del hamiltoniano cuando el sistema químico es alterado mediante una perturbación temporal periódica. En este caso el punto fijo inestable en cuya vecindad se calculan los objetos geométricos de interés para la aplicación de la TST, se transforma en una órbita periódica del mismo periodo que la perturbación. Esto ha permitido la simulación de la reactividad en presencia de un pulso láser. Conocer el efecto de esta perturbación posibilita el control de la reactividad química. Además de obtener los objetos geométricos que rigen la dinámica en una cierta vecindad de la órbita periódica y que son la clave de la TST, se ha estudiado el efecto de los parámetros del pulso sobre la reactividad en el espacio de fases global así como sobre el flujo reactivo que atraviesa la superficie divisoria que separa reactivos de productos. Así, se ha puesto de manifiesto, que la amplitud del pulso es el parámetro más influyente sobre la reactividad química, pudiendo producir la aparición de flujos reactivos a energías inferiores a las de aparición del sistema aislado y el aumento del flujo reactivo a valores constantes de energía inicial. ABSTRACT We have studied the isomerization reaction LiNC/LiCN isolated and perturbed by a laser pulse. Transition State theory (TST) is the main tool we have used. The basis of this theory is knowing the dynamics close to a fixed point of the potential energy surface. It is possible to calculate kinetic magnitudes by knowing the dynamics in a neighbourhood of the fixed point. TST was first formulated in the 30's and there were 2 points of view, one thermodynamical by Eyring (Eyr38) and another dynamical one by Wigner (Wig38). The latter one has grown lately due to the growth of the dynamical systems leading to a geometrical view of the TST. This is the basis of the work shown in this thesis. As the TST has one main handicap: the high computational cost, one of the main goals of this work is to find an efficient method. We have adapted a methodology developed in the field of celestial mechanics (Jor99). The result: an efficient, fast and accurate algorithm that allows us to obtain the geometric objects that lead the dynamics close to the fixed point. Flux across the dividing surface, density of states and reaction rate coefficient have been calculated and compared with previous statistical results, (Mül07), leading to the conclusion that the method is accurate and good enough. We have widen the methodology to include a time dependent perturbation. If the perturbation is periodic in time, the fixed point becomes a periodic orbit whose period is the same as the period of the perturbation. This way we have been able to simulate the isomerization reaction when the system has been perturbed by a laser pulse. By knowing the effect of that perturbation we will be able to control the chemical reactivity. We have also studied the effect of the parameters on the global phase space dynamics and on the flux across the dividing surface. It has been prove that amplitude is the most influent parameter on the reaction dynamics. Increasing amplitude leads to greater fluxes and to some flux at energies it would not if the systems would not have been perturbed.

Sistema de control para optimizar la eficiencia de un equipo fotovoltaico de bombeo directo accionado por un motor de inducción

El objetivo principal está recogido en el título de la Tesis. Ampliando éste para hacerlo más explícito, puede decirse que se trata de “desarrollar un sistema de control para que una instalación fotovoltaica de bombeo directo con una bomba centrífuga accionada por un motor de inducción trabaje de la forma más eficiente posible”. Para lograr ese propósito se establecieron los siguientes objetivos específicos: 1. Diseñar y construir un prototipo de instalación fotovoltaica de bombeo directo que utilice principalmente elementos de bajo coste y alta fiabilidad. Para cumplir esos requisitos la instalación consta de un generador fotovoltaico con módulos de silicio monocristalino, una bomba centrífuga accionada por un motor de inducción y un inversor que controla vectorialmente el motor. Los módulos de silicio monocristalino, el motor asíncrono y la bomba centrífuga son, en sus respectivas categorías, los elementos más robustos y fiables que existen, pudiendo ser adquiridos, instalados e incluso reparados (el motor y la bomba) por personas con una mínima formación técnica en casi cualquier lugar del mundo. El inversor no es tan fiable ni fácil de reparar. Ahora bien, para optimizar la potencia que entrega el generador y tener algún tipo de control sobre el motor se necesita al menos un convertidor electrónico. Por tanto, la inclusión del inversor en el sistema no reduce su fiabilidad ni supone un aumento del coste. La exigencia de que el inversor pueda realizar el control vectorial del motor responde a la necesidad de optimizar tanto la operación del conjunto motor-bomba como la del generador fotovoltaico. Como más adelante se indica, lograr esa optimización es otro de los objetivos que se plantea. 2. Reducir al mínimo el número de elementos de medida y control que necesita el sistema para su operación (sensorless control). Con ello se persigue aumentar la robustez y fiabilidad del sistema y reducir sus operaciones de mantenimiento, buscando que sea lo más económico posible. Para ello se deben evitar todas las medidas que pudieran ser redundantes, tomando datos sólo de las variables eléctricas que no pueden obtenerse de otra forma (tensión e intensidad en corriente continua y dos intensidades en corriente alterna) y estimando la velocidad del rotor (en vez de medirla con un encoder u otro dispositivo equivalente). 3. Estudiar posibles formas de mejorar el diseño y la eficiencia de estas instalaciones. Se trata de establecer criterios para seleccionar los dispositivos mas eficientes o con mejor respuesta, de buscar las condiciones para la operación óptima, de corregir problemas de desacoplo entre subsistemas, etc. Mediante el análisis de cada una de las partes de las que consta la instalación se plantearán estrategias para minimizar pérdidas, pautas que permitan identificar los elementos más óptimos y procedimientos de control para que la operación del sistema pueda alcanzar la mayor eficiente posible. 4. Implementar un modelo de simulación del sistema sobre el que ensayar las estrategias de control que sean susceptibles de llevar a la práctica. Para modelar el generador fotovoltaico se requiere un conjunto de parámetros que es necesario estimar previamente a partir de datos obtenidos de los catálogos de los módulos a utilizar o mediante ensayos. Igual sucede con los parámetros para modelar el motor. Como se pretende que el motor trabaje siempre con la máxima eficiencia será necesario realizar su control vectorial, por lo que el modelo que se implemente debe ser también vectorial. Ahora bien, en el modelo vectorial estándar que normalmente se utiliza en los esquemas de control se consideran nulas las pérdidas en el hierro, por lo que sólo se podrá utilizar ese modelo para evaluar la eficiencia del motor si previamente se modifica para que incluya el efecto de dichas pérdidas. 5. Desarrollar un procedimiento de control para que el inversor consiga que el motor trabaje con mínimas pérdidas y a la vez el generador entregue la máxima potencia. Tal como se ha mencionado en el primer objetivo, se trata de establecer un procedimiento de control que determine las señales de consigna más convenientes para que el inversor pueda imponer en cada momento al motor las corrientes de estator para las que sus pérdidas son mínimas. Al mismo tiempo el procedimiento de control debe ser capaz de variar las señales de consigna que recibe el inversor para que éste pueda hacer que el motor demande más o menos potencia al generador fotovoltaico. Actuando de esa forma se puede lograr que el generador fotovoltaico trabaje entregando la máxima potencia. El procedimiento de control desarrollado se implementará en un DSP encargado de generar las señales de referencia que el inversor debe imponer al motor para que trabaje con mínimas pérdidas y a la vez el generador fotovoltaico entregue la máxima potencia.

Levantamiento batimétrico en el islote Montaña Clara (Islas Canarias).

Los objetivos del presente Proyecto Fin de Carrera, serán el levantamiento cartográfico en detalle de la zona de trabajo propuesta, mediante la utilización del Buque Oceanográfico Emma Bardán y su sistema de ecosonda multihaz de alta resolución. Los trabajos asociados serán:  El reconocimiento batimétrico sistemático del Islote Montaña Clara con recubrimiento al 100% (ecosonda Multihaz EM-3002D).  Toma de datos sistemáticos de perfiles de velocidad del sonido en el agua. Los objetivos específicos del presente Proyecto Fin de Carrera son proporcionar documentación cartográfica que contribuya a la investigación del fondo marino de la plataforma continental del Islote Montaña Clara, tanto por métodos indirectos como directos. Las fases de trabajo son:  Procesar y tratar los datos obtenidos en la campaña, para la obtención de mapas,cartas e informes.  Elaborar una carta batimétrica basada en datos multihaz con cobertura al 100% del área propuesta.  Realizar el tratamiento informático necesario para la integración de los datos en un Sistema de Información Geográfica (S.I.G.).

Estudio numérico de flujo turbulento cargado con partículas sólidas a través de canales y tuberías de sección variable

La presente Tesis Doctoral tiene como objetivo el estudio de flujo turbulento cargado con partículas sólidas a través de canales y tuberías de sección constante usando un enfoque Euleriano-Lagrangiano. El campo de flujo de la fase de transporte (aire) se resuelve usando simulación de grandes escalas (LES), implementada en un programa de volúmenes finitos mientras que las ecuaciones gobernantes de la fase dispersa son resueltas por medio de un algoritmo de seguimiento Lagrangiano de partículas que ha sido desarrollado y acoplado al programa que resuelve el flujo. Se estudia de manera sistemática y progresiva la interacción fluido→partícula (one-way coupling), a través de diferentes configuraciones geométricas en coordenadas cartesianas (canales de sección constante y variable) y en coordenadas cilíndricas (tuberías de sección constante y sección variable) abarcando diferentes números de Reynolds y diferentes tamaños de partículas; todos los resultados obtenidos han sido comparados con datos publicados previamente. El estudio de flujo multifásico a través de, tuberías de sección variable, ha sido abordada en otras investigaciones mayoritariamente de forma experimental o mediante simulación usando modelos de turbulencia menos complejos y no mediante LES. El patrón de flujo que se verifica en una tubería con expansión es muy complejo y dicha configuración geométrica se halla en múltiples aplicaciones industriales que involucran el transporte de partículas sólidas, por ello es de gran interés su estudio. Como hecho innovador, en esta tesis no solo se resuelven las estadísticas de velocidad del fluido y las partículas en tuberías con diferentes tamaños de expansión y diferentes regímenes de flujo sino que se caracteriza, usando diversas formulaciones del número de Stokes y el parámetro de arrastre, el ingreso y acumulación de partículas dentro de la zona de recirculación, obteniéndose resultados coincidentes con datos experimentales. ABSTRACT The objective of this Thesis research is to study the turbulent flow laden with solid particles through channels and pipes with using Eulerian-Lagrangian approach. The flow field of the transport phase (air ) is solved using large eddy simulation ( LES ) implemented in a program of finite volume while the governing equations of the dispersed phase are resolved by means of a particle Lagrangian tracking algorithm which was developed and coupled to principal program flow solver . We studied systematically and progressively the fluid interaction → particle ( one- way coupling ) , through different geometric configurations in Cartesian coordinates ( channel with constant and variable section) and in cylindrical coordinates ( pipes with constant section and variable section ) covering different Reynolds numbers and different particle sizes, all results have been compared with previously published data . The study of multiphase flow through, pipes with variable section has been addressed in other investigations predominantly experimentally or by simulation using less complex models and no turbulence by LES. The flow pattern is verified in a pipe expansion is very complex and this geometry is found in many industrial applications involving the transport of solid particles, so it is of great interest to study. As an innovator fact , in this Thesis not only finds fluid velocity statistics and particles with different sizes of pipe expansion and different flow regimes but characterized, using various formulations of the Stokes number and the drag parameter are resolved, the entry and accumulation of particles within the recirculation zone , matching results obtained with experimental data.

Enfoque al ahorro de energía en la conducción de ferrocarriles utilizando el método múltiple de restricción de la optimización

El ahorro de la energía consumida no es un tema menor en el sector ferroviario. Las infraestructuras ferroviarias y los planes operativos implicados en distintos sistemas como el metropolitano, trenes locales, larga distancia y alta velocidad, deben ser realizados teniendo en cuenta las múltiples alternativas para alcanzar un sistema de optimización en el ahorro de energía. Para conseguir alcanzar la mejor solución final, las alternativas deben recoger y trabajar con el máximo de variables asociadas al movimiento del tren y a los planes de operación. El método presentado aquí usa una solución semianalítica que conlleva una discretización del problema, para posteriormente aplicar el método de los Multiplicadores de Lagrange para resolver la optimización de las n-tuplas de velocidad. Con esta alternativa propuesta, es posible incluir todos los detalles sobre las operaciones del tren, como restricciones en el horario o tipos de frenado.

Análisis del efecto de las actividades físicas grupales en pacientes con daño cerebral adquirido en fase subaguda

Introducción. El daño cerebral adquirido (DCA) hace referencia a cualquier tipo de lesión no degenerativa que se produce en el cerebro. Las actividades físicas grupales (AFG) se presentan como un tratamiento efectivo para la mejora de la capacidad funcional. Objetivo. Analizar la eficacia de un programa de AFG en personas con DCA en fase subaguda para su integración física en la comunidad. Pacientes y métodos. Treinta y tres pacientes con DCA, con una edad de 33,18 ± 10,39 años, participaron en un programa de AFG (talleres de circuito, equilibrio simple, equilibrio dual, desplazamiento dual y actividades físico-deportivas) de 10 semanas. Al comenzar y concluir el programa se evaluaron las variables de velocidad (prueba de velocidad de la marcha en 10 metros), resistencia (prueba de marcha de seis minutos), equilibrio dinámico (Step Test), capacidad funcional(Timed Up & Go), escala de percepción de seguridad (Activities-specific Balance Confidence Scale) y Physical Activity and Disability Survey (PADS), el promedio por hora de la intensidad de la actividad y el número de pasos fuera del centro de rehabilitación (usando monitores de actividad física). La prueba t para muestras relacionadas se utilizó para evaluar las diferencias en las variables. Resultados. Se hallaron diferencias significativas (p ≤ 0,05) en las variables de velocidad, resistencia, equilibrio, capacidad funcional, percepción de seguridad, percepción de realización de actividad general (pregunta 3 del PADS) y número de pasos. Conclusión. Los programas de AFG mejoran las capacidades físicas, percepción de seguridad, realización de actividad en general y número de pasos, lo que puede conllevar una mayor participación en la comunidad.

Arquitectura de Radios Wake-up para redes de sensores inalámbricas basada en FPGA

En este artículo se muestra la implementación de una Wake-up radio para nodos de redes de sensores inalámbricas basada en FPGAS de ultra bajo consumo. El objetivo principal es evaluar la utilización de dispositivos de lógica programable para realizar el procesamiento de los mensajes y explotar su velocidad de, flexibilidad y bajo consumo comparado con las implementaciones más tradicionales basadas en ASIC o microcontroladores.

Algoritmo para el análisis del Comportamiento de fresas periféricas en condiciones variables

La simulación de procesos de mecanizado supone hoy en día una herramienta de gran interés para predecir el comportamiento de la herramienta en las operaciones de corte y por lo tanto, la posibilidad de optimizar estas operaciones, permitiendo alcanzar una mayor productividad en los procesos de fabricación. Los algoritmos que actualmente se están utilizando para la predicción de fuerzas de corte son variados y su eficiencia diferente. La mayoría de los algoritmos desarrollados se centran en determinar la evolución de las fuerzas de corte en cada vuelta de la herramienta sin variar las condiciones de mecanizado en este intervalo. En este sentido, se ha desarrollado un algoritmo para fresado periférico basado en el espesor de viruta medio, que permite simplificar el algoritmo de estimación de fuerzas y ser ejecutado con una mayor velocidad, manteniendo el mismo nivel de precisión en la estimación. Con este nuevo modelo es posible realizar la estimación de fuerzas en fresado no solo cuando las condiciones de corte son uniformes sino también cuando se producen cambios en la configuración del corte.

La productividad en la carga, la resistencia de la roca y la energía del explosivo

Twenty production blasts in two open pit mines were monitored, in rocks with medium to very high strength. Three different blasting agents (ANFO, watergel and emulsion blend) were used, with powder factors ranging between 0.88 and 1.45 kg/m3 . Excavators were front loaders and rope shovels. Mechanical properties of the rock, blasting characteristics and mucking rates were carefully measured. A model for the calculation of the productivity of excavators is developed thereof, in which the production rate results as a product of an ideal, maximum, productivity rate times an operating efficiency. The maximum rate is a function of the dipper capacity and the efficiency is a function of rock density, strength, and explosive energy concentration in the rock. The model is statistically significant and explains up to 92 % of the variance of the production rate measurements. RESUMEN Se han analizado veinte voladuras de producción en dos minas a cielo abierto, en rocas con resistencia entre media y muy alta. Se usaron tres agentes de voladura diferentes (ANFO, hidrogel y emulsión mezcla), con consumos específicos entre 0.88 y 1.45 kg/m3 . Las excavadoras eran cargadoras frontales y palas de cables. Se midió cuidadosamente las propiedades de la roca, las características de las voladuras y las velocidades de carga, a partir de las cuales se desarrolla un modelo para el cálculo de la productividad de las excavadoras en el que ésta es función de una productividad máxima o ideal multiplicada por una eficacia operativa. La producción máxima es función de la capacidad del cazo y la eficiencia es función de la densidad y resistencia de la roca, y de la concentración energética del explosivo en la roca. El modelo tiene significación estadística y explica un 92 % de la varianza de las medidas de velocidad.

High speed railroad bridges. Conception, bases of design, typological possibilities and construction methods

Diseño conceptual de puentes de alta velocidad ferroviarios. Railroad bridges, in general, and those for high speed railways, in particular, demand very special conditions. The traffic loads are much higher than for road bridges. Loads due to braking and acceleration determine, due to their magnitude, the structural layout. Because of the speed of the vehicles there are specific dynamic effects which need to be considered. In order to ensure passenger comfort, compatible with speeds of up to 350 km/h, it is necessary to meet very demanding conditions with respect to stiffness, displacements and dynamic behavior. In this paper these conditions are briefly described and different typological possibilities to satisfy them are presented as well as the main construction methods applicable to this kind of bridges.

Influencia del espesor de la junta de mortero en muros de ladrillo sometidos a esfuerzos de compresión

En este trabajo se analiza la influencia del espesor de la junta de mortero en la resistencia a compresión de la fábrica de ladrillo cerámico bajo esfuerzos de compresión uniaxial. Además se estudia el progresivo agrietamiento del material durante el proceso de carga. Se han ensayado nueve muros y nueve pilares, combinando tres espesores de junta de mortero relizando medidas con ultrasonidos, además se han realizado medidas con extensómetros y esclerómetro. Los resultados obtenidos muestran que a menor espesor de junta mayor carga de rotura y mayor valor de la velocidad de pulso ultrasónico. La medida con ultrasonidos muestra una buena correlación con las medidas extensométricas y permite una eficaz detección del agrietamiento interior del material durante el proceso de rotura

Herbigacion en cítricos a través de sistemas de riego por goteo: factores que influyen en el comportamiento y distribución de herbicidas en el bulbo húmedo

Existe una serie de características que diferencian malherbologicamente las plantaciones con sistemas de riego por goteo de las tradicionales. Asimismo, el comportamiento de los herbicidas en el suelo tiene unas determinadas peculiaridades inherentes al tipo de riego, que hay que tener en cuenta a la hora de utilizarlos correctamente en estos sistemas. El objetivo fundamental de este trabajo ha sido estudiar como influyen algunos factores en ese comportamiento. Para ello, se ha realizado una serie de ensayos de herbigación en campo, con diferentes productos, diferentes momentos y tiempos de incorporación a lo largo del riego, fraccionamientos de dosis, secuencias de incorporaciones o bajo diferentes intensidades de riego. Los trabajos de campo se completan en laboratorio con el análisis de las muestras de suelo, tomadas en diferentes prospecciones de los bulbos húmedos, mediante bioensayos con plantas indicadoras. Paralelamente a los trabajos de herbigación, se realizan secuencias de ensayos de herbicidas en aplicación tradicional, con productos remanentes en pulverización superficial y programas a base de herbicidas exclusivamente foliares, estudiando su evolución a lo largo de varios años. De todas estas experiencias se han podido extraer conclusiones sobre la influencia de distintos factores en el comportamiento de herbicidas en parcelas con sistemas de riego por goteo, como son las características de los productos, fraccionamiento de dosis, momento de la incorporación a lo largo del riego o la época del cmo/intensidad de riegos. Así, hay productos, como el oxifluorfén, que aplicados en herbigación, muestran una escasa movilidad en el bulbo húmedo, quedando retenidos en la parte central del mismo, mientras otros tienen un desplazamiento constante, más o menos rápido, hacia la periferia del bulbo, como sucede con terbacilo, terbumetona+terbitilazina, norflurazona y, en menor medida, con azafenidin. En estos casos, para un determinado herbicida, la velocidad de desplazamiento, y por lo tanto su efecto en el control de hierbas, depende en gran medida, de la intensidad de riegos. La solubilidad y el coeficiente de adsorción Koc, son características de los herbicidas que, en gran medida, determinan ese comportamiento en el bulbo húmedo. Excepto para los productos muy poco móviles, la dosificación en herbigación ha de ser fraccionada, siendo la intensidad de riegos uno de los factores clave para optimizar la frecuencia de incorporaciones. Junto a esa intensidad de riegos, el momento y tiempo de la incorporación determinan también los resultados. Mientras el desplazamiento de algunos productos depende casi exclusivamente del riego de incorporación, como en el oxifluorfen y azafenidin, en otros casos se ve también muy influenciado por la intensidad dé los riegos posteriores, como sucede especialmente con los más móviles: terbacilo, la mezcla terbutilazina + terbumetona o norflurazona. En las aplicaciones tradicionales de herbicidas remanentes, puede comprobarse como los resultados obtenidos con algunos productos, en parcelas con riego por goteo, son muy diferentes a las de riego tradicional Así, terbacilo o norflurazona, reducen extraordinariamente su persistencia activa en la zona central del bulbo, especialmente en épocas con elevadas intensidades de riego, mientras la alargan fuera de esas zonas de influencia de los goteros. En base a las conclusiones de estos trabajos se apunta una serie de posibles mejoras para optimizar los programas de control de hierbas en plantaciones de cítricos con riego por goteo, que incluyen la herbigación, con combinación de productos de diferentes comportamientos en el bulbo húmedo, la aplicación tradicional de herbicidas remanentes, programas con foliares y combinaciones de estas técnicas.

High fidelity simulations of failure in fiber-reinforced composites

Los ensayos virtuales de materiales compuestos han aparecido como un nuevo concepto dentro de la industria aeroespacial, y disponen de un vasto potencial para reducir los enormes costes de certificación y desarrollo asociados con las tediosas campañas experimentales, que incluyen un gran número de paneles, subcomponentes y componentes. El objetivo de los ensayos virtuales es sustituir algunos ensayos por simulaciones computacionales con alta fidelidad. Esta tesis es una contribución a la aproximación multiescala desarrollada en el Instituto IMDEA Materiales para predecir el comportamiento mecánico de un laminado de material compuesto dadas las propiedades de la lámina y la intercara. La mecánica de daño continuo (CDM) formula el daño intralaminar a nivel constitutivo de material. El modelo de daño intralaminar se combina con elementos cohesivos para representar daño interlaminar. Se desarrolló e implementó un modelo de daño continuo, y se aplicó a configuraciones simples de ensayos en laminados: impactos de baja y alta velocidad, ensayos de tracción, tests a cortadura. El análisis del método y la correlación con experimentos sugiere que los métodos son razonablemente adecuados para los test de impacto, pero insuficientes para el resto de ensayos. Para superar estas limitaciones de CDM, se ha mejorado la aproximación discreta de elementos finitos enriqueciendo la cinemática para incluir discontinuidades embebidas: el método extendido de los elementos finitos (X-FEM). Se adaptó X-FEM para un esquema explícito de integración temporal. El método es capaz de representar cualitativamente los mecanismos de fallo detallados en laminados. Sin embargo, los resultados muestran inconsistencias en la formulación que producen resultados cuantitativos erróneos. Por último, se ha revisado el método tradicional de X-FEM, y se ha desarrollado un nuevo método para superar sus limitaciones: el método cohesivo X-FEM estable. Las propiedades del nuevo método se estudiaron en detalle, y se concluyó que el método es robusto para implementación en códigos explícitos dinámicos escalables, resultando una nueva herramienta útil para la simulación de daño en composites. Virtual testing of composite materials has emerged as a new concept within the aerospace industry. It presents a very large potential to reduce the large certification costs and the long development times associated with the experimental campaigns, involving the testing of a large number of panels, sub-components and components. The aim of virtual testing is to replace some experimental tests by high-fidelity numerical simulations. This work is a contribution to the multiscale approach developed in Institute IMDEA Materials to predict the mechanical behavior of a composite laminate from the properties of the ply and the interply. Continuum Damage Mechanics (CDM) formulates intraply damage at the the material constitutive level. Intraply CDM is combined with cohesive elements to model interply damage. A CDM model was developed, implemented, and applied to simple mechanical tests of laminates: low and high velocity impact, tension of coupons, and shear deformation. The analysis of the results and the comparison with experiments indicated that the performance was reasonably good for the impact tests, but insuficient in the other cases. To overcome the limitations of CDM, the kinematics of the discrete finite element approximation was enhanced to include mesh embedded discontinuities, the eXtended Finite Element Method (X-FEM). The X-FEM was adapted to an explicit time integration scheme and was able to reproduce qualitatively the physical failure mechanisms in a composite laminate. However, the results revealed an inconsistency in the formulation that leads to erroneous quantitative results. Finally, the traditional X-FEM was reviewed, and a new method was developed to overcome its limitations, the stable cohesive X-FEM. The properties of the new method were studied in detail, and it was demonstrated that the new method was robust and can be implemented in a explicit finite element formulation, providing a new tool for damage simulation in composite materials.

Análisis del contenido en frecuencias de vibraciones por voladura

Este proyecto tiene por objeto el análisis de las diferentes frecuencias y tipos de onda presentes en los registros de vibraciones de las voladuras. Se han utilizado 108 registros de 4 barrenos sencillos disparados individualmente y una voladura de producción. Se estudian así mismo la evolución de las frecuencias (en relación a la distancia y a la geología particular de cresta y plaza) y la ley de atenuación del terreno que relaciona la evolución velocidad de partícula en función de la distancia escalada. Para la voladura de producción se analiza si las vibraciones cumplen la normativa vigente según la norma UNE 22.381.93 “Control de Vibraciones Producidas por Voladuras”. Abstract The main objective of this project is to analyze the different types and frequencies of seismic waves found in blasting logs. 108 logs from 4 single shots and a standard production blasting have been used. The evolution of frequencies (regarding distance and different geologies of crest and toe) and the attenuation law which allows us to obtain peak particle velocity regarding distance are also studied. For the production blasting, data is compared to the limits set by current regulations following the standard UNE 22.381.93 “Control de Vibraciones Producidas por Voladuras”.

Mejora y optimización de un código de simulación hidrodinámico con transporte de radiación para el estudio de sistemas en condiciones de alta densidad de energía

Se presentan las mejoras introducidas en un código de transporte de radiación acoplada a la hidrodinámica llamado ARWEN para el estudio de sistemas en el rango de física de alta densidad de energía (High Energy Density Physics). Los desarrollos introducidos se basan en las siguientes áreas: ít>,~ Ecuaciones de estado: se desarrolla una nueva metodología mediante la cual es posible ajustar los resultados de un modelo simple de ecuaciones de estado como QEOS a datos experimentales y resultados de AIMD. Esta metodología tiene carácter general para poder ser aplicada en gran cantidad de materuales de interés y amplia la flexibilidad de ajuste de los métodos de los que ha partido como base este trabajo. En segundo lugar, se ha desarrollado una librería para la gestión de tablas de datos de ecuaciones de estado que también incluye la gestión de tablas con datos de opacidades y de ionización. Esta nueva librería extiende las capacidades de la anterior al tener llamadas más específicas que aceleran los cálculos, y posibilidad de uso de varias tablas a la vez. Solver de difusión: se ha desarrollado un nuevo paquete para resolver la ecuación de difusión que se aplicará a la conducción de calor dentro del plasma. El método anterior no podía ser ejecutado en paralelo y producía resultados dependientes de la resolución de la malla, mientras que este método es paralelizable y además obtiene una solución con mejor convergencia, lo que supone una solución que no depende del refinamiento del mallado. Revisión del paquete de radiación: en primer lugar se ha realizado una revisión de la implementación del modelo de radiación descubriendo varios errores que han sido depurados. También se ha incluido la nueva librería de gestión de tablas de opacidades que permiten la obtención de las propiedades ópticas del plasma en multigrupos de energía. Por otra parte se ha extendido el cálculo de los coeficientes de transporte al esquema multimaterial que ha introducido David Portillo García en el paquete hidrodinámico del código de simulación. Por último se ha revisado el esquema de resolución del transporte y se ha modificado para hacerlo paralelizable. • Se ha implementado un paquete de trazado de rayos para deposición láser que extiende la utilidad del anterior al ser en 3D y poder utilizar entonces diferentes configuraciones. • Una vez realizadas todas estas tareas se ha aplicado el código ARWEN al estudio de la astrofísica de laboratorio simulando los experimentos llevados a cabo en la instalación PALS por Chantal Stehlé acerca de ondas de choque radiativas. Se han comparado los resultados experimentales frente a las predicciones del código ARWEN obteniéndose una gran concordancia en la velocidad de la onda de choque generada y en las dimensiones del precursor. El código de simulación sobre el que se ha trabajado, junto con los desarrollos aportados por otros investigadores durante la realización de esta tesis, ha permitido participar en colaboraciones con laboratorios de Francia o Japón y se han producido resultados científicos publicados basados en el trabajo descrito en esta tesis. ABSTRACT Improvements in radiation hydrodynamic code ARWEN for the study of systems in the range of physics high energy density (High Energy Density Physics) are presented. The developments introduced are based on the following áreas: • Equations of state: a new methodology was developed to adjust the results of a simple Equation of State model like QEOS to experimental data and results of AIMD. This methodology can be applied to a large amount of materials and it increases the flexibility and range of the previous methods used as basis for this work. Also a new computer library has been developed to manage data tables of thermodynamic properties as well as includes the management of opacity and ionization data tables. This new library extends the capabilities of the previous one with more specific routines, and the possibility of using múltiple tables for several materials. • Diffusion solver: a new package has been developed to solve the diffusion equation applied to the heat conduction of the plasma. The previous method is not parallelizable and it produced mesh dependent results, while this new package can be executed in parallel and achieves a more converged solution that does not depend on the refinement of the mesh. • Radiation package: the check of the radiation model rose several bugs in the implementation that had been removed. The new computer library for EOS managing includes capabilities to store opacity tables for multigroups of energy. Moreover the transport coefficients calculations have been extended for the new multimaterial hydrodynamic package developed by David Portillo García. Also the solving methodology for the transport equation has been modified to make the code run in parallel. • A new ray tracing package has been introduced to extend the previous one to 3D. Once all these tasks has been implemented, the ARWEN code has been applied to study laboratory astrophysics systems. Simulations have been done in order to reproduce the results of the experiments carried out in PALS facility by Chantal Stehlé in radiative shock production. Experimental results are in cióse agreement to the ARWEN estimations of the speed of the shock wave and the length of the precursor. The simulation code used in this thesis, including the work done in ARWEN by other colleagues at the time of this research, allowed the collaboration with other research institution in France and Japan and some of the results presented in this thesis have been published in scientific journals.