Laboratorio 5 - Calculo Por Elementos finitos-FIM

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 23

UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA MECÁNICA

LABORATORIO Nº5: OFFSSHORE TURBINA EOLICA

CURSO: CALCULO POR ELEMENTOS FINITOS- MC516


SECCION: “A” GRUPO: ‘3’
PROFESOR: ABREGÚ LEANDRO EDWIN ASENCIÓN

ALUMNOS: VIGURIA BALDEON LUIS GUSTAVO 20151139B


BUENDIA CHIPANA CARLOS ALBERTO 20156222I
GARCIA VILLAFRANCA DIEGO FLORENCIO 20152094B
GARCIA PILLACA ALEX JESUS 20167515I
CASTILLO HUARCAYA BRYAN CARLOS 20130045I
MOLINA MOSCOSO DENNIS AHMED 20112050D
CASAZOLA NINANTAY RAY BRIAN 20132650G

2019-1
QUINTA PRÁCTICA CALIFICADA

PROBLEMA

¿Es posible determinar la geometría de un monopile para una turbina eólica


offshore sometida al viento corriente y ola marina?

𝛒𝐚𝐢𝐫𝐞 1.225 kg/m3


𝛒𝐚𝐠𝐮𝐚 1027 kg/m3
V viento 17.4 m/s
E 206.8 GPa
D 6m
𝛔𝐩𝐞𝐫𝐦𝐢𝐬𝐢𝐛𝐥𝐞 540 Pa
SOLUCIÓN DEL PROBLEMA

1. MODELADO DEL CUERPO REAL

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

Calculamos la fuerza que ejerce el aire

𝒒𝒂𝒊𝒓𝒆 = 𝟎. 𝟓 ∗ 𝝆𝒂𝒊𝒓𝒆 ∗ 𝒗𝟐
𝑭𝒂𝒊𝒓𝒆 = 𝟏𝟎𝟎 ∗ 𝒒𝒂𝒊𝒓𝒆
𝑞𝑎𝑖𝑟𝑒 = 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 𝑚

Reemplazaremos las cargas distribuidas por sus fuerzas equivalentes luego


procederemos a calcular por elementos finitos sus deformaciones, esfuerzos, fuerza
cortante, momento flector.

Matriz de esfuerzos
0 0
−𝐹𝑎𝑖𝑟𝑒 −1065.75
𝐹=[ ] = [−1582.5569 ∗ 6]
−𝐹𝑚𝑎𝑟𝑒𝑎
𝑅1 𝑅1
2. GRADOS DE LIBERTAD NODALES (Vector Desplazamiento)

A través del grafico se muestran los grados de libertad


nodales globales.

Se puede considerar en el nodo 4 como punto fijo pues


están situados bajo tierra por ende el matriz
desplazamiento será:

𝑄1 𝑄1
𝑄2 𝑄2
𝑄=[ ]=[ ]
𝑄3 𝑄3
𝑄4 0

3. MATRIZ DE RIGIDEZ

Y la matriz de rigidez local está dado por:


𝟏𝟐 𝟔 ∗ 𝒍𝒆 −𝟏𝟐 𝟔 ∗ 𝒍𝒆
𝑬∗𝑰 𝟔 ∗ 𝒍𝒆 𝟒 ∗ 𝒍𝟐𝒆 −𝟔 ∗ 𝒍𝒆 𝟐 ∗ 𝒍𝟐𝒆
𝑲𝒆 = ( 𝟑 ) ∗
𝒍𝒆 −𝟏𝟐 −𝟔 ∗ 𝒍𝒆 𝟏𝟐 −𝟔 ∗ 𝒍𝒆
[𝟔 ∗ 𝒍𝒆 𝟐 ∗ 𝒍𝟐𝒆 −𝟔 ∗ 𝒍𝒆 𝟒 ∗ 𝒍𝟐𝒆 ]
Reemplazando se obtendrá:
12 6 ∗ 110 −12 6 ∗ 110
206.8 ∗ 109 ∗ 64 6 ∗ 110 4 ∗ 1102 −6 ∗ 110 2 ∗ 1102
𝐾𝑒 = ( )∗[ ]
1103 ∗ 64 −12 −6 ∗ 110 12 −6 ∗ 110
6 ∗ 110 2 ∗ 1102 −6 ∗ 110 4 ∗ 1102
4. ECUACIONES DE RIGIDEZ Y CONDICIONES DE CONTORNO

La ecuación de rigidez está determinada por la siguiente ecuación:

𝑭𝒊 = 𝑲𝒊𝒇 ∗ 𝑸𝒇

Lo que con nuestros valores calculados tenemos:

0 12 6 ∗ 110 −12 6 ∗ 110 𝑄1


−1065.75 4 6 ∗ 110 4 ∗ 1102 −6 ∗ 110 2 ∗ 1102 𝑄2
[−1582.5569 ∗ 6] = 2427.69 ∗ 6 ∗ [ ]∗[ ]
12 −6 ∗ 110 12 −6 ∗ 110 𝑄3
𝑅1 6 ∗ 110 2 ∗ 1102 −6 ∗ 110 4 ∗ 1102 0

0 = 2427.69 ∗ 64 ∗ (12 ∗ 𝑄1 + 660 ∗ 𝑄2 − 12 ∗ 𝑄3 )

−1065.75 = 2427.69 ∗ 64 ∗ (660 ∗ 𝑄1 + 48400 ∗ 𝑄2 − 660 ∗ 𝑄3 )

−1582.5569 ∗ 6 = 2427.69 ∗ 64 ∗ (−12 ∗ 𝑄1 − 660 ∗ 𝑄2 + 12 ∗ 𝑄3 )

𝑅1 = 2427.69 ∗ 64 ∗ (660 ∗ 𝑄1 + 24200 ∗ 𝑄2 − 660 ∗ 𝑄3 )

Resolviendo la matriz podremos obtener los valores de 𝑄2 , 𝑄3 , 𝑄4 , 𝑄5 , 𝑄6 , 𝑅1

𝑄1 = 0.5432 𝑚𝑚

𝑄2 = 0.3234 𝑚𝑚

𝑄3 = 0.1273 𝑚𝑚

𝑅1 = 6324.7 𝑁

5. DISEÑO DE ESFUERZOS FUERZAS Y MOMENTO FLECTOR

Para calcular el esfuerzo cortante máximo aplicaremos la siguiente ecuación

6∗𝐸∗𝐼
𝜏 =𝛼∗( ) ∗ (2 ∗ 𝑄1 + 𝑙𝑒 ∗ 𝑄2 − 2 ∗ 𝑄3 + 𝑙𝑒 ∗ 𝑄4 )
𝐴 ∗ 𝑙𝑒3

𝜏 = 425.432 𝑃𝑎

𝜏 < 𝜎𝑝𝑒𝑟𝑚𝑖𝑠𝑖𝑏𝑙𝑒

Para hallar la fuerza cortante se usará:

6∗𝐸∗𝐼
𝑉=( ) ∗ (2 ∗ 𝑄1 + 𝑙𝑒 ∗ 𝑄2 − 2 ∗ 𝑄3 + 𝑙𝑒 ∗ 𝑄4 )
𝑙𝑒3

𝑉 = 6704.23 𝑁
6. RESULTADOS

Finalmente, los resultados son mostrados en la siguiente tabla:

𝑸𝟏 = 𝟎. 𝟓𝟒𝟑𝟐 𝒎𝒎
𝑸𝟐 = 𝟎. 𝟑𝟐𝟑𝟒 𝒎𝒎
𝑸𝟑 = 𝟎. 𝟏𝟐𝟕𝟑 𝒎𝒎
𝑹𝟏 = 𝟔𝟑𝟐𝟒. 𝟕 𝑵
𝝉 = 𝟒𝟐𝟓. 𝟒𝟑𝟐 𝑷𝒂
𝑽 = 𝟔𝟕𝟎𝟒. 𝟐𝟑 𝑵
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

Ingrese Número de Elementos Finitos: 4

Espesor de las alas (mm):13

Espesor del alma (mm):25

Longitud de las alas (mm):100

Ingrese Longitud de la Viga (mm):3000

Módulo de Elasticidad(N/mm2):2.1e5

Ingrese Peso Específico(N/mm3):76.518e-6

Carga Distribuida Externa(N/mm):-5

MOMENTOS DE INERCIA

1.0e+008 *

0.33741516666667 2.06373183333333 2.06373183333333

0.33741516666667

MATRIZ DE RIGIDEZ ESTRUCTURAL K

1.0e+017 *

Columns 1 through 4

0.00000113371496 0.00042514311000 -0.00000113371496 0.00042514311000

0.00042514311000 0.21257155500000 -0.00042514311000 0.10628577750000

-0.00000113371496 -0.00042514311000 0.00000806785392 0.00217515900000

0.00042514311000 0.10628577750000 0.00217515900000 1.51272261000000

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.00001386827792 0 -0.00000693413896 0.00260030211000

0 2.60030211000000 -0.00260030211000 0.65007552750000

-0.00000693413896 -0.00260030211000 0.00000806785392 -0.00217515900000

0.00260030211000 0.65007552750000 -0.00217515900000 1.51272261000000

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

FUERZAS DEBIDO A LA CARGA DISTRIBUIDA

-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)

Ingrese punto genérico a analizar: 50

z=-1

1.0e-005 *
CONCLUSIONES

 se puede apreciar que la fuerza ejercida por la marea depende de la sección de


la planta eólica por lo que si se conoce su esfuerzo permisible y el factor de
seguridad se puede determinar el diámetro máximo que debe tener para que no
colapse con el tiempo
 el diseño de la sección transversal influye mucho en el cálculo de la fuerza
ejercida por la marea pues permite desplazar la fuerza en sentido contrario
anulándola entre sí
 las deformaciones no son significantes con respecto al tamaño del molino por lo
que no presenta ningún riesgo de ruptura

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{

 DESIGN AND FINITE ELEMENT ANALYSIS OF AN OCEAN CURRENT


TURBINE BLADE, Nicholas S. Asseff, The College Engineering and Computer
Science in Partial Fulfillment of the Requirements for the Degree of Master of
Science, August 200
 http://data.irestal.com/files/files/2012030204152933979.pdf
 https://www.empresascarbone.com/pdf/ficha-tecnica-del-acero-inoxidable.pdf
ANEXOS
Turbina eólica Offshore:

Tipos de base de turbinas:


Partes de la turbina eolica:

Partes del soporte de la turbina:


Teniendo en cuenta el viento y las olas marinas, como se muestra en la figura calcularemos la
geometria del monopile.

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

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:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy