Numerical study of semi-elliptical cracks in the critical position of pipe elbow

semi-elliptical crack angle orientation was studied at the critical position to estimate the critical angle. In the last part, the failure assessment diagram (FAD) was used to show the critical crack depth ratios at a critical position and a critical angle.

semi-elliptical crack angle orientation was studied at the critical position to estimate the critical angle. In the last part, the failure assessment diagram (FAD) was used to show the critical crack depth ratios at a critical position and a critical angle.

INTRODUCTION
ipelines are very important tool to transport hydrocarbons with high safety properties over long distances [1][2][3]. Indeed, the piping systems including the elbows are subjected to severe factors which may menace the integrity and safety of the pipe [4][5]. Due to the importance of elbows in the petroleum industries including their geometry and location, they are considered to be as critical parts in piping systems [6][7]. Due to their geometry, stress amplification occurs which promote failure and leak. Moreover, their stiffness is less in comparison with straight pipes having the same crosssection and material properties, inducing significantly higher stresses and deformations which may lead to failure. The corrosion phenomenon of elbow is one of the major damage problems in the energy transportation due to unsuitable factors such as fluid nature, operating conditions and the material of pipeline [8][9][10][11]. This situation may affect the energy transport significantly and increase the possibility of corrosion to occur. Many authors have studied the interaction between fluid and solid or solid-solid for two types of damage such as erosion and corrosion of steel pipelines [12][13][14][15]. El-Gammal et al. [16] have presented the flow accelerated corrosion of a 90-degree elbow by computational fluid dynamics (CFD) simulation. They have proved that the wear maximum value at elbow extrados is approximately 37% less than the maximum value along intrados. Muthanna et al. [17] have investigated the erosion corrosion in the internal surface of pipe elbow. They have found that the small sand particles and the attack of fluid flow on critical positions were responsible for this serious problem. Tian et al. [18] have examined the effect of flow velocity on the corrosion behavior of AZ91D magnesium alloy at pipe elbow to analyze the interaction between fluid and internal wall of elbow in the critical zones. Their results have showed that the corrosion rate augment with the increasing of flow velocity. The erosion-corrosion phenomenon occurs in difficult situations where the interaction is a complex field which leads to make a degradation for wall thickness of pipeline steels [19][20][21]. Barros et al. [22] have studied the repairing of a real corrosion defect by composite sleeve at welded joints. Moreover, many researchers [23][24][25][26] have studied the circumferential stress in the pipe elbow by analytical models, numerical methods and/or experimental tests. They proved that the maximal circumferential stress located in the intrados section of elbow. In addition, other authors [27][28][29][30][31] have used the Failure Assessment Diagram (FAD) to classify the grade of safety of pipelines or elbows. This work aimed to study the safety of a pipe elbow under the effect of gas pressure with change of radius bending, crack position, crack angle orientation, and semi elliptical crack depth ratios. To evaluate the influence of the flow in the internal wall of the critical position of a pipe elbow, a semi-elliptic crack was used. The results of the stress intensity factor (mode I) can give the solutions to avoid corrosion phenomena due to the nature of liquid and the material quality of pipe elbow. This study was conducted to improve the durability and the performance of elbow steel.

NUMERICAL STUDY AND SOFTWARE ANALYSIS
orrosion and erosion induce the most common defects occurring to pipeline. An example for a corrosion case of a pipeline is shown in Fig. 1. We have performed a numerical study in order to understand, relief and avoid the occurrence of such corrosion problems as much as possible. In this numerical study, the interaction between elbow and natural gas analyses was carried out by ANSYS software [32]. 90° pipe elbow made of API X52 steel was investigated and the pipe and its dimensions are shown in Fig. 1. Tab. 1 and 2 present the chemical composition and mechanical properties of API 5L X52 pipeline steel, respectively. Tab. 3 shows the parameters of the natural gas flow properties. The geometry of the elbow is shown in Fig. 1b. Dimensions are as follow: internal radius R i =298.45 mm, wall thickness t = 12.7 mm and the length L=1000 mm. is the elbow bending radius radius.

RESULTS AND DISCUSSIONS
ue to the complexities of carrying out an experimental study, we have used a numerical software in order to understand the interaction between elbow and natural gas. In particular, we have used ANSYS Workbench software with the computational fluid dynamics Software (CFD section). A steady state gas flow analysis and pressure distributions were carried out. These pressures across the elbow were obtained by the importation of the other section called APDL (Ansys Parametric Design Language) where the Fig. 2 illustrates the methodology of numerical steps. First of all, elbow radius was applied progressively in order to show the stress. At stabilization, the study of stress distributions along elbow curvature under gas pressure of 7MPa was conducted in order to get the critical position defined as the position of maximum stress. At this critical position, a semi-elliptical crack was introduced and its orientation  effect was studied in order to get the critical crack angle. This will create different semi-elliptical crack depths (a/t) in the range of 0.1 to 0.8; where "a" is the crack depth and "t" is the pipe elbow thickness. Failure assessment diagram (FAD) curve was used to evaluate the critical crack depth ratios at a critical position and a critical angle.

Effect of elbow radius
In order to find out the adequacy of different finite element types and the sensitivity of our models to different parameters, different steps were studied by extracting their critical positions and applying boundary conditions with suitable mesh to get better results. Fig. 3 illustrates the elbow geometry, meshing and boundaries conditions of our model. Pressure distributions were applied in the internal surface of pipe elbow. In order to analyze the effect of the cracked elbow on the safety of piping system at critical zones, the radius bending was taken into our consideration as it plays an important role in the security and assessment for pipeline failure where the critical zones are affected by the interaction fluid/ structure including elbow geometry. Figs. (4)(5)(6) show the distribution of pressure, velocity and Von-Mises stress along the pipe elbows, respectively. The Effect of elbow radius was investigated in order to obtain the radius value at which the stress is stabilized. The change in the elbow radius leads to change the behavior flow including the maximum stress and therefore to security and defect assessment of pipeline as shown in Fig. 7.  The effect of elbow radius was also examined according to the maximal values of Von-Mises stress as it shown in Fig. 7. Von Mises stress evolution can be divided into four zones as indicated in Tab. 4. The admissible elbow radius refers to the value of Yield strength of API X52 pipeline steel which is equal to 410 MPa and the traditional safety factor taken as F s = 2. Any bending radius leading to a maximum stress greater than 410 MPa divided by the safety factor is not acceptable for industrial service. As a real case (Fig. 1a), a value of 500 mm of elbow radius was taken in our investigation. For this case: the stress amplification due to the presence of elbow is k t = 1.2.

Critical position along elbow
ANSYS simulation offers a good way to increase the pipeline efficiency and to prevent the occurrence of corrosion. The presence of high stress in elbow can be considered as a crucial problem depending on elbow angle. Critical zone presented as a red zone in Fig. 6 that is resulting from the internal natural gas pressure. Fig. 8 presents the corresponding distribution of different stresses (Von-Mises, normal stress and principal stress). Fig. 8 indicates that the maximum value of Von-Mises stress is 198.96 MPa. Normal stress  xx and principal stress   have the same value equal to 219.56 MPa. The maximum stress criterion apply to the 3 different stresses gives a maximum stress localization in the same node with an angle α=72°. In the next part, a semi-elliptical crack is introduced with in different angular position along the intrados of elbow (=5° to = 85°); is defined in Fig. 9. The semi-elliptical crack has a relative depth ratio a/t=0.5 and a/c=0.5. In the next figures, the stress intensity factor in mode I as a function of crack tip angleis presented. The original point of cracking at critical zone is one of among important reasons responsible for pipe elbow corrosion issue. Fig. 10 presents the stress intensity factor distributions in a semi-elliptical surface crack for different angular positions along the elbow intrados surface. Maximum mode I stress intensity factor KI,max is localized at =75° and is given by the following formulae [33][34][35][36] : " max " is the maximum applied stress, "a" is the crack depth and "F" is the geometric function which depends of crack length, aspect ratio, as well as the type of applied load.

Effect of angle crack orientation on SIFs
Crack is not naturally oriented perpendicular to normal stress therefore mixed mode of loading (KI, KII, KIII) is applied. In addition, elbow curvatures induce mixed mode of loading however crack orientation is  = 0°. Stress distributions at semi-elliptical cracks tip are illustrated in Fig. 11

Effect of crack depth ratio on SIF
As can be seen from SIF results, the crack is submitted to mixed mode of loading. Although the opening mode (I) is the most significant mode, still modes II and III are not negligible. Mode II is induced by longitudinal stress and mode III by pipe and elbow curvatures. Opening mode is at maximum for  =90°. This angle point can be named as the critical angle. For this angle point, KII is negligible and KIII is less than KI of about one forth. Therefore, failure criterion is given in function of equivalent stress intensity factor which is given by: where  is the Poisson's ratio. Keq for crack orientation angle = 90° is higher than KI, this indicates again that pure mode I is not obtained. Fig. 13 presents the equivalent stress intensity factor with the orientation  of semi-elliptical crack in different angular position It shows that orientation of 90° has more efficiency for elbow pipeline. The influence of relative depth ratios (a/t) changing from 0.1 to 0.8 on stress intensity factors (K I , K II , K III ) is shown in Fig. 14 for a/t = 0.5. Figure 13: Equivalent stress intensity factor versus relative crack orientation for different angular position (Elbow: internal radius Ri = 298.45 mm, wall thickness t = 12.7 mm, elbow radiusmm, P = 7 MPa, a/t =0.5).
Due to the concentration of energy in the critical zones in the elbow, stress intensity factors in mode K I and K eq increase with relative crack depth ratios (a/t) as seen in Fig. 13 and 14.
Figure14: Stress intensity factor in mode KI versus angular position for different relative crack depth ratios (a/t) The maximum equivalent stress intensity factor K eq was obtained for  = 90° while K II is approximately negligible. Semielliptical crack with different relative depths was examined at a critical position ( = 75°) and a critical crack angle (θ = 90°). A failure criterion for mixed mode of loading was used: Failure condition is given by C=1, K I , K II and K III are applied in elbow at critical zone for mode I, II and III K IC , K IIC and KIIIC are the fracture toughness of three modes equal to 116.6, 116.6 and 67.27 MPa √ m respectively for API X52 pipe steel. In steel; And where the equation of calculation K IIIC is given by as follow where c and c are critical shear and normal stresses. Fig. 15 gives an evolution of fracture condition C for different crack orientations under the same internal pressure and for a relative crack depth a/t =0.5. Fracture condition is fulfilled for crack orientation in the value range of θ from 32.52° to87.66° for semi-elliptical crack front position (

PARTICULAR FAILURE ASSESSMENT DIAGRAM (FAD)USE FOR ELBOW CRACK DEFECT
n order to analyze the behavior of elbows failure the task of Failure Assessment Diagram was focused on the determination of stress intensity factors which was examined depending the critical crack depth ratios at a critical position and a critical angle. In the case that the stress intensity factors was not enough to show the critical zones along the propagation of the cracks in the elbow, the FAD becomes a necessary step to predict the crack zones depending on three zones: the security zone, the safety zone, and the failure zone. The assessment of elbow defect harmfulness is made using Failure Assessment Diagram (FAD). FAD presents defect loading conditions in a graph where the non-dimensional crack driving force kr is plotted versus non-dimensional applied stress Lr. Due to a mixed mode loading, the non-dimensional crack driving force kr is defined as the ratio of maximum equivalent stress intensity factor, K eq,max , to the fracture toughness of material in mode K IC . kr=Keq/KIc (6) Here in this particular FAD Non-dimensional load Lr is defined as the ratio of maximum circumferential stress max and flow stress  o :  Lr=θ θ,max/0 with  0 = y + ul /2 (8)  y is the yield stress, ult the ultimate stress. In FAD, failure is given by the following condition; if the assessment point of coordinate * * [ , ] r r L k is under the failure curve given by the following equation , = , ) where the subscript c indicates critical condition , the structure is safe. If the assessment point is above the curve, failure occurs. Fig. 16 represents the Failure Assessment Diagram for elbow and straight pipe for different value of relative depth crack ratio (a/t) in the range 0.1 -0.8. Values of fracture toughness take into account the constraint effect [29,37]. Constraint is defined by T stress [38] and determined by the Stress Difference Method [39,40]. Material Failure Master Curve (MFMC) expressed the dependence of fracture toughness with effective critical constrain , . For pipe steel, MFMC has been determined by Eqn. 9 where it obeys a linear relationship where a is equal to -0.069√m and KIC,0is the reference fracture toughness for Tef,c=0 andKIC,0=77.28 MPa√m. For a straight pipe, a value of K IC =116.6Mpa√m was obtained and corresponds to a low constrain which has generally encountered in a pipe submitted to internal pressure. For elbow, the value of fracture toughness was reduced to KIC,0=95 Mpa√m due to a higher constrain. Evolutions of assessment points for different relative crack depth ratios are reported in FAD for straight pipe and elbow with same pressure, diameter and thickness in Fig. 16. It is noted that due to stress amplification and higher constrain, the assessment points for the same defect size is located with higher ordinates. Critical defect sizes are defined when evolution curves cross the failure curve. For elbow, the relative critical defect depth is (a⁄t=0.28), which is higher for a straight pipe (a⁄t=0.66). Assessment points for elbow and below failure are located in the brittle zone of the FAD. This is proving the use of a linear behavior for computing SIF. This assumption is not valid for a straight pipe but it is given here for comparison.

CONCLUSION
he presence of elbow in piping system leads to activate a locally stresses. For the studied pipe (internal radius Ri = 298.45 mm, wall thickness t = 12.7 mm) and for internal pressure P = 7 MPa, stress amplification factor kt was found to be 1.2. The maximum stress was located at a curvature angle α=72°. This localized an eventual defect which is considered as a semi-0elliptical surface defect. This defect is submitted to mixed mode of loading due to longitudinal and circumferential pipe stress distribution and elbow and pipe curvature. The mixed mode consists of superposition of mode I, II and III. Therefore, an equivalent stress intensity factor was computed and a mixed mode failure criterion was used to determine critical semi elliptical crack orientation. The assessment of defects in elbow can be made with a modified Failure Assessment Diagram where the non-dimensional crack driving force refers to the equivalent stress intensity factor and the value of fracture toughness which take into account high constrain in elbow. Severe stress conditions of a defect in elbow locate its assessment point in the brittle part of the FAD.