Based on the unequal-length crack interference method above, the simulation research of unequal-length crack interference was carried out, and the characteristic quantity fluid pressure load P in the process of crack propagation, the peak value of von Mises stress distribution field in the crack propagation path, the difference value ∆Bx between the peak value of magnetic induction component at the crack, and the value at its symmetric position were compared and analyzed. On this basis, the crack interaction ratio factor was introduced to study the influence of crack spacing on crack interference and crack size on crack interference.
3.1. Interference Phenomenon of Unequal Length Cracks
For the problem of ferromagnetic pipeline welds with double cracks, in order to study the interference effect of unequal-length double cracks, a model of unequal-length double cracks in pipeline welds was established, as shown in
Figure 2a. The double crack consists of a crack in the outer wall of the pipeline weld and a crack in the inside of the pipeline weld. In order to compare the interference effect, a single crack model on the outer wall of the pipeline weld was established at the same position, as shown in
Figure 2b, and another single crack model inside the pipeline weld was established, as shown in
Figure 2c. Three models of cracks are set in
Figure 2, with the yellow lines representing cracks. The two unequal-length cracks in
Figure 2a are located at the same position and have the same length as the single cracks in
Figure 2b,c, that is, the initial length of the outer wall crack is
lo =
lo′ = 2 mm, and the initial length of the internal crack is
li =
li′ = 4 mm, and the initial distance between the crack tips of two unequal-length cracks is
s = 2 mm. In
Figure 2a, the crack tips are indicated by
T1,
T2, and
T3; in
Figure 2b, the crack tips are indicated by
T1′; and in
Figure 2c, the crack tips are indicated by
T2′ and
T3′.
EN represents the finite element,
ENT1 represents the number of finite elements extended by the crack tip
T1, and so on. When meshing, the discrete size of the unit on the crack propagation path was 0.25 mm.
The von Mises stress during the propagation process of the outer wall single crack (
p = 20.7534 MPa), the internal single crack (
p = 21.3534 MPa), and the unequal-length double cracks (
p = 18.3534 MPa) in the loading step before the crack propagation were extracted. Their distribution cloud maps are shown in
Figure 3a,
Figure 4a and
Figure 5a, and the von Mises stress distribution fields of their corresponding crack propagation path are shown in
Figure 3b,
Figure 4b and
Figure 5b. It can be seen from
Figure 3 that the peak value of the von Mises stress distribution field of the crack propagation path is located at 2 mm from the outer wall, which is at the crack tip
T1′ of the single crack on the outer wall. Additionally, it can be seen from
Figure 4 that the two peaks of the von Mises stress distribution field of the crack propagation path are located at 4 mm and 8 mm from the outer wall, which are the crack tips
T2′ and
T3′ of the internal single crack. In
Figure 5, it appears that the three peaks of the von Mises stress distribution field of the crack propagation path are located at 2 mm, 4 mm, and 8 mm from the outer wall, which are the double crack tips
T1,
T2, and
T3. From
Figure 3,
Figure 4 and
Figure 5, it can be seen that the place where the von Mises stress is the largest is at the crack tip; the distance between the peak value of the von Mises stress distribution field and the internal or outer wall of the pipeline is consistent with the distance from the crack tip to the internal or outer wall of the pipeline. Therefore, the position of the single crack or double crack tips can be judged from the peak position of the von Mises stress distribution field.
Figure 6 shows the required pressure load for the crack tips corresponding to the outer wall single crack, internal single crack, and double crack to extend by one finite element (
EN), respectively. It can be seen from the figure that due to the interference effect of the internal crack, the pressure load required for the outer wall crack tip
T1 to extend by 1
EN in the double cracks is reduced by 21.0534 − 18.6534 = 2.4 MPa compared to the pressure load required for the outer wall single crack tip
T1′ to extend by 1
EN. Similarly, due to the interference effect of the outer crack, the pressure load required for the internal crack tip
T2 to extend by 1
EN is reduced by 21.6534 − 18.8715 = 2.7819 MPa compared to the pressure load required for the internal single crack tip
T2′ to extend by 1
EN. Additionally, the pressure load required for the crack tip
T3 to extend by 1
EN in the double cracks is reduced by 21.8775 − 19.0113 = 2.8662 MPa compared to the pressure load required for the internal single crack tip
T3′ to extend by 1
EN. Judging from the pressure load
P required for crack propagation, the interference effect between the double cracks intensifies the crack propagation process compared with a single-crack case.
The von Mises stress distribution field of the crack propagation paths for the outer wall single crack, internal single crack, and unequal-length double cracks under the same pressure load
p = 18.3534 MPa without any propagation is shown in
Figure 7. It can be seen from the figure that under the same pressure load
p = 18.3534 MPa, due to the interference effect of the internal crack, the peak von Mises stress at the outer wall crack tip
T1 of the double cracks increases by 851.96 − 651.05 = 200.91 MPa compared to the peak von Mises stress at the outer wall single crack tip
T1′. Because of the interference effect of the outer wall crack on the internal crack, the von Mises stress peak value of the crack tip
T2 in the double cracks increases by 817.03 − 691.5 = 125.53 MPa compared to the peak von Mises stress at the internal single crack tip
T2′. The peak von Mises stress at the crack tip
T3 increases by 782.97–693.4 = 89.57 MPa compared to the peak von Mises stress at the internal single crack tip
T3′. Judging from the von Mises stress at the crack tip, compared with the single crack, the interference effect between the double cracks intensifies the crack propagation process.
The component curves of magnetic induction intensity for the outer wall single crack, internal single crack, and unequal-length double cracks under the same pressure load
p = 18.3534 MPa without any propagation were extracted and are shown in
Figure 8. Since the numerical model only establishes the relevant cracks on the right side, and the previous research results show that the magnetic induction intensity curve is symmetrical when there is no crack in the pipeline weld, the characteristic quantity ∆
Bx is defined as the difference between the peak value of the magnetic induction intensity component at the crack and its symmetrical position [
24,
25]. In this study, there are, in total, three ∆
Bx characteristic quantities, with ∆
Bxop representing the outer wall single crack, ∆
Bxip for the internal single crack, and ∆
Bxdp for the double cracks. In
Figure 8, ∆
Bxop = 0.0065 T, ∆
Bxip = 0.0025 T, ∆
Bxdp = 0.0246 T. The interference effect between the double cracks leads to the superposition of the leakage magnetic field of the outer wall crack with that of the internal crack. Furthermore, it leads to the detected peak value of the magnetic induction intensity component of the double cracks increasing by ∆
Bxdp − ∆
Bxop = 0.0181 T compared to the peak value of the magnetic induction intensity component of the single crack in the outer wall. So, the detected peak value of the magnetic induction intensity component of the double cracks increases by ∆
Bxdp − ∆
Bxip = 0.0221 T compared to the peak value of the magnetic induction intensity component of the internal single crack. Judging from the peak value of the magnetic induction intensity component, compared with the single crack, the interference effect between the double cracks intensifies the crack propagation process. It can be seen that when the double cracks interfere with each other, the cracks are more likely to propagate toward each other.
Since the peak position of the von Mises stress distribution field is at the crack tip, the position and number of cracks can be judged. Extracting
s = 2 mm during the crack propagation process of the double crack in the
Figure 2a above, the von Mises stress distribution field of the crack tip
T1 extending by one, two, three, and four finite elements, that is, the von Mises stress distribution fields of
ENT1 = 1,
ENT1 = 2,
ENT1 = 3, and
ENT1 = 4, are depicted in
Figure 9a,
Figure 10a,
Figure 11a and
Figure 12a. In
Figure 9b,
Figure 10b,
Figure 11b and
Figure 12b, the length of the outer wall crack is expressed by
loe, and the length of the internal crack is expressed by
lie. The distance between the crack tip
T1 and the crack tip
T2 is denoted by
se. Grid reconstructions are shown in
Figure 9c,
Figure 10c,
Figure 11c and
Figure 12c. It is seen from
Figure 9,
Figure 10,
Figure 11 and
Figure 12 that as the crack propagates, the peak position of the von Mises stress distribution field also changes, and there is a one-to-one correspondence. For example, in
Figure 12b, when
ENT1 = 4 and
ENT2 = 3, the distance between the crack tip
T1 and the crack tip
T2 is a finite element,
se = 0.25 mm, while in
Figure 12a, the corresponding peak position coordinates of the von Mises stress distribution field at this time are (15.752, 1741.84) and (16.002, 1856.01), and the distance between the two peaks is 16.002–15.752 = 0.25 mm. The von Mises stress distribution fields of
ENT1 = 1,
ENT1 = 2,
ENT1 = 3, and
ENT1 = 4 are summarized in
Figure 13. This figure shows that when the crack tip begins to propagate, the von Mises stress value at the crack tip increases as the crack propagation process intensifies. From
Figure 13, the crack propagation process and position can be monitored more clearly.
Figure 14a–d show the magnetic induction intensity nephograms when extracting crack tip
T1 extending by one, two, three, and four finite elements, and
Figure 14e shows the magnetic induction intensity component curves when extracting crack tip
T1 extending by one, two, three, and four finite elements. With the increase in the crack propagation process, the greater the peak value of the von Mises stress at the crack tip, the greater the deformation at the crack position, which affects the magnetic field distribution of the pipeline weld during the crack propagation process, so the peak value of the detected magnetic induction intensity component shows an increasing trend.
3.2. Effect of Crack Spacing on Crack Interference
In order to analyze the influence of crack spacing on crack propagation, finite element models with initial crack spacings
s = 2 mm, 4 mm, 6 mm, and 8 mm were established, respectively.
Figure 15a–c, respectively, show the model diagrams of the crack tip
T1 extending by four finite elements when
s = 4 mm,
s = 6 mm, and
s = 8 mm. When
s = 2 mm, the model of crack tip
T1 extending by four finite elements is depicted in
Figure 12b.
For different distances, the fluid pressure required is calculated when the double crack tip
T1 extends by one, two, three, and four finite elements. That is, at distances
s = 2 mm, 4 mm, 6 mm, and 8 mm, respectively, the required fluid pressures are given in
Table 2 and
Figure 16 when
ENT1 = 1,
ENT1 = 2,
ENT1 = 3, and
ENT1 = 4. In
Table 2 and
Figure 16, to compare the influence of crack spacing on crack interference, the fluid pressures required for the expansion of a single outer wall crack
T1′ extending to one, two, three, and four finite elements are also listed. The interference effects at different spacing from the required pressures can also be found in
Table 2 and
Figure 16.
When the double-crack distances are s = 2 mm, 4 mm, 6 mm, and 8 mm, respectively, under the same pressure load without any propagation (taking
p = 18.3534 MPa as an example), the von Mises stress distribution field is extracted, which is shown in
Figure 17. To compare the interference effect with various crack spacing, the von Mises stress distribution field of the outer wall single crack is also listed.
Figure 17a–c show the von Mises stress distribution field at
s = 4 mm,
s = 6 mm, and
s = 8 mm, respectively. Additionally, the von Mises stress distribution field at
s = 2 mm can be viewed from
Figure 7.
Under the same pressure load without any propagation (taking
p = 18.3534 MPa as an example), the magnetic induction intensity component curves of double cracks
s = 2 mm, 4 mm, 6 mm, and 8 mm are shown in
Figure 18. To compare the interference effect with various crack spacing, the magnetic induction intensity component curve of the outer wall single crack is also listed (
Figure 8 shows the case for
s = 2 mm already).
Based on
Table 2 and
Figure 16, the required fluid pressure increases with the initial spacing
s of the double cracks. From
Figure 7 and
Figure 17, under the same pressure load, with the increase in the different initial distances between the double cracks, the peak value of the von Mises stress at the crack decreases. It can be seen from
Figure 18 that with the increase in the different initial spacing
s of the double cracks, the difference between the peak values of the detected magnetic induction intensity components of the double cracks and the values of their symmetrical positions decreases. In order to describe the effect of crack spacing on crack interference more conveniently, the crack interaction scale factors
,
, and
are introduced.
is the ratio of the pressure required for the outer wall single crack to the pressure required for the double cracks when extending by the same number of finite elements;
is defined as the ratio of the peak value of von Mises stress at the crack tip
T1 of the double cracks to the peak value of the von Mises stress at the crack tip
T1′ of the outer wall single crack under the same pressure load;
is defined as the ratio of the peak value difference of the magnetic induction intensity component of the double crack to the peak value difference of the magnetic induction intensity component of the single crack in the outer wall under the same pressure load. The expressions are:
where
Po is the pressure required for the outer wall single crack, in units of MPa;
Pd is the pressure required for the double cracks, in units of MPa.
where Mises
T1 is the peak von Mises stress of the crack tip
T1 of double cracks, in units of MPa; Mises
T1′ is the peak von Mises stress of the crack tip
T1′ of the outer wall single crack, in units of MPa.
where ∆
Bxdp represents the peak value difference of the magnetic induction intensity component of the double cracks, in units of T; ∆
Bxop is the peak value difference of the magnetic induction intensity component of the outer wall single crack, in units of T.
Table 3 is obtained from
Table 2;
Table 4 is obtained from
Figure 7 and
Figure 17. According to the extraction method of the characteristic quantity ∆
Bx value in
Figure 8, the value of ∆
Bxdp when
s = 2 mm, the value of ∆
Bxop of the outer wall single crack, and the value of ∆
Bxdp when s = 4 mm, 6 mm, and 8 mm are extracted from
Figure 18 and summarized in
Table 5.
It can be seen from
Table 3 that when
T1 and
T1′ are extended by the same number of finite elements, the value of
decreases as the crack spacing increases. That is, the gap between the pressure required for the outer wall single crack and the pressure required for the double cracks decreases. Taking
ENT1′ = 3 and
ENT1 = 3 for example, when
s changes from 2 mm to 4 mm to 6 mm to 8 mm,
changes from 1.1623 to 1.0918 to 1.0283 to 1.0036. It can be seen from
Table 4 that under the same pressure load, as the crack spacing increases, the value of
decreases. That is, the difference between the peak von Mises stress of the crack tip
T1 of the double crack and the peak value of the von Mises stress of the crack tip
T1′ of the outer wall single crack decreases. It can be seen from
Table 5 that under the same pressure load, as the crack spacing increases, the value of
decreases. That is, the peak difference between the magnetic induction intensity component of the double crack and the magnetic induction intensity component peak difference of the outer wall single crack decreases. As the crack spacing increases,
,
, and
decrease, and the interference effect between the double cracks becomes smaller and smaller. When the crack spacing
s = 8 mm (double cracks model: the initial length of the outer wall crack
lo = 2 mm, the initial length of the internal crack
li = 4 mm), the values of
,
, and
are close to 1, and the interference effect of the double cracks is negligible. The interference trend of the double cracks gradually weakens with the increasing tip distance and finally tends to the situation of the single crack. Thus, when conducting a safety assessment for multicrack oil and gas pipelines, it is possible to directly simplify the multicrack treatment to the single crack treatment without considering the interaction between cracks.
3.3. Effect of Crack Size on Crack Interference
Through the analysis in
Section 3.1, compared with the single-crack case, the interference effect between double cracks intensifies the crack propagation process. Three characteristic values are defined as the pressure
P required for crack propagation, the peak value of von Mises stress, and the difference ∆
Bx between the peak value of the magnetic induction intensity component at the crack and the value of its symmetrical position. These characteristic values can be used to measure the progress of crack propagation. According to the analysis in
Section 3.2, the double cracks’ interference effect weakens as the crack tip distance increases. To further study the interference effect of the unequal-length double cracks, two sets of the above three characteristic quantities are established and extracted. For the first set, the initial crack length in the outer wall is
lo =
lo′= 2 mm, the initial length of the internal crack is
li =
li′ = 6 mm, and the initial distance between the crack tips of two unequal-length cracks is
s = 2 mm. For the second numerical example, the initial length of the outer wall crack is
lo =
lo′ = 2 mm, the initial length of the internal crack is
li =
li′ = 8 mm, and the initial distance between the crack tips of two unequal-length cracks is
s = 2 mm. The model diagrams of the first calculation example are shown in
Figure 19a–c, and the model diagrams of the second calculation examples are shown in
Figure 19d–f. By comparing the crack interaction scale factors (
,
,
), the interference effect of large-size cracks on small-size cracks and the interference effect of small-size cracks on large-size cracks are studied. Equation (4) is obtained from Formula (1), that is, when extending by the same number of finite elements, the pressure ratio required for the internal single cracks and the double cracks is expressed as:
where
Pi is the pressure required for the internal single crack, in units of MPa, and
Pd is the pressure required for the double cracks, in units of MPa.
The pressure
P required to extract the characteristic quantity crack propagation process is listed in
Table 6 and
Table 7.
Table 8 is obtained from Formulas (1) and (4) and
Table 6 and
Table 7.
By comparing the interference scale factors
and
in
Table 8, it can be seen from the first set of calculation examples (
lo =
lo′
= 2 mm,
s = 2 mm,
li =
li′
= 6 mm) that when the crack tip extends by one finite element (
ENT1 = 1,
ENT1′ = 1 or
ENT2 = 1,
ENT2′ = 1 or
ENT3 = 1,
ENT3′ = 1), due to the interference of the 6 mm crack on the 2 mm crack,
(
ENT1 = 1,
ENT1′ = 1), namely, the extension pressure
Po required for the 2 mm single-crack case is 1.3197 times as large as the expansion pressure
Pd required for the 2 mm-and-6 mm double-crack case. Owing to the interference of the 2 mm crack on the 6 mm crack,
(
ENT2 = 1,
ENT2′ = 1) and
(
ENT3 = 1,
ENT3′ = 1), and the required propagation pressure
Pi for the 6 mm single-crack case is 1.1742 times (
ENT2 = 1,
ENT2′ = 1) and 1.1805 times (
ENT3 = 1,
ENT3′ = 1) as large as
Pd required for the 2 mm-and-6 mm double-crack case. Therefore, the change range of the crack tip propagation pressure of the 2 mm crack under the interference of the 6 mm crack is larger than that of the crack tip propagation pressure in the opposite case. In other words, the proportional factor
; the interference effect of 6 mm crack on 2 mm crack is more severe than that of the 2 mm crack on the 6 mm crack. This phenomenon can still be observed when the crack tip extends by two finite elements.
For the second numerical example with lo = lo′ = 2 mm, s = 2 mm, and li = li′ = 8 mm, when the crack tip extends by one finite element (ENT1 = 1, ENT1′ = 1 or ENT2 = 1, ENT2′ = 1 or ENT3 = 1, ENT3′ = 1), (ENT1 = 1, ENT1′ = 1), (ENT2 = 1, ENT2′ = 1), and (ENT3 = 1, ENT3′ = 1); the interference effect of 8 mm crack on 2 mm crack is more severe than that of 2 mm crack on 8 mm crack. Based on the two numerical examples, it is found that the interference factor () of 8 mm crack on the 2 mm crack is times as large as the interference factor () of the 6 mm crack on the 2 mm crack. The interference factor (ENT2 = 1, ENT2′ = 1) and (ENT3 = 1, ENT3′ = 1) of the 2 mm crack on 8 mm crack is (ENT2 = 1, ENT2′ = 1) times and (ENT3 = 1, ENT3′ = 1) times as large as the interference factor (ENT2 = 1, ENT2′ = 1) and (ENT3 = 1, ENT3′ = 1) of the 2 mm crack on the 6 mm crack. The interference effect of the 2 mm crack on the 6 mm crack and 8 mm crack is not much different. The interference effect of large-size cracks such as 8 mm and 6 mm on small-size cracks such as 2 mm is significantly greater than the influence of the small-size crack such as 2 mm on large-size cracks such as 8 mm and 6 mm. The influence of the 8 mm crack on the 2 mm crack is much greater than that of the 6 mm crack on the 2 mm crack.
In Formula (5), the von Mises stress ratio is defined as, under the same pressure load, the ratio of the von Mises stress peak value at crack tip
T2 of the double crack to that at crack tip
T2′ of the internal single crack. In Formula (6), the von Mises stress ratio is defined as, under the same pressure load, the ratio of the von Mises stress peak value at crack tip
T3 of the double crack to that at crack tip
T3′ of the internal single crack.
Under the same pressure load,
p = 15.6534 MPa is taken as the first example and
p = 13.2534 MPa is taken as the second example. The von Mises stress distribution field is shown in
Figure 20 and
Figure 21. The von Mises stress peak values are extracted from
Figure 20 and
Figure 21. The results are listed in
Table 9 and
Table 10 by applying Formulas (2), (5) and (6).
In
Table 9 and
Table 10, for the first example with
lo =
lo′
= 2 mm,
s = 2 mm, and
li =
li′
= 6 mm, under the same pressure load
p = 15.6534 MPa, the interference scale factors including
,
, and
are compared. Owing to the interference of the 6 mm crack on the 2 mm crack,
1.3897, i.e., the von Mises stress peak value at crack tip
T1 for the double-crack case (2 mm and 6 mm, respectively) is 1.3897 times as large as the peak value at the crack tip
T1′ when only the single 2 mm crack exists. Owing to the interference of the 2 mm crack on the 6 mm crack,
1.2028 and
1.1254, i.e., the von Mises stress peak value at the crack tip
T2 for the double crack case is 1.2028 times as large as the one at crack tip
T2′ when only the single 6 mm crack exists. The von Mises stress peak value at the crack tip
T3 for the double crack case is 1.1254 times as large as the one at crack tip
T3’ when only the single 6 mm crack exists. Therefore, the change amplitude of the von Mises stress peak at the crack tip under the interference of the 6 mm crack to the 2 mm crack is larger than that of the other way around. The scaling factor
, and the interference effect of 6 mm crack to 2 mm crack is stronger than that of the 2 mm crack to 6 mm crack. This phenomenon is more obviously observed in the second example, where
lo =
lo′
= 2 mm,
s = 2 mm, and
li =
li′
= 8 mm. Under the same pressure load
p = 13.2534 MPa,
1.4138,
1.2130, and
1.1262, that is, the proportional factor
, and the interference effect of the 8 mm crack to 2 mm crack is more severe than the other way around.