Analysis of Inclined Cracks in Thin-Walled Circular Tube under Mixed-Mode I + II Fracture

This paper provides a study on mixed-mode fracture mechanics in thin-walled tube which is subjected to tension, shear and torsion loading. This type of loading causes an inclined crack to develop and generate a mixture of normal and shear stresses ahead of a crack tip. The stress state ahead of a crack tip is frequently based on mixed-mode type of interactions which designate the amplitude of the crack tip stresses. The analytical expressions for the stress intensity factors for mixed-mode I + II approach are presented. The Paris law for mixed-modes I + II has been discussed. Mixed-mode fracture mechanics is used with theoretical models to predict the path of crack growth when an inclined crack is subjected to a combination of mode I and mode II deformations. The torque at which crack propagation can be expected has been determined. The numerical calculations have been carried out by using MATLAB code. The results are good and could be useful for companies working with thin-walled circular tubes.


Introduction
Production and use of steel structures were found to be accompanied by failures and breakdowns of structures such as rails, vessels, large tanks, boilers, bridges and many more. Defects in structures may be the main cause of failures and breakdowns. Some problems in thin-walled circular tubes occur due to sudden brittle fractures in welded structures.
The stress field at the tip of the cracks should be studied in order to ensure safety of structures. A fracture in real metals is carried out by a process of deformation and local failure at the front of a slowly or rapidly propagating crack. Cracks in brittle structures are often governed by linear elastic fracture mechanics. Most of the previous studies on structures are focused on the mechanisms of crack initiation and propagation of Mode I cracks [1][2][3][4][5].
The most common understanding of Mode II failure is that shear mode crack growth is caused by a linkage of many small tensile cracks [6,7]. It is believed that at a mixed-mode I + II fracture, the crack grows in a direction that maximizes one of the stress intensity factors or , rather than in a direction that maximizes both stress intensity factors and simultaneously. The stress intensity factor is mostly used as a parameter to evaluate the state of stress near the crack tip according to Linear Elastic Fracture Mechanics (LEFM) approach. The inclined propagating cracks are analyzed for the direction of propagation of cracks. Numerical methods have many advantages to solving the fracture problems [8][9][10][11] and this is the approach of this study.

Mathematical Modeling for the Mixed-Mode I + II
In this study, a central inclined crack of length 2a in a thin-walled circular tube as shown in Figure 1a, is considered. When the tube is loaded biaxially, flaws will be subjected to mixed loading as shown in Figure 1b. The stresses are and = along the y and x directions, respectively. Mode I will result in opening of the crack faces while Mode II will cause a sliding motion of the crack faces along the length of the crack. The stress intensity factors and for the tensile and shear mode [12] are derived using Equations (1) and (2) as where is the hoop stress on the tube and the tensile stress and the shear stress acting on the crack are given in the form The longitudinal and hoop stresses in the circular tube are obtained from the equilibrium and are given as respectively. Equations (4) and (5) can be used to find = . The stress intensity factors for Mode I and Mode II in Equations (1) and (2) The tube is now loaded by a torque T and an axial tensile force F. In this case, it is assumed that r >> t so that the curvature of the tube may be neglected when the stress intensity factor is determined. The normal stress is given by the axial tensile force, F in the longitudinal direction and the shear stress is given by the torque T, as shown in Figure 2a. In Figure 2b, a thin-walled circular tube is subjected to a twisting moment T. The mean radius is given by , a typical small area of length and thickness t transmitting an increment of force = ( ) are assumed [14]. The moment of the incremental force about the axis of the circular tube equals the applied torque T. Therefore [14] The two-dimensional Cartesian stress transformation equations [15,16] are given as For the purpose of this study, the stress gives a normal stress ∅ across the crack which then loads the crack in Mode I. The maximum circumferential stress criterion [18], is based on the assumption that the crack will grow in the direction such that the Model I stress intensity factor is maximized. During a biaxial loading cycle, it is crucial to determine the fatigue crack growth direction. Applying the maximum tangential stress criterion, the crack extension direction can be given b [18] = 2 1 4 4 ± 1 4 4 + 8 where is the crack propagation direction. The equivalent stress intensity factors can be given in the form [17,19,20] ( ) = + + 4 + 1 .
Fracture is expected when = if LEFM can be used, therefore the value of the torque can be found as where is the critical crack length. Fracture toughness is a measure of the ability of a material to resist the growth of a pre-existing crack or flaw, that is, it is a critical value of the stress intensity factor at the time of an unstable crack propagation, that is, = .

Law of Fatigue Crack Propagation
In the 1960's, Paris [21] postulated that the range of stress intensity factor might characterize sub-critical crack growth under fatigue loading in the same way that ∆ characterized fast fracture. To analyze the mixed-mode fatigue crack, the Paris law [21] should be redefined by replacing ∆ with an equivalent stress intensity factor which takes into account both modes of fracture. Thus, where Equations (14)- (16) are used, C is the Paris coefficient and m is the slope stress ratio. The equivalent stress intensity factors in Equations (14)- (16) are selected for the calculation of the fatigue crack growth rate considering relative variation of the and .

Numerical Experimental Procedures
For this study, the crack propagation is the main parameter for the fatigue damage. A thinwalled circular tube is made from Ti alloy (Ti-6Al-4V Grade 5) [22] with of 107 MPa√ , yield strength of 1200 MPa, elastic modulus E = 119 GPa, the coefficient of Paris law of C = 6.0 × 10 the internal pressure p = 250 MPa producing a circumferential hoop stress of 360 MPa and the Poisson's ratio = 0.37. The tube has surface flaw oriented at angle ∅ = 60 to the uniform tensile hoop stress. The tube was regularly subjected to Nondestructive Testing inspection. The inspection showed that the tube contained an angled flaw of length 2a (a = 1.0 mm). The inner radius is 45 mm and the outside radius is 50 mm, while the wall-thickness of the tube is = − , and length l = 5 m.

Calculations and Results
The results were obtained using the Matlab script. Figure 3a shows the and values that have been estimated theoretically and it can be noticed that values are predominant over values. Therefore the crack is influenced more by stress intensity factor in Mode I. The torque from Equation (17) has been calculated and the value is 145.2 kNm at = 7 mm. The predominant Mode loading caused a crack path deviation by 0.83 when Equation (13) was used. The critical crack length was obtained to be 7.1987 × 10 m using Equation (18). In Figure 4c the equivalent stresses of Equations (14) and (15) are plotted against the crack length and it can be seen that there is a slight significance between them. In Figure 4d, it can be confirmed that the ( ) is predominant over both ( ) and ( ), while ( ) is predominant over ( ). All for A, B and C behave the same way as seen in Figure 4d. It is interesting that in Figure   4c, da/dN B dominates both da/dN A and da/dN C respectively.

Conclusions
Analysis of an inclined crack on thin-walled circular tube under mixed-mode I + II fracture has been studied in the present work. The crack propagation is the main parameter for the fatigue damage. The flaw is inclined at an angle to the direction of the hoop stress. So the hoop and longitudinal stresses are the most critical stresses. The numerical calculation of K for Mode I and Mode II were carried out. The results were obtained using the Matlab script.
It can be noticed that the values obtained for are predominant over those obtained for the . The major results of this study are summarized in the plots. The torque at which a crack propagation can be expected has been found and looks reasonable for this study. The angle for the crack extension direction has been determined. The results obtained in this study are good and could be useful for companies working with thin-walled circular tubes.