Armaduras Planas
Armaduras Planas
FACULTAD DE INGENIERIA
MECANICA
CALCULO POR ELEMENTOS FINITOS
TERCERA PRACTICA: ARMADURAS PLANAS
ALUMNO:
CODIGO:
PROFESOR:
ING. RONALD CUEVA PACHECO
CURSO SECCIN:
MC516 D
2012-II
20092563A
NDICE
2. Solucin................................................................................................4
4. Vector Carga........................................................................................6
5. Matriz de Rigidez.................................................................................7
7. Esfuerzos .........................................................................................14
8. Resultados16
9. Diagrama de Flujo.............................................................................17
10. Uso de Matlab
(DIGITAL).........................................................................................19
11. Conclusiones23
2
Datos:
E =3.1x105 N/mm2
D= 150mm
Cargas a considerar (seccin D)
PA=5000N , PB=4000N , PE=2000N
L=1500 mm
2. SOLUCION
MODELADO DEL CUERPO REAL
Se tienen 6 elementos con 5 nodos y 10 grados de libertad. Las coordenadas para los nodos
son:
Nodo
1
2
3
4
5
X(mm)
0
1500
3000
1500
0
Y(mm)
0
0
0
1500
1500
Elemento
1
2
3
4
5
6
Nodos
(1) (2)
1
2
2
3
3
4
4
2
4
1
4
5
1
1
3
5
7
7
7
GDL
2 3 4
2 3 4
4 5 6
6 7 8
8 3 4
8 1 2
8 9 10
Le (mm)
1500.00
1500.00
2121.321
1500.00
2121.32
1500.00
Ae (
1963.5
1963.5
1963.5
1963.5
1963.5
1963.5
Ee (N/
Donde
, pues la viga esta empotrada en los nodos 1 y 5,
los dems desplazamientos son incgnitas que tendremos que calcular.
4. VECTOR CARGA
F1 R1
F 2 R2
F3 0
F 4 2000
F 5 5000
F
N
F 6 4000
F7 0
F8 0
F 9 R3
F10 R 4
5. MATRIZ DE RIGIDEZ
A continuacin pasamos a calcular la matriz de Rigidez Global, que est determinada
por la siguiente ecuacin:
( ) *
Respecto a
Respecto a (X, Y):
Elemento 1:
x2 x1
1
L
y y1
m1 2
0
L
l1
K1=405790*[
(traccin simple)
donde
( ) [
Resulta:
Elemento 2:
x3 x2
1
L
y y2
m2 3
0
L
l2
K2=405790*[
Elemento 3:
l3
x 4 x3
m3
L 2
y 4 y3
L 2
2
2
K3=286937*[
Elemento 4:
x4 x2
0
L
y y2
m4 4
1
L
l4
K4=405790*[
Elemento5:
x4 x1
1
2
L 2
y y1
m5 4
1
2
L 2
l5
K5=286937*[
10
Elemento 6:
x 4 x1
1
L
y y1
m1 4
0
L
l1
K6=405790*[
17530087
0
35060174.16
0
0
0
0
6197822
0
0
17530089
0
k
0
0
0
0
6197822 6197822
0
0
0
6197822
6197822 6197822
0
0
0
0
0
0
0
0
11
0
0
17530089
0
23727909
6197822
6197822
6197822
0
0
6197822
6197822
6197822
6197822
0
0
0
0
0
0
0
0
0
6197822
0
0
6197822 6197822 6197822
0
0
.
6197822
6197822 6197822
0
0
6197822
2992531
0
17530087 0
6197822
0
29925731 17530087 0
0
17530087 17530087 17530087 0
0
0
0
0
0
R1 K K1 K 2 K 3 K 4 K 5 K 6
R2
23727909 6197822 17530089
6197822 6197822
0
0
17530087
0
35060174.16
2000
0
0
0
5000
0
17530089
F
=k 0
4000
0
0
0
0
6197822 6197822
0
0
6197822 6197822
0
0
0
0
R3
0
0
0
R 4
17530089
6197822
0
0
0
0
23727909 6197822
6197822 6197822
0
6197822 6197822
6197822 6197822 6197822
0
6197822 6197822
6197822 6197822
0
0
0
0
0
0
0
6197822
0
0
6197822 6197822
0
0
.
6197822 6197822
0
0
2992531
0
17530087 0
0
29925731 17530087 0
17530087 17530087 17530087 0
0
0
0
0
12
Q3 0.005367mm
Q4 0.031421mm
Q5 0.010734mm
Q6 0.054296mm
Q7 0.009022mm
Q8 0.027428mm
R1 20816N
R2 11408N
R3 15816N
R4 0
13
7. ESFUERZOS
Para calcular los valores de los esfuerzos por elemento, aplicamos la siguiente
ecuacin:
Q1'
'
Q
E
e l m l m 2'
Q3
L e
'
Q4
Y obtenemos lo siguiente:
0
3.1x10 5
0.8319N / mm 2
1 0 1 0
1
0.005367
1500
0.031421
0.005367
0.031421
5
3.1x10
0.8319N / mm 2
1 0 1 0
2
0.010734
1500
0.054296
14
0.010734
0.054296
3.1x10 5
0.55124N / mm 2
1 / 2 1 / 2 1 / 2 1 / 2
0
.
009022
2828
.
427
0.027428
0.005367
0.031421
5
3.1x10
0.618915N / mm 2
0 1 0 1
0.009022
2121.321
0.027428
0
0
3.1x10 5
1.42652N / mm 2
1 / 2 1 / 2 1 / 2 1 / 2
0
.
009022
2121
.
321
0.027428
0
5
0
3.1x10
1.398N / mm 2
1 0 1 0
0.009022
1500
0.027428
15
8. RESULTADOS
Finalmente los resultados son mostrados en la siguiente tabla:
R1 6711 .2 N
R2 6665 .97 N
R3 18854 .2 N
R4 0 N
1 0.00062 N / mm 2
2 0.007233 N / mm 2
3 0.55124 N / mm 2
4 0.618915 N / mm 2
5 1.42652 N / mm 2
6 1.398 N / mm 2
16
9. DIAGRAMA DE FLUJO
INICIO
Calcula lee, l y m de
cada elemento
Para i=1:n
17
Para i=1:n
Calculo de esfuerzos:
s(1)=E/le(i)*[-l(i) -m(i) l(i)
m(i)]*Q(1:4)
Imprime Reacciones,
desplazamientos y
esfuerzos
FIN
18
%Ingreso de Datos:
E=input('Ingrese el modulo de Elasticidad (N/mm^2,MPa): ');
d=input('Ingrese diametro de la seccion de cada elemento (mm): ');
PA=input('Ingrese magnitud de fuerza PA (N): ');
PB=input('Ingrese magnitud de fuerza PB (N): ');
PC=input('Ingrese magnitud de fuerza PC (N): ');
A=(pi*d^2)/4;
19
le(5)=sqrt((x(4)-x(1))^2+(y(4)-y(1))^2);%Elemento 5
le(6)=sqrt((x(5)-x(4))^2+(y(5)-y(4))^2);%Elemento 6
%Calculo de Cosenos directores:
l(1)=(x(2)-x(1))/le(1);m(1)=(y(2)-y(1))/le(1);%Elemento 1
l(2)=(x(3)-x(2))/le(2);m(2)=(y(3)-y(2))/le(2);%Elemento 2
l(3)=(x(4)-x(3))/le(3);m(3)=(y(4)-y(3))/le(3);%Elemento 3
l(4)=(x(4)-x(2))/le(4);m(4)=(y(4)-y(2))/le(4);%Elemento 4
l(5)=(x(4)-x(1))/le(5);m(5)=(y(4)-y(1))/le(5);%Elemento 5
l(6)=(x(4)-x(5))/le(6);m(6)=(y(4)-y(5))/le(6);%Elemento 6
fprintf('\nTabla de Datos: \n')
fprintf('\nElemento
Long. Elemento
\n')
for i=1:6
fprintf('%4d %18.4f %16.4f %16.4f\n',i,le(i),l(i),m(i))
end
20
R1\n')
R2\n')
fprintf('\n %d',f(3:8)')
fprintf('\n
21
R9\n')
fprintf('
R10\n')
%d',Q)
R1:
%d\n',R1)
fprintf('
R2:
%d\n',R2)
fprintf('
R9:
%d\n',R9)
fprintf('
R10: %d\n',R10)
%Distribucion de Esfuerzos:
Js=[];
for i=1:6
J=E/le(i)*[-l(i) -m(i) l(i) m(i)];
Js=[Js;J];
end
s(1)=Js(1,:)*Q(1:4);
s(2)=Js(2,:)*Q(3:6);
s(3)=Js(3,:)*Q(5:8);
s(4)=Js(4,:)*[Q(3:4);Q(7:8)];
s(5)=Js(5,:)*[Q(1:2);Q(7:8)];
s(6)=Js(6,:)*[Q(9:10);Q(7:8)];
fprintf('\nTabla de Esfuerzos: \n')
fprintf('\nElemento
Esfuerzo (N/mm^2,MPa)\n')
for i=1:6
fprintf('%4d %18.4f \n',i,s(i))
22
11. CONCLUSIONES
23
24
25