La mecánica de fracturas es el campo de la mecánica que se ocupa del estudio de la propagación de grietas en materiales. Utiliza métodos de mecánica analítica de sólidos para calcular la fuerza impulsora sobre una grieta y los de mecánica experimental de sólidos para caracterizar la resistencia del material a la fractura .
En teoría, la tensión delante de la punta de una grieta afilada se vuelve infinita y no puede usarse para describir el estado alrededor de una grieta. La mecánica de fractura se utiliza para caracterizar las cargas en una grieta, generalmente usando un solo parámetro para describir el estado de carga completo en la punta de la grieta. Se han desarrollado varios parámetros diferentes. Cuando la zona plástica en la punta de la grieta es pequeña en relación con la longitud de la grieta, el estado de tensión en la punta de la grieta es el resultado de fuerzas elásticas dentro del material y se denomina mecánica de fractura elástica lineal ( LEFM ) y se puede caracterizar utilizando la tensión. factor de intensidad . Aunque la carga sobre una grieta puede ser arbitraria, en 1957 G. Irwin descubrió que cualquier estado podía reducirse a una combinación de tres factores de intensidad de tensión independientes:
Cuando el tamaño de la zona plástica en la punta de la grieta es demasiado grande, se puede utilizar la mecánica de fractura elástico-plástica con parámetros como la integral J o el desplazamiento de apertura de la punta de la grieta .
El parámetro caracterizador describe el estado de la punta de la grieta que luego puede relacionarse con las condiciones experimentales para garantizar la similitud . El crecimiento de grietas ocurre cuando los parámetros generalmente exceden ciertos valores críticos. La corrosión puede causar que una grieta crezca lentamente cuando se excede el umbral de intensidad de la tensión de corrosión bajo tensión . De manera similar, pequeños defectos pueden provocar el crecimiento de grietas cuando se someten a cargas cíclicas. Conocida como fatiga , se descubrió que para grietas largas, la tasa de crecimiento se rige en gran medida por el rango de intensidad de la tensión experimentada por la grieta debido a la carga aplicada. Se producirá una fractura rápida cuando la intensidad de la tensión exceda la tenacidad a la fractura del material. La predicción del crecimiento de grietas es el núcleo de la disciplina de diseño mecánico con tolerancia al daño .
Los procesos de fabricación, procesamiento, mecanizado y conformado de materiales pueden introducir fallas en un componente mecánico terminado. En todas las estructuras metálicas se encuentran defectos interiores y superficiales derivados del proceso de fabricación. No todos estos defectos son inestables en condiciones de servicio. La mecánica de fracturas es el análisis de fallas para descubrir aquellas que son seguras (es decir, que no crecen) y aquellas que pueden propagarse como grietas y causar fallas en la estructura defectuosa. A pesar de estos defectos inherentes, es posible lograr mediante el análisis de tolerancia al daño el funcionamiento seguro de una estructura. La mecánica de fracturas como tema de estudio crítico apenas existe desde hace un siglo y, por lo tanto, es relativamente nueva. [1] [2]
La mecánica de fracturas debe intentar proporcionar respuestas cuantitativas a las siguientes preguntas: [2]
La mecánica de fracturas fue desarrollada durante la Primera Guerra Mundial por el ingeniero aeronáutico inglés AA Griffith (de ahí el término grieta de Griffith ) para explicar la falla de materiales frágiles. [5] El trabajo de Griffith estuvo motivado por dos hechos contradictorios:
Se necesitaba una teoría para conciliar estas observaciones contradictorias. Además, los experimentos con fibras de vidrio que realizó el propio Griffith sugirieron que la tensión de fractura aumenta a medida que disminuye el diámetro de la fibra. Por lo tanto, la resistencia a la tracción uniaxial, que se había utilizado ampliamente para predecir la falla del material antes de Griffith, no podía ser una propiedad del material independiente de la muestra. Griffith sugirió que la baja resistencia a la fractura observada en los experimentos, así como la dependencia del tamaño de la resistencia, se debía a la presencia de defectos microscópicos en el material a granel.
Para verificar la hipótesis del defecto, Griffith introdujo un defecto artificial en sus muestras de vidrio experimentales. El defecto artificial tenía la forma de una grieta superficial mucho más grande que otros defectos de una muestra. Los experimentos demostraron que el producto de la raíz cuadrada de la longitud del defecto ( ) y la tensión de fractura ( ) era casi constante, lo que se expresa mediante la ecuación:
Una explicación de esta relación en términos de la teoría de la elasticidad lineal es problemática. La teoría de la elasticidad lineal predice que la tensión (y por tanto la deformación) en la punta de un defecto agudo en un material elástico lineal es infinita. Para evitar ese problema, Griffith desarrolló un enfoque termodinámico para explicar la relación que observó.
El crecimiento de una grieta, la extensión de las superficies a ambos lados de la grieta, requiere un aumento en la energía superficial . Griffith encontró una expresión para la constante en términos de la energía superficial de la grieta resolviendo el problema de elasticidad de una grieta finita en una placa elástica. Brevemente, el enfoque fue:
donde es el módulo de Young del material y es la densidad de energía superficial del material. Suponiendo y dando una excelente concordancia entre la tensión de fractura predicha por Griffith y los resultados experimentales para el vidrio.
Para el caso simple de una placa rectangular delgada con una grieta perpendicular a la carga, la tasa de liberación de energía, , se convierte en:
donde es la tensión aplicada, es la mitad de la longitud de la grieta y es el módulo de Young , que para el caso de deformación plana debe dividirse por el factor de rigidez de la placa . La tasa de liberación de energía de deformación puede entenderse físicamente como: la tasa a la que la energía es absorbida por el crecimiento de la grieta .
Sin embargo, también tenemos eso:
Si ≥ , este es el criterio por el cual la grieta comenzará a propagarse.
Para materiales altamente deformados antes de la propagación de la grieta, la formulación de la mecánica de fractura elástica lineal ya no es aplicable y es necesario un modelo adaptado para describir el campo de tensión y desplazamiento cerca de la punta de la grieta, como en la fractura de materiales blandos .
El trabajo de Griffith fue en gran medida ignorado por la comunidad de ingenieros hasta principios de la década de 1950. Las razones de esto parecen ser (a) en los materiales estructurales reales, el nivel de energía necesario para causar la fractura es órdenes de magnitud mayor que la energía superficial correspondiente, y (b) en los materiales estructurales siempre hay algunas deformaciones inelásticas alrededor de la grieta. frente que haría que la suposición de un medio elástico lineal con tensiones infinitas en la punta de la grieta sea muy poco realista. [6]
La teoría de Griffith coincide excelentemente con los datos experimentales de materiales frágiles como el vidrio. Para materiales dúctiles como el acero , aunque la relación aún se mantiene, la energía superficial ( γ ) predicha por la teoría de Griffith suele ser irrealmente alta. Un grupo que trabajó bajo la dirección de GR Irwin [7] en el Laboratorio de Investigación Naval (NRL) de EE. UU. durante la Segunda Guerra Mundial se dio cuenta de que la plasticidad debe desempeñar un papel importante en la fractura de materiales dúctiles.
En materiales dúctiles (e incluso en materiales que parecen frágiles [8] ), se desarrolla una zona plástica en la punta de la grieta. A medida que aumenta la carga aplicada , la zona plástica aumenta de tamaño hasta que la grieta crece y el material elásticamente deformado detrás de la punta de la grieta se descarga. El ciclo de carga y descarga de plástico cerca de la punta de la grieta conduce a la disipación de energía en forma de calor . Por tanto, hay que añadir un término disipativo a la relación de balance de energía ideada por Griffith para materiales frágiles. En términos físicos, se necesita energía adicional para el crecimiento de grietas en materiales dúctiles en comparación con materiales frágiles.
La estrategia de Irwin fue dividir la energía en dos partes:
Entonces la energía total es:
donde es la energía superficial y es la disipación plástica (y la disipación de otras fuentes) por unidad de área de crecimiento de grieta.
La versión modificada del criterio energético de Griffith se puede escribir como
Para materiales frágiles como el vidrio, domina el término de energía superficial y . Para materiales dúctiles como el acero, domina el término de disipación plástica y . Para polímeros cercanos a la temperatura de transición vítrea , tenemos valores intermedios de entre 2 y 1000 .
Otro logro significativo de Irwin y sus colegas fue encontrar un método para calcular la cantidad de energía disponible para la fractura en términos de tensiones asintóticas y campos de desplazamiento alrededor del frente de una grieta en un sólido elástico lineal. [7] Esta expresión asintótica para el campo de tensión en la carga del modo I está relacionada con el factor de intensidad de tensión siguiente: [9]
donde están las tensiones de Cauchy , es la distancia desde la punta de la grieta, es el ángulo con respecto al plano de la grieta y son funciones que dependen de la geometría de la grieta y las condiciones de carga. Irwin llamó a esta cantidad factor de intensidad del estrés. Dado que la cantidad no tiene dimensiones, el factor de intensidad de tensión se puede expresar en unidades de .
La intensidad de la tensión reemplazó la tasa de liberación de energía de deformación y un término llamado tenacidad a la fractura reemplazó a la energía de debilidad de la superficie. Ambos términos están simplemente relacionados con los términos energéticos que utilizó Griffith:
y
donde es la intensidad del estrés modal , la tenacidad a la fractura y es el índice de Poisson.
La fractura ocurre cuando . Para el caso especial de deformación por deformación plana, se convierte y se considera una propiedad del material. El subíndice surge debido a las diferentes formas de cargar un material para permitir que se propague una grieta . Se refiere a la carga denominada "modo " en contraposición al modo o :
La expresión será diferente para geometrías distintas a la placa infinita con fisura central, como se analiza en el artículo sobre el factor de intensidad de tensión. En consecuencia, es necesario introducir un factor de corrección adimensional , , para caracterizar la geometría. Este factor de corrección, también denominado a menudo factor de forma geométrica , viene dado por series determinadas empíricamente y tiene en cuenta el tipo y la geometría de la grieta o muesca. Así tenemos:
donde es una función de la longitud y el ancho de la grieta de la lámina dada, para una lámina de ancho finito que contiene una grieta en todo el espesor de longitud , por:
Irwin fue el primero en observar que si el tamaño de la zona plástica alrededor de una grieta es pequeño en comparación con el tamaño de la grieta, la energía requerida para hacer crecer la grieta no dependerá críticamente del estado de tensión (la zona plástica) en la punta de la grieta. [6] En otras palabras, se puede utilizar una solución puramente elástica para calcular la cantidad de energía disponible para la fractura.
La tasa de liberación de energía para el crecimiento de la grieta o la tasa de liberación de energía de deformación se puede calcular como el cambio en la energía de deformación elástica por unidad de área de crecimiento de la grieta, es decir,
donde U es la energía elástica del sistema y a es la longitud de la grieta. La carga P o el desplazamiento u son constantes al evaluar las expresiones anteriores.
Irwin demostró que para un modo I crack (modo de apertura), la tasa de liberación de energía de deformación y el factor de intensidad de tensión están relacionados por:
donde E es el módulo de Young , ν es la relación de Poisson y K I es el factor de intensidad de tensión en el modo I. Irwin también demostró que la tasa de liberación de energía de deformación de una grieta plana en un cuerpo elástico lineal se puede expresar en términos del modo Factores de intensidad de tensión I, modo II (modo deslizante) y modo III (modo desgarro) para las condiciones de carga más generales.
A continuación, Irwin adoptó el supuesto adicional de que el tamaño y la forma de la zona de disipación de energía permanecen aproximadamente constantes durante la fractura frágil. Esta suposición sugiere que la energía necesaria para crear una superficie de fractura unitaria es una constante que depende únicamente del material. Esta nueva propiedad del material recibió el nombre de tenacidad a la fractura y se denominó G Ic . Hoy en día, es el factor crítico de intensidad de tensión K Ic , encontrado en la condición de deformación plana, el que se acepta como la propiedad definitoria en la mecánica de fractura elástica lineal.
En teoría, la tensión en la punta de la grieta, donde el radio es casi cero, tendería a infinito. Esto se consideraría una singularidad de tensión, lo cual no es posible en aplicaciones del mundo real. Por esta razón, en estudios numéricos en el campo de la mecánica de fracturas, a menudo es apropiado representar las grietas como muescas de punta redonda , con una región de concentración de tensiones dependiente de la geometría que reemplaza la singularidad de la punta de la grieta. [9] En realidad, se ha descubierto que la concentración de tensión en la punta de una grieta dentro de materiales reales tiene un valor finito pero mayor que la tensión nominal aplicada a la muestra.
Sin embargo, debe existir algún tipo de mecanismo o propiedad del material que impida que dicha grieta se propague espontáneamente. La suposición es que la deformación plástica en la punta de la grieta efectivamente la embota. Esta deformación depende principalmente de la tensión aplicada en la dirección aplicable (en la mayoría de los casos, esta es la dirección y de un sistema de coordenadas cartesiano regular), la longitud de la grieta y la geometría de la muestra. [10] Para estimar cómo esta zona de deformación plástica se extendía desde la punta de la grieta, Irwin equiparó el límite elástico del material con las tensiones de campo lejano de la dirección y a lo largo de la grieta (dirección x) y resolvió el radio efectivo. A partir de esta relación, y suponiendo que la grieta está cargada hasta el factor de intensidad de tensión crítica, Irwin desarrolló la siguiente expresión para el radio idealizado de la zona de deformación plástica en la punta de la grieta:
Los modelos de materiales ideales han demostrado que esta zona de plasticidad está centrada en la punta de la grieta. [11] Esta ecuación proporciona el radio ideal aproximado de la deformación de la zona plástica más allá de la punta de la grieta, lo cual es útil para muchos científicos estructurales porque brinda una buena estimación de cómo se comporta el material cuando se somete a tensión. En la ecuación anterior, los parámetros del factor de intensidad de tensión e indicador de tenacidad del material, y el límite elástico , son importantes porque ilustran muchas cosas sobre el material y sus propiedades, así como sobre el tamaño de la zona plástica. Por ejemplo, si es alto, entonces se puede deducir que el material es tenaz, y si es bajo, se sabe que el material es más dúctil. La relación de estos dos parámetros es importante para el radio de la zona plástica. Por ejemplo, si es pequeño, entonces la proporción al cuadrado de a es grande, lo que da como resultado un radio plástico mayor. Esto implica que el material puede deformarse plásticamente y, por tanto, es resistente. [10] Esta estimación del tamaño de la zona plástica más allá de la punta de la grieta se puede utilizar para analizar con mayor precisión cómo se comportará un material en presencia de una grieta.
El mismo proceso descrito anteriormente para la carga de un solo evento también se aplica a la carga cíclica. Si hay una grieta en una muestra que sufre cargas cíclicas, la muestra se deformará plásticamente en la punta de la grieta y retrasará el crecimiento de la grieta. En caso de sobrecarga o excursión, este modelo cambia ligeramente para adaptarse al aumento repentino de la tensión que experimentó el material anteriormente. Con una carga suficientemente alta (sobrecarga), la grieta crece fuera de la zona plástica que la contenía y deja atrás la bolsa de la deformación plástica original. Ahora, suponiendo que la tensión de sobrecarga no es lo suficientemente alta como para fracturar completamente la muestra, la grieta sufrirá una mayor deformación plástica alrededor de la nueva punta de la grieta, ampliando la zona de tensiones plásticas residuales. Este proceso endurece y prolonga aún más la vida útil del material porque la nueva zona plástica es más grande de lo que sería en las condiciones de tensión habituales. Esto permite que el material pase por más ciclos de carga. Esta idea se puede ilustrar mejor con el gráfico del aluminio con una grieta central sometida a eventos de sobrecarga. [12]
Pero a los investigadores del NRL les surgió un problema porque los materiales navales, como por ejemplo el acero para barcos, no son perfectamente elásticos sino que sufren una deformación plástica significativa en la punta de una grieta. Una suposición básica en la mecánica de fractura elástica lineal de Irwin es la fluencia a pequeña escala, la condición de que el tamaño de la zona plástica sea pequeño en comparación con la longitud de la grieta. Sin embargo, esta suposición es bastante restrictiva para ciertos tipos de fallas en aceros estructurales, aunque dichos aceros pueden ser propensos a fracturarse por fragilidad, lo que ha llevado a una serie de fallas catastróficas.
La mecánica de fractura elástica lineal tiene un uso práctico limitado para aceros estructurales y las pruebas de tenacidad a la fractura pueden ser costosas.
La mayoría de los materiales de ingeniería muestran cierto comportamiento elástico e inelástico no lineal en condiciones de operación que involucran grandes cargas. [ cita necesaria ] En tales materiales, los supuestos de la mecánica de fractura elástica lineal pueden no ser válidos, es decir,
Por lo tanto, se necesita una teoría más general del crecimiento de grietas para materiales elástico-plásticos que pueda explicar:
Históricamente, el primer parámetro para la determinación de la tenacidad a la fractura en la región elastoplástica era el desplazamiento de apertura de la punta de la grieta (CTOD) o "abertura en el vértice de la grieta". Este parámetro fue determinado por Wells durante los estudios de aceros estructurales, que debido a su alta tenacidad no podían caracterizarse con el modelo de mecánica de fractura elástica lineal. Observó que, antes de que ocurriera la fractura, las paredes de la grieta se estaban saliendo [ se necesita aclaración ] y que la punta de la grieta, después de la fractura, variaba de aguda a redondeada debido a la deformación plástica. Además, el redondeo de la punta de la grieta fue más pronunciado en aceros con tenacidad superior.
Hay varias definiciones alternativas de CTOD. En las dos definiciones más comunes, CTOD es el desplazamiento en la punta original de la grieta y la intersección de 90 grados. Rice sugirió la última definición y se usa comúnmente para inferir CTOD en modelos de elementos finitos de los mismos. Tenga en cuenta que estas dos definiciones son equivalentes si la punta de la grieta se embota en un semicírculo.
La mayoría de las mediciones de laboratorio de CTOD se han realizado en muestras con bordes agrietados cargados con flexión de tres puntos. Los primeros experimentos utilizaron un calibre plano en forma de paleta que se insertaba en la grieta; Cuando se abrió la grieta, el medidor de paletas giró y se envió una señal electrónica a un trazador xy. Sin embargo, este método era inexacto porque era difícil alcanzar la punta de la grieta con el medidor de paleta. Hoy en día, se mide el desplazamiento V en la boca de la grieta y el CTOD se infiere suponiendo que las mitades de la muestra son rígidas y giran alrededor de un punto de articulación (la punta de la grieta).
Un primer intento en la dirección de la mecánica de fractura elástico-plástica fue la curva de resistencia a la extensión de grietas de Irwin , la curva de resistencia al crecimiento de grietas o la curva R. Esta curva reconoce el hecho de que la resistencia a la fractura aumenta con el tamaño de la grieta en materiales elástico-plásticos. La curva R es una gráfica de la tasa total de disipación de energía en función del tamaño de la grieta y puede usarse para examinar los procesos de crecimiento lento y estable de grieta y fractura inestable. Sin embargo, la curva R no se utilizó ampliamente en aplicaciones hasta principios de la década de 1970. Las razones principales parecen ser que la curva R depende de la geometría de la muestra y la fuerza impulsora de la grieta puede ser difícil de calcular. [6]
A mediados de la década de 1960, James R. Rice (entonces en la Universidad de Brown ) y GP Cherepanov desarrollaron de forma independiente una nueva medida de tenacidad para describir el caso en el que hay suficiente deformación en la punta de la grieta como para que la pieza ya no obedezca a la aproximación lineal-elástica. El análisis de Rice, que supone una deformación elástica no lineal (o teoría de la deformación monótona plástica ) delante de la punta de la grieta, se denomina integral J. [13] Este análisis se limita a situaciones en las que la deformación plástica en la punta de la grieta no se extiende hasta el borde más alejado de la parte cargada. También exige que el comportamiento elástico no lineal supuesto del material sea una aproximación razonable en forma y magnitud a la respuesta de carga del material real. El parámetro de falla elástico-plástica se denomina J Ic y se convierte convencionalmente a K Ic usando la siguiente ecuación. Tenga en cuenta también que el enfoque integral J se reduce a la teoría de Griffith para el comportamiento lineal-elástico.
La definición matemática de J-integral es la siguiente:
dónde
Desde que los ingenieros se acostumbraron a utilizar K Ic para caracterizar la tenacidad a la fractura, se ha utilizado una relación para reducir J Ic a esta:
Cuando una región importante alrededor de la punta de una grieta ha sufrido deformación plástica, se pueden utilizar otros enfoques para determinar la posibilidad de una mayor extensión de la grieta y la dirección de crecimiento y ramificación de la grieta. Una técnica sencilla que se incorpora fácilmente a los cálculos numéricos es el método del modelo de zona cohesiva , que se basa en conceptos propuestos independientemente por Barenblatt [14] y Dugdale [15] a principios de los años 1960. La relación entre los modelos de Dugdale-Barenblatt y la teoría de Griffith fue discutida por primera vez por Willis en 1967. [16] Rice demostró la equivalencia de los dos enfoques en el contexto de la fractura frágil en 1968. [13]
Supongamos que un material tiene un límite elástico y una tenacidad a la fractura en el modo I. Según la mecánica de fractura, el material fallará bajo tensión . Según la plasticidad, el material cederá cuando . Estas curvas se cruzan cuando . Este valor de se denomina tamaño de defecto de transición y depende de las propiedades del material de la estructura. Cuando la falla se rige por la fluencia plástica, y cuando la falla se rige por la mecánica de fractura. El valor de para aleaciones de ingeniería es de 100 mm y para cerámicas es de 0,001 mm. [ cita necesaria ] Si asumimos que los procesos de fabricación pueden dar lugar a fallas del orden de micrómetros , entonces, se puede ver que las cerámicas tienen más probabilidades de fallar por fractura, mientras que las aleaciones de ingeniería fallarían por deformación plástica.
El análisis de fractura del hormigón es parte de la mecánica de fractura que estudia la propagación de grietas y los modos de falla relacionados en el hormigón . [17] Como se utiliza ampliamente en la construcción, el análisis de fracturas y los modos de refuerzo son una parte importante del estudio del hormigón, y los diferentes hormigones se caracterizan en parte por sus propiedades de fractura. [18] Las fracturas comunes incluyen las fracturas en forma de cono que se forman alrededor de los anclajes bajo resistencia a la tracción.
Bažant (1983) propuso un modelo de banda de fisuras para materiales como el hormigón cuya naturaleza homogénea cambia aleatoriamente en un cierto rango. [17] También observó que en concreto simple, el efecto del tamaño tiene una fuerte influencia en el factor de intensidad de tensión crítica , [19] y propuso la relación
= / √(1+{ / }), [19] [20]
donde = factor de intensidad de tensión, = resistencia a la tracción, = tamaño de la muestra, = tamaño máximo del agregado y = una constante empírica.
La Mecánica de Fractura Atomística (AFM) es un campo relativamente nuevo que estudia el comportamiento y las propiedades de los materiales a escala atómica cuando se someten a fractura. Integra conceptos de la mecánica de fracturas con simulaciones atomísticas para comprender cómo las grietas inician, se propagan e interactúan con la microestructura de los materiales. Mediante el uso de técnicas como las simulaciones de dinámica molecular (MD), AFM puede proporcionar información sobre los mecanismos fundamentales de formación y crecimiento de grietas, el papel de los enlaces atómicos y la influencia de los defectos e impurezas de los materiales en el comportamiento de las fracturas. [21]