Empirical Solution of Stress Intensity Factors for the Inclined Inner Surface Crack of Pipe under External Pressure and Axial Compression

Abstract

Based on fracture mechanics theory, a finite element method was used to determine the stress intensity factors of the inclined crack on the inner surface of the pipe under axial compression load and external pressure. The effects of different influencing factors on the stress intensity factor along the crack front considering crack closure were systematically explored, which were different to those under internal pressure. The effects of high aspect ratio on K_II, the crack inclination asymmetry caused by curvature and the effects of the friction coefficient on the stress intensity factors of the pipe with an inclined inner surface crack under axial compression load and external pressure were explored in this paper. To be fit for defect assessment, the solutions for stress intensity factors K_II and K_III were derived, and new correction factors f_θ and f_μ were proposed in the empirical solutions to accommodate the crack inclination asymmetry and the friction coefficient, respectively.

1. Introduction

Pipes are widely used in various fields of petroleum, chemical industry, and natural gas. In engineering practice, there are a large number of shells subjected to external pressure loads such as submarines, aerospace simulators, deep-sea pipelines, and buried pipelines. Due to the limitations of the manufacturing process, the structures may have tiny defects such as cracks. Lin and Smith [1] found that arbitrarily shaped cracks become semi-elliptical after a certain period of expansion. The stress intensity factor is an important fracture parameter for judging crack propagation and failure in defect assessment. Therefore, since the last century, scholars have conducted many studies on this subject.

In the early years, Raju and Newman [2] established a stress intensity factor empirical solution for mode I internal and external surface cracks of cylindrical vessels through FEA. The formula gave the influence coefficients based on different shell and crack sizes. Navid and Fenner [3] calculated the stress intensity factors of a thick-walled cylinder with internal pressure for different crack sizes based on the boundary integral equation method. Afterward, various numerical techniques were developed to calculate the stress intensity factor of pipe surface cracks such as the line spring model [4] and the weight function method [5,6,7]. Kamaya [8] combined the finite element alternating method and the finite element analysis to calculate the stress intensity factors of shell surface cracks. Wallbrink [9] proposed a semi-analytical method of conformal transformation to predict the stress intensity factors of circumferential surface cracks. However, most of these studies were limited to simple loading and the crack direction was axial or circumferential. In fact, in most cases, the loading conditions of the pipelines are complex and the crack direction is inclined.

The research on mixed mode cracks first started from a flat plate. Murakami [10] analyzed K_Ⅰ, K_Ⅱ, K_Ⅲ by the body force method, and provided a formula for the maximum value K_θmax at the crack front according to the projected area of the crack in the direction of the maximum principal stress. Zeng and Dai [11] proposed a simplified analytical model of an inclined surface crack under biaxial stress, and proposed a closed solution for K_Ⅰ and K_Ⅲ at the deepest point of the crack front. Some scholars have also used FEA [12,13], the numerical influence function method [14], and the experimental method [15] to study the mixed mode cracks of flat plates. The study of mixed mode cracks in shell can be divided into two types. First, the crack direction is determined, usually axial or circumferential, and the mixed mode crack is formed by applying a complex far-field load on the structure. Shahani and Habibi [16] considered the stress intensity factor of cylinders with a circumferential crack under the action of axial force, bending moment, and torque. Predan et al. [17] studied the stress intensity factor for circumferential semi-elliptical surface cracks in a hollow cylinder subjected to pure torsion. Ramezani et al. [18] demonstrated the empirical solution of the stress intensity factor of cylinder surface crack under pure torsion. Second, the mixed mode cracks were formed by changing the inclination angle of the surface cracks. For example, Li et al. [19,20] studied the stress intensity factor of the inclined crack in the pipeline under the far field tension and tension-bending, and provided an empirical formula to calculate the complete value of the stress intensity factor at the crack front through the influence coefficient under different aspect ratios. Li and Mao [21] calculated the stress intensity factor of the inclined surface crack of the outer wall of the heat exchanger by numerical simulation.

In the research on cracks, it is inevitable to encounter the problem of complete or partial closure of the cracks. For simplicity, more papers have chosen to ignore the surface contact of the cracks. However, the contact closure of the crack face significantly changes the stress–strain distribution at the crack front, thereby affecting the fracture behaviors of the structure. In fact, since the contact area is time-varying, the contact problem is a typical nonlinear problem with boundary conditions. In static fracture mechanics, the contact friction of the crack face has an obvious influence on the stress intensity factor. At present, there have been few studies on mixed mode cracks on the friction surfaces under compressive loading, which have mainly focused on flat plates [22,23,24]. Bowie and Freese [25] proposed a crack closure technique to correct the solution of overlapping crack surfaces, but their method does not consider the sliding between crack surface, and is only suitable for the case of a large friction coefficient. Liu and Tan [26] used the boundary element method to study the effect of the interaction between the friction and crack surface on the stress intensity factor. Hammouda et al. [27,28] used finite element analysis to study the effect of crack surface friction and crack inclination on the stress intensity factor of the central cracked plate under unidirectional compressive load. In addition, Dorogoy and Banks-Sills [29] investigated the effects of loading angle and friction coefficient on the stress intensity factor and the crack length of Brazilian disc cracks under concentrated loading using a finite difference solution.

To the best of the authors’ knowledge, researchers have not focused on the surface crack stress intensity factor of shell structures subjected to compressive loads, however, in engineering practice, fracture damage occurs in shell structures subjected to simultaneous external and axial pressures. In this paper, the stress intensity factors of the shell subjected to simultaneous external and axial pressures were evaluated by using the three-dimensional finite element analysis method, and the effects of the relative depth of the surface cracks, aspect ratio, crack inclination, and friction coefficient on the stress intensity factor of the crack fronts were investigated, which were different from those under internal pressure. The empirical solutions of the mode II and mode III stress intensity factors along the crack front are given using numerical methods, respectively.

2. Finite Element Model of Pipe Crack

2.1. Meshing of the Pipe with Inclined Crack on the Inner Surface

For the study of mixed mode cracks, the most widely used and accurate method is FEA, which involves the extended finite element method (XFEM) and contour integral. Based on the static loading model in this paper, the fluctuation of results caused by inaccurate meshing can be better avoided by using contour integrals. The commercial code ABAQUS [30] was used for FEA. The material used in this paper was TA2 with a modulus of elasticity of 101,901 MPa and a Poisson’s ratio of 0.348. The pipe size was fixed at R_i = 100 mm, R₀ = 110 mm, t = 10 mm. In order to avoid the influence of the boundary effect on the stress intensity factor, the pipe length should be six times that of R_i and 20 times that of the crack half-length c, l = 600 mm [31,32]. The crack size is described by dimensionless parameters: the relative crack depth (a/t) and the ratio of long and minor axes of semi-elliptical cracks (a/c). The crack inclination angle is θ. The pipe model is shown in Figure 1a. The external pressure of the pipe P₀ = 100 MPa and the axial pressure p = 525 MPa. The boundary conditions were that one end of the pipe was completely fixed and the other end restricted its circumferential constraint, as shown in Figure 1b.

The current general semi-elliptical crack meshing method was used by the authors [31,33]. The six-node triangular prism element C3D6 was used to sweep along the path around the crack tip to simulate the stress–strain singularity at the crack front. C3D6 belongs to a fully integrated element with linear interpolation in all directions, the number of nodes is 6, and the number of integration points is 6. The eight-node hexahedral element C3D8R was used in the nearby area. C3D8R belongs to a linear reduced integral element with eight nodes and only one integral point exists in the center of the cell, which is equivalent to a constant stress cell. The reduced integration element can effectively avoid stress discontinuity problems. The meshes at the crack were subdivided, the mesh size was controlled within 1 mm, and the mesh size of the rest zone was controlled within 10 mm to reduce the calculation time. The mesh division is shown in Figure 2.

2.2. Penalty Function Method for Contact Problem

Since the pipe is not only subjected to external pressure but also to axial compressive load, contact should be set between the crack surfaces. Generally speaking, the crack surfaces in the contact state have three characteristics:

• (1). The contact surfaces do not penetrate or overlap each other;

• (2). The contact surfaces are able to transmit normal pressure and tangential friction;

• (3). The contact surfaces generally do not transmit normal tension and can be separated freely.

In ABAQUS, the penalty function method, Lagrange multiplier method, and static friction–kinetic friction index decay method can be used to solve the contact problem. In this paper, a penalty function was used to characterize the tangential friction force between the crack surfaces. The penalty function method requires the normal and tangential friction coefficients, which is similar to setting a “spring” between the contact surfaces, and the “spring” works only when the contact surfaces are closed. In this way, the normal contact pressure can be expressed as:

(1)$P_n=\{\begin{array}{l}0 u_n>0\\ k_n{u}_n u_n\le 0\end{array}$

According to Coulomb’s law of friction, the frictional stress at the contact surface can be written as:

(2)$\begin{array}{l}{\tau }_s=k_s{u}_s-\mu P_n\\ {\tau }_t=k_t{u}_t-\mu P_n\end{array}$

where t, n, s represent the tangential, normal and sub-normal directions, respectively, as shown in Figure 3. The advantage of the penalty function is that it does not increase the degree of freedom of the problem. The disadvantage is that when the friction coefficient is too large, the convergence will become more difficult. The values of the friction coefficient μ in this paper were 0, 0.2, and 0.4, respectively.

2.3. Mesh Independence Verification

Raju and Newman [2] first proposed a formula to calculate the stress intensity factor of a mode I crack in a pipe, and Chun-Qing Li [19] extended this formula to calculate the mixed-mode stress intensity factor of an inclined crack, as follows:

(3)$\mathit{K}={\sigma }_0\sqrt{\pi \frac{a}{Q}}\mathit{F}\left(\frac{a}{t},\frac{a}{c},\frac{t}{R_i},\xi ,\theta \right)$

where σ₀ = PR/t is the average hoop stress of the pipe; K and F are the stress intensity factor and influence coefficient function under different cracking modes, respectively; Q is the shape factor obtained from the second type of elliptic integral; and the empirical formula is given by Shiratori et al. [34].

(4)$\mathit{K}=\left\{\begin{array}{l}{K}_I\hfill \\ K_{II}\hfill \\ K_{III}\hfill \end{array}\right\}$

(5)$\mathit{F}\left(\frac{a}{t},\frac{a}{c},\frac{t}{R_i},\xi ,\theta \right)=\left\{\begin{array}{l}{F}_I\hfill \\ F_{II}\hfill \\ F_{III}\hfill \end{array}\right\}$

(6)$\begin{array}{ll}Q=1+1.464{\left(\frac{a}{c}\right)}^{1.65}\hfill & \frac{a}{c}⩽1\hfill \\ Q=\left[1+1.464{\left(\frac{c}{a}\right)}^{1.65}\right]{\left(\frac{a}{c}\right)}^2\hfill & \frac{a}{c}>1\hfill \end{array}$

In order to unify the evaluation indicators, it is necessary to normalize the stress intensity factor and the position along the crack front. In this paper, the normalization factor K₀ was defined as [18]:

(7)$K_0=\sigma \sqrt{\pi \cdot a_0} a_0=1$

where σ = PR/2t is the axial compressive stress. The position along the crack front can be normalized according to the number of equal meshes. The advantage of this normalization method is that K₀ becomes a constant, which can truly reflect the variation trend of K.

Figure 4 shows the normalized stress intensity factors at the crack front of three kinds of mesh sizes at different angles. The mesh size is shown in Table 1. It can be seen that the results of the first two mesh sizes were basically the same. A comparison with the normalized K_II and K_III results with Chun-Qing Li [19] after changing the load using mesh2 is given in Figure 5. The average differences were 8.98% and 2.55%, respectively, which demonstrates the accuracy of the results.

2.4. Comparison of SIFs between Crack Opening and Crack Closing

In order to compare the change in the stress intensity factor for crack opening and crack closing, two sets of cracks with different angles were selected in this paper. The contact surface friction coefficient μ = 0. The pipe load for crack opening was changed to the internal pressure and axial tension, and the results of the comparison are shown in Figure 6. To facilitate comparison of the differences, the normalized stress intensity factors under external pressure and axial compression were expressed symmetrically about the transverse coordinate. It can be seen that the absolute values of the normalized stress intensity factors at crack closure were obviously higher for both the mode II crack and mode III crack. For this phenomenon, the authors used the stress field calculation method for the plastic zone of the crack tip mentioned by Shlyannikov and Tumanov [35] to convert the three stress states obtained from the finite element calculations to give the sub-normal stresses and tangential stresses at the crack front in these two cases, corresponding to mode II and mode III cracks, respectively, as shown in Figure 7a,b, and it can be seen that τ_s and τ_t at the crack front of the pipe under axial compression and external pressure were higher than those under axial tension and internal pressure, which is the fundamental reason for the difference in SIFs. This difference also proves that it is necessary to study the stress intensity factor of cracks under axial compression load and external pressure.

3. Results and Discussion

The stress intensity factor for mode I crack is 0 under the compression load [36,37]. However, due to the crack inclination and external pressure, the K_II and K_III generated by the interaction of tangential and frictional forces at the crack front cannot be neglected. In this paper, the SIFs of pipe surface crack fronts were obtained by FEA under different influencing factors, which included the relative depth (a/t), the aspect ratio (a/c), the crack inclination angle (θ), and the friction coefficient (μ). The specific calculation scheme is shown in Table 2. A total of 336 finite element models were calculated in this paper.

3.1. Effects of Crack Geometry and Inclination Angle

Figure 8 shows the variation in the mode II and mode III SIFs with relative depth. The long axes of the semi-elliptical cracks were fixed as c = 5 mm and c = 10 mm, respectively, the friction coefficient μ was 0. For the convenience of presentation, only the comparison results of the maximum value along the crack front are given, and the rules were consistent for the remaining positions. It can be seen that the absolute value of the mode II and mode III stress intensity factors increased with the deepening of the crack, and the increase in the absolute values was more obvious when the crack size was larger. The surface crack was most dangerous at θ = 45° and became safer with the offset of crack inclination to the axial and circumferential directions.

The variation in the mode II and mode III SIFs with the aspect ratio is given in Figure 9, and the minor axis of the semi-elliptical crack was fixed as a = 4 mm and a = 6 mm, respectively, the friction coefficient μ was 0. For the mode II cracks, when a/c < 1, the absolute value of the normalized K_II increased with the aspect ratio, and when a/c > 1, the absolute value of the normalized K_II decreased with the increase in the aspect ratio. The K_II values of the crack fronts for different aspect ratios at two depths are given in Figure 9e,f. For semi-elliptical cracks with lower aspect ratios, the maximum value of K_II always appeared at the surface point of the crack front. However, for semi-elliptical cracks with higher aspect ratios, the maximum value of K_II appeared near the surface point, which led to a smaller value at the surface point. These variations are consistent with the results of Yang et al. [38] for high aspect ratio cracks. For mode III cracks, the absolute value of K_III decreased monotonically with the increase in aspect ratio, with the most pronounced decrease at θ = 45°.

The variation in the mode II and mode III SIFs with the crack inclination angle is given in Figure 10, the friction coefficient μ was 0. It can be seen that in all cases, the absolute values of K_II and K_III were symmetrical about the deepest point. The maximum value of the stress intensity factor always appeared at the surface point or the deepest point at θ = 45°, and the SIF value was smaller when the crack inclination was closer to the axial and circumferential directions, which was due to the more obvious closure effect under compression load, and the tangential component parallel to the crack surface became smaller. It should be noted that for small-size cracks, K_II and K_III were almost identical when θ was symmetric about 45° (θ = 30°, 60° and θ = 15°, 75°), and the pipe surface cracks at this time could be approximated as flat plate cracks, which has been verified in the study of mixed mode cracks in flat plates and pipes [13,20]. However, for large size pipe surface cracks, the difference caused by the angle was more obvious. The crack surface had a different curvature on the inner wall surface of the pipe at different inclination angles, which made the stress intensity factors no longer consistent for θ = 30°, 60° and θ = 15°, 75°. As shown in Figure 10e,f, this difference due to curvature was more pronounced in spherical shells [39].

3.2. Effect of Friction Coefficient on Contact Surface

The variation in the mode II and mode III SIFs with the friction coefficient is given in Figure 11. Due to the symmetry of the curves, only the results for the half-length along the crack front are shown in the figure. The absolute values of K_II and K_III decreased with an increasing friction coefficient for different crack sizes and crack inclination angles. At θ = 45°, the friction coefficient had the greatest effect on the SIFs, and this effect also increased as the crack size increased. For the mode II cracks, the friction coefficient had the greatest effect on the SIF at the surface point, and this effect decreased as the location of the crack front moved forward, reaching a minimum at the deepest point. In contrast to mode II cracks, the friction coefficient significantly changed the SIF at the deepest point of mode III, but had a dropping effect on the surface point.

4. Empirical Solution of SIFs for the Pipe with Inclined Inner Surface Cracks under External Pressure and Axial Compression

Although FEA is one of the most effective methods for calculating SIFs, it is difficult to apply in engineering practice due to the complexity of modeling and time-consuming calculations. In this paper, a new form of empirical solution was proposed by least-squares fitting based on the FEA results to give the influence coefficients at different a/t with the normalized crack front position ξ as the basis function. A new correction factor for the inclination angle, f_θ, was proposed based on the effect caused by the curvature of the pipe above-mentioned in the effect of the crack inclination angle. The friction coefficient influence coefficient, f_μ, was proposed because the pipe is subjected to external pressure and axial pressure. The influence coefficients are shown in Table 3 and Table 4, and the fitted equations are shown in Equations (8)–(16).

(8)$\mathit{K}=\sigma \sqrt{\pi \frac{a}{Q}}\mathit{F}\left(\frac{a}{t},\frac{a}{c},\xi ,f_{\theta },f_{\mu }\right)$

where σ is the far-field compressive stress, $\mathit{K}=\left\{\begin{array}{l}{K}_{II}\hfill \\ K_{III}\hfill \end{array}\right\}$,$\mathit{F}\left(\frac{a}{t},\frac{a}{c},\xi ,f_{\theta },f_{\mu }\right)=\left\{\begin{array}{l}{F}_{II}\hfill \\ F_{III}\hfill \end{array}\right\}$.

Influence coefficient functions for K_II:

(9)$F_{II}=(H_1+H_2\xi +H_3{\xi }^2+H_4{\xi }^3)f_{\theta }{f}_{\mu }$

(10)$H_{i(i=1,2,3,4)}=h_1+h_2(\frac{a}{c})+h_3{(\frac{a}{c})}^2+h_4{(\frac{a}{c})}^3$

(11)$f_{\theta }=\left\{1+\frac{ac}{600}\left[\sin \left(\theta -45\right)\right]\right\}\sin 2\theta$

(12)$f_{\mu }=\left(1+0.575\mu \right){\left(ac\right)}^{-0.505\mu }$

Influence coefficient functions for K_III:

(13)$F_{III}=(H_1+H_2\xi +H_3{\xi }^2)f_{\theta }{f}_{\mu }$

(14)$H_{i(i=1,2,3)}=h_1+h_2(\frac{a}{c})+h_3{(\frac{a}{c})}^2+h_4{(\frac{a}{c})}^3$

(15)$f_{\theta }=\left\{1+\frac{ac}{1500}\left[\sin \left(45-\theta \right)\right]\right\}\sin 2\theta$

(16)$f_{\mu }=\left(1+0.281\mu \right){\left(ac\right)}^{-0.511\mu }$

The fitting results for θ = 45° and μ = 0 at the three sizes are given in Figure 12a,b. The results of the fit at different angles for the same crack size are given in Figure 12c,d, where μ = 0. The fitting results for the same crack size with different angles and different friction coefficients are given in Figure 12e,f; due to the symmetry of the curves, only the results for the half-length of the crack front are shown in the figure. It can be seen that the fitting results achieved good agreement with the FEA results, which can be applied for the defect assessment.

5. Conclusions

In this paper, the SIFs of the inclined crack on the inner surface of the pipe under combined axial pressure load and external pressure were investigated by the finite element method, and the main conclusions are as follows.

(1) By comparing the SIF values along the crack fronts, it can be found that K_II and K_III for the closed surface cracks under external and axial pressure were higher than K_II and K_III for the open surface cracks under internal pressure and axial tension.

(2) The effects of the relative depth and aspect ratio on the SIFs of the inclined inner crack were analyzed. The results showed that the larger the a and c, the easier it is for the crack to expand, and when a/c >1, K_IImax does not appear at the surface point of the crack, but near the surface point.

(3) The effects of the crack inclination angle on the SIFs of the inclined inner crack were analyzed. The results showed that the SIFs at the crack front were the largest at θ = 45°, and their values decreased with an inclination angle toward the axial and circumferential directions. The larger the crack size, the more obvious the asymmetry of the SIFs about θ = 45°.

(4) The effects of the friction coefficient on the SIFs of surface crack were analyzed. The larger the friction factor, the smaller the SIFs along the crack front. The friction coefficient had the greatest effect on the surface point of mode II cracks and the deepest point of mode III cracks, and the greatest effect on cracks with θ = 45°.

(5) Based on the above results, a new solution for stress intensity factors K_II and K_III were proposed, and the corresponding coefficients for different crack sizes, an angle correction factor f_θ, and a friction correction factor f_μ are given.

Author Contributions

X.-M.Y.: Conception, methodology, software, writing—manuscript, and designed the research. Y.-C.Z., Q.P. and L.-Z.J.: Software. T.-H.M.: Investigation. X.-H.H.: Funding acquisition. C.-Y.Z.: Formal analysis, conception and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors gratefully acknowledge the financial support of the National Natural Science Foundation of China (51975271).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

a	Depth of a semi-elliptical surface crack
c	Half-length of a semi-elliptical surface crack
FI	Influence coefficient functions for KI
FII	Influence coefficient functions for KII
FIII	Influence coefficient functions for KIII
fθ	Angle effect correction factor
fμ	Frictional influence correction factor
H1, H2, H3, H4	Sub-curve-fitting functions for FI, FII, FIII
h1, h2, h3, h4	Sub-curve fitting constants for H1, H2, H3, H4
KI	Stress intensity factors for Modes I crack
KII	Stress intensity factors for Modes II crack
KIII	Stress intensity factors for Modes III crack
K0	Normalizing factor for the stress intensity factors
kn	Main normal stiffness of crack front
ks	Sub-normal stiffness of crack front
kt	Tangential stiffness of crack front
l	Length of a pipe
P	Axial pressure of a pipe
P0	External pressure of a pipe
Pn	Main normal pressure of crack front
Q	Shape factor for elliptical crack
Ri,R0	Internal and external radius of a pipe
t	Wall thickness of a pipe
un	Main normal displacement of crack front
us	Sub-normal displacement of crack front
ut	Tangential displacement of crack front
θ	Inclination angle of a surface crack
ξ	Normalized position along the crack front
σ	Far-field compressive stress
τs	Crack front sub-normal shear stress
τt	Crack front tangential shear stress
μ	Friction coefficient of the cracked surface

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. (a) Pressure pipe with an inclined crack on the inner surface. (b) Loads and boundary conditions.

Figure 2. Meshing of the pipe with the inclined internal surface crack:(a) global view of the cracked pipe; (b) equal divisions at the crack front; (c) hexahedron elements around the crack front; (d) wedge-shaped elements at the crack front.

Figure 3. Local coordinate system of the quarter-elliptic crack front.

Figure 4. Effect of the mesh density on the normalized SIF at different angles, a/t = 0.2, a/c = 0.2. (a) Normalized KII, (b) normalized KIII.

Figure 5. Comparison of the stress intensity factors with mixed modes, a/t = 0.2, a/c = 0.4, θ = 45°. (a) Normalized KII, (b) normalized KIII [19].

Figure 6. Comparison of SIFs under tension and compression loads with different angles, μ = 0, a/t = 0.2, a/c = 0.8. (a) Normalized KII, (b) normalized KIII.

Figure 7. Comparison of shear stress along the crack fronts under tension and compression loads. (a) Shear stress along the crack front of mode II. (b) Shear stress along the crack front of mode III.

Figure 8. Effect of the relative depth on the normalized SIFs, μ = 0. (a) Normalized KII (c = 5 mm, a/c = 0.4, 0.8, 1.6); (b) normalized KIII (c = 5 mm, a/c = 0.4, 0.8, 1.6); (c) normalized KII (c = 10 mm, a/c = 0.2, 0.4, 0.8); (d) normalized KIII (c = 10 mm, a/c = 0.2, 0.4, 0.8).

Figure 9. Effect of the aspect ratio on the normalized SIFs (μ = 0). (a) Normalized KII (a = 4 mm, a/t = 0.4); (b) normalized KIII (a = 4 mm, a/t = 0.4); (c) normalized KII (a = 6 mm, a/t = 0.6); (d) normalized KIII (a = 6 mm, a/t = 0.6); (e) normalized KII near the crack surface (a = 4 mm, a/t = 0.4); (f) normalized KII near the crack surface (a = 6 mm, a/t = 0.6).

Figure 10. Effect of the inclination angle on the normalized SIFs, μ = 0. (a) Normalized KII (a/t = 0.2, a/c = 1.6); (b) normalized KIII (a/t = 0.2, a/c = 1.6); (c) normalized KII (a/t = 0.4, a/c = 0.8); (d) normalized KIII (a/t = 0.4, a/c = 0.8); (e) normalized KII (a/t = 0.6, a/c = 0.4); (f) normalized KIII (a/t = 0.6, a/c = 0.4).

Figure 11. Effect of the coefficient of friction on normalized SIFs. (a) Normalized KII (a/t = 0.2, a/c = 1.6); (b) normalized KIII (a/t = 0.2, a/c = 1.6); (c) normalized KII (a/t = 0.6, a/c = 0.4); (d) normalized KIII (a/t = 0.6, a/c = 0.4).

Figure 12. Comparison of SIFs fitting results with (a,b) different crack sizes (a/t = 0.2, 0.4, 0.6, a/c = 1.6, 0.8, 0.4, θ = 45°, μ = 0); (c,d) different inclination angles (a/t = 0.6, a/c = 0.4, μ = 0); (e,f) different friction coefficients (a/t = 0.4, a/c = 0.8, μ = 0, 0.2, 0.4).

Table 1

Crack front meshing.

Division Strategy	Mesh1	Mesh2	Mesh3
Wedge element size (mm)	0.05	0.05	0.1
Hexahedral element layers	8	6	3
Equal divisions at the crack front	180	90	45

Table 2

Finite element calculation scheme.

Influencing Factors
a/t	0.2	0.4	0.6	0.8
a/c	0.2	0.4	0.8	1.6
θ(°)	0	15	30	45	60	75	90
μ	0	0.2	0.4

Table 3

The sub-curve fitting constants in Equation (10) h_i.

	a/t = 0.2				a/t = 0.4
Constants	H 1	H 2	H 3	H 4	H 1	H 2	H 3	H 4
h 1	0.161	−0.323	0.159	−0.117	−0.121	1.404	−3.037	1.865
h 2	−2.105	10.376	−18.905	12.616	−0.753	3.320	−7.455	5.953
h 3	1.795	−14.836	34.186	−22.826	0.285	−6.828	21.492	−15.793
h 4	−0.492	4.623	−11.073	7.399	0.020	1.850	−6.686	5.024
	a/t = 0.6				a/t = 0.8
Constants	H 1	H 2	H 3	H 4	H 1	H 2	H 3	H 4
h 1	0.338	0.317	−2.580	1.678	0.385	1.134	−5.297	3.484
h 2	−2.987	9.044	−10.429	7.042	−3.309	7.181	−2.875	2.012
h 3	3.129	−14.121	24.981	−16.777	3.648	−13.768	20.741	−13.951
h 4	−0.999	4.641	−8.408	5.654	−1.206	4.955	−8.080	5.433

Table 4

The sub-curve fitting constants in Equation (14) h_i.

	a/t = 0.2			a/t = 0.4
Constants	H 1	H 2	H 3	H 1	H 2	H 3
h 1	−0.448	−0.582	0.587	−0.497	−1.179	1.171
h 2	0.599	−1.998	2.023	0.727	1.049	−0.907
h 3	−0.322	1.544	−1.596	−0.491	−2.190	1.959
h 4	0.053	−0.313	0.335	0.116	1.001	−0.910
	a/t = 0.6			a/t = 0.8
Constants	H 1	H 2	H 3	H 1	H 2	H 3
h 1	−0.451	−0.187	0.210	−0.375	−1.135	1.159
h 2	0.494	−3.521	3.478	0.154	−0.353	0.306
h 3	−0.214	3.819	−3.781	0.200	0.413	−0.373
h 4	0.021	−1.192	1.181	−0.126	−0.058	0.045

References

1. Lin, X.B.; Smith, R.A. Shape Growth Simulation of Surface Cracks in Tension Fatigued Round Bars. Int. J. Fatigue; 1997; 19, pp. 461-469. [DOI: https://dx.doi.org/10.1016/S0142-1123(97)00012-1]

2. Raju, I.S.; Newman, J.C., Jr. Stress-Intensity Factors for Internal and External Surface Cracks in Cylindrical Vessels. J. Press. Vessel Technol.; 1982; 104, pp. 293-298. [DOI: https://dx.doi.org/10.1115/1.3264220]

3. Navid, H.; Fenner, R.T.; Nadiri, F.; Webster, G.A. Stress Intensity Factors for Internal and External Cracks in Pressurised Thick-Walled Cylinders. Int. J. Press. Vessel Pip.; 1985; 18, pp. 241-254. [DOI: https://dx.doi.org/10.1016/0308-0161(85)90013-4]

4. Kumar, V.; German, M.D.; Schumacher, B.I. Analysis of Elastic Surface Cracks in Cylinders Using the Line-Spring Model and Shell Finite Element Method. J. Press. Vessel Technol.; 1985; 107, pp. 403-411. [DOI: https://dx.doi.org/10.1115/1.3264474]

5. Zheng, X.J.; Kiciak, A.; Glinka, G. Weight Functions and Stress Intensity Factors for Internal Surface Semi-Elliptical Crack in Thick-Walled Cylinder. Eng. Fract. Mech.; 1997; 58, pp. 207-221. [DOI: https://dx.doi.org/10.1016/S0013-7944(97)00083-0]

6. Jones, I.S.; Rothwell, G. Reference Stress Intensity Factors with Application to Weight Functions for Internal Circumferential Cracks in Cylinders. Eng. Fract. Mech.; 2001; 68, pp. 435-454. [DOI: https://dx.doi.org/10.1016/S0013-7944(00)00111-9]

7. Eshraghi, I.; Soltani, N. Stress Intensity Factor Calculation for Internal Circumferential Cracks in Functionally Graded Cylinders Using the Weight Function Approach. Eng. Fract. Mech.; 2015; 134, pp. 1-19. [DOI: https://dx.doi.org/10.1016/j.engfracmech.2014.12.007]

8. Kamaya, M.; Nishioka, T. Analysis of Surface Crack in Cylinder by Finite Element Alternating Method. J. Press. Vessel Technol.; 2005; 127, pp. 165-172. [DOI: https://dx.doi.org/10.1115/1.1903003]

9. Wallbrink, C.D.; Peng, D.; Jones, R. Assessment of Partly Circumferential Cracks in Pipes. Int. J. Fract.; 2005; 133, pp. 167-181. [DOI: https://dx.doi.org/10.1007/s10704-005-0628-0]

10. Murakami, Y. Analysis of Stress Intensity Factors of Modes I, II and III for Inclined Surface Cracks of Arbitrary Shape. Eng. Fract. Mech.; 1985; 22, pp. 101-114. [DOI: https://dx.doi.org/10.1016/0013-7944(85)90163-8]

11. Zhao-Jing, Z.; Shu-Ho, D. Stress Intensity Factors for an Inclined Surface Crack under Biaxial Stress State. Eng. Fract. Mech.; 1994; 47, pp. 281-289. [DOI: https://dx.doi.org/10.1016/0013-7944(94)90228-3]

12. Ayhan, A.O. Mixed Mode Stress Intensity Factors for Deflected and Inclined Corner Cracks in Finite-Thickness Plates. Int. J. Fatigue; 2007; 29, pp. 305-317. [DOI: https://dx.doi.org/10.1016/j.ijfatigue.2006.03.006]

13. Şahin, H.; Ayhan, A.O. Three-Dimensional Mixed Mode Stress Intensity Factors for Inclined Elliptical Surface Cracks in Plates under Uniform Tensile Load. Procedia Struct. Integr.; 2019; 21, pp. 38-45. [DOI: https://dx.doi.org/10.1016/j.prostr.2019.12.084]

14. Phyo Myat, K.; Osawa, N.; Gadallah, R.; Tanaka, S. Accurate and Efficient Method for Analyzing Mixed-Mode SIFs for Inclined Surface Cracks in Semi-Infinite Bodies by Using Numerical Influence Function Method. Theor. Appl. Fract. Mech.; 2020; 106, 102471. [DOI: https://dx.doi.org/10.1016/j.tafmec.2019.102471]

15. Bhat, M.A.; Shaikh, A.A. Effect of Specimen Parameters on Mixed-Mode I/II Stress Intensity Factors for Additive Manufactured Slant Edge Crack Plate. Mater. Today Proc.; 2021; 44, pp. 4305-4308. [DOI: https://dx.doi.org/10.1016/j.matpr.2020.10.549]

16. Shahani, A.R.; Habibi, S.E. Stress Intensity Factors in a Hollow Cylinder Containing a Circumferential Semi-Elliptical Crack Subjected to Combined Loading. Int. J. Fatigue; 2007; 29, pp. 128-140. [DOI: https://dx.doi.org/10.1016/j.ijfatigue.2006.01.017]

17. Predan, J.; Močilnik, V.; Gubeljak, N. Stress Intensity Factors for Circumferential Semi-Elliptical Surface Cracks in a Hollow Cylinder Subjected to Pure Torsion. Eng. Fract. Mech.; 2013; 105, pp. 152-168. [DOI: https://dx.doi.org/10.1016/j.engfracmech.2013.03.033]

18. Ramezani, M.K.; Purbolaksono, J.; Andriyana, A.; Ramesh, S.; Putra, I.S. Empirical Solutions for Stress Intensity Factors of a Surface Crack in a Solid Cylinder under Pure Torsion. Eng. Fract. Mech.; 2018; 193, pp. 122-136. [DOI: https://dx.doi.org/10.1016/j.engfracmech.2018.02.015]

19. Li, C.-Q.; Fu, G.; Yang, W. Stress Intensity Factors for Inclined External Surface Cracks in Pressurised Pipes. Eng. Fract. Mech.; 2016; 165, pp. 72-86. [DOI: https://dx.doi.org/10.1016/j.engfracmech.2016.08.014]

20. Fu, G.; Yang, W.; Li, C.-Q. Stress Intensity Factors for Mixed Mode Fracture Induced by Inclined Cracks in Pipes under Axial Tension and Bending. Theor. Appl. Fract. Mech.; 2017; 89, pp. 100-109. [DOI: https://dx.doi.org/10.1016/j.tafmec.2017.02.001]

21. Li, D.; Mao, Z. Experimental and Numerical Simulations on Compound Stress Intensity Factor of Semi-Elliptical Cracks on the Exchanger Outer Walls with Inclined Angles. Alex. Eng. J.; 2022; 61, pp. 5065-5072. [DOI: https://dx.doi.org/10.1016/j.aej.2021.09.030]

22. Woo, C.W.; Cheung, Y.K.; Chen, Y.Z.; Wang, Y.H. A Simple Model for the Contact Problem of a Finite Cracked Plate in Bending. Eng. Fract. Mech.; 1988; 29, pp. 227-231. [DOI: https://dx.doi.org/10.1016/0013-7944(88)90049-5]

23. Beghini, M.; Bertini, L. Effective Stress Intensity Factor and Contact Stress for a Partially Closed Griffith Crack in Bending. Eng. Fract. Mech.; 1996; 54, pp. 667-678. [DOI: https://dx.doi.org/10.1016/0013-7944(95)00234-0]

24. Thiagarajan, S.; Alwar, R.S. Influence of Crack Closure in the Case of an Angled Crack. Eng. Fract. Mech.; 1986; 24, pp. 533-537. [DOI: https://dx.doi.org/10.1016/0013-7944(86)90227-4]

25. Bowie, O.L.; Freese, C.E. On the “Overlapping” Problem in Crack Analysis. Eng. Fract. Mech.; 1976; 8, pp. 373-379. [DOI: https://dx.doi.org/10.1016/0013-7944(76)90017-5]

26. Liu, S.B.; Tan, C.L. Two-Dimensional Boundary Element Contact Mechanics Analysis of Angled Crack Problems. Eng. Fract. Mech.; 1992; 42, pp. 273-288. [DOI: https://dx.doi.org/10.1016/0013-7944(92)90218-4]

27. Hammouda, M.M.I.; Fayed, A.S.; Sallam, H.E.M. Mode II Stress Intensity Factors for Central Slant Cracks with Frictional Surfaces in Uniaxially Compressed Plates. Int. J. Fatigue; 2002; 24, pp. 1213-1222. [DOI: https://dx.doi.org/10.1016/S0142-1123(02)00048-8]

28. Hammouda, M.M.I.; Fayed, A.S.; Sallam, H.E.M. Stress Intensity Factors of a Shortly Kinked Slant Central Crack with Frictional Surfaces in Uniaxially Loaded Plates. Int. J. Fatigue; 2003; 25, pp. 283-298. [DOI: https://dx.doi.org/10.1016/S0142-1123(02)00146-9]

29. Dorogoy, A.; Banks-Sills, L. Effect of Crack Face Contact and Friction on Brazilian Disk Specimens—A Finite Difference Solution. Eng. Fract. Mech.; 2005; 72, pp. 2758-2773. [DOI: https://dx.doi.org/10.1016/j.engfracmech.2005.05.005]

30. DASSAULT SYSTEMES. 2021; Available online: https://www.3ds.com/support/hardware-and-software/simulia-system-information/abaqus-2021/ (accessed on 19 December 2022).

31. Zhang, T.; Lu, K.; Katsuyama, J.; Li, Y. Stress Intensity Factor Solutions for Surface Cracks with Large Aspect Ratios in Cylinders and Plates. Int. J. Press. Vessel Pip.; 2021; 189, 104262. [DOI: https://dx.doi.org/10.1016/j.ijpvp.2020.104262]

32. Randeniya, C.; Robert, D.J.; Fu, G.; Li, C.Q. The Effect of Corrosion Patch Geometry on Stress Intensity Factors for External Surface Cracks in Cast Iron Water Mains. Proceedings of the Fourth International Conference on Sustainable Construction Materials and Technologies; Las Vegas, NV, USA, 7–11 August 2016.

33. Teh, S.; Andriyana, A.; Ramesh, S.; Putra, I.S.; Kadarno, P.; Purbolaksono, J. Tetrahedral Meshing for a Slanted Semi-Elliptical Surface Crack at a Solid Cylinder. Eng. Fract. Mech.; 2021; 241, 107400. [DOI: https://dx.doi.org/10.1016/j.engfracmech.2020.107400]

34. Shiratori, M.; Miyoshi, T.; Sakai, Y.; Zhang, G.R. Analysis of Stress Intensity Factors for Surface Cracks Subjected to Arbitrarily Distributed Surface Stresses. Trans. Jpn. Soc. Mech. Eng.; 1987; 53, pp. 1651-1657. [DOI: https://dx.doi.org/10.1299/kikaia.53.1651]

35. Shlyannikov, V.N.; Tumanov, A.V. An Inclined Surface Crack Subject to Biaxial Loading. Int. J. Solids Struct.; 2011; 48, pp. 1778-1790. [DOI: https://dx.doi.org/10.1016/j.ijsolstr.2011.02.024]

36. Cordes, R.D.; Joseph, P.F. Crack Surface Contact of Surface and Internal Cracks in a Plate with Residual Stresses. Int. J. Fract.; 1994; 66, pp. 1-17. [DOI: https://dx.doi.org/10.1007/BF00012632]

37. Zhu, F.; Liu, H.; Yao, L.; Mei, G. Stress Field Analysis of an Infinite Plate with a Central Closed Inclined Crack under Uniaxial Compression. Theor. Appl. Fract. Mech.; 2021; 116, 103111. [DOI: https://dx.doi.org/10.1016/j.tafmec.2021.103111]

38. Yang, S.T.; Ni, Y.L.; Li, C.Q. Weight Function Method to Determine Stress Intensity Factor for Semi-Elliptical Crack with High Aspect Ratio in Cylindrical Vessels. Eng. Fract. Mech.; 2013; 109, pp. 138-149. [DOI: https://dx.doi.org/10.1016/j.engfracmech.2013.05.014]

39. Ayhan, A.O.; Kurt, E. Three-Dimensional Mixed-Mode Stress Intensity Factors for Deflected External Surface Cracks in Thin and Midsize-Thick-Walled Spherical Pressure Vessels. Int. J. Press. Vessel Pip.; 2022; 195, 104596. [DOI: https://dx.doi.org/10.1016/j.ijpvp.2021.104596]