(Tag: Visual edit) |
|||
Line 408: | Line 408: | ||
En la figura se muestra un ensayo ADCB de una probeta de dimensiones aproximadas a las utilizadas por Mollón et al [2]. Las propiedades del material corresponden al AS4/8552 de Hexcel Composites. La probeta tiene una longitud total de 150 mm, la grieta es de 50 mm, y el espesor de las partes superior e inferior es 1 y 5 mm, respectivamente. Se impone un desplazamiento de 5 mm a la parte superior. Puede apreciarse la inclinación de la probeta debida a la asimetría de los desplazamientos de apertura. | En la figura se muestra un ensayo ADCB de una probeta de dimensiones aproximadas a las utilizadas por Mollón et al [2]. Las propiedades del material corresponden al AS4/8552 de Hexcel Composites. La probeta tiene una longitud total de 150 mm, la grieta es de 50 mm, y el espesor de las partes superior e inferior es 1 y 5 mm, respectivamente. Se impone un desplazamiento de 5 mm a la parte superior. Puede apreciarse la inclinación de la probeta debida a la asimetría de los desplazamientos de apertura. | ||
− | [[File:Mujika_et_al_2017a-image5-c.png|centre| | + | [[File:Mujika_et_al_2017a-image5-c.png|centre|623x623px]] |
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> |
Los laminados compuesto que se utilizan en aplicaciones estructurales consisten en muchos casos en laminados constituidos por capas. En cada una de ellas el refuerzo es unidireccional y las propiedades mecánicas del laminado dependen de los materiales de las capas, de la orientación del refuerzo en cada una de ellas, así como del orden en que se apilan las mismas. La ausencia de refuerzo en la dirección perpendicular al laminado está relacionada con la fractura interlaminar, que es una de las limitaciones en la utilización de estos materiales.
Por ejemplo, un impacto de baja energía puede provocar una pequeña deslaminación entre capas. Si una de las partes deslaminadas está sometida a compresión, puede generarse pandeo local que provoca un avance de grieta que puede derivar en el fallo catastrófico del componente.
Para caracterizar la fractura interlaminar, es habitual determinar experimientalmente la tasa crítica de liberación de energía , fundamentalmente en modos I, II y modo mixto I/II. Se realizan también ensayos en modo III, que no son objetod del presente trabajo. Las probetas utilizadas tienen geometría de tira y existen diferentes ensayos para determinar . Se describirán únicamente los más utilizados. Para determinar GI se utiliza el ensayo de Viga en Doble Voladizo (Double Cantilever Beam, DCB). Para determinar el modo II se utiliza el ensayo de Flexión con Entalla Final (End Notched Flexure, ENF). Para determinar el modo mixto se utiliza el ensayo de flexión en modo mixto (Mixed Mode Bending, MMB) que es una combinación de los dos métodos anteriores.
En los ensayos DCB, ENF y MMB se utilizan probetas que tienen una grieta inicial generada en el proceso de fabricación. Esta grieta está situada en el plano central del laminado, por lo que la configuración de probeta es simétrica en todos los casos. En el ensayo de viga asimétrica en doble voladizo (Asymmetric Double Cantilever Beam, ADCB) los brazos de grieta son de distintos espesores y se produce modo mixto I/II. En la mayoría de los casos los ensayos se realizan en probetas unidireccionales, con la fibra orientada en la longitud de la probeta. Sin embargo, las grietas interlaminares que se producen en servicio ocurren en circunstancias muy distintas a las de los ensayos: distintos espesores de los brazos de grieta, piezas con geometría de placa y distintas orientaciones de la interlámina.
En el presente estudio se aborda la descomposición de modos de probetas con geometría de tira, en la que los brazos de grieta pueden ser de distintos espesores y también de material diferente. Se limita al caso en que el material es ortótropo.
Se analiza un sistema en equilibrio sometido a un sistema de fuerzas independientes , siendo los desplazamientos de los puntos de aplicación de las fuerzas en la misma dirección. El trabajo correspondiente a una variación elemental de los desplazamientos viene dado por:
|
(1) |
Habiéndose utilizado el convenio de índices mudo para suprimir el sumatorio de la ecuación (1).
En un avance infinitesimal de grieta, el balance de energía viene dado por:
|
(2) |
Donde representa la energía interna, la energía liberada por unidad de área y el área de la nueva superficie de grieta creada. Si el ancho b es uniforme donde es la longitud del avance de grieta.
Según las ecuaciones (1) y (2), la energía interna viene dada por:
|
(3) |
Aun cuando el avance de grieta es un proceso irreversible, se asume que la energía interna es una función de estado, siendo las variables de estado los desplazamientos independientes y la longitud de grieta . Por lo tanto, también puede expresarse como:
|
(4) |
Identificando términos en las ecuaciones (3) y (4) resulta:
|
(5) |
Con lo que se obtienen el primer teorema de Castigliano y la tasa de liberación de energía como liberación de la energía de deformación.
Por otra parte, se define la energía de deformaicón complementaria :
|
(6) |
Diferenciando la ecuación (6) y sustituyendo en el balance de energía de la ecuación (3) se obtiene:
|
(7) |
Asumiendo que es función de las variables de estado y :
|
(8) |
Identificando términos en las ecuaciones (7) y (8) se tiene:
|
(9) |
Se obtienen así el primer teorema de Engesser y la tasa de liberación de energía como el incremento de energía de deformación complementaria. En sistemas elásticos lineales el teorema de Engesser se convierte en el teorema de Castigliano, por lo que se denomina a veces teorema de Engesser-Castigliano. Para utilizar este teorema se debe expresar en función de las fuerzas aplicadas.
En el caso particular de un sistema elástico lineal sometido a una única carga con un desplazamiento asociado , siendo se tiene:
|
(10) |
Donde representa la flexibilidad del sistema. Utilizando las ecuaciones (5) y (9) resulta en ambos casos que:
|
(11) |
Se ha suprimido el subíndice o indicativo de la derivación parcial, ya que la flexibilidad depende úncamente de la longitud de grieta.
Rice definió la integral J como una integral independiente del camino en una curva cerrada que incluía la punta de grieta:
|
(12) |
Donde es el camino cerrado de integración en sentido anthihorario, la densidad de energía de deformación, el vector tensión y el vector desplazamiento. Se demuestra además que en el caso de comportamiento elástico.
Tomando en las cercanías de la punta de grieta es nulo por lo que la ecuación (12) queda:
|
(13) |
Siendo la tensión normal de apertura, el desplazamiento relativo asociado, la tensión cortante y el desplazamiento relativo asociado. Por lo tanto, la primera integral de la ecuación representa la contribución del modo I y la segunda integral incluye la contribución de los modos II y III, aunque este último no se considera en el presente trabajo.
Las tensiones y desplazamientos citados corresponden a una zona próxima a la punta de grieta, denominada zona cohesiva. Esta zona está relacionada con diversos factores: plastificación debido a que las tensiones alcanzan el límite elástico, puenteo de fibras en la zona de grieta o aparición de microgrietas. La relación entre las tensiones y los desplazamientos mencionados anteriormente en dicha zona se denomina ley cohesiva y se puede obtener de la siguiente manera: asumiendo que es función de los desplazamientos independientes y , puede escribirse:
|
(14) |
Asumiendo ahora que es diferencial exacta:
|
(15) |
Identificando términos en las ecuaciones (14) y (15) resulta:
|
(16) |
En caso de modos puros, únicamente existe un desplazamiento relativo como variable de estado.
Los desplazamientos relativos y son la suma de los desplazamiento correspondiente a puntos homólogos de cada brazo de grieta, entendiendo como puntos homólogos aquellos que están sometidos a la misma tensión. Por lo tanto, en el caso de que los desplazamientos de los puntos homólogos sean distintos, se podrán descomponer en una parte simétrica que realiza trabajo y una parte antisimétrica que no realiza trabajo y está asociada a un movimiento de sólido rígido.
En esta sección se utiliza un procedimiento basado en el estudio de Williams [], determinando directamente la energía complementaria antes y después del avance de grieta sin modificar la carga y utilizando la ecuación (9) para determinar .
Se considera un elemento diferencial de longitud antes y después del avance de grieta. El brazo superior se denomina 1 y el brazo inferior 2. Ambos brazos pueden ser de material y espersor diferente. Considerando únicamente el efecto de la flexión, la energía complementaria antes del avance de grieta viene dada por:
|
(17) |
Donde es la rigidez a flexión de toda la sección. Tras el avance de grieta, la energía complementaria es:
|
(18) |
Donde y son las rigideces a flexión de las partes 1 y 2, respectivamente. De acuerdo a la ecuación (9) la tasa crítica de liberación de energía es:
|
(19) |
Donde y son las relaciones entre rigideces a flexión. Aplicando el principio de superposición los momentos flectores de cada brazo de grieta se descomponen en los correspondientes a los modos de fractura I y II:
|
(20) |
Sustituyendo la ecuación (20) en la ecuación (19) se puede escribir en función de los momentos en modo I y II y por lo tanto se obtienen las tasas de liberación de energía correspondientes a cada modo:
|
(21) |
Por otra parte, los términos que contienen productos de momentos en modos I y II deben ser nulos dado que los modos son independientes, es decir:
|
(22) |
Se supone ahora que los momentos de cada modo en ambos brazos de grieta están relacionados según dos parámetros y , es decir:
|
(23) |
Sustituyendo la ecuación (23) en la ecuación (22):
|
(24) |
Por otra parte, combinando las ecuaciones (20) y (23):
|
(25) |
La determinación de y se trata en la próxima sección. Conocidos estos parámetros, de las ecuaciones (25) se pueden obtener los momentos en los distintos modos:
|
(26) |
Combinando las ecuaciones (26) con las ecuaciones (21) y (23) pueden determinarse y .
El análisis de Williams está particularizado al caso de único material y distintos espesores. Este investigador estableció que para modo puro I se debía verificar que los momentos fueran opuestos, es decir . Por otra parte, para modo puro II estableció la condición de que las curvaturas fueran iguales, con lo que en el caso actual . Utilizando estos valores en el caso de igual espesor, la ecuación (24) se satisface idénticamente.
No obstante en el ensayo ADCB, en el que el material es único pero los espesores distintos, los momentos son iguales pero diferentes estudios muestran que existe combinación de modos. Por lo tanto, los valores utilizados en la descomposición de Williams son válidos únicamente en el caso de igual espesor y un único material.
Teniendo en cuenta la ecuación (13) en la que G se obtiene en función de las tensiones y desplazamientos relativos de la zona cohesiva, se puede afirmar que:
En consecuencia, se pueden probar dos criterios distintos para determinar los parámetros y :
Si la hipótesis es correcta, los valores de y obtenidos en ambos casos deben ser iguales.
No obstante, los desplazamientos normales y tangenciales deben ser calculados teniendo en cuenta todo el sistema, por lo que no pueden relacionarse directamente con los momentos en la punta de grieta. Ello significa que el análisis debe realizarse para cada configuración de ensayo.
Como primera aproximación, puede suponerse que las deformaciones normales correspondientes a modo II son iguales y opuestas en puntos homólogos. Posteriormente, se utiliza la ecuación (24) como segunda condición.
En la figura se muestra un ensayo ADCB de una probeta de dimensiones aproximadas a las utilizadas por Mollón et al [2]. Las propiedades del material corresponden al AS4/8552 de Hexcel Composites. La probeta tiene una longitud total de 150 mm, la grieta es de 50 mm, y el espesor de las partes superior e inferior es 1 y 5 mm, respectivamente. Se impone un desplazamiento de 5 mm a la parte superior. Puede apreciarse la inclinación de la probeta debida a la asimetría de los desplazamientos de apertura.
Mollón et al. proponen una fórmula para obtener la relación de modos para distintas relaciones de espesores. Dicha fórmula está ajustada a partir de resultados numéricos obtenidos con el Método de Dos Pasos desarrollado por los autores. Los resultados obtenidos a partir de la condición de igualdad de deformaciones mencionada anteriormente concuerdan razonablemente con los obtenidos con la fórmula propuesta.
La descomposición de modos en probetas de fractura interlaminar con geometría de tira asimétricas está relacionada con los desplazamientos relativos en las cercanías de la punta de grieta. El método de descomposición de Williams sólo es válido en el caso de probetas de un mismo material con la grieta centrada. Las relaciones de modos obtenidas en el caso de ADCB mediante una primera aproximación basada en deformaciones normales concuerdan razonablemente con resultados obtenidos mediante el MEF.
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] J.G. Williams, International Journal of Fracture 36, pág. 101-119 (1988).
[2] V.Mollón, J. Bonhomme, J. Viña, A. Argüelles, Polymer Testing 29, pág. 766-770 (2010). doi:10.1016/j.polymertesting.2010.04.001
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.020
Licence: Other
Are you one of the authors of this document?