Abstract

During the Wenchuan Earthquake, with a magnitude of 5.12, collapses and rockfall hazards persisted for a long time after the initial investigations carried out by research fellow S. M. He and his team at the scene of the disaster in October 2008. It is possible that additional incidents of rockfalls in large quantities may continue in the same areas over the next ten to fifteen years. Furthermore, in the vast mountainous region of western China, the topographic relief is evident, and earthquakes occur frequently. Therefore, it is difficult to effectively defend against rockfall hazards. When designing protective structures, the key issue is the analysis of the mechanical response mechanism of the soil cushion layer of the upper cushion when subjected to the impact of rockfall. As such, a theoretical method was used to perform such an analysis. The cavity expansion and energy conservation model were adopted. Analytical solutions for the impact force and penetration depth were then derived. Furthermore, the impact force and penetration depth of rockfall were studied with the LS-DYNA software to obtain values for the impact forces and the penetration depth. Finally, the reliability of the theoretical method was evaluated using the cavity expansion, energy conservation, numerical simulation, Hertz, Japanese, Swiss, Australian, B. S. Guan, tunnel manual, and subgrade methods based on an engineering model. The results show that the cavity expansion and the energy conservation methods yielded consistent results. Meanwhile, the cavity expansion and the energy conservation methods also yielded consistent results with the numerical simulation, Japanese (obtained by laboratory experiment), Swiss (obtained by laboratory experiment), and Australian (obtained by field experiment) methods. The relevant methods and conclusions shall therefore be applied to the design of rockfall protection structure in future investigations.

1. Introduction

As a key aspect in the engineering of highway routes, the technology for building tunnels has gradually improved over time [13]. Shed tunnel protection structures are special tunnel structures. These structures and others such as open cut tunnels and passive protective structures in hazardous areas are designed primarily to protect against the impact forces of rockfall. In order to avoid immediate impact load on protected structures, a soil layer of a certain thickness is generally used [4]. During the rockfall impact on the soil cushion layer, the time interval is short and the energy conversion is complicated. In addition, elastoplastic deformation of the soil cushion layer occurs. Therefore, the impact process of the rockfall is a complex impact dynamics problem. In the design of protective structures, the key issue is the analysis of the mechanical response mechanism of the soil cushion layer under the impact of rockfall.

Many scholars and research institutes are directing increasing attention to rockfall disaster research. The classic Hertz contact theory and Thornton elastoplastic impact theory have been applied widely [58]. Johnson described in his book “Contact Mechanics” the details of the mechanical properties of an object under the effects of static contact, sliding, rolling, and impact [9]. Labiouse et al. presented a semiempirical and semitheoretical method to analyze rockfall impact forces based on impact experiments (hereinafter referred to as the Swiss method) [10]. Kishi and Ikeda studied the mechanical response characteristics of a retaining wall under rockfall impact based on field experiments and numerical simulations [11]. Heidenreich and Labiouse studied the mechanical characteristics of a slope with a loose structure under the impact of rockfall using small-scale model experiments [12]. Calvetti and Prisco analyzed the mechanical response characteristics of gravel soil under rockfall impact based on experiments and numerical simulations [13]. Based on the field test method and the theorem of impulse, Pichler and Hillmich developed a calculation method for impact force and penetration depth associated with rockfall (hereinafter referred to as the Australian method) [7]. Mougin et al. examined the mechanical response of rockfall on a concrete plate under direct impact conditions based on model experiments [14]. Delhomme et al. adopted model experiments and numerical simulations to analyze the energy dissipation and seismic mitigation effects of the new type of shed tunnel structures [15]. The Japanese Road Association developed a semiempirical and semitheoretical calculation method for rockfall impact force using model experiments (hereinafter referred to as the Japanese method) [16]. Calvetti and Prisco analyzed the design method for shed tunnel structures based on the uncoupled calculation method [17]. Ye and Chen performed a comparative study on rockfall impact force values by employing the Japanese, Swiss, and Hertz methods in conjunction with an engineering model [18].

Ye et al. analyzed the rockfall impact forces under oblique impact on the soil cushion layer [19]. Qi et al. set up an analysis model for rockfall impact forces based on a theoretical approach and studied the amplification coefficient of rockfall impact forces [20]. Yuan et al. analyzed the influence of the weight, shape, falling height, and impact angle of rockfall on the impact force based on experiments [21]. Sun et al. attempted to add discarded tires to the soil cushion layer and investigated the impact resistance of the new cushion structure via an experimental approach [22]. Zhang et al. investigated the mechanical response mechanism of oil/gas pipe under rockfall impact using a numerical simulation method [23]. Hu et al. discovered the law of jumping range of rockfall based on three influential factors: shape, weight, and the release altitude of rockfall using an experimental method [24]. Lam et al. studied the overturning stability of L-shaped rigid barriers subjected to rockfall impacts [25]. Basharat et al. analyzed the effects of volume and topographic parameters on rockfall travel distance based on a case study from the NW Himalayas, Pakistan [26].

Zhu et al. analyzed particle size and the thickness of the gravel cushion layer on the restitution coefficient of rockfall in rockfall impact procession [27]. Luo et al. experimentally studied the rockfall impact force of the frame shed tunnel structure [28]. Asteriou and Tsiambaos analyzed the effect of rockfall impact velocity, block mass, and hardness on rockfall restitution coefficients [29]. Toe et al. adopted a metamodel hazard mitigation strategy to improve the effectiveness of rockfall protection barriers [30]. Wang et al. adopted an experimental approach to analyzing the mechanical response law of an EPS cushion layer under rockfall impact [31]. Bhatti combined experimental and numerical simulation methods to study the energy dissipation property of the EPS cushion layer under the impact of rockfall [32]. Bi et al. adopted a numerical simulation method to analyze the mechanical response mechanism of a bumper plate under the rockfall impact [33]. Yu et al. studied the maximum impact force of rockfall using an experimental approach [34]. Zhang et al. devised an estimation model for calculating the maximum impact force of rockfall based on contact theory and verified this model experimentally [35]. Yan et al. utilized numerical simulations to analyze the mechanical response mechanism of reinforced concrete slabs under the impact of rockfall [36].

Prof. B. S. Guan proposed a calculation method for impact forces and penetration depth based on experiments. Further, he considered the effect on the thickness of the soil cushion layer (hereinafter referred to as the B. S. Guan method) [37]. Using the theory of momentum, the impact force and penetration depth are given in the Technical Manual for Railway Engineering Design Tunnel (Revised Edition) (hereinafter referred to as the tunnel manual method). In addition, in Code for Design of Highway Subgrade of China, the rockfall impact force and penetration depth are given based on the work-energy principle (hereinafter referred to as the subgrade method).

The Japanese, Swiss (based on indoor experiments), and Australian methods (based on field tests) are three classical methods for the design of protection structures. However, these three methods are all semitheoretical and semiempirical. Moreover, few scholars have focused on theoretical derivations to obtain the penetration depth and impact force. Meanwhile, some researchers have indicated that the value of the impact force is small in China [18].

Based on the dynamic cavity expansion model and the energy conservation model, a theoretical mechanical system of rockfall impacting the soil cushion layer is established and analytic solutions of rockfall impact force and penetration depth are obtained. In addition, the LS-DYNA software is used to analyze the change laws of impact force and penetration depth during rockfall impact processes. Finally, this study will compare and analyze the calculation results for penetration depth and impact forces calculated based on the cavity expansion, energy conservation, numerical simulation, Japanese, Swiss, Australian, Hertz, B. S. Guan, tunnel manual, and subgrade methods.

2. Conventional Classical Calculation Theories

2.1. Hertz Method

In 1985, Hertz derived expressions for the maximum deformation and maximum impact based on elastic theory.where is the rockfall impact force, is the elasticity modulus of soil cushion layer, is the radius of the rockfall, is the mass of the rockfall, is the initial velocity of rockfall, and is the penetration depth of the rockfall impact.

2.2. Swiss Method

In 1996, a Swiss scholar, Labiouse, developed the relevant empirical method for rockfall impact forces based on rockfall experiments.where is the deformation modulus of soil cushion layer, is the quality of rockfall, and is the falling height of the rockfall.

2.3. Japanese Method

In 2000, based on rockfall experiments and combined with Hertz elastic theory, the Japanese Road Association presented the relevant semiempirical and semitheoretical calculation method of rockfall impact forces.where is the gravitational acceleration, is the lame constants, and is the Poisson’s ratio of soil cushion layer.

2.4. Australian Method

In 2005, Australian scholar Pichler proposed the semiempirical and semitheoretical calculation method of rockfall impact force and penetration depth combined with the rockfall field test.where is the impact time of the rockfall, is the diameter of the rockfall, and is the indentation resistance of target materials.

2.5. B. S. Guan Method

In 1996, the scholar B. S. Guan established the empirical method of impact force and penetration depth based on laboratory experiments. This method considers the influence of the thickness of the soil cushion layer on impact force.where is the correction coefficient of the rockfall impact force, is the coefficient of rockfall penetration depth, is the acceleration of the rockfall impact, and is the thickness of soil cushion layer.

2.6. Method of Tunnel Manual

The method of tunnel manual refers to an approach recommended in Technical Manual for Railway Engineering Design Tunnel (Revised Edition), and in essence, it is also an approximation method of the theorem of momentum.where is the reciprocating velocity of compression waves in the soil cushion layer, is the specific gravity of the soil cushion layer, is the cross-sectional area of an equivalent sphere for rockfall, is the inner friction angle of soil cushion layer, and ρ0 is the density before deformation of spherical cavity microbody.

2.7. Subgrade Method

Code for Design of Highway Subgrade provides a computational formula for rockfall impact force and penetration depth based on the work-energy principle.

3. Theoretical Analysis Based on Dynamic Cavity Expansion Model

The spherical expansion model is widely applied to penetration mechanical model analysis. In this section, the impact of rockfall on a soil cushion layer is simplified by considering the normal impact of a rigid ball on a semi-infinite soil cushion layer.

By referring to Figure 1, the normal stress of point A on the soil cushion layer σA could be calculated using of the cavity expansion sphere at the same position A, so a series of spherical cavities originating from the bottom of the sphere could be used to simulate the process of crater formation due to rockfall impact on the soil cushion layer. The radius of these cavity spheres gradually increases from 0 to rq, thus forming the so-called dynamic spherical cavity expansion model.

The radius of cavity expansion sphere iswhere is the radius of the spherical cavity, is the x coordinate value of point A, and is the acute angle between the normal line and the vertical line at point A.

During the normal impact of rockfall on the soil cushion layer at a velocity of , the velocity of the dynamic spherical cavity expansion iswhere is the velocity of spherical cavity expansion, is the expansion velocity of spherical cavity, and is the velocity of rockfall impact.

And the acceleration of the dynamic cavity expansion iswhere is the acceleration of spherical cavity expansion, is the expansion acceleration of spherical cavity, and is the acceleration of a rockfall impact.

3.1. Mechanical Analysis of Microunit with Spherical Coordinates

As shown in Figure 2, the Euler coordinate of spherical coordinates is employed to perform mechanical analysis on the selected microbody.

By referring to Figure 2, a mechanical analysis of the microbody selected yields equation (7) according to the theorem of momentum:where is the radial stress of the spherical cavity microbody, is the hoop stress of the spherical cavity microbody, is the distance from the coordinate origin to the mass point of the spherical cavity microbody, is the moving velocity of the mass point of the spherical cavity microbody, and is the opening angle of spherical cavity microbody.

The velocity of the mass point of the microbody in this coordinate iswhere is the displacement of the mass point of the spherical cavity microbody and is the moving time of the mass point of the spherical cavity microbody.

And the acceleration of the mass point could be expressed by formula (8):

For the smaller , could be expressed as approximately . Substituting this simplified relation into formula (11) and dividing both sides of the simplified equation by r2dr dθ2, if the high order trace is omitted, considering equations (11) and (13), we obtain

Figure 3 represents the force and displacement of a spherical cavity microelement.

Combined with the law of conservation of mass, we obtainwhere is the density after deformation of spherical cavity microbody.

After the partial derivative of equation (15), we obtain

Taking the integral of equation (16) and assuming the medium of the plastic zone is incompressible, we obtain

Further, combined with the boundary conditions of , , we obtain the displacement of the microbody, as follows:where is the radius of spherical cavity.

Combined with equations (18) and (12), we obtain an expression for the velocity of the microbody, as follows:

According to the continuous medium assumption, the expression of for the velocity shall be applied to the medium of both elastic and plastic zones.

3.2. Mechanical Analysis of the Elastic Zone

After rockfall impact on the soil cushion layer, the elastic zone formed within this layer has the following relationships:(1)Physical equation(2)Geometric equationwhere is the radial strain of spherical cavity microbody and is the hoop strain of spherical cavity microbody.

Substituting the expression for the displacement of the microbody s into equation (21) and using a Taylor series to expand, only the first two terms of the expansion are substituted into equation (20) on the premise of ensuring the effectiveness of the calculation results. We then obtain

Combining equations (22), (19), and (14), we obtain

Considering the boundary conditions r ⟶ ∞ and σr = 0 and taking the integral of equation (23), we obtain

3.3. Mechanical Analysis of Plastic Zone

For the ideal elastic-plastic soil cushion layer in the case of shear failure of the soil cushion layer, combined with Thornton’s theory, equation (25) is established:

Combined with the dynamic strength theory proposed by Richart, can be determined using equation (26):where is the shear strength of soil cushion layer, is the undrained shear strength of soil cushion layer, and is the shear strength coefficient of soil cushion layer.

Considering equations (22), (19), and (14), we obtain

Taking the boundary conditions and into account, we obtainwhere is the radial stress on the spherical cavity surface.

Using the boundary conditions of and according to the assumption of the soil cushion layer being in the elastic-plastic zone and substituting the boundary conditions into equation (24) and (28), respectively, we can obtain an expression for the normal stress on the surface of cavity .where is the plastic region radius of spherical cavity.

Combining equations (22) and (25), we obtain

Substituting equation (30) into (29), we obtain

The process of rockfall impact on cratering can be imitated through a series of infinitesimal cavity expansion originating from the bottom of the sphere. As shown in Figure 4, for the horizontal infinitesimal cavity that lies at point , the cavity radius gradually increases from 0 to , and then forms the final infinitesimal cavity . The formation process for the infinitesimal cavity is also the same. Further, we perform an integration to analyze the entire process for the impact of rockfall on Earth mass.

From Figure 4, the velocity and acceleration of the dynamic cavity expansion arewhere is the acute angle between the tangent at the spherical cavity surface through point with x axis.

The value of cos φ is determined using relevant geometrical relations. A tangent line of the cavity is made through point to obtain the geometrical relation expression as follows:where is the x coordinate value of point on spherical cavity and is the y coordinate value of point on spherical cavity.

Since L1 is the tangent line through point , we obtainwhere is the slope of the tangent line L1.

According to the trigonometric function relation,and we can combine equations (33)–(35) to obtain

3.4. Impact Force

Equation (31) is the expression of the normal stress on the surface of cavity , over which the integral is performed so as to obtain the impact force of the rockfall.where is the radius of the horizontal infinitesimal cavity at the upper surface of soil cushion layer, and is the width value of infinitesimal cavity.

is calculated according to the equation below:

Combined equations (31), (32), (37), and (38), and given that when rq=R and Lmax ≤ R, the impact force of rockfall on the soil cushion layer shall be calculated according to equation (39):where is the rockfall impact force.

Generally, the velocity of rockfall for the calculation is 4–24 m/s. For the convenience of calculation, we omit and in equation (31), so equation (31) can be rewritten asand equation (39) for the calculation of the impact force can be simplified intowhere is the rockfall impact force.

During the rockfall impact process of the rolling stone, when Lmax ≥ R, the value of can be substituted into that of in equation (37) and (41).

3.5. Penetration Depth

Combined with Newton’s second law, we obtain the relation between the impact of velocity of rockfall and penetration depth:

During the impact process, if the velocity and acceleration of the rockfall are not considered for the influence on penetration depth, and when R, we multiply the rockfall mass M on both ends of equation (42) and obtain the equation for the impact force as follows:where is the maximum penetration depth of the rockfall by the cavity method.

That is,

Combined with the initial conditions of and L = 0 and taking the integral of equation (44), we obtain

From equation (45), we can obtain the penetration depth of rockfall at any time when the velocity is , during the period in which impact velocity decreases from to 0. When  = 0, we obtain the maximum penetration depth of the rockfall Lmax.

If the calculation result satisfies Lmax ≥ R, in equation (44), we set L = R.

If the influence of rockfall velocity and acceleration during the impact process is considered and when Lmax ≤ R, we obtain

Substituting equations (39) into (47), we obtain

Considering the initial conditions of and , when the calculation result for the penetration depth does not satisfy Lmax ≥ R, in equation (45), we set L = R. Using an analytical method, we can obtain the relation between L and .

Using equation (48), we obtain

We can substitute the values of L and obtained using the analytical method into equation (49), to calculate the value of the impact force F.

4. Theoretical Derivation Based on Energy Conservation Model

In this section, the impact of rockfall on the soil cushion layer is simplified into the normal impact of a rigid ball on a semi-infinite soil cushion layer, and further derivation is based on the principle of conservation of energy.

During the process of rockfall impact on the soil cushion layer, there are four main forms of energy:(1)Initial kinetic energy of rockfall at the moment of contact with the cushion layer(2)Plastic deformation and friction generate forms of energy dissipation(3)Microgravitational potential energy of rockfall during the impact process

4.1. Surface Stress Analysis of Cylindrical Cavity

After rockfall impact on the cushion layer, we take the cross section of a cylindrical cavity with a radius of rq, and the plastic zone and elastic zone are shown in Figure 5.

We then consider the microbody as shown in Figure 6 during mechanical analysis. The value of displacement of the microbody is positive in the outward direction, and the values of the stress and and strain and are positive with pressure.

So, we obtain the equation below:where is the radial stress of the cylindrical cavity microbody, is the hoop stress of the cylindrical cavity microbody, and is the opening angle of cylindrical cavity microbody.

Further, in the elastic zone, the geometric relation is given by equation (51):where is the radial strain of cylindrical cavity microbody and is the hoop strain of cylindrical cavity microbody.

For the medium of the elastic zone, according to Hooke’s law and the reference Cavity Expansion Methods in Geomechanics [38], the elastic relation expression is given by equation (52):

Further, according to the law of the conservation of mass, we obtainwhere is the distance from the coordinate origin to the mass point of the cylindrical cavity microbody, is the displacement of the mass point of the cylindrical cavity microbody, is the density before deformation of cylindrical cavity microbody, and is the density after deformation of cylindrical cavity microbody.

Here, we consider that only the medium of the elastic zone is compressed, so we obtain

Considering the boundary conditions  =   =  and taking the integral of equation (65), we obtain the expression for the displacement of a microbody:where is the radius of the cylindrical cavity.

Substituting equation (55) in equation (51), and given that the radius of the plastic zone of the soil cushion layer is far less than the radius of elastic zone , we obtain

Combining equations (56) and (52), we obtain

Substituting equation (57) in (50) and considering the boundary conditions ⟶ ∞ and  = 0, we obtain

Assuming the soil cushion layer is an ideal elastic-plastic body, we obtain the failure criterion according to equation (59):

Substituting equation (59) into (50), we obtain

Further, consider the boundary conditions  = ,  = , we obtainwhere is the radial stress on the cylindrical cavity surface.

According to the boundary conditions  =  and  =  for the assumption of a continuous medium in the elastic-plastic zone and combined with the yield criterion of equation (59), we obtain t equation (62) for the surface stress of a cylindrical cavity:where is the plastic region radius of cylindrical cavity.

4.2. Energy Dissipation due to Crater Formation

Using the slices method, the energy dissipation due to crater formation is equal to the sum of energy works of the surface stress of each cylindrical cavity in the stripe microbody, as shown in Figure 7.where and are expressed by the equations below:where is the x coordinate value of point on the rockfall surface and is the y coordinate value of point on the rockfall surface.

Substitute equation (64) into (63), we obtain

Using the above equations, we can obtain the final energy dissipation due to crater formation of the rockfall:where is the energy dissipated due to the deformation of thin layer of soil cushion layer, is the total energy dissipated due to the deformation of soil cushion layer, and is the maximum penetration depth of the rockfall by energy method.

4.3. Energy Dissipation due to Friction

Combined with Figures 7 and 8, the slices method is adopted to analyze the internal energy dissipation due to the friction of rockfall and soil cushion layer.

In Figure 9, the friction of thin layer of microsoil cushion layer and rockfall Ff iswhere is the friction between the thin layer of soil cushion layer and the rockfall, is the friction coefficient between soil cushion layer and the rockfall, and is the acute angle between the inclined plane of the thin soil layer with the vertical line.

Combined with Figure 7, when the thin layer of soil cushion layer has the displacement relative to rockfall, the energy work of friction iswhere is the horizontal movement distance of the thin soil layer, is the thickness of the soil layer and is the energy dissipated due to friction between the rockfall and the soil cushion layer.

Further, the expression for the side profile of rockfall is

Combining equations (62), (68), and (69), we obtain

In equation (70), is the distance of movement of a thin layer of soil cushion relative to the rockfall. This can be determined using equation (71):

To determine as represented in Figure 7, a tangent line L2 is made through the contact point of the thin layer of soil cushion and rockfall and the included angle of L2. The horizontal axis is , so we obtainwhere is the slope of the tangent line L2 and is the acute angle between the inclined plane of the thin soil layer with the horizontal line.

Further, based on the trigonometric function relation,and combining the equations (72) and (73), we obtain

Combining equations (71) and (72), we obtain

Combining equations (74) and (75), we obtain

According to Figure 8, we can write the equation (77):and we can state that

According to the above equations, we obtain the final expression of equation (79):

Therefore, when rockfall has the maximum penetration depth of , the energy dissipation due to friction iswhere is the total energy dissipated due to friction between the rockfall and the soil cushion layer.

4.4. Energy Dissipation of Shed Tunnel Structure

Under the impact effect of rockfall load, the roof of a shed tunnel structure will have slight deflection deformation due to rockfall impact, and therefore, energy is dissipated. The deflection of the roof of the shed tunnel under an impact force Pf is [39]where is the roof deflection of the shed tunnel structure, is the Given rockfall impact load, is the elasticity modulus of the roof of the shed tunnel structure, is the inertia moment of the roof of the shed tunnel structure, and is the span of the shed tunnel structure.

According to structural mechanics and equation (81), we obtain the relevant energy of deformation as

Further, can be determined using equation (83):where is the yield strength of soil cushion layer and is the total energy dissipated due to the deflection deformation of the roof of the shed tunnel.

4.5. Penetration Depth

The initial kinetic energy of rockfall during the impact process is ; the energy of deformation due to crater shall be calculated using equation (66), the energy of friction using equation (80), and the energy of deflection deformation of the roof of the shed tunnel using equation (83). The work done by gravity during rockfall impact process is ; therefore, according to the law of energy conservation, we obtain:

There is only one unknown parameter in equation (76), so we obtain the solution using equation (84). The value of the impact force F is obtained based on .

4.6. Impact Force

According to the principle of work and energy, the relation expression during the contact of the rockfall and the soil cushion layer iswhere is the average impact force of rockfall.

The value determined using equation (85) is the average of the impact force. However, the rockfall impact process is an impulse type, Yang et al. [37], and with the vibration curve of the impact force over time during the impact process that simulates the process of changes of impact force approximately in the form of a sine function, we obtainwhere is the maximum impact force of the rockfall, is the angular frequency, and is the time.

Since the impulse due to actual impact force is equal to that of the average impact force of rockfall, we obtain equation (87):

Combining the equations (86) and (87), we obtain

Further, we haveand the time duration of rockfall impact is

5. LS-DYA Numerical Simulation Analysis

In order to evaluate the theoretical calculation result and explore the detailed process of rockfall impacting on the soil cushion layer, LS-DYA software was used to analyze the impact process. A typical working condition is taken as an example: size of the soil cushion layer is 6 m (length) × 6 m (width) × 3 m (height); the radius of the rockfall is 0.5 m; and impact velocity of rockfall is 20 m/s. Face-to-face contact is set between the rockfall and the soil cushion layer, and grid refinement is conducted within the range of 1.5 m × 1.5 m in the center of the soil cushion layer. All the nodes in the bottom of the cushion are restricted in three directions: 0.06 s is set as the total time of the analysis step and 200 steps are chosen as the total calculation steps. The calculation parameters are shown in Tables 1 and 2.

Figure 10 shows the results of the LS-DYNA numerical model and the numerical simulation results at 0.02 s, 0.04 s, and 0.06 s. Figure 11 shows the variation curves of impact force and penetration depth of the rockfall with time. After rockfall impact on soil cushion layer, the impact force of the rockfall reaches a maximum value of 1559801 N at 0.004 s (14/200 steps), then reduces gradually, and finally reaches a value of 0 at 0.44 s (147/200 steps). However, the curve shows a stronger vibration between 0.018 s and 0.026 s. The penetration depth of the rockfall gradually increases in the range 0∼0.035 s and reaches a maximum value 0.35 m at 0.035 s (117/200 step), but approximately 9% spring back quantity is generated between 0.035 s and 0.05 s, forming the ultimate penetration depth of 0.312 m.

6. Engineering Example

6.1. Influence of Rockfall Velocity

We have investigated the effects of the collisions of a rockfall block, 0.5 m in radius, with the soil cushion layer when the impact speed of the block is 4, 6, 8, 10, 12, 16, 20, and 24 m/s.

6.1.1. “Impact Velocity-Penetration Depth” Relation Analysis

Figure 12 shows the plot of “impact velocity vs penetration depth.” With an increase in the rockfall impact velocity, the results of these eight methods all reveal a trend of linear growth. The calculation results for the rockfall penetration depth obtained using the cavity expansion, energy conservation, numerical simulation, Australian, and B. S. Guan methods all demonstrated good consistency. Comparatively, the Hertz method has smaller results, while the results from the tunnel manual method and subgrade method are obviously larger.

6.1.2. “Impact Velocity-Contact Radius” Relation Analysis

Figure 9 shows a plot of “impact velocity vs contact radius.” The “contact radius” in this study is defined as follows: when a rock falls on the earth, there is a certain contact area between them. The contact area is calculated and converted into an equivalent circle, and the radius of the circle is the contact radius. From Figure 9, the variation trend for the results of the eight methods is basically consistent with that in Figure 12. In addition, the calculation results for the contact radius based on the cavity expansion, energy conservation, numerical simulation, Australia, and B. S. Guan methods reveal better consistency.

6.1.3. “Impact Velocity vs Impact Force” Relation Analysis

Figure 13 represents the plot of “Impact velocity vs impact force.” With the increase in rockfall impact velocity, the calculation results for the 10 methods all reveal a linear growth trend. The calculation results for the cavity expansion, energy conservation, numerical simulation, Japanese, Swiss, and Australian methods demonstrate better consistency. In contrast, the calculation results for the B. S. Guan, tunnel manual, and subgrade methods are smaller overall, while the calculation result for the Hertz method is slightly larger.

6.1.4. “Penetration Depth-Impact Force” Relation Analysis

Figure 14 shows the curve of “penetration depth vs impact force.” The difference between the curves for the eight methods is evident. The calculation results based on the cavity expansion, energy conservation, numerical simulation, and Australian methods demonstrate better consistency. In comparison, the B. S. Guan method exhibited a slightly smaller impact force with an appropriate penetration depth, the tunnel manual method, and subgrade method exhibited weaker impact forces with larger penetration depth, while the Hertz Method resulted in a greater force with a smaller penetration depth.

6.1.5. “Contact Radius vs Impact Force” Relation Analysis

Figure 15 shows the curve for the “contact radius vs impact force.” According to the figure, the difference between the plots for the eight methods is basically the same as that of Figure 14. The calculation results for the cavity expansion, energy conservation, numerical simulation, and Australian methods demonstrated better consistency and moderate positions.

6.2. Influence of Rockfall Radius

The impact velocity of 20 m/s and rockfall radii of 0.2 m, 0.3 m, 0.4 m, 0.5 m, 0.6 m, 0.7 m, and 0.8 m were considered to analyze the impact effect of rockfall on the soil cushion layer.

6.2.1. “Rockfall Radius vs Penetration Depth” Relation Analysis

Figure 16 shows the curve of “rockfall radius vs penetration depth.” The calculation results for the eight methods all reveal linear growth trends. The calculation results for the penetration depth based on the cavity expansion, energy conservation, numerical simulation, Australian, and B. S. Guan methods demonstrated better consistency. In contrast, the results for the Hertz method are smaller, while that of the tunnel manual method and subgrade method are obviously larger.

6.2.2. “Rockfall Radius vs Impact Force” Relation Analysis

Figure 17 represents the curve of “rockfall radius vs impact force.” According to Figure 17, with the increase in the rockfall radius, the calculation results for the 10 methods all reveal a linear growth trend. The calculation results for the cavity expansion, energy conservation, numerical simulation, Japanese, Swiss, and Australian methods demonstrated better consistency, while the calculation results for the B. S. Guan, tunnel manual, and subgrade methods are all smaller. In addition, the calculation result for the Hertz method is slightly larger.

6.2.3. “Rockfall Radius vs Contact Radius” Relation Analysis

Figure 18 shows the curve of “rolling stone radius vs contact radius.” As shown in this figure, the variation trend of the results for the 8 methods is basically consistent with that shown in Figure 7. Moreover, the calculation results for the contact radius obtained using the cavity expansion, energy conservation, Australian, numerical simulation, and B. S. Guan methods are relatively consistent.

As observed from the preceding analysis, the calculation results for the rockfall impact force and penetration depth based on the cavity expansion, energy conservation, numerical simulation, Japanese, Swiss, and Australian methods reveal better consistency, which explains why the cavity expansion method and energy conservation method have better accuracy. In comparison, the penetration depth obtained for the B. S. Guan method is accurate, but the impact force is slightly smaller. The calculation results for the impact force obtained using the tunnel manual method and subgrade method are smaller, while the calculation result for the penetration depth is larger. In addition, the calculation results for the penetration depth based on the Hertz method are smaller, while that of the impact force is larger.

6.3. Energy Dissipation Analysis
6.3.1. Analysis on Impact Velocity Influence on Energy Dissipation of Rockfall

We analyzed the energy dissipation of the rockfall impact with the rockfall radius of 0.5 m and impact velocity of 4, 6, 8, 10, 12, 16, 20, and 24 m/s.

Figure 19 represents the “influence of impact velocity on energy dissipation”, from which it can be seen that with an increase in the impact velocity, the energy dissipation due to crater formation and due to friction increase rapidly. The increase rate and tendency of the two curves are basically consistent. Furthermore, the gross energy dissipation due to friction is always higher than that due to crater formation, and the energy dissipation of the shed tunnel structure is small overall. Figure 19 also shows the proportion of three kinds of energy dissipations under impact velocity, from which it can be seen that the proportion of energy dissipation due to crater formation gradually rises from approximately 25% to approximately 40% and then becomes stable. Energy dissipation due to friction decreases from approximately 70% to 60% and becomes stable. The energy dissipation of the shed tunnel structure also decreased from approximately 5% to approximately 0%. Therefore, the energy dissipation of the shed tunnel structure could be ignored.

6.3.2. Analysis on the Influence of the Rockfall Radius on Energy Dissipation

Futher, the impact velocity of 20 m/s and rockfall radii of 0.2 m, 0.3 m, 0.4 m, 0.5 m, 0.6 m, 0.7 m, and 0.8 m were considered to analyze the energy dissipation of the rockfall impact.

Figure 20 represents the “Influence of rockfall radius on energy dissipation.” Based on this figure, it can be seen that the energy dissipation due to crater formation and friction increases rapidly with an increase in the rockfall radius. In addition, the increasing tendency of the two curves is consistent. Furthermore, energy dissipation due to friction is always higher than that due to crater formation. Figure 20 also shows the proportion of energy dissipation for different rockfall radii. It can be seen that energy dissipation due to crater formation is maintained at approximately 40% and that friction is approximately 60%, which further indicates that the energy dissipation of a conventional shed tunnel structure can be ignored.

7. Conclusion

In this study, the impact force and penetration depth of rockfall were analyzed theoretically. The main conclusions are as follows:(1)This investigation began with an analysis of the mechanical response mechanism of rockfall impact on the soil cushion layer. Based on the theoretical derivation and the combined cavity expansion model and energy conservation model, respectively, two analytical solutions for the penetration depth and impact forces were obtained. Further, the LS-DYNA software was employed to study the variation law of penetration depth and impact force.(2)For the impact force of the rockfall, either with the change of impact velocity or its radius, the calculation results for the cavity expansion, energy conservation, and numerical simulation methods demonstrated better consistency and coincided well with the results for the Japanese, Swiss, and Australian methods. The results of the Japanese method and the Swiss method were obtained based on laboratory tests. The results for the Australian method were acquired based on field tests.(3)For the penetration depth and contact radius of the rockfall either with the change of impact velocity or its radius, the calculation results for the cavity expansion, energy conservation, numerical simulation, and B. S. Guan methods exhibited better consistency and coincided well with the results for the Australia method.(4)For the relation between the “penetration depth-impact force” and the “contact radius-impact force,” the calculation results for the cavity expansion, energy conservation, numerical simulation, and Australian methods exhibited better consistency. The penetration depth obtained using the Hertz Method was smaller, and the value for the impact force was greater. In comparison, the tunnel manual method and ubgrade method showed inverse results. The penetration depth obtained using the B. S. Guan method is accurate, but the impact force is slightly smaller. Thus, this indicates that the cavity expansion method and energy conservation method have better accuracy.(5)When the impact velocity or radius of the rockfall increases, both the deformation energy consumption and friction energy consumption increase in the form of a parabola. Further, the friction energy consumption is a bit larger than the deformation energy consumption. The impact velocity of the rockfall will influence the energy dissipation ratio for various energies, while the radius of the rockfall has a slight influence. A normal shed hole structure rarely influences the energy dissipation process.

Notations

:Normal stress at point A of soil cushion layer
:Normal stress at point A of spherical cavity
:Radius of the spherical cavity
xA:The x coordinate value of point A
α:The acute angle between the normal line and the vertical line at point A
:The velocity of rockfall impact
:The acceleration of a rockfall impact
:The velocity of spherical cavity expansion
:The acceleration of spherical cavity expansion
:Radius of spherical cavity
:Radius of the cylindrical cavity
:The expansion velocity of spherical cavity
:The expansion acceleration of spherical cavity
θ:The opening angle of spherical cavity microbody
:The opening angle of cylindrical cavity microbody
:The density before deformation of spherical cavity microbody
:The density before deformation of cylindrical cavity microbody
ρ:The density after deformation of spherical cavity microbody
:The density after deformation of cylindrical cavity microbody
r:Distance from the coordinate origin to the mass point of the spherical cavity microbody
:Distance from the coordinate origin to the mass point of the cylindrical cavity microbody
s:The displacement of the mass point of the spherical cavity microbody
:The displacement of the mass point of the cylindrical cavity microbody
:Moving velocity of the mass point of the spherical cavity microbody
t:Moving time of the mass point of the spherical cavity microbody
σr:The radial stress of the spherical cavity microbody
:The radial stress of the cylindrical cavity microbody
σθ:The hoop stress of the spherical cavity microbody
:The hoop stress of the cylindrical cavity microbody
εr:The radial strain of spherical cavity microbody
:The radial strain of cylindrical cavity microbody
εθ:The hoop strain of spherical cavity microbody
:The hoop strain of cylindrical cavity microbody
:The plastic region radius of spherical cavity
:The plastic region radius of cylindrical cavity
υ:Poisson's ratio of soil cushion layer
E:Elasticity modulus of soil cushion layer
τ0:The shear strength of soil cushion layer
τ:The undrained shear strength of soil cushion layer
K:Shear strength coefficient of soil cushion layer
:Radial stress on the spherical cavity surface
:Radial stress on the cylindrical cavity surface
φ:The acute angle between the tangent at the spherical cavity surface through point M with x axis
x:x coordinate value of point on spherical cavity
y:y coordinate value of point on spherical cavity
:x coordinate value of point on the rockfall surface
:y coordinate value of point on the rockfall surface
:The slope of the tangent line L1
:The slope of the tangent line L2
R1:Radius of the horizontal infinitesimal cavity at the upper surface of soil cushion layer
l:The width value of infinitesimal cavity
R:The radius of the rockfall block
D:The diameter of the rockfall block
M:Mass of the rockfall block
L:The penetration depth of the rockfall impact
Lmax:Maximum penetration depth of the rockfall by cavity method
:Maximum penetration depth of the rockfall by energy method
F, F1, F2:Rockfall impact force
:Average impact force of rockfall
:The maximum impact force of the rockfall
:Yield strength of soil cushion layer
:The horizontal movement distance of the thin soil layer
β:The acute angle between the inclined plane of the thin soil layer with the vertical line
dz:Thickness of the soil layer
χ:The acute angle between the inclined plane of the thin soil layer with the horizontal line
:The energy dissipated due to the deformation of thin layer of soil cushion layer
:The total energy dissipated due to the deformation of soil cushion layer
Ff:The friction between the thin layer of soil cushion layer and the rockfall
μ:The friction coefficient between soil cushion layer and the rockfall
:The energy dissipated due to friction between the rockfall and the soil cushion layer
:The total energy dissipated due to friction between the rockfall and the soil cushion layer
:The energy dissipated due to the deflection deformation of the roof of the shed tunnel structure
δf:The roof deflection of the shed tunnel structure
Pf:Given rockfall impact load
E1:Elasticity modulus of the roof of the shed tunnel structure
I:The inertia moment of the roof of the shed tunnel structure
:Span of the shed tunnel structure
ω0:Angular frequency
t1:Time
Q:Quality of rockfall block
ME:The deformation modulus of soil cushion layer
H:Falling height of the rockfall block
:The initial velocity of the rockfall block
λ:Lame constants
:Volume of the rockfall block
Rc:Indentation resistance of target materials
h:The thickness of soil cushion layer
tw:The impact time of the rockfall
c:The reciprocating velocity of compression waves in the soil cushion layer
γ:The specific gravity of the soil cushion layer
ψ:The cross-sectional area of an equivalent sphere for rockfall
ζ:The correction coefficient of the rockfall impact force
κ:The coefficient of rockfall penetration depth
:Gravitational acceleration
a:The acceleration of the rockfall impact
:Inner friction angle of soil cushion layer.

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 there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was made possible by the Sichuan Provincial Traffic Construction Science and Technology Project (2015A1-3) and Fundamental Research Funds for the Central Universities (no. 310826151048). The two sources of financial support are highly appreciated.