Análisis por Elementos Finitos
Clase 5
Ing. Eduardo Orcés P.
Junio 15/2011
Ing. Eduardo Orcés
2011-I
TEMAS
Objetivos de aprendizaje
Análisis de Vigas
Programación en MATLAB / CALFEM
Ing. Eduardo Orcés 2011-I
Vigas
Son elementos estructurales esbeltos que
soportan cargas transversales y momentos.
Sus esfuerzos internos son de flexión.
Ing Eduardo Orcés 2011-I
Vigas
Objetivos de aprendizaje:
o Entender el comportamiento estructural de las vigas.
o Estar en capacidad de crear modelos estructurales
apropiados para vigas sometidas a flexión.
o Aplicar el método de EF para analizar vigas usando
Matlab/Calfem.
Ing Eduardo Orcés 2011-I
Vigas
• Ecuaciones de Mecánica de Sólidos:
Ing Eduardo Orcés 2011-I
Vigas
• Modelo estructural de una viga:
Ing Eduardo Orcés 2011-I
• Asumiendo que la carga se aplica solamente en
los nodos:
• Se escoge un polinomio de interpolación cúbico,
el cual es la solución exacta de la ecuación
diferencial anterior:
Ing Eduardo Orcés 2011-I
• Expresando las ecuaciones anteriores en forma matricial:
• El desplazamiento v(x) dentro del elemento está dado por
la siguiente expresión:
Ing Eduardo Orcés 2011-I
• Los Ni(x) son los polinomios osculatorios de Hermite.
• Notar que,
• Lo cual es requerido para que el movimiento de cuerpo
rígido sea representado por el modo de deformación
asumido.
Ing Eduardo Orcés 2011-I
• Fuerzas Cortantes y Momentos Flectores en los extremos
del elemento viga:
Ing Eduardo Orcés 2011-I
• Expresando en forma matricial los resultados anteriores
obtenemos la matriz de rigidez del elemento viga.
Ing Eduardo Orcés 2011-I
• Interpretación de los términos de la matriz de rigidez:
El término kij se puede interpretar como la fuerza (o
momento) requerida en la dirección i para producir un
desplazamiento ( o rotación) unitario en la dirección j.
Ing Eduardo Orcés 2011-I
Método de la Energía.
• Partimos de las funciones de forma:
• Para calcular la deformación en el elemento:
Ing Eduardo Orcés 2011-I
• Curvatura de la viga:
• Donde B es la matriz de gradientes, dada por:
Ing Eduardo Orcés 2011-I
• Calculamos la Energía de deformación elástica:
Ing Eduardo Orcés 2011-I
• El término entre paréntesis da una expresión
general para la matriz de rigidez del elemento
viga:
Ing Eduardo Orcés 2011-I
• Realizando las operaciones matriciales indicadas
obtenemos la matriz k(e) del elemento viga.
Ing Eduardo Orcés 2011-I
Ejemplo: Usando el Método Directo de la Rigidez,
determinar desplazamientos y fuerzas nodales,
así como los Diagramas de MF y FC.
• Solución: Matriz de rigidez de la estructura:
Ing Eduardo Orcés 2011-I
C.B.:
• Eliminando 3ª, 5ª y 6ª filas y columnas:
Ing Eduardo Orcés 2011-I
• Resolviendo este sistema de ecuaciones, se obtienen los
desplazamientos nodales („GDL activos‟).
• Substituyendo en F = KQ, se obtienen las fuerzas
nodales globales (externas):
Ing Eduardo Orcés 2011-I
• Fuerzas y Momentos nodales elementales:
Ing Eduardo Orcés 2011-I
Cargas Distribuidas
• Hallamos las cargas nodales equivalentes
aplicando el principio de trabajo virtual.
L L L
p( x)v( x)dx p( x)N(xv) qdx p( x)N(x)dx q
0 0 0
L
T T
P q
E PE p( x)N(x)dx
0
Ing Eduardo Orcés 2011-I
• Por ejemplo, para p(x) = p = constante:
1
PET p 2 3
3 2
1, L 3
2 2
, -2 3
3 2, L 3 2
Ld
0
pl pl 2 pl pl 2 x
, , , -
2 12 2 12 L
Ing Eduardo Orcés 2011-I
Ejercicio 5.1: Hallar los desplazamientos y
fuerzas nodales en la viga en voladizo sometida
a la carga uniforme mostrada. También se
muestran las cargas nodales equivalentes.
Ing Eduardo Orcés 2011-I
Ejemplo: (a) Hallar desplazamientos y fuerzas nodales
globales.
(b) Hallar fuerzas nodales elementales.
(c) Determinar la deflexión en el punto medio del
elemento 2.
Ing Eduardo Orcés 2011-I
• Solución: E = 200 GPa, I = 4 x 106 mm4
• Cargas nodales equivalentes:
• Rigideces elementales:
Ing Eduardo Orcés 2011-I
• Ecuaciones de equilibrio global:
• C.B.:
• GDL activos:
Ing Eduardo Orcés 2011-I
• Fuerzas nodales globales (“reacciones”): Substituyendo
(B) en (A).
• Fuerzas nodales elementales:
Ing Eduardo Orcés 2011-I
• Usando conectividad para relacionar los GDL globales
y locales:
Ing Eduardo Orcés 2011-I
• Diagramas de Fuerza Cortante y Momento Flector:
Ing Eduardo Orcés 2011-I
• Desplazamiento en el punto medio del elemento 2:
Sumamos el desplazamiento causado por las „cargas
nodales equivalentes‟ y por las „cargas de empotramiento
perfecto‟:
q1 0
q2 Q4 - wL4
v x/L 0.5
N q x/L 0.5
vEP N1 N 2 N 3 N 4
q3 0 384EI
q4 Q6 x / L 0. 5
8.93 10 5 m 3.906 10-5 m
0.0128 mm
Ing Eduardo Orcés 2011-I
Ejercicio 5.2: Hallar las deflexiones y reacciones de la
viga mostrada, cuyo extremo derecho está apoyado en
un resorte (apoyo flexible). Datos: P = 50 kN, k = 200
kN/m, L = 3 m, E = 210 GPa, I = 2 x 10-4 m4.
Ing Eduardo Orcés 2011-I
• Solución:
• Usamos 2 elementos viga y 1 elemento resorte.
• La matriz de rigidez del elemento resorte es:
• La matriz de rigidez global de la estructura es:
Ing Eduardo Orcés 2011-I
• Para simplificar la notación, se ha usado:
• Condiciones de borde y carga:
• Obtenemos los GDL activos eliminando las filas y columnas
1ª, 2ª, 3ª y 7ª:
Ing Eduardo Orcés 2011-I
• Resolviendo este sistema de ecuaciones, obtenemos:
• Resolviendo las ecuaciones globales de equilibrio, obtenemos
las fuerzas nodales de reacción:
• Chequeo:
Ing Eduardo Orcés 2011-I
Programa en Matlab / Calfem
% Progama 'Viga2'
% Analisis de una viga de 2 tramos, con carga distribuida
% sobre uno de ellos.
%
%----- Conectividad ---------------------------------------------
conec1=[1 1 2 3 4 5 6];
conec2=[2 4 5 6 7 8 9];
%----- Matriz de rigidez K, vector de cargas P -------------------
K=zeros(9); P=zeros(9,1);
%----- Matrices de rigidez de los elementos ---------------------
E=2.0e11; A=1.0e-4; I=4.0e-6; ep=[E A I]; % N, m
ex1=[0 1]; ey1=[0 0]; ex2=[1 2]; ey2=[0 0];
eq=[0 -12e3];
Ke1=beam2e(ex1,ey1,ep)
[Ke2,Pe2]=beam2e(ex2,ey2,ep,eq)
%----- Ensamblaje de K(e) en K -----------------------------------
K=assem(conec1,K,Ke1);
[K,P]=assem(conec2,K,Ke2,P,Pe2);
Ing Eduardo Orcés 2011-I
Programa en Matlab / Calfem
%----- Resolucion del sistema de ecuaciones --------------------
cb=[1 0; 2 0; 3 0; 5 0; 8 0];
[D,R]=solveq(K,P,cb);
%----- Fuerzas en los elementos -------------------------------
ed1=extract(conec1,D);
ed2=extract(conec2,D);
[es1,edi1,eci1]=beam2s(ex1,ey1,ep,ed1,[0 0],20);
[es2,edi2,eci2]=beam2s(ex2,ey2,ep,ed2,eq,20);
%----- Resultados ----------------------------------------------------
Desplazamientos= D
Reacciones= R
es1;
es2;
Ing Eduardo Orcés 2011-I
Programa en Matlab / Calfem
%----- Resolucion del sistema de ecuaciones --------------------
cb=[1 0; 2 0; 3 0; 5 0; 8 0];
[D,R]=solveq(K,P,cb);
%----- Fuerzas en los elementos -------------------------------
ed1=extract(conec1,D);
ed2=extract(conec2,D);
[es1,edi1,eci1]=beam2s(ex1,ey1,ep,ed1,[0 0],20);
[es2,edi2,eci2]=beam2s(ex2,ey2,ep,ed2,eq,20);
%----- Resultados ----------------------------------------------------
Desplazamientos= D
Reacciones= R
es1;
es2;
Ing Eduardo Orcés 2011-I
Programa en Matlab / Calfem
%----- Diagramas de Desplazamientos / FC /MF ------
figure(1)
plot(eci1,edi1(:,2),'-r',1+eci2,edi2(:,2),'-r')
title('Desplazamientos'),xlabel('m'),ylabel('m')
grid;
figure(2)
plot(eci1,-es1(:,2),'-r',1+eci2,-es2(:,2),'-r')
title('Fuerzas Cortantes'),xlabel('m'),ylabel('N')
grid;
figure(3)
plot(eci1,es1(:,3),'-r',1+eci2,es2(:,3),'-r')
title('Momentos Flectores'),xlabel('m'),ylabel('N.m')
grid;
Ing Eduardo Orcés 2011-I
Tareas
Leerlas siguientes secciones del texto:
- Cap. 8, Vigas y Marcos, Secs. 8.1 a 8.5
Resolver problemas 1, 2 y 3 del Deber #3
Resolver ejercicios 5.1 y 5.2 de la clase 5.
Ing. Eduardo Orcés 2011-I
