Mixed Mode Cohesive Law of Adhesive Joints

Latest revision as of 21:35, 26 November 2022

1 Introducción

El uso de materiales compuestos y la utilización cada vez más generalizada de adhesivos estructurales frente a otros procesos de unión tradicional implica la necesidad de un conocimiento específico del comportamiento de dichas uniones ante solicitaciones que provocan la separación de las piezas unidas.

El empleo de adhesivos estructurales en este tipo de uniones, en vez de los sistemas tradicionales, como soldadura, tornillos o remaches, no sólo implican una reducción de peso por la eliminación de estos elementos sino que supone una mejora en las características estructurales [1].

Sin embargo, el tipo de exposición que comúnmente sufren las uniones adhesivas pueden afectar a la durabilidad de las mismas [2,3]. Por tanto, la evaluación de la calidad de la unión adhesiva y la obtención de parámetros que caractericen su comportamiento mecánico ante la fractura es de vital importancia, tanto para el diseño como para el control de calidad del proceso productivo de elementos estructurales.

Las uniones adhesivas han sido intensamente investigadas durante los últimos años, dando lugar a diversas aproximaciones analíticas [4-5,67] y propuestas de ensayo para la caracterización mecánica de estas uniones.[8-910]

El modelo de zona cohesiva (CZM) es uno de los métodos disponibles para la investigación de la propagación de la grieta mediante el método de elementos finitos.

CZM fue presentado a principios de los sesenta por Barenblatt [11] y Dugdale [12] y describe los procesos locales de fractura cerca de las puntas de la grieta como un fenómeno gradual donde la separación tiene lugar a través de una zona cohesiva. Es ampliamente utilizado en los paquetes de elementos finitos comerciales debido a la aplicabilidad para analizar la fractura de diferentes materiales bajo diferentes condiciones de carga. Con CZM la fractura se modeliza como un proceso en el que una superficie en el material forma primero una zona cohesiva denominada zona de proceso de fractura (FPZ). Esta FPZ se separa posteriormente en dos superficies de grieta. En la zona cohesiva, las superficies de las grietas se mantienen unidas por tensiones cohesivas.

La idea general es que el proceso de fractura se describe por una relación de tensión - desplazamiento relativo local en la zona de proceso de fractura (FPZ). Esta relación es una ley constitutiva del material, llamada ley cohesiva.

Debido a que la ley cohesiva juega un papel importante en la simulación del comportamiento de la fractura de los materiales los métodos para determinar la ley cohesiva de uniones adhesivas han sido desarrollados por diversos autores [13-,14,15,16].

Sin embargo, estos métodos requieren el uso de equipo externos para determinar los desplazamientos relativos en la punta de la grieta, como pueden ser la Correlación Digital de Imagen (DIC) o el Transformador Diferencial De Voltaje Lineal (LVDT) [13-1415]. Estos métodos de medición incorporan dificultades experimentales asociadas, relacionadas con resultados inexactos en el caso de desplazamientos relativos de puntas de grieta muy pequeños.

En este trabajo se presenta un nuevo método para la determinación de la ley cohesiva en modo mixto de uniones adhesivas de pequeño espesor.

El método propuesto se basa en el ensayo MMB (Mixed Mode Bending test) [17].

La reducción de datos se realiza mediante la combinación de dos método recientemente presentados para modo I y modo II [18, 19]. La combinación de estos métodos basados en teoría de vigas, permite determinar la longitud de grieta equivalente para cada punto del ensayo mediante la variación de flexibilidad.

Utilizando únicamente el registro la carga y desplazamiento del punto de carga del ensayo, se determinan la tasa de liberación de energía J y el desplazamiento relativo en la punta de grieta Δ, en función del avance de grieta equivalente, asumiendo que el desarrollo de la Zona de Proceso de Fractura (FPZ) es análogo a este avance.

Finalmente, representando J en función de Δ ¸'Ø��? la ley cohesiva τ = τ (Δ) se determina mediante derivación numérica.

2 Metodología

La tasa de Liberación de Energía (J) se determina mediante la integral J, integral de contorno independiente de contorno adoptado para evaluarla, presentada por Rice [20]

$J=\int _{C}^{}\left(Wdy-T{\frac {\partial u}{\partial x}}\,dC\right)$	(1)
Donde C es el camino de integración en sentido contrario a las agujas del reloj, W es la densidad de energía de deformación, T el vector tensión y u el vector de desplazamiento. Al elegir c cerca de la punta de la grieta, T es nulo [20,21] y la Eq. (1) se convierte en:
$J=\int _{C}^{}Wdy=\int _{0}^{\Delta }{\mathit {\boldsymbol {\tau }}}\,d\Delta$	(2)

Donde τ es tensión cohesiva y Δ es el desplazamiento total en la punta de la grieta

Se utiliza el ensayo de Mixed Mode Bending MMB para analizar la fractura interlaminar en modo mixto.

El ensayo MMB es una combinación del ensayo ENF y DCB (Figura 1). Consecuentemente, el sistema de carga y desplazamiento correspondiente al ensayo MMB se puede representar como la superposición de modo I y modo II puro. Las fuerzas correspondientes al modo I y modo II son:

(3)

El desplazamiento de la palanca

(4)

De manera que la flexibilidad de total es

(5)

donde ${\textstyle {C}_{I}={\frac {{\delta }_{I}}{{P}_{I}}}\,\quad \quad \,\,\,{C}_{II}=}$ ${\frac {{\delta }_{II}}{{P}_{II}}}\,\,$

Sustituyendo las expresiones obtenidas en trabajos previos para la flexibilidad en Modo I [22] y Modo II [23] e igualando la expresión a la flexibilidad experimental, obtendremos la grieta equivalente a_e mediante un método iterativo.

Captura de pantalla 2022-11-26 202809.png

Figura 1. Configuración de ensayo MMB. Superposición de un ensayo ENF y DCB

Arratibel et al 2021a-picture-x0000 s1073.svg

Figura 2. Definición del sistema equivalente para cada punto del ensayo MMB.

Asumiendo que el desarrollo de la zona de proceso de fallo es análogo al avance de grieta equivalente determinado por la variación de la flexibilidad, se determinan las rotaciones en los puntos de aplicación de carga y se obtiene la tasa de liberación de energía [14]

${\textstyle {J}_{I}={\frac {{2P}_{I}}{w}}\left[{\theta }_{AI}\right]\,}$ ${\textstyle {J}_{II}=}$ ${\frac {{P}_{II}}{2w}}\left[{\theta }_{AII}-{2\theta }_{CII}+{\theta }_{BII}\right]$

(6)

El desplazamiento relativo total en la punta de la grieta se determina como la norma euclídea del desplazamiento relativo de apertura (modo I) [18] y tangencial (modo II) [19], determinados en función del avance de grieta equivalente.

(7)

Según la Eq. (2) si la relación entre J y Δ es conocida, suponiendo que J es una diferencial exacta, la ley cohesiva se puede obtener de acuerdo a:

(8)

3 Resultados

Los datos de carga desplazamiento que se han utilizado para realizar la validación numérica del método propuesto se han obtenido de la publicación de Turón et al. [24] y se muestran en la figura 3. Las curvas corresponden a ensayos virtuales MMB para ratios de modo mixto de %20, %50 y %80.

Figura 3. Curvas Carga-Desplazamiento de ensayos virtuales MMB para ratios de modo mixto de %20, %50 y %80 publicados por Turón et al [24].

Las muestras simuladas son de 150 mm de largo, 20 mm de ancho, con dos brazos de 1,55 mm de espesor, con una longitud inicial de grieta de 35 mm. Se utiliza una rigidez inicial de k_p = 10⁶N/mm³. En la tabla 1 se muestran las propiedades del material

Tabla 1: Propiedades del material

E₁₁ (GPa)	E₂₂=E₃₃ (GPa)	G₁₂=G₁₃ (GPa)	G₂₃ (GPa)	ν₁₂=ν₁₃	ν₂₃
120.0	10.5	5.25	3.48	0.3	0.5
σ_n (MPa)	τ_t (MPa)	J_Ic N/mm	J_IIc N/mm	η
30	58.9	0.260	1.002	2

Para poder obtener los ratios de modo mixto deseados, la longitud de aplicación de carga c se ajusta según [25]. Los ratios de modo mixto B de %20, %50 y %80 se obtienen ajustando la distancia c de la palanca en 149.05 mm, 63.18 mm y 41.06 mm, respectivamente.

Los datos de las leyes cohesivas de entrada, que se muestran en la siguiente tabla, se obtienen sustituyendo los valores correspondientes a los ratios de modo mixto fijos de %20, %50 y %80 en las expresiones definidas en [8] donde se aplica el criterio de Benzeggagh y Kenane siendo η es el parámetro de interacción de modo mixto y B el ratio de modo mixto.

Tabla 2: Leyes cohesivas de entrada correspondientes a ratios de modo mixto de %20, %50 y %80

B	J_c (N/mm)	Tensión (MPa)	Δ₀ (mm)	Δ_f (mm)
%20	0.2897	31.67	3.10^-5	0.0183
%50	0.4455	39.27	4.10^-5	0.0227
%80	0.7349	50.44	5.10^-5	0.0291

Arratibel et al 2021a 8713 arratibel fig4.png

Figura 4. Leyes cohesivas que describen las ecuaciones constitutivas del Modelo de zona cohesiva.

De acuerdo al procedimiento propuesto, la Tasa de Liberación de Energía J se obtienen para cada punto del ensayo sustituyendo en la Eq. (6) las rotaciones en los puntos de aplicación de carga determinados en función del correspondiente avance de grieta equivalente determinado utilizando únicamente datos de carga y desplazamiento de la máquina.

En las siguientes figuras se muestra la evolución de J con el desplazamiento del punto de carga δ. La J determinada mediante el método propuesto correspondiente a cada modo se compara con los resultados esperados correspondientes a ratios de modo mixto fijos de %20, %50 y %80.

(a)

(b)

(c)

Figura 5. Evolución de la tasa de liberación de energía correspondiente a cada modo de fractura para ratios de modo mixto de a)%20,b) %50 y c)%80.

El Desplazamiento Relativo total en la punta de grieta Δ se obtienen para cada punto del ensayo sustituyendo en la Eq. (7) el desplazamiento relativo de apertura Δ_n [18] y el tangencial Δ_t [19], determinados en función del avance de grieta equivalente.

La Fig 4. Se muestra la curva J-Δ que corresponde a ratios de modo mixto de %20, %50 y %80

Figura 6. Evolución J-Δ para ratio de modo mixto de %20, %50 y %80.

Según la Eq. (2) si la relación entre J y Δ es conocida, se obtienen las leyes cohesivas mediante derivación numérica Eq. (8).

Figura 7. Ley cohesiva para ratio de modo mixto de %20, %50 y %80.

4 Conclusiones

En el presente trabajo se propone un método para la determinación de la ley cohesiva en modo mixto 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 (J) en función del desplazamiento relativo de la punta de grieta (Δ).

La Curva de Resistencia se determina obteniendo J 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 relativo 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 relativo de la punta de la grieta determinados, la ley cohesiva experimental se determina por diferenciación numérica.

El método propuesto proporciona una manera sencilla de obtener la ley cohesiva 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] T.A Barnes, I.R. Pashby, Journal of Materials Processing Technology. 99, pág. 72-79 (2000). https://doi.org/10.1016/S0924-0136(99)00361-1

[2] M. Fernando, WW. Harjoprayitbi, AJ. Kinloch, International Journal of Adhesion and Adhesives 16, 2, pág. 113-119, (1996). https://doi.org/10.1016/0143-7496(96)89799-5

[3] J.Jumel, MK. Budzik, N Ben Salem, MER Shanahan, International Journal Of Solids an Structures, 50, 2, pág. 297-309, (2013). https://doi.org/10.1016/j.ijsolstr.2012.08.028

[4] G. Fernlund, JK. Spelt, Composites Science and Technology, 50, 4, pag. 441–449, (1994). https://doi.org/10.1016/0266-3538(94)90052-3

[5] LFM da Silva, PJC das Neves, RD Adams, JK Spelt, International Journal of Adhesion and Adhesives, 29, 3, pág. 319-330, (2009). https://doi.org/10.1016/j.ijadhadh.2008.06.005

[6] RQ. Rodríguez, W. Portilho de Paiva, P. Sollero, MRB. Rodrigues, E. Lima de Albuquerque, International Journal of Adhesion and Adhesives, 37, pág. 26–36,( 2012) https://doi.org/10.1016/j.ijadhadh.2012.01.009

[7] BRK Blackman, AJ Kinloch, M Paraschi. Engineering Fracture Mechanics, 72, pág. 877-897, (2005) https://doi.org/10.1016/j.engfracmech.2004.08.007

[8] WS. Kim, JJ. Lee, Engineering Fracture Mechanics, 76, 12, pág. 1785–1799, (2009) https://doi.org/10.1016/j.engfracmech.2009.03.014

[9] MV. Fernández, MFSF de Moura, LFM. da Silva, AT. Marques, Composite Structures, 95, pág. 222–226, (2013) https://doi.org/10.1016/j.compstruct.2012.07.031

[10] H. Osnes, D. McGeorge, Composites Part B: Engineering, 40, pág. 29–40, (2009) https://doi.org/10.1016/j.compositesb.2008.07.002

[11] GI. Barenblatt, Advances in Applied Mechanics, 7, pág. 55-129, (1962) https://doi.org/10.1016/S0065-2156(08)70121-2

[12] DS. Dugdale, Journal of the Mechanics and Physics of Solids, 8, 2, pág.100–104, (1960) https://doi.org/10.1016/0022-5096(60)90013-2

[13 ] K. Leffler, KS. Alfredsson, U. Stigh, International Journal of Solids and Structures, 44, pág. 530- 545, (2007) https://doi.org/10.1016/j.ijsolstr.2006.04.036

[14] C. Sarrado, A. Turon, J. Costa, J. Renart, Composites Part A: Applied Science and Manufacturing, 90, pág. 234–242 (2016) https://doi.org/10.1016/j.compositesa.2016.07.004

[15] T. Andersson, U. Stigh, International Journal of Solids and Structures, 41, 2, pág. 413–434, (2004) https://doi.org/10.1016/j.ijsolstr.2003.09.039

[16] G. Ji, Z. Ouyang, G. Li, S. Ibekwe, SS. Pang, International Journal of Solids and Structures, 47, 18-19, pág. 2445-2458, (2010) https://doi.org/10.1016/j.ijsolstr.2010.05.006

[17] JR. Reeder, JH Crews, AIAA Journal, 28, 7, pág. 1270-1276 (1998)

[18] A. Arrese, A. Boyano, J. de Gracia, F. Mujika, Composites Science and Technology, 152, pág. 76- 84, (2017) https://doi.org/10.1016/j.compscitech.2017.09.012

[19] A. Arrese, N. Insausti, F. Mujika, M. Perez-Galmés, J. Renart, Composites Science and Technology, 170, pág. 42–50, (2019) https://doi.org/10.1016/j.compscitech.2018.11.031

[20] JR. Rice, Journal of Applied Mechanics, 35, 2, pág. 379–386, (1968) https://doi.org/doi:10.1115/1.3601206

[21] Z. Suo, G. Bao, B. Fan, Journal of the Mechanics and Physics of Solids, 40, 1, pág. 1–16, (1992) https://doi.org/10.1016/0022-5096(92)90198-B

[22] J. de Gracia, A. Boyano, A. Arrese, F. Mujika, Engineering Fracture Mechanics, 135, pág. 274-285, (2015) http://dx.doi.org/10.1016/j.engfracmech.2015.01.016

[23] A. Arrese, N. Carbajal, G. Vargas, F. Mujika, Engineering Fracture Mechanics, 77, 1, pág. 51-70, (2010) https://doi.org/10.1016/j.engfracmech.2009.09.008

[24] A. Turon, PP. Camanho, J. Costa, J. Renart, Composite Structures, 92, 8, pág. 1857–1864, (2010) https://doi.org/10.1016/j.compstruct.2010.01.012

[25] N. Blanco, A. Turon, J. Costa, Composites Science and Technology, 66, 10, pág.1256–1258, (2006) https://doi.org/10.1016/j.compscitech.2005.10.028

Latest revision as of 21:35, 26 November 2022

1 Introducción

2 Metodología

Donde τ es tensión cohesiva y Δ es el desplazamiento total en la punta de la grieta

3 Resultados

4 Conclusiones

Agradecimientos

Referencias

Document information

Document Score

Share this document

Keywords

claim authorship

@@ Line 1: / Line 1: @@
-== Abstract ==
+==1 Introducción ==
-This work presents an experimental method to determine the cohesive law of thin adhesive joints in mixed mode.
+El uso de materiales compuestos y la utilización cada vez más generalizada de adhesivos estructurales frente a otros procesos de unión tradicional implica la necesidad de un conocimiento específico del comportamiento de dichas uniones ante solicitaciones que provocan la separación de las piezas unidas.
-The proposed method is based on the MMB (Mixed Mode Bending test).
-The Energy release Rate and Crack Tip Relative Displacement are determined based on a compliance variation data reduction method by processing the global load displacement curve, without monitoring the crack length and the crack tip relative displacement during the test.
-Representing the energy release rate versus relative displacement the cohesive law is determined by numerical differentiation.
-The method is validated numerically by virtual tests including the cohesive zone model.
-The results reveal that the proposed data reduction scheme is suitable to obtain the mixed mode cohesive law using only the load and displacement data obtained from the testing machine, without any external displacement measurement technique and without any assumption of the form of the cohesive law.
-== Full document ==
+El empleo de adhesivos estructurales en este tipo de uniones, en vez de los sistemas tradicionales, como soldadura, tornillos o remaches, no sólo implican una reducción de peso por la eliminación de estos elementos sino que supone una mejora en las características estructurales [<span id='cite-1'></span>[[#1|1]]].
-<pdf>Media:Draft_Content_635036536native-articles-304.pdf</pdf>
+Sin embargo, el tipo de exposición que comúnmente sufren las uniones adhesivas pueden afectar a la durabilidad de las mismas [<span id='cite-2'></span>[[#2|2]],<span id='cite-3'></span>[[#3|3]]]. Por tanto, la evaluación de la calidad de la unión adhesiva y la obtención de parámetros que caractericen su comportamiento mecánico ante la fractura es de vital importancia, tanto para el diseño como para el control de calidad del proceso productivo de elementos estructurales.
+<span id='_Ref413929177'></span>Las uniones adhesivas han sido intensamente investigadas durante los últimos años, dando lugar a diversas aproximaciones analíticas [<span id='cite-4'></span>[[#4|4]]-<span id='cite-5'></span>[[#5|5]],<span id='cite-6'></span>[[#6|6]]<span id='cite-7'></span>[[#7|7]]] y propuestas de ensayo para la caracterización mecánica de estas uniones.[<span id='cite-8'></span>[[#8|8]]-<span id='cite-9'></span>[[#9|9]]<span id='cite-10'></span>[[#10|10]]]
+El modelo de zona cohesiva (CZM) es uno de los métodos disponibles para la investigación de la propagación de la grieta mediante el método de elementos finitos.
+CZM fue presentado a principios de los sesenta por Barenblatt [<span id='cite-11'></span>[[#11|11]]] y Dugdale [<span id='cite-12'></span>[[#12|12]]] y describe los procesos locales de fractura cerca de las puntas de la grieta como un fenómeno gradual donde la separación tiene lugar a través de una zona cohesiva. Es ampliamente utilizado en los paquetes de elementos finitos comerciales debido a la aplicabilidad para analizar la fractura de diferentes materiales bajo diferentes condiciones de carga. Con CZM la fractura se modeliza como un proceso en el que una superficie en el material forma primero una zona cohesiva denominada zona de proceso de fractura (FPZ). Esta FPZ se separa posteriormente en dos superficies de grieta. En la zona cohesiva, las superficies de las grietas se mantienen unidas por tensiones cohesivas.
+La idea general es que el proceso de fractura se describe por una relación de tensión - desplazamiento relativo local en la zona de proceso de fractura (FPZ). Esta relación es una ley constitutiva del material, llamada ley cohesiva.
+<span id='_Ref476748718'></span><span id='_Ref477195745'></span><span id='_Ref477195653'></span><span id='_Ref477195755'></span><span id='_Ref477195659'></span>Debido a que la ley cohesiva juega un papel importante en la simulación del comportamiento de la fractura de los materiales  los métodos para determinar la ley cohesiva de uniones adhesivas han sido desarrollados por diversos autores [<span id='cite-13'></span>[[#13|13]]-,<span id='cite-14'></span>[[#14|14]],<span id='cite-15'></span>[[#15|15]],<span id='cite-16'></span>[[#16|16]]].
+Sin embargo, estos métodos requieren el uso de equipo externos para determinar los desplazamientos relativos en la punta de la grieta, como pueden ser la Correlación Digital de Imagen (DIC) o el Transformador Diferencial De Voltaje Lineal (LVDT) [<span id='cite-_Ref477195745'></span>[[#_Ref477195745|13]]-<span id='cite-_Ref477195653'></span>[[#_Ref477195653|14]]<span id='cite-_Ref477195755'></span>[[#_Ref477195755|15]]]. Estos métodos de medición incorporan dificultades experimentales asociadas, relacionadas con resultados inexactos en el caso de desplazamientos relativos de puntas de grieta muy pequeños.
+En este trabajo se presenta un nuevo método para la determinación de la ley cohesiva en modo mixto de uniones adhesivas de pequeño espesor.
+El método propuesto se basa en el ensayo MMB (Mixed Mode Bending test) [<span id='cite-17'></span>[[#17|17]]].
+<span id='_Ref2606677'></span><span id='_Ref2606714'></span>La reducción de datos se realiza mediante la combinación de dos método recientemente presentados para modo I y modo II [<span id='cite-18'></span>[[#18|18]], <span id='cite-19'></span>[[#19|19]]]. La combinación de estos métodos basados en teoría de vigas, permite determinar la longitud de grieta equivalente para cada punto del ensayo mediante la variación de flexibilidad.
+Utilizando únicamente el registro la carga y desplazamiento del punto de carga del ensayo, se determinan la tasa de liberación de energía ''J'' y el desplazamiento relativo en la punta de grieta ''&#x0394;'', en función del avance de grieta equivalente, asumiendo que el desarrollo de la Zona de Proceso de Fractura (FPZ) es análogo a este avance.
+Finalmente, representando J en función de &#x0394;&#x00a0;¸'Ø��&#x00a0;la ley cohesiva &#x03c4;&#x00a0;=&#x00a0;&#x03c4; (&#x0394;) se determina mediante derivación numérica.
+==2 Metodología ==
+<span id='_Ref512346843'></span>La tasa de Liberación de Energía (''J'') se determina mediante la integral ''J'', integral de contorno independiente de contorno adoptado para evaluarla, presentada por Rice [<span id='cite-20'></span>[[#20|20]]]
+{| style="width: 100%;border-collapse: collapse;"
+|-
+|  style="text-align: center;"|<math>J=\int_{C}^{}\left( Wdy-T\frac{\partial u}{\partial x}\, dC\right)</math>
+|  style="text-align: right;"|(1)
+|-
+|  colspan='2' |<span id='_Ref487194928'></span>Donde '''C''' es el camino de integración en sentido contrario a las agujas del reloj, '''W''' es la densidad de energía de deformación, '''T''' el vector tensión y '''u''' el vector de desplazamiento. Al elegir '''c''' cerca de la punta de la grieta, '''T''' es nulo [<span id='cite-_Ref512346843'></span>[[#_Ref512346843|20]],<span id='cite-21'></span>[[#21|21]]] y la Eq. (1) se convierte en:
+|-
+|  style="text-align: center;"|<math>J=\int_{C}^{}Wdy=\int_{0}^{\Delta }\mathit{\boldsymbol{\tau }}\, d\Delta</math>
+|  style="text-align: center;"|(2)
+|}
+==Donde &#x03c4;&#x00a0;es tensión cohesiva y &#x0394; es el desplazamiento total en la punta de la grieta ==
+Se utiliza el ensayo de ''Mixed Mode Bending'' MMB para analizar la fractura interlaminar en modo mixto.
+El ensayo MMB  es una combinación del ensayo ENF y DCB (Figura 1). Consecuentemente, el sistema de carga y desplazamiento correspondiente al ensayo MMB se puede representar como la superposición de modo I y modo II puro. Las fuerzas correspondientes al modo I y modo II son:
+{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="text-align: center;"|<math>{\quad \quad P}_{I}=\frac{3c-L}{4L}P\quad \quad \, \, {P}_{II}=\frac{c+L}{L}P</math>
+|  style="text-align: right;"|(3)
+|}
+El desplazamiento de la palanca
+{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="text-align: center;"|<math>\delta ={\delta }_{I}\left[ \frac{3c-L}{4L}\right] +{\delta }_{II}\left[ \frac{c+L}{L}\right]</math>
+|  style="text-align: right;"|(4)
+|}
+De manera que la flexibilidad de total es
+{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="text-align: center;"|<math display="inline">C=\frac{\delta }{P}={C}_{I}{\left[ \frac{3c-L}{4L}\right] }^{2}+</math><math>{C}_{II}{\left[ \frac{c+L}{L}\right] }^{2}</math>
+|  style="text-align: right;"|(5)
+|}
+donde  <math display="inline">{C}_{I}=\frac{{\delta }_{I}}{{P}_{I}}\, \quad \quad \, \, \, {C}_{II}=</math><math>\frac{{\delta }_{II}}{{P}_{II}}\, \,</math>
+Sustituyendo las expresiones obtenidas en trabajos previos para la flexibilidad en Modo I [<span id='cite-22'></span>[[#22|22]]] y Modo II [<span id='cite-23'></span>[[#23|23]]] e igualando la expresión a la flexibilidad experimental, obtendremos la grieta equivalente ''a<sub>e</sub>'' mediante un método iterativo.
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">[[File:Captura de pantalla 2022-11-26 202809.png|centre|713x713px]]</span></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 de ensayo MMB. Superposición de un ensayo ENF y DCB </span></div>
+[[Image:Arratibel_et_al_2021a-picture-x0000_s1073.svg|center|600x600px]]
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">'''Figura 2. '''Definición del sistema equivalente para cada punto del ensayo MMB.</span></div>
+Asumiendo que el desarrollo de la zona de proceso de fallo es análogo al avance de grieta equivalente determinado por la variación de la flexibilidad, se determinan las rotaciones en los puntos de aplicación de carga y se obtiene la tasa de liberación de energía [<span id='cite-_Ref477195653'></span>[[#_Ref477195653|14]]]
+{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="text-align: center;"|<math>J={J}_{I}+{J}_{II}</math>
+<math display="inline">{J}_{I}=\frac{{2P}_{I}}{w}\left[ {\theta }_{AI}\right] \,</math>              <math display="inline">{J}_{II}=</math><math>\frac{{P}_{II}}{2w}\left[ {\theta }_{AII}-{2\theta }_{CII}+{\theta }_{BII}\right]</math>
+|  style="text-align: right;"|(6)
+|}
+El desplazamiento relativo total en la punta de la grieta se determina como la norma euclídea del desplazamiento relativo de apertura (modo I) [<span id='cite-_Ref2606677'></span>[[#_Ref2606677|18]]] y tangencial (modo II) [<span id='cite-_Ref2606714'></span>[[#_Ref2606714|19]]], determinados en función del avance de grieta equivalente.
+{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="text-align: center;"|<math>\Delta =\sqrt{{{\Delta }_{n}}^{2}+{{\Delta }_{t}}^{2}}</math>
+|  style="text-align: right;"|(7)
+|}
+Según la Eq. (2) si la relación entre ''J'' ''y'' ''&#x0394;'' es conocida, suponiendo que ''J'' es una diferencial exacta, la ley cohesiva se puede obtener de acuerdo a:
+{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="text-align: center;"|<math>\tau \left( \Delta \right) =\frac{\partial J}{\partial \Delta }\quad \quad \quad</math>
+|  style="text-align: right;"|(8)
+|}
+==3 Resultados ==
+<span id='_Ref185029'></span>Los datos de carga desplazamiento que se han utilizado para realizar la validación numérica del método propuesto se han obtenido de la publicación de Turón et al. [<span id='cite-24'></span>[[#24|24]]]  y se muestran en la figura 3.  Las curvas corresponden a ensayos virtuales MMB para ratios de modo mixto de %20, %50 y %80.
+[[Image:Arratibel_et_al_2021a_9053_arratibel_fig3.png|center|587x587px]]
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">'''Figura 3'''. Curvas Carga-Desplazamiento de ensayos virtuales MMB para ratios de modo mixto de %20, %50 y %80 publicados por Turón et al [</span><span id='cite-_Ref185029'></span>[[#_Ref185029|<span style="text-align: center; font-size: 75%;">24</span>]]<span style="text-align: center; font-size: 75%;">].</span></div>
+Las muestras simuladas son de 150 mm de largo, 20 mm de ancho, con dos brazos de 1,55 mm de espesor, con una longitud inicial de grieta de 35 mm. Se utiliza una rigidez inicial de k<sub>p</sub> = 10<sup>6 </sup>N/mm<sup>3</sup>. En la tabla 1 se muestran las propiedades del material
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">Tabla 1: Propiedades del material</span></div>
+{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="border-top: 1pt solid black;text-align: center;"|''E<sub>11</sub>''
+(GPa)
+|  style="border-top: 1pt solid black;text-align: center;"|''E<sub>22</sub>=E<sub>33</sub>''
+(GPa)
+|  style="border-top: 1pt solid black;text-align: center;"|''G<sub>12</sub>=G<sub>13</sub>''
+(GPa)
+|  style="border-top: 1pt solid black;text-align: center;"|''G<sub>23</sub>''
+(GPa)
+|  style="border-top: 1pt solid black;text-align: center;"|''ν<sub>12</sub>=ν<sub>13</sub>''
+|  style="border-top: 1pt solid black;text-align: center;"|''ν<sub>23</sub>''
+|-
+|  style="border-bottom: 1pt solid black;text-align: center;"|120.0
+|  style="border-bottom: 1pt solid black;text-align: center;"|10.5
+|  style="border-bottom: 1pt solid black;text-align: center;"|5.25
+|  style="border-bottom: 1pt solid black;text-align: center;"|3.48
+|  style="border-bottom: 1pt solid black;text-align: center;"|0.3
+|  style="border-bottom: 1pt solid black;text-align: center;"|0.5
+|-
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''&#x03c3;<sub>n</sub>''
+(MPa)
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''&#x03c4;<sub>t</sub>''
+(MPa)
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''J<sub>Ic</sub>''
+N/mm
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''J<sub>IIc</sub>''
+N/mm
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''&#x03b7;''
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|
+|-
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|30
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|58.9
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|0.260
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|1.002
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|2
+|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|
+|}
+Para poder obtener los ratios de modo mixto deseados, la longitud de aplicación de carga ''c'' se ajusta según [<span id='cite-25'></span>[[#25|25]]]. Los ratios de modo mixto ''B'' de %20, %50 y %80 se obtienen ajustando la distancia ''c'' de la palanca en 149.05 mm, 63.18 mm y 41.06 mm, respectivamente.
+Los datos de las leyes cohesivas de entrada, que se muestran en la siguiente tabla, se obtienen sustituyendo los valores correspondientes a los ratios de modo mixto fijos de %20, %50 y %80 en las expresiones definidas en [8] donde se aplica  el criterio de Benzeggagh y Kenane siendo ''&#x03b7;'' es el parámetro de interacción de modo mixto y ''B'' el ratio de modo mixto.
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">Tabla 2: Leyes cohesivas de entrada correspondientes  a ratios de modo mixto de %20, %50 y %80</span></div>
+{| style="width: 84%;margin: 1em auto 0.1em auto;border-collapse: collapse;"
+|-
+|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">''B''</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">''J''<sub>c</sub></span>
+<span style="text-align: center; font-size: 75%;">(N/mm)</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Tensión </span>
+<span style="text-align: center; font-size: 75%;">(MPa)</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x0394;''<sub>0</sub></span>
+<span style="text-align: center; font-size: 75%;">(mm)</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x0394;<sub>f</sub>''</span>
+<span style="text-align: center; font-size: 75%;">(mm)</span>
+|-
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">%20</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2897 </span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">31.67</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.10<sup>-5</sup></span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0183</span>
+|-
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">%50</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.4455 </span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">39.27</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.10<sup>-5</sup></span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0227</span>
+|-
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">%80</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.7349 </span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">50.44</span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.10<sup>-5</sup></span>
+|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.0291</span>
+|}
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+[[Image:Arratibel_et_al_2021a_8713_arratibel_fig4.png|center]]
+</div>
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">'''Figura 4. '''Leyes cohesivas que describen las ecuaciones constitutivas del Modelo de zona cohesiva.</span></div>
+De acuerdo al procedimiento propuesto, la Tasa de Liberación de Energía ''J'' se obtienen para cada punto del ensayo sustituyendo en la Eq. (6) las rotaciones en los puntos de aplicación de carga determinados en función del correspondiente avance de grieta equivalente determinado utilizando únicamente datos de carga y desplazamiento de la máquina.
+En las siguientes figuras se muestra la evolución de ''J'' con el desplazamiento del punto de carga &#x03b4;. La ''J'' determinada mediante el método propuesto correspondiente a cada modo se compara con los resultados esperados correspondientes a ratios de modo mixto fijos de %20, %50 y %80.
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+[[File:Arratibel_et_al_2021a_9272_arratibel_fig5a.png|center|500px]]</div>
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">(a)</span></div>
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+[[File:Arratibel_et_al_2021a_1268_arratibel_fig5b.png|center|500px]]
+</div>
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">(b)</span></div>
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+[[File:Arratibel_et_al_2021a_3470_arratibel_fig5c.png|center|500px]]
+ </div>
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">(c)</span></div>
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">'''Figura 5.''' Evolución de la tasa de liberación de energía correspondiente a cada modo de fractura para ratios de modo mixto de a)%20,b) %50 y c)%80.</span></div>
+El Desplazamiento Relativo total en la punta de grieta ''&#x0394;'' se obtienen para cada punto del ensayo sustituyendo en la Eq. (7) el desplazamiento relativo de apertura ''&#x0394;<sub>n</sub>'' [<span id='cite-_Ref2606677'></span>[[#_Ref2606677|18]]] y el tangencial ''&#x0394;<sub>t</sub>'' [<span id='cite-_Ref2606714'></span>[[#_Ref2606714|19]]], determinados en función del avance de grieta equivalente.
+La Fig 4. Se muestra la curva J-&#x0394; que corresponde a ratios de modo mixto de %20, %50 y %80
+[[File:Arratibel_et_al_2021a_5309_arratibel_fig6.png|center|601x601px]]
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">'''Figura 6.''' Evolución J-&#x0394;&#x00a0;para ratio de modo mixto de %20, %50 y %80.</span></div>
+Según la Eq. (2) si la relación entre ''J'' ''y'' ''&#x0394;'' es conocida, se obtienen las leyes cohesivas mediante derivación numérica Eq. (8).
+[[File:Arratibel_et_al_2021a_3381_arratibel_fig7.png|center|606x606px]]
+<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
+<span style="text-align: center; font-size: 75%;">'''Figura 7.''' Ley cohesiva para ratio de modo mixto de %20, %50 y %80.</span></div>
+==4 Conclusiones==
+En el presente trabajo se propone un método para la determinación de la ley cohesiva en modo mixto 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 (''J'') en función del desplazamiento relativo de la punta de grieta (''&#x0394;'').
+La Curva de Resistencia se determina obteniendo J 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 relativo 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 relativo de la punta de la grieta determinados, la ley cohesiva experimental se determina por diferenciación numérica.
+El método propuesto proporciona una manera sencilla de obtener la ley cohesiva 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]]] T.A Barnes, I.R. Pashby, ''Journal of Materials Processing Technology''. '''99,''' pág. 72-79 (2000). [https://doi.org/10.1016/S0924-0136 https://doi.org/10.1016/S0924-0136](99)00361-1
+<div id="2"></div>
+[[#cite-2|[2]]] M. Fernando, WW. Harjoprayitbi, AJ. Kinloch, ''International Journal of Adhesion and Adhesives '' '''16''', 2, pág. 113-119, (1996).  [https://doi.org/10.1016/0143-7496 https://doi.org/10.1016/0143-7496](96)89799-5
+<div id="3"></div>
+[[#cite-3|[3]]] J.Jumel, MK. Budzik, N Ben Salem, MER Shanahan, ''International Journal Of Solids an Structures'', '''50''', 2, pág. 297-309, (2013). [https://doi.org/10.1016/j.ijsolstr.2012.08.028 https://doi.org/10.1016/j.ijsolstr.2012.08.028]
+<div id="4"></div>
+<span style="text-align: center; font-size: 75%;">[[#cite-4|[4]]] G. Fernlund, JK. Spelt, ''Composites Science and Technology'', '''50''', 4, pag. 441–449, (1994).</span> <span style="text-align: center; font-size: 75%;">[https://doi.org/10.1016/0266-3538 https://doi.org/10.1016/0266-3538](94)90052-3</span>
+<div id="5"></div>
+[[#cite-5|[5]]] LFM da Silva, PJC das Neves, RD Adams, JK Spelt, ''International Journal of Adhesion and Adhesives'', '''29''', 3, pág. 319-330, (2009). [https://doi.org/10.1016/j.ijadhadh.2008.06.005 https://doi.org/10.1016/j.ijadhadh.2008.06.005]
+<div id="6"></div>
+[[#cite-6|[6]]] RQ. Rodríguez, W. Portilho de Paiva, P. Sollero, MRB. Rodrigues, E. Lima de Albuquerque, ''International Journal of Adhesion and Adhesives'', '''37''', pág. 26–36,( 2012) [https://doi.org/10.1016/j.ijadhadh.2012.01.009 https://doi.org/10.1016/j.ijadhadh.2012.01.009]
+<div id="7"></div>
+[[#cite-7|[7]]] BRK  Blackman, AJ Kinloch, M Paraschi. ''Engineering Fracture Mechanics'', '''72''', pág. 877-897, (2005) [https://doi.org/10.1016/j.engfracmech.2004.08.007 https://doi.org/10.1016/j.engfracmech.2004.08.007]
+<div id="8"></div>
+[[#cite-8|[8]]] WS. Kim, JJ. Lee, ''Engineering Fracture Mechanics'', '''76''', 12, pág. 1785–1799, (2009) [https://doi.org/10.1016/j.engfracmech.2009.03.014 https://doi.org/10.1016/j.engfracmech.2009.03.014]
+<div id="9"></div>
+[[#cite-9|[9]]] MV. Fernández, MFSF de Moura, LFM. da Silva, AT. Marques, ''Composite Structures'', '''95''', pág. 222–226, (2013) [https://doi.org/10.1016/j.compstruct.2012.07.031 https://doi.org/10.1016/j.compstruct.2012.07.031]
+<div id="10"></div>
+[[#cite-10|[10]]] H. Osnes, D. McGeorge, ''Composites Part B: Engineering'', '''40''', pág. 29–40, (2009) [https://doi.org/10.1016/j.compositesb.2008.07.002 https://doi.org/10.1016/j.compositesb.2008.07.002]
+<div id="11"></div>
+[[#cite-11|[11]]] GI. Barenblatt, ''Advances in Applied Mechanics,'' 7, pág. 55-129, (1962) [https://doi.org/10.1016/S0065-2156 https://doi.org/10.1016/S0065-2156](08)70121-2
+<div id="12"></div>
+<span style="text-align: center; font-size: 75%;">[[#cite-12|[12]]] DS. Dugdale, ''Journal of the Mechanics and Physics of Solids'', '''8''', 2, pág.100–104, (1960)</span> <span style="text-align: center; font-size: 75%;">[https://doi.org/10.1016/0022-5096 https://doi.org/10.1016/0022-5096](60)90013-2</span>
+<div id="13"></div>
+<span style="text-align: center; font-size: 75%;">[[#cite-13|[13]] ] K. Leffler, KS. Alfredsson, U. Stigh, ''International Journal of Solids and Structures'', '''44''', pág. 530- 545, (2007)</span> <span style="text-align: center; font-size: 75%;">[https://doi.org/10.1016/j.ijsolstr.2006.04.036 https://doi.org/10.1016/j.ijsolstr.2006.04.036]</span>
+<div id="14"></div>
+[[#cite-14|[14]]] C. Sarrado, A. Turon, J. Costa, J. Renart, ''Composites Part A: Applied Science and Manufacturing'', '''90''', pág. 234–242 (2016) [https://doi.org/10.1016/j.compositesa.2016.07.004 https://doi.org/10.1016/j.compositesa.2016.07.004]
+<div id="15"></div>
+[[#cite-15|[15]]]  T. Andersson, U. Stigh'', International Journal of Solids and Structures'', '''41''', 2, pág. 413–434, (2004) [https://doi.org/10.1016/j.ijsolstr.2003.09.039 https://doi.org/10.1016/j.ijsolstr.2003.09.039]
+<div id="16"></div>
+[[#cite-16|[16]]] G. Ji, Z. Ouyang, G. Li, S. Ibekwe, SS. Pang, ''International Journal of Solids and Structures'', '''47''', 18-19, pág. 2445-2458, (2010) [https://doi.org/10.1016/j.ijsolstr.2010.05.006 https://doi.org/10.1016/j.ijsolstr.2010.05.006]
+<div id="17"></div>
+[[#cite-17|[17]]] JR. Reeder, JH Crews, ''AIAA Journal, '' '''28''', 7, pág. 1270-1276 (1998)
+<div id="18"></div>
+[[#cite-18|[18]]] A. Arrese, A. Boyano, J. de Gracia, F. Mujika, ''Composites Science and Technology,'' '''152''', pág. 76- 84, (2017) [https://doi.org/10.1016/j.compscitech.2017.09.012 https://doi.org/10.1016/j.compscitech.2017.09.012]
+<div id="19"></div>
+[[#cite-19|[19]]] A. Arrese, N. Insausti, F. Mujika, M. Perez-Galmés, J. Renart, ''Composites Science and Technology,'' '''170''', pág. 42–50, (2019) [https://doi.org/10.1016/j.compscitech.2018.11.031 https://doi.org/10.1016/j.compscitech.2018.11.031]
+<div id="20"></div>
+<span style="text-align: center; font-size: 75%;">[[#cite-20|[20]]] JR. Rice, ''Journal of Applied Mechanics'', '''35''', 2, pág. 379–386, (1968)</span> h<span style="text-align: center; font-size: 75%;">ttps://doi.org/doi:10.1115/1.3601206</span>
+<div id="21"></div>
+[[#cite-21|[21]]] Z. Suo, G. Bao, B. Fan, ''Journal of the Mechanics and Physics of Solids'', '''40''', 1, pág. 1–16, (1992) [https://doi.org/10.1016/0022-5096 https://doi.org/10.1016/0022-5096](92)90198-B
+<div id="22"></div>
+[[#cite-22|[22]]] '' ''J. de Gracia, A. Boyano, A. Arrese, F. Mujika, ''Engineering Fracture Mechanics'', '''135''', pág. 274-285, (2015) [http://dx.doi.org/10.1016/j.engfracmech.2015.01.016 http://dx.doi.org/10.1016/j.engfracmech.2015.01.016]
+<div id="23"></div>
+[[#cite-23|[23]]] '' ''A. Arrese, N. Carbajal, G. Vargas, F. Mujika'', Engineering Fracture Mechanics, 77, 1, pág. 51-70, (2010)'' [https://doi.org/10.1016/j.engfracmech.2009.09.008 https://doi.org/10.1016/j.engfracmech.2009.09.008]
+<div id="24"></div>
+[[#cite-24|[24]]] A. Turon, PP. Camanho, J. Costa, J. Renart, ''Composite Structures'',''' 92''', 8, pág. 1857–1864, (2010) [https://doi.org/10.1016/j.compstruct.2010.01.012 https://doi.org/10.1016/j.compstruct.2010.01.012]
+<div id="25"></div>
+[[#cite-25|[25]]]  N. Blanco, A. Turon, J. Costa, ''Composites Science and Technology, '' '''66''', 10, pág.1256–1258, (2006) [https://doi.org/10.1016/j.compscitech.2005.10.028 https://doi.org/10.1016/j.compscitech.2005.10.028]