UNIVERSIDAD NACIONAL DE

INGENIERÍA
FACULTAD DE INGENIERÍA MECÁNICA

CÁLCULO POR ELEMENTOS FINITOSi

INFORME Nº1
TRACCIÓN UNIDIMENSIONAL

Estudiantes:
Pinto Quispe, José Manuel 20222141D
Custodio Medrano, Yerik Joel 20220225F
Contreras Cornejo, Emanuel Contreras 20222150C

Profesor:
Ampuero Suárez Miguel Ángel

Lima, 2025
Índice general

1. Problema 2

2. Calculos 3
2.1. Analisis para 3 elementos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2. Analisis para 4 elementos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

3. Codigo en Matlab 7
3.1. Compilación para 3 elementos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3.2. Compilación para 4 elementos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

4. Conclusiones 10

5. Recomendaciones 11

1
Problema

2
Calculos
2.1. Analisis para 3 elementos

P = 10, 000 N
E = 3 × 105 MPa
3
ρ = 7,84 × 10−5 N/mm
b1 = 833,3 mm, b2 = 500 mm, b3 = 166,7 mm

No. (1) (2) GDL L (mm) A (mm²)

1 1 2 1, 2 400 125000
2 2 3 1, 2 400 75000
3 3 4 2, 3 400 25000

Cuadro 2.1: Tabla de L y A

 
0
Q2 
Q=
Q3 


Q4

Obteniendo la matriz de fuerzas:

F11 = 1960 + R1 , F1 = 1960 + R1

F21 = 1960
F22 = 1176 F2 = 3136
F32 = 1176
F33 = 392 F3 = 1568

3
 F43 = 392 + 10000 F5 = 10392
 
1960 + R1
 3136 
F =
 1568 
 N
10392

Hallando K:
   
−1 1 0 0 0 0 0 0
1 −1 0 0 + 56,25 × 106 × 0 1 −1 0
K = 93,75 × 106 × 
 

0 0 0 0 0 −1 1 0
0 0 0 0 0 0 0 0
 
0 0 0 0
0 0 0 0
+ 18,75 × 106 × 


0 0 1 −1
0 0 −1 1
 
93,75 −93,75 0 0
−93,75 150 −56,25 0 
K = 106 ×  
 0 −56,25 75 −18,75
0 0 −18,75 18,75

Calcularemos Q y R1:
F = KQ
    
3136 150 −56,25 0 Q2
 1568  = 106 × −56,25 75 −18,75 Q3 
10392 0 −18,75 18,75 Q4
   
Q2 161,024
Q3  = 10−6 × 373,646
Q4 427,886
 
0

6

Q2 
1960 + R1 = 10 × 93,75 −93,75 0 0  
Q3 
Q4
R1 = −17056
Por equilibrio: R1 = P + f, f = p × Am × t
1000
R1 = 10000 + 7,84 × 10−5 × ( ) × 1200 × 150 = 17056
2

Finalizando con los esfuerzos:

3 × 105
 
0 3
σ1 = ( ) × (−1 1) × 10−6 = 0,1207 N/mm
400 161,024

3 × 105
 
161,024 3
σ2 = ( ) × (−1 1) × 10−6 = 0,1594 N/mm
400 373,646
3 × 105
 
373,646 3
σ3 = ( ) × (−1 1) × 10−6 = 0,4156 N/mm
400 427,886

4
2.2. Analisis para 4 elementos

P = 10, 000 N
E = 3 × 105 MPa
3
ρ = 7,84 × 10−5 N/mm
b1 = 675 mm, b2 = 625 mm, b3 = 375 mm, b4 = 125 mm

No. (1) (2) GDL L (mm) A (mm²)

1 1 2 1, 2 300 131250
2 2 3 1, 2 300 93750
3 3 4 2, 3 300 56250
4 4 5 3, 4 300 18750

Cuadro 2.2: Tabla de L y A


0
Q2 
 
Q3 
Q= 
Q4 
Q5

Obteniendo la matriz de fuerzas:

F11 = 1543,5 + R1 , F1 = 1543,5 + R1

F21 = 1543,5
F22 = 1102,5 F2 = 2646
F32 = 1102,5
F33 = 661,5 F3 = 1764
F43 = 661,5
F44 = 220,5 F4 = 882
F54 = 220,5 + 10000 F5 = 10220,5

5
  
1543,5 + R1

 2646 

F =
 1764 
 N
 882 
10220,5

Hallando K:
   
−1 1 0 0 0 0 0 0 0 0
 1 −1 0 0 0 0 1 −1 0 0
K = 131, 25 × 106 ×  6
   
0 0  + 93, 75 × 10 × 0
0 0 0  −1 1 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
   
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
+ 56, 25 × 106 ×  6
   
0 0 1 −1 0
 + 18, 75 × 10 × 0
 0 0 0 0
0 0 −1 1 0 0 0 0 1 −1
0 0 0 0 0 0 0 0 −1 1
 
131,25 −131,25 0 0 0
−131,25 225 −93,75 0 0 
K = 106 × 
 
 0 −93,75 150 −56,25 0  
 0 0 −56,25 75 −18,75
0 0 0 −18,75 18,75

Calcularemos Q y R1:
F = KQ
    
2646 225 −93,75 0 0 Q2
 1764 
  = 106 × 
−93,75 150 −56,25 0  Q3 
  
 882   0 −56,25 75 −18,75 Q4 
10220,5 0 0 −18,75 18,75 Q5
   
Q2 118,19
Q3  −6
255,433
  = 10 ×  
Q4  452,811
Q5 997,904
 
0
Q2 
1543,5 + R1 = 106 × 131,25 −131,25 0 0 0 


 
Q3 

Q4 
Q5
R1 = −17056
Por equilibrio:R1 = P + f, f = p × Am × t
1000
R1 = 10000 + 7,84 × 10−5 × ( ) × 1200 × 150 = 17056
2

Finalizando con los esfuerzos:

3 × 105
 
0 3
σ1 = ( ) × (−1 1) × 10−6 = 0,11819 N/mm
300 118,19

3 × 105
 
118,19 3
σ2 = ( ) × (−1 1) × 10−6 = 0,137243 N/mm
300 255,433
3 × 105
 
255,433 3
σ3 = ( ) × (−1 1) × 10−6 = 0,197378 N/mm
300 452,811
3 × 105
 
452,811 3
σ4 = ( ) × (−1 1) × 10−6 = 0,545093 N/mm
300 997,904

6
Codigo en Matlab
El código que se utilizó para el laboratorio fue el siguiente:

clear all;
clc;
E=3*(105 );
p = 7,84 ∗ (10( − 5));
t = 150;
w1 = 1000;
w2 = 0;
L = 1200;
P = 10000;
N = input(′ IN GRESAR LA CAN T IDAD DE ELEM EN T OS CON LOS QU E SE T RABAJARA :′
);
W 1 = zeros(1, N + 1);
W = zeros(1, N + 1);
f or i = 1 : N + 1
W 1(1, i) = (w2/2) + (i − 1) ∗ (w1 − w2)/(2 ∗ N );
end
f or i = 1 : N + 1
W (1, i) = W 1(1, N + 2 − i);
end
W = 2 ∗ W;
B = zeros(1, N );
f orı = 1 : N
B(1, i) = (W (1, i) + W (1, i + 1))/2;
end
A = zeros(1, 3);
f or i = 1 : N
A(1, i) = B(1, i) ∗ t;
end
LV = zeros(1, N );
f or i = 1 : N
LV (1, i) = L/N ;
end
R1 = 0;
f = zeros(N, 2);
f or i = 1 : N
f or j = 1 : 2
f (i, j) = (A(1, i) ∗ LV (1, i) ∗ p)/2;
end
end
F = zeros(1, N + 1);
f or i = 1 : N − 1
F (1, i + 1) = f (i, 2) + f (i + 1, 1);
end
F (1, 1) = f (1, 1);
F (1, N + 1) = f (N, 2) + P ;
k = zeros(N + 1, N + 1);
K = cell(1, N );
KK = zeros(N + 1, N + 1);
f or i = 1 : N

7
KK(i, i) = 1;
KK(i + 1, i) = −1;
KK(i + 1, i + 1) = 1;
KK(i, i + 1) = −1;
K1, i = KK;
KK = zeros(N + 1, N + 1);
end
f or i = 1 : N
K1, i = ((A(1, i) ∗ E)/LV...
(1, i)) ∗ K1, i;
end
f or i = 1 : N
k = k + K1, i;
end
F P RIM A = zeros(1, N );
f or i = 1 : N
F P RIM A(1, i) = F (1, i + 1);
end
kP RIM A = zeros(N, N );
f or i = 1 : N
f or j = 1 : N
kP RIM A(i, j) = k(i + 1, j + 1);
end
end
QP RIM A = ((kP RIM A( − 1)) ∗ (F P RIM A′ ));
Q = zeros(1, N + 1);
f or i = 1 : N
Q(1, i + 1) = QP RIM A(i, 1);
end
R1 = −f (1, 1) + k(1, 2) ∗ Q(1, 2);
tetha = zeros(1, N );
QP RIM A2 = zeros(1, 2);
f or i = 1 : N
QP RIM A2(1, 1) = Q(1, i);
QP RIM A2(1, 2) = Q(1, i + 1);
tetha(1, i) = (E/LV (1, i)) ∗ [−1, 1] ∗ (QP RIM A2′ );
end
disp(′ W ′ );
disp(W );
disp(′ El vector base :′ );
disp(B);
disp(′ El vector carga total :′ );
disp(F );
disp(′ El vector desplazamiento total :′ );
disp(Q);
disp(′ El vector esf uerzo :′ );
disp(tetha);

8
3.1. Compilación para 3 elementos

3.2. Compilación para 4 elementos

9
Conclusiones
Resultados para el primer caso:

Cálculos efectuados Resultados MATLAB %error

Q2(mm) 0,000161024 0,000161 0,01490461049
Q3(mm) 0,000373646 0,0003736 0,01231111801
Q4(mm) 0,000927886 0,0009279 0,001508806039
σ1 (N/mm2) 0,1207 0,1208 0,08285004143
σ2 (N/mm2) 0,1594 0,1595 0,06273525721
σ3 (N/mm2) 0,4156 0,4157 0,02406159769

Resultados para el segundo caso:

Cálculos efectuados Resultados MATLAB %error

Q2(mm) 0,00011819 0,0001182 0,008460952703
Q3(mm) 0,000255433 0,0002554 0,0129192391
Q4(mm) 0,000452811 0,0004528 0,002429269607
Q5(mm) 0,000997904 0,0009979 0,000400840161
σ1 (N/mm2) 0,11819 0,1182 0,008460952703
σ2 (N/mm2) 0,137243 0,1372 0,0313312883
σ3 (N/mm2) 0,197378 0,1974 0,01114612571
σ4 (N/mm2) 0,545093 0,5451 0,001284184534

Se puede concluir que en este laboratorio se llevó acabo de manera exitosa apreciándose que existe
poco error en la comparación de los resultados de desplazamiento y esfuerzos; el error se debe en mayor
medida al redondeo del programa MATLAB ya que las reacciones (R1) resultaron ser idénticas.

El laboratorio permitió comprender la relación directa entre las fuerzas aplicadas y los desplazamientos
resultantes en un elemento sometido a tracción.

10
Recomendaciones
Se recomendarı́a hacer uso del programa Ansys para simular los ensayos y tener una experiencia más
visible.
Verificar cuidadosamente las dimensiones de la probeta antes del ensayo para asegurar la exactitud
en el cálculo de esfuerzos y deformaciones.
Se obtendrian error

11