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.

Los Métodos de las Características y Elementos Finitos en propagación de ondas

Se presentan a continuación algunas consideraciones acerca de la propagación de ondas en medios unidimensionales elásticos y viscoplásticos. El objetivo básico es la comparación entre diferentes modelizaciones del sistema establecidas siguiendo las técnicas de Elementos Finitos. Como criterio comparativo se utilizan resultados obtenidos mediante el método de las características, que, para el caso descrito proporciona la solución más aproximada. La justificación del trabajo radica en que las técnicas de Elementos Finitos son fácilmente generalizables a problemas bi o tridimensionales,mientras que las de características pierden rápidamente la ventaja de su simplicidad.

Representación de isostáticas para la mejora del mallado de elementos finitos en sólidos elásticos planos

Uno de los problemas primordiales en el cálculo por elementos finitos ha sido la obtención del mallado óptimo tal que se minimice el error obtenido, pudiendo distinguirse los siguientes procedimientos: - Aumento del número de nudos de la malla, fundamentalmente en las zonas del modelo donde aparece un error mayor. - Incrementando el grado de los polinomios de interpolación en los elementos donde el modelo presenta un error mayor. - Una combinación entre el primer y el segundo procedimiento. Según los trabajos realizados en la tesis doctoral de D. Rubén Martínez Marín1, se llega a la conclusión de que, tras medir el error del mallado por dos procedimientos distintos; los nudos de la malla óptima se sitúan a lo largo de las líneas isostáticas. Lo destacable de este resultado es que se obtiene sin variar el número de nudos iniciales, y sin incrementar el grado de los polinomios de interpolación; es decir, únicamente buscando la posición óptima de los nudos. Así, en el presente documento se plantea la realización de dos cálculos por elementos finitos; uno con un mallado convencional formado por elementos rectangulares, y otro con un mallado isostático, y la comparación de su error. Los dos mallados tendrán un número similar de nudos. Como modelo se utiliza una viga en ménsula de 6 m de longitud y 2 m de canto con una carga puntual vertical en su extremo. Todos los algoritmos utilizados se encuentran programados en MATLAB. El presente documento se estructura en las siguientes partes: - Capítulo 1.- Descripción de los trabajos. Donde se realiza un resumen de los trabajos realizados en la creación del presente documento. - Capítulo 2.- Trabajos previos. En el que se resumen los trabajos realizados por otros autores antecedentes del presente documento. - Capítulo 3.- Fundamentos teóricos. Donde se explican las bases teóricas que se van a aplicar en la creación del algoritmo y en su análisis. - Capítulo 4.- Descripción del algoritmo implementado en este trabajo. En este capítulo se analiza la estructura del algoritmo empleado. Incluye diagramas de proceso del programa base y de las principales subrutinas. - Capítulo 5.- Resultados y discusión. Donde se realiza la comparación del error del mallado convencional y del mallado isostático; por un lado comparando las flechas obtenidas en el extremo de la viga en voladizo con el valor exacto de la flecha, y por otro lado utilizando el Error Cuadrático Medio de las tensiones medias. Se termina con un análisis crítico de los resultados. - Capítulo 6.- Conclusiones y futuras líneas de investigación.

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.).

Curso básico de programación del método de los elementos finitos

Estas notas que se publican a continuación corresponden a un curso de postgrado impartido durante el primer semestre del año 1983. El interés mostrado por los asistentes a dicho curso nos ha animado a escribir un resumen de las clases. Este libro supone un conocimiento te6rico de las ideas básicas del método de los elementos finitos. No obstante en una primera lección se resumen y ordenan aquellos aspectos mas importantes, que serán utilizados en lecciones sucesivas. En estas se desarrolla un programa de computador muy sencillo -sin complicaciones informáticas que obscurezcan la simplicidad del método- y se analiza de un modo detallado -en forma de organigramas y listados comentados- las distintas rutinas en lenguaje FORTRAN de este programa. Asimismo, y respetando el carácter elemental de la exposición se abren algunas posibilidades de ampliación y nuevos desarrollos del método. Algunos ejercicios y ejemplos al final de cada capítulo se espera permitan clarificar los puntos mas conflictivos del método. Finalmente se reúne en un apéndice, los distintos programas que se han mostrado en las sucesivas lecciones y que con objeto de que puedan ser procesados en microcomputadores se han traducido al lenguaje BASIC. Creemos y la experiencia del curso así nos la ha confirmado, que el método de elementos finitos se debe enseñar y aprender mediante la praxis y presentar los sucesivos desarrollos del método de un modo motivado como solución a problemas numéricos e informáticos que aparecen en su desarrollo. Si las lecciones que aquí se presentan permiten transmitir mejor estas ideas, los autores se sentirán más que recompensados por el trabajo que ha supuesto dar a luz a esta publicación.

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.

Problemas de elementos finitos : curso monográfico : 26-30 de mayo de 1980

Curso monográfico realizado del 26 al 30 de mayo de 1980, en el que se resuelven una serie de problemas de elementos finitos. Los temas tratados son: 1. Problemas generales. Ancho de banda. Diseño de mallas. 2.-Términos numéricos de ecuaciones diferenciales. 3.-Funciones de forma. 4 y 5 - Matrices de rigidez y cargas equivalentes. 6.- Aplicaciones. 7.- Ejercicios de cálculo numérico. 8.-Ejercicios de programación

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.

Optimización de mallas bidimensionales en elementos finitos

La solución al problema de encontrar la malla óptima en Elementos Finitos (EF), con un determinado número de grados de libertad, presenta un indudable interés en la aplicación del método. En la actualidad, el problema se plantea en términos de un proceso que permite obtener una mejor malla de elementos finitos a partir de una inicial. La nueva malla se diseña matemáticamente (remallado) de forma que el error del método sea lo más uniforme posible en todo el dominio de cálculo. Sin embargo, esta técnica de indudable interés y aplicación, al aumentar el número de grados de libertad (gdl) de la aproximación, no permite deducir de un modo directo el problema de la malla óptima condicionada a un número fijo de gdl. Con la solución de este problema se podrán deducir algunos criterios y recomendaciones para el diseño de una malla de elementos finitos, que exigirá, en general, en un proceso de remallado, modificaciones menores. Para problemas unidimensionales (barras y pilares simples), se pueden encontrar soluciones analíticas. Para problemas 2-D más complicados (tensión y deformación plana), se han utilizado métodos numéricos para obtener la malla óptima. Existen varios criterios de optimización, aquí se utiliza el del mínimo de la energía potencial total (EPT). Algunos ejemplos ilustrativos del método de optimización se presentan, indicándose algunas conclusiones.

Mallas isostáticas bidimensionales en elementos finitos

Se plantea el problema de conseguir un método que, partiendo de una malla inicial en elementos finitos, genere una malla próxima a la óptima, conservando el mismo número de grados de libertad y, por consiguiente, sin penalizar los tiempos de análisis por ordenador. La bondad de una malla se mide por un funcional determinado (Energía Potencial Total, Error Cuadrático Medio, etc.) La Técnica de gradiente descendente, aplicada al funcional de Energía Potencial Total (1) permite, a partir de una inicial dada, obtener una malla mejorada. No obstante, este método requiere un esfuerzo de computación muy grande. Como consecuencia de la aplicación del método del gradiente descendente a numerosos casos, se ha observado que la geometría que adopta la malla mejorada, se aproxima a la de una malla tal que sus nudos se sitúen sobre las líneas isostáticas (envolventes de las tensiones principales) y generen elementos regulares (de lados iguales) o cuasirregulares. La conclusión más importante es, precisamente, que a partir de una malla inicial, razonablemente regular, se puede realizar un único análisis y definir las isostáticas correspondientes a este modelo. Ajustando una malla, con el mismo número de nudos que la inicial, a las líneas isostáticas con nudos situados de forma regular, se obtiene otra que está próxima a la óptima. A la malla así obtenida, se la denomina isostática isométrica. Esta conclusión se ha probado que es válida para los diferentes funcionales que se utilizan para evaluar la bondad de la malla en elementos finitos.

Una familia de elementos finitos clase C.1. Aplicación a la flexión de placas de Kirchhoff

Se hace una pequeña introducción y después un estudio sobre las posibilidades y limitaciones en análisis de placas delgadas de elementos simples polinómicos de clase C.1. Se expone una familia jerárquica de dichos elementos, que se aplica a varios casos particulares. En base a estos se deducen algunas conclusiones, especialmente en lo que se refiere a eficacia computacional. Al final se proponen trabajos a realizar a partir de los datos existentes.

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.

Métodos numéricos y estadísticos de caracterización de uniones tubulares soldadas para su aplicación en modelos de elementos finitos de estructuras de vehículos

El análisis de estructuras mediante modelos de elementos finitos representa una de las metodologías más utilizadas y aceptadas en la industria moderna. Para el análisis de estructuras tubulares de grandes dimensiones similares a las sobrestructuras de autobuses y autocares, los elementos de tipo viga son comúnmente utilizados y recomendados debido a que permiten obtener resultados satisfactorios con recursos computacionales reducidos. No obstante, los elementos de tipo viga presentan importante desventaja ya que las uniones modeladas presentan un comportamiento infinitamente rígido, esto determina un comportamiento mas rígido en las estructuras modeladas lo que se traduce en fuentes de error para las simulaciones estructurales (hasta un 60%). Mediante el modelado de uniones tubulares utilizando elementos de tipo área o volumen, se pueden obtener modelos más realistas, ya que las características topológicas de la unión propiamente dicha pueden ser reproducidas con un mayor nivel de detalle. Evitándose de esta manera los inconvenientes de los elementos de tipo viga. A pesar de esto, la modelización de estructuras tubulares de grandes dimensiones con elementos de tipo área o volumen representa una alternativa poco atractiva debido a la complejidad del proceso de modelados y al gran número de elementos resultantes lo que implica la necesidad de grandes recursos computacionales. El principal objetivo del trabajo de investigación presentado, fue el de obtener un nuevo tipo de elemento capaz de proporcionar estimaciones más exactas en el comportamiento de las uniones modeladas, al mismo tiempo manteniendo la simplicidad del procesos de modelado propio de los elementos de tipo viga regular. Con el fin de alcanzar los objetivos planteados, fueron realizadas diferentes metodologías e investigaciones. En base a las investigaciones realizadas, se obtuvo un modelo de unión viga alternativa en el cual se introdujeron un total seis elementos elásticos al nivel de la unión mediante los cuales es posible adaptar el comportamiento local de la misma. Adicionalmente, para la estimación de las rigideces correspondientes a los elementos elásticos se desarrollaron dos metodologías, una primera basada en la caracterización del comportamiento estático de uniones simples y una segunda basada en la caracterización del comportamiento dinámico a través de análisis modales. Las mejoras obtenidas mediante la implementación del modelo de unión alternativa fueron analizadas mediante simulaciones y validación experimental en una estructura tubular compleja representativa de sobrestructuras de autobuses y autocares. En base a los análisis comparativos realizados con la uniones simples modeladas y los experimentos de validación, se determinó que las uniones modeladas con elementos de tipo viga son entre un 5-60% más rígidas que uniones equivalentes modeladas con elementos área o volumen. También se determinó que las uniones área y volumen modeladas son entre un 5 a un 10% mas rígidas en comparación a uniones reales fabricadas. En los análisis realizados en la estructura tubular compleja, se obtuvieron mejoras importantes mediante la implementación del modelo de unión alternativa, las estimaciones del modelo viga se mejoraron desde un 49% hasta aproximadamente un 14%. ABSTRACT The analysis of structures with finite elements models represents one of the most utilized an accepted technique in the modern industry. For the analysis of large tubular structures similar to buses and coaches upper structures, beam type elements are utilized and recommended due to the fact that these elements provide satisfactory results at relatively reduced computational performances. However, the beam type elements have a main disadvantage determined by the fact that the modeled joints have an infinite rigid behavior, this shortcoming determines a stiffer behavior of the modeled structures which translates into error sources for the structural simulations (up to 60%). By modeling tubular junctions with shell and volume elements, more realistic models can be obtained, because the topological characteristics of the junction at the joint level can be reproduced more accurately. This way, the shortcoming that the beam type elements present can be solved. Despite this fact, modeling large tubular structures with shell or volume type elements represents an unattractive alternative due to the complexity of the modeling process and the large number of elements that result which imply the necessity of vast computational performances. The main objective of the research presented in this thesis was to develop a new beam type element that would be able to provide more accurate estimations for the local behavior of the modeled junctions at the same time maintaining the simplicity of the modeling process the regular beam type elements have. In order to reach the established objectives of the research activities, a series of different methodologies and investigations have been necessary. From these investigations an alternative beam T-junction model was obtained, in which a total of six elastic elements at the joint level were introduced, the elastic elements allowed us to adapt the local behavior of the modeled junctions. Additionally, for the estimation of the stiffness values corresponding to the elastic elements two methodologies were developed, one based on the T-junction’s static behavior and a second one based on the T-junction’s dynamic behavior by means of modal analysis. The improvements achieved throughout the implementation of this alternative T-junction model were analyzed though mechanical validation in a complex tubular structures that had a representative configuration for buses and coaches upper structures. From the comparative analyses of the finite element modeled T-junctions and mechanical experimental analysis, was determined that the beam type modeled T-junctions have a stiffer behavior compared to equivalent shell and volume modeled T-junctions with average differences ranging from 5-60% based on the profile configurations. It was also determined that the shell and volume models have a stiffer behavior compared to real T-junctions varying from 5 to 10% depending on the profile configurations. Based on the analysis of the complex tubular structure, significant improvements were obtained by the implementation of the alternative beam T-junction model, the model estimations were improved from a 49% to approximately 14%.