m (Marherna moved page Review 246308283647 to Mujika et al 2023a) |
|||
(9 intermediate revisions by 2 users not shown) | |||
Line 325: | Line 325: | ||
{| style="text-align: center; margin:auto;" | {| style="text-align: center; margin:auto;" | ||
|- | |- | ||
− | | | + | | <span id='Eq21'></span> <math>C'=a{\int }_{L_x}NN'dx+d{\int }_{L_x}MM'dx+</math><math>s{\int }_{L_x}VV{}'dx+b\left({\int }_{L_x}N'Mdx+\right. </math><math>\left. {\int }_{L_x}NM'dx\right)+</math> |
+ | <br><math>+(aN_{HT}+bM_{HT}){\int }_{L_x}N'dx+(bN_{HT}+dM_{HT}){\int }_{L_x}M'dx</math> | ||
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (21) | | style="width: 5px;text-align: right;white-space: nowrap;" | (21) | ||
|} | |} | ||
− | |||
===3.4. Tasa critica de liberación de energía=== | ===3.4. Tasa critica de liberación de energía=== | ||
Line 367: | Line 367: | ||
''W'' es la energía de deformación por unidad de volumen. Γ es una curva que rodea la punta de grieta, evaluándose la integral en sentido anti horario comenzando por la zona inferior de la zona cohesiva y continuando por la zona superior, siguiendo la curva Γ. <math display="inline">\overset{\rightarrow}{T}</math> es el vector tensión definido según la normal externa a la curva Γ, <math display="inline">\overset{\rightarrow}{\delta }</math> es el vector desplazamiento y ''ds'' es un elemento de longitud de arco. De acuerdo a la Figura 1, ''dy'' = 0. | ''W'' es la energía de deformación por unidad de volumen. Γ es una curva que rodea la punta de grieta, evaluándose la integral en sentido anti horario comenzando por la zona inferior de la zona cohesiva y continuando por la zona superior, siguiendo la curva Γ. <math display="inline">\overset{\rightarrow}{T}</math> es el vector tensión definido según la normal externa a la curva Γ, <math display="inline">\overset{\rightarrow}{\delta }</math> es el vector desplazamiento y ''ds'' es un elemento de longitud de arco. De acuerdo a la Figura 1, ''dy'' = 0. | ||
− | <div class="center" style="width: auto; margin-left: | + | <div class="center" style="width: auto; margin-left: 200px; margin-right: auto;"> |
[[File:Review_246308283647_9191_mat_com_mujika1.png|400px]]</div> | [[File:Review_246308283647_9191_mat_com_mujika1.png|400px]]</div> | ||
Line 549: | Line 549: | ||
2. Determinar las fuerzas y desplazamientos en la punta de grieta para determinar las tasas de liberación de energía de cada modo. | 2. Determinar las fuerzas y desplazamientos en la punta de grieta para determinar las tasas de liberación de energía de cada modo. | ||
− | <div id="_GoBack" class="center" style="width: auto; margin-left: | + | <div id="_GoBack" class="center" style="width: auto; margin-left: 20px; margin-right: auto;"> |
− | [[ | + | [[File:Review_246308283647_5389_mat_com_mujika2.png|372px]] </div> |
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
Line 559: | Line 559: | ||
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;"> | ||
− | Figura | + | Figura 3. Ensayo de Flexión con Entalla Final Asimétrico (AENF).</div> |
==Bibliografía== | ==Bibliografía== |
ABSTRACT: The fundamentals of new analytical models based on energy and coenergy of deformation are explained to analytically address asymmetric interlaminar fracture tests. After reviewing some concepts of laminated beam theory, the energy release rate as a function of coenergy is determined from a global point of view, including in an unprecedented way the contribution of hygrothermal stresses. On the other hand, after the determination of the J integral in the cohesive zone, the conditions required for the asymmetric fracture tests to occur in pure mode in linear elastic fracture are explained, based on the forces and displacements of the crack tip. Furthermore, knowing these forces and displacements, it is possible to determine the energy release rate that corresponds to each fracture mode.
Keywords: Fractura interlaminar, Modo I, Modo II, Combinación de modos I-II
La fractura interlaminar en materiales compuestos se analiza habitualmente mediante los ensayos de Viga en Doble Voladizo (Double Cantilever Beam, DCB) y de Flexión con Entalla Final (End Notched Flexure, ENF) en modos I y II, respectivamente. Estas mismas configuraciones se utilizan también para caracterizar uniones adhesivas, con brazos de grieta simétricos para que la fractura se produzca en modo puro. En caso de que exista asimetría geométrica o material, se puede producir combinación de modos I/II. En este trabajo, se abordan las bases de nuevos planteamientos analíticos para los ensayos asimétricos ADCB y AENF.
Además de la flexibilidad de la probeta y de la tasa crítica de liberación de energía total, se determinan las tasas de energía correspondientes a modo I y modo II, con un planteamiento local en la punta de grieta. Dicho planteamiento es equivalente a la Técnica de Cierre Virtual de Grieta (Virtual Crack Closure Technique, (VCCT), utilizado cuando el problema se aborda numéricamente con el Método de los Elementos Finitos.
Se supone un cuerpo sometido a cargas generalizadas concentradas Fi. El desplazamiento generalizado del punto de aplicación de dicha fuerza en su dirección es δi. En un pequeño incremento de grieta, el trabajo realizado en un desplazamiento infinitesimal dδi es:
|
(1) |
donde dW es el trabajo realizado por las fuerzas aplicadas; dU es el cambio de energía de deformación; G es la energía necesaria para el avance de grieta por unidad de área; w es el ancho de la grieta; y da es el avance diferencial de grieta. El trabajo diferencial realizado por las fuerzas externas Fi y sus respectivos desplazamientos δi, utilizando la convención de índices repetidos, es . Así la Ec. (1) se puede escribir como:
|
(2) |
Asumiendo que la energía de deformación es una función de estado que depende de los desplazamientos δi y de la longitud de grieta a:
|
(3) |
Identificando términos en las Ecs. (2) y (3) resulta:
|
(4) |
|
(5) |
La Ec. (4) es el primer teorema de Castigliano y la Ec. (5) es la Tasa de Liberación de Energía de Deformación. La energía complementaria o coenergía se define como:
|
(6) |
El nombre de coenergía se utiliza habitualmente en el caso de fuerzas magnéticas [1], pero también se ha utilizado en Mecánica [2]. En el presente estudio, se adopta para la energía de deformación complementaria en el caso mecánico, utilizando la letra C para denominarla. Desde un punto de vista termodinámico, corresponde a la energía libre de Gibbs. Derivando la Ec. (6) y sustituyendo la Ec. (2) resulta:
|
(7) |
Asumiendo que la coenergía es una función de estado de las fuerzas generalizadas Fi y la longitud de grieta a:
|
(8) |
Identificando términos en las Ecs. (7) y (8) resulta:
|
(9) |
|
(10) |
La Ec. (9) es el teorema de Engesser-Castigliano y la Ec. (10) es la tasa de liberación de energía de deformación. En el presente estudio, se utiliza la coenergía para considerar las fuerzas generalizadas como variables de estado.
Suponiendo un laminado de N láminas, la tensión normal en la lámina k es:
|
(11) |
Donde son las tensiones higrotérmicas de la lámina k. Dado que la fuerza normal N y el momento flector M son la resultante y el momento resultante de las fuerzas normales, resulta:
|
(12) |
Donde el subíndice G se refiere a Global y HT hace referencia a la fuerza y momento higrotérmico. Los coeficientes de rigidez son:
|
(13) |
Inviertiendo la relación se obtiene:
|
(14) |
Las tensiones normales vienen dadas por:
|
(15) |
Suponiendo comportamiento elástico lineal, y teniendo en cuenta las deformaciones normales higrotérmicas, la coenergía por unidad de volumen viene dada por:
|
(16) |
Integrando en el área de la sección, la coenergía por unidad de longitud C1 es:
|
(17) |
La coenergía debida a las tensiones cortantes es:
|
(18) |
donde s es la flexibilidad equivalente de cortadura. En el caso de un material de una única capa ortótropo es . La coenergía de toda la viga que corresponde a tensiones normales y cortantes es:
|
(19) |
Como se ha visto en la Ed. (9), según el teorema de Engesser-Castigliano, el desplazamiento generalizado δk del punto de aplicación de la fuerza generalizada Fk es:
|
(20) |
En la Ec. (20), las primas indican derivadas respecto a Fk. Dado que las fuerzas y momentos higrotérmicos no dependen de las fuerzas aplicadas Fk, la Ec. (20) puede expresarse como:
|
(21) |
Según la Ec. (10), la tasa crítica de liberación de energía es:
|
(22) |
En el caso de la Ec. (22) es necesario tener en cuenta que los límites de integración pueden depender de la longitud de grieta, por lo que es necesario determinar la coenergía antes de derivarla respecto de a. Otra alternativa, es utilizar la fórmula de Leibniz relativa a diferenciación de integrales respecto a un parámetro.
En este apartado se realiza una explicación de la integral J definida por Rice [3]. Se aplica a la zona cohesiva, ya que será la herramienta básica para la descomposición de modos. La zona cohesiva corresponde a la zona de dañada cercana a la punta de grieta. En fractura elástica lineal la integral J es equivalente a la tasa de liberación de energía G. Dicha integral se define como:
|
(23) |
W es la energía de deformación por unidad de volumen. Γ es una curva que rodea la punta de grieta, evaluándose la integral en sentido anti horario comenzando por la zona inferior de la zona cohesiva y continuando por la zona superior, siguiendo la curva Γ. es el vector tensión definido según la normal externa a la curva Γ, es el vector desplazamiento y ds es un elemento de longitud de arco. De acuerdo a la Figura 1, dy = 0.
Según la Figura 1:
|
(24) |
La curva Γ se divide Γ en ΓI y ΓII, que corresponden a las caras inferior y superior de la zona cohesiva, respectivamente, como se muestra en la Figura 1. Siendo las tensiones iguales en ambas curvas, teniendo en cuenta la ecuación 3 resulta:
|
(25) |
Sustituyendo los resultados de la Ec. (25) la integral J es:
|
(26) |
Los desplazamientos relativos normal y tangencial en la Ec.(26) se definen como:
|
(27) |
Dado que los desplazamientos relativos son función de x, resulta:
|
(28) |
Además, los límites de la integral de la Ec.(26) se modifican según:
|
(29) |
siendo Δn y Δt los valores de los desplazamientos relativos en la punta de grieta. Además, siendo y , combinando las Ecs. (26), (28) y (29) resulta:
|
(30) |
La deducción de la Ec.(30) no depende de la longitud de la zona cohesiva. El cambio de variable permite además interpretar la integral J de dos formas: según la Ec. (26), es una integral extendida a lo largo de la zona cohesiva. Según la Ec. (30) puede considerarse como el trabajo por unidad de área realizado por las tensiones normal y cortante de la punta de grieta, cuando el desplazamiento relativo varía desde cero hasta su valor final.
La primera integral de la Ec. (30) es la contribución de modo I y la segunda integral es la contribución que corresponde a los modos II y III. Considerando sólo modos I y II, si la segunda integral es nula, se tiene modo puro I, siendo y si la primera integral es nula, se tiene modo puro II, siendo .
Se realiza el desarrollo para modo I. Suponiendo una ley cohesiva que corresponde a una rotura frágil, se supone una ley cohesiva lineal, en la cual cuando el desplazamiento relativo del punto de grieta alcanza su valor final, el vector tensión tiende a 0 y la grieta avanza da. Siendo Fn la fuerza normal en la punta de grieta, resulta.
|
(31) |
Análogamente, en modo II
|
(32) |
Según las Ecs. (31) y (32), conociendo las fuerzas y desplazamientos de la punta de grieta, se pueden obtener las tasas de liberación de energía que corresponden a cada modo. Además, proporcionan criterios para determinar las condiciones de modo puro:
La fuerza nula en la punta de grieta indica que no es necesario impedir ningún desplazamiento de la misma dirección cercano a la punta de grieta. De ahí la simultaneidad de condiciones de fuerzas y derivadas de desplazamiento, ya que dichas derivadas indican desplazamientos cercanos a la punta de grieta.
Las ecuaciones desarrolladas en los apartados anteriores pueden aplicarse para realizar modelos matemáticos del comportamiento mecánico de probetas de fractura interlaminar. En particular, se pueden analizar los ensayos de Viga en Doble Voladizo Asimétrica (Asymmetric Double Cantielver Beam, ADCB) y de Flexión con Entalla Final Asimétrica (Asymmetric End Notched Flexure, ENF). Dado que en ambos ensayos puede existir combinación de modos, como se muestra en las Figuras 2 y 3, el problema general puede descomponerse en otros 2:
1. Hallar las condiciones que se deben cumplir para que los ensayos se realicen en modo puro.
2. Determinar las fuerzas y desplazamientos en la punta de grieta para determinar las tasas de liberación de energía de cada modo.
[1] | O.B. Mawardi (1957). On the concept of coenergy. Journal of the Franklin Institute, 264(4), 313–332. DOI: 10.1016/0016-0032(57)90190-4 |
[2] | Seeler, K. A. (2014). System Dynamics. An Introduction for Mechanical Engineers, Springer. DOI: 10.1007/978-1-4614-9152-1 |
[3] | J.R. Rice (1968). A path independent integral and the approximate analysis of strain concentration by notches and cracks, Journal of Applied Mechanics, 35(2), 379–386. DOI: 10.1115/1.3601206 |
Published on 27/01/24
Accepted on 18/12/23
Submitted on 19/05/23
Volume 08 - COMUNICACIONES MATCOMP21 (2022) Y MATCOMP23 (2023), Issue Núm. 3 - Caracterización - Daño, Fractura y Fatiga, 2024
Licence: Other
Are you one of the authors of this document?