The Specificity of Dynamic Resistance of Existing Bridge Structures in Mining Areas

The article presents the results of a dynamic resistance analysis of selected existing bridge structures. The numerical calculations performed with the use of THA (time history analysis) pertained to road and rail bridges of steel and concrete structures—not paraseismically designed. The purpose of the work was to determine the dynamic resistance of objects to the influence of paraseismic impacts resulting from mining operations. The specificity of performing calculations in this area of engineering issues is presented. Real mining tremor signals, registered in the vicinity of the objects were used in the numerical analyses. The correct determination of the dynamic resistance of the object can be used as a basis for performing repair/strengthening of existing objects, which in the future will be loaded with additional paraseismic influences.


Introduction
Bridge structures are an important element of infrastructure in urban areas. Through the use of bridge structures, it becomes possible to overcome obstacles, such as rivers or valleys. Often, such constructions are located in areas, where soil accelerations of paraseismic origin occur. In such cases, these existing objects should be subject to an estimation of dynamic resistance if they are not paraseismically designed in advance. Their value should be sufficient to safely transfer additional dynamic influences of, e.g., mining origin. There are a number of buildings erected before mining operations begin, and that is why they are usually not designed with the possible paraseismic impact properly taken into account. Additionally, during use, the bridge may be damaged or its condition will deteriorate over time (e.g., due to the impact of various atmospheric conditions). Hence, it becomes necessary to estimate the current resistance of the real bridge object and compare it with the predicted dynamic impact resulting from human activity in the vicinity of the bridge location.
Intensive human activity in industrial areas, manifested in mining exploitation of deposits, including hard coal, is a source of static and dynamic subsoil impacts on buildings. These impacts can cause damage to objects due to a lack of adequate protection and poor technical condition of these structures. Some regions in Poland, e.g., Upper Silesia, are constantly exposed to mining tremors with peak acceleration values of up to 2 m/s 2 . This value is not significant compared to those recorded during seismic excitations [1]; however, the number of paraseismic impacts recorded during the year is disproportionately higher than seismic. We can observe approximately 1000 paraseismic excitations per year on 40 km 2 area square. The more information could be find on [2].
The problem of design and structure analysis in the field of dynamic impact is complex [3][4][5][6] and still uncertain [7]. It is worth noting that the greatest simplifications usually relate to dynamic load [8]. It is often treated as static without its time-spectral characteristics. Strict dynamic analyses, e.g., using the load history method (THA) are not common. Most often, dynamic resistance assessment is carried out using simplified methods with the use of values characterizing the dynamic properties of the structure, such as frequencies and forms of natural vibrations. The results of numerical analyses can be verified with the results of tests on real objects [9][10][11][12]. The premise for conducting these type of tests, resulting from the nature of the moving load and from the point of view of the randomness of the paraseismic load, may be insufficient. In such a situation, a realistic assessment of the dynamic characteristics of the existing object [13][14][15] should be completed using numerical models.
In terms of construction, bridges are diverse. Thus, it becomes necessary to conduct a detailed analysis of dynamic response to the impact of mining tremors using the response history method with the use of kinematic excitations. These excitations should be a record of the course of soil acceleration change per unit of time in three orthogonal directions. It is important, in this case, that the accelerometer location is in the area of the foundation of the bridge object which is undergoing analysis. This assumption is related to the characteristics of tremors, which depend, among others, on the mining activities being carried out and the parameters of the soil lying in the area, where the bridge is located.
The subject of the described analysis of dynamic resistance to the impact of mining tremors are representative existing objects located in Poland. They were selected from all atypical objects located in the area subject to the impact of the planned exploitation of hard coal by mines in Upper Silesia. The selection was made in such a way as to present the specifics of the procedure for determining dynamic resistance. Not without significance were: the age of the object, the universality of the structure occurrence in the analysed area and its purpose. Selected existing bridges were marked as: A-reinforced concrete road viaduct about 45 m long, over a railway line.
B-steel railway viaduct over a river, about 27 m long. C-reinforced concrete road viaduct about 56 m long, over a railway line.
D-steel railway viaduct over a road, about 56 m long.

Data on Mining Tremors Adopted in Dynamic Analyses
One of the first important stages of the analysis of dynamic resistance of existing bridges is the collection of information on the paraseismic activity of the mining area, where the object is located. This is related to the conditions of the selected method of analysis of the load-bearing structure resistance. Due to the unusual nature of this analysis, it was first and foremost necessary to obtain accelerograms registered at measuring stations existing in the area of the location of the analysed objects. Obtained accelerograms came from two sites located on the surface of the land. These places are located in the basements of a residential building and an outbuilding. The system of recording vibrations generated by mining tremors has been operating since August 2010. It is also known that in the previous 10 years there were no tremors with energy greater than 1.0 MJ.
The accelerograms contained data of eight actual tremors with a time step dt = 0.0026 s, recorded in three orthogonal directions. While analysing the received signals, attention was paid to: maximum horizontal acceleration, resultant acceleration at the measurement site after filtration in the 0-10 Hz band as well as distance between the object and the epicenter. The most intense recorded tremors from the entire period of monitoring the ground surface vibration were selected (Table 1), i.e., from August 2010. The data does not include measurements, as well as estimates of peak accelerations in the epicentral zone.
Generally, when assessing the usefulness of the compiled data for the numerical analysis of objects, it should be stated that the tremors are of low intensity, and therefore, only three out of eight were selected for further calculations. These were designated as Z107, Z137 and P6 and with the maximum actual horizontal acceleration of 201, 148 and 155 mm/s 2 . Graphical presentation of time-varying acceleration has been compiled in Figs. 1-2 which correspond to the order number of tremors.
The listed accelerograms (Figs. 1, 2, 3) were selected from those available as kinematic excitations of the numerical model. Adoption of the optimal signal will finally make it possible to achieve the lowest dynamic resistance of the analysed object, and thus ensure the safety of the structure. The optimal signal should be considered the one which has the highest intensity and time-spectral properties that harmonize with the modal characteristics of the analysed object. Fulfilling these demands guarantees an extreme dynamic response of the object. From a practical point of view, the most important element of the analysis is the correct selection of signals, due to the time-spectral properties that resonate with the viaduct model. The use of signals recorded in the area, where the analysed object is located makes the occurrence of tremors with similar characteristics more realistic. In turn, this makes the procedure of determining dynamic resistance more realistic.
Each signal was subject to an independent assessment as to its usefulness due to: • acceleration peak values, • duration, • spectral characteristics using the Gabor transform.
In general, the Gabor transformation [16] can be defined as the representation of a non-stationary signal in the form of the sum of base functions. They arise from the assumed prototype function as a result of its movement along the time and frequency axis, by multiplying complex harmonic signals. The time-frequency Gabor transformation of the continuous signal x(t) is defined by the formula: where g m;n t ð Þ ¼ g t À m Á Dt ð Þ e j2mDf Át ; Dt Á Df 1-shifting in time domain by m Á Dt and frequency domain by n Á Df the so-called function of the synthesis window g(t) with energy equal to one, Dt, Df-preset shift in time and frequency domain, respectively.
The synthesis window should be characterized by a concentration of energy in both domains. Practically, this means that it is compact, i.e., has a narrow base. This ensures good resolution of the signal synthesis and the accuracy of results. Represented in the form of changes in the spectral and time properties of the tested signal (Fig. 3). Usually a Gaussian window is used, because it can be narrowed and expanded depending on the needs and computational conditions, primarily from the period of discretization of the analysed signal.
Gabor transformation decomposition coefficients are calculated from the signal using the so-called analysis window function c(t): This function (2) is associated with the synthesis function (1), because it is biorthogonal to it (complex conjugate function). This means that it is simultaneously orthogonal in time and frequency domain, which is expressed by the formula: X m;n g m;n t ð Þc Ã m;n t 0 When analysing one signal, the same analysis window is used, which moves along the signal in the time domain. In this way, the values of the coefficients are obtained in different points on the plane defined by the time and frequency axes.
The Gabor transform modules, i.e., the resulting timefrequency representations of Gabor, i.e., the equivalents of spectral power density, but calculated simultaneously in both time and frequency domains, are determined by the formula: By analogy with the Fourier transform, function (4) can be interpreted as spectral power density.   (1) and (2) also have a discrete version. It was implemented using a script written in the Matlab environment. Chebyshev polynomials showing energy concentration around the central frequency and minor errors related to the so-called side leaves adopted as a synthesis window. An example of a solution to this task is a map in the frequency and time structure, which is presented in Fig. 4. The dominant frequencies in the Z107 signal in the Y direction (placed on the map as an overview) are marked in a dark colour. Significant signal frequencies in this direction are located between 0.75 and 1.00 s of their duration and take values between 4 and 8 Hz. Similar analyses were carried out for all signals included in the study.
According to the current prognosis regarding the impact of mining tremors on surface objects, acceleration values from 110 to 230 mm/s 2 were predicted (by the mining plant) depending on the location of the object and the mining area. In this study, 270 mm/s 2 is assumed as the extreme boundary acceleration forecasted for the analysed objects. This approach allows to safely determine the resistance of an object assuming possible differences in the prognosis.

Description of the Road Viaduct-A
The viaduct (Fig. 5) was built in 1952. It is located in the area of Upper Silesia. It enables passage over a twotrack railway line. The building has a monolithic reinforced concrete supporting structure in the form of a three-span frame with theoretical spans of 13.75-18.55-13.40 m. The viaduct crosses the railway tracks at an angle of 79.5°. This is also the slant of the span and supports relative to the span. The superstructure of the viaduct platform is a reinforced concrete slab, based on five continuous support beams. Beam spacing is 1.44 m. Their cross-section is 1.05 m high and 0.55 m wide, which is increased in the support area on the frame pillars. On both sides of the road, cantilever pavement slabs, with a projection of 1.80 m beyond the face of the extreme support beams, were made. The platform width is 9.92 m and its length is 46.35 m. The viaduct is supported by two monolithic frame pillars with five posts each, founded on a foundation footing. Abutments were also made as monolithic, full-walled, with wings located parallel to the road. The beams were based on abutments via slip layers. The facility was not designed with mining loads taken into account. In addition, mining tremors were also not considered. Currently, the object has undergone renovation together with the strengthening of support columns. The article analyses the object's prestrengthening resistance.

Technical Description of the Railway Viaduct-B
The facility is located on the railway line above the stream and road. This line has two independent tracks, and therefore, the object consists of two independent viaducts: the older one, erected in 1946 under one track, located on the west side, the newer one from 1973, under the second track, located on the east side.
Each of the viaducts has two spans and is structurally independent of the other (Fig. 6). The choice of the older bridge structure was dictated by the greater span length and   Viaduct girders were based on a pillar with fixed hull bearings. On the abutments there are single-direction sliding bearings, of a double-shaft design. The eastern viaduct has lateral expansion joints on the pillar and abutments.
The facility was not designed with taking possible paraseismic loads into account.

Description of the Road Viaduct-C
The supporting structure of the road viaduct over the railway line consists of three freely supported spans with a grid and plate structure. The theoretical lengths of the spans are 14.10 m, 17.00 m and 13.75 m. The total length of the supporting structure is 46.9 m. The total length of the viaduct, including the abutments is 55.95 m. Abutments and pillars were made as reinforced concrete, monolithic. The grate of each span has four load-bearing beams with a cross section of 0.

Technical Description of the Railway Viaduct-D
The steel railway viaduct over the two-lane road was constructed as a single-span structure with a theoretical span of 22.7 m. The total length of the bridge together with the abutments is 56.5 m. There is one railway track on the bridge structure, the axis of which crosses the support axis at the abutment at an angle 49°. In the cross section, the bridge is a system of two I-ribbed girders with a height of 2.0 m each. The upper and lower flange of the welded plate girder are 600 9 25 mm rectangles in cross section. A steel web is located between them with a cross section of 15 9 1950 mm, which local stability is ensured by a system of welded vertical ribs spaced every 1.5 m. A steel beam pavement structure with a 1.6 m reach has been attached to the outside of the main girders. There is a reinforced concrete abutment with transverse dimensions of 5.4 9 2.0 m at the bearing location and 7.4 9 2.0 m at direct foundation. There are fixed bearings on one of the abutments and one multi-directional sliding on the other. This is where the supporting system in the form of two plate girders is based.

Assumptions and Data for Models and Numerical Calculations of Selected Bridges
The analysis of this dynamic resistance was carried out in a unique computational situation assuming that the paraseismic load is distinct. The following assumptions resulted from the analysis: • analogy between unique computational situations with paraseismic and seismic loading, • taking into account the almost constant variable load value in the combination, • forgoing the effects of wind and snow on the analysed object.
The unique combination according to Eurocode accepts in the analysed case: where G k,j -unfavourable values of permanent loads on the structure, w 2,i -nearly constant value variable load factor, Q k,i -characteristic value of the variable load, A d -the calculated value of the unique load generated by mining tremors.
The basic premise for using in the analyses only of the values of almost constant variable (utility) loads results from the fact that the probability of paraseismic load of extreme intensity, i.e., from regional tremor in the mining area, coinciding with the maximum possible load of the viaduct resulting from road or rail traffic is very low. The estimated period of occurrence of a significant (causing significant response of the object) excitation is up to several years. In addition, the duration of this interaction with extreme intensity is no more than a few seconds (which is confirmed by Figs. 1-2), and the accumulation time of extremely intense car or rail traffic on the viaduct span certainly is not longer than a few minutes. This fact results from field observations of traffic on the analysed objects. This approach seems to be fully justified, and therefore, it has been implemented. As a result, the coefficients of almost constant values of major mobile interactions were adopted as w 2,2 = 0.
Permanent, variable and movable loads were taken from the bridge design and supplemented with standard values. The masses distributed in each model were calculated in accordance with the principles of conducting paraseismic analyses, i.e., they result from the characteristic values of permanent object loads and the fact that the almost constant (long-term) part of the moving load is zero.
Considering the above, it was assumed that independent model calculations were necessary: 1. Static analysis including permanent and variable loads (including moving loads) in the basic combination. 2. Modal analysis taking into account the real rigidity of structural elements and masses attributable to these elements. 3. Dynamic analysis (THA) allowing for real signals excluding constant, variable and moving loads.
In the numerical analyses, spatial structural calculation models were used by means of the Finite Element Method (FEM) implemented in ABAQUS software. Each bridge model reflects the geometry and the load-bearing system using the structural elements listed in the technical description. The elastic work of materials (steel and concrete) was adopted in all models. The boundary conditions at the point, where the supporting elements were supported on the abutment or pillar took into account the actual bearings used (freedom of movement was defined). The presented cases have the abutment supported on the foundations piles. The dynamic signal is applied to the supports as kinematic displacement over time. This simplification causes not consider the ground stiffness and damping.

Model and Data for Numerical Calculations-A
The spatial numerical model used in the calculations was based on the FEM method. It consisted of 14,000 shell and bar finite elements with a side dimension of approx. 0.2 m. It was assumed that the columns in the foundation footing were fixed and the beams rested on the abutments on the contact layer with a friction coefficient of 0.5. The model view is shown in Fig. 7. Permanent, variable and movable loads were taken from the bridge design and supplemented with standard values. The masses distributed in the model were calculated in accordance with the principles of conducting paraseismic analyses: • the distributed mass of the main beam with the associated part of the slab and topcoats is 13 460 kg/lm, • the distributed pavement mass with surface is 1 820 kg/ lm.
Data characterizing of the construction materials was taken from the building design: The values of internal quantities in the critical crosssections of the main beam from self-load and useful load were taken from the design calculations. In addition, the design load-bearing capacity of these cross-sections due to bending and shear calculated according to standards was also taken from the design calculations. They are associated with the following loads on the main beams in static analysis:

Model and Data for Numerical Calculations-B
In the numerical calculations, the spatial viaduct calculation model using FEM was used. It included structural elements of the object, i.e., main beams, transverse and wind bracing as well as stiffening ribs in the scope of elastic work of the material. The model uses square shell finite elements (side dimension 0.15 m) in plate girder and bar in concentrations. The total number of finite elements (shell and bar) was over 5000. The view of the model is shown in Fig. 8. The calculations took into account the viaduct loads, among others the distributed mass of the main beam with the associated part of secondary elements and track structure as 1550 kg/lm.
The following loads were applied to the main beams of the object:  The view of the model is shown in Fig. 9.

Model and Data for Numerical Calculations-D
The spatial calculation model of the steel viaduct was used in the numerical calculations. Bar and shell finite elements were used, mainly quadrangular with dimensions of about 0.25 9 0.25 m. This provides for a total amount of about 7 thousand elements. The view of the model is shown in Fig. 10. In the numerical model, in addition to the masses of load-bearing elements, substitute masses were introduced from service walkways, cantilever extensions of the main girders and tracks with equipment and the track bed. The masses have the following values for: • pavements-4 kN/m, they are arranged in the plane of the webs of the main girders along the upper ledge, differences in values resulting from the existence of steel pipelines made under the pavement are also taken into account, • railway tracks-1.06 kN/m2, which is evenly distributed on the surface of the viaduct orthotropic plate, it includes tracks, fenders, wooden sleepers with washers, a layer of crushed stone and shotcrete.

Presentation of the Results of Numerical Calculations
Calculations for each model object were carried out in stages. In the first (static) stage, the created model was subjected to constant and variable loads, including mobile interaction, in the basic combination. It is possible to omit this step in the dynamic resistance analysis using the available viaduct design calculations. According to the authors, this approach may generate differences in the level of determined resistance due to the fact that different calculation models are used. Thus, it is more advantageous to use the ''own'' model, and the results of static analyses obtained in this way can be referred to those contained in the viaduct design project. The second stage was a modal analysis taking into account the real rigidity of structural elements and masses attributable to these elements. The common procedure at this stage is the conversion of constant loads to masses distributed over structural elements. Performing the aforementioned postulates makes it possible to draw conclusions from the results of viaduct analyses regarding the level of projected effort (safety margin), and also, what is important, the optimal selection of the soil acceleration signal. The last stage was to conduct a full dynamic analysis (THA) with a kinematic load dependent on the actual signal. At this stage, constant, variable and movable loads were omitted. Numerical dynamic calculations of viaduct models were carried out in accordance with the following assumptions: • the actual stiffnesses of the structural elements together with the masses attributable to them (e.g., road surface) were taken into account, • basic forms of natural vibrations and corresponding frequency values were sought, • material damping with the first natural frequency were taken into account, • the models were only subjected to dynamic loads, excluding constant, variable and movable loads, • in the numerical calculations, time analysis (THA) was used, i.e., integration of the equation of motion by the Newmark method [14,17,18], • kinematic excitations were applied, i.e., real accelerograms recorded for tremors, respectively, Z107, Z137 and P6, • independent cases of excitation loads were simultaneously horizontal components acting on the X and Y directions, then only the vertical component and the combined action of the three components, • on the basis of dynamic fields of normal and cutting stresses in the critical cross-sections of the viaduct's structural elements, the dynamic cross-sectional sizes and support reactions were calculated, • finally, from the obtained results, extreme cases were selected, which formed the basis for estimating the dynamic resistance of the object.
In the following table, the results of the dynamic analysis of the A and B bridge models are described in detail, in the case of the other models (C and D) the algorithm was identical, which is why the brief versions of the results are presented.

Analysis of the Characteristics and Dynamic
Response of the Model-A The results of the first stage, i.e., static analysis, are the values of cross-sectional numbers resulting from the bridge design calculations, without taking into account paraseismic influences. They are summarized in Table 2. It takes into account extreme bending moments and lateral forces in the critical cross-sections of the bridge, i.e., the support zone (above the columns) and the span (in the middle of the span). These results correlate with the values included in the building design of the viaduct A.
The numerical solution of the model in the field of modal analysis are the values of frequencies and forms of natural vibrations. Figure 11 presents three forms that determine the dynamic response of the model to kinematic excitations. The first form of natural vibrations of the presented system is displacement across the axis of the object. The second is the wave along the longitudinal axis of the object, the third figure is associated with torsion. The presented values and forms enable the proper choice of the loading direction by kinematic excitation. Having a signal with a dominant frequency around 6 Hz recorded in a given direction, it must be entered (in the example shown) in the direction transverse to the axis of the object. This will cause a coincidence of the natural frequency with the frequency forced by the seismic tremor, and, consequently, extreme dynamic response of the object. In specific conditions, the introduction of kinematic excitation of supports in a given direction does not have to correlate with the axial arrangement of the object, and can be rotated relative to it so that it causes an extreme response. Figure 12 shows the time-varying course of horizontal displacements of the pillar head relative to its base, caused by the action of horizontal components of excitation by the Z137 tremor signal. The displacement value at t = 0 s represents the response to static load (from the masses themselves). It is worth noting that the relative values of the calculated displacements in the Y direction (along the object axis) are 8 times greater than those in the X direction (across the object axis). The Z137 signal itself (Fig. 2) has twice the intensity in the Y direction than in X direction. Figure 13 shows maps of extreme normal stresses associated with the impact on the model of horizontal kinematic excitations resulting from the Z137 signal. The presentation was made in the structural elements of the platform slab and beams at the time corresponding to the largest displacement from Fig. 12.
Generally, it can be stated that the effect of the applied kinematic excitations on the object model in the form of displacements is insignificant. For example, the peak values of the displacement components along and across the axis of the viaduct crossover, above the pillar, do not exceed u y = 24 lm and u x = 35 lm. Due to the high rigidity of the supporting structure, they are accompanied by large cross-sectional sizes.

Analysis of the Characteristics and Dynamic Response of the Model-B
The static calculations gave the values of internal quantities in the critical cross-sections of the main beam from selfload and functional loads, as well as the design loadbearing capacity of these cross-sections due to bending and shearing calculated according to standards. The figures of cross-sectional values are presented in Table 3. Figure 14 presents three frequencies and forms of natural vibrations characterizing the numerical model of object B, determining the dynamic response in kinematic excitation. The first form is a bending form in the X direction, the second in the Z direction, and the third is a torsional form along the Y axis. The presented solution is consistent with the results of tests on actual objects [9,10].
Performing numerical calculations in the field of time analysis with the use of selected signals allowed for obtaining time-varying values of support reactions. Figure 15 shows the course of three translational reactions registered on the support, caused by the action of both horizontal components of excitation with the Z107 signal.
The presented results have been limited to dynamic load, and the model response is significant in the Y direction, i.e., along the axis of the object. This means that the bridge support parameters can determine the dynamic resistance of the object. THA (Time History Analysis) enables tracking changes in stress values in structural loadbearing elements. Figure 16 shows maps of extreme normal stresses associated with the impact on the model of horizontal and vertical kinematic excitations resulting from the use of the Z137 signal. Stress distribution maps refer to a specific time moment. In the case, where the vertical component is taken into account, stress extremes can be obtained at a completely different time moment than when this component is not taken into account. The example illustrated in Fig. 16 clearly shows the change in location of the greatest effort of the bridge structural elements.

Analysis of the Characteristics and Dynamic
Response of the Model-C At this point, the authors decided to present only the results of modal analysis, while all stages of model analysis were carried out in a manner analogous to previously presented cases. Figure 17 presents three forms of natural vibrations of the C model together with the corresponding frequency values. Due to the rigidity of the structure, the first form of

Analysis of the Characteristics and Dynamic Response of the Model-D
The dynamic characteristics of the steel viaduct model are depicted in Fig. 18, which shows three frequencies and forms of free vibrations. The design of the viaduct determines subsequent forms, where the first is associated with the bending form in the horizontal direction transverse to the axis of the object. The second is the bending form in the vertical direction and the third is the torsional form along the longitudinal axis of the object. It is worth noting that in the solutions of all the analysed viaduct models, regardless of the material from which they are made, the natural frequency of the bending form in the vertical direction was 5.3-6.3 Hz. In the subsequent part of the article, the method of determining the dynamic resistance of the viaduct is presented with the examples of models A and B.

Specification for Determining the Dynamic Resistance of Existing Bridge
The resistance of the supporting structure of the viaducts to the impact of mining tremors was determined by performing an exceptional form of calculation. The loads were, therefore, set according to the unique combination, in the form known as EC (Eurocode).   Taking this combination into account (5), together with the simultaneous application of cross-sectional forces and cross-sectional load capacity of the critical structural elements (due to bending and shear), the structure resistance to the impact of mining tremors can be estimated according to the following relationship (6): where • R d , R M,d , R V,d -design load capacity of the critical section, in general: due to bending moments and lateral forces, respectively, • E d,stat , E d,dyn -values of computational effects impact, • M d,stat , M d,dyn ,-design values of bending moments in critical sections, respectively: static and dynamic, • V d,stat , V d,dyn -as above but with respect to lateral forces, • a max , a 0,max -values of peak accelerations caused by mining tremors, respectively: kinematic excitations and those resulting from measurements and characterizing the most intense estimated mining tremor.
The dynamic resistance method defined above was used in the analysis of the viaduct models presented.

Estimation of Dynamic Resistance of the Viaduct-A
The assessment of the viaduct supporting structure resistance using bending moments and lateral forces was supplemented by checking the resistance of the bearing beams' support on the abutments due to slip. To this end, the total friction force transmitted by the beam support was estimated: • the reaction over the abutment of a single beam according to the design is 139.14 kN, • total friction force (10 supports) H t = 0.5Á10Á139.14 = 696.0 kN.
The results of calculations according to formula (6) are summarized in Table 4. It includes design values and those obtained from numerical calculations using three mining tremor signals.
Based on the analysis of the results of numerical calculations, it can be assumed that the resistance of the viaduct expressed by the extreme values of the resulting horizontal accelerations of the earth's surface is: • maximum acceleration value after filtration in the 0-10 Hz band about 550 mm/s 2 , • peak acceleration actual value about 700 mm/s 2 .
The presented results allow for the assumption that the dominant parameter determining the multitude of dynamic resistance is the object's ability to transfer horizontal forces through abutments. This approach was already suggested by the results presented in Fig. 12, where the displacement response dominated, recorded along the axis of the object.
The values of the estimated resistances of the A viaduct model are greater than the forecast that was 270 mm/s 2 . It should be concluded that, in this case, the anticipated dynamic impacts of mining origin will not threaten the safety the bridge structure.

Estimation of Dynamic Resistance of the Viaduct-B
The resistance of the supporting structure of the railway viaduct to the impact of mining tremors was determined by performing an exceptional form of calculation. Loads were set according to the unique combination (5) with simultaneous application of cross-sectional forces and load-bearing capacity of critical sections, due to bending and shearing. Therefore, it is possible to estimate the resistance of the structure to the impact of mining tremors according to the relationship (6). The assessment of the viaduct load-bearing structure resistance using bending moments and lateral forces was supplemented by checking the load-bearing system resistance to horizontal forces. For this purpose, the computational values have been determined: • braking forces H 1 = 208.1 kN, • forces from side impact of rolling stock H 2 = 150 kN.
The results of calculations according to formula (6) are summarized in Table 5. The calculation values of the loadbearing capacity of girders of the viaduct were taken from   Based on the analysis of the results of numerical calculations, it can be assumed that the resistance of the viaduct expressed by extreme values of horizontal and accidental accelerations of the Earth's surface is: • maximum acceleration value after band filtration 0-10 Hz ok. 650 mm/s 2 , • peak value of actual acceleration about 900 mm/s 2 .
As in the case of the viaduct A, the horizontal resistance of the railway viaduct B calculated numerically using the time analysis method (THA) is determined by the horizontal force, but induced in a direction perpendicular to the track axis. The resistances obtained with regard to bending and shearing of the main beam are not less than 2 m/s 2 of the horizontal acceleration of soil due to paraseismic tremor. Such resistance values are significantly higher than those forecasted for this area at 270 mm/s 2 .

Estimation of Dynamic Resistance of the Viaduct-C
The dynamic resistance of flyover C was carried out analogously to those mentioned above. Based on the results of the numerical analysis of the viaduct subjected to paraseismic interactions, it can be stated that: • there is no threat to the facility's safety related to exceeding the bending strength of the main beams and bearing capacity of the vertical load provided that mining tremors acting on the object will be characterized by vertical soil acceleration not greater than 2.6 m/ s 2 , • The critical limit state is the possibility of slip (travel) of the sliding (roller) bearings in the reverse direction, i.e., perpendicular to the main axis of the viaduct, resistance would be characterized by acceleration of soil in this direction equal to 1.4 m/s 2 and would signify he viaduct to be in good technical condition, • considering the poor condition of the support bearings (which results in a higher coefficient of friction), the actual viaduct resistance can be estimated at no more than 0.9 m/s 2 .
Bearing in mind the above statements, it should be assumed that the resistance of the viaduct in question expressed by extreme values of horizontal and accidental accelerations of the earth's surface is, after filtration, in the band of 0-10 Hz ok. 900 mm/s 2 .

Estimation of Dynamic Resistance of the Viaduct-D
As a result of the load capacity analysis based on the numerical model of the viaduct D, the limit state of load bearing capacity of fixed and unidirectional sliding bearings was selected due to the horizontal reaction components induced by them, caused by paraseismic loads. The supporting structure of the viaduct due to vertical variable load-movement is also subjected to horizontal forces from side impacts and associated with accelerating or braking trains. The first operate transversely to the viaduct axis, and the others along it. The side impact force is not normally located along the viaduct span. It indicates that possible impact can occur in each of the two support zones or in the span. Their full load-bearing capacity in this direction was used in the analysis. In the analysed viaduct there is a clear coupling of vertical and horizontal vibrations resulting in the appearance of horizontal components of support reactions caused by the vertical component of paraseismic kinematic excitation. This phenomenon is unusual in building constructions. In this case, it is caused by the asymmetry of the structure due to the location of the bearings on axes not perpendicular to the main axis of the span. This promotes its torsional vibrations around this axis. From the point of view of the viaduct's resistance to mining excitations, the situation is unfavorable. This is due to the fact that the vertical components and horizontal excitations jointly affect the critical elements of the supporting structure-in the case of the viaduct D-the structural bearings. Based on the results of numerical calculations, considering a slight margin of error due to the inaccuracies of the analysis, it can be assumed that the resistance of the viaduct D expressed in extreme values of horizontal and accidental accelerations of the Earth's surface, after filtration, in the 0-10 Hz band is about 850 mm/s 2 .

Summary
The above analyses of resistance of railway viaducts (B and D), as well as road viaducts (A and C), were made using reliable calculation models of load-bearing structures of both types of existing objects. The load history method with kinematic excitations was used, which were real accelerograms recorded in the region close to the location of the objects. The final assessments were made with caution and, therefore, were done conservatively. Worthy of note is the fact that the final dynamic resistance of both existing bridges was determined by the combined horizontal force, which suggests paying close attention to the methods of foundation of viaducts on abutments and columns. The general technical condition of the object as well as the condition of the support and expansion joints are not without significance. Unfortunately, these elements negatively affect the value of dynamic resistance (by lowering it) and at the same time they are often difficult to take into account.
However, it should be borne in mind and the authors clearly emphasize that the analysis carried out in this respect is unusual, and what is more, there are no results available for measuring the viaduct response exposed to mining tremors. No such measurements were carried out in any region of our country. Therefore, there is no relevant experimental base that could be used as a source of information on the behaviour of viaducts during mining tremors. Therefore, there is also no useful numerical data for comparing and possibly validating numerical models used in the analysis of various bridge structures threatened by the impact of mining tremors. So far, the authors have not received any information about possible damage to existing bridge structures resulting from paraseismic interactions. Most likely the resistance of existing bridge structures is higher than that resulting from the conducted analyses. This finding requires broader field tests in the future.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.