UNIVERSIDAD NACIONAL DE INGENIERIA

FACULTAD DE INGENIERIA MECANICA

LABORATOIO N°1

TRACCION
CURSO: CALCULO POR ELEMENTOS FINITOS MC516-E

ALUMNO: QUISPE TICONA, MANUEL ALONSO 20170100K

DOCENTE: ING. RONALD CUEVA PACHECO

LIMA-PERU

2020-1
 INDICE
CAP. 1 ENUNCIADO........................................................................................3

1. ENUNCIADO:......................................................................................3

2. MODELADO DEL CUERPO REAL:....................................................4

3. GRADOS DE LIBERTAD NODALES:.................................................5

4. VECTOR CARGA:..............................................................................5

5. MATRIZ DE RIGIDEZ:........................................................................7

6. EUACION DE RIGIDEZ Y CONDICIONES DE CONTORNO:...........7

7. ESFUERZOS......................................................................................8

8. RESULTADOS....................................................................................9

CAP. 2 DIAGRAMA DE FLUJO.....................................................................10

CAP. 3 DIGITACION Y EJECUCION.............................................................11

CAP 4. CONCLUSIONES DEL TRABAJO REALIZADO...............................12

 CAP. 1 ENUNCIADO

1. ENUNCIADO:
En la figura se muestra una placa triangular simétrica; con espesor constante
igual a 120 mm; está sometida a su peso propio y a una fuerza concentrada, tal
como se indica.

Datos:

L = 1500 mm

h = 1000 mm

t = 120 mm (espesor)

P = 50000 N

ɣ = 8 gr-f/cm3 = 78.45x10-6 N/mm3

E = 2 x 105 N/mm2

HALLAR:

 La distribución de esfuerzos a lo largo de la barra; modelándola con

cuatro (04) elementos finitos (unidimensionales).
 La fuerza de reacción en el apoyo.
 2. MODELADO DEL CUERPO REAL:
Para resolver el problema, se considerará 4 elementos finitos
(unidimensionales), entonces tomaremos las alturas de 750,250,250,250
para estos cálculos.

obtenemos los anchos de cada

elemento unidimensional:

b1= 750 mm

b2= 417 mm

b3= 250 mm

b4= 83 mm

nodos GDL le b Ae
e primer segundo
1 2 (mm) (mm) (mm2)
nodos nodo
1 1 2 Q1 Q2 750 750 90000
2 2 3 Q2 Q3 250 417 50040
3 3 4 Q3 Q4 250 250 30000
4 4 5 Q4 Q5 250 83 9960

3. GRADOS DE LIBERTAD NODALES:

Luego obtendremos el vector desplazamiento

Q1

[]
Q2
Qi= Q3 mm
Q4
Q5

Donde Q1=0 ya que la placa esta empotrada y los demás desplazamientos

son incógnitas

4. VECTOR CARGA:

Analizando las fuerzas en cada elemento finito (en N):

ɣ ( AxL)1
F 11= + R1=2647.688+ R1
2

ɣ ( AxL )1
F 12= + P=52647.688
2

ɣ ( AxL )2
F 22= =490.7047
2

ɣ ( AxL )2
F 23= =490.7047
2
 ɣ ( AxL )3
F 33= =294.1875
2

ɣ ( AxL )3
F 34= =294.1875
2

ɣ ( AxL )4
F 44= =97.67025
2

ɣ ( AxL )4
F 45= =97.67025
2

Ahora analizaremos las fuerzas para todo el cuerpo (en N):

F 1=F11=2647.687+ R1

F 2=F12 + F 22=53138.392

F 3=F23 + F 33=784.892

F 4=F 34 + F 44 =391.857

F 5=F54=97.6702

Entonces el vector carga será:

2647.687+ R 1

[ ]
53138.392
F i= 784.892
391.857
97.6702

5. MATRIZ DE RIGIDEZ:

Ahora veremos la matriz de rigidez global que está determinado en la

siguiente ecuación
 1 −1 0 0 0 0 0 0 0 0 0 0 0 0 0

K ij =
AE
( )
l
(1)

[ −1 1 0 0 0
0
0
0
] [
0 0 0 0+
0 0 0 0
0 0 0 0
AE
( )
l
(2) 0 1 −1 0 0
0 −1 1 0 0 +
0 0
0 0
0 0 0
0 0 0
] [
AE
( )
l
(3) 0
0
0
0
0 0 0 0

]
0 1 −1 0 +
0 −1 1 0
0 0 0 0
AE
l(

Reemplazando los valores obtenidos anteriormente:

240 −240 0 0 0
5
K ij =10 ∗ 0
0
0
[
−240 640.32 −400.32
−400.32 640.32
0
0
−240
0
0
−24 0
0
0
319.68 −79.68
−79.68 79.68
N
mm
]
6. EUACION DE RIGIDEZ Y CONDICIONES DE CONTORNO:

De la ecuación de rigidez está determinada por:

F i=K ij∗Q ij

Remplazando nuestros valores en la ecuación anterior:

2647.687+ R1 240 −240 0 0 0 0

[ ] [
53138.392

391.857
97.6702
5
−240 640.32 −400.32
784.892 =10 ∗ 0
0
0
−400.32 640.32
0
0
−240
0
0
−240
0
0
Q2
Q3
319.68 −79.68 Q4
−79.68 79.68 Q5
][ ]
Resolviendo las ecuaciones:

Q 1=0

Q2=2.267∗10−3 mm

Q 3=2.299∗10−3 mm

Q4 =2.319∗10−3 mm
Q 5=2.332∗10−3 mm

R1=−5 7060.5 N

7. ESFUERZOS
Se aplicará la siguiente ecuación:

E e q
e
σ =
x
L ( )
[ −1 1 ] 1
q2 []
Reemplazando, se obtiene lo siguiente:

2∗10 5
σ1= ( ) [ −1 1 ] 0 ∗10−3=0.6045 N 2
[ ]
750 2.267 mm

2∗10 5
σ 2= ( ) [−1 1 ] 2.267 ∗10−3=0.0256 N 2
[ ]
250 2.299 mm

2∗105
σ3= ( )[−1 1 ] 2.299 ∗10−3 =0.016 N 2
[ ]
250 2.319 mm

2∗105
σ 4= ( )[ −1 1 ] 2.319 ∗10−3=0.0104 N 2
[ ]
250 2.332 mm

8. RESULTADOS
R1=−57060.5 N

N
σ 1 =0.6045
m m2

N
σ 2=0.0256
mm2

N
σ 3 =0.0 16
m m2

N
σ 4=0.0 104
m m2
CAP. 2 DIAGRAMA DE FLUJO
 CAP. 3 DIGITACION Y EJECUCION
clc
clear all
format short
b0=1000
bn=0
t=120
h=1500
n=4
E=200000
y=0.00007845
P=50000

le=zeros(n,1);
ho=zeros(n,1);
bo=zeros(n,1);
b=zeros(n,1);
a=zeros(n,1);
Fe=zeros(n+1,1);
bo(1)=b0;
ho(1)=h;

for i=1:n
if n>=i
le(i)=input('Ingrese la longitud del elemento finito(mm)= ');
b(i)=0.5*((bo(i)*(ho(i)-le(i)))/ho(i)+bo(i));
a(i)=b(i)*t;
ho(i+1)=ho(i)-le(i);
bo(i+1)=(bo(i)*(ho(i)-le(i)))/ho(i);
else
le(i)=ho(i);
b(i)=(bn+bo(i))/2;
a(i)=b(i)*t;
end
end

disp('bases(mm)=')
disp(b')
disp('longitudes(mm)=')
disp(le')
disp('Areas(mm^2)=')
disp(a')
for i=1:n
Fe(i)=y*a(i)*le(i)/2;
end

for i=1:n+1
if i==1
F(i)=Fe(i);
elseif i==n+1
F(i)=Fe(i-1);
else
F(i)=Fe(i-1)+Fe(i);
end
end

F(2)=F(2)+P;
disp('el vector fuerza(N)= ')
disp(F')

k=zeros(n+1);
for i=1:n
x=zeros(n+1);
x(i,i)=1;x(i+1,i)=-1;x(i,i+1)=-1;x(i+1,i+1)=1;
k=k+((a(i)*E)/(le(i)))*x;
end

disp('la matriz rigidez es (N/mm)')

disp(k')
Q=inv(k(2:n+1,2:n+1))*((F(2:n+1))');
Q=[0;Q];
disp('los desplazamientos son(mm)=')
disp(Q')
for i=1:n
e(i)=(E/(le(i)))*[-1 1]*[Q(i);Q(i+1)];
end

disp('los esfuerzos son(N/mm^2)= ')

disp(e');
R1=k*Q-F;
disp('la fuerza de reaccion(en N)=')
disp(R1(1,1)')
CAP 4. CONCLUSIONES DEL TRABAJO REALIZADO

 Se puede ver que, si se usan más nodos, se podría obtener valores

más exactos así el error seria aún menor.
 Se observa que las deformaciones son pequeñas, y salieron positivas
en el sentido asumido así que los esfuerzos son de tracción.
 Se recomendaría hacer otro tipo de mallado en igual cantidad de
particiones para observar si se presenta otro error.
 También se obtuvo error por los cálculos al aproximar.
 Si se hace el cálculo de reacción por estática, la reacción saldría
57060.5 igual que nuestra respuesta, entonces el método si es muy
efectivo en este caso.
 Al usar el programa Matlab mobile para usar el código, la respuesta
es muy similar a la obtenida de forma manual, se obtiene algo de
error en unas variables ya que se usa el format short.