Implementación del método simplet en un programa de elementos finitos basados en volumenes de control

En este trabajo se presenta la implantación del método SIMPLET en un programa de elementos finitos basados en volúmenes de control. Este método toma en cuanta las variaciones de temperatura en la determinación del campo de presiones en problemas de convección libre, con el próposito de acelerar la convergencia y disminuir los tiempos de computación. El método SIMPLET, que originalmente ha sido propuesto en base al método de diferencias finitas con mallas desplazadas, se implementó en un programa de elementos finitos basados en volúmenes de control con mallas no desplazadas. Se resolvió un problema de convección libre en una cavidad cuadrada y los resultados obtenidos en términos del número de iteracciones y tiempo de computación se compararon se compararon con los resultados del método SIMPLE. Los resultados muestran que en este tipo de problemas el método SIMPLET es más rápido que le SIMPLE cuando el número de Rayleigh es bajo, mientras que para números de Rayleigh elevados, el desempeño de ambos métodos es similar.

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.

Interpolación de Birkhoff a partir de una formulación variacional asociada al método de los elementos finitos

El problema de la interpolación polinómica de Birkhoff, como es conocido, no siempre admite solución. Se presenta un método como alternativa a este tipo de problemas para la obtención de interpolantes con unas condiciones de continuidad determinadas y en el cual el criterio de aproximación es el de minimización de un cierto funcional real, utilizando el método de los elementos finitos. Se analiza la eficiencia del método en diversos ejemplos de aplicación 1-0 y 2-D.

Un estudio sobre la exactitud del método de los elementos finitos : aplicación a la barra recta de sección variable bajo esfuerzos axiles

El objetivo de este trabajo se centra en la formulación clásica del método de los elementos finitos y dentro de ella, únicamente en un aspecto todavía no bien estudiado: su exactitud. En opinión del autor, la complejidad de este estudio, en el estado actual de conocimientos, motiva que aquel no sea susceptible de un planteamiento general. Por esta razón, ha parecido conveniente la consideración de un caso estructural muy simple -la columna de sección variable-, correspondiente a un problema en elementos finitos - de clase C0, que, sin embargo, permite, por un lado, su extensión a situaciones estructurales análogas (problemas de torsión, membranas de revolución, bandas finitas, etc.), sin apenas modificación conceptual, y, por otra parte, aportar indicaciones sobre las posibilidades de llevar a cabo un análisis paralelo en elementos y estructuras más complejos, bien en mayor dimensión (estructuras 2-D y 3-D), bien en requerimientos de continuidad más elevada (estructuras de flexión, tipo c1 , etc.).

Una aplicación del método de los elementos finitos al trazado de carreteras

El método de los elementos finitos ha encontrado numerosas aplicaciones en la resolución de problemas de tipo estructural. No obstante sus posibilidades son mucho más amplias, como se pone de manifiesto en la aplicación que en el siguiente artículo se hace de dicho método al trazado de carreteras, de interés en casos específicos, como el tratamiento de problemas de optimización.

Modelización de la dispersión de contaminantes en la bahía de Santander, mediante el método de los elementos finitos

Existen diferentes métodos para el estudio de la calidad del agua en los estuarios, ríos, lagos, en zonas afectadas por los vertidos de sustancias contaminantes; que van desde los métodos experimentales tradicionales, hasta los más recientes modelos matemáticos. El modelo desarrollado permite resolver la ecuación de dispersión 2—D mediante técnicas en elementos finitos y por lo tanto obtener la evolución espacio-temporal de la concentración de constituyente. Se ha observado su exactitud en diferentes situaciones prácticas, y en particular se ha aplicado a la Bahía de Santander, analizándose el efecto de diferentes tipos de vertidos. El modelo es susceptible de utilizarse conjuntamente con otro modelo hidrodinámico capaz de simular la evolución del campo de velocidades.

Modelo numérico basado en el método de los elementos finitos para el estudio de la hidrodinámica de bahías y estuarios

Se ha desarrollado un modelo implícito no lineal 2-D en EF para la resolución de las ecuaciones de aguas poco profundas. La discretización espacial se ha realizado por medio de elementos lagrangianos isoparamétricos. Se ha aplicado la integración numérica de Simpson para obtener las matrices elementales, y para la integración temporal se han utilizado diferentes esquemas en diferencias finitas, comprobándose el modelo con diferentes ejemplos.

Aplicación del método de los elementos finitos a problemas de interpolación

Como es conocido, no siempre admite solución el problema de interpolación polinómica de Birkhoff. En este trabajo se presenta un método como alternativa a este tipo de problemas para la obtención de interpolantes con unas determinadas conclusiones de continuidad y en el cual el criterio de aproximación es el de minimización de un cierto funcional real utilizando el método de los elementos finitos. Se describe el método empleado, así como diversos ejemplos l-D y la extensión a problemas 2-D.

Aplicación de método de los elementos finitos a los problemas de amplificación dinámica de suelos

En el análisis de estructuras situadas en un emplazamiento de riesgo sísmico, es necesario conocer las característiias de un posible movimiento sísmico actuante, en la cota de cimentación de la estructura. La importancia del papel que juega el suelo circundante a la estructura en la modificación de las características del movimientos sísmico (principalmente en cuanto a amplitudes y contenido frecuencial), es bien conocida, denominándose amplificación sísmica a este fenómeno. En el caso de estructuras importantes, como Centrales Nucleares, presas, etc., o terrenos muy blandos, la presencia de la propia estructura también afecta de forma significativa al movimiento, dando lu gar al fenómeno de la interacción suelo-estructura. En esta comunicación se tratará únicamente del problema de la amplificación de las ondas sísmicas causantes del movimiento del suelo, y del análisis de la respuesta del mismo en ausencia de la estructura. Este movimiento del suelo denominado "movimiento del campo libre"(free field motion) puede, bien ser utilizado directamente en el análisis de la estructura, despreciando los efectos de la interacción, o emplearse como "input" en ciertas técnicas de análisis de interacción suelo-estructura.

Modelos de elementos finitos con aplicación tridimensional de la analogía fluido-fibra, para predecir el comportamiento mecánico de piezas estructurales de madera bajo la presencia de los principales tipos de nudos

Los nudos son los defectos que más disminuyen la resistencia de piezas de madera en la escala estructural, al ocasionar no solo una discontinuidad material, sino también la desviación de las fibras que se encuentran a su alrededor. En la década de los 80 se introdujo la teoría de la analogía fluido-fibra, como un método que aproximaba adecuadamente todas estas desviaciones. Sin embargo en su aplicación tridimensional, nunca se consideraron las diferencias geométricas en el sentido perpendicular al eje longitudinal de las piezas estructurales, lo cual imposibilitaba la simulación numérica de algunos de los principales tipos de nudos, y disminuía la precisión obtenida en aquellos nudos en los que la modelización sí era viable. En este trabajo se propone un modelo programado en lenguaje paramétrico de un software de elementos finitos que, bajo una formulación en tres dimensiones más general, permitirá estudiar de forma automatizada el comportamiento estructural de la madera bajo la influencia de los principales tipos de nudos, a partir de la geometría visible de los mismos y la posición de la médula en la pieza, y el cual ha sido contrastado experimentalmente, simulando de forma muy precisa el comportamiento mecánico de vigas sometidas a ensayos de flexión a cuatro puntos. Knots are the defects that most reduce the strength of lumber at the structural level, by causing not only a material discontinuity but also the deviation of the fibers that surround them. In the 80's it was introduced the theory of the flow-grain analogy as a method to approximating adequately these deviations. However, in three-dimensional applications, geometrical differences in the direction perpendicular to the longitudinal axis of the structural specimens were never considered before, which prevented the numerical simulation of some of the main types of knots, and decreased the achieved precision in those kind of knots where modeling itself was possible. This paper purposes a parametric model programmed in a finite element software, in the way that with a more general three-dimensional formulation, an automated study of the structural behavior of timber under the influence of the main types of knots is allowed by only knowing the visible geometry of such defects, and the position of the pith. Furthermore that has been confirmed experimentally obtaining very accurately simulations of the mechanical behavior of beams under four points bending test.

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.

Amplificación sísmica: una aproximación por Elementos Finitos

En las páginas que siguen se presenta el estudio de uno de los más típicos problemas de dinámica estructural, cual es la obtención de la respuesta de una estructura excitada por un movimiento de la base. Este es un caso muy frecuente en ingeniería sísmica, donde el objeto del estudio puede ser el edificio (sometido a un movimiento en la cimentación) o un estrato de terreno sobre fondo rígido. Al objeto de facilitar un soporte intuitivo a la exposicióri, ésta se organiza en base al segundo de los casos citados (estrate en base rígida). La aproximación escogida, elementos finitos, pone de relieve una vez más la potencia y generalidad del método en lo que respecta a la formulación del sistema de equilibrio. La discusión se centra en un aspecto concreto del método: la elección de funciones de forma.

El método de las ecuaciones integrales singulares (B.I.E.M.)

El Método de Ecuaciones Integrales Singulares es una alternativa a los métodos dominantes, como el conocido Método de los Elementos Finitos. La idea esencial es la combinación de las clásicas relaciones recíprocas con la filosofía de la discretización del Método de los Elementos Finitos.

Propuesta de Configuración y de Método de Inspección de Uniones Mixtas Mediante Pernos Conectores

La utilización conjunta de los soportes de acero con los forjados de hormigón armado aúna las ventajas constructivas y arquitectónicas de ambos sistemas de ejecución. El principal inconveniente que se presenta es la complejidad de sus nudos. Las uniones estructurales viga-pilar mediante pernos tipo Nelson o Köco aparecen sólo ocasionalmente en el mundo de la edificación. Los casos en los que esto ocurre normalmente corresponden a obras de especial singularidad, a uniones secundarias o a anclajes de reducida responsabilidad estructural. La presente Tesis Doctoral tiene el objetivo principal de aportar una alternativa a las uniones convencionales entre losas o vigas de hormigón armado con pilares metálicos, utilizando para ello los pernos conectores. En este trabajo se establece un campo de aplicación concreto de las estructuras de edificación en el que la utilización de uniones mediante pernos conectores resulta óptima, dada la posibilidad de industrialización de estos nudos, su sencillez y la alta resistencia de los pernos en la transmisión del cortante. Se ha procedido al estudio teórico-analítico de la resistencia de estas uniones y a su posterior comparación experimental. Se presentan las diferencias obtenidas entre los resultados experimentales y los previstos por las formulaciones teóricas, debidas principalmente a las particularidades de las configuraciones de la unión y al efecto del rozamiento. El excelente comportamiento mecánico obtenido en los resultados experimentales induce a presagiar que la configuración de unión que se propone, junto con el método de cálculo aportado, puede llegar a sustituir a las soluciones convencionales utilizadas hasta el momento en las estructuras mixtas de edificación. Los pernos conectores aportan grandes ventajas al campo de la construcción aunque sus soldaduras, pese a la gran seguridad que las caracteriza, son difícilmente comprobables mediante ensayos no destructivos. Los ensayos destructivos no son en absoluto recomendables para las configuraciones de unión estudiadas en esta Tesis. Esto se debe principalmente al reducido número de pernos que configuran las uniones y a la gran responsabilidad que cada perno asume respecto a la seguridad global de la estructura. Los métodos convencionales de inspección presentan grandes dificultades e inconvenientes para comprobar las soldaduras de los pernos, por lo que se han investigado nuevas posibilidades de inspección y control. Esta Tesis plantea un método alternativo para la inspección de las soldaduras de los pernos conectores. El método desarrollado se basa en la acústica y consiste en el cálculo analítico de los modos propios de vibración del perno y su posterior comparación con las frecuencias naturales obtenidas a partir de una medición experimental. La investigación analítica comprende numerosas simulaciones por elementos finitos en las que se varían las condiciones de contorno para simular los defectos de soldadura. La investigación experimental consiste en la medición acústica de los pernos conectores y su posterior ensayo mediante métodos destructivos y no destructivos. Finalmente se han comparado los resultados analíticos con los obtenidos en el laboratorio y se ha obtenido una clara correlación entre las frecuencias calculadas y la calidad de las soldaduras. De este modo se puede averiguar el estado interno de una soldadura a partir de medición de su frecuencia natural de vibración. Se propone un método de inspección acústica que permite comprobar la calidad de las soldaduras de modo sencillo, económico y de fácil ejecución en obra. En conclusión, esta Tesis presenta una nueva configuración de unión mixta mediante pernos conectores y un completo método de calculo y comprobación. Finalmente se ha desarrollado un método no destructivo para la inspección de los pernos que permite garantizar la calidad de sus soldaduras en obra.

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.