Study on Instability Mechanism of Extraction Structure under Undercut Space Based on Thin Plate Theory in Block Caving Method

Xia, Zhi-Yuan; Tan, Zhuo-Ying

doi:https://doi.org/10.1155/2021/5548213

Shock and Vibration

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Advanced Computational Methods for Mitigating Shock and Vibration Hazards in Deep Mines

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 5548213 | https://doi.org/10.1155/2021/5548213

Study on Instability Mechanism of Extraction Structure under Undercut Space Based on Thin Plate Theory in Block Caving Method

Zhi-Yuan Xia^1,2,3and Zhuo-Ying Tan^1,2,3

Academic Editor: Junfei Zhang

Received03 Feb 2021

Revised16 Feb 2021

Accepted22 Feb 2021

Published28 Feb 2021

Abstract

The instability of extraction structure under the undercut space in the block caving stope presents specific characteristics: rib spalling and floor heave in ore-loading roadway and collapse of major apexes. In order to study the stress and displacement evolution law of extraction structure under undercut space and reveal the instability mechanism of extraction structure, the numerical simulation model of block caving stope was established using the finite difference software FLAC^3D. According to the postundercutting strategy in Tongkuangyu Mine in China, extraction structure was formed first in the simulation process, and then the undercut level was divided into eight units for excavation step by step. The stress and displacement of extraction structure after each step of undercutting were monitored and analyzed. Based on the thin plate theory, the mechanism of stress change and deflection deformation of extraction structure was revealed. The research results show that, under the action of high horizontal tectonic stress and vertical stress, the extraction structure under undercut space produces vertical upward bending deformation after undercutting during the block caving. The tension stress concentration gradually appears in the side wall of the ore-loading roadway and the tip of the major apexes; with the increase of the undercutting area, the degree of tensile stress concentration gradually becomes strong; when the tensile strength of the rock mass in extraction structure is exceeded, extraction structure presents instability. It is necessary to make the overlying ore collapse on extraction structure as soon as possible after undercutting, which is beneficial to release the tension stress in the extraction structure under undercutting space.

1. Introduction

As a large-scale, low-cost, and high-efficiency underground mining method, block caving mining method is the preferred method for deep ore resources when conditions permit [1–3]. The extraction structure of block caving mining method refers to a series of roadways and excavations between the undercut level and the production level [4, 5]. The extraction structure undertakes the task of ore drawing, and all ores need to be transported out of the stope through the extraction structure. Therefore, ensuring the safety and stability of the extraction structure is one of the key factors for the successful application of the block caving method. Due to the particularity of block caving mining, the service life of extraction structure is long, and thus extraction structure is in a complex and changing high stress environment, which makes it difficult to maintain its stability, so that it has high risk of instability.

Many scholars have studied the theories and methods related to the stability of the extraction structure and ground pressure behavior in block caving method. Trueman et al. [6–8] used numerical simulation, theoretical analysis, on-site monitoring, and other means to study the influence of undercut strategy, tectonic stress, extraction structure height, advancing undercut front retention time, and other factors on the extraction structure of block caving method, and control measures were put forward. Rojas et al. [9, 10] studied the causes of the failure of the extraction structure in the ore caving stage by monitoring the stress, deformation, crack propagation, and microseismic events of the ore-loading roadway. Xia et al. [11, 12] studied the instability mechanism of the extraction structure under the conditions of multiple production levels in one panel and poor undercutting in block caving method by using three-dimensional numerical simulation method. The research on the force of ore bulk on the extraction structure under the undercut space has gradually attracted the attention of mining scholars. Pierce [13–15] used limit equilibrium method, empirical formula, and numerical simulation to study the force of ore bulk on the extraction structure in the ore drawing stage, so as to guide the maintenance of the extraction structure in the ore drawing process. Castro et al. [16, 17] studied the variation of vertical stress of extraction structure in stope under three kinds of ore drawing conditions by using laboratory physical model; the research shows that the vertical stress was mainly affected by the area of ore drawing area and the distance from the ore drawing advancing front.

In the actual production of block caving mine, under the condition of large high horizontal tectonic stress, with the advance of undercutting, the extraction structure under the undercut space presented instability in some areas. However, previous studies have not really revealed the root cause of stress and deformation of the extraction structure under undercut space in the block caving. Based on the engineering background of Tongkuangyu Mine, a typical block caving mine in China, the evolution characteristics of the stress and displacement of the extraction structure under the undercut space in the process undercutting of block caving were simulated by using the three-dimensional finite difference software FLAC^3D, the stress of the extraction structure under the undercut space was analyzed based on the thin plate theory, and the instability mechanism of the extraction structure under the undercut pulling space was revealed. It is of great significance to guide the stability maintenance of the extraction structure in block caving mines.

2. Instability Characteristics of Extraction Structure under Undercut Space in Tongkuangyu Mine

Postundercutting strategy is adopted in Tongkuangyu Mine; the excavation of the extraction structure precedes the blasting of undercut level. The disadvantage of the postundercutting strategy is that stress concentration is formed easily in the extraction structure in front of an advancing undercut front, which has been confirmed by many scholars [6]. The in-situ stress of Tongkuangyu Mine was measured by stress relief method; the maximum principal stress was in the horizontal direction, which was almost consistent with the strike of ore body. According to the investigation on-site, it was found that under the influence of high horizontal tectonic stress, with the increase of undercut area, the extraction structure under undercut space presented partial instability; the instability characteristics presented rib spalling and floor heave in ore-loading roadway and collapse of major apexes. As shown in Figure 1, the most easily observed features at the production level were the deformation and failure of the ore-loading roadway.

(a)

(b)

3. Stress Distribution and Evolution Law of Extraction Structure under Undercut Space

3.1. Model Construction and Simulation Process

In this study, a three-dimensional finite difference software FLAC^3D was used to establish the numerical model, and the model range was taken from the actual mining range of Tongkuangyu Mine. The model was 385 m long, 250 m wide, and 260 m high, with 989824 units in total. The overall structure of the model is shown in Figure 2(a). There are 8 ore-crossing veins and 31 ore-loading roadways and draw-bells in production level, and the burial depth of production level is 550 m. The draw-bell is 13 m long, 10 m high, and 6.4 m wide at the bottom and 11 m wide at the top. The size of ore-loading roadways and the ore drawing roadways is 3.8 m × 3.2 m (width × height). The internal structure diagram and roadway name of the model are shown in Figure 2(b).

In the simulation process, the extraction structure was formed firstly, and then undercut level was divided into 7 units for excavation step by step. The width of the undercutting unit advanced in each step was 30 m, which was the distance between the two adjacent ore-loading roadways. The stress and displacement of extraction structure under undercut space were monitored after the excavation of each unit. The excavation sequence of undercut unit is shown in Figure 2(c).

The rock mass of Tongkuangyu Mine was sampled on-site and then made into standard rock samples for laboratory test. The measured rock parameters were reduced and weakened to obtain the parameters of engineering rock mass which could be used for numerical simulation. The Mohr–Coulomb failure criterion was used in the model calculation, but it could not effectively describe the influence of joints, fracture, and structural planes on the strength of rock mass; however, Hoek–Brown criterion could accurately characterize the failure of rock mass materials with a large number of joints and fracture. The cohesion and internal friction angle in the equivalent Mohr–Coulomb failure criterion were converted by Hoek–Brown criterion, so that the Mohr–Coulomb failure criterion could be better applied to the engineering rock mass in the numerical model. The physical and mechanical parameters required for simulation are shown in Table 1. The vertical stress was applied on the top of the model and the in-situ stress was applied inside the model. The horizontal movement is limited on the side of the model, and the vertical movement was limited on the bottom of the model.

The characteristic equations of in-situ stress are shown as follows:

3.2. Distribution Characteristics of Stress and Displacement of Extraction Structure under Undercut Space

After three-dimensional numerical simulation, the characteristic cutting planes of stress and displacement contours after each excavation step were selected for comparative analysis, and the stress concentration area and large deformation area after each excavation step could be found, thus the location of the extraction structure prone to ground pressure failure could be obtained.

Evolution characteristics contours of minimum principal stress in XOZ cutting plane of extraction structure are shown in Figure 3. In the contours, the negative value represents the compressive stress, and the positive value represents the tensile stress. The red area is the high stress concentration area of the minimum principal stress, which is positive, so it represents tension stress. Figure 3(a) shows the stress state of the extraction structure when the undercut strike length is 90 m. It can be seen that the side wall of ore-loading roadways and the tip of major apexes under the undercut space are in the tension stress concentration area, and the maximum tension stress value reaches 4.82 MPa. With the continuous advancement of undercutting, the strike length of undercut gradually reaches 210 m, and the stress state of the extraction structure is shown in Figure 3(b) at this time. It can be seen that the concentration of tension stress in the side wall of ore-loading roadways and the tip of major apexes under the undercut space increases, and the maximum tensile stress value reaches 5.07 MPa. Therefore, with the increase of undercut area, it is easy to reach the tension failure condition at the side wall of ore-loading roadways and the tip of major apexes under the undercut space, which leads to rib spalling of ore-loading roadways and the collapse of the tip of major apexes.

(a)

(b)

Evolution characteristics contours of vertical displacement in XOZ cutting plane of extraction structure are shown in Figure 4. In the contours, the negative value represents the vertical downward displacement, and the positive value represents the vertical upward displacement. The red position is the area with large vertical displacement. Figure 4(a) shows the vertical displacement state of the extraction structure when the undercut strike length is 90 m. It can be seen that the whole extraction structure under the undercut space has a vertical upward displacement, and the vertical upward displacement of major apexes tip in the middle of the stope is the largest, reaching 98 mm. With the continuous advancement of undercutting, the strike length of undercut gradually reaches 210 m, and the vertical displacement state of the extraction structure is shown in Figure 4(b) at this time. It can be seen that the vertical upward displacement of the side wall of the ore-loading roadways and the tip of major apexes under the undercut space increases, and the maximum vertical displacement reaches 116 mm. Therefore, with the increase of the undercut area, the vertical displacement of the extraction structure under the undercut space increases.

(a)

(b)

3.3. Evolution Law of Stress and Displacement of Extraction Structure under Undercut Space

The stress and displacement monitoring points were arranged at the side wall of the ore-loading roadways and the tip of major apexes under the undercut space. The variation law of the maximum tensile stress and maximum vertical displacement with the undercut strike length is shown in Figure 5. It can be seen from Figure 5(a) that with the increase of the advancing strike length of the undercut, the tension stress concentration of the major apexes under the undercut space becomes intense, the tension stress of the side wall of the ore-loading roadways increases in a oscillatory manner, and they are all in a state of high tension stress concentration. It can be seen from Figure 5(b) that the maximum vertical displacement of the major apexes and the side wall of the ore-loading roadways under the undercut space increases with the increasing of the advancing strike length of the undercut.

(a)

(b)

4. Analysis of Instability Mechanism of Extraction Structure under Undercut Space Based on Thin Plate Theory

4.1. Establishment of Equivalent Model of Extraction Structure Thin Plate

In order to reveal the mechanism of the bending deformation and tension stress increase of the extraction structure under the undercut space, based on the thin plate theory of elasticity, an equivalent model of the extraction structure under the lateral and longitudinal compression of Winkler elastic support was proposed in this section. The extraction structure rock mass under the undercut space was equivalent to a rectangular plate with reduced elastic modulus, which was simply supported on four sides, as shown in Figure 6. Due to the bottom of the thin plate equivalent model is connected with the rock mass floor, it can be assumed that the displacement at any point on the floor surface is proportional to the stress at that point, that is, the rock mass floor is composed of a series of spring elements. The lateral compressive stress is perpendicular to the plate surface, and the longitudinal pressure is the in-plane pressure of the thin plate; there are the two horizontal principal stresses and . The x-axis direction of the model is the inclination of extraction structure in the stope, and the length is set as .The direction of y-axis is the strike of extraction structure in the stope, and the length is set as . The direction perpendicular to the extraction structure of stope is z-axis, and the thickness of thin plate is .

Because the height of the extraction structure is about 16 m, the minimum width of the undercut range of the stope in Tongkuangyu Mine is 80 m. According to the definition of the thin plate, the ratio of the thickness of the thin plate to the minimum size of the middle surface is in the range of 1/80∼1/5, which is the thin plate. It is obvious that the model is in the range of the thin plate definition. When the thin plate is bent, the middle surface becomes a curved surface, and the z-direction displacement of each point on the middle surface is called deflection (). When , it is called small deflection bending. According to the field investigation results of extraction structure instability, the extraction structure deformation meets the condition of small deflection bending.

The extraction structure rock mass under the undercut space is approximately equivalent to the elastic thin plate after the reduction of elastic modulus. Due to the extraction structure thin plate is an incomplete thin plate, it would cause incompatible deformation and stress inside. However, the causes and influencing factors of the extraction structure bending deformation and tension stress are analyzed in this section only, so the extraction structure elastic thin plate is equivalent to a homogeneous thin plate. The elastic modulus of the equivalent elastic thin plate should be reduced according to the excavation rate of the roadway. According to statistics, the excavation rate of the roadway in the extraction structure rock space of Tongkuangyu Mine is about 40%, thus the elastic modulus of the equivalent elastic thin plate can be calculated according to the following formula:

In formula (2), is the equivalent elastic modulus of the extraction structural thin plate, is the elastic modulus of complete rock mass, and is the excavation rate of extraction structure rock mass.

4.2. Solution of Deflection and Tension Stress of Extraction Structure under Undercut Space

According to the static equilibrium condition of the equivalent thin plate differential body, the small deflection bending differential equation can be derived as follows:where is the bending stiffness of the extraction structure thin plate, is the deflection of the extraction structural plate, is the Laplacian operator, , , and are the internal forces of thin plate, is the initial lateral stress of thin plate model of extraction structure, is the foundation modulus of floor under extraction structure, is the thickness of extraction structure thin plate, and is Poisson’s ratio of thin rock mass of extraction structure.

Due to the thin plate of the extraction structure is approximately regarded as satisfying the boundary condition of simply supported four edges, the deflection of the extraction structure can be calculated by the Navier method, and the deflection expression is assumed to be a double trigonometric series:

According to the longitudinal load of the thin plate, the internal force of the middle plane is synthesized as follows:

By combining formula (6) with formula (3), the following formula can be given:

In order to obtain in the deflection expression, we need to expand in the form of double trigonometric series:

In order to find a suitable coefficient of this trigonometric series, multiply on both sides of equation (8) to find the integral from 0 to b, when , , and when , .

Therefore, the following formula can be given:

Multiply on both sides of formula (8) and find the integral from 0 to b,

When and are odd numbers, a suitable coefficient is obtained:

Therefore, the double trigonometric series expansion of can be obtained:

By combining formula (12), deflection expression (5), and formula (7), the following formula can be given:

In order to simplify the calculation, according to the actual deflection convergence degree of the extraction structure under the undercut space, when and are taken as 1, the accuracy requirement of deflection can be basically met. Therefore, formula (13) can be written as

Due to cannot be 0, the following formula can be given:

From formula (15), it can be obtained that

The deflection formula of the equivalent thin plate model of the extraction structure can be obtained:

According to the physical equation of thin plate, three main stress components are expressed as deflection:

By combining formula (17) and deflection expression (18), the stress of each point of the extraction structure thin plate model can be obtained:

When and , the deflection of the extraction structure thin plate model is the largest:

When , , and , at this point, the bending moment of the thin plate is the largest, and the tension stress is the largest at the upper part of the thin plate:

As the strike length of stope is greater than the dip length of stope , according to equation (21), it can be concluded that , thus is the maximum tensile stress of thin plate.

4.3. Variation Law of Deflection and Tension Stress of Extraction Structure under Undercut Space

In order to study the variation law of the maximum deflection and maximum tension stress of the equivalent model of the extraction structure thin plate with the advancing of the undercutting and the initial lateral stress, according to the geological data of Tongkuangyu Mine, the calculation parameters in the formula and are assigned; the calculation parameters are shown in Table 2.

According to the actual process of extraction drawing in Tongkuangyu Mine, when the horizontal inclined length of undercut reached a certain value, it remained unchanged, and the horizontal strike length of undercut increased continuously until the end of ore drawing in the middle section. Therefore, the variation of the maximum deflection and the maximum tension stress of the thin plate equivalent model with the increase of the undercut area can be simplified as the variation with the increase of the horizontal strike length of the undercut. After assignment, the variation rules of the maximum tension stress and the maximum deflection with the increase of the horizontal strike length under different initial stresses are shown in Figure 7.

(a)

(b)

It can be seen from Figure 7(a) that the maximum tension stress of thin plate increases with the increase of strike length, but the growth rate slows down and finally tends to be stable. The decrease of initial lateral stress weakens the maximum tension stress concentration but does not affect the variation trend of maximum tension stress with strike length. Therefore, with the increase of the undercut area, the maximum tension stress of the extraction structure rock mass under the undercut space would increase. When the tensile strength of the rock mass is exceeded, extraction structure presents instability. To make overlying ore cave on the extraction structure as soon as possible, the initial lateral stress in the vertical direction can be reduced, thus the tension stress of the extraction structure below the undercut space can be released to a certain extent, which can be conducive to the stability of the extraction structure. It can be seen from Figure 7(b) that the maximum vertical displacement of thin plate increases with the increase of strike length, but the growth rate slows down and finally tends to be stable. The decrease of initial lateral stress weakens the maximum vertical displacement, but does not affect the variation trend of maximum vertical displacement with strike length. Therefore, with the increase of the undercut area, the maximum vertical displacement of the extraction structure rock mass under the undercut space would increase. To make overlying ore cave on the extraction structure as soon as possible, the vertical displacement of the extraction structure below the undercut space can be released to a certain extent.

5. Analysis of Instability Mechanism of Extraction Structure under Undercut Space

The extraction structure under the undercut space is similar to a weak “thin rock stratum plate” at the bottom of the stope. Under the influence of high horizontal tectonic stress and vertical stress, the excavation space is squeezed, and the extraction structure would produce upward bending deformation. In order to resist the deflection deformation, the extraction structure gradually presents tension state. With the increasing of the undercut space, the concentration degree of tension stress becomes higher and higher. When the tension strength of the rock mass is exceeded, extraction structure presents instability. The tension stress is easy to concentrate on the side wall of the ore-loading roadways and the tip of major apexes under the undercut space. Therefore, with the increase of the undercut area and the insufficient support strength of the extraction structure, the ore-loading roadways would show spalling, and the major apexes would gradually collapse. Due to the whole extraction structure would have vertical upward bending deformation, if the mine only strengthens the support of side wall and roof and does not pay attention to the maintenance of the floor of the production roadways, the floor of the production roadway would heave. At the same time, inducing the overlying strata to collapse on the extraction structure and carrying out orderly and balanced ore drawing can reduce the vertical deflection and tension stress concentration of the extraction structure under the undercut space, so as to ensure the stability of the extraction structure. The combined support form of bolt mesh cable shotcreting and floor concrete reverse arch could effectively control the loose deformation of the surrounding rock of the ore drawing roadway and improve the overall strength of the extraction structure and increase the stability of the extraction structure.

6. Conclusions

The evolution characteristics of the stress and displacement of the extraction structure under the undercut space in the process undercutting of block caving were simulated using FLAC^3D based on the engineering background of Tongkuangyu Mine, the stress of the extraction structure under the undercut space was analyzed based on the thin plate theory, and the instability mechanism of the extraction structure under the undercut pulling space was revealed. The conclusions are as follows:(1)Under the action of high horizontal tectonic stress and vertical stress, the extraction structure under undercut space produces vertical upward bending deformation after undercutting during the block caving;(2)The tension stress concentration gradually appears in the side wall of the ore-loading roadway and the tip of the major apexes; with the increase of the undercutting area, the degree of tensile stress concentration gradually becomes strong; when the tensile strength of the rock mass in extraction structure is exceeded, extraction structure presents instability;(3)It is necessary to make the overlying ore collapse on extraction structure as soon as possible after undercutting, which is beneficial to release the tension stress in the extraction structure under undercutting space.

Data Availability

The 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 have no conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 51574015).

References

J. P. Hurtado, N. Díaz, E. I. Acuña, and J. Fernández, “Shock losses characterization of ventilation circuits for block caving production levels,” Tunnelling and Underground Space Technology, vol. 41, p. 88, 2014.
View at: Publisher Site | Google Scholar
K.-S. woo, E. Eberhardt, D. Elmo, and D. Stead, “Empirical investigation and characterization of surface subsidence related to block cave mining,” International Journal of Rock Mechanics and Mining Sciences, vol. 61, p. 31, 2013.
View at: Publisher Site | Google Scholar
R. Trueman, R. Castro, and A. Halim, “Study of multiple draw-zone interaction in block caving mines by means of a large 3D physical model,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 7, pp. 1044–1051, 2008.
View at: Publisher Site | Google Scholar
V. Sánchez, R. L. Castro, and S. Palma, “Gravity flow characterization of fine granular material for Block Caving,” International Journal of Rock Mechanics and Mining Sciences, vol. 114, p. 24, 2019.
View at: Publisher Site | Google Scholar
R. Rafiee, M. Ataei, R. KhalooKakaie, S. E. Jalali, F. Sereshki, and M. Noroozi, “Numerical modeling of influence parameters in cavabililty of rock mass in block caving mines,” International Journal of Rock Mechanics and Mining Sciences, vol. 105, p. 22, 2018.
View at: Publisher Site | Google Scholar
R. Trueman, M. Pierce, and R. Wattimena, “Quantifying stresses and support requirements in the undercut and production level drifts of block and panel caving mines,” International Journal of Rock Mechanics and Mining Sciences, vol. 39, no. 5, p. 617, 2002.
View at: Publisher Site | Google Scholar
R. J. Butcher, “Design rules for avoiding draw horizon damage in deep level block caves,” Journal of the Southern African Institute of Mining and Metallurgy, vol. 99, no. 2, pp. 151–155, 1999.
View at: Google Scholar
P. J. Bartlett, Planning, Implementation, Operation and Monitoring of a Cave Mining Method with Coarse Fragmentation with Reference to Cave Mining at Premier, University of Pretoria, Pretoria, South Africa, 1998.
E. Rojas, R. Molina, A. Bonani et al., “The pre-undercut caving method at the el teniente mine, codelco Chile,” in Massmin 2000 proceedings, pp. 261–266, Brisbane, Australia, November 2000.
View at: Google Scholar
G. Viegas, K. Bosman, D. Angus et al., “Mapping cave front growth utilising the collective behaviour of seismicity and velocity fields,” in Proceedings of the Caving 2018, Fourth International Symposium on Block and Sublevel Caving, pp. 577–588, Vancouver, Canada, October 2018.
View at: Google Scholar
Z. Xia, Z. Tan, and Y. Miao, “Damage evolution mechanism of extraction structure during mining gently dipped orebody by block caving method,” Geotechnical and Geological Engineering, vol. 38, no. 4, 2020.
View at: Publisher Site | Google Scholar
Z. Xia, Z. Tan, Q. Pei, and J. Wang, “Ground pressure damage evolution mechanism of extraction level excavations induced by poor undercutting in block caving method,” Geotechnical and Geological Engineering, vol. 37, no. 5, pp. 4461–4472, 2019.
View at: Publisher Site | Google Scholar
M. E. Pierce, “Forecasting vulnerability of deep extraction level excavations to draw-induced cave loads,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 11, no. 3, pp. 527–534, 2019.
View at: Publisher Site | Google Scholar
H. Sahupala, C. Brannon, S. Annavarapu et al., “Recovery of extraction pillars in the deep ore zone (DOZ) block cave, PT freeport Indonesia,” in Proceedings of the Massmin 2008, 5th International Conference & Exhibition on Mass Mining, pp. 191–202, Lulea, Sweden, June 2008.
View at: Google Scholar
M. Orellana, C. Cifuentes, and J. Díaz, “Caving experiences in esmeralda sector, el teniente mine,” in Proceedings of the 3rd International Symposium On Block And Sublevel Caving, pp. 78–90, Santiago, Chile, June 2014.
View at: Google Scholar
R. Castro, R. Gomez, M. Pierce et al., “Experimental quantification of vertical stresses during gravity flow in block caving,” International Journal of Rock Mechanics and Mining Sciences, vol. 127, 2020.
View at: Publisher Site | Google Scholar
R. Castro, Study of the Mechanisms of Gravity Flow for Block Caving, The University of Queensland, Brisbane, Australia, 2006.

Copyright

Copyright © 2021 Zhi-Yuan Xia and Zhuo-Ying Tan. 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.

PDF Download Citation

Download other formats

Order printed copies