Temperature Effects on Critical Energy Release Rate for Aluminum and Titanium Alloys

: This work investigates temperature’s effect on the critical energy release rate using damage mechanics material models and the element deletion method. The energy release rate describes the decrease in total potential energy per increase in crack surface area. The critical energy release rate is widely used as the failure criterion for various elastic and plastic materials. In real-life scenarios, fractures may occur at different temperatures. The temperature dependency of the critical energy release rate for aluminum 2024-T351 and titanium Ti-6Al-4V is studied in this work. We utilized test-data-based advanced material models of these two alloys, considering the strain rate, temperature, and state of stress for plasticity and failure. These material models are used to simulate a three-dimensional fracture specimen to find the critical energy release rate at different temperatures. A new method to calculate the critical energy release rate with the element deletion method is introduced and verified with the virtual crack opening method. This method enables the calculation of the energy release rate in a classical damage mechanics simulation for dynamic cack propagation. The simulation result indicates that the critical energy release rate increases with rising temperatures for these alloys.


Introduction
Fracture mechanics is an important topic in engineering sciences, especially for automotive and aerospace applications [1][2][3][4].The history of fracture mechanics dates back 100 years.Inglis first proposed the analytical solution of a stress field around a crack tip [5].However, his derivation was based on small-strain theory, which leads to an unrealistic stress singularity at the crack tip [6].Irwin defines the stress intensity factor [7], but the concept is limited to small-strain, linear elastic material and brittle fractures.James Rice [8] cleverly avoids the crack tip singularity by defining a path-independent J-integral around the crack tip.However, this idea is not easily extended to three-dimensional fractures, mixed-mode fractures, or fractures with large-scale plasticity, which is a common case for metallic materials [9,10].The energy release rate concept was first introduced by Griffith based on a glass experiment [11].Griffith's energy release rate states that the dissipated energy is converted to an increment in the crack surface energy; however, this statement is not totally supported by experiments.For example, the surface energy of glass measured from the contact angle is 0.0834 J/m 2 [12].The critical energy release rate of glass is around 7 J/m 2 [13].The modern energy release rate is defined as the amount of total potential energy decrease per increase in the fracture surface area [14].In principle, the energy release rate concept can be extended to large-strain theory and treating materials with non-linearity in elastic and plastic deformation.The energy release rate is more general, as it hints at energy conservation.In practice, the energy release rate value can be calculated using either numerical or experimental methods.There are several ways to calculate energy release rates.
The critical energy release rate can be calculated from the force-displacement curve measured from a fracture specimen during a loading and unloading process [14][15][16].This method is also applied in experiments for composite [17] and biological tissue material [18][19][20][21], etc.The energy release rate can also be measured with the virtual crack opening method and J-integral.Notice that these three methods can be applied in both experiments and simulations.The node release method can be used to calculate the energy release rate in finite element simulation, as demonstrated by Wang et al. [6,9].In this method, the crack opening is represented by the separation of two previously connected nodes at the crack tip.The connection force proportionally drops to zero while the displacement is measured.The force-displacement data are used to calculate the energy release rate.
In this study, a method is developed to calculate the energy release rate with the element deletion method (element erosion method), and its verification is presented.Notice that the element deletion method has already been widely adopted to model failure in finite element simulation [22][23][24][25].Examples include biological material [26], metal manufacturing processes [27,28], civil engineering applications [29], etc.However, the element deletion method has not been used to calculate the energy release rate.We first verify this method with the virtual crack opening method [14,15] to demonstrate its feasibility.Then, we use this method to simulate fracture tests of realistic materials and calculate the critical energy release rate.Based on this, we found the relation between critical energy release and temperature.
The primary motivation of this research is to bridge the gap between fracture mechanics and damage mechanics.Notice that element deletion is the default method for modeling failure in damage mechanics simulations with a finite element analysis.In principle, we should be able to use accurate damage mechanics material models to study fracture behavior.Damage mechanics describe material deformation and failure in terms of stress-strain relations and the strain at failure.Fracture mechanics describe material failure with material fracture toughness parameters.For example, the critical energy release rate, critical J-integral, critical stress intensity factor, critical Crack Tip Opening Displacement (CTOD), and critical Crack Tip Opening Angle (CTOA) are used as material model parameters in fracture mechanics.One material cannot have two independent failure criteria.Theoretically, fracture mechanics material parameters can be computed from the damage mechanics material model using finite element simulation.
Realistic damage mechanics material models were used in the simulation.Plasticity and failure material models of a 0.5-inch aluminum 2024-T351 plate and a 0.5-inch titanium-6Al-4V plate were previously developed by the Federal Aviation Administration to study aircraft engine containment [1,2,30,31].In these material models, the plasticity depends on temperature and the strain rate; the failure depends on the strain rate, temperature, and state of the stress.These models are based on realistic test data, and they represent realistic material behavior under different temperatures, strain rates, and states of stress.We utilized these material models in the simulation for the calculation of the critical energy release rate.

Theroy of Energy Release Rate
The energy release rate is defined as the amount of energy dissipated during crack propagation to increase the crack length.The energy release rate in two dimensions is defined in Equation (1) [9,14], Symmetry 2024, 16, 142 where a is the incremented crack length, U represents the strain energy in Equation ( 2), and V is the potential of the external force in Equation (3).
In Equation ( 2), W stands for the strain energy per unit area, and A 0 is the area of a body.In Equation (3), T i describes the external traction force on boundary Γ t , and u i measures the displacement of the boundary Γ t .Classical fracture mechanics view the energy release rate as a variable that can be calculated under any loading condition.The crack starts to extend when the energy release rate exceeds the critical energy release rate described in Equation (4).
Classical fracture mechanics treat the critical energy release rate as a constant in material model parameters.In this study, we found that the critical energy release rate depends on temperature.

Simulation Material Model
The MAT 224 material subroutine in LS-DYNA is a fully tabulated elastic-thermalvisco-plastic deformation and failure model.The linear elasticity part uses Young's modulus and Poisson's ratio as material model parameters.The thermal-visco-plasticity part features a fully tabulated input.The strain rate sensitivity table, LCK1, contains the stress vs. plastic strain curves of different strain rates.The temperature sensitivity table, LCKT, contains the stress vs. plastic strain curves for different temperatures.Each curve in the table is defined by discrete data points, and they represent the stress-strain curve under certain constant strain rates and constant temperatures.During the simulation process, the material subroutine automatically interpolates between input data points to obtain the stress value under any strain, strain rate, and temperature at the timestep.The failure criterion is primarily based on plastic strain.A 3D tabulated failure surface, LCF, defines the plastic failure strain as a function of stress triaxiality and Lode parameters.Similar to the plasticity model, the failure subroutine interpolates between the input data points to obtain the failure strain at any strain rate, temperature, and state of the stress.The tabulated curve LCG scales the plastic failure strain with respect to different plastic strain rates.The tabulated curve LCH relates plastic failure strain to temperature.The tabulated curve LCI compensates for the element size effect, making coarser mesh and finer mesh fail under the same conditions.The MAT 224 material subroutine has semi-thermal-mechanical coupling built in.The heat generation is calculated inside the material subroutine and influences the material behavior, but the heat cannot be dissipated if only the mechanical solver is used.The parameter settings for the material model are based on fitting the test data of a 0.5-inch aluminum 2024-T351 plate and a 0.5-inch Ti-6Al-4V plate with an inverse analysis [1,31].In a quasistatic test, the heat is completely dissipated into the environment, so the temperature is almost constant.To simulate such an isothermal condition, the Taylor-Quinney Coefficient BETA, or the ratio to convert the dissipated energy into heat, needs to be set to a very small number.The unit system for the material model is mm, second, ton, Newton, and Kelvin.

Critical Energy Release Rate Calculated from Element Deletion Method
We introduced a new method to calculate the critical energy release rate with the element deletion method.The method is later verified against the virtual crack opening method [14,15].The material model of aluminum 2024-T3 is used in the verification simulation for both methods.In the simulation, the rectangular sample is 24 mm in height, 12.5 mm in width, and 1 mm in thickness and carries a line crack.The simulation model for the element deletion method is shown in Figure 1.The element size is uniformly 1 mm.The velocity-controlled boundary condition is adopted.The strain rate of the simulation is roughly around 8.33(1/s).
simulation for both methods.In the simulation, the rectangular sample is 24 mm in heig 12.5 mm in width, and 1 mm in thickness and carries a line crack.The simulation mo for the element deletion method is shown in Figure 1.The element size is uniformly 1 m The velocity-controlled boundary condition is adopted.The strain rate of the simulat is roughly around 8.33 (1/s).

Critical Energy Release Rate from Element Deletion Method
In the element deletion method, an element is deleted if it satisfies the failure crite of the material model.The deleted element will be set to zero stiffness and will be exclud from the calculation in future time steps.The strain energy of the deleted element, ∆ can be directly calculated by the finite element solver.For example, the crack tip elem (brown) is the target element used to calculate the energy release in Figure 1b,c.In DYNA, the strain energy of the deleted element can be output by assigning each elem a separate part ID and outputting the internal energy of the deleted part.The crit energy release rate is calculated by Equation ( 5),

𝐺 = ∆𝑈 2𝑧∆𝑎
where ∆ is the maximum internal energy of the element right before the element fa ∆ is the incremented crack length due to the element failing, and  is the thickness the specimen.This equation is another version that states the energy dissipated for increasing crack area.

Element Deletion Method Result
In the simulation, an element is deleted if it meets the failure criteria, as shown Figure 1.The first failed element, in theory, is behind the crack tip, which is the dar brown element shown in Figure 1c.Because we used the classical element deletion meth without any special treatment of the crack tip element, the "sharp" crack tip will fi become "blunt" by deleting the yellow element.At this stage, the sharp crack tip in Fig 1b becomes the blunt crack tip in Figure 1c.Further crack extensions will always hav blunt crack tip.This is the limitation of the element deletion method; the crack geome cannot be precise.This will be less of a problem if an ultra-fine element size is chosen.
Here, the yellow element is to the left of the original crack tip, while the bro element is our target element used for the energy release rate calculation.The crack

Critical Energy Release Rate from Element Deletion Method
In the element deletion method, an element is deleted if it satisfies the failure criteria of the material model.The deleted element will be set to zero stiffness and will be excluded from the calculation in future time steps.The strain energy of the deleted element, ∆U, can be directly calculated by the finite element solver.For example, the crack tip element (brown) is the target element used to calculate the energy release in Figure 1b,c.In LS-DYNA, the strain energy of the deleted element can be output by assigning each element a separate part ID and outputting the internal energy of the deleted part.The critical energy release rate is calculated by Equation ( 5), where ∆U is the maximum internal energy of the element right before the element fails, ∆a is the incremented crack length due to the element failing, and z is the thickness of the specimen.This equation is another version that states the energy dissipated for an increasing crack area.

Element Deletion Method Result
In the simulation, an element is deleted if it meets the failure criteria, as shown in Figure 1.The first failed element, in theory, is behind the crack tip, which is the darker brown element shown in Figure 1c.Because we used the classical element deletion method without any special treatment of the crack tip element, the "sharp" crack tip will first become "blunt" by deleting the yellow element.At this stage, the sharp crack tip in Figure 1b becomes the blunt crack tip in Figure 1c.Further crack extensions will always have a blunt crack tip.This is the limitation of the element deletion method; the crack geometry cannot be precise.This will be less of a problem if an ultra-fine element size is chosen.
Here, the yellow element is to the left of the original crack tip, while the brown element is our target element used for the energy release rate calculation.The crack tip element (brown) has a maximum internal energy of 61.0690N•mm in Figure 2. It is considered that the strain energy decreased with an increase in the crack surface area.The new crack surface area caused by the element deletion is 2 mm 2 .The energy release rate is calculated as 30.5345N/mm by Equation (5).Notice that the symmetry in the geometry and boundary condition makes the kinetic energy approximately zero, as shown in Figure 2.
element (brown) has a maximum internal energy of 61.0690 N • mm in Figure 2. It is considered that the strain energy decreased with an increase in the crack surface area.The new crack surface area caused by the element deletion is 2 mm .The energy release rate is calculated as 30.5345N/mm by Equation (5).Notice that the symmetry in the geometry and boundary condition makes the kinetic energy approximately zero, as shown in Figure 2.

Critical Energy Release Rate from Virtual Crack Opening Method
The virtual crack opening method is a conventional method for finding the critical energy release rate [14][15][16].It includes displacement control and force control methods.The displacement control method is applied here as it is closer to the physical test condition.The simulation model for the virtual crack opening method is shown in Figure 3. Notice that the simulation model for the element deletion method is shown in Figure 1.An identical mesh is used for the simulation of both methods.The element size is a uniform 1 mm .First, the specimen with an initial crack length a is pulled apart right before the crack starts to extend, as shown in Figure 3a.The simulation results of the crosssectional force ( ) vs. the grip displacement ( ) are recorded.Second, the geometry of the finite element model is modified to have an initial crack size of  + ∆.The modified specimen is pulled apart right before the crack starts to extend, as shown in Figure 3b.The simulation results of the cross-sectional force ( ) vs. the grip displacement ( ) are recorded.The critical energy release rate can be calculated with Equation (6).

Critical Energy Release Rate from Virtual Crack Opening Method
The virtual crack opening method is a conventional method for finding the critical energy release rate [14][15][16].It includes displacement control and force control methods.The displacement control method is applied here as it is closer to the physical test condition.The simulation model for the virtual crack opening method is shown in Figure 3. Notice that the simulation model for the element deletion method is shown in Figure 1.An identical mesh is used for the simulation of both methods.The element size is a uniform 1 mm.First, the specimen with an initial crack length a is pulled apart right before the crack starts to extend, as shown in Figure 3a.The simulation results of the cross-sectional force (F 1 ) vs. the grip displacement (S 1 ) are recorded.Second, the geometry of the finite element model is modified to have an initial crack size of a + ∆a.The modified specimen is pulled apart right before the crack starts to extend, as shown in Figure 3b.The simulation results of the cross-sectional force (F 2 ) vs. the grip displacement (S 2 ) are recorded.The critical energy release rate can be calculated with Equation (6).
ry 2024, 16, x FOR PEER REVIEW 5 of 14 element (brown) has a maximum internal energy of 61.0690 N • mm in Figure 2. It is considered that the strain energy decreased with an increase in the crack surface area.The new crack surface area caused by the element deletion is 2 mm .The energy release rate is calculated as 30.5345N/mm by Equation (5).Notice that the symmetry in the geometry and boundary condition makes the kinetic energy approximately zero, as shown in Figure 2.

Critical Energy Release Rate from Virtual Crack Opening Method
The virtual crack opening method is a conventional method for finding the critical energy release rate [14][15][16].It includes displacement control and force control methods.The displacement control method is applied here as it is closer to the physical test condition.The simulation model for the virtual crack opening method is shown in Figure 3. Notice that the simulation model for the element deletion method is shown in Figure 1.An identical mesh is used for the simulation of both methods.The element size is a uniform 1 mm .First, the specimen with an initial crack length a is pulled apart right before the crack starts to extend, as shown in Figure 3a.The simulation results of the crosssectional force ( ) vs. the grip displacement ( ) are recorded.Second, the geometry of the finite element model is modified to have an initial crack size of  + ∆.The modified specimen is pulled apart right before the crack starts to extend, as shown in Figure 3b.The simulation results of the cross-sectional force ( ) vs. the grip displacement ( ) are recorded.The critical energy release rate can be calculated with Equation (6).The force vs. displacement curve for the model with the original and incremented crack length is shown in Figure 4.The integration of the force-displacement curve gives the potential energy.The decrease in the potential energy can thus be calculated.The numerical integration is carried out by using the trapezoidal rule in MATLAB, and the value is 62.2991 N•mm.Since the incremented area for the crack surface is 2 mm 2 , the critical energy release rate ends up being 31.1496N/mm.The crack tip is at the center of the yellow region.This is a side view of a three-dimensional specimen.The crack tip front is perpendicular to the two-dimensional surface and goes inside the surface.

Virtual Crack Opening Method Result
The force vs. displacement curve for the model with the original and incremented crack length is shown in Figure 4.The integration of the force-displacement curve gives the potential energy.The decrease in the potential energy can thus be calculated.The numerical integration is carried out by using the trapezoidal rule in MATLAB, and the value is 62.2991 N • mm.Since the incremented area for the crack surface is 2 mm , the critical energy release rate ends up being 31.1496N/mm.

Critical Energy Release Rate Comparison and Verification
The critical energy release rate values calculated by both methods are almost the same, as summarized in Table 1.This result verified that the element deletion method can be used to calculate the critical energy release rate.

Method
Virtual Crack Opening Element Deletion  (N/mm) 31.149630.5345

Critical Energy Release Rate Comparison of Verification
The previous section, Section 2.3, demonstrated the feasibility of using the element deletion method to calculate the critical energy release rate.In this section, we will implement this method for dynamic crack propagation.The 3D thin plate still has dimensions of 24 mm in length, 12.5 mm in width, and 1 mm in thickness.This simulation model has a line crack of 2.5 mm as shown in Figure 5.A line crack is modeled by not letting the adjacent elements share nodes.

Critical Energy Release Rate Comparison and Verification
The critical energy release rate values calculated by both methods are almost the same, as summarized in Table 1.This result verified that the element deletion method can be used to calculate the critical energy release rate.

Critical Energy Release Rate Comparison of Verification
The previous section, Section 2.3, demonstrated the feasibility of using the element deletion method to calculate the critical energy release rate.In this section, we will implement this method for dynamic crack propagation.The 3D thin plate still has dimensions of 24 mm in length, 12.5 mm in width, and 1 mm in thickness.This simulation model has a line crack of 2.5 mm as shown in Figure 5.A line crack is modeled by not letting the adjacent elements share nodes.
A fully integrated solid mesh with an element of 0.125 mm is used in the simulation.A smaller element size is used for better accuracy.The aspect ratio of the elements is close to one, representing a perfect cube, as shown in Figure 5.The geometry is the same as for the mesh used in Section 2.3, Figures 1 and 3, except for the crack length and mesh size.
The velocity-specified boundary condition is applied to the top and bottom nodes to pull the specimen apart.The velocity on each end has the same value but in opposite directions.The loading velocity is 100 mm/s, which roughly represents a strain rate of 8.33(1/s).The Taylor-Quinney Coefficient used to convert the dissipated plastic energy into heat is set to 1 × 10 6 to make the temperature constant in the simulation.We chose four constant temperature cases to study the temperature effect: 186 K, 300 K, 420 K, and 537 K.
At the crack extension moment, the maximum internal energy of the first deleted element is recorded and used as the dissipated energy, ∆U.The newly generated crack surface area is recorded as ∆A.The critical energy release rate is calculated as G = ∆U /∆A for all elements on the crack extension path.It is worthwhile to mention that the kinetic energy in the center is almost zero due to the symmetry in the geometry and boundary conditions.A fully integrated solid mesh with an element of 0.125 mm is used in the simulation.A smaller element size is used for better accuracy.The aspect ratio of the elements is close to one, representing a perfect cube, as shown in Figure 5.The geometry is the same as for the mesh used in Section 2.3, Figures 1 and 3, except for the crack length and mesh size.
The velocity-specified boundary condition is applied to the top and bottom nodes to pull the specimen apart.The velocity on each end has the same value but in opposite directions.The loading velocity is 100 mm/s, which roughly represents a strain rate of 8.33 (1/s).The Taylor-Quinney Coefficient used to convert the dissipated plastic energy into heat is set to 1 × 10 to make the temperature constant in the simulation.We chose four constant temperature cases to study the temperature effect: 186 K, 300 K, 420 K, and 537 K.
At the crack extension moment, the maximum internal energy of the first deleted element is recorded and used as the dissipated energy, ∆.The newly generated crack surface area is recorded as ∆ .The critical energy release rate is calculated as  = ∆ /∆ for all elements on the crack extension path.It is worthwhile to mention that the kinetic energy in the center is almost zero due to the symmetry in the geometry and boundary conditions.
The specimen is pulled apart until the center elements fail and are deleted one by one until the crack goes through the entire specimen.Each element on the crack extension path is assigned an individual part ID in LS-DYNA, so the part energy file MATSUM can be used to export the strain energy time history of the element, as shown in Figure 6.Since this is a symmetry problem, we chose to output the energy for the element to the bottom of the symmetry plane, as shown in the section view in Figure 6 The 300 K temperature case is chosen to study: 1.The critical energy release rate vs. the crack length of a dynamic crack.2. The crack speed vs. the incremented crack length of a dynamic crack.The specimen is pulled apart until the center elements fail and are deleted one by one until the crack goes through the entire specimen.Each element on the crack extension path is assigned an individual part ID in LS-DYNA, so the part energy file MATSUM can be used to export the strain energy time history of the element, as shown in Figure 6.Since this is a symmetry problem, we chose to output the energy for the element to the bottom of the symmetry plane, as shown in the section view in Figure 6 The 300 K temperature case is chosen to study: 1.The critical energy release rate vs. the crack length of a dynamic crack.2. The crack speed vs. the incremented crack length of a dynamic crack.

Temperature Dependency of Critical Energy Release Rate
Table 2 records the average critical energy release rate under different temperatures for aluminum 2024-T351 and titanium-6Al-4V.It is seen that the critical energy release

Temperature Dependency of Critical Energy Release Rate
Table 2 records the average critical energy release rate under different temperatures for aluminum 2024-T351 and titanium-6Al-4V.It is seen that the critical energy release rate increases with an increase in temperature.This observation is in accordance with the common understanding that metal is easier to fracture at low temperatures than at high temperatures.For aluminum 2024-T351, the critical energy release rate has no significant difference between 186 K and 420 K.The martial critical energy release rate increases dramatically from 420 K to 537 K.For titanium-6Al-4V, the critical energy release rate increases almost linearly with respect to temperature, as shown in Table 2.

Dynamic Crack Size Dependency of Critical Energy Release Rate
In classical fracture mechanics theory, the critical energy release rate is treated as a constant of material model parameters, regardless of the fracture length in fracture propagation.To verify this theory, the critical energy release rate vs. incremented crack size is shown in Figure 7. Notice that the relationship between the crack size and critical energy release rate comes from one simulation of a dynamically extending crack, where the crack size naturally grows from small to large.For aluminum 2024-T351 and titanium-6Al-4V at a temperature of 300 K, the simulation result shows that the critical energy release rate is approximately constant regardless of the fracture extension length.The critical energy release rate calculated from the elements on the surface of the specimen is higher than the value calculated from the center elements.A comparison is shown in Table 2.A similar pattern is also seen for the critical energy release rate at other temperatures as well.This may be caused by the center element and the edge element failing in different states of stress (triaxiality and lode parameter).In the material model we used, different states of stress correspond to different failure strains, and thus, different strain energies are deleted at the moment of failure.

Crack Speed in Dynamic Crack Popogation
The crack speed vs. incremental crack size is shown in Figure 8.The simulation is for aluminum 2024-T351 and titanium-6Al-4V at a temperature of 300 K.The crack speed of aluminum 2024-T351 increases from below 1000 mm/s to around 2000 mm/s .A large fluctuation is seen when the material is about to fail.This is similar to the prediction in classical fracture mechanics.The titanium-6Al-4V simulation shows an unstable crack speed.

Crack Speed in Dynamic Crack Popogation
The crack speed vs. incremental crack size is shown in Figure 8.The simulation is for aluminum 2024-T351 and titanium-6Al-4V at a temperature of 300 K.The crack speed of aluminum 2024-T351 increases from below 1000 mm/s to around 2000 mm/s.A large fluc-Symmetry 2024, 16, 142 9 of 13 tuation is seen when the material is about to fail.This is similar to the prediction in classical fracture mechanics.The titanium-6Al-4V simulation shows an unstable crack speed.

Crack Speed in Dynamic Crack Popogation
The crack speed vs. incremental crack size is shown in Figure 8.The simulation is for aluminum 2024-T351 and titanium-6Al-4V at a temperature of 300 K.The crack speed of aluminum 2024-T351 increases from below 1000 mm/s to around 2000 mm/s .A large fluctuation is seen when the material is about to fail.This is similar to the prediction in classical fracture mechanics.The titanium-6Al-4V simulation shows an unstable crack speed.

Crack Tip Profile at Critical Time
Figure 9 demonstrates the critical moment right before element deletion at 300 K as an example.The highlighted elements are the first deleted elements to extend the crack length.The high ductility of the metal alloy causes the excessive deformation of the crack tip element before failure.The elements ahead of the crack edge were damaged first before the crack propagated.This is caused by the limitations of the element deletion method.The cases of 186 K, 420 K, and 537 K have similar configurations.

Crack Front Profile during Propagation
Figure 10 shows the crack's front profile during the propagation process.The crack grows along the positive X-axis, or to the right.The crack's front profile is symmetrical along the Z-axis, or normal to the paper.This is because of the symmetry in the geometry.During the process, the crack's front profile grows self-similarly, which is consistent with classic fracture mechanics theory [7,8,11].For example, the crack's front profiles in Figure 10b,h repeat at later stages in Figure 10d,j separately.We selected and recorded 12 stages in the same time interval.

Crack Front Profile during Propagation
Figure 10 shows the crack's front profile during the propagation process.The crack grows along the positive X-axis, or to the right.The crack's front profile is symmetrical along the Z-axis, or normal to the paper.This is because of the symmetry in the geometry.During the process, the crack's front profile grows self-similarly, which is consistent with classic fracture mechanics theory [7,8,11].For example, the crack's front profiles in Figure 10b,h repeat at later stages in Figure 10d,j separately.We selected and recorded 12 stages in the same time interval.moment right after the deletion of the first element.(c) The moment crack become "blunt".(d) The moment right after the target element is deleted.

Crack Front Profile during Propagation
Figure 10 shows the crack's front profile during the propagation process.The crack grows along the positive X-axis, or to the right.The crack's front profile is symmetrical along the Z-axis, or normal to the paper.This is because of the symmetry in the geometry.During the process, the crack's front profile grows self-similarly, which is consistent with classic fracture mechanics theory [7,8,11].For example, the crack's front profiles in Figure 10b,h repeat at later stages in Figure 10d,j separately.We selected and recorded 12 stages in the same time interval.

Initial Crack Size Dependency of Critical Energy Release Rate
Temperature's effect on the critical energy release rate of a specimen with different initial crack sizes is shown in Figure 11.Aluminum 2024-T351 and titanium-6Al-4V show that the initial crack length does not affect the critical energy release rate under a specific temperature.This is in agreement with classical fracture mechanics theory.This means that the critical energy release rate can be used as a constant material model parameter for predicting the extension of crack with different crack lengths.. Notice that these conclusions come from damage mechanics simulations with an advanced material model, not from fracture mechanics theory.

Initial Crack Size Dependency of Critical Energy Release Rate
Temperature's effect on the critical energy release rate of a specimen with different initial crack sizes is shown in Figure 11.Aluminum 2024-T351 and titanium-6Al-4V show that the initial crack length does not affect the critical energy release rate under a specific temperature.This is in agreement with classical fracture mechanics theory.This means that the critical energy release rate can be used as a constant material model parameter for predicting the extension of crack with different crack lengths.. Notice that these conclusions come from damage mechanics simulations with an advanced material model, not from fracture mechanics theory.

Discussion
Test data for the critical energy release rate under high and low temperatures are not found in the literature for aluminum 2024-T351 and titanium-6Al-4V.The roomtemperature fracture test data in the literature [32] only report the critical stress intensity factor.Although the critical energy release rate can be converted from a critical stress intensity factor, the conversation is based on a linear elastic material, which is not the case for aluminum and titanium.If we use the conversion equation regardless, the aluminum

Discussion
Test data for the critical energy release rate under high and low temperatures are not found in the literature for aluminum 2024-T351 and titanium-6Al-4V.The room-temperature fracture test data in the literature [32] only report the critical stress intensity factor.Although the critical energy release rate can be converted from a critical stress intensity factor, the conversation is based on a linear elastic material, which is not the case for aluminum and titanium.If we use the conversion equation regardless, the aluminum 2024-T3 alloy compact tension specimens at low strain rates give critical energy release rates ranging from 20N/mm to 100N/mm with a 2 to 6mm crack extension [32].This test for aluminum 2024-T3 is conducted in different environments (exfoliation corrosion solution, air, and sodium chloride) at low strain rates.A different result is found for aluminum 2024-T351, with critical energy release rates ranging from 24N/mm to 32N/mm using compact tension specimens with different geometries [33].Our aluminum result is 17.9N/mm.The titanium-6Al-4V critical energy release rate is reported to be around 37N/mm in a test on compact tension specimens made by additive manufacturing [34].Our result is 24N/mm.Besides the conversion error above, the material properties of the aluminum and titanium alloys may be quite different with respect to different heat treatments, plate thicknesses, manufacturing conditions, plate stocks, etc. [31].So, a comparison with a similar alloy can only serve as a reference.
An experiment on a long glass specimen [35] and laminate composite [36] shows that the critical energy release rate increases with temperature.This agrees with our result.The crack speed for aluminum 2024-T351 agrees with classical dynamic fracture mechanics theory [37].However, the crack velocity oscillations for titanium-6Al-4V need further investigation.Interestingly, a similar study reported that the crack speed oscillates when it reaches a critical value on a PMMA plate [38].
This work introduced a new method to calculate the critical energy release rate in finite element analysis.The simulation procedure is identical to a standard finite element analysis with the element deletion method.The simulation result can be converted into a critical energy release rate based on the suggested method.With this method, the damage mechanics simulation based on plasticity and failure material models can be used to calculate the critical energy release rate.The method is validated with the virtual crack opening method, and it can be used in dynamic fracture propagation.Unlike the virtual crack opening method, this method only needs a single simulation.It is worthwhile to mention that the key reason for using this method to calculate the critical energy release rate is to have an advanced material model which is correlated with experiment data on deformation and failure.
In summary, our simulation result demonstrates that the initial crack size does not affect the critical energy release rate.Therefore, the critical energy release rate fracture criterion is still valid in damage mechanics simulations with the use of plasticity and failure models as the material models.Ideally, the parameters used in fracture criteria should not be related to the specimen shape, crack length, crack mode, etc.For example, the stress intensity factor requires a manual distinction between crack modes, which limits its use in general three-dimensional problems.This work verified that the critical energy release rate can be used as a fracture criterion.

Conclusions and Future Work
The goal and major contribution of this work is to demonstrate that the critical energy release rate increases with rising temperatures for aluminum and titanium alloys.Our numerical results demonstrate this phenomenon by simulating a fracture test with a testdata-based material model.The specific contributions are summarized as follows.(1) We introduce a method to calculate the critical energy release rate with the element deletion method and verify its feasibility.(2) We demonstrate that the element deletion method is able to calculate the critical energy release rate in dynamic crack propagation.(3) We discussed the dependency of the critical energy release rate with respect to different temperatures, crack extension lengths, crack initial lengths, and crack speeds for the titanium and aluminum alloys.(4) We verified a self-similar crack profile in a dynamic crack propagation simulation.

Figure 1 .
Figure 1.The scheme of the element deletion method.(a) Initial configuration; (b) the configura right before element deletion.(c) The configuration before the crack is incremented along the id propagation direction.The yellow element at the crack surface is first deleted.(d) The bro element ahead of the crack tip is deleted next.

Figure 1 .
Figure 1.The scheme of the element deletion method.(a) Initial configuration; (b) the configuration right before element deletion.(c) The configuration before the crack is incremented along the ideal propagation direction.The yellow element at the crack surface is first deleted.(d) The brown element ahead of the crack tip is deleted next.

Figure 2 .
Figure 2. Energy vs. displacement response by element deletion method.

Figure 2 .
Figure 2. Energy vs. displacement response by element deletion method.

Figure 2 .
Figure 2. Energy vs. displacement response by element deletion method.

Figure 3 .
Figure 3. Virtual crack opening method.(a) Initial crack; (b) incremented crack.The crack tip is at the center of the yellow region.This is a side view of a three-dimensional specimen.The crack tip front is perpendicular to the two-dimensional surface and goes inside the surface.

Figure 3 .
Figure 3. Virtual crack opening method.(a) Initial crack; (b) incremented crack.The crack tip is at the center of the yellow region.This is a side view of a three-dimensional specimen.The crack tip front is perpendicular to the two-dimensional surface and goes inside the surface.

Figure 4 .
Figure 4. Force vs. displacement response by virtual crack opening method.

Figure 4 .
Figure 4. Force vs. displacement response by virtual crack opening method.

Symmetry 2024 , 14 Figure 5 .
Figure 5.The geometry and mesh of the finite element model with a perfect crack.(Left): Geometry.(Right): Zoomed-in mesh.

Figure 5 .
Figure 5.The geometry and mesh of the finite element model with a perfect crack.(Left): Geometry.(Right): Zoomed-in mesh.

Symmetry 2024 , 14 Figure 6 .
Figure 6.Finite element model for crack speed study.(Left): Whole model.(Right): Zoomed-in mesh.The colored element on the crack path is used in the calculation of energy release rate.

Figure 6 .
Figure 6.Finite element model for crack speed study.(Left): Whole model.(Right): Zoomed-in mesh.The colored element on the crack path is used in the calculation of energy release rate.

Symmetry 2024 ,
16,  x FOR PEER REVIEW 9 of 14 states of stress correspond to different failure strains, and thus, different strain energies are deleted at the moment of failure.

3. 4 .Figure 9 .
Figure9demonstrates the critical moment right before element deletion at 300 K as an example.The highlighted elements are the first deleted elements to extend the crack length.The high ductility of the metal alloy causes the excessive deformation of the crack tip element before failure.The elements ahead of the crack edge were damaged first before the crack propagated.This is caused by the limitations of the element deletion method.The cases of 186 K, 420 K, and 537 K have similar configurations.Symmetry 2024, 16, x FOR PEER REVIEW 10 of 14

Figure 9 .
Figure 9. Crack edge profile before the crack prorogation of aluminum 2024-T351 at 300 K.The green colored elements are the target elements.(a) The moment crack is about to extend.(b) The moment right after the deletion of the first element.(c) The moment crack become "blunt".(d) The moment right after the target element is deleted.

Figure 10 .
Figure 10.Dynamic crack front profile of aluminum 2024-T351 at a temperature of 300 K.The time frame increases (a-l).The element with pink color is the internal element of the original specimen.The element with white color is the new crack surface.

Figure 10 .
Figure 10.Dynamic crack front profile of aluminum 2024-T351 at a temperature of 300 K.The time frame increases (a-l).The element with pink color is the internal element of the original specimen.The element with white color is the new crack surface.

Table 1 .
Critical energy release rate by virtual crack opening and element deletion methods.

Table 1 .
Critical energy release rate by virtual crack opening and element deletion methods.

Table 2 .
Simulation critical energy release rate vs. temperature.