Facultad de ingeniería

mecánica
2019-2

ARMADURAS EN 3D

 Alumno:
LA MADRID VILLAVERDE RICARDO 20151178H
 Profesor: Msc. Ing. Abregu Leandro Edwin Asencion
 Sección: A

1
 INDICE

Enunciado del problema…………………………………………………………….(3)

Solución…………………….………………………………………………………...(4)

Método de cálculo…………………………………………………………………...(4)

Etapa de cálculos………………..…………………………………………………..(7)

Programa de Matlab……………………………...………………………………..(10)

Conclusiones…………………………………………………………………….....(14)

Bibliografía…………………………………………………………………………..(15)

2
 1. ENUNCIADO DEL PROBLEMA
Para la armadura de tres dimensiones que se muestra en la figura adjunta, que
está sometida a una carga de 889.644 N en el nodo 2, determinar:
a) El vector de desplazamientos nodales
b) El vector fuerza
c) Los esfuerzos en los elementos de la armadura

Fig. 1. Representación gráfica para el enunciado del problema

Datos del material:

𝐸 = 73084,43 𝑁/𝑚𝑚2

𝐴 = 1006.45 𝑚𝑚2

3
2. SOLUCIÓN
Se elaboran los cuadros de conectivdad y de coordenadas Nodales.

Tabla 1. Coordenadas Nodales

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

1 0 0 76.4

2 152.8 0 0

3 0 0 -76.4

4 0 152.4 0

Tabla 2. Cuadro de Conectividad

e NODOS GDL L e (mm) E e (N/mm 2 )

1 1 2 1 2 3 4 5 6 170.38 7.31*104
1
1
2 1 3 1 2 3 7 8 9 152.4 7.31*104
2
1
3 1 4 1 2 3 10 11 12 170.38 7.31*104
3
1
4 2 4 4 5 6 10 11 12 215.52 7.31*104
1
5 2 3 4 5 6 7 8 9 170.38 7.31*104
1
6
3 4 7 8 9 10 11 12 170.38 7.31*104
1

3. MÉTODO DE CÁLCULO

La matriz de transformación para un elemento recto en 3D está representado por la

ecuación 1.
𝑙 𝑚 𝑛 0 0 0
𝐿=[ ] (1)
0 0 0 𝑙 𝑚 𝑛
𝑙 , m, n son los cosenos directores:

4
Cálculo de la longitud de la barra:

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

3.2. Transformación de coordenadas

Desplazamientos nodales:
𝑞 ′ = 𝐿𝑞 (3)
Según el problema planteado se tiene cuatro nodos y al tratarse de una armadura
tridimensional se tendrá 12 grados de libertad.

Q=[𝑞1 𝑞2 𝑞3 … … … … … … 𝑞24 ] (4)

Cálculo de la matriz de rigidez local:

𝐾𝑒 = 𝐿𝑇 𝑘′𝑒 𝐿 (5)

𝒍𝟐 𝒍𝒎 𝒍𝒏 −𝒍𝟐 −𝒍𝒎 −𝒍𝒏

𝟐
𝒍𝒎 𝒎 𝒎𝒏 −𝒍𝒎 −𝒎𝟐 −𝒎𝒏
𝑲𝒆 = 𝒆𝒍 𝒆 𝒍𝒏𝟐
𝑬 𝑨 𝒎𝒏 𝒏𝟐 −𝒍𝒏 −𝒎𝒏 −𝒏𝟐 (6)
𝒆 −𝒍 −𝒍𝒎 −𝒍𝒏 𝒍𝟐 𝒍𝒎 𝒍𝒏
𝟐 𝟐
−𝒍𝒎 −𝒎 −𝒎𝒏 𝒍𝒎 𝒎 𝒎𝒏
[ −𝒍𝒏 −𝒎𝒏 −𝒏𝟐 𝒍𝒏 𝒎𝒏 𝒏𝟐 ]

Cálculo de los cosenos directores

Dadas las coordenadas de cada nodo, ver tabla 1, los cosenos directores respecto a
cada eje coordenado se calculan según:

𝒙𝟏 − 𝒙𝒐 𝒚𝟏 − 𝒚 𝒐 𝒛𝟏 − 𝒛𝒐
𝒍= 𝒎= 𝒍= (𝟕)
𝒍𝒆 𝒍𝒆 𝒍𝒆

5
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

1 0.894 0.000 -0.447

1
2 0.000 0.000 -1.000
1
3 0.000 0.894 -0.447
1
4 -0.707 0.707 0.000
1

5 -0.895 0.000 -0.447

1
6 0.000 0.895 0.447
1

Cálculo del vector fuerza

La formula que permite determinar la fuerza resultante en cada nodo es de la forma:

𝐹 =𝐾∗𝑄 (8)

6
Vector desplazamiento

Fig. 2. Representación gráfica de las deformaciones en los nodos

Según la figura 2, se tiene:

𝑄 = [𝑄1 𝑄2 𝑄3 𝑄4 𝑄5 𝑄6 𝑄7 𝑄8 𝑄9 𝑄10 𝑄11 𝑄12]𝑇 (9)

Las condiciones de contorno establecidas por el problema son:

𝑄1 = 𝑄3 = 𝑄7 = 𝑄10 = 𝑄11 = 0 (10)

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)

7
vector fuerza

Fig. 3. Representación gráfica de las fuerzas actuantes sobre la

armadura

el vector fuerza para la armadura es:

𝐹 = [𝐹1 𝐹2 𝐹3 𝐹4 𝐹5 𝐹6 𝐹7 𝐹8 𝐹9 𝐹10 𝐹11 𝐹12]𝑇

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

𝐹 = [𝑅1 𝑅2 𝑅3 0 − 889.644 0 𝑅7 𝑅8 𝑅9 𝑅10 𝑅11 𝑅12]𝑇

De la figura 2, las reacciones en determinadas direcciones tendrán valor de cero debido

a que no está restringidas a ningún apoyo rígido. Donde resulta:

𝐹 = [𝑅1 0 𝑅3 0 − 889.644 0 𝑅7 0 0 𝑅10 𝑅11 0]𝑇 (13)

8
 matriz de rigidez Global

Debido a la complejidad de la armadura, los valores de la matriz de rigidez se calcularon

con software (Matlab o Excel).

se aprecia que es una matriz simétrica, estás matrices son las que gobiernan el
comportamiento de cualquier tipo de armadura.

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 ]

Cálculo de los desplazamientos.

De (8), (12), (13) y de la matriz de rigidez global, eliminando convenientemente filas y
columnas correspondientes a los desplazamientos nulos, se tiene:

0 28779.69 0.00 0.00 0.00 0.00 0.00 14389.84 𝑄2

0 0.00 71782.08 −14220.20 −5.03 0.00 −14384.81 0.00 𝑄4
−889.644 0.00 −14220.20 14220.20 0.00 0.00 0.00 0.00 𝑄5
0 = 0.00 −5.03 0.00 14384.18 0.00 −7189.26 0.00 𝑄6
0 0.00 0.00 0.00 0.00 28782.21 14384.81 −14384.81 𝑄8
0 0.00 −14384.81 0.00 −7189.26 14384.81 54599.33 −7189.26 𝑄9
[ 0 ] [14389.84 0.00 0.00 0.00 −14384.81 −7189.26 14384.18 ] [𝑄12]

Los valores de los desplazamientos son:

NODO X (x103mm) Y (x103mm) Z (x103mm)

1 0 0.86 0
2 -16.83 -79.4 -2.77
3 0 1.91 -5.53
4 0 0 -1.71

9
Las reacciones nodales son:

NODO X (N) Y (N) Z (N)

1 444.822 0 0
3 444.822 0 0
4 889.644 889.644 0

Cálculo de los esfuerzos.

Utilizando la ecuación (11), se obtienen los esfuerzos en cada elemento finito:

ELEMENTO ESFUERZO (N/mm2)

1 -0.49
2 0.22
3 0.00
4 1.25
5 -0.50
6 0.00

PROGRAMA PRINCIPAL

clear all

clc

E=200000;

A=pi*900*ones(1,16)'; % Area de las secciones

% Ingresa coordenada de los nodos

nodeCoordinates=[0 0 0;

-35 0 0 ;

-400 400 0 ;

400 400 0 ;

205 205 980;

-205 205 980;

-205 -205 980;

205 -205 980;];

10
% Ingresa la conectividad de los nodos

elementNodes=[1 2;1 4;1 5;1 8;2 3; 2 7; 2 8; 3 4; 3 6; 3 7; 4 5; 4 6; 5

8; 5 6; 6 7; 7 8];

numberElements=size(elementNodes,1);

numberNodes=size(nodeCoordinates,1);

xx=nodeCoordinates(:,1);

yy=nodeCoordinates(:,2);

zz=nodeCoordinates(:,3);

% Para la estructura:

% displacements: vector de desplazamiente

% force : vector de fuerza

% stiffness: matriz de rigidez

% GDof: número global de grados de libertad

GDof=3*numberNodes;

U=zeros(GDof,1);

force=zeros(GDof,1);

% Aplicando la carga en los nodos 7 y 8.

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

% Salida de desplazamientos y reacciones

outputDisplacementsReactions(displacements,stiffness,...

GDof,prescribedDof)

% Esfuerzo de los elementos

stresses3Dtruss(numberElements,elementNodes,nodeCoordinates,...

displacements,E)

FUNCIONES AUXILIARES

Función que calcula la rigidez

11
function
[stiffness]=formStiffness3Dtruss(GDof,numberElements,elementNodes,numb
erNodes,nodeCoordinates,E,A);

stiffness=zeros(GDof);

% Calcula la matriz de rigidez del sistema

for e=1:numberElements;

% elementDof: Grado de libertad del elemento

indice=elementNodes(e,:) ;

elementDof=[3*indice(1)-2 3*indice(1)-1 3*indice(1)...

3*indice(2)-2 3*indice(2)-1 3*indice(2)] ;

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

L = sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1) +...

(z2-z1)*(z2-z1));

disp('L');

disp(L);

CXx = (x2-x1)/L;CYx = (y2-y1)/L;CZx = (z2-z1)/L;

disp(e)

disp(':')

disp(CXx)

disp(CYx)

disp(CZx)

T = [CXx*CXx CXx*CYx CXx*CZx ; CYx*CXx CYx*CYx CYx*CZx ; ...

CZx*CXx CZx*CYx CZx*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

Función que calcula desplazamientos

12
function displacements=solution(GDof,prescribedDof,stiffness,force)

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)

%MUESTRA LOS DESPLAZAMIENTOS Y REACCIONES EN FORMA TABULAR

%GDof=total numero de grados de libertad del problema

%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

Función que calcula el esfuerzo

function
stresses3Dtruss(numberElements,elementNodes,nodeCoordinates,displaceme
nts,E)

% Esfuerzo en un elemento 3D

fprintf('Esfuerzos en los elementos barra\n')

ff=zeros(numberElements,6);

13
for e=1:numberElements;

% elementDof: Grados de libertad de un elemento

indice=elementNodes(e,:) ;

elementDof=[3*indice(1)-2 3*indice(1)-1 3*indice(1) 3*indice(2)-2

3*indice(2)-1 3*indice(2)] ;

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

L = sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1) +...

(z2-z1)*(z2-z1));

CXx = (x2-x1)/L;CYx = (y2-y1)/L;CZx = (z2-z1)/L;

u=displacements(elementDof);

member_stress(e)=E/L*[-CXx -CYx -CZx CXx CYx CZx]*u;

fprintf('%3d %12.8f\n',e, member_stress(e));

end

5. CONCLUSIONES
 Los elementos finitos 2 y 4 están en compresión, ya que se obtuvieron
esfuerzos positivos.
 Los elementos finitos 1 y 5 están en tracción, ya que se obtuvieron
esfuerzos negativos.
 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.
 Para consideraciones de diseño, es de gran interés saber qué elementos
están más propensos a sufrir fallas. Los valores de esfuerzos que se

14
 determinaron nos permiten predecir el comportamiento bajo las
condiciones de diseño. Conociendo dichos valores se puede saber en qué
parte de la armadura se deben de usar materiales más resistentes.
6. BIBLIOGRAFIA
 “Análisis de elementos finitos”. Saeed Moaveni.
 “Introducción a los elementos finitos en ingeniería”. Ashok D. Belegundu,
Tiruphati R. Chandrupatla.
 Cuaderno profesor Cueva Ronald.

15