Modeling of Hydrogen-Charged Notched Tensile Tests of an X70 Pipeline Steel with a Hydrogen-Informed Gurson Model

Hydrogen can degrade the mechanical properties of steel components, which is commonly referred to as “hydrogen embrittlement” (HE). Quantifying the effect of HE on the structural integrity of components and structures remains challenging. The authors investigated an X70 pipeline steel through uncharged and hydrogen-charged (notched) tensile tests. This paper presents a combination of experimental results and numerical simulations using a micro-mechanics-inspired damage model. Four specimen geometries and three hydrogen concentrations (including uncharged) were targeted, which allowed for the construction of a fracture locus that depended on the stress triaxiality and hydrogen concentration. The multi-physical finite element model includes hydrogen diffusion and damage on the basis of the complete Gurson model. Hydrogen-Assisted degradation was implemented through an acceleration of the void nucleation process, as supported by experimental observations. The damage parameters were determined through inverse analysis, and the numerical results were in good agreement with the experimental data. The presented model couples micro-mechanical with macro-mechanical results and makes it possible to evaluate the damage evolution during hydrogen-charged mechanical tests. In particular, the well-known ductility loss due to hydrogen was captured well in the form of embrittlement indices for the different geometries and hydrogen concentrations. The limitations of the damage model regarding the stress state are discussed in this paper.


Introduction
Hydrogen gas is an important energy carrier in the transition towards a low carbon economy [1]. Metal infrastructure has been considered for the storage and transport of hydrogen gas. However, the presence of hydrogen atoms in the steel can cause a loss in mechanical properties, such as ductility, fracture toughness, or fatigue resistance [2,3]. This phenomenon, known as 'Hydrogen Embrittlement' (HE), has to be analyzed carefully to avoid structural integrity issues leading to accidents.
Predicting fracture is of crucial importance when assessing structural integrity [4]. Continuum computational methods can be employed as an alternative to experimental testing [5,6]. These methods can describe damage development [7,8], replicate experimental test results, and support predictions of the structural behavior of structures with a complex geometry and/or loading conditions [9][10][11]. Additionally, considering the complexity of performing hydrogen-charged mechanical tests, a computational method for describing HE would be desirable. However, the accurate modelling of HE is difficult due to its complicated nature [12]. acceleration in the void coalescence process, which was also based on the HELP mechanism. Finally, Lin et al. [33] showed the potential of an adapted Gurson model to include failure by decohesion according to the HEDE mechanism. It should be noted that Yu et al. [31], Depraetere et al. [32] and Lin et al. [33] only considered the development of the model, and did not apply the model yet to compare with experimental tests. Asadipoor et al. [20], on the other hand, employed a Gurson model for modeling hydrogen-charged flat strip tensile tests. However, only one geometry was modeled, and instead of making the Gurson parameters dependent on the hydrogen concentration, all Gurson parameters were simply recalibrated for the hydrogen-charged specimens. In other words, this approach is similar to the fracture strain approaches due to its phenomenological nature. Furthermore, since only a single geometry is investigated, the transferability of the obtained Gurson parameters to other geometries is uncertain.
In conclusion, all of the aforementioned studies concerning the continuum modeling of the plasticity-based mechanisms for accurately predicting HE lack at least one of the following three aspects. Firstly, some studies do not include hydrogen diffusion. Secondly, most studies consider a single geometry for determining the degradation parameters, thereby effectively assuming that HE is independent of the stress state. Thirdly, some studies make use of a phenomenological approach such as the fracture strain model, which complicates the connection with the physics observed during experiments.
The current work reports the experimental test data of an API 5L [34] X70 pipeline steel using four different geometries and three hydrogen concentrations (including uncharged specimens). Based on the macro-mechanical data and evidence of the failure mechanisms, a hydrogen-informed Gurson model was fitted. Finally, the hydrogen-affected fracture locus was analyzed to investigate the fracture process over a range of positive stress triaxialities.

Materials
The current study used the experimental test data of an API 5L [34] X70 pipeline steel. Test specimens originated from a pipeline that has been in use for approximately 30 years and features a banded ferrite-pearlite micro-structure with some degree of plastic anisotropy [28]. The chemical composition is given in Table 1. For a more detailed microstructural characterization of the material, the reader is referred to Cauwels et al. [35]. Four different tensile test geometries were extracted from the pipe along the longitudinal direction (also corresponding to the rolling direction of the steel), as shown in Figure 1. Both smooth (R∞) and double-notched (R6, R2, and R1.2) specimens were tested and covered a broad range of stress triaxialities. The specific levels of stress triaxiality associated with these geometries are discussed in Section 4. The use of two notches allowed for the post-mortem micro-scopic investigation of two separate deformation states: maximum force (considering the unfailed notch) and complete failure (considering the failed notch).

Experimental Methods
Prior to tensile testing, hydrogen charging was performed on a subset of the samples by submersion for 6 h in an electrolytic cell with a current density of 0.8 mA/cm 2 . Two different electrolytes were used, which resulted in two different nominal hydrogen concentrations C L,0 after charging. The first electrolyte was 0.1 M NaOH with 1 g/L NH 4 SCN, which resulted in C L,0 = 0.36 ± 0.03 wppm. The second electrolyte was 0.5 M H 2 SO 4 with 1 g/L CH 4 N 2 S, which resulted in C L,0 = 1.09 ± 0.33 wppm. These con-centrations were measured with hot extraction (G8 Galileo) and refer to the diffusible hydrogen concentrations.
After charging, the specimens were tensile tested at a fixed displacement rate. This rate was adapted for each geometry, such that the global strain rate in each specimen is equal to 2.5 × 10 −4 s −1 . The reference length for the strain rate is the total 'reduced area section length', which is denoted as l 0 for each geometry on Figure 1. Note that the local strain rate in the double-notched specimens effectively doubles after maximum force, since all deformation will concentrate in the 'weakest' notch. The time between extracting the specimen from the electrolyte and starting the tensile test was approximately 7 min. During the test, the lateral contractions of the specimen were monitored using two perpendicularly positioned cameras. This technique is also described in [28,36]. After post-processing, the diameters in both the transversal-to-rolling T direction D T and the through-thickness S direction D S were extracted. Consequently, the minimal instantaneous cross-sectional area was calculated as A = πD T D S /4. The advantage of this technique is that the strain values obtained are local strains at the minimal cross-section, rather than global strains averaged over the gauge length of an extensometer. Additionally, this enables the extraction of post-necking true stress versus true strain behavior [37]. Using the measured force F and the initial minimal cross-sectional area A 0 , the true stress σ t and the true strain t are calculated as: Table 2 presents an overview of the experimental test data used for the calibration of the damage model parameters. Note that only a limited number of successful tests could be performed in the high triaxiality regime, since such a stress state leads to fracture of the micro-structural bands oriented along the S-direction, thus resulting in the occurrence of splits on the fracture surface that are accompanied by a sudden force drop [36]. The occurrence of a split is a stochastic event, and the cohesive strength of the interface between different layers is not included in the continuum HE model. Additionally, splitting is associated with a sudden drop of the stress triaxiality in the center. Therefore, these tests were not considered for the calibration of the model parameters and, thus, are not included in Table 2. In total, 14 uncharged tensile tests, 12 tensile tests with C L,0 = 0.36 wppm, and 19 tests with C L,0 = 1.09 wppm were successful. Table 2. Numbers of tests for each geometry and nominal hydrogen concentration C L,0 , together with the average maximal normalized load (F/A 0 ) max , the average critical strain c , and the average embrittlement index EI . When applicable, the standard deviation is also reported.

Macro-Mechanical Responses
The resulting normalized load and true stress versus true strain curves are displayed in Figures 2 and 3, respectively. The figures show that for all four geometries, the presence of hydrogen did not affect the strength. However, hydrogen did affect the ductility, with a larger hydrogen concentration leading to reduced fracture strains.
A point of interest is the true strain at the maximum true stress, since it is associated with the macro-crack initiation at the center of the specimen [6]. After this point, the loadcarrying capacity rapidly drops, which can also be distinguished by a so-called 'knee-point' in the normalized load versus the plastic strain curve [38] (Figure 2). This point is further referred to as the critical strain c and is marked with a star in Figures 2 and 3.  An embrittlement index is commonly used for reflecting the amount of hydrogenassisted degradation: Data of the observed maximum normalized load, critical strain, and embrittlement index are reported in Table 2. Note that, due to the limited number of tests for R1.2, it was impossible to determine a standard deviation for that specimen geometry. Overall, increasing embrittlement can be observed for increasing nominal hydrogen concentrations C L,0 and increasing stress triaxiality, notwithstanding the exception of R1.2 with C L,0 = 0.36 wppm. In particular, hydrogen-charged smooth round bars (R∞) show significantly lower degradation compared to the notched round bars, as determined by one-tailed independent t-tests (p < 0.015 for all four combinations of R2 or R6 and for C L,0 = 0.36 wppm or C L,0 = 1.09 wppm).

Failure Mechanisms
The failure mechanisms of the tested API 5L X70 pipeline steel without and with hydrogen are described in another study [28] and will be summarized here. Illustrative fracture surfaces of the R2 and R∞ geometry with C L,0 = 0 wppm and C L,0 = 1.09 wppm are provided in Figure 4. The uncharged tensile tests feature regular ductile fracture by micro-void coalescence. In the hydrogen-charged specimens, a mix of ductile dimples and quasi-cleavage was observed. The presence of MnS inclusions was observed on both uncharged and hydrogen-charged fracture surfaces (Figure 4e-g), thereby implying a significant role in the fracture process. From the fracture surfaces ( Figure 4), it is clear that fractures in the lower stress triaxiality regime were increasingly dominated by shear failure. Moreover, the presence of hydrogen appeared to enhance the domination of shear-driven failure.
The damage processes were investigated in more detail using interrupted (notched) tensile tests and X-ray micro-CT [28] (The previous work [28] reported a hydrogen concentration of 0.89 wppm which corresponds to C L,0 = 1.09 wppm in the present work. The difference is explained by a different method for measuring the hydrogen concentration, since the previous work used hot extraction at 300°C, whereas the current work used hot extraction at 900°C). Stress triaxiality levels ranging between 0.4 and 1.3 (averaged over the entire load-deformation trajectory) were tested using various geometries. The fracture processes in uncharged and hydrogen-charged specimens were observed to be similar, in the sense that they both consisted of void nucleation occurring in planes along the micro-structural bands, followed by void growth and void coalescence. However, two major differences were observed, which could have contributed to the reported ductility loss. First, the presence of hydrogen caused void nucleation to occur at much lower strains. While the minimal true strain at which voids were observed in uncharged smooth round bar samples was around 0.7, voids in hydrogen-charged samples could already be observed at strains of around 0.1. Second, hydrogen enhances lateral void growth, which reduces the cross-sectional ligament. However, this was not associated with an acceleration in the void volume growth. Both the accelerated void nucleation, as well as the enhanced lateral void growth, were present for all stress triaxiality levels tested. The reader is referred to Ref. [28] for more information.

Hydrogen-Informed Gurson Model
The hydrogen-informed Gurson model is employed for modeling hydrogen-charged tensile tests. In this way, the ductile damage processes such as void nucleation, void growth, and void coalescence can be accelerated depending on the local lattice hydrogen concentration C L , as postulated by the HELP and HESIV mechanisms. The diffusion of hydrogen is typically modeled using the derivation by Krom et al. [39]: which includes the hydrogen diffusion towards the hydrostatic stress gradient ∇σ h and the division between the lattice hydrogen concentration C L and trapped hydrogen concentration C T (both in wppm). The lattice diffusion coefficient D L (in mm 2 /s) determines the rate of diffusion. Equation (4) also includes the equivalent plastic strain p and the partial molar volume of hydrogen (V H = 2 × 10 −6 m 3 /mol). The equilibrium between C L and C T was postulated by Oriani [40] and is dependent on the trap binding energy E B (in kJ/mol): with N L and N T being the density of metal atoms and trapping sites, respectively, β being the number of interstitial sites per metal atom, and α being the number of sites per trap. A finite element implementation in Abaqus using user subroutines is provided in Depraetere et al. [32].
The evolution of damage is described by the complete Gurson model (CGM) [41], where the void volume fraction f is the damage indicator. Each element has an initial assigned void volume fraction f 0 and is completely fractured when reaching the final void volume fraction f f . The yield criterion is given by the following: with σ being the stress tensor,σ being the flow stress, f * being the effective void volume fraction, σ e being the von Mises stress, and q 1 and q 2 being empirical constants. Note that the effective void volume fraction f * is equal to the real void volume fraction f before coalescence, and it is accelerated upon coalescence.
The evolution of f is governed by the sum of void nucleation and void growth. Void nucleation is modeled using a strain-controlled approach based on a Weibull distribution: with k > 0 being the shape parameter and N > 0 being the scale parameter of the Weibull distribution, f N being the volume fraction of void nucleating particles, and p being the equivalent plastic strain. Note that this is based on the well-known strain-controlled approach by Chu and Needleman [42], who used a Gaussian distribution. However, since a Weibull distribution produces strictly positive values, we regard this as a better choice for modeling void nucleation at reduced plastic strains [43]. For k > 1, the scale parameter N will approximate the mean value of the void nucleation distribution (both values differing by less than 10%) and serve a similar role as the mean nucleating strain N in the approach by Chu and Needleman [42]. The void growth rate is derived from the volume conservation using the trace of the plastic strain rate tensor˙ p kk [41]. Void coalescence is based on the limit load criterion by Thomason [41].
Considering the desire for a micro-mechanics-inspired model for HE predictions, the Gurson parameters should be adapted to reflect hydrogen-assisted degradation based on the experimental observations of the failure mechanisms. An important aspect is that void nucleation in hydrogen-charged samples occurred at significantly reduced plastic strains. Accordingly, the mean void nucleation parameter N was adapted for the presence of hydrogen. Data from several studies suggest that the degradation of the mechanical properties increases rapidly with increasing hydrogen concentration, after which saturation occurs [3,44]. For this reason, a degradation function was proposed that is exponentially decaying with increasing hydrogen concentration: where N,0 equals the scale parameter without any hydrogen, and ξ 1 and ξ 2 are degradation parameters describing HE. It can be easily derived that N,0 − ξ 1 represents the lowest possible value of N , while ξ 2 controls the speed of the degradation. Notably, Equation (8) only includes the lattice hydrogen concentration C L , which disregards degradation due to the trapped hydrogen concentration C T . This was done due to the lack of consensus as to whether numerical degradation models for hydrogen embrittlement should be based on C L , on C T , or on the sum of both [32,45]. Using Oriani's equilibrium (Equation (5)), it can be derived that C L typically predominates over C T [45]. Accordingly, the current study assumes that C L is mainly responsible for hydrogen-related degradation. Another micro-mechanical aspect observed in the experiments was the hydrogen enhanced lateral void growth. Since Gurson-Based models assume spherical void growth, void shape changes are neglected [5]. However, void shape effects may play a prominent role in the low stress triaxiality regime [6,46]. While extensions for including void shape and orientation changes exist [5,6,47], these models were deemed outside the scope of the present work.
Note that, when making the parameters N dependent on the local lattice hydrogen concentration C L , the nucleating distribution (Equation (7)) will change upon deformation due to hydrogen diffusion. For this reason, Equation (7) features the absolute nucleated voids f nucleation rather than the more common nucleation rateḟ nucleation . In this way, the number of voids nucleated at any point in time is a direct result of the instantaneous nucleating distribution. Additionally, it is specified thatḟ nucleation ≥ 0 such that nucleated voids cannot disappear.

Finite Element (FE) Models
Finite element (FE) models were constructed in Abaqus (2019) for simulating the (notched) uncharged and hydrogen-charged tensile tests ( Figure 5). A quarter-circular crosssection was modeled with appropriate symmetry boundary conditions. Three-dimensional linear brick elements were employed (C3D8T in Abaqus). The mesh at the expected fracture zone consisted of cubic elements with a length of 0.2 mm. Note that only a single notch was modeled in Figure 5b, whereas the experiments featured double-notched round bars (Figure 1). To best represent the localized deformation in a single notch in the experiments, the initial displacement rate applied in the notched geometries resulted in a local strain rate of 1.25 × 10 −4 s −1 and was doubled when the maximum force was reached. For the smooth round bar, a constant strain rate of 2.5 × 10 −4 s −1 was used.

R1.2 R2
⌀10 ⌀6 Given the predominance of C L over C T according to Equation (5), the specimens were given an initial lattice hydrogen concentration C L,0 corresponding to the nominal diffusible hydrogen concentration after charging (Table 2). At the edge of the specimen, the hydrogen flux J was assumed to be zero, such that hydrogen remained inside the specimen. It should be noted that this is an idealization, and that setting the edge at zero is an alternative boundary condition to represent the contact with the atmosphere during the mechanical test [32,43]. However, a comparison of these two extreme boundary conditions resulted in only a minor difference (+/−10%) of the macro-mechanical response for the employed strain rate and diffusion coefficient. In other words, the hydrogen concentration at the center, where failure is expected to start, was only marginally affected by the hydrogen concentration at the edge for the performed tests. The effect of the adopted boundary condition is described in Appendix A.
The hydrogen-related properties employed in the simulations can be found in Table 3. The diffusion coefficient D L was determined experimentally by permeation testing using a Devanathan-Stachursky setup, while the other parameters were adopted from literature [45,48]. The flow stress behavior of the API 5L X70 steel was determined as the true stress versus the true plastic strain obtained from a smooth round bar tensile test [43] and is presented in Figure 6. Note that this is an approximation, but it is frequently used for modeling ductile failure [49]. Table 3. Hydrogen-Related material properties used for the FE simulations.

Calibration Procedure
The hydrogen-informed Gurson model contains several parameters that require calibration. The empirical constants q 1 and q 2 were determined using the tabulated values as a function of the yield strength and hardening coefficient from the unit cell approach by Faleskog et al. [50]. The remaining parameters ( f 0 , f N , f f , k, N,0 , ξ 1 , and ξ 2 ) were determined through inverse analysis by comparing the simulation output with experimental data, as explained below. This is a common approach for calibrating ductile damage models [8,[51][52][53][54].
An objective function is defined based on the root mean squared error (RMSE) between the experimental and the numerical normalized load F/A 0 at equidistant levels of the true strain ln(A 0 /A) for all four geometries. Note that the FE model does not take the observed plastic anisotropy into account, but an engineering approach is followed by comparing the plastic strains.
To properly reflect the adapted fracture behavior in the presence of hydrogen, the parameters f 0 , f N , f f , k, and N,0 are first determined by minimizing the objective function using uncharged experimental data. Thereafter, using these calibrated parameters, the degradation parameters ξ 1 and ξ 2 are calibrated to best represent the hydrogen-charged tensile tests.

Agreement with Experiments
After minimizing the objective function, the calibrated model parameters were obtained (Table 4). Figure 7 shows the comparison between the performed experiments and the predicted simulations. Note that only the filled markers were employed for computing the objective function, since the open markers correspond to a stress-strain state after a delamination occurrence. In general, an excellent agreement was found for all four geometries and three hydrogen concentrations. Only at higher hydrogen concentrations did the model overestimate the fracture strains in the lower stress triaxiality regime. Table 4. Calibrated damage parameters for the investigated API 5L X70 steel. Similarly to the experiments, the 'kneepoint' can be extracted from simulations as the strain at the maximum true stress c . Using Equation (3), it can be transformed into an embrittlement index. Figure 8 compares the experimentally obtained embrittlement indices for each notch and geometry with the embrittlement indices from the simulations. In general, the model was capable of capturing the general trends of the hydrogen-assisted degradation quite well. Exceptions to this were the significantly overestimated fracture strains for the combination of the lowest considered stress triaxiality levels (R6, R∞) with the highest-considered hydrogen concentration (C L,0 = 1.09 wppm). We put forward two hypotheses to explain the poor agreement for the low triaxiality tests at higher hydrogen concentrations. A first hypothesis is the lack of an internal shearing coalescence criterion. The micromechanisms of void coalescence can be categorized into internal necking occurring at higher stress triaxialities and internal shearing at lower stress triaxialities [5,6,55,56]. Regarding the employed numerical damage model, Gurson-based models predicted the internal necking mechanism, and thus can yield poor numerical predictions when internal necking is not the dominant fracture mechanism [57][58][59][60]. Extensions, including the internal shearing coalescence, have been developed [47,61,62], but these typically require additional specific shear-dominated tests for calibration. Moreover, when modeling smooth round bar tensile tests with both coalescence mechanisms, it has been reported that internal shearing coalescence dominates [47,53]. Additionally, it is recognized that the presence of hydrogen can even change the mode from internal necking to internal shearing [29]. Micro-mechanical investigation of the fracture surfaces confirmed that failure was predominantly shear-driven in the low stress triaxiality regime and appeared to be enhanced by hydrogen (Figure 4). Accordingly, neglecting this increased shear-driven failure may lead to underestimated numerical predictions of the amount of embrittlement.

R1.2 R2 R6 R
An additional hypothesis for the poor agreement is the omission of void shape effects in the hydrogen-informed Gurson model. Hydrogen enhances void growth perpendicular to the loading axis, hereby decreasing the remaining ligament area compared to the spherical void growth. Furthermore, it has been reported that void shape effects may be particularly important in the low stress triaxiality regime [6,46]. Therefore, a model that does not account for these void shape effects may lead to underestimated numerical predictions of the amount of embrittlement in this regime.
Both aspects should be addressed if the aim is to accurately model fracture in the low stress triaxiality regime.
The optimal degradation function of the mean nucleating strain N (Equation (8)) is displayed in Figure 9. The mean nucleating strains for the nominal hydrogen concentrations tested were 0.31 and 0.1 for C L,0 = 0.36 wppm and C L,0 = 1.09 wppm respectively. Remarkably, N for C L,0 = 0 wppm and for C L,0 = 1.09 wppm agreed well with the minimal strain at which voids were observed (around 0.7 and 0.1, respectively) using the X-ray micro-CT method that was described above [28]. This suggests a strong connection between the experiments and the hydrogen-informed Gurson model from a micromechanical point of view, since the model appeared to quantitatively capture the void nucleation mechanism.

Fracture Locus
An alternative method to investigate the effect of hydrogen and stress triaxiality on the fracture behavior of the API 5L X70 steel is through the fracture locus. From the FE simulations with the calibrated parameters (Table 4), an average triaxiality η avg can be extracted from the central element: where i corresponds to the strain at which the central element coalesces according to Thomason's limit criterion. The average triaxiality can be used to construct the fracture locus c = f (η avg , C L,0 ). Such a fracture locus is relevant, since it can be used as an input for other HE model predictions [12,19,21]. Note that, after the critical strain c is reached, crack propagation occurs and the stress fields are disturbed. For this reason, the fracture locus was constructed using c , which is thought to be a better indicator of ductility than the fracture strain f [6,38]. Additionally, it is well established that the fracture behavior in the lower stress triaxiality regime is affected not only by the stress triaxiality, but also by the Lode angle parameter that can be derived from the third stress invariant [63]. However, since the current study features only axisymmetric tensile experiments, which all have an identical Lode angle parameter, this effect was not investigated. Figure 10 shows the resulting fracture locus. Note that the range of stress triaxiality values tested was roughly between 0.4 and 1.5. A common equation for expressing the fracture locus is as follows [27,[64][65][66]: The following was the best fit obtained using the averaged experimental data of each geometry from the uncharged tests and was added to Figure 10 To include the effect of hydrogen in the fracture locus, a reduced fracture locus c,H was considered. Two forms of degradation are presented in Table 5, which applied a subtraction and a multiplier to Equation (11). Both forms were fitted on the experimentally obtained c , which was averaged for each geometry and hydrogen concentration ( Figure 10). The goodness of both fits was evaluated through the mean squared error MSE, which is found in Table 5.  From Figure 10 and the MSE, it is clear that a subtraction of the fracture locus better represented the performed experiments than a multiplication. A multiplication of the fracture locus corresponds to the EI being independent of stress triaxiality, which was not the case for these performed experiments (Table 2), as well as for other experimental studies [24][25][26][27]. Remarkably, a multiplication has been commonly assumed for degrading the fracture locus in the presence of hydrogen, without any experimental evidence [21][22][23]. It must be noted that Yu et al. [31] studied hydrogen degradation through unit cell simulations and derived a hydrogen-degraded failure locus that was best represented by a multiplication. However, they only considered internal necking failure and triaxiality values greater than 1. Moreover, their conclusions might have been significantly impacted by their assumptions regarding the modeling of degradation. We advise that particular attention must be given to the best representation of the triaxiality dependence of the hydrogen-affected fracture strains of a specific material.

Conclusions
The effect of hydrogen on the stress-triaxiality dependent fracture behavior of an API 5L X70 steel was studied from a combined experimental and numerical viewpoint. For this purpose, four different specimen geometries were tested at three different hydrogen concentrations (0, 0.36, and 1.09 wppm).
Hydrogen-assisted degradation appeared to increase with increasing stress triaxiality. Based on the plasticity-dominated HE mechanisms HELP and HESIV, a hydrogen-informed Gurson model was proposed for modeling ductile damage. The model parameters were fitted based on macro-mechanical deformations, as well as micro-mechanical evidence of the failure mechanisms. By reducing the void nucleating strain with increasing hydrogen concentration, good agreement between the experiments and the micro-mechanics inspired model predictions was obtained. Furthermore, the calibrated void nucleating strain agreed well with the experimental observations.
The predictions of the model significantly overestimated the fracture strains of the low stress triaxiality specimens. This was attributed to a lack in the hydrogen-informed Gurson model of either an internal shearing coalescence criterion or void shape effects capturing enhanced lateral void growth. While these issues should be addressed in future work, the presented model captured the experimental trends well within the higher stress triaxiality regime.
Finally, the effect of hydrogen on the fracture locus was investigated. It was shown that the effect of hydrogen on the triaxiality-dependent fracture strains was better represented by a subtraction, rather than a multiplication, of the uncharged fracture locus.