Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 8049318, 12 pages
https://doi.org/10.1155/2018/8049318

Research Article

Interaction between a Macrocrack and a Cluster of Microcracks by Muskhelishvili’s Complex Potential Method

Xu Li, Xiaotao Li, Hongda Yang, and Xiaoyu Jiang

Applied Mechanics and Structure Safety Key Laboratory of Sichuan Province, School of Mechanics and Engineering, Southwest Jiaotong University, Chengdu 610031, China

Correspondence should be addressed to Xiaoyu Jiang; xiaoyujiang8@sina.com

Received 22 July 2018; Accepted 17 September 2018; Published 9 October 2018

Academic Editor: Paolo Lonetti

Copyright © 2018 Xu Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The interaction between a macrocrack and a cluster of microcracks has been investigated based on Muskhelishvili’s complex potential method. A step-by-step subproblem procedure is used to satisfy the stress boundary conditions on each crack surface. The interactions between a cluster of microcracks and a macrocrack and the interaction among microcracks are analyzed. Three damage configurations as chained, reverse-chained, and randomly distributed microcracks have been designed to simulate the damage around the macrocrack tip. The solution of an infinite elastic plane containing a macrocrack and a cluster of microcracks is presented for the plane subjected to a uniform tensile load. The stress intensity factor (SIF) at the macrocrack tip and the microcrack tips is obtained. The results show that the inclination angle of the microcrack and the distance between the macrocrack and microcracks have a great influence on SIF. When the inclination angle is small, the SIF at microcrack tips may be larger than other inclination angles. These results are helpful to analyze the fracture or damage behaviors of materials.

1. Introduction

Many experimental investigations demonstrate that profuse microcrack propagations can be observed during the evolution of fracture process in the high-stress concentration regions, as well at the tip of the macrocrack [1–4]. The microcrack has a great effect on SIF at the macrocrack tip. Thus, the interaction between a macrocrack and multiple microcracks is a topical problem in elastic theory and fracture mechanics.

Many results concerning the interaction of macrocrack and arrays of microcracks were obtained using various numerical methods such as serious expansions, boundary element methods, and collocation methods. Crack propagation mechanisms can be developed if the basic question whether the microcracks enhance or suppress the macrocrack propagation is answered. The pseudotractions method has been developed by many studies [5–9] to investigate the interaction problem between main crack and microcracks. The core idea of this method was that a superposition of a number of subproblems each contains one crack loaded by unknown tractions. Kachanov [10–12] proposed the interaction matrix between cracks in the theory. He also replaced the actual traction by average traction. Based on small parameter method, Tamuzs et al. [13, 14] investigated the interaction of a macrocrack with a cluster of cracks; Petrova et al. [15–17] effectively applied it in thermal fracture of a functionally graded/homogeneous bimaterials with series of cracks. Hoagland [18] and Ismail [19] studied the process of the nucleation and interaction of microcracks near a macrocrack tip by using the complex potential method and the theory of dislocation. Yan Xiangqiao [20, 21] presented a numerical approach for the interaction problem of macrocrack with microcrack in an infinite elastic plane. Based on the principle of superposition and a flattening method, Xia Xiaozhou et al. [22] analyzed the interaction between cracks and effect of microcrack zone on main crack tip. Rafiee [23] and Seelig [24] investigated the interaction between cracks and branching of fast running cracks.

However, microcracks are mainly arrays in parallel and collinear in previous studies. In fact, the formation of microcracks is random, and randomly distributed microcracks in front of a macrocrack tip are usually used to model the material damage in some materials such as metals [25–27]. In this paper, three kinds of damage configurations as chained, reverse-chained, and randomly distributed microcracks have been designed to simulate the damage around the macrocrack tip.

The outline of this paper is as follows. In Section 2, a high efficient method based on Muskhelishvili’s complex potential method is presented to study the interaction between a finite macrocrack and a cluster of microcracks and the interaction among microcracks. A step-by-step subproblem procedure is used to satisfy the boundary conditions on each crack surface. The solution of an infinite plane that contains a macrocrack and arbitrary located and oriented microcracks is presented. In Section 3, numerical examples are employed to illustrate the accuracy of the proposed method. The influence of three kinds of arrangements of microcracks on SIF in the macrocrack tip is analyzed. In Section 4, some conclusions are drawn.

2. Formulation

2.1. Problem Description and Solution Scheme

The problem in terms of an infinite elastic plane containing a macrocrack and randomly distributed microcracks under remote uniform traction T is considered, as shown in Figure 1. In order to a convenient statement, the microcrack from 1 to is called microcrack 1 to microcrack k in order. In Figure 1, oxy is the global coordinate system with the origin at the center of the macrocrack. is the local coordinate system with the origin at the center of microcrack k.

Figure 1: An infinite plane containing a macrocrack and multiple microcracks under mode I load.

The problem in Figure 1 can be replaced by a series of subproblems, as shown in Figure 2. Each subproblem is independent.

(a) Subproblem (i)

(b) Subproblem (ii)

(c) Subproblem (iii)

Figure 2: The transformation of the original problem into a series of subproblems.

With reference to Figure 2, is the distance between the macrocrack tip and the center of microcrack k; is the inclination angle to the x-axis of microcrack k; is the direction with respect to the macrocrack tip. In addition, a is the half-length of the macrocrack, and is the half-length of microcrack k.

The subproblems are as follows: In subproblem (i), an infinite plane containing a finite macrocrack, without the microcrack, subjected to uniaxial tension, refer to Figure 2(a). In subproblem (ii), an infinite plane containing an arbitrarily oriented microcrack, without the macrocrack, subjected to the tractions obtained in subproblem (i) but opposite in sign. There are n microcracks corresponding to n subproblems in subproblem (ii), refer to Figure 2(b). In subproblem (iii), an infinite plane containing a finite macrocrack, without the microcrack, subjected to the tractions obtained in subproblem (ii) but opposite in sign, refer to Figure 2(c). Subproblem (iv) and subproblem (v) may be obtained through a similar process of subproblem (ii) and subproblem (iii) and so on.

It should be noted that the original problem is under remote uniform traction T, as shown in Figure 1. The stress boundary condition on each crack surface is traction free. The stress boundary condition is satisfied by the step-by-step subproblems. For example, the stresses along the line which is the location of microcrack j in subproblem (i) are obtained. Then the opposite stresses are applied to the crack surface of microcrack j in subproblem (ii). At this time, the stresses on the crack surface of microcrack j are zero. But the stresses on the crack surface of microcrack j are not zero because of the affection of other microcracks in subproblem (ii). And the nonzero value of stress on the crack surface of microcrack j is a small order of magnitude. The nonzero value of stress on the crack surface of microcrack j is added to the stress from subproblem (iii), and then the resultant stress (opposite in sign) is applied to the microcrack j in subproblem (iv), and so on. In the end, the stress boundary condition of traction free on each crack surface is satisfied.

2.2. Mathematical Formulation of Subproblems

For plane elasticity, the stress fields of a single crack in the center of an infinite plane can be expressed in terms of two complex potentials and as [28]

where ; denotes the normal stress in the horizontal direction; denotes the normal stress in the vertical direction; denotes the shear stress; for plane strain, for plane stress; is shear modulus of the material; is Poisson's ratio; and are the first derivative and the second derivative of z, respectively.

2.2.1. Subproblem (i): The Macrocrack

In subproblem (i), the corresponding stress functions may be given by [28]

where T is the applied loading as shown in Figure 1. Substituting (2) into (1), the well-known stress fields near the macrocrack tip may be obtained

For the nonexistent microcrack k, substituting the position parameters , , , and of the microcrack into (3), the stresses , , and along the position of the microcrack in the global coordinate system are obtained, the superscript (i) on the stress and denotes the results in subproblem (i), and the subscript (k) denotes microcrack k.

By coordinate transformation, the normal and tangential stresses along the position of microcrack k in the local coordinate system may be obtained by

where the subscripts N and T denote the normal and tangential direction of the crack surface, respectively.

2.2.2. Subproblem (ii): The Microcracks

In subproblem (ii), the crack surface is subjected to a nonuniform distributed load that is equal in magnitude and opposite in sign to the stresses obtained in subproblem (i). To solve subproblem (ii), it is necessary to introduce the problem shown in Figure 3. In the coordinate system , a pair of concentrated forces acts at a point of the crack surface in an infinite elastic plane, where P is in the normal direction and Q is in the tangential direction.

Figure 3: A crack subjected to concentrated forces. (a) ; (b) .

Then, introduce the mapping functions (see (5)), where is the inverse function of z, is another complex variable, , , and a is the half-length of the crack. The crack surface in the is transformed into a unit circle in the transformed plane, and is corresponding to .

Therefore, the stress functions of the concentrated forces applied to the upper and lower boundary of the crack will have the following forms [28]:

With (5), the point in is transformed to the point in the plane, and the basic equation (6) may be rewritten in terms of the transformed and .

Thus, in subproblem (ii), when a pair of concentrated forces is applied at the edge of the microcrack j in the coordinate system , where

P and Q are equal in magnitude and opposite in sign to the nonuniform stresses and obtained in subproblem (i) along the line which is the presumed location of the microcrack j.

The stresses (, , , , , and ) along the position of the macrocrack and the microcrack can be calculated on the basis of (1), (5), (6), and (7). The superscript (ii) on the stress and denotes the results in subproblem (ii), and the subscript F denotes the results under concentrated forces F.

Then, when the nonuniform stresses act on the crack surface of the microcrack j in an infinite plane, the stress fields of the macrocrack and microcrack in the local coordinate system are obtained by

For the macrocrack, the subscript (j) denotes the effect of microcrack j on the macrocrack. For the microcrack, the subscript (kj) denotes the effect of microcrack j on microcrack k. At this time, the stress on the crack surface of microcrack j is zero because of the superposition of subproblem (i) and subproblem (ii).

Finally, transform the local coordinate system to the coordinate system of macrocrack oxy and microcrack , the solution (, , , and ) may be obtained for a nonuniform distributed load acting on the crack surface of microcrack j in an infinite elastic plane. In the same way, we may use the formulas to derive the rest solutions in subproblem (ii).

There are n microcracks corresponding to n subproblems in subproblem (ii), the whole traction components of the macrocrack and microcracks are given by

It is shown from (10) that the stress on the crack surface of each microcrack is not zero because of the affection of other microcracks in subproblem (ii). But the nonzero value of stress on the crack surface of each microcrack is a small order of magnitude compared with the stresses and obtained in subproblem (i).

2.2.3. Subproblem (iii): The Macrocrack

It is needed to study the problem of an infinite plane containing a macrocrack and without the microcracks. The crack surface is subjected to a nonuniform distributed load that is equal in magnitude and opposite in sign to the stresses in (9). It is used to solve the problem that the same method as described in subproblem (ii). When the concentrated forces F acting at the upper and lower surface of the macrocrack, the stress fields , , and of the nonexist microcrack can be obtained; the expression of F is

The stress components along the position of the microcrack are

Finally, the normal and shear stresses and in the coordinate system are obtained by the coordinate transformation from oxy to . In the end, the results and may be obtained for microcracks in subproblem (iii).

2.2.4. Subproblem (iv): The Microcracks

and are added with (10) to obtain the load on microcrack surface as and in subproblem (iv). The results , , , and may be obtained. The results are similar to (9) and (10).

In the following section, the numerical results for the interaction effect between the macrocrack and an array of microcracks and the interaction effect among microcracks have been presented.

By the aforementioned method, the interactions between the cracks are reflected in every subproblem, for example, in subproblem (i), the macrocrack effect on the microcrack; in subproblem (ii), the microcrack effect on the macrocrack and the interaction effect among the microcracks; in subproblem (iii), the macrocrack effect on the microcrack once again. Besides, the more the decomposition of the subproblems is, the higher the calculation accuracy can be obtained.

2.3. The Solution of KI and KII

As mentioned above, the traction-free boundary condition on the crack surface is satisfied by a step-by-step subproblem procedure. The results of each subproblem satisfy the fundamental equations of linear elastic mechanics independently.

In subproblem (i), the load on the macrocrack surface being and , and the SIFs at the tip of the macrocrack are

In the step-by-step procedure, the loads applied on the macrocrack surface are and in subproblem (i), and in subproblem (iii), and in subproblem (v), etc. The SIFs at the tip of the macrocrack may be written as

The loads applied on each microcrack surface are and in subproblem (ii), and in subproblem (iv), and in subproblem (vi), etc. The SIFs at the tip of the microcrack may be written as

3. Results and Discussion

In order to understand the macrocrack propagating and the microcracks getting together, it is designed to simulate the damage at the macrocrack tip that three distributions of microcracks as chained arranged microcracks, reverse-chained arranged microcracks and randomly oriented microcracks. At the first, the correctness of the step-by-step procedure should be verified.

3.1. Computational Accuracy Verification

An infinite plane containing two collinear cracks with the same length under mode I load is considered in Figure 4(a), where . Crack AB may be considered as the macrocrack in Section 2. The variation of normalized SIF at crack tip B versus distance is depicted in Figure 4(b). The results obtained by the proposed theoretical method are in good agreement with the corresponding theoretical results [29, 30] and FEM results. The finite element model is established by ANSYS 12.0. It is observed that the normalized SIF at the crack tip B increase with the decrease of the distance between cracks.

(a)

(b)

Figure 4: (a) Two collinear cracks under tensile load; (b) versus for under tensile load.

3.2. Interaction between a Macrocrack and a Cluster of Chained Arranged Microcracks

It is considered that an infinite plane containing a macrocrack of length and a cluster of chained arranged microcracks of length as shown in Figure 5. The uniform tensile load T is applied at infinity perpendicular to the macrocrack. The length of each microcrack is taken as . The fourteen microcracks are symmetrical about the x-axis. Some configurations of microcrack 1 to microcrack 7 (as in Figure 2) are given in Tables 1–3.

Table 1: The direction angle of each microcrack.

20
25
15
20
15
10
5

Table 2: The distance between the macrocrack and each microcrack.

2
5
8

3
6
9

4
7
10

5
8
11

Table 3: The inclination angle of each microcrack.

10
5
-5
2.5
-2.5
2.5
-2.5

20
10
-10
5
-5
5
-5

30
15
-15
7.5
-7.5
7.5
-7.5

40
20
-20
10
-10
10
-10

50
25
-25
12.5
-12.5
12.5
-12.5

Figure 5: An infinite plane with a macrocrack and fourteen microcracks.

The nondimensional parameters for the macrocrack and each microcrack are taken as

In Figure 6, (+1) and (-1) represent the right and left crack tips respectively, where may be considered as the distance parameter between the macrocrack and each microcrack, and may be considered as the inclination parameter of the cluster of microcracks from the presupposed relation of each crack angle in Table 3.

(a)

(b)

Figure 6: Normalized SIF at the right macrocrack tip versus inclination parameter for different distance parameter of between the cluster of microcracks and the macrocrack as in Tables 2-3. (a) ; (b) .

Due to the influence of microcracks, all values of at the macrocrack tip are greater than 1, all values of at the macrocrack tip are not equal to zero. In Figure 6(a), when is lower than 20 degrees, the strong amplification effects (which means the normalized SIF is larger than 1, the interaction effect between cracks enlarges the SIF at the crack tip) are observed at the macrocrack tip. The SIF decreases with the increase of inclination parameter and distance parameter . On the contrary, shielding effect means the normalized SIF is less than 1, and the interaction effect between cracks decreases the SIF at the crack tip. In Figure 6(b), the normalized is small relative to the normalized and its influence is very small.

Due to the symmetry, only the normalized SIF over the x-axis is shown in Figure 7, where . In Figure 7(a), when is between -30 and 20 degrees, the value of SIF is near or over 1, that means the microcracks act an amplification effect each other, and this is due to the interaction effect among microcracks. When is more than 20 degrees, a change from the amplification to the shielding occurs.

(a)

(b)

Figure 7: Normalized SIF at the left tip of each microcrack versus the inclination angle of each microcrack, and 1 to 7 represent the microcracks in Figure 5. (a) ; (b) .

The results in Figure 7(b) show the SIF begins to decline from the peak when is more than 20 degrees. Moreover, the normalized and of the microcracks have the same order of magnitude; it means that the microcracks are mixed type crack of mode I and mode II.

3.3. Interaction between a Macrocrack and a Cluster of Reverse-Chained Arranged Microcracks

An infinite plane containing a macrocrack of length and a cluster of reverse-chained arranged microcracks of length is considered in Figure 8. The tensile stress T is applied at infinity perpendicular to the macrocrack. The length of each microcrack is taken as . Some configurations of microcracks (as in Figure 2) are given in Tables 4–6.

Table 4: The direction angle of each microcrack.

40
30
20
10
15
10
5

Table 5: The distance between the macrocrack and each microcrack.

2
5
8

3
6
9

4
7
10

5
8
11

Table 6: The inclination angle of each microcrack.

2.5
-2.5
2.5
-2.5
5
-5
10

5
-5
5
-5
10
-10
20

7.5
-7.5
7.5
-7.5
15
-15
30

10
-10
10
-10
20
-20
40

12.5
-12.5
12.5
-12.5
25
-25
50

Figure 8: An infinite plane with a macrocrack and fourteen microcracks.

The results in Figure 9 show that the effect of distance on SIF at the macrocrack tip has a similar behavior to Figure 6; it will be getting weaker as the distance between microcrack and macrocrack increases. However, the value of SIF is larger than that in Figure 6, due to the more concentrated microcracks in front of the macrocrack tip.

(a)

(b)

Figure 9: Normalized SIF at the right macrocrack tip versus inclination parameter for different distance parameter of between the cluster of microcracks and the macrocrack as in Table 5. (a) ; (b) .

The normalized SIF at the microcrack tip over the x-axis is shown in Figure 10, where . The greater the number of microcracks at the macrocrack tip, the bigger the value of normalized . The microcracks are mixed type crack of mode I and mode II because of the normalized and of the microcracks being the same order of magnitude.

(a)

(b)

Figure 10: Normalized SIF at the left tip of each microcrack versus the inclination angle of each microcrack, and 1 to 7 represent the microcracks in Figure 9. (a) ; (b) .

3.4. Interaction between a Macrocrack and a Cluster of Randomly Distributed Microcracks

In this section, consider an infinite elastic plane containing a macrocrack and nine randomly distributed microcracks, subjected to uniform tensile load in direction y. Each microcrack is assigned a random orientation and a prospective microcrack site. Figure 11 shows the randomly distributed microcracks in front of the macrocrack tip. Some configurations of microcracks (as in Figure 2) are given in Tables 7–9. The length of each microcrack is taken as . Four cases of inclination angles of microcrack are considered in Table 9.

Table 7: The direction angle of each microcrack.

80
60
10
30
0
-45
-20
-70
-90

Table 8: The distance between the macrocrack and each microcrack.

3.5
2.5
4.5
5.5
3.5
2.5
3.5
4.5
5.5

Table 9: The inclination angle of each microcrack.

case

1
60
0
45
-20
-80
-10
-70
60
30

2
0
70
-60
90
10
-45
-30
-20
80

3
-45
20
0
70
-30
90
-80
10
60

4
90
45
-70
0
20
30
-60
80
-10

Figure 11: Nine randomly distributed microcracks in front of the macrocrack tip.

With reference to (16), Tables 10 and 11 show the normalized SIF of macrocrack and microcracks; these SIFs are , , and .

Table 10: The mode-I SIF at the right macrocrack tip and left microcrack tips 1 to 9.

case

1
0.956
0.06
1.89
0.15
1.29
0.06
1.47
0.03
0
1.77

2
0.995
1.63
0.06
0
0.08
1.52
0.03
0.58
0.64
0

3
0.999
0.11
0.76
1.36
0.08
0.70
0.07
0.08
1.63
0

4
0.971
0.05
0.02
0
1.37
1.23
0.83
0
0
0.79

Table 11: The mode-II SIF at the right macrocrack tip and left microcrack tips 1 to 9.

case

1
0.122
0.11
-0.02
0.51
-0.32
0.006
-0.46
-0.06
0.30
0.43

2
-0.015
0.26
0.06
-0.23
-0.09
0.13
-0.38
-0.71
-0.63
-0.18

3
-0.053
-0.70
0.72
-0.05
0.11
-0.74
-0.08
-0.03
0.38
0.35

4
0.086
-0.16
0.34
-0.06
-0.09
0.56
0.85
-0.16
-0.13
-0.55

The results show the microcracks are mainly mixed type crack of mode I and mode II. Compared to the two chained arranged microcracks as discussed above, the macrocrack was affected little by randomly distributed microcracks; the value of the normalized even decreased (shielding, ). In other words, the concentrated distributed microcracks have greater effect on macrocrack than randomly distributed microcracks.

However, the interaction among microcracks does not decreased obviously in the four cases. The microcrack always acts an amplification effect() at about degrees; the microcrack acts a shielding effect() at about and degrees.

Moreover, it may be predicted that the microcrack is more likely to be propagated when the inclination angle is between -10 and 10 degrees when the load increases or is subjected to fatigue load.

4. Conclusions

In this paper, a high efficient method based on Muskhelishvili’s complex potential method is presented to study the interaction between a finite macrocrack and a cluster of microcracks. A step-by-step subproblem procedure is used to satisfy the stress boundary conditions on each crack surface. The stress intensity factor (SIF) of the macrocrack and the microcracks is obtained. With the comparison of these three kinds of damage configurations around the macrocrack tip, the interactions between a cluster of microcracks and a macrocrack and the interaction among microcracks are analyzed. Some conclusions can be summarized in the following:

The chained and reverse-chained arranged microcracks increase SIF of macrocrack at about degrees, while the randomly distributed microcracks have little effects on the value of SIF of macrocrack. In other words, the concentrated distributed microcracks have greater effect on macrocrack than randomly distributed microcracks.

The effect of microcracks on SIF at the macrocrack tip will be getting weaker as the distance between microcracks and macrocrack increases.

For the chained and reverse-chained arranged microcracks, the SIF at the microcrack tip may be increased when the inclination angle of each microcrack is at about degrees. When the inclination angle at about degrees, the shielding effect appears gradually.

For the randomly distributed microcracks, the SIF at the microcrack tip may be increased when the inclination angle of each microcrack is at about degrees.

The microcracks are mixed type crack of mode I and mode II for the three distributions of microcracks.

As a future work on the subject, some guidelines may be given as follows.

The proposed theoretical method based on the Muskhelishvili’s complex potential theory is a simple tool to solve crack problems, even for the complex cases such as a macrocrack interacting with a cluster of nine microcracks. It is worth noting that the proposed theoretical method is also appropriate to obtain the linear fracture mechanics solution of the pure mode II and mixed mode I/II crack problems, and the free traction conditions are also satisfied.

The results in the paper can provide some useful information to predict the fracture or damage behaviors in the macrocrack tip for the materials containing microcracks. What is more, different damage configurations may play an important role for crack shielding or enhancement. The SIF of the crack is affected by many factors, such as the number of microcracks, the distance between macrocrack and microcrack, and the orientation and direction of the microcrack. Thus, the analysis of different damage configurations at a macrocrack tip is important to materials design.

Data Availability

All the data used in the manuscript are available upon request by contact with the corresponding author (E-mail: xiaoyujiang8@sina.com).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work is supported by the National Natural Science Foundation of China (11472230).

References

Y. Gao, C. Zhang, X. Xiong, Z. Zheng, and M. Zhu, “Intergranular corrosion susceptibility of a novel Super304H stainless steel,” Engineering Failure Analysis, vol. 24, pp. 26–32, 2012.

K. C. Hwang, S. Lee, and H. C. Lee, “Effects of alloying elements on microstructure and fracture properties of cast high speed steel rolls part II. Fracture behavior,” Materials Science & Engineering, vol. 254, no. 1-2, pp. 296–304, 1998.

J. A. Ronevich, B. P. Somerday, and C. W. San Marchi, “Effects of microstructure banding on hydrogen assisted fatigue crack growth in X65 pipeline steels,” International Journal of Fatigue, vol. 82, pp. 497–504, 2016.

J. Song, C. Fan, H. Ma, L. Liang, and Y. Wei, “Crack deflection occurs by constrained microcracking in nacre,” Acta Mechanica Sinica, vol. 34, no. 1, pp. 143–150, 2018.

S.-X. Gong and H. Horii, “General solution to the problem of microcracks near the tip of a main crack,” Journal of the Mechanics and Physics of Solids, vol. 37, no. 1, pp. 27–46, 1989.

S. X. Gong and S. A. Meguid, “A general solution to the antiplane problem of an arbitrarily located elliptical hole near the tip of a main crack,” International Journal of Solids and Structures, vol. 28, no. 2, pp. 249–263, 1991.

S. X. Gong and S. A. Meguid, “Microdefect interacting with a main crack. A general treatment,” International Journal of Mechanical Sciences, vol. 34, no. 12, pp. 933–945, 1992.

M. Hori and S. Nemat-Nasser, “Interacting microcracks near the tip in the process zone of a macrocrack,” Journal of the Mechanics and Physics of Solids, vol. 35, no. 5, pp. 601–629, 1987.

W. Tian and Y. Chen, “Further investigation of interaction between interface macrocrack and parallel microcracks in bimaterial anisotropic solids,” Acta Mechanica Sinica, vol. 15, no. 3, pp. 262–264, 1999.

M. Kachanov, “A simple technique of stress analysis in elastic solids with many cracks,” International Journal of Fracture, vol. 28, no. 1, pp. R11–R19, 1985.

M. Kachanov, “Elastic solids with many cracks: A simple method of analysis,” International Journal of Solids and Structures, vol. 23, no. 1, pp. 23–43, 1987.

M. Kachanov, “Effective elastic properties of cracked solids: critical review of some basic concepts,” Applied Mechanics Reviews, vol. 45, no. 8, pp. 304–335, 1992.

V. Tamuzs, N. Romalis, and V. Petrova, Fracture of Solids with Microdefects, 2000.

V. P. Tamuzs and V. E. Petrova, “On macrocrack-microdefect interaction,” International Applied Mechanics, vol. 38, no. 10, pp. 1157–1177, 2002.

V. Petrova and S. Schmauder, “Thermal fracture of a functionally graded/homogeneous bimaterial with system of cracks,” Theoretical and Applied Fracture Mechanics, vol. 55, no. 2, pp. 148–157, 2011.

V. Petrova and S. Schmauder, “Mathematical modelling and thermal stress intensity factors evaluation for an interface crack in the presence of a system of cracks in functionally graded/homogeneous bimaterials,” Computational Materials Science, vol. 52, no. 1, pp. 171–177, 2012.

V. Petrova, V. Tamuzs, and V. Petrova, “A survey of macro- microcrack interaction problems,” Applied Mechanics Reviews, vol. 53, no. 5, pp. 117–146, 2000.

R. G. Hoagland and J. D. Embury, “A Treatment of Inelastic Deformation Around a Crack Tip due to Microcracking,” Journal of the American Ceramic Society, vol. 63, no. 7-8, pp. 404–410, 2010.

I. Demir and H. M. Zbib, “A mesoscopic model for inelastic deformation and damage,” International Journal of Engineering Science, vol. 39, no. 14, pp. 1597–1615, 2001.

X. Yan, “Effective numerical approach for interaction of macrocrack with mircocracks,” Chinese Journal of Theoretical & Applied Mechanics, vol. 38, no. 1, pp. 119–124, 2006.

X. Q. Yan, “Microcrack interacting with a finite main crack,” Journal of Harbin Institute of Technology, vol. 40, no. 5, pp. 421–430, 2006.

X. Z. Xia, Q. Zhang, P. Z. Qiao, and L. I. Li-Juan, “Interaction between cracks and effect of microcrack zone on main crack tip,” Applied Mathematics and Mechanics, vol. 31, no. 1, pp. 67–76, 2010.

S. Rafoee, D. Gross, and T. Seelig, “The influence of microcrack nucleation on dynamic crack growth - a numerical study,” Engineering Fracture Mechanics, vol. 71, no. 4, pp. 849–857, 2004.

T. Seelig and D. Gross, “On the interaction and branching of fast running cracks - A numerical investigation,” Journal of the Mechanics and Physics of Solids, vol. 47, no. 4, pp. 935–952, 1999.

A. Brencich and A. Carpinteri, “Interaction of a main crack with ordered distributions of microcracks: a numerical technique by displacement discontinuity boundary elements,” International Journal of Fracture, vol. 76, no. 4, pp. 373–389, 1986.

X. Huang and B. L. Karihaloo, “Interaction of penny-shaped cracks with a half-plane crack,” International Journal of Solids and Structures, vol. 30, no. 15, pp. 2117–2139, 1993.

J.-P. Laures and M. Kachanov, “Three-dimensional interactions of a crack front with arrays of penny-shaped microcracks,” International Journal of Fracture, vol. 48, no. 4, pp. 255–279, 1991.

N. I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticity, Springer, Netherlands, 1975.

M. Kachanov, “Elastic solids with many cracks and related problems,” Advances in Applied Mechanics, vol. 30, no. 114, pp. 259–445, 1993.

Y. Murakami and L. M. Keer, “Stress intensity factors handbook,” Journal of Applied Mechanics, vol. 60, no. 4, article no. 1063, 1987.