Scribd
Cargar
31 vistas15 páginas

MC516 Lab2

El informe presenta un análisis de tracción simple unidimensional considerando el efecto de temperatura, utilizando el método de elementos finitos. Se realizan cálculos para dos configuraciones diferentes de elementos (3 y 4 elementos), incluyendo la obtención de matrices de fuerzas y de rigidez. Finalmente, se discuten los resultados y se ofrecen conclusiones y recomendaciones sobre el análisis realizado.

Cargado por

joce2005pinto
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
31 vistas15 páginas

MC516 Lab2

El informe presenta un análisis de tracción simple unidimensional considerando el efecto de temperatura, utilizando el método de elementos finitos. Se realizan cálculos para dos configuraciones diferentes de elementos (3 y 4 elementos), incluyendo la obtención de matrices de fuerzas y de rigidez. Finalmente, se discuten los resultados y se ofrecen conclusiones y recomendaciones sobre el análisis realizado.

Cargado por

joce2005pinto
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

UNIVERSIDAD NACIONAL DE

INGENIERÍA
FACULTAD DE INGENIERÍA MECÁNICA

CÁLCULO POR ELEMENTOS FINITOS


INFORME Nº2
TRACCIÓN SIMPLE UNIDIMENSIONAL – 1D CON
EFECTO DE TEMPERATURA

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

Profesor:
Miguel Ángel Ampuero Suárez
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 8
3.1. Compilación para 3 elementos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.2. Compilación para 4 elementos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

4. Conclusiones 12

5. Conclusiones 13

6. Recomendaciones 14

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
δT = 5 C
C=1
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 2, 3 400 75000
3 3 4 3, 4 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

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

Obteniendo la matriz debido al efecto de la temperatura:


   
1 5 −6 −1 −1
θ = 3 · 10 · 25000 · 11 · 10 · s = 2062500
1 1
   
2 5 −6 −1 −1
θ = 3 · 10 · 75000 · 11 · 10 · s = 1237500
1 1
   
−1 −1
θ3 = 3 · 105 · 25000 · 11 · 10−6 · s = 412500
1 1

θ1 = −2062500 N θ2 = 825000 N θ3 = 825000 N θ4 = 412500 N


Siendo la suma de vectores anteriores el vector carga global:
   
F1 −2060540 + R1
F2   828136 
F= F3  = 
  N
826568 
F4 422892
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
    
828136 150 −56,25 0 Q2
826568 = 106 × −56,25 75 −18,75 Q3 
422892 0 −18,75 18,75 Q4
   
Q2 22161,02
Q3  = 10−6 × 44373,65
Q4 66927,89
 
0
Q2
−2060540 + R1 = 106 × 93,75 −93,75 0 0 
   
Q3 
Q4

4
R1 = −17055, 625
Finalizando con los esfuerzos:
3 × 105
 
0 2
σ1 = ( ) × (−1 1) × 10−6 − 3 × 105 × 11 × 10−6 × 5 = 0,120765 N/mm
400 22161,02

3 × 105
 
22161,02 2
σ2 = ( ) × (−1 1) × 10−6 − 3 × 105 × 11 × 10−6 × 5 = 0,1594725 N/mm
400 44373,65
3 × 105
 
44373,65 2
σ3 = ( ) × (−1 1) × 10−6 − 3 × 105 × 11 × 10−6 × 5 = 0,41568 N/mm
400 66927,89

2.2. Analisis para 4 elementos

P = 10, 000 N
E = 3 × 105 MPa
3
ρ = 7,84 × 10−5 N/mm
δT = 5 C
C=1
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 2, 3 300 93750
3 3 4 3, 4 300 56250
4 4 5 4, 5 300 18750

Cuadro 2.2: Tabla de L y A




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

5
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
 
1543,5 + R1

 2646 

F =
 1764 
 N
 882 
10220,5

Obteniendo la matriz debido al efecto de la temperatura:


   
−1 −1
θ1 = 3 · 105 · 131250 · 11 · 10−6 · s = 2165625
1 1
   
−1 −1
θ2 = 3 · 105 · 93750 · 11 · 10−6 · s = 1546875
1 1
   
−1 −1
θ3 = 3 · 105 · 56250 · 11 · 10−6 · s = 928125
1 1
   
−1 −1
θ4 = 3 · 105 · 18750 · 11 · 10−6 · s = 309375
1 1

θ1 = −2165625 N θ2 = 618750 N θ3 = 618750 N θ4 = 618750 N θ5 = 309375 N

Siendo la suma de vectores anteriores el vector carga global:


   
F1 −2164081,5 + R1
F2   621396 
   
F = F3  = 
   620514 N

F4   619632 
F5 319595,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

6
 
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
    
621396 225 −93,75 0 0 Q2
 620514  6
 −93,75 150 −56,25 0  Q3 
 619632  = 10 ×  0
    
−56,25 75 −18,75 Q4 
319595,5 0 0 −18,75 18,75 Q5
   
Q2 16618,19
Q3 
  = 10−6 × 33255,43
 
Q4  49952,81
Q5 66997,9
 
0
 Q2 
 
−2164081,5 + R1 = 106 × 131,25 −131,25

0 0 0 Q3 

Q4 
Q5
R1 = −17055,9375
Finalizando con los esfuerzos:
3 × 105
 
0 2
σ1 = ( ) × (−1 1) × 10−6 − 3 × 105 × 11 × 10−6 × 5 = 0,11819 N/mm
400 16618,19

3 × 105
 
16618,19 2
σ2 = ( ) × (−1 1) × 10−6 − 3 × 105 × 11 × 10−6 × 5 = 0,13724 N/mm
400 33255,43
3 × 105
 
33255,43 2
σ3 = ( ) × (−1 1) × 10−6 − 3 × 105 × 11 × 10−6 × 5 = 0,19738 N/mm
400 49952,81
3 × 105
 
49952,81 2
σ4 = ( ) × (−1 1) × 10−6 − 3 × 105 × 11 × 10−6 × 5 = 0,54509 N/mm
400 66997,9

7
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;
Alpha = 11 ∗ (10( − 6));
delT = 264;
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;
T heT = zeros(N, 2);
f ori = 1 : N
f orj = 1 : 2
T heT (i, j) = (A(1, i) ∗ Alpha ∗ E ∗ delT ) ∗ (−1)( j);
end
end
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

8
F (1, i + 1) = f (i, 2) + f (i + 1, 1);
end
F (1, 1) = f (1, 1) + T heT (1, 1);
F (1, N + 1) = f (N, 2) + P + T heT (N, 2);
k = zeros(N + 1, N + 1);
K = cell(1, N );
KK = zeros(N + 1, N + 1);
f or i = 1 : N
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′ ) − E ∗ Alpha ∗ delT ;
end
disp(′ W ′ );
disp(W );
disp(′ Elvectorbase :′ );
disp(B);
disp(′ Elvectorarea :′ );
disp(A);
disp(′ Elvectordef ormacionporcalor :′ );
disp(T heT );
disp(′ Elvectorcargatotal :′ );
disp(F );
disp(′ Elvectordesplazamientototal :′ );
disp(Q);
disp(′ Elvectoresf uerzo :′ );

9
disp(tetha);
disp(′ ElvalordeR1 :′ );
disp(R1);

3.1. Compilación para 3 elementos

10
3.2. Compilación para 4 elementos

11
Conclusiones
Resultados para el primer caso:

Cálculos efectuados Resultados MATLAB %error


Q2(mm) 0,02216102 0,0222 0,1758944309
Q3(mm) 0,04437365 0,0444 0,0593820882
Q4(mm) 0,06692789 0,0669 0,0416717156
σ1 (N/mm2) 0.120765 0,1208 0,028981907
σ2 (N/mm2) 0,1594725 0,1595 0,0227714984
σ3 (N/mm2) 0,41568 0,4157 0,0048113934
R1(N) 17055,625 17060 0,02564478312

Resultados para el segundo caso:

Cálculos efectuados Resultados MATLAB %error


Q2(mm) 0,01661819 0,0166 0,1094583706
Q3(mm) 0,03325543 0,0333 0,1340232257
Q4(mm) 0,04995281 0,05 0,09446916
Q5(mm) 0,0669979 0,067 0,0031344266
σ1 (N/mm2) 0,11819 0,1182 0,0084609527
σ2 (N/mm2) 0,13724 0,1373 0,0437190324
σ3 (N/mm2) 0,19738 0,1974 0,0101327389
σ4 (N/mm2) 0,54509 0,5451 0,0018345594
R1(N) 17055.9375 17060 0.0238130129

12
Conclusiones
Se puede conlcuir en este laboratorio que se llevó acabo que el codigo de MATLAB funciona de
manera correcta, pues como se puede ver en el cuadro comparativo vemos que lo teórico con respecto
a lo programado no hay errores pronunciados pues a las justas llega al 0,1 por ciento de error como
maximo. Se analizo los desplazamientos globalos,esfuerzos y la fuerza que se genera al estar empotrado
en la parte de arriba sienod todos los valores muy cercanos siendo despreciable el error. Los errores que
pueden darse son debido a que manualmente no se llegan a trabajar con tantos decimales como si lo
haria el MATLAB haciendolo una herramienta muy útil al momento de querer precisión y rapidez.

El laboratorio fue muy similar al primero, con la diferencia que aqui se consideró la temperatura y
eso influyó más en los resultados de los desplazamientos y no tanto en los esfuerzos.

Se observa también que existe una mayor diferencia en los cálculos efectuados de R1 con respecto a
la obtenida por MATLAB.

13
Recomendaciones
Se recomienda tener cuidado en la parte de programación con respecto a la parte del manejo de
matrices que puede ser la más engorriosa. Aparte que en este caso no se puede hacer simulación con

ANSYS debido a que este no puede simular deformaciones por calor.

14

También podría gustarte