Determination of the flexural and shear moduli in curved section specimens by means of bending tests



<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''DETERMINACIÓN DE LOS MÓDULOS DE FLEXIÓN Y CORTADURA EN PROBETAS DE SECCIÓN CURVA MEDIANTE ENSAYOS DE FLEXIÓN'''</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
I. Adarraga <sup>(1)</sup> *, F. Mujika <sup>(1)</sup>, F. de Caso<sup> (2)</sup></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
(1) GRUPO MECÁNICA DE MATERIALES. Departamento de Ingeniería Mecánica. Escuela de Ingeniería de Gipuzkoa (UPV/EHU)</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
(2) University of Miami, Florida, USA.</div>

==1. INTRODUCCIÓN==

En los últimos años se están sustituyendo los postes de madera y hormigón por otros fabricados de poliéster reforzado con fibra de vidrio (PRFV). Las ventajas de este material frente a los tradicionales son su bajo peso, facilidad para el transporte y resistencia a la corrosión, al fuego y a la intemperie.

Debido a la gran variabilidad en las propiedades mecánicas aportadas por los diferentes fabricantes, puede ser de utilidad definir un procedimiento de ensayo de flexión para caracterizar el material utilizando probetas obtenidas a partir del poste. Dado que estas probetas tienen sección curva, es preciso analizar la influencia de la cortadura según varía el ancho y longitud de la probeta.

En este estudio se analiza el efecto de la cortadura en el cálculo los módulos de flexión y cortadura mediante ensayos de flexión.

<span id='_Ref100070158'></span>

==2. CÁLCULO DE LOS MÓDULOS DE FLEXIÓN Y CORTADURA MEDIANTE ENSAYOS DE FLEXIÓN ==

==='''2.1.''' INTRODUCCIÓN===

Mujika [<span id='cite-1'></span>[[#1|1]]] propuso un procedimiento para determinar los módulos de flexión y cortadura del material a partir de la curva carga-desplazamiento en el ensayo de flexión de tres puntos. Considerando la condición de simetría del ensayo se ha modelado como la media viga empotrada de la figura 1.

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

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 1. Modelo del ensayo de tres puntos y sección transversal</div>

De acuerdo con el teorema de Enguesser-Castigliano [<span id='cite-2'></span>[[#2|2]]], el desplazamiento en el punto de aplicación de la carga ''P'' es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math>\delta =\frac{P{L}^{3}}{3{E}_{f}{I}_{z}}+\chi \frac{PL}{GA}</math>
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|(1)
|}


Donde ''M'' es el momento flector, ''V'' es la fuerza cortante, ''E<sub>f</sub>'' es el módulo de flexión, ''G'' es el módulo de cortadura, ''I<sub>z</sub>'' es el momento de inercia de la sección respecto al eje ''z'' que para por su centro de gravedad, ''A'' es el área de la sección transversal y ''&#x03c7;'' es el factor de cortadura. Reordenando los términos resulta:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math>{E}_{f}={E}_{0}\left[ 1+\chi \frac{3{E}_{f}{I}_{z}}{GA}{\left( \frac{1}{L}\right) }^{2}\right]</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref100144974'></span>(2)
|}


Donde <math display="inline">{E}_{0}=\frac{m{L}^{3}}{3{I}_{z}}</math> y <math display="inline">m=</math><math>\frac{P}{\delta }</math> es la pendiente de la curva carga-desplazamiento.

Reordenando la ecuación <span id='cite-_Ref100144974'></span>[[#_Ref100144974|(2)]] resulta:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math>{E}_{0}^{-1}={E}_{f}^{-1}\left[ 1+\chi \frac{3{E}_{f}{I}_{z}}{GA}{\left( \frac{1}{L}\right) }^{2}\right]</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|(3)
|}


Esta ecuación es equivalente a una ecuación lineal <math display="inline">y={A}_{r}+</math><math>{B}_{r}x.</math> Donde:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math display="inline">y={E}_{0}^{-1}</math>, <math display="inline">x={\left( \frac{1}{L}\right) }^{2}</math>, <math display="inline">{A}_{r}=</math><math>{E}_{f}^{-1}</math>, <math display="inline">{B}_{r}=\chi \frac{3{I}_{z}}{A}{G}^{-1}</math>


|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|

<span id='_Ref96946959'></span>(4)


|}


Los coeficientes ''A<sub>r</sub>'' y ''B<sub>r</sub>'' se obtienen por regresión lineal a partir de las ''n'' curvas carga-desplazamiento obtenidas variando'' ''la separación entre apoyos. Conocidos los coeficientes ''A<sub>r</sub>'' y ''B<sub>r</sub>'' se tienen ''E<sub>f</sub>''y ''G'':

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math display="inline">{E}_{f}={{A}_{r}}^{-1}</math>, <math display="inline">G=</math><math>\chi \frac{3{I}_{z}}{A}{{B}_{r}}^{-1}</math>
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref100146883'></span>(5)
|}


Estas ecuaciones son válidas para cualquier tipo de sección, pero previamente hay que calcular el factor de cortante de la sección transversal.

==='''2.2.''' VALORES ESTÁTICOS DE LA SECCIÓN TRANSVERSAL===

En la Fig. 2 se muestra la sección transversal de las probetas objeto de este análisis, donde ''R''<sub>1</sub> y ''R''<sub>2</sub> son los radios exterior e interior respectivamente y ''&#x03c6;''(-''&#x03b8;,&#x03b8;'') el ángulo.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Adarraga_385429259-image2-c.png|186px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 2. Sección transversal</div>

El momento de inercia respecto al eje ''z'' que pasa por el centro de gravedad de la sección es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math>{I}_{z}={R}_{1}^{4}\left( 1-{f}^{2}\right) \left[ \frac{1}{8}\left[ 2\theta +sin\left( 2\theta \right) \right] \left( 1+\right. \right. </math><math>\left. \left. {f}^{2}\right) -\frac{4}{9}\frac{{sin}^{2}\theta }{\theta }{R}_{1}{\left( 1+\frac{\, {f}^{2}}{1+f}\right) }^{2}\right]</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|(6)
|}


Siendo <math display="inline">f={\frac{{R}_{2}}{{R}_{1}}}</math>. La posición del centro de gravedad de la sección, ''y<sub>G</sub>'', es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math>{y}_{G}=\frac{2}{3}\frac{sin\theta }{\theta }\frac{{R}_{1}^{3}-{R}_{2}^{3}}{{R}_{1}^{2}-{R}_{2}^{2}}=</math><math>\frac{2}{3}\frac{sin\theta }{\theta }{R}_{1}\left( 1+\frac{\, {f}^{2}}{1+f}\right)</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref96683171'></span>(7)
|}


Y el área de la sección transversal:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math>A={R}_{1}^{2}\left( 1-{f}^{2}\right)</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref96683179'></span>(8)
|}


==='''2.3.''' FACTOR DE CORTANTE===

El factor de cortante de una sección se calcula a partir de la energía de deformación de cortadura de la viga y su expresión es [<span id='cite-3'></span>[[#3|3]]]:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math>\chi =\frac{A}{{T}^{2}}\int_{A}^{}{\tau }^{2}dA</math>
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref97128431'></span>(9)
|}


Donde ''T'' es el esfuerzo cortante, ''A'' el área de la sección transversal y ''&#x03c4;'' la tensión de cortadura. Descomponiendo la tensión cortante en sus componentes radial y tangencial:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math>{\tau }^{2}={\tau }_{rx}^{2}+{\tau }_{\varphi x}^{2}</math>
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref97128401'></span>(10)
|}


Sustituyendo la ecuación <span id='cite-_Ref97128401'></span>[[#_Ref97128401|(10)]] en <span id='cite-_Ref97128431'></span>[[#_Ref97128431|(9)]],

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math>\chi =\frac{A}{{T}^{2}}\left( \int_{A}^{}{\tau }_{rx}^{2}dA+\int_{A}^{}{\tau }_{\varphi x}^{2}dA\right)</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref97014451'></span>(11)
|}


El factor de cortante se puede descomponer como la suma de los factores de cortadura tangencial y radial.

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math>\chi ={\chi }_{r}+{\chi }_{\varphi }</math>
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref97014446'></span>(12)
|}


La componente tangencial es:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math>{\chi }_{\varphi }=\frac{A}{{T}^{2}}\int_{A}^{}{\tau }_{\varphi x}^{2}dA=\frac{A}{{T}^{2}}\int_{R2}^{R1}rdr\int_{-\theta }^{\theta }{\tau }_{\varphi x}^{2}d\varphi</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref97128543'></span>(13)
|}


Y la componente radial:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math>{\chi }_{r}=\frac{A}{{T}^{2}}\int_{A}^{}{\tau }_{rx}^{2}dA=\frac{A}{{T}^{2}}\int_{R2}^{R1}rdr\int_{-\theta }^{\theta }{\tau }_{rx}^{2}d\varphi</math> 
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref96685162'></span>(14)
|}


==='''2.4.''' ANÁLISIS DE LOS RESULTADOS===

A continuación, se ha analizado la influencia del radio exterior, el espesor y el ángulo en el factor de cortante. Para ángulos pequeños la componente <math display="inline">{\chi }_{\varphi }</math> es despreciable y el factor de cortante tiende a 1,2 como en una sección rectangular. Para 180º la componente <math display="inline">{\chi }_{r}</math> es despreciable y el factor toma el valor 2 como en un tubo hueco delgado.

En la Figura 5 se muestra la influencia de las componentes radial y tangencial del factor de cortante en una probeta de 39mm de radio exterior y 3mm de espesor. El ángulo a partir del cual el efecto del espesor de la probeta puede despreciarse oscila entre 25º y 40º en función de la geometría de la probeta, pero hasta ese valor su influencia debe tenerse en cuenta.

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

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 5 Componentes radial y tangencial del factor de cortante </div>

En la figura 6 se muestra la variación del factor de cortante en una sección curva de 6 mm de espesor, para radios exteriores que varían desde 30 mm hasta 165 mm y ángulos de corte desde 5º hasta 180º.

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

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 6 Influencia del radio y del ángulo de corte </div>

En el gráfico se observa que entre 5º y 45º el factor de cortante aumenta de forma pronunciada y ese efecto aumenta con el radio. Para un radio de 165mm y 20º el factor de cortante es 20. A partir de 90º el valor se aproxima a 2 en todos los casos.

En la figura 7 se analiza la variación del factor de cortante en una sección de 165 mm de radio exterior a la que se va variando el espesor desde 6mm hasta 15mm y ángulos de corte desde 5º hasta 180º. En este caso el factor de cortante disminuye cuando aumenta el espesor, para un mismo ángulo.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Adarraga_385429259-image5.png|474px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig. 7 Influencia del espesor y del ángulo de corte para un radio de 165 mm</div>

==='''2.5.''' COMPARACIÓN CON RESULTADOS NUMÉRICOS===

Con el objetivo de validar numéricamente la expresión analítica del factor de cortante, se ha realizado un análisis por Elementos Finitos de un ensayo de flexión de 3 puntos de una viga de sección curva de poliéster reforzado con fibra de vidrio. En la simulación del ensayo las propiedades del material son datos de entrada y, teniendo en cuenta la condición de simetría del ensayo de tres puntos, se ha modelado como media viga empotrada. No se han incluido los rodillos de apoyo ni de carga, por tanto, no hay que considerar la reducción de la luz entre apoyos que sucedería en un ensayo real. Para evitar los efectos de las deformaciones locales, se ha impuesto un desplazamiento vertical en todos los puntos de la sección del extremo libre.

En el análisis por Elementos Finitos como datos de entrada se han supuesto las propiedades mecánicas que figuran en la Tabla 1 [<span id='cite-4'></span>[[#4|4]]]:

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 1: propiedades del material introducidas en el análisis EF.</div>

{| style="width: 85%;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''E''<sub>1</sub> (GPa)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''E''<sub>2</sub>, ''E''<sub>3</sub> (GPa)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''G''<sub>12</sub>, ''G''<sub>13 </sub>(GPa)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''G''<sub>23</sub> (GPa)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|&#x00a0;&#x03bd;<sub>12</sub>,&#x03bd;<sub>13</sub>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|&#x03bd;<sub>23</sub>
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|39,9
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|8,3
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|3,6
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|3,7
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0,27
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0,11
|}


Para el análisis EF se han utilizado elementos hexaédricos de 8 nodos incompatibles (C3D8I) en el programa ABAQUS ESTÁNDAR [<span id='cite-5'></span>[[#5|5]]]. A partir de las expresiones de <span id='cite-_Ref96946959'></span>[[#_Ref96946959|(4)]] por regresión lineal se obtienen los parámetros'' A<sub>r</sub> ''y'' B<sub>r</sub>''. Como el valor de ''G'' es conocido, de la ecuación <span id='cite-_Ref100146883'></span>[[#_Ref100146883|(5)]] se obtiene el factor de cortante:

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="vertical-align: top;margin:auto;width: 100%;"
|-
| <math>\chi =\frac{{B}_{r}A}{3{I}_{z}G}</math>
|}
|  style="width: 5px;text-align: right;white-space: nowrap;"|(15)
|}


Para realizar la validación numérica, se ha analizado 5 vigas de sección curva de 65mm de radio exterior y 6mm de espesor, con ángulos ''&#x03b8;'' de 5º, 10º, 20º, 30º y 45º.

En la tabla 2 se muestra una comparativa entre los resultados numéricos y analíticos. En todos los casos el error es menor que el 1% excepto para un ángulo de 45º en el que el error es del 2%.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 2: Factores de cortante. ''R''<sub>2</sub>=65mm; ''H''=6mm</div>

{| style="width: 53%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''&#x03b8;'' (º)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|&#x03c7;<sub>analítico</sub>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|&#x03c7;<sub>numérico</sub>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|Error (%)
|-
|  style="text-align: center;"|5
|  style="text-align: center;"|1,20
|  style="text-align: center;"|1,19
|  style="text-align: center;vertical-align: top;"|0,9 %
|-
|  style="text-align: center;"|10
|  style="text-align: center;"|1,24
|  style="text-align: center;"|1,25
|  style="text-align: center;vertical-align: top;"|-0,4 %
|-
|  style="text-align: center;"|20
|  style="text-align: center;"|3,88
|  style="text-align: center;"|3,84
|  style="text-align: center;vertical-align: top;"|1,0 %
|-
|  style="text-align: center;"|30
|  style="text-align: center;"|7,94
|  style="text-align: center;"|7,96
|  style="text-align: center;vertical-align: top;"|-0,2 %
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|45
|  style="border-bottom: 1pt solid black;text-align: center;"|6,52
|  style="border-bottom: 1pt solid black;text-align: center;"|6,66
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-2,1 %
|}


A continuación se ha analizado una viga de 165mm de radio exterior y ''&#x03b8;'' de 10º, para espesores de 2, 4 y 6 mm. En la Tabla 3 se muestran los factores de cortante analíticos y numéricos obtenidos con los diferentes espesores.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 3: Factores de cortante. ''R''<sub>2</sub>=165mm; ''&#x03b8;'' =10º</div>

{| style="width: 55%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''H''(mm)
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|&#x03c7;<sub>analítico</sub>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|&#x03c7;<sub>numérico</sub>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|Error (%)
|-
|  style="text-align: center;"|2
|  style="text-align: center;"|55,55
|  style="text-align: center;"|55,28
|  style="text-align: center;"|0,5 %
|-
|  style="text-align: center;"|4
|  style="text-align: center;"|12,69
|  style="text-align: center;"|12,7
|  style="text-align: center;"|-0,1 %
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|6
|  style="border-bottom: 1pt solid black;text-align: center;"|4,20
|  style="border-bottom: 1pt solid black;text-align: center;"|4,16
|  style="border-bottom: 1pt solid black;text-align: center;"|-1,0 %
|}


En los tres casos analizados el error es menor que el 1%, por lo tanto, se da por validada la aproximación analítica del factor de cortante.

==3. DETERMINACIÓN DE LOS MÓDULOS DE FLEXIÓN Y CORTADURA EN PROBETAS DE SECCIÓN CURVA Y ANCHO CONSTANTE==

A la hora de extraer las probetas, éstas normalmente se cortarían de diferentes partes del poste y con una anchura constante. Teniendo en cuenta que estos postes tienen diferentes radios en cabeza y en base se han analizado probetas obtenidas en la cabeza y en la base de dos postes de diferentes alturas. En esta sección se van a calcular los módulos de flexión y cortadura mediante el procedimiento descrito en el apartado (<span id='cite-_Ref100070158'></span>[[#_Ref100070158|2]]), en probetas de 15mm de anchura [<span id='cite-6'></span>[[#6|6]]] y tomando como factor de cortante el obtenido mediante la aproximación analítica.

En la tabla 4 se muestran los resultados de la simulación del ensayo de 3 puntos. Los módulos de flexión y cortadura se han comparado con los valores introducidos en el programa como datos de entrada.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Tabla 4: Módulos de flexión y cortadura</div>

{| style="width: 80%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
|  colspan='2'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|POSTE 1
|  colspan='2'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|POSTE 2
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Altura (m)
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|14
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Radio exterior en base (mm)
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|165
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|65
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Radio exterior en cabeza (mm)
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|39
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|39
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|Espesor (mm)
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|10
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|&#x03c7;<sub>Base</sub>
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1,200
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1,297
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|&#x03c7;<sub>Cabeza</sub>
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1,207
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1,434
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|BASE
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Numérico
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Error
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Numérico
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Error
|-
|  style="border-top: 1pt solid black;vertical-align: top;"|''E<sub>f</sub>'' (MPa)
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|39983
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|0,2%
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|40239
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|0,8%
|-
|  style="vertical-align: top;"|''G'' (MPa)
|  style="text-align: center;vertical-align: top;"|3576
|  style="text-align: center;vertical-align: top;"|-0,7%
|  style="text-align: center;vertical-align: top;"|3569
|  style="text-align: center;vertical-align: top;"|-0,9%
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|CABEZA
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Numérico
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Error
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Numérico
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Error
|-
|  style="border-top: 1pt solid black;vertical-align: top;"|''E<sub>f</sub>'' (MPa)
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|39953
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|0,1%
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|40268
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|0,9%
|-
|  style="vertical-align: top;"|''G'' (MPa)
|  style="text-align: center;vertical-align: top;"|3571
|  style="text-align: center;vertical-align: top;"|-0,4%
|  style="text-align: center;vertical-align: top;"|3575
|  style="text-align: center;vertical-align: top;"|-0,7%
|}


En todos los casos el error cometido tanto en el módulo de flexión como en el de cortadura es menos del 1%. Si en la probeta obtenida de la parte superior del poste 2, se hubiera supuesto un factor de cortante de 1.2 correspondiente a una sección rectangular, se obtendría un módulo de cortadura de 2981 MPa, con un error de -17,2%.

==4. CONCLUSIONES==

Se ha realizado una aproximación analítica para determinar el factor de cortante en probetas de sección curva obtenidas de postes de poliéster reforzado con fibra de vidrio.

Para validar la aproximación analítica, se ha realizado una simulación por elementos finitos del ensayo de 3 puntos con probetas con diferentes geometrías para obtener el valor del factor de cortante numérico. La aproximación analítica ha sido validada numéricamente.

Finalmente, se han calculado los módulos de flexión y cortadura del material en probetas con diferentes geometrías mediante el ensayo de flexión de 3 puntos, tomando como factor de cortante el obtenido mediante la aproximación analítica. En todos los casos el error ha sido inferior al 1%.

Dependiendo de la geometría de la probeta la influencia del cortante puede ser muy grande y no tener en cuenta este efecto puede originar errores considerables en el cálculo del módulo de cortadura.

==5. AGRADECIMIENTOS==

Los autores agradecen a la Universidad del País Vasco (UPV/EHU) la financiación del Grupo de Investigación Mecánica de Materiales GIU20/060 en la convocatoria de 2020.

==6. REFERENCIAS==

<div id="1"></div>
[[#cite-1|1]]. Mujika F. On the effect of shear and local deformation in three-point bending tests. Polym Test. 2007; 26(7):869-77.

<div id="2"></div>
[[#cite-2|2]]. Oden JT, Ripperger EA. Mechanics of Elastic Structures. 2nd ed. Washington: Hemisphere Publishing Corporation; 1981.

<div id="3"></div>
[[#cite-3|3]]. Gere JM, Timoshenko SP. Mechanics of materials. 4th ed. Pws; 1997

<div id="4"></div>
[[#cite-4|4]]. Morgado T, Silvestre N, Correia JR. Simulation of fire resistance behavior of pultruded GFRP beams – Part I: Models description and kinematic issues. Compos Struct. 2018;187:269-80.

<div id="5"></div>
[[#cite-5|5]]. ABAQUS version 6.12: ABAQUS user's manual, SIMULIA World Headquarters, 166 Valley Street, Providence, RI 02909, USA; 2012

<div id="6"></div>
[[#cite-6|6]] ISO 14125 Fibre-reinforced plastic composites - Determination of flexural properties