^{1}

^{2}

^{2}

^{3}

^{2}

^{2}

^{2}

^{1}

^{2}

^{3}

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.

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 [

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 [

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 [

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 [

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.

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

For the purpose to solve the homogeneous equation of equation (

The vibration of the thin plate is the harmonic oscillator with time. Therefore,

With (^{4} = ^{2}/

In order to decouple the physical quantities in equation (_{x}, _{y}, _{xy}, _{x}, _{y}, _{x}, and _{y} are the bending moment, torque, transverse shear force, and transverse total shear force, respectively. If we set vector _{x}, _{x}], then the dual equation of the Hamilton system could be expressed as follows [

Then,

Based on the symplectic geometry method, the solutions of (_{1} and _{2} are characteristic values in the _{1} and _{2} are characteristic values in the _{1}, _{2}, _{1}, and _{2} meet the following rules:

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.

The deflection equation of the thin plate rock mass with both sides simply supported (S-S) along the

The boundary conditions of the thin plate rock mass with both free sides (F-F) along the

With (

Then, we insert (_{1}−_{1}^{2}, _{2}−_{1}^{2}, _{1}^{3}−(2−_{1}^{2}_{1}), _{2}^{3}−(2−_{1}^{2}^{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

We set _{2} = 1and insert it into (

Then,

As for the thin plate rock mass with four-sided fixed supported (C-C-C-C), the boundary conditions C-C along the

With (

And the deflection equation could be expressed as follows:

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

^{2}^{2} = 0; ^{2}^{2} = 0.

With (

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

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

Thin plate rock mass | Frequency equations | Deflection equations |
---|---|---|

_{1} | ||

_{1}−_{1}^{2}, _{2}−_{1}^{2}, _{1}^{3}−(2−_{1}^{2}_{1}), _{2}^{3}−(2−_{1}^{2}_{2,} | ||

_{1} | ||

The solutions’ form of the nonhomogeneous control equation (

We insert (

With (

Then,

The orthogonality of the deflection equations may be as follows [

We multiply both sides of (_{m}(

We set

Then,

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

If the external dynamic load

Then,

Therefore, the solutions of control equation (see (

Taking granite in a mine in Sichuan, the parameter values of the thin plate rock mass were as follows: ^{3}, ^{9} Pa,

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 ^{st} order, 3^{rd} order, 8^{th} order, and 10^{th} 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 _{n}/_{0} of the 3^{rd} order, 8^{th} order, and 10^{th} order was about 6.4%, 4.5%, and 0.5% of the 1^{st} order.

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

Parameter | 1^{st} order | 2^{nd} order | 3^{rd} order | 4^{th} order | 5^{th} order | 6^{th} order | 7^{th} order | 8^{th} order | 9^{th} order | 10^{th} order |
---|---|---|---|---|---|---|---|---|---|---|

_{1} | 1.047 | 1.047 | 1.047 | 2.094 | 2.094 | 1.047 | 2.094 | 3.142 | 2.094 | 1.047 |

_{1} | ≈0 | 0.716 | 1.469 | ≈0 | 0.75 | 2.268 | 1.578 | ≈0 | 2.382 | 3.105 |

_{2} | 1.481 | 1.645 | 2.086 | 2.961 | 3.056 | 2.709 | 3.356 | 4.443 | 3.801 | 3.44 |

1.047 | 1.269 | 1.804 | 2.094 | 2.225 | 2.498 | 2.622 | 3.142 | 3.172 | 3.277 |

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

Order | Parameter | |||
---|---|---|---|---|

_{n} | _{n}/_{0} | |||

^{st} | ^{4} | ^{−5} | ||

2^{nd} order | 148 | 4.536 .510^{3} | −4.642 .610^{−6} | −6.914 .910^{−12} |

^{rd} | ^{3} | − | −^{−7} | |

4^{th} order | 404 | 2.197 .110^{4} | 9.059 .010^{−10} | 1.020 .010^{−16} |

5^{th} order | 456 | 8.646 .610^{3} | 9.876 .810^{−20} | 2.505 .510^{−26} |

6^{th} order | 575 | 1.208 .210^{3} | 2.881 .810^{−6} | 4.148 .110^{−12} |

7^{th} order | 633 | 3.684 .610^{3} | −9.813 .810^{−15} | −4.208 .210^{−21} |

^{th} | ^{4} | ^{−7} | ||

9^{th} order | 926 | 2.082 .010^{3} | 6.097 .010^{−21} | 3.163 .110^{−27} |

^{th} | ^{3} | − | −^{−8} |

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 ^{st} order, 5^{th} order, and 6^{th} order were effective vibration modes, and the natural frequencies of which were 308 rad/s, 975 rad/s, and 1309 rad/s, respectively (Table _{n}/_{0} of the 5^{th} order and 6^{th} order were about 9.4% and 4.9% of the 1^{st} order.

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

Parameter | 1^{st} order | 2^{nd} order | 3^{rd} order | 4^{th} order | 5^{th} order | 6^{th} order | 7^{th} order | 8^{th} order | 9^{th} order | 10^{th} order |
---|---|---|---|---|---|---|---|---|---|---|

_{1} | 1.443 | 1.301 | 2.573 | 2.483 | 1.227 | 3.643 | 2.404 | 3.589 | 2.346 | 2.346 |

_{2} | 2.145 | 3.247 | 2.949 | 3.765 | 4.437 | 3.895 | 4.802 | 4.514 | 5.930 | 5.930 |

_{1} | 1.122 | 2.104 | 1.019 | 2.001 | 3.015 | 0.975 | 2.940 | 1.935 | 3.904 | 3.851 |

_{2} | 2.329 | 2.795 | 3.778 | 4.042 | 3.479 | 5.243 | 4.494 | 5.433 | 4.248 | 5.083 |

1.828 | 2.474 | 2.767 | 3.189 | 3.255 | 3.771 | 3.797 | 4.078 | 4.079 | 4.509 |

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

Order | Parameter | |||
---|---|---|---|---|

_{n} | _{n}/_{0} | |||

1^{st} order | 308 | 5.787 .710^{3} | 13.589 | 7.624 .610^{−6} |

5^{th} order | 975 | 1.267 .210^{4} | 8.503 | 6.883 .810^{−7} |

6^{th} order | 1309 | 2.587 .510^{4} | 12.565 | 3.71 .710^{−7} |

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 ^{st}, 5^{th}, 6^{th}, and 10^{th} 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 _{n}/_{0} of the 5^{th} order, 6^{th} order, and 10^{th} order were about 10.6%, 2.8%, and 1.3% of the 1^{st} order.

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} | 1.047 | 2.094 | 1.047 | 2.094 | 1.047 | 3.142 | 2.094 | 3.142 | 1.0472 | 3.142 |

_{1} | 1.183 | 1.052 | 2.127 | 2.031 | 3.025 | 0.992 | 2.960 | 1.959 | 3.908 | 2.896 |

_{2} | 1.896 | 3.143 | 2.592 | 3.591 | 3.368 | 4.552 | 4.187 | 4.855 | 4.180 | 5.303 |

1.580 | 2.343 | 2.371 | 2.917 | 3.201 | 3.294 | 3.626 | 3.702 | 4.046 | 4.273 |

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

Order | Parameter | |||
---|---|---|---|---|

_{n} | _{n}/_{0} | |||

1^{st} order | 230 | 1.407 .410^{3} | 7.171 | 2.511 .510^{−5} |

5^{th} order | 943 | 9.980 .910^{2} | 2.540 | 2.699 9610^{−6} |

6^{th} order | 999 | 8.898 .810^{3} | 6.171 | 6.943 .910^{−7} |

10^{th} order | 1680 | 1.856 .810^{3} | 1.045 | 3.352 .310^{−7} |

According to the deflection equations of the thin plate rock mass under different boundary conditions shown in Table

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 ^{st} order was mainly distributed in the middle of the two short sides and four concentration regions of the two long sides, that of the 5^{th} order in the middle of the two short sides and the region along the long central axis of the thin plate (Figure ^{th} order in the region along the short center axis of the thin plate (Figure ^{th} order in any region (Figure

The maximum shear of the 1^{st} order mode (S-C-S-C).

The maximum shear of the 5^{th} order mode (S-C-S-C).

The maximum shear of the 6^{th} order mode (S-C-S-C).

The maximum shear of the 10^{th} 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

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

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 ^{st} order was mainly distributed in the middle of the four sides, that of the 5^{th} order in the middle of the two short sides and the region along the long central axis of the thin plate (Figure ^{th} order in the middle of the two long sides and the region along the short central axis of the thin plate (Figure

The maximum shear of 1^{st} order mode (C-C-C-C).

The maximum shear of the 5^{th} order mode (C-C-C-C).

The maximum shear of the 6^{th} 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

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

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 ^{st} order was mainly distributed in the middle of the two short sides, that of the 3^{rd} order at the four corners of the thin plate (Figure ^{th} order in the six concentration regions of the two short sides of the thin plate (Figure ^{th} 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

The maximum shear of the 1^{st} order mode (S-F-S-F).

The maximum shear of the 3^{rd} order mode (S-F-S-F).

The maximum shear of the 8^{th} order mode (S-F-S-F).

The maximum shear of the 10^{th} 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

According to Table ^{st} order (main frequency), following by the 5^{th}, 6^{th}, and 10^{th} order. From Figures ^{st} 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 5^{th} and 6^{th} 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 ^{st} order (main frequency), following by the 5^{th} and 6^{th} order. From Figures ^{st} order was mainly distributed in the middle of the four sides where plate cracks occurred firstly and that of the 5^{th} and 6^{th} 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 ^{st} order (main frequency), following by the 3^{rd}, 8^{th}, and 10^{th} order. From Figures ^{st} 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 3^{rd}, 8^{th}, and 10^{th} 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 _{n}/_{0} of 1^{st} order is the largest (Tables

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.

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.

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.

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.

The data used to support the findings of this study are available from the corresponding author upon request.

The authors declare that they have no conflicts of interest.

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).