Análisis de un accidente de rotura de tubos de GV mediante la utilización del SGI de Zorita

En el año 2008 la Universidad Politécnica de Madrid y la empresa Gas Natural Fenosa firmaron un acuerdo por el que se creaba el Aula José Cabrera en el Departamento de Ingeniería Nuclear de la UPM. Dicho aula cuenta con el simulador gráfico interactivo de la central nuclear José Cabrera, que es un simulador de alcance total de una central nuclear PWR de un lazo. El objetivo de la ponencia es demostrar la gran aplicación didáctica que tiene dicho aula. El simulador es una herramienta de uso interactivo para trabajo individual o en grupo con los alumnos. Dentro de la asignatura de “Fiabilidad y Análisis del Riesgo” del Máster Ciencia y Tecnología Nuclear de la ETSII-UPM, se propuso la realización de un árbol de eventos para un accidente propuesto por los alumnos. El trabajo que se presenta en esta ponencia ha consistido en el análisis del accidente de rotura de tubos en el generador de vapor usando el simulador gráfico interactivo.

Protocolo de recuperación y reentrenamiento de un jugador de fútbol profesional con una rotura muscular de sóleo.

En el fútbol profesional, la exigencia a la que están sometidos los deportistas y su perseverancia para aumentar el rendimiento ha incrementado considerablemente la incidencia de lesiones. En concreto, son especialmente frecuentes las lesiones musculares. Este trabajo tiene como objetivo el estudio y análisis del protocolo aplicado a un jugador de fútbol profesional de la Primera División Española quien sufre una lesión muscular durante la temporada 2011/2012. Desde el momento en que se produce la lesión se inicia un trabajo paralelo de tratamiento, rehabilitación y readaptación del deportista a la competición. En este proceso interviene un grupo de profesionales especializados en diferentes áreas interrelacionadas dentro de las ciencias del deporte que componen un equipo de trabajo multidisciplinar, destacando la figura del preparador físico en su papel de readaptador. Se estudian las técnicas empleadas para la recuperación funcional y reentrenamiento al esfuerzo antes de la reincorporación progresiva a los entrenamientos, analizando con especial atención los métodos y protocolos empleados por el readaptador en todas las fases de la temporada, tanto a nivel preventivo como terapéutico.

Simulación del proceso de rotura por tracción de alambres de acero perlítico con un modelo de fisura cohesiva

La fractura de metales dúctiles como el acero suele explicarse a partir de la hipótesis de nucleación, crecimiento y coalescencia de microhuecos. A partir de esta teoría, se han desarrollado diversos modelos numéricos, entre los que el modelo de Gurson y sus variantes son los más extendidos. Dichos modelos reproducen matemáticamente el fenómeno físico de crecimiento de huecos resultando en un desarrollo progresivo del daño en el interior del material durante un ensayo de tracción. En estos modelos, el daño comienza a desarrollarse en fases muy tempranas del ensayo, incluso anteriores a la carga máxima. Ensayos realizados por los autores parecen indicar, sin embargo, que en el caso de barras de acero eutéctico empleado en la fabricación de alambres de pretensado, el daño originado en el interior del material como consecuencia del crecimiento de microhuecos sólo es apreciable en un estado muy avanzado del ensayo, momentos antes de producirse la rotura. Además, desde hace décadas se conoce que la triaxialidad de tensiones tiene una fuerte influencia sobre la rotura de los materiales. En este trabajo se presenta un modelo de rotura para elementos de acero sometidos a tracción, basado en un comportamiento cohesivo del material y que contempla el valor de la triaxialidad de tensiones, diferente en cada punto de la sección crítica de rotura. The fracture of ductile materials, such as steel, is usually explained with the theory of nucleation, growth and coalescence of microvoids. Based on this theory, many numerical models have been developed, with a special mention to Gurson-type models. These models simulate mathematically the physical growth of microvoids, leading to a progressive development of the internal damage that takes place during a tensile test. In these models, the damage starts to develop in very early stages of the test. Tests carried out by the authors seem to point out that, in the case of eutectoid steel bars used for manufacturing prestressing steel wires, the internal damage that takes place as a result of the growth of microvoids is only noticeable in a very advanced state of the test. In addition to this, it is known that the stress triaxiality has a strong influence over the fracture of ductile materials. This work presents a fracture model for steel specimens in a tensile test, based on a cohesive behaviour and taking into account the effect of stress triaxiality, which is different in every point of the fracture plane.

Simulación de terremotos por rotura de falla usando el método de los elementos de contorno

Se presenta un método para la simulación de terremotos, basado en la modelización conjunta de la fuente del sismo,emplazamiento de la estructura y trayectoria intermedia. Esta aproximación tiene la ventaja, respecto a los métodos clásicos, de hacer innecesarias muchas de las hipótesis obligadas en estos, fundamentalmente en lo que respecta al tipo y características del tren de ondas incidente. Por otra parte, al englobar la fuente en la discretización se facilita la cuantificación del terremoto, según las ideas más actuales (basadas precisamente en las características físicas de la falla). Naturalmente, el método presenta una parte negativa obvia: la complejidad y extensión del sistema a modelizar. Es por ello que el método de los elementos de contorno, al discretizar tan sólo el contorno, parece la solución adecuada. La comunicación se centra en la formulación del método, poniendo de relieve la facilidad con que la aproximación se incluye en un programa usual de ordenador como una simple subrutina más.

Simulación de la rotura de barras de acero perlítico sometidas a tracción mediante un modelo de fisura cohesiva. Validación experimental

Las probetas cilíndricas fabricadas con materiales metálicos de elevada ductilidad, como el aluminio o el cobre, sometidas a tracción suelen presentar una rotura comúnmente denominada rotura en copa y cono, debido a su geometría. Este tipo de rotura se reproduce numéricamente con éxito mediante el modelo de Gurson-Tvergaard- Needleman, cuya formulación matemática se basa en el fenómeno físico de nucleación, crecimiento y coalescencia de microhuecos. A diferencia de dichos materiales, las barras de acero perlítico, material con una ductilidad apreciable, presentan un frente de rotura plano que no puede simularse correctamente con los modelos antes mencionados, apareciendo una región interior de daño que, en principio, también puede atribuirse a un fenómeno de nucleación y crecimiento de microhuecos, mientras que en el exterior aparece una zona cuya micrografía permite asociar su rotura a un mecanismo de clivaje. En trabajos anteriores los autores han presentado un elemento de intercara cohesivo dependiente de la triaxialidad de tensiones que, incorporado a un código de elementos finitos, permite reproducir de forma razonable el daño que se desarrolla en la región interior mencionada. En este trabajo se presentan los resultados de una campaña experimental que permite validar el modelo desarrollado. Para ello, se ensayan probetas de diferentes diámetros y se comparan los resultados con los obtenidos numéricamente, empleando tres bases extensométricas diferentes en cada uno de los diámetros. Los resultados numéricos se ajustan razonablemente bien a los obtenidos experimentalmente.The cylindrical specimens made of high-ductility metallic materials, such as aluminium and copper, usually fail showing a fracture surface commonly known as cup-cone fracture because of its shape. This type of fracture is successfully reproduced using the Gurson-Tvergaard-Needleman model, which is based on the physical process of nucleation, growth and coalescence of microvoids. Unlike these materials, pearlitic steel bars, which are considerably ductile, show a flat fracture surface that cannot be correctly reproduced with the aforementioned models. In this flat fracture surface, a dark region can be observed in the centre of the specimen, which is the result of a process of nucleation and growth of microvoids, while in the rest of the fracture surface a different region can be identified, which a micrographic study reveals to be the result of a process of cleavage. In previous works, the authors presented a triaxiality-dependent cohesive interface element that, implemented in a finite element code, can reproduce in a reasonably accurate manner the damage that takes place in the dark region mentioned before. The results of an experimental campaign designed to validate the model are presented in this paper. For it, different diameter specimens are tested and these results are compared to those obtained with the numerical models, using three different initial lengths for the strain. Numerical results agree reasonably well with those obtained experimentally.

Estudio mediante modelización numérica del efecto del estado tensional inicial del terreno en la carga de rotura de un ensayo trapdoor

Uno de los fenómenos que tiene más relevancia en la estabilidad de excavaciones subterráneas es el denominado efecto arco. Su estudio se puede llevar a cabo mediante el ensayo ?trapdoor?, conocido principalmente a partir de su presentación por Terzaghi en la primera ICSMFE en 1936 [1], aunque ya con alguna referencia en el siglo XIX. En este trabajo se lleva a cabo un estudio del ensayo "trapdoor", mediante modelización numérica, con el objetivo de analizar el efecto que tiene sobre la formación del arco estable el estado tensional inicial del terreno. Para ello se ha realizado un análisis paramétrico con el programa de diferencias finitas FLAC de ITASCA. La modelización numérica posibilita este tipo de estudios al poder fijar y variar cómodamente todos los parámetros que definen el ensayo. Pero, además, permite un exhaustivo seguimiento del proceso de rotura, por lo que se facilita el análisis de los resultados obtenidos. En un ensayo ?trapdoor? la carga última es función de la geometría del problema y de las propiedades del material, particularmente de su dilatancia. No obstante, los resultados de este trabajo muestran que el estado tensional del terreno determina, conjuntamente con el resto de parámetros del material, el proceso de rotura. Así, para situaciones con una relación entre la tensión horizontal y vertical superior a 0.5, se forma un primer arco estable de dimensiones reducidas que hace que se presente un mínimo en la carga de rotura, lo cual no ocurre en caso contrario. En consecuencia, se justifica que el estado tensional inicial del terreno incide claramente en la presión necesaria para estabilizar una excavación subterránea.

Simulación de la rotura de barras de acero perlítico sometidas a tracción mediante un modelo de fisura cohesiva. Validación experimental

Las probetas cilíndricas fabricadas con materiales metálicos de elevada ductilidad, como el aluminio o el cobre, sometidas a tracción suelen presentar una rotura comúnmente denominada rotura en copa y cono, debido a su geometría. Este tipo de rotura se reproduce numéricamente con éxito mediante el modelo de Gurson-Tvergaard- Needleman, cuya formulación matemática se basa en el fenómeno físico de nucleación, crecimiento y coalescencia de microhuecos. A diferencia de dichos materiales, las barras de acero perlítico, material con una ductilidad apreciable, presentan un frente de rotura plano que no puede simularse correctamente con los modelos antes mencionados, apareciendo una región interior de daño que, en principio, también puede atribuirse a un fenómeno de nucleación y crecimiento de microhuecos, mientras que en el exterior aparece una zona cuya micrografía permite asociar su rotura a un mecanismo de clivaje. En trabajos anteriores los autores han presentado un elemento de intercara cohesivo dependiente de la triaxialidad de tensiones que, incorporado a un código de elementos finitos, permite reproducir de forma razonable el daño que se desarrolla en la región interior mencionada. En este trabajo se presentan los resultados de una campaña experimental que permite validar el modelo desarrollado. Para ello, se ensayan probetas de diferentes diámetros y se comparan los resultados con los obtenidos numéricamente, empleando tres bases extensométricas diferentes en cada uno de los diámetros. Los resultados numéricos se ajustan razonablemente bien a los obtenidos experimentalmente.The cylindrical specimens made of high-ductility metallic materials, such as aluminium and copper, usually fail showing a fracture surface commonly known as cup-cone fracture because of its shape. This type of fracture is successfully reproduced using the Gurson-Tvergaard-Needleman model, which is based on the physical process of nucleation, growth and coalescence of microvoids. Unlike these materials, pearlitic steel bars, which are considerably ductile, show a flat fracture surface that cannot be correctly reproduced with the aforementioned models. In this flat fracture surface, a dark region can be observed in the centre of the specimen, which is the result of a process of nucleation and growth of microvoids, while in the rest of the fracture surface a different region can be identified, which a micrographic study reveals to be the result of a process of cleavage. In previous works, the authors presented a triaxiality-dependent cohesive interface element that, implemented in a finite element code, can reproduce in a reasonably accurate manner the damage that takes place in the dark region mentioned before. The results of an experimental campaign designed to validate the model are presented in this paper. For it, different diameter specimens are tested and these results are compared to those obtained with the numerical models, using three different initial lengths for the strain. Numerical results agree reasonably well with those obtained experimentally.

Procedimiento de diseño de protecciones tipo repié para evitar la rotura por deslizamiento en masa de presas de escollera sometidas a percolaciones extremas

Síntesis del procedimiento de diseño formulado en la tesis doctoral de Rafael Morán. La demanda de seguridad en infraestructuras críticas ha aumentado en los últimos años, especialmente en los países con mayor nivel de desarrollo económico y social. En estos países, el control del estado de la seguridad de estas infraestructuras y las medidas para disminuir los riesgos asociados en caso de avería están adquiriendo una gran relevancia, lo que se ha plasmado en unas normativas más exigentes. Las presas están dentro de este tipo de infraestructuras tanto por su interés estratégico como por los daños que pueden llegar a ocasionar en caso de un mal funcionamiento o de su eventual rotura. La tecnología de protecciones de presas surge como respuesta a esta demanda de seguridad y, por este motivo, ha experimentado un gran avance en las últimas décadas. Dentro de este artículo se va a describir la protección tipo repié y el procedimiento de diseño propuesto para evitar el deslizamiento en masa de presas con espaldón de aguas abajo de escollera, cuando se produce una filtración anormalmente elevada debido a un fallo en la impermeabilidad en el conjunto presa-cimiento. El artículo describe brevemente el estado del arte de la tecnología y de sus aplicaciones así como de la innovación conseguida en este campo a partir de las investigaciones desarrolladas en la Universidad Politécnica de Madrid (UPM), en colaboración con el Centro de Estudios Hidrográficos del CEDEX y el Centro Internacional de Métodos Numéricos en la Ingeniería (CIMNE), en el marco del proyecto de investigación del Plan Nacional denominado EDAMS.

Un método de comprobación de secciones en losas de hormigón armado con familias de armaduras en direcciones arbitrarias. Parte I: Lajas y membranas

Las armaduras en estructuras bidimensionales de hormigón (losas y láminas) se suelen disponer en dos direcciones, típicamente ortogonales. Sin embargo, a veces, particularmente en zonas en las que las tensiones principales son elevadas, se disponen más de dos familias de armaduras y si la geometría del contorno de la estructura no es regular o no es rectangular es preciso colocar familias de armaduras formando ángulos oblicuos entre sí. En general, las direcciones de las tensiones principales en un punto de una estructura bidimensional no coinciden con las de las armaduras, lo que implica una incertidumbre acerca del trabajo de éstas. Esta problemática que aparece en el diseño usual de las estructuras de hormigón armado y pretensado, no suele estar recogida en la mayoría de las instrucciones. En particular, el tratamiento que presenta la norma española HE acerca del armado de las estructuras y elementos distintos de los monodimensionales, es decir, de la viga, es muy escaso. Este trabajo, que se ha dividido en dos partes, presenta un tratamiento unificado de comprobación de las armaduras en estructuras bidimensionales. En esta primera parte se recoge su aplicación a estructuras, tipos laja y membrana, sometidas a esfuerzos de extensión, es decir, axiles y rasantes, contenidos en su plano medio en el caso de lajas o en su plano tangente a la superficie media en el punto de comprobación, si se trata de una membrana. Como es usual, los esfuerzos, que se determinan a partir de un cálculo elástico y lineal, se mayoran mediante los pertinentes coeficientes de seguridad para obtener los llamados esfuerzos de cálculo. En este articulo, las armaduras en el punto en el que se comprueba la estructura se disponen con la máxima generalidad, es decir, una o varias familias formando ángulos arbitrarios en planta, y colocadas bien en el plano medio o simétricamente en planos paralelos equidistantes del anterior y separados de las caras superior e inferior de la estructura por los mismos recubrimientos. La segunda parte de este trabajo, que representa una extensión de la metodología al caso general de flexión-extensión, es objeto de una siguiente publicación. La metodología en este trabajo tiene en cuenta las ecuaciones, dadas por la elasticidad, de equilibrio, compatibilidad y constitutivas entre los esfuerzos conocidos y las tensiones y deformaciones en ambos materiales, hormigón y acero. Naturalmente, la ecuación constitutiva del hormigón no considera su resistencia a tracción, y por concreción se utiliza la conocida parábola rectángulo con posibilidad de rama descendente. Para el acero se supone para la relación tensiones-deformaciones un diagrama bilineal, es decir, se tiene en cuenta el posible endurecimiento. El cálculo, que se lleva a cabo mediante un simple programa de computador, permite obtener en pocos segundos las curvas de las tensiones y de las deformaciones en cada una de las familias de barras, así como de las tensiones principales en el hormigón en función del factor de amplificación de los esfuerzos. De esta forma se deduce el nivel de seguridad que se alcanza en un punto de la estructura de hormigón armado.

Comparativa de protocolos de recuperación de rotura de ligamento cruzado anterior de la rodilla en el fútbol profesional y amateur

La lesión dentro del deporte es un proceso complejo con un gran impacto a nivel tanto socio-económico como en la salud y el rendimiento del deportista. Debido a la amplitud de este campo de investigación el presente trabajo se centra en el estudio de una lesión, denominada de gravedad, la rotura del ligamento cruzado anterior de la rodilla. La elección de esta lesión, como mayor exponente de repercusión de una lesión y en específico en el fútbol, deporte de mayor transcendencia en nuestra sociedad, hace este caso un magnífico ejemplo para apreciar la complejidad del proceso de recuperación de un deportista. Esto se hará analizando el proceso de recuperación de dicha lesión desde el momento que se produce hasta la vuelta del deportista a la competición, centrando la atención en la figura del readaptador, como profesional de la actividad física, y las técnicas de recuperación funcional, reentrenamiento al esfuerzo y trabajo preventivo que se utilizan. Para ello se lleva a cabo el análisis de un caso real a nivel profesional, en la Primera División Española, durante la temporada 2011/2012 y además se valora la situación actual del proceso de rehabilitación y readaptación a nivel amateur mediante la comparación con un caso real de categoría juvenil.

Aplicación de la regresión logística en la predicción empírica de fenómenos complejos en obras subterráneas : squeezing y rotura de pilares de carbón

La heterogeneidad del medio geológico introduce en el proyecto de obra subterránea un alto grado de incertidumbre que debe ser debidamente gestionado a fin de reducir los riesgos asociados, que son fundamentalmente de tipo geotécnico. Entre los principales problemas a los que se enfrenta la Mecánica de Rocas moderna en el ámbito de la construcción subterránea, se encuentran la fluencia de roca en túneles (squeezing) y la rotura de pilares de carbón. Es ampliamente conocido que su aparición causa importantes perjuicios en el coste y la seguridad de los proyectos por lo que su estudio, ha estado tradicionalmente vinculado a la predicción de su ocurrencia. Entre las soluciones existentes para la determinación de estos problemas se encuentran las que se basan en métodos analíticos y numéricos. Estas metodologías son capaces de proporcionar un alto nivel de representatividad respecto del comportamiento geotécnico real, sin embargo, su utilización solo es posible cuando se dispone de una suficiente caracterización geotécnica y por tanto de una detallada definición de los parámetros que alimentan los complejos modelos constitutivos y criterios de rotura que los fenómenos estudiados requieren. Como es lógico, este nivel de definición solo es posible cuando se alcanzan etapas avanzadas de proyecto, incluso durante la propia construcción, a fin de calibrar adecuadamente los parámetros introducidos en los modelos, lo que supone una limitación de uso en etapas iniciales, cuando su predicción tiene verdadero sentido. Por su parte, los métodos empíricos permiten proporcionar soluciones a estos complejos problemas de un modo sencillo, con una baja parametrización y, dado su eminente enfoque observacional, de gran fiabilidad cuando se implementan sobre condiciones de contorno similares a las originales. La sencillez y escasez de los parámetros utilizados permiten a estas metodologías ser utilizadas desde las fases preliminares del proyecto, ya que estos constituyen en general, información habitual de fácil y económica adquisición. Este aspecto permite por tanto incorporar la predicción desde el principio del proceso de diseño, anticipando el riesgo en origen. En esta tesis doctoral, se presenta una nueva metodología empírica que sirve para proporcionar predicciones para la ocurrencia de squeezing y el fallo de pilares de carbón basada en una extensa recopilación de información de casos reales de túneles y minas en las que ambos fenómenos fueron evaluados. Esta información, recogida de referencias bibliográficas de prestigio, ha permitido recopilar una de las más extensas bases de datos existentes hasta la fecha relativa a estos fenómenos, lo que supone en sí mismo una importante contribución sobre el estado del arte. Con toda esta información, y con la ayuda de la teoría de clasificadores estadísticos, se ha implementado sobre las bases de datos un clasificador lineal de tipo regresión logística que permite hacer predicciones sobre la ocurrencia de ambos fenómenos en términos de probabilidad, y por tanto ponderar la incertidumbre asociada a la heterogeneidad incorporada por el medio geológico. Este aspecto del desarrollo es el verdadero valor añadido proporcionado por la tesis y la principal ventaja de la solución propuesta respecto de otras metodologías empíricas. Esta capacidad de ponderación probabilística permite al clasificador constituir una solución muy interesante como metodología para la evaluación de riesgo geotécnico y la toma de decisiones. De hecho, y como ejercicio de validación práctica, se ha implementado la solución desarrollada en un modelo coste-beneficio asociado a la optimización del diseño de pilares involucrados en una de mina “virtual” explotada por tajos largos. La capacidad del clasificador para cuantificar la probabilidad de fallo del diseño, junto con una adecuada cuantificación de las consecuencias de ese fallo, ha permitido definir una ley de riesgo que se ha incorporado al balance de costes y beneficios, que es capaz, a partir del redimensionamiento iterativo del sistema de pilares y de la propia configuración de la mina, maximizar el resultado económico del proyecto minero bajo unas condiciones de seguridad aceptables, fijadas de antemano. Geological media variability introduces to the subterranean project a high grade of uncertainty that should be properly managed with the aim to reduce the associated risks, which are mainly geotechnical. Among the major problems facing the modern Rock Mechanics in the field of underground construction are both, the rock squeezing while tunneling and the failure of coal pillars. Given their harmfulness to the cost and safety of the projects, their study has been traditionally linked to the determination of its occurrence. Among the existing solutions for the determination of these problems are those that are based on analytical and numerical methods. Those methodologies allow providing a high level of reliability of the geotechnical behavior, and therefore a detailed definition of the parameters that feed the complex constitutive models and failure criteria that require the studied phenomena. Obviously, this level of definition is only possible when advanced stages of the project are achieved and even during construction in order to properly calibrate the parameters entered in the models, which suppose a limited use in early stages, when the prediction has true sense. Meanwhile, empirical methods provide solutions to these complex problems in a simple way, with low parameterization and, given his observational scope, with highly reliability when implemented on similar conditions to the original context. The simplicity and scarcity of the parameters used allow these methodologies be applied in the early stages of the project, since that information should be commonly easy and cheaply to get. This aspect can therefore incorporate the prediction from the beginning of the design process, anticipating the risk beforehand. This thesis, based on the extensive data collection of case histories of tunnels and underground mines, presents a novel empirical approach used to provide predictions for the occurrence of both, squeezing and coal pillars failures. The information has been collected from prestigious references, providing one of the largest databases to date concerning phenomena, a fact which provides an important contribution to the state of the art. With all this information, and with the aid of the theory of statistical classifiers, it has been implemented on both databases, a type linear logistic regression classifier that allows predictions about the occurrence of these phenomena in terms of probability, and therefore weighting the uncertainty associated with geological variability. This aspect of the development is the real added value provided by the thesis and the main advantage of the proposed solution over other empirical methodologies. This probabilistic weighting capacity, allows being the classifier a very interesting methodology for the evaluation of geotechnical risk and decision making. In fact, in order to provide a practical validation, we have implemented the developed solution within a cost-benefit analysis associated with the optimization of the design of coal pillar systems involved in a "virtual" longwall mine. The ability of the classifier to quantify the probability of failure of the design along with proper quantification of the consequences of that failure, has allowed defining a risk law which is introduced into the cost-benefits model, which is able, from iterative resizing of the pillar system and the configuration of the mine, maximize the economic performance of the mining project under acceptable safety conditions established beforehand.

Modelización de la rotura y fluidificación en geomateriales y geoestructuras

La rotura de las geoestructuras puede ser causada por cambios en las tensiones efectivas debidos bien por cargas externas (terremotos, por ejemplo), bien por variación de las presiones intersticiales (lluvia), o cambios en la geometría (erosión), así como por una disminución de las propiedades resistentes de los materiales (meteorización, ataque químico, etc). El caso particular de los deslizamientos es interesante, existiendo diversas clasificaciones que tienen en cuenta su forma, velocidad de propagación, etc. Dos de estos casos son los deslizamientos propiamente dichos y los flujos. En el primer caso, la deformación se concentra en zonas de pequeño espesor, que se idealiza como una superficie (superficie de rotura). La cinemática de esta rotura se puede considerar como el movimiento relativo de dos masas cuyas deformaciones no son grandes. Este mecanismo está usualmente asociado a materiales sobreconsolidados que presentan reblandecimiento. Los flujos se producen generalmente en materiales de baja densidad y estructura metaestable, con tendencia a compactar, de forma que se generan presiones intersticiales que aumentan el ángulo de rozamiento movilizado, pudiéndose llegar en algunos casos a la licuefacción. Este mecanismo de rotura se conoce como rotura difusa, y no ha sido tan estudiado como el de localización a pesar de se trata frecuentemente de roturas de tipo catastrófico. De hecho, el suelo pasa de un estado sólido a un estado fluidificado, con una gran movilidad. Es importante para el ingeniero predecir tanto el comportamiento de las geoestructuras bajo las cargas de cálculo como las condiciones en las que se producirá la rotura. De esta manera, en algunos casos, se podrán reforzar las zonas más débiles aumentando así su seguridad. En otros casos, no se podrá realizar este refuerzo (grandes deslizamientos como avalanchas, lahares, etc), pero sí se podrán conocer las consecuencias de la rotura: velocidad de propagación, alcance, espesores, etc. La modelización de estos problemas es compleja, ya que aparecen dificultades en los modelos matemáticos, constitutivos o reológicos y numéricos. Dado que en los geomateriales aparece una interacción fuerte entre el esqueleto sólido y los fluidos intersticiales, esto debe ser tenido en cuenta en los modelos matemáticos. En este trabajo se describirán, pues, el desarrollo y aplicación de técnicas avanzadas de modelización; matemática, constitutiva/reológica y numérica. Se dedicará especial atención a los problemas de transición entre suelos y suelos fluidificados, que hoy en día se estudian en una gran mayoría de los casos con modelos diferentes. Así por ejemplo, se emplean modelos constitutivos para el comportamiento previo a la rotura, y reológicos para los materiales fluidificados. En lo que respecta a los modelos matemáticos, existen formulaciones nuevas en velocidades (o desplazamientos), tensiones, y presiones de aire y agua intersticial, de los que se pueden obtener modelos simplificados integrados en profundidad para deslizamientos rápidos. Respecto de los modelos constitutivos, es interesante la teoría de la Plasticidad Generalizada (modelo básico de Pastor-Zienkiewicz y extensiones a suelos no saturados). Se estudiará la extensión de estos modelos elastoplásticos a la viscoplasticidad (Perzyna y consistente), explorando la posibilidad de emplearlos tanto antes como después de la rotura. Finalmente, en lo que a modelos numéricos se refiere, se describirá la implementación de los modelos matemáticos y constitutivos desarrollados en (i) modelos clásicos de elementos finitos, como el GeHoMadrid desarrollado en los grupo investigador M2i al que pertenece el autor de este trabajo, (ii) Métodos de tipo Taylor Galerkin, y (iii) métodos sin malla como el SPH y el Material Point Model. Estos modelos se aplicarán, principalmente a (i) Licuefacción de estructuras cimentadas en el fondo marino (ii) presas de residuos mineros (iii) deslizamientos rápidos de laderas.

Diseño de un Procesador de Membranas Paralelo para los Sistemas P de Transición

La computación con membranas surge como una alternativa a la computación tradicional. Dentro de este campo se sitúan los denominados Sistemas P de Transición que se basan en la existencia de regiones que contienen recursos y reglas que hacen evolucionar a dichos recursos para poder llevar a cada una de las regiones a una nueva situación denominada configuración. La sucesión de las diferentes configuraciones conforman la computación. En este campo, el Grupo de Computación Natural de la Universidad Politécnica de Madrid lleva a cabo numerosas investigaciones al amparo de las cuales se han publicado numerosos artículos y realizado varias tesis doctorales. Las principales vías de investigación han sido, hasta el momento, el estudio del modelo teórico sobre el que se definen los Sistemas P, el estudio de los algoritmos que se utilizan para la aplicación de las reglas de evolución en las regiones, el diseño de nuevas arquitecturas que mejoren las comunicaciones entre las diferentes membranas (regiones) que componen el sistema y la implantación de estos sistemas en dispositivos hardware que pudiesen definir futuras máquinas basadas en este modelo. Dentro de este último campo, es decir, dentro del objetivo de construir finalmente máquinas que puedan llevar a cabo la funcionalidad de la computación con Sistemas P, la presente tesis doctoral se centra en el diseño de dos procesadores paralelos que, aplicando variantes de algoritmos existentes, favorezcan el crecimiento en el nivel de intra-paralelismo a la hora de aplicar las reglas. El diseño y creación de ambos procesadores presentan novedosas aportaciones al entorno de investigación de los Sistemas P de Transición en tanto en cuanto se utilizan conceptos que aunque previamente definidos de manera teórica, no habían sido introducidos en el hardware diseñado para estos sistemas. Así, los dos procesadores mantienen las siguientes características: - Presentan un alto rendimiento en la fase de aplicación de reglas, manteniendo por otro lado una flexibilidad y escalabilidad medias que son dependientes de la tecnología final sobre la que se sinteticen dichos procesadores. - Presentan un alto nivel de intra-paralelismo en las regiones al permitir la aplicación simultánea de reglas. - Tienen carácter universal en tanto en cuanto no depende del carácter de las reglas que componen el Sistema P. - Tienen un comportamiento indeterminista que es inherente a la propia naturaleza de estos sistemas. El primero de los circuitos utiliza el conjunto potencia del conjunto de reglas de aplicación así como el concepto de máxima aplicabilidad para favorecer el intra-paralelismo y el segundo incluye, además, el concepto de dominio de aplicabilidad para determinar el conjunto de reglas que son aplicables en cada momento con los recursos existentes. Ambos procesadores se diseñan y se prueban mediante herramientas de diseño electrónico y se preparan para ser sintetizados sobre FPGAs. ABSTRACT Membrane computing appears as an alternative to traditional computing. P Systems are placed inside this field and they are based upon the existence of regions called “membranes” that contain resources and rules that describe how the resources may vary to take each of these regions to a new situation called "configuration". Successive configurations conform computation. Inside this field, the Natural Computing Group of the Universidad Politécnica of Madrid develops a large number of works and researches that provide a lot of papers and some doctoral theses. Main research lines have been, by the moment, the study of the theoretical model over which Transition P Systems are defined, the study of the algorithms that are used for the evolution rules application in the regions, the design of new architectures that may improve communication among the different membranes (regions) that compose the whole system and the implementation of such systems over hardware devices that may define machines based upon this new model. Within this last research field, this is, within the objective of finally building machines that may accomplish the functionality of computation with P Systems, the present thesis is centered on the design of two parallel processors that, applying several variants of some known algorithms, improve the level of the internal parallelism at the evolution rule application phase. Design and creation of both processors present innovations to the field of Transition P Systems research because they use concepts that, even being known before, were never used for circuits that implement the applying phase of evolution rules. So, both processors present the following characteristics: - They present a very high performance during the application rule phase, keeping, on the other hand, a level of flexibility and scalability that, even known it is not very high, it seems to be acceptable. - They present a very high level of internal parallelism inside the regions, allowing several rule to be applied at the same time. - They present a universal character meaning this that they are not dependent upon the active rules that compose the P System. - They have a non-deterministic behavior that is inherent to this systems nature. The first processor uses the concept of "power set of the application rule set" and the concept of "maximal application" number to improve parallelism, and the second one includes, besides the previous ones, the concept of "applicability domain" to determine the set of rules that may be applied in each moment with the existing resources.. Both processors are designed and tested with the design software by Altera Corporation and they are ready to be synthetized over FPGAs.

La estabilidad del frente de excavación de un túnel en terrenos heterogéneos o con un criterio de rotura no-lineal

El frente de un túnel puede colapsar si la presión aplicada sobre el es inferior a un valor limite denominado presión “critica” o “de colapso”. En este trabajo se desarrolla y presenta un mecanismo de rotura rotacional generado punto a punto para el cálculo de la presión de colapso del frente de túneles excavados en terrenos estratificados o en materiales que siguen un criterio de rotura nolineal. La solución propuesta es una solución de contorno superior en el marco del Análisis Límite y supone una generalización del mecanismo de rotura mas reciente existente en la bibliografía. La presencia de un terreno estratificado o con un criterio de rotura no-lineal implica una variabilidad espacial de las propiedades resistentes. Debido a esto, se generaliza el mecanismo desarrollado por Mollon et al. (2011b) para suelos, de tal forma que se puedan considerar valores locales del ángulo de rozamiento y de la cohesión. Además, la estratificación del terreno permite una rotura parcial del frente, por lo que se implementa esta posibilidad en el mecanismo, siendo la primera solución que emplea un mecanismo de rotura que se ajusta a la estratigrafía del terreno. Por otro lado, la presencia de un material con un criterio de rotura no-lineal exige introducir en el modelo, como variable de estudio, el estado tensional en el frente, el cual se somete al mismo proceso de optimización que las variables geométricas del mecanismo. Se emplea un modelo numérico 3D para validar las predicciones del mecanismo de Análisis Limite, demostrando que proporciona, con un esfuerzo computacional significativamente reducido, buenas predicciones de la presión critica, del tipo de rotura (global o parcial) en terrenos estratificados y de la geometría de fallo. El mecanismo validado se utiliza para realizar diferentes estudios paramétricos sobre la influencia de la estratigrafía en la presión de colapso. Igualmente, se emplea para elaborar cuadros de diseño de la presión de colapso para túneles ejecutados con tuneladora en macizos rocosos de mala calidad y para analizar la influencia en la estabilidad del frente del método constructivo. Asimismo, se lleva a cabo un estudio de fiabilidad de la estabilidad del frente de un túnel excavado en un macizo rocoso altamente fracturado. A partir de el se analiza como afectan las diferentes hipótesis acerca de los tipos de distribución y de las estructuras de correlación a los resultados de fiabilidad. Se investiga también la sensibilidad de los índices de fiabilidad a los cambios en las variables aleatorias, identificando las mas relevantes para el diseño. Por ultimo, se lleva a cabo un estudio experimental mediante un modelo de laboratorio a escala reducida. El modelo representa medio túnel, lo cual permite registrar el movimiento del material mediante una técnica de correlación de imágenes fotográficas. El ensayo se realiza con una arena seca y se controla por deformaciones mediante un pistón que simula el frente. Los resultados obtenidos se comparan con las estimaciones de la solución de Análisis Límite, obteniéndose un ajuste razonable, de acuerdo a la literatura, tanto en la geometría de rotura como en la presión de colapso. A tunnel face may collapse if the applied support pressure is lower than a limit value called the ‘critical’ or ‘collapse’ pressure. In this work, an advanced rotational failure mechanism generated ‘‘point-by-point” is developed to compute the collapse pressure for tunnel faces in layered (or stratified) grounds or in materials that follow a non-linear failure criterion. The proposed solution is an upper bound solution in the framework of limit analysis which extends the most advanced face failure mechanism in the literature. The excavation of the tunnel in a layered ground or in materials with a non-linear failure criterion may lead to a spatial variability of the strength properties. Because of this, the rotational mechanism recently proposed by Mollon et al. (2011b) for Mohr-Coulomb soils is generalized so that it can consider local values of the friction angle and of the cohesion. For layered soils, the mechanism needs to be extended to consider the possibility for partial collapse. The proposed methodology is the first solution with a partial collapse mechanism that can fit to the stratification. Similarly, the use of a nonlinear failure criterion introduces the need to introduce new parameters in the optimization problem to consider the distribution of normal stresses along the failure surface. A 3D numerical model is employed to validate the predictions of the limit analysis mechanism, demonstrating that it provides, with a significantly reduced computational effort, good predictions of critical pressure, of the type of collapse (global or partial) in layered soils, and of its geometry. The mechanism is then employed to conduct parametric studies of the influence of several geometrical and mechanical parameters on face stability of tunnels in layered soils. Similarly, the methodology has been further employed to develop simple design charts that provide the face collapse pressure of tunnels driven by TBM in low quality rock masses and to study the influence of the construction method. Finally, a reliability analysis of the stability of a tunnel face driven in a highly fractured rock mass is performed. The objective is to analyze how different assumptions about distributions types and correlation structures affect the reliability results. In addition, the sensitivity of the reliability index to changes in the random variables is studied, identifying the most relevant variables for engineering design. Finally, an experimental study is carried out using a small-scale laboratory model. The problem is modeled in half, cutting through the tunnel axis vertically, so that displacements of soil particles can be recorded by a digital image correlation technique. The tests were performed with dry sand and displacements are controlled by a piston that supports the soil. The results of the model are compared with the predictions of the Limit Analysis mechanism. A reasonable agreement, according to literature, is obtained between the shapes of the failure surfaces and between the collapse pressures observed in the model tests and computed with the analytical solution.

Reducción de psicofármacos en aguas residuales después de la filtración con membranas

Se analiza la presencia de diferentes psicofármacos en el agua residual y los rendimientos de eliminación obtenidos en una depuradora convencional. Los psicofármacos prácticamente no se eliminan en el tratamiento biológico. Sin embargo, se ha obtenido una biodegradación considerable debido al tipo de compuesto y de la edad del fango. Además se estudian las retenciones de estos compuestos empleando tratamientos terciarios convencionales y plantas piloto de ultrafiltración y nanofiltración. Las eliminaciones obtenidas no difieren demasiado de un tipo de membrana a otro a pesar del diferente tamaño de paso entre ellas