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.

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.

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.

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.

Aplicação do método da resistência directa a vigas contínuas de aço

Dissertação para a obtenção do grau de Mestre em Engenharia Civil

Análise pushover com elementos shell em paredes estruturais

Dissertação para obtenção do Grau de Mestre em Engenharia Civil - Perfil Estruturas

Determinação da pressão de colapso em tubos API 5L X56 aplicando o método R6

O presente trabalho tem como objetivo determinar a pressão de colapso de tubos API 5L X56 de parede fina com trincas longitudinais internas e externas submetidas a pressão interna. Foi utilizado o método R6 conforme a norma BS7910:1999 e a Análise por Elementos Finitos. Foram feitos ensaios de tração e de caracterização da tenacidade do material utilizado por integral J, assim como o modelamento por Elementos Finitos do tubo e a aplicação dos “Failure Assessment Diagram” FAD para determinar a pressão de colapso. Verificou-se que para trincas internas a pressão de colapso é levemente maior que em trincas externas. Isto é mais acentuado para trincas com uma profundidade maior que a metade da parede do tubo. As predições da pressão de colapso em função do tamanho da trinca do método R6 em comparação com alguns resultados experimentais mostraram uma precisão satisfatória, porem em alguns casos não conservadora.

Análise estática e dinâmica, linear e não-linear geométrica, através de elementos hexaédricos de oito nós com um ponto de integração

Para a análise estática e dinâmica, linear e não-linear de placas, cascas e vigas, implementa-se neste trabalho o elemento hexaédrico com integração reduzida, livre de travamento volumétrico e travamento de cisalhamento e que não apresenta modos espúrios. Na formulação do elemento, utiliza-se apenas um ponto de integração. Desta forma, a matriz de rigidez é dada de forma explícita e o tempo computacional é significativamente reduzido, especialmente em análise não-linear. Os modos espúrios são suprimidos através de um procedimento de estabilização que não exige parâmetros especificados pelo usuário. Para evitar o travamento de cisalhamento, desenvolve-se o vetor de deformações num sistema co-rotacional e remove-se certos termos não constantes nas componentes de deformações de cisalhamento. O travamento volumétrico é resolvido fazendo-se com que a parte dilatacional (esférica) da matriz gradiente seja avaliada apenas no ponto central do elemento. Como a eliminação do travamento de cisalhamento depende de uma abordagem no sistema local, emprega-se um procedimento co-rotacional para obter o incremento de deformação no sistema local e atualizar os vetores de tensões e forças internas na análise não-linear Para a solução das equações de equilíbrio na análise estática, utilizam-se métodos diretos baseados na eliminação de Gauss ou métodos iterativos de Gradientes Conjugados Precondicionado elemento-por-elemento (EBE). Para a análise dinâmica, as equações de equilíbrio são integradas através do método explícito de Taylor-Galerkin ou do método implícito de Newmark. Através de exemplos numéricos demonstra-se a eficiência e o potencial do elemento tridimensional na análise de casca, placas e vigas submetidas a grandes deslocamentos e grande rotações. Os resultados são comparados com trabalhos que utilizam elementos clássicos de placa e casca.

El método del mecano y análisis isogeométrico

[EN]En este trabajo presentamos avances del método del mecano para la generación de mallas de tetraedros y la parametrización volumétrica de sólidos. El método combina varios procedimientos: una transformación desde la frontera del mecano (espacio paramétrico sencillo) a la superficie del sólido (espacio físico), un algoritmo de refinamiento local de mallas de tetraedros, y un desenredo y suavizado simultáneo de mallas. La base del método consiste en la definición de una transformación volumétrica biyectiva entre el mecano y el sólido. Se mostrarán resultados usando elementos finitos adaptativos en varios problemas medioambientales : campos de viento, contaminación atmósférica, etc…

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.

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