Análisis de una Placa triangular de espesor
constante, mediante el método de
elementos finitos, con tres elementos de
diferentes dimensiones.
Practica Nro-1
Calculo por elementos finitos
MC-516
¨ CÁLCULO POR ELEMENTOS FINITOS-MC 516¨
¨ PRÁCTICA NRO-1¨
APELLIDOS Y NOMBRES: RAMIREZ CASTAÑEDA, EDWARD
CÓDIGO: 20150170C
SECCIÓN: ¨F¨
DOCENTE: CUEVA PACHECO RONALD
FECHA: 30-09-2019
2019-II
MC-516 1
MC-516 2
1. Enunciado del problema .-Cuerpo real
Dada la siguiente placa triangular, de espesor constante
t 150 mm
Datos.-
Cargas externas
P 60000 N
ᵞ (peso específico) 8 gr-f/cm^3
Material:
E 3.*10^5 N/mm^2
Calcular:
a) .- Calcular los esfuerzos en determinadas secciones de la barra de
sección rectangular.
b) .-La reaccionen los apoyos.
MC-516 3
2. Modelado del cuerpo Real.-
El modelado del cuerpo es con una malla de rectángulos de diferentes
dimensiones. Se considerará 3 Elementos Finitos
Para facilitar los cálculos se tendrán longitudes de 500, 500,1000 mm
Los elementos para el análisis son rectángulos de un ancho de la base (b)
MC-516 4
Para cada elemento se calculará el ancho de su base, que vendría a ser el
promedio
b1 750 mm
=(500+1000)/2
b2 =(500+250)/2 375 mm
b3 =(250+0)/2 125 mm
Cuadro de conectividad
GRADOS DE
NODOS L(LONGITUD) A(AREA)
ELEMENTO LIBERTAD
mm mm^2
(1) (2) 1 2
1 1 2 1 2 1000 112500
2 2 3 2 3 500 56250
3 3 4 3 4 500 18750
Se procedió a calcular el área de cada elemento finito
Área =bixt
MC-516 5
Elemento finito maestro: en posición vertical
3. Grados de libertad Nodales.- (Vector Desplazamiento)
3.1.- Local.-
MC-516 6
3.2.-Global.-
Donde Q1 es cero por que esta
empotrada en la base, los demás
desplazamientos se tendremos
que calcularlos.
0
𝑄2
𝑄 = [𝑄3]
𝑄4
4. Vector carga.-
4.1.-Local.-
MC-516 7
4.2.-Global.-
Analizando las fuerzas involucradas en cada elemento finito (Global)
4i.-Componentes del vector de carga.- Local:
MC-516 8
Dónde: PE es el la acción del peso específico
𝐹 𝑒 𝑆=1,2 = 𝑃𝑥𝑙𝑒/2 − 𝑃𝐸(𝐴1𝑥𝐿1)/2 − 𝑃𝑆
Para nuestros 3 elementos
𝐹11 = 𝑅1 − 𝑃𝐸(𝐴1𝑥𝐿1)/2
𝐹1 2 = −𝑃𝐸(𝐴1𝑥𝑙1)/2 − 𝑃(60𝐾𝑁)
𝐹 2 2 = −𝑃𝐸(𝐴2𝑥𝑙2)/2
𝐹 2 3 = −𝑃𝐸(𝐴2𝑥𝑙2)/2
𝐹 3 3 = −𝑃𝐸(𝐴3𝑥𝑙3)/2
𝐹 3 4 = −𝑃𝐸(𝐴3𝑥𝑙3)/2
Donde:
𝐹1 𝐹11 𝑅1 − 𝑓(𝐴1𝑥𝑙1)/2
1 2 −𝑓(𝐴1𝑥𝑙1)/2 − 𝑃(60𝐾𝑁) − 𝑓(𝐴2𝑥𝑙2)/2
𝐹2 𝐹 +𝐹 2
𝐹 = [ ] = 22 3 =
𝐹3 𝐹 3+𝐹 3 −𝑓(𝐴2𝑥𝑙2)/2 − 𝑓(𝐴3𝑥𝑙3)/2
𝐹4 [ 𝐹 3
] [ −𝑓(𝐴3𝑥𝑙3)/2 ]
4
5. Matrices de rigidez
5.1.-Local.-
1 −1 1
𝑅 𝑒 𝑟𝑠 = (𝐸𝑥𝐴/𝑙)𝑒 𝑥 [ ]
−1 1 2
5.2.-Global.-
∈
𝐾𝑖𝑗 = ∑ [𝑅 𝑒 𝑟𝑠 ] 𝑠=𝑖
𝑒=1 𝑟=𝑗
Desarrollo:
1 −1 0 0 0 0 0 0
𝐾𝑖𝑗 = (𝐸𝑥𝐴/𝑙)1 𝑥 [−1 1 0 0] + (𝐸𝑥𝐴/𝑙)2 𝑥 [ 0 1 −1 0]
0 0 0 0 0 −1 1 0
0 0 0 0 0 0 0 0
0 0 0 0
+(𝐸𝑥𝐴/𝑙) 𝑥 [0
3 0 0 0 ]
0 0 1 −1
0 0 −1 1
MC-516 9
Kij =
[33750000 -33750000 0 0
-33750000 67500000 -33750000 0
0 -33750000 45000000 -11250000
0 0 -11250000 11250000]
6. Ecuación de Rigidez y condiciones de contorno.-
𝐹𝑖 = 𝐾𝑖𝑗 . 𝑄
Desarrollo:
𝑅1 − 𝑓(𝐴1𝑥𝑙1)/2 𝐾11 𝐾12 𝐾13 𝐾14 0
−𝑓(𝐴1𝑥𝑙1)/2 − 𝑃(60𝐾𝑁) − 𝑓(𝐴2𝑥𝑙2)/2 𝐾21 𝐾22 𝐾23 𝐾24 𝑄2
=[ ] . [ ] … (𝛼)
−𝑓(𝐴2𝑥𝑙2)/2 − 𝑓(𝐴3𝑥𝑙3)/2 𝐾31 𝐾32 𝐾33 𝐾34 𝑄3
[ −𝑓(𝐴3𝑥𝑙3)/2 ] 𝐾41 𝐾32 𝐾33 𝐾44 𝑄4
Los esfuerzos de contacto Superficial no existen, en este caso particular.
Todas las incógnitas, tanto en el primer miembro (R1) como en el segundo
miembro (Q2, Q3, Q4); deben resolverse, mediante un artificio, tomamos las
submatriz
−𝑓(𝐴1𝑥𝑙1)/2 − 𝑃(60𝐾𝑁) − 𝑓(𝐴2𝑥𝑙2)/2 𝐾22 𝐾23 𝐾24 𝑄2
[ −𝑓(𝐴2𝑥𝑙2)/2 − 𝑓(𝐴3𝑥𝑙3)/2 ] = [𝐾32 𝐾33 𝐾34 ] . [𝑄3]
−𝑓(𝐴3𝑥𝑙3)/2 𝐾32 𝐾33 𝐾44 𝑄4
Resolviendo:
Q2 =-0.0020 mm
Q3 =-0.0021 mm
Q4 =-0.0021 mm
MC-516 10
Con ello reemplazamos en la ecuación (α) y obtenemos F1
F1 = (6.7357).10^4 N
Donde despejando
R1=F1+ (PE*A1*l1/2)
R1 =7.177210^4 N (Respt. b)
7. Esfuerzos (En cada elemento finito)
En el método numérico ya desarrollado:
𝑄𝑖
𝛿 𝑒 = 𝑒𝑠𝑓 𝑒 = (𝐸/𝑙)𝑒 . [−1 1] [ ]
𝑄𝑖+1
El esfuerzo en cada elemento:
esf1 =-0.5987 [N/mm^2]
esf2 =-0.0327 [N/mm^2]
esf3 =-0.0196 [N/mm^2]
(Respt. a)
El signo de los esfuerzos, según el elemento maestro escogido, son esfuerzos
de compresión.
MC-516 11
Uso de software matlab
Diagrama de Flujo
MC-516 12
% Elementos finitos-Practica 1-Placa triangular-Código en MATLAB
clc
clear all
%t=150 mm
t=input('Ingrese el espesor de la placa en milímetros =')
disp('t')
%P= 60000 N
P=input ('Carga externa de la placa en Newton =')
disp('P')
%peso específico= 8 gr-f/cm^3
PE1=input('Ingrese el peso específico de la placa en gr-f/cm^3=')
PE=PE1*9.81*10^-6
disp('PE')
disp('PE esta en Newtons')
%E = 3*10^5 N/mm^2
E=input('Ingrese Modulo de la elasticidad del material en N/mm^2 =')
B=input('Ingrese la base de la placa triangular en milímetros =')
H=input('Ingrese la altura de la placa triangular en milímetros =')
l1=input('Ingrese longitud del primer elemento finito en milímetros =')
l2=input('Ingrese longitud del segundo elemento finito en milímetros =')
l3=input('Ingrese longitud del tercer elemento finito en milímetros =')
x1=(H-l1)*B/H;
x2=l3*B/H;
b1=(B+x1)/2;
b2=(x1+x2)/2;
b3=(x2+0)/2;
A1=b1*t; A2=b2*t; A3=b3*t
Q1=0; % empotrado en el pizo no hay desplazamiento
% desplazamientos Q2;Q3;Q4;
%F1= R1-f(A1*l1)/2
F2=-PE*A1*l1/2 -P -PE*A2*l2/2
F3=-PE*A2*l2/2-PE*A3*l3/2
F4=-PE*A3*l3/2
%Matrices de rigidez
Kij=(E*A1/l1)*[1 -1 0 0;-1 1 0 0;0 0 0 0;0 0 0 0]+(E*A2/l2)*[0 0 0 0;0 1 -1 0;0 -1
1 0;0 0 0 0]+(E*A3/l3)*[0 0 0 0;0 0 0 0;0 0 1 -1;0 0 -1 1]
SubF=[F2;F3;F4]
SubK=Kij([2,3,4],[2,3,4])
SubQ=(inv(SubK))*SubF
Q2=SubQ(1,:)
Q3=SubQ(2,:)
Q4=SubQ(3,:)
F=Kij*[Q1;Q2;Q3;Q4]
F1=F(1,:)
R1=F1+(PE*A1*l1/2)
esf1 = (E/l1)*[-1 1]*[Q1;Q2]
esf2 = (E/l2)*[-1 1]*[Q2;Q3]
esf3 = (E/l3)*[-1 1]*[Q3;Q4]
MC-516 13
EN MATLAB
MC-516 14
Con la ayuda de MATLAB realizamos varias pruebas variando el peso específico
del material, hasta que este sea cero.
PE[gr-f/cm^3] esf1 [N/mm^2] esf2[N/mm^2] esf3[N/mm^2] R1 [N]
8 -5.987E-01 -3.270E-02 -1.960E-02 7.177E+04
7 -0.5906 -0.0286 -0.0172 7.03E+04
6 -0.5824 -0.0245 -0.0147 68829
5 -0.5742 -0.0204 -0.0123 6.74E+04
4 -0.566 -0.0164 -0.0098 6.59E+04
3 -0.5579 -0.0123 -0.0074 6.44E+04
2 -0.5497 -0.0082 -0.0049 62943.0000
1 -0.5415 -0.0041 -0.0025 6.15E+04
0 -0.5333 0 0 6.00E+04
PE(peso específico) Vs Esfuerzos
0.000E+00
0 1 2 3 4 5 6 7 8 9
-1.000E-01
-2.000E-01
-3.000E-01
-4.000E-01
-5.000E-01
-6.000E-01
-7.000E-01
PE [gr-f/cm^3]
esf1 [N/mm^2] esf2[N/mm^2] esf3[N/mm^2]
MC-516 15
PE(peso específico) vs R1 (reaccion)
7.400E+04
7.200E+04
7.000E+04
R1[Newton]
6.800E+04
6.600E+04
6.400E+04
6.200E+04
6.000E+04
5.800E+04
0 1 2 3 4 5 6 7 8 9
PE[gr-f/cm^3]
Del modelado en SolidWorks
0.744
N/mm2
0.595
N/mm2
0.446
N/mm2
0.297
N/mm2
0.148
N/mm2
0
MC-516 16
CONCLUSIONES
Aplicando MEF (0.5987N/mm2) con 3 elementos aún hay un gran margen de error
comparando con el esfuerzo máximo hallado en SolidWorks (0.744 N/mm2 ) cerca
de un 20%, recordemos que para la ingeniería el valor aproximado debe estar cerca del
5 a 3% de error (a medida que aumentemos la cantidad de elementos este valor se
acerca más al del software).
Con ayuda del MATLAB, tabulamos una pequeña tabla PE (peso específico) vs
Esfuerzos en cada elemento, con esta podemos afirmar que el peso del cuerpo influye
de manera casi insignificante en los cálculos, por qué el error está cerca del 5%, muy
admisible en ingeniería.
PE[gr-f/cm^3] % error-esf1 % error-esf2 % error-esf3 %R1
8 1.353E+00 1.254E+01 1.224E+01 2.051E+00
7 1.388E+00 1.434E+01 1.453E+01 2.092E+00
6 1.408E+00 1.673E+01 1.633E+01 2.137E+00
5 1.428E+00 1.961E+01 2.033E+01 2.185E+00
4 1.431E+00 2.500E+01 2.449E+01 2.234E+00
3 1.470E+00 3.333E+01 3.378E+01 2.284E+00
2 1.492E+00 5.000E+01 4.898E+01 2.339E+00
1 1.514E+00 1.000E+02 1.000E+02 2.393E+00
Acerca del problema de la exactitud lo solucionamos con ayuda del MATLAB
recordemos que al operar a mano hay pérdida de cifras decimales, centésimas,
milésimas que al final influyen en el resultado final, estos problemas se
solucionan usando Matlab con una precisión casi exacta de 1.86*10^-14
aproximadamente 0.
El valor de los esfuerzos y la reacción varia de manera lineal con respecto al
peso específico, en solo caso de la placa es triangular.
MC-516 17
