An Extended Finite Element Method (XFEM) Study on the Elastic T-Stress Evaluations for a Notch in a Pipe Steel Exposed to Internal Pressure

: The work investigates the importance of the K-T approach in the modelling of pressure cracked structures. T-stress is the constant in the second term of the Williams expression; it is often negligible, but recent literature has shown that there are cases where T-stress plays the role of opening the crack, also T-stress improves elastic modeling at the point of crack. In this research study, the most important effects of the T-stress are collected and analyzed. A numerical analysis was carried out by the extended ﬁnite element method (X-FEM) to analyze T-stress in an arc with external notch under internal pressure. The different stress method (SDM) is employed to calculate T-stress. Moreover, the inﬂuence of the geometry of the notch on the biaxiality is also examined. The biaxiality gave us a view on the initiation of the crack. The results are extended with a comparison to previous literature to validate the promising investigations.


Introduction
In the field of fracture mechanics especially linear elastics, the vicinity of notch tip is often symbolized by singular stress entities. Their resistance is measured through stress intensity factor (SIF). The SIF is mainly depended on the distance r from the tip of notch. In this field, parameter T is introduced to enrich the parameter K (SIF) to make the model better in the elastic stress field; this is the K-T approach [1].
In fracture mechanics body of knowledge, it is established that the same Stress Intensity Factor (SIF) is required for two cracks to propagate in the same way. Experiment [2] have shown that two plates with the same SIF and different crack length a 1 > a 2 show different the propagation of the cracks. The results have shown that the propagation speed of a 2 is higher than that of a 1 . The study concluded that the first term of asymptotic development is not sufficient to predict crack behavior. Therefore, it is necessary to increase the order. The first term asymptotic development is the SIF that determines the initiation and propagation of the crack. In addition, the second term is constant and controls the stability of the propagation direction. It is the transverse component symbolized by T.
Several studies have shown the significance of the T-Stress, and its influence on different parameters of mechanics. Jayadevan [3] highlighted that plastic zone is manipulated by the variation of T-stress. It means, the plastic area escalates with the increase in absolute value of T-stress and changes its form. Sobotka et al. [4] demonstrated the alteration in the plastic wake with many T-stresses which depends on plastic wakes height (HPW).
T-Stress has been an essential contributing factor in the stability of the direction of the crack propagation. For instance, T negative gives a stable direction, and for T positive it is unstable [5]. Fayed et al. [6] explored the impact of T-stress on propagating crack direction by Maximum Tangential Stress (MTS). The principle of MTS is that the crack propagates in the trend of maximum tangential stress. They obtained that the directions of the overall crack no coincide with the initial direction of the crack. Many studies have concluded the behavior by the fact that the tangential stress is affected by the T-Stress. MTS becomes generalized maximum tangential stress (GMTS), which considers the constraint T in the expression of stress. In the same context Shahani [1] studied the result on the initiating angle of propagating crack of the stress T. The study has shown a negative T value declines the angle of crack initiation, and a positive T value enhances it. Nejati et al. [7] gauged the relationship between T-stress and material properties. Chen et al. [8] has shown that Graded Poisson's ratio affects the T-stress. Additionally, Toshio et al. [9] concluded that the Poison's ratio influences the T-stress on a three-dimensional edge-cracked plate.
Other important research that has shown the influences of the T-stress includes: Zhang et al. [10], which used numerical manifold method (NMM) is employed to calculate the T-stress for two-dimensional functionally graded material (FGM) having numerous cracks. Noritaka et al. [11] has resulted the T-stress might open micro-branches in the mist region. For the bending and tension load, Hancock [12] determined that T decreases with escalating crack length. Matvienko [13] explored the influence of T-stress in problems of the elastic and the elastic-plastic fracture mechanics.
Conventionally, the T-stress is often calculated at the crack tip, which is certainly not the case in this research. The research evaluates the T-stress at the tip of notch through extended finite element method (X-FEM).
The finite element method FE method is limited by the simple cases, as well as the presence of a singularity greatly degrades the convergence of the FEM. Belytschko and Blacken in 1999 added discontinuously enriching function in finite element approximation by respecting boundary conditions. Later, Moës et al. [14] developed the technique and called it as extended finite element method which is abbreviated as X-FEM. The efficiency of the X-FEM is well endorsed in the literature. To simulate the propagation of cracks in porous media, Wang et al. [15] integrated embedded discrete fracture method (EDFM) with X-FEM simulating fracture associated fluid and solid mechanics. Shu et al. [16] investigated the fatigue growth of 3-D multiple cracks by X-FEM. The implementation of X-FEM for composites resulted successfully [17][18][19]. X-FEM is also used for the calculation and analysis of failure mechanics parameters. Fakkoussi et al. [20] calculated stress intensity factor for mode one by X-FEM. Llavori et al. [21] studied the problems of contact fatigue by X-FEM.
The research study presents the use of X-FEM to calculate T-stress in the notch tip for an arc of the pipe of steel P264GH. Further, it demonstrates the benefits of using the X-FEM approach to compute K-T at the notch point in an arc under pression and revealing the possibility of detecting the crack initiation.
The remaining article is organized as, Section 2 talks about the K-T approach, X-FEM, and the geometry used. Section 3 deals with the numerical result obtained, compared against FEM conclusions, along with the discussion. Finally, the conclusion is in the other section.

K-T Approach
Stress Intensity Factor has been an essential parameter in the field of linear fracture mechanics, widely used for crack evaluation. SIF measures the strength of the singularity, and integrates various parameters such as load, geometry, and shape of crack. M. Hadj Meliani [22] suggested that a thin structure such as a thin pipe with a longitudinal crack, it is very difficult to characterize the stress field by a single parameter. For linear elasticity, the enrichment of the SIF with the T-stress is required to model the notch tip. In literature, studies show that T explains how geometry influences tenacity (K IC ) [22]. T-stress helps in approximating the level of stress at a crack or tip of the notch. The possibility of constructing the K(T) curve numerically has given the opportunity to predict the loading of a crack initiation [23,24]. Including T-stress in calculations, it improves the prediction of propagating crack under the control of mixed loading. The K-T was additionally calculated for through-wall-cracked pipes under various pressure conditions by three dimension-3D FE [25]. The importance of the T-Stress is established in many works.
• T-stress enhances the possibility of crack opening stresses in the context of small crack [26]; Taking σ xx (r, θ) = σ cr , for the crack propagation, the first term tends to K 1 √ 2πr f xx (θ) to a because K 1 = S yy √ πa, and if a tends to zero; the first term becomes negligible intheface to T at the crack point. lim In this case, the T-stress cannot be ignored, T play the role of crack opening, and therefore the importance of T varies with the size of the cracks.
• T-stress influences the plastic zone. Jayadevan [3] highlighted that plastic zone is affected by the variation of T-stress. The plastic area escalates with the upsurge in the absolute value of T-stress and changes its form. Sobotka et al. [4] demonstrated the deviation in the plastic wake with various T-stress which depends on plastic wakes height (hpw). • Propagation's direction: Fayed et al. [6] analyzed the impact of T-stress on the propagation's direction by Maximum Tangential Stress (MTS). The principle of MTS is in the crack propagation direction which is in line with that of maximum tangential stress. MTS only considers the term singularity. Therefore, the direction of crack does not coincide with the initial directions of the crack. This behavior shows that the tangential stress is affected by the T-Stress. MTS becomes generalized maximum tangential stress (GMTS), which considers the T-stress in the expression of stress. • T-stress has an impact on crack initiation angle. Shahani [1] analyzed the consequence of the T-stress on the angle of initiation of crack propagation. The study has shown that a negative T value declines the angle of crack initiation, and a positive T value enhances it. • T-Stress is of prime significance when ensuring the stability of the direction of crack propagation, such as T negative gives a stable direction, and for T positive it is unstable [5].
T-Stress could be computed through numerous techniques. Weight Function Method has shown its efficiency in several problems cracking-related such as edge-cracked rectangular plate, circular disk [27]. Kfouri [28] developed a technique for evaluating the T-stress. The method uses the attributes of the path-independent J-integral and is called the Esheby-Integral method. The stress different method (SDM) has been proposed by Yang [29]. The idea of this method is the errors of the numerical values of σ 11 and σ 22 near a crack point progress with r in the same way, and the variation must effectively eliminate errors.
Biaxiality is a parameter that relates the SIF and T-stress:

Extended Finite Elements
One of the uses of the FE method is the study of crack propagation but is limited for simple cases. If the mesh size does not conform to the crack, the FE method does not treat the propagation, and the presence of a singularity degrades the convergence of the FE method.
The solution is to add enrichment function to the FE approximation (see Figure 1); this is the extended finite element method.
where: Mathematics 2021, 9, x FOR PEER REVIEW edge-cracked rectangular plate, circular disk [27]. Kfouri [28] developed a techniq evaluating the T-stress. The method uses the attributes of the path-independent J-in and is called the Esheby-Integral method. The stress different method (SDM) has been proposed by Yang [29]. The idea method is the errors of the numerical values of σ11 and σ22 near a crack point pr with r in the same way, and the variation must effectively eliminate errors.
Biaxiality is a parameter that relates the SIF and T-stress:

Extended Finite Elements
One of the uses of the FE method is the study of crack propagation but is limit simple cases. If the mesh size does not conform to the crack, the FE method does no the propagation, and the presence of a singularity degrades the convergence of t method.
The solution is to add enrichment function to the FE approximation (see Figu this is the extended finite element method.  Step of enrichment methods. The X-FEM requires operation to confirm the enrichment status of a knot acco to its position in reference to the crack and to evaluate the functions H(x) and F(x position (r, θ) in relation to the notch point is calculated herein to know if x is ab below the crack. These operations are carried out using the level sets method. The nique for describing crack is known as level set method. In the X-FEM, it determin location of the crack and the crack point, and the position to apply discontinuo richment and enrichment to the crack point figure1. Most importantly, the level se vides an instant result that helps track crack propagation, i.e., as the crack propa enrichment at the crack front becomes discontinuous enrichment, and nodes (n riched) become enriched. There are two-level functions, and the first describes the surface (φ), second gives the crack front (ψ) [30].
The X-FEM method was applied in several studies. Yousheng Xie et al. [31 implemented the X-FEM method in the study of propagating crack in mixed mod The X-FEM requires operation to confirm the enrichment status of a knot according to its position in reference to the crack and to evaluate the functions H(x) and F(x). The position (r, θ) in relation to the notch point is calculated herein to know if x is above or below the crack. These operations are carried out using the level sets method. The technique for describing crack is known as level set method. In the X-FEM, it determines the location of the crack and the crack point, and the position to apply discontinuous enrichment and enrichment to the crack point Figure 1. Most importantly, the level set provides an instant result that helps track crack propagation, i.e., as the crack propagates, enrichment at the crack front becomes discontinuous enrichment, and nodes (not enriched) become enriched. There are two-level functions, and the first describes the crack surface (ϕ), second gives the crack front (ψ) [30].
The X-FEM method was applied in several studies. Yousheng Xie et al. [31] have implemented the X-FEM method in the study of propagating crack in mixed mode, and evaluated the crack initiation angle. Reference [32] has shown the performance of the method. The study applied X-FEM to calculate SIF for 3D crack propagation problems for a Compact Tension C-T specimen. X-FEM was also implemented in the analysis of bi-material interfaces, calculating service life and fatigue resistance [33]. An integration between the X-FEM and embedded discrete fracture method (EDFM) is established for simulation of the process of fluid fracture propagation in porous media [15]. Savenkov et al. [34] employed the X-FEM to represent the central surface of the crack. The application of X-FEM for composite models is also carried out supporting current investigation [35,36].
X-FEM was used to predict components failure from a different form of notch [37]. Patria et al. [38] adopted X-FEM to study the mechanical attributes and fracture behavior of (Reinforced Polymeric Composites) RPC materials with single edge notch three-point bending.

Geometry
An arc of pipe containing a notch under pression was numerically analyzed using X-FEM in ABAQUS software. The material used is a steel P264GH. The arc characterized by an inner radius Ri = 219.55 mm and thickness t = 6.1 mm. More details on geometry, the shape of the notch, and boundary requirements used are illustrated in Figure 2 and Table 1.
a Compact Tension C-T specimen. X-FEM was also implemented in the analysis of bi-material interfaces, calculating service life and fatigue resistance [33]. An integration between the X-FEM and embedded discrete fracture method (EDFM) is established for simulation of the process of fluid fracture propagation in porous media [15]. Savenkov et al. [34] employed the X-FEM to represent the central surface of the crack. The application of X-FEM for composite models is also carried out supporting current investigation [35,36].
X-FEM was used to predict components failure from a different form of notch [37]. Patria et al. [38] adopted X-FEM to study the mechanical attributes and fracture behavior of (Reinforced Polymeric Composites) RPC materials with single edge notch three-point bending.

Geometry
An arc of pipe containing a notch under pression was numerically analyzed using X-FEM in ABAQUS software. The material used is a steel P264GH. The arc characterized by an inner radius Ri = 219.55 mm and thickness t = 6.1 mm. More details on geometry, the shape of the notch, and boundary requirements used are illustrated in Figure 2 and Table 1.  The mechanical attributes of the material used are presented in Table 2, and the chemical composition of the material are included in Table 3.    The mechanical attributes of the material used are presented in Table 2, and the chemical composition of the material are included in Table 3.

Results
This section presented the results of SIF and T-stress given by X-FEM via an user element UEL subroutine, the calculation was executed by ABAQUS software, we used quadratic element C3D20 in the mesh Figure 3, with a size of 0.5 mm. The number of elements is 246117. The number of nodes is 363256.

Results
This section presented the results of SIF and T-stress given by X-FEM via an user element UEL subroutine, the calculation was executed by ABAQUS software, we used quadratic element C3D20 in the mesh Figure 3, with a size of 0.5 mm. The number of elements is 246117. The number of nodes is 363256. The T-stress is calculated through stress different method, and normalization is done for effect of the T-stress relative to the stress intensity factor by a parameter dimensionless termed biaxiality.

= T√πa
K is the value of stress intensity factor, and a is the notch length. Figure 4 illustrates the difference of the SIF in mode 1 as a function of r, for a/t = 0.2, by extended finite element. The elastic SIF distribution at the notch tip decreases with distance from the tip of notch, for a/t = 0.2 the maximum value of SIF is 21 MPa√m at notch tip, i.e., = 0 (see Figure 4). Near the notch tip SIF decreases rapidly to 6 MPa√m at r = 0.3 mm, then its variation becomes slower. Figure 5 shows the variation of stress σ xx, σ yy and T. T-stress increases with increasing r up to r = 0.43, and after that it starts to stabilize. The numerical cal- The T-stress is calculated through stress different method, and normalization is done for effect of the T-stress relative to the stress intensity factor by a parameter dimensionless termed biaxiality. β = T √ πa K K is the value of stress intensity factor, and a is the notch length. Figure 4 illustrates the difference of the SIF in mode 1 as a function of r, for a/t = 0.2, by extended finite element.

Results
This section presented the results of SIF and T-stress given by X-FEM via an user element UEL subroutine, the calculation was executed by ABAQUS software, we used quadratic element C3D20 in the mesh Figure 3, with a size of 0.5 mm. The number of elements is 246117. The number of nodes is 363256. The T-stress is calculated through stress different method, and normalization is done for effect of the T-stress relative to the stress intensity factor by a parameter dimensionless termed biaxiality.

= T√πa
K is the value of stress intensity factor, and a is the notch length. Figure 4 illustrates the difference of the SIF in mode 1 as a function of r, for a/t = 0.2, by extended finite element. The elastic SIF distribution at the notch tip decreases with distance from the tip of notch, for a/t = 0.2 the maximum value of SIF is 21 MPa√m at notch tip, i.e., = 0 (see Figure 4). Near the notch tip SIF decreases rapidly to 6 MPa√m at r = 0.3 mm, then its variation becomes slower. Figure 5 shows the variation of stress σ xx, σ yy and T. T-stress increases with increasing r up to r = 0.43, and after that it starts to stabilize. The numerical cal-  The elastic SIF distribution at the notch tip decreases with distance from the tip of notch, for a/t = 0.2 the maximum value of SIF is 21 MPa √ m at notch tip, i.e., r = 0 (see Figure 4). Near the notch tip SIF decreases rapidly to 6 MPa √ m at r = 0.3 mm, then its variation becomes slower. Figure 5 shows the variation of stress σ xx, σ yy and T. T-stress increases with increasing r up to r = 0.43, and after that it starts to stabilize. The numerical calculation of the stress σ xx, σ yy by X-FEM is executed by ABAQUS software. Figure 6 gives the Von Mises stress obtained by the Abaqus software. culation of the stress σ xx, σ yy by X-FEM is executed by ABAQUS software. Figure  6 gives the Von Mises stress obtained by the Abaqus software.  The result obtained of biaxiality are compared with H. Moustabchir [39], who calculated the biaxiality through the finite element method, for the same geometry which is used in this study. Figure 7 gives the variation of biaxiality as a function of a/t, by X-FEM and FEM.
The biaxiality levels up with the increase a/t, The results are identical with the H. Moustabchir [39] recommendations. Kim and Paulino [40] has also got the same variation. By X-FEM the biaxiality varies from −0.21 to 0.76 for a/t = 0.2 and 0.8, respectively, and it goes from a positive to a negative value at a/t = 0.3. The difference between the results given by X-FEM and FEM is 0.013. culation of the stress σ xx, σ yy by X-FEM is executed by ABAQUS software. Figure  6 gives the Von Mises stress obtained by the Abaqus software.  The result obtained of biaxiality are compared with H. Moustabchir [39], who calculated the biaxiality through the finite element method, for the same geometry which is used in this study. Figure 7 gives the variation of biaxiality as a function of a/t, by X-FEM and FEM.
The biaxiality levels up with the increase a/t, The results are identical with the H. Moustabchir [39] recommendations. Kim and Paulino [40] has also got the same variation. By X-FEM the biaxiality varies from −0.21 to 0.76 for a/t = 0.2 and 0.8, respectively, and it goes from a positive to a negative value at a/t = 0.3. The difference between the results given by X-FEM and FEM is 0.013. The result obtained of biaxiality are compared with H. Moustabchir [39], who calculated the biaxiality through the finite element method, for the same geometry which is used in this study. Figure 7 gives the variation of biaxiality as a function of a/t, by X-FEM and FEM.  The biaxiality levels up with the increase a/t, The results are identical with the H. Moustabchir [39] recommendations. Kim and Paulino [40] has also got the same variation. By X-FEM the biaxiality varies from −0.21 to 0.76 for a/t = 0.2 and 0.8, respectively, and it goes from a positive to a negative value at a/t = 0.3. The difference between the results given by X-FEM and FEM is 0.013.

Discussion
SIF measures the strength of the singularity, which explains the rapid variation of the SIF obtained for r = [0, 0.3]. The more approach is made towards the notch tip, the more the stress concentration increases, and therefore, SIF increases. Moustabchir et al. [39] used the volumetric approach to calculate SIF, in the same condition that this research studied. Moustabchir et al. obtained for mode I at tip notch K 1 = 21.6 MPa √ m, which differs from our result by 0.5. The T-stress can be analyzed by the SDM. The biaxiality increases with the rise in a/t, and the notch depth affects the value of T-stress. The same variation was obtained in other investigations on various materials. Bouchard et al. [41] have shown that T increases with increasing depth for a mono silicon. In [42], Sherry et al. obtained an increase in T in absolute value with the variation of the crack size over the width of a plate. In addition, Ayatollahi et al. [43] obtained for mode I, an increase in T-stress as a function of the depth of the crack for a single edge notched. If Which is not the case for this study β max = 0.78, so T = σ, which implies that T has no influence on the notch in our case and our condition.
Many studies have resulted that the influence of the T-stress is remarkable and significant when T is negative [1,3,6], however, in this study for r < 0.43, negative T-stress causes an increase in the plastic zone [27]. This will cause a crack to initiate. In the presence of a crack, negative T can change the direction of propagating crack and decreases the crack growth initiation angle [31].
The K-T approach is an integration between the SIF and the T-stress to improve modelling the elastic stress at the point of the crack. The importance of T-stress has been highlighted, and it takes the place of short crack opening stress. Besides, the importance of the cooperation of SIF and T-stress, such as the K(T) curve was elaborated which gives a prediction of the stress of crack initiation [44]. Neggaz et al. [44] studied the influences of the reinforcements in the structure of composites, with the aim of reducing constraints at notch-tip. Moreover, authors evaluated the effective stress intensity factors in the regard of propagating crack in thin and thick panels. Therefore, an Extended Finite Element Method (XFEM) is novel and improved technique on the elastic T-stress evaluations for a notch in a pipe steel exposed to internal pressure.

Conclusions
Three-dimensional Extended Finite Element (X-FEM) analysis is applied to evaluate the stress intensity factor and the T-stress for an arc of pipe with external notch under internal pressure. The results are presented below:

•
To study the influence of geometry and notch size on the T-stress, authors have approached the biaxiality as a function of a/t. The evaluation endorsed that biaxiality β increases with the increasing a/t which is in accord with president results.

•
The integration of biaxiality allowed us to determine the state of the crack initiation, and we can say that the pipe is safe in the used conditions.

•
The SIF alone does not characterize the behavior of notches. T-stress is obtained by Stress Difference Method (SDM) along with the notch in mode I. SDM is an efficient and simple method to calculate the fracture parameters. • With ABAQUS-based investigations, the numerical results achieved by X-FEM are in good agreement with the Moustabchir result [39]. The implementation of X-FEM in the presence of a notch corrected the problems of the standard finite element method. The advantage of the X-FEM is that the mesh is independent of the notch.
For more precision, the future objective is to calculate the parameters of fracture mechanics by iso-geometrical analysis.

Data Availability Statement:
The data that support the findings of this study are available on request from the corresponding author.