Research Article  Open Access
Effect of the Preexisting Fissure with Different Fillings in PMMA on BlastInduced Crack Propagation
Abstract
In order to study the dynamic crack propagation law in fissured rock under the different fillings, a borehole with 7 mm diameter was processed in the center of a polymethyl methacrylate (PMMA) specimen. The preexisting fissure with different angles (θ = 0°, 45°, and 90°) and different distances (L = 20, 30, 40, 50, and 60 mm) was prefabricated around the borehole. Air, soil, and water were employed as fillings in the fissure, respectively. The experiment of explosive loading was carried out by a single detonator, and the dynamic crack propagation process of the experimental specimens was simulated by nonlinear dynamics software AUTODYN. The results show that the blastinduced cracks are the most favorable and unfavorable to propagate when θ = 0° and θ = 45°, respectively. The length of the farend wing crack decreases with the increase of the distance L, and the length of the farend wing crack in the airfilled specimens is larger than those in soilfilled and waterfilled specimens. The damagepressure curve of the farend wing crack initiation point shows “S”type change, and the damagepressure curve shows two obvious damage evolution processes of initial nonlinear and later linear stages. With the increase of the angle, the distance from the borehole to the crack initiation point decreases and the compressive stress wave peak value should increase, but the tensile force peak value decreases. Meanwhile, the relationships between pressure and average velocity of the initiation point and L, θ, and fillings are established, respectively. The numerical simulation agrees with the experimental results well. It can be seen that the fillings types, angle, and distance have a mutual restraint relationship with the reflected and absorbed stress wave energy. The phenomenon of crack propagation under different fillings can be explained well from the viewpoint of discontinuity degree and stress wave energy, which reveals the general law of blastinduced crack propagation.
1. Introduction
A large number of randomly distributed joints and fissures make the rock mass to show properties of the discontinuity, anisotropy, and inhomogeneity. For the fissured rock, when the stress waves travel to the fissure surface, the stress waves can produce the phenomena of reflection, refraction, diffraction, and transmission [1–3], which will influence the propagation direction and length of blastinduced cracks. Meanwhile, the loss of explosive gas in fissures causes the rock mass far away from borehole to decrease significantly in the explosive energy utilization rate, which will reduce the scope and effect of the blasting fragmentation, and not only affect the contour shape of the tunnel or mining slope and increase the explosive consumption greatly but also produce unfavorable effects on safety of construction, support, and postoperation, and so on as well as ultimately result in lowering production efficiency and increasing production cost.
Song and Kim [4] analyzed the mechanical properties of crack initiation, propagation, and bifurcation in the explosion process of rock fracture by numerical simulation. Zhu et al. [5–7] used AUTODYN to simulate the blastinduced crack propagation and concluded that the fissure hindered crack propagation seriously and fissure width and fillings also have an important influence on crack propagation. Ma et al. [8] also utilized LSDYNA to simulate the crack propagation in fissured rock mass and obtained the rule that the fissure greatly impeded crack propagation as well and the rock near the fissure or free side was more seriously broken caused by the reflected tensile stress waves [9]. Wang and Konietzky [10] analyzed the dynamic fracture process of the fissured rock mass in combination with LSDYNA and UDEC and concluded that different fissure types could cause different crack propagation laws. Sharafisafa et al. [11] studied the effect of fissure types on blastinduced crack propagation and concluded that the horizontal fissure was more favorable for crack propagation than inclined and vertical fissures. Wu and Wong [12] investigated initiation stress, initiation angle, and propagation law of wing crack caused by preexisting fissure end at different angles under uniaxial compression by using the numerical method. Li and Wong [13] revealed that the preexisting fissures of different angles had an important effect on wing crack initiation position, propagation angles, and propagation priority by AUTODYN. Zhuang et al. [14] researched the difference of wing crack propagation law between unfilled and filled preexisting fissures in rocklike brittle materials under uniaxial compression in the crack initiation stress and initiation angle. Yang et al. [15] studied the crack propagation when explosion stress waves with oblique incidence propagate to preexisting fissures and considered the propagation law of wing crack had an important relationship with initial compressive stress. Hu et al. [16] and Yang et al. [17, 18] both researched blastinduced crack propagation when stress wave traveled to preexisting fissure at different incident angles in PMMA and concluded the different angles would cause the different propagation mechanism of the farend wing crack.
Although many researchers have studied largely the blastinduced crack propagation law in fissured rock mass, they mostly studied the influence of single variable on crack propagation, such as stress wave incident angles [14, 16, 17] and fillings [5]. So they lack to consider synthetically the influence of filling material types and the angle and distance of fissure to borehole on the blastinduced crack propagation. Firstly, this paper mainly studies the influence of the preexisting fissure with the different filling materials (air, soil, and water) on the blastinduced crack propagation and the crack propagation law at the end of preexisting fissure by experimental studies, and then, the crack propagation law in the experiment will be further revealed by the corresponding numerical simulation through AUTODYN. The research of this paper is conducive to improving the energy utilization efficiency of rock blasting and the effect of rock fragmentation and the safety level of fissured rock slope and tunnel, which can provide theoretical basis for improving blasting design and optimizing blasting parameters.
2. Explosive Experiment
2.1. Experimental Scheme
PMMA and rock exhibit an anisotropic behavior and brittle characteristic under explosive loading [19, 20], and thus, PMMA is chosen as an experimental material in our study. With a length 400 mm × width 400 mm × height 5 mm of the PMMA plate as the experimental specimen, a 7 mm diameter hole is located in the center of the specimen, and the preexisting fissure size is length 60 mm × width 2 mm × height 5 mm, as shown in Figure 1. θ is the angle between the line from the borehole center to the proximal end of preexisting fissure and axis in the longitudinal direction of the preexisting fissure, and θ is called the angle between borehole and preexisting fissure for short, L is the distance from the borehole center to the proximal end of preexisting fissure. According to the fillings types, experiments are divided into three groups, namely, air group, soil group, and water group; each group is divided into three subgroups with θ as variable, namely, 0° subgroup, 45° subgroup, and 90° subgroup; each subgroup has five PMMA specimens with L as variable, and L is 20, 30, 40, 50, and 60 mm, respectively, as shown in Table 1.

2.2. Experimental Procedures
The experimental procedures are as follows:(i)The laser lathe was applied to incise the PMMA specimen and prefabricate the borehole and preexisting fissure, which requires walls of the borehole and preexisting fissure to be smooth and perpendicular to the specimen surface; Sellotape was pasted tightly on the backside of preexisting fissure in each soilfilled specimen and waterfilled specimen; gaskets were used to form a holder for underlaying the specimen via its four corners to the same height on the ground, which aimed at preventing the bottom of the detonator from touching the ground to change the main charge zone position.(ii)To the airfilled group, ① a #8 instantaneous electric detonator was fixed with the same batch, the similar resistance values in the center of the borehole and the detonator must be perpendicular to the specimen surface, and all the main charge zones of the detonator were required to align the thickness of the specimen in the same position; ② the specimen was put on the holder; ③ a thin wood plate that had a slightly bigger hole than the borehole with the same size of the specimen was covered above the specimen to prevent the detonator debris from scratching in the explosion; and ④ the detonator was initiated, and the specimen was recycled after explosion.(iii)To the soilfilled group, firstly soil was filled in preexisting fissure completely, which requires soil to be consistent with the thickness of the specimen so that it does not invade into the rear paste area of Sellotape, and then, the steps ①, ②, ③, and ④ were followed in the airfilled group.(iv)To the waterfilled group, firstly the steps ① and ② were followed in the airfilled group; then, a syringe was used to fill water in the preexisting fissure completely, and it was ensured that no bubbles exist in the fissure; and lastly, the steps ③ and ④ were followed in the airfilled group.
2.3. Experimental Results and Discussion
The crack propagation effects of the experimental specimens after explosion are shown in Figures 2–4.
2.3.1. The Influence of Fillings, Angle, and Distance on Crack Propagation
The mechanical properties of the fillings affect the discontinuity degree of the specimen, which was measured by the wave impedance. The wave impedances of four media are shown in Table 2 [21–24], where ρ is the density, c is the longitudinal wave velocity, η is the wave impedance, and Δη is the wave impedance difference value. According to the rock wave impedance matching viewpoint [25], the difference value of wave impedance between the PMMA and fillings is inversely proportional to the discontinuity degree of the specimen and is proportional to the magnitude of the blocked energy of stress wave, which means that the bigger the wave impedance difference value, the smaller the discontinuity degree and the more the blocked energy. Compared with wave impedances of PMMA, soil, and water, it can be considered that stress waves are completely blocked in the airfilled specimen because of its minimum wave impedance [26]; soilfilled and waterfilled specimens have a bigger discontinuity degree, which is conducive to stress wave travel. For the angle θ and the distance L, it is generally believed that the larger the angle and the smaller the distance of the three kinds of filling specimens, the larger the reflected stress wave energy. But the smaller the distance in the soilfilled and waterfilled specimens, the greater the absorbed stress wave energy. It can be seen that the fillings types, angle, and distance have a mutual restraint relationship with the reflected and absorbed stress wave energy. From the experimental results, it is because that the three variables lead to different crack propagation results, and the analysis should also consider the influence of three variables and the mutual restraint relationship of stress wave energy on blastinduced crack propagation so as to make a more comprehensive explanation.

2.3.2. The Influence of Fillings and Angle on the Total Number N of BlastInduced Cracks
Figure 5 shows the relationship between the total number N of blastinduced cracks exceeded 3 mm in length and fillings. Due to the minimum discontinuity degree of the airfilled specimen, the preexisting fissure completely blocks the explosive energy, and the reflected stress wave energy is larger, which is favorable for the formation of cracks. The soilfilled and waterfilled specimens have the bigger discontinuity degree, and a large amount of compressive stress wave energy is absorbed, so N in the airfilled specimens is substantially larger than that in the other two groups.
(a)
(b)
(c)
As shown in Figure 6, since absorbed stress wave energy at 45° is larger than that at 0°, and the reflected stress waves mostly travel far away from explosion source. When θ is 90°, the reflected stress waves mostly travel to explosion source, so N in the airfilled and the other specimens decreases initially and increases afterwards as the angle changes. Compared with the airfilled specimens, due to the bigger discontinuity degree of waterfilled specimens, the obtained lower energy of the reflected stress waves restrains the crack propagation, so N in the waterfilled specimens first increases and then decreases.
(a)
(b)
(c)
2.3.3. The Wing Crack Analysis
In the actual blasting project, the wing crack propagation in the fissured rock mass seriously affects the rock fragment range and the later engineering support, so it is necessary to analyze the farend wing crack. In the case of two farend wing cracks in the 0° group, for the convenience of studying, the numbers show the average values in Figures 7 and 8.
(a)
(b)
(c)
(a)
(b)
(c)
In Figure 7, the farend wing crack length in the whole specimens shows a decreasing trend, and the reason is that the tensionshear stress at the crack tip decreases with the increase of L, which causes the decrease in length. Because of the minimum discontinuity degree of the airfilled specimen, the farend wing crack length is larger than that in the soilfilled and waterfilled specimens; and the discontinuity degree of the soilfilled and waterfilled specimens is relatively bigger and closer, so both the farend wing cracks are shorter and more similar in length.
In Figure 8, when the angle θ increases, the loss of energy in the stress wave travel process decreases as the distance from the borehole center to the preexisting fissure farend decreases, so the farend wing crack in all the specimens increases as the angle increases when L is 20 to 30 mm. When L is 40 to 60 mm, the farend wing crack first decreases and then increases. The main reasons are as follows: the increase of L leads to the increase of energy loss in the stress wave travel process, compressive stress wave energy is further absorbed by fillings of soil and water, and relatively smaller energy loss in travel is caused because of the closer distance from the borehole center to the preexisting fissure farend in the 90° group than that in the 45° group.
3. The Numerical Simulation Analysis of Crack Propagation under the Different Fillings
In order to further reveal the blastinduced crack propagation law in the fissured rock mass under the explosive loads, the nonlinear dynamics software AUTODYN is used for numerical simulation according to the corresponding parameters of the experimental specimen, involving materials including PMMA, fillings (air, soil, and water), and explosive.
3.1. Material Model
3.1.1. PMMA
To better cooperate with the tensile crack softening failure model, the linear equation of state (EOS) is applied in PMMA [27]:where is the explosive pressure, is the volume modulus of the material, and is the ratio of the current density to the initial density during the explosion.
The tensile crack softening model is applied as the PMMA failure model, and the detailed description of this model is in [16, 27, 28]. Generally, the relevant physical and mechanical parameters of PMMA under the dynamic loads are difficult to measure, which is related to the material properties, experimental conditions, temperature, and so on. So the author revised some parameters according to the parameters listed in [13], as shown in Tables 3 and 4 [29, 30].


3.1.2. Fillings
(a)AirFor air, the EOS is described by the ideal gas:where is the pressure; is the adiabatic exponent, 1.4; is the density, 1.25 kg/m^{3}; and is the internal energy.(b)SoilThe soil is described by the shock EOS and the Drucker–Prager strength model. The concrete parameters are shown in Table 5 [31].(c)WaterA polynomial EOS is applied to water:For (compression),For (tension),where in the above two formulas is the pressure, is the density, is the internal energy, and are the constants, and the specific parameters are shown in Table 6.


3.1.3. Explosive
The explosive is ammonium nitratefuel oil (ANFO) explosive, and the EOS is described by JWL:where is the pressure, is the initial internal energy of the explosive product, is the relative volume, and are the constants. The material parameters of ANFO are listed in Table 7.

3.2. Material Location and Grid Mesh
According to the experimental specimen, the 2D model and noncoupling charge structure are applied. The diameters of the borehole and explosive are 7 mm and 6 mm, respectively. A variety of grid sizes were applied to simulate the crack propagation, and we found that too big or too small grid sizes would lead to unsatisfactory or nonideal crack results compared with the experimental crack results. Meanwhile, we also did the convergence computation through all kinds of grid sizes. Based on the above previous works before simulation, we chose 1.5 mm as the grid size at last and the total number elements were 78655. The mesh method is sweep, the boundary condition is the free boundary, and the calculation time is 100 μs. The location of the materials and the mesh are shown in Figure 9.
(a)
(b)
3.3. Simulation Results
The numerical simulation results are shown in Figures 10–12. Due to the high explosive pressure on the borehole wall, a crushed zone develops near the borehole and a fractured zone also develops because of the radial cracks generation after the formation of the crushed zone. Blastinduced cracks almost cannot get over the preexisting fissure to continue to propagate, and the wing cracks appear at both ends of the preexisting fissure. It should be noted that the farend wing cracks propagate horizontally at θ = 0°, while the farend wing cracks propagate above obliquely at θ = 45° and θ = 90°. Additionally, the crack is produced on the other side of the preexisting fissure in the airfilled specimen of L = 30 mm at θ = 45°. All these phenomena appearing in numerical simulation agree basically with the experimental phenomena.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
3.4. The BlastingInduced Crack Propagation Mechanism
Because the explosion of the detonator in experiment has less gas and the calculation time is relatively shorter in numerical simulation, it is considered that the stress wave is the most important factor in blastinduced crack initiation and propagation. From the experimental and simulative results, it can be concluded that the crack propagation process can be divided into two stages through the time points that whether the reflected stress wave formed in the preexisting fissure has common interaction with the compressive stress wave or not. The first is the phase of the compressive stress wave, and the second is the phase of the common interaction between compressive stress wave and reflected stress wave.
3.4.1. The Radial Crack Propagation Mechanism around the Borehole
After detonator or explosive detonation, the shock wave compresses on the borehole wall intensely to develop the crushed zone around the borehole because of the shear stress. The shock wave attenuates the stress wave during the process of the crushed zone formation. The stress wave compresses radially the media outside the crushed zone to cause the tangential stress, and once the tangential stress reaches the dynamic tensile strength of the medium, it will lead to radial crack initiation and propagation. Then, the reflected stress wave and the compressive stress wave act together to promote the crack propagation again.
Considering the different positions between preexisting fissure and borehole, the reflected stress wave that travels back to the borehole is different, and radial crack propagation mechanism is also different, so it is necessary to discuss the influence of the distance, angle, and fillings on crack propagation.(1)Distance L
In numerical simulation, take the airfilled specimens at θ = 90° as an example and choose the element as the gauss point at 10 mm below the preexisting fissure center from L = 20 to 60 mm; the relationship between the peak pressure and L is shown in Figure 13. It shows the law that the peak pressure decreases when L increases.(2)Angle θ
When θ = 0°, the compressive stress wave travels parallely to the preexisting fissure, and the reflected stress wave at the proximal end almost does not act on crack propagation, so the crack is mostly caused by the compressive stress wave; when θ = 45°, the compressive stress wave travels obliquely to the preexisting fissure, and the reflected stress wave travels far from the media around the borehole, in which another small part of the reflected stress wave has limited effect on the propagation and initiation of crack; when θ = 90°, the compressive stress wave travels vertically to the preexisting fissure, and the more reflected stress wave can help to crack propagation. Therefore, the increase of the angle θ is beneficial to the crack propagation.(3)Fillings
When θ = 90° and L = 50 mm, the element is also selected as the gauss point at 10 mm below the preexisting fissure center of the three fillings specimens. As L is the same, theoretically compressive stress wave peak pressure should also be basically the same. To reflect the influence of fillings on the reflected stress wave, take the peak pressure of stress wave pressure as a reference, as shown in Figure 14. The peak value of the reflected stress wave in the airfilled specimen is maximum, which shows that the airfilled specimens are more conducive to the crack propagation than the soilfilled and waterfilled specimens; comparing the peak value of the reflected stress wave with Δη (the difference value of wave impedance between the three kinds of fillings and PMMA), it can be seen that the peak pressure value increases with the increase of Δη, which further indicates that it is valid to use the difference value of wave impedance to measure the discontinuity degree of the specimen and adopt the discontinuity degree to explain the phenomena of crack propagation.
3.4.2. The Wing Crack Propagation Mechanism
In numerical simulation, the damagetime curve is obtained by setting the initiation point of the farend wing crack as the gauss point. Taking the airfilled specimen at θ = 45° and L = 20 mm as the example, the damage value of this specimen is fitted by Slogistic function, as shown in Figure 15, and the fitting function iswhere is the time (s) and is the damage. By fitting the curve, it can be seen that the damage value is the “S” curve changing with time. According to the characteristics of change in the fitting curve, damage evolution is divided into three stages: OA stage is the early stage of damage, and the damage relatively increases more slowly with time; AB stage is the middle stage of damage, and the damage value increases rapidly in the short time; BC stage is the late stage of damage because the damage value of B point is close to 1, so in a very short time the damage value achieves 1, but this stage curve changes slowly, and it can be speculated that the damage growth rate decreases. Overall, the damagetime experiences the evolution process that OA stage keeps an upward trend but in a less significant way, AB stage shows an obvious soar trend, and then BC stage rises moderately.
The dynamic tensile damage evolution process is reflected fully in Figure 16; namely, there are two growth stages of significant difference in the tensile damage evolution: one is the AB phase of nonlinear fast growth, the other is the BC phase of linear growth, and the transition of the two stages has an obvious critical characteristic. When the initiation point of the farend wing crack is subjected to tensile force, damage occurs, so initial damage point A occurs at the tension stage, but the tensile force at point A is not the peak value. Once the damage occurs, the damage will experience a rapid nonlinear growth to the turning point B, whose corresponding tensile force is the peak value; the damage increases linearly when the tensile force decreases, and finally, the damage value reaches 1 when the medium becomes totally fractured.
Figure 17 shows the typical pressuretime history curves at different angles of the initiation point of the farend wing crack in the airfilled specimen at L = 20 mm. The compressive stress wave acts on the initiation point before the reflected stress wave and the diffracted stress wave in the preexisting fissure, so firstly the initiation point is subjected to the pressure and then the tensile force. With the increase of the angle, the distance from the borehole to the crack initiation point decreases and the compressive stress wave peak value increases, but the tensile force peak value decreases. The reason is that the reflected stress wave has a serious impact on the subsequent compressive stress wave, which results in reduction of the tensile force peak when the phenomenon of two kinds of stress wave superimposition appears. Additionally, it is also found that peak pressure decreases faster than the tensile force peak at the farend initiation point, and the whole pressuretime decreases when the angle increases.
The effects of distance, angle, and fillings on the initial damage and fracture time of the initiation point of the farend wing crack are discussed below. Generally, the distance that the stress wave travels to the far end of the preexisting fissure also increases as the distance L increases, which causes the delay of initial damage and fracture time of the farend wing crack initiation point. If the angle θ increases, the initial damage and fracture time of the initiation point is earlier, and the distance that the stress wave travels to the far end of the preexisting fissure decreases as the angle θ increases. When the fillings change, it can be seen that the initial damage and fracture time of the farend wing crack initiation point in the airfilled specimen is the earliest and the soilfilled specimen is the latest, and the reason is that the smaller the Δη value is, the larger the discontinuity degree is and the more the energy absorption of the compressive stress wave is, so the initial damage and fracture time is different.
The compressive stress wave travels to the preexisting fissure, and the pressures are firstly generated at the initiation points of the both ends. The pressure is gradually converted to the tensile force as the reflected stress wave generated in the prefabricated crack continues to pressure on initiation points. As time goes on, the initiation point begins to damage, and the damage value increases rapidly, and finally the damage value reaches 1. During the whole process, the stress concentration occurs at the crack initiation point. When the stress strength exceeds the dynamic tensile strength of PMMA, the wing cracks begin to initiate and propagate. In the crack propagation process, the shear stress is produced by the shearslip caused by the uneven stress at the crack tip, which indicates that the wing crack propagates in a tensileshear mode. When the stress strength is lower than the dynamic initiation strength, the wing crack arrests.
Through the analysis of the physical experiment and numerical simulation, for most of the specimens, the wing crack propagation direction exists in two situations: when the distance L is relatively smaller (L = 20–30 mm), the blastinduced crack outside the crushed zone directly propagates through the proximal end of preexisting fissure; when the distance L is relatively larger (L = 40–60 mm), the proximalend wing crack propagates toward the borehole. However, farend wing cracks propagate away from the borehole in all the specimens. Meanwhile, it is found that the wing cracks propagation paths are almost parallel to the connection line between borehole center and two ends of the preexisting fissure, which indicates that the wing crack propagation path is directly related to the incident direction of the stress wave.
3.5. The Influence of Fillings and Angles on Total Number N of Cracks
In Figure 18, the total number N of cracks in the airfilled specimens is almost larger than that in the solifilled and waterfilled specimens due to the minimal discontinuity degree of the airfilled specimen. The total number N of cracks in the solifilled and waterfilled specimens is almost close due to the close discontinuity degree of the two kinds of specimens. This indicates the numerical simulation results are in accordance with the experimental results.
(a)
(b)
(c)
In Figure 19, when L = 20 mm and 40 mm, the total number N increases first and then decreases when the angle increases; when L = 30 mm, the total number N decreases first and then increases; when L = 50 mm, the total number N decreases when the angle increases; when L = 60 mm, the total number N decreases first and then increases when the angle increases, which indicates that blastinduced cracks show a complex phenomenon of mutual restraint under the variables of the fillings, angle θ, and distance L.
(a)
(b)
(c)
3.6. The FarEnd Wing Crack Propagation Law
3.6.1. Analysis on Length of FarEnd Wing Crack
In Figure 20, the length of the farend wing crack gradually decreases with the increase of the distance L. Because of the minimum discontinuity degree in the airfilled specimens, the length of the farend wing crack is larger than that of the soilfilled and waterfilled specimens. Therefore, the law of the farend wing crack propagation in numerical simulation and physical experiment is more consistent.
(a)
(b)
(c)
It is well known that the reason for the decrease in the farend wing crack length when the distance L increases is due to the decrease in the explosion pressure at the crack tip. To specify this reason, during the simulation, the peak pressure value at the different distance L is obtained by setting the initiation point of the farend wing crack as the gauss point, as shown in Figure 21. It can be seen that the peak pressure value of the initiation point does decrease with the increase of the distance L, and the peak pressure value of the initiation point of the farend wing crack in the airfilled specimens is larger than that of the soilfilled and waterfilled specimens, which explains the phenomena that the farend wing crack length increases when the distance L increases; the wing crack length in the airfilled specimens is much larger than that in the soilfilled and waterfilled specimens.
(a)
(b)
(c)
In Figure 22, the stress wave energy that travels to the farend of the preexisting fissure increases when angle θ increases, and the more energy is obtained when the farend wing crack initiates, which contributes to the wing crack propagation. Therefore, the length of the farend wing crack increases as a whole when the angle θ increases, which is consistent with the experimental results.
(a)
(b)
(c)
3.6.2. Analysis on Average Velocity of the FarEnd Wing Crack Propagation
Figure 23 shows the relationship between the average velocity of farend wing crack propagation () and the fillings. The trend of in the airfilled specimen is different at different angles, which is related to the stress of the crack tip caused by fillings and distance L. As a whole, the discontinuity degree of soilfilled and waterfilled specimens is more similar, and the pressures of the initiation point of the farend wing crack are close, so the trends of are consistent in the two kinds of specimens.
(a)
(b)
(c)
Figure 24 shows the variation trend of average crack propagation velocity () at different angles. In general, with the increase of the angle, of different fillings specimens shows the downtrend; when θ = 0°, is larger than that at θ = 45° and 90°, which indicates that θ = 0° is conducive to crack propagation and the validity of the literature [8]. When L = 20–40 mm, of the different fillings specimens decreases when the angle θ increases; when L = 50 mm, first increases and then decreases when the angle increases, and is inconsistent when L = 60 mm.
(a)
(b)
(c)
4. Conclusions
(1)Compared with the experimental results, the numerical simulation can reflect better the entire variation law of blastinduced crack propagation that the total number N of cracks and the length of the farend wing crack in the airfilled specimen are greater than that in the soilfilled and waterfilled specimens.(2)The wing crack is a kind of mixed mode crack by propagating in the tensileshear mode. The similar propagation path of the farend wing crack is related to the direction of the compressive stress wave propagation at the same angle. The length of the farend wing crack in airfilled specimens is larger than that in the soilfilled and waterfilled specimens, which shows that the damage range caused by explosion is the largest under this condition.(3)The damagetime curve of the initiation point of the farend wing crack presents “S”type change, and the damagepressure curve of the initiation point can fully reflect the law of the dynamic tensile damage evolution process.(4)The specimen discontinuity degree is measured by the wave impedance difference value between the three fillings and PMMA, and the greater the difference value, the smaller the discontinuity degree and the greater the reflected wave energy. It is generally believed that the larger the angle and the smaller the distance, the greater the reflected wave energy should be in the three kinds of filling specimens, but meanwhile the greater the transmitted stress wave energy also is. So there is a mutual restraint relationship between the distance L, angle θ, and fillings types and reflected and transmitted stress wave energy that affects the effect of blastinduced crack propagation. The phenomenon of crack propagation under different fillings can be explained well from the viewpoint of discontinuity degree and stress wave energy, which reveals the general law of blastinduced crack propagation.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the support of NSAF (no. U1530140) and the Beijing Institute of Applied Physics and Computational Mathematics through the contract of HG7017105.
References
 Y. Li, Z. Zhu, B. Li, J. Deng, and H. Xie, “Study on the transmission and reflection of stress waves across joints,” International Journal of Rock Mechanics and Mining Sciences, vol. 48, no. 3, pp. 364–371, 2011. View at: Publisher Site  Google Scholar
 L. F. Fan and L. N. Y. Wong, “Stress wave transmission across a filled joint with different loading/unloading behavior,” International Journal of Rock Mechanics and Mining Sciences, vol. 60, pp. 227–234, 2013. View at: Publisher Site  Google Scholar
 X. Huang, S. Qi, A. Williams, Y. Zou, and B. Zheng, “Numerical simulation of stress wave propagating through filled joints by particle model,” International Journal of Solids and Structures, vol. 6970, pp. 23–33, 2015. View at: Publisher Site  Google Scholar
 J. Song and K. Kim, “Micromechanical modeling of the dynamic fracture process during rock blasting,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 33, no. 4, pp. 387–394, 1996. View at: Publisher Site  Google Scholar
 Z. Zhu, B. Mohanty, and H. Xie, “Numerical investigation of blastinginduced crack initiation and propagation in rocks,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 3, pp. 412–424, 2007. View at: Publisher Site  Google Scholar
 Z. Zhu, H. Xie, and B. Mohanty, “Numerical investigation of blastinginduced damage in cylindrical rocks,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 2, pp. 111–121, 2008. View at: Publisher Site  Google Scholar
 Z. Zhu, “Numerical prediction of crater blasting and bench blasting,” International Journal of Rock Mechanics and Mining Sciences, vol. 46, no. 6, pp. 1088–1096, 2009. View at: Publisher Site  Google Scholar
 G. W. Ma and X. M. An, “Numerical simulation of blastinginduced rock fractures,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 6, pp. 966–975, 2008. View at: Publisher Site  Google Scholar
 M. Bendezu, C. Romanel, and D. Roehl, “Finite element analysis of blastinduced fracture propagation in hard rocks,” Computers and Structures, vol. 182, pp. 1–13, 2017. View at: Publisher Site  Google Scholar
 Z. L. Wang and H. Konietzky, “Modelling of blastinduced fractures in jointed rock masses,” Engineering Fracture Mechanics, vol. 76, no. 12, pp. 1945–1955, 2009. View at: Publisher Site  Google Scholar
 M. Sharafisafa, Z. Aliabadian, R. Alizadeh, and M. Ali, “Distinct element modelling of fracture plan control in continuum and jointed rock mass in presplitting method of surface mining,” International Journal of Mining Science and Technology, vol. 24, no. 6, pp. 871–881, 2014. View at: Publisher Site  Google Scholar
 Z. Wu and L. N. Y. Wong, “Frictional crack initiation and propagation analysis using the numerical manifold method,” Computers and Geotechnics, vol. 39, pp. 38–53, 2012. View at: Publisher Site  Google Scholar
 H. Li and L. N. Y. Wong, “Influence of flaw inclination angle and loading condition on crack initiation and propagation,” International Journal of Solids and Structures, vol. 49, no. 18, pp. 2482–2499, 2012. View at: Publisher Site  Google Scholar
 X. Zhuang, J. Chun, and H. Zhu, “A comparative study on unfilled and filled crack propagation for rocklike brittle material,” Theoretical and Applied Fracture Mechanics, vol. 72, pp. 110–120, 2014. View at: Publisher Site  Google Scholar
 L. Yang, R. Yang, G. Qu, and Y. Zhang, “Caustic study on blastinduced wing crack behaviours in dynamicstatic superimposed stress field,” International Journal of Mining Science and Technology, vol. 24, no. 4, pp. 417–423, 2014. View at: Publisher Site  Google Scholar
 R. Hu, Z. Zhu, J. Xie, and D. Xiao, “Numerical study on crack propagation by using softening model under blasting,” Advances in Materials Science and Engineering, vol. 2015, Article ID 108580, 9 pages, 2015. View at: Publisher Site  Google Scholar
 R. Yang, Y. Wang, and C. Ding, “Laboratory study of wave propagation due to explosion in a jointed medium,” International Journal of Rock Mechanics and Mining Sciences, vol. 81, pp. 70–78, 2016. View at: Publisher Site  Google Scholar
 R. Yang, P. Xu, Z. Yue, and C. Chen, “Dynamic fracture analysis of crack–defect interaction for mode I running crack using digital dynamic caustics method,” Engineering Fracture Mechanics, vol. 161, pp. 63–75, 2016. View at: Publisher Site  Google Scholar
 H. P. Rossmanith, A. Daehnke, R. E. K. Nasmillner, N. Kouzniak, M. Ohtsu, and K. Uenishi, “Fracture mechanics applications to drilling and blasting,” Fatigue and Fracture of Engineering Materials and Structures, vol. 20, no. 11, pp. 1617–1636, 1997. View at: Publisher Site  Google Scholar
 D. Rittel and A. Brill, “Dynamic flow and failure of confined polymethylmethacrylate,” Journal of the Mechanics and Physics of Solids, vol. 56, no. 4, pp. 1401–1416, 2008. View at: Publisher Site  Google Scholar
 X. Song, Study on the Stratification Cracking Effect of the Layered Rock Slope Caused by Blasting and Its Influence on the Stability of Layered Slope, Southwest Jiaotong University, Chengdu, China, 2007.
 R. Yang, Z. Yue, T.S. Xiao et al., “Dynamic caustics experiment on crack propagation of jointed medium with controlled blasting,” Chinese Journal of Rock Mechanics and Engineering, vol. 27, no. 2, pp. 244–250, 2008. View at: Google Scholar
 P. Zhang, Y. Li, X. Wu et al., “Locating system of pipeline inspection gauge based on acoustic detecting,” Modern Scientific Instruments, vol. 1, pp. 45–47, 2011. View at: Google Scholar
 W. Wang, Experimental Study on the Acoustics Characteristic of Layer Rock and Rock Damage under Explosion Load, China University of Geosciences (Beijing), Beijing, China, 2009.
 X. Guo and J. Zhang, Blasting Engineering, China Communications Press, Beijing, China, 2007.
 Y. Guo, H. Meng, Y. Qi et al., “Study on mechanisms of the transmission of blasting vibration in the presplit crack,” Chinese Journal of Scientific Instrument, vol. 31, no. 4, pp. 17–20, 2010. View at: Google Scholar
 Century Dynamics Inc., AUTODYN Theory Manual, Century Dynamics Inc., Concord, CA, USA, 2003.
 Z. Tu and Y. Lu, “Modifications of RHT material model for improved numerical simulation of dynamic response of concrete,” International Journal of Impact Engineering, vol. 37, no. 10, pp. 1072–1082, 2010. View at: Publisher Site  Google Scholar
 N. Murphy, M. Ali, and A. Ivankovic, “Dynamic crack bifurcation in PMMA,” Engineering Fracture Mechanics, vol. 73, no. 16, pp. 2569–2587, 2006. View at: Publisher Site  Google Scholar
 Z. W. Yue, L. Y. Yang, and Y. B. Wang, “Experimental study of crack propagation in polymethyl methacrylate material with double holes under the directional controlled blasting,” Fatigue and Fracture of Engineering Materials and Structures, vol. 36, no. 8, pp. 827–833, 2013. View at: Publisher Site  Google Scholar
 B. Luccioni, D. Ambrosini, G. Nurick, and I. Snyman, “Craters produced by underground explosions,” Computers and Structures, vol. 87, no. 2122, pp. 1366–1373, 20s09. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Xin Yang 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.