Abstract
Accurately measuring a crack length is a crucial aspect of experimental fracture tests. In this work, we present an innovative application of the A* (A-star) shortest path methodology to track different shapes of cracks from numerical simulations. This approach is highly efficient, significantly improving the speed and accuracy of crack length measurements. Furthermore, we introduce a modified weight cost function that follows the crack path in the damage field, enhancing the accuracy of our method. The effectiveness of the proposed procedure is shown by fabricating damage fields with different geometry and good agreement when compared to the exact values. In addition, we evaluate a time-dependent crack propagation case, achieving high accuracy. We present all features and steps of the procedure to showcase its efficacy in accurately measuring the length of a crack path. Finally, we validate our method using a phase-field fracture framework and compare it with the compliance technique. The results show that the proposed method is applicable in finite element analyses with recovering accurate results.
Keywords
Phase Field; Crack length; A-Star; Shortest-Path
1 INTRODUCTION
The fracture of materials has been of concern in the industry and academia over the years. In response to this issue, experimental studies have been conducted to better understand this phenomenon and develop strategies to mitigate catastrophic failures.
The introduction and subsequent widespread utilization of material parameters and their correlation, such as fracture toughness, , and crack propagation rate, (obtained from the Paris curve), have made the investigation of crack behavior simpler and the prediction of potential mechanical failure of materials. These parameters are closely linked to crack length, for which precise measurement is crucial.
Numerous crack length measurement methods are available for data obtained from experimental tests. Examples include the compliance technique (Saxena and Hudak, 1978Saxena, A. and Hudak, S. J., (1978), Review and extension of compliance information for common crack growth specimens. International Journal of Fracture, 14(5):453–468.; Kanters et al., 2015Kanters, M.J.W., Stolk, J., and Govaert, L.E., (2015). Direct comparison of the compliance method with optical tracking of fatigue crack propagation in polymers. Polymer Testing, 46:98–107.), the direct length measurement using an optical microscope (Lee et al., 2015Lee, S.Y., Huang, E.W., Woo, W., Yoon, C., Chae, H., and Yoon, S.G., (2015), Dynamic strain evolution around a crack tip under steady- and overloaded-fatigue conditions. Metals, 5(4):2109–2118.), and the direct current potential drop (Bär, 2020Bär, J., (2020) Crack detection and crack length measurement with the DC potential drop method–possibilities, challenges, and new developments. Applied Sciences (Switzerland), 10(23):1–14.). These methods rely on force and displacement measurements, which require accurate experimental equipment. However, alternative methods have been developed to address this challenge that requires no experimental output. These methods include thermoelasticity (Díaz et al., 2004Díaz, F. A., Patterson, E. A., Tomlinson R. A., and Yates J. R. (2004) Measuring stress intensity factors during fatigue crack growth using thermoelasticity. Fatigue and Fracture of Engineering Materials and Structures, 27(7):571–583.), photoelasticity (Jobin et al., 2020Jobin, T. M., Khaderi S. N., and Ramji M. (2020). Experimental evaluation of the strain intensity factor at the inclusion tip using digital photoelasticity. Optics and Lasers in Engineering, 126(June 2019):105855.), and digital image correlation (DIC) (Mokhtarishirazabad et al., 2018Mokhtarishirazabad, M., Lopez-Crespo, P., and Zanganeh, M. (2018), Stress intensity factor monitoring under cyclic loading by digital image correlation. Fatigue and Fracture of Engineering Materials and Structures, 41(10):2162–2171.).
Researchers have developed new methods and improvements for measuring crack length over the years. In a study conducted by Creel et al. (2009)Creel J.A., Stover S.M., Martin R.B., Fyhrie D.P., Hazelwood S.J., and Gibeling J.C.(2009). Compliance calibration for fracture testing of anisotropic biological materials. Journal of the Mechanical Behavior of Biomedical Materials, 2(5):571–578., compliance parameters for fracture testing of an anisotropic biological material were calibrated. Fünfschilling et al. (2010)Fünfschilling, S., Fett, T., Oberacker, R., Hoffmann, M.J., Ozcoban H., Jelitto, H., Schneider, G.A., and Kruzic, J.J. (2010) R curves from compliance and optical cracklength measurements. Journal of the American Ceramic Society, 93(9):2814–2821. compared crack length measurements obtained using optical and compliance methods when applied to R-curves. In more recent years, Yuan et al. (2021)Yuan, Y., Ge, Z., Su, X, Guo X., Suo T., Liu Y., and Yu Q. (2021) Crack length measurement using convolutional neural networks and image processing. Sensors, 21(17):1–16. utilized artificial intelligence based on convolution neural networks (CNN) and digital image processing to monitor crack length. Similarly, Farahani et al. (2022)Farahani, B.V., Direito F., Sousa P.J., Tavares P.J., Infante V., and Moreira P.P.M.G. (2022), Crack tip monitoring by multiscale optical experimental techniques. International Journal of Fatigue, 155(October 2021):106610. combined electronic speckle pattern interferometry (ESPI) with digital image correlation (DIC) to obtain the displacement field, calculate the stress intensity factor (SIF), and monitor the crack opening evolution.
Recently, phase field-based methods have been developed for modeling fracture phenomena. They use a continuous damage variable, usually denoted by with range from for undamaged material and for fully fractured material, to represent physical crack discontinuity. Many publications have reported good qualitative and quantitative agreement with experimental tests. However, measuring crack length in numerical simulations using a predetermined threshold damage remains challenging. Klinsmann et al. (2015)Klinsmann, M., Rosato, D., Kamlah, M., and McMeeking, R.M. (2015), An assessment of the phase field formulation for crack growth. Computer Methods in Applied Mechanics and Engineering, 294:313–33., Zhang et al. (2018)Zhang, P., Hu, X., Wang, X., and Yao, W. (2018), An iteration scheme for phase field model for cohesive fracture and its implementation in Abaqus. Engineering Fracture Mechanics, 204(October):268–287., Carrara et al. (2020)Carrara, P., Ambati, M., Alessi, R., and De Lorenzis, L. (2020) A framework to model the fatigue behavior of brittle materials based on a variational phase-field approach. Computer Methods in Applied Mechanics and Engineering, 361:112731., and Liu et al. (2020)Liu, J., Liang, X., Xue Y., Fu, Y., Yao, K., and Dou, F. (2020) Investigation on crack initiation and propagation in hydraulic fracturing of bedded shale by hybrid phase-field modeling. Theoretical and Applied Fracture Mechanics, 108(May):102651. used a damage threshold of and emphasize that the uncertainty in measuring the crack length contributes significantly to the discrepancy of results. This method becomes inaccurate when dealing with general shape cracks under mixed-mode loading (Zhang et al., 2017Zhang, X., Sloan, S.W., Vignes, C., and Sheng, D. (2017), A modification of the phase-field model for mixed mode crack propagation in rock-like materials. Computer Methods in Applied Mechanics and Engineering, 322:123–136.) or dynamic crack propagation cases with crack branching (Bleyer et al., 2017Bleyer J., Roux-Langlois, C., and Molinari, J.F. (2017) Dynamic crack propagation with a variational phase-field model: limiting speed, crack branching and velocity-toughening mechanisms. International Journal of Fracture, 204(1):79–100.). For accurate numerical analysis, better numerical methodologies are required to follow experimental tests and evaluate crack propagation (ASTM E647, 2016ASTM. E647-15e1 - Standard Test Method for Measurement of Fatigue Crack Growth Rates 1. ASTM Book of Standards, 03(July):1–49, 2016.) and fracture toughness (ASTM E399, 2017ASTM. E399 17 - Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIc of. Annual Book of ASTM Standards, pages 1–33, 2017.). Recently, Gehri et al. (2020)Gehri, N., Falcón J.M.-, and Kaufmann W. (2020). Automated crack detection and measurement based on digital image correlation. Construction and Building Materials, 256:119383. introduced innovative techniques that use digital image correlation (DIC) and image processing to measure cracks experimentally, resulting in improvements in multiple crack detection. Nonetheless, it is important to note that the use of image filters can substantially increase computational costs, which is dependent on the analysis periodicity and the total number of image pixels.
A-Star () is a popular shortest-path method with numerous applications, including map routing (Wang et al., 2014Wang, H., Zhou J., Zheng G., and Liang, Y. (2014). Has: Hierarchical a-star algorithm for big map navigation in special areas. In 2014 5th International Conference on Digital Home, pages 222–225.), autonomous vehicle path planning (Erke et al, 2020Erke, S., Bin D., Yiming, N., Qi, Z., Liang, X., and Dawei, Z. (2020). An improved A-Star based path planning algorithm for autonomous land vehicles. International Journal of Advanced Robotic Systems, 17(5):1–13.), movement of military troops (Dawid and Pokonieczny, 2021Dawid, W. and Pokonieczny, K. (2021). Methodology of using terrain passability maps for planning the movement of troops and navigation of unmanned ground vehicles. Sensors, 21:4682, 07.), and commercial video games (Lawande et al, 2022Lawande, S.R., Jasmine, G., Anbarasi, J., and Izhar L.I. (2022). A Systematic Review and Analysis of Intelligence-Based Pathfinding Algorithms in the Field of Video Games. Applied Sciences (Switzerland), 12(11).). This method was introduced by Peter Hart, Nils Nilsson, and Bertram Raphael (Hart et al, 1968Hart, P.E., Nilsson, N.J., and Raphael, B. (1968). A formal basis for the heuristic determination of minimum cost paths. IEEE Transactions on Systems Science and Cybernetics, 4(2):100–107.) and can be seen as an extension of Dijkstra's algorithm (Dijkstra, 1959Dijkstra, E. W. (1959). A note on two problems in connexion with graphs. Numerische Mathematik, 1(1):269–271.). Considering heuristics in leads to improved performance and efficiency in finding the shortest path with low cost and fast implementation.
This work proposes a novel application of the pathfinding algorithm to numerically measure the crack length using the finite element method (FEM) and a phase field fracture model. The key concept is to adapt the method, typically associated with graph analysis, to the finite element mesh and use the damage scalar field as a guide to track the crack path. For this purpose, the vertices, edges, and faces of graphs are associated, respectively, to the nodes, edges, and elements of the finite element meshes to make possible the adequate operation of the method.
This paper is structured as follows. Firstly, we provide a concise introduction to the phase field fracture model and the equations system to be solved. Subsequently, we present the typical pathfinding method, followed by the novel weighted cost function and the procedure for measuring the crack length. Finally, we illustrate the effectiveness of our approach by presenting the path length results for artificially fabricated damage fields and a crack resulting from fracture simulation with the phase field methodology.
2 PHASE FIELD MODELLING
Phase field methodology utilizes the continuous variable (see Figure 1) to describe an interface problem. It was initially introduced to solidification (Boettinger et al., 2002Boettinger, W J, Warren, J A, Beckermann, C, and Karma, (2002), A. phase -field simulation of solidification.) and material phase transformation (Penrose and Fife, 1990Penrose, O. and Fife, P.C. (1990). Thermodynamically consistent models of phasefield type for the kinetics of phase transitions Oliver PENROSE. 43:44–62, 1990.). More recently, this methodology has found applications in material modelling for fracture analysis, fatigue (Boldrini et al, 2016Boldrini, J. L., Barros de Moraes, E. A., Chiarelli, L. R., Fumes, F. G., and Bittencourt, M. L. (2016). A non-isothermal thermodynamically consistent phase field framework for structural damage and fatigue. Computer Methods in Applied Mechanics and Engineering, 312:395–427.; Alessi et al, 2018Alessi, R., Vidoli, S., and De Lorenzis, L. (2018) A phenomenological approach to fatigue with a variational phase-field model: The one-dimensional case. Engineering Fracture Mechanics, 190:53–73. and Haveroth et al, 2020Haveroth, G.A., Vale, M.G., Bittencourt, M.L., and Boldrini, J.L.(2020). A non-isothermal thermodynamically consistent phase field model for damage, fracture, and fatigue evolutions in elasto-plastic materials. Computer Methods in Applied Mechanics and Engineering, 364:112962.), anisotropy (Petrini et al, 2021Petrini, A.L.E.R., Boldrini J.L., and Bittencourt, M.L. (2021). A thermodynamically consistent phase field framework for anisotropic damage propagation. Latin American Journal of Solids and Structures, 18(1):1–17.), and hydrogen embrittlement (Martínez-Pañeda et al, 2018Martínez-Pañeda, E., Golahmar, A., and Niordson, C.F. (2018). A phase field formulation for hydrogen assisted cracking. Computer Methods in Applied Mechanics and Engineering, 342(March 2020):742–761.).
A phase field model with thermodynamic consistency is derived using the balance equations of continuum mechanics, resulting in a general and robust framework. In this work, the formulation is given according to Boldrini et al. (2016)Boldrini, J. L., Barros de Moraes, E. A., Chiarelli, L. R., Fumes, F. G., and Bittencourt, M. L. (2016). A non-isothermal thermodynamically consistent phase field framework for structural damage and fatigue. Computer Methods in Applied Mechanics and Engineering, 312:395–427. and Haveroth et al. (2020)Haveroth, G.A., Vale, M.G., Bittencourt, M.L., and Boldrini, J.L.(2020). A non-isothermal thermodynamically consistent phase field model for damage, fracture, and fatigue evolutions in elasto-plastic materials. Computer Methods in Applied Mechanics and Engineering, 364:112962..
The conservation of mass is given by
where is the volumetric mass density and is velocity vector field given by time derivative of the displacement field, , both represented in the Eulerian description.
The principle of virtual power (PVP) is applied to obtain the following equilibrium equation at the macro scale for an arbitrary part for the body at time :
where is the Cauchy stress tensor field, is the body load vector field and is the normal vector field to the boundary .
The PVP is also used at the microscale resulting in the following equilibrium equation:
where is the microstress vector field, is the microforce scalar field, and is the density of external body microforce.
The balance of energy, or the first principle of thermodynamics, is given by
where is the specific internal energy density, is the heat flux vector field, and is the specific heat source density.
Finally, the second law of thermodynamics, stated by the generalized Clausius-Duhem inequality, is
where is the specific entropy density, is the entropy flux due to microscopic evolution, is the absolute temperature, is the specific entropy production, and is the entropy production due to microscopic evolution.
To obtain thermodynamically consistent relations, the constitutive properties of the material are expressed in terms of the free-energy potential, , and the pseudo-potential of dissipation, . Using the balance of energy of Eq. 4, the entropy inequality, Eq. 5, is rewritten in terms of the specific free-energy by
The free energy potential depends on the internal variables as
Moreover, the terms , , and are split in their reversible (non-dissipative) and irreversible (dissipative) parts and represented as and respectively. In addition, it is assumed that and .
For the sake of brevity, some manipulations from the previous equations will be omitted here. Further details can be found in Boldrini et al. (2016)Boldrini, J. L., Barros de Moraes, E. A., Chiarelli, L. R., Fumes, F. G., and Bittencourt, M. L. (2016). A non-isothermal thermodynamically consistent phase field framework for structural damage and fatigue. Computer Methods in Applied Mechanics and Engineering, 312:395–427..
Reversible terms are chosen such that, for any admissible process, there is no entropy increase. Considering irreversible terms, and after some simplifications, we have the following equations:
To satisfy the previous inequality, we define the pseudo-potential of dissipation by
with
If , and is continuous and convex with respect to the variables and , Eq. 8 is satisfied.
Collecting the previous equations and after some simplifications, the system of governing equations is expressed by
This system of equations is very general. Changing the free energy and the pseudo-potential , different thermodynamically consistent models that account for fracture phenomena can be obtained for many materials.
Let us consider the particular case of an isothermal, isotropic, and linear elastic material with damage.
The volumetric free-energy density for a material subjected to damage is given by the sum of the elastic energy density and the energy density related to damage as
where a nearly-incompressible material is such that
The elastic strain energy density of the damaged material is obtained from the respective energy of the undamaged state by
where is the degradation function and is the elastic strain energy density of the undamaged material. The degradation function may depend on several variables. In this work, we assume that
where is a material dependent real constant.
Moreover, we split the elastic strain energy density, as proposed by Amor et al. (2009)Amor, H., Marigo, J., and Maurini, C. (2009). Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments. Journal of the Mechanics and Physics of Solids, 57(8):1209–1229., into its volumetric and deviatoric (volume-preserving) contributions such that
where
and
where and are the shear and the bulk modulus of the material, respectively. The operator is defined as .
The crack energy density is the sum of the energy density associated with the transition layers and the total bulk energy density. Therefore,
where is a parameter related to the width of the fracture layers and is the Griffith fracture energy parameter.
The non-plastic pseudo-potential of dissipation is defined by
where controls the damage change rate, see Boldrini et al. (2016)Boldrini, J. L., Barros de Moraes, E. A., Chiarelli, L. R., Fumes, F. G., and Bittencourt, M. L. (2016). A non-isothermal thermodynamically consistent phase field framework for structural damage and fatigue. Computer Methods in Applied Mechanics and Engineering, 312:395–427..
A quasi-static process is here considered for the motion equation. General thermodynamically consistent equations are expressed in the following way:
with
where is a constant that controls the damage growth rate and a small constant to avoid numerical singularity.
3 DETERMINATION OF CRACK PATH
This section describes the proposed procedure for tracking crack paths using a modified version of the method.
The based methods require prior knowledge of the starting and ending points of the path. While the crack nucleation position is generally known from experimental tests, finding the respective initial node of the mesh can be challenging in complex geometries. However, it can be estimated using the damage phase field variable and the maximum von Mises stress. Once is found, it remains constant along the procedure to determine the crack path.
Similarly, the crack tip is chosen to be located at the point with the maximum local von Mises stress. Its position must be updated at every load step of the simulation to track the crack path.
3.1 Weighted Method
The algorithm is a commonly used informed search approach in pathfinding problems. The algorithm aims to find the path from a starting point (crack nucleation site) to a goal point (crack tip) with the minimum value of the cost function for moving between and , which are both known positions.
The path is defined as a set of consecutive nodes, denoted by , in the finite element mesh, where is a positive integer. The path connects and with the minimum cost by selecting the nodes with the least accumulated cost.
The shortest path between points and are described in terms of mesh nodes denoted by the set of parent nodes which start at point . Variables with superscript denote values computed at node , while those variables with superscript represent the corresponding values for a parent node.
A cost function is typically used at each node and expressed as
Term represents the accumulated distances between the previous parent nodes and the distance between node and its parent node ; term is an estimated heuristic cost that approximates the distance between node and the final goal point (Tortora and Francese, 2003Tortora, G. and Francese, R., (2003). Graph Algorithms.).
Specifically, can be expressed as , where $ is the Euclidean norm. Similarly, .
Data for each considered node is allocated in the structure called OPEN, which stores the node number and its parent node number. The cost function is recursively calculated using Algorithm 1, while the cost function is calculated using the Euclidean distance mentioned before.
Algorithm 1: Recursive function forInitialize: | |
Input: , , , | |
Output: | |
1: | Identify the parent node of the actual node in |
2: | Calculate: |
3: | if Current node is the start node, then |
4: | Return to the main program |
5: | end if |
6: | Update: |
7: | Call this function |
The method involves tracking nodes with the minimum cost from point to in order to find the shortest path between them. At each iteration, the program allocates information for the neighboring nodes , , , for the current node also in the OPEN structure. These neighboring nodes can be chosen from the adjacent nodes (i.e., nodes connected by an element edge) or the nodes that belong to the neighboring elements (see Figure 2). When searching for adjacent nodes in a quadrilateral element, the path length will depend on both the element edges and the crack directions.
Finite element mesh with a central node (black) and its neighboring nodes connected by an edge (red) and the nodes of the neighboring elements (blue).
Node is appended to the CLOSED matrix, which contains the nodes already analyzed. The OPEN matrix must be sorted in ascending order of the cost function before the new iteration. The loop criterion is achieved when the analyzed node is the goal node. The procedure is summarized in the Algorithm 2.
The nodes in CLOSED matrix are sequenced backward for building the nodal path sequence. For example, starting from the goal node (the first position) then the parent node provides the next node in the series, and so forth, up to the start node.
Algorithm 2: WeightedInput: Mesh, damage field | |
Output: Nodal path | |
1: | Allocate the start node at |
2: | while Node Goal Node, do |
3: | for Every neighbor node by edge or by element, do |
4: | if node is not in or not in , then |
5: | Calculate , , and |
6: | Store the node information in |
7: | end if |
8: | end for |
9: | Store Node in |
10: | Remove Node From |
11: | Sort in ascending order of cost function |
12: | end while |
The previous procedure is efficient in searching for the shortest path in a finite element mesh. However, it does not guarantee that the path will pass through the crack. To enforce this condition, the cost functions and are weighted by the damage parameter, , as
where is the weight function, which reduces the cost functions where the nodal damage, , approaches , as indicated by the considered degradation function
where .
There are different weight functions, , proposed in the literature. In this case, it is crucial to ensure that the total cost decreases as the damage increases, forcing the path to pass through the most damaged nodes. Consequently, the weighted cost function no longer searches for the globally shortest path. It should be noted that the classical method is recovered when or .
3.2 Approximation of the crack length
We propose a simple and effective method for calculating the length of the crack path based on a line integral.
The parametrization of the curve in terms of the local coordinate from one node with coordinates to another node with coordinates is given by
where
In the finite element mesh, the scalar field is interpolated using the one-dimensional linear nodal shape functions expressed in the local coordinates. Eq. 27 is then rewritten for an element that contains two consecutive nodes of the path as
Since is independent of , Eq.29 is rewritten as
where is the curve length and is the nodal damage weighted by the nodal shape function.
Finally, the crack length is given by
For elements with , i.e., almost fully cracked, $ is considered .
4 RESULTS
This section presents fabricated solution examples to illustrate the capabilities and features of the proposed method. Additionally, a time-dependent case that utilizes the phase-field crack solution is presented to validate the effectiveness of the method.
4.1 Fabricated damage fields
To illustrate the influence of the weight function and validate the proposed method, the following fabricated damage field is considered:
where and are the coordinates of the nodes at a critical damage state.
We analyze the two crack configurations of Figure 3 to test the effectiveness of the proposed method. The first configuration, (Figure 3a), represents a curved crack configuration, while the second configuration (Figure 3b), depicts a branched crack configuration. Circular shapes were chosen for the crack configurations motivated by the fact that approximating circular geometry using straight lines is particularly challenging. Consequently, obtaining good results for these cases ensures that simpler cases can also be accurately approximated.
Two-dimensional domains: (a) and (b) , illustrating the starting points (blue) and goal points (red) associated with each domain.
The critical damage nodes coordinates are defined as
where is the radius of the circumference and is a small increment.
The starting and goal points coordinates () are, respectively, and for . The geometry represents a straight crack followed by branching and described by two paths. The starting point coordinate is and the two goal points coordinates are and .The mesh domain is defined in the interval for both and coordinates and discretized using equally spaced elements in each direction.
The resulting paths for two different choices of are presented in Figure 4 and 5 using Algorithms 1 and 2. Firstly, the usual method is recovered for and the shortest path using Eq. 24 is obtained and illustrated in Figs. 4b and 5b. For the weighted method using Eq. 26, the path pass through the nodes with the largest values of as shown in Figs. 4c and 5c.
(a) Scalar field for the curve and the obtained paths associated to the weigh function where (b) and (c) .
(a) Scalar field for the curve and the obtained paths associated to the weigh function where (b) and (c) .
The proposed method showed qualitative agreement compared to the original shape. However, a quantitative analysis must be considered to validate the methodology. For this purpose, the method presented in Section 3.2 is applied, and the crack length is obtained for different mesh refinements () and different neighbor search options (by edge or by element).
The results are summarized in Table 1 and the accuracy depends on the path shape and on the choice of searching for neighbor nodes (by edge or by element). Evidently, a circular shape in a mesh with squared elements will present lower accuracy when compared to a straight line. In this case, searching for neighbors by node provides more possible nodes and directions improving the path length result.
4.2 Crack under propagation
The proposed method has a practical application in tracking the path of a stable crack over time. Specifically, the crack path is measured at regular intervals during the simulation, between the initial time and final time . This technique is commonly used in fatigue crack propagation tests to determine the propagation velocity of a stable crack under a given stress intensity factor. The resulting data is represented by the Paris-Erdogan curve, which provides valuable material properties, including the threshold stress intensity factor and stable propagation parameters.
The considered square domain with side is discretized into structured meshes with elements. Again, the damage scalar field is given by Eq. 32 and updated every time step. The simulation time ranges from to with time step .
A straight crack pattern is here considered. The nodal positions are constant for every time step and the critical damaged nodal coordinates are given by the curve
with
where is the time-dependent crack tip position function, is the starting crack tip coordinate in the -axis and is the width of the domain in -direction. The starting crack parameters are , , .
The crack length is measured using three meshes with different refinements () and compared with the exact value given by Eq. 35 and shown in Figure 6.
The proposed method can closely follow the damage field and the crack tip, even with a coarser mesh. Evidently, the accuracy of the result depends on the mesh refinement. For the coarsest mesh, Figs. 6a-6b shows that the crack length evolution behaves as a step function where the damage rate is lower compared to Fig. 6c. Therefore, in the crack propagation simulation, the result accuracy is dependent on the measurement frequency (), the crack speed (), and the mesh refinement () in such a way that .
The average percentage error is calculated for the three mesh refinements with values of , and , respectively.
To avoid dependence on the sampling frequency, the measurements are counted only when the crack advances to a new node. Figure 7 shows the results for . The differences are more pronounced in the beginning but more accurate when measured. The average percentage error for this case is .
4.3 Comparison with the compliance method
The resulting crack length obtained from the method is compared with the crack length using the compliance technique. The compliance technique is often used in experimental tests since it is simpler to obtain the crack length when compared with other more accurate methods. The length of a straight crack can be indirectly measured using the crack mouth opening displacement (CMOD) and the applied force. For a standard specimen geometry, the crack length results from a calibrated function of the CMOD, the applied force, , and the effective stiffness, . Also, the geometric parameters and are, respectively, the width (measured from the center of the fixation holes) and the thickness. The C(T) specimen is used with dimensions shown in Figure 8 according to the ASTM E399 standard (ASTM E399, 2017ASTM. E399 17 - Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIc of. Annual Book of ASTM Standards, pages 1–33, 2017.).
A fracture toughness test (mode I opening test) in the specimen under displacement control is simulated using the phase-field model stated in Eq. 21. The CMOD and the reaction force are provided from the finite element analysis, which are inputs for the crack length functional given by
with
where is the compliance crack length, is a dimensionless variable, and is the secant compliance given by
Eq. 37 is valid for crack lengths of (ASTM E399, 2017ASTM. E399 17 - Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIc of. Annual Book of ASTM Standards, pages 1–33, 2017.).
For the phase-field simulations, the parameters are , , , , and . The analysis is done under the elastic regime, without plastic deformations. In this study, we have discretized the geometry utilizing over 75000 linear quadrilateral elements which have been strategically concentrated along the expected crack path.
The crack length is computed using the proposed method every 20 time steps and also with the compliance technique every time step with . The mesh is unstructured and refined in the expected crack path with an average element width of . The goal node is located at the maximum von Mises stress and a physical crack is summed when .
The crack length measurements for the proposed and the compliance methods are shown in Figure 9. Notably, the weighted approach demonstrates remarkable agreement with the compliance method across the entire crack propagation.
5 CONCLUSIONS
Accurate measurement of crack paths is crucial for validating numerical simulations and improving the reliability of results. In this work, we introduced a novel application of the method for tracking cracks in numerical simulations. Unlike conventional methods that only measure the direct path from the crack nucleation site to the crack tip, our approach can track the entire crack path, regardless of its shape. The method is efficient because it searches the nodes closer to the crack path.
We illustrated the effectiveness of the proposed method by presenting fabricated solution cases for both curved and branched crack patterns. Our method was compared to the traditional approach, and the results showed excellent qualitative and quantitative agreement. We also presented a time-dependent fabricated case that highlighted the influence of the mesh, crack velocity, and sampling frequency of the results.
Furthermore, we applied our method to crack propagation simulations using the phase-field methodology and obtained an excellent agreement with the compliance technique. This improved accuracy allows for more reliable results and refined parameters.
In future studies, we plan to extend the proposed method to the experimental field using digital image correlation (DIC) methodology, which bridges experimental and numerical analyses.
Acknowledgements:
The authors would like to acknowledge CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior) for the scholarship under grant 88882.435204/2019-01. J.A. Avila is a Serra Hunter Fellow and a CNPq fellow.
References
- Alessi, R., Vidoli, S., and De Lorenzis, L. (2018) A phenomenological approach to fatigue with a variational phase-field model: The one-dimensional case. Engineering Fracture Mechanics, 190:53–73.
- Amor, H., Marigo, J., and Maurini, C. (2009). Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments. Journal of the Mechanics and Physics of Solids, 57(8):1209–1229.
- ASTM. E399 17 - Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIc of. Annual Book of ASTM Standards, pages 1–33, 2017.
- ASTM. E647-15e1 - Standard Test Method for Measurement of Fatigue Crack Growth Rates 1. ASTM Book of Standards, 03(July):1–49, 2016.
- Bär, J., (2020) Crack detection and crack length measurement with the DC potential drop method–possibilities, challenges, and new developments. Applied Sciences (Switzerland), 10(23):1–14.
- Bleyer J., Roux-Langlois, C., and Molinari, J.F. (2017) Dynamic crack propagation with a variational phase-field model: limiting speed, crack branching and velocity-toughening mechanisms. International Journal of Fracture, 204(1):79–100.
- Boettinger, W J, Warren, J A, Beckermann, C, and Karma, (2002), A. phase -field simulation of solidification.
- Boldrini, J. L., Barros de Moraes, E. A., Chiarelli, L. R., Fumes, F. G., and Bittencourt, M. L. (2016). A non-isothermal thermodynamically consistent phase field framework for structural damage and fatigue. Computer Methods in Applied Mechanics and Engineering, 312:395–427.
- Carrara, P., Ambati, M., Alessi, R., and De Lorenzis, L. (2020) A framework to model the fatigue behavior of brittle materials based on a variational phase-field approach. Computer Methods in Applied Mechanics and Engineering, 361:112731.
- Creel J.A., Stover S.M., Martin R.B., Fyhrie D.P., Hazelwood S.J., and Gibeling J.C.(2009). Compliance calibration for fracture testing of anisotropic biological materials. Journal of the Mechanical Behavior of Biomedical Materials, 2(5):571–578.
- Dawid, W. and Pokonieczny, K. (2021). Methodology of using terrain passability maps for planning the movement of troops and navigation of unmanned ground vehicles. Sensors, 21:4682, 07.
- Díaz, F. A., Patterson, E. A., Tomlinson R. A., and Yates J. R. (2004) Measuring stress intensity factors during fatigue crack growth using thermoelasticity. Fatigue and Fracture of Engineering Materials and Structures, 27(7):571–583.
- Dijkstra, E. W. (1959). A note on two problems in connexion with graphs. Numerische Mathematik, 1(1):269–271.
- Erke, S., Bin D., Yiming, N., Qi, Z., Liang, X., and Dawei, Z. (2020). An improved A-Star based path planning algorithm for autonomous land vehicles. International Journal of Advanced Robotic Systems, 17(5):1–13.
- Farahani, B.V., Direito F., Sousa P.J., Tavares P.J., Infante V., and Moreira P.P.M.G. (2022), Crack tip monitoring by multiscale optical experimental techniques. International Journal of Fatigue, 155(October 2021):106610.
- Fünfschilling, S., Fett, T., Oberacker, R., Hoffmann, M.J., Ozcoban H., Jelitto, H., Schneider, G.A., and Kruzic, J.J. (2010) R curves from compliance and optical cracklength measurements. Journal of the American Ceramic Society, 93(9):2814–2821.
- Gehri, N., Falcón J.M.-, and Kaufmann W. (2020). Automated crack detection and measurement based on digital image correlation. Construction and Building Materials, 256:119383.
- Hart, P.E., Nilsson, N.J., and Raphael, B. (1968). A formal basis for the heuristic determination of minimum cost paths. IEEE Transactions on Systems Science and Cybernetics, 4(2):100–107.
- Haveroth, G.A., Vale, M.G., Bittencourt, M.L., and Boldrini, J.L.(2020). A non-isothermal thermodynamically consistent phase field model for damage, fracture, and fatigue evolutions in elasto-plastic materials. Computer Methods in Applied Mechanics and Engineering, 364:112962.
- Jobin, T. M., Khaderi S. N., and Ramji M. (2020). Experimental evaluation of the strain intensity factor at the inclusion tip using digital photoelasticity. Optics and Lasers in Engineering, 126(June 2019):105855.
- Kanters, M.J.W., Stolk, J., and Govaert, L.E., (2015). Direct comparison of the compliance method with optical tracking of fatigue crack propagation in polymers. Polymer Testing, 46:98–107.
- Klinsmann, M., Rosato, D., Kamlah, M., and McMeeking, R.M. (2015), An assessment of the phase field formulation for crack growth. Computer Methods in Applied Mechanics and Engineering, 294:313–33.
- Lawande, S.R., Jasmine, G., Anbarasi, J., and Izhar L.I. (2022). A Systematic Review and Analysis of Intelligence-Based Pathfinding Algorithms in the Field of Video Games. Applied Sciences (Switzerland), 12(11).
- Lee, S.Y., Huang, E.W., Woo, W., Yoon, C., Chae, H., and Yoon, S.G., (2015), Dynamic strain evolution around a crack tip under steady- and overloaded-fatigue conditions. Metals, 5(4):2109–2118.
- Liu, J., Liang, X., Xue Y., Fu, Y., Yao, K., and Dou, F. (2020) Investigation on crack initiation and propagation in hydraulic fracturing of bedded shale by hybrid phase-field modeling. Theoretical and Applied Fracture Mechanics, 108(May):102651.
- Martínez-Pañeda, E., Golahmar, A., and Niordson, C.F. (2018). A phase field formulation for hydrogen assisted cracking. Computer Methods in Applied Mechanics and Engineering, 342(March 2020):742–761.
- Mokhtarishirazabad, M., Lopez-Crespo, P., and Zanganeh, M. (2018), Stress intensity factor monitoring under cyclic loading by digital image correlation. Fatigue and Fracture of Engineering Materials and Structures, 41(10):2162–2171.
- Penrose, O. and Fife, P.C. (1990). Thermodynamically consistent models of phasefield type for the kinetics of phase transitions Oliver PENROSE. 43:44–62, 1990.
- Petrini, A.L.E.R., Boldrini J.L., and Bittencourt, M.L. (2021). A thermodynamically consistent phase field framework for anisotropic damage propagation. Latin American Journal of Solids and Structures, 18(1):1–17.
- Saxena, A. and Hudak, S. J., (1978), Review and extension of compliance information for common crack growth specimens. International Journal of Fracture, 14(5):453–468.
- Tortora, G. and Francese, R., (2003). Graph Algorithms.
- Wang, H., Zhou J., Zheng G., and Liang, Y. (2014). Has: Hierarchical a-star algorithm for big map navigation in special areas. In 2014 5th International Conference on Digital Home, pages 222–225.
- Yuan, Y., Ge, Z., Su, X, Guo X., Suo T., Liu Y., and Yu Q. (2021) Crack length measurement using convolutional neural networks and image processing. Sensors, 21(17):1–16.
- Zhang, P., Hu, X., Wang, X., and Yao, W. (2018), An iteration scheme for phase field model for cohesive fracture and its implementation in Abaqus. Engineering Fracture Mechanics, 204(October):268–287.
- Zhang, X., Sloan, S.W., Vignes, C., and Sheng, D. (2017), A modification of the phase-field model for mixed mode crack propagation in rock-like materials. Computer Methods in Applied Mechanics and Engineering, 322:123–136.
Edited by
Publication Dates
-
Publication in this collection
17 July 2023 -
Date of issue
2023
History
-
Received
15 Mar 2023 -
Reviewed
24 Mar 2023 -
Accepted
30 May 2023 -
Published
05 June 2023