Una singularidad que puede aparecer en la mayoría de los problemas de fractura tridimensional bajo comportamiento elástico es la singularidad de esquina o de borde libre, localizada en la intersección del frente de grieta con una frontera libre. Sin embargo, su efecto es habitualmente ignorado. Por otro lado, el método de los elementos finitos extendido, XFEM, es una técnica que permite el modelado numérico eficiente de problemas de fractura, para lo que incorpora la geometría de la grieta a través de funciones de enriquecimiento dentro de un modelo de elementos finitos, cuya malla ya no necesita adaptarse a la geometría de la fisura. Sin embargo, cuando se aplica al estudio de problemas de fractura que presentan singularidad de esquina, esta no se captura adecuadamente. La razón es que las funciones de enriquecimiento del XFEM solo describen la singularidad típicamente asociada al frente de grieta. Por tanto, para poder modelar la singularidad del borde libre es necesario su introducción en el enriquecimiento. En este trabajo presentamos un nuevo conjunto de funciones de enriquecimiento basadas en armónicos esféricos que consiguen introducir el comportamiento asociado a la singularidad de borde libre.
One type of singularity that may appear in a three-dimensional fracture problem under elastic behavior is the free corner singularity, which occurs at the intersection of the crack with a free boundary and whose effect is usually ignored. The extended finite element method (XFEM) is a technique that allows the efficient numerical modeling of fracture problems, by using enrichment functions within a finite element model that incorporate the geometry and effect of the crack. However, when applied to the study of a problem with corner singularity, the singular behavior is not properly captured. The reason is that the usual enrichment in the XFEM only describes the typical crack front singularity. Hence, in order to include the effect of the free border singularity, the enrichment has to be modified. In this work, we present a new set of functions which relies on spherical harmonics and is able to capture the behavior of the free border singularity.
El estudio de la singularidad que aparece en la intersección del frente de grieta con una superficie libre bajo comportamiento elástico, conocida como singularidad de esquina, de borde libre o de vertex, es de gran interés. Sin embargo, a pesar de diversos estudios previos [1], aún no se han logrado resultados teóricos sobre el estado tensional existente. Además, las únicas formulaciones disponibles están limitadas a casos particulares [2–5].
Debido a la falta de resultados analíticos o experimentales que permitan el desarrollo de un modelo aproximado, la consideración de este efecto se suele realizar mediante métodos numéricos [6–11]. Aun así, debido a las complicaciones asociadas al modelado y a la discretización mediante elementos finitos, lo habitual es ignorar su efecto, ya que su zona de influencia es pequeña en relación con la zona de influencia de la singularidad del frente de grieta [6,7,12]. Sin embargo, su presencia modifica el comportamiento en fatiga [9,12–14] y, por tanto, debe tenerse en cuenta para un estudio adecuado.
Por otro lado, el método de los elementos finitos extendido –XFEM– es una técnica numérica que simplifica el estudio de problemas de mecánica de la fractura [15–17]. No obstante, cuando se utiliza el XFEM para estudiar una grieta tridimensional que presenta singularidades de esquina, las funciones de enriquecimiento usuales solo describen la singularidad típica del frente de grieta y no permiten capturar correctamente otras singularidades. Para describir mejor el efecto de esquina, en este trabajo proponemos añadir un conjunto de nuevas funciones de enriquecimiento que puedan representar la singularidad de borde libre.
Habitualmente, para definir las funciones de enriquecimiento, utilizadas en un problema, se parte de una expresión analítica conocida del campo de desplazamientos o de tensiones. Sin embargo, no se conoce ninguna expresión analítica para el comportamiento asociado a la singularidad de borde libre, por lo que no se dispone de una expresión a partir de la cual se puedan construir las funciones de enriquecimiento. Por tanto, se ha intentado capturar el efecto de esta singularidad mediante un enriquecimiento genérico.
2Singularidad de borde libreLa singularidad de borde libre, que aparece en la intersección entre un frente de grieta tridimensional y una superficie sin restricciones, es un tipo de comportamiento singular propio de la mecánica de la fractura tridimensional y presenta un comportamiento que no se corresponde con la singularidad clásica de los campos de Williams presente en el frente de grieta, tal y como se acepta en mecánica de la fractura. Este problema y su influencia en la propagación por fatiga han motivado que desde los años setenta del siglo pasado se esté intentando, sin éxito, llegar a una expresión analítica para su aproximación [2–8,18–20]. Sin embargo, se acepta que, en las cercanías del vertex, el campo de tensiones muestra un comportamiento esférico y, así, se le asocia un sistema local de referencia esférico con coordenadas (ρ, θ, ϕ) centrado en el vertex (fig. 1), en contraposición al sistema local cilíndrico, de coordenadas (r, θ, z), usado en la descripción habitual de la grieta de mecánica de la fractura según los campos de Williams y definido de forma que el eje z coincida con el frente de grieta.
Designando σijbl al campo de tensiones para la singularidad de borde libre, este presenta la siguiente forma funcional [7]:
En esta relación, el término fij(θ, ϕ) corresponde a una función angular esférica sin expresión analítica conocida. Además, fij(θ, ϕ) y λ dependen de ν, es decir, del coeficiente de Poisson. Kbl se corresponde con el concepto intuitivo del factor de intensidad de tensiones, FIT, en dicho punto, aunque existen diferencias desde un punto de vista estricto. Hay que destacar que en problemas de modo i, Kbl es tal que, si el frente de grieta es recto e intersecta perpendicularmente con la superficie libre, el FIT en modo i para la zona donde hay vertex disminuye bruscamente. Este comportamiento es uno de los efectos más documentados de dicha singularidad [7,8].
La dependencia de λ con ν, para una grieta perpendicular a la superficie libre, se estudia en [4], obteniéndose los resultados de la tabla 1.
3Enriquecimiento para la singularidad de borde libreA continuación se propone un enriquecimiento para el modelado de la singularidad de borde libre basado en armónicos esféricos. El objetivo es introducir un enriquecimiento genérico para fenómenos cuya influencia se extienda esféricamente y aplicarlo al caso de la singularidad de borde libre.
3.1Breve introducción a los armónicos esféricosSe conoce como armónicos esféricos a la parte angular de un conjunto de soluciones de la ecuación de Laplace en coordenadas esféricas. Además, son funciones asociadas a los polinomios de Legendre y pueden generarse directamente a partir de ellos mediante recurrencias [21–23]. Las funciones normalizadas de armónicos esféricos denominados mediante Ylm(θ, ϕ) son:
siendo Plm(cosθ) los polinomios asociados de Legendre definidos por:
La representación habitual de las 4 primeras bandas de los armónicos esféricos reales se puede apreciar en la figura 2.
Los armónicos esféricos presentan relaciones de completitud, ortogonalidad y normalidad, constituyendo una base completa para las funciones con dependencia angular esférica. Por tanto, toda función que presenta solo dependencia angular del sistema de referencia esférico puede ser expresada en términos de armónicos esféricos, por ejemplo:
donde alm es el coeficiente real asociado con el armónico esférico Ylm(θ, ϕ) y que se debe determinar para cada problema.3.2Armónicos esféricos como enriquecimiento para la singularidad de borde libre
Si consideramos la ecuación (1) y la relación de la ecuación (4), entonces es posible utilizar los armónicos esféricos para intentar reproducir la parte angular del campo de desplazamientos de la singularidad. Consecuentemente, es posible introducir un nuevo conjunto de funciones de enriquecimiento relacionadas con el campo de desplazamientos, como:
Así se puede expresar el desplazamiento en la aproximación de XFEM para simular la existencia de grietas, incluyendo la singularidad de borde libre o esquina, por:
donde los subconjuntos J y K contienen los nodos enriquecidos de la aproximacion clásica de XFEM, V contiene los nodos enriquecidos con armónicos esféricos, y ai, bi y cilm son los correspondientes grados de libertad de cada enriquecimiento. Los grados de libertad ui se corresponden con los desplazamientos nodales. La selección del conjunto V se realiza tomando los elementos que incluyan el punto vertex y cuyos nodos pertenezcan al subconjunto K. El número de términos de armónicos esféricos añadidos se controla con L.
Hay algunos aspectos delicados para la construcción y el enriquecimiento del subconjunto V. La primera dificultad es la elección del orden de singularidad, siendo λ el indicado en la tabla 1. La segunda complicación es que la influencia de la singularidad de borde libre únicamente se aprecia en una región pequeña y cercana al punto de singularidad, cuyo tamaño es desconocido e implica limitaciones respecto al tamaño de elemento.
4Verificación numérica4.1Modelo geométrico y problema de referenciaPara la verificación se considera un caso tridimensional en modo i. La geometría corresponde a una placa rectangular de espesor t, anchura w y altura 2h (fig. 3), con una grieta de longitud a. La placa tiene las siguientes dimensiones: a/t=0,5, w/t=1,5 y h/t=2,5. El coeficiente de Poisson es 0,3 y el módulo de Young es 20.000 unidades. Este problema se utiliza para estudiar la singularidad de borde libre en [24], tomando el valor de λ=0,5. El cálculo de los FIT de un problema similar aparece en [25] y, además, es empleado como problema de referencia en [15] para la verificación del XFEM. Para el estudio se considera un sistema coordenado cartesiano centrado en el frente de grieta y localizado en el plano medio del espesor del modelo geométrico.
4.2Modelos de elementos finitos
Se han construido 2 modelos. El primero, realizado en ABAQUS, utiliza una malla estándar de elementos finitos (fig. 4) y gracias a las simetrías, presentes en el plano de grieta y en el plano medio, solo ha sido necesario considerar un cuarto de la placa. La malla presenta 2 regiones: una que está adaptada a la singularidad y a la morfología de la grieta, concentrada hacia el borde libre, y la otra que permite capturar el comportamiento general del modelo. Se han utilizado elementos tipo hexaedro y prisma de orden cuadrático, con 20 nodos por elemento, siendo el tamaño mínimo de elemento del orden de 0,002a. Este modelo se utiliza para obtener una solución de referencia con la que comparar los demás resultados.
El segundo modelo de elementos finitos lineales se ha creado para la aproximación de XFEM y utiliza una malla regular de hexaedros (fig. 5). Solo se ha aplicado la simetría que aparece en el plano medio, analizando una secuencia de mallas para observar el efecto del enriquecimiento, con un tamaño de elemento respecto a la longitud de grieta de 0,125a, 0,083a y 0,062a.
4.3Resultados para los factores de intensidad de tensiones
Los resultados obtenidos con los 2 modelos de elementos finitos se comparan con la solución que aparece en [24]. Se realiza un análisis usando 3 tamaños de elemento en la malla de XFEM. El tamaño de elemento se elige observando la solución de referencia y la obtenida mediante ABAQUS. En estas soluciones la región de singularidad de esquina es alrededor del 3% de la longitud del frente de grieta, t. Se construyen 3 mallas: una con un tamaño de elemento más grande que la región de singularidad (0,125a), otra con un tamaño de elemento similar (0,083a) y otra, la tercera, con un tamaño de elemento inferior (0,062a).
Este problema presenta una configuración de modo i puro. Se puede, por tanto, utilizar la integral J para extraer el FIT. Los resultados de los FIT para diversos órdenes de armónicos esféricos, seleccionados mediante L, se presentan en las figuras 6, 7 y 8. En estas gráficas se puede apreciar que el comportamiento de la singularidad de borde libre en los fit se captura mejor si se usan algunos términos de enriquecimiento mediante armónicos esféricos.
En la zona de la intersección de la grieta con el borde libre, se detecta una caída brusca en los FIT para la solución de elementos finitos muy refinada que se toma como referencia. Sin embargo, el FIT obtenido mediante XFEM sin armónicos esféricos no muestra tal comportamiento, al igual que si solo se introduce el primer término de armónicos esféricos. De hecho, ambos valores coinciden en las gráficas.
A medida que se añaden más términos, correspondientes a armónicos esféricos, esta caída brusca también aparece, como se puede observar a partir de L=1. Más aún, a partir de L=3 el efecto de la mejora se estabiliza. Por tanto, únicamente es necesaria la introducción de un número finito de términos para mejorar la solución. Este resultado está de acuerdo con la experiencia en otros ámbitos tratados con desarrollos en armónicos esféricos [21–23].
Se observa, también, que el enriquecimiento con armónicos esféricos introduce comportamientos no deseados. Este comportamiento se aprecia en el aumento incorrecto del FIT justo antes de la zona afectada por la singularidad de borde libre, z/t≈0,40−0,45, siendo mayor para tamaños de elemento muy grandes y disminuyendo conforme se refina. Se trata de un efecto de la zona de transición, es decir, de los elementos no totalmente enriquecidos, fenómeno que también se encuentra para el enriquecimiento habitual de XFEM.
Hay que destacar que, dado que el enriquecimiento con armónicos esféricos simplemente constituye una base para efectos esféricos, se puede utilizar para introducir cualquier efecto local con comportamiento esférico, para lo cual, habitualmente, solo es necesario considerar un número finito de términos de enriquecimiento. Sin embargo, hay que considerar los problemas en la zona de transición. Existe también la posibilidad de que se introduzca una contaminación en la solución para mallas no suficientemente refinadas si la zona enriquecida es mucho mayor que la región afectada por la singularidad. Estos efectos pueden suponer limitaciones desde un punto de vista práctico.
5ConclusionesEn este trabajo se ha propuesto un nuevo enriquecimiento capaz de modelar y capturar parte del comportamiento de la singularidad de esquina. El comportamiento en el factor de intensidad de tensiones en la región cercana a la singularidad de esquina muestra una mejora respecto al enriquecimiento estándar de XFEM. Más aún, para esta mejora en los resultados solo se necesitan unos pocos términos del nuevo enriquecimiento. Hay que añadir que, al ser la aplicación de una base genérica para efectos esféricos, también constituye una opción natural para la simulación de otros efectos esféricos, como inclusiones, mediante una aproximación de XFEM.
Sin embargo, también se puede observar que aparece un comportamiento no deseado en la zona de transición entre la región dominada por la singularidad del frente de la grieta y la zona dominada por la singularidad de esquina. Dado que la región exacta de influencia del borde libre y el orden de la singularidad no se conocen con exactitud, es necesario seguir investigando para lograr una mejor comprensión de la singularidad en el punto de vertex. No obstante, el nuevo enriquecimiento también puede ser utilizado como una herramienta con este objetivo: puede emplearse para aumentar la información sobre la parte angular –comparando los grados de libertad asociados a cada término de enriquecimiento– y como función de prueba para realizar estudios sobre los diferentes órdenes de singularidad, permitiendo hacer comparaciones y análisis de la influencia de cada expresión, lo que será objeto de estudio en trabajos futuros.