Effect of Geo-Material on Dynamic Response of Tunnel Subjected to Surface Explosion

: Prime materials involved in a problem such as underground structures are concrete, reinforcement steel, and geo-material surrounding the tunnel. Among these three materials, concrete and steel are manufactured materials and their properties can be controlled up to a certain extent. However, geo-material is a naturally occurring material whose constitutive properties vary from region to region, making it highly unpredictable. Findings from one study cannot be applied to other geotechnical problems directly, especially in the case of tunnels subjected to surface explosions. The blast wave generated has to travel through the geo-material before it interacts with the tunnel. As the shock wave propagates radially, its characteristics are likely to be altered by the geo-material. Limited study has been carried out considering this problem. In the present study, the effect of various types of geo-material on the blast response of tunnels subjected to surface explosions is investigated. Finite element analysis has been carried out using LS-DYNA ® , wherein the problem has been modeled using the multi-material arbitrary Lagrangian–Eulerian (MM-ALE) method. Materials with ﬂuid behavior such as air, explosives, and soil are modeled using ALE formulation. Other materials including tunnel lining, reinforcement steel, and rock are modeled using Lagrangian formulation. Blast loading is simulated using the Jones–Wilkins–Lee (JWL) equation of state. Geo-materials considered for the comparative study are sandy loam, saturated clayey soil, sandstone, and granite. Vertical displacement measured at the crown of the tunnel is used to determine the response of the tunnel. Sandy loam soil, being a highly compressible soil, exhibits non-linear and ﬂuid-like behavior under high-strain loading such as explosions. Tunnels undergo extreme deformation in the case of sandy loam soil and clayey soil compared to rock cases. Further, the effect of saturation in sandy loam on tunnel stability is studied. It is observed that with the increase in saturation of soil, more blast energy is transmitted to the structure, which results in higher deformation. Lastly, the effect of the weathering of rock on the tunnel’s response is investigated in the case of sandstone and granite. It was observed that weathering in rock led to more displacement of tunnel crown when compared to intact rock.


Introduction
The prime materials involved in a problem such as underground structures are concrete, reinforcement steel, and geo-material surrounding the tunnel. Among these three materials, concrete and steel are manufactured materials, and the properties which define their behavior can be controlled up to a certain extent. However, geo-material is a naturally occurring material whose constitutive properties vary from region to region, making it highly unpredictable. Findings from one study cannot be applied to other geotechnical problems directly, especially in the case of a tunnel subjected to an external explosion. In the surface explosion, the blast wave generated has to travel through the geo-material before it interacts with the tunnel. As the shock wave propagates radially, its characteristics are likely to be altered by the geo-material. The wave may undergo attenuation or amplification; however, the degree of change depends on many factors such as type of Geotechnics 2022, 2, FOR PEER REVIEW 3 A soil cover of 6 m from the crown of the tunnel is assumed. Considering the symmetry of the problem, the model is cut along two symmetrical planes, and symmetrical boundary conditions are applied to the nodes to constrain their lateral movement. A nominal mesh size of 150 mm is chosen from a convergence study where three sizes were considered: 300 mm, 150 mm, and 100 mm. A 2 m high column of air is modeled on top of the soil to accommodate expulsion of ejecta and formation of craters in the soil. The explosive is modeled using the volume fraction geometry method and is located at the interface of soil and air, over the crown of tunnel. Non-reflecting boundary conditions are applied on the outer sides of the soil domain to reduce the reflection of blast waves, which can cause interferences with the results.
Two cases of geo-material formation have been considered in the study. In the first case, the entire geo-material domain is considered to be made of sandy loam soil. In the second case, the geo-material domain is considered to be made of two layers; the top layer is made of sandy loam, which is overlayed on top of rock strata, as shown in Figure 1. The thickness of topsoil in the second case is taken as 5 m. The entire problem has been modeled using the multi-material arbitrary Lagrangian-Eulerian (MM-ALE) method where tunnel, rock, and rebar are modeled using Lagrangian element formulation. Air, soil, and explosive are modeled using ALE element formulation. Lagrangian and ALE elements are coupled using the fluid-solid coupling method. All the materials are modeled using solid elements with reduced integration except for reinforcement. The rebar steel is modeled using beam elements and they are coupled with concrete using a constrained Lagrangian in the solid option.

Geo-Materials
A soil and foam material model has been used to simulate the highly compressible nature of sandy loam soil. This model is capable of simulating fluid-like behavior under low yield, which makes it suitable for soft soils such as one considered in the present study. The model is based on a pressure-dependent failure surface incorporating volumetric deformation. The deviatoric behavior of the model is defined by a pressure-dependent perfectly plastic yield function, (1) A soil cover of 6 m from the crown of the tunnel is assumed. Considering the symmetry of the problem, the model is cut along two symmetrical planes, and symmetrical boundary conditions are applied to the nodes to constrain their lateral movement. A nominal mesh size of 150 mm is chosen from a convergence study where three sizes were considered: 300 mm, 150 mm, and 100 mm. A 2 m high column of air is modeled on top of the soil to accommodate expulsion of ejecta and formation of craters in the soil. The explosive is modeled using the volume fraction geometry method and is located at the interface of soil and air, over the crown of tunnel. Non-reflecting boundary conditions are applied on the outer sides of the soil domain to reduce the reflection of blast waves, which can cause interferences with the results.
Two cases of geo-material formation have been considered in the study. In the first case, the entire geo-material domain is considered to be made of sandy loam soil. In the second case, the geo-material domain is considered to be made of two layers; the top layer is made of sandy loam, which is overlayed on top of rock strata, as shown in Figure 1. The thickness of topsoil in the second case is taken as 5 m. The entire problem has been modeled using the multi-material arbitrary Lagrangian-Eulerian (MM-ALE) method where tunnel, rock, and rebar are modeled using Lagrangian element formulation. Air, soil, and explosive are modeled using ALE element formulation. Lagrangian and ALE elements are coupled using the fluid-solid coupling method. All the materials are modeled using solid elements with reduced integration except for reinforcement. The rebar steel is modeled using beam elements and they are coupled with concrete using a constrained Lagrangian in the solid option.

Geo-Materials
A soil and foam material model has been used to simulate the highly compressible nature of sandy loam soil. This model is capable of simulating fluid-like behavior under low yield, which makes it suitable for soft soils such as one considered in the present study. The model is based on a pressure-dependent failure surface incorporating volumetric deformation. The deviatoric behavior of the model is defined by a pressure-dependent perfectly plastic yield function, where J 2 is the second deviatoric invariant; a 0 , a 1 , and a 2 are constants for plastic yield function; and p is mean stress. The compaction and crushing of soil are defined using a hydrostatic compression curve which needs to be input in tabulated form. Compressibility of soil decreases as its degree of saturation is increased. Hence, in this study, the hydrostatic compression curve is used as the measure of the degree of saturation in the soil considered.
The properties used to define the behavior of sandy loam soil and saturated clayey soil are given in Table 1. The tabulated data of hydrostatic compressibility are given in Table 2. Three levels of saturation in soil have been considered wherein the degree of saturation increases from sandy loam sample SL1 to sandy loam sample SL3. Sandy loam sample SL2 has a saturation level approximately in the mid-range between SL1 and SL3, and has been taken as a standard sandy loam soil sample for comparison with other geo-materials. Rock has been modeled using the Johnson-Holmquist concrete (JHC-2) model. This model was initially developed to simulate concrete under high-strain-rate loading and large deformation. However, since both concrete and rocks are brittle materials with similar fracture mechanisms, this JHC-2 is also used for rock simulation under extreme loads. It uses a three-stage polynomial equation of state to determine the pressure-volumetric strain of material. The three phases considered here are the elastic phase, plastic phase, and compaction phase. It incorporates a damage failure model which is controlled by defining two damage factors. The yield surface of the model is expressed as: where σ * is the dimensionless equivalent stress and defined as: where σ is the actual equivalent stress; f c is the quasi-static yield strength; A, B, D, N, and C are the material constitutive model parameters; p * is the dimensionless pressure and is the ratio of actual pressure to quasi-static yield strength; and ε * is the ratio of strain rate to reference strain rate. The value of damage factor D defines the damage accumulation as a percentage of the cohesive strength of the material and varies from 0 to 1, wherein D = 1 indicates full damage. Material parameters of all the rocks considered in the study are listed in Table 3.

Air
Air has been modeled using the null material model, where its density is taken as 1.29 kg/m 3 . This model allows incorporation of equation of states without calculating deviatoric stress. Since the model has zero shear stiffness, it requires appropriate hourglass control to avoid energy loss. For the present study, the hourglass coefficient was taken 0.1. The air has been assumed as an ideal gas and its behavior is controlled using a linear polynomial equation of state, which is governed by the following equation: Since the air is assumed as an ideal gas here, the gamma-law equation of state has been applied. Based on this law, the pressure (P) equation can be simplified by substituting the coefficients C 0 , C 1 , C 2 , C 3 , and C 6 with zero, and coefficients C 4 and C 5 with (γ−1). Thus, the expression becomes Here, γ (ratio of specific heat) = 1.4, ρ 0 (initial density of the air) = 1 MPa, and E (internal energy per initial volume) = 0.25 MPa.

Explosive
For the present study, Tri-nitro-toluene (TNT) is taken as the explosive. The behavior of this explosive is defined using high explosive burn material model and the Jones-Wilkins-Lee equation of state. The pressure generated at any given time by an explosive element is Geotechnics 2022, 2

640
where F is a burn factor that controls the chemical energy generated from detonation; P EOS is the pressure estimated from the JWL equation of state; V is the relative volume; and E is internal energy per unit volume. The value of the burn factor, F, depends on the beta setting. There are two burn options: beta burn, wherein detonation occurs due to volumetric compression; and programmed burn, wherein the detonation is controlled by the INITIAL_DETONATION option, which defines the time and location of detonation. The equation of state is not required when the beta option is set to beta burn. For the study, the beta option was set to default, which incorporated both beta burn and programmed burn options. The burn factor, F, was taken as the maximum of the two: where F 1 is for the beta burn option and F 2 is for the programmed burn option, Here, t is the current time and t l is a lighting time which dictates the initialization phase of denotation simulation. For each explosive element, the distance from its center to the detonation point is divided by the detonation velocity, D, which gives the lighting time. V CJ is the Chapman-Jouguet relative volume; υ e is the element volume; and A e,max is the maximum element area. On fulfilment of criteria for the beta option, blast energy is released based on the pressure calculated by the JWL equation of state, which is The parameters required to model the explosive and its detonation are listed in Table 4, where ρ is density; υ D is the velocity of detonation; P cut is the Chapman-Jouguet pressure; and ω, A, B, R 1 and R 2 are JWL parameters. Table 4. Material properties for TNT [6].

Tunnel
The tunnel is a composite made of concrete and reinforcement steel bars. These materials are modeled separately as two parts in order to simulate the elastic-plastic nature of reinforced concrete. Reinforcement steel is modelled using a plastic kinematic model with material properties density = 7850 kg/m 3 , Young's modulus = 210 GPa, Poisson's ratio = 0.29, and yield and stress = 500 MPa. Concrete is modeled using the Winfrith model, which is based on the shear failure surface and is expressed as Strain softening of concrete is simulated in terms of tensile cracking, which is controlled by defining unconfined tensile strength and crack width at which the tensile strength of concrete will be zero across the crack.
The formulation for crack width is based on the results of experiments conducted by Wittmann in 1988 [27], where G F is specific fracture energy and f t ' is unconfined tensile strength. The material parameters are given in Table 5.

Validation of Numerical Scheme
The FE modeling approach developed is validated based on crater size. Baker et al. [28] developed an empirical formula based on a series of experimental tests, which can predict the size of the crater formed in geo-materials due to explosions. It is expressed as Here, W is the explosive charge weight; d is the buried depth of explosive; R is the apparent radius of the crater; σ is the stress factor; and K is the gravitational factor. For validation purposes, a series of numerical simulations were conducted with varying explosive weight and depth of buried explosion. The results were compared with the predicted values of the scaled radius of the crater using the above empirical formula for the corresponding set of explosive charge and depth of burial, as shown in Figure 2. The results agree well with the empirical formula, thus validating the FE modeling scheme.

Validation of Numerical Scheme
The FE modeling approach developed is validated based on crater size. Baker et al. [28] developed an empirical formula based on a series of experimental tests, which can predict the size of the crater formed in geo-materials due to explosions. It is expressed as Here, W is the explosive charge weight; d is the buried depth of explosive; R is the apparent radius of the crater; σ is the stress factor; and K is the gravitational factor. For validation purposes, a series of numerical simulations were conducted with varying explosive weight and depth of buried explosion. The results were compared with the predicted values of the scaled radius of the crater using the above empirical formula for the corresponding set of explosive charge and depth of burial, as shown in Figure 2. The results agree well with the empirical formula, thus validating the FE modeling scheme. Ground motion is measured using displacement, peak particle velocity (PPV), and peak particle acceleration (PPA), whereas ground shock is often measured using soil pressure. Fundamentals of Protective Design for Conventional Weapons (TM5-855-1) [29] is a technical manual published by the United States Army that offers an equation to predict soil pressure (P, in MPa) in sandy loam soil for a given charge weight (W, in kg) and standoff distance of explosives (R, in m), Ground motion is measured using displacement, peak particle velocity (PPV), and peak particle acceleration (PPA), whereas ground shock is often measured using soil pressure. Fundamentals of Protective Design for Conventional Weapons (TM5-855-1) [29] is a technical manual published by the United States Army that offers an equation to predict soil pressure (P, in MPa) in sandy loam soil for a given charge weight (W, in kg) and stand-off distance of explosives (R, in m), Here, f is the coupling factor, which is 1 for a buried charge. To simulate a fully buried explosion, a 50 kg TNT spherical charge buried 5 m below the earth's surface is modeled in a 20 m × 20 m × 20 m soil domain. For the equivalent stand-off distance, R, the soil pressure predicted using the aforementioned empirical equation is plotted with the soil pressure acquired via FE simulation. Figure 3 shows strong agreement between the results, especially at larger stand-off distances.
Geotechnics 2022, 2, FOR PEER REVIEW 9 Figure 3. Comparison of numerically determined soil pressure using the current numerical technique with pressure computed for a buried explosion using the TM5-855-1 empirical equation [29].
Additional validation is performed using ground-motion parameters, displacement, and PPA. Soheyli et al. [30] tested a small-scaled reinforced concrete tunnel with a boxshaped cross section (1 m × 1 m) and a wall thickness of 100 mm. A side-buried explosion of a 1.69 kg explosive (TNT-equivalent) situated 4 m from the tunnel wall was detonated and its effect on the tunnel was measured in terms of displacement and acceleration. This experimental setup is computationally replicated, and displacement and acceleration along the length of the wall facing the explosion are measured. Figure 4a,b shows a comparison between numerically derived and experimental data.

Results and Discussions
The effect of a surface explosion of yield 200 kg TNT on the stability of tunnel is studied for all the geo-materials considered in the study. The crater formed on the topsoil layer made of sandy loam soil can be seen in Figure 5. When the entire geo-material domain consists of sandy loam soil, the apparent depth of the crater formed is 4.3 m. In the second case, the tunnel is surrounded by rock mass with a cover consisting of two differ- Additional validation is performed using ground-motion parameters, displacement, and PPA. Soheyli et al. [30] tested a small-scaled reinforced concrete tunnel with a boxshaped cross section (1 m × 1 m) and a wall thickness of 100 mm. A side-buried explosion of a 1.69 kg explosive (TNT-equivalent) situated 4 m from the tunnel wall was detonated and its effect on the tunnel was measured in terms of displacement and acceleration. This experimental setup is computationally replicated, and displacement and acceleration along the length of the wall facing the explosion are measured. Figure 4a,b shows a comparison between numerically derived and experimental data. Additional validation is performed using ground-motion parameters, displacement, and PPA. Soheyli et al. [30] tested a small-scaled reinforced concrete tunnel with a boxshaped cross section (1 m × 1 m) and a wall thickness of 100 mm. A side-buried explosion of a 1.69 kg explosive (TNT-equivalent) situated 4 m from the tunnel wall was detonated and its effect on the tunnel was measured in terms of displacement and acceleration. This experimental setup is computationally replicated, and displacement and acceleration along the length of the wall facing the explosion are measured. Figure 4a,b shows a comparison between numerically derived and experimental data.

Results and Discussions
The effect of a surface explosion of yield 200 kg TNT on the stability of tunnel is studied for all the geo-materials considered in the study. The crater formed on the topsoil

Results and Discussions
The effect of a surface explosion of yield 200 kg TNT on the stability of tunnel is studied for all the geo-materials considered in the study. The crater formed on the topsoil layer made of sandy loam soil can be seen in Figure 5. When the entire geo-material domain consists of sandy loam soil, the apparent depth of the crater formed is 4.3 m. In the second case, the tunnel is surrounded by rock mass with a cover consisting of two different geo-materials, i.e., topsoil of 5 m thick sandy loam overlaying 1 m thick rock mass. For this case, the apparent depth of the crater is 3.45 m and it is found to be consistent for all types of rock considered. The two predominant causes behind crater formation during the surface explosion is shock-wave generation and rapid gas expansion. Both cause compression and erosion of geo-material in a radially outward direction. Compressive strength and toughness of sandy loam soil are significantly less than that of all rock types. The crater formation in the rock case stays constricted in the topsoil layer and is not extended into the rock layer. Hence, the crater formed in the first case is larger than that of the second case.
Geotechnics 2022, 2, FOR PEER REVIEW 10 and toughness of sandy loam soil are significantly less than that of all rock types. The crater formation in the rock case stays constricted in the topsoil layer and is not extended into the rock layer. Hence, the crater formed in the first case is larger than that of the second case. The blast-wave propagation and crater formation in sandy loam is compared with that of saturated clayey soil. The crater formed in sandy loam soil is wider and the side slope of the crater is flatter in comparison to the crater formed in clayey soil, as evident in Figure 6. The apparent size of the crater in sandy loam soil is relatively big, since the soil has larger air voids and is highly compressible in nature. Further, the shear strength of soil plays an important role in the development of craters as well. An increase in shear strength increases the side slope and depth-to-width ratio of the crater. The shear strength of clayey soil is considerably higher than that of sandy loam soil and, hence, a steeper side slope crater is observed. The clayey soil behaves comparatively stiffer and is less compressible. Since dissipation of blast energy due to compression is small, the tunnel roof suffers full collapse in the case of clayey soil.
It can be seen in Figure 7 that the deflection of tunnel lining is by far the highest for sandy loam soil compared to rocks. If the 1 m thick band of geo-material on top of the tunnel lining is considered as a beam, then the downward displacement measured at the tunnel crown can be taken as the center point deflection of the beam. The downward pressure exerted by the shock wave and gas expansion along with the surcharge of topsoil can be considered as the load acting on the beam, which is consistent for all cases of geo-material. Deflection of a beam is a function of beam stiffness which, in turn, depends on the strength characteristics of the material, especially the modulus strength. The smaller the modulus strength of the material, the higher the deflection of the beam will be. Revisiting the material characteristics of the geo-material, it can be noticed that the shear modulus of all rock types considered is almost 9-10 times greater than the shear modulus of sandy loam soil. Thus, even though the attenuation capacity of sandy loam soil is higher than The blast-wave propagation and crater formation in sandy loam is compared with that of saturated clayey soil. The crater formed in sandy loam soil is wider and the side slope of the crater is flatter in comparison to the crater formed in clayey soil, as evident in Figure 6. The apparent size of the crater in sandy loam soil is relatively big, since the soil has larger air voids and is highly compressible in nature. Further, the shear strength of soil plays an important role in the development of craters as well. An increase in shear strength increases the side slope and depth-to-width ratio of the crater. The shear strength of clayey soil is considerably higher than that of sandy loam soil and, hence, a steeper side slope crater is observed. The clayey soil behaves comparatively stiffer and is less compressible. Since dissipation of blast energy due to compression is small, the tunnel roof suffers full collapse in the case of clayey soil.
It can be seen in Figure 7 that the deflection of tunnel lining is by far the highest for sandy loam soil compared to rocks. If the 1 m thick band of geo-material on top of the tunnel lining is considered as a beam, then the downward displacement measured at the tunnel crown can be taken as the center point deflection of the beam. The downward pressure exerted by the shock wave and gas expansion along with the surcharge of topsoil can be considered as the load acting on the beam, which is consistent for all cases of geomaterial. Deflection of a beam is a function of beam stiffness which, in turn, depends on the strength characteristics of the material, especially the modulus strength. The smaller the modulus strength of the material, the higher the deflection of the beam will be. Revisiting the material characteristics of the geo-material, it can be noticed that the shear modulus of all rock types considered is almost 9-10 times greater than the shear modulus of sandy loam soil. Thus, even though the attenuation capacity of sandy loam soil is higher than that of rock, it fails to mitigate the shock energy in the case of shallowly buried structures.
Geotechnics 2022, 2, FOR PEER REVIEW 11 for both sandstone and granite. Better visualization of the damage sustained by the tunnel due to the surface explosion is presented in Figure 8.  Geotechnics 2022, 2, FOR PEER REVIEW 11 for both sandstone and granite. Better visualization of the damage sustained by the tunnel due to the surface explosion is presented in Figure 8.  Negligible difference is observed between the deflections measured in sandstone and granite cases. The effect of the weathering of rock has also been briefly investigated. Although, compared to sandy loam soil, the deflections for weathered rock cases are much less, an increase in deflection of tunnel crown is observed when compared to intact rock for both sandstone and granite. Better visualization of the damage sustained by the tunnel due to the surface explosion is presented in Figure 8.  The effect of water saturation on blast response of the tunnel has been studied and the displacement-time histories for each case of saturation level are compared in Figure  9. Maximum displacement is observed for case SL3 with the highest water saturation and least compressibility. Soil is a complex heterogeneous material, made up of a skeleton of weathered solid particles with pores filled with air and water. Initial theories assumed that soil particles did not undergo deformation; however, it was later verified that soil solids deformed under loading. Since water is known to be incompressible, the load is shared by the soil grain skeleton and water. Under a slow rate of loading, the water and air get expelled as the soil compresses. These assumptions worked well under static and quasi-static loading. However, under high-strain-rate loading such as blasts, water does not get sufficient time to drain out during the loading period. The soil particles undergo deformations, and water and air get trapped in the soil skeleton. Water, being incompressible in nature, offers resistance against loading and does not allow compression of the soil The effect of water saturation on blast response of the tunnel has been studied and the displacement-time histories for each case of saturation level are compared in Figure 9. Maximum displacement is observed for case SL3 with the highest water saturation and least compressibility. Soil is a complex heterogeneous material, made up of a skeleton of weathered solid particles with pores filled with air and water. Initial theories assumed that soil particles did not undergo deformation; however, it was later verified that soil solids deformed under loading. Since water is known to be incompressible, the load is shared by the soil grain skeleton and water. Under a slow rate of loading, the water and air get expelled as the soil compresses. These assumptions worked well under static and quasi-static loading. However, under high-strain-rate loading such as blasts, water does not get sufficient time to drain out during the loading period. The soil particles undergo deformations, and water and air get trapped in the soil skeleton. Water, being incompressible in nature, offers resistance against loading and does not allow compression of the soil matrix. Thus, when the water saturation in the soil increases, the compressibility of the soil decreases. matrix. Thus, when the water saturation in the soil increases, the compressibility of the soil decreases.
When the soil with a lower saturation level is loaded with a blast wave, compression of soil mass occurs due to the collapse of voids and the crushing of soil solids. The blast energy is partially utilized during compression of soil, which leads to attenuation of blast waves before it reaches the tunnel lining. However, when the saturation level is increased, the voids are filled with water, leading to decreased compressibility and attenuation. Thus, a large amount of impulse is transferred to the tunnel lining, which explains the highest displacement observed for soil case SL3.

Conclusions
The influence of geo-materials on the dynamic response of shallow-buried tunnels under surface explosions is investigated. Four geo-materials have been considered for the comparison: sandy loam soil, clayey soil, sandstone, and granite. An FE model of the problem is developed using the MM-ALE approach with help of software LS-DYNA ® . Two consecutive aspects of surface explosions are studied, which are crater formation and shock-wave generation. The outcomes of the study are presented as follows: 1. The apparent depth of the crater is found to be higher in the case of sandy loam when compared to clayey soil and other rock types. On comparing the shape of craters between sandy loam and clayey soil, it was noticed that while sandy loam has a larger crater with a wider opening, the crater formed in clayey soil had a larger depth-towidth ratio with steeper side slopes; 2. The clayey soil behaved stiffer under the loading and did not undergo much compression compared to sandy loam soil. A large amount of blast energy was transferred to the tunnel structure. The tunnel sustains significant damage in the case of sandy loam soil; however, full collapse of tunnel roof is seen with saturated clayey soil. Sandy loam soil aided in dissipation of blast energy, resulting in a smaller degree of damage on the tunnel. Thus, a soil stratum composed of granular soil is more preferable than fine-grained soil for a loading condition as described in the study; 3. The deflection measured in the case of rocks is very small. The effect of weathering on the dynamic response of tunnels under blast load, although minute in comparison to that of sandy loam soil, is also observed. Higher deflection of tunnel crown is measured in the cases of weathered counterparts of granite and sandstone when compared to the intact rocks. The damage can be minimized by adopting adequate tunnel support systems for tunnels constructed in highly weathered rocks; When the soil with a lower saturation level is loaded with a blast wave, compression of soil mass occurs due to the collapse of voids and the crushing of soil solids. The blast energy is partially utilized during compression of soil, which leads to attenuation of blast waves before it reaches the tunnel lining. However, when the saturation level is increased, the voids are filled with water, leading to decreased compressibility and attenuation. Thus, a large amount of impulse is transferred to the tunnel lining, which explains the highest displacement observed for soil case SL3.

Conclusions
The influence of geo-materials on the dynamic response of shallow-buried tunnels under surface explosions is investigated. Four geo-materials have been considered for the comparison: sandy loam soil, clayey soil, sandstone, and granite. An FE model of the problem is developed using the MM-ALE approach with help of software LS-DYNA ® . Two consecutive aspects of surface explosions are studied, which are crater formation and shock-wave generation. The outcomes of the study are presented as follows: 1.
The apparent depth of the crater is found to be higher in the case of sandy loam when compared to clayey soil and other rock types. On comparing the shape of craters between sandy loam and clayey soil, it was noticed that while sandy loam has a larger crater with a wider opening, the crater formed in clayey soil had a larger depth-to-width ratio with steeper side slopes; 2.
The clayey soil behaved stiffer under the loading and did not undergo much compression compared to sandy loam soil. A large amount of blast energy was transferred to the tunnel structure. The tunnel sustains significant damage in the case of sandy loam soil; however, full collapse of tunnel roof is seen with saturated clayey soil. Sandy loam soil aided in dissipation of blast energy, resulting in a smaller degree of damage on the tunnel. Thus, a soil stratum composed of granular soil is more preferable than fine-grained soil for a loading condition as described in the study; 3.
The deflection measured in the case of rocks is very small. The effect of weathering on the dynamic response of tunnels under blast load, although minute in comparison to that of sandy loam soil, is also observed. Higher deflection of tunnel crown is measured in the cases of weathered counterparts of granite and sandstone when compared to the intact rocks. The damage can be minimized by adopting adequate tunnel support systems for tunnels constructed in highly weathered rocks; 4.
The study concluded with a brief investigation of the effect of saturation. More impulse was transferred to the structure as the degree of saturation in soil was increased. The tunnel roof deflected the most in the case of sandy loam soil sample SL3 with the highest water saturation, implying that the presence of a high degree of moisture in soil can be detrimental to the stability of the tunnel. Prior to construction of tunnels under saturated soil conditions, the water table in the concerned region should be lowered by adopting a suitable dewatering method.
This study can be extended further by considering other types of geo-materials and different influential factors related to geo-materials.