CALCULO DE ESTRUCTURAS MEDIANTE EL EMPLEO DE

MODELOS DE ELEMENTOS FINITOS

Lara Menéndez, Rafael

Báculos, S.A. C/ Juan de Mariana, 4 28045 – Madrid, Teléfonos: 91 467 05 20 / 91
527 92 58; Fax: 91 467 35 25; e-mail: basa@retemail.es

RESUMEN

Se pretende explicar el procedimiento de cálculo de estructuras mediante

modelo de elementos finitos a través del empleo de varios ejemplos de soportes de
diseño de alumbrado ejecutados por BACULOS, S.A.

Se expondrá la forma de trabajar desde el diseño del modelo en 3D y creación

del modelo de elementos finitos hasta la obtención de resultados tanto de forma
numérica como visual en pantalla.

Los programas de calculo de estructuras mediante elementos finitos

proporcionan una potente herramienta de trabajo. Se emplean para la optimización
de resultados mediante la creación de un modelo sólido y el posterior análisis
elástico – tensional para un estado de cargas que simule su implantación en
servicio.

INTRODUCCION

Con la llegada de los ordenadores se ha facilitado en gran medida la solución de

los sistemas de ecuaciones lineales, lo que ha representado un retorno a los
procedimientos tradicionales de cálculo. Estos métodos matriciales implican un gran
número de operaciones repetitivas particularmente adecuadas para el cálculo

1
automático y han sido formuladas con el fin de obtener las mayores ventajas de los
ordenadores.

El cálculo de ciertas estructuras se complica si su forma geométrica y/o cargas

adoptan formas irregulares en cuyo caso es de aplicación el método de los
elementos finitos consistente en sustituir la estructura por un modelo matemático de
un cierto numero de elementos independientes, denominados elementos finitos o
discretos, unidos entre si a través de un número finito de puntos o nodos

1. LOS ELEMENTOS FINITOS:

1.1 Breve introducción histórica:

• 1950 Aplicación de métodos matriciales en análisis de estructuras.

• 1950 Introducción de los desplazamientos y método de matriz de
rigidez para estructuras aeroespaciales complejas.
• 1960 Nacimiento del nombre “Elementos Finitos”.
• 1970 Aplicaciones de modelos de Elementos Finitos principalmente en
el campo aeroespacial y la industria del automóvil.
• 1977 Los Modelos de Elementos Finitos son la herramienta numérica
para análisis estructural.

1.2 Terminología. Conceptos previos:

- CLOUGH: “Los elementos finitos no son simple trozos segregados de la

estructura original, sino modelos específicos elásticos obligados a
deformarse de tal manera que tienda a mantenerse toda la continuidad del
conjunto”
- Las estructuras complejas se subdividen en puntos discretos (NODOS) los
cuales están interconectados entre si mediante elementos.
- Cada NODO tiene 6 posibilidades independientes de movimiento, 3
traslaciones y 3 rotaciones.
- Estas posibilidades de movimiento independientes se denominan grados
de libertad (DOF).
- Los desplazamientos de los nodos son compatibles con la deformación del
elemento.

2
 - Las fuerzas aplicadas a la estructura, se introducen como un sistema de
cargas equivalentes de cargas sobre los nodos de tal manera que se
simula las acciones reales que van a sufrir.
- Las condiciones de apoyo de la estructura, se introducen como
restricciones en giro y desplazamiento sobre los nodos.

1.3 ¿Cómo se crea un modelo de elementos finitos? Discretización:

A cada elemento finito estudiado aisladamente se le aplica la teoría clásica de

cálculo, lo que es posible dada la forma sencilla elegida (triángulos, rectángulos,…)
estableciéndose las condiciones de contorno y equilibrio a través de los nodos. Esta
discretización de la estructura permite su resolución planteando un sistema de
ecuaciones lineales y se aplica fácilmente a cualquier estructura por complicadas
que sean sus características geométricas y condiciones de carga. Para llevar a cabo
esto, se necesitan un gran número de operaciones matemáticas, que debido a
naturaleza repetitiva se adaptan con total normalidad a la programación numérica y
a su resolución con ordenador.

Fundamentalmente, la aplicación del método consiste en plantear para cada

elemento finito la matriz de rigidez que relaciona las fuerzas con las deformaciones.
A continuación se procede al ensamblaje de la matriz de rigidez total para toda la
estructura.

Los elementos finitos se enlazan entre ellos solamente mediante los nodos, lo
que significa que las condiciones de contorno de la estructura no se satisfacen a lo
largo de los lados. No obstante, a veces, la elección adecuada de un modelo de
deformación para los elementos finitos puede satisfacer esta condición en algunos
elementos o en todos.

En general, dado que el método de cálculo mediante elementos finitos es un

procedimiento aproximado, la precisión requerida aumenta directamente con el
número de elementos empleados. Evidentemente, un mayor número de elementos
necesita mayor tiempo de cálculo por lo que aumentaría el costo.

Resulta ventajoso en algunos casos utilizar una discretización gradual a fin de

conseguir un estudio mas detallado en aquellas zonas donde se tiene previsto que
puedan concentrarse tensiones. Esta discretización resulta eficaz y rebaja el tiempo
de cálculo sin perder precisión. No es posible concretar el número de elementos que
se requieren para obtener en cada caso la solución más satisfactoria, ya que va a

3
depender de la estructura objeto de estudio. La elección de la subdivisión más
conveniente ha de realizarse de acuerdo con la experiencia, basándose, si es
posible, en resultados obtenidos mediante ensayos. En cualquier caso, se deben
poner en practica soluciones que hagan intervenir mallas de dife rentes dimensiones
con el fin de asegurar la convergencia de los resultados.

Las cargas externas aplicadas a la estructura se sustituyen por sistemas de

fuerzas equivalentes concentradas en los nodos. De existir cargas puntuales ha de
distribuirse la malla de manera que los nodos coincidan con sus puntos de
aplicación.

2. PROGRAMA DE CALCULO DE ESTRUCTURAS MEDIANTE EL EMPLEO DE

MODELOS DE ELEMENTOS FINITOS.

2.1 Introducción:

Mediante programas de elemento finitos se pueden calcular todo tipo de

estructuras. Los programas de elementos finitos ofrecen resultados exactos frente a
resultados obtenidos mediante cálculos manuales.

Para la resolución de estructuras por el método de cálculo mediante elementos

finitos se precisa el planteamiento de las matrices de rigidez de todos los elementos
discretos y ensamblarlos adecuadamente con el fin de formar la matriz de rigidez [K]
de la totalidad de la estructura.

2.2 Etapas para el cálculo con ordenador:

Se distinguen tres bloques bien diferenciados que componen el proceso de

calculo de estructuras mediante el empleo de elementos finitos. Se enumeran:

• PREPROCESADOR: El objetivo principal del preprocesador es

conseguir la estructura mallada mediante nodos y elementos. Estos,
son los puntos de control de la estructura mallada a través de los
cuales se introducen datos y se obtienen resultados.
• SOLVER: Parte interna del proceso. Resolución de sistemas de
ecuaciones a partir de los datos introducidos.
• POSTPROCESADOR: Análisis de resultados.

4
 En el esquema adjunto se resumen las etapas principales en base a la
programación con ordenador cuando se emplea el método de cálculo de estructuras
mediante elementos finitos.

CREACIÓN GEOMETRÍA
Nodos y elementos

PREPROCESADOR
ENTRADA DE DATOS
Propiedades del material, cargas
y condiciones de contorno

FORMACIÓN DE LAS MATRICES DE RIGIDEZ

DE LOS ELEMENTOS FINITOS

ENSAMBLAJE DE LA MATRIZ DE
RIGIDEZ DE LA
TOTALIDAD DE LA ESTRUCTURA
SOLVER

APLICACION CONDICIONES DE CONTORNO

CALCULO DE DESPLAZAMIENTOS

CALCULO DE TENSIONES

ANÁLISIS E
IMPRESION DE RESULTADOS
POSTPROCESADOR

Figura A: Esquema Proceso de Cálculo mediante ordenador.

2.2.1 Creación de geometría y entrada de datos:

En la entrada de datos se define la geometría de la estructura, las propiedades

del material, las cargas que ha de soportar y los apoyos que la sustentan.

A modo de ejemplo se detalla a continuación la que pudiera ser una secuencia

idónea de introducción de datos a la hora de crear un modelo de elementos finitos.

1 – Crear o importar geometría,

2 – Especificar material y propiedades de cada elemento que
componen el modelo,

5
 3 – Mallado de la geometría mediante elementos y nodos.
4 – Aplicar condiciones de contorno o restricciones de movimientos en
los nodos según correspondan,
5 – Introducir sistemas de cargas equivalentes que simulen la puesta
en servicio de la estructura.

2.2.2 Matriz de rigidez: Obtención de la matriz de rigidez en cada elemento.

En esencia el método de los elementos finitos se resuelve cuando en cada

elemento se plantea la ecuación que liga las deformaciones de los nodos con las
fuerzas en ella aplicadas a través de la matriz de rigidez. El procedimiento que se
sigue consiste en suponer unas leyes, generalmente polinómicas, que definen los
recorridos de los puntos del elemento finito en función de las coordenadas
asignadas a los nodos. A partir de aquí, se deducen para los nodos, aplicando el
Teorema de los Trabajos Virtuales, las fuerzas equivalentes al estado tensional
creado como consecuencia de los recorridos asignados a los puntos del elemento.
Finalmente se relacionan fuerzas en nodos y recorridos de estos, con lo que se
define la matriz de rigidez.

2.2.3 Ensamblaje de la matriz: Obtención de la matriz de rigidez global.

Tras la obtención de las matrices de rigidez de todos los elementos se procede

a su ensamblaje a fin de calcular la matriz de rigidez [K] de la totalidad de la
estructura.

Dado que cada nodo puede tener hasta 6 grados de libertad, la estructura
tendrá un número total de grados de libertad igual a 6 veces el número total de
nodos del modelo. Con todo esto cabe esperar que la matriz de rigidez [K] ha de ser
una matriz cuadrada de orden 6 veces el número de nodos que componen la
estructura.

2.2.4 Obtención de los desplazamientos:

Consiste en la resolución de la ecuación {F} = [K] · {δ} empleándose para ello

como procedimientos mas adecuados el método de eliminación de Gauss, el de
Cholesky y el iterativo de Ga uss – Seidel. Esta es la fase que ocupa la mayor parte
del tiempo de cálculo.

2.2.5 Obtención de tensiones:

6
 A partir de la deducción de los desplazamientos en los nodos se obtienen las
tensiones o fuerzas de sección para cada elemento finito.

2.2.6 Análisis y presentación de resultados:

Una vez resueltas todas las ecuaciones y calculadas todas las incógnitas, los
resultados se presentan tanto en listados como de forma visual en pantalla
empleando colores para señalar zonas críticas y puntos donde los esfuerzos son
mayores.

2.3 Tipos de Elementos:

Los principales tipos de elementos que podemos usar para crear modelos se
enumeran a continuación:

Elemento Plano (MEMBRANE):

- 2 grados de libertad cada nodo,
- Aplicación: estructuras planas,
- Ejemplos: triángulo de 3 nodos, cuadrilátero de 4 nodos.

Figura B: Ejemplo de elementos con forma triangular y cuadrada.

Elemento Placa (PLATE BENDING):

- 3 grados de libertad cada nodo (1 desplazamiento + 2 giros),
- Aplicación: chapas donde las posibles fuerzas actúan siempre
perpendicularmente al plano.

Elemento Cascara (SHELL):

- 5 grados de libertad cada nodo (3 desplazamientos + 2 giros),

7
 - Aplicación: todo tipo de estructuras laminares, espesor< 0,1 su
longitud menor,
- Precaución: singularidad de 6 grados de libertad en nodos con
todos sus elementos en el mismo plano.

Elemento Sólido (SOLID):

- 3 grados de libertad cada nodo (3 desplazamientos),
- Aplicación: estructuras de paredes gruesas,
- Ejemplos: tetraedro, cuñas, hexaedro.

Elementos tipo barra:

- ROD: 2 grados de libertad cada nodo. Elemento que trabaja tracción
– compresión + torsión.
- BAR: 6 grados de libertad cada nodo (3 desplazamientos + 3 giros).
- BEAM: Igual que BAR pero con posibilidad de sección variable.

Otros tipos de elementos:

- SPRING: Elemento que permite rigidez y amortiguamiento.
- GAP: permite definir dos rigideces distintas según trabaja a
compresión o tracción.
- MASS: Define una masa concentrada en un nodo,
- RIGID: Elemento totalmente rígido,
- SLIDE LINE: Elementos para definir contactos.

3. EJEMPLOS DE APLICACIÓN.

3.1 Presentación de modelos:

Se presentan cuatro tipos de columnas soporte de proyectores o luminarias.

- TIPO 1: Columna de diseño especial para soportar dos luminarias.

Su característica fundamental es la de tener sección constante a lo
largo de toda su longitud. Altura útil 12 metros.
- TIPO 2: Columna inclinada para soportar proyectores que
iluminarán instalaciones deportivas.
Características principales: Columna formada por dos mástiles
inclinados unidos entre sí mediante tubos.

8
 - TIPO 3: Columna CA-14 de 35 metros de altura de dimensiones
normalizadas para alumbrar grandes áreas. Es de utilidad para
estudiar partes concretas y mejorar el funcionamiento del conjunto
estructural.
- TIPO 4: Columna invertida para alumbrado. Su principal
característica es que presenta la sección de menor tamaño en su
base, teniendo su mayor sección en la punta.

3.2 Presentación de resultados:

A través de las figuras que forman parte del presente informe se muestran los
pasos seguidos para generar el modelo y resultados obtenidos a los cálculos
realizados:

• Figura-1: Geometría de la columna de sección cuadrada formada a partir

de puntos, líneas y superficies. Es el punto de partida del estudio.
• Figura-2: Detalle geometría zona inferior de columna TIPO 1.
• Figura-3: Mallado columna TIPO 1. Se realiza en placa base un mallado
formado por elementos triangulares tipo placa, el resto de la columna con
elementos cuadrados de igual dimensión. Las zonas de unión se realizan
con mallado manual mediante elementos triangulares. Se observa
mediante diferentes colores los espesores de los materiales que forman
parte del modelo.
• Figura-4: Condiciones de contorno en columna TIPO 1. Se anula giro y
desplazamiento en zona de taladros simulando la acción de los pernos de
anclaje.
• Figura-5: Aplicación de cargas en columna TIPO 1. Para este caso, y
debido a la geometría del modelo, se aplican las cargas en dos direcciones
ortogonales. Con los resultados se decide cual de las dos opciones de
aplicación de carga condiciona el dimensionamiento.
• Figura-6: Estado de tensiones Von Misses (MPa) en columna TIPO 1.
Mediante colores se refleja los puntos y zonas críticas, que en el caso
particular se localizan en unión de fuste con placa de anclaje y
alrededores.
• Figura-7: Geometría Pórtico. Geometría compuesta por dos fustes de
sección cuadrada variable a lo largo de su longitud, unidos entre si
mediante tubo redondo estructural.

9
• Figura-8: Detalle mallado zona inferior PORTICO. El mallado de la placa
base se realiza automático con elementos tipo placa de forma triangular.
• Figura-9: Mallado zona superior PORTICO. Mediante colores se
representan los distintos espesores que forman parte de la estructura.
• Figura-10: Mallado zona superior PORTICO.
• Figura-11: Aplicación de cargas en PORTICO. Las cargas que se aplican
son viento sobre fuste, viento sobre proyectores y peso de proyectores. El
peso propio de la estructura lo calcula internamente el programa.
• Figura-12: Estado de tensiones Von Misses (MPa) en zona inferior de
fuste.
• Figura-13: Estado de tensiones Von Misses (MPa) en placa de anclaje de
fuste de PORTICO
• Figura-14: Deformada de PORTICO según cargas aplicadas. La figura
muestra como se deforma (de forma exagerada) el modelo cuando se le
aplican las cargas.
• Figura-15: Geometría de columna tipo CA-14 en parte inferior donde se
localiza el registro.
• Figura-16: Detalle mallado de placa de anclaje, zona inferior de fuste y
registro para columna CA-14.
• Figura-17: Esquema de aplicación de cargas. Las cargas se aplican en
dirección perpendicular al plano que contiene el registro inferior por ser
ésta la orientación más débil de la columna.
• Figura-18: Estado de tensiones Von Misses (MPa) para estado de cargas
descrito en zona inferior de fuste a barlovento.
• Figura-19: Estado de tensiones Von Misses (MPa) para estado de cargas
descrito en zona inferior de fuste a sotavento.
• Figura-20: Estado de tensiones Von Misses (MPa) en placa de anclaje de
columna CA-14.
• Figura-21: Vista de estado de tensiones para estado de carga aplicado.
Vista general en empotramiento de tramos de la columna.
• Figura-22: Geometría de columna INVERTIDA.
• Figura-23: Detalle geometría zona inferior columna INVERTIDA.
• Figura-24: Mallado general de columna tipo INVERTIDA. El mallado del
fuste se realiza de forma automática mediante elementos tipo placa con
forma cuadrada de igual dimensión.
• Figura-25: Detalle mallado zona inferior columna INVERTIDA. El mallado
en la placa base se realiza mediante elementos con forma triangular.
• Figura-26: Esquema de aplicación de cargas equivalente.

10
• Figura-27: Tensiones Von Misses (MPa) en zona inferior de fuste para
cargas aplicadas.
• Figura-28: Tensiones Von Misses (MPa) en placa de anclaje para cargas
aplicadas.

11
Figura 1: Geometría columna sección cuadrada.
13
Figura 2: Detalle geometría zona inferior en columna sección cuadrada.
14
Figura 3: Detalle mallado zona inferior columna sección cuadrada.
15
Figura 4: Aplicación de condiciones de conto rno en columna sección cuadrada.
16
Figura 5: Aplicación de cargas en columna de sección cuadrada.
17
Figura 6: Estado de tensiones en columna de sección cuadrada.
18
Figura 7: Aspecto general geometría pórtico.
19
Figura 8: Detalle mallado zona inferior de pórtico.
20
Figura 9: Detalle mallado zona superior de pórtico.
21
Figura 10: Detalle mallado superior de pórtico.
22
Figura 11: Aplicación de cargas sobre pórtico.
23
Figura 12: Estado de tensiones en zona inferior de pórtico.
24
Figura 13: Detalle estado de tensiones en placa de anclaje de pórtico.
25
Figura 14: Deformada de pórtico según cargas aplicadas.
26
Figura 15: Detalle geometría zona inferior de columna tipo CA-14.
27
Figura 16: Detalle mallado zona inferior de columna tipo CA-14.
28
Figura 17: Aplicación de cargas en modelo de columna tipo CA-14.
29
Figura 18: Estado de tensiones en parte inferior de columna tipo CA-14.
30
Figura 19: Estado de tensiones en columna tipo CA-14.
31
Figura 20: Estado de tensiones en placa de anclaje de columna tipo CA-14.
32
Figura 21: Estado de tensiones en columna tipo CA-14.
33
Figura 22: Geometría de modelo tipo columna Invertida.
34
Figura 23: Detalle geometría placa de anclaje y registro en columna Invertida.
35
Figura 24: Aspecto general mallado columna Invertida.
36
Figura 25: Esquema aplicación cargas equivalentes en columna Invertida.
37
Figura 26: Estado de tensiones en zona inferior de columna Invertida.
38
Figura 27: Estado de tensiones en placa de anclaje de columna invertida.
39