Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
MODIFICACIÓN DE CONDICIONES INICIALES A UNA BARRA CON
ÁREA VARIABLE
Se presentan las condiciones iniciales del problema, donde mediante el código MATLAB se
obtienen variables como esfuerzos y desplazamientos, luego se realizarán modificaciones para
ver las diferencias entre unos y otros, estas modificaciones pueden ser de densidad, morfología
de la barra, módulo de Young, etc.
Esquema de la barra:
El esquema 1 corresponde al problema original y el 2 a el modelado para una placa con dos
elementos.
Código MATLAB:
1 clc
2 clear all
3
4 %%%%asignación de variables
5 P=100; %lb En el centro
6 t=1; %in espesor de placa
7 E=60e6; %psi modulo de Young ESTOS SON LOS
8 rho=0.2836; %lb/in^3 Densidad VALORES QUE SE IRÁN
9 L=24; %longitud de la barra
10 B=12; %ancho inicial MODIFICANDO.
11 b=12; %ancho final
12
13 e=2; %elementos
14 n=e+1; %nodos.
15 l=L/e; %diferencial de longitud
16
17 %vector de anchos=#de nodos x=[x1 x2 ... xn]=[B x2 ... b]
18 x=zeros(1,n);
19 for i=1:n
20 x(i)=B-(B-b)*(i-1)*l/L;
21 end
22 x';
1
Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
23 %%vector de areas= x.*t= espesor * ancho
24
25 for i=1:e
26 A(i)=((x(i)+x(i+1))/2)*t;
27 end
28 A';
29
30 %matriz de rigidez global 3x3
31 K=zeros(n,n);
32 %matriz rigidez de cada elemento
33 for i=1:e
34 k1=E*A(i)/l*[1 -1; -1 1]; %matriz de rigidez de cada elemento
35 for j=1:2
36 for p=1:2
37 K(j+(i-1),p+(i-1))=K(j+(i-1),p+(i-1))+k1(j,p);
38 end
39 end
40 end
41 K;
42
43 %vector de fuerzas global
44 F=zeros(n,1);
45 for i=1:e
46 f=A(i)*l*rho/2*[1; 1]; %vector de fuerzas de cada elemento
47 for j=1:2
48 F(j+(i-1))=F(j+(i-1))+f(j);
49 end
50 end
51 F %imprimir F
52
53 %localización de la fuerza
54 Px=e/2+1; %%%fuerza puntual localizada en la mitad
55 F(Px)=F(Px)+P;
56 F;
57
58 %matriz reducida
59 Kr=K;
60 Kr(1,:)=[]; %elimina fila 1 del elemento Kr
61 Kr(:,1)=[]; %elimina columna 1 del elemento Fr
62 Kr;
63
64 %Vector F reducido
65 Fr=F;
66 Fr(1,:)=[]; %elimina fila 1 del elemento Fr
67 Fr;
68 %
69 %desplazamientos
70 U=zeros(n,1);
71 Ur=inv(Kr)*Fr;
72 U(2:n)=Ur;
73 U
74 %
2
Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
75 %%%Esfuerzos
76 for i=1:e
77 u=[U(i); U(i+1)];
78 d(i)=E/l*[-1 1]*u;
79 end
80 d
81
82 %vector de reaciones en nodos
83 R=K*U-F;
84 R1=R(1);
85 sum(F);
86
87 %discretización
88 xx=0:l:L;
89
90 figure
91 subplot(2,1,1);
92 plot(xx,U,'o-') %desplazamientos
93 xlabel('Longitud de barra (in)')
94 ylabel('Desplazamiento (in)')
95
96 subplot(2,1,2);
97 plot(d,'o-') %esfuerzos
98 xlabel('Elemento')
99 ylabel('Esfuerzo (psi)')
PROBLEMA INICIAL:
• Carga aplicada (P)=100 (lb)
• Espesor de placa(t)=1(in)
• Módulo de Young(E)=30*106(psi)
• Densidad (ρ)=0.2836(lb/in3)
• Longitud de barra(L)= 24 (in)
• Ancho inicial (B)=6(in)
• Ancho final(b)=3(in)
Obteniendo el vector fuerza y desplazamiento
8.93 0
𝐹 = (15.31) [𝑙𝑏] ; 𝑢 = (0.9272 ∗ 10−5 ) [𝑖𝑛]
6.38 0.9953 ∗ 10−5
𝑙𝑏 𝑙𝑏
𝜎1 = 23.18 ( 2
) ; 𝜎2 = 1.70 ( 2 )
𝑖𝑛 𝑖𝑛
3
Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
Con sus respectivos gráficos
Del gráfico 1 se tiene: un patrón de desplazamiento que muestra cómo la placa se deforma en
respuesta a la carga aplicada y las condiciones de contorno. Los nodos cercanos a la carga
puntual tienen un mayor desplazamiento en comparación con los nodos lejanos.
Del gráfico 2: se identifica como responden los esfuerzos a medida que se avanza a través de los
elementos, siendo el esfuerzo máximo percibido en el elemento 1.
4
Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
VARIACIÓN 1: AUMENTO DE DENSIDAD Y ESPESOR.
Teóricamente se tiene: si se aumenta la densidad se obtiene un aumento de inercia (barra más
resistente a vibración y deformación) ya que aumenta la masa y con esto la rigidez, conllevando
a una distribución de esfuerzos mas uniforme. El aumento de espesor hace que se tenga una
mayor sección transversal, disminuyendo deformaciones y concentraciones de esfuerzo.
Se verá mediante el programa de MATLAB si se cumple lo que dice la teoría. Por ejemplo, si
se utiliza un acero con densidad de 0.5672 (lb/in3) y se ensancha el espesor de la placa a 2
(in), el resto de las variables se mantiene igual. Se tiene:
35.73 0
𝐹 = (61.26) [𝑙𝑏] ; 𝑢 = (0.7115 ∗ 10−5 ) [𝑖𝑛]
25.52 0.8477 ∗ 10−5
𝑙𝑏 𝑙𝑏
𝜎1 = 17.8 ( 2
) ; 𝜎2 = 3.4 ( 2 )
𝑖𝑛 𝑖𝑛
Observando lo obtenido se tiene: aumentó la fuerza resultante F (no coindice con teoría) lo
cual se puede deber a algún error o factor no considerado en el modelo, sin embargo, los valores
de σ1 y σ2 disminuye lo que lo hace una estructura más eficiente en cuanto a distribución de
esfuerzos, además los valores de desplazamiento disminuyen lo que es coherente con la teoría
que hace referencia a que se tendrá una barra más rígida.
5
Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
VARIACIÓN 2: AUMENTO DE CARGA APLICADA Y LARGO TOTAL DE
LA BARRA.
Teóricamente: el aumento de longitud implica mayor distancia entre los puntos de apoyo
(condiciones de contorno), aumentando la deformación de la barra y los esfuerzos, además el
aumento de la carga P presentaría una distribución de esfuerzos mas alta en la estructura.
Manteniendo todas las condiciones originales del problema, modificándose P=200 (lb) y L=48
(in). Se obtiene:
17.87 0
𝐹 = (30.63) [𝑙𝑏] ; 𝑢 = (0.37 ∗ 10−4 ) [𝑖𝑛]
12.76 0.39 ∗ 10−4
𝑙𝑏 𝑙𝑏
𝜎1 = 46.36 ( 2
) ; 𝜎2 = 3.4 ( 2 )
𝑖𝑛 𝑖𝑛
Observaciones: se comprueba todo lo preestablecido por la teoría, los cambios aplicados en la
variación 2 conllevan a una barra más flexible y experimentará una mayor concentración de
esfuerzos. aun así, se está hablando de desplazamientos casi imperceptibles al ojo humano de
orden 10-4, por ende, no habría complicaciones de derrumbe del sistema aun aumentando el
doble la carga aplicada.
6
Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
VARIACIÓN 3: MODIFICACIÓN DEL MODULO DE ELASTICIDAD Y ÁREA
CONSTANTE.
Teóricamente se tiene: El módulo de elasticidad (E) se refiere a la capacidad de un material
para resistir la deformación elástica bajo cierta carga aplicada, un (E) más alto conlleva a
menores desplazamientos respecto a la situación inicial, ya que la barra es más rígida y resiste
mejor a la carga aplicada. Por otra parte, mantener constante la el área transversal (A) significa
que la barra podrá resistir la carga con menos concentración de esfuerzos.
Los valores modificados quedan como E=60*106 psi, B=12 y b=12.
20.42 0
𝐹 = (40.83) [𝑙𝑏] ; 𝑢 = (0.27 ∗ 10−5 ) [𝑖𝑛]
20.42 0.30 ∗ 10−4
𝑙𝑏 𝑙𝑏
𝜎1 = 13.43 ( 2
) ; 𝜎2 = 1.7 ( 2 )
𝑖𝑛 𝑖𝑛
Observaciones: nuevamente los valores teóricos coinciden con los experimentales, volviendo una
barra más rígida y resistente a la deformación, así como una distribución de esfuerzo mucho
mejor sobre todo comparando el esfuerzo 𝜎1 en la condición inicial que era 23.18 (lb/in2), se
redujo considerablemente para este caso.
7
Elementos Finitos Tarea 06 Ignacio Muñoz Baeza
CONCLUSIONES Y COMENTARIOS
Las respectivas variaciones aplicadas en cada uno de los casos influyen notoriamente en los
esfuerzos y desplazamientos de la barra, los resultados experimentales en su mayoría poseen
concordancia con la teoría, sin embargo, puede haber factores externos no considerados en el
modelo como por ejemplo la temperatura ambiente, el sitio donde se va a trabajar, si la barra
puede estar cerca de elementos corrosivos, humedad, etc. por lo tanto siempre es necesario en
el mundo real garantizar esos entornos y con eso la integridad de la barra, pero aun así
mediante el uso de Elementos Finitos con MATLAB se puede hacer una muy buena
aproximación y simulación de lo que se tendrá en el análisis estructural.
La elección de la modificación a hacer a la barra dependerá de factores como los requisitos de
rendimiento, presupuesto, disponibilidad de materiales, durabilidad y mantenimiento, factores
de diseño y la misma viabilidad del diseño y fabricación, todo depende de las necesidades
específicas de cada proyecto.
Finalmente y siguiendo con lo último, se considera la variación numero 3 (modificación de
módulo de elasticidad y área constante) como la mejor y más viable, el aumento del módulo
de elasticidad genera una estructura más rígida y eficiente en términos de distribución de
esfuerzos, mantener el área constante permite controlar los costos de materiales y fabricación
de la barra (menos complejidad de diseño), además se tiene mejor comportamiento ante la
carga lo cual garantiza la seguridad y eficiencia estructural.
