Response spectrum method for spatial seismic ground motion

For seismic analysis of wide structures (for example, airports, stadiums, bridges, etc.), it is important to consider the wave seismic motion of the ground with spatial variability. The differential seismic motion is described by the vector field of kinematic parameters at each point of the base. The differential model of seismic excitation is used for structures on a flexible foundation or discrete supports, if the short-length waves (compared to the size of the foundation) are dominated in seismic waves spectrum. Such excitation is accompanied by angular displacements due to uneven translational ground motions, and it is not yet possible to consider using standard seismic analysis. This paper presents a method of differential seismic analysis with the corresponding modification of the response spectrum (RSM) method. The proposed method requires three accelerograms at each reference point. These accelerograms can be obtained during detailed seismic zoning. To illustrate the method, a model of a stadium cover is considered, for which it is important to apply a spatial differential seismic excitation, and some properties of the excitation are described and analyzed.


Introduction
The data from seismic arrays can show significant spatial variability of the seismic ground motion within small sites [1,2]. The spatial variability can be considered in the earthquake-resistant design on the base of differential models of ground motion, with estimating seismic loads and obtaining seismic response by Response Spectrum Method (RSM) [3][4][5][6].
This paper presents a new modal approach to the seismic analysis of extended structures, when the spatial nature of the seismic excitation is taken into account. It includes a description of spatial differential motion for structural seismic analysis and calculation of seismic response of a structure using the RSM.

Differential seismic excitation
The differential seismic excitation is transferred to a structure as kinematic excitation of the ground base. The ground motion at the support points can be described three translational displacements. For structure with support points the ground motion at the th point , , is described by the vectors of absolute translational displacements, velocities, and acceleration , , , , , , , , , , , , 1, . . . , . The vector of translational accelerations at the th point is equal to: Peak translational accelerations at the th point are equal to = max . For the modal RSM the vector and directional cosines of the seismic excitation are considered as constant: The directions are found from maximum of the dynamic structural response and known as critical directions of the seismic excitation. The directions of the seismic loads applied to the structure are determined solely from its dynamic properties.
Spatial kinematic parameters at the different points is expressed in terms of the same parameters at the anchor point, for example, at the first point with displacements and accelerations . According to Eq. (2), at the first point the accelerations are expressed in terms of the peak translational acceleration = and directional cosines : , , = . For the spatial kinematic field of accelerations, scalar matrix functions of coordinates are introduced. At the th point with coordinates = , , : Then three-dimensional vector of accelerations at the th support point is written via peak acceleration at the anchor point (with regard to = ⋅ ) = = . The 3 dimensional vector of accelerations is equal to: Amplification factors (AFs) can be different at different points of the spatial kinematic field. The ground motion at the th support point is described by three-dimensional vector of spectral AFs Ω = . At the anchor point Ω = , where Ω is the natural frequency.
Functions Θ are introduced for describing of spatial spectral variation of seismic excitation ( -the number of the accelerations vector component , = 1,2,3; -the number of the support point, = 1, . . . , ). At the th -support point with = , , the AFs are equal:

Equations of motion
Consider a structure with degrees of freedom under multidimensional kinematic excitation transferred to the structure through support points. Denote -the vector of absolute generalized displacements, -the vector of relative generalized displacements, -the vector of transport motion and: The transport motion is defined as displacements caused by the displacements of all support points. Equations of absolute motion of the structure are: where × , × and × -the mass, damping and stiffness matrices, × , × the stiffness and damping matrices of support elements, and -the vectors of displacements and velocities of the ground at support points. From the equilibrium of the structure at the support points ( -the vector of transport velocities): Substituting Eq. (4) in Eq. (5) and considering Eq. (6) we get the equation of the relative motion: For the case of small damping + ≫ + , therefore from Eq. (6) + = . Denote the flexibility matrix = , then the transport motion vector is = − , and the equation of relative motion Eq. (7) has the form: Denote the matrix of the transport inertial coefficients × = − , then Eq. (6) takes the ordinary form of the relative motion equation:

Structural seismic response
For solution of Eq. (8) is used modal analysis with orthogonal transformation = , where -the vector of normal coordinates, -the matrix of modes , = 1,2, . . . , . Eq. (8) about normal coordinates is: where is the matrix of squares of natural frequencies, = − -the vector of transport seismic forces. Considering Eq.

Research and testing
To test the method, we perform calculations of many types of structures, both very simple [4,5] and real complex models of structures under various types of earthquakes.
A good example for studying the effect of spatial seismic variability on the dynamic response is a stadium cover model in Fig. 1, which has 8 main supports arranged in a circle with a diameter of 250 m. The range of natural frequencies taken for seismic analysis is 0.13-1.45 Hz, in natural periods 0.7-7.7 s. The ground accelerations at each support are different. To determine the differential ground motion, we have records of two weak natural earthquakes recorded at one of industrial sites [7] at points with distances between them of 150-200 m. A soil is composed of sand and clay layers. A shear waves velocity at near-surface is 230-300 m/s. A seismic event No. 10 is local, short and more high-frequency; an event No. 22 is regional, long and more lowfrequency. In Fig. 2 and 3 are shown the amplification factors (AF) of both earthquakes depending on natural periods for one of horizontal components (direction ).  Table 1 shows the standard deviations (SD) of accelerations for both events for each of the three directions at points 54, 57, and 58. These points are located on a flat site at the vertices of an almost equilateral triangle with sides of about 175 m. Both events demonstrate a significant variability of the acceleration field, more expected for a short-wave event No. 10. In Table 1 the difference between the maximum and minimum SD reaches almost 27 % for event No. 10 and almost 32 % for event No. 22; it indicates the need to apply seismic analysis based on a differential model of ground motion.
Usually, an earthquake engineering of many wide structures is assumed identical ground motion at all support points and using only one accelerogram. Comparing the seismic loads calculated in the traditional way with the loads obtained by the proposed method, we can make brief conclusions: spatial seismic impact leads to a different distribution of seismic loads and internal forces, to a significant increase of seismic loads on some structural elements, to other mechanisms of damage accumulation and destruction.

Conclusions
The paper presents modified RSM that is applied for differential seismic analysis of structures. This method works reliably for random unpredictable wave processes at the base, considering the worst possible situations, for example, the motion of different supports in opposite phases.
Currently, the proposed method is being tested on real large complex objects, a dependence between structural response and variability of spatial seismic ground motion is investigated, new software is developed and put into practice.
The results of this work can be applied in new seismic standards.