Next Article in Journal
Additive Manufacturing of CrFeNiTi Multi-Principal Element Alloys
Previous Article in Journal
Hydrotalcite-Modified Clinoptilolite as the Catalyst for Selective Catalytic Reduction of NO with Ammonia (NH3-SCR)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Studying the Characteristics of Chaos and Fractals of Construction Rocks under Different Loading Velocities

1
Reserch Center for Wind Engineering and Engineering Vibration, Guangzhou University, Guangzhou 510006, China
2
School of Physics and Electronic Science, Qiannan Normal University for Nationalities, Duyun City 558000, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(22), 7890; https://doi.org/10.3390/ma15227890
Submission received: 14 October 2022 / Revised: 31 October 2022 / Accepted: 4 November 2022 / Published: 8 November 2022

Abstract

:

Highlights

  • Rock materials exhibit obvious strain rate effects at different loading velocities.
  • The damage evolution process of rock materials contains chaotic features.
  • Logical mapping is used to analyze rock material damage–strain information sequences.
  • The fractal dimension of rock sample crack is positively correlated with load velocity.

Abstract

Rock is a widely used construction material; its mechanical properties change due to the influence of different load speed. In this study, the split Hopkinson pressure bar (SHPB) was used to test the dynamic properties of rock samples by loading four different pressures (0.05, 0.08, 0.14, and 0.23 MPa). The peak stress of the sample increases from 82.19 to 284.16 MPa, and the particle size of the sample debris decreases from 46.57 to 18.34 mm as the impact pressure increases from 0.05 to 0.23 MPa. As a chaos method in nonlinear dynamics, it is introduced into the quantitative evaluation of the sample at four loading pressures, which is then calculated. The damage evolution process of the sample under four loading pressures is calculated, and the chaotic characteristics contained in the process are analyzed. Based on the logistic mapping, the increase in the load velocity can delay the entry of the damage variable into the period-doubling bifurcation and chaotic states. Finally, the fractal dimension of the rock crack at the corresponding time under different load speeds is calculated, and the results showed that the increase in the load velocity can increase the uniformity of the crack distribution.

1. Introduction

Construction rock, as a kind of widely distributed natural geological material, has been widely used in engineering construction [1,2]. Rock materials often contain defects due to long-term geological effects, which show significant ambiguity, nonlinearity and anisotropy on the macroscopic scale [3,4,5,6]. It is an effective method to study the variation characteristics of rock materials on dynamic loading by using chaos theory in nonlinear dynamics.
As a typical form of motion of nonlinear systems, chaotic behavior is characterized by sensitivity to initial values and has very complex fractal and self-similar structures [7,8]. Chaos discrimination based on the Lyapunov index is used to observe information sequences of rock materials and is considered to be an effective method to analyze the failure properties of construction rocks [9]. As a typical representative, Liu [10], based on the Logistics equation, observed the damage evolution sequence in the process of rock loading failure and analyzed the chaotic characteristics of the rock material deformation and failure processes under static loading conditions, and obtained good results. Meanwhile, in fractal theory, Sun et al. [11] provided quantitative guidance for analyzing the stability of a rock mass during excavation engineering by fractal characterization, and Ning et al. [12] studied the fractal characteristics of fatigue failure of coal rock samples. Both studies successfully analyzed the characteristics of rock failure.
Construction rock not only accounts for the action of static load, but it also bears the effect of a wide range of dynamic loads [13,14]. However, the rock material will exhibit significant strain rate effects on dynamic loading, and the properties of the rock will also be significantly different from those on the static loading [15,16,17,18,19,20]. Sunita et al. [18] studied the dynamic loading properties of three kinds of rocks (compact basalt, metadolerite and granite) by SHPB, electron microscope and X-ray equipment. It was concluded that the dynamic peak strength increased by 4.14–7.52 times in the range of the test strain rate. Meanwhile, Millon et al. [21] filled the gap of information that exists regarding dynamic behavior of sedimentary rocks under strain rates between 100 and 5.2 × 102 s−1. Gong et al. [20] analyzed the mechanical properties and failure behavior of rock-like materials in the range of 2.7 to 4.0 × 105 MPa/s, where an increase in peak strength and elastic modulus of rock together with a decrease in peak strain and dynamic failure duration in the post-peak period are shown.
Xia et al. [22] described the development and research progress of dynamic compression, dynamic Brazil test and dynamic notched semi-circular bend (three rock dynamic testing methods recommended by the International Society for Rock Mechanics) and summarized the different properties of rock materials under dynamic loading. Zhang et al. [23] summarized the research work in this field, and some popular semi-empirical rate-dependent equations for the enhancement of dynamic mechanical properties were given. The above research work explored the properties of the rock materials on dynamic loading conditions. It is very meaningful research work, and the theory of chaotic fractals can provide an effective exploration path for this research work.
In this paper, the split Hopkinson pressure bar (SHPB) test equipment is used to analyze the characteristics of rock materials at different loading speeds. The evolution law of chaos contained in the damage–strain relation is analyzed. At the same time, the variation of dimension damage of rock samples at different loading speeds is studied by fractal theory.

2. SHPB Test

2.1. Sample Preparation

Construction rock samples were drilled from Xia’men mining area, China. The samples were Φ80 × 40 mm cylinders that exhibited scratches smaller than 0.025 mm with less than 0.025 mm deviation in the flatness. Meanwhile, the samples were cleaned and dried in a natural environment, as shown in Figure 1. These samples were tested in four series (DT-A/B/C/D), depending on the test conditions.

2.2. Experimental Setup and Test Method

The split Hopkinson pressure bar (SHPB) test device with a diameter of 100 mm was used in this study, as shown in Figure 2 (bullet length: 800 mm, incident rod length: 4100 mm, transmission rod length: 3100 mm). Butyl rubber wafers (diameter: 30 mm, thickness: 1.5 mm) were used as waveform shapers. The purpose of waveform shaping is to make the waveform approximate to the sinusoidal curve shape and to lengthen the rising time of the waveform, so as to help the sample achieve stress balance as soon as possible during loading. In the long-term use of the equipment (SHPB), this specification waveform shaper was found to meet the test requirements very well.
Before each test, the impact load is carried out to determine the interference degree of noise to the circuit (strain gauge and Wheatstone bridge), and to remove high frequency noise by adjusting the low-pass filter in the dynamic strain gauge to ensure credibility of the experimental data.
In this test, four different loading pressures (0.05, 0.08, 0.14, and 0.23 MPa) were set up to produce different stress wave amplitudes. The samples were tested at each impact speed, and multiple tests were performed under each operating condition to omit outliers and to ensure good test repeatability. A high-speed camera system was used to record the destruction process of the samples.
The screening tool used in the test is shown in Figure 2b. The apertures were 2.0, 5.0, 10.0, and 20 mm, and the sample was sieved into 0.0–2.0, 2.0–5.0, 5.0–10.0, 10.0–20.0, and 20.0–80.0 mm (80 mm is the diameter before the sample is broken) groups of particle size. The five groups of particles were numbered according to particle size from small to large (I = 1, 2, 3, 4, 5). The weighing instrument was a high-precision electronic scale (Model: ACS-S) with a range of 1 kg and an accuracy of 0.1 g.

3. Results

3.1. Stress–Strain Curves

The stress–strain curve is an important expression of information reflecting the deformation and failure process of rock under external load. It is of great significance to accurately understand the physical and mechanical properties of rock materials and to guide engineering practice.
The rock sample exhibited a significant strain rate effect with the increase in load velocity, as shown in Figure 3 and Figure 4 and Table 1. With the impact pressure increases from 0.05 to 0.08, 0.14, and 0.23 MPa, the peak stress (mean) increases from 82.19 to 163.18, 229.62, and 284.16 MPa, respectively, increases of 98.54%, 179.38%, 245.74%, and the process shows the law of power function ( y = 30300.56 x 0.00473 + 29816.3 ), This is because the sample is input with more energy per unit time at high-speed loading conditions, and the energy dissipation rate of rock materials is higher than that at low-speed loading conditions, such that the crack along the weakest path at a low strain rate evolved along the maximum energy dissipation path at a high strain rate. The cracks in the rocks continued to change from intergranular to transgranular, and the peak intensity of the sample on the macro scale increased, showing a significant strain rate effect [24,25,26].

3.2. The Sample Debris after Loading

After the sample is damaged by the impact load, the relationship between the number and the size of blocks with the loading speed is one of the basic problems of rock dynamics. The fracture state of rock samples is different at different loading speeds, and the fragmentation degree of broken pieces often has certain mass distribution characteristics. In order to characterize this relationship, the average block size is used as a measure:
x ¯ = x i η i η i
where x ¯ is the average size of the broken sample (mm); x i is the average size of debris remaining on each stage of the screen (mm); η i is the percentage of the mass of the debris when the average size is x i , %.
The broken rock samples are collected one by one, sieved and weighed one by one. The typical distribution of the sample debris after screening and separation at different impact speeds is shown in Figure 5a. The percentage of the mass fraction of the debris is gradually changed with the increase in the impact pressure (the average particle diameter is 46.57 mm at 0.05 MPa; 38.96 mm at 0.08 MPa; 30.83 mm at 0.14 MPa; and 18.34 mm at 0.23 MPa). The variation of the particle size of the rock sample is shown in Figure 5b and Table 2.
As the loading speed increases, the average particle size of the rock debris decreases gradually, and the degree of fracture is more severe. However, the cracks are basically carried out parallel to the loading direction, and the broken rock debris mostly exhibits a prism shape without significant change. The number of small-sized rocks as debris increases sharply when the loading speed increases, mainly because the cracks are more densely developed and the friction between the prisms is more severe.

4. Discussion

4.1. The Chaotic Characteristics of Rock Sample Damage Evolution

4.1.1. Rock Damage

Liu [10] obtained the rock damage evolution law based on the Lemaitre strain equivalent hypothesis [26] and Weibull distribution (Equation (2)).
D = k 1 + e a r ε ,   0 < k 1
In Equation (2), D is the damage variable ( D = 1 A A , A is the initial cross-sectional area of the non-destructive material, A is the effective bearing area (net area) after the material is damaged); ε is the strain; a is the parameter reflecting the degree of initial damage, r is the internal growth rate of damage evolution, and k is the maximum damage degree.
Typical rock sample stress–strain and damage–strain curves at different loading speeds are shown in Figure 3a–d. As the impact pressure increases from 0.05 to 0.08, 0.14, and 0.23 MPa, the stress value of the damage initiation time of the sample increases from 36.1 to 75.0, 119.4, and 203.3 MPa, and the ratio of the value to the corresponding peak intensity increases from 42.3% to 45.1%, 50.2%, and 73.8%, respectively. On the contrary, the corresponding strain at this time was reduced from 0.00262 to 0.00168, 0.00079, and 0.00039, respectively. The higher the loading speed, the shorter the destruction process of the sample. The damage value of the sample is increased from 0.589 to 0.715, 0.972, and 0.996, as shown in Figure 6.
This is because the sample inputs more energy per unit time as the impact velocity increases, the crack that is along the weakest path at the low speed load is converted along the maximum energy-consuming path, a large number of tiny cracks are not fully expanded, and macroscopic cracks caused by transgranular fracture are formed. The appearance of macroscopic cracks causes the material to be in a state of failure. Therefore, the strain corresponding to the peak stress moment of the sample will be smaller at high-speed loading. Meanwhile, the phenomenon of the damage–strain curve of the sample gradually changes from a concave curve into a convex curve, due to the damage surge phase changes from one peak to the peak before.

4.1.2. Chaotic Characteristics

As a heterogeneous natural material at the mesoscale, the dynamic response of rock under high-speed impact is a nonlinear behavior. Based on the analysis method of chaos theory, the assumption that rock materials are isotropic is the premise for subsequent analyses. The deformation and failure process of rock at the external loading can be obtained from the sequence information of the experiment, so as to study the characteristics of the rock damage. Deriving the damage evolution Equation (2), the differential form of the damage evolution equation can be obtained [10]:
d D d ε = r D ( 1 D k )
The above equation can be written as the following iterative expression:
D n + 1 = μ D n ( 1 D n ) = f ( D n ) ,   n = 0 , 1 , 2 , 3 , μ = r Δ ε + 1 ,   D n = r Δ ε ( r Δ ε + 1 ) k D n
The mapping f is a logistic mapping reflecting the evolution of the generalized damage variable D n . Since the function f ( D n ) is nonlinear, and the inverse map D n = f ( D n + 1 ) cannot be defined singly, such nonlinear maps are irreversible [27]. For the control parameters μ [ 0 , 4 ] , x n + 1 = f ( x n ) [ 0 , 1 ] is the invariant interval of the logistic map. In this interval, as μ increases, the logistic map undergoes a period-doubling bifurcation and finally enters a chaotic state, as shown in Figure 7.
As from Figure 3, based on the experimental data, the period-doubling bifurcation points A1~4( μ [ 3 , 3.5699 ) ) of the generalized damage variable at different loading speeds are obtained, and the stress value B1~4 corresponding to point A1~4 is also obtained. At the same time, the starting point C1~4( μ [ 3.5699 , 4 ) ) and the corresponding stress values D1~4 of the generalized damage variable entering the chaotic state are obtained.
From Figure 8a, the stress value corresponding to the generalized damage variable at the period-doubling bifurcation points to the stress peak increases with the increase in the loading speed; when the impact pressure is 0.05 MPa, the ratio is 87.8%. As the impact pressure is increased to 0.08, 0.14, and 0.23 MPa, respectively, the value becomes 95.4%, 97.5%, and 97.7% (Figure 8a, black fold line). The ratio of the stress value to the peak stress corresponding to the chaotic starting point of the generalized damage variable has a similar law (Figure 8a red line). Increases in the impact velocity can delay the generalized damage variable enter the periods of period-doubling bifurcation and chaotic states.

4.2. Effect of Loading Velocity on Chaotic Characteristics

The damage–strain cures of the sample from the initial damage time to the start of the period-doubling bifurcation at different loading speeds were treated as dimensionless; the data were fitted separately, as shown in Figure 8b. It can be found that the difference in the loading speed has little influence on the curve form of the stage. The damage–strain curve shows linear growth, and the evolution law (trend) shows similarity. The difference in the fitted linear slope is less than 7%.
This is because the voids inside the rock material are tightly closed at this stage. Usually only a small amount of plastic deformation occurs, which is the initial stage of the damage process. The increase in the load speed only affects the duration of the process and does not affect the damage mode, and the generalized damage variable attracts a fixed point D * = μ 1 μ ( lim n + f n ( D 0 ) = μ 1 μ ) in the logistic equation. Due to the strain rate effect of the rock material, the difference in duration makes the stress of the sample different at the end of the stage.

4.3. The Fractal Characteristics of Rock Sample Damage Evolution

4.3.1. The Damage Fractal Dimension

The damage material is (equivalent to) the porous media, and the damage evolution process is equivalent to the micropore evolution process in the material. The fractal dimension of the rock crack reflects the characteristics of the crack and the complexity of the section, which directly affect the judgment of the properties of the rock materials. The Hausdorff dimension (Equation (5)) and the box dimension method are used to calculate the damage situation corresponding to the period-doubling bifurcation and the chaotic states while using the fractal geometry covering method to quantitatively investigate the statistical self-similarity of the rock crack [27,28,29,30,31].
D H = lim δ 0 ln N ( δ ) ln ( 1 / δ )
where δ is the metric indicating the measure; N ( δ ) represents the number of targets covered by the scale; D H is the covered dimension
A FASTCAM SA-Z (200K-C-32 GB) high-speed camera system (Japan) was used to record the failure pattern of the samples during the loading process; photos are shown in Figure 9. Without obstructions on the sides of the samples in the SHPB apparatus, the failure mode of the half-circumference surface was randomly analyzed by the high-speed camera. To depict the development of rock cracks at different impact velocities, the high-speed camera system was used to photograph the sample failure process. Resultant images were processed by binary image processing technology [32], as shown in Figure 9.
The damage evolution process of rock sample is analyzed by the fractal geometry method. The fractal dimension values of the sample at the moment of period-doubling bifurcation and chaotic states are calculated, the crack binary images corresponding to the above moments are the measurement objects, with log (1/L) as the abscissa and log N(L) as the ordinate, and the log (1/L) and log N (L) curves are drawn, as shown in Figure 9a–h. As can be seen from the figure, there is a linear correlation between log (1/L) and log N(L), and the slope is the fractal dimension of each stress state.

4.3.2. The Variation of Fractal Dimension

At different loading speeds, the sample damage evolution process has fractal characteristics and statistical self-similarity at the moment of period-doubling bifurcation and chaotic states because log (1/L) and log N (L) exhibit a linear correlation.
Rock materials have initial damage, and their damage fractal dimension is generally greater than 1. When the damage is evenly distributed absolutely, the fractal dimension is 3. When the damage distribution is more localized, the value is closer to 1; on the contrary, when the damage distribution is less localized, the value is closer to 3. Therefore, the fractal dimension of the damage reflects the degree of localization of the damage distribution.
As shown in Figure 10, the damage fractal dimension value of the same sample is usually 0.26–0.36 times larger at the chaos initial moment than the initial moment of the period-doubling bifurcation, which indicates that the crack of the sample has a more discrete distribution during the whole loading process. The fractal dimension increases with the increase in impact pressure, from 2.07 at 0.05 MPa to 2.44 at 0.23 MPa, increasing by 17.87% at the start time of the period-doubling bifurcation state. The fractal dimension at the beginning of the chaotic state also shows the same increasing trend with the increase in impact pressure, which increases from 2.33 at 0.05 MPa to 2.70 at 0.23 MPa, and increases by 15.88%. This indicates that the increase in impact velocity obviously improves the homogenization degree of crack development in different stages of rock sample damage and failure process [33,34,35].

5. Conclusions

In this paper, the dynamic failure properties of the rock samples under different impact velocities were analyzed by the SHPB test system and sizing tests, based on chaos and fractal theory, the damage process and characteristics of the sample were analyzed, and the following conclusions were drawn:
  • Construction rock exhibits significant strain rate effects under dynamic loading conditions. As the impact pressure increases from 0.05 to 0.23 MPa, the peak stress of the sample increases from 82.19 to 284.16 MPa (an increase of 245.74%), and the increase process is in accordance with the power function law; the average particle size of the sample debris is reduced from 46.57 to 18.34 mm (reduced by 60.62%).
  • The damage–strain evolution information of the construction rock shows chaotic characteristics. The increase in the load velocity can delay the evolution of the generalized damage variable into the bifurcation state and the chaotic state.
  • The damage–strain evolution law (trend) of the sample from the initial moment of damage to the start of the period of bifurcation at different loading pressures shows similarity, and the slope of the fitted straight line is less than 7%.
  • Under dynamic loading conditions, the damage evolution process of the samples has fractal characteristics at the start of the bifurcation and chaotic periods, and the damage cracks have statistical self-similarity. Meanwhile, the increase in load velocity significantly increases the degree of homogenization of crack development in the process of damage.

Author Contributions

Conceptualization, N.W.; methodology, N.W.; formal analysis, N.W.; investigation, N.W.; resources, J.F.; data curation, C.X.; writing—original draft preparation, N.W.; writing—review and editing, N.W.; supervision, J.F.; funding acquisition, J.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China Youth Fund (grant no. 52108448).

Data Availability Statement

The introduction data supporting this manuscript are from previously reported studies and datasets, which have been cited. The processed data are available from the corresponding author upon request. The raw test data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they do not have any commercial or associated interests that might represent a conflict of interest in connection with the work submitted.

References

  1. Li, D.X.; Wang, E.Y.; Kong, X.G.; Jia, H.S.; Wang, D.M.; Muhammad, A. Damage precursor of construction rocks under uniaxial cyclic loading tests analyzed by acoustic emission. Constr. Build. Mater. 2019, 206, 169–178. [Google Scholar]
  2. Gu, Q.H.; Ma, Q.; Tan, Y.L.; Jia, Z.C.; Zhao, Z.H.; Huang, D.M. Acoustic emission characteristics and damage model of cement mortar under uniaxial compression. Constr. Build. Mater. 2019, 213, 377–385. [Google Scholar] [CrossRef]
  3. Sirdesai, N.N.; Gupta, T.; Singh, T.N.; Ranjith, P.G. Studying the acoustic emission response of an Indian monumental sandstone under varying temperatures and strains. Constr. Build. Mater. 2018, 168, 346–361. [Google Scholar] [CrossRef]
  4. Sun, B.; Chen, R.; Ping, Y.; Zhu, Z.; Wu, N.; He, Y. Dynamic Response of Rock-like Materials Based on SHPB Pulse Waveform Characteristics. Materials 2022, 15, 210. [Google Scholar] [CrossRef]
  5. Sun, B.; Chen, R.; Ping, Y.; Zhu, Z.; Wu, N.; Shi, Z. Research on Dynamic Strength and Inertia Effect of Concrete Materials Based on Large-Diameter Split Hopkinson Pressure Bar Test. Materials 2022, 15, 2995. [Google Scholar] [CrossRef]
  6. Sun, B.; Chen, R.; Ping, Y.; Zhu, Z.; Wu, N. Study on Axial Tensile Strain Rate Effect on Concrete Based on Experimental Investigation and Numerical Simulation. Materials 2022, 15, 5164. [Google Scholar] [CrossRef]
  7. Hirata, T. Fractal structure of spatial distribution earthquake-The two point correlation function. Geophys. J. R. Astron. Soc. 1980, 62, 303–320. [Google Scholar]
  8. Chen, Y.F.; Lu, J.H.; Zhou, C.B. Predicition of slope displacement by using Lyapunov exponent improved technique. Chin. J. Rock Mech. Eng. 2001, 20, 671–675. (In Chinese) [Google Scholar]
  9. Wang, L.G.; Song, Y.; Miao, X.X. Lyapunov exponent of chaos feature in the processing of deformation failure for coal floor. Chin. J. Geotech. Eng. 2002, 24, 356–359. (In Chinese) [Google Scholar]
  10. Liu, D.Q. Study on Rock Damage Constitutive Model and Chaos Characteristics in the Process of Deformation and Failure. Ph.D. Thesis, China University of Mining and Technology, Beijing, China, 2014. [Google Scholar]
  11. Sun, H.; Liu, X.L.; Zhu, J.B. Correlational fractal characterisation of stress and acoustic emission during coal and rock failure under multilevel dynamic loading. Int. J. Rock Mech. Min. 2019, 117, 1–10. [Google Scholar] [CrossRef]
  12. ZZhu, F.; Wu, X.; Zhou, M.; Sabri, M.M.S.; Huang, J. Intelligent Design of Building Materials: Development of an AI-Based Method for Cement-Slag Concrete Design. Materials 2022, 15, 3833. [Google Scholar] [CrossRef] [PubMed]
  13. Geng, J.S.; Sun, Q.; Zhang, Y.C.; Cao, L.W.; Zhang, W.Q. Studying the dynamic damage failure of concrete based on acoustic emission. Constr. Build. Mater. 2017, 149, 9–16. [Google Scholar] [CrossRef]
  14. Wu, N.; Zhu, Z.D.; Zhu, C.; Luo, Z.H. Dynamic Behavior of Rock Joint Under Different Impact Loads. KSCE J. Civ. Eng. 2019, 23, 541–548. [Google Scholar] [CrossRef]
  15. Fakhimi, A.; Azhdari, P.; Kimberley, J. Physical and numerical evaluation of rock strength in Split Hopkinson Pressure Bar testing. Comput. Geotech. 2018, 102, 1–11. [Google Scholar] [CrossRef]
  16. Huang, S.; Chen, R.; Xia, K.W. Quantification of dynamic tensile parameters of rocks using a modified Kolsky tension bar apparatus. J. Rock Mech. Geotech. Eng. 2010, 20, 162–168. [Google Scholar] [CrossRef] [Green Version]
  17. Feng, P.; Dai, F.; Liu, Y.; Xu, N.; Fan, P. Effects of coupled static and dynamic strain rates on the mechanical behaviors of rock-like specimens containing preexisting fissures under uniaxial compression. Can. Geotech. J. 2018, 55, 640–652. [Google Scholar] [CrossRef] [Green Version]
  18. Mishra, S.; Chakraborty, T.; Matsagar, V.; Loukus, J.; Bekkala, B. High strain rate charactterization of deccan trap rocks using SHPB device. J. Mater. Civ. Eng. 2017, 30, 04018059. [Google Scholar] [CrossRef]
  19. Li, H.M.; Li, H.G.; Gao, B.B.; Jiang, D.J.; Feng, J.F. Study of acoustic emission and mechanical characteristics of coal samples under different loading rates. Shock. Vib. 2015, 2015, 458519. [Google Scholar] [CrossRef] [Green Version]
  20. Gong, F.Q.; Ye, H.; Luo, Y. The effect of high loading rate on the behaviour and mechanical properties of coal-rock combined body. Shock. Vib. 2018, 2018, 4374530. [Google Scholar] [CrossRef] [Green Version]
  21. Millon, O.; Ruiz-Ripoll, M.L.; Hoerth, T. Analysis of the Behavior of Sedimentary Rocks Under Impact Loading. Rock Mech. Rock Eng. 2016, 49, 4257–4272. [Google Scholar] [CrossRef]
  22. 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] [Green Version]
  23. Zhang, Q.B.; Zhao, J. A Review of Dynamic Experimental Techniques and Mechanical Behaviour of Rock Materials. Rock Mech. Rock Eng. 2014, 47, 1411–1478. [Google Scholar] [CrossRef] [Green Version]
  24. Wu, N.; Zhu, Z.D.; Zhou, Y.J.; Gao, S.H. A Comparative Study on Rock Properties in Split and Compression Dynamic Test. Shock. Vib. 2018, 2018, 2861537. [Google Scholar]
  25. Hokka, M.; Black, J.; Tkalich, D.; Fourmeau, M.; Kane, A.; Hoang, H.; Li, C.; Chen, W.W.; Kuokkala, V.T. Effects of strain rate and confining pressure on the compressive behavior of Kuru granite. Int. J. Impact Eng. 2016, 91, 183–193. [Google Scholar] [CrossRef]
  26. Lemaitre, J. A continuous damage mechanics model for ductile materials. J. Eng. Mater. Technol. 1985, 107, 83–89. [Google Scholar] [CrossRef]
  27. Gabriel, J.P.; Saucy, F.; Bersier, L.F. Paradoxes in the logistic equation. Ecol. Model. 2005, 185, 147–151. [Google Scholar] [CrossRef] [Green Version]
  28. Chen, Y.S. Studies on Damage Evolution Characteristics in Non-Interpenetrated Material under Axial Compressive Condition. Master’s Thesis, Xi’an University of Technology, Xi’an, China, 2002. [Google Scholar]
  29. Balka, R.; Buczolich, Z.; Elekes, M. A new fractal dimension: The topological Hausdorff dimension. Adv. Math. 2015, 274, 881–927. [Google Scholar] [CrossRef] [Green Version]
  30. Fernández-Martínez, M.; Sánchez-Granero, M.A. How to calculate the Hausdorff dimension using fractal structures. Appl. Math. Comput. 2015, 264, 116–131. [Google Scholar] [CrossRef]
  31. Croydon, D.A. The Hausdorff dimension of a class of random self-similar fractal trees. Adv. Appl. Probab. 2007, 39, 708–730. [Google Scholar] [CrossRef] [Green Version]
  32. Luo, Z.H.; Zhu, Z.D.; Ruan, H.N.; Shi, C. Extraction of microcracks in rock images based on heuristic graph searching and application. Comput. Geosci.-UK 2015, 85, 22–35. [Google Scholar] [CrossRef]
  33. Tavares, L.M. Chapter 1 Breakage of Single Particles: Quasi-Static. Handb. Powder Technol. 2007, 12, 3–68. [Google Scholar]
  34. Tavares, L.M.; King, R.P. Modeling of particle fracture by repeated impacts using continuum damage mechanics. Powder Technol. 2002, 123, 138–146. [Google Scholar] [CrossRef]
  35. King, R.P.; Bourgeois, F. Measurement of fracture energy during single-particle fracture. Miner. Eng. 1993, 6, 353–367. [Google Scholar] [CrossRef]
Figure 1. The rock samples.
Figure 1. The rock samples.
Materials 15 07890 g001
Figure 2. Test equipment. (a) SHPB apparatus system. (b) Screening device.
Figure 2. Test equipment. (a) SHPB apparatus system. (b) Screening device.
Materials 15 07890 g002
Figure 3. The stress/damage–strain relationship curves of samples at different loading conditions: (a) 0.05 MPa; (b) 0.08 MPa; (c) 0.14 MPa; (d) 0.23 MPa.
Figure 3. The stress/damage–strain relationship curves of samples at different loading conditions: (a) 0.05 MPa; (b) 0.08 MPa; (c) 0.14 MPa; (d) 0.23 MPa.
Materials 15 07890 g003aMaterials 15 07890 g003b
Figure 4. The relationship between peak stress and impact pressure.
Figure 4. The relationship between peak stress and impact pressure.
Materials 15 07890 g004
Figure 5. The sample debris after loading: (a) photo of the typical distribution of the sample debris; (b) figure of variation of the particle size.
Figure 5. The sample debris after loading: (a) photo of the typical distribution of the sample debris; (b) figure of variation of the particle size.
Materials 15 07890 g005aMaterials 15 07890 g005b
Figure 6. The effects of different loading speeds on stress–strain of samples.
Figure 6. The effects of different loading speeds on stress–strain of samples.
Materials 15 07890 g006
Figure 7. Solution graph of logistic equation.
Figure 7. Solution graph of logistic equation.
Materials 15 07890 g007
Figure 8. The chaotic characteristics of rock damage evolution: (a) influence of loading speed on chaotic characteristics of samples; (b) similarity of damage evolution.
Figure 8. The chaotic characteristics of rock damage evolution: (a) influence of loading speed on chaotic characteristics of samples; (b) similarity of damage evolution.
Materials 15 07890 g008
Figure 9. The fractal characteristics of rock damage at different loadings: (a,c,e,g) starting time of period-doubling bifurcation state of sample; (b,d,f,h) starting time of chaotic state of sample.
Figure 9. The fractal characteristics of rock damage at different loadings: (a,c,e,g) starting time of period-doubling bifurcation state of sample; (b,d,f,h) starting time of chaotic state of sample.
Materials 15 07890 g009aMaterials 15 07890 g009bMaterials 15 07890 g009cMaterials 15 07890 g009d
Figure 10. The fractal value at different loading conditions.
Figure 10. The fractal value at different loading conditions.
Materials 15 07890 g010
Table 1. The stress value of rock sample damage evolution.
Table 1. The stress value of rock sample damage evolution.
Sample NumberImpact PressureImpact Velocity
(m/s)
Peak Strength (MPa)Stress Values Corresponding to Different Starting Times (MPa)
(MPa)DamagePeriod-Doubling BifurcationChaotic
DT-A-10.053.5281.1036.171.279.4
DT-A-20.053.4280.3132.367.577.4
DT-A-30.053.4685.1737.673.382.3
DT-B-10.086.65167.5172.0159.8165.4
DT-B-20.086.91159.9178.3145.7155.7
DT-B-30.086.67162.1175.7163.8159.9
DT-C-10.149.71237.02119.4231.7236.3
DT-C-20.149.77221.52123.1227.5225.7
DT-C-30.149.90230.33117.9237.9227.8
DT-D-10.233.61275.10203.3266.6273.4
DT-D-20.233.67287.38207.9273.7284.9
DT-D-30.233.55289.99200.3269.9288.1
Table 2. The quality percentage of each grade of sample debris.
Table 2. The quality percentage of each grade of sample debris.
Sample NumberImpact Pressure
(MPa)
mvi/%m/gAverage Particle Diameter (mm)
12345
DT-A-10.050.93%1.31%2.24%4.29%91.23%53646.48
DT-A-20.051.01%1.11%2.04%4.63%91.11%54046.45
DT-A-30.051.29%1.10%2.39%6.99%90.99%54446.77
DT-B-10.083.00%2.81%5.24%16.67%72.28%53439.16
DT-B-20.083.35%2.97%5.02%15.43%73.23%53839.45
DT-B-30.083.54%3.54%7.08%15.27%70.58%53738.27
DT-C-10.146.62%6.17%10.17%24.36%52.68%53931.04
DT-C-20.147.63%5.97%11.31%25.08%50.47%54330.13
DT-C-30.146.31%6.77%12.32%20.61%53.99%53631.31
DT-D-10.2313.10%10.41%17.31%31.65%23.53%53218.31
DT-D-20.2313.95%10.52%20.80%31.07%23.64%53518.55
DT-D-30.2314.10%11.81%17.71%34.04%22.34%54318.16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, N.; Fu, J.; Xiong, C. Studying the Characteristics of Chaos and Fractals of Construction Rocks under Different Loading Velocities. Materials 2022, 15, 7890. https://doi.org/10.3390/ma15227890

AMA Style

Wu N, Fu J, Xiong C. Studying the Characteristics of Chaos and Fractals of Construction Rocks under Different Loading Velocities. Materials. 2022; 15(22):7890. https://doi.org/10.3390/ma15227890

Chicago/Turabian Style

Wu, Nan, Jiyang Fu, and Chao Xiong. 2022. "Studying the Characteristics of Chaos and Fractals of Construction Rocks under Different Loading Velocities" Materials 15, no. 22: 7890. https://doi.org/10.3390/ma15227890

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