Elementos Finitos Elementos Resorte Barra y Viga
Elementos Finitos Elementos Resorte Barra y Viga
Elementos Finitos Elementos Resorte Barra y Viga
ESTRUCTURAS III
Autores:
-2007-
Estructuras III Introducción a los Elementos Finitos
1-Introducción
Si bien fue originalmente desarrollado para el análisis de estructuras, con este método se
pueden representar entre otros, los siguientes fenómenos físicos:
2- Concepto
[K] .{δ}={F}
Donde :
Página 1/26
Estructuras III Introducción a los Elementos Finitos
Página 2/26
Estructuras III Introducción a los Elementos Finitos
Modelado
Geométrico
Modelado
de
Elementos Finitos
Definición del
Ambiente
Análisis
Corroboración de Resultados
• Modelado Geométrico: Creación del modelo matemático del objeto o del conjunto.
Reproducción del sólido en forma precisa y de la geometría de la superficie.
• Definición del Ambiente: Aplicar las cargas y las condiciones de borde para simular
el ambiente de la operación.
Página 3/26
Estructuras III Introducción a los Elementos Finitos
4- Método Directo
Página 4/26
Estructuras III Introducción a los Elementos Finitos
4-1-Elemento Resorte
• Dos Nodos: i, j
• Constante Elástica del Resorte: k
El elemento tiene dos grados de libertad, en el sentido longitudinal del elemento, cualquier
desplazamiento de los nodos en el sentido normal al elemento no generara esfuerzos internos:
Nodo i:
Nodo j:
Expresado matricialmente:
O bien, k.u =f
Donde:
• k: es la matriz rigidez
• u: es la vector desplazamiento
• f: es la vector de fuerzas internas
Página 5/26
Estructuras III Introducción a los Elementos Finitos
4-2-Sistema de Resortes
Para el Elemento 1:
Para el Elemento 2:
Nodo 1:
Nodo 2:
Nodo 3:
Desarrollando con los valores de las fuerzas internas en función de la rigidez de cada
elemento:
Página 6/26
Estructuras III Introducción a los Elementos Finitos
De forma Matricial:
O bien, K.U=F
Reemplazando:
Se reduce a:
Página 7/26
Estructuras III Introducción a los Elementos Finitos
Y,
Resolviendo,
Como conclusión, para un sistema de “n” nodos, el método de elementos finitos permite
generar “n” ecuaciones, las cuales deberán tener “n” incógnitas para ser un sistema definido.
Las incógnitas podrán ser parte del vector desplazamiento o ser parte del vector fuerza. Cada
nodo deberá tener su desplazamiento o su fuerza actuante como condición de borde impuesta.
Este sistema permite, cómo veremos más adelante resolver sistemas isoestáticos e hiperestáticos
sin necesidad de cambiar el método.
Página 8/26
Estructuras III Introducción a los Elementos Finitos
• Dos Nodos: i, j
• Modulo de Elasticidad E
• Área de la Sección Transversal A
• Longitud del Elemento L
El elemento tiene dos grados de libertad, en el sentido longitudinal del elemento, cualquier
desplazamiento de los nodos en el sentido normal al elemento no generara esfuerzos internos:
Por lo tanto,
O bien,
Página 9/26
Estructuras III Introducción a los Elementos Finitos
Para la resolución de este sistema se procede de la misma manera que en el elemento resorte.
Local Global
x,y X,Y
ui´,vi´ ui,vi
1 grado de Libertad 2 grados de Libertad
Donde:
Página 10/26
Estructuras III Introducción a los Elementos Finitos
O bien,
A modo de ejemplo, se puede decir que para un sistema de elemento barra con dos nodos
tenemos que:
O bien,
Con,
Página 11/26
Estructuras III Introducción a los Elementos Finitos
O bien,
Obtenemos:
La cual es una matriz simétrica de 4X4. Escrita de manera explicita tenemos que:
4-5-Elemento Viga
Página 12/26
Estructuras III Introducción a los Elementos Finitos
Se considera una viga en el plano. Esta toma esfuerzos de Corte, Axiles y Momentos,
todas consideradas en el plano. Cada Nodo posee tres Grados de Libertad (u, v, q).
Un elemento que toma estas cargas, tiene asociado para el calculo a E, J, l y A.
• Dos Nodos: i, j
• Modulo de Elasticidad E
• Área de la Sección Transversal A
• Longitud del Elemento L
• Momento de Inercia I
Para crear la Matriz Rigidez se suponen casos con desplazamientos unitarios, que luego
mediante Superposición se ensamblan y dan forma a dicha matriz.
Se adoptan giros en sentido horario y desplazamientos positivos.
¾ Se supone ui=1
uj, vi, vj, θi, θj =0
Página 13/26
Estructuras III Introducción a los Elementos Finitos
Aplicando la Ley de Hooke, tal como se hace con elemento barra, tenemos que:
EA
Hi = ⋅ ui
l
Por lo tanto,
EA
Hi =
l
Realizando un equilibrio de fuerzas,
EA
Hj = −
l
¾ Se supone vi =1
ui, uj, vj, θi, θj = 0
Se puede demostrar calculando por Método de las Fuerzas que para un Desplazamiento
Transversal en el extremo i, los esfuerzos en el sistema son:
12EJ 6 EJ
Vi = Mi =
l3 l2
12EJ 6 EJ
Vj = − Mj =
l3 l2
Hi = 0 Hj =0
¾ Se supone θi=1
Página 14/26
Estructuras III Introducción a los Elementos Finitos
6 EJ 4 EJ
Vi = − 2 Mi =
l l
6 EJ 2 EJ
Vj = Mj =
l2 l
Hi = 0 Hj =0
Procediendo de forma análoga para los desplazamientos del Nodo j , obtendremos los restantes
coeficientes de la Matriz Rigidez del Elemento.
5- Ejemplos Prácticos
Página 15/26
Estructuras III Introducción a los Elementos Finitos
5-1-Elemento Resorte
Datos:
Se pide:
Solución
Página 16/26
Estructuras III Introducción a los Elementos Finitos
o bien,
Página 17/26
Estructuras III Introducción a los Elementos Finitos
donde,
Página 18/26
Estructuras III Introducción a los Elementos Finitos
Se pide:
Solución:
Página 19/26
Estructuras III Introducción a los Elementos Finitos
5-3-Elemento Barra
Sean dos barras de igual Longitud, igual Modulo de Elasticidad y el Área de una
es dos veces la de la otra:
Se pide:
Solución:
Elemento 1
Elemento 2
Página 20/26
Estructuras III Introducción a los Elementos Finitos
y reemplazando, se obtiene:
Por lo tanto,
Sean dos barras idénticas que poseen el mismo Modulo de Elasticidad, la misma
Área Transversal y la misma Longitud:
Página 21/26
Estructuras III Introducción a los Elementos Finitos
Se pide:
Solución:
Estas matrices no pueden ser ensambladas juntas, ya que cada una esta en
diferentes Sistemas Coordinados. Por esta razón se debe trabajar con un Sistema
Coordenado Global.
Trabajaremos con cada elemento por separado:
Elemento 1:
Elemento 2:
Página 22/26
Estructuras III Introducción a los Elementos Finitos
Resolviendo:
5-5-Elemento Viga
Página 23/26
Estructuras III Introducción a los Elementos Finitos
Se pide:
Nota: Para este ejemplo, se tomo como convención que el momento antihorario es
positivo. No solo como condición de borde, sino en el análisis matricial, y es por eso
que la matriz rigidez difiere en signos con la de la explicación de la pagina 15.
Solución:
Página 24/26
Estructuras III Introducción a los Elementos Finitos
Reemplazamos y obtenemos:
Resolviendo tenemos,
Página 25/26
Estructuras III Introducción a los Elementos Finitos
6- Referencias
Página 26/26