(Tag: Visual edit) |
|||
(3 intermediate revisions by the same user not shown) | |||
Line 1: | Line 1: | ||
− | + | ==1 Introducción == | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
La tenacidad a fractura de un material se determina mediante la tasa crítica de liberación de energía, Gc. Para que ocurra el avance de grieta, la tasa de liberación de energía debe superar ese valor crítico. Así, uno de los objetivos principales de los ensayos de fractura es calcular la tasa crítica de liberación de energía. Se supone que se requiere una propagación estable de grieta para poder representar adecuadamente las curvas R [1]. | La tenacidad a fractura de un material se determina mediante la tasa crítica de liberación de energía, Gc. Para que ocurra el avance de grieta, la tasa de liberación de energía debe superar ese valor crítico. Así, uno de los objetivos principales de los ensayos de fractura es calcular la tasa crítica de liberación de energía. Se supone que se requiere una propagación estable de grieta para poder representar adecuadamente las curvas R [1]. | ||
Line 76: | Line 7: | ||
En este trabajo se estudia el comportamiento del crecimiento de grieta en modo mixto I/II mediante la configuración de ensayo ENFR. Se define la condición de estabilidad de grieta en base a la derivada de una nueva tasa de liberación de energía propuesta en este trabajo. Analizando los datos experimentales durante toda la propagación, se comprueba la validez de la tasa de liberación de energía propuesta, y se estudia la estabilidad de grieta. Por último, se resumen algunas recomendaciones para realizar ensayos ENFR en condiciones óptimas | En este trabajo se estudia el comportamiento del crecimiento de grieta en modo mixto I/II mediante la configuración de ensayo ENFR. Se define la condición de estabilidad de grieta en base a la derivada de una nueva tasa de liberación de energía propuesta en este trabajo. Analizando los datos experimentales durante toda la propagación, se comprueba la validez de la tasa de liberación de energía propuesta, y se estudia la estabilidad de grieta. Por último, se resumen algunas recomendaciones para realizar ensayos ENFR en condiciones óptimas | ||
− | == | + | ==2 Descripción del ensayo ENFR== |
La configuración de modo mixto ENFR está basada en el ensayo ENF. Se trata de una probeta con grieta final centrada en el espesor. Con el fin de conseguir un modo mixto, se introduce un rodillo metálico entre ambos brazos de la grieta. El modo II se consigue mediante la aplicación de carga externa, y el modo I o de apertura, por la inserción del rodillo. En este trabajo el modelo análitico incluye las tres posiciones del rodillo mostradas en la <span id='cite-_Ref481564898'></span>[[#_Ref481564898|Figura 1]]. | La configuración de modo mixto ENFR está basada en el ensayo ENF. Se trata de una probeta con grieta final centrada en el espesor. Con el fin de conseguir un modo mixto, se introduce un rodillo metálico entre ambos brazos de la grieta. El modo II se consigue mediante la aplicación de carga externa, y el modo I o de apertura, por la inserción del rodillo. En este trabajo el modelo análitico incluye las tres posiciones del rodillo mostradas en la <span id='cite-_Ref481564898'></span>[[#_Ref481564898|Figura 1]]. | ||
− | + | [[File:Boyano_et_al_2017a-image5-c.png|centre|600x600px]] | |
− | [[ | + | |
<div id="_Ref481564898" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481564898" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | <span style="text-align: center; font-size: 75%;">'''Figura 1 '''Descripción de la configuración de ensayo</span></div> | + | <span style="text-align: center; font-size: 75%;">'''Figura 1 '''Descripción de la configuración de ensayo.</span></div> |
Dónde ''a'' es la longitud de grieta, 2L la distancia entre los apoyos, y c es la distancia de la colocación del rodillo desde el soporte. Cuando se introduce el rodillo, el punto de aplicación de carga sufre un desplazamiento inicial hacia arriba, <math display="inline">{\delta }_0</math> . Al comenzar el ensayo el desplazamiento pasar por cero cuando la carga tiene un valor positivo denominado <math display="inline">P_0</math> , tal y como se muestra en el gráfico de la <span id='cite-_Ref481565380'></span>[[#_Ref481565380|Figura 2]]. | Dónde ''a'' es la longitud de grieta, 2L la distancia entre los apoyos, y c es la distancia de la colocación del rodillo desde el soporte. Cuando se introduce el rodillo, el punto de aplicación de carga sufre un desplazamiento inicial hacia arriba, <math display="inline">{\delta }_0</math> . Al comenzar el ensayo el desplazamiento pasar por cero cuando la carga tiene un valor positivo denominado <math display="inline">P_0</math> , tal y como se muestra en el gráfico de la <span id='cite-_Ref481565380'></span>[[#_Ref481565380|Figura 2]]. | ||
Line 91: | Line 21: | ||
<div id="_Ref481565380" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481565380" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | <span style="text-align: center; font-size: 75%;">'''Figura 2 '''Gráfica teórica de carga-desplazamiento</span></div> | + | <span style="text-align: center; font-size: 75%;">'''Figura 2 '''Gráfica teórica de carga-desplazamiento.</span></div> |
Por tanto, la flexibilidad del ensayo se puede expresar como: | Por tanto, la flexibilidad del ensayo se puede expresar como: | ||
− | <math display="inline">C=\frac{\delta }{P-P_0}=\frac{-{\delta }_0}{P_0}=\frac{\delta -{\delta }_0}{P}</math> () | + | <math display="inline">C=\frac{\delta }{P-P_0}=\frac{-{\delta }_0}{P_0}=\frac{\delta -{\delta }_0}{P}</math>(1) |
Obteniendo los valores de [[Image:Boyano_et_al_2017a-image10-c.png|30px]] y [[Image:Boyano_et_al_2017a-image11-c.png|30px]] para las tres posiciones del ensayo [6], la flexibilidad es: | Obteniendo los valores de [[Image:Boyano_et_al_2017a-image10-c.png|30px]] y [[Image:Boyano_et_al_2017a-image11-c.png|30px]] para las tres posiciones del ensayo [6], la flexibilidad es: | ||
− | [[Image:Boyano_et_al_2017a-image12.png|200px]] () | + | [[Image:Boyano_et_al_2017a-image12.png|200px]] (2) |
Dónde R es el radio del rodillo, c es la distancia desde el apoyo, siendo c0 cuando está hacia el exterior, ci cuando está hacia el interior y c cuando está en la vertical del apoyo. P es la carga aplicada, Ef es el módulo de flexión, b es la anchura de la probeta y 2h es la altura total de la misma. | Dónde R es el radio del rodillo, c es la distancia desde el apoyo, siendo c0 cuando está hacia el exterior, ci cuando está hacia el interior y c cuando está en la vertical del apoyo. P es la carga aplicada, Ef es el módulo de flexión, b es la anchura de la probeta y 2h es la altura total de la misma. | ||
Line 105: | Line 35: | ||
La tasa de liberación de energía para el ensayo ENFR, se puede expresar como [6]: | La tasa de liberación de energía para el ensayo ENFR, se puede expresar como [6]: | ||
− | [[Image:Boyano_et_al_2017a-image13.png|350px]] () | + | [[Image:Boyano_et_al_2017a-image13.png|350px]] (3) |
Los términos en negro corresponden a GI que es la tasa de liberación de energía debida al modo I de fractura. Los términos en magenta corresponden a GII, es decir al modo II de fractura. Esta aproximación análitica ha sido validada mediante elementos finitos utilizando el método de TSEP (Two Step Extension Procedure), [6]. | Los términos en negro corresponden a GI que es la tasa de liberación de energía debida al modo I de fractura. Los términos en magenta corresponden a GII, es decir al modo II de fractura. Esta aproximación análitica ha sido validada mediante elementos finitos utilizando el método de TSEP (Two Step Extension Procedure), [6]. | ||
− | == | + | ==3 Tasa de liberación de energía equivalente== |
La propagación de la grieta sólo se produce si la energía disponible para una extensión de la grieta, G, es suficiente para proporcionar toda la energía que se requiere para crecimiento de la grieta [1]. La energía necesaria para el crecimiento de grieta se denomina tasa de liberación de energía crítica Gc. En el caso de modo mixto I/II, la tasa de liberación de energía es G= GI+ GII, pero GC ≠ GIc+ GIIc. El criterio lineal es adecuado para la representación de la propagación de grieta en la configuración ENFR, donde la relación de modo varía durante el ensayo [7]. Por tanto, la condición para la propagación de las grietas se puede definir como: | La propagación de la grieta sólo se produce si la energía disponible para una extensión de la grieta, G, es suficiente para proporcionar toda la energía que se requiere para crecimiento de la grieta [1]. La energía necesaria para el crecimiento de grieta se denomina tasa de liberación de energía crítica Gc. En el caso de modo mixto I/II, la tasa de liberación de energía es G= GI+ GII, pero GC ≠ GIc+ GIIc. El criterio lineal es adecuado para la representación de la propagación de grieta en la configuración ENFR, donde la relación de modo varía durante el ensayo [7]. Por tanto, la condición para la propagación de las grietas se puede definir como: | ||
Line 116: | Line 46: | ||
{| | {| | ||
|- | |- | ||
− | | <math display="inline"></math> | + | |(4)<math display="inline"></math> |
| <math display="inline">G_IG_{IIc}+G_{II}G_{Ic}\geq G_{Ic}G_{IIc}</math> | | <math display="inline">G_IG_{IIc}+G_{II}G_{Ic}\geq G_{Ic}G_{IIc}</math> | ||
|} | |} | ||
− | |||
El término de la izquierda es la energía disponible para una extensión de grieta y el término de la derecha es el valor crítico de esa energía. Por consiguiente, se puede definir una tasa de liberación de energía equivalente Geq, y una tasa crítica de liberación de energía equivalente Geqc de la siguiente forma: | El término de la izquierda es la energía disponible para una extensión de grieta y el término de la derecha es el valor crítico de esa energía. Por consiguiente, se puede definir una tasa de liberación de energía equivalente Geq, y una tasa crítica de liberación de energía equivalente Geqc de la siguiente forma: | ||
Line 126: | Line 55: | ||
G_{eq}={\left(G_IG_{IIc}+G_{II}G_{Ic}\right)}^{1/2}\\ | G_{eq}={\left(G_IG_{IIc}+G_{II}G_{Ic}\right)}^{1/2}\\ | ||
G_{eqc}={\left(G_{Ic}G_{IIc}\right)}^{1/2} | G_{eqc}={\left(G_{Ic}G_{IIc}\right)}^{1/2} | ||
− | \end{array}</math> () | + | \end{array}</math> (5) |
En consecuencia, se puede decir que la grieta crecerá cuando se cumpla la siguiente condición: | En consecuencia, se puede decir que la grieta crecerá cuando se cumpla la siguiente condición: | ||
− | <math>G_{eq}\geq G_{eqc}</math> () | + | <math>G_{eq}\geq G_{eqc}</math> (6) |
− | == | + | ==4 Estabilidad de la propagación de grieta== |
==='''4.1''' Definición de estabilidad=== | ==='''4.1''' Definición de estabilidad=== | ||
Line 138: | Line 67: | ||
Una definición general de estabilidad es que un sistema es estable si un cambio finito en los parámetros de entrada no causa un cambio infinito en los valores de salida [8]. Con respecto a ensayos de fractura realizados en en máquinas de ensayo, los parámetros de entrada son la carga o el desplazamiento aplicado y el valor de salida es la longitud de la grieta. Así, un ensayo de fractura es estable durante el crecimiento de grieta cuando un cambio infinitesimal de la carga o el desplazamiento no causa un cambio infinito en la longitud de la grieta. En el presente trabajo, la condición de estabilidad se aplica a la tasa de liberación de energía equivalente Geq de la siguiente forma: | Una definición general de estabilidad es que un sistema es estable si un cambio finito en los parámetros de entrada no causa un cambio infinito en los valores de salida [8]. Con respecto a ensayos de fractura realizados en en máquinas de ensayo, los parámetros de entrada son la carga o el desplazamiento aplicado y el valor de salida es la longitud de la grieta. Así, un ensayo de fractura es estable durante el crecimiento de grieta cuando un cambio infinitesimal de la carga o el desplazamiento no causa un cambio infinito en la longitud de la grieta. En el presente trabajo, la condición de estabilidad se aplica a la tasa de liberación de energía equivalente Geq de la siguiente forma: | ||
− | <math>\frac{dG_{eq}}{da}\leq 0</math> () | + | <math>\frac{dG_{eq}}{da}\leq 0</math> (7) |
Si el resultado de la derivada es positivo, el crecimiento de grieta es inestable debido a que la energía liberada es mayor que la necesaria para crear una nueva superficie. Si es negativo, se ha de aplicar trabajo externo para mantener el crecimiento de grieta. | Si el resultado de la derivada es positivo, el crecimiento de grieta es inestable debido a que la energía liberada es mayor que la necesaria para crear una nueva superficie. Si es negativo, se ha de aplicar trabajo externo para mantener el crecimiento de grieta. | ||
Line 145: | Line 74: | ||
Los ensayos llevados a cabo en máquinas de ensayo son considerados de desplazamiento controlado [1]. La condición de estabilidad en desplazamiento controlado se puede expresar como: | Los ensayos llevados a cabo en máquinas de ensayo son considerados de desplazamiento controlado [1]. La condición de estabilidad en desplazamiento controlado se puede expresar como: | ||
+ | [[File:Captura de pantalla 2022-09-20 160317.png|left|479x479px]] | ||
+ | (8) | ||
− | < | + | <nowiki> </nowiki>El material empleado en el estudio analítico de la estabilidad, es un composite unidireccional T300/F593 de matriz epoxi reforzados con fibra de carbono. En la <span id='cite-_Ref481566397'></span>[[#_Ref481566397|Tabla 1]] se presentan algunas características de dicho material, entre las que se encuentran los valores de GIc y GIIc obtenidos del ajuste lineal realizado previamente para comprobar la validez del criterio de fallo lineal [7]. El valor de Geqc que se presenta en la <span id='cite-_Ref481566397'></span>[[#_Ref481566397|Tabla 1]], está calculado sustituyendo esos valores en la ecuación 5. |
− | + | ||
− | El material empleado en el estudio analítico de la estabilidad, es un composite unidireccional T300/F593 de matriz epoxi reforzados con fibra de carbono. En la <span id='cite-_Ref481566397'></span>[[#_Ref481566397|Tabla 1]] se presentan algunas características de dicho material, entre las que se encuentran los valores de GIc y GIIc obtenidos del ajuste lineal realizado previamente para comprobar la validez del criterio de fallo lineal [7]. El valor de Geqc que se presenta en la <span id='cite-_Ref481566397'></span>[[#_Ref481566397|Tabla 1]], está calculado sustituyendo esos valores en la ecuación 5. | + | |
<div id="_Ref481566397" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481566397" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
Line 183: | Line 112: | ||
En la <span id='cite-_Ref481496627'></span>[[#_Ref481496627|Figura 3]], se muestran los resultados de aplicar la condición de la ecuación de estabilidad a las diferentes condiciones de ensayo. | En la <span id='cite-_Ref481496627'></span>[[#_Ref481496627|Figura 3]], se muestran los resultados de aplicar la condición de la ecuación de estabilidad a las diferentes condiciones de ensayo. | ||
− | + | [[File:Boyano_et_al_2017a-image20-c.png|centre|600x600px]] | |
− | [[ | + | |
<div id="_Ref481496627" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481496627" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | <span style="text-align: center; font-size: 75%;">'''Figura 3 '''Estabilidad teórica de la propagación de grieta en desplazamiento controlado</span></div> | + | <span style="text-align: center; font-size: 75%;">'''Figura 3 '''Estabilidad teórica de la propagación de grieta en desplazamiento controlado.</span></div> |
Excepto el caso del ensayo ENF, todas las condiciones de ensayo estudiadas son estables. En el caso de ENF el límite de estabilidad se encuentra en a=42 que es el valor correspondiente a a/L=0,7 [5]. Es más, todas las curvas tienden hacia un valor negativo por lo que parece que son propensas a la estabilidad. | Excepto el caso del ensayo ENF, todas las condiciones de ensayo estudiadas son estables. En el caso de ENF el límite de estabilidad se encuentra en a=42 que es el valor correspondiente a a/L=0,7 [5]. Es más, todas las curvas tienden hacia un valor negativo por lo que parece que son propensas a la estabilidad. | ||
− | == | + | ==5 Evaluación experimental== |
==='''5.1''' Material y aparatos de ensayo=== | ==='''5.1''' Material y aparatos de ensayo=== | ||
Line 243: | Line 171: | ||
En la <span id='cite-_Ref481566908'></span>[[#_Ref481566908|Figura 4]] se representan los valores experimentales punto a punto de la Geq propuesta en este trabajo y definida en la ecuación 8. | En la <span id='cite-_Ref481566908'></span>[[#_Ref481566908|Figura 4]] se representan los valores experimentales punto a punto de la Geq propuesta en este trabajo y definida en la ecuación 8. | ||
− | + | [[File:Boyano_et_al_2017a-image21-c.png|centre|600x600px]] | |
− | [[ | + | |
<div id="_Ref481566908" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481566908" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | <span style="text-align: center; font-size: 75%;">'''Figura 4 '''Valores experimentales de la tasa de liberación de energía equivalente</span></div> | + | <span style="text-align: center; font-size: 75%;">'''Figura 4 '''Valores experimentales de la tasa de liberación de energía equivalente.</span></div> |
Todas las curvas de la tasa equivalente de liberación de energía tienden a un valor constante. El valor promedio se encuentra en el intervalo de 505 y 555 (J/m2) durante toda la propagación. El valor constante mencionado concuerda con el valor teórico de GCeq=550 (J/m2) que se muestra en la <span id='cite-_Ref481566397'></span>[[#_Ref481566397|Tabla 1]]. En cuanto a la tasa equivalente de liberación de energía, la desviación máxima se da en el caso de los ensayos 1 y 2, y el error relativo máximo en este caso es del 8,5 %. Por lo tanto, la aproximación propuesta de la tasa equivalente de liberación de energía se considera válida para definir la condición de propagación de grieta en el ensayo ENFR para el material estudiado. | Todas las curvas de la tasa equivalente de liberación de energía tienden a un valor constante. El valor promedio se encuentra en el intervalo de 505 y 555 (J/m2) durante toda la propagación. El valor constante mencionado concuerda con el valor teórico de GCeq=550 (J/m2) que se muestra en la <span id='cite-_Ref481566397'></span>[[#_Ref481566397|Tabla 1]]. En cuanto a la tasa equivalente de liberación de energía, la desviación máxima se da en el caso de los ensayos 1 y 2, y el error relativo máximo en este caso es del 8,5 %. Por lo tanto, la aproximación propuesta de la tasa equivalente de liberación de energía se considera válida para definir la condición de propagación de grieta en el ensayo ENFR para el material estudiado. | ||
Line 256: | Line 183: | ||
En la <span id='cite-_Ref481566947'></span>[[#_Ref481566947|Figura 5]] se presentan las curvas de carga-desplazamiento de los cinco ensayos estudiados. El marcador circular indica el punto de iniciación de crecimiento de grieta. | En la <span id='cite-_Ref481566947'></span>[[#_Ref481566947|Figura 5]] se presentan las curvas de carga-desplazamiento de los cinco ensayos estudiados. El marcador circular indica el punto de iniciación de crecimiento de grieta. | ||
− | + | [[File:Captura de pantalla 2022-09-20 160527.png|centre|600x600px]] | |
− | [[ | + | |
− | + | ||
<div id="_Ref481566947" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481566947" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | <span style="text-align: center; font-size: 75%;">'''Figura 5 '''Curvas experimentales de carga-desplazamiento</span></div> | + | <span style="text-align: center; font-size: 75%;">'''Figura 5 '''Curvas experimentales de carga-desplazamiento.</span></div> |
Cuando la relación inicial de modo es cercana al modo II puro, el caso de los ensayos 1 y 2, al comenzar la propagación de grieta, la carga continúa aumentando. Eso significa que se está propagando de una manera estable. La carga una vez alcanzado su valor máximo desciende súbitamente. Además, aparecen saltos en dos curvas. Por consiguiente, los ensayos 1 y 2 se consideran estables durante la mayor parte de la propagación, pero en la etapa final se vuelven inestables. En el resto de los casos, no hay saltos o caídas de carga repentinas durante la propagación de grietas. Se puede concluir entonces, que todos ellos son estables | Cuando la relación inicial de modo es cercana al modo II puro, el caso de los ensayos 1 y 2, al comenzar la propagación de grieta, la carga continúa aumentando. Eso significa que se está propagando de una manera estable. La carga una vez alcanzado su valor máximo desciende súbitamente. Además, aparecen saltos en dos curvas. Por consiguiente, los ensayos 1 y 2 se consideran estables durante la mayor parte de la propagación, pero en la etapa final se vuelven inestables. En el resto de los casos, no hay saltos o caídas de carga repentinas durante la propagación de grietas. Se puede concluir entonces, que todos ellos son estables | ||
La relación de la longitud de grieta a lo largo del desplazamiento se puede observar en la <span id='cite-_Ref481566974'></span>[[#_Ref481566974|Figura 6]]. | La relación de la longitud de grieta a lo largo del desplazamiento se puede observar en la <span id='cite-_Ref481566974'></span>[[#_Ref481566974|Figura 6]]. | ||
− | + | [[File:Boyano_et_al_2017a-image23-c.png|centre|600x600px]] | |
− | [[ | + | |
<div id="_Ref481566974" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481566974" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | <span style="text-align: center; font-size: 75%;">'''Figura 6 '''Longitud de grieta en función del desplazamiento</span></div> | + | <span style="text-align: center; font-size: 75%;">'''Figura 6 '''Longitud de grieta en función del desplazamiento.</span></div> |
De acuerdo a la <span id='cite-_Ref481566974'></span>[[#_Ref481566974|Figura 6]] , en las curvas correspondientes a los ensayos 1 y 2, en la primera etapa de la propagación, un cambio finito en el desplazamiento causa un cambio finito en la longitud de la grieta y por lo tanto puede considerarse estable, de acuerdo a la definición de estabilidad de la sección 3.1. Sin embargo, en la fase final de la propagación, un pequeño cambio en el desplazamiento causa un cambio mucho mayor en longitud de la grieta, y hay incluso saltos en la curva. Por lo tanto, se vuelven inestables. Para el resto de las pruebas, la propagación de la grieta comienza antes, y la pendiente de la curva es similar, lo que significa que se necesita un cambio finito en el desplazamiento para tener un cambio finito en la longitud de la grieta. En consecuencia, todos los ensayos son estables según esta definición. Esta es la misma conclusión que se extrae mediante el análisis de curvas de carga-desplazamiento. | De acuerdo a la <span id='cite-_Ref481566974'></span>[[#_Ref481566974|Figura 6]] , en las curvas correspondientes a los ensayos 1 y 2, en la primera etapa de la propagación, un cambio finito en el desplazamiento causa un cambio finito en la longitud de la grieta y por lo tanto puede considerarse estable, de acuerdo a la definición de estabilidad de la sección 3.1. Sin embargo, en la fase final de la propagación, un pequeño cambio en el desplazamiento causa un cambio mucho mayor en longitud de la grieta, y hay incluso saltos en la curva. Por lo tanto, se vuelven inestables. Para el resto de las pruebas, la propagación de la grieta comienza antes, y la pendiente de la curva es similar, lo que significa que se necesita un cambio finito en el desplazamiento para tener un cambio finito en la longitud de la grieta. En consecuencia, todos los ensayos son estables según esta definición. Esta es la misma conclusión que se extrae mediante el análisis de curvas de carga-desplazamiento. | ||
− | En la <span id='cite-_Ref481566989'></span>[[#_Ref481566989|Figura 7]] se puede observar el valor experimental de la ecuación 5 . | + | En la <span id='cite-_Ref481566989'></span>[[#_Ref481566989|Figura 7]] se puede observar el valor experimental de la ecuación 5. |
− | + | [[File:Boyano_et_al_2017a-image24-c.png|centre|600x600px]] | |
− | [[ | + | |
<div id="_Ref481566989" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div id="_Ref481566989" class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | <span style="text-align: center; font-size: 75%;">'''Figura 7 '''Estudio experimental de la condición de estabilidad</span></div> | + | <span style="text-align: center; font-size: 75%;">'''Figura 7 '''Estudio experimental de la condición de estabilidad.</span></div> |
De acuerdo a la <span id='cite-_Ref481566989'></span>[[#_Ref481566989|Figura 7]], en el caso de los ensayos 1 y 2, la derivada de Geq durante la propagación tiene un valor cercano a cero. Aún así, todos los ensayos cumplen la condición de estabilidad ya que la derivada es siempre negativa. En los casos de ensayos con relación inicial de modo baja, el valor negativo de la derivada es mayor lo que significa que el ensayo es más estable. | De acuerdo a la <span id='cite-_Ref481566989'></span>[[#_Ref481566989|Figura 7]], en el caso de los ensayos 1 y 2, la derivada de Geq durante la propagación tiene un valor cercano a cero. Aún así, todos los ensayos cumplen la condición de estabilidad ya que la derivada es siempre negativa. En los casos de ensayos con relación inicial de modo baja, el valor negativo de la derivada es mayor lo que significa que el ensayo es más estable. | ||
− | == | + | ==6 Recomendaciones para realizar ensayos== |
Al realizar ensayos, para una geometría de probeta dada, son tres los parámetros que se han de definir: longitud inicial de grieta, radio del rodillo y posición del rodillo. Diferentes combinaciones de estos tres parámetros, provocarán diferentes ratios de modo. Por otra lado, se ha de tener en cuenta la estabilidad de la propagación de la grieta. | Al realizar ensayos, para una geometría de probeta dada, son tres los parámetros que se han de definir: longitud inicial de grieta, radio del rodillo y posición del rodillo. Diferentes combinaciones de estos tres parámetros, provocarán diferentes ratios de modo. Por otra lado, se ha de tener en cuenta la estabilidad de la propagación de la grieta. | ||
Line 296: | Line 219: | ||
Entre las combinaciones de los valores propuestos se elegirán las que mejor cumplan el criterio de fallo lineal [7] siendo éstas las que tienen una relación de modo inicial entre 65 % y 75 %. Este hecho se puede constatar, observando la <span id='cite-_Ref481566908'></span>[[#_Ref481566908|Figura 4]] que muestra valores experimentales de Geq. | Entre las combinaciones de los valores propuestos se elegirán las que mejor cumplan el criterio de fallo lineal [7] siendo éstas las que tienen una relación de modo inicial entre 65 % y 75 %. Este hecho se puede constatar, observando la <span id='cite-_Ref481566908'></span>[[#_Ref481566908|Figura 4]] que muestra valores experimentales de Geq. | ||
− | == | + | ==7 Conclusiones== |
Se ha propuesto una tasa de liberación de energía equivalente basada en el cumplimiento del criterio de fallo lineal. Esta tasa equivalente tiene en cuenta la interacción entre ambos modos de fractura. Los resultados experimentales obtenidos de probetas unidireccionales de carbono/epoxy F593/T300 muestran un buen acuerdo con la propuesta. Además, durante la propagación de la grieta el valor crítico de Geq, es decir Geqc, alcanza un valor constante y similar para todas las condiciones de ensayo estudiadas, por lo que puede considerarse una propiedad del material. | Se ha propuesto una tasa de liberación de energía equivalente basada en el cumplimiento del criterio de fallo lineal. Esta tasa equivalente tiene en cuenta la interacción entre ambos modos de fractura. Los resultados experimentales obtenidos de probetas unidireccionales de carbono/epoxy F593/T300 muestran un buen acuerdo con la propuesta. Además, durante la propagación de la grieta el valor crítico de Geq, es decir Geqc, alcanza un valor constante y similar para todas las condiciones de ensayo estudiadas, por lo que puede considerarse una propiedad del material. |
La tenacidad a fractura de un material se determina mediante la tasa crítica de liberación de energía, Gc. Para que ocurra el avance de grieta, la tasa de liberación de energía debe superar ese valor crítico. Así, uno de los objetivos principales de los ensayos de fractura es calcular la tasa crítica de liberación de energía. Se supone que se requiere una propagación estable de grieta para poder representar adecuadamente las curvas R [1].
En las últimas décadas, se han desarrollado diferentes configuraciones para analizar el modo mixto I/II de fractura. Sin embargo, en la norma no se define ninguna condición de estabilidad. Se menciona que para casos de relación de modo alta la propagación suele ser inestable y que eso impide conocer los valores de propagación de la tenacidad a fractura. También se comenta en la norma, que las grietas iniciales más largas, tienen mayor tendencia a una propagación estable [2]. En cuanto a los ensayos de modos puros I y II, la configuración DCB se considera estable siempre [3], y la configuración ENF se considera estable cuando se cumple que a/L>0,7 [4,5].
En este trabajo se estudia el comportamiento del crecimiento de grieta en modo mixto I/II mediante la configuración de ensayo ENFR. Se define la condición de estabilidad de grieta en base a la derivada de una nueva tasa de liberación de energía propuesta en este trabajo. Analizando los datos experimentales durante toda la propagación, se comprueba la validez de la tasa de liberación de energía propuesta, y se estudia la estabilidad de grieta. Por último, se resumen algunas recomendaciones para realizar ensayos ENFR en condiciones óptimas
La configuración de modo mixto ENFR está basada en el ensayo ENF. Se trata de una probeta con grieta final centrada en el espesor. Con el fin de conseguir un modo mixto, se introduce un rodillo metálico entre ambos brazos de la grieta. El modo II se consigue mediante la aplicación de carga externa, y el modo I o de apertura, por la inserción del rodillo. En este trabajo el modelo análitico incluye las tres posiciones del rodillo mostradas en la Figura 1.
Dónde a es la longitud de grieta, 2L la distancia entre los apoyos, y c es la distancia de la colocación del rodillo desde el soporte. Cuando se introduce el rodillo, el punto de aplicación de carga sufre un desplazamiento inicial hacia arriba, . Al comenzar el ensayo el desplazamiento pasar por cero cuando la carga tiene un valor positivo denominado , tal y como se muestra en el gráfico de la Figura 2.
Por tanto, la flexibilidad del ensayo se puede expresar como:
(1)
Obteniendo los valores de y
para las tres posiciones del ensayo [6], la flexibilidad es:
Dónde R es el radio del rodillo, c es la distancia desde el apoyo, siendo c0 cuando está hacia el exterior, ci cuando está hacia el interior y c cuando está en la vertical del apoyo. P es la carga aplicada, Ef es el módulo de flexión, b es la anchura de la probeta y 2h es la altura total de la misma.
La tasa de liberación de energía para el ensayo ENFR, se puede expresar como [6]:
Los términos en negro corresponden a GI que es la tasa de liberación de energía debida al modo I de fractura. Los términos en magenta corresponden a GII, es decir al modo II de fractura. Esta aproximación análitica ha sido validada mediante elementos finitos utilizando el método de TSEP (Two Step Extension Procedure), [6].
La propagación de la grieta sólo se produce si la energía disponible para una extensión de la grieta, G, es suficiente para proporcionar toda la energía que se requiere para crecimiento de la grieta [1]. La energía necesaria para el crecimiento de grieta se denomina tasa de liberación de energía crítica Gc. En el caso de modo mixto I/II, la tasa de liberación de energía es G= GI+ GII, pero GC ≠ GIc+ GIIc. El criterio lineal es adecuado para la representación de la propagación de grieta en la configuración ENFR, donde la relación de modo varía durante el ensayo [7]. Por tanto, la condición para la propagación de las grietas se puede definir como:
(4) |
El término de la izquierda es la energía disponible para una extensión de grieta y el término de la derecha es el valor crítico de esa energía. Por consiguiente, se puede definir una tasa de liberación de energía equivalente Geq, y una tasa crítica de liberación de energía equivalente Geqc de la siguiente forma:
(5)
En consecuencia, se puede decir que la grieta crecerá cuando se cumpla la siguiente condición:
(6)
Una definición general de estabilidad es que un sistema es estable si un cambio finito en los parámetros de entrada no causa un cambio infinito en los valores de salida [8]. Con respecto a ensayos de fractura realizados en en máquinas de ensayo, los parámetros de entrada son la carga o el desplazamiento aplicado y el valor de salida es la longitud de la grieta. Así, un ensayo de fractura es estable durante el crecimiento de grieta cuando un cambio infinitesimal de la carga o el desplazamiento no causa un cambio infinito en la longitud de la grieta. En el presente trabajo, la condición de estabilidad se aplica a la tasa de liberación de energía equivalente Geq de la siguiente forma:
(7)
Si el resultado de la derivada es positivo, el crecimiento de grieta es inestable debido a que la energía liberada es mayor que la necesaria para crear una nueva superficie. Si es negativo, se ha de aplicar trabajo externo para mantener el crecimiento de grieta.
Los ensayos llevados a cabo en máquinas de ensayo son considerados de desplazamiento controlado [1]. La condición de estabilidad en desplazamiento controlado se puede expresar como:
(8)
El material empleado en el estudio analítico de la estabilidad, es un composite unidireccional T300/F593 de matriz epoxi reforzados con fibra de carbono. En la Tabla 1 se presentan algunas características de dicho material, entre las que se encuentran los valores de GIc y GIIc obtenidos del ajuste lineal realizado previamente para comprobar la validez del criterio de fallo lineal [7]. El valor de Geqc que se presenta en la Tabla 1, está calculado sustituyendo esos valores en la ecuación 5.
Propiedad | Valor |
E1f1 (MPa) | 107,4 |
E22 (MPa) | 10,6 |
G12 (MPa) | 4,3 |
GIc (J/m2) | 263 |
GIIc (J/m2) | 1152 |
Geqc (J/m2) | 550 |
Las dimensiones de probeta son: longitud 200 mm, espesor (2h) 3 mm y anchura (b) 15 mm. Se considera una luz (2L) de 120 mm.
Para determinar la condición de estabilidad es necesario conocer la carga aplicada para cada longitud de grieta. Se ha elegido el punto crítico o punto de iniciación para realizar el estudio de la estabilidad. La nomenclatura utilizada en este trabajo, es: ai-Rj-ck. Siendo ai la longitud de grieta inicial, Rj el radio del rodillo y ck el parámetro que define la posición del mismo. Por ejemplo, a40-R0,5-c10 hace referencia a un ensayo en el que se ha introducido un rodillo de 0,5 mm de radio a 10 mm de distancia del soporte hacia el interior de la grieta siendo la grieta inicial de 40mm.
En la Figura 3, se muestran los resultados de aplicar la condición de la ecuación de estabilidad a las diferentes condiciones de ensayo.
Excepto el caso del ensayo ENF, todas las condiciones de ensayo estudiadas son estables. En el caso de ENF el límite de estabilidad se encuentra en a=42 que es el valor correspondiente a a/L=0,7 [5]. Es más, todas las curvas tienden hacia un valor negativo por lo que parece que son propensas a la estabilidad.
Las probetas utilizadas para el estudio experimental de la propagación de grieta, son composites unidireccionales T300/F593 de matriz epoxi reforzados con fibra de carbono. Se han fabricado laminados unidireccionales de 16 capas con un 55% de porcentaje de fibra. La grieta inicial se ha creado colocando una lámina de teflón en mitad del proceso de apilamiento. Se han determinado los módulos Ef y GLT’ en cada probeta, siendo los resultados obtenidos:
Ef = 107.4 (±1.4) GPa
GLT’ = 4.3 (±0.4) GPa
En cuanto a las dimensiones de las probetas, la longitud es 200 mm. Los valores nominales de la anchura (b) y espesor (2h) son 15 mm y 3 mm respectivamente. Todos los ensayos se realizan con una luz (2L) de 120 mm en una máquina de ensayos MTS Insight 100KN equipada con una célula de carga de 5KN con control de desplazamiento.
Se han llevado a cabo diferentes ensayos de modo mixto I/II con el objetivo de analizar la propagación de grieta. En la configuración ENFR la longitud de grieta se determina sin utilizar ningún medio óptico, en los puntos en los que se conocen la carga aplicada y el desplazamiento experimental, aplicando el método basado en la variación de la flexibilidad experimental [7].
Debido a que la relación de modo varía durante el ensayo, se utiliza la relación de modo en el punto de iniciación como parámetro para definir el modo mixto del ensayo. En la Tabla 2 se presenta un resumen de las condiciones de ensayo estudiadas.
ID | NOMENCLATURA | RELACIÓN DE MODO GII/G INICIAL (%) |
1. | a40-R0,5-c0 | 96 |
2. | a31-R0,5-c0 | 90 |
3. | a45-R1,5-c0 | 74 |
4. | a43-R0,9-c8 | 66 |
5. | a38-R0,7-c8 | 59 |
En la Figura 4 se representan los valores experimentales punto a punto de la Geq propuesta en este trabajo y definida en la ecuación 8.
Todas las curvas de la tasa equivalente de liberación de energía tienden a un valor constante. El valor promedio se encuentra en el intervalo de 505 y 555 (J/m2) durante toda la propagación. El valor constante mencionado concuerda con el valor teórico de GCeq=550 (J/m2) que se muestra en la Tabla 1. En cuanto a la tasa equivalente de liberación de energía, la desviación máxima se da en el caso de los ensayos 1 y 2, y el error relativo máximo en este caso es del 8,5 %. Por lo tanto, la aproximación propuesta de la tasa equivalente de liberación de energía se considera válida para definir la condición de propagación de grieta en el ensayo ENFR para el material estudiado.
Durante el trabajo experimental se ha considerado que el crecimiento inestable de la grieta se produce cuando hay un salto en la curva de carga-desplazamiento.
En la Figura 5 se presentan las curvas de carga-desplazamiento de los cinco ensayos estudiados. El marcador circular indica el punto de iniciación de crecimiento de grieta.
Cuando la relación inicial de modo es cercana al modo II puro, el caso de los ensayos 1 y 2, al comenzar la propagación de grieta, la carga continúa aumentando. Eso significa que se está propagando de una manera estable. La carga una vez alcanzado su valor máximo desciende súbitamente. Además, aparecen saltos en dos curvas. Por consiguiente, los ensayos 1 y 2 se consideran estables durante la mayor parte de la propagación, pero en la etapa final se vuelven inestables. En el resto de los casos, no hay saltos o caídas de carga repentinas durante la propagación de grietas. Se puede concluir entonces, que todos ellos son estables
La relación de la longitud de grieta a lo largo del desplazamiento se puede observar en la Figura 6.
De acuerdo a la Figura 6 , en las curvas correspondientes a los ensayos 1 y 2, en la primera etapa de la propagación, un cambio finito en el desplazamiento causa un cambio finito en la longitud de la grieta y por lo tanto puede considerarse estable, de acuerdo a la definición de estabilidad de la sección 3.1. Sin embargo, en la fase final de la propagación, un pequeño cambio en el desplazamiento causa un cambio mucho mayor en longitud de la grieta, y hay incluso saltos en la curva. Por lo tanto, se vuelven inestables. Para el resto de las pruebas, la propagación de la grieta comienza antes, y la pendiente de la curva es similar, lo que significa que se necesita un cambio finito en el desplazamiento para tener un cambio finito en la longitud de la grieta. En consecuencia, todos los ensayos son estables según esta definición. Esta es la misma conclusión que se extrae mediante el análisis de curvas de carga-desplazamiento.
En la Figura 7 se puede observar el valor experimental de la ecuación 5.
De acuerdo a la Figura 7, en el caso de los ensayos 1 y 2, la derivada de Geq durante la propagación tiene un valor cercano a cero. Aún así, todos los ensayos cumplen la condición de estabilidad ya que la derivada es siempre negativa. En los casos de ensayos con relación inicial de modo baja, el valor negativo de la derivada es mayor lo que significa que el ensayo es más estable.
Al realizar ensayos, para una geometría de probeta dada, son tres los parámetros que se han de definir: longitud inicial de grieta, radio del rodillo y posición del rodillo. Diferentes combinaciones de estos tres parámetros, provocarán diferentes ratios de modo. Por otra lado, se ha de tener en cuenta la estabilidad de la propagación de la grieta.
En consecuencia, los rangos propuestos de los parámetros para realizar ensayos ENFR son:
Entre las combinaciones de los valores propuestos se elegirán las que mejor cumplan el criterio de fallo lineal [7] siendo éstas las que tienen una relación de modo inicial entre 65 % y 75 %. Este hecho se puede constatar, observando la Figura 4 que muestra valores experimentales de Geq.
Se ha propuesto una tasa de liberación de energía equivalente basada en el cumplimiento del criterio de fallo lineal. Esta tasa equivalente tiene en cuenta la interacción entre ambos modos de fractura. Los resultados experimentales obtenidos de probetas unidireccionales de carbono/epoxy F593/T300 muestran un buen acuerdo con la propuesta. Además, durante la propagación de la grieta el valor crítico de Geq, es decir Geqc, alcanza un valor constante y similar para todas las condiciones de ensayo estudiadas, por lo que puede considerarse una propiedad del material.
Se ha analizado la estabilidad de la propagación de grieta en la configuración ENFR mediante la derivada de Geq. Los resultados experimentales han confirmado que en ensayos a desplazamiento constante todos los casos estudiados han resultado estables.
Los autores agradecen a la Universidad del País Vasco (UPV/EHU) la financiación del Grupo de Investigación Mecánica de Materiales GIU 16/51 en la convocatoria de 2016.
[1] Anderson TL. Fracture mechanics :fundamentals and applications. 3rd ed. : CRC Press, 2005.
[2] ASTM D6671-01. Standard test method for mixed mode I-mode II interlaminar fracture toughness of unidirectional fiber-reinforced polymer matrix composites 2006.
[3] ASTM D5528 – 01 Standard Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites. 2001.
[4] ASTM D7905/D7905M-14 Standard Test Method for Determination of the Mode II Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites. 2014.
[5] Carlsson LA, Gillespie JW, JR, Pipes RB. On the Analysis and Design of the End Notched Flexure (ENF) Specimen for Mode II Testing. Journal of Composite Materials 1986;20:594-604.
[6] Boyano A, Mollón V, Bonhomme J, De Gracia J, Arrese A, Mujika F. Analytical and numerical approach of an End Notched Flexure test configuration with an inserted roller for promoting mixed mode I/II. Eng Fract Mech 2015;143:63-79.
[7] Boyano A, De Gracia J, Arrese A, Mujika F. Experimental assessment of an End Notched Flexure test configuration with an inserted roller for analyzing mixed-mode I/II fracture toughness. Eng Fract Mech 2016;163:462-75.
[8] Bazant Z, Cedolin L. Stability of structures: elastic, inelastic, fracture, and damage theories. Oxford University Press. Inc, New York 1991.
Published on 19/10/17
Accepted on 19/10/17
Submitted on 19/10/17
Volume 01 - Comunicaciones Matcomp17 (2017), Issue Núm. 1 - Comportamiento en Servicio de los Materiales Compuestos, 2017
DOI: 10.23967/r.matcomp.2017.10.018
Licence: Other
Are you one of the authors of this document?