Figura 1: Gráfico típico de la tasa de crecimiento de grietas versus el rango de intensidad de tensión. La ecuación de París-Erdogan se ajusta a la región lineal central del Régimen B.
Se utiliza una ecuación de crecimiento de grietas para calcular el tamaño de una grieta por fatiga que crece debido a cargas cíclicas. El crecimiento de una grieta por fatiga puede provocar fallas catastróficas, particularmente en el caso de aeronaves. Cuando muchas grietas por fatiga en crecimiento interactúan entre sí, se conoce como daño por fatiga generalizado . Se puede utilizar una ecuación de crecimiento de grietas para garantizar la seguridad, tanto en la fase de diseño como durante la operación, prediciendo el tamaño de las grietas. En estructuras críticas, las cargas se pueden registrar y utilizar para predecir el tamaño de las grietas para garantizar que se realice el mantenimiento o el retiro antes de que falle alguna de las grietas. Los factores de seguridad se utilizan para reducir la vida de fatiga prevista a una vida de fatiga de servicio debido a la sensibilidad de la vida de fatiga al tamaño y la forma de los defectos que inician grietas y la variabilidad entre la carga supuesta y la carga real experimentada por un componente.
La vida a fatiga se puede dividir en un período de iniciación y un período de crecimiento de la grieta. [1] Las ecuaciones de crecimiento de grietas se utilizan para predecir el tamaño de la grieta a partir de un defecto inicial determinado y generalmente se basan en datos experimentales obtenidos de pruebas de fatiga de amplitud constante .
donde es la longitud de la grieta y es el crecimiento de la grieta por fatiga para un solo ciclo de carga . Se han desarrollado una variedad de ecuaciones de crecimiento de grietas similares a la ecuación de París-Erdogan para incluir factores que afectan la tasa de crecimiento de grietas, como la relación de tensión, las sobrecargas y los efectos del historial de carga.
El rango de intensidad de tensión se puede calcular a partir de la intensidad de tensión máxima y mínima para un ciclo.
Se utiliza un factor de geometría para relacionar la tensión del campo lejano con la intensidad de la tensión de la punta de la grieta usando
.
Existen referencias estándar que contienen los factores de geometría para muchas configuraciones diferentes. [3] [4] [5]
Historia de las ecuaciones de propagación de grietas.
A lo largo de los años se han propuesto muchas ecuaciones de propagación de grietas para mejorar la precisión de la predicción e incorporar una variedad de efectos. Los trabajos de Head, [6] Frost y Dugdale, [7] McEvily e Illg, [8] y Liu [9] sobre el comportamiento de crecimiento de grietas por fatiga sentaron las bases de este tema. La forma general de estas ecuaciones de propagación de grietas se puede expresar como
donde, la longitud de la grieta se denota por , el número de ciclos de carga aplicados está dado por , el rango de tensión por y los parámetros del material por . Para configuraciones simétricas, la longitud de la grieta desde el eje de simetría se define como y es la mitad de la longitud total de la grieta .
Las ecuaciones de crecimiento de grietas de la forma no son una verdadera ecuación diferencial ya que no modelan el proceso de crecimiento de grietas de manera continua durante todo el ciclo de carga. Como tal, se requieren algoritmos de identificación o recuento de ciclos separados, como el algoritmo de recuento de flujo de lluvia comúnmente utilizado , para identificar los valores máximo y mínimo en un ciclo. Aunque se desarrolló para los métodos de tensión/deformación-vida, el recuento de flujo de lluvia también ha demostrado funcionar para el crecimiento de grietas. [10] También se han desarrollado un pequeño número de ecuaciones derivadas verdaderas de crecimiento de grietas por fatiga. [11] [12]
Factores que afectan la tasa de crecimiento de las grietas.
Regímenes
La Figura 1 muestra un gráfico típico de la tasa de crecimiento de la grieta en función de la intensidad de la tensión alterna o la fuerza impulsora de la punta de la grieta representada en escalas logarítmicas. El comportamiento de la tasa de crecimiento de las grietas con respecto a la intensidad de la tensión alterna se puede explicar en diferentes regímenes (ver figura 1) de la siguiente manera
Régimen A: A tasas de crecimiento bajas, las variaciones en la microestructura , la tensión media (o relación de carga) y el medio ambiente tienen efectos significativos en las tasas de propagación de grietas. Se observa que en relaciones de carga bajas la tasa de crecimiento es más sensible a la microestructura y en materiales de baja resistencia es más sensible a la relación de carga. [13]
Régimen B: En el rango medio de tasas de crecimiento, las variaciones en la microestructura, la tensión media (o relación de carga), el espesor y el entorno no tienen efectos significativos sobre las tasas de propagación de grietas.
Régimen C: A altas tasas de crecimiento, la propagación de grietas es muy sensible a las variaciones en la microestructura, la tensión media (o relación de carga) y el espesor. Los efectos ambientales tienen relativamente menos influencia.
Efecto de la relación de estrés
Los ciclos con una relación de tensión más alta tienen una mayor tasa de crecimiento de grietas. [14] Este efecto a menudo se explica utilizando el concepto de cierre de grietas , que describe la observación de que las caras de las grietas pueden permanecer en contacto entre sí con cargas superiores a cero. Esto reduce el rango del factor de intensidad de tensión efectivo y la tasa de crecimiento de grietas por fatiga. [15]
Efectos de secuencia
Una ecuación da la tasa de crecimiento para un solo ciclo, pero cuando la carga no es de amplitud constante, los cambios en la carga pueden conducir a aumentos o disminuciones temporales en la tasa de crecimiento. Se han desarrollado ecuaciones adicionales para tratar algunos de estos casos. La tasa de crecimiento se retrasa cuando ocurre una sobrecarga en una secuencia de carga. Estas cargas generan una zona plástica que puede retrasar la tasa de crecimiento. Dos ecuaciones notables para modelar los retrasos que ocurren mientras la grieta crece a través de la región de sobrecarga son: [16]
El modelo Wheeler (1972)
con
donde es la zona plástica correspondiente al i-ésimo ciclo que ocurre después de la sobrecarga y es la distancia entre la grieta y la extensión de la zona plástica en la sobrecarga.
El modelo Willenborg
Ecuaciones de crecimiento de grietas
Ecuación de umbral
Para predecir la tasa de crecimiento de grietas en la región cercana al umbral, se ha utilizado la siguiente relación [17]
Ecuación París-Erdoğan
Para predecir la tasa de crecimiento de grietas en el régimen intermedio, se utiliza la ecuación de París-Erdoğan [2]
ecuación de Forman
En 1967, Forman propuso la siguiente relación para tener en cuenta el aumento de las tasas de crecimiento debido a la relación de tensiones y al aproximarse a la tenacidad a la fractura [18]
Ecuación de McEvily-Groeger
McEvily y Groeger [19] propusieron la siguiente relación de ley de potencia que considera los efectos de valores altos y bajos de
.
ecuación nasgro
La ecuación de NASGRO se utiliza en los programas de crecimiento de crack AFGROW, FASTRAN y NASGRO. [20] Es una ecuación general que cubre la tasa de crecimiento más baja cerca del umbral y la tasa de crecimiento aumentada que se acerca a la tenacidad a la fractura , además de permitir el efecto de tensión media al incluir la relación de tensión . La ecuación de NASGRO es
donde , , , , y son los coeficientes de la ecuación .
donde es la tensión de flujo, es el módulo de Young y es una constante que normalmente se encuentra en el rango de 0,1 a 0,5.
Ecuación de Walker
Para tener en cuenta el efecto de la relación de tensiones, Walker sugirió una forma modificada de la ecuación de París-Erdogan [22]
donde, es un parámetro del material que representa la influencia de la relación de tensiones en la tasa de crecimiento de las grietas por fatiga. Normalmente, toma un valor cercano a , pero puede variar entre . En general, se supone que la porción de compresión del ciclo de carga no tiene ningún efecto sobre el crecimiento de la grieta considerando lo que da. Esto puede explicarse físicamente considerando que la grieta se cierra con carga cero y no se comporta como una grieta bajo cargas de compresión. En materiales muy dúctiles como el acero Man-Ten, la carga de compresión contribuye al crecimiento de grietas según . [23]
ecuación de elber
Elber modificó la ecuación de Paris-Erdogan para permitir el cierre de grietas con la introducción del nivel de intensidad de la tensión de apertura en el que se produce el contacto. Por debajo de este nivel no hay movimiento en la punta de la grieta y, por tanto, no hay crecimiento. Este efecto se ha utilizado para explicar el efecto de la relación de tensiones y la mayor tasa de crecimiento observada con grietas cortas. La ecuación de Elber es [16]
donde, y son parámetros materiales. Con base en diferentes mecanismos de avance de grietas y protección de puntas de grietas en metales, cerámicas e intermetálicos , se observa que la tasa de crecimiento de grietas por fatiga en metales depende significativamente del término, en cerámicas de y los intermetálicos tienen una dependencia casi similar de los términos y .
Predicción de la vida por fatiga.
Programas de computador
Existen muchos programas informáticos que implementan ecuaciones de crecimiento de crack, como Nasgro , [24] AFGROW y Fastran . Además, también existen programas que implementan un enfoque probabilístico del crecimiento de grietas que calculan la probabilidad de falla a lo largo de la vida útil de un componente. [25] [26]
Los programas de crecimiento de grietas hacen crecer una grieta desde un tamaño de falla inicial hasta que excede la tenacidad a la fractura de un material y falla. Debido a que la tenacidad a la fractura depende de las condiciones límite, la tenacidad a la fractura puede cambiar desde condiciones de deformación plana para una grieta de superficie semicircular a condiciones de tensión plana para una grieta pasante. La tenacidad a la fractura en condiciones de tensión plana suele ser dos veces mayor que en condiciones de deformación plana. Sin embargo, debido a la rápida tasa de crecimiento de una grieta cerca del final de su vida, las variaciones en la tenacidad a la fractura no alteran significativamente la vida de un componente.
Los programas de crecimiento de crack suelen ofrecer una opción de:
Métodos de conteo cíclico para extraer los extremos del ciclo.
Factores geométricos que seleccionan la forma de la grieta y la carga aplicada.
ecuación de crecimiento de grietas
modelos de aceleración/retraso
Propiedades del material como límite elástico y tenacidad a la fractura.
Solucion analitica
El factor de intensidad del estrés está dado por
donde es la tensión de tracción uniforme aplicada que actúa sobre la muestra en la dirección perpendicular al plano de la grieta, es la longitud de la grieta y es un parámetro adimensional que depende de la geometría de la muestra. La intensidad del estrés alternante se vuelve
donde es el rango de amplitud de la tensión cíclica.
Suponiendo que el tamaño inicial de la grieta es , el tamaño crítico de la grieta antes de que falle la muestra se puede calcular usando como
La ecuación anterior es de naturaleza implícita y puede resolverse numéricamente si es necesario.
Caso I
Porque el cierre de grietas tiene un efecto insignificante sobre la tasa de crecimiento de las grietas [27] y la ecuación de París-Erdogan se puede utilizar para calcular la vida a fatiga de una muestra antes de que alcance el tamaño crítico de la grieta como
Modelo de crecimiento de crack con valor constante de 𝛽 y R = 0
Figura 2: Representación geométrica de la muestra de ensayo de tensión con fisura central
Para el modelo de crecimiento de grietas de Griffith-Irwin o grieta central de longitud en una lámina infinita como se muestra en la figura 2, tenemos y es independiente de la longitud de la grieta. Además, se puede considerar que es independiente de la longitud de la grieta. Al asumir la integral anterior se simplifica a
Al integrar la expresión anterior para los casos y , el número total de ciclos de carga viene dado por
Ahora bien, para que un tamaño de grieta crítico sea muy grande en comparación con el tamaño de grieta inicial, se obtendrá
Las expresiones analíticas anteriores para el número total de ciclos de carga hasta la fractura se obtienen asumiendo . Para los casos en los que depende del tamaño de la grieta, como las geometrías de tensión de muesca de borde único (SENT) y tensión de fisura central (CCT), se puede utilizar la integración numérica para calcular .
Caso II
Porque el fenómeno de cierre de grietas tiene un efecto en la tasa de crecimiento de las grietas y podemos invocar la ecuación de Walker para calcular la vida a fatiga de una muestra antes de que alcance el tamaño crítico de la grieta como
Cálculo numérico
Figura 3: Representación esquemática del proceso de predicción de la vida por fatiga [28]
Este esquema es útil cuando depende del tamaño de la grieta . Se considera que el tamaño inicial de la grieta es . El factor de intensidad de tensión en el tamaño de grieta actual se calcula utilizando la tensión máxima aplicada como
Si es menor que la tenacidad a la fractura , la grieta no ha alcanzado su tamaño crítico y la simulación continúa con el tamaño de grieta actual para calcular la intensidad de la tensión alterna como
Ahora, sustituyendo el factor de intensidad de tensión en la ecuación de París-Erdogan, el incremento en el tamaño de la grieta se calcula como
¿Dónde está el tamaño del paso del ciclo? El nuevo tamaño de grieta se vuelve
donde índice se refiere al paso de iteración actual. El nuevo tamaño de grieta se utiliza para calcular la intensidad de la tensión en la tensión máxima aplicada para la siguiente iteración. Este proceso iterativo continúa hasta
Una vez que se cumple este criterio de falla, se detiene la simulación.
La representación esquemática del proceso de predicción de la vida por fatiga se muestra en la figura 3.
Ejemplo
Figura 4: Representación geométrica de la muestra de prueba de tensión de muesca de borde único
El factor de intensidad de tensión en una muestra SENT (ver figura 4) bajo crecimiento de grietas por fatiga viene dado por [5]
Para el cálculo se consideran los siguientes parámetros
mm, mm, mm, , ,
MPa, , .
La longitud crítica de la grieta, , se puede calcular cuando :
Resolviendo la ecuación anterior, la longitud crítica de la grieta se obtiene como .
Ahora bien, invocando la ecuación París-Erdogan se obtiene
Por integración numérica de la expresión anterior, el número total de ciclos de carga hasta la falla se obtiene como .
Referencias
^ Schijve, J. (enero de 1979). "Cuatro conferencias sobre el crecimiento de grietas por fatiga". Ingeniería Mecánica de Fracturas . 11 (1): 169–181. doi :10.1016/0013-7944(79)90039-0. ISSN 0013-7944.
^ ab París, ordenador personal; Erdogan, F. (1963). "Un análisis crítico de las leyes de propagación de grietas". Revista de Ingeniería Básica . 18 (4): 528–534. doi :10.1115/1.3656900..
^ Murakami, Y.; Aoki, S. (1987). Manual de factores de intensidad del estrés . Pérgamo, Oxford.
^ Rooke, DP; Cartwright, DJ (1976). Compendio de Factores de Intensidad del Estrés . Oficina de papelería de Su Majestad, Londres.
^ ab Tada, Hiroshi; París, Paul C.; Irwin, George R. (1 de enero de 2000). Manual de análisis de tensión de grietas (tercera ed.). Three Park Avenue Nueva York, NY 10016-5990: ASME. doi : 10.1115/1.801535. ISBN0791801535.{{cite book}}: CS1 maint: location (link)
^ Head, AK (septiembre de 1953). "El crecimiento de las grietas por fatiga". Revista filosófica y revista científica de Londres, Edimburgo y Dublín . 44 (356): 925–938. doi :10.1080/14786440908521062. ISSN 1941-5982.
^ Escarcha, NE; Dugdale, DS (enero de 1958). "La propagación de grietas por fatiga en probetas de chapa". Revista de Mecánica y Física de Sólidos . 6 (2): 92-110. Código bibliográfico : 1958JMPSo...6...92F. doi :10.1016/0022-5096(58)90018-8. ISSN 0022-5096.
^ McEvily, Arthur J.; Ilg, Walter (1960). "Un método para predecir la tasa de propagación de grietas por fatiga". Simposio sobre fatiga de estructuras de aeronaves . ASTM Internacional. págs. 112–112–8. doi :10.1520/stp45927s. ISBN9780803165793.
^ Liu, HW (1961). "Propagación de grietas en láminas de metal delgadas bajo carga repetida". Revista de Ingeniería Básica . 83 (1): 23–31. doi : 10.1115/1.3658886. hdl : 2142/111864 . ISSN 0021-9223.
^ Dividir, R.; Seetharam, SA; Bhaskaran, TA (1984). "Recuento de ciclos para el análisis del crecimiento de grietas por fatiga". Revista Internacional de Fatiga . 6 (3): 147-156. doi :10.1016/0142-1123(84)90032-X.
^ Pommier, S.; Risbet, M. (2005). "Ecuaciones derivadas del tiempo para el crecimiento de grietas por fatiga en metales en modo I". Revista Internacional de Fatiga . 27 (10-12): 1297-1306. doi :10.1016/j.ijfatigue.2005.06.034.
^ Lu, Zizi; Liu, Yongming (2010). "Análisis del crecimiento de grietas por fatiga a pequeña escala". Revista Internacional de Fatiga . 32 (8): 1306-1321. doi :10.1016/j.ijfatigue.2010.01.010.
^ Ritchie, RO (1977). "Propagación de grietas por fatiga casi umbral en acero de ultra alta resistencia: influencia de la relación de carga y la resistencia cíclica". Revista de Tecnología y Materiales de Ingeniería . 99 (3): 195-204. doi : 10.1115/1.3443519. ISSN 0094-4289. S2CID 136642892.
^ Maddox, SJ (1975). "El efecto de la tensión media en la propagación de grietas por fatiga: una revisión de la literatura". Revista Internacional de Fractura . 1 (3).
^ Elber, W. (1971), "La importancia del cierre de grietas por fatiga", Tolerancia al daño en estructuras de aeronaves , ASTM International, págs. 230-242, doi :10.1520/stp26680s, ISBN9780803100312
^ ab Suresh, S. (2004). Fatiga de Materiales . Prensa de la Universidad de Cambridge. ISBN978-0-521-57046-6.
^ Allen, RJ; Stand, GS; Jutla, T. (marzo de 1988). "Una revisión de la caracterización del crecimiento de grietas por fatiga mediante la mecánica de fractura elástica lineal (LEFM). Parte II: documentos de asesoramiento y aplicaciones dentro de los estándares nacionales". Fatiga y fractura de estructuras y materiales de ingeniería . 11 (2): 71–108. doi :10.1111/j.1460-2695.1988.tb01162.x. ISSN 8756-758X.
^ Forman, RG; Kearney, VE; Engle, RM (1967). "Análisis numérico de la propagación de grietas en estructuras cargadas cíclicamente". Revista de Ingeniería Básica . 89 (3): 459–463. doi : 10.1115/1.3609637. ISSN 0021-9223.
^ McEvily, AJ; Groeger, J. (1978), "En el umbral del crecimiento de grietas por fatiga", Avances en la investigación sobre la resistencia y la fractura de materiales , Elsevier, págs. 1293–1298, doi :10.1016/b978-0-08-022140-3.50087 -2, ISBN9780080221403
^ Forman, RG; Shivakumar, V.; Cardenal, JW; Williams, LC; McKeighan, ordenador personal (2005). "Base de datos de crecimiento de grietas por fatiga para análisis de tolerancia a daños" (PDF) . FAA . Consultado el 6 de julio de 2019 .
^ ab Ritchie, RO (1 de noviembre de 1999). "Mecanismos de propagación de grietas por fatiga en sólidos dúctiles y frágiles". Revista Internacional de Fractura . 100 (1): 55–83. doi :10.1023/A:1018655917051. ISSN 1573-2673. S2CID 13991702.
^ Walker, K. (1970), "El efecto de la relación de tensiones durante la propagación de grietas y la fatiga para el aluminio 2024-T3 y 7075-T6", Efectos del entorno y el historial de cargas complejas sobre la vida útil de la fatiga , ASTM International, págs. 1-14 , doi :10.1520/stp32032s, ISBN9780803100329
^ Dowling, Norman E. (2012). Comportamiento mecánico de materiales: métodos de ingeniería para deformación, fractura y fatiga . Pearson. ISBN978-0131395060. OCLC 1055566537.
^ "Software de crecimiento de grietas por fatiga y mecánica de fracturas NASGRO®". 26 de septiembre de 2016 . Consultado el 14 de julio de 2019 .
^ "Actualización del programa informático de probabilidad de fractura (PROF) para el análisis de riesgo de aeronaves envejecidas. Volumen 1: Modificaciones y guía del usuario" . Consultado el 14 de julio de 2019 .
^ "Software de evaluación de confiabilidad y mecánica de fracturas DARWIN". 14 de octubre de 2016 . Consultado el 14 de julio de 2019 .
^ Zehnder, Alan T. (2012). Mecánica de fracturas . Apuntes de conferencias sobre mecánica aplicada y computacional. vol. 62. Dordrecht: Springer Países Bajos. doi :10.1007/978-94-007-2595-9. ISBN9789400725942.
^ "Crecimiento de grietas por fatiga" . Consultado el 6 de julio de 2019 .
enlaces externos
Forman, RG; Shivakumar, V.; Cardenal, JW; Williams, LC; McKeighan, ordenador personal (2005). "Base de datos de crecimiento de grietas por fatiga para análisis de tolerancia a daños" (PDF) . FAA . Consultado el 6 de julio de 2019 .
Gallagher, JP; Giessler, FJ; Berens, AP; Engle, Jr, JM "Manual de diseño tolerante a daños de la USAF: Directrices para el análisis y diseño de estructuras de aeronaves tolerantes a daños. Revisión B". Archivado desde el original el 9 de julio de 2019 . Consultado el 9 de julio de 2019 .
"Manual de evaluación de tolerancia a daños, volumen I: introducción, mecánica de fracturas, propagación de grietas por fatiga" (PDF) . Administración Federal de Aviación. 1993 . Consultado el 16 de julio de 2019 .