Elementos de contorno adaptables

Como es bien sabido, en el método de los elementos finitos se suele hablar de dos tipos de convergencia. La primera, o convergencia h, se refiere a la mejora del resultado que se obtiene refinando la malla. Debido a la correspondencia elemento-variables nodales-funciones de interpolación, ello implica un ajuste progresivo de los resultados en aquellas zonas donde se produce el refinamiento. Se trata del método más usado cuando, de forma pragmática, se desea tener una idea de la convergencia de los resultados. Su principal inconveniente radica en el hecho que cada refinamiento exige el cálculo de matrices de rigidez diferentes de las anteriores, de modo que la información debe ser rehecha en cada caso y, por tanto, los costes son elevados. El segundo método analiza la convergencia p, o refinamiento de la aproximación mediante el incremento del grado del polinomio definido sobre cada elemento. Se trata de abandonar la idea de asociar a cada nodo el valor físico de la variable correspondiente en la aproximación típica: u ~ a1Ø1 + a2Ø2 + a3Ø3+ … + anØn; donde las funciones Ø son unidad en el nodo correspondiente y cero en el resto. Por el contrario, se vuelve a la idea original de Ritz, semejante al de un desarrollo en la serie de Fourier, donde las funciones Ø están definidas globalmente y los coeficientes de ponderación no tienen por qué presentar un significado físico concreto. Evidentemente la vuelta no es total; se siguen manteniendo elementos y dentro de cada uno de ellos se establece una jerarquía de funciones Øi. Con esta situación intermedia entre la globalidad absoluta de Ritz y la correspondencia absoluta de la discretización con las variables se consigue, por un lado, mantener una versatilidad suficiente para el ajuste por trozos y, por otro, refinar la aproximación de forma inteligente ya que, al igual que sucede en una serie de Fourier, cada término que se añade produce un efecto menor, lo que posibilita el truncamiento cuando se alcanza un determinado nivel de precisión. Además, puesto que cada Ø tiene un soporte perfectamente definido desde un principio, cada etapa del refinamiento aprovecha todos los cálculos anteriores y sólo se necesita evaluar los nuevos términos de la matriz de rigidez. La primera idea fue propuesta por Zienckiewicz et al.(1970), y posteriormente han desarrollado el método Szabo et al.(1978), Babuska (1975,1978), Peano (1978)etc. El proceso operativo incluye así: a)Establecimiento de una malla amplia sobre el dominio a analizar; b)Definición de una jerarquía de funciones de interpolación dentro de cada elemento; c)Establecimiento de un "indicador" de las zonas que precisen la adición de nuevas funciones jerarquizadas; d)Establecimiento de un "estimador a posteriori" que evalúe el error cometido y precise el momento en que pueda ser detenido el proceso. Un método que sigue los pasos anteriores se denomina autoadaptable y, como se puede comprender, resulta interesantísimo para problemas no triviales. En este artículo, se contempla la posibilidad de extender las ideas anteriores al método de los elementos de contorno.

Elementos de contorno adaptables

Como es bien sabido, en el método de los elementos finitos se suele hablar de dos tipos de convergencia. La primera, o convergencia h, se refiere a la mejora del resultado que se obtiene refinando la malla. Debido a la correspondencia elemento-variables nodales-funciones de interpolación, ello implica un ajuste progresivo de los resultados en aquellas zonas donde se produce el refinamiento. Se trata del método más usado cuando, de forma pragmática, se desea tener una idea de la convergencia de los resultados. Su principal inconveniente radica en el hecho que cada refinamiento exige el cálculo de matrices de rigidez diferentes de las anteriores, de modo que la información debe ser rehecha en cada caso y, por tanto, los costes son elevados. El segundo método analiza la convergencia p, o refinamiento de la aproximación mediante el incremento del grado del polinomio definido sobre cada elemento. Se trata de abandonar la idea de asociar a cada nodo el valor físico de la variable correspondiente en la aproximación típica: u ~ a1Ø1 + a2Ø2 + a3Ø3+ … + anØn; donde las funciones Ø son unidad en el nodo correspondiente y cero en el resto. Por el contrario, se vuelve a la idea original de Ritz, semejante al de un desarrollo en la serie de Fourier, donde las funciones Ø están definidas globalmente y los coeficientes de ponderación no tienen por qué presentar un significado físico concreto. Evidentemente la vuelta no es total; se siguen manteniendo elementos y dentro de cada uno de ellos se establece una jerarquía de funciones Øi. Con esta situación intermedia entre la globalidad absoluta de Ritz y la correspondencia absoluta de la discretización con las variables se consigue, por un lado, mantener una versatilidad suficiente para el ajuste por trozos y, por otro, refinar la aproximación de forma inteligente ya que, al igual que sucede en una serie de Fourier, cada término que se añade produce un efecto menor, lo que posibilita el truncamiento cuando se alcanza un determinado nivel de precisión. Además, puesto que cada Ø tiene un soporte perfectamente definido desde un principio, cada etapa del refinamiento aprovecha todos los cálculos anteriores y sólo se necesita evaluar los nuevos términos de la matriz de rigidez. La primera idea fue propuesta por Zienckiewicz et al.(1970), y posteriormente han desarrollado el método Szabo et al.(1978), Babuska (1975,1978), Peano (1978)etc. El proceso operativo incluye así: a)Establecimiento de una malla amplia sobre el dominio a analizar; b)Definición de una jerarquía de funciones de interpolación dentro de cada elemento; c)Establecimiento de un "indicador" de las zonas que precisen la adición de nuevas funciones jerarquizadas; d)Establecimiento de un "estimador a posteriori" que evalúe el error cometido y precise el momento en que pueda ser detenido el proceso. Un método que sigue los pasos anteriores se denomina autoadaptable y, como se puede comprender, resulta interesantísimo para problemas no triviales. En este artículo, se contempla la posibilidad de extender las ideas anteriores al método de los elementos de contorno.

Extreme temperature events analyzed with Fast Fourier Transform

Extreme weather and climate events have received increased attention in the last few years, due to the often large loss of agriculture business and exponentially increasing costs associated with them and insurance planning. This increased attention raises the question as to whether extreme weather and climate events are truly increasing, whether this is only a perceived increase exacerbated by enhanced media coverage, or both. There are a number of ways extreme climate events can be defined, such as extreme daily temperatures, extreme daily rainfall amounts, and large areas experiencing unusually warm monthly temperatures, among others. In this study, we will focus our attention in frost and heatstroke events measuring it as the number of days under 0 ºC and number of days with daily maximum over 30ºC monthly respectively. We have studied the trends in these extreme events applying a Fast Fourier Transform to the series to clarify the tendency. Lack of long-term climate data suitable for analysis of extremes is the single biggest obstacle to quantifying whether extreme events have changed over the twentieth century, including high temporal and spatial resolution observations of temperatures. However, several series have been grouped in different ways: chosen the longest series independently, by provinces, by main watersheds and altitude. On the other hand, synthetic series generated by Luna and Balairón (AEMet) were also analyzed. The results obtained by different pooling data are discussed concluding the difficulties to assess the extreme events tendencies and high regional variation in the trends.

Formulación tridimensional del método de los elementos de contorno con interpolación parabólica

El método de los elementos de contorno ha suscitado un interés creciente en los últimos diez años, presentándose como una herramienta útil para la resolución de problemas de la ingeniería estructural , modelados matemáticamente por sistemas de ecuaciones diferenciales en derivadas parciales. A través de una formulación integral en el contorno del dominio, donde se define el problema, se consigue la resolución de éste, mediante la sola discretización de aquél, en contraposición a los métodos de dominio, que necesitan la completa discretización del mismo (Elementos Finitos ). Estos años han servido, por un lado para cimentar el método anterior - en cuanto a sus bases matemáticas se refiere, y por otro para definir los campos de la Física donde su utilización podría ser más efectiva. Puede decirse que han sido lo que la década de los 60 para elementos finitos : la preparación inicial para el desarrollo espectacular hoy alcanzado. De hecho, el método de los elementos finítos es capaz de abordar eficazmente los grandes retos que plantea el cálculo estructural en la actualidad, mientras que los elementos de contorno no han alcanzado aún dicho estadio. Esta Tesis pretende establecer el camino definitivo para la resolución de este tipo de problemas en el caso elástico tridimensional. Para ello, se ha desarrollado la formulación pertinente tanto en su fase, puramente matemática, como numérica para medios tridimensionales heterogéneos, así como un programa de grandes posibilidades que permiten atacar cualquier caso prácticamente sin limitaciones de tamaño, geometría ó condiciones de contorno en el marco de la Elasticidad Tridimensional. El trabajo se compone de dos partes claramente diferenciadas, el estudio y planteamiento del método, así como la descripción de las soluciones adoptadas para la resolución de los múltiples problemas numéricos que plantea no solo el método en sí, sino la implementación de un programa de tales características , y de un anexo que comprende la descripción detallada de cada uno de los apartados que componen el programa. Por último y paralelamente se ha analizado también el caso axisimétrico, problema tridimensional con características muy especificas, desarrollandose su formulación e implementando asimismo un programa que demuestra la exactitud de dicha formulación.

Evaluación de integrales singulares y cuasisingulares en el método de los elementos de contorno bidimensional mediante transformaciones no lineales

Se presentan en este trabajo distintos niétodos, basados en transformaciones de coordenadas no-lineales, para la evaluación de las integrales singulares y cuasisingulares que aparecen en el Método Directo de los Elementos de Contorno. Se detecta un error inherente a algunas de las transformaciones propuestas y finalmente se sugieren dos nuevas transformaciones que mejoran las actualmente disponibles.

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.

Un nuevo procedimiento numérico para la evaluación de integrales singulares de Valor Principal de Cauchy en métodos de contorno

Es indudable que las estrategias para evaluar numéricamente en el computador integrales singulares en el sentido de Valor Principal de Cauchy (VPC) constituyen uno de los aspectos clave en la precisión y confiabilidad del Método de los Elementos de Contorno. Así pues, en este trabajo se presenta un nuevo procedimiento numérico orientado a la evaluación de integrales singulares VPC y de aquellas que contengan otro tipo de singularidades, basado en una transformación bi-cúbica de coordenadas. Se desarrollan las ideas y detalles de la formulación propuesta, obteniéndose expresiones sencillas y atractivas para la cuadratura numérica y rápidamente incorporables a los códigos MEC ya existentes. Se presentan algunos ejemplos numéricos que avalan la estabilidad y precisión de los nuevos algoritmos.

Aplicaciones de la versión p-adaptable del método de contorno en el análisis de medios continuos

Este trabajo está dedicado a la presentación y discusión de los resultados obtenidos al aplicar la versión p-adaptable del Método de los Elementos de Contorno (MEC) a problemas de elastostática bidimensional. Se describen brevemente algunos criterios básicos inherentes al desarrollo e implementación de la versión p-adaptable del MEC. Se presentan y discuten algunos ejemplos ilustrativos con singularidades producidas por cambios repentinos en la geometría y las condiciones de contorno, los cuales demuestran la potencia y veqsatilidad de la técnica p-adaptable propuesta.

Simulación de excavaciones mediante un método de acoplamiento elementos finitos – elementos de contorno

Cuando se modelan sistemas físicos no lineales de extensión infinita, como las excavaciones, se hace necesario simular adecuadamente tanto la solución en el infinito como la no linealidad. El método de elementos finitos es una herramienta efectiva para representar la no linealidad. Sin embargo, el tratamiento del campo infinito truncando el dominio es bastante cuestionable. Por otro lado, el método de elementos de contorno es adecuado para simular el comportamiento en el infinito sin truncamientos. Por combinación de ambos métodos, se puede obtener un uso adecuado de las ventajas de cada uno. En este trabajo se proponen diversas posibilidades de acoplamiento entre los dos métodos. Se desarrollan algoritmos de acoplamiento basados en una descomposición de dominios y se comparan con los esquemas más tradicionales de acoplamiento.

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.

Un nuevo enfoque en elastostática: elementos de contorno p-autoadaptables

Este trabajo describe la aplicación del Método de los Elementos de Contorno en su versión p-adaptable a problemas de elastostática lineal y su implementación en micro-ordenadores. Se demuestra que la técnica aquí propuesta es sumamente ventajosa en el análisis de medios continuos, consiguiéndose velocidades de convergencia mayores que las que se obtienen mediante el refinamiento tipo h. Se incluyen ejemplos ilustrativos.

Aplicación del Método de los Elementos de Contorno al problema termoelástico

El problema termoelástico general admite una formulación lineal y desacoplada en los casos en que las deformaciones y los incrementos de temperatura son pequeños y las variaciones de temperatura son lentas. En general, la resolución por vía analítica del sistema de ecuaciones diferenciales que gobierna el problema elástico lineal presenta grandes dificultades obligando al empleo de métodos aproximados, entre los que se encuentra el método de los elementos de contorno. Desde esta perspectiva se realiza una formulación por la que se aplica el método de los elementos de contorno en el caso termoelástico estacionario y, como aplicación de lo expuesto anteriormente, se desarrolla un programa de ordenador para resolver problemas termoelásticos estacionarios por el método de los elementos de contorno.

Elementos de contorno p-adaptables en problemas tridimensionales

Se presentan en este artículo los resultados obtenidos con la aplicación de la filosofía de p-convergencia a problemas tridimensionales controlados por la ecuación de Laplace. Se pretende con ello sintetizar las ventajas inherentes a la discretización de contorno y a la jerarquización de las funciones de interpolación.

Problemas con singularidades de contorno. Aplicación del Método de los Elementos de Contorno

Esta comunicación tiene como fundamento el demostrar la importancia de las singularidades geométricas para la resolución de sistemas que impliquen condiciones peculiares de contorno. No es necesario a estas alturas, romper ninguna lanza a favor de los métodos numéricos de discretización del contorno para el tratamiento de problemas frente a los métodos de discretización del dominio. Con todo, el simple hecho de ser un procedimiento numérico conlleva ciertas peculiaridades o problemas de difícil soslayo, pero también pueden presentar serias ventajas basadas en los fundamentos del método. En esta comunicación se pretende hacer constancia no solo de la importancia de este método, palpablemente demostrada en las referencias consignadas, sino algunas de las ventajas que puede proporcionar la resolución de problemas que impliquen potencial, como por ejemplo la transmisión de calor en régimen permanente, y que han sido demostradas para problemas de filtraciones por ejemplo, puestos de manifiesto en este artículo.

La extracción de soluciones óptimas en la versión p-adaptable del Método de los Elementos de Contorno

En este articulo se resumen las principales ideas relacionadas con la resolución de problemas elípticos mediante fórmulas de representación. El uso de una familia de funciones interpolantes jerarquizadas permite el establecimiento de un sistema de resolución autoadaptable a un nivel de exactitud prefijado. Se incluye también una comparación descriptiva con el método de los elementos finitos. La extracción de una mejor solución sin refinar la malla se obtiene en el Método de los Elementos de Contorno, gracias a la aplicación de la fórmula de representación para puntos de contorno. Ello permite diseñar una estrategia de indicadores y estimadores que mejora la eficacia de intentos anteriores y permite controlar el desarrollo de la solución tanto en las versiones p como en la h o mixtos.