Next Article in Journal
Sensitivity of Axial Velocity at the Air Gap Entrance to Flow Rate Distribution at Stator Radial Ventilation Ducts of Air-Cooled Turbo-Generator with Single-Channel Ventilation
Previous Article in Journal
Assessing the Impact of DFIG Synthetic Inertia Provision on Power System Small-Signal Stability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fracture Evaluation and Dynamic Stress Concentration of Granite Specimens Containing Elliptic Cavity under Dynamic Loading

1
School of Resources and Safety Engineering, Central South University, Changsha 410083, China
2
State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, China
*
Authors to whom correspondence should be addressed.
Energies 2019, 12(18), 3441; https://doi.org/10.3390/en12183441
Submission received: 24 June 2019 / Revised: 18 August 2019 / Accepted: 23 August 2019 / Published: 6 September 2019

Abstract

:
The Split-Hopkinson pressure bar (SHPB) was used to determine the fracture characteristics of a long bar rock specimen with an elliptical cavity under different axial ratios and dip angles. A high speed camera was applied to record the fracturing process of the granite specimen around the cavity. The experimental results showed that the fracture characteristics around the elliptical cavity were closely related to the axial ratio and dip angle. A three-dimensional numerical model was established using LS-DYNA to quantitatively analyze the dynamic stress state around the cavity. The numerical results indicate that the dip angle and axial ratio of the elliptical cavity significantly affected the dynamic stress concentration factor (DSCF), then affected the rock failure. The location of higher DSCF led to a higher possibility of spalling failure. The maximum DSCF remarkably decreased with a decreasing dip angle and increased the axial ratio. In the dynamic loading propagate process, the stress concentration distribution around the cavity formed by a compression stress wave had a certain damaging effect on the destruction of rock around the cavity, and the stress concentration generated by the tensile stress wave was the main factor of the rock fracture, which was most notable in the peak area of the stress concentration.

1. Introduction

The construction of cavities widely exists in civil and mining engineering. Under external disturbances such as blasting and earthquakes, the stress waves that spread to cavities will produce a series of reflection and transmission phenomena. The stress wave arrives at different diffraction degrees around the cavity and strengthens the local stress, which influences the stability of the engineering structure.
The Hopkinson bar is commonly used for testing rock dynamic properties [1]. Based on the Hopkinson bar, many scholars have studied the failure characteristics of rock and rock-like materials [2,3,4]. Regarding rock specimens with a cavity, Li et al. [5] conducted an experimental study on the dynamic failure characteristics of a rock specimen containing an elliptic hole based on the Split-Hopkinson pressure bar (SHPB). In their experiment, the dynamic compressive strength, failure mode, and crack propagation characteristics of a rock specimen with a hole under different dynamic loads were analyzed. Tao et al. [6,7] carried out a static and dynamic combination test of the prefabricated circular hole specimen by using an SHPB device and compared the results with theoretical calculations. They concluded that the dynamic stress concentration area was the first region damaged in the specimen.
Theoretically, the stress distribution around holes under dynamic incident loading is commonly studied as a diffraction or scattering of elastic waves. Pao et al. [8] systematically studied the basic principle of wave diffraction and the main directions of the related theory to determine the diffraction in various shapes of holes under steady-state and transient wave incidence conditions. Liu et al. [9] analyzed the scattering of SH waves by a rectangular cavity. Yi et al. [10,11] obtained the dynamic response of a circular lining tunnel as function of a cylindrical P-wave. Tao et al. investigated the dynamic stress concentration factor (DSCF) around a circular hole subjected to transient waves by using the complex variable method [12].
Additionally, some scholars have conducted thorough research on the cavity diffraction question from aspects of numerical simulation. Lin et al. [13] proposed a simple and easy numerical solution technique for the wave diffraction problem. She and Guo [14] studied the stress concentration factors under tensile stress in the elliptical hole wall of isotropic materials by the finite element method. Zhou et al. [15] verified that the dynamic stress state around the elliptical hole had a relationship with the wave number. Deng et al. [16] found that the joint position of rock mass had a great influence on the tunnel damage by using UDEC. Du et al. [17] analyzed specimen shape and cross-section effects on rocks properties. More recently, Fang et al. [18] analyzed the influence of a viscoelastic interface on the dynamic stress around the tunnel by means of numerical examples. Kung et al. [19] concluded that the influence of tunnel radius on the earthquake stress increment was stronger than that of the tunnel depth by using ABAQUS. Li et al. [20], combined with PFC-2D, analyzed the dynamic stress concentration and the energy evolution process surrounding the hole affected by the explosion stress wave. Additionally, Wang et al. [21,22] reported that the hole defects of brittle materials had a great influence on dynamic fracturing, and mechanized excavation face.
Theoretical research, experiments, and numerical simulations of the holes above showed that there are more studies on circular holes and that the various theoretical methods are mature. However, studies on the effect of the dynamic stress state formed by the diffraction of the stress wave through elliptical holes on the defective rock are rare. In practical engineering, defects exist in small cracks and large chambers. The problem of defective rocks can be described more comprehensively in terms of the shape of elliptical holes under different axial ratios and dip angles. In contrast to circular holes, the change of the axial ratio and dip angle of an elliptical hole will greatly affect the stress concentration distribution around the hole when the stress wave passes through the hole. Therefore, it is of great practical importance to study the stress concentration and failure characteristics of an elliptical-hole rock based on different axial ratios and dip angles under the influence of a stress wave.
In the present study, based on the SHPB experimental platform, spallation failure tests for long bar granite specimens with different elliptical holes were conducted. The fracture location and process were recorded by a high speed camera. Furthermore, numerical simulations were employed to discuss the dynamic stress distribution of the compressed and tensile stress waves propagating to the elliptical hole under the different dip angles and axial ratios.

2. Experimental Process and Result Analysis

2.1. Material Identification and Experimental Program

Granite samples were chosen of a relatively homogeneous granite from a quarry mine in Dingzi Town, Changsha, Hunan Province, China. Granite, which is plutonic acidic igneous rock, is generated from the condensation of magma. The dominant minerals identified were feldspar, quartz, mica, and a small amount of dark minerals. The complexity and diversity of the mineral component make the rock fracture uniquely inhomogeneous and the study of the rock fracture more difficult. To provide guidance for the microscopic analysis of the rock fracture, the mineralogical characteristics of the granite samples were studied before the spall tests.
After identification, the rock sample was named as gray-white biotite granite, as shown in Figure 1a and the mineral content is shown in Figure 1b. Among them, quartz accounted for about 50% of the total composition, which occurred as an irregular granular form with flaws.
The modified SHPB experiment apparatus was adopted in this study, and the details have been described in our previous studies [23,24,25]. The wavelength of the incident wave was about 72 mm, and wave period was about 240 μs. The rock specimen was cuboid with a cross-section of 35 mm × 35 mm and a length of 1200 mm; the density was 2694 kg/m3, the Poisson’s ratio was 0.19, the elasticity modulus was 45.3 GPa, P-wave velocity was 4104 m/s, uni-axial tensile strength was 9.4 MPa, and the uni-axial compression strength was 128 MPa. Six surfaces of the cuboid were carefully polished and non-parallelism of its surface were within 0.02 mm. An elliptical hole with a different axial ratio of the long axis to the short axis and dip angles was then drilled in the center of the lateral surface. Keeping the long axis of the elliptical hole as 10 mm, and the short axis as 2 mm, 4 mm, 7 mm, and 10 mm, respectively, the corresponding axial ratios were m = 5, m = 2.5, m = 1.43, and m = 1. The dip angle of holes is the angle between the long axis direction and the loading direction, and there are five kinds of dip angles of a specimen with an elliptical hole under each axial ratio, which are β = 0°, β = 30°, β = 45°, β = 60°, and β = 90°. The specimen was in contact with the input bar, and the schematic views of the modified SHPB device are shown in Figure 2.
Due to the existence of holes, the incident compressive stress wave will produce a diffraction phenomenon, leading to a dynamic stress concentration, thus causing the non-uniform stress distribution surrounding the hole. When the incident compressive stress wave passes through the hole and continues to spread toward the right free surface of the specimen, it is reflected as the tensile stress wave. When the reflected tensile stress wave spreads from right to left, it will also generate diffraction after encountering the hole, which will lead to dynamic stress concentration generated by a tensile wave in the periphery of the hole. The wavelength of the incident wave is considered in the design process, which avoids encountering incident and reflective stress waves near the holes. Given the length of the specimen is two times that of the wavelength, the reflective stress wave has no effect on the incident wave when a whole half-sine incident stress wave passes through a hole.

2.2. Experimental Results and Analyses

In the experimental process, the constant gas pressure of 0.5 MPa was applied for the impact, and the whole rock damage process around the hole caused by the dynamic loading was observed using a high-speed camera employed to record the failure process. The resolution ratio of the high-speed camera adopted in this test was 128 × 128, and the videography frequency was 180,000 fps. Herein, the arrival time of the incident stress wave at the left end of the specimen was defined as the zero point of the time record, based on the geometry and wave velocity of the specimen. The complete failure processes are presented in Figure 3, where the left end is the impact incidence.
In addition, Figure 4 shows the rock fracture around the hole that changed with time at different axial ratios when β = 90°. When t = 139 μs, the incident compressive stress wave almost arrived at the left side of the hole. The incident wave completely passed through the hole at about t = 378 μs. When t = 533 μs and t = 589 μs, the reflective tensile wave spread through the hole. Up until the t = 828 μs moment, the reflective stress wave completely passed through the hole.
The high-speed camera images showed that the dip angle β = 90°. When the incident compressive stress wave passed through the hole, there was some rock debris ejection at the upper and lower positions of the holes, whereas such phenomena did not appear at other positions. This shows that powder-like particles pop up when the axial ratio is small. As the axial ratio increases, the number of ejected debris and debris size increase, but at this time only as small debris without cracks. When the reflective tensile stress wave gradually returns and re-incidents through the hole, a crack gradually forms at the upper and lower positions. Moreover, the crack slowly stretches along the upper and lower positions. Figure 4 shows that there were some differences in the time of the crack initiation and elongation around the hole under different axial ratios. The crack initiation at an axial ratio of 2.5 occurred earlier than that related to the other two axial ratios. At t = 533 μs, signs of a crack appeared in the upper part of the hole. At t = 589 μs, the crack notably formed and the width of the crack was greater than that of the other two axial ratios. Similarly, in comparison with the cracks of m = 1.43 and m = 1, signs of cracking appeared in the upper and lower parts of the hole at t = 589 μs and m = 1.43, which were especially notable in the upper part. At the same time, only a tiny crack appeared in the upper part of the hole related to m = 1; the signs of cracks in the lower part were not stronger. At t = 828 μs, the crack around the hole had completely formed for each axial ratio. In contrast, the width of the crack was significantly different, that is, the width of the crack with a high axial ratio was greater than that of the crack with a low axial ratio.
The results show that the specimen with the defective elliptic hole had different degrees of fracture around the hole under compressive and tensile stress waves, but that both failure positions were in the defective hole, and that the shape of the fracture changed with the variation in the axial ratio and dip angle of the hole. Figure 5 shows a comparison of the failure modes of partially prefabricated holes under different axial ratios and dip angles. When the dip angle β = 90°, the failure point around the holes under different axial ratios was at the top and bottom of the hole, and the damage position was symmetrical with the short axis, that is, basically on a straight line with the regular destruction form. As the dip angle decreased to β = 45° and β = 30°, the failure point appeared in the upper right and lower left corners of the hole. At this time, the damage position was no longer symmetrical with the short axis, but was close to the central hole, whereas the failure surface was still regular, that is, perpendicular to the spread direction of the stress wave.
When the dip angle β = 0°, the fracture position and destruction form were similar to that of β = 90°. Based on the damage form of all holes in the picture, we concluded that apart from the fracture around the hole, the other locations ere intact and only some of the holes had uneven pits. Compared with the failure form of the hole, the dip angle of the hole had a certain influence on the failure of the rock. Taking the center of the hole as the reference point, the angle of the failure point will fluctuate with the change of the dip angle.
Figure 4 and Figure 5 show that the ejection of rock debris caused by the incident compressive stress wave and the crack initiation and propagation caused by the reflective tensile stress wave all occurred in the upper and lower positions of the hole. Similar phenomena were not observed in other places. The rock debris ejected by the incident compressive stress wave through the hole led to the initial damage of the defective hole in a certain area, but the initial damage was not enough to break the defective specimen. When the tensile stress wave passed through the hole, it formed different types of cracks; at the same time, the crack formation around the hole under the high axial ratio was faster than that under the low axial ratio. In other words, when the same reflective tensile wave spread to the defect, the net tensile stress of the same magnitude was enough to cause an initial crack around the hole with a high axial ratio.
Based on this, the local stress concentration based on the same compressive and tensile stress waves was greater under a high axial ratio than that under a low axial ratio around the hole. The dynamic damage around the hole under a high axial ratio was stronger than that under a low axial ratio. Therefore, there was a larger amount and size of rock debris cuttings around the hole under the high axial ratio than that under the low axial ratio and that the crack initiation occurred earlier than that under low axial ratio.

3. Numerical Material Model and Validation

The experiment qualitatively explained that the failure point was related to the stress concentration, however, the form of destruction was not enough to absolutely ensure that the stress concentration distribution as stress wave passing through the hole could promote the fracture of the rock sample, or that the damage position had a correlation with the enhancement of local stress. In order to explain more meticulously that rock damage is closely related to axial ratio and dip angle, the finite element program LS-DYNA was employed to simulate the spalling process of rocks with a different elliptical hole, which has been proven to be suitable for rock simulation [26,27].
The model and test system were performed at a ratio of 1:1 and the experimental model is shown in Figure 6a, and the mesh around of the elliptical hole is shown in Figure 6b. The contact among the entities adopted single-sided automatic contact and the boundary was set as a free boundary.
Additionally, the elastic material model was used for the SHPB elastic bars. The continuous surface cap model (CSCM) was used to model the rock bar, which has been widely used for rocks [7,27]. Therefore, we used CSCM and the rock mechanic parameters to calculate the related simulation [27,28,29]. By comparing the results with the experiment, we aimed to verify the correctness of the material model. The damage patterns around the holes using numerical simulations and the experimental fracture are shown in Figure 7. The stress-time curve of the incident wave under the numerical simulation and experimental fracture are shown in Figure 8.
The above results indicate that the plastic deformation of the numerical simulation and experimental fracture around the hole were generally consistent. The period, amplitude, and variation trend were close to each other. In summary, the material model is suitable for simulation analysis of hard rock.

4. Numerical Simulation Results and Analyses

Based on the validated material model and its parameters, the damage around the elliptical holes was simulated under different axial ratios (m = 5, m = 2.5, m = 1.43, and m = 1) and different dip angles (β = 0°, β = 30°, β = 45°, β = 60°, and β = 90°). For each axial ratio, there were five groups of simulations under different dip angles except for a circular hole with m = 1. The failure mechanism of elliptical holes under different axial ratios and different dip angles were further analyzed.

4.1. Stress Evaluation and Failure Characteristics around the Hole

We defined the moment as the record zero point when the incident wave entered the specimen, based on the material parameters and geometry of the model, which was about 146 μs, and the incident wave spread from the left side of the hole to arrive at the hole. At approximately 300 μs, the incident wave spread through the hole. At about 386 μs, the incident wave completely passed through the hole. At about 438 μs, the reflective tensile wave spread via the reflection of the free surface to arrive at the right side of the hole. At about 678 μs, the reflective tensile wave completely passed through the hole.
With an increase in the axial ratio, the time of the failure of the hole moved ahead of schedule. Figure 9 shows that the larger the axial ratio, the greater the high effective strain produced by the compressive stress wave. When m = 5, the strain region that appeared at the upper and lower positions of the elliptical holes was larger than that of other axial ratios. These phenomena indicate that the axial ratio affected the dynamic stress concentration.
Figure 10 shows the change in the plastic zone around the hole with time when the axial ratio m = 2.5 and the dip angle β = 0°, β = 30°, β = 45°, β = 60°, and β = 90.
Figure 10 shows that the larger the dip angle, the greater the strain produced by the compressive stress wave passing the hole. The position of the damage is consistent with the position of the effective strain, which shows that the change of the dip angle has great influence on the rock failure around the hole.
Figure 9 and Figure 10 show that the effective strain around the hole due to the spread of the stress wave will cause a different magnitude of strain of the rock around the hole, so the position will also be different. With the arrival of the reflective tensile wave, the effective tensile strain gradually increased and finally reached a maximum value. Subsequently, intense strain occurred around the hole, causing plastic deformation that led to rock failure. Based on the damage position and changing region of the effective strain, the strain area was subjected to dynamic damage due to the stress concentration. When the tensile stress arrived, the first damage occurred in the damaged area, showing that the first-damage position is the weakest part of the rock.

4.2. Dynamic Stress Concentration Factor (DSCF) around the Hole

A stress field was gradually produced around the holes due to the scattering or diffraction of the incident wave, and the dynamic stress was concentrated during the dynamic loading process. The magnitude of the dynamic stress concentration can be characterized by the dynamic stress concentration factor (DSCF), which is determined by normalizing the hoop stress by the magnitude of the radial stress of the incident wave at the same point in the medium with no opening. [10]
k DSCF = σ θ σ r
where kDSCF is the dynamic stress concentration factor; kMDSCF denotes the maximum dynamic stress concentration factor; σ θ is the hoop stress value around the hole; and σ r is the axial stress propagated by stress waves when the medium has no holes.
When calculating the distribution of the DSCF around the hole, the calculation model is normally simplified as a plane problem to facilitate the analysis of the description of the phenomenon and results. The established plane coordinates and corresponding simplified model are shown in Figure 11.
The vector parameters r θ , r θ θ , and r are introduced here. The parameter | r θ | is defined as the distance between the center of the hole and any point of the surrounding rock; | r θ θ | is defined as the distance from the center of the hole to the hole boundary; and | r | is defined as the vertical distance from any point of the surrounding rock to the boundary of the hole. The relationship between these three parameters is as follows:
r θ = r θ θ + r

4.2.1. DSCF at the Same Dip Angle and Different Axial Ratios

Figure 12 and Figure 13 present the change of the hoop dynamic stress concentration factor with different axial ratios and the same dip angle around the holes under compressive and tensile stress waves when r = 0, i.e., r θ = r θ θ , respectively. Figure 14 shows the change of the maximum dynamic stress concentration factor with the axial ratio at the same dip angle.
Figure 12 and Figure 13 show the DSCF around the hole. The results indicate that the axial ratio was the main factor affecting the stress concentration distribution under the same dip angle, but that minimum stress occurred near θ = 0°.
Figure 14 shows the maximum DSCF value with varying axial ratios at the same dip angle. When the incident wave was reflected as tensile stress wave, the change in the maximum stress concentration was close to that of the compressive stress wave. The kDSCF first increased, and then decreased with the increase of the axial ratio.

4.2.2. DSCF at the Same Axial Ratio and Different Dip Angles

The previous section showed that the axial ratio had a significant effect on the DSCF around the hole. Under the same axial ratio, this section compares and analyzes the influence of the DSCF for each dip angle. Figure 15 and Figure 16 show the change of the DSCF for the same axial ratio and different dip angles around the holes under compressive and tensile stress waves, respectively, when r = 0 (i.e., r θ = r θ θ ).
Figure 17 shows the change in the maximum DSCF with the dip angle at the same axial ratio. Based on the change trend, the angles of the peak stress appeared to first decrease, and then increase with an increase in the dip angle of the hole.
Based on Figure 15, Figure 16 and Figure 17, the dynamic stress concentration factors for each axial ratio can be described in detail as follows:
First, when the axial ratio m = 1, the hole was circular, the stress distribution was not correlated with the dip angle β, and the maximum stress appeared at 90°. At this point, the change in the stress concentration under compressive and tensile stress waves was the same, but the maximum DSCF largely differed, which showed that kCMDSCF was 2.22 times higher than that of kTMDSCF. When the axial ratio m ≠ 1, the elliptical hole under the same axial ratio increased with the dip angle and the stress distribution changed from an initial y-axis symmetry to a central symmetry. When β = 90°, the distribution was symmetrical with the y-axis and the compressive stress wave had a larger range of stress concentration than the tensile stress wave.
Second, in the case of m = 1.43 and 2.5, kCMDSCF under the compressive stress wave resembled linear growth with a gradually increasing dip angle and the linear rate under the high axial ratio was higher than that under a lower axial ratio. The kTMDSCF under the tensile stress wave considerably fluctuated with the increasing dip angle and the whole trend of m = 1.43 steadily rose. However, in the case of m = 2.5, kTMDSCF was different before and after β = 60°. It resembles linear growth before β = 60° and the growth rate was large, while it decreased after β = 60° and the growth rate was small.
Third, when m = 5, the growth rate of kCMDSCF under the compressive stress wave notably varied before and after β = 60° with an increase in the dip angle of the hole. In the early stage, it showed rapid growth but turned flat in the later stage, with a small growth range. Under the function of the tensile stress wave, kTMDSCF rapidly decreased, reaching the maximum value at ~45°. The maximum DSCF increased rapidly before 45°, with a growth rate larger than that of other axial ratios, while it decreased after 45°, finally reaching the minimum value when β = 90°. At this point, the stress concentration was smaller than that under a small axial ratio when compared with the curve of the maximum stress concentration for the same axial ratio that changed with the dip angle under the compressive stress wave. When the compressive wave passes through the hole, it will cause initial damage, and the results indicated that the factors causing the decrease in the maximum DSCF related to the initial damage. In comparison with the change in the maximum DSCF of the three types of axial ratios corresponding to the dip angle of β = 90°, under the tensile stress wave, the result coincided with the situation in the previous section, that is, the larger the axial ratio, the smaller the tensile stress concentration. Therefore, the effect of compressive stress wave is one of the reasons for this situation. In addition, in the same axial ratios, with the increase in dip angle, the void area perpendicular to the direction of stress wave gradually increased, and the weak surface of the defect rock was larger and larger, so the increase in the weak surface also promotes the occurrence of destruction.

5. Conclusions

This study indicates that rock failure is affected by the dip angle and axial ratio of an elliptical hole. Under the same axial ratio, the change of the compressive and tensile stress concentrations varied with the increase in the dip angle of the holes. The compressive stress concentration showed an increasing trend, while the tensile stress concentration showed a decrease in the high axial ratio. The numerical results indicate that the stress concentration produced by the compressive stress wave passing through the hole promotes damage to the rock surrounding the hole. However, the damage is insufficient to solely cause the fractured rock of the holes. When the tensile stress wave arrives, the position of the dynamic tensile stress concentration coincides with the previous incident compressive stress wave.

Author Contributions

This paper is a result of the collaboration of all co-authors. M.T. conceived and designed the study; A.M. conducted the experimental, and drafted the manuscript; K.P. provided experimental guidance; Y.W. and K.D. revised the paper. All authors have read and approved the final manuscript.

Funding

The research presented in this paper was funded by the National Natural Science Foundation of China (Grant Nos. 41630642 and 11772357).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Xia, K.W.; Yao, W. Dynamic rock tests using split Hopkinson (Kolsky) bar system—A review. J. Rock Mech. Geotech. Eng. 2015, 7, 27–59. [Google Scholar] [CrossRef]
  2. Gu, H.L.; Tao, M.; Li, X.B.; Li, Q.Y.; Cao, W.Z.; Wang, F.; Zaretsky, E. , Dynamic response and failure mechanism of fractured coal under different soaking times. Theor Appl. Fract. Mech. 2018, 98, 112–122. [Google Scholar] [CrossRef]
  3. Zhu, J.B.; Liao, Z.Y.; Tang, C.A. Numerical SHPB tests of rocks under combined static and dynamic loading conditions with application to dynamic behavior of rocks under in situ stresses. Rock Mech. Rock Eng. 2016, 49, 3935–3946. [Google Scholar] [CrossRef]
  4. Liao, Z.Y.; Zhu, J.B.; Xia, K.W.; Tang, C.A. Determination of dynamic compressive and tensile behavior of rocks from numerical tests of split Hopkinson pressure and tension bars. Rock Mech. Rock Eng. 2016, 49, 3917–3934. [Google Scholar] [CrossRef]
  5. Li, D.Y.; Cheng, T.J.; Zhou, T.; Li, X.B. Experimental study of the dynamic strength and fracturing characteristics of marble specimens with a single hole under impact loading. Chin. J. Rock Mech. Eng. 2015, 34, 249–260. [Google Scholar]
  6. Tao, M.; Ma, A.; Cao, W.; Li, X.; Gong, F. Dynamic response of pre-stressed rock with a circular cavity subject to transient loading. Int. J. Rock Mech. Min. Sci. 2017, 99, 1–8. [Google Scholar] [CrossRef]
  7. Tao, M.; Zhao, H.; Li, X.; Li, X.; Du, K. Failure characteristics and stress distribution of pre-stressed rock specimen with circular cavity subjected to dynamic loading. Tunn. Undergr. Space Technol. 2018, 81, 1–15. [Google Scholar] [CrossRef]
  8. Pao, Y.H.; Mow, C.C.; Achenbach, J.D. Diffraction of Elastic Waves and Dynamic Stress Concentrations. J. Appl. Mech. 1973, 40, 872. [Google Scholar] [CrossRef]
  9. Liu, Q.J.; Zhang, C.; Todorovska, M.I. Scattering of SH waves by a shallow rectangular cavity in an elastic half space. Soil. Dynam. Earthq. Eng. 2016, 90, 147–157. [Google Scholar] [CrossRef]
  10. Yi, C.P.; Zhang, P.; Johansson, D.; Nyberg, U. Dynamic response of a circular lined tunnel with an imperfect interface subjected to cylindrical P-waves. Comput. Geotech. 2014, 55, 165–171. [Google Scholar] [CrossRef]
  11. Yi, C.P.; Lu, W.B.; Zhang, P.; Johansson, D.; Nyberg, U. Effect of imperfect interface on the dynamic response of a circular lined tunnel impacted by plane P-waves. Tunn. Undergr. Space. Technol. 2016, 51, 68–74. [Google Scholar] [CrossRef]
  12. Tao, M.; Li, Z.W.; Cao, W.Z.; Li, X.B.; Wu, C.Q. Stress redistribution of dynamic loading incident with arbitrary waveform through a circular cavity. Int. J. Numer. Anal. Methods Geomech. 2019, 43, 1279–1299. [Google Scholar] [CrossRef]
  13. Lin, G.; Li, Z.Y.; Li, J.B. A substructure replacement technique for the numerical solution of wave scattering problem. Soil. Dynam Earthq. Eng. 2018, 111, 87–97. [Google Scholar] [CrossRef]
  14. She, C.M.; Guo, W.L. Three-dimensional stress concentrations at elliptic holes in elastic isotropic plates subjected to tensile stress. Int. J. Fatigue 2007, 29, 330–335. [Google Scholar] [CrossRef]
  15. Zhou, X.L.; Wang, J.H.; Jiang, L.F. Dynamic response of a pair of elliptic tunnels embedded in a poroelastic medium. J. Sound Vib. 2009, 325, 816–834. [Google Scholar] [CrossRef]
  16. Deng, X.F.; Zhu, J.B.; Chen, S.G.; Zhao, Z.Y.; Zhou, Y.X.; Zhao, J. Numerical study on tunnel damage subject to blast-induced shock wave in jointed rock masses. Tunn. Undergr. Space Technol. 2014, 43, 88–100. [Google Scholar] [CrossRef]
  17. Du, K.; Su, R.; Tao, M.; Yang, C.Z.; Momeni, A.; Wang, S.F. Specimen shape and cross-section effects on the mechanical properties of rocks under uniaxial compressive stress. Bull Eng. Geol. Environ. 2019. [Google Scholar] [CrossRef]
  18. Fang, X.Q.; Jin, H.X. Dynamic response of a non-circular lined tunnel with visco-elastic imperfect interface in the saturated poroelastic medium. Comput. Geotech. 2017, 83, 98–105. [Google Scholar] [CrossRef]
  19. Kung, C.L.; Wang, T.T.; Chen, C.H.; Huang, T.H. Response of a Circular Tunnel Through Rock to a Harmonic Rayleigh Wave. Rock Mech. Rock Eng. 2018, 51, 547–559. [Google Scholar] [CrossRef]
  20. Li, X.B.; Li, C.J.; Cao, W.Z.; Tao, M. Dynamic stress concentration and energy evolution of deep-buried tunnels under blasting loads. Int. J. Rock Mech Min. Sci. 2018, 104, 131–146. [Google Scholar] [CrossRef]
  21. Wang, Y.B.; Yang, R.S.; Zhao, G.F. Influence of empty hole on crack running in PMMA plate under dynamic loading. Polym. Test. 2017, 58, 70–85. [Google Scholar] [CrossRef]
  22. Wang, P.F.; Li, Y.G.; Liu, R.H.; Shi, Y.J. Effect of forced-to-exhaust ratio of air volume on dust control of wall-attached swirling ventilation for mechanized excavation face. Tunn. Undergr. Space Technol. 2019, 90, 194–207. [Google Scholar] [CrossRef]
  23. Li, X.B.; Zhou, Z.L.; Lok, T.S.; Hong, L.; Yin, T.B. Innovative testing technique of rock subjected to coupled static and dynamic loads. Int. J. Rock Mech Min. Sci. 2008, 45, 739–748. [Google Scholar] [CrossRef]
  24. Liu, Z.H.; Du, C.L.; Jiang, H.X.; Liu, K. Analysis of Roadheader for Breaking Rock Containing Holes under Confining Pressures. Energies 2017, 10, 1154. [Google Scholar] [CrossRef]
  25. Peng, K.; Zhou, J.Q.; Zou, Q.L.; Zhang, J.; Wu, F. Effects of stress lower limit during cyclic loading and unloading on deformation characteristics of sandstones. Constr. Build. Mater. 2019, 217, 202–215. [Google Scholar] [CrossRef]
  26. Guo, J.S.; Ma, L.Q.; Wang, Y.; Wang, F.T. Hanging Wall Pressure Relief Mechanism of Horizontal Section Top-Coal Caving Face and Its Application—A Case Study of the Urumqi Coalfield, China. Energies 2017, 10, 1371. [Google Scholar] [CrossRef]
  27. Tao, M.; Li, X.B.; Li, D.Y. Rock failure induced by dynamic unloading under 3D stress state. Theor. Appl. Fract. Mech. 2013, 65, 47–54. [Google Scholar] [CrossRef]
  28. Tao, M.; Li, X.B.; Wu, C.Q. Characteristics of the unloading process of rocks under high initial stress. Comput. Geotech. 2012, 45, 83–92. [Google Scholar] [CrossRef]
  29. Peng, K.; Liu, Z.P.; Zou, Q.L.; Zhang, Z.Y.; Zhou, J.Q. Static and dynamic mechanical properties of granite from various burial depths. Rock Mech. Rock Eng. 2019, 1–22. [Google Scholar] [CrossRef]
Figure 1. Granite mineral identification map and mineral composition map. (a) Sheet identification map (pad plasterboard). (b) Sample mineral composition map.
Figure 1. Granite mineral identification map and mineral composition map. (a) Sheet identification map (pad plasterboard). (b) Sample mineral composition map.
Energies 12 03441 g001
Figure 2. The actual and schematic views of modified SHPB device test system. (a) Actual view; (b) Schematic view.
Figure 2. The actual and schematic views of modified SHPB device test system. (a) Actual view; (b) Schematic view.
Energies 12 03441 g002aEnergies 12 03441 g002b
Figure 3. The complete rock spalling failure process.
Figure 3. The complete rock spalling failure process.
Energies 12 03441 g003
Figure 4. Fracture around the hole under different axial ratios at the same time (β = 90°).
Figure 4. Fracture around the hole under different axial ratios at the same time (β = 90°).
Energies 12 03441 g004
Figure 5. Final failure modes of partial holes under different axial ratios and dip angles.
Figure 5. Final failure modes of partial holes under different axial ratios and dip angles.
Energies 12 03441 g005
Figure 6. The model geometry and corresponding elliptical hole mesh. (a) Model geometry, (b) Elliptical hole mesh.
Figure 6. The model geometry and corresponding elliptical hole mesh. (a) Model geometry, (b) Elliptical hole mesh.
Energies 12 03441 g006
Figure 7. Comparison between the numerical simulation and experimental fracture around the hole (The green part of the figure is the plastic deformation, and the blue part is the elastic deformation).
Figure 7. Comparison between the numerical simulation and experimental fracture around the hole (The green part of the figure is the plastic deformation, and the blue part is the elastic deformation).
Energies 12 03441 g007
Figure 8. Comparison of the experimental and simulation incident wave patterns.
Figure 8. Comparison of the experimental and simulation incident wave patterns.
Energies 12 03441 g008
Figure 9. Distribution evolution diagram of the plastic variation of holes with different axial ratios but the same dip angle (β = 90°).
Figure 9. Distribution evolution diagram of the plastic variation of holes with different axial ratios but the same dip angle (β = 90°).
Energies 12 03441 g009aEnergies 12 03441 g009b
Figure 10. Distribution evolution diagram of the plastic strain variation of holes with different dip angles and the same axial ratio (m = 2.5).
Figure 10. Distribution evolution diagram of the plastic strain variation of holes with different dip angles and the same axial ratio (m = 2.5).
Energies 12 03441 g010aEnergies 12 03441 g010b
Figure 11. Simplified coordinate model of the incident waves and prefabricated ellipse holes.
Figure 11. Simplified coordinate model of the incident waves and prefabricated ellipse holes.
Energies 12 03441 g011
Figure 12. The DSCF with different axial ratios around the holes under a compressive stress wave ( r = 0).
Figure 12. The DSCF with different axial ratios around the holes under a compressive stress wave ( r = 0).
Energies 12 03441 g012
Figure 13. The DSCF with different axial ratios around the holes under a tensile stress wave ( r = 0).
Figure 13. The DSCF with different axial ratios around the holes under a tensile stress wave ( r = 0).
Energies 12 03441 g013
Figure 14. The maximum DSCF with varying axial ratio at the same dip angle. (a) Compressive stress waves; (b) Tensile stress waves.
Figure 14. The maximum DSCF with varying axial ratio at the same dip angle. (a) Compressive stress waves; (b) Tensile stress waves.
Energies 12 03441 g014
Figure 15. The DSCF at the same axial ratio and different dip angles under the compressive stress wave ( r = 0).
Figure 15. The DSCF at the same axial ratio and different dip angles under the compressive stress wave ( r = 0).
Energies 12 03441 g015aEnergies 12 03441 g015b
Figure 16. The DSCF at the same axial ratio and different dip angles under the tensile stress wave ( r = 0).
Figure 16. The DSCF at the same axial ratio and different dip angles under the tensile stress wave ( r = 0).
Energies 12 03441 g016
Figure 17. The maximum DSCF at the dip angle and axial ratio. (a) Compressive stress waves, (b) Tensile stress waves.
Figure 17. The maximum DSCF at the dip angle and axial ratio. (a) Compressive stress waves, (b) Tensile stress waves.
Energies 12 03441 g017

Share and Cite

MDPI and ACS Style

Tao, M.; Ma, A.; Peng, K.; Wang, Y.; Du, K. Fracture Evaluation and Dynamic Stress Concentration of Granite Specimens Containing Elliptic Cavity under Dynamic Loading. Energies 2019, 12, 3441. https://doi.org/10.3390/en12183441

AMA Style

Tao M, Ma A, Peng K, Wang Y, Du K. Fracture Evaluation and Dynamic Stress Concentration of Granite Specimens Containing Elliptic Cavity under Dynamic Loading. Energies. 2019; 12(18):3441. https://doi.org/10.3390/en12183441

Chicago/Turabian Style

Tao, Ming, Ao Ma, Kang Peng, Yiqing Wang, and Kun Du. 2019. "Fracture Evaluation and Dynamic Stress Concentration of Granite Specimens Containing Elliptic Cavity under Dynamic Loading" Energies 12, no. 18: 3441. https://doi.org/10.3390/en12183441

APA Style

Tao, M., Ma, A., Peng, K., Wang, Y., & Du, K. (2019). Fracture Evaluation and Dynamic Stress Concentration of Granite Specimens Containing Elliptic Cavity under Dynamic Loading. Energies, 12(18), 3441. https://doi.org/10.3390/en12183441

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop