Laboratorio 5 - Calculo Por Elementos finitos-FIM
Laboratorio 5 - Calculo Por Elementos finitos-FIM
Laboratorio 5 - Calculo Por Elementos finitos-FIM
2019-1
QUINTA PRÁCTICA CALIFICADA
PROBLEMA
Según la figura que plantea el problema, esta se puede aproximar a una viga con 2
cargas distribuidas. En la cual podemos apreciar un punto fijo y 3 puntos libres por lo
cual procederemos a hacer los cálculos
𝒒𝒂𝒊𝒓𝒆 = 𝟎. 𝟓 ∗ 𝝆𝒂𝒊𝒓𝒆 ∗ 𝒗𝟐
𝑭𝒂𝒊𝒓𝒆 = 𝟏𝟎𝟎 ∗ 𝒒𝒂𝒊𝒓𝒆
𝑞𝑎𝑖𝑟𝑒 = 10.6575 𝑁/ 𝑚
𝐹𝑎𝑖𝑟𝑒 = 1065.75 𝑁
Calculamos la fuerza que ejerce la marea
𝑭𝒎𝒂𝒓𝒆𝒂 = 𝜸𝒂𝒈𝒖𝒂 ∗ 𝑨 ∗ 𝒉𝒄
𝐴 = 10 ∗ 6 ∗ 𝜋
ℎ𝑐 = 5 𝑚
𝐹𝑚𝑎𝑟𝑒𝑎 = 1582.5569 ∗ 6 𝑁
Calculamos la posición donde ejerce la fuerza de la marea
𝑰𝒄
𝒍𝒑 = 𝒍𝒄 +
𝒍𝒄 ∗ 𝑨
𝑙𝑐 = 5 𝑚
𝐴 = 10 ∗ 𝐷 ∗ 𝜋
1
𝐼𝑐 = ∗ 𝐷 ∗ 103
12
𝑙𝑝 = 7.6526 𝑚
Matriz de esfuerzos
0 0
−𝐹𝑎𝑖𝑟𝑒 −1065.75
𝐹=[ ] = [−1582.5569 ∗ 6]
−𝐹𝑚𝑎𝑟𝑒𝑎
𝑅1 𝑅1
2. GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)
𝑄1 𝑄1
𝑄2 𝑄2
𝑄=[ ]=[ ]
𝑄3 𝑄3
𝑄4 0
3. MATRIZ DE RIGIDEZ
𝑭𝒊 = 𝑲𝒊𝒇 ∗ 𝑸𝒇
𝑄1 = 0.5432 𝑚𝑚
𝑄2 = 0.3234 𝑚𝑚
𝑄3 = 0.1273 𝑚𝑚
𝑅1 = 6324.7 𝑁
6∗𝐸∗𝐼
𝜏 =𝛼∗( ) ∗ (2 ∗ 𝑄1 + 𝑙𝑒 ∗ 𝑄2 − 2 ∗ 𝑄3 + 𝑙𝑒 ∗ 𝑄4 )
𝐴 ∗ 𝑙𝑒3
𝜏 = 425.432 𝑃𝑎
𝜏 < 𝜎𝑝𝑒𝑟𝑚𝑖𝑠𝑖𝑏𝑙𝑒
6∗𝐸∗𝐼
𝑉=( ) ∗ (2 ∗ 𝑄1 + 𝑙𝑒 ∗ 𝑄2 − 2 ∗ 𝑄3 + 𝑙𝑒 ∗ 𝑄4 )
𝑙𝑒3
𝑉 = 6704.23 𝑁
6. RESULTADOS
𝑸𝟏 = 𝟎. 𝟓𝟒𝟑𝟐 𝒎𝒎
𝑸𝟐 = 𝟎. 𝟑𝟐𝟑𝟒 𝒎𝒎
𝑸𝟑 = 𝟎. 𝟏𝟐𝟕𝟑 𝒎𝒎
𝑹𝟏 = 𝟔𝟑𝟐𝟒. 𝟕 𝑵
𝝉 = 𝟒𝟐𝟓. 𝟒𝟑𝟐 𝑷𝒂
𝑽 = 𝟔𝟕𝟎𝟒. 𝟐𝟑 𝑵
DIAGRAMA DE FLUJO Y FUNCION EN MATLAB
DIGITACIÓN POR MATLAB
PROGRAMA EN MATLAB
clc;
format long;
n=input('Ingrese Número de Elementos Finitos:');
e1=input('Espesor de las alas(mm):');
e2=input('Espesor del alma(mm):');
l1=input('Longitud de las alas(mm):');
L=input('Ingrese Longitud de la Viga(mm):');
E=input('Modulo de Elasticidad(N/mm2):');
yp=input('Ingrese Peso Específico(N/mm3):');
pe=input('Carga Distribuida Externa(N/mm):');
disp('MOMENTOS DE INERCIA')
for i=1:(n/2)
d(i)=(4*(L*(i-1)/n)/15+100+4*(L*i/n)/15+100)/2;
I(i)=2*(l1*e1^3/12)+e2*(d(i)-2*e1)^3/12+2*(d(i)/2-e1/2)^2*l1*e1;
end
for i=((n/2)+1):n
d(i)=(900-4*(L*(i-1)/n)/15+900-4*(L*i/n)/15)/2;
I(i)=2*(l1*e1^3/12)+e2*(d(i)-2*e1)^3/12+2*(d(i)/2-e1/2)^2*l1*e1;
end
disp(I)
disp('MATRIZ DE RIGIDEZ ESTRUCTURAL K')
k=zeros(2*(n+1),2*(n+1));
for i=1:n
l=L/n;
ke(:,:,i)=E*I(i)/l*[12 6*l -12 6*l;
6*l 4*l*l -6*l 2*l*l;
-12 -6*l 12 -6*l;
6*l 2*l*l -6*l 4*l*l];
gl1=2*i-1;
gl2=gl1+1;
gl3=2*(i+1)-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
k(gl,gl)=k(gl,gl)+ke(:,:,i);
end
disp(k)
disp('FUERZAS DEBIDO AL PESO DEL MATERIAL')
for i=1:n
A(i)=l1*e1*2+e2*(d(i)-2*e1);
p(i)=-yp*A(i);
end
w=zeros(1,2*(n+1));
for i=1:n
l=L/n;
we(:,:,i)=[p(i)*l/2 p(i)*l^2/12 p(i)*l/2 -p(i)*l^2/12];
gl1=2*i-1;
gl2=gl1+1;
gl3=2*(i+1)-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
w(1,gl)=w(1,gl)+we(:,:,i);
end
wt=w';
disp(wt)
disp('FUERZAS DEBIDO A LA CARGA DISTRIBUIDA')
c=zeros(1,2*(n+1));
for i=2:3
l=L/n;
ce(:,:,i)=[pe*l/2 pe*l^2/12 pe*l/2 -pe*l^2/12];
gl1=2*i-1;
gl2=gl1+1;
gl3=2*(i+1)-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
c(1,gl)=c(1,gl)+ce(:,:,i);
end
ct=c';
disp(ct)
disp('FUERZA TOTAL')
f=ct+wt;
disp(f)
disp('DESPLAZAMIENTOS')
disp('Q=')
kf=k(3:8,3:8);
ff=f(3:8,1);
qf=inv(kf)*ff;
Q=[0;0;qf;0;0];
disp(Q)
disp('ESFUERZOS LONGITUDINALES(N/mm2)')
y=input('Ingrese punto genérico a analizar:');
z=-1;
es1=zeros(n,1);
for i=1:n
gl1=i*2-1;
gl2=gl1+1;
gl3=(i+1)*2-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
q=Q(gl);
es1(i)=-E*y/l^2*[6*z (3*z-1)*l -6*z (3*z+1)*l]*q;
end
disp('z=-1')
disp(es1)
z=1;
es2=zeros(n,1);
for i=1:n
gl1=i*2-1;
gl2=gl1+1;
gl3=(i+1)*2-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
q=Q(gl);
es2(i)=-E*y/l^2*[6*z (3*z-1)*l -6*z (3*z+1)*l]*q;
end
disp('z=1')
disp(es2)
z=0;
es0=zeros(n,1);
for i=1:n
gl1=i*2-1;
gl2=gl1+1;
gl3=(i+1)*2-1;
gl4=gl3+1;
gl=[gl1 gl2 gl3 gl4];
q=Q(gl);
es0(i)=-E*y/l^2*[6*z (3*z-1)*l -6*z (3*z+1)*l]*q;
end
disp('z=0')
disp(es0)
EJECUCIÓN DEL PROGRAMA
Módulo de Elasticidad(N/mm2):2.1e5
MOMENTOS DE INERCIA
1.0e+008 *
0.33741516666667
1.0e+017 *
Columns 1 through 4
0 0 -0.00000693413896 -0.00260030211000
0 0 0.00260030211000 0.65007552750000
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0
Columns 5 through 8
0 0 0 0
0 0 0 0
-0.00000693413896 0.00260030211000 0 0
-0.00260030211000 0.65007552750000 0 0
0 0 -0.00000113371496 -0.00042514311000
0 0 0.00042514311000 0.10628577750000
Columns 9 through 10
0 0
0 0
0 0
0 0
0 0
0 0
-0.00000113371496 0.00042514311000
-0.00042514311000 0.10628577750000
0.00000113371496 -0.00042514311000
-0.00042514311000 0.21257155500000
FUERZAS DEBIDO AL PESO DEL MATERIAL
1.0e+004 *
-0.01994250375000
-2.49281296875000
-0.05423213250000
-1.79339062500000
-0.06857925750000
-0.05423213250000
1.79339062500000
-0.01994250375000
2.49281296875000
-1875
-234375
-3750
-1875
234375
0
FUERZAS DEBIDO A LA CARGA DISTRIBUIDA
-1875
-234375
-3750
-1875
234375
FUERZA TOTAL
1.0e+005 *
-0.00199425037500
-0.24928129687500
-0.02417321325000
-2.52308906250000
-0.04435792575000
-0.02417321325000
2.52308906250000
-0.00199425037500
0.24928129687500
DESPLAZAMIENTOS
Q= 1.0e-006 *
-0.06887011717571
-0.00007462641465
-0.10005353989168
-0.00000000000000
-0.06887011717571
0.00007462641465
ESFUERZOS LONGITUDINALES(N/mm2)
z=-1
1.0e-005 *
CONCLUSIONES
BIBLIOGRAFÍA
https://slideplayer.es/slide/1727186/
http://cybertesis.uach.cl/tesis/uach/2014/bmfciv473d/doc/bmfciv473d.pdf
https://es.slideshare.net/rosamarisolcrisantoflores9/accion-del-viento-en-las-
estructuras{
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies: