Abstract:
Los polímeros se han convertido en el material por excelencia en la industria de la distribución de gas. Por ello, el desarrollo de nuevas resinas para la fabricación de tuberías plásticas presurizadas con fluidos ha reavivado el interés por la comprensión y modelado del fenómeno de fractura conocido como Propagación Rápida de Fisuras (RCP).
El RCP se caracteriza por el crecimiento de una grieta frágil en estado cuasi-estacionario, capaz de propagarse a velocidades sónicas por todo el material, y producir consecuencias catastróficas. Este tipo de fractura tiene lugar bajo determinadas condiciones críticas de presión y temperatura, y suele iniciarse por la acción de fuentes externas de daño, como un impacto, o por el crecimiento lento de fisuras (SCG). Debido a la complejidad que entraña la prevención de la iniciación de grietas desde el diseño, el enfoque está puesto en que dichas condiciones críticas no se alcancen nunca. Para ello se han diseñado dos pruebas, ahora estandarizadas por la Organización Internacional de Normalización (ISO). La primera de ellas se trata de un costoso ensayo a escala real (FST), mientras que el segundo es un ensayo a escala de laboratorio (S4T) diseñado inicialmente con fines de investigación por P. S. Leevers y colaboradores. Ambas pruebas intentan simular las condiciones de operación para determinar la presión crítica, por debajo de la cual, el crecimiento rápido es insostenible y la grieta se detiene. Sin embargo, la correlación entre sus resultados no es directa y todavía hoy no se comprenden bien las causas de esta discordancia. Además, aún no se ha desarrollado un método para medir la resistencia a la fractura del material de la tubería. Esta propiedad debe determinarse mediante un análisis profundo de los datos obtenidos a través de la versión instrumentada de estos experimentos.
Este trabajo describe el estudio teórico del RCP y el desarrollo de una novedosa herramienta computacional ideada para simularlo, y ayudar a la interpretación de los resultados de la prueba S4T. El modelo propuesto recoge los tres aspectos físicos principales del fenómeno de una manera sencilla: (--) la deformación de las paredes de la tubería, (----) el escape transitorio de gas a través de la grieta abierta, (----) el comportamiento de fractura del material, y sus efectos acoplados.
Aprovechando su naturaleza cuasi-estacionaria, el análisis se centra en una porción reducida de la tubería, un volumen de control que encierra la punta de la grieta y se traslada con ella a la misma velocidad. Este volumen se divide en secciones anulares contiguas cuya deformación, causada por la fractura y la fuga de fluido, se representa con la flexión y extensión axial de una viga curva. Las condiciones de transición entre dos anillos adyacentes se consiguen introduciendo una restricción elástica traslacional cuya rigidez varía axialmente junto con el perfil de distribución de la presión. La magnitud de la presión, la tasa de crecimiento y la geometría de la fisura medidos en un ensayo S4T instrumentado, junto con las propiedades del material y de la tubería, representan los parámetros de entrada.
El modelo mecánico utilizado para simular la fractura se deriva de un arco con una restricción traslacional externa en la ubicación de la grieta. El problema de contorno, sus ecuaciones diferenciales y las condiciones de transición se obtienen con la aplicación del Cálculo de Variaciones. El método de Ritz se implementa para obtener resultados numéricos del campo de desplazamiento. Este último se utiliza para realizar un balance de energía dentro del volumen de control que permite la cuantificación de la fuerza impulsora para la propagación. Esta magnitud puede utilizarse para caracterizar la resistencia a la fractura del polímero, siempre que los datos de entrada correspondan a las condiciones críticas.
Description:
Polymers have become the material par excellence in the industry of gas distribution. Therefore, the development of new resins for the manufacture of fluid-pressurized plastic pipes has revived the interest in the understanding and modeling of the fracture phenomenon known as Rapid Crack Propagation (RCP).
The RCP is characterized by the growth of a fragile crack in a quasi-stationary state, capable of propagating at sonic speeds throughout the material, thus producing catastrophic consequences. It takes place under certain critical conditions of pressure and temperature and is usually initiated by the action of external sources of damage, such as impacts, or Slow Crack Growth (SCG). Due to the complexity involved in the prevention of crack initiation from design, the focus is on preventing that such critical conditions are never reached. Two tests have been designed for this purpose, now standardized by the International Organization for Standardization (ISO). The first one is an expensive full-scale test (FST) while the second is a laboratory-scale test (S4 Small Scale Steady State Test), initially designed for research purposes by P. S. Leevers and collaborators. Both tests attempt to simulate the operating conditions to determine the critical pressure, below which, the rapid growth is untenable and the crack is arrested. However, the correlation between their results is not direct and still today, the causes of this discordance are not well understood. Furthermore, a method to measure the fracture toughness of the pipe material has not been developed yet. This property has to be determined by a profound analysis of the data obtained through the instrumented version of these experiments.
This paper describes the theoretical study of RCP and a novel computational tool developed to simulate it and to assist the interpretation of the S4 test results. It compiles the three physical aspects of the complex phenomenon in a simple manner: (𝑖) the deformation of the pipe walls, (𝑖𝑖) the transient scape of gas through the opened crack, (𝑖𝑖𝑖) the material’s fracture behavior, and their coupled effects. Taking advantage of the quasi-stationary nature of the phenomenon, the analysis focuses on a reduced portion of the pipe, a control volume that encloses the crack tip and moves along with it. This volume is divided into contiguous ring sections whose deformation caused by the fracture and fluid leakage is modeled as a curved beam. The transition conditions between two adjacent rings are achieved by introducing an elastic translational restraint whose stiffness varies axially together with the pressure distribution profile. Both the spring stiffness, the magnitude of the pressure, and the growth rate are input parameters obtained from an S4 test.
The mechanical model used to simulate the fracture process is derived from an arc with an internal translational restraint in the crack location. The boundary value problem and the transition conditions are obtained with the application of the Calculus of Variations. The Ritz method is used to obtain numerical results of the displacement field. The latter is used to carry out an energy balance within the fracture control volume that allows the quantification of the crack driving force. This magnitude can be used to characterize the fracture toughness of the polymer, provided that the input data corresponds to the critical conditions.