UNIVERSIDAD NACIONAL
DE INGENIERA
Facultad de Ingeniera
Mecnica
CALCULO POR
ELEMENTOS FINITOS
2017-2
LABORATORIO N5
FLEXION DE VIGAS
Ing. Edwin Ascencin Abreg Leandro
Alumno: Jerson Arnold Suarez Becerra 20144019E
1 Calculo por Elementos Finitos
NDICE
Enunciado del Problema 3
Solucin 3
Flujograma 7
Cdigo de Matlab 8
Bibliografa 13
Anexos 13
2 Calculo por Elementos Finitos
ENUNCIADO DEL PROBLEMA
La viga mostrada en la figura es una brida ancha W310x52 con un rea de
seccin transversal de 6650 mm2 y una profundidad de 317 mm. El segundo
momento de la zona es 118.6x106mm4. El haz se somete a una carga
uniformemente distribuida de 25000 N / m. El mdulo de elasticidad de la viga
es E = 200GPa. Determine el desplazamiento vertical de un nodo 3 y las
rotaciones en los nodos 2 y 3.
Solucin
SOLUCION
1. MODELADO Y GDL
QJ Q1 Q2 Q3 Q4 Q5 Q6
Desplazamientos y torsiones de los nodos 1, 2 y 3
2. CARGA APLICADA (NODAL EQUIVALENTE)
P=25000N/m
p * le 25000 * 7.5
F1' F3' F5' 62500 N
3 3
p * (le) 2 25000 * (7.5)
F4' F6' 117187.5 N m
12 12
3 Calculo por Elementos Finitos
F1' 62500N F3' 62500 N F5' 62500 N
F4' 117187.5N m F6' 117187.5 N m
3. MATRICES DE RIGIDEZ (LOCALES)
El factor numrico:
E * I 200 *10 9 *118.6 *10 6
1.97666 *10 6 N / m
le 3
12
12 6le 12 6le 12 6le 12 6le
6le 4le 2 6le 2le 2 6le 4le 2 6le 2le 2
Kij1 1.97666 *106 * K 2 1.97666 *106 *
12 6le 12 6le ij 12 6le 12 6le
6le 2le
2
6le 4le 2 6le 2le
2
6le 4le 2
4. MATRIZ DE RIGIDEZ (GLOBAL)
12 6 12 6 0 0
6 4 6 2 0 0
12 6 24 0 12 6
K ij 1.97666 *10 6 *
6 2 0 8 6 2
0 0 12 6 12 6
0 0 6 2 6 4
5. VECTOR CARGA
F1' 62500N F3' 62500 N F5' 62500 N
F4' 117187.5N m F6' 117187.5 N m
R5 0
R2 R4 0
R6 0
R1 R3
4 Calculo por Elementos Finitos
Resulta:
Fi R1 62500 R2 R3 62500 117187.5 625000 117187.5
6. CONDICIONES DE CONTORNO (GDL)
Q j 0 0 0 Q4 Q5 0
Variables de contorno que necesitan ser hallados
6. ECUACION DE RIGIDEZ
Fi K ij * QJ
7.1 RESOLVIENDO:
117187.5 8 6 Q4
625000 1.97666 *10 * 6 12 * Q
6
5
Q4 0.0434rad
Q5 0.04808mm
7.2 RELACION EN LOS APOYOS
El vector carga se puede escribir:
Fi Ri Fi '
En la ecuacin de rigidez
Ri Kij * Q j Fi '
8. FLECHA
V e H r' * qr
1
l
l
V e * 2 3 3 q1 e 1 2 3 q2 2 3 3 q3 e 1 2 3 q4
4 2 2
Hallando la deformacin en el punto final (3)
Flecha en el punto medio de (3) en
0
1
V 2 * 2q1 e 1q2 2q3 e 1q4
l l
4 2 2
Dnde:
5 Calculo por Elementos Finitos
q1 Q3 0
q2 Q4 0.0434rad
q3 Q5 0.04808mm
q4 Q6 0
Resolviendo:
V 2 0.926mm
9. PENDIENTE
dV e 3 3
l
3
l
1 q1 e 1 2 3 2 q2 2 1 q3 e 1 2 3 2 q4
d 4
8 4 8
Flecha en el punto medio de (2) en
0
dV e
0.00694
d
10. ESFUERZO
E* y
xe 2 6q1 (3 1)l2 q2 6q3 3 1le q4
le
200 *109 * 0.04808
1
x 2
(1)5(0) 0 1(5)( 0.0434)
7.5
x 91.5065Gpa
1
200 *109 *118.6 6(0.04808)(0) (3(0.04808) 1)( 2.5)( 0.0434)
x2
7.52 6(0.04808)( 0.04808) 3(0.04808) 1(2.5)(0)
x2 45.0026Gpa
11. RESULTADOS
DESPLAZAMIENTO Y TORSION ESFUERZO
Q1 0 1x 91.5065Gpa
Q2 0 x2 45.0026Gpa
Q3 0
Q4 0.0434rad
Q5 0.04808mm
Q6 0
DIAGRAMA DE FLUJO
6 Calculo por Elementos Finitos
INICIO
Leer
datos
de
Para i=1:4
Calcula la matriz de
rigidez de cada elemento
y tambin la global.
Calcula
desplazamientos,
reacciones
Para i=1:4
Calcula
esfuerzos para
e=-1,1
Si
ES1<=ES2
Emax=E Emax=ES
S2 1
Imprime esfuerzos y
reacciones.
CDIGO DE MATLAB
7 Calculo por Elementos Finitos
clc;
format long;
n=input('Ingrese Numero de Elementos Finitos:');
e1=input('Espesor de las alas(mm):');
e2=input('Espesor del alma(mm):');
l1=input('Longitud de las alas(mm):');
L=input('Ingrese Longitud de la Viga(mm):');
E=input('Modulo de Elasticidad(N/mm2):');
yp=input('Ingrese Peso Especifico(N/mm3):');
pe=input('Carga Distribuida Externa(N/mm):');
disp('MOMENTOS DE INERCIA')
for i=1:(n/2)
d(i)=(4*(L*(i-1)/n)/15+100+4*(L*i/n)/15+100)/2;
I(i)=2*(l1*e1^3/12)+e2*(d(i)-2*e1)^3/12+2*(d(i)/2-e1/2)^2*l1*e1;
end
for i=((n/2)+1):n
d(i)=(900-4*(L*(i-1)/n)/15+900-4*(L*i/n)/15)/2;
I(i)=2*(l1*e1^3/12)+e2*(d(i)-2*e1)^3/12+2*(d(i)/2-e1/2)^2*l1*e1;
end
disp(I)
disp('MATRIZ DE RIGIDEZ ESTRUCTURAL K')
k=zeros(2*(n+1),2*(n+1));
for i=1:n
l=L/n;
ke(:,:,i)=E*I(i)/l*[12 6*l -12 6*l;
6*l 4*l*l -6*l 2*l*l;
-12 -6*l 12 -6*l;
6*l 2*l*l -6*l 4*l*l];
gl1=2*i-1;
gl2=gl1+1;
gl3=2*(i+1)-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
k(gl,gl)=k(gl,gl)+ke(:,:,i);
end
disp(k)
disp('FUERZAS DEBIDO AL PESO DEL MATERIAL')
for i=1:n
A(i)=l1*e1*2+e2*(d(i)-2*e1);
p(i)=-yp*A(i);
end
w=zeros(1,2*(n+1));
for i=1:n
l=L/n;
we(:,:,i)=[p(i)*l/2 p(i)*l^2/12 p(i)*l/2 -p(i)*l^2/12];
gl1=2*i-1;
gl2=gl1+1;
gl3=2*(i+1)-1;
gl4=gl3+1;
8 Calculo por Elementos Finitos
gl=[gl1 gl2 gl3 gl4];
w(1,gl)=w(1,gl)+we(:,:,i);
end
wt=w';
disp(wt)
disp('FUERZAS DEBIDO A LA CARGA DISTRIBUIDA')
c=zeros(1,2*(n+1));
for i=2:3
l=L/n;
ce(:,:,i)=[pe*l/2 pe*l^2/12 pe*l/2 -pe*l^2/12];
gl1=2*i-1;
gl2=gl1+1;
gl3=2*(i+1)-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
c(1,gl)=c(1,gl)+ce(:,:,i);
end
ct=c';
disp(ct)
disp('FUERZA TOTAL')
f=ct+wt;
disp(f)
disp('DESPLAZAMIENTOS')
disp('Q=')
kf=k(3:8,3:8);
ff=f(3:8,1);
qf=inv(kf)*ff;
Q=[0;0;qf;0;0];
disp(Q)
disp('ESFUERZOS LONGITUDINALES(N/mm2)')
y=input('Ingrese punto generico a analizar:');
z=-1;
es1=zeros(n,1);
for i=1:n
gl1=i*2-1;
gl2=gl1+1;
gl3=(i+1)*2-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
q=Q(gl);
es1(i)=-E*y/l^2*[6*z (3*z-1)*l -6*z (3*z+1)*l]*q;
end
disp('z=-1')
disp(es1)
z=1;
es2=zeros(n,1);
for i=1:n
gl1=i*2-1;
gl2=gl1+1;
9 Calculo por Elementos Finitos
gl3=(i+1)*2-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
q=Q(gl);
es2(i)=-E*y/l^2*[6*z (3*z-1)*l -6*z (3*z+1)*l]*q;
end
disp('z=1')
disp(es2)
z=0;
es0=zeros(n,1);
for i=1:n
gl1=i*2-1;
gl2=gl1+1;
gl3=(i+1)*2-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
q=Q(gl);
es0(i)=-E*y/l^2*[6*z (3*z-1)*l -6*z (3*z+1)*l]*q;
end
disp('z=0')
disp(es0)
EJECUCIN DEL PROGRAMA
Ingrese Nmero de Elementos Finitos: 4
Espesor de las alas (mm):13
Espesor del alma (mm):25
Longitud de las alas (mm):100
Ingrese Longitud de la Viga (mm):3000
Modulo de Elasticidad(N/mm2):2.1e5
Ingrese Peso Especifico(N/mm3):76.518e-6
Carga Distribuida Externa(N/mm):-5
MOMENTOS DE INERCIA
1.0e+008 *
0.33741516666667 2.06373183333333 2.06373183333333
0.33741516666667
MATRIZ DE RIGIDEZ ESTRUCTURAL K
1.0e+017 *
Columns 1 through 4
0.00000113371496 0.00042514311000 -0.00000113371496
0.00042514311000
10 Calculo por Elementos Finitos
0.00042514311000 0.21257155500000 -0.00042514311000
0.10628577750000
-0.00000113371496 -0.00042514311000 0.00000806785392
0.00217515900000
0.00042514311000 0.10628577750000 0.00217515900000
1.51272261000000
0 0 -0.00000693413896 -0.00260030211000
0 0 0.00260030211000 0.65007552750000
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
Columns 5 through 8
0 0 0 0
0 0 0 0
-0.00000693413896 0.00260030211000 0 0
-0.00260030211000 0.65007552750000 0 0
0.00001386827792 0 -0.00000693413896 0.00260030211000
0 2.60030211000000 -0.00260030211000 0.65007552750000
-0.00000693413896 -0.00260030211000 0.00000806785392 -
0.00217515900000
0.00260030211000 0.65007552750000 -0.00217515900000
1.51272261000000
0 0 -0.00000113371496 -0.00042514311000
0 0 0.00042514311000 0.10628577750000
Columns 9 through 10
0 0
0 0
0 0
0 0
0 0
0 0
-0.00000113371496 0.00042514311000
-0.00042514311000 0.10628577750000
0.00000113371496 -0.00042514311000
-0.00042514311000 0.21257155500000
FUERZAS DEBIDO AL PESO DEL MATERIAL
1.0e+004 *
-0.01994250375000
-2.49281296875000
-0.05423213250000
-1.79339062500000
11 Calculo por Elementos Finitos
-0.06857925750000
0
-0.05423213250000
1.79339062500000
-0.01994250375000
2.49281296875000
FUERZAS DEBIDO A LA CARGA DISTRIBUIDA
0
0
-1875
-234375
-3750
0
-1875
234375
0
0
FUERZA TOTAL
1.0e+005 *
-0.00199425037500
-0.24928129687500
-0.02417321325000
-2.52308906250000
-0.04435792575000
0
-0.02417321325000
2.52308906250000
-0.00199425037500
0.24928129687500
DESPLAZAMIENTOS
Q=
1.0e-006 *
0
0
-0.06887011717571
-0.00007462641465
-0.10005353989168
-0.00000000000000
-0.06887011717571
0.00007462641465
0
0
ESFUERZOS LONGITUDINALES(N/mm2)
Ingrese punto generico a analizar: 50
12 Calculo por Elementos Finitos
z=-1
1.0e-005 *
0.56239135133738
-0.06865358764229
-0.14030037338828
-0.35343739030681
z=1
1.0e-005 *
-0.35343739030681
-0.14030037338828
-0.06865358764229
0.56239135133738
z=0
1.0e-005 *
0.10447698051528
-0.10447698051528
-0.10447698051528
0.10447698051528
BIBLIOGRAFA
[1] Introduccin al estudio del elemento finito en ingeniera. Segunda
Edicin. Tirupathi. R. Chandrupatla y Ashok D. Belegundu.
[2] Finite Element Analysis. Tercera Edicin. Saeed Moaveni.
[3] Apuntes de clase del Ing. Cueva y de Ing. Abreg.
ANEXOS
a. Simulacin del problema con el programa ANSYS 18.1
13 Calculo por Elementos Finitos
