Due to the importance of the land displacement caused by shallow tunnels’ drillings in urban areas, many studies have been done on this subject by different researchers, since some decades ago. Pack (1969) did field studies on some excavated tunnels in different soils and showed that ground-level profiles’ subsidence statistical distribution is a normal curve (Gaussian) [1] . Then Atkinson and Potts [2] (1977), Clough and Schmidt [3] (1981), Mill and colleagues (1983) [4] and Leach (1985) [5] offered different relationships to estimate subsidence caused by tunneling operations. Pitilakis and colleagues in 2014 studied the seismic behavior of circular tunnels due to the interaction of ground structures; their main focus was on the tunnel responses, and their priority was tunnel’s behaviors. The results of this study show that in all analyzed cases, increased dynamic response of the tunnel in terms of deformation and lining in shallow tunnels is more significant [6] .

Alielahi and colleagues (2014) studied the use of boundary element method on time-domain for seismic response of buried cylindrical holes in soil. As a result of this research; the numerical results suggested that in a random SV wave with buried hole, horizontal component at the top of the cavity is reduced in comparison with displacement of the empty earth [7] .

Sika and partners examined the surface ground-domain under affections of buried and underground cavities and the ground-level domain was their main effort and activity. The results showed that the numerical methods suggest that the presence of buried craters on the earth has influence on the vacant land’s reinforce and this effect is not negligible [8] . Liu and colleagues in 2013 focused on the seismic behavior of buried tunnels with treating P and SV relay waves. The numerical results show that due to the effects of lined tunnels buried in the ground, linear shear modulus ratio of the average thickness is linear and has significant amounts [9] .

But most studies have been stated on the horizontal surface and few studies have done about the effects of the earth’s tilt on the subsidence which has been caused by tunneling operations. Heydari and his colleagues (2003) expressed that the land slope in the tunnel subsidence is negligible [10] . Ghazavi and Karimpour (2006) examined the slope of the surface subsidence caused by tunneling operations an announced that the land slope is not affected by this issue [11] .

Hasanloo Rad and colleagues in 2014 analyzed the surface slope effects on vertical displacement of tunnel’s crown and floor by ABAQUS 3D Tunnel software help, and mentioned that with increase in surface slope, displacement of tunnel’s crown and floor would have increase and the displacement of the tunnel’s walls at the bottom and top of the slope would encounter decrease [12] . Moosiv and Bayat in 2015 used software ABAQUS 3D Tunnel, and examined the surface slope on the surface subsidence and horizontal displacement of tunnel surrounding and they announced that any increase in ground-level slope will cause subsidence in surface level at the top of the tunnel and horizontal relocation with 6 m from the center of the tunnel at the bottom and top of the slope will be decreased or increased, respectively [13] . Rostami et al. [14] (2016) focus on investigating Effect of Tunnel Gate Shapes on ground settlements and acceleration they concluded that the rectangular tunnels had more effect in case of other tunnel shapes on ground responses. Rostami et al. [15] (2016) evaluated the effect of tunnel maintenances on surface and concluded that an umbrella method had significant effect on surface. Rostami et al. [16] (2016) evaluated effect of tunnels with different shapes on ground acceleration. Mousivan and Askari Ziarati [17] (2014) focused on depth of tunnel in different conditions with plaxis software. Also Alielahi et al. [18] (2015) worked on seismic ground amplification by unlined tunnels and concluded that The effect of the tunnel on the seismic ground surface response is gradually decreased or becomes insignificant with increasing the buried depth of the tunnel. Rostami et al. (2016) investigated Effect of Tunnel Gate Shapes with Similar Cross Section on Inserted Forces on Its Coverage and Soil Surface Settlement and conclude that the forces and settlements on rectangular tunnel shapes more than other tunnel shapes [19] . All the above issues are conducted in static loading state, and despite the expanded studies, evaluating the effects of dynamic loading is less considered. In this study, the effects of frequency and surface slopes on the displacements caused by tunneling under dynamic load will be discussed, using the software ABAQUS.

Earth has a dust layer according to the specifications set in Table 1 [1] . Underground water level is lower than 50 meters from the ground.

In modeling with finite element method, choosing the model’s behavior is very important, because the stress-strain analysis is highly dependent on the behavioral model. To explain the soil behavior, the Hardening Soil Model has been used. The hardening model is a dough model which merges the shear hardening and volume hardening by using the united law to cap submission and using the in-united issue to non-diversion submission [13] (Figure 1).

Main features of this model can be summarized as follows:

・ Stress is dependent to stiffness according to the 3^rd power (hard-tension power dependency)

・ The relationship between the strain and stress contribute to the deviant stress in the tri-axial tests

・ Segregation of initial loading in contrast with of reloading process

・ Failure according to Mohr-Coulomb criteria

・ Material parameters

There is a tunnel in the above land with a diameter of 8.85 m at a distance from the surface and the specifications are mentioned in Table 2 and Table 3. Also the thickness of the tunnel in the form of Lining has been defined in Table 4.

In this research, modeling has been done by using finite element software ABAQUS. Meshing is in average kind with 15-points, and in sensitive areas like the tunnel’s walls finer mesh has been considered (Figure 2 and Figure 3). Dynamic loads applied to the soil mass, are sinusoidal harmonic in range of 100 KN with different frequencies of 2, 4, 8 and 16 Hz and are applied to the bottom of the mass (Figure 4).

In the amplitude above, A is domain, w is angular velocity and f is frequency. The amount of the initial phase (Φ) is considered as zero. Standard boundaries and energy-absorbing barriers have been used in this model. Table 5 shows studies’ programs; in order to analyze the effects of frequency and slopes, 20 models have been considered.

By analyzing the models, land surface subsidence curves changes are as forms 5 to 9 (Figures 5-9).

In static mode with increasing slope, vertical subsidence on top of the tunnel will have increase too, and this subsidence would decrease in the left side of the tunnel and would increase for the right side. In dynamic modes, with slope increase in equal frequencies, vertical subsidence in tunnel’s summit part will increase too, this subsidence would encounter decrease in left side of the tunnel and increase in right side; these changes are due to the decrease or increase of the overloads because of the transformations in surface slopes. In general we can say that the general view of vertical displacement changes in static mode and in dynamic mode is equal. According to the above charts f = 2 in all cases the maximum amount of subsidence is in summit and around the tunnel, this amount of frequency can be considered as dominant frequency and surface slope changes has no effect on the amount of dominant frequency.

By analyzing models, land surface subsidence curves’ changes can be seen as Figures 20-24.

In static mode with increasing slopes, upright subsidence will increase on top of the tunnel and decrease on the left and again increase in the right side of the tunnel. In dynamic state in same frequencies, with increase in the surface slope, subsidence will increase on top of the tunnel and also this subsidence will decrease on the left and increase on the right side of the tunnel; these changes are due to the increase or decrease of the overhead because of the slope changes. In general we can say that the general state of changes of vertical subsidence of surface in both static and dynamic states is equal. According to the above graphs f = 2, in all cases the maximum amount of subsidence is in the top and around the tunnel and this value of frequency can be expressed as the dominant frequency, and changes of surface slope have no effect on the dominant frequency.

In this study, the effects of frequency and slopes on the surface subsidence and horizontal displacement (a distance of 6 km of tunnel X) caused by tunneling under dynamic loads have been considered and summaries of the results are as follows

1) In the static mode vertical subsidence on top of the tunnel will increase with increase in surface slope, and this subsidence would decrease in the left side of the tunnel and increase in the right side.

2) Horizontal displacement of the surface beside tunnel is under affection of surface slope and this surface slope would cause decrease and increase of horizontal displacement for down side of the slope and up side of that, respectively.

3) In dynamic state with increase in the surface slope at same frequency, vertical subsidence on top of the tunnels would increase, and also this subsidence will decrease on the left side and increase on the right side of the tunnel. Horizontal displacement of the ground on the left and right side of the tunnel would decrease in the surface slope up to 5˚ and would increase in angles greater than 5˚.

4) Only one frequency in all cases (except when the slope surface is 5˚) causes the maximum displacement (f = 2) which is called dominant frequency. In the case where the surface slope is 5˚C, the dominant frequency is (f = 4), in other words it can be said that the dominant frequency change is effective on slopes. It is suggested to do complementary studies on the dominant frequency studies.

5) With increasing the depth, the amount of subsidence because of the tunnel excavations would reduce, it can be concluded that the depth is inversely related to the subsidence.

Arash Rostami,Hamid Alielahi,Mehdi Zare,Kaveh Haghighi, (2016) Frequency and Surface Slope’s Effects on the Surface Displacement by Drilling Shallow and Deep Tunnels under Dynamic Loads. Open Journal of Marine Science,06,353-370. doi: 10.4236/ojms.2016.63030