Modelación numérica de deslizamientos de ladera en embalses mediante el método de partículas y elementos finitos (PFEM)

Se presentan resultados obtenidos mediante la aplicación del MétododePartículasyElementosFinitos (PFEM) en la simulación dedeslizamientosdeladeraenembalses. Es un fenómeno complejo, por la interacción entre el material deslizado, la masa de agua del embalse, yel conjunto formado por el vaso y la presa. PFEM es un esquema numérico original con el que se ha afrontado con éxito la resolución de problemas de interacción fluido-estructura. Combina un enfoque Lagrangiano con la resolución de las ecuaciones deelementosfinitosmediante la generación de una malla, que se actualiza en cada paso de tiempo. Se presentan resultados de casos de validación en los que se han reproducido ensayos en laboratorio existentes en la bibliografía técnica. Se muestran también otros cálculos más complejos, sobre la cartografía a escala real de un embalse, donde se aprecia el fenómeno de generación de la ola, su propagación por el embalse y la afección a la presa. Por último, se ha modelado el deslizamiento ocurrido en 1958 en la bahía de Lituya (Alaska), enel que la caída de 90 millones de toneladas de roca produjo una ola que alcanzó una sobreelevación máxima de 524 m en la ladera opuesta. Los resultados permiten afirmar que PFEM puede ser una herramienta útil enel análisis de riesgos frente a este tipo de fenómenos, ofreciendo una buena aproximación de las afecciones potenciales

Modelización numérica del comportamiento mecánico de un fieltro de polipropileno

En el presente trabajo se presenta un modelo del continuo para un fieltro denso. Igualando la densidad de potencia de un elemento del continuo a la densidad de potencia mecánica que actúa sobre el conjunto de las fibras se obtiene una expresión del tensor de tensiones en la configuración de referencia. El modelo se completa mediante la inclusión de un modelo de daño para modelar fenomenológicamente los mecanismos de extracción (pull-out) y rotura de las fibras. Se ha implementado el modelo como una subrutina de material de usuario para un código de elementos finitos (ABAQUSExplicit), formulado en grandes deformaciones. Los resultados obtenidos se han comparado con experimentos realizados sobre un fieltro comercial (geotextil) de fibras de polipropileno y muestran que el modelo es capaz de reproducir el comportamiento del material hasta la localización del daño y pérdida de capacidad portante del mismo.

Simulación numérica del cruce de trenes: combinación de un método de reducción del modelo con una implementación de los elementos de contorno

En este artículo se presentan los resultados obtenidos con un método aproximado para calcular las historias temporales de los coeficientes aerodinámicos correspondientes a las fuerzas y los momentos a que están sometidos los coches de cabeza durante el cruce de dos trenes de alta velocidad al aire libre.

Estructuras de hormigón armado bajo carga dinámica severa. Parte III: ejemplo de simulación numérica

Publicamos la tercera parte y última de este artículo sobre estructuras de hormigón armado bajo carga dinámica severa en el que, los autores, nos han ido exponiendo diferentes ejemplos. En esta ocasión, se pone patente que con el desarrollo de los ordenadores durante las últimas décadas se ha dado la posibilidad de utilizar el método de los elementos finitos (FEM: Finite Elements Method) para estudiar los efectos de una carga dinámica severa, como en el caso de las explosiones. A continuación se describe el uso del método para este tipo de situaciones de carga.

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

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

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

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

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

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

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

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

Simulación numérica por el Método de los Elementos de Contorno. Utilización de elementos isoparamétricos parabólicos

Se presenta en esta comunicación el tratamiento de problemas de potencial en sistemas bidimensionales, haciendo uso de la discretización de su contorno o frontera mediante elementos parabólicos tanto en geometría como en las variables de campo. Se estudian las ventajas frente al uso de elementos isoparamétricos lineales dentro de la teoría del potencial. Se presenta también un estudio sobre las zonas singulares a que dan lugar los elementos parabólicos degenerados = This paper presents a B.I.E.M. for potential theory, using in the discretization a completely isoparametric parabolic formulation; that is, the field variable, its first derivative and the boundary domain are interpolated using second orden piecewise polinomic. Several results are presented and comparison is mode with other simpler formulations. Also treated is the posibility of modelling singular behavior by moving the midside mode of selected elements.

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.

Modelación física y numérica de aliviaderos con cajeros altamente convergentes

Presentación de las conclusiones del proyecto de investigación ALCON. La publicación del Reglamento Técnico sobre seguridad de Presas y Embalses en el año 1996 ha generado la necesidad de adaptar muchas de las presas en explotación en España para cumplir esta nueva normativa, debiendo ampliar gran parte de ellas su capacidad de desagüe y también aumentar los resguardos. Muchas de las presas españolas se encuentran aún en proceso de adaptación a los nuevos criterios debido en ocasiones a que la ampliación del aliviadero existente es muy costosa o plantea dificultades técnicas importantes. Por esta razón, en los últimos años se han desarrollado nuevas tipologías de aliviaderos, capaces de satisfacer las necesidades planteadas. Una de estas soluciones innovadoras son los aliviaderos con cajeros altamente convergentes. Esta tipología, adecuada para presas de gravedad de hormigón, tiene la particularidad de que la longitud del vertedero es sensiblemente mayor que el ancho del cuenco amortiguador. Para conducir el caudal vertido al cuenco, se disponen unos cajeros con fuerte convergencia. Si bien este tipo de aliviaderos puede plantearse en presas de nueva construcción, la ampliación de capacidad de presas existentes mediante esta tipología resulta sencilla.Desde el punto de vista de la explotación hidroeléctrica, la ampliación de capacidad conseguida con esta tipología permite aumentar la cota del nivel máximo normal (NMN) gracias a la reducción de la altura de lámina vertiente, manteniendo los resguardos exigidos por la normativa. Ello se traduce en un incremento considerable de la capacidad de generación hidroeléctrica, aumentando tanto la altura del salto como el volumen embalsado. Los aliviaderos con cajeros altamente convergentes ya han sido utilizados en varias presas tanto en España como en el resto del mundo. Sin embargo no existen referencias con criterios de diseño de carácter general para esta tipología, por lo que se hace necesario a la hora de proyectar un aliviadero de este tipo realizar ensayos en modelo físico específicos en cada caso. Con el objetivo desarrollar una metodología para el cálculo hidráulico y proyecto de este tipo de aliviaderos se está llevando a cabo el proyecto ALCON, incluido en el Plan Nacional de Investigación Científica del Ministerio de Economía y Competitividad, en el que participan ALATEC, UPM Y CIMNE. En el proyecto se combinan la modelación física y numérica para alcanzar los objetivos planteados.

Simulación numérica de un problema de contaminación atmosférica

La ecuación en derivadas parciales de advección difusión con reacción química es la base de los modelos de dispersión de contaminantes en la atmósfera, y los diferentes métodos numéricos empleados para su resolución han sido objeto de amplios estudios a lo largo de su desarrollo. En esta Tesis se presenta la implementación de un nuevo método conservativo para la resolución de la parte advectiva de la ecuación en derivadas parciales que modela la dispersión de contaminantes dentro del modelo mesoescalar de transporte químico CHIMERE. Este método está basado en una técnica de volúmenes finitos junto con una interpolación racional. La ventaja de este método es la conservación exacta de la masa transportada debido al empleo de la ley de conservación de masas. Para ello emplea una formulación de flujo basado en el cálculo de la integral ponderada dentro de cada celda definida para la discretización del espacio en el método de volúmenes finitos. Los resultados numéricos obtenidos en las simulaciones realizadas (implementando el modelo conservativo para la advección en el modelo CHIMERE) se han comparado con los datos observados de concentración de contaminantes registrados en la red de estaciones de seguimiento y medición distribuidas por la Península Ibérica. Los datos estadísticos de medición del error, la media normalizada y la media absoluta normalizada del error, presentan valores que están dentro de los rangos propuestos por la EPA para considerar el modelo preciso. Además, se introduce un nuevo método para resolver la parte advectivadifusiva de la ecuación en derivadas parciales que modeliza la dispersión de contaminantes en la atmósfera. Se ha empleado un método de diferencias finitas de alto orden para resolver la parte difusiva de la ecuación de transporte de contaminantes junto con el método racional conservativo para la parte advectiva en una y dos dimensiones. Los resultados obtenidos de la aplicación del método a diferentes situaciones incluyendo casos académicos y reales han sido comparados con la solución analítica de la ecuación de advección-difusión, demostrando que el nuevo método proporciona un resultado preciso para aproximar la solución. Por último, se ha desarrollado un modelo completo que contempla los fenómenos advectivo y difusivo con reacción química, usando los métodos anteriores junto con una técnica de diferenciación regresiva (BDF2). Esta técnica consiste en un método implícito multipaso de diferenciación regresiva de segundo orden, que nos permite resolver los problemas rígidos típicos de la química atmosférica, modelizados a través de sistemas de ecuaciones diferenciales ordinarias. Este método hace uso de la técnica iterativa Gauss- Seidel para obtener la solución de la parte implícita de la fórmula BDF2. El empleo de la técnica de Gauss-Seidel en lugar de otras técnicas comúnmente empleadas, como la iteración por el método de Newton, nos proporciona rapidez de cálculo y bajo consumo de memoria, ideal para obtener modelos operativos para la resolución de la cinética química atmosférica. ABSTRACT Extensive research has been performed to solve the atmospheric chemicaladvection- diffusion equation and different numerical methods have been proposed. This Thesis presents the implementation of an exactly conservative method for the advection equation in the European scale Eulerian chemistry transport model CHIMERE based on a rational interpolation and a finite volume algorithm. The advantage of the method is that the cell-integrated average is predicted via a flux formulation, thus the mass is exactly conserved. Numerical results are compared with a set of observation registered at some monitoring sites in Spain. The mean normalized bias and the mean normalized absolute error present values that are inside the range to consider an accurate model performance. In addition, it has been introduced a new method to solve the advectiondiffusion equation. It is based on a high-order accurate finite difference method to solve de diffusion equation together with a rational interpolation and a finite volume to solve the advection equation in one dimension and two dimensions. Numerical results obtained from solving several problems include academic and real atmospheric problems have been compared with the analytical solution of the advection-diffusion equation, showing that the new method give an efficient algorithm for solving such problems. Finally, a complete model has been developed to solve the atmospheric chemical-advection-diffusion equation, adding the conservative method for the advection equation, the high-order finite difference method for the diffusion equation and a second-order backward differentiation formula (BDF2) to solve the atmospheric chemical kinetics. The BDF2 is an implicit, second order multistep backward differentiation formula used to solve the stiff systems of ordinary differential equations (ODEs) from atmospheric chemistry. The Gauss-Seidel iteration is used for approximately solving the implicitly defined BDF solution, giving a faster tool than the more commonly used iterative modified Newton technique. This method implies low start-up costs and a low memory demand due to the use of Gauss-Seidel iteration.

Integración numérica de la ecuación DNLS

Interacciones no lineales de ondas de Alfven existen tanto para plasmas en el espacio como en laboratorios. En ingeniería aeroespacial amarras electrodinámicas espaciales ("tethers") generan emisión de ondas de Alfven en estructuras denominadas "Alas de Alfven". La ecuación Derivada no lineal de Schrödinger (DNLS) posee la capacidad de describir la propagation de ondas de Alfven de amplitud finita circularmente polarizadas tanto para plasmas fríos como calientes. En esta investigación, dicha ecuacion es solucionada numéricamente por medio de tecnicas espectrales para las derivadas espaciales y un esquema de Runge-Kutta de 4to orden para evaluar el avance en el tiempo. Se considera la ecuacion DNLS sin efectos difusivos, sin embargo se mantienen el termino lineal y lineal y el dispersivos. Se ha trabajado con dos condiciones iniciales: 1 - Una onda, 2 - Tres ondas cerca de resonancia (k1 = k2 + k3), la primera onda (correspondiente a k1) está excitada y las otras dos amortiguadas. En el caso de una única onda los resultados numéricos verifican las condiciones analíticas de estabilidad modular, adems se ha encontrado que el tiempo en el que se produce la inestabilidad y la forma en que evoluciona el sistema depende, para un mismo número de onda, de la amplitud inicial. Con tres ondas se realiza un estudio numérico tanto para plasmas frios como calientes encontrándose que aparece redistribución de energía en un gran número de nodos tanto para ondas polarizadas hacia la izquierda como hacia la derecha.

Avances en investigación aplicada mediante modelación física y numérica en el diseño de la ingeniería de presas

En la actualidad, y en todo el mundo, hay en desarrollo un muy importante número de proyectos de obras hidráulicas de diversa naturaleza (presas, canales, desaladoras, tanques de tormentas, centrales hidroeléctricas, obras de saneamiento, etc.), donde España es, en muchas ocasiones, el marco de referencia. En este artículo, se presentan algunas de las principales investigaciones en curso en el campo de la modelación física y numérica de la ingeniería de presas, con el objetivo de mejorar el conocimiento de los fenómenos hidráulicos que intervienen en su gestión y desarrollar nuevas herramientas de diseño que permitan dar solución a problemas hidráulicos complejos.

Simulación numérica computacional de fenómenos dinámicos

En este proyecto se trata la simulación numérica de un fenómeno dinámico, basado en el comportamiento de una onda transmitida a lo largo de una cuerda elástica de un instrumento musical, cuyos extremos se encuentran anclados. El fenómeno físico, se desarrolla utilizando una ecuación en derivadas parciales hiperbólicas con variables espacial y temporal, acompañada por unas condiciones de contorno tipo Dirichlet en los extremos y por más condiciones iniciales que dan comienzo al proceso. Posteriormente se han generado algoritmos para el método numérico empleado (Diferencias finitas centrales y progresivas) y la programación del problema aproximado con su consistencia, estabilidad y convergencia, obteniéndose unos resultados acordes con la solución analítica del problema matemático. La programación y salida de resultados se ha realizado con Visual Studio 8.0. y la programación de objetos con Visual Basic .Net In this project the topic is the numerical simulation of a dynamic phenomenon, based on the behavior of a transmitted wave along an elastic string of a musical instrument, whose ends are anchored. The physical phenomenon is developed using a hyperbolic partial differential equation with spatial and temporal variables, accompanied by a Dirichlet boundary conditions at the ends and more initial conditions that start the process. Subsequently generated algorithms for the numerical method used (central and forward finite differences) and the programming of the approximate problem with consistency, stability and convergence, yielding results in line with the analytical solution of the mathematical problem. Programming and output results has been made with Visual Studio 8.0. and object programming with Visual Basic. Net