Estudo do fenômeno de instabilidade de estruturas cilíndricas de paredes finas submetidas ao carregamento de pressão externa

Thin walled cylindrical shells are widely used in many areas of industry, including civil, mechanical, nuclear, marine, petroleum and aerospace engineering. The wide application of thin cylindrical shells and the importance of instability phenomenon are the motivation basis to this study, since these factors have a great importance in engineering projects. It is presented a detailed study about the instability of cylindrical shells based on theoretical calculation, which results are compared with finite elements method calculation. The loading and boundary conditions analyzed are based on the most common types verified in real engineering projects and refer respectively to lateral (external) pressure and cylinders with simply supported edges. The calculation based on the finite elements method was executed with ANSYS 13.0 software. The results obtained with this calculation are in good agreement with the analytical theory presented in the technical note NACA No 1341 (BATDORF, 1947) considering a wide range of applicability. On the other hand, the analytical method presented in the book Theory of Elastic Stability (TIMOSHENKO; GERE, 1936) has a very restrict applicability and has presented considerable deviations in a great sort of the analyzed cases

Controle ativo de vibrações em vigas utilizando controle feedforward e elementos espectrais

Pós-graduação em Engenharia Mecânica - FEB

Implementação de aparato experimental para medição de instabilidades tipo Roll Waves em fluidos não- newtonianos

Pós-graduação em Engenharia Mecânica - FEIS

Cargas explosivas sobre pilares y forjados modelados con elementos tipo viga y lámina

Los fenómenos de impacto y explosión sobre estructuras de hormigón tienen efectos en muchos casos catastróficos a pesar de su reducida probabilidad. Las estructuras de hormigón no suelen estar diseñadas para resistir este tipo de solicitaciones dinámicas. El análisis numérico mediante elementos finitos con integración explícita permite una aproximación suficiente a los efectos de la onda explosiva sobre pilares y forjados de estructuras reticuladas de hormigón. Los materiales recientemente implementados en LS-Dyna para hormigón como el CSCM [1], para elementos de continuo 3D, y la formulación que proporciona la debida compatibilidad con los elementos viga de acero dispuestos de forma segregada, permite estudiar de forma realista modelos detallados de pilares y forjados. Pero las limitaciones computacionales hacen inviable emplear estos métodos en estructuras completas. Como alternativa es posible usar modelos de elementos estructurales de vigas y láminas para el análisis de estas estructuras. Sin embargo es necesario un adecuado ajuste de parámetros y propiedades en estos modelos. Este trabajo muestra un método con en el que obtener modelos de elementos estructurales, elementos viga y lámina, usando modelos de material [2] adecuados para ellos, junto a un procedimiento para incluir la armadura de forma adecuada. Utilizando este método es posible representar con suficiente aproximación el comportamiento de modelos detallados realistas de forjados y pilares de estructuras reticuladas de hormigón frente a acciones explosivas, posibilitando el análisis de una estructura completa frente a explosión.

Fenómenos de resonancia asociados al transporte de mercancías líquidas por ferrocarril

Los vagones de mercancías se unen entre sí mediante acoplamientos UIC estándar. Aunque cada plataforma está equipada con un tensor de enganche en ambos extremos, sólo uno de ellos se utiliza, mientras que el otro descansa apoyado en un gancho. Debido a los efectos de oleaje del líquido contenido en las cisternas, este tensor en reposo puede entrar en resonancia con la frecuencia natural de algunos fluidos, llegando a oscilar de forma apreciable. En ocasiones llega a saltar fuera de su alojamiento, provocando ciertas molestias durante la circulación, aunque sin entrañar ningún riesgo para la seguridad de la marcha. Con el propósito de simular la influencia de diversos factores sobre el comportamiento dinámico asociado a estos fenómenos de oleaje, se ha desarrollado un modelo de la interacción fluido-vehículo, empleado técnicas de modelado de sistemas multicuerpo y de elementos finitos.

Análisis estructural 3D de los invernaderos multitunel del sudeste de España. Elementos verticales

Este trabajo analiza la capacidad estructural de los elementos verticales que componen el pórtico resistente más limitativo de la variante de invernadero multitúnel más utilizada en el sudeste de España. Se estudia la influencia de la longitud y anchura del invernadero en la distribución de tensiones que se produce en dichos elementos. Se ha realizado un modelo tridimensional de elementos finitos con el programa ROBOT MILLENNIUM en su versión 20.1 en el que se ha representado la estructura de barras que conforma el invernadero. Se han analizado los resultados de tensiones normales máximas de 16 estructuras con geometría de túnel similar, pero de diferente anchura y longitud, de tal forma que la superficie en planta oscila entre los 600 m 2 y los 38.400 m 2 . El cálculo de las acciones existentes sobre los invernaderos se ha realizado según la norma UNE-EN 13031-1. Los resultados numéricos obtenidos han sido estudiados para determinar la influencia de las dimensiones del invernadero en la distribución de tensiones. Aunque se observan significativas reducciones en las tensiones con la adición de túneles, no se puede generalizar la influencia del aumento de las dimensiones del invernadero en los cambios tensionales del mismo. Los resultados obtenidos se consideran de gran utilidad para la optimización de las estructuras de los invernaderos multitúnel así como para garantizar la seguridad mecánica de las mismas

Análisis estructural 3D de los invernaderos multitunel del sudeste de España. Elementos horizontales

Este trabajo analiza la capacidad estructural de los Este trabajo analiza la capacidad estructural de los elementos horizontales que componen el pórtico resistente más limitativo de la variante de invernadero multitúnel más utilizada en el sudeste de España. Se analiza la influencia de la anchura y longitud del invernadero en la distribución de tensiones que se producen en dichos elementos. Se ha realizado un modelo tridimensional de elementos finitos con el programa ROBOT MILLENNIUM en su versión 20.1 en el que se ha representado la estructura de barras que conforma el invernadero. Se han analizado los resultados de tensiones normales máximas de 16 estructuras con geometría de túnel similar, pero de diferente anchura y longitud, de tal forma que la superficie en planta oscila entre los 600 m 2 y los 38.400 m 2 .El cálculo de las acciones existentes sobre los invernaderos se ha realizado según la norma UNE-EN 13031-1. El aumento de las dimensiones del invernadero conlleva importantes variaciones en las tensiones máximas registradas por canalones. Sin embargo, en arcos y tirantes de cercha, las tensiones permanecen constantes. Los resultados obtenidos se consideran de gran utilidad para la optimización de las estructuras de los invernaderos multitúnel así como para garantizar la seguridad mecánica de las mismas.

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.

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

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.

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.

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.

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.