Informe4 MC516 Pumarrumi
Informe4 MC516 Pumarrumi
Informe4 MC516 Pumarrumi
LABORATORIO N°4
ARMADURA EN EL ESPACIO
CURSO:
CÁLCULO POR ELEMENTOS FINITOS
SECCIÓN:
E
ESTUDIANTE:
PUMARRUMI ESCOBAR JASSELY MILAGROS 20184116A
DOCENTE:
CUEVA PACHECO RONALD
LIMA, 2021
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
ÍNDICE
2. SOLUCIÓN...........................................................................................................3
3. DIAGRAMA DE FLUJO........................................................................................6
4. CÓDIGO EN MATLAB..........................................................................................7
5. RESULTADOS EN MATLAB..............................................................................10
6. CONCLUSIONES.............................................................................................. 14
7. BIBLIOGRAFIA…............................................................................................... 14
1
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
CUARTA PRÁCTICA
En la figura se muestra una viga con dos apoyos simples, conformada mediante una
armadura en 3D; con las cargas de servicio indicadas.
HALLAR:
DATOS:
• Dimensiones:
o L = 600 mm
2
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
o b = 500 mm
o Sección de barras: ∅ = 50 𝑚𝑚
• Material:
o E=2.1*105 MPa
• Cargas:
o PA = 10000 N
o PB = 8000 N
▪ Ángulo de inclinación 𝛽 = 70
2. SOLUCIÓN:
2.1 MODELADO DEL CUERPO REAL
3
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
Nodos GDL
e
(1) (2) 1 2 3 4 5 6
1 1 2 27 28 29 1 2 3 600 1963.495409
2 2 5 1 2 3 8 9 10 346.41016 1963.495409
3 1 5 27 28 29 8 9 10 692.82032 1963.495409
4 2 3 1 2 3 4 5 6 600 1963.495409
5 3 5 4 5 6 8 9 10 692.82032 1963.495409
6 5 6 8 9 10 11 12 13 600 1963.495409
7 3 6 4 5 6 11 12 13 346.41016 1963.495409
8 3 4 4 5 6 30 7 31 600 1963.495409
9 4 6 30 7 31 11 12 13 692.82032 1963.495409
10 7 8 32 33 34 14 15 16 600 1963.495409
11 8 11 14 15 16 21 22 23 346.41016 1963.495409
12 7 11 32 33 34 21 22 23 692.82032 1963.495409
13 8 9 14 15 16 17 18 19 600 1963.495409
14 9 11 17 18 19 21 22 23 692.82032 1963.495409
15 11 12 21 22 23 24 25 26 600 1963.495409
16 9 12 17 18 19 24 25 26 346.41016 1963.495409
17 9 10 17 18 19 35 20 36 600 1963.495409
18 10 12 35 20 36 24 25 26 692.82032 1963.495409
19 1 7 27 28 29 32 33 34 500 1963.495409
20 1 8 27 28 29 14 15 16 781.02497 1963.495409
21 2 8 1 2 3 14 15 16 500 1963.495409
22 3 8 4 5 6 14 15 16 781.02497 1963.495409
23 3 9 4 5 6 17 18 19 500 1963.495409
24 4 9 30 7 31 17 18 19 781.02497 1963.495409
25 4 10 30 7 31 35 20 36 500 1963.495409
26 5 11 8 9 10 21 22 23 500 1963.495409
27 6 11 11 12 13 21 22 23 781.02497 1963.495409
28 6 12 11 12 13 24 25 26 500 1963.495409
29 1 11 27 28 29 21 22 23 854.40037 1963.495409
30 2 11 1 2 3 21 22 23 608.27625 1963.495409
31 3 11 4 5 6 21 22 23 854.40037 1963.495409
32 3 12 4 5 6 24 25 26 608.27625 1963.495409
33 4 12 30 7 31 24 25 26 854.40037 1963.495409
4
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
𝐹𝑖 = 𝐾𝑖∫ 𝑄∫
5
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
𝐸 𝑄
𝜎𝑒 = ( ) [−1 1] [ 𝑖 ]
𝑙 𝑒 𝑄𝑖+1
3. DIAGRAMA DE FLUJO
INICIO
LECTURA DE DATOS
CALCULO DE ÁREAS DE
ELEMENTOS FINITOS
CALCULO
LA MATRIZ DE RIGIDEZ GLOBAL
ECUACIÓN DE RIGIDEZ Y
CONDICIONES DE CONTORNO
CALCULO DE ESFUERZOS
𝐸 𝑄
𝜎𝑒 = ( ) [−1 1] [ 𝑖 ] − (𝐸𝛼)𝑒 ∆𝑇
𝑙 𝑒 𝑄𝑖+1
IMPRESIÓN DE RESULTADOS
FIN
6
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
4. CÓDIGO EN MATLAB
clc; clear;
format
%Datos de entrada (Modificables)
numelementos=33;
El=input ('Ingrese el modulo de elasticidad o Young MPa:');
E(1:numelementos)=El;
pi=3.14159265359;
b=input('Ingrese valor de b de la armadura en mm:');
l=input('Ingrese valor de L de la armadura en mm:');
d=input ('Ingrese sección de la barra(diametro)mm:');
PA=input('Ingrese valor de fuerza A en N:');
PB=input('Ingrese valor de fuerza B en N:');
beta=input('Ingrese angulo beta en ° :');
area=(pi*d*d/4)*ones(1,33);
nodos=[b 0 0 ; b l 0 ; b 2*l 0 ; b 3*l 0 ; b l -346.41016 ; b 2*l -
346.41016 ; 0 0 0 ; 0 l 0 ; 0 2*l 0 ; 0 3*l 0 ; 0 l -346.41016 ; 0 2*l -
346.41016];
UnionNodos=[1 2;2 5;1 5;2 3;3 5;5 6;3 6;3 4;4 6;7 8;8 11;7 11;8 9;9 11;11
12;9 12;9 10;10 12;1 7;1 8;2 8;3 8;3 9;4 9;4 10;5 11;6 11;6 12;1 11;2 11;3
11;3 12;4 12];
Desplazamientos=[0 0 0 1 1 1 1 1 1 0 1 0 1 1 1 1 1 1 0 0 0 1 1 1 1 1 1 0 1
0 1 1 1 1 1 1];
Fuerzas=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 -PA 0 PB*sind(beta) -PB*cosd(beta) 0
0 0 0 0 0 0 0 0 0 0 0 0 0 -PA 0 PB*sind(beta) -PB*cosd(beta)];
L = zeros(1,numelementos);
Cx = zeros(1,numelementos);
Cy = zeros(1,numelementos);
Cz = zeros(1,numelementos);
LAMDA = zeros(6,6);
for i = 1:numelementos
indice = UnionNodos(i,:);
P1 = [nodos(indice(1),1) nodos(indice(1),2) nodos(indice(1),3)];
P2 = [nodos(indice(2),1) nodos(indice(2),2) nodos(indice(2),3)];
L(i) = norm(P1-P2);
Cx(i) = (P2(1) - P1(1))/ L(i);
Cy(i) = (P2(2) - P1(2))/ L(i);
Cz(i) = (P2(3) - P1(3))/ L(i);
lamda = [Cx(i)^2 Cx(i)*Cy(i) Cx(i)*Cz(i) ; Cy(i)*Cx(i) Cy(i)^2
Cy(i)*Cz(i) ;...
Cz(i)*Cx(i) Cz(i)*Cy(i) Cz(i)^2];
LAMDA(:,:,i) = [lamda -lamda ; -lamda lamda];
end
for i = 1:numelementos
7
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
C(j(2),j(1),i) = B(2,1,i);
C(j(2),j(2),i) = B(2,2,i);
end
A
S = 3*size(nodos,1); %Dimensiones de la matriz global
m = cell(S/3,S/3);
for i = 1:size(nodos,1)
for j = 1:size(nodos,1)
v = find(Desplazamientos==0);
MGR = MG;
MGR(v,:) = 0;
MGR(:,v) = 0;
indicefil = zeros(1,S);
indicecol = zeros(1,S);
for i = 1:S
if MGR(i,:) == 0
indicefil(i) = i;
end
if MGR(:,i) == 0
indicecol(i) = i;
end
end
MGR(indicefil~=0,:) = []; %Eliminar filas y columnas de ceros para tener
la matriz global reducida
MGR(:,indicecol~=0) = []
Fuerzas(indicefil~=0) = []; %Eliminar filas y columnas de ceros de las
fuerzas
%Calcular los desplazamientos nodales en ambas direcciones
d = MGR\Fuerzas';
dfinal = zeros(S,1);
k = 1;
8
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
% <-Diagrama 3D-->
d3 = reshape(dfinal,[3,size(nodos,1)])';
NodosDesp = nodos + d3;
d4 = reshape(Desplazamientos,[3,size(nodos,1)])';
for j = 1:numelementos
indice = UnionNodos(j,:);
line([nodos(indice(1),1) nodos(indice(2),1)],...
[nodos(indice(1),2) nodos(indice(2),2)],...
[nodos(indice(1),3) nodos(indice(2),3)],...
'LineWidth',1.5,'Color','k');
hold on
line([NodosDesp(indice(1),1) NodosDesp(indice(2),1)],...
[NodosDesp(indice(1),2) NodosDesp(indice(2),2)],...
[NodosDesp(indice(1),3) NodosDesp(indice(2),3)],...
'LineWidth',1,'Color','b');
end
for i = 1:size(nodos,1)
plot3(nodos(i,1),nodos(i,2),nodos(i,3),'ro','MarkerSize',6,'MarkerFaceColo
r','r');
hold on
end
grid on
NodosDespx = NodosDesp(:,1);
NodosDespy = NodosDesp(:,2);
NodosDespz = NodosDesp(:,3);
[fil,col] = find(d4 ~= 0);
plot3(NodosDespx(fil),NodosDespy(fil),NodosDespz(fil),'og','MarkerSize',6,
'MarkerFaceColor','g');
view(45,45);
axis equal
9
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
DATOS INGRESADOS
10
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
11
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
ESFUERZOS
REACCIONES
12
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
DEFORMADA EN LA ARMADURA
13
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM
6. CONCLUSIONES
• Para la elaboración de los cálculos, el software Matlab fue de mucha utilidad, debido
a que en el informe hubo cálculos con matrices de orden 33 que resultaban
operativos.
• En los nodos establecido 1,7,4 y 10 extremos las deformaciones obtenidas son cero,
hay 10 barras en las cuales no hay deformación.
• También están los desplazamientos pequeños, del orden de los milímetros, que son
efecto únicamente de las deformaciones por tensión o compresión de las barras que
componen la armadura.
• Hay barras que presentan valor negativo en los esfuerzos, lo cual significa que al
asumir las cargas son de sentido contrario, tener en cuenta que no se ha considerado
el peso propio de cada viga.
• Para reducir el número de ecuaciones y cálculos en las operaciones para las
reacciones, las condiciones de contorno son útiles.
7. BIBLIOGRAFIA
14