3700x
000434
2024-01-26
Estructura

Elementos finitos en RFEM

Elementos en 1D

Para elementos de barra, se asume que la sección permanece plana cuando se deforma. Los 'Elementos en 1D' se utilizan para representar vigas, cerchas, nervios, cables y acoplamientos rígidos. Un elemento de barra 1D tiene un total de doce grados de libertad − seis al inicio y seis al final del elemento, cada uno para los desplazamientos (ux, uy, uz ) y giros (φx, φy, φz ). Al calcular los datos estructurales linealmente, la tracción, compresión y torsión se expresan como funciones lineales del eje x de la barra, independientemente de la flexión y el cortante. Estas funciones se aproximan mediante un polinomio de tercer orden en x, incluyendo la influencia de las tensiones tangenciales de los esfuerzos cortantes Vy y Vz. La matriz de rigidez KL (12, 12) describe el comportamiento lineal de los elementos 1D. La interacción mutua entre el esfuerzo axil y la flexión en el caso de problemas geométricamente no lineales se expresa en la matriz de rigidez KNL (12, 12). Para obtener más información, consulte [1] y [2].

Para cálculos según el análisis de grandes deformaciones, se recomienda utilizar un Refinamiento de malla de líneas para que los resultados de la barra se puedan determinar con la precisión correspondiente.

Elementos en 2D

Para las superficies, se crean Elementos en 2D. Donde no se pueden usar elementos cuadrangulares, el generador de malla inserta elementos triangulares.

Los grados de libertad de los elementos cuadrangulares y triangulares en los nudos de esquina son los mismos que para los elementos en 1D: grados de libertad del desplazamiento (ux, uy, uz) y giro (φx, φy, φz). Esto asegura la compatibilidad entre los elementos 1D y 2D en los nudos. Los parámetros se definen en el sistema de coordenadas local plano de los elementos y se convierten en el sistema de coordenadas global al componer la matriz de rigidez global.

La imagen Elementos cuadrangulares para láminas muestra la aproximación para elementos planos de láminas. Se basa en la teoría de Mindlin/Reissner. Para asegurar el acoplamiento directo con los elementos de la barra, se elige una aproximación cuadrada en el plano de la lámina (ux, uy ). Al eliminar los nudos intermedios, se crea un elemento de cuatro nudos con grados de libertad de giro adicionales en las direcciones x, y y z. Por lo tanto, es posible un acoplamiento directo con elementos de viga para elementos de diafragma. Basado en una interpolación combinada de desplazamientos transversales, giros de la sección y deformaciones por cortante transversal, los métodos de Bathe y Dvorkin [3] presenta los elementos MITC ( interpolación mixta de componentes tensoriales ): MITC3+ para triángulos, MITC4 para cuadrángulos.

En la actualidad, se consideran los elementos de barra al resolver directamente la ecuación diferencial según el análisis de segundo orden. Los efectos de alabeo no se pueden considerar cuando se usa la torsión de Saint Venant.

Los cálculos de las membranas se basan en los principios de Bergan [4] [5] [6]. Para elementos triangulares, por ejemplo, las funciones básicas se subdividen en tres deformaciones de cuerpo rígido, tres condiciones de deformación constante y tres gradientes lineales especiales de tensión y deformación. El campo de deformación dentro de un elemento es cuadrático y el campo de tensiones es lineal. La matriz de rigidez del elemento KL se transforma posteriormente en nueve parámetros colectivos de los tiposux, uy, φz. Los componentes de esta matriz se insertan junto con los componentes para flexión y cortante en la matriz de rigidez general (18, 18). Esta matriz es el resultado del concepto de Lynn/Dhillon. De esta manera, se aplican las llamadas placas de Mindlin (placas con una deformación por cortante distinta) y se analizan según Timoshenko. De esta manera, RFEM encuentra la solución correcta tanto para placas gruesas como finas (placas según el método de Navier).

En caso de problemas geométricamente no lineales, no es posible dividir la condición tensión-deformación en un estado plano y en flexión/cortante mencionada anteriormente. Las influencias mutuas de estos estados se consideran en la matriz KNL. RFEM utiliza una forma bastante simple pero eficaz de la matriz KNL que se basa en los enfoques deZienkiewicz [7]. Se aplica la componente cuadrada ε2 del tensor de deformación de Green/Lagrange ε = ε1 + ε2. Se suponen una distribución linear de uz(x, y) de la condición de tensiones planas y distribuciones lineares de ux(x, y) y uy(x, y) de la interacción con flexión. Esta suposición es posible debido a que el efecto principal de la interacción depende de la primera derivación de la ecuación diferencial y porque la influencia de componentes de un orden superior disminuye con la división en elementos más pequeños. La exactitud de este procedimiento ha sido probada en varios análisis numéricos.

Información

El espesor de los elementos para elementos de lámina debe ser significativamente menor que las dimensiones. Si no es así, se recomienda modelar los objetos como sólidos. Además, con respecto a los elementos de la lámina, evite introducir tensiones de torsión de forma selectiva: el grado de libertad de giro de la superficie normal reacciona con mucha sensibilidad.

Elementos en 3D

Los Elementos en 3D se utilizan para los sólidos. Los siguientes tipos de elementos están implementados en RFEM: tetraedro, pentaedro (prisma, pirámide) y hexaedro. Para información detallada sobre elementos y matrices aplicados, consulte los 'Elementos finitos en 3D con grados de libertad al giro' de Sevcik [8] para más detalles. Puede solicitarle esta documentación a Dlubal Software.

Generalmente, se deben considerar todos los grados de libertad al giro como críticos para sólidos. Como la deformación de un sólido se determina únicamente a partir de los vectores de desplazamiento, el giro de un nudo de malla (por ejemplo, debido a una torsión introducida de forma singular) no afecta a la deformación dentro del sólido.


Referencias
  1. Vladimír Kolář et al. Bemessung von zwei- und dreidimensionalen Strukturen mit FEM. Springer-Verlag, New York / Wien, 1975. Kapitel 1 (1D-Element) und 6 (Variationsprinzip)
  2. Vladimír Kolář und Ivan Němec. Finite Element Analysis of Structures. United Nations Development Program, Economic Com. for Europe, Workshop on CAD Techniques, Prague - Geneva, 1984.
  3. Eduardo N. Dvorkin und Klaus-Jürgen Bathe. A continuum mechanics based four-node shell element for general non-linear analysis. Engineering Computations, 1, 1984.
  4. P. G. P.rgan. P.nite Elements Based on Energy Orthogonal Functions. P.ternational Journal for Numerical Methods in Engineering, 15, 1980.
  5. P. G. P.rgan und M. P. P.gård. P.nite Elements With Increased Freedom in Choosing Shape Functions. P.ternational Journal for Numerical Methods in Engineering, 20, 1984.
  6. P. G. P.rgan und Carlos A. P.lippa. P.Triangular Membrane Element With Rotational Degrees of Freedom. P.mputer Methods in Applied Mechanics and Engineering, 50, 1985.
  7. Olgierd Cecil Zienkiewicz. The Finite Element Method in Engineering Science. Mc Graw-Hill, London 3. Auflage, 1979. Chapter 18 - 19 (Nonlinear Problems).
  8. I. Sevčík. I. Finite Elements with Rotational Degrees of Freedom. I.M Consulting s.r.o, Brno.
Sección original