Mechanism and Method of Testing Fracture Toughness and Impact Absorbed Energy of Ductile Metals by Spherical Indentation Tests

To address the problem of conventional approaches for mechanical property determination requiring destructive sampling, which may be unsuitable for in-service structures, the authors proposed a method for determining the quasi-static fracture toughness and impact absorbed energy of ductile metals from spherical indentation tests (SITs). The stress status and damage mechanism of SIT, mode I fracture, Charpy impact tests, and related tests were first investigated through finite element (FE) calculations and scanning electron microscopy (SEM) observations, respectively. It was found that the damage mechanism of SITs is different from that of mode I fractures, while mode I fractures and Charpy impact tests share the same damage mechanism. Considering the difference between SIT and mode I fractures, uniaxial tension and pure shear were introduced to correlate SIT with mode I fractures. Based on this, the widely used critical indentation energy (CIE) model for fracture toughness determination using SITs was modified. The quasi-static fracture toughness determined from the modified CIE model was used to evaluate the impact absorbed energy using the dynamic fracture toughness and energy for crack initiation. The effectiveness of the newly proposed method was verified through experiments on four types of steels: Q345R, SA508-3, 18MnMoNbR, and S30408.


Introduction
Fracture toughness and impact absorbed energy are two crucial mechanical parameters that reflect a material's resistance to fracture.Accurate measurements of these two parameters are of great significance in evaluating the structural integrity of materials and ensuring equipment safety in long-term service.Conventional fracture tests and Charpy impact tests require destructive sampling, and thus cannot be used on in-service structures.
A method for estimating the quasi-static fracture toughness of materials from SITs was first proposed by Byun et al. [6], in which the work done by the average contact pressure in the direction of the impress was defined as the indentation energy necessary to fracture W IEF .Then, the fracture toughness K IC in N/mm 3/2 of the sample material can be determined by applying the Griffith theory [7] where h and h C in mm are the indentation depth and critical indentation depth, respectively; a C in mm is the critical contact radius (the contact radius a is the radius of the projection area of the indentation, and can be given as (2Rh − h 2 ) 0.5 , where R is the radius of the indenter), and P in N is the external load ( P = Sh , where S in N/mm is the slope of the loading line).
where E is the Young's modulus of the specimen material.W F in kJ/m 2 is the energy of formation for a unit area of a crack and can be approximated as W IEF in kJ/m 2 .
To overcome the drawback of Byun's method, which requires an additional destructive tensile test, Haggag [8] empirically judged the critical indentation depth by considering the difference between brittle and ductile fractures (named Haggag toughness method, HTM).The calculation formula for fracture toughness was modified by considering the lower platform of the fracture toughness at low temperatures, as shown in Eq. (3): The HTM is applicable to a wide range of temperatures and does not require any additional destructive tests.However, the judgment of the critical indentation depth is empirical, and its effectiveness largely depends on the specimen materials [9].Additionally, the HTM does not provide any explanation of the mechanism for obtaining the fracture toughness from SITs.
Using continuum damage mechanics (CDM), the mechanism of the fracture toughness calculation from SITs was investigated by Lee [10], who developed the critical indentation energy (CIE) model by finite element (FE) analysis and microscopic observation.The critical indentation energy is calculated as follows: where the term on the left represents the energy required for the formation of two crack surfaces.In Lee's definition, the value of W f is only half of that defined by Haggag or Byun.
According to Griffith's theory, the fracture toughness K IC can be calculated using Eqs.( 4) and ( 2): (1) The CIE model, which introduced CDM to obtain the critical indentation depth, has been widely applied in previous studies [11][12][13].However, the critical damage variable D* used in the previous CIE model was determined from the analysis of a moving crack tip, in which the specimen was considered to be an elastic and perfectly plastic material (i.e., the hardening behavior was neglected) [14].Additionally, the difference between the damage at high stress triaxiality and low stress triaxiality was neglected [15].
In contrast to numerous fracture toughness calculation models, few studies have evaluated the impact absorbed energy by SITs.The existing researches have primarily focused on correlating the impact absorbed energy to the quasi-static and dynamic fracture toughness [16].Barsom et al. [17] conducted a large number of impact tests and conventional fracture tests, and proposed empirical formulae for the impact absorbed energy (KV 2 ), static fracture toughness (K IC ), and dynamic fracture toughness (K ID ).When J-integral theory was proposed, Hübner [18] established the relationship between fracture toughness and impact absorbed energy (J and KV 2 ) from the perspective of energy.Combining experimental and numerical approaches, Smith [19] contrasted the fracture behavior of pre-cracked fracture specimens tested quasi-statically with the Charpy V-notch specimens tested both quasi-statically and dynamically.They proposed an energy scaling model for relating the impact of the absorbed energy to the fracture toughness, which is suitable for both static and dynamic loading.Obtaining the impact absorbed energy from an empirical model is simple and reliable, but an empirical model often requires a large number of destructive tests for the determination of the fracture toughness, and these tests are more complex than impact tests.
To summarize, it can be seen from the aforementioned research that there are two or three methods for obtaining the fracture toughness of materials using SITs, but that the calculation accuracy of these models needs to be improved [20].Studies on the fracture toughness and absorbed impact energy are relatively common.However, most testing methods require destructive sampling, which may not be suitable for inservice structures.
The authors of this study focused on proposing a method for determining the quasi-static fracture toughness and impact absorbed energy of ductile metals from SITs.The main technical route used in this (5) study is shown in Figure 1.The mechanism analysis for specimens under the SITs, mode I fracture tests, and Charpy impact tests were carried out with microobservation and finite element (FE) calculations.The similarities of these damage mechanisms are discussed through comparison.Subsequently, based on the existing theory [21,22] for mechanism correlation, a correlation analysis of the three damage mechanisms was conducted.Based on the correlation analysis of the quasi-static fracture toughness, dynamic fracture toughness, and energy for crack initiation, calculation methods for the fracture toughness and impact absorbed energy by SITs were proposed.Finally, the reliability of the proposed method was verified through experiments on four types of steel.

Materials
Four types of steels commonly used for pressure vessels, Q345R, SA508-3, 18MnMoNbR, and S30408, were used in quasi-static tensile tests, dynamic tensile tests, mode I fracture tests, SITs, and Charpy impact tests.The chemical composition of each material was tested using the Bruker Tasman Q4 spectrograph and the results are provided in Table 1 with the corresponding standards.Metallographic observations of Q345R, SA508-3, 18MnMoNbR, and S30408 are shown in Figure 2. It was found that the average grain size did not exceed 50 μm for all four metals, which indicates that the indenter can cover at least three grains even in the 1st cycle.Therefore, the effect of the indentation size was not considered.

Quasi-Static and Dynamic Tensile Tests
Quasi-static and dynamic tensile tests were conducted to obtain the tensile properties of the experimental materials under different strain rates.The quasi-static tensile specimens were machined into rod-shaped specimens with a diameter of 10 mm, a gauge length of 50 mm, and a total length of 194 mm in accordance with the Chinese standard [26].The specific geometries of the tensile specimens are shown in Figure 3. Tensile tests were performed on an AG-IC universal testing machine at a loading speed of 1 mm/min until breakage, and the displacement was measured using an extensometer.
Dynamic tensile specimens were machined into plate specimens with a thickness of 3 mm, a gauge length of 40 mm, and a width of 12.5 mm in accordance with the Chinese standard [26].A detailed configuration of the dynamic tensile specimens is shown in Figure 4.In accordance with the strain rate of the Charpy impact test (10 s -1 ), strain-controlled (set tensile speed at 0.4 m/s) tests were conducted on HTM-5020.
In the high-strain-rate tests, the strain of a specimen was determined using digital image correlation (DIC) instead of the common extensometer measurements.Before the tests, speckles were sprayed onto the gauge  section of the specimens, and the speckle quality was calibrated statically.The speckle patterns in the tensile tests were captured using a super-high-speed digital camera (APX-RS) at a speed of 10 6 s -1 .Subsequently, the plastic strain distribution on the specimen surface was obtained using a commercial DIC system VIC-2D 6 (Correlated Solutions) and the captured images.

Spherical Indentation Tests
The specimens were cut into blocks with the dimensions 10 × 10 × 50 mm 3 and polished to a smooth surface.SITs were carried out on the stress-strain microprobe system B4000, as shown in Figure 5, using a tungsten carbide indenter with a diameter of 0.76 mm (Young's modulus and Poisson's ratio of the indenter were 710 GPa and 0.23, respectively).Considering the quasi-static characteristics of the test, the loading and unloading rates were set to 0.1 mm/min.The load-depth curve obtained from the 8-cycle spherical indentation test is shown in Figure 6, whereby the maximum indentation depth h max was set to 0.24R.

Mode I Fracture Tests
Mode I fracture tests were conducted to determine the fracture toughness of the experimental materials.A  standard compact tensile (CT) specimen with a thickness of 0.45 inches and a ratio of nominal width to nominal thickness of 2 was used, as shown in Figure 7. Precracking was performed on a high-frequency fatigue testing machine with the stress ratio and frequency set to 0.1 and 100 Hz, respectively.The fatigue crack length on the surface was recorded every 20000 cycles during pre-cracking while reversing the specimen to ensure that the crack lengths on both sides of the specimen were essentially the same.To strengthen the plane strain state of the fatigue crack tip, a side groove with a depth of 1.5 mm and an angle of 90° was machined onto the surface of the specimen after pre-cracking.CT tests were also performed on a stress-strain microprobe (SSM) system B4000.Preformed cracks propagate slowly and steadily in ductile metals; thus, the plane strain fracture toughness K IC cannot be ascertained directly.According to the relevant standard [27], the unloading compliance method can be used on a single specimen to determine the crack initiation resistance curve of ductile materials, that is, the J-R resistance curve.The intersection point of the structural line with a 0.2 mm offset and the J-R curve was used as the engineering fracture toughness, that is, J IC in kJ/m 2 .Then, K IC was calculated as

Charpy Impact Tests
Charpy impact tests were conducted to obtain the crack initiation energy and the impact absorbed energy.The test was performed on an RKP450 impact testing machine with a hammer in accordance with the requirements of the Chinese standard [28].To visually observe the processes of crack initiation and extension, the test was carried out on a self-designed tooling that stops the hammer after crack initiation and protects the sample from fracture.The configurations of the Charpy (6) impact test specimen and the tooling for crack arrest are shown in Figures 8 and 9, respectively.In Figure 9(a), the red box represents a cushion block, the green box represents the specimen, and the blue box represents the hammer.It has been reported that ductile metal materials crack when the hammer moves approximately 1.8 mm after impact along the end face of the specimen [29].To ensure the effectiveness of the test, the crack arrest tests were carried out with cushion blocks of different thicknesses (11.7 mm, 11.8 mm, and 11.9 mm).

FE Simulation
Considering that the damage to ductile materials is closely related to their stress state, the SIT, mode I fracture, and Charpy impact tests were first numerically investigated using ABAQUS 6.14.The purpose of these simulations was to analyze the stress status of the specimens.Thus, a damage criterion was not developed in the calculations, and cracks were not allowed to propagate.

FE Modelling of Mode I Fracture Tests
Considering the characteristics of the CT specimen (mode I fracture), an FE model was developed.As a contrast to the SEM observations described in Section 4, the material properties of the specimens were characterized using a previously performed   quasi-static tensile test on Q345R.In this calculation, the contact between the pin and the specimen was regarded as a surface-to-surface contact, with the friction coefficient set to 0.2 [30].To ensure calculation accuracy and reduce cost, a 8-node linear brick with reduced integration (C3D8R in ABAQUS) was used.The meshes around the pre-cracked region (set to 2 mm × 2 mm) were gradually refined to 0.05 mm, while those in the distal region were set to around 0.4 mm, as shown in Figure 10.Similar to the fracture toughness test, the pin on the left of the model was fixed, and the pin on the other side was displacement-controlled with a maximum displacement of 3 mm.
A comparison of the load P-displacement q curve from the mode I fracture tests and the numerical simulations with different mesh densities is provided in Figure 11.The grid independence of the model was verified by refining the mesh size at the pre-crack tip to 0.15 mm, 0.10 mm and 0.05 mm.The initial part of the P-q curves from the FE calculations coincided well with the initial parts of those from the experiment.As the elastoplastic FE calculation did not consider the crack propagation in the CT specimen, the FE calculated Pq curves gradually deviated from their experimental counterparts after the crack began to propagate.

FE Modelling of Charpy Impact Tests
Considering the characteristics of the Charpy impact specimen, a quarter model, as shown in Figure 12, was created to determine the stress distribution and the variation in stress triaxiality σ Tri at the V-notch tip.In this simulation, the Johnson-Cook model [31], shown in Eq. ( 7), was chosen for describing the mechanical properties of the specimen (the temperature was not considered in this simulation).
where σ t is the flow stress, ε p is the equivalent strain, ε * represents the plastic strain rate, ε r * is usually taken as 1 s -1 and is the reference plastic strain rate.A, B, C, and n are material parameters that can be determined by fitting the data from the quasi-static and dynamic tensile tests of Q345R (A = 332.53,B = 1015.68,C = 0.0128, and n = 0.56).
8-node linear brick elements with reduced integration (C3D8R in ABAQUS) were used.The meshes around the V-notch and hammer tips were refined to 0.1 mm for calculation accuracy, and the overall mesh sizes of the specimens were 0.5 mm.According to the impact test, the initial velocity of the hammer was set to 5.2 m/s in the direction of the y-axis.The support was regarded as a rigid body, and its movement was fixed.

FE Modelling of SITs
Considering the characteristics of the SITs, an axisymmetric model was established for analyzing the (7)  stress state, as shown in Figure 13.The indenter was considered as a rigid body whose diameter was the same as of that used in the experiment, while the mechanical properties of the specimen were determined by quasistatic tensile tests of Q345R.
The element type of the specimen was a 4-node bilinear axisymmetric quadrilateral with reduced integration (CAX4R).The contact settings of the SITs were similar to those in the FE model of the mode I fracture test, where the indenter and specimen surfaces were defined as the master and slave surfaces, respectively.According to a previous study [21], the minimum mesh size in the region near the indenter was 0.002 mm, and the maximum mesh size in the edge region was 0.018 mm.The FE calculation of the SIT involved displacement control, and the maximum depth was set in accordance with the SITs (i.e., h max = 0.24R).

Calculation Method for Fracture Toughness
Based on SITs

Stress State Analysis 4.1.1 Mode I Fracture Tests
The von Mises stress distribution of the pre-crack tip with a loading line displacement q of 1.5 mm (as the FE calculated P-q curves deviated from their experimental counterpart for q > 1.5 mm) on the specimen is shown in Figure 14(a).The maximum von Mises stress appears on the front edge of the pre-crack tip, and the stress concentration zone surrounding the crack has the shape of a pincer, which is the same as the theoretically determined plane strain plastic zone shape [7].
To better understand the stress state of the CT specimen, the stress triaxiality σ Tri of the region before of the pre-crack tip, which was determined using Eq. ( 8), is shown in Figure 14(b): where σ m and σ eq are the hydrostatic and von Mises equivalent stresses, respectively.It can be observed that the stress triaxiality was relatively high near the bottom of the pre-crack in the specimen, which is regarded as a typical high stress triaxiality region (σ Tri ≥ 0.4).According to relevant studies [15,32], failure of the high-stress triaxiality region is always caused by void nucleation, growth, and coalescence, which is called void accumulation fracture.

Spherical Indentation Tests
Similar to the spherical indentation tests, the simulation process consisted of eight times equal depth loading and unloading.The von Mises stress and shear stress (8) σ Tri = σ m σ eq , distributions at the end of the 8th loading (h = 0.24R) around the indentation area are presented in Figures 15  and 16, respectively.It can be seen that the distributions of the von Mises stress and the shear stress are basically the same and approximately in the shape of a spherical crown, where the maximum values are located in the 45°-60° range inclined to the y-axis.Owing to friction, the plastic flow at the indentation bottom was suppressed, leading to a relatively small stress at the indentation center.
The stress triaxiality σ Tri in the area around the indentation, as shown in Figure 17, is almost entirely negative.The minimum stress triaxiality was less than − 1.6 at the bottom of the indentation, while the maximum was located at the edge of the indentation with − 0.2.This indicates that the region around the indentation is a typical low-stress triaxiality region (σ Tri ≤ 0), where the fracture always occurs due to shear stress-induced damage [15,32].

Damage Mechanism Analysis
All microscopic observations were performed using the same material (Q345R) to ensure consistency.The micromorphology of the indentation section surface is shown in Figure 18, in which several wedge-shaped voids are visible around the border of the residual indentation.This is similar to the shear damage explained by dislocation pile-up [33].It was found that many wedge-shaped voids are located in the 'wing region' (i.e., the 45°-60° range inclined to the y-axis), while almost no voids were observed at the indentation bottom.The damage distribution is consistent with the shear stres s concentration presented in Section 4.1.2.Thus, it can be concluded that shear stress is the main reason for the initiation and localized distribution of voids in the indentation test, and that the damage mechanism in SITs is the same as that in the pure shear test.As reported in a previous study [21], the mode I fracture surface exhibits a typical dimple feature, which is similar to the tensile fracture.The micromorphology near the crack tip is shown in Figure 19, with several spherical or ellipsoidal voids visible in the region under mode I loading.To analyze the relationship between the fracture mode and the stress status, a comparison of the stress triaxiality-load deflection curves at the stress concentration zone in the FE model of SITs and the mode I fracture (which corresponds to the damage concentration area in microscopic observations) is provided in Figure 20.σ Tri increased significantly with h at the beginning of loading, but finally stabilized at approximately − 1.4.Combined with the microstructure of the residual indentation, it was proved that the inner damage of the material was caused by shear stress, and that the material more easily produced shear failure in the low stress triaxiality region.As shown in Figure 20(b), σ Tri varies between 0.8 and 1.2, indicating that the damage in the mode I specimen is similar to that in the uniaxial tensile tests, originating from the nucleation, growth, and coalescence of voids in the high-stress triaxiality region [15,32].
The mechanisms of tensile and shear fractures are essentially different.However, the initiation and growth of voids and microcracks lead to a decrease in the effective bearing area, and thus, cause a reduction in the effective elastic modulus (Young's modulus is used to characterize tensile damage, and the shear modulus is used to characterize shear damage).According to continuum damage mechanics, the tensile damage variable D I and the shear damage variable D II can be determined using Eqs.( 9) and (10), respectively [34]: where E eff , G eff and G are the effective Young's, effective shear, and shear moduli of the undamaged material, respectively.
Shi et al. [22] proposed a critical damage strain energy release rate Y * , as shown in Eqs. ( 11) and ( 12), for the tensile and shear forces, respectively is a material constant that is independent of the stress state (i.e., where v is Poisson's ratio and can be considered as a material constant for each type of metal (e.g., 0.3 for steels used in this study).D I * is the critical tensile damage variable and D II * is the critical shear damage variable of the material.
According to von Mises yield criterion, there is:  substituting Eq. ( 11), and Eq. ( 12), into Eq.( 13), yields Eq. ( 14): According to a mechanism analysis, the critical damage variables under the mode I fracture test and the SITs correspond to D I * and D II * , respectively.

Modified CIE Model
The previous CIE model creatively introduced CDM into the indentation test, but the model has some shortcomings, as mentioned in Section 1. First, the critical damage variable D* used in the previous CIE model was derived from Andersson's numerical analysis [35], in which the object was assumed to be an elastic perfectly plastic material (i.e., the effect of work hardening was neglected).In addition, D* has not yet been experimentally verified.Moreover, the difference between the damage mechanisms in high stress triaxiality and low stress triaxiality regions was neglected in the model.
To obtain the critical damage variable of the materials, cyclic loading tensile tests were performed, as mentioned in Section 2.2.The true stress-strain curve of Q345R is shown in Figure 21.Considering the effect of necking on the determination of the stressstrain curve, only the data before obvious necking were used.In addition, only the first 40% of the unloading curve data were used to avoid the effect of back-stress relaxation [36].
The variation in D I with the true strain is shown in Figure 22.The damage variable of the four materials initially increased with an increase in the true strain and (13) .
then tended to remain constant.As shown in Figure 22, the value of the critical damage variable D I * ranged between 0.20 and 0.23, which is consistent with D I * (approximately 0.22) as measured through tension tests by predecessors [34].To obtain a more conservative fracture toughness calculation for engineering applications, a conservative value of D I * = 0.2, was used for all materials in this study.According to Eq. ( 14), the critical shear damage variable of all materials, D II * was calculated to be 0.25.Because the damage around the indentation was almost entirely caused by shear stress, the critical indentation depth h c should be determined by D II * .To determine E eff from SITs, the Pharr-Oliver formula was used in the previous CIE model [37].However, it is important to note that the plastic deformation of the specimen before unloading was not considered in the Pharr-Oliver formula, that is, Eq. ( 15), which leads to an evident reduction in the effective Young's modulus at a small indentation depth [37,38]: where v and v ind are the Poisson's ratios of the specimen and spherical indenter, respectively; L in N/mm is the slope of the unloading curve in the load-displacement curve obtained from the SITs, and E ind is the Young's modulus of the spherical indenter.
To consider the plastic deformation of the material during SITs, a modified effective Young's modulus formula was proposed in our previous study [30], as follows: (15)  R 0 and h r in mm are the radius of the residual indentation and secondary loading depth, respectively, and are obtained using Eqs.( 17) and ( 18), respectively: where h max and h p are the maximum and residual indentation depths, respectively, as shown in Figure 23.
It was found from experiments that the relationship between lnE eff and lnh can be described by a linear function, as shown in Figure 24.( 16) Similar to the previous CIE model, the critical effective Young's modulus E eff * can be calculated using Eq. ( 19), where D II * and E are known: The critical indentation depth h c can then be determined using the fitting curve shown in Figure 24.The fracture toughness can be calculated using Eq. ( 5), with h C .An illustration of the calculation processes of the previous and modified CIE models is shown in (19)

Stress State Analysis of Impact Test
According to the displacement and the drop speed of the hammer during the crack arrest test, the Mises stress and the stress triaxiality distribution of the impact specimen when h reached 1.8 mm (initiation of crack propagation) are shown in Figure 26.The stress concentration zone can be seen in the vicinity of the V-notch, which is in the shape of a bowknot.The stress gradually decreased from the middle to both sides, with the maximum located at 1 mm in front of the bottom of the V-notch.
From Figure 26(b), it was found that the stress triaxiality around the stress concentration zone was greater than 0.4, which is regarded as a typical high-stress triaxiality region, indicating that the main factors resulting in to the specimen were void nucleation, growth, and coalescence [15,32].

Damage Mechanism Analysis
The SEM observations of the section surfaces of the crack-arrest specimens are shown in Figure 27.Voids and microcracks can be observed in front of the crack tip, which is located in the high-stress triaxiality region mentioned in the numerical analysis.According to micro-mechanics, the initiation of voids can be explained by the fact that the local stress in the damage region could effectively offset the bonding force between the particles and the matrix.Moreover, the volume of voids increases rapidly, particularly in the high-stress triaxiality region.With the development of voids, adjacent voids coalesce, leading to the initiation and propagation of microcracks.This damage morphology was similar to that of the CT specimen, as shown in Figure 19.An observation of the fracture surface of the impact test specimen (Q345R), shown in Figure 28, shows obvious dimple features,  which are similar to those of the CT specimen [21].The variation in σ Tri with the displacement of the hammer, shown in Figure 29, also proved the impact damage formed in the high stress triaxiality region, similar to the CT specimen shown in Figure 20(b).Therefore, it can be concluded that the impact specimen shares the same tensile damage mechanism with the CT specimen.

Correlation between Dynamic and Quasi-Static Fracture Toughness
High strain rate loading affects the yield strength of the specimens, which is known as the strain rate effect.Therefore, the variation in yield strength can reflect the strain rate effect of the materials to a certain extent.Based on this, Chaouadi et al. [39] proposed a method for characterizing the relationship between the dynamic and quasi-static fracture toughness (J IC to J ID ) using yield strength: where σ YD in MPa is the dynamic yield strength, σ YS in MPa is the quasi-static yield strength, and α is the proportion coefficient.This method was verified by Pan et al. [16] through a large number of experiments on Q345R, from which the constant α was determined as being 0.28 (with a (20) ratio of J IC to J ID of 0.44 determined from quasi-static and dynamic J-R curves).The engineering stressstrain curves from the dynamic tensile tests mentioned in Section 2 are presented in Figure 30, with the corresponding yield strengths shown in Table 2.The average ratio of σ YD to σ YS was 1.258.Compared with the ratio in Ref. [34], the ratio of J IC to J ID determined from the tests in this study was 0.443 (the difference was less than 1%).

Correlation between Dynamic Fracture Toughness and Impact Absorbed Energy
Due to the complexity of dynamic fractures, the solution of most problems still has to be based on numerical analysis or empirical formulae, with a few exceptions, such as an infinite crack body and a plane crack in elastic-perfectly plastic material.However, the finite element method cannot be used to obtain the dynamic fracture parameters of a material.Therefore, the empirical correlation method is currently the main approach for analyzing dynamic parameters [40,41].
Based on the impact tests conducted by Chaouadi [29], the location of crack initiation can be determined from the load-displacement curve.As shown in Figure 31, crack initiation occurred at a load between the general yield (i.e., F gy ) and the maximum load (i.e., F max ).The initiation load can be written as The impact absorbed energy KV 2 is the area under the load-displacement curve, which can be obtained by integrating over the entire curve.The energy-to-crack initiation KV 2i is the area of the purple zone, which can also be determined by the curvilinear integral.
The experimental data were analyzed from the viewpoint of energy, and it was found that the ratio of KV 2i to KV 2 is η, which is approximately equal to a constant, as shown in Eq. ( 22): The Charpy impact tests shown in Section 2.5 yield an η of around 7%, which coincided well with those provided in previous studies [16,29].Thus, the authors attempted to determine the relationship between the dynamic fracture toughness J ID and the impact absorbed energy KV 2 by correlating J ID and KV 2i .The specific analysis is as follows.
Comparing the damage mechanisms mentioned in Sections 4.2 and 5.2, it was concluded that the V-notch tip of the impact specimen and the crack tip of the mode I fracture specimen were observed in the high stress triaxiality region.Furthermore, they shared the same fracture mechanism, in which the damage was caused by the nucleation, growth, and coalescence of voids under tensile stress fields.Therefore, the Charpy impact test can be correlated with the mode I fracture test, which helps to prove that there is a certain quantitative relationship between the dynamic fracture toughness and the crack initiation energy, and also explains that it is rational to calculate the impact absorbed energy by SITs.
As shown in Table 3, the quasi-static fracture toughness J IC was obtained by the mode I fracture test, and J ID was calculated using Eq. ( 20).The crack initiation energy was determined using the load-displacement curve from (21)   the Charpy impact test.The KV 2i -J ID plot and fitting curve are shown in Figure 32.The crack initiation energy KV 2i can be calculated using the fitting formula shown in Eq. ( 23): In combination with Eq. ( 22), the impact absorbed energy KV 2 of a material can be obtained.The value of η is discussed in Section 6.

Fracture Toughness by Modified CIE Model
According to the SITs mentioned in Section 2.3, the effective Young's modulus was calculated using the modified formula in Eq. ( 16), and the Pharr-Oliver formula in Eq. (15) (i.e., E eff and E PO ), and the E PO and E eff from the 1st indentation cycle were assumed to be the E of the undamaged material.Comparisons of E PO and E eff from the 1st indentation cycle with those from the tensile test (E T ) are presented in Table 4. From the comparison, it was found that the error between E PO and E T is up to 30%.However, the error between E eff and E T is less than (23) KV 2i = 0.00036J 2 ID − 0.3127J ID + 79.749.
6%, which proves the validity of the modified effective Young's modulus model.The load-displacement curves obtained from the spherical indentation tests are shown in Figure 33.The loading curve is approximately a straight line, which is the result of a superposition of the two nonlinear behaviors in the SITs.One of these is the nonlinearity of the loaddisplacement curve caused by the geometric nonlinearity of the spherical indenter, and the other is the nonlinearity of the work hardening of the materials.Combined with the load-displacement curve shown in Figure 33, the fracture toughness was calculated using the previous CIE model and the modified CIE model (i.e., K ICP and K ICM , as calculated from Figure 25).A comparison of K ICP and K ICM, with the average results of the mode I fracture test (K ICT ) in Section 2.5 is shown in Table 5.
By comparison, it was found that the fracture toughness calculated using the previous CIE model was higher than that of the mode I fracture test.The minimum error of the previous model was more than 20%, and the error for S30408 was up to 51.82%.The larger error is mainly because the effective Young's modulus obtained in the previous CIE model is not sufficiently accurate [35].Additionally, the value of the critical damage variable from the previous model was larger (illustrated in Section 4.3), which also led to a large error in the fracture toughness.Compared with the previous CIE model, the modified CIE model result was closer to the experimental results.Except for the error between K ICM and K ICT for 18MnMoNbR, which is approximately 15%, that of other materials is approximately 5%.In addition, the fracture toughness calculated by the modified CIE model is conservative, which is desirable in engineering applications.To address the problem of conventional approaches for mechanical property determination requiring destructive sampling, which may be unsuitable for in-service structures, a method for determining the quasi-static fracture toughness and impact absorbed energy from SITs is proposed in this study.The stress status and damage mechanism of specimens under SITs, mode I fracture tests, and Charpy impact tests were first investigated using FE calculations and SEM observations.Subsequently, the widely used CIE model for quasi-static fracture toughness determination from the SITs was modified.The quasistatic fracture toughness, as determined by the CIE model, was employed for the assessment of impact absorbed energy, in conjunction with dynamic fracture toughness and crack initiation energy.The effectiveness of the newly proposed method was verified through experiments on four types of steel: Q345R, SA508-3, 18MnMoNbR, and S30408.The following conclusions were drawn from the investigations mentioned above.
(1) The damage in the mode I fracture originates from void nucleation, growth, and coalescence in the high stress triaxiality region, which is similar to the results from uniaxial tensile tests.However, the indentation damage originates from a dislocation pile-up in a low stress triaxiality region, which is similar to that from pure shear tests.Therefore, the SIT and mode I fractures can be correlated by determining the relationship between the uniaxial tensile and pure shear tests.(2) From cyclic-loading tensile tests on the four metal materials, it was found that the critical tensile damage variable varied between 0.20 and 0.23 (conservatively set to 0.2).Correspondingly, the critical shear damage variable, which is applicable to SITs, was 0.25.(3) Compared with the previous CIE model, the modified CIE model considers the damage mechanism and is more effective in the quasi-static fracture toughness calculation (with the maximum error decreasing from 51.82% to 17.32%).(4) When the impact absorbed energy is calculated from the SITs using the material's own crack initiation energy ratio (i.e., η = η i ), the maximum error is approximately 25%.Enhanced calculation accuracy can be achieved upon the development of a more comprehensive database for crack-initiation energy ratios.In this case, SIT can be considered as an alternative method for evaluating the impact of absorbed energy.

Figure 1
Figure 1 Main technical route

Figure 7 Figure 8
Figure 7 Configuration of the CT specimen (in mm)

Figure 9
Figure 9 Configuration of the crack arrest tooling (in mm)

Figure 11 Figure 12
Figure 10 Illustration of the FE model of the mode I fracture tests

Figure 13 Figure 14
Figure 13 Illustration of the FE modeling of SITs

Figure 15 Figure 16 Figure 17 Figure 18
Figure 15 Mises stress distribution around indentation

Figure 19
Figure 19Observation of section surface for mode I specimen[21]

Figure 20
Figure 20 Comparison of σ Tri from FE calculations on SITs and the mode I fracture at damage zone

Figure 26 Figure 27 Figure 28
Figure 26 Mises stress and stress triaxiality distribution of specimen under impact test

Figure 29 Figure 30
Figure 29 Variation of σ Tri with displacement of hammer x (region of crack tip)

Figure 32
Figure 32 KV 2i -J ID plot from tests and the fitting curve

Table 1
Chemical composition of the material and corresponding standard values

Table 2
Dynamic and static yield strengths of each material [29]re31Determination of the onset of crack initiation using the load-displacement curve[29]

Table 3
Fracture toughness and impact crack initiation energy from tests

Table 4
Comparison of E PO and E eff from the 1st indentation cycle