Análisis de sensibilidad paramétrica sobre un modelo computacional XFEM para la propagación de grietas en una probeta CT de acero de fase dual

Abstract

The implementation of the extended finite element method (XFEM) has allowed modeling without increasing computational costs and minimizing the use of excessively fine meshes, discontinuities, complex geometric singularities, crack propagation processes, among others. Being a recent method, the effects of the input parameters on the response variables of the model are unknown. For this reason, the sensitivity analysis is relevant to determine the effect of various parameters involved in the modeling of crack propagation. The objective of this work is to carry out a parametric sensitivity analysis for an XFEM computational model of the propagation of cracks in a test piece for fracture toughness of the compact tension (CT) type made of DP600 steel, minimizing the number of nodes to 1000, for which five parameters are selected that vary in two conditions each. The selected parameters are the model dimensions, finite element size, material behavior, evolution, and damage tolerance in the material. As response variables, the Force vs. Displacement diagram and the crack trajectory were used.

For statistical analysis, the Pareto graph, main effect graphs, and an analysis of variance were used, to determine the parameter with the most significant effect in the models. To quantify the error of the two model responses for the experimental data, a relative error was used between the strength of the computational model for the experimental results and a mean squared error (RMSE), totaling the error in only one for the Force vs. Displacement, in the case of two-dimensional models the average RMS errors were 46.52 %, while the three-dimensional models were closer to the experimental data with average RMS errors of 2.20 %. In the case of crack growth, RMS errors below 10 % were obtained. Through the statistical analyzes carried out, a high statistical significance was established in the parameter of the model type, being the main one for the decrease in the error calculated for experimental cases.

La implementación del método de elementos finitos extendidos (XFEM) ha permitido modelar sin aumentar costos computacionales y minimizando el uso de mallas excesivamente finas las discontinuidades, singularidades geométricas complejas, procesos de propagación de grietas, entre otros. Al ser un método reciente, se desconoce los efectos de los parámetros de entrada sobre las variables de respuesta del modelo. Por tal motivo, el análisis de sensibilidad toma relevancia para determinar el efecto de varios parámetros que intervienen en el modelado de la propagación de grietas. El objetivo de este trabajo es realizar un análisis de sensibilidad paramétrica para un modelo computacional por XFEM de la propagación de grietas en una probeta para tenacidad a la fractura del tipo compact tension (CT) fabricada en acero DP600 minimizando la cantidad de nodos a 1000, para el cual se seleccionó cinco parámetros que varían en dos condiciones cada uno. Los parámetros seleccionados son las dimensiones del modelo, tamaño del elemento finito, comportamiento del material, evolución y tolerancia del daño en el material. Como variables de respuesta se usó el diagrama de Fuerza Vs. Desplazamiento y la trayectoria de la grieta.

Para análisis estadadistico se utilizó el gráfico de Pareto, gráficas de efectos principales y un análisis de varianza, con el fin de determinar el parámetro con el efecto más significativo en los modelos. Para cuantificar el error de las dos respuestas del modelo, se utilizó un error relativo entre la fuerza del modelo computacional respecto a los resultados experimentales y un error cuadrático medio (RMSE) totalizando el error en uno solo para el diagrama de Fuerza vs. Desplazamiento, en el caso de los modelos de dos dimensiones los errores RMS en promedio fue del 46.52 %, mientras que los modelos en tres dimensiones fue más cercanos a los datos experimentales con errores RMS en promedio de 2.20 %. En el caso del crecimiento de grietas, se obtuvo errores RMS por debajo del 10 %. A través de los análisis estadísticos realizados, se estableció significancia estadadística alta en el parámetro del tipo modelo, siendo principal para la disminución del error calculado respecto a casos experimentales.