UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA MECÁNICA

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

1. ENUNCIADO DEL PROBLEMA...........................................................................2

2. SOLUCIÓN...........................................................................................................3

2.1. MODELADO DEL CUERPO REAL.................................................................3

2.2. GRADOS DE LIBERTAD NODALES.............................................................4

2.3. VECTOR CARGA...........................................................................................5

2.4. MATRIZ DE RIGIDEZ....................................................................................5

2.5. ECUACIÓN DE RIGIDEZ Y CONDICIÓN DE CONTORNO..........................5

2.6. ESFUERZOS DE CADA ELEMENTO...........................................................6

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

1. ENUNCIADO DEL PROBLEMA

En la figura se muestra una viga con dos apoyos simples, conformada mediante una
armadura en 3D; con las cargas de servicio indicadas.

La unidad rígida en el espacio es el tetraedro; todos los espacios de la armadura en 3D

deben ser tetraedros:

HALLAR:

• El esfuerzo en cada elemento finito.

• Las reacciones en los apoyos.

• La deformada de la armadura (gráfico).

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

Figura No1. Vista isométrica de la armadura espacial.

A continuación, el siguiente cuadro muestra la ubicación espacial de cada uno de los 12

nodos establecidos para la armadura en el espacio.

NODO X(mm) Y(mm) Z(mm)

1 500 0 0
2 500 600 0
3 500 1200 0
4 500 1800 0
5 500 600 -346.41
6 500 1200 -346.41
7 0 0 0
8 0 600 0
9 0 1200 0
10 0 1800 0
11 0 600 -346.41
12 0 1200 -346.41

Tabla No1. Coordenadas geométricas de los nodos de la armadura en el espacio.

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

Tabla No2. Cuadro de conectividad de la armadura en el espacio.

2.2. GRADOS DE LIBERTAD NODALES(GDL) (VECTOR DESPLAZAMIENTO)

Debido al gran número de grados de libertad de la armadura espacial el vector
desplazamiento no podrá ser representado en su totalidad. El vector desplazamiento
estará dado por:
𝑄1
𝑄2
𝑄𝑗 = 𝑄3 𝑚𝑚
⋮
[𝑄36 ]
Donde 𝑄27 = 𝑄28 = 𝑄29 = 𝑄30 = ⋯ = 𝑄36 = 0, ya que la armadura está fijada en los
nodos 1 y 7 en su totalidad; y en los nodos 4 y 10 en las direcciones ′𝑥′ y ′𝑧′. Por lo tanto,
los demás desplazamientos deberán ser calculados.

4
 INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM

2.3. VECTOR CARGA

El vector carga estará dado por las cargas puntuales y las fuerzas de reacción pues no
se considera el peso de las barras ni variación de temperatura.
𝐹1
𝐹2
𝐹𝑖 = 𝐹3 𝑁
⋮
[𝐹36 ]
2.4. MATRIZ DE RIGIDEZ LOCALES
Sea la matriz de rigidez local para un elemento cualquiera en el espacio:

−𝑙 2 + 𝑙𝑚 +𝑙𝑛 + −𝑙 2 −𝑙𝑚 −𝑙𝑛

−𝑙𝑚 + 𝑚2 +𝑚𝑛 −𝑙𝑚 −𝑚2 −𝑚𝑛
𝑒′
𝐸𝐴 𝑒 +𝑙𝑛 −𝑚𝑛 +𝑛2 −𝑙𝑛 −𝑚𝑛 −𝑛2
𝑘𝑟𝑠 =( )
𝑙𝑒 −𝑙 2 −𝑙𝑚 −𝑙𝑛 +𝑙 2 𝑙𝑚 𝑙𝑛
2
−𝑙𝑚 −𝑚 −𝑚𝑛 +𝑙𝑚 𝑛2 𝑚𝑛
[ −𝑙𝑛 −𝑚𝑛 −𝑛2 +𝑙𝑛 𝑚𝑛 𝑛2 ]
Donde los cosenos directores en función de las coordenadas nodales serán:
𝑥2 − 𝑥1
𝑙=
𝑙𝑒
𝑦2 − 𝑦1
𝑚=
𝑙𝑒
𝑧2 − 𝑧1
𝑛=
𝑙𝑒

𝑙𝑒 = √(𝑥2 − 𝑥1 )2 + (𝑦2 − 𝑦1 )2 + (𝑧2 − 𝑧1 )2

2.5. MATRIZ DE RIGIDEZ GLOBAL

A continuación, pasamos a calcular la Matriz de Rigidez Global para la armadura,
que está determinada por la siguiente ecuación:
33
𝑒′
𝐾𝑖𝑗 = ∑ 𝑘𝑟𝑠
𝑒=1

2.6. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO

La ecuación de rigidez está determinada por la siguiente ecuación:

𝐹𝑖 = 𝐾𝑖∫ 𝑄∫

5
 INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM

2.7. ESFUERZO DE CADA ELEMENTO

Para calcular los valores de los esfuerzos por elemento, aplicamos la siguiente ecuación:

𝐸 𝑄
𝜎𝑒 = ( ) [−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 LAS REACCIONES

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)];

%Construcción de las matrices de rigidez de cada barra.

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

k = (E.*area)./L; %Constante AE/L para cada elemento

A = zeros(6,6); %Inicializamos las matrices de rigidez de cada elemento

%Ensamble de la matriz global de rigidez

for i = 1:numelementos

%Calculo de la matriz de rigidez de cada elemento (AE/L)*k_local

A(:,:,i) = k(i)*LAMDA(:,:,i);

7
INFORME N°4 DE ARMADURA EN EL ESPACIO UNI-FIM

%Convertimos el arreglo A en celdas pero las matrices de 6x6 de A

las dividimos en 4 paquetes de 3x3
j = UnionNodos(i,:); %Par de índice de cada elemento
B(:,:,i) = mat2cell(A(:,:,i),[3 3],[3 3]);

%Asignamos cada paquete en los índices correspondientes de la matriz

global de rigidez
C(j(1),j(1),i) = B(1,1,i);
C(j(1),j(2),i) = B(1,2,i);

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)

%En cada vuelta recogemos todos los elementos con el mismo

índice(i,j), los superponemos y sumamos entre sí.
clear x
x(:,:,:) = cell2mat(reshape(C(i,j,:),1,[],numelementos));
m(i,j) = {sum(x,3)};

%Si en ese índice (i,j) no se asignó ningú paquete, metemos un

paquete de ceros de 2x2
if size(m{i,j}) == [0 0]
m(i,j) = {zeros(3,3)};
end
end
end
MG = cell2mat(m) %Convertimos la matriz global en un arreglo numérico

%Reducir la matriz global

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

for i = 1:length(Desplazamientos) %Insertar los desplazamientos calculados

en el vector de desplazamientos original
if Desplazamientos(i) == 0
dfinal(i,1) = 0;
else
dfinal(i,1) = d(k);
k = k+1;
end
end
% <-Resultados- >
d2 = mat2cell(dfinal,3*ones(1,size(nodos,1)),1); %Dividimos dfinal en
paquetes de 3x1
Esfuerzos = zeros(1,numelementos);
Flocal = zeros(numelementos,6);
j = 1;
for i = 1:numelementos
indice = UnionNodos(i,:);
Esfuerzos(i) = (E(i)./L(i)) * [-Cx(i) -Cy(i) -Cz(i) Cx(i) Cy(i) Cz(i)]
* [d2{indice(1,1)} ; d2{indice(1,2)}];
Flocal(i,:) = A(:,:,i)*[d2{indice(1,1)} ; d2{indice(1,2)}];
j = j + 2;
end
dfinal %Imprimir resultados
disp('LOS ESFUERZOS (MPa):');
disp(Esfuerzos);
Reacciones=MG*dfinal %Reacciones globales
Flocal

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

5. RESULTADOS OBTENIDOS EN MATLAB:

DATOS INGRESADOS

Ingrese el modulo de Young en MPa:2.1*10^5

Ingrese el valor de b de la armadura en mm :500

Ingrese el valor de L de la armadura en mm :600

Ingrese el valor de la fuerza A en N :10000

Ingrese el valor de la fuerza B en N :8000

Ingrese el valor del ángulo beta en °: 70

RESULTADOS OBTENIDOS DE LOS 33 ELEMENTOS FINITOS

MATRIZ DE RIGIDEZ LOCAL

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

Para realizar esto, la armadura fue simulada en el software ANSYS.

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

“Introducción a los elementos finitos en ingeniería”. Ashok D. Belegundu, Tiruphati R.

Chandrupatla.

“Análisis de elementos finitos”. Saeed Moaveni.

Teoría brindada en las clases del Ing. Cueva Ronald en la sección E.

14