Impact of High Energy Mining-Induced Seismic Shocks from Different Mining Activity Regions on a Multiple-Support Road Viaduct

In this paper, the dynamic responses of a large-scale multiple-support road viaduct to mining-induced seismic events registered in two regions of mining activity were compared. The regions differ in geological structure, which results in discrepancies in the dominant frequency content. Spatial variation of ground motion causing the kinematic excitation non-uniformity was accounted for in the dynamic analyses of this large-scale structure. Non-uniform mining-induced kinematic excitation models were proposed, with respect to the specificity of mining origin quakes. The dynamic performance of the viaduct was determined using three different methods of calculation: the time history analysis, the response spectrum analysis, and the multiple support response spectrum analysis. Both the uniform and non-uniform kinematic excitation models were adopted for the dynamic performance assessment. The research revealed that the dynamic response of some members of the structure, determined using the non-uniform excitation model, was significantly greater than that obtained for the uniform one. Hence, in the dynamic analysis of multiple-support structures under mining-induced events, the effect of spatial variation of ground motion should be considered. The study pointed out that the commonly used response spectrum analysis may lead to the underestimation of the dynamic response of large-scale multiple-support structures. Instead, the multiple support response spectrum method, which takes into account the non-uniformity of ground motion, is recommended as a conservative approximation. This method provides a safe upper estimation of the full-dynamic analysis results of large-scale structures under mining-induced tremors. Finally, the research indicated that the dynamic performance of a structure strongly depends on the frequency range attributed to a specific mining region. The dynamic performance of identical engineering structures under tremors of similar maximal amplitudes may differ significantly due to discrepancies in frequency contents of shocks occurring in various mining regions.


Introduction
The growing demand for fuel and mineral extraction has become a fundamental challenge to the sustainable development of mining clusters. In the next decades, scientists will have to find a balance between the economic benefits and the environmental damage caused by mining. Mining works have a significant impact not only on the environment but on civil surface infrastructure as well. The reduction of the environmental damage caused to land and air resources by mines [1][2][3] must go hand-in-hand with the development of infrastructure protection strategies [4]. The extraction of coal, natural gas, oil, or copper is the main generator of additional static (ground deformation) [5][6][7] and dynamic (shocks of mining origin) loads, which have a substantial impact on engineering structures [8]. In recent years, the evaluation of risk resulting from mining-induced tremors for surface infrastructure became a task of intensive studies, since energies and amplitudes of vibration caused by high-energy mining-induced events are comparable to those resulting from small earthquakes. So far the vast majority of the studies on the influences of mining shocks on surface structures have been focused on residential buildings [9][10][11][12][13]. The recognition of mining-related effects on untypical engineering structures is still insufficient. In particular, the effects of mining shocks on large-scale structures are not well researched.
To assess the impact of shocks triggered by mining activities on engineering objects, dynamic analyses are carried out using various calculation methods. On account of similarities between energies and amplitudes of vibrations caused by high-energy mining-related tremors and small earthquakes, methods applied in seismic engineering can also be utilized for the evaluation of dynamic responses of structures to mining tremors. Due to the relatively simple application and in accordance with international standards and national codes recommendations [14][15][16], scientists and engineers often use simplified dynamic analysis methods, in particular, response spectrum analysis (RSA) [17][18][19]. This approach allows for the estimation of the stress level in structure elements without advanced numerical calculations and leads to correct results for typical infrastructure facilities, such as residential, industrial, or warehouse buildings. However, mining activity areas are often highly urbanized, with many atypical structures, the design process of which requires more advanced calculation tools. The development of computational techniques opens new perspectives in the field of modeling for example large-scale engineering structures.
In the case of the dynamic analysis of large-scale objects, the spatial variation of earthquake ground motion (SVEGM) must be taken into consideration [20]. In such cases, the non-uniformity of kinematic excitation, both in terms of amplitude and frequency, may play a key role in the dynamic response level. Most of the simplified methods like the response spectrum analysis mentioned above, do not take into account this excitation complexity. These methods are based on the assumption that the ground motion under the structure's supports is uniform. Hence, in the case of large-scale structures, with dimensions comparable to the length of the seismic wave, such simplified analysis may lead to significant errors, for example underestimation of the dynamic response to ground motion [21,22].
Typical large-scale objects, sensitive to the above-mentioned excitation non-uniformity, are called multiple-support structures [23]. The dynamic analysis of such objects, taking into account the non-uniform ground motion, is usually carried out using full-time dynamic analysis (time history analysis-THA). The vast majority of such studies address objects localized in areas of seismic activity. The non-uniform kinematic excitations of seismic origin are mostly considered in scientific research analyzing long bridges and footbridges [24][25][26][27], footbridges [28], tunnels [29], dams [30,31], pipelines [32], and large-dimensional cooling towers [33]. There are only a few works on the dynamic response of multiple-support structures, subjected to non-uniform kinematic excitation triggered by mining-induced shocks [34,35].
In addition to THA, there is another method meeting the requirements of the SVEGM effect used for the assessment of the dynamic performance of large-dimensional structures, called the multiple support response spectrum (MSRS) analysis [36]. The MSRS method is based on the RSA, but additionally takes into account the non-uniformity of ground motion. So far MSRS has been utilized mainly in the seismic analysis [37,38], but not in studies concerning the dynamic analysis of structures exposed to non-uniform excitation of mining origin. The lack of such investigations has become a genesis of the present study on the numerical analysis of the dynamic response of the multiple-support structure subjected to mining-related seismic events. Both methods, the full-time dynamic analysis (THA) and the approximate multiple support response spectrum (MSRS), require the assumption of the non-uniform kinematic excitation models. Such models are widely accessible in literature but, being appropriate for natural seismic cases, they do not address characteristics of mining-induced excitations. This notion leads to an urgent need to create models of non-uniform Energies 2020, 13, 4045 3 of 25 kinematic excitation dedicated strictly for mining-related tremors and suitable for utilizing in the above-mentioned methods.
Mining-related seismicity is a significant problem in regions where numerous underground mines are located. For example, an urgent need to protect both existing and newly designed engineering structures against mining-induced shocks exists in two regions of strong mining activity in Poland: the Upper Silesian Coal Basin and the Legnica-Glogow Copper District. The regions differ in geological structure, which results in differences in parameters characterizing mining-triggered tremors, like energy, the frequency contents, and the lengths of shocks.
The main objective of this paper is to assess and compare the dynamic performance of a multiple-support road viaduct subjected to the mining-induced seismic events that occurred in two different regions of mining activity. The main criterion of the selection of these events was the difference in their frequency contents. The dynamic performance of similar structures located in these two particular mining regions may differ in the dynamic response level, due to the dissimilarities of the dominant frequency ranges of both shocks, even if their peak ground accelerations are equal.
To accomplish the above-mentioned objective the acceleration-time histories registered in both regions were applied as the input kinematic excitations. The stress levels in the structural members of the viaduct under the action of both mining tremors were compared. The dynamic performance of the viaduct subjected to the mining-triggered excitations was determined using three different methods of calculation: the time history analysis (THA), the response spectrum analysis (RSA), and the multiple support response spectrum (MSRS) analysis. Both the uniform and non-uniform kinematic excitation models were adopted for the dynamic performance assessment.
Comparison of results of these simultaneous computations allowed for evaluation of the impact of the non-uniformity of excitation on the structure's dynamic response level, as well as assessing the applicability of both simplified approaches (RSA and MSRS) for proper, conservative estimation of stress levels for both excitation models. This was the second purpose of the research.
Investigations on the dynamic performance of large-dimensional structures subjected to spatially varying natural seismic excitation can be found in numerous recent studies. However, the problem of non-uniformity of kinematic excitation is, to the best of the authors' knowledge, rarely if ever addressed in the context of mining-induced events. The novelty of this research lies in:

•
Proposing the models of non-uniform kinematic excitation of mining origin; • Adopting the multiple-support response spectrum method for assessing the dynamic performance of a multiple-support structure under mining-triggered shock; and • Performing comprehensive research on the dynamic performance of a multiple-support structure exposed to non-uniform kinematic excitations of mining origin.
The above-mentioned aspects have not been highlighted in the available literature, and make this study innovative in the fields of civil surface infrastructure protection strategies and sustainable resources extraction development.

Non-Uniformity of Kinematic Excitation of Multiple-Support Structures
In the case of long multiple-support structures, the dimensions of which are comparable with a seismic wavelength, the problem of the spatial variation of earthquake ground motion (SVEGM) becomes an issue [20]. Generally, authors suggest that incorporating the effect of non-uniformity of seismic excitation decreases the dynamic response of a structure, due to a reduction of average amplitudes of ground motion. On the other hand, some authors mention that occurrence of pseudo-static effects, resulting from differences in excitation in particular points of foundations, may lead to an increased response in some structural members [20,28].
In practice, four non-uniformity reasons are most commonly taken into account during dynamic analyses: the wave passage effect, the incoherence effect, the attenuation effect, and the site effect [20,23]. The wave passage effect results from the difference in the arrival times of the wavefront at different locations, so the subsequent ground points localized in the direction of wave propagation repeat the same oscillations with a time delay. The time lag is directly related to the distances between the structure's supports and the seismic shock wave velocity, which strongly depends on the ground stiffness. The effect of incoherence, arising from the scattering of waves in the heterogeneous ground, can be taken into account by means of the coherency function, based on the power and cross-power spectrum density. The function parameters should be defined based on experimental studies. The attenuation effect describes the amplitude decay with increasing distance from the vibration source. This effect strongly depends on the soil conditions in the object localization area and the energy of the earthquake. The relation between the amplitude reduction and the distance from the source point can take different forms for different seismic regions [10,39]. The site effect depends directly on local ground conditions. The complicated system of heterogeneous soil layers results in the modifications of the seismic wave velocity, amplitude, and frequency domain, and in consequence, in various ground motion which causes the complex dynamic response of the object.
In the paper, the problem of non-uniform kinematic excitation was solved using the classical equation of motion for a multi-degrees of freedom system. The classic formula describing the ground excitation can be written as [20]: In Equation (1) M, C, K mean the mass, damping, and stiffness matrices, respectively, with indices s referring to un-supported degrees of freedom, g-supported ones, whereas sg and gs representing the coupling matrices associated with both types of degrees of freedom. X s is the total displacement vector of the structure points. X g and P mean the ground displacement vector and the reaction vector at the supports, respectively.
The total displacement of the system is a sum of the quasi-static X qs s and the dynamic X d s components (Equation (2)). The quasi-static displacements of the non-supported degrees of freedom are obtained from the equation of motion based solely on the stiffness matrix and the supports' motion, with the exclusion of the dynamic effects (Equation (3)): Substituting the quasi-static displacements (Equation (3)) into the equation of motion (Equation (1)) leads to Equation (4), from which the dynamic displacements X d s can be obtained: In the case of kinematic excitation, the term comprising damping components (C sg + C ss ·R) is significantly smaller than the mass damping, and as such can be omitted [15,20]. After the above-mentioned simplification, Equation (4) takes the form: Analyzing Equations (3) and (5) one can notice that the total displacements of a structure depend solely on the ground accelerations. It should be noted, that the vector .. X g contains accelerations of all supports of a structure, and the time histories of ground accelerations are usually registered by a seismological station at one control point only. While analyzing multiple-support large-scale structures the variation of earthquake ground motion should be additionally accounted for. Hence, a model of kinematic excitation has to be implemented. On the basis of accelerations registered at one control point and the adopted kinematic excitation model, accelerations of all supports of a structure can be determined. In the model of non-uniform mining-induced kinematic excitation adopted in the present study, three sources of non-uniformity, discussed earlier in this section, were taken into consideration: the wave passage effect, the incoherence effect, and the attenuation effect. The site effect was neglected in this case due to the high homogeneity of ground conditions at the analyzed location. The first shock analyzed in this study was registered by the seismic station located in the Upper Silesian Coal Basin (USCB). The shock duration was over 6 s, with the strongest intensity phase lasting the first 3 s. The energy of the shock was 1 × 10 7 J, which puts it in the category of high-energy events [10]. The maximum horizontal peak ground acceleration (PGA) of the event reached 0.35 m/s 2 (in the WE direction), whereas in the NS and vertical directions PGA equaled 0.28 m/s 2 and 0.12 m/s 2 , respectively. For the purposes of this study, all registered data was scaled up proportionally, so that PGA in the horizontal direction WE was raised to 1.7 m/s 2 . This value corresponds to the PGA of a mining-induced seismic shock registered in the USCB region in September 2015 [4], classified as a high energy event of 1 × 10 9 J. The acceleration-time histories of the shock in three directions are presented in Figure 1, whereas the frequency spectra are shown in Figure 2.
Energies 2020, 13, x FOR PEER REVIEW 5 of 26 a model of kinematic excitation has to be implemented. On the basis of accelerations registered at one control point and the adopted kinematic excitation model, accelerations of all supports of a structure can be determined. In the model of non-uniform mining-induced kinematic excitation adopted in the present study, three sources of non-uniformity, discussed earlier in this section, were taken into consideration: the wave passage effect, the incoherence effect, and the attenuation effect. The site effect was neglected in this case due to the high homogeneity of ground conditions at the analyzed location. The first shock analyzed in this study was registered by the seismic station located in the Upper Silesian Coal Basin (USCB). The shock duration was over 6 s, with the strongest intensity phase lasting the first 3 s. The energy of the shock was 1 × 10 7 J, which puts it in the category of high-energy events [10]. The maximum horizontal peak ground acceleration (PGA) of the event reached 0.35 m/s 2 (in the WE direction), whereas in the NS and vertical directions PGA equaled 0.28 m/s 2 and 0.12 m/s 2 , respectively. For the purposes of this study, all registered data was scaled up proportionally, so that PGA in the horizontal direction WE was raised to 1.7 m/s 2 . This value corresponds to the PGA of a mining-induced seismic shock registered in the USCB region in September 2015 [4], classified as a high energy event of 1×10 9 J. The acceleration-time histories of the shock in three directions are presented in Figure 1, whereas the frequency spectra are shown in Figure 2.   a model of kinematic excitation has to be implemented. On the basis of accelerations registered at one control point and the adopted kinematic excitation model, accelerations of all supports of a structure can be determined. In the model of non-uniform mining-induced kinematic excitation adopted in the present study, three sources of non-uniformity, discussed earlier in this section, were taken into consideration: the wave passage effect, the incoherence effect, and the attenuation effect. The site effect was neglected in this case due to the high homogeneity of ground conditions at the analyzed location. The first shock analyzed in this study was registered by the seismic station located in the Upper Silesian Coal Basin (USCB). The shock duration was over 6 s, with the strongest intensity phase lasting the first 3 s. The energy of the shock was 1 × 10 7 J, which puts it in the category of high-energy events [10]. The maximum horizontal peak ground acceleration (PGA) of the event reached 0.35 m/s 2 (in the WE direction), whereas in the NS and vertical directions PGA equaled 0.28 m/s 2 and 0.12 m/s 2 , respectively. For the purposes of this study, all registered data was scaled up proportionally, so that PGA in the horizontal direction WE was raised to 1.7 m/s 2 . This value corresponds to the PGA of a mining-induced seismic shock registered in the USCB region in September 2015 [4], classified as a high energy event of 1×10 9 J. The acceleration-time histories of the shock in three directions are presented in Figure 1, whereas the frequency spectra are shown in Figure 2.   the tremor as a high-energy event [10]) was five times larger than that of the shock recorded in the USCB region. Again, for the purposes of this study, all registered data was scaled up proportionally, so that the PGA in the horizontal WE direction was raised to 1.7 m/s 2 . This value corresponds to the maximum PGA of a real mining-induced seismic shock registered in the LGCD region in January 2003 [4]. The components of the shock acceleration and its frequency spectra are presented in Figures 3  and 4, respectively.

Data on the
Energies 2020, 13, x FOR PEER REVIEW 6 of 26 (classifying the tremor as a high-energy event [10]) was five times larger than that of the shock recorded in the USCB region. Again, for the purposes of this study, all registered data was scaled up proportionally, so that the PGA in the horizontal WE direction was raised to 1.7 m/s 2 . This value corresponds to the maximum PGA of a real mining-induced seismic shock registered in the LGCD region in January 2003 [4]. The components of the shock acceleration and its frequency spectra are presented in Figures 3 and 4, respectively.  Scaling up the maximum values of PGA to the same value for both analyzed regions enabled focusing on the discrepancies in the frequency spectra only, and not taking into consideration differences in the shock intensities. It is also worth noticing that the value of 1.7 m/s 2 corresponds to maximal recorded PGA of mining tremors in both regions in the last two decades [4]. The dynamic responses of the bridge to the tremors with such a large value of maximum PGA may be considered as substantial and crucial results for the assessment of the mining tremors impact on the structure.

The Comparison of Frequency Spectra of the Selected Shocks
The shocks' frequency contents differ for various mining regions in Poland. Typically mininginduced shocks are registered with sensors covering the frequency range up to 30 Hz. Dominant frequency ranges of mining tremors registered in the Legnica-Glogow Copper District generally consist of higher frequencies than those registered in the Upper Silesian Coal Basin [10].
The dominant frequencies of the selected shock from the USCB region in all directions are located in the range from 1.6 to 4.8 Hz, with a noticeable peak at 3.5 Hz (see Figure 2). The frequencies of the mining shock registered in the LGCD region show maxima at the dominant frequencies from 5 to 10 Hz, with a visible peak at 7 Hz. Additionally, in the direction WE the amplified acceleration occurs at the frequency range from 17 to 22 Hz with a clear peak at 20 Hz.
Hence, it is clearly noticeable that the tremors from both analyzed mining regions differ significantly, as far as the dominant frequency ranges are concerned. (classifying the tremor as a high-energy event [10]) was five times larger than that of the shock recorded in the USCB region. Again, for the purposes of this study, all registered data was scaled up proportionally, so that the PGA in the horizontal WE direction was raised to 1.7 m/s 2 . This value corresponds to the maximum PGA of a real mining-induced seismic shock registered in the LGCD region in January 2003 [4]. The components of the shock acceleration and its frequency spectra are presented in Figures 3 and 4, respectively.  Scaling up the maximum values of PGA to the same value for both analyzed regions enabled focusing on the discrepancies in the frequency spectra only, and not taking into consideration differences in the shock intensities. It is also worth noticing that the value of 1.7 m/s 2 corresponds to maximal recorded PGA of mining tremors in both regions in the last two decades [4]. The dynamic responses of the bridge to the tremors with such a large value of maximum PGA may be considered as substantial and crucial results for the assessment of the mining tremors impact on the structure.

The Comparison of Frequency Spectra of the Selected Shocks
The shocks' frequency contents differ for various mining regions in Poland. Typically mininginduced shocks are registered with sensors covering the frequency range up to 30 Hz. Dominant frequency ranges of mining tremors registered in the Legnica-Glogow Copper District generally consist of higher frequencies than those registered in the Upper Silesian Coal Basin [10].
The dominant frequencies of the selected shock from the USCB region in all directions are located in the range from 1.6 to 4.8 Hz, with a noticeable peak at 3.5 Hz (see Figure 2). The frequencies of the mining shock registered in the LGCD region show maxima at the dominant frequencies from 5 to 10 Hz, with a visible peak at 7 Hz. Additionally, in the direction WE the amplified acceleration occurs at the frequency range from 17 to 22 Hz with a clear peak at 20 Hz.
Hence, it is clearly noticeable that the tremors from both analyzed mining regions differ significantly, as far as the dominant frequency ranges are concerned. Scaling up the maximum values of PGA to the same value for both analyzed regions enabled focusing on the discrepancies in the frequency spectra only, and not taking into consideration differences in the shock intensities. It is also worth noticing that the value of 1.7 m/s 2 corresponds to maximal recorded PGA of mining tremors in both regions in the last two decades [4]. The dynamic responses of the bridge to the tremors with such a large value of maximum PGA may be considered as substantial and crucial results for the assessment of the mining tremors impact on the structure.

The Comparison of Frequency Spectra of the Selected Shocks
The shocks' frequency contents differ for various mining regions in Poland. Typically mining-induced shocks are registered with sensors covering the frequency range up to 30 Hz. Dominant frequency ranges of mining tremors registered in the Legnica-Glogow Copper District generally consist of higher frequencies than those registered in the Upper Silesian Coal Basin [10].
The dominant frequencies of the selected shock from the USCB region in all directions are located in the range from 1.6 to 4.8 Hz, with a noticeable peak at 3.5 Hz (see Figure 2). The frequencies of the mining shock registered in the LGCD region show maxima at the dominant frequencies from 5 to 10 Hz, with a visible peak at 7 Hz. Additionally, in the direction WE the amplified acceleration occurs at the frequency range from 17 to 22 Hz with a clear peak at 20 Hz.
Hence, it is clearly noticeable that the tremors from both analyzed mining regions differ significantly, as far as the dominant frequency ranges are concerned. The attenuation effect describes the reduction in ground motion intensity with increasing distance from a vibration source. This phenomenon is more rapid in the case of mining shocks compared to seismic events [10]. Hence, the reduction of tremors' amplitudes is non-negligible and may play an important role in the dynamic response of long structures. The attenuation function adopted in this study is based on the empirical dependence of the level of vibration amplitudes on distance from the epicenter, elaborated by seismologists for the analyzed mining regions [10,39]: where a(r) is the acceleration at a distance r; a e is the acceleration in the epicenter zone; r is the distance from the epicenter (km); and H(r) is the empirical function describing the decrease of amplitudes resulting from geometric damping outside the epicenter zone.
The empirical function H(r) takes the following forms: for the Upper Silesian Coal Basin [39]: and for the Legnica-Glogow Copper District [10]: The graphs of H(r) functions for the Upper Silesian Coal Basin and the Legnica-Glogow Copper District are presented in Figure 5. The attenuation effect describes the reduction in ground motion intensity with increasing distance from a vibration source. This phenomenon is more rapid in the case of mining shocks compared to seismic events [10]. Hence, the reduction of tremors' amplitudes is non-negligible and may play an important role in the dynamic response of long structures. The attenuation function adopted in this study is based on the empirical dependence of the level of vibration amplitudes on distance from the epicenter, elaborated by seismologists for the analyzed mining regions [10,39]: where a(r) is the acceleration at a distance r; ae is the acceleration in the epicenter zone; r is the distance from the epicenter (km); and H(r) is the empirical function describing the decrease of amplitudes resulting from geometric damping outside the epicenter zone.
The graphs of H(r) functions for the Upper Silesian Coal Basin and the Legnica-Glogow Copper District are presented in Figure 5.

Calculation Methods of the Dynamic Response of a Structure to Kinematic Excitation
In the present study, the complex dynamic performance of the viaduct subjected to the mininginduced excitation was determined using three different methods of calculation of the dynamic response to kinematic excitation: the time history analysis (THA), the response spectrum analysis (RSA), and the multiple support response spectrum (MSRS) analysis.

Time History Analysis (THA)
The THA algorithm [40] is based on the direct integration of equations of motion at each time increment. Due to its high accuracy, the THA method is suggested in EC standards [14] for both linear and nonlinear dynamic analyses. For the THA algorithm kinematic excitation can be defined by the vector of the supports' accelerations. Both the uniform and non-uniform kinematic excitation models can be adopted for the dynamic performance assessment of a structure using the THA. In this study, the mining shocks recorded in the USCB and LGCD regions (see Figures 1 and 3, respectively) were the sources of accelerations of the viaduct supports.

Calculation Methods of the Dynamic Response of a Structure to Kinematic Excitation
In the present study, the complex dynamic performance of the viaduct subjected to the mining-induced excitation was determined using three different methods of calculation of the dynamic response to kinematic excitation: the time history analysis (THA), the response spectrum analysis (RSA), and the multiple support response spectrum (MSRS) analysis.

Time History Analysis (THA)
The THA algorithm [40] is based on the direct integration of equations of motion at each time increment. Due to its high accuracy, the THA method is suggested in EC standards [14] for both linear and nonlinear dynamic analyses. For the THA algorithm kinematic excitation can be defined by the vector of the supports' accelerations. Both the uniform and non-uniform kinematic excitation models can be adopted for the dynamic performance assessment of a structure using the THA. In this study, the mining shocks recorded in the USCB and LGCD regions (see Figures 1 and 3, respectively) were the sources of accelerations of the viaduct supports.

Response Spectrum Analysis (RSA) with Spectral Curves for the Selected Mining Regions
The RSA is a simplified method of determining the impact of vibration on a structure [40], based on the assumption that any structure movement can be represented by a superposition of its modal responses. Due to this assumption, the application of the superposition principle limits this method to linear problems. The RSA leads to the upper estimation of the structure response (displacements) only, without providing it as a function of time. The method allows for estimating the dynamic behavior of a structure under uniform kinematic excitation only.
The RSA requires adopting spectral curves (response spectrum functions) representing the peak response of a simple harmonic oscillator subjected to a base motion. In the EC standards [14], for seismic analysis purposes, the spectral functions are stochastically determined for typical seismic shocks, taking into account the local soil conditions at a structure's location. However, for the purposes of mining-induced seismicity, response spectrum curves related to different mining regions have been developed on the basis of representative tremors and local soil conditions.
In the present study, the functions of the local response spectra have been adopted for the discussed regions of mining activity. In both cases, the damping ratio of 5% was assumed. The following response spectrum functions [10,41] were applied ( Figure 6): 14] Hz f > 7.14 Hz The RSA is a simplified method of determining the impact of vibration on a structure [40], based on the assumption that any structure movement can be represented by a superposition of its modal responses. Due to this assumption, the application of the superposition principle limits this method to linear problems. The RSA leads to the upper estimation of the structure response (displacements) only, without providing it as a function of time. The method allows for estimating the dynamic behavior of a structure under uniform kinematic excitation only.
The RSA requires adopting spectral curves (response spectrum functions) representing the peak response of a simple harmonic oscillator subjected to a base motion. In the EC standards [14], for seismic analysis purposes, the spectral functions are stochastically determined for typical seismic shocks, taking into account the local soil conditions at a structure's location. However, for the purposes of mining-induced seismicity, response spectrum curves related to different mining regions have been developed on the basis of representative tremors and local soil conditions.
In the present study, the functions of the local response spectra have been adopted for the discussed regions of mining activity. In both cases, the damping ratio of 5% was assumed. The following response spectrum functions [10,41] were applied ( Figure 6): LGCD a It is worth mentioning that the applied response spectrum curves were developed on the basis of many years of research on mining-induced seismicity, and they correlate with the representative tremors and local soil conditions occurring in the USBC and LGCD regions.

Multiple Support Response Spectrum (MSRS) Analysis
The novelty of the present study lies in the implementation of the multiple support response spectrum (MSRS) method for the estimation of the dynamic performance of the viaduct under nonuniform mining-induced excitation. The MSRS method is generally based on the RSA algorithm but, contrary to the classical RSA, the non-uniformity of kinematic excitation may be taken into consideration [36].
The i-th modal response Determined with the MSRS method is the sum of dynamic and quasi- It is worth mentioning that the applied response spectrum curves were developed on the basis of many years of research on mining-induced seismicity, and they correlate with the representative tremors and local soil conditions occurring in the USBC and LGCD regions.

Multiple Support Response Spectrum (MSRS) Analysis
The novelty of the present study lies in the implementation of the multiple support response spectrum (MSRS) method for the estimation of the dynamic performance of the viaduct under Energies 2020, 13, 4045 9 of 25 non-uniform mining-induced excitation. The MSRS method is generally based on the RSA algorithm but, contrary to the classical RSA, the non-uniformity of kinematic excitation may be taken into consideration [36].
The i-th modal response Determined with the MSRS method is the sum of dynamic and quasi-static components, defined by the formula: where a k = q k ·R is the response quantity of interest for unit displacement applied to the k-th support, R = K −1 ss ·K sg , and q k is the influence vector. The parameter b ki can be presented as b ki = q k ·Φ i ·β i , where β i = (Φ T ·M·R)/(Φ T ·M·Φ) and Φ is the matrix of modal shapes. The maximum support displacement is denoted as u gkmax , and D ki is the displacement response spectrum function.
The maximum response of the structure may be calculated using the Complete Quadratic Calculation (CQC) algorithm [36,40]. The total response of the structure can be estimated using Equation (12), or in the simpler version of Equation (13): In the above formulas B = a i ·u gi and Φ BD = b ki ·D ki , whereas l uu = ρ u i u j , l uz = ρ u i s ki , l zz = ρ s ki s lj are the correlation matrices: l uu is the correlation matrix between the displacements at the supports, l uz is the correlation matrix between the supports displacements and modal displacements and l zz is the correlation matrix between the modal displacements. The determination of these matrices plays a key role in the MSRS algorithm. The elements of the correlation matrices are given by: In the above formulae σ u is a mean square value of the ground displacements (Equation (17)) and σ s is a mean square value of the normalized structure's response (Equation (18)): The h function is the frequency response function of the i-th modal equation, while the S u i u j and S ..   u g and the coherency function γ, and can be determined from the following formulae: To conduct the dynamic analysis of a structure using the MSRS method, the algorithm parameters must be defined. The power spectral density function S .. u g can be determined based on the response spectral functions, using formula (22) taken from the literature [36]: where p 0 = 2· ln(2.8·ω·τ/2π), τ is the duration time of excitation and ω f f and Θ are stochastic parameters, most often adopted as: ω f f = 0.705 and Θ = 3.0.

Comparison of the Introduced Calculational Methods of the Dynamic Behavior of Structures
The general comparison of the calculational methods presented in this section can be summarized as follows: • The THA is more precise than methods based on the modal analysis (RSA and MSRS). It allows for finding the time-histories of structure response measures (stresses, strains, displacements) at any element. It also enables incorporating different problem nonlinearities like complex material behavior, large displacements or strains, or unilateral contact conditions. • The methods based on system eigenmodes are much more cost-effective than direct time integration of all the degrees of freedom of the system in the dynamic analysis performed with THA.

•
Taking into account the non-uniformity of ground motion, both in THA or by adopting MSRS method may play an important role in the analysis of the dynamic response of large-scale structures to seismic shocks. However, one should bear in mind that the correct definition of the non-uniform kinematic excitation model requires the knowledge of local soil conditions and in situ tests. Such a model for the current analysis is presented in Section 3.2.

Structural Layout and Material Data of the Viaduct
The analyzed object is an existing road viaduct designated for traffic and pedestrian communication, localized in Chrzanow, Southern Poland. The location is a coal-mining region, so the viaduct can be subjected to mining shocks.
The viaduct is created as an arch bridge (Figure 7). The load-bearing system contains a continuous girder and two arches. The girder consists of two five-span pre-stressed concrete beams integrated with a concrete slab (Figure 8c The viaduct is created as an arch bridge (Figure 7). The load-bearing system contains a continuous girder and two arches. The girder consists of two five-span pre-stressed concrete beams integrated with a concrete slab (Figure 8c), additionally linked by crossbars in the middle span and over supports (Figure 8a). The span lengths are 18.5 m, 23 m, 80 m, 23 m, and 18.5 m. The longest span is suspended to the arches by two sets of 14 hangers located on both sides of the girder in equal spacings of 4 m (see Figure 8b). Each beam is pre-stressed by 15 tendons with an area of 150 mm 2 .   The girder is supported on two concrete abutments and four double pillars acting as intermediate supports, and founded on pile foundations. Due to the deck geometry (the longitudinal and the transverse decline), the height of the pillars differs from 4.2 m on the left side to 5.5 m on the right side. Loads from the superstructure to the supports are transferred by a system of pot bearings. Three kinds of pot bearings are used: fixed, unidirectional, and multidirectional. Over the third support, a pair of fixed and unidirectional bearings is used. Over other supports pairs of unidirectional and multidirectional bearings are applied in a manner ensuring structure stability, yet providing freedom of some displacements to minimize unnecessary forces due to over-constraints.
Concrete classes used for the viaduct's structural elements are C50/60 (arches), C40/50 (girder), and C30/37 (pillars, abutments). Prestressing tendons are made of steel grade Y1860, while steel grade The girder is supported on two concrete abutments and four double pillars acting as intermediate supports, and founded on pile foundations. Due to the deck geometry (the longitudinal and the transverse decline), the height of the pillars differs from 4.2 m on the left side to 5.5 m on the right side. Loads from the superstructure to the supports are transferred by a system of pot bearings. Three kinds of pot bearings are used: fixed, unidirectional, and multidirectional. Over the third support, a pair of fixed and unidirectional bearings is used. Over other supports pairs of unidirectional and multidirectional bearings are applied in a manner ensuring structure stability, yet providing freedom of some displacements to minimize unnecessary forces due to over-constraints.

The Numerical Model of the Viaduct
The finite element model of the viaduct (Figure 9) was created with the ANSYS Workbench software [42]. Numerous types of finite elements were used: eight-node brick elements (for modeling beams, crossbars, arches, and pillars), continuum shell elements (for the slab), and truss elements (for the hangers). The total number of elements in the model amounted to 530,000. The density of the mesh was determined based on convergence analysis, with the natural frequency values serving as the convergence criterion. In order to account for the frictional forces acting in sliding pot bearings, the horizontal springs of constant stiffnesses k = 4e7 N/m were introduced [43]. Due to efficiency reasons rigid structural elements, like foundation footings or abutments, were neglected in the numerical model, since they have the negligible influence on the dynamic characteristics of the object. Moreover, the soil to structure interaction was also neglected in the model. The material parameters applied to the numerical model correspond to the type of planned dynamic analysis. The response spectrum analysis necessitates the linearization of both material and geometric behavior. Therefore linear elastic material characteristics are used in all kinds of dynamic analysis. The pre-compression of the girders is incorporated in the model by the modification of its Young modulus according to Equation (23) [44]: where Ecm/Es are Young moduli of concrete/steel and Ac, As, A are cross-section areas of concrete, steel, and the whole element. The material data used in the numerical model is summarized in Table 1.  The material parameters applied to the numerical model correspond to the type of planned dynamic analysis. The response spectrum analysis necessitates the linearization of both material and geometric behavior. Therefore linear elastic material characteristics are used in all kinds of dynamic analysis. The pre-compression of the girders is incorporated in the model by the modification of its Young modulus according to Equation (23) [44]:

Results
where E cm /E s are Young moduli of concrete/steel and A c , A s , A are cross-section areas of concrete, steel, and the whole element. The material data used in the numerical model is summarized in Table 1.

The Natural Frequencies and Modes of Vibration of the Viaduct
During the first stage of the dynamic analysis, the natural frequencies and modes of vibration of the viaduct were calculated. The numerically evaluated natural frequencies as well as the modes of vibration were verified by the in situ experiment. In Figure 10 the mounting of accelerometers on the viaduct is illustrated. The experimental modal assessments of the viaduct, presented in detail in [45], was based on the operational modal analysis (OMA) techniques. Stochastic subspace identification (SSI) [46] was implemented to obtain natural frequencies of the viaduct and the time domain decomposition (TDD) [47] was utilized for mode shapes evaluation. The comparison of the natural frequencies of the numerical model with those obtained based on the in situ experiment, presented in Table 2, shows high conformity, with discrepancies below 10%. To verify the compliance of the numerical mode shapes of the viaduct (see Figure 11) with the The experimental modal assessments of the viaduct, presented in detail in [45], was based on the operational modal analysis (OMA) techniques. Stochastic subspace identification (SSI) [46] was implemented to obtain natural frequencies of the viaduct and the time domain decomposition (TDD) [47] was utilized for mode shapes evaluation. The comparison of the natural frequencies of the numerical model with those obtained based on the in situ experiment, presented in Table 2, shows high conformity, with discrepancies below 10%. To verify the compliance of the numerical mode shapes of the viaduct (see Figure 11) with the experimental results, the modal assurance criterion (MAC) [48] was used. Over 95% of the obtained MAC matrix diagonal values, related to the same modes (MAC i,i ), were greater than 0.8 (see Table 2), whereas the off-diagonal values (related to different modes) were smaller than 0.2, which indicates a high level of compatibility of the proposed numerical model with the real structure [48]. For the dynamic analysis, the Rayleigh model of mass and stiffness proportional damping was applied [40]. The damping coefficients α = 0.2815 and β = 0.0028 were used. The parameters were calculated based on the damping ratios of 5% and the frequencies corresponding to the first and the second vertical modes of vibration.

The Proposed Model of Non-Uniform Mining-induces Excitation Used for the Time History Analysis
In the present study, the viaduct dynamic performance was numerically assessed using the THA, under the action of both the uniform and non-uniform excitation models. In the case of the uniform excitation model, the time-acceleration histories of the tremor recorded in three directions in the USCB region (see Figure 1) or in the LGCD region (see Figure 3) were applied to all structure supports: WE shock component in the horizontal direction along the viaduct axis, NS component in the transverse horizontal direction, and Z component in the vertical direction. In the case of the model of non-uniform mining-induced excitations, the original records were assigned only to the first support. Signals applied to the consecutive supports were modified, to incorporate the phenomena discussed in Section 2.1 (i.e., the effect of incoherence, the wave passage effect, and the attenuation effect).
In order to account for the local soil conditions, the in situ experiment was carried out by the authors of this study. The ground accelerations under the viaduct were induced by a seismic vibrator (see Figure 12) and recorded at different points of the field. The obtained data was the basis for determining the apparent wave velocity and coefficients of the coherence function. The comprehensive description of the performed experiments, as well as the details of determination of For the dynamic analysis, the Rayleigh model of mass and stiffness proportional damping was applied [40]. The damping coefficients α = 0.2815 and β = 0.0028 were used. The parameters were calculated based on the damping ratios of 5% and the frequencies corresponding to the first and the second vertical modes of vibration.

The Proposed Model of Non-Uniform Mining-Induces Excitation Used for the Time History Analysis
In the present study, the viaduct dynamic performance was numerically assessed using the THA, under the action of both the uniform and non-uniform excitation models. In the case of the uniform excitation model, the time-acceleration histories of the tremor recorded in three directions in the USCB region (see Figure 1) or in the LGCD region (see Figure 3) were applied to all structure supports: WE shock component in the horizontal direction along the viaduct axis, NS component in the transverse horizontal direction, and Z component in the vertical direction. In the case of the model of non-uniform mining-induced excitations, the original records were assigned only to the first support. Signals applied to the consecutive supports were modified, to incorporate the phenomena discussed in Section 2.1 (i.e., the effect of incoherence, the wave passage effect, and the attenuation effect).
In order to account for the local soil conditions, the in situ experiment was carried out by the authors of this study. The ground accelerations under the viaduct were induced by a seismic vibrator (see Figure 12) and recorded at different points of the field. The obtained data was the basis for determining the apparent wave velocity and coefficients of the coherence function. The comprehensive description of the performed experiments, as well as the details of determination of the above-mentioned parameters, are presented in [45].

The Effect of Incoherence
The signals applied to the consecutive supports of the viaduct were firstly modified to incorporate the phenomenon of incoherence, by means of the algorithm described in papers [49,50], in which the authors developed a method of conditional stochastic simulation of seismic wave propagation using the spatiotemporal correlation function in the time domain. The algorithm is used to generate acceleration-time histories at various points of ground motion random field. The process of signal generation is based on the reference event specified at one location, wave velocity, and the adopted coherence function: where σ is the standard deviation of the field of the ground motion; di is the distance between the first and the i-th support; α is the scale parameter depending on local geological and topographical conditions and specifying the degree of correlation between the points of the field; v is the mean apparent seismic wave velocity. The parameters of the coherence function (24) σ = 0.964 and α = 10.64, determined based on the above mentioned in situ test, coincide well with the values most commonly used by other authors [50,51]. The apparent wave velocity in the ground, v = 243 m/s, was determined based on the cross-correlation function of accelerations recorded at points located along the direction of wave propagation.

The Wave Passage Effect
In the model of non-uniform excitation, it is assumed that subsequent points of the ground in the direction of wave propagation repeat the same motions with a time delay dependent on wave velocity. Having the experimentally determined velocity v = 243 m/s and the distances between the supports of the viaduct the time delay was calculated and implemented to modify the accelerationtime histories at subsequent supports of the viaduct.

The Attenuation Effect
The attenuation effect represents the amplitude decay with increasing distance from the vibration source. As mentioned in Section 2.2.4, the attenuation function, Equation (6), was adopted in this study, with empirical formulae elaborated for the discussed mining regions: Equation (7) for USCB and Equation (8) for LGCD. According to these equations, the signals at the consecutive

The Effect of Incoherence
The signals applied to the consecutive supports of the viaduct were firstly modified to incorporate the phenomenon of incoherence, by means of the algorithm described in papers [49,50], in which the authors developed a method of conditional stochastic simulation of seismic wave propagation using the spatiotemporal correlation function in the time domain. The algorithm is used to generate acceleration-time histories at various points of ground motion random field. The process of signal generation is based on the reference event specified at one location, wave velocity, and the adopted coherence function: where σ is the standard deviation of the field of the ground motion; d i is the distance between the first and the i-th support; α is the scale parameter depending on local geological and topographical conditions and specifying the degree of correlation between the points of the field; v is the mean apparent seismic wave velocity. The parameters of the coherence function (24) σ = 0.964 and α = 10.64, determined based on the above mentioned in situ test, coincide well with the values most commonly used by other authors [50,51]. The apparent wave velocity in the ground, v = 243 m/s, was determined based on the cross-correlation function of accelerations recorded at points located along the direction of wave propagation.

The Wave Passage Effect
In the model of non-uniform excitation, it is assumed that subsequent points of the ground in the direction of wave propagation repeat the same motions with a time delay dependent on wave velocity. Having the experimentally determined velocity v = 243 m/s and the distances between the supports of the viaduct the time delay was calculated and implemented to modify the acceleration-time histories at subsequent supports of the viaduct.

The Attenuation Effect
The attenuation effect represents the amplitude decay with increasing distance from the vibration source. As mentioned in Section 2.2.4, the attenuation function, Equation (6), was adopted in this study, with empirical formulae elaborated for the discussed mining regions: Equation (7) for USCB and Equation (8) for LGCD. According to these equations, the signals at the consecutive supports were scaled using the reduction factor, appropriate for the analyzed region (Table 3). Table 3. Reduction factors of the amplitude accelerations for the consecutive supports of the viaduct in the USCB and the LGCB regions. Applying the above-discussed modifications to the recorded tremors (presented in Figures 1  and 3) results in the acceleration-time histories of the subsequent viaduct supports in three directions. For example, the horizontal acceleration component WE for both analyzed mining regions is presented in Figure 13. It is worth pointing out that two analyzed mining activity areas, i.e., the USCB and the LGCD region, differ significantly in geological structure. This results in discrepancies in parameters characterizing mining-induced tremors, like energy, the lengths of shocks, and, especially, the frequency contents.

Supports
Energies 2020, 13, x FOR PEER REVIEW 17 of 26 Table 3. Reduction factors of the amplitude accelerations for the consecutive supports of the viaduct in the USCB and the LGCB regions. Applying the above-discussed modifications to the recorded tremors (presented in Figures 1 and  3) results in the acceleration-time histories of the subsequent viaduct supports in three directions. For example, the horizontal acceleration component WE for both analyzed mining regions is presented in Figure 13. It is worth pointing out that two analyzed mining activity areas, i.e., the USCB and the LGCD region, differ significantly in geological structure. This results in discrepancies in parameters characterizing mining-induced tremors, like energy, the lengths of shocks, and, especially, the frequency contents.

Adoption of the Multiple-Support Response Spectrum (MSRS) Method for the Assessment of the Dynamic Performance of the Viaduct under Mining-Induced Excitations
Since the MSRS method takes into account the non-uniformity effects, the response spectrum functions have to be determined separately for each support of the viaduct. The original local spectral functions developed for both the USCB and the LGCD regions, presented in Section 2.3.2 (Equations (9) and (10)), were applied in three directions at the first support of the viaduct only. The spectral functions applied to the consecutive supports were modified to incorporate the phenomena of incoherence and attenuation. The original spectral functions have been modified using the coherency function with parameters presented in Section 3.2.1 and the relationship described by Equation (22). Additionally, the reduction factors determined for the USCB and the LGCB regions (see Table 3) were applied to the consecutive supports' accelerations, to include the attenuation effect. The modified response spectrum functions applied to all supports in the MSRS method for both mining regions are presented in Figure 14.
Energies 2020, 13, x FOR PEER REVIEW 18 of 26 functions applied to the consecutive supports were modified to incorporate the phenomena of incoherence and attenuation. The original spectral functions have been modified using the coherency function with parameters presented in Section 3.2.1 and the relationship described by Equation (22). Additionally, the reduction factors determined for the USCB and the LGCB regions (see Table 3) were applied to the consecutive supports' accelerations, to include the attenuation effect. The modified response spectrum functions applied to all supports in the MSRS method for both mining regions are presented in Figure 14.

Comparative Analysis of the Dynamic Responses of the Viaduct to Mining-Induced Shocks
The dynamic performance of the viaduct under the mining-induced excitations was determined using different calculational methods: the time history analysis under both uniform and non-uniform excitations, the response spectrum analysis, and the multiple support response spectrum analysis.
Comparison of the results of these simultaneous computations allowed for the evaluation of the impact of the non-uniformity of excitation on the structure's dynamic response level, as well as assessing the applicability of the simplified approaches (RSA and MSRS) for proper, conservative estimation of stress levels for both excitation models.
The dynamic response of the viaduct was determined in all structure members. The response histories were recorded and assessed in over 200 elements. Due to the different predicted behavior of the elements located in the span and the support zones, the detailed analysis is presented for six representative elements placed in those zones ( Figure 15). Elements W1-W3 were situated in the span zones, whereas elements W4-W6 were located above the supports.

Comparative Analysis of the Dynamic Responses of the Viaduct to Mining-Induced Shocks
The dynamic performance of the viaduct under the mining-induced excitations was determined using different calculational methods: the time history analysis under both uniform and non-uniform excitations, the response spectrum analysis, and the multiple support response spectrum analysis.
Comparison of the results of these simultaneous computations allowed for the evaluation of the impact of the non-uniformity of excitation on the structure's dynamic response level, as well as assessing the applicability of the simplified approaches (RSA and MSRS) for proper, conservative estimation of stress levels for both excitation models.
The dynamic response of the viaduct was determined in all structure members. The response histories were recorded and assessed in over 200 elements. Due to the different predicted behavior of the elements located in the span and the support zones, the detailed analysis is presented for six representative elements placed in those zones ( Figure 15). Elements W1-W3 were situated in the span zones, whereas elements W4-W6 were located above the supports. The dynamic responses of the viaduct, in terms of extremal principal stresses (maximal or minimal), calculated at the elements W1-W6 for shocks from both mining activity regions are compared in Figures 16-19. The results of the THA with the uniform and non-uniform excitation models are marked by the black and red solid lines, respectively. The results of the RSA and MSRS analysis are presented by the dashed and solid grey lines, respectively.  The dynamic responses of the viaduct, in terms of extremal principal stresses (maximal or minimal), calculated at the elements W1-W6 for shocks from both mining activity regions are compared in Figures 16-19. The results of the THA with the uniform and non-uniform excitation models are marked by the black and red solid lines, respectively. The results of the RSA and MSRS analysis are presented by the dashed and solid grey lines, respectively.  The comparison of the principal stresses for the USCB tremor, determined by various calculational methods at the span zone elements W1-W3, is shown in Figure 16.
In the span zones, the stresses obtained by the THA under non-uniform excitation are smaller than those determined under uniform excitation. Both approximate methods, RSA and MSRS, lead in this case to a conservative upper estimation of the acceleration-time histories given by the THA. However, the stresses determined with the MSRS are smaller than those acquired using the RSA and closer to the results of the THA. Hence, the MSRS method provides a more accurate solution, better estimating the full-dynamic analysis results.
In Figure 17 the principal stresses at the support zone elements W4-W6 for the USCB tremor are compared. In the vicinity of the supports, the stresses obtained by the THA with non-uniform excitation are up to 15% greater than those obtained with uniform excitation. It can be noticed that under uniform excitation the RSA provides safe (upper) stress estimation of the THA. However, this method underestimates the results obtained by the THA under non-uniform excitation. Although the underestimation is not very significant (about 10%), the results given by the RSA cannot be considered a correct solution to the problem. The comparison of results from the RSA and MSRS shows that the incorporation of non-uniformity effects in the MSRS method results in an increased stress level in the support zones. The key observation is that the dynamic response obtained from the MSRS is greater than that from the THA. This makes the MSRS a conservative method, providing a safe estimation of the full-dynamic analysis results.
The comprehensive analysis of the results obtained for the LGCD region, presented in Figures 18  and 19 leads to quite similar observations.
In the span zones taking into account the tremor non-uniformity leads to a reduction of stress levels obtained from the THA (see Figure 18). The RSA gives a safe upper estimation of the THA results for both excitation models. The MSRS seems to give a conservative estimation of THA calculated stresses only when the excitation nonuniformity is accounted for. However, in some span zone elements, MSRS analysis generates stress levels lower than those obtained from the THA under uniform excitation.
Contrary to the span areas, in the support zones accounting for the shock nonuniformity in THA results in an increase of the stress level up to 30% (see Figure 19). What is important, for elements W4-W6, RSA underestimates the peak stress values obtained by THA under non-uniform excitation by at least 10%. Hence, the RSA method, being nonconservative, cannot be recommended as an approximate method for determining the dynamic response of structures that undergo non-uniform excitations.

Assessment of the Dynamic Performance of the Viaduct Subjected to the Mining-Induced Events from Different Mining Activity Regions
The level of the dynamic response of a structure to a kinematic excitation depends not only on the shock amplitude level, but on the frequencies of excitation as well. The dominant frequency contents are different for various mining regions. In particular, the dominant frequency ranges of mining-induced tremors registered in the Legnica-Glogow Copper District consist of higher frequencies than those registered in the Upper Silesian Coal Basin [10].
To assess the dynamic performance of the viaduct under the shocks registered in both considered mining regions (USCB and LGCD), the maximum stress levels obtained by THA in all analyzed elements for the non-uniform excitation model were collated in Table 4. The comparison of maximal principal stresses at point W4, obtained for the mining-induced non-uniform shocks from both regions, based on THA and MSRS methods is provided in Figure 20. The acceleration-time histories registered in the USCB and the LGCD regions (see Figures 1  and 3, respectively) show the same PGA value, since they were both scaled up to 1.7 m/s 2 . However, the frequency spectra of these tremors differ significantly. The dominant frequencies of the shock from Energies 2020, 13, 4045 22 of 25 the USCB region were located in the range from 1.6 to 4.8 Hz (see Figure 2), whereas the amplitudes of the shock registered in the LGCD region show maxima at frequencies from 5 to 10 Hz (see Figure 4).
It is clearly visible from Table 4 and Figure 20 that the dynamic response of the viaduct to the mining-induced event is 4 to 5 times greater for the shock from the USCB than the LGCD region, despite both having the same PGA. It indicates that the viaduct dynamic performance is strongly dependent on the dominant frequency range of mining tremors attributed to each region. The dominant frequency range of the shock from the USCB region included the first four natural frequencies of the viaduct (see Table 2). This resulted in the resonance effect which tended to increase the dynamic response. One can conclude that the dynamic performance of identical engineering structures located in various mining regions under tremors of similar maximal amplitudes may differ significantly, due to dissimilar shock frequency contents. This assessment was the main objective of this research.

Conclusions
In this study, the dynamic performance of a multiple-support road viaduct subjected to mining-induced kinematic excitations registered in various mining activity regions was assessed. Three methods of calculation were used: the time history analysis (THA), the response spectrum analysis (RSA), and the multiple support response spectrum (MSRS) analysis. The uniform and non-uniform kinematic excitation models were adopted for dynamic performance evaluation. The following conclusions can be formulated based on the analyses performed: • The dynamic response of some parts of the structure, like the support zones, obtained using the non-uniform excitation model was meaningfully greater than that calculated for the uniform excitation. The research proved that some members of a structure may be strongly affected by pseudo-static effects resulting from the non-uniformity of kinematic excitation. Hence, in the assessment of the dynamic response of large multiple-support structures to mining-induced events, the effect of spatial variation of ground motion seems to be vital and should be accounted for.

•
The commonly used RSA method may lead to an underestimation of the dynamic response of a structure to non-uniform mining-induced kinematic excitation. Hence, as nonconservative, the RSA cannot be recommended as an appropriate method in determining the dynamic response of large-scale structures that undergo non-uniform mining-induced excitation.

•
The performed investigation demonstrated that the dynamic response obtained from the multiple support response spectrum analysis was greater than the response received from the THA analyses. The MSRS method, which takes into account the non-uniformity of ground motion, is a conservative approximation and provides a safe upper estimation of the full-dynamic analysis results of large-scale structures under mining-induced tremors.

•
In the dynamic analyses of multiple-support structures subjected to spatially varying ground motion, an adequate model of non-uniform mining-induced kinematic excitation should be applied, taking into account features characterizing events arising in mining zones. In the proposed models the main aspects of ground motion non-uniformity, specific for quakes of mining origin, were accounted for. The parameters of the models, like the velocity of seismic wave propagation, coherence, and attenuation functions, should be obtained experimentally for the local soil condition in mining activity zones.

•
Mining activity regions usually differ in geological structure and topography, which results in discrepancies in frequency contents of mining-induced shocks in these zones. The presented research indicated that the dynamic performance of a structure strongly depends on a frequency range attributed to the specific mining region. The dynamic performance of identical engineering structures located in various mining regions under tremors of similar maximal amplitudes may differ significantly, due to the dissimilarity of frequency contents.
The last conclusion is of great importance for large-scale infrastructure erected in mining activity zones. If the lower natural frequencies of a structure fall into a frequency range typical for shocks Energies 2020, 13, 4045 23 of 25 in a given mining region, the amplification of the structural dynamic response may appear, due to the resonance phenomenon. For this reason, the predictive modal analysis of such structures should be conducted. The analysis methods presented in this research may be useful at the design stage of large-scale multiple-support structures situated in mining activity zones, regarding infrastructure protection strategies and minimizing the risk of the resonance phenomenon appearance.