(Created page with "== Abstract == A novel method is presented for the determination of mode II cohesive law for the characterization of delamination in thin adhesive joints The proposed method...")
 
 
(6 intermediate revisions by 2 users not shown)
Line 1: Line 1:
== Abstract ==
+
==1 Introducción ==
  
A novel method is presented for the determination of mode II cohesive law for the characterization of delamination in thin adhesive joints
+
Se presenta un nuevo modelo para la determinación de la ley cohesiva en modo II de uniones adhesivas de pequeño espesor. La tasa de liberación de energía G en función del avance de grieta &#x0394;a se determina mediante el método de reducción de datos denominado Beam Theory including Bending Rotation effects (BTBR) [<span id='cite-1'></span>[[#1|1]]]. Mediante este método, la longitud de grieta para cada punto del ensayo se determina mediante la variación de flexibilidad. El desplazamiento relativo de cortadura en la punta de grieta &#x0394; se determina en función del avance de grieta equivalente, asumiendo que el desarrollo de la Zona de Proceso de Fractura (FPZ) es análogo al avance de grieta. Finalmente, se representa G en función de &#x0394;&#x00a0;y se determina la ley cohesiva &#x03c4;&#x00a0;= &#x03c4;&#x00a0;(&#x0394;) mediante derivación numérica..
The proposed method is based on the three point end notched flexure test. The data reduction method used is the Beam Theory including Bending Rotation effects method (BTBR), which includes the effect of local deformations, shear and bending rotations. By this method, the initial crack length, the apparent crack length during the crack advance, and the mode II resistance curve are determined.
+
The relative shear displacement at the crack tip is determined as a function of the apparent crack advance that occurs when the generation of the fracture zone progresses. In order to determine the shear displacement an analytical analytical model based on the classical theory of beams is used.
+
Finally, the fracture toughness versus the shear relative displacement is represented and the cohesive law is determined by numerical derivation.
+
  
== Full document ==
+
==2 Metodología ==
<pdf>Media:Draft_Content_550729540native-articles-177.pdf</pdf>
+
 
 +
La Fig. 1 muestra la configuración del ensayo de flexión en tres puntos con entalla final (ENF) de acuerdo a BTBR
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
[[Image:Arrese_et_al_2018a-image5-c.png|534px]] </div>
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 1.''' Configuración del ensayo ENF de acuerdo a BTBR.</span></div>
 +
 
 +
La integral-J es equivalente a la tasa de liberación de energía G para un cuerpo elástico. [<span id='cite-2'></span>[[#2|2]]]
 +
 
 +
 
 +
{|
 +
|-
 +
| [[Image:Arrese_et_al_2018a-image6-c.png|330px]]
 +
| [[Image:Arrese_et_al_2018a-image7-c.png|center|90px]]
 +
|}
 +
(1)
 +
 
 +
En la aproximación propuesta se asume que el desarrollo de la FPZ es análogo al avance de grieta equivalente &#x0394;''a'', relacionado con la variación de la flexibilidad.
 +
[[File:Arrese_et_al_2018a-image8-c.png|centre|600x600px]]
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 2. '''Zona de proceso de Fractura.</span></div>
 +
 
 +
De acuerdo a la teoría clásica de vigas &#x0394; puede expresarse de la siguiente forma:
 +
 
 +
[[Image:Arrese_et_al_2018a-image9-c.png|456px]] (2)
 +
 
 +
donde
 +
 
 +
[[Image:Arrese_et_al_2018a-image10-c.png|240px]]
 +
 
 +
La relación G=G(&#x0394;) se aproxima a una función logística de la siguiente forma:
 +
 
 +
[[Image:Arrese_et_al_2018a-image11-c.png|354px]] (3)
 +
 
 +
F(&#x0394;) se determina mediante un ajuste lineal por regiones de los datos experimentales. Finalmente la ley cohesiva se obtiene median derivación numérica.
 +
 
 +
[[Image:Arrese_et_al_2018a-image12-c.png|336px]] (4)
 +
 
 +
==3 Procedimiento Experimental ==
 +
 
 +
==='''3.1''' Material ===
 +
 
 +
Para la fabricación de las probetas se han utilizado probetas unidireccionales fabricadas a partir de prepregs de carbono/epoxi T6T/F593 proporcionados por Hexcel Composites con 16 capas unidireccionales que han sido  delaminadas completamente.  Los adherentes fruto de la delaminación completa han sido unidas por el adhesivo Loctite 9644 Hysol.
 +
 
 +
Las dimensiones nominales de la probeta son 3 mm de espesor y 16mm de ancho. La grieta inicial se genera mediante un inserto de teflón.
 +
 
 +
Los ensayos se han realizado en una maquina de ensayo universal MTS-Insight 10 con célula de carga de 10kN
 +
 
 +
==='''3.2'''  Ensayos Preliminares===
 +
 
 +
Las propiedades elásticas E<sub>f</sub> and G<sub>13</sub>, se han determinado mediante el procedimiento propuesto por Mujika [<span id='cite-3'></span>[[#3|3]]], siendo el modulo de flexión 107GPa y el modulo de cortadura 3 GPa.
 +
 
 +
El desplazamiento de la probeta (&#x03b4;<sub>spec</sub>) se determina de las curvas carga desplazamiento. El desplazamiento experimental (&#x03b4;<sub>exp</sub>) es la suma  del desplazamiento de la probeta y el desplazamiento debido a la flexibilidad del sistema.
 +
 
 +
[[Image:Arrese_et_al_2018a-image13-c.png|414px]]
 +
 
 +
La flexibilidad del sistema se determina experimentalmente siendo su valor promedio C<sub>s</sub> = 1/15•10<sup>-3</sup> mm/N.
 +
 
 +
==='''3.3''' Resultados===
 +
 
 +
En la Figura 3 se muestra la curva experimental carga desplazamiento correspondiente a la probeta S1.
 +
[[File:Arrese_et_al_2018a-image14-c.png|centre|600x600px]]
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 3.''' Curva carga-desplazamiento.</span></div>
 +
 
 +
Una vez determinada la longitud de grieta para cada punto del ensayo, la longitud de grieta inicial se determina como promedio de los valores en el rango de carga donde la longitud de grieta se mantiene constante, tal y como se aprecia en la Fig. 4. Por lo tanto, la iniciación del daño local está relacionada con el primer punto donde el avance de grieta es distinto de cero.
 +
[[File:Arrese_et_al_2018a-image15-c.png|centre|600x600px]]
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 4.''' Evolución de la longitud de grieta. </span></div>
 +
 
 +
Teniendo la longitud de de grieta inicial y la longitud de grieta para cada punto del ensayo, es posible obtener la curva-R, la cual representa la evolución de G respecto al avance de grieta, Fig. 5.
 +
[[File:Arrese_et_al_2018a-image16-c.png|centre|522x522px]]
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 5.''' Curva R.</span></div>
 +
 
 +
De acuerdo a la  Eq. (2) el desplazamiento de cortadura en la punta de grieta puede ser determinado para cada punto del ensayo sustituyendo el correspondiente avance de grieta. La Fig 6. muestra la curva  G¸Ç_ëž&#x0394; que corresponde a la probeta S1.
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;"> [[Image:Arrese_et_al_2018a-image17-c.png|522px]] </span></div>
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 6.''' Curva.</span>G-&#x0394;<span style="text-align: center; font-size: 75%;">.</span></div>
 +
 
 +
Para evitar el ruido excesivo de los datos experimentales, los datos los datos G-&#x0394; se aproximan a la función logística de la Eq (3). La función F (&#x0394;) se ajusta a los datos experimentales mediante un  ajuste lineal por regiones. Fig. 7 muestra los datos experimentales obtenidos de la Eq. (3) y la curva ajustada. Con el fin de evaluar la precisión del ajuste, los datos experimentales y la función logística aproximada se representan en la Fig. 8.
 +
[[File:Arrese_et_al_2018a-image18-c.png|centre|600x600px]]
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 7.''' La función F.(&#x0394;)</span></div>
 +
[[File:Arrese_et_al_2018a-image19-c.png|centre|600x600px]]
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 8.''' (G-G<sub>ini</sub>)/(G<sub>ss</sub>-G<sub>ini</sub>) experimental y aproximada. </span></div>
 +
 
 +
Reemplazando los valores de F (&#x0394;) y sus derivadas obtenidas por diferenciación numérica en la Eq. (4), se obtiene la ley cohesiva.
 +
 
 +
La Fig. 9 muestra la ley cohesiva de la probeta S1
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
[[Image:Arrese_et_al_2018a-image20-c.png|600px]] </div>
 +
 
 +
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 +
<span style="text-align: center; font-size: 75%;">'''Figura 9.''' Ley cohesiva de la probeta S1.</span></div>
 +
 
 +
==4 Conclusiones==
 +
 
 +
En el presente trabajo se propone un nuevo método para la determinación de la ley cohesiva en modo II para uniones adhesivas de pequeño espesor.
 +
 
 +
La aproximación propuesta, se basa en la correlación entre la tasa de liberación de energía desarrollada (G) en función del desplazamiento de cortadura de la punta de grieta (&#x0394;).
 +
 
 +
La Curva de Resistencia se determina obteniendo G en función del avance de grieta cuya determinación se basa en la variación de la flexibilidad de la probeta debido al avance del daño.
 +
 
 +
El desplazamiento de cortadura de la punta de la grieta se determina en función del avance de grieta equivalente suponiendo que el desarrollo de la Zona de Proceso de Fractura es análogo a un avance grieta.
 +
 
 +
Finalmente,  representando la tasa de liberación de energía y el desplazamiento de cortadura de la grieta determinados, y mediante un ajuste a una función logística de la curva G-&#x0394;&#x00a0;y la ley cohesiva experimental se determina por diferenciación numérica.
 +
 
 +
Mediante el método propuesto se ha evaluado una unión realizada con LOCTITE HYSOL 9644. La resistencia cohesiva obtenida es del orden de 34 MPa,
 +
 
 +
El método propuesto proporciona una manera sencilla de obtener la ley cohesiva en modo II utilizando únicamente datos de desplazamiento y carga proporcionados por la máquina de ensayo universal, sin necesidad de ninguna medición de desplazamiento externo y sin asumir ninguna forma de la ley cohesiva.
 +
 
 +
==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 GIU 16/51 en la convocatoria de 2016.
 +
 
 +
==Referencias==
 +
 
 +
<div id="1"></div>
 +
[[#cite-1|1]]<span style="text-align: center; font-size: 75%;"> [?] A. Arrese, N. Carbajal, G. Vargas, F. Mujika F. A new method for determining mode II R-curve by the End-Notched flexure test. Engineering Fracture Mechanics, 77, 2010, p.77-20. </span>
 +
 
 +
<div id="2"></div>
 +
[[#cite-2|2]]<span style="text-align: center; font-size: 75%;"> [?] Anderson TL. Fracture Mechanics. Fundamentals and applications. 2005</span>
 +
 
 +
<div id="3"></div>
 +
[[#cite-3|3]]<span style="text-align: center; font-size: 75%;"> [?] Mujika F. On the effect of shear and local deformation in three point bending tests. Polymer testing 2007;26:869-77</span>

Latest revision as of 18:32, 20 October 2022

1 Introducción

Se presenta un nuevo modelo para la determinación de la ley cohesiva en modo II de uniones adhesivas de pequeño espesor. La tasa de liberación de energía G en función del avance de grieta Δa se determina mediante el método de reducción de datos denominado Beam Theory including Bending Rotation effects (BTBR) [1]. Mediante este método, la longitud de grieta para cada punto del ensayo se determina mediante la variación de flexibilidad. El desplazamiento relativo de cortadura en la punta de grieta Δ se determina en función del avance de grieta equivalente, asumiendo que el desarrollo de la Zona de Proceso de Fractura (FPZ) es análogo al avance de grieta. Finalmente, se representa G en función de Δ y se determina la ley cohesiva τ = τ (Δ) mediante derivación numérica..

2 Metodología

La Fig. 1 muestra la configuración del ensayo de flexión en tres puntos con entalla final (ENF) de acuerdo a BTBR

Arrese et al 2018a-image5-c.png
Figura 1. Configuración del ensayo ENF de acuerdo a BTBR.

La integral-J es equivalente a la tasa de liberación de energía G para un cuerpo elástico. [2]


Arrese et al 2018a-image6-c.png
Arrese et al 2018a-image7-c.png

(1)

En la aproximación propuesta se asume que el desarrollo de la FPZ es análogo al avance de grieta equivalente Δa, relacionado con la variación de la flexibilidad.

Arrese et al 2018a-image8-c.png
Figura 2. Zona de proceso de Fractura.

De acuerdo a la teoría clásica de vigas Δ puede expresarse de la siguiente forma:

Arrese et al 2018a-image9-c.png (2)

donde

Arrese et al 2018a-image10-c.png

La relación G=G(Δ) se aproxima a una función logística de la siguiente forma:

Arrese et al 2018a-image11-c.png (3)

F(Δ) se determina mediante un ajuste lineal por regiones de los datos experimentales. Finalmente la ley cohesiva se obtiene median derivación numérica.

Arrese et al 2018a-image12-c.png (4)

3 Procedimiento Experimental

3.1 Material

Para la fabricación de las probetas se han utilizado probetas unidireccionales fabricadas a partir de prepregs de carbono/epoxi T6T/F593 proporcionados por Hexcel Composites con 16 capas unidireccionales que han sido delaminadas completamente. Los adherentes fruto de la delaminación completa han sido unidas por el adhesivo Loctite 9644 Hysol.

Las dimensiones nominales de la probeta son 3 mm de espesor y 16mm de ancho. La grieta inicial se genera mediante un inserto de teflón.

Los ensayos se han realizado en una maquina de ensayo universal MTS-Insight 10 con célula de carga de 10kN

3.2 Ensayos Preliminares

Las propiedades elásticas Ef and G13, se han determinado mediante el procedimiento propuesto por Mujika [3], siendo el modulo de flexión 107GPa y el modulo de cortadura 3 GPa.

El desplazamiento de la probeta (δspec) se determina de las curvas carga desplazamiento. El desplazamiento experimental (δexp) es la suma del desplazamiento de la probeta y el desplazamiento debido a la flexibilidad del sistema.

Arrese et al 2018a-image13-c.png

La flexibilidad del sistema se determina experimentalmente siendo su valor promedio Cs = 1/15•10-3 mm/N.

3.3 Resultados

En la Figura 3 se muestra la curva experimental carga desplazamiento correspondiente a la probeta S1.

Arrese et al 2018a-image14-c.png
Figura 3. Curva carga-desplazamiento.

Una vez determinada la longitud de grieta para cada punto del ensayo, la longitud de grieta inicial se determina como promedio de los valores en el rango de carga donde la longitud de grieta se mantiene constante, tal y como se aprecia en la Fig. 4. Por lo tanto, la iniciación del daño local está relacionada con el primer punto donde el avance de grieta es distinto de cero.

Arrese et al 2018a-image15-c.png
Figura 4. Evolución de la longitud de grieta.

Teniendo la longitud de de grieta inicial y la longitud de grieta para cada punto del ensayo, es posible obtener la curva-R, la cual representa la evolución de G respecto al avance de grieta, Fig. 5.

Arrese et al 2018a-image16-c.png
Figura 5. Curva R.

De acuerdo a la Eq. (2) el desplazamiento de cortadura en la punta de grieta puede ser determinado para cada punto del ensayo sustituyendo el correspondiente avance de grieta. La Fig 6. muestra la curva G¸Ç_ëž?Δ que corresponde a la probeta S1.

Arrese et al 2018a-image17-c.png
Figura 6. Curva.G-Δ.

Para evitar el ruido excesivo de los datos experimentales, los datos los datos G-Δ se aproximan a la función logística de la Eq (3). La función F (Δ) se ajusta a los datos experimentales mediante un ajuste lineal por regiones. Fig. 7 muestra los datos experimentales obtenidos de la Eq. (3) y la curva ajustada. Con el fin de evaluar la precisión del ajuste, los datos experimentales y la función logística aproximada se representan en la Fig. 8.

Arrese et al 2018a-image18-c.png
Figura 7. La función F.(Δ)
Arrese et al 2018a-image19-c.png
Figura 8. (G-Gini)/(Gss-Gini) experimental y aproximada.

Reemplazando los valores de F (Δ) y sus derivadas obtenidas por diferenciación numérica en la Eq. (4), se obtiene la ley cohesiva.

La Fig. 9 muestra la ley cohesiva de la probeta S1

Arrese et al 2018a-image20-c.png
Figura 9. Ley cohesiva de la probeta S1.

4 Conclusiones

En el presente trabajo se propone un nuevo método para la determinación de la ley cohesiva en modo II para uniones adhesivas de pequeño espesor.

La aproximación propuesta, se basa en la correlación entre la tasa de liberación de energía desarrollada (G) en función del desplazamiento de cortadura de la punta de grieta (Δ).

La Curva de Resistencia se determina obteniendo G en función del avance de grieta cuya determinación se basa en la variación de la flexibilidad de la probeta debido al avance del daño.

El desplazamiento de cortadura de la punta de la grieta se determina en función del avance de grieta equivalente suponiendo que el desarrollo de la Zona de Proceso de Fractura es análogo a un avance grieta.

Finalmente, representando la tasa de liberación de energía y el desplazamiento de cortadura de la grieta determinados, y mediante un ajuste a una función logística de la curva G-Δ y la ley cohesiva experimental se determina por diferenciación numérica.

Mediante el método propuesto se ha evaluado una unión realizada con LOCTITE HYSOL 9644. La resistencia cohesiva obtenida es del orden de 34 MPa,

El método propuesto proporciona una manera sencilla de obtener la ley cohesiva en modo II utilizando únicamente datos de desplazamiento y carga proporcionados por la máquina de ensayo universal, sin necesidad de ninguna medición de desplazamiento externo y sin asumir ninguna forma de la ley cohesiva.

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 GIU 16/51 en la convocatoria de 2016.

Referencias

1 [?] A. Arrese, N. Carbajal, G. Vargas, F. Mujika F. A new method for determining mode II R-curve by the End-Notched flexure test. Engineering Fracture Mechanics, 77, 2010, p.77-20.

2 [?] Anderson TL. Fracture Mechanics. Fundamentals and applications. 2005

3 [?] Mujika F. On the effect of shear and local deformation in three point bending tests. Polymer testing 2007;26:869-77

Back to Top
GET PDF

Document information

Published on 14/10/18
Accepted on 14/10/18
Submitted on 14/10/18

Volume 02 - Comunicaciones Matcomp17 (2018), Issue Núm. 4 - Comportamiento en servicio de los materiales compuestos (2), 2018
DOI: 10.23967/r.matcomp.2018.10.016
Licence: Other

Document Score

0

Views 2
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?