Informe 4 - Armadura 3D
Informe 4 - Armadura 3D
Informe 4 - Armadura 3D
INGENIERÍA
FACULTAD DE INGENIERÍA MECÁNICA
SECCION:
“A”
PROFESOR:
ING. ABREGU LEANDRO, EDWIN ASENCION
ALUMNOS:
AYALA ROJAS, CARLOS CHRISTIAN
20151333C
2019 – II
INDICE
Solución…………………….………………………………………………………...(4)
Método de cálculo………………………………………………………………..….(4)
Etapa de cálculos………………..………………………………………………….(7)
Programa de Matlab……………………………...…………………………………(9)
Conclusiones……………………………………………………………………...(14)
Bibliografía………………………………………………………………………...(15)
𝐸 = 73084,43 𝑁/𝑚𝑚2
𝐴 = 1006.45 𝑚𝑚2
2. SOLUCIÓN
Se elaboran los cuadros de conectivdad y de coordenadas Nodales.
1 0 0 914.4
2 1828.8 0 0
3 0 0 -914.4
4 0 1828.8 0
𝒙𝟏 − 𝒙𝒐 𝒚𝟏 − 𝒚 𝒐 𝒛𝟏 − 𝒛𝒐
𝒍= 𝒎= 𝒍= (𝟕)
𝒍𝒆 𝒍𝒆 𝒍𝒆
Usando la ecuación 7, se obtienen los valores de los cosenos directores.
Tabla 3. Valores de los Cosenos Directores
COSENOS DIRECTORES
Elemento l m n
𝐹 =𝐾∗𝑄 (8)
Vector desplazamiento
Esfuerzos nodales
El cálculo de los esfuerzos en los nodos se realiza utilizando la siguiente ecuación:
𝑄1
𝑄2
𝑄3
𝜎 𝑒 = (𝐸/𝑙 𝑒 )[−𝑙 −𝑚 −𝑛 𝑙 𝑚 𝑛] 𝑄4 (11)
𝑄5
[𝑄6]
4. ETAPA DE CÁLCULOS
Se presentan los valores incógnita y los valores conocidos.
Vector desplazamiento
De (9) y (10) se tiene
𝑄 = [0 𝑄2 0 𝑄4 𝑄5 𝑄6 0 𝑄8 𝑄9 0 0 𝑄12]𝑇 (12)
Vector fuerza
El vector fuerza para la armadura es:
Por las condiciones de contorno y de la figura 3, el vector fuerza queda limitado solo a
las fuerzas de reacción y a la carga aplicada en el nodo 2.
Se tendrá entonces
28779.67 0.00 −14389.84 −28779.67 0.00 14389.84 0.00 0.00 0.00 0.00 0.00 0.00
0.00 28779.69 −14389.84 0.00 0.00 0.00 0.00 0.00 0.00 0.00 −28779.69 14389.84
−14389.84 −14389.84 54610.65 14389.84 0.00 −7194.92 0.00 0.00 −40220.81 0.00 14389.84 −7194.92
−28779.67 0.00 14389.84 71782.08 −14220.20 −5.03 −28782.21 0.00 −14384.81 −14220.20 14220.20 0.00
0.00 0.00 0.00 −14220.20 14220.20 0.00 0.00 0.00 0.00 14220.20 −14220.20 0.00
14389.84 0.00 −7194.92 −5.03 0.00 14384.18 −14384.81 0.00 −7189.26 0.00 0.00 0.00
𝐾= 𝑁/𝑚𝑚
0.00 0.00 0.00 −28782.21 0.00 −14384.81 28782.21 0.00 14384.81 0.00 0.00 0.00
0.00 0.00 0.00 0.00 0.00 0.00 0.00 28782.21 14384.81 0.00 −28782.21 −14384.81
0.00 0.00 −40220.81 −14384.81 0.00 −7189.26 14384.81 14384.81 54599.33 0.00 −14384.81 −7189.26
0.00 0.00 0.00 −14220.20 14220.20 0.00 0.00 0.00 0.00 14220.20 −14220.20 0.00
0.00 −28779.69 14389.84 14220.20 −14220.20 0.00 0.00 −28782.21 −14384.81 −14220.20 71782.10 −5.04
[ 0.00 14389.84 −7194.92 0.00 0.00 0.00 0.00 −14384.81 −7189.26 0.00 −5.04 14384.18 ]
PROGRAMA en MATLAB
clear all
clc
E=200000;
nodeCoordinates=[0 0 0;
-35 0 0 ;
-400 400 0 ;
400 400 0 ;
numberElements=size(elementNodes,1);
numberNodes=size(nodeCoordinates,1);
xx=nodeCoordinates(:,1);
yy=nodeCoordinates(:,2);
zz=nodeCoordinates(:,3);
% Para la estructura:
GDof=3*numberNodes;
U=zeros(GDof,1);
force=zeros(GDof,1);
force(21)=-4000;
force(24)=-2000;
% Matriz de rigidez
[stiffness]=formStiffness3Dtruss(GDof,numberElements,elementNodes,numb
erNodes,nodeCoordinates,E,A);
% Condiciones de contorno
prescribedDof=[3 4 6 7 8 9 11 12]';
% Solución
displacements=solution(GDof,prescribedDof,stiffness,force);
outputDisplacementsReactions(displacements,stiffness,...
GDof,prescribedDof)
stresses3Dtruss(numberElements,elementNodes,nodeCoordinates,...
displacements,E)
FUNCIONES AUXILIARES
stiffness=zeros(GDof);
for e=1:numberElements;
indice=elementNodes(e,:) ;
x1=nodeCoordinates(indice(1),1);
y1=nodeCoordinates(indice(1),2);
z1=nodeCoordinates(indice(1),3);
x2=nodeCoordinates(indice(2),1);
y2=nodeCoordinates(indice(2),2);
z2=nodeCoordinates(indice (2),3);
(z2-z1)*(z2-z1));
disp('L');
disp(L);
disp(e)
disp(':')
disp(CXx)
disp(CYx)
disp(CZx)
disp('T')
disp(T)
E*A(e)/L*[T -T ; -T T]
stiffness(elementDof,elementDof)=...
stiffness(elementDof,elementDof)+E*A(e)/L*[T -T ; -T T];
end
activeDof=setdiff([1:GDof]',[prescribedDof]);
stiffness(activeDof,activeDof)
force(activeDof)
U=stiffness(activeDof,activeDof)\force(activeDof);
displacements=zeros(GDof,1);
displacements(activeDof)=U;
end
Función que calcula las reacciones
function
outputDisplacementsReactions(displacements,stiffness,GDof,prescribedDo
f)
%dezplazamientos
disp('Displacements')
%displacements=displacements1;
jj=1:GDof; %format
A=[jj' displacements];
disp(A)
%REACCIONES
F=stiffness*displacements;
reactions=F(prescribedDof);
disp('reactions')
B=[prescribedDof reactions];
disp(B)
end
% Esfuerzo en un elemento 3D
ff=zeros(numberElements,6);
for e=1:numberElements;
indice=elementNodes(e,:) ;
x1=nodeCoordinates(indice(1),1);
y1=nodeCoordinates(indice(1),2);
z1=nodeCoordinates(indice(1),3);
x2=nodeCoordinates(indice(2),1);
y2=nodeCoordinates(indice(2),2);
z2=nodeCoordinates(indice(2),3);
(z2-z1)*(z2-z1));
u=displacements(elementDof);
end
5. CONCLUSIONES
Los elementos finitos 1 y 5 están en compresión, ya que se obtuvieron esfuerzos
negativos.
Los elementos finitos 2 y 4 están en tracción,
Los elementos finitos que están sometidos a valores de esfuerzo bajos (valores
prácticamente nulos) son los elementos finitos 3 y 6.
El elemento finito 4 es el que está sometido a mayor esfuerzo. Esto debido a que
la carga de 889.644 N está actuando en su extremo más vulnerable mientras su
otro extremo está empotrado en el plano de la carga.
Los siguientes cuadros evidencian la gran utilidad del software (en nuestro caso
ANSYS 18.2) para la simulación y predicción de fallas en las armaduras.
Conociendo dichos valores se puede saber en qué parte de la armadura se
deben de usar materiales más resistentes.
Desplazamientos en la armadura
MATLAB ANSYS 18.2 ERROR
NODOS
X [10^-3 mm] Y [10^-3 mm] Z [10^-3 mm] X [10^-3 mm]Y [10^-3 mm] Z [10^-3 mm] X [mm] Y [mm] Z [mm]
1 0 5.4415 0 0 -2.5323 0 0 7.9738 0
2 -6.63 -31.1266 -1.0084 -6.1509 -29.005 -1.0097 -0.4791 -2.1216 0.0013
3 0 5.4415 -2.1769 0 3.539 -2.0194 0 1.9025 -0.1575
4 0 0 -1.0084 0 0 5.0616 0 0 -6.07
Reacciones en la armadura
MATLAB ANSYS 18.2 ERROR
NODOS
X [N] Y [N] Z [N] X [N] Y [N] Z [N] X [N] Y [N] Z [N]
1 444.82 0 0 444.82 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0
3 444.82 0 0 444.82 0 0 0 0 0
4 -889.64 889.64 0 -889.64 889.64 0 0 0 0
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies: