Facultad de

CÁLCULO POR ELEMENTOS FINITOS

ingeniería
mecánica
2018-1

ARMADURAS EN 3D

 Alumno:
TARAZONA PINEDO, Wilmer Billi 20141371j
 Profesor: Msc. Ing. Abregu Leandro Edwin Asencion
 Sección: A

1
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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 914.4

2 1828.8 0 0

3 0 0 -914.4

4 0 1828.8 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 2044.66 7.31*104
1
1
2 1 3 1 2 3 7 8 9 1828.8 7.31*104
2
1
3 1 4 1 2 3 10 11 12 2044.66 7.31*104
3
1
4 2 4 4 5 6 10 11 12 2586.314 7.31*104
1
5 2 3 4 5 6 7 8 9 2044.66 7.31*104
1
6
3 4 7 8 9 10 11 12 2044.66 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 POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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. ESTAPA 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
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

% 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
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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
 CÁLCULO POR ELEMENTOS FINITOS

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 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.
 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
determinaron nos permiten predecir el comportamiento bajo las

14
 CÁLCULO POR ELEMENTOS FINITOS

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