Análisis experimental y computacional de modelos termofluidomecánicos para el diseño de ingeniería con metales líquidos

Es importante disponer de una herramienta con la cual diseñar dispositivos de uso industrial y comercial que trabajen con metales líquidos (fuentes de neutrones de alta intensidad, núcleos de sistemas de transmutación nuclear, reactores de fisión de nueva generación, instalaciones de irradiación de materiales o reactores de fusión nuclear). Los códigos CFD (Computational Fluid Dynamics) son una de esas herramientas, y la manera de llevar a cabo su validación es la simulación de experimentos existentes. La turbulencia y la presencia de dos o más fases, son los dos principales problemas a los que tiene que hacer frente un código CFD. La mayoría de los modelos de turbulencia presentes en los códigos CFD se basan en considerar la proporcionalidad directa entre el transporte de cantidad de movimiento turbulento y el transporte turbulento de calor. Precisamente, el coeficiente de difusión del calor turbulento, se asume que sea proporcional a la viscosidad turbulenta a través de una constante empírica, llamada número de Prandtl turbulento. El valor de este número, en los códigos comerciales está entre 0,9 y 0,85 dependiendo del modelo de turbulencia, lo cual significa que en los códigos se asume que el transporte turbulento tanto de cantidad de movimiento como de calor, son prácticamente equivalentes. Esta asunción no es cierta en los flujos de metales líquidos, donde se demuestra que la transmisión de calor por turbulencia es pequeña frente a la transmisión de calor molecular. La solución pasa por aumentar el número de Prandtl turbulento, o abandonar la analogía de Reynolds, en el tratamiento de la turbulencia. Por otro lado, en los metales líquidos la capa límite térmica es más ancha que la de velocidad, y las funciones de pared incluidas en los códigos no satisfacen adecuadamente los flujos turbulentos de los fluidos con bajo número de Prantdl (los metales líquidos). Sí serían adecuados, si el mallado es tal, que la celda más cercana a la pared, está dentro de la subcapa laminar, en la cual la propiedad dominante es la conductividad molecular. En la simulación de flujo multifase los códigos se encuentran con una serie de dificultades, que en el caso de que las densidades de los fluidos que intervienen sean muy diferentes entre sí (como ocurre con los metales líquidos y los gases), serán aún mayores. La modelización de la interfase gas metal líquido, así como el encontrar una correlación válida para los coeficientes de resistencia y sustentación para el movimiento de las burbujas en el seno del metal líquido, son dos de los principales retos en la simulación de este tipo de flujos. Las dificultades no se limitan sólo a la simulación mediante CFD, las medidas experimentales de velocidad de las burbujas y del metal líquido también son complicadas. Hay parámetros que no se pueden definir bien: la trayectoria y la forma de las burbujas entre ellos. En el campo de aplicación industrial de los metales líquidos, los altos valores de los coeficientes de expansión volumétrica y de conductividad térmica hacen que estos fluidos sean muy atractivos en la refrigeración por convección libre en dispositivos de alta densidad de potencia. Tomando como base uno de los diseños de ADS (Accelerator Driven System), y teniendo en cuenta la dificultad que conlleva el uso de múltiples modelos físicos, los cálculos realizados muestran cómo, en caso de fallo eléctrico, la operación de la instalación puede continuar de forma segura. Para la validación de los códigos CFD en su uso como herramienta de diseño, uno de los fenómenos donde cuantitativamente más dificultades encuentran los códigos es en los que aparecen en la modelización de las superficies libres. Un buen ajuste de los modelos multifase y de turbulencia es imprescindible en este tipo de simulaciones. Efectivamente, en la instalación de irradiación de materiales IFMIF, la formación de ondas en la superficie libre del flujo de Litio, es un fenómeno que hay que tratar de evitar, y además se requiere predecir las temperaturas, para ver si hay peligro de ebullición del metal líquido. La simulación llevada a cabo se enfoca al análisis termohidráulico. Variando la velocidad de inyección de Litio desde 10 hasta 20 m/s, se comprueba que las temperaturas máximas quedan alejadas del punto de ebullición del Litio, debido al aumento de presión producido por la fuerza centrífuga. Una de las cuestiones más críticas que se presentan en las fuentes de neutrones sería la refrigeración de la ventana metálica sobre la que incide el haz de protones. La simulación de experimentos como MEGAPIE y TS-1, permite la “visualización” de recirculación en el flujo, de los puntos de estancamiento, de los puntos calientes, etc, y da una fotografía de las zonas críticas del diseño.

Desarrollo de un procedimiento para el cálculo de la fuerza de fricción en un contacto EHD. Validación experimental del procedimiento

Con objeto de minimizar la fricción y desgaste, la industria diseña sus máquinas para que estas funcionen en régimen elastohidrodinámico durante el mayor tiempo posible. Científicos de todos los tiempos han dedicado sus vidas al estudio de la lubricación en régimen de película gruesa. Con especial interés, Ingenieros del último siglo y medio han realizado multitud de experimentos y, con los medios disponibles en su momento, han analizado y tabulado resultados, con objeto de disponer de herramientas que permitieran calcular la fuerza de fricción en contactos EHD. En el presente trabajo de tesis se ha elaborado un modelo propio de lubricación en régimen elastohidrodinámico. Este modelo permite la resolución del problema EHD por dos vías. En primer lugar la resolución analítica del problema EHD permite obtener de forma rápida valores aproximados de coeficiente de fricción y espesor de película, comprobar las hipótesis de cada modelo, y estudiar la influencia de los parámetros. En segundo lugar es posible resolver el problema EHD de forma numérica, realizándose un cálculo totalmente paramétrico de coeficiente de fricción y de la película lubricante. El cálculo se realiza de forma interactiva y a través de una herramienta informática. El objetivo del trabajo es disponer de una herramienta que aplicada a cualquier contacto tribológico (engranajes, rodamientos, etc), permita estudiar su eficiencia, predecir flash de temperatura, posibilidad de fallo por desgaste o gripado. Para el ajuste del modelo ha sido necesaria la realización de unos experimentos previos de ajuste de los equipos de ensayo. Los equipos de ensayo de laboratorio nos ofrecen unos resultados de coeficiente de fricción. Estos valores son introducidos en los modelos de la teoría elastohidrodinámica clásica disponibles en la literatura. Los modelos nos permiten estimar valores teóricos de coeficiente de fricción y espesor de película, que son contrastados con los valores de coeficiente de fricción obtenidos experimentalmente. Por repetición de experiencias ha sido posible ajustar los parámetros de material y lubricante que definen los modelos, hasta ajustar los resultados calculados con los obtenidos por vía experimental. Este ajuste ha permitido comprobar los cálculos y validar las hipótesis de los regímenes y condiciones de operación de los equipos de ensayo. Finalmente se han empleado los parámetros de material y lubricante obtenidos, para formar el modelo EHD propio. Introduciendo en el modelo los resultados obtenidos, es posible calcular el espesor de película y el coeficiente de fricción, para cualquier condición de velocidad, temperatura y carga. Para la elaboración de este modelo se ha partido de las ecuaciones de estado de la hidrodinámica (Ley de Continuidad de Reynols). Este modelo predice el comportamiento ideal en régimen hidrodinámico. Para aproximarnos al régimen elastohidrodinámico no es posible considerar a los cuerpos en contacto como sólidos rígidos que mantienen su geometría inalterada, y por tanto se introducen en este modelo las leyes que consideran la deformación elástica de los materiales. Por último es necesario tener en cuenta que sólo cuando el fluido lubricante está sometido a bajas velocidades de deformación y las condiciones de carga no son elevadas, es posible considerar su comportamiento ideal newtoniano, según el cual la tensión tangencial que ofrece es proporcional al gradiente de velocidad al que se encuentra sometido. Para el ajuste del modelo y adecuada predicción de la película lubricante es necesario introducir modelos de comportamiento del lubricante más complejos que tienen en cuenta las propiedades viscoelásticas de este, tales como Newton/Barus, Tensión Cortante Límite y Carreau. Para ello se han tomado medidas reológicas de los fluidos a distintas presiones y temperaturas. Se ha demostrado que los modelos isotermos no predicen adecuadamente el coeficiente de fricción en condiciones de funcionamiento con deslizamientos elevados, debido al calentamiento del aceite por compresión y cizalla. Es estos casos es necesario emplear un modelo termoelastohidrodinámico.

Neutronics design and radiological aspects of HiPER facilities

En este trabajo se han cubierto diferentes asuntos del diseño neutrónico de los aspectos radiológicos de las dos instalaciones del proyecto HiPER. El proyecto HiPER es un proyecto europeo concebido en el marco del programa ESFRI (European Scientific Facilities Research Infrastructure). Está destinado al desarrollo de la energía de fusión nuclear inercial mediante el uso de láseres y el esquema iluminación directa. Consecuentemente, se trata de una instalación con fines exclusivamente civiles. Se divide en dos fases, correspondientes con dos instalaciones: HiPER Engineering y HiPER Reactor. La instalación HiPER Engineering desarrollará las tecnologías implicadas en la ignición de alta repetición de cápsulas de DT por iluminación directa. El HiPER Reactor será una planta demostradora que produzca electricidad haciendo uso de las tecnologías desarrolladas durante la fase HiPER Engineering. El HiPER Engineering se centrará en las tecnologías relevantes para las igniciones a alta repetición de cápsulas de DT usando la iluminación directa. El principal esfuerzo de desarrollo tecnológico se hará en todos los asuntos directamente relacionados con la ignición: láseres, óptica, inyector, y fabricación masiva de cápsulas entre otros. Se espera una producción de entre 5200 MJ/año y 120000 MJ/año dependiendo del éxito de la instalación. Comparado con la energía esperada en NIF, 1200 MJ/año, se trata de un reto y un paso más allá en la protección radiológica. En este trabajo se ha concebido una instalación preliminar. Se ha evaluado desde el punto de vista de la protección radiológica, siendo las personas y la óptica el objeto de protección de este estudio. Se ha establecido una zonificación durante la operación y durante el mantenimiento de la instalación. Además, se ha llevado a cabo una evaluación de la selección de materiales para la cámara de reacción desde el punto de vista de gestión de residuos radiactivos. El acero T91 se ha seleccionado por, siendo un acero comercial, presentar el mismo comportamiento que el acero de baja activación EUROFER97 al evaluarse como residuo con el nivel de irradiación de HiPER Engineering. Teniendo en cuenta los resultados obtenidos para la instalación preliminar y las modificaciones de la instalación motivadas en otros campos, se ha propuesto una instalación avanzada también en este trabajo. Un análisis más profundo de los aspectos radiológicos, así como una evaluación completa de la gestión de todos los residuos radiactivos generados en la instalación se ha llevado a cabo. La protección radiológica se ha incrementado respecto de la instalación preliminar, y todos los residuos pueden gestionarse en un plazo de 30 sin recurrir al enterramiento de residuos. El HiPER Reactor sera una planta demostradora que produzca electricidad basada en las tecnologías de ignición desarrolladas durante la fase HiPER Engineering. El esfuerzo de desarrollo tecnológico se llevará a cabo en los sistemas relacionados con la generación de electricidad en condiciones económicas: manto reproductor de tritio, ciclos de potencia, vida y mantenimiento de componentes, o sistemas de recuperación de tritio entre otros. En este trabajo la principal contribución a HiPER Reactor está relacionada con el diseño de la cámara de reacción y sus extensiones en la planta. La cámara de reacción es la isla nuclear más importante de la planta, donde la mayoría de las reacciones nucleares tienen lugar. Alberga la primera pared, el manto reproductor de tritio y la vasija de vacío. Todo el trabajo realizado aquí ha pivotado en torno al manto reproductor de tritio y sus interacciones con el resto de componentes de la planta. Tras una revisión profunda de la bibliografía de los diseños recientes de cámaras de reacción con características similares a HiPER Reactor, se ha propuesto y justificado un esquema tecnológico innovador para el manto reproductor de tritio. El material fértil selecconado es el eutéctico 15.7 at.% Litio – 84.3 at.% Plomo, LiPb, evitando el uso de berilio como multiplicador neutrónico mientras se garantiza el ajuste online de la tasa de reproducción de tritio mediante el ajuste en el enriquecimiento en 6Li. Aunque se podría haber elegido Litio purom el LiPb evita problemas relacionados con la reactividad química. El precio a pagar es un reto materializado como inventario radiactivo de Z alto en el lazo de LiPb que debe controlarse. El material estructural seleccionado es el acero de baja activación EUROFER97, que estará en contacto directo con le LiPb fluyendo a alta velocidad. En este esquema tecnológico, el LiPb asegurará la autosuficiente de tritio de la planta mientras el propio LiPb extrae del manto el calor sobre él depositado por los neutrones. Este esquema recibe el nombre de manto de Litio-Plomo auto-refrigerado (SCLL por sus siglas en inglés). Respecto de los conceptos SCLL previos, es destacable que nos e requieren componentes del SiC, puesto que no hay campos magnéticos en la cámara de reacción. Consecuentemente, el manto SCLL propuesto para HiPER presenta riesgo tecnológicos moderados, similares a otros dispositivos de fusión magnética, como el HCLL, e incluso inferiores a los del DCLL, puesto que no se require SiC. Los retos que se deben afrontar son el control del inventario de Z alto así como las tasas de corrosión derivadas de la interacción del LiPb con el EUROFE97. En este trabajo se abordan ambos aspectos, y se presentan los respectivos análisis, junto con otros aspectos neutrónicos y de activación, tales como la protección de la vasija de vacío por parte del material fértil para garantizar la resoldabilidad de por vida en la cara externa de la vasija. También se propone y se estudio un ciclo de potencia de Brayton de Helio para dos configuraciones diferentes de refrigeración del sistema primera pared-manto reproductor. Las principales conclusiones de estos estudios son: i) el inventario de Z alto puede controlarse y es comparable al que se encuentra en dispositivos de fusión similares, ii)la vasija de vacío requiere una mayor protección frente a la radiación neutrónica y iii) las tasas de corrosión son demasiado altas y la temperatura media de salida del LiPb es demasiado baja. Tiendo en cuenta estos resultados juntos con otras consideraciones relacionadas con el mantenimiento de componentes y la viabilidad constructiva, se ha propuesto una evolución de la cámara de reacción. Las evoluciones más destacables son la introducción de un reflector neutrónico de grafito, la modificación de la configuración de la óptica final, la forma y el tamaño de la cámara de vacío y una nueva subdivisión modular del manto. Se ha evaluado desde el punto de vista neutrónico, y su análisis y posterior evolución queda fuera del objeto de este trabajo. Los códigos utilizados en este trabajo son: CATIA para la generación de geometrías 3D complejas MCAM para la traducción de archivos de CATIA a formato de input de MCNP MCNP para el transporte de la radiación (neutrones y gammas) y sus respuestas asociadas ACAB para la evolución del inventario isotópico y sus respuestas asociadas MC2ACAB para acoplar MCNP y ACAB para el cómputo de dosis en parada usando la metodología R2S basada en celda. Moritz para visualizar los reultados de MCNP FLUENT para llevar a cabo cálculos de fluido-dinámica Para llevar a cabo este trabajo, han sido necesarias unas destrezas computacionales. Las más relevantes utilizadas son: generación de geometrás 3D complejas y transmisión a MCNP, diferentes tñecnica de reducción de varianza como importancia por celdas y weight windows basado en malla, metodología Rigorous-two-Steps basada en celdas para el cálculo de dosis en parada y la modificación del código ACAB para el cálculos con múltiples espectros en la misma simulación. Como resumen, la contribución de este trabajo al proyecto HiPER son dos diseños conceptuales de instalación: una para HiPER Engineering y otra para HiPER Reactor. La primera se ha estudio en profundidad desde el punto de vista de protección radiológica y gestión de residuos, mientras que la segunda se ha estudiado desde el punto de vista de operación: seguridad, comportamiento, vida y mantenimiento de componentes y eficiencia del ciclo de potencia.

Estimación del error de discretización con principios variacionales multicampo: II. Flexión de placas

Desarrollos recientes para encajar dentro de un marco variacional la llamada Formulación Libre sugieren la posibilidad de introducir un nuevo tipo de estimador de error para cálculos por elementos finitos. Este estimador se basa en una diferencia entre ciertos funcionales multicampo, que toman el mismo valor para la solución exacta del problema. En el presente artículo, dividido en dos partes, se introduce la formulación del estimador para problemas de elasticidad y de flexión de placas según las hipótesis clásicas de Kirchhoff. Se presentan también algunos ejemplos para dar idea de los comportamientos numéricos observados.

Algunos aspectos de la plantación mecánica

El autor analiza en este interesante artículo las ventajas de la adopción de los modernos sistemas de plantación mecánica en sustitución del método tradicional de la plantación a mano y los mayores rendimientos que naturalmente se obtienen.

Recolección mecánica del cacahuete

El cacahuete es una de las plantas hortícolas de mayor riqueza alimenticia y que tiene grandes posibilidades de expansión en nuestro país. Su contenido en grasa varía del 47 al 50 %; el de hidratos de carbono, alrededor del 18 %, y el de proteínas, alrededor del 26 %. El resto hasta 100 (5-8 %) es el tanto por ciento de humedad del cacahuete sin tostar. Su valor energético es considerable, produciendo 1 kg de cacahuetes alrededor de las 6.300 kilocalorías. La producción mundial de cacahuetes es de unos 17 millones de toneladas, siendo los principales productores: India (6 millones de t), China (2,5 millones), Nigeria (1,4 millones), EE.UU. (1,1 millones) y Senegal (1 millón). España encabeza la lista de los productores europeos, con cerca de 15.000 t, producidas en su mayor parte en la región valenciana, habiéndose sembrado cantidades mínimas en Extremadura y Andalucía [1] *. En la región de Valencia ocupa el cacahuete un puesto más entre los cultivos hortícolas de verano; se siembra en mayo-junio y se recoge en octubre-noviembre. Es una especie que fructifica de forma característica.- enterrando el ovario fecundado en la zona superior de la parte subterránea de la planta y formando ahí las vainas (fig. 1). Por tanto, en cierto modo ha de asemejarte su cultivo y recolección a los de las patatas y otros productos subterráneos. El sistema de cultivo es ancestral, tal como se realizaba hace doscientos años, cuando fue introducido en la región valenciana; así también la recolección, que consta de las siguientes operaciones manuales: 1 2.° Arranque: Tirando hacia arriba de la planta firme, pero suavemente, para disminuir al máximo las pérdidas en el suelo. Sacudida de la tierra y cambio de posición durante el secado en la hilera. 3.° Trilla o sacudida enérgica sobre un cesto para la separación de las vainas (flg. 2). 4.° Cribado y limpieza, variable según las condiciones del producto. Las pérdidas son muy altas (pueden llegar hasta el 50 %), especialmente en el arranque, que, según las condiciones del terreno, se hace muy difícil (en estos casos se ayuda el operario de una azada). El rendimiento de estas labores es muy bajo, sobre todo el del arranque, labor ésta que incluye además gran porcentaje de tiempo muerto, debido a ser altamente fatigosa. Se calculan 182 h TH/ha * * para el conjunto de arranque y trilla manuales.

Estructura histológica de la piel de tomate (Lycopersicon esculentum) en relación con su resistencia mecánica.

La resistencia de la piel es un factor determinante de la resistencia mecánica de los frutos de tomate. Los ensayos de punción se han mostrado muy útiles para el conocimiento de las diferencias entre las variedades en ese aspecto. El estudio histológico de las zonas epidérmica e hipodérmica de este mismo material muestra la importancia de la cutícula y de la forma de las células epidérmicas en relación a la resistencia a la punción.

Ensayos de recolección mecánica de tomates para industria en Extremadura

En los procesos de mecanización de los cultivos hortícolas, la mecanización de la recolección es un punto clave porque generalmente es la operación más difícil de mecanizar y, en cierto modo, condiciona la mecanización de las operaciones anteriores y las técnicas de cultivo que se utilicen (material vegetal, marco de cultivo, fertilización, riegos, etc)

Propiedades físicas de las variedades de tomate para la recolección mecánica

La utilización de máquinas para la recolección y manipulación de los productos agrícolas, aún los más frágiles (y perecederos) es un hecho. Sin embargo, se está lejos de llegar a una perfección en los resultados. Las relaciones entre los elementos que constituyen las máquinas, y los productos a manejar son desconocidas: Por un lado, falta el análisis de los elementos de los sistemas de recolección Por otro lado, se desconocen las causas intrínsecas o biológicas (histológicas, fisiológicas, bioquímicas) de la resistencia de los frutos a las distintas agresiones. Siendo el tomate uno de los frutos (de tipo "perecedero" o frágil) en los que más se ha avanzado en lo que respecta a recolección mecánica, sin haberse resuelto totalmente el problema, se cuenta con una experiencia, y un material que constituye una buena base de partida. La amplitud de la variabilidad de los materiales cultiva dos actualmente, es otro punto a favor del estudio del tomate. Aunque actualmente toda esta experiencia se basa en tomate para industria, la mayor parte de las conclusiones de estos estudios serán de inmediata aplicación en la mejora de variedades de tomate para fresca En éstos, ya existe actualmente la necesidad de técnicas de medida, así como de datos básicos, de las características físico-mecánicas de los frutos. En los estudios que se vienen realizando sobre cuajado de los frutos de tomate en invernadero son también necesarias estas técnicas y datos. El propósito de estos trabajos es conseguir técnicas para la diferenciación de variedades de tomate, respecto a las propiedades físicas que influyen en su resistencia a la manipulación mecánica. Dentro de ello, se intenta encontrar las relaciones causales entre las diversas características que puedan contribuir a un mayor - conocimiento de las propiedades mecánicas de estos frutos. La correlación de los ensayos de laboratorio con los sistemas de determinación de la calidad del trabajo - de recolección llevado a cabo por las máquinas contribuirá al desarrollo final de una óptima solución de mecanización. El estudio de los agentes o causas reales de los daños en los procesos de manipulación se inició en 1955, en que Halsey realizó un estudio preliminar sobre los - daños producidos a los frutos de tomate de mercado, in dicando ya que los daños eran producidos principalmente por impactos. A partir de los años 60 se han realizado numerosos estudios sobre daños producidos a los frutos de tomate en la recolección mecánica, dirigiéndose especialmente a los daños en el transporte a granel (aunque también a los producidos en la máquina cosechadora) (O'Brien - et al. 1965, id 1968, id. 1971, O'Brien 1974, Angle et al. 1974). Las características que confieren resistencia a los frutos en estos procesos suponen un complejo de factor e s , entre los que sobresalen la elasticidad del fruto - entero y la resistencia de la piel. En este proceso del transporte los estudios se basan en las características de vibración de los frutos, que dependen a su vez de los factores indicados (elasticidad, etc.).

Efectos dinámicos debidos al tráfico de ferrocarril sobre la infraestructura de vía y las estructuras

La vía tradicional sobre balasto sigue siendo una selección para las líneas de alta velocidad a pesar de los problemas técnicos y la prestación del funcionamiento. El problema de la vía sobre balasto es el proceso continuo del deterioro de éste debido a las cargas asociadas al tráfico ferroviario. En consecuencia es imprescindible un mantenimiento continuado para mantener un alineamiento adecuado de la vía. Por eso se surge la necesidad de comprender mejor el mecanismo involucrado en el deterioro de la vía y los factores claves que rigen su progresión a lo largo de ciclos de carga con el fin de reducir los costos del mantenimiento de la vía y mejorar el diseño de las nuevas vías. La presente tesis intenta por un lado desarrollar los modelos más adecuados y eficientes del vehículo y de la vía para los cálculos de los efectos dinámicos debido al tráfico de ferrocarril sobre la infraestructura de la vía sobre balasto, y por otro evaluar estos efectos dinámicos sobre el deterioro de la vía sobre balasto a largo plazo, empleando un adecuado modelo de predicción del deterioro de la misma. Se incluye en el trabajo una recopilación del estado del arte en lo referente a la dinámica de la vía, a la modelización del vehículo, de la vía y de la interacción entre ambos. También se hace un repaso al deterioro de la vía y los factores que influyen en su proceso. Para la primera línea de investigación de esta tesis, se han desarrollado los diferentes modelos del vehículo y de la vía y la modelización de la interacción entre ambos para los cálculos dinámicos en dos y tres dimensiones. En la interacción vehículo-vía, se ha empleado la formulación de contacto nodo-superficie para establecer la identificación de las superficies en contacto y el método de los multiplicadores de Lagrange para imponer las restricciones de contacto. El modelo de interacción se ha contrastado con los casos reportados en la literatura. Teniendo en cuenta el contacto no lineal entre rueda-carril y los perfiles de irregularidades distribuidas de la vía, se han evaluado y comparado los efectos dinámicos sobre el sistema vehículo-vía en la interacción de ambos, para distintas velocidades de circulación del vehículo, en los aspectos como la vibración del vehículo, fuerza de contacto, fuerza transmitida en los railpads, la vibración del carril. También se hace un estudio de la influencia de las propiedades de los componentes de la vía en la respuesta dinámica del sistema vehículo-vía. Se ha desarrollado el modelo del asiento de la vía que consiste en la implementación del modelo de acumulación de Bochum y del modelo de hipoplasticidad en la subrutina del usuario \UMAT" del programa ABAQUS. La implementación numérica ha sido comprobado al comparar los resultados de las simulaciones numéricas con los reportados en la literatura. Se ha evaluado la calidad geométrica de la vía sobre balasto de los tramos de estudio con datos reales de la auscultación proporcionados por ADIF (2012). Se ha propuesto una metodología de simulación, empleando el modelo de asiento, para reproducir el deterioro de la geometría de la vía. Se usan los perfiles de la nivelación longitudinal de la auscultación como perfiles de irregularidades iniciales de la vía en las simulaciones numéricas. También se evalúa la influencia de la velocidad de circulación sobre el deterioro de la vía. The traditional ballast track structures are still being used in high speed railways lines with success, however technical problems or performance features have led to ballast track solution in some cases. The considerable maintenance work is needed for ballasted tracks due to the track deterioration. Therefore it is very important to understand the mechanism of track deterioration and to predict the track settlement or track irregularity growth rate in order to reduce track maintenance costs and enable new track structures to be designed. This thesis attempts to develop the most adequate and efficient models for calculation of dynamic track load effects on railways track infrastructure, and to evaluate these dynamic effects on the track settlement, using a track settlement prediction model, which consists of the vehicle/track dynamic model previously selected and a track settlement law. A revision of the state of the knowledge regarding the track dynamics, the modelling of the vehicle, the track and the interaction between them is included. An overview related to the track deterioration and the factors influencing the track settlement is also done. For the first research of this thesis, the different models of vehicle, track and the modelling of the interaction between both have been developed. In the vehicle-track interaction, the node-surface contact formulation to establish the identification of the surfaces in contact and the Lagrange multipliers method to enforce contact constraint are used. The interaction model has been verified by contrast with some benchmarks reported in the literature. Considering the nonlinear contact between wheel-rail and the track irregularities, the dynamic effects on the vehicle-track system have been evaluated and compared, for different speeds of the vehicle, in aspects as vehicle vibration, contact force, force transmitted in railpads, rail vibration. A study of the influence of the properties of the track components on the the dynamic response of the vehicle-track system has been done. The track settlement model is developed that consist of the Bochum accumulation model and the hipoplasticity model in the user subroutine \UMAT" of the program ABAQUS. The numerical implementation has been verified by comparing the numerical results with those reported in the literature. The geometric quality of the ballast track has been evaluated with real data of auscultation provided by ADIF (2012). The simulation methodology has been proposed, using the settlement model for the ballast material, to reproduce the deterioration of the track geometry. The profiles of the longitudinal level of the auscultation is used as initial profiles of the track irregularities in the numerical simulation. The influence of the running speed on the track deterioration is also investigated.

Elementos Finitos no Lineales en Mecánica de Sólidos y Estructuras. Introducción básica y formulación del Método de los Elementos Finitos : Curso

Con el método de elementos finitos se ha desarrollado un procedimiento para el tratamiento de ecuaciones diferenciales y su resolución numérica. En este capítulo se va a realizar la aplicación a problemas de potencial regidos por la ecuación de Laplace. Se trata de una situación relativamente sencilla pero que permite comentar la mayoría de las características del método de elementos finitos en desplazamientos y sirve como introducción a situaciones más complicadas. Además, esta ecuación rige problemas de interés como la filtración, la torsión de cilindros, la distribución de termperaturas en régimen estacionario, etc. En lo que sigue la exposición se limita a problemas planos. Tras un primer ejemplo representativo del tema se realiza la formulación débil del problema de potencial. Acto seguido se introduce el concepto de discretización con toda su parafernalia de elementos, nudos, variables nodales, funciones de interpolación etc. El capítulo finaliza con la mención de algunos temas más especializados que el lector estudioso deberá proseguir en otro tipo de textos.

Nuevo procedimiento para determinar la curva tensión deformación plástica a partir del ensayo de compresión diametral

El objetivo de la presente comunicación es presentar un nuevo procedimiento que permite determinar la curva tensión deformación plástica de un material a partir de un ensayo de compresión diametral. El método combina los resultados experimentales del ensayo de compresión diametral, cálculos numéricos por elementos finitos y un algoritmo iterativo original para ajustar los resultados numéricos a los resultados experimentales. El procedimiento se ha aplicado a una aleación de circonio, sometida a tres temperaturas 20°C, 135"C y 300°C con diferentes niveles de hidruración. Se ha obtenido una curva tensión deformación plástica para cada configuración que ajusta los valores experimentales a los numéricos de forma prácticamente perfecta.

Obtención de la energía de fractura a distintas temperaturas en una aleación de circonio a partir del ensayo de compresión diametral

En este trabajo se presenta un procedimiento original para obtener la energía de fractura a partir de ensayos de compresión diametral realizados en tubos de pared delgada. El método propuesto combina resultados experimentales, cálculos numéricos por elementos finitos, en los que el proceso de rotura se ha modelizado con la teoría de la fisura cohesiva, y técnicas de optimización no lineal. Los parámetros del modelo cohesivo se han ajustado de forma automática de modo que la diferencia entre los resultados numéricos y los registros experimentales sea mínima. La técnica se ha aplicado para determinar la energía de fractura de una aleación comercial de circonio a distintas temperaturas 20°C, 135°C y 300°C, cargada con diferentes niveles de hidrógeno, desde 0 hasta 1200 ppm.

Una metodología para el autoaprendizaje y la autoevaluación en Mecánica Básica

En la Escuela de Minas de la Universidad Politécnica de Madrid (UPM) se ha llevado a cabo el proceso de adaptación a los principios de Bolonia. Esto ha implicado cambios en la metodología de la enseñanza. En este artículo se describe una nueva metodología aplicada a un curso introductorio de Mecánica ubicado en el primer semestre del grado de dos diferentes títulos de ingeniería. Los diferentes resultados se presentan mediante un índice que permite evaluar el resultado global del proceso de enseñanza-aprendizaje alcanzado con esta metodología

Capacidades y potencialidades endógenas, de los centros y laboratorios de investigación y desarrollo tecnológico, en las áreas de materiales y mecánica de Ciudad Guayana, Venezuela

El Estado Bolívar cuenta con una superficie de 238.000 Kms2 se encuentra ubicado al sur-este de Venezuela y su capital es Ciudad Bolívar. La región se ha caracterizado a nivel mundial por poseer una gran riqueza que se manifiesta en las reservas de los principales recursos minerales tales como oro, diamante, manganeso, hierro y bauxita, estos dos últimos han permitido la instalación de empresas básicas para el desarrollo de la industria siderúrgica y del aluminio, ubicados en Ciudad Guayana, Además cuenta con el parque industrial metalmecánico más grande del país concentrando un gran numero de Pequeñas y Medianas Industrias (PYMIS). En tal sentido el sector productivo de Guayana no se ha desarrollado integralmente, debido a que no ha orientando sus esfuerzos, en desarrollar e innovar en nuevos productos o mejoras de procesos debido a la poca capacidad y baja inversión en Investigación y Desarrollo Tecnológico (I&DT), para transformar nuestras materias primas y generar productos competitivo con un alto valor agregado. El presente proyecto se planteo como objetivo principal Diseñar Estrategias para Desarrollar las Capacidades y Potencialidades Endógenas, de los Centros y Laboratorios de Investigación y Desarrollo, para la Vinculación y Fortalecimiento de las PYMIS de Manufacturas en las Áreas de Materiales y Mecánica.