Shock and Vibration
Volume 2021 (2021), Article ID 9998558, 15 pages
https://doi.org/10.1155/2021/9998558

Research Article

The Dynamic Response and Failure Model of Thin Plate Rock Mass under Impact Load

Feng Li,1,2 Xinhui Dong,² Yue Wang,³ Hanwu Liu,² Chuang Chen,² and Xin Zhao²

¹State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China
²School of Emergency Management and Safety Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China
³School of Safety Science and Engineering, Xinjiang Institute of Engineering, Xinjiang, Urumqi 830023, China

Correspondence should be addressed to Feng Li; lifengcumtb@126.com

Received 10 March 2021; Revised 25 March 2021; Accepted 31 March 2021; Published 20 April 2021

Academic Editor: Gong-Da Wang

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

Abstract.

The layered rock mass widely exists in mining, construction, transportation, and water conservancy projects, and the damage phenomena of plate crack and spalling often occurs in the process of coal and rock dynamic disaster in deep mining. Therefore, the rock mass nearby excavation surface is usually considered to be composed of layers of thin plate rock mass to reveal the damage and failure mechanism of rock mass. In the whole dynamic process of mining and coal and rock dynamic disaster, rock mass would bear the dynamic disturbance from mine earthquake, and at present, the mechanical characteristics of rock mass are mainly studied under static load, while dynamic mechanical response characteristics and the mechanisms of dynamic damage, failure, and disaster-causing are still unclear. This study mainly focused on the dynamic response characteristic and failure mechanism of rock mass based on a rectangular thin plate model. The frequency equations and deflection equations of the thin plate rock mass with different boundary conditions (S-F-S-F, S-C-S-C, and C-C-C-C) were established under free vibration by the thin plate model and the dual equation of the Hamilton system, and the deflection equations under impact load were derived based on the Duhamel integral. And then, the effective vibration modes of the thin plate rock mass with different boundary conditions and their natural frequencies were obtained by Newton’s iterative method. Based on the third-strength theory and the numerical simulation results by LS-DYNA, the maximum shear of the effective vibration modes and the processes of damage and failure under impact load were analyzed. The research results showed that the initial position of damage and failure may be determined by effective vibration mode with the lowest frequency; the develop tendency of which by the combined actions of other effective vibration modes and the effective vibration modes with lower frequency could have greater influence on the process of damage and failure of the thin plate rock mass, which are beneficial to revealing the mechanism of coal and rock dynamic disaster.

1. Introduction

Thin plate theory has been widely used and achieved good results in the research of stope movement and the failure mechanism of stope roof. The sedimentary process and the arrangement of mineral particles during the formation of coal measure rock mass make it have obvious bedding structure and show layered structure macroscopically [1–4]. This layered rock mass widely exists in mining, construction, transportation, and water conservancy projects, such as slate [5], shale [6], sandstone [7], and phyllite [8]. During the excavation of underground engineering, a series of layered rock mass failure is a key factor for major geological disasters. For example, the water inrush accident in Chaoyang Tunnel of Luo Tuoshan and Gui Nan High-Speed Railway and the surface collapse in Ma Jiliang Mine of Datong were all caused by the instability of layered rock structure. A group of dominant parallel bedding exists in layered rock mass generally, which are usually regarded as transversely isotropic thin plates in mechanics. Therefore, many scholars have introduced the thin plate theory into the calculation of stope movement and put forward the theory of key strata [9] and stope thin plate pressure [10], the fracture law of basic roof strata structure [11–13], and the limit analysis method of strata movement [14]. Based on the thin plate theory, the “O-X” and “K” fracture patterns of stope roof were deduced in overlayer strata movement [15–20], and the specific location of crack bifurcation and formation process of spatial “X” crack pattern were obtained [21, 22]. Currently, this theory is widely used in theoretical research, experimental research, and numerical simulation to reveal the laws of plate rock failure and stope movement.

Thin plate rock mass model has been used widely to reveal the mechanism and process of coal and rock dynamic disaster in deep mining where coal and gas outburst, rockburst, rockfall, and chamber instability are prone to occur and usually accompanied by the damage phenomena of plate crack and spalling in stress concentration areas [23–28] (Figure 1).

Figure 1: The damage phenomena of plate crack in rock mass.

On the one hand, the excavation of coal and rock mass in heading face of deep mine would lead to stress redistribution in front of heading face, and the distribution of stress concentration state will present layered phenomenon. And the tangential stress increases and the radial stress decreases in the surrounding rock of the heading face. Under the loading, the cracks in the rock expand along the maximum principal stress, which can form a plate failure surface roughly parallel to the heading face [29]. On the other hand, the elastic compressive strain energy stored in the coal and rock mass with high stress in deep area would be released and reflected on the excavation surface to form tensile stress wave because of the unloading of dynamic load caused by excavation. When the tensile stress exceeds the tensile strength of coal and rock mass, the failure surface with plate crack will be formed and roughly parallel to the excavation surface. And the rapid propagation of crack will split the surrounding rock into several rock plates and form thin plate rock finally [30]. Currently, based on the thin plate model, the mechanical behaviors of plate cracking and spalling failure of rock mass were studied deeply and widely [31–34].

In the whole dynamic process of mining, and during the beginning, formation, development, and termination of coal and rock dynamic disaster, rock mass will bear the dynamic disturbance from mine earthquake, such as blasting, roof and floor breaking, and rock instability. At present, the mechanical characteristics of rock mass are mainly studied under static load, while the research studies on dynamic mechanical response characteristics are only in the stage of experimental exploration, and the mechanisms of dynamic damage, failure, and disaster-causing are still unclear [35–40], because the vibration of rock mass under external force has a significant impact on its mechanical properties [35, 37, 41].

This study mainly focused on the dynamic response characteristic and failure mechanism of rock mass based on a rectangular thin plate model. The Hamilton dual system was introduced to research the dynamic response characteristics of thin plate rock mass under external forced vibration. Besides, according to the characteristics of uniform impact load, the effective vibration mode solution was studied by using Duhamel integral principle. And the process and mechanism of dynamic damage and failure of thin plate rock mass were analyzed, which are beneficial to revealing the mechanism of coal and rock dynamic disaster.

2. Free Vibration Model of the Thin Plate Rock Mass

2.1. Control Equation of Free Vibration of the Thin Plate Rock Mass

In this study, the assumptions were as follows: the rectangular plates of thin rock mass are transversely isotropic in mechanics; the density of the thin plate was ρ; the thickness was h; the size was a × b; the elastic modulus was E; and Poisson’s ratio was υ. Based on the transient equilibrium conditions of the internal mechanics of the thin plate, the differential equation of thin plate rock mass can be derived under forced vibration [42]:

where (x, y, t) is the deflection, D is the bending stiffness, and q(x, y, t) is external dynamic load.

For the purpose to solve the homogeneous equation of equation (1), we set q(x, y, t) = 0 and insert it into (1), and the free vibration differential equation of the thin plate could be obtained:

The vibration of the thin plate is the harmonic oscillator with time. Therefore, (x, y, t) could be expressed as follows:

where ω is the natural frequency and W(x, y) is the deflection mode function.

With (2),

where k⁴ = ρhω²/D.

2.2. Hamilton Dual Vibration Equation of the Thin Plate Rock Mass

In order to decouple the physical quantities in equation (4), it was necessary to introduce the Hamilton dual equations. If we set , according to the mechanical analysis results of the thin plate, the relationship between the physical quantities may be as follows [43, 44]:

where Mx, My, Mxy, Qx, Qy, Vx, and Vy are the bending moment, torque, transverse shear force, and transverse total shear force, respectively. If we set vector  = [W, θ, −Vx, Mx], then the dual equation of the Hamilton system could be expressed as follows [45]:

Then,

2.3. Solution of the Hamilton Dual Equations

Based on the symplectic geometry method, the solutions of (6) could be obtained by the separation variable method [43]. And if W(x) and W(y) were the deflection modes along the x and y directions, respectively, their specific forms may be as follows [46, 47]:

where α₁ and α₂ are characteristic values in the x direction; β₁ and β₂ are characteristic values in the y direction. The parameters α₁, α₂, β₁, and β₂ meet the following rules:

, , , .

3. Deflection Equation of the Thin Plate Rock Mass under Free Vibration

It could be considered that the rock mass nearby excavation surface is composed of layers of thin plate rock mass in mining. In order to deduce the deflection equation of thin plate rock mass under free vibration, its boundary conditions need to be introduced, which is fixed (C) on four sides before failure, with the failure of rock mass near boundary caused by external dynamic load, these thin plate rock mass will gradually become simply supported (S) or free (F) state. Therefore, it is necessary to analyze the deflection equations of thin plate rock mass under different boundary conditions.

3.1. Deflection Equation of the Thin Plate Rock Mass S-F-S-F under Free Vibration

The deflection equation of the thin plate rock mass with both sides simply supported (S-S) along the x direction could be expressed as follows [46–48]:

The boundary conditions of the thin plate rock mass with both free sides (F-F) along the y direction may be as follows:

With (8) and (9),

Then, we insert (9) into (12) and (13):

where j = −α₁−υβ₁², f = α₂−υβ₁², m = −α₁³−(2−υ) β₁²α₁), n = α₂³−(2−υ) β₁²α².

Because the above formula has nonzero solutions, the determinant of the coefficient matrix on the left was zero. Therefore, the frequency equation along the y direction could be obtained:

We set C₂ = 1and insert it into (14):

A2 = −f/j, B2 = −nk2/m, C2 = 1, D2 = k3.

Then,

3.2. Deflection Equation of the Thin Plate Rock Mass C-C-C-C under Free Vibration

As for the thin plate rock mass with four-sided fixed supported (C-C-C-C), the boundary conditions C-C along the y direction may be as follows:

W(x, 0) = 0, ∂W(x, 0)/∂y = 0; W(x, b) = 0, ∂W(x, b)/∂y = 0.

With (9), the frequency equation along the y direction could be obtained similarly:

And the deflection equation could be expressed as follows:

where

, .

3.3. Deflection Equation of the Thin Plate Rock Mass S-C-S-C under Free Vibration

As for the thin plate rock mass with two sides fixed and other two sides simply supported (S-C-S-C), the boundary conditions S-S along the x direction may be as follows:

W(0, y) = 0, ∂²W(0, y)/∂x² = 0; W(a, y) = 0, ∂²W(a, y)/∂x² = 0.

With (8), the frequency equation could be derived similarly:

And the deflection equation could be expressed as follows:

Deflection and frequency equations of the thin plate rock mass under different boundary conditions are shown in Table 1.

Table 1: Deflection and frequency equations of the thin plate rock mass under different boundary conditions.

Thin plate rock mass
Frequency equations
Deflection equations W(x, y)

x direction:
x direction: W(x) = sin β₁x

y direction:
y direction:

j = −α₁−υβ₁², f = α₂−υβ₁², m = −α₁³−(2−υ) β₁²α₁), n = α₂³−(2−υ) β₁²α2,

x direction:
x direction:

y direction:
y direction:

,

x direction:
x direction: W(x) = sin β₁x

y direction:
y direction:

4. Deflection Equation of the Thin Plate Rock Mass under Forced Vibration

4.1. Solution of the Control Equation

The solutions’ form of the nonhomogeneous control equation (1) can be expressed as follows:

We insert (21) into (1):

With (4),

Then,

The orthogonality of the deflection equations may be as follows [48]:

We multiply both sides of (24) by Wm(x, y) and do integral over the thin plate:

We set

Then,

Based on the Duhamel integral, the solutions can be expressed as follows:

If the external dynamic load q(x, y, t) is assumed to be impact load, it may be as follows:

Then,

Therefore, the solutions of control equation (see (1)) of the thin plate rock mass under impact load can be expressed as follows:

where .

5. Effective Vibration Modes of the Thin Plate Rock Mass

Taking granite in a mine in Sichuan, the parameter values of the thin plate rock mass were as follows: h = 0.06 m, ρ = 2800 kg/m³, E = 72 × 10⁹ Pa, υ = 0.3, a × b = 3 m × 3.6 m.

5.1. Effective Vibration Modes of the Thin Plate Rock Mass S-F-S-F

Based on the frequency equation of the thin plate rock mass S-F-S-F, the Newton iteration method was introduced to solve the effective vibration modes. The calculated results of the first 10 order vibration parameters under free vibration are shown in Table 2, and the calculated results under impact load are shown in Table 3. According to Table 3, the 1st order, 3rd order, 8th order, and 10th order were effective vibration modes, and the natural frequencies of which were 101 rad/s, 300 rad/s, 909 rad/s, and 989 rad/s, respectively, and we could obtain that the coefficient An/q₀ of the 3rd order, 8th order, and 10th order was about 6.4%, 4.5%, and 0.5% of the 1st order.

Table 2: The first 10 order vibration parameters of the thin plate rock mass under free vibration (S-F-S-F).

Parameter
1st order
2nd order
3rd order
4th order
5th order
6th order
7th order
8th order
9th order
10th order

β₁
1.047
1.047
1.047
2.094
2.094
1.047
2.094
3.142
2.094
1.047

α₁
≈0
0.716
1.469
≈0
0.75
2.268
1.578
≈0
2.382
3.105

α₂
1.481
1.645
2.086
2.961
3.056
2.709
3.356
4.443
3.801
3.44

K
1.047
1.269
1.804
2.094
2.225
2.498
2.622
3.142
3.172
3.277

Table 3: The first 10 order vibration parameter values of the thin plate rock mass under impact load (S-F-S-F).

Order
Parameter

Mn

An/q₀

1storder
101
2.798 × 10⁴
38.097
1.348 × 10⁻⁵

2nd order
148
4.536 .510³
−4.642 .610⁻⁶
−6.914 .910⁻¹²

3rdorder
300
1.736 × 10³
−0.451
−8.663 × 10⁻⁷

4th order
404
2.197 .110⁴
9.059 .010⁻¹⁰
1.020 .010⁻¹⁶

5th order
456
8.646 .610³
9.876 .810⁻²⁰
2.505 .510⁻²⁶

6th order
575
1.208 .210³
2.881 .810⁻⁶
4.148 .110⁻¹²

7th order
633
3.684 .610³
−9.813 .810⁻¹⁵
−4.208 .210⁻²¹

8thorder
909
1.786 × 10⁴
9.968
6.140 × 10⁻⁷

9th order
926
2.082 .010³
6.097 .010⁻²¹
3.163 .110⁻²⁷

10thorder
989
1.052 × 10³
−0.071
−6.797 × 10⁻⁸

5.2. Effective Vibration Modes of the Thin Plate Rock Mass C-C-C-C

Based on the frequency equation of the thin plate rock mass with four-side fixed support C-C-C-C, the effective vibration modes were also solved by Newton’s iterative method. The calculated results of the first 10 order vibration parameters under free vibration are shown in Table 4. Similarly, from the analysis results of the first 10 order vibration parameter of the thin plate rock mass C-C-C-C under impact load, we could obtain that 1st order, 5th order, and 6th order were effective vibration modes, and the natural frequencies of which were 308 rad/s, 975 rad/s, and 1309 rad/s, respectively (Table 5). And the results showed that the coefficient An/q₀ of the 5th order and 6th order were about 9.4% and 4.9% of the 1st order.

Table 4: The first 10 vibration parameters of the thin plate rock mass under free vibration (C-C-C-C).

Parameter
1st order
2nd order
3rd order
4th order
5th order
6th order
7th order
8th order
9th order
10th order

β₁
1.443
1.301
2.573
2.483
1.227
3.643
2.404
3.589
2.346
2.346

β₂
2.145
3.247
2.949
3.765
4.437
3.895
4.802
4.514
5.930
5.930

α₁
1.122
2.104
1.019
2.001
3.015
0.975
2.940
1.935
3.904
3.851

α₂
2.329
2.795
3.778
4.042
3.479
5.243
4.494
5.433
4.248
5.083

k
1.828
2.474
2.767
3.189
3.255
3.771
3.797
4.078
4.079
4.509

Table 5: The effective vibration modes of the thin plate rock mass under impact load (C-C-C-C).

Order
Parameter

Mn

An/q₀

1st order
308
5.787 .710³
13.589
7.624 .610⁻⁶

5th order
975
1.267 .210⁴
8.503
6.883 .810⁻⁷

6th order
1309
2.587 .510⁴
12.565
3.71 .710⁻⁷

5.3. Effective Vibration Modes of the Thin Plate Rock Mass S-C-S-C

The calculated results of the first 10 order vibration parameters of the thin plate rock mass S-C-S-C under free vibration are shown in Table 6. Similarly, according to the first 10 order vibration parameter values of the thin plate rock mass S-C-S-C under impact load, it could be obtained that the 1st, 5th, 6th, and 10th order were effective vibration modes, and the natural frequencies of which were 230 rad/s, 943 rad/s, 999 rad/s, and 1680 rad/s, respectively (Table 7). And the results showed that the coefficient An/q₀ of the 5th order, 6th order, and 10th order were about 10.6%, 2.8%, and 1.3% of the 1st order.

Table 6: The first 10 order vibration parameters of the thin plate rock mass under free vibration (S-C-S-C).

Parameter
1
2
3
4
5
6
7
8
9
10

β₁
1.047
2.094
1.047
2.094
1.047
3.142
2.094
3.142
1.0472
3.142

α₁
1.183
1.052
2.127
2.031
3.025
0.992
2.960
1.959
3.908
2.896

α₂
1.896
3.143
2.592
3.591
3.368
4.552
4.187
4.855
4.180
5.303

K
1.580
2.343
2.371
2.917
3.201
3.294
3.626
3.702
4.046
4.273

Table 7: The effective vibration modes of the thin plate rock mass under impact load (S-C-S-C).

Order
Parameter

Mn

An/q₀

1st order
230
1.407 .410³
7.171
2.511 .510⁻⁵

5th order
943
9.980 .910²
2.540
2.699 9610⁻⁶

6th order
999
8.898 .810³
6.171
6.943 .910⁻⁷

10th order
1680
1.856 .810³
1.045
3.352 .310⁻⁷

6. Dynamic Damage and Failure of the Thin Plate Rock Mass

According to the deflection equations of the thin plate rock mass under different boundary conditions shown in Table 1, it could be obtained that the solutions’ form of the thin plate rock mass S-C-S-C was the simplest and that of the thin plate rock mass S-F-S-F was the most complex. And then, we would analyze the process of dynamic damage and failure of thin plate rock mass in order of the complexity of their solutions’ form.

6.1. Dynamic Damage and Failure Process of the Thin Plate Rock Mass S-C-S-C

Based on the third-strength theory, the maximum shear of effective vibration modes of the thin plate rock mass S-C-S-C under impact load is shown in Figures 2–5. It can be seen from Figure 2 that the maximum shear of the 1st order was mainly distributed in the middle of the two short sides and four concentration regions of the two long sides, that of the 5th order in the middle of the two short sides and the region along the long central axis of the thin plate (Figure 3), that of the 6th order in the region along the short center axis of the thin plate (Figure 4), and that of the 10th order in any region (Figure 5).

Figure 2: The maximum shear of the 1st order mode (S-C-S-C).

Figure 3: The maximum shear of the 5th order mode (S-C-S-C).

Figure 4: The maximum shear of the 6th order mode (S-C-S-C).

Figure 5: The maximum shear of the 10th order mode (S-C-S-C).

Based on LS-DYNA, the PLASTIC_KINEMATIC material model was used to simulate and analyze the process of damage and failure of the thin plate rock mass S-C-S-C under impact load. The long side of the model is horizontal, and the short side is vertical. The model size is 3 m × 3.6 m, the Cowper–Symonds model is used to describe the yield process of coal and rock, and the impact load acts uniformly on the surface of coal and rock mass. The results could be shown as follows (Figure 6): firstly, plate cracks occurred in the middle of the two short sides and the four concentration regions of the two long sides (Figures 6(a) and 6(b)), which was consistent with Figure 2, then the develop tendency of plate cracks was along the long and short central axes toward the center of the thin plate rock mass (Figures 6(c) and 6(d)), which was consistent with Figures 3 and 4, and last, plate cracks would be distributed in any region (Figure 6(d)), which was consistent with Figure 5.

(a)

(b)

(c)

(d)

Figure 6: Dynamic damage and failure of thin plate rock mass S-C-S-C under impact load.

6.2. Dynamic Damage and Failure Process of the Thin Plate Rock Mass C-C-C-C

Based on the third-strength theory, the maximum shear of effective vibration modes of the thin plate rock mass C-C-C-C under impact load is shown in Figures 7–9. It can be seen from Figure 7 that the maximum shear of the 1st order was mainly distributed in the middle of the four sides, that of the 5th order in the middle of the two short sides and the region along the long central axis of the thin plate (Figure 8), and that of the 6th order in the middle of the two long sides and the region along the short central axis of the thin plate (Figure 9).

Figure 7: The maximum shear of 1st order mode (C-C-C-C).

Figure 8: The maximum shear of the 5th order mode (C-C-C-C).

Figure 9: The maximum shear of the 6th order mode (C-C-C-C).

Similarly, based on LS-DYNA, the PLASTIC_KINEMATIC material model was used to simulate and analyze the process of damage and failure of the thin plate rock mass C-C-C-C under impact load, which could be shown as follows (Figure 10): firstly, plate cracks occurred in the middle of the two short sides and two long sides and these cracks would propagate rapidly along the boundary of the thin plate rock mass(Figures 10(a) and 10(b)), which was consistent with Figure 7, then plate cracks occurred at the central of the thin plate and the main develop tendency of these cracks was outward along the long central axis (Figures 10(c) and 10(d)), which was consistent with Figure 8, and these cracks also tended to expand outward along the direction of short side (Figure 10(d)), which was consistent with Figure 9.

(a)

(b)

(c)

(d)

Figure 10: Dynamic damage and failure of thin plate rock mass C-C-C-C under impact load.

6.3. Dynamic Damage and Failure Process of the Thin Plate Rock Mass S-F-S-F

Similarly, based on the third-strength theory, the maximum shear of effective vibration modes of the thin plate rock mass S-F-S-F under impact load is shown in Figures 11–14. It can be seen from Figure 11 that the maximum shear of the 1st order was mainly distributed in the middle of the two short sides, that of the 3rd order at the four corners of the thin plate (Figure 12), that of the 8th order in the six concentration regions of the two short sides of the thin plate (Figure 13), and that of the 10th order at the four corners of the thin plate, in the region along the long center axis and in the four concentration regions of the two long sides (Figure 14).

Figure 11: The maximum shear of the 1st order mode (S-F-S-F).

Figure 12: The maximum shear of the 3rd order mode (S-F-S-F).

Figure 13: The maximum shear of the 8th order mode (S-F-S-F).

Figure 14: The maximum shear of the 10th order mode (S-F-S-F).

It could be analyzed that the process of damage and failure of the thin plate rock mass S-F-S-F under impact load may be as follows: firstly, plate cracks occurred in the middle of the two short sides and the four corner regions (Figures 11 and 12), and then these cracks would propagate rapidly along the two short sides (Figure 13) and also tended to expand inward along the direction of long side (Figures 13 and 14).

7. Discussion on Dynamic Failure Mechanism of the Thin Plate Rock Mass

According to Table 7, the vibration modes of the thin plate rock mass S-C-S-C under impact load were dominated by the 1st order (main frequency), following by the 5th, 6th, and 10th order. From Figures 2–6, the maximum shear of the 1st order was mainly distributed in the middle of the two short sides and four concentration regions of the two long sides where plate cracks occurred firstly and that of the 5th and 6th order in the regions along the long and central axes of the thin plate rock mass where the plate cracks tended to expand. Similarly, according to Table 5, the vibration modes of the thin plate rock mass C-C-C-C under impact load were dominated by the 1st order (main frequency), following by the 5th and 6th order. From Figures 7–10, the maximum shear of the 1st order was mainly distributed in the middle of the four sides where plate cracks occurred firstly and that of the 5th and 6th order in the regions along the long and central axes of the thin plate rock mass where the plate cracks tended to expand.

According to Table 3, the vibration modes of the thin plate rock mass S-F-S-F under impact load were dominated by the 1st order (main frequency), following by the 3rd, 8th, and 10th order. From Figures 11–14, it could be concluded by the maximum shear of the 1st order that plate cracks occurred firstly in the middle of the two short sides and the four corner regions, and these cracks would propagate rapidly along the two short sides and also tended to expand inward along the direction of long side based on that of the 3rd, 8th, and 10th order.

According to the above discussion, it could be concluded that the initial position of damage and failure of the thin plate rock mass may be determined by effective vibration mode with the lowest frequency, and the develop tendency of damage and failure may be determined by the combined actions of other effective vibration modes. Furthermore, we could obtain that the coefficient An/q₀ of 1st order is the largest (Tables 3, 5 and 7), and it could be concluded that the effective vibration modes with lower frequency could have greater influence on the process of damage and failure of the thin plate rock mass.

8. Conclusions

(1)Based on the thin plate model and the dual equation of the Hamilton system, the frequency equations and deflection equations of the thin plate rock mass with different boundary conditions (S-F-S-F, S-C-S-C, and C-C-C-C) were established under free vibration, and the deflection equations under impact load were derived by the orthogonality of vibration modes and Duhamel integral.(2)Taking granite in a mine in Si Chuan, the effective vibration modes of the thin plate rock mass with different boundary conditions and their natural frequencies were obtained by the frequency equations and Newton’s iterative method.(3)Based on the third-strength theory, the maximum shear of the effective vibration modes was obtained, and the processes of damage and failure of the thin plate rock mass under impact load were analyzed with the numerical simulation results by LS-DYNA.(4)Based on the analysis of damage and failure of the thin plate rock mass under impact load, we could concluded that the initial position of damage and failure may be determined by effective vibration mode with the lowest frequency, the develop tendency of which by the combined actions of other effective vibration modes. In addition, the effective vibration modes with lower frequency could have greater influence on the process of damage and failure of the thin plate rock mass.

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 study was financially supported by the National Natural Science Foundation of China (51804311 and 52064046), China National University of Mining and Technology (Beijing), State Key Laboratory of Deep Geotechnical Mechanics and Underground Engineering, Open Fund Project (SKLGDUEK1818), and the Fundamental Research Funds for the Central Universities (2020YJSAQ13).

References

Y. S. Liu, C. H. He, H. L. Fu et al., “Study on dynamic and static mechanical properties and failure mode of layered slate,” Journal of Railway Science and Engineering, vol. 17, no. 11, pp. 2789–2797, 2020.

P. T. Wang, M. F. Cai, and F. H. Ren, “Anisotropy and directionality of tensile behaviours of a jointed rock mass subjected to numerical Brazilian tests,” Tunnelling and Underground Space Technology, vol. 73, pp. 139–153, 2018.

D.-P. Xu, X.-T. Feng, D.-F. Chen, and Q.-X. Fan, “Constitutive representation and damage degree index for the layered rock mass excavation response in underground openings,” Tunnelling and Underground Space Technology, vol. 64, pp. 133–145, 2017.

G. Zhang, W. Shen, K. Cui et al., “Degradation behavior and mechanism of slate under alternating conditions of freeze-thaw and wet-dry,” Journal of Central South University (Science and Technology), vol. 50, no. 6, pp. 1392–1402, 2019.

P. Zhang, C. Yang, H. Wang et al., “Stress-strain characteristics and anisotropy energy of shale under uniaxial compression,” Rock and Soil Mechanics, vol. 39, no. 6, pp. 2106–2114, 2018.

Y. Liu, S. Wang, S. Yan et al., “Properties and failure mechanism of layered sandstone based on acoustic emission experiments,” Journal of Central South University (Science and Technology), vol. 50, no. 6, pp. 1419–1427, 2019.

P. Zhou, Y. Shen, J. Zhao et al., “Research on disaster-induced mechanism of tunnels with steeply dipping phyllite strata based on an improved ubiquitous-joint constitutive model,” Chinese Journal of Rock Mechanics and Engineering, vol. 38, no. 9, pp. 1870–1883, 2019.

M. Qian, X. Miao, and J. Xu, Theory of Keystratum in Ground control, China University of Mining and Technology Press, Xuzhou, China, 2003.

X. Jia, Rock Mechanics and Strata Control, China University of Mining and Technology Press, Xuzhou, China, 2010.

Z. Song, Practical Mine Pressure Control, China Universityof Mining and Technology Press, Xuzhou, China, 1988.

J. Jiang, “The fracturing rules of main roof strata as a plant structure,” Journal of Shandong Mining Institute, vol. 1, pp. 51–58, 1988.

J. Zuo, M. Yu, S. Hu et al., “Experimental investigation on fracture mode of different thick rock strata,” Journal of Mining and Strata Control Engineering, vol. 1, no. 1, Article ID 013007, 2019.

E. Hoek, “Progressive caving induced by mining an inclined corebody,” Transactions of the Institution of Mining and Metallurgy, vol. 83, pp. 133–139, 1974.

S. Xie, D. Chen, Y. Sun et al., “Analysis on thin plate model of basic roof at elastic foundation boundary(I) : first breaking,” Journal of China Coal Society, vol. 41, no. 6, pp. 1360–1368, 2016.

F. He, J. Zhao, and Z. Yao, The Control Theory of Mining field, Metallurgical Industry Press, Beijing, China, 2009.

D. Chen, F. He, S. Xie et al., “First fracture of the thin plate of main roof with three sides elastic foundation boundary and one side coal pillar,” Journal of China Coal Society, vol. 42, no. 10, pp. 2528–2536, 2017.

H. Tu, S. Tu, C. Fang et al., “Study on the deformation and fracture feature of steep inclined coal seam roof based on the theory of thin plates,” Journal of Mining & Safety Engineering, vol. 31, no. 1, pp. 49–54, 2014.

X. Wang and M. Gao, “Mechanical model of fracture mechanism of stope roof for working face with variable length,” Journal of China University of Mining & Technology, vol. 44, no. 1, pp. 36–45, 2015.

X. Li, F. Gao, and W. Zhong, “Analysis of fracturingmechanism of stope Roof based on plate model,” Journal of Mining& Safety Engineering, vol. 25, no. 2, pp. 180–183, 2008.

H. Wu, “Structural models of longwall stope and their algorithm,” Ground Pressure & Strata Control, vol. 1, pp. 10–13, 1997.

H. Pu, Y. Huang, and R. Chen, “Mechanical analysis for X-O type fracture morphology of stope roof,” Journal of China University of Mining & Technology, vol. 40, no. 6, pp. 835–840, 2011.

F. Feng, X. B. Li, D. Y. Li et al., “Mechanism and control strategy of buckling rock bursts of orthotropic slab,” Chinese Journal of Geotechnical Engineering, vol. 39, no. 7, pp. 1302–1311, 2017.

L. Yuan, “Control of coal and gas outbursts in Huainan mines in China: a review,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 8, pp. 559–567, 2016.

W. D. Ortlepp, “The behaviour of tunnels at great depth under large static and dynamic pressures,” Tunnelling and Underground Space Technology, vol. 16, no. 1, pp. 41–48, 2001.

Y. Potvin, J. Hadjigeorgiou, and D. Stacey, Challenges in Deep and High Stress mining, Australian Center for Geomechanics, Nedlands, Australia, 2007.

W. Blake and D. G. F. Hedley, Rockbursts: Case Studies from North American Hard-Rock Mines, Society for Mining, Metallurgy & Exploration, Englewood, CO, USA, 2003.

X.-c. Zhang and J.-q. Wang, “Research on the mechanism and prevention of rockburst at the Yinxin gold mine,” Journal of China University of Mining and Technology, vol. 17, no. 4, pp. 541–545, 2007.

D. Li, C. C. Li, and X. LI, “Influence of sample height-to-width ratios on failure mode for rectangular prism samples of hard rock loaded in uniaxial compression,” Rock Mechanics and Rock Engineering, vol. 44, no. 3, pp. 253–267, 2011.

M. Tao, X. Li, and D. Li, “Rock failure induced by dynamic unloading under 3D stress state,” Theoretical and Applied Fracture Mechanics, vol. 65, pp. 47–54, 2013.

X. Li, K. Du, and D. Li, “True triaxial strength and failure modes of cubic rock specimens with unloading the minor principal stress,” Rock Mechanics and Rock Engineering, vol. 48, no. 6, pp. 2185–2196, 2015.

Z.-S. Hou, Q.-M. Gong, and Z.-H. Sun, “Primary failure types and their failure mechanics of deep buried and intact marble at Jinping II hydropower station,” Chinese Journal of Rock Mechanics and Engineering, vol. 30, no. 4, pp. 727–732, 2011.

S.-Y. Wu, Q.-M. Gong, G. Wang et al., “Experimental study of slabbing failure for deep-buried marble at Jinping II hydropower station and its influence of TBM excavation,” Chinese Journal of Rock Mechanics and Engineering, vol. 29, no. 6, pp. 1089–1095, 2010.

F. Li, M. X. Bi, F. C. Zhang et al., “Dynamic response characteristics and damage mechanisms of coal-rock masses under P-waves,” Journal of China Coal Society, vol. 43, no. 5, pp. 1272–1280, 2018.

G. Su, L. Hu, X. Feng et al., “True triaxial experimental study of rockburst process under low frequency cyclic disturbance load combined with static load,” Chinese Journal of Rock Mechanics and Engineering, vol. 35, no. 7, pp. 1309–1322, 2016.

F. Cai, Z. Liu, and Y. Luo, “Propagation and attenuation characteristics of stress waves generated by explosion in high-gas coal-beds,” Journal of China Coal Society, vol. 39, no. 1, pp. 110–114, 2014.

W. Chang, S. Fan, and L. Zhang, “A model based on explosive stress wave and tectonic coal zone which gestate dangerous state of coal and gas outburst,” Journal of China Coal Society, vol. 39, no. 11, pp. 2226–2231, 2014.

L. Dou, X. He, and T. Ren, “Mechanism of coal-gas dynamic disasters caused by the superposition of static and dynamic loads and its control technology,” Journal of China University of Mining & Technology, vol. 47, no. 1, pp. 48–59, 2018.

M. Luo, C. Fan, Li Sheng et al., “Failure criteria of the geological dynamic system of coal and gas outburst,” Journal of China Coal Society, vol. 47, no. 1, pp. 137–144, 2018.

F. Li, F. Sha, M. Zhang, Z. He, X. Dong, and B. Ren, “Dynamic failure model of thin coal and rock mass under uniform impact load,” Engineering, vol. 12, no. 10, pp. 699–714, 2020.

F. Li, M. Bi, and J. Tian, “The dynamic damage mechanisms and failure modes of coal-rock masses under the action of high order P-waves,” Shock and Vibration, vol. 2018, Article ID 5386123, 11 pages, 2018.

Z. Y. Fang, Analytical Solution for the Forced Vibration of Orthotropic Rectangular Thin Plates, Da Lian University of Technology, Dalian, China, 2011.

W. A. Yao and W. X. Zhong, Symplectic Elastic Mechanics, Higher Education Press, Beijing, China, 2002.

Y. Xing and B. Liu, “New exact solutions for free vibrations of rectangular thin plates by symplectic dual method,” Acta Mechanica Sinica, vol. 25, no. 2, pp. 265–270, 2009.

S. Y. Bao and Z. C. Deng, “A general solution of free vibration for rectangular thin plates in Hamilton systems,” Journal of Dynamics and Control, vol. 3, no. 2, pp. 10–16, 2005.

S. Y. Bao and Z. C. Deng, “High-precision approximate analytical solutions for free bending vibrations of thin plates and an improvement,” Applied Mathematics and Mechanics, vol. 37, no. 11, pp. 1169–1180, 2016.

R. Li, B. Wang, and P. Li, “Hamiltonian system-based benchmark bending solutions of rectangular thin plates with a corner point‐supported,” International Journal of Mechanical Sciences, vol. 85, pp. 212–218, 2014.

G. Z. Yang and C. S. Lin, “Forced vibration analysis of rectangular thin plates with mixed boundary,” Journal of Chifeng University, vol. 21, no. 4, pp. 16-17, 2005.

M. Zhang, L. Cui, W. Hu, J. Du, Z. Zhang, and C. Wang, “Acoustic emission experimental research of the damage characteristics of raw coal under different loading and unloading rates,” Shock and Vibration, vol. 2020, no. 7, Article ID 9063929, 15 pages, 2020.