Sobre la versión p-adaptable del método de los elementos de contorno

En este artículo 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.

Aplicación del Método de los Elementos de Contorno a problemas dinámicos en medios incomprensibles

Los métodos numéricos, como es sabido, son en muchas ocasiones altamente sensibles a las situaciones singulares que se presentan en diversos problemas de ingeniería, pudiendo dar lugar a errores, locales o globales, que invalidan los resultados proporcionados por un modelo que, en principio, se supone bien planteado. Tal es el caso de la aplicación del Método de los Elementos de Contorno a problemas dinámicos en medios incompresibles, donde el Coeficiente de Poisson se aproxima a su valor límite de 0,5. En estas circunstancias, según el modelo clásico de propagación de ondas, la velocidad de propagación de las ondas longitudinales (ondas P) tiende a infinito, lo que hace inoperantes las expresiones explícitas de la solución fundamental (o función de ponderación) tal como aparecen en la literatura técnica. El recurso de aproximar el valor de dicho Coeficiente como 0,499 ... conduce en ocasiones a resultados insuficientemente exactos cuando no manifiestamente erróneos, especialmente cuando se trabaja con simple precisión, como consecuencia de los errores de truncamiento en expresiones que, en el límite, serían de la forma: (cero)x(infinito). En el presente artículo se muestran los resultados analíticos de las diversas formas indeterminadas que aparecen en la aplicación del método para el caso aludido. Los resultados son de aplicación a diversos problemas de fluidos, así como a problemas de interacción suelo-estructura en terrenos arcillosos. Se presenta un ejemplo de aplicación que muestra la convergencia de la solución obtenida por el procedimiento propuesto, en contraste con la no-convergencia a que conduce la aproximación del Coeficiente de Poisson como 0,499…

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 jerarquizados

Se presenta en este trabajo un nuevo enfoque del M.E.C., según el cual el refinamiento de la solución se logra mediante incremento del número de funciones de forma, en contraposición al tradicional de refinamiento de la discretización. De esta forma se consigue una convergencia más rápida, y de forma automática. Asimismo la discretización del contorno se realiza de una forma natural, viniendo los puntos que la definen impuestos por cambios de geometría o condiciones de contorno solamente. Como conclusión, cabe señalar que la idea de jerarquías de funciones de ínterpolación permite discretizar el contorno en elementos naturales de acuerdo con las discontinuidades previsibles en las variables del problema. Se puede mejorar de una forma independiente y automática la aproximación tanto en diferentes elementos como en las variables esenciales y naturales del problema. Se puede reducir el cálculo de integrales (la parte mas costosa del M.E.C.) al mínimo imprescindible. Se produce un mejor condicionamiento de la matriz y se debilitan los acoplamientos sucesivos. Toda la información precisa es aprovechada al máximo al pasar de una aproximación a la siguiente, lo que facilita el uso de soluciones iterativas.

La versión p-adaptable del Método de los Elementos Finitos y del Método de los Elementos de Contorno: aplicaciones en elasticidad bidimensional

Este trabajo discute la bondad de las técnicas numéricas p-autoadaptables, comparando el Método de los Elementos Finitos (MEF) y el Método de los Elementos de Contorno(MEC). Se presenta un breve resumen de las herramientas matemáticas necesarias para gobernar el proceso de refinamiento en ambos métodos. Finalmente, se presenta un ejemplo ilustrativo de relevancia práctica en ingeniería, el cual pone de manifiesto la potencia y versatilidad de las técnicas p-adaptables frente a situaciones reales.

El Método de los Elementos de Contorno en problemas dinámicos

En esta nota se presentan algunas posibilidades de aplicación de las técnicas de contorno a los problemas dinámicos. El desarrollo hace especial hincapié en la situación estacionaria puesto que es el área en la que tenemos más experiencia. Ello no significa ninguna limitación del método ya que el uso de transformaciones integrales es una clara posibilidad de obtener soluciones transitorias. Uno de los ejemplos presentados se refiere al estado estacionario de una laja sometida a tracción cíclica, con y sin fisuras. El siguiente está relacionado con el cálculo de las impedancias del suelo, necesarias para los estudios de interacción terreno-estructura y, finalmente, se presenta el efecto de ondas incidentes sobre cimientos rígidos.

El Método de los Elementos de Contorno en termoelasticidad axisimétrica

El método de elementos de contorno se aplica con éxito a problemas en los que sea decisiva la reducción de la dimensionalidad del dominio de integración. Por ello es interesante su aplicación a cuerpos tridimensionales simétricos respecto a un eje, donde la reducción es doble y la discretización puede limitarse a un contorno monodimensional. En este artículo se presenta el método aplicado al caso de la termoelasticidad lineal. Como es bien sabido en dicho problema aparecen fuerzas de volumen que, sin embargo, pueden reducirse con cierta facilidad al contorno haciendo el método muy fructífero.

Los elementos mixtos en la teoría de elementos de contorno

En este trabajo se expone la formulación del B.I.E.M. (en problemas de potencial)con elementos mixtos, que representan una interpolación lineal en la función de campo y constante en su derivada. El objetivo primordial de dicha formulación, es el soslayar los problemas que se presentan con los planteamientos anteriores, cuando existen puntos singulares. Los ejemplos que se incluyen, resueltos mediante un programa desarrollado en un miniordenador, confirman que el método propuesto consigue mejores resultados, sin ninguna complicación adicional sobre la formulación general, y con un ahorro significativo en cuanto al número de incógnitas y ecuaciones que se han de resolver = In this paper the mixed B.I.E.M. for bidimensional potential problems is presented. The interpolation of the field variable is done by piecewise constant one. The main idea is to swep out the solution the parasitic disturbances introduced near singular points by a finite value although an important byproduct is the possibility of conexion with domain methods, as F.E.M., in which the same kind of interpolation is worked. Results are very good, as shown by the exemples, and also interesting is the reduction in computer time comparatively to the classical B.l.E.M approach.

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.

Método de los Elementos de Contorno en algunos problemas de interacción suelo-estructura

Los problemas del comportamiento sismico de estructuras masivas de gran responsabilidad y edificios de gran altura asi como el clásico problema del cimiento de las maquinas vibrantes, hacen que el estudio de la interacci6n suelo- estructura adquiera una gran actualidad. Estos problemas, que implican formas y propiedades complicadas, suponen siempre la necesidad de utilizar un modelo numerico del medio considerado. Aqui se emplea el metodo de los elementos de contorno que dadas sus caracteristicas resulta una alternativa muy sugestiva para modelar el suelo y que hace posible el estudio de problemas tridimensionales a un precio razonable. Se introduce un tipo de elementos para problenas bidimensionales , que incluye una singularidad de tipo logaritmico en uno de sus extremos.Se muestran distribuciones de tensiones obtenidas con este tipo de elementos.

El método de los elementos de contorno en los problemas de mecánica del suelo

La preparación de estas notas ha llevado, al más veterano de los autores, a rememorar sus primeros tanteos con los métodos numéricos. Tratando de desarrollar su tesis doctoral sobre efectos dinámicos en puentes de ferrocarril, descubrió, en 1968, en la biblioteca del Laboratorio de Transporte (donde el profesor ]iménez Salas era Director) las Actas de la reunión ASTM en las que Quilan y Sung proponían la asimilación del comportamiento dinámico del semiespacio elástico a un sistema con un grado de libertad. Además de incorporar estos resultados a un modelo de puente para tener en cuenta los fenómenos de interacción dinámica terreno-estructura dicho autor entró en contacto con algunos miembros del equipo de investigación del Prof. ]iménez Salas que, por entonces, estaba explorando la posibilidad de aplicación del ordenador y los métodos numéricos necesarios para tratar los problemas más difíciles de Mecánica de los Medios Continuos. De hecho fue ese grupo quien contribuyó a introducir en España el método de los elementos finitos en la ingeniería civil, pero además, y en relación directa con el título de este artículo fue el propio ]iménez Salas quién inició la línea de trabajo de lo que mas tarde se ha llamado Método Indirecto de Elementos de Contorno que luego fue seguida por otros miembros de su grupo. En aquélla época poco podía sospechar el autor precitado que iba a dedicar una parte sustancial de su vida al desarrollo de ese método numérico en su versión Directa y mucho menos que gran parte de la motivación vendría del problema de interacción dinámica terreno-estructura, una de cuyas primeras soluciones había obtenido en la mencionada visita al Laboratorio de Transporte. En efecto los autores trataban en 1975 de encontrar un procedimiento que les permitiera afrontar el estudio de la interacción en túneles sometidos a carga sísmica y tropezaron, al utilizar el método de elementos finitos, con el problema de las reflexiones de ondas en los contornos artificiales creados al truncar la malla de cálculo. Deseando evitar el uso contornos absorbentes y otros recursos similares se exploró la posibilidad de soluciones fundamentales que incorporasen el comportamiento en el infmito y, fruto de ello, fueron los primeros trabajos que introdujeron el Método Directo de los Elementos de Contorno en España en problemas estáticos. La extensión a teoría del potencial, dinámica en el dominio de la frecuencia, plasticidad, etc tuvo lugar inmediatamente siendo en la mayoría de los casos los problemas típicos de mecánica del suelo los que motivaron y justifican el esfuerzo realizado. Un campo apasionante, el de la poroelasticidad ha dado lugar a nuevas contribuciones y también se han escrito libros de diverso calado que describen las posibilidades del método para dar contestación a preguntas de gran importancia técnica. Los autores quieren poner de manifiesto que la redacción de este trabajo, debe considerarse no solo como la muestra de algunos resultados de aplicación a problemas prácticos, sino también como un homenaje y reconocimiento explícito a la labor precursora del Prof. ]iménez Salas y a su espíritu de permanente curiosidad por el conocimiento científico.

Guía de introducción al Método de los Elementos de Contorno

Entre la impresionante floración de procedimientos de cálculo, provocada por la aplicación intensiva del ordenador, el llamado Método de los Elementos de Contorno (Boundary Element Method o Boundary Integral Equation Method) parece afianzarse como una alternativa útil al omnipresente Método de los Elementos Finitos que ya ha sido incorporado, como una herramienta de trabajo más, al cotidiano quehacer de la ingeniería. En España, tras unos intentos precursores que se señalan en el texto, la actividad más acusada en su desarrollo y mejora se ha centrado alrededor del Departamento que dirige uno de los autores. Después de la tesis doctoral de J. Domínguez en 1977 que introdujo en España la técnica del llamado "método directo", se han producido numerosas aportaciones en forma de artículos o tesis de investigación que han permitido alcanzar un nivel de conocimientos notable. En esta obrita se pretende transmitir parte de la experiencia adquirida, siquiera sea a nivel elemental y en un campo limitado de aplicación. La filosofía es semejante a la del pequeño libro de Hinton y Owen "A simple guide to finite elements" (Pineridge Press, 1980) que tanta aceptación ha tenido entre los principiantes. El libro se articula alrededor de un sólo tema, la solución del problema de Laplace, y se limitan los desarrollos matemáticos al mínimo imprescindible para el fácil seguimiento de áquel. Tras unos capítulos iniciales de motivación y centrado se desarrolla la técnica para problemas planos, tridimensionales y axisimétricos, limitando los razonamientos a los elementos más sencillos de variación constante o lineal. Finalmente, se incluye un capítulo descriptivo donde se avizoran temas que pueden provocar un futuro interés del estudioso. Para completar la información se ha añadido un apéndice en el que se recoge un pequeño programa para microordenador, con el objetivo de que se contemple la sencillez de programación para el caso plano. El programa es mejorable en muchos aspectos pero creemos que, con ello, mantiene un nivel de legibilidad adecuado para que el lector ensaye sobre él las modificaciones que se indican en los ejercicios al final del capítulo y justamente la provocación de ese aprendizaje es nuestro objetivo final.

Sobre el establecimiento de estimadores e indicadores para los métodos de contorno p-adaptables en problemas de potencial

El establecimiento de métodos numéricos capaces de garantizar una cota de exactitud previamente establecida y de autoadaptar la discretización inicial hasta conseguirlo es una de las tendencias más atractivas en la actualidad. En artículos previos se ha mostrado cómo se puede conseguir ésto con el método de los elementos de contorno y aquí se desarrollan nuevas ideas, inspiradas en tratamientos paralelos con el método de los elementos finitos, para fijar los criterios de adaptaci6n y estimación.

Nueua relacion, en curioso romance, en que se declaran las copiosas lluvias, y grandes inundaciones, que ha padecido la ciudad de Sevilla, y sus contorno, los barrios, arrabales, y casas : y las copiosas limosnas que han repartido a los pobres el Cabildo de la Santa Metropolitana Iglesia, el señor Arzobispo, el Cabildo de la ciudad ... este año de 1736

Fecha tomada del tit

El papel del ferrocarril de contorno en el paisaje urbano de la zona sur de Madrid

El papel del ferrocarril de contorno en el paisaje urbano de la zona sur de Madrid