 Open Access
 Authors : Ankit Gupta , Shashikant Srivastava ,
 Paper ID : IJERTV9IS060100
 Volume & Issue : Volume 09, Issue 06 (June 2020)
 Published (First Online): 08062020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Seismic Analysis of Underground Structures
Ankit Gupta * Shashikant Srivastava*
Rama University, Kanpur
Abstract: This paper extends to other types of structure the simplified methodology proposed by Constantopoulos et al. (1979) for the seismic design of tunnels. As a practical example, a large structure of reinforced concrete, of box shape and totally embedded in soil, is analyzed. The dynamic pressures acting on walls, roof and floor, due to body and surface waves, are considered in the analyses. A set of seismic load combination hypotheses are proposed to account for the different polarization planes of the seismic waves. The influence of neighbouring buildings can be taken into account considering the new soil stress states that they produce.
Keywords: Seismic waves, Earthquakes forces, Static analysis, dynamic loading

INTRODUCTION
Under seismic conditions, underground structures respond to different seismic waves propagating through soil media. One way to analyze this problem involves three dimensional finite element analyses using appropriate transmitting boundaries, which allow accurate modelling of the structure and the surrounding soil. These full numerical analyses should be carried out in time domain to account properly for soil nonlinear behaviour. However this methodology is expensive when utilized as a design tool of the structure, because it demands long computer time.
Since the pioneering paper of Yeh (1974), in which a simplified analysis for buried pipes subjected to seismic loads was proposed, some attempts have been made to apply those theoretical concepts to the structural design of heavier and more complicated structures, such as underground tunnels and galleries in freefield condition (Constantopoulos et al. (1979 and 1980), Christiano et al. (1983) and Navarro & Samart{n (1988)).
In their methodology, Constantopoulos et al. (1979) assumed that tunnels move in a manner similar to that of the surrounding soil, and thus soiltunnel interaction may be neglected, at least in its traditional sense, because the most tunnel vibration energy is radiated away by soil elastic waves. This method requires only static analyses and accounts very roughly for the difference between tunnel and soil rigidities. They divided the seismic tunnel analysis into two parts: transverse and longitudinal analyses. In the first, a tunnel crosssection is analyzed, considering the soil pressures against tunnel walls, roof and floor caused by the different seismic waves. These pressures are obtained from the stresses in soil in freefield condition, taking into account the changes in the stress distribution around the tunnel originated by its physical presence on soil mass. The types of seismic waves considered in this analysis are: shear and compression waves propagating vertically, and Rayleigh wave
components (distorsional and dilatational) travelling horizontally and 1939
perpendicularly to the longitudinal axis of the tunnel. Love waves are not considered in this part of the analysis because they are polarized in horizontal planes parallel to the ground surface.
For longitudinal analysis, surface waves propagating parallel to the tunnel axis are considered and the tunnel can be modeled as an elastic beam connected by springs to the far field (Winkler model). Navarro and Samartin (1988) provided an analytical solution for this problem. They considered the tunnel bending analyses, as suggested by Constantopoulos et al. (1979), analyzing the effects on the tunnel caused by the distorsional component of Rayleigh and Love waves, but they incorporated a "pushpull" analysis to account for the influence of the dilatational component of Rayleigh waves, improving that suggested for Constantopoulos et al. (1979). To check the modified analysis, Navarro and Samartin (1988) compared the values so obtained with those recorded in actual earthquakes, observing a high degree of accuracy.
In this paper, the methodology proposed by Cons tantopoulos et al. (1979) is extended to other types of buried structures. As a practical example, a reinforced concrete large box (12.512.512.0 meters), protecting diesel tanks, is analyzed seismically. It is founded on intact rock (shear wave velocity v2 = 2,960 m/s and density p2 = 2.71 tons/m3) and the walls and the roof are in contact with a granular backfill (v, = 300m/s and P1 = 2.06 tons/m3). The rock material behaviour was supposed elastic, and the variation of the dynamic properties of the granular backfill versus shear strain were taken into account (Hardin and Dmevich (1972)). A structural transverse section is shown in Figure 1.
The seismic excitation consists of two statistically independent accelerograms with a broad frequency content. They are defined at rock/backfill interface level, with a peak acceleration of 0.2g, in both horizontal and vertical directions, and match the response spectra given in the the Regulatory Guide 1.60. 1.00

SEISMIC ANALYSIS

Preliminary SSI studies
To check the importance of soilstructure interaction (SSI) phenomena in these problems, finite element analyses were carried out, in the frequency domain, using the wellknown numerical tool FLUSH (Lysmer et al. (1975) which takes into account threedimensional effects in an approximate manner. The results confirmed that SSI effects were not very important for the case of tunnels (Constantopoulos et al. (1979) and Navarro (1992)), nor for the problem under study.
1.0 10.50 1.00
1.00
0.50
6.50 BACKFILL
1.50
ENFORCED CONCRE
Figure 1. Geometrical definition of the analyzed structure.
Acceleration (*g)
Acceleration (*g)
Acceleration( g)
Acceleration( g)
Point A Point B
Frequency (Hz) Frequency (Hz)
F'igure 2.Horizontal respouse Spectra at points A and B (continous line) and at its level in freefieldcondition
(dotted line). Damping ratio 5%
In this figure, the horizontal response spectra at structure points A and B (see Figure 1 for situation), corresponding to horizontal excitation, are compared with those obtained in freefield conditions at the same levels, Maximum spectral ordinates appear for the natural frequency of the soil layer. Differences of about 10% are observed: that for point A is greater than that of freefield, whereas for point B the opposite occurs. For vertical excitation these differences are much less.

Seismic environment
The determination of a seismic input for buried 1940 structures is very complex due to the randomness of direction and the magnitude of the seismic motion, so simplifying and conservative hypotheses should be used to account for the uncertainties of the seismic event. The seismic waves considered are: vertically propagating body (shear and compression), and horizontally propagating Rayleigh waves, all polarized in vertical planes, and Love waves polarized in horizontal planes parallel to the ground surface. The last two waves, normally known as surface waves, have a dispersive character when they propagate
through a layered halfspace, ie., wave phase velocity is a function of the frequency of the wave component considered. For the case of Love waves, results of the dispersion equation are shown graphically in Figure 3, for the simple layered halfspace considered. For the fundamental frequency range of an earthquake (15 Hz), Love wave components propagate to a phase velocity very near to the shear waves velocity of the halfspace. A similar conclusion is achieved when treating Rayleigh wves in the same soil system (Ewing et al. (1957)). Fundamental modes of such waves, in the simple layered halfspace, show that no changes in soil stresses are observed through layer depth. This means that in freefield conditions, the maximum soil stresses at layer points, due to surface waves.
The soil stresses caused by vertically propagating body waves may be calculated by means of the computer program SHAKE (Schnabel et al. (1975)) or by using the simplified method suggested by Seed and Idriss (1971) to evaluate soil shear stresses caused by the vertical propagation of shear waves. The author has checked that
the approximate expression given by Seed and Idriss (1971) can also be used for vertically propagating compression waves, considering the same coefficients but introducing the vertical acceleration at ground level.

Seismic soil pressures
Once the soil stresses in freefield, caused by each basic wave type, are known, soil pressure distributions can be computed using a static finite element mesh in which the structure and a part of the surrounding soil are modeled, as suggested by Constantopoulos et al. (1979). In this manner, pressure concentration effects at Structure corners, caused by the structure and soil stiffness differences, can be taken into account. Thus the normal and shear stress distribution, acting on the vertical wall zone and at the middle structure width, due to shear, compression and Rayleigh waves are obtained (Figure 4).
This figure presents some ponits of interest. First of all, the pressure distributions do not follow the Mononobe Okabe theory (Seed and Whitman (1970)) of seismic earth thrust predictions. This is due to the wall flexibility; the theory is not applicable to this type of structure. The shear stress values due to vertical propagation of shear wave are much higher than those obtained from either Rayleigh or compression waves. Note that the first reach a maximum of about 10.5 tons/m2 whereas the second values between 0 and 2
tons/m2. Although the shear stress values would be limited by soilwall friction phenomena, the axial forces at the wall/floor level will be governed by shear wave propagation. The normal stresses acting on the wall due to Rayleigh and compression waves are very similar, and in practice can be assumed as constant along the wall. Those caused by shear wave present a minimum value about the middle of the wall and maximum values.

Influence of neighbouring buildings
When massive building are founded near the considered structure, this methodology does not necessarily lead to conservative design values. Recent works (GÃ©mez MassÃ©6 and Atalla (1984) and Navarro (1992)) on this last subject prove that the influence of buildings can greatly alter the stress states generated in soil mass regarding those obtained in freefield situation.
So, for the case of a building (8,000 tons) founded at surface (Figure 7(a)) on a soil system of the same dynamic properties and seismic input as said before, Figure 8(a) shows the stress state, generated into the soil, at different distances from the building edge and in the middle of the layer depth, for horizontal excitation. For approximately 30 meters from the building edge, the influence of the building has disappeared. For vertical excitation.the influence decreases much faster.

CONCLUSIONS
This paper outlines a simplified methodology for the design of large structures other than tunnels and galleries. It is an extension of that proposed by Constantopoulos et al. (1979).
Some recommendations are made about the choice of Rayleigh and Love wave velocities to be used in the analysis. This greatly simplifies the consideration of the effects of the seismic waves.
Particular attention should be given when massive buildings are near the buried structure because then some of the results of the simplified methodology may not be conservative.
REFERENCES

Constantopoulos, IV. et al. 1979. Dynamic analysis of tunnels. Proc. 3rd. Int. Conf. Num. Meth. Geomech..

Constantopoulos, I.V. et al. 1980. Seismic analysis of buried tunnels. Proc. 7th. World Con. Earthquake Engng.

Christiano, P.P. et al. 1983. Response of buried tunnels to earthquake excitation. Trans. 7th Int. Conf. Struct. Mech. Reactor Tech., Vol. K: 287294.

Ewing, W. M. et al. 1957. Elastic waves in layered media. PrenticeHall.

G6Ã©mezMass6, A. & I. Atalla 1984. Finite element versus simplified methods in the seismic analysis of underground structures. Earthquake Engng. Struct. Dyn. 12: 347367.

Hardin, B.O. & V.P. Drnevich 1972. Shear modulus and damping in soils: design equations and curves. J.

_. Soil Mech. Found. Div. ASCE 98: 603624.

Lysmer, J. et al. 1975. FLUSH, a computer program for approximate 3D analysis of soilstructure interaction problems. Rep. EERC 7530.

Navarro, C. & Samartfn, A. 1988. Simplified longitudinal analysis of buried tunnels. Software for Engng. Workstations 4: 310.

Navarro, C. 1992. Effect of adjoining structures on seismic response of tunnels. Int. J. Num. Anal. Meth. Geomech. (in press).

Newmark, N.M. et al. 1973. Seismic design spectra for Nuclear Power Plants. J. Power Div. ASCE 99.

Schnabel, P.B. et al. 1975. SHAKE A computer pro gram for earthquake response analysis of horizontally latered sites. Rep. EERC 7212.

Seed, H.B. & Idriss 1971. Simplified procedure for evaluating soil liquefaction potential. J. Soil Mech. Found. Div., ASCE 97:12491273.

Sedd, H.B. & Whitman, R.V. 1970. Design of earth retaining structures for dynamic loads. Lateral stresses in the ground & design of earth retaining structures. ASCE.

USNRC. Regulatory Guide 1.60. Design response spectra for seismic design of Nuclear Power Plants.

USNRC. Regulatory Guide 1.92. Combining modal 1944 res ponses and spatial components in seismic response analysis. .

Yeh, G. C. K. 1974. Seismic analysis of slender buried beams. Bul. Seis. Soc. America 64: 15511562. _ Wolf, J.P. 1985. Dynamic soilstructure interaction. PrenticeHall.
