UNIVERSIDAD

NACIONAL DE
INGENIERA
Facultad De Ingeniera Mecnica
Laboratorio: N5 (Flexin)
Curso:

MC 516-A Calculo por elementos

finitos

Profesor:

Cueva Pacheco, Ronald

Estudiante: Arroyo Cndor, Jean Marco

Codigo:

20102678K

2013-II

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

INDICE
Enunciado del Problema....................................................................2
Solucin (Modelado de la viga)..........................................................3
Matriz de Rigidez de los Elementos...................................................3
Matriz de Rigidez Local.............................................5
Fuerza Total Sometida a la viga........................................................7
Esfuerzos Longitudinales...................................................................8
Diagrama de Flujo..............................................................................9
Uso de Matlab....................................................................................10
Ejecucin del Programa.....................................................................12
Conclusiones................................................................................... 15

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

QUINTA PRCTICA CALIFICADA

(FLEXIN)

ENUNCIADO DEL PROBLEMA:

Modelar la viga mostrada con 4 elementos finitos (por
lo menos), y calcular en ellos los esfuerzos debido a la flexin de la misma.

100 mm
13 mm

d (variable)

25 mm

A-A

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

Material:
Acero estructural A-36
E=2.1x10^5 N/mm2
=7.8 gr-f/cm3
-Esfuerzos.- En cada elemento finito de la viga; en un punto genrico ( , y):
Ey
e = 2 [ 6 q1 + ( 3 1 )e q 2 6 q3 + ( 3 + 1 )e q 4 ]
e
6 EI
V
e
max
= = 3 [ 2q1 + e q 2 2q3 + e q 4 ]
A
A e
Donde y es la distancia del punto genrico a la fibra neutra.

1. MODELADO DE LA VIGA
Hacemos el modelado de la viga, en 4 elementos finitos:

1
2

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Q1

Calculo por elementos Finitos

Flexion

Q3

1
Q2

Q5

Q4

Q9

Q7

Q10

Q8
Q6

2. MATRIZ DE RIGIDEZ DE LOS ELEMENTOS

100 mm

Para el elemento finito 1:

13 mm

100 x13 3 25 x(200 13 13) 3 100 x133

200 13 2
+
+
+ 2 x(
) x100 x13
12
12
12
2
2
101224550
I1 =
mm4
3
I1 =

25 mm

200 mm

Matriz de Rigidez Local:

k1 =

101224550 12
) 4500
3

12
750

4500

(2.1x10 5 ) x(

4500
2250000
4500
1125000

12
4500
12
4500

4500
1125000

4500

2250000

Para el elemento finito 2:

100 mm
13 mm

25 mm

I2 =

619119550
100 x133 25 x (400 13 13) 3 100 x133
400 13 2
mm4
+
+
+ 2 x(
) x100 x13 I 2 =
3
12
12
12
2
2

Matriz de Rigidez Local:

5

400 mm

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

k2 =

619119550 12
) 4500
3

12
750

4500

( 2.1x10 5 ) x(

Calculo por elementos Finitos

Flexion

4500
2250000
4500
1125000

12
4500
12
4500

4500
1125000

4500

2250000

Para el elemento finito 3:

100 mm
13 mm

25 mm

I3 =

400 mm

619119550
100 x133 25 x (400 13 13) 3 100 x133
400 13 2
mm4
+
+
+ 2 x(
) x100 x13 I 3 =
3
12
12
12
2
2

Matriz de Rigidez Local:

k3 =

619119550 12
) 4500
3

12
750

4500

( 2.1x10 5 ) x (

4500
2250000
4500
1125000

12
4500
12
4500

4500
1125000

4500

2250000

Para el elemento finito 4:

100 mm
13 mm

25 mm

I4 =

101224550
100 x133 25 x (200 13 13) 3 100 x133
200 13 2
mm4
+
+
+ 2 x(
) x100 x13 I 4 =
3
12
12
12
2
2

3. MATRIZ DE RIGIDEZ LOCAL:

200 mm

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

101224550 12
( 2.1x10 5 ) x (
) 4500
3

k4 =
12
750

4500

4500
2250000
4500
1125000

12
4500
12
4500

4500
1125000

4500

2250000

Hallamos las fuerzas que es sometida la viga debido al peso del material:
= 7.8 gr f / cm 3 = 76.518 x10 6 N / mm 3
p1 = p 4 = A1 = A4 = 76.518 x(100 x13 x 2 + 25 x( 200 26)) = 0.5318001 N/mm
p 2 = p 3 = A2 = A3 = 76.518 x (100 x13 x 2 + 25 x (400 26)) = 0.9143901 N/mm
0.5318001x 750
W1 =
2

0.9143901x750
W2 =
2

0.5318001x750 2
12

0.5318001x750
2

0.5318001x750 2

12

0.9143901x 750 2
12

0.9143901x750
2

0.9143901x750 2

12

0.9143901x750
W3 =
2

0.5318001x750
W4 =
2

0.9143901x750 2
12

0.9143901x750
2

0.9143901x750 2

12

0.5318001x750 2
12

0.5318001x750
2

0.5318001x 750 2

12

199.4250375
24928.1296875

542.321325

17933.90625
685.792575
W =

542.321325

17933.90625
199.4250375

24928.1296875

Hallamos las fuerzas debido a la carga distribuida:

5 x750
P2 =
2

5 x750
P3 =
2

5 x750 2
12

5 x750
2

5 x750 2
12

5 x 750 2
12

5 x750
2

5 x750 2

12

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

1875

234375
3750
P =

1875

234375

4. LA FUERZA TOTAL A LA QUE ES SOMETIDA LA VIGA ES:

199.4250375
24928.1296875

2417.321325

252308.90625
4435.792575
F =

2417.321325

252308.90625
199.4250375

24928.1296875

Como los desplazamientos Q1, Q2, Q9 y Q10, quedan restringidos a cero, necesitamos
encontrar Q3, Q4, Q5, Q6, Q7 y Q8.
KQ = F

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

2175159000
8067853.92
2175159000 1512722610000

6934138.96 2600302110

650075527500
2600302110

0
0

0
0

Calculo por elementos Finitos

Flexion

6934138.96

2600302110

2600302110
13868277.92
0
6934138.96
2600302110

650075527500
0
26003021100000
2600302110
650075527500

Q3 2417.321325
Q

4 252308.90625
Q5 4435.792575

=
0

Q6
Q7 2417.321325
252308.90625
Q8
Obtenemos:
0

6.88701171757 x10 8

11
7.46264146538 x10
1.00053539892 x10 7
Q =

6.88701171757 x10 8

11
7.46264146538 x10

5. LOS ESFUERZOS LONGITUDINALES:

Para un punto genrico (z,y), donde z[-1,1]

e = (

Ey
)[ 6 zq1 + (3 z 1)l e q 2 6 zq 3 + (3 z + 1)l e q 4 ]
l e2

Para y=50 mm
Para z=-1

0
0

6934138.96
2600302110

2600302110 650075527500
8067853.92
2175159000

2175159000 1512722610000

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

5.62391351347 x10 6

6.8665358764 x10 7
=
1.40300373388 x10 6

6
3.53437390307 x10

Para z=1
3.53437390307 x10 6

1.40300373388 x10 6
=
6.8665358764 x10 7

6
5.62391351347 x10

Para z=0
1.04476980515 x10 6

1.04476980515 x10 6
=
1.04476980515 x10 6

6
1.04476980515 x10

1
0

Calculo por elementos Finitos

Flexion

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

6. DIAGRAMA DE FLUJO DEL PROGRAMA:

1
1

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

INICI
O
Leer
datos de
entrada
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=ES1

Emax=ES2

Imprime esfuerzos y
reacciones.

1
2

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

7. USANDO MATLAB
PROGRAMA EN MATLAB
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));
1
3

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

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;
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];
1
4

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

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;
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
1
5

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

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
0
0
0
0
0
0

0.00042514311000 -0.00000113371496 0.00042514311000

0.21257155500000 -0.00042514311000 0.10628577750000
-0.00042514311000 0.00000806785392 0.00217515900000
0.10628577750000 0.00217515900000 1.51272261000000
0 -0.00000693413896 -0.00260030211000
0 0.00260030211000 0.65007552750000
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

1
6

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

Columns 9 through 10
0
0
0
0
0
0
-0.00000113371496
-0.00042514311000
0.00000113371496
-0.00042514311000

0
0
0
0
0
0
0.00042514311000
0.10628577750000
-0.00042514311000
0.21257155500000

FUERZAS DEBIDO AL PESO DEL MATERIAL

1.0e+004 *
-0.01994250375000
-2.49281296875000
-0.05423213250000
-1.79339062500000
-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 *
1
7

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

-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
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
8

Universidad Nacional de Ingeniera

Facultad de Ingeniera Mecnica

Calculo por elementos Finitos

Flexion

1.0e-005 *
0.10447698051528
-0.10447698051528
-0.10447698051528
0.10447698051528
8. CONCLUSIONES:

Las matrices que se analizan en estos sistemas son de orden muy elevado, por
tal razn es necesario utilizar un lenguaje de programacin que nos permita
manejar las variables con mayor flexibilidad y poder generalizar el mtodo de
anlisis.

El anlisis de viga de seccin variable es la generalizacin del anlisis de una

viga se seccin constante.

El vector desplazamiento es desarrollado en base a la conectividad de los

elementos, por ello es importante manejar una tabla de conectividad ordenada y
secuencial.

En este caso de viga se seccin variable era de esperarse que cada elemento
tuviera 4 grados de libertad.

Cada elemento de la viga est sujeto a fuerzas y un momento; las fuerzas que
pueden ser de compresin o tensin directa mientras los momentos son de
flexin.

Como es propio de la viga, en este caso todas las cargas son aplicadas en los
nodos, adems los clculos se realizan despreciando la friccin en los nodos.

1
9