Design of a small sustainable wind blade




<span id='_GoBack'></span>

=1. Introducción=

Para conseguir una descarbonización del planeta, una condición necesaria es pasar de una economía lineal a una economía circular. Uno de los objetivos es que los productos tengan una mayor vida útil. Y cuando éstos lleguen al fin de su vida, los residuos generados sean mínimos, buscando una solución para que estos residuos pasen a ser materia prima de otros nuevos productos.  Enfocándonos en la industria eólica, actualmente de un aerogenerador se recicla un importante porcentaje del mismo, cercano al 90%. Las palas de los aerogeneradores están fabricadas de composites formados por fibras (vidrio y carbono) y resinas termoestables (poliéster, viniléster o epoxi), que hacen que al final de su vida útil su reciclado sea costoso. Existen varias tecnologías de reciclaje con diferentes grados de desarrollo, como pueden ser trituración, pirolisis o solvólisis. Cada una de ellas cuenta con ventajas e inconvenientes, pero se está trabajando en obtener solución para poder alcanzar el 100% de reciclado de los aerogeneradores.

En este contexto, el trabajo que aquí se presenta pone el foco de la economía circular desde el diseño de las palas del aerogenerador, pensando en qué hacer con ellas al final de su vida útil, antes incluso de fabricarlas. Es decir, incluir el concepto de sostenibilidad desde el diseño inicial de la pala. Para ello se ha diseñado una pequeña pala de aerogenerador de pequeña potencia, sustituyendo la resina termoestable, material que nos dificulta el reciclado al final de la vida útil del aerogenerador, por una novedosa resina termoplástica. Al fabricar composites con resinas termoplásticas, se han conseguido recuperar las fibras y resinas mediante procesos más sencillos tecnológicamente y utilizando menores recursos energéticos a pequeña escala (1).

Esta modificación en el diseño de una pequeña pala eólica sostenible, hace que el volumen de residuo al final de la vida útil disminuya significativamente, y se recuperen la mayor parte de las materias primas usadas en la fabricación de la pala y puedan ser usado en la fabricación de nuevos composites.

=2. Diseño de la pala=

El primer paso ha sido realizar el diseño aerodinámico de una pala de 2m de longitud. Para ello se ha utilizado el software QBlade (2) para definir la pala a partir del perfil aerodinámico SG6043, desarrollado por Selig y Giguere (SG). Es un perfil que tiene un 5,1 % de altura máxima de la línea media, situado al 49,7 % del borde de ataque y con un espesor relativo del 10,0 %, situado al 32,1 % del borde de ataque. La raíz de la pala se ha definido de sección cilíndrica. Desde la raíz de la pala se ha interpolado el perfil SG6043 con otros dos perfiles intermedios hasta llegar a la zona aerodinámica de la pala, a partir del cual se ha usado únicamente el perfil SG definido. La <span id='cite-_Ref193698746'></span>[[#_Ref193698746|Figura  '''1''']] muestra los cuatro perfiles usados para el diseño de la pala.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Carnicero_712546528-image1-c.jpg|384px]] </div>

<div id="_Ref193698746" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figura  1 – Perfiles usados en el diseño de la pala a partir de perfil SG6043 (verde)</div>

Posteriormente se ha usado el software SolidWorks (3) para el diseño CAD de la pala, usando los perfiles de las diferentes secciones obtenidos con Qblade.

[[Image:Draft_Carnicero_712546528-image2-c.png|600px]]

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figura  2 – Pala diseñada con SolidWorks</div>

El siguiente paso ha sido diseñar los diferentes tejidos que va a llevar la pala. Para ello se ha utilizado el software ANSYS Workbench (4). Se decide diseñar la pala con fibra de carbono. Se usarán tejidos unidireccionales TFX – 322 g/m<sup>2,</sup> así como tejidos biaxiales a +- 45. La resina a utilizar es una resina termoplástica llamada Akelite (5), desarrollada y patentada por ITCP-CSIC. Puesto que en el momento no se disponía de las propiedades mecánicas provenientes de ensayo se usaron los valores estimados que se muestran en las Tablas 1 y 2.

{| style="width: 56%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Propiedad
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Valor
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Densidad [g/cm<sup>3</sup>]
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|1.49
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|E<sub>11</sub> [GPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|110.83
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|E<sub>22</sub> [GPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|5.07
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|E<sub>33</sub> [GPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|5.07
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|v<sub>12</sub>  
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.31
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|v<sub>13</sub>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0.31
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|G<sub>12</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|3344
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|G<sub>23</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|3005
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|G<sub>13</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|3344
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 1. Propiedades mecánicas estimadas para una lámina UD CFRP TFX – 322 g/m<sup>2</sup>.</div>

{| style="width: 56%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Propiedad
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Valor
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|X<sub>T</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|1979
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|Y<sub>T</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|26
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|Z<sub>T</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|26
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|X<sub>C</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|893
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|Y<sub>C</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|139
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|Z<sub>C</sub> [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|139
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|Shear XY [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|100
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|Shear YZ [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|50
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|Shear XZ [MPa]
|  style="border: 1pt solid black;text-align: center;vertical-align: bottom;"|100
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 2. Resistencias de una lámina UD CFRP TFX – 322 g/m<sup>2</sup>.</div>

Se diseña el laminado a partir de tres zonas de igual longitud, buscando mantener la simetría en cada zona y la continuidad del mayor número de capas posible. Se considera un espesor de capa de 0.3 mm. El laminado propuesto para la pala de CFRP/Akelite se muestra en la <span id='cite-_Ref159863682'></span>[[#_Ref159863682|Figura  3]]. Las capas ±45<sup>o</sup> se modelan usando dos capas UD del material con propiedades mostradas en las Tablas 2 y 3, de manera que se obtengan unas propiedades equivalentes al tejido biaxial ± 45 BC-50K-300 - 305 g/m<sup>2</sup>.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Carnicero_712546528-image3.png|456px]]  [[Image:Draft_Carnicero_712546528-image4-c.png|384px]] </div>

<span id='_Ref159863682'></span><div id="_Ref159863676" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figura  3 – Laminado para las palas de CFRP/Akelite</div>

==2.1. Casos de carga==

Se evalúan las cargas definidas en UNE-EN 61400-2 (6), según metodología de cargas simplificadas. Se utilizan los parámetros de diseño mostrados en la <span id='cite-_Ref166750018'></span>[[#_Ref166750018|Tabla 3]].

{| style="width: 82%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Parámetro 
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Símbolo
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Valor
|-
|  style="border: 1pt solid black;vertical-align: top;"|Masa pala
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''mB [kg]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|5.7
|-
|  style="border: 1pt solid black;vertical-align: top;"|Distancia entre CdG pala y CdG rotor
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''Rcog [m]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0,86
|-
|  style="border: 1pt solid black;vertical-align: top;"|Velocidad del rotor de diseño
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''n,design [rpm]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|250
|-
|  style="border: 1pt solid black;vertical-align: top;"|Velocidad del rotor máxima
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''n,max [rpm]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|350
|-
|  style="border: 1pt solid black;vertical-align: top;"|Par rotor de diseño
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''Qdesign (Nm)''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|164
|-
|  style="border: 1pt solid black;vertical-align: top;"|Número de palas
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''B (-)''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|3
|-
|  style="border: 1pt solid black;vertical-align: top;"|Gravedad
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''g (m/s<sup>2</sup>)''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|9,81
|-
|  style="border: 1pt solid black;vertical-align: top;"|Radio rotor
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''R  (m)''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|2,135
|-
|  style="border: 1pt solid black;vertical-align: top;"|Velocidad del viento de diseño
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''V,design (m/s)''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|10,5
|-
|  style="border: 1pt solid black;vertical-align: top;"|Distancia centro rotor a eje orientación
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''Lrt [m]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0,26
|-
|  style="border: 1pt solid black;vertical-align: top;"|Distancia centro rotor y primer cojinete
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''Lrb [m]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0,28
|-
|  style="border: 1pt solid black;vertical-align: top;"|Densidad del aire
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''Ρ [kg/m<sup>3</sup>]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|1,225
|-
|  style="border: 1pt solid black;vertical-align: top;"|Área de pala proyectada en plano perpendicular
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''Aproj,B [m<sup>2</sup>]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0,35
|-
|  style="border: 1pt solid black;vertical-align: top;"|Coeficiente de empuje
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''CT [-]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0,5
|-
|  style="border: 1pt solid black;vertical-align: top;"|Velocidad del viento media anual
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|''Vave [m/s]''
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|3,94<span id="fnc-1"></span><span style="text-align: center; font-size: 75%;">[[#fn-1|<sup>1</sup>]]</span>
|}


<div id="_Ref166750018" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 3. Parámetros de diseño de la pala</div>

Los casos de carga más desfavorables se recogen en la <span id='cite-_Ref192671839'></span>[[#_Ref192671839|Tabla 4]], sombreados en gris. El caso que produce un mayor momento flector es el C (error de orientación). A partir de éste se propone el ensayo con dos cargas de 500 N situadas a 0,66m y 1,44m de la raíz, que produce un momento flector y unas solicitaciones similares a los casos de carga más desfavorables.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;vertical-align: bottom;"|
|  colspan='6'  style="border: 1pt solid black;text-align: center;"|'''Cargas'''
|-
|  style="border: 1pt solid black;text-align: center;"|'''Hipótesis'''
|  style="border: 1pt solid black;text-align: center;"|'''F<sub>x </sub>[N]'''
|  style="border: 1pt solid black;text-align: center;"|'''F<sub>y</sub> [N]'''
|  style="border: 1pt solid black;text-align: center;"|'''F<sub>z </sub>[N]'''
|  style="border: 1pt solid black;text-align: center;"|'''M<sub>x </sub>[N·m]'''
|  style="border: 1pt solid black;text-align: center;"|'''M<sub>y</sub> [N·m]'''
|  style="border: 1pt solid black;text-align: center;"|'''M<sub>z</sub> [N·m]'''
|-
|  style="border: 1pt solid black;text-align: center;"|A: Operación normal
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|6720,73
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-150,8607
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-291,0045
|-
|  style="border: 1pt solid black;text-align: center;"|B: Orientación
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-174,880
|-
|  style="border: 1pt solid black;text-align: center;"|C: Error orientación
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-904,0493
|-
|  style="border: 1pt solid black;text-align: center;"|D: empuje máximo
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-476,7770
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|-
|  style="border: 1pt solid black;text-align: center;"|E: Velocidad rotación máxima
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|6586,319
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|-
|  style="border: 1pt solid black;text-align: center;"|F: cortocircuito
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-194,76
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|-
|  style="border: 1pt solid black;text-align: center;"|G: desconexión
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-125,62
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|-
|  style="border: 1pt solid black;text-align: center;"|H: carga extrema de viento
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|0
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|-679,8018
|-
|  style="border: 1pt solid black;text-align: center;"|I: Inmovilización
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|-
|  style="border: 1pt solid black;text-align: center;"|J: Transporte
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|N/A
|-
|  style="border: 1pt solid black;text-align: center;"|
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|-
|  style="border: 1pt solid black;text-align: center;"|'''Cargas de ensayo*'''
|  style="border: 1pt solid black;text-align: center;"|-
|  style="border: 1pt solid black;text-align: center;"|-1000
|  style="border: 1pt solid black;text-align: center;"|-
|  style="border: 1pt solid black;text-align: center;"|-
|  style="border: 1pt solid black;text-align: center;"|-
|  style="border: 1pt solid black;text-align: center;"|-
|}


<span style="text-align: center; font-size: 75%;">*2 cargas de 500 N situadas en dos secciones situadas a 0,66m y 1,44m de la longitud de la pala.</span>

<span id='_Ref192671839'></span><span id='_Ref166750109'></span><div id="_Ref166750103" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 4. Casos de carga</div>

==2.2. Modelo de elementos finitos==

Se crea un modelo de elementos finitos usando ANSYS Workbench, definiendo los laminados de la <span id='cite-_Ref159863682'></span>[[#_Ref159863682|Figura  '''3''']] usando el módulo ACP. Se añade la masa del adhesivo (0.7 kg) repartido en las aristas de unión entre la cara superior y la inferior (o extrados e intrados). La pala resultante tiene una masa de 4.23 kg, ligeramente inferior a la considerada en el cálculo de las cargas. Se utilizan 7410 elementos Shell cuadráticos, de tamaño medio 10 mm y verificando la calidad del mallado. Únicamente se tienen 11 elementos con calidad inferior a 0.7 y 23 con aspecto ratio superior a 2. El mallado se muestra en la <span id='cite-_Ref192672715'></span>[[#_Ref192672715|Figura  4]]. La pala se empotra por la raíz y se aplican las cargas definidas en la <span id='cite-_Ref192671839'></span>[[#_Ref192671839|Tabla 4]].

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;width: 51%;"|[[Image:Draft_Carnicero_712546528-image5.png|276px]] 
|  style="text-align: center;vertical-align: top;width: 48%;"|[[Image:Draft_Carnicero_712546528-image6.png|276px]] 
|}


<div id="_Ref192672715" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figura  4 – Imagen general mallado</div>

==2.3. Resultados de los casos de carga==

Se evalúan las deflexiones en cada una de las tres direcciones y los criterios de fallo de Hashin y Puck.

{| style="width: 82%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;"|'''Hipótesis'''
|  style="border: 1pt solid black;text-align: center;"|'''δx [mm]'''
|  style="border: 1pt solid black;text-align: center;"|'''δy [mm]'''
|  style="border: 1pt solid black;text-align: center;"|'''δz [mm]'''
|  style="border: 1pt solid black;text-align: center;"|'''SF (Hashin)'''
|  style="border: 1pt solid black;text-align: center;"|'''SF (Puck)'''
|-
|  style="border: 1pt solid black;text-align: center;"|'''A: Operación normal'''
|  style="border: 1pt solid black;text-align: center;"|-0.65
|  style="border: 1pt solid black;text-align: center;"|-47.25
|  style="border: 1pt solid black;text-align: center;"|9.71
|  style="border: 1pt solid black;text-align: center;"|4.73
|  style="border: 1pt solid black;text-align: center;"|4.73
|-
|  style="border: 1pt solid black;text-align: center;"|'''C: Error orientación'''
|  style="border: 1pt solid black;text-align: center;"|1
|  style="border: 1pt solid black;text-align: center;"|-139.03
|  style="border: 1pt solid black;text-align: center;"|21.63
|  style="border: 1pt solid black;text-align: center;"|1.90
|  style="border: 1pt solid black;text-align: center;"|1.90
|-
|  style="border: 1pt solid black;text-align: center;"|'''D: empuje máximo'''
|  style="border: 1pt solid black;text-align: center;"|0.36
|  style="border: 1pt solid black;text-align: center;"|-27.317
|  style="border: 1pt solid black;text-align: center;"|4.88
|  style="border: 1pt solid black;text-align: center;"|4.47
|  style="border: 1pt solid black;text-align: center;"|4.47
|-
|  style="border: 1pt solid black;text-align: center;"|'''J: Cargas de ensayo*'''
|  style="border: 1pt solid black;text-align: center;"|-1.85
|  style="border: 1pt solid black;text-align: center;"|-136.17
|  style="border: 1pt solid black;text-align: center;"|22.21
|  style="border: 1pt solid black;text-align: center;"|1.49
|  style="border: 1pt solid black;text-align: center;"|1.49
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 5 – Resultados de deflexión (δ) y factor de seguridad (SF) para cada caso de carga.</div>

Se verifica que la pala cumple los requisitos de solicitaciones estáticas y que el ensayo propuesto reproduce las condiciones más desfavorables.

==2.4. Análisis modal==

Para el análisis modal se impiden todos los grados de libertad en la raíz de la pala. Previamente al cálculo definitivo, se analizó la sensibilidad al mallado y se concluyó que con el mismo usado para el análisis estructural se obtienen resultados con suficiente precisión. Se calculan los 10 primeros modos usando el método de Lanczos.

{| style="width: 67%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Tamaño elemento (mm)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Frecuencia (Hz)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Modo vibración
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;vertical-align: top;"|10
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|12.90
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Flapwise
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|31.32
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Edgewise
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 6 – Resultados análisis modal</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|  colspan='2'  style="vertical-align: top;width: 100%;"|[[Image:Draft_Carnicero_712546528-image7.png|390px]] 
|-
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Carnicero_712546528-image8.png|318px]] 
|  style="vertical-align: top;width: 50%;"|[[Image:Draft_Carnicero_712546528-image9.png|318px]] 
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figura  5 – Formas de los modos de vibración: Flapwise (arriba), Edgewise (abajo)</div>

=3. Fabricación y Ensayo de la pala=

Se ha fabricado una pala de fibra de carbono y resina termoplástica Akelite, con el laminado de fibras de carbono definido en la <span id='cite-_Ref159863682'></span>[[#_Ref159863682|Figura  3]] del apartado anterior. Se ha realizado el proceso de infusión de las dos valvas de la pala y posteriormente se han pegado dichas valvas mediante adhesivo estructural. La zona cercana a la raíz de la pala contiene 13 capas de fibra de carbono con diferentes direcciones de 0, 45 y 90º, obteniendo un espesor de 3.54mm. La zona intermedia central de la pala está formada por 9 tejidos de fibra de carbono, con un espesor total de 2.34mm. Por último, la zona de la punta de la pala lleva 5 tejidos de fibra de carbono, y un espesor de 1.74mm.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Carnicero_712546528-image10-c.jpeg|600px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figura  6 – Pala de fibra de carbono fabricada con resina termoplástica Akelite</div>

El último paso ha sido el ensayo de la pala. Para ello se ha realizado el ensayo de propiedades de la pala según la norma IEC 61400-23, midiendo el peso, el centro de gravedad y las frecuencias naturales.

{| style="width: 85%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Peso

(kg)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Centro de gravedad (mm)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Frecuencia Flapwise (Hz)
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Frecuencia Edgewise (Hz)
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|3.48
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|807
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|12.50
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|24.14
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 7 – Resultados del ensayo de propiedades</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Carnicero_712546528-image11-c.jpeg|342px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figura  7 – Medida de las frecuencias naturales de la pala CFRP + Akelite </div>

=4.    Conclusiones=

En el presente artículo se ha mostrado el diseño de una pequeña pala eólica sostenible, definiendo desde la fase de diseño el uso de una resina termoplástica para el uso en la fabricación de palas eólicas.  Posteriormente se han hecho las simulaciones para calcular las cargas estáticas de ensayo y el análisis modal de las frecuencias naturales de la pala. El siguiente paso ha sido la fabricación de la pala según el diseño previamente establecido. Y finalmente se ha realizado el ensayo de propiedades de la pala para comprobar las posibles desviaciones del modelo.

Comparando los valores de las frecuencias naturales de la pala, se puede observar que la desviación de la frecuencia de flapwise está muy cercana a lo previsto en el diseño, con una desviación de 3.1%; mientras que para la frecuencia en edgewise el desvío ha sido mayor, con una desviación de 22.9%.  Obtenidos estos valores, se entiende que es necesario ajustar el modelo para disminuir este error en la frecuencia en edgewise, antes de realizar el ensayo estático de la pala, ya que las cargas de diseño han de ser ajustadas. Por otro lado, la desviación del centro de gravedad de la pala respecto al valor de diseño ha sido del 6.16%. Los resultados obtenidos son bastante prometedores, pero se necesita seguir dando pasos en el desarrollo de esta novedosa resina termoplástica Akelite para implementarla en la fabricación de palas eólicas de mayor envergadura.

=5.    Bibliografía=

1. ''Lessons learned from the sustainable recycling process of a wind blade made with a new thermoplastic resin. ''Carnicero R., Cano L., Cruz I. s.l. : IOP, Journal of Physics: Conference Series, 2024.

2. Marten, David. ''QBlade: A modern tool for the aero elastic simulation of wind turbines. ''2020.

3. [https://www.solidworks.com/. https://www.solidworks.com/.] [En línea]

4. [https://www.ansys.com/products/ansys-workbench. https://www.ansys.com/products/ansys-workbench.] [En línea]

5. Verdejo R., López-Manchado M.A. ''ES2906674B2. ''España, 2021.

6. AENOR. UNE-EN 61400-2: Aerogeneradores. Parte 2: Aerogeneradores pequeños.

----

<span id="fn-1"></span>([[#fnc-1|<sup>1</sup>]])  Se utilizan datos de ["https://globalwindatlas.info/en/ "https://globalwindatlas.info/en/]area/Spain/ para emplazamiento en Soria