Научная статья на тему 'DYNAMIC INTERACTION OF HIGH-SPEED TRAINS WITH SPAN STRUCTURES AND FLEXIBLE SUPPORT'

DYNAMIC INTERACTION OF HIGH-SPEED TRAINS WITH SPAN STRUCTURES AND FLEXIBLE SUPPORT Текст научной статьи по специальности «Строительство и архитектура»

CC BY
82
16
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Magazine of Civil Engineering
Scopus
ВАК
RSCI
ESCI
Ключевые слова
RAILWAY BRIDGE / TRAIN

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Smirnov V.N., Shestakova E.B., Chizhov S.V., Antonyuk A.A., Lediaev A.P.

To ensure reliable and safe operation of the bridge structure throughout the life cycle, it is necessary to analyze and take into account many important factors, including the interaction of the main load-bearing structures. The presented work has new theoretical information, which gives the main provisions in the field of design of artificial structures for high-speed railroads. To analyze the system of a multi-span split bridge design using flexible intermediate supports of a flyover type, the Newton method (algorithm) is used. The basic data on the interaction of span structures and the design of flexible supports that are not taken into account in the design of facilities are not specified or regulated in the basic normative space, either by domestic JV standards or by foreign European EN standards, including national normative bases of the CIS countries. Harmonic analysis of the recording of the interaction of high-speed rolling stock and the joint operation of the main bridge structures of man-made structures is necessary for the design of high-speed railroads of transport infrastructure, especially in conditions of high-speed rolling stock. The article proposes a methodology for taking into account the interaction of the elements of the "bridge-train" system and determines the directions for further research to take into account the joint work and optimization of the basic designs of modern bridges and rolling stock in the region of high train speeds.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ДИНАМИЧЕСКОЕ ВЗАИМОДЕЙСТВИЕ ВЫСОКОСКОРОСТНЫХ ПОЕЗДОВ С ПРОЛЕТНЫМИ СТРОЕНИЯМИ И ГИБКИМИ ОПОРАМИ

Для обеспечения надежной и безопасной работы мостового сооружения на протяжении всего жизненного цикла необходим анализ и учет многих важных факторов, в том числе и учет взаимодействия основных несущих конструкций. Представленная работа обладает новыми сведениями теоретического характера, дающими основные положения в области проектирования искусственных сооружений для скоростных железнодорожных магистралей. Для анализа системы многопролетной разрезной конструкции моста с применением гибких промежуточных опор эстакадного типа используется метод (алгоритм) Ньютона. Представлены основные данные о взаимодействии пролетных строений и конструкции гибких опор, которые не учитываются при проектировании сооружений, не оговариваются и не регламентированы в базовом нормативном пространстве, как отечественными нормами СП, так и зарубежными европейскими нормами EN, включая национальные нормативные базы стран СНГ. Гармонический анализ учета взаимодействия высокоскоростного подвижного состава и совместной работы основных мостовых конструкций искусственных сооружений необходим для проектирования скоростных железнодорожных магистралей транспортной инфраструктуры, особенно в условиях высокоскоростного движения подвижного состава. В статье предложена методика учета взаимодействия элементов системы «мост-поезд» и определены направления дальнейших исследований для учета совместной работы и оптимизации основных конструкций современных мостов и подвижного состава в области высоких скоростей движения поездов.

Текст научной работы на тему «DYNAMIC INTERACTION OF HIGH-SPEED TRAINS WITH SPAN STRUCTURES AND FLEXIBLE SUPPORT»

doi: 10.18720/MCE.76.11

Dynamic interaction of high-speed trains with span structures and

flexible support

Динамическое взаимодействие высокоскоростных поездов с пролетными строениями и гибкими опорами

V.N. Smirnov, E.B. Shestakova, S.V. Chizhov, A.A. Antonyuk, L.A. Lediaev, I.A. Indeykin ,

Petersburg State Transport University, St. Petersburg, Russia

E.S. Evtukov,

Saint-Petersburg State University of Architecture and Civil Engineering, St. Petersburg, Russia

Д-р техн. наук, заведующий кафедры В.Н. Смирнов,

канд. техн. наук, доцент Е.Б. Шестакова, канд. техн. наук, доцент С.В. Чижов, аспирант А.А. Антонюк, д-р техн. наук, профессор А.П. Ледяев, д-р техн. наук, заведующий кафедры А.В. Индейкин,

Петербургский государственный университет путей сообщения Императора Александра I, г. Санкт-Петербург, Россия д-р техн. наук, профессор С.А. Евтюков, Санкт-Петербургский государственный архитектурно- строительный университет, г. Санкт-Петербург, Россия

Key words: railway bridge; train; supports; high- Ключевые слова: железнодорожный мост; speed railway mainline (HSRM) поезд; опоры; высокоскоростная железная

магистраль (ВСМ)

Abstract. To ensure reliable and safe operation of the bridge structure throughout the life cycle, it is necessary to analyze and take into account many important factors, including the interaction of the main load-bearing structures. The presented work has new theoretical information, which gives the main provisions in the field of design of artificial structures for high-speed railroads. To analyze the system of a multi-span split bridge design using flexible intermediate supports of a flyover type, the Newton method (algorithm) is used. The basic data on the interaction of span structures and the design of flexible supports that are not taken into account in the design of facilities are not specified or regulated in the basic normative space, either by domestic JV standards or by foreign European EN standards, including national normative bases of the CIS countries. Harmonic analysis of the recording of the interaction of high-speed rolling stock and the joint operation of the main bridge structures of man-made structures is necessary for the design of high-speed railroads of transport infrastructure, especially in conditions of high-speed rolling stock. The article proposes a methodology for taking into account the interaction of the elements of the "bridge-train" system and determines the directions for further research to take into account the joint work and optimization of the basic designs of modern bridges and rolling stock in the region of high train speeds.

Аннотация. Для обеспечения надежной и безопасной работы мостового сооружения на протяжении всего жизненного цикла необходим анализ и учет многих важных факторов, в том числе и учет взаимодействия основных несущих конструкций. Представленная работа обладает новыми сведениями теоретического характера, дающими основные положения в области проектирования искусственных сооружений для скоростных железнодорожных магистралей. Для анализа системы многопролетной разрезной конструкции моста с применением гибких промежуточных опор эстакадного типа используется метод (алгоритм) Ньютона. Представлены основные данные о взаимодействии пролетных строений и конструкции гибких опор, которые не учитываются при проектировании сооружений, не оговариваются и не регламентированы в базовом нормативном пространстве, как отечественными нормами СП, так и зарубежными европейскими нормами EN, включая национальные нормативные базы стран СНГ. Гармонический анализ учета взаимодействия высокоскоростного подвижного состава и совместной работы основных мостовых конструкций искусственных сооружений необходим для проектирования скоростных железнодорожных магистралей транспортной инфраструктуры, особенно в условиях

высокоскоростного движения подвижного состава. В статье предложена методика учета взаимодействия элементов системы «мост-поезд» и определены направления дальнейших исследований для учета совместной работы и оптимизации основных конструкций современных мостов и подвижного состава в области высоких скоростей движения поездов.

Introduction

According to preliminary forecast, the total length of high-speed mainlines (HSRM) will reach 60 thousand kilometers until 2020 in the world. Nowadays, China possesses the densest and developed network of high-speed railway mainlines in the world, holds first place in the world - 40% of total length of HSRM, will hold this position until 2020 increasing the length till 30 thousand kilometers, and 38 thousand kilometers until 2025. As a rule, more than half of the total length HSRM consists of elevates, 48-80% tunnels and bridges.

The normative documents used in the design of bridges operating in the Russian Federation [1, 2] and the countries of the near and far abroad [2, 3] formulate only general basic requirements for the corresponding projected structure that is part of the railway network. The above normative documents do not take into account and do not reflect the effect of the interaction of high-speed rolling stock and basic structures of bridges (towers and span structures), these norms do not apply at all to the design of bridges intended for high-speed traffic - HSRM. The main drawback is that the regulatory framework is not updated for the design of high-speed railways.

In compliance with the program for creating national standards, leading project Institute Gyprostroymost (The largest organization of total bridge engineering design in Russian Federation), there was prepared first redaction of Project for set of rules "Artificial constructions of high -speed railway lines. The rules of projecting and building" [5, 6].

To develop project-specific design codes (PSDC) for designing and building the infrastructure of the Moscow to Kazan high-speed long-distance railway line in order to create an up-to-date regulatory document for designing bridges on high-speed long-distance railway lines [7, 8].

In the works dedicated to analysis of bridge construction works of HSRM, as a rule, the dynamical reaction of the bridge is not evaluated on longitudinal train action that is necessary for constructions with (flexible) supports, characterized for high-speed mainlines.

Dynamic interaction of train loading on bridges usually either brings to analysis of span fluctuations or is considered as a interaction of spans and trains, as in 60-70s years [9, 10]; 80-90s years [11,12]; and 2000s [13-15].

Method

Influence of solid characteristics support of various types of railway bridges while constructing "bridge-train" models is not practically taken into account.

With the development of high-speed railways, the dynamic behaviour of trains and bridges has been studied more thoroughly. However, it is difficult to find papers in the scientific, normative and technical literatures about the lateral response of high-speed trains travel over long viaducts with pier.

The dynamic interaction between high-speed train and bridge is studied by theoretical analysis and field experiment. The main purpose of this paper is to develop a simple-model moving wheel/rail contact element, so that the sticking, sliding, and separation modes of the wheel/rail contact can be appropriately simulated. The three-dimensional (3D) contact finite element analysis for a realistic wheel and rail was used to accurately model the wheel/rail contact stiffness [16].

A computational model of train-bridge system with 24 m-span PC box girders are simulated. The dynamic responses of the bridge such as dynamic deflections, lateral amplitudes, lateral and vertical accelerations, lateral pier amplitudes, and the vehicle responses such as derail factors, offload factors, wheel/rail forces and car-body accelerations are calculated [17, 18]. Experimental and theoretical studies have been performed to determine the dynamic behavior of bridges crossed by the Korean high-speed train (KHST) [19, 20].

In the article of Russian authors deals with peculiar features of dynamic interaction of high-speed train loading and beam superstructures of bridges on the basis of numerical experiment. It also presents dependence of the quan tity of dynamic influence of trains on bridges from the speed of their movement, dynamic characteristics of the superstructures ^ore than 10 basic types and their lengths from 2.55 m to

90 m), as well as the results of pilot studies of the operation of steel superstructures at 250 km/h speed of movement of the high-speed train "Sapsan" [21].

In the article of Spain authors influence of pier height on the response of train and bridge is also studied. Continuous bridges, straight and constant section deck viaducts with variable height and tapered piers are the structures which have been considered. The height of the tested viaduct ranges between 60 and 120 m [22-25].

To solve the problem of oscillations of a railway beam bridge in its plane, taking into account the work of the bridge, the following assumptions were made.

1. A multi-span bridge is viewed as a system with a finite number of degrees of freedom.

2. The bases of the supports are perfectly elastic (the movements of the foundation of the support and the reaction of the base are connected by a linear dependence).

3. The bridge cloth is modeled by an elastic bar fixed to the elastic at the ends of the structure and lying on the base, elastic in the longitudinal direction. The last assumption is based, first, on the fact that, with terminal fasteners and elastic gaskets, it is possible for the track to work elastically at relatively large longitudinal displacements of the rails, and secondly because the values of the experimental relative displacements of the rails on the bridges at longitudinal effects of trains are very small.

4. When simulating a temporary mobile longitudinal load, the latter is assumed to be a one-dimensional system with one degree of freedom, which simplifies the solution of the problem, making it possible to obtain reliable results.

5. Dynamic impact on the bridge along the track axis is realized with a change in the time of traction forces or braking of the train f„. For example, the traction force of a locomotive can be expressed as

fn (t) = F0(1 - e ); (1.1)

and the intensity of the brake load

r(t) = z(1 - e ). (1.2)

Here Fo, т - the maximum values of the relevant factors (for example, for the freight train т = 0.1 q, where q - the uniformly distributed vertical train load), у, д - coefficients, t- time.

It is assumed that at the time of the longitudinal load on the structure the train is on the bridge, and the transmission of the horizontal longitudinal force in the case of split beam bridges is carried out through the fixed support parts of the beams and the path. Thus, between the rolling stock and the rail way of the bridge, a frictional coupling that varies in time is assumed. Between the beams of span structures and the supports in the places where the beams are supported on the movable bearing parts, friction bonds are also provided.

Under the assumed assumptions, the calculation scheme of, for example, the four-span bridge will have the form shown in Figure 1. In accordance with the design scheme, the bridge is an elastic system with elastic connections between track rails and span structures and with frictional connections in the moving support parts.

Equations of free vibrations of a bridge as a conservative system can be obtained by making Lagrange equations of the second kind, which have the form

d dz dnw

+ = 0 (13)

dt dqj dqj

where qj - displacement in the direction of the j-th generalized coordinate;

Kw - potential energy of deformation of the bridge;

tw - kinetic energy of the structure oscillations.

To determine tw it is necessary to establish the dependence of the crossings of the sections of the track u(x) from generalized coordinate's qj. Vertical displacements of the centers and longitudinal displacements of the ends of span beams, vertical and longitudinal displacements of the top of supports,

as well as displacement along the axis of the path of points located in the middle part of each intermediate support are accepted as generalized coordinates.

We will determine the desired displacements of the points of the track as a function of the displacement of the top of the beams of the span structures. To compose the differential equation of motion of the elastic rod (on the elastic in the direction along the bridge base), which simulates the rail track on the bridge, consider the element of the rod with the stiffness EAp of length dx. It can be seen from the figure that

dN =r(u - v)dx;

where

dN = r[u(x) - v(x)] (1.4)

dx

Here r[u(x) - v(x)] is the linear resistance of the track to the shear; r - coefficient of proportionality, characterizing the elastic properties of the path in the direction along the axis of the bridge (longitudinal modulus of elasticity of the under-rail base). Taking into account that, according to Hooke's law, the elementary displacement of a bar from a pair of rails with rigidity EAp is determined by expression

Ndx

du =

EAp

du N

dx EAp d 2u dN

dx2 dx EAp

(1.5)

(1.6)

Equating the expressions for dN/dx, obtained by (1.4) and (1.6), we obtain the required differential equation:

d 2u r

dx2 EAp

(u - v) (1.7)

where, f=

We represent the resulting equation in the following form

u/(x) - Y2 u(x) = - Y2 v(x); (1.8)

r

EAp

The boundary conditions for equation (1.7) are obtained by writing down the relations for the normal stresses at the ends of the rail track:

(n\ R-u(0) ,r) - R-u(l)

a(0) = ; °(L) =-T^ . (1.9)

ap aP

Representing the equation (1.9) in the form

o(0) = E-d— , o(L) = E, (1.10)

dxx=0 dxx=L

we find the following parameters

1

du ! R ■ u(о) du I _ - ^ • u(L)

u (0 ) du, _- R ■ u( dxu=0 = ; dx|x=L = EAp ■ (1/M)

The solution of equation (1.7) taking into account condition (1.11) has the form

L

ufe)=-y2 -iG(x' v(x)dx; (1.12)

0

where G(x,$ - the Green's function, taking into account the boundary conditions (1.11), has the form

I x fA{z\U(x) if 0 < x G(x,^)=j (s) () s . (1.13)

[b(^)-F(x) if x < L

Here

A(ç) = al ■ ¿Ä+ a2 ■ cA y 2, ; (1.14)

£(£) = b ■ iÄ+ b2 ■ cAy£, ; (1.14)

Where

Y- sHyL + H ■ cHyL H ■ SHYL + Y- CHYL ^ =-; =--; (1.15)

1 D D

h=R/(EAp);

bx = h/D; b 2 = y/D, D = y(y2 + H2)sHyL + 2y2H.eHyL ;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

U(x) = hshY x +y chY x;

V(x) = (h.ch y L +y shy L) shYx - (h.shY L +Y chY L)chyx. As regards the second factor of the integrand in equation (1.12), which represents the displacement of the points of the top of the i-th span structure v,(x), it is determined on the basis of the following considerations. We will assume that the span structure is bent along the half-wave of the sinusoid, that is, the vertical displacement of the point of the i-th span structure with abscissa x is expressed by the dependence

Ttx

y(x) =ymax Sin —.

It is clear that the longitudinal displacement of the lower belt of the movable end of the beam from its deflection y will amount to is 2qe. Moving the movable end of the beam along the upper fiber will be

2^ - cpHB = (— e--y = — (e )y;

n

* = ly.

where e, Hb - the eccentricity and height of the span of the span structure, respectively.

Longitudinal displacements of the points of the top of the beam of the i-th span during its deflection are given by

ViB (x) = (A + B cos n x);

li

x E { 0, 1i }. (1.16)

where:

A + B = piHBi; 0-17)

A -B = ( -piHBs. By solving the system of equations (1.17), we find

A = (1e1, B = - p,ei +vHbi • (1.18)

Substituting equation (1.18) into (1.16), we obtain the current displacement of the top of the i-th beam from its sagging along the sinusoid:

*b K K b

vi (x) = — y (ei - (ei - HBi )cos- x) = vi (x)yt • (1.19)

Accordingly, the longitudinal displacement of the beam bottom is determined by the expression

*h *b/ x j j K

V1 = vi (x) -(HBi cosyx (1.20)

li

The expression for the total longitudinal displacement of the top of the i-th beam (taking into account its deformation both under the action of longitudinal forces and the compliance of the supports in the vertical direction) will have the form

right left n H

Vi (x) = vB (x)(y - (y + ) + yH + x; (1.21)

where yright, yleft - vertical displacement of the right and left ends of the beam of the i-th span, respectively; y n,y H - longitudinal displacement, respectively, of the movable and fixed end of the beam

of the i-th span in the level of the centers of gravity of the section; y - vertical deflection of the middle of the beam of the i-th span.

Now, in the case of a multi-span bridge or viaduct (see, for example, Fig. 2), the solution of equation (1.7) can be written in the form

1

L1 L2 l3 Lk

— u(%)= J(x)dx + JG(x,%)v2(x)dx + JG(x,^)v3(x)dx +... JG(x,%)vk(x)dx (1.22)

y 0 ¿1 ¿2 ¿t-1

(0 < x < k) , k - number of spans.

Figure 1 shows

Vi (x) = vB (x)U - |] + f x (0 < x < Li);

v2(x) = vB(x^ -^J-h q4 + ^x L < x < L2); V3 (x) = vB (xqn - ^^ j + q9+ x (L2 < x < L3).

vB « = ?

el -(ei - Hm)cos j(x - L-1)

(1.23)

Considering, as an example, a four-span viaduct as a system with seventeen degrees of freedom (see Figure 1), we can, based on the above, write an expression for the displacement of the S point of the railroad track of the bridge in the form

1 L L х

--2 ) = q J g(- {x)dx + q2 J G(x, )—dx

У 0 0 l1

—1 qs J G(- Zs kB (x )dx —1 qs J G(- ZS v f (x)dx

2 о 2 A

L Г x — L ^ L2 L2 x — L -q4 JG(x,Z5) 1--1 dx + q6 JG(x,^)v(x)dx + q7 JG(x,Z5)-1 dx —

L2 V L2 ./4 4 L

1 l2 1 l3 l3 f x _

1 qs J G(x, Vf (x)dx — 1 qs J G(- Zs К (x)dx + q9 J G(- ZS 1 1 —

x—L

lз L3 ^ _ J^ l3

11J g(- zs К(x)dx + q12 J g(- zs) 7 2 dx — 1 q13 J g(- zs К(x)dx

l2 l2 l3 2 l2

^ l4 l4 ^ 4

— 1 q13 J G(- Zs Vf (x)dx + q14 J G(- Zs )

dx +

l3 l3 V 4 У

l4 l 4 X — L

+ q16 J G(- Zs К (x )dx+qn J G(- ZS ) -—lL dx ■

where Lm, Li - the abscissa of the initial and final reference sections of the i-th span structure, respectively (Figure 1).

We divide the rail whip with a length equal to the length of the viaduct into m sections and we will consider its displacements at m + 1 boundary points of these sections (Figure 1). We introduce the vector {u} of displacements of the rail track at these points

{u(Z)} ={u(Zl=0) U(Z2) ■■■ u(Zm+l)}={Ul U2 ■■■ Um+1}; (125)

then the equation (1.24) for a system with n degrees of freedom in the matrix form is written in the form

и = У Z G(x, Z) ai(x)d(x)q (1.26)

where q = {qi q2 ... qn} - column of generalized coordinates of the system; и - a column of offsets of track points within the structure; a,(x) - i-th matrix-string, given for the calculation scheme in Figure 2.

Finally, the equation 1.24 for a system with n degrees of freedom in the matrix form is written in the

form

и = Ф■ q; (1.27)

where Ф - the matrix of coefficients j determined from expression

k

ф?} = —У2 ■Z ZS )■ f (x)dx (1.28)

1=1

(S=1, 2 ..., m+1; i=1, 2 ..., k; j=1, 2 ..., n). where G(x,Zs) - function of the effect of the movements of the base on the displacements of the rails of

the bridge web, resiliently fixed at the ends of the bridge (Green's function); Zs - abscissa of the s-th point of the track;

fj- a function that determines the longitudinal displacement of the top of the i-th span when the qj of the system is moved;

m - number of sections of railroad lashing on the bridge;

k - number of spans;

n - number of degrees of freedom of the system. The matrix O can be represented in the form

O =

P( P2

(1)

Pi P2

(2)

m(m+1) m(m+1) P1 P2

.p

(1)

.P

(2)

.P

( m+1) ]

P0)

pf

■p(:+1)

(1.29)

The coefficients js ) for the example of the four-span bridge (Fig. 2) have the form

p2() = -y2

Pi(1) =-y2 J G(x,^ y, (x)dx ;

0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

L 1

J G(x, ^ )vf (x) — xdx ;

0 A

(1.30)

and so on in accordance with equation (1.24).

Returning to the Lagrange equations (1.36), we note that the potential strain energy of the railway bridge can be expressed by the formula

n = n + n + ku+K ■ ;

w r ec b pier'

(1.31)

where, nr, , nb, n ler - the potential energy of deformation, respectively, of the rails of the track, elastic connections between rails and beams of span structures, beams of span structures, supports. The kinetic energy of the oscillations of the bridge can be expressed by the dependence

=Tb +Tpter ; (1.32)

where Tb, Tpier - kinetic energy of oscillations, respectively, of span structures and supports. In the matrix form, expression (1.31) can be represented in the form

nw=(1/2)qTn q; (133)

where n - matrix of coefficients under generalized coordinates q of the system. The expression (1.33) can be written in the matrix form

Tw=(1/2)q*TT q*; (134)

where T - matrix of coefficients for the vector q*.

Substituting the expressions (1.33) and (1.34) into the Lagrange equations (1.36), we obtain a system of ordinary differential equations of the form

Tq**+nq=0. For q = z sinwt, instead of (1.35), we have

Tz - (1/w2)nz =0.

(1.35)

(1.36)

It is known that the shapes and the corresponding frequencies of free oscillations described by the system of equations (1.36) are obtained by analyzing the eigenvectors z and the corresponding frequencies of their own problem

(T - A n)z = 0;

(1.37)

where in х= \/ , w - frequency of free oscillations.

2

w

On the basis of the above methodology, it is possible to obtain a solution to the problem of free oscillations of a beam rail bridge structure in its plane, taking into account the influence of the track.

Analysis of forced oscillations of bridges with longitudinal effects of a train

In non-stationary (transient) modes of train traffic movement, characteristic for real operating conditions, there are forced oscillations of the structure in its plane.

The oscillations of the train-bridge system from the longitudinal effects of the train under the assumptions made earlier will be described by a system of nonlinear differential equations of the form

Mp** + ]g fn ¡t{ i sign [ (v+p*) - 0q]}+(M-MB) g fn sign(v+p)=O;

m

m (2.1)

T-1 T q** + T 1 B q* + T 1 nq +T 1 Q = T-1 w —Bg fn It{ I sign [(v+p*)-0 q*]}.

m

where q, q*, q** - the column of the generalized coordinates, velocities and accelerations of the bridge, respectively, as a system with n degrees of freedom (see Figure 1);

T, B, and P - the matrix of the coefficients of inertia, resistance and rigidity, respectively;

Q - column of generalized frictional forces in the movable supporting parts of span beams;

I - column of the form {0,7 1 1... 0,7} m + 1, reflecting the discreteness of the mobile load;

v - the speed of the moving load relative to the fixed path;

p - the longitudinal displacement of the moving load band caused by the bridge oscillations; g - the acceleration due to gravity;

M, Mb - the mass of the temporary mobile load is respectively common and on the bridge;

fn - function describing the impact on the bridge on the side of the train longitudinal load (for example, braking);

0 - the matrix characterizing the elasticity of the under-rail base along the axis of the path, formed from expression (1.29);

w - column of the form {1 1 1 ... 1}n.

As a result of the calculation, the movements of the points of the structure along the directions of the generalized coordinates q are determined, the values of the forces acting on the structure and applied along the directions of the generalized coordinates are determined for any instant of time by the formula

s = nq; (2.2)

and the longitudinal forces at the m + 1 point of the track on the bridge according to formula

Nm+1 = RoQq; (2.3)

where

EAP

Ro =-p- K. (2.4)

L / m

where EAp - the linear rigidity of rail track sections on a bridge of length L;

K = [-1 1

1 -2 1

........ 1 -2 1 (2.5)

1 -1]m+1

Thus, the proposed method of dynamic calculation of a beam rail bridge on the longitudinal effects of a train includes the following steps.

1. We construct the matrices T, n, 0 (according to 1.29),

2. A matrix of coefficients of resistance B is formed.

3. The columns of external loads are determined, as well as the generalized friction forces Q in the moving support parts.

4. The system of non-linear differential equations (2.1), which describes the oscillations of the interacting "train-bridge" system, is solved numerically, as a result of which the displacements along directions of the generalized coordinates qj are determined.

5. According to the expression (2.2), the forces acting on the structure along the directions of the generalized coordinates are determined for the longitudinal action of the rolling stock (braking or starting from the place).

6. According to the formula (2.3), axial dynamic forces are sought in the rails of the track on the bridge, which arise when braking or starting from the place of temporary loading.

7. The values of the dynamic coefficients to the longitudinal loads on the supports are determined

dyn

by the formula 1 + j =

q

q,

st

Where qJdyn, q}st - the dynamic and static displacement of the support, respectively, in the

st

direction of the corresponding generalized coordinate (static displacements q. are obtained, for

example, by assuming a coefficient y = 0 for the function fP determined by (2.1), which ensures a "quasistatic" effect of the train load).

According to the methodology outlined above, a dynamic calculation of a four-span railway bridge with reinforced concrete beams of span structures of 27.6 m in length, intermediate supports of 17.85 m in height and ballast rides (see Figure 1) is performed when the train brakes with a vertical load equal to 80 kN / m of the way. The module of longitudinal elasticity of the under-rail base was adopted U = 0 (the work of the rail track for longitudinal actions is not taken into account), as well as U = 4.5 MPa and U = 26 MPa (in the case of bridges with riding on wooden cross bars with crutches and riding on reinforced concrete sleepers with crushed stone ballast and separate fasteners).

The resistance of the path along the ends of the bridge was not taken into account (in Figure 1 R = 0). Integration of the system equations (2.1) was carried out numerically by the Runge-Kutta method using a special program. Some results of the calculations are given below. The role of the connections between the rail and the sub-rail base on the dynamic reaction of supports is evaluated. In Figure 2 is a graph of the dynamic movements of the top of the support along the axis of the path under the longitudinal action of the rolling stock. The vibrational nature of the motion of the support is seen, and in the absence of elastic bonds between the rails and beams of the span structures (U = 0), the support oscillations are realized at a frequency of about 3.5 Hertz (circular frequency w = 22.5 1 / s), and in the presence of such bonds (U = 26 MPa) the frequency increases to 9.7 Hertz. The values of the dynamic displacements of the support, calculated with the account of the work of the bridge cloth for longitudinal forces, are several times smaller than the corresponding displacements found by the traditional method -without taking into account the work of the path on the bridge (at U = 0).

In Fig. 2, it is shown that the largest values of the amplitude of the oscillations of the top of the support are observed in the step-like form of the longitudinal action function, when the resistivity of the train increases almost instantaneously (from fn = 0.01 to fn = 0.1), which is theoretically possible and sometimes adopted in dynamic calculations. However, the practical implementation of the fn transition from the minimum to the regulated by the norms is stretched in time, which is confirmed experimentally.

In order to reflect this circumstance in the computational model, the coefficient of friction fn in exponential calculations varied exponentially (in practice, changing discrete values of fn were determined at different instants of time).

Due to the gradual increase in the value of fn in time, the dynamic displacements of the top of the supports with longitudinal actions of the rolling stock are significantly reduced in comparison with the case of instantaneous changes in the values of fn.

As can be seen from the figure, with static application of the braking force, the amount of displacement of the top of the support (without taking into account the path) was 1.2 mm. The dynamic nature of the process of interaction between the train and the bridge caused an increase in movement by 10-12% (Figure 1). This circumstance gives grounds to assume in the future (with the accumulation of a pachet and experimental material) the possibility of introducing corresponding dynamic coefficients to the longitudinal loads on the supports, the values of which depend on the dynamic properties of the system and the nature of the interaction of the temporary load and the structure. The recommendations of the norms for assigning the values of the dynamic coefficient for the supports appear to be poorly grounded.

Result and Discussion

While passing of train by bridges sprung crew experience the dynamic interaction specified with rebound deflection of spans and other irregularities on roads. Especially, adverse conditions are created while passing of rolling-stock by multi-span bridge with equal spans. Spans can experience in this case essential vertical fluctuations, especially in case of one-type crew. On these occasions dynamic processes, leaked into system "train-bridge", in substantial measure may depend on solid, inertial and dissipative characteristic intermediate support which are projected for moderns elevates in the form of eased constructions of reduced rigidity. To define dynamic forces in supports as elements of the system is becoming more controversial task.

Calculating scheme of the system "bridge-train" is shown on the Figure 1. Elevate is presented in the view of system of discrete lumps, each of them has 2 levels of freedom. Every crew in case of this dynamic task also owns 2 levels of freedom, having an opportunity of vertical moving (bumping) and turning (galloping). Foundations are accepted absolutely hard.

Figure 1. Calculating scheme of the system "bridge-train"

By accepted assumptions differential equations of system's fluctuations "bridge-train" may be considered as:

Г

<

ъ =

M 3 q 3+A q 3

с

¿33 + ^

- M5 q5 S35 - M7 q7 ô31 - M9 q9S39 -

m MS " П MS "

- Z M_(z 1 - g¿31(г)- Z M_(z2 - g)S32(i) + SiFxsi§n Si q3+q6- qi2

i 2 i 2 v

qio =

v

mio q 10 +^10 qi0 r 10,10 m8 q8 s10,8 - m12 q12 s10 ,12- m14 q 14 14

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

j

14 14 10,14 \

m16 q 16 ¿10

v / v y

Here Zi - absolute movement of body i-th crew on point 1 (first by passing in a run), defined by rigidity spring and irregularities of passage.

qj - movement of j-th point of the bridge by directions of summing coordinations.

¿k - solitary movement of the point "k" of the bridge from force, given on point /;

Мэ - mass i-th crew;

Mi -th discrete lump of construction;

Fi- friction force of movable bearing parts, put on i-th support, from permanent and temporary loading;

s31, i = r33isin mi1 +

r

55

2

S32i = S„ sin 7TU^ +

s

55

'33

2

A1 vt

u, = — = —

(x 2 vt\ I 2 I 2 )

v - speed of train's running;

S - coefficient of passing from vertical movement in the middle of spans to longitudinal movement of construction on joint hinge's level of support part.

The solution of nonlinear differential equation system provides to get the image of dynamic interaction elements of the system "bridge-train" on conditions of high-speed running and to reveal the influence of support parameters on construction's dynamical work.

Conc/usions

1. To ensure reliable and safe operation of a bridge structure under conditions of high-speed movement within the limits of long-term operation throughout the life cycle, it is necessary to take into account the joint operation of the rolling stock and the main load-bearing bridge structures: span structures and supports.

2. Harmonic analysis of the recording of the interaction of high-speed rolling stock and the joint operation of the main bridge structures of man-made structures is necessary from the point of view of optimal design of high-speed railroads of transport infrastructure, especially in conditions of high-speed rolling stock movement.

3. As a result of using the above methodology for calculating the interaction of the elements of the "bridge-train" system, the dependence of the dynamic characteristics of bridge structures on the inertial and rigid parameters of the supports is revealed, which makes it possible to design the supports of the flyover with the provision of their minimum possible massiveness (and, consequently, material capacity). This significantly reduces the material consumption of the support and the laboriousness of its erection and, as a consequence, increases the economy of the structure, ensuring operational reliability and durability.

4. The proposed method for the bridge design makes it possible to determine the stress-strain state of the elements of the structure (including the rail unshackled path) with longitudinal influences to avoid the appearance of unacceptable efforts in the rails and at the same time to ensure the design of non-material-intensive economical supports, avoiding excessive reserves of massive railroad bridge supports, and with optimal geometric parameters and ensuring the necessary operational parameters of the structure for reliability and durability awns. Thus, by varying the geometric parameters of the supports, the designer can evaluate the dynamic response under the influence of the railway train load under different initial data and choose the optimal design solution for the bridge structure taking into account economic indicators.

11 11

5. The results of the calculation procedure presented in this article can be used in the future for practical applications and for a deeper study:

- It is recommended to update the current standards taking into account the interaction of the elements of the "bridge-train" system and the criterion of the economics of the support of bridge structures;

- Development of recommendations and an album of technical solutions for new railroad bridge supports operated at high speeds of railway transport;

- It is recommended to apply this method in software tools Revit and SOFiSTiK or others to expand the basic kit for modeling the dynamic effects on a bridge structure from high-speed trains.

References

1. Stroitelnyye pravila SP 35.13330.2011. Mosty i truby [Construction rules SP 35.13330.2011. Bridges and pipes]. Moscow, 2011. (rus)

2. Stroitelnyye normy i pravila SNiP 2.05.03-84*. Mosty i truby [Construction Norms and Rules SNiP 2.05.03-84*. Bridges and pipes]. Gosstroy Rossii. Moscow, GUP TsPP, 1998. 214 p. (rus)

3. EN 1990. Eurocode: Basis of structural design. Annex A2 application for bridges. European Committee for Standardization. 2002.

4. EN 13848-5. Railway applications - track - track geometry quality - Part 5: geometric quality assessment. European Committee for Standardization. 2009.

5. Ofitsialnyy sayt Proyektnogo Instituta Giprostroymost [Official site of the Project Institute Giprostroymost.]. [Online]. URL: http://www.giprostroymost.ru/news/events/432-vsm. (date of reference: 31.05.2017). (rus)

6. Ivanchenko I.I. Dynamic interaction of bridges and highspeed trains. Mechanics of Solids. 2011. No. 46(455).

7. Benin A.V. Dyachenko L.K., Smirnov V.N. Osobennosti proyektirovaniya i stroitelstva mostov vysokoskorostnoy zheleznodorozhnoy magistrali «Moskva - Kazan» [Features of designing and building bridges of the high-speed railway "Moscow - Kazan"]. Izvestiya Peterburgskogo universiteta putey soobshcheniya. 2015. No. 4(45). Pp. 15-20. (rus)

8. Wu Y.S., Yang Y.B., Yau J.D. Three-dimensional analysis of train-rail-bridge interaction problems. Vehicle Syst Dyn. 2001. No. 36(1). Pp. 1-35.

9. Doronin F.A., Indeykin A.V. Issledovaniye dinamiki proletnykh stroyeniy mosta s uchetom yego sovmestnoy raboty s oporami [Investigation of the dynamics of span structures of a bridge, taking into account its joint work with supports]. Byulleten rezultatov nauchnykh issledovaniy. 2012. No. 4(3). Pp. 131-136. (rus)

10. Kozmin Yu.G. Vzaimodeystviye balochnykh proletnykh stroyeniy zheleznodorozhnykh mostov i podvizhnogo sostava pri vysokikh skorostyakh dvizheniya [Interaction of beam span structures of railway bridges and rolling stock at high speeds]. Doc. Thesis: 05.23.17. LIIZhT, 1973. 356 p. (rus)

11. Zhang N., Tian Y., Xia H. A train-bridge dynamic interaction analysis method and its experimental validation. Engineering. 2016. Vol. 2(4). Pp. 528-536.

12. Bondar N.G. Kozmin Yu.G., Roytburd Z.G. et al. Vzaimodeystviye zheleznodorozhnykh mostov s podvizhnym sostavom [Interaction of railway bridges with rolling stock]. Moscow: Transport, 1984. 272 p.

13. Fridkin V.M., Arkhipenko Yu.V. Iteratsionnyy metod resheniya uravneniy dvizheniya yedinoy dinamicheskoy sistemy «mostovoye sooruzheniye podvizhnoy sostav» [Iterative method for solving the equations of motion of a single dynamic system "bridge structure rolling stock"]. Nauchnyye trudy OAO TsNIIS. Moscow, OAO TsNIIS, 2006. № 230. Pp. 32-39. (rus)

14. Chen Z.W., Zhai W.M., Cai C.B., Sun Y. Safety threshold of high-speed railway pier settlement based on train-track-

1.

3.

4.

7.

Литература

Строительные правила СП 35.13330.2011. Мосты и трубы. Актуализированная редакция Строительные нормы и правила СНиП 2.05.03-84*. (утв. Приказом Минрегиона РФ от 28.12.2010 N 822). M., 2011.

Строительные нормы и правила СНиП 2.05.03-84*. Мосты и трубы/ Госстрой России. М., ГУП ЦПП, 1998. 214 с.

EN 1990. Eurocode: Basis of structural design. Annex A2 application for bridges. European Committee for Standardization; 2002.

EN 13848-5. Railway applications - track - track geometry quality - Part 5: geometric quality assessment. European Committee for Standardization; 2009.

Официальный сайт проектного Института Гипростроймост. [Электронный ресурс]. URL: http://www.giprostroymost.ru/news/events/432-vsm. (дата обращения: 31.05.2017)

Ivanchenko I.I. Dynamic interaction of bridges and highspeed trains // Mechanics of Solids. 2011. № 46(455).

Бенин А.В. Дьяченко Л.К., Смирнов В.Н. Особенности проектирования и строительства мостов высокоскоростной железнодорожной магистрали «Москва - Казань» // Известия Петербургского университета путей сообщения. 2015. № 4(45). С. 1520.

8.

9.

Wu Y.S., Yang Y.B., Yau J.D. Three-dimensional analysis of train-rail-bridge interaction problems // Vehicle Syst Dyn. 2001. № 36(1). Pp. 1-35.

Доронин Ф.А., Индейкин А.В. Исследование динамики пролетных строений моста с учетом его совместной работы с опорами // Бюллетень результатов научных исследований. 2012. № 4(3). С. 131-136.

10. Козьмин Ю.Г. Взаимодействие балочных пролетных строений железнодорожных мостов и подвижного состава при высоких скоростях движения. Дис. докт. техн. наук: 05.23.17. ЛИИЖТ, 1973. 356 с.

11. Zhang N., Tian Y., Xia H. A train-bridge dynamic interaction analysis method and its experimental validation // Engineering. 2016. Vol. 2(4). Pp. 528-536.

12. Бондарь Н.Г., Козьмин Ю.Г., Ройтбурд З.Г. и др. Взаимодействие железнодорожных мостов с подвижным составом. М.: Транспорт, 1984. 272 с.

13. Фридкин В.М., Архипенко Ю.В. Итерационный метод решения уравнений движения единой динамической системы «мостовое сооружение подвижной состав» // Научные труды ОАО ЦНИИС. М.: ОАО ЦНИИС, 2006. № 230. С. 32-39.

14. Chen Z.W., Zhai W.M., Cai C.B., Sun Y. Safety threshold of high-speed railway pier settlement based on train-trackbridge dynamic interaction // Science China Technological Sciences. 2015. Vol. 58(2). Pp. 202-210.

15. Смирнов В.Н. Динамическая работа мостов высокоскоростных железнодорожных магистралей при продольных воздействиях поездной нагрузки: моногр. СПб. : Петербург. гос. ун-т путей сообщения, 2013. 65 с.

16. Ju S.H. A frictional contact finite element for wheel/rail

bridge dynamic interaction. Science China Technological Sciences. 2015. Vol. 58(2). Pp. 202-210.

15. Smirnov V.N. Dinamicheskaya rabota mostov vysokoskorostnykh zheleznodorozhnykh magistraley pri prodolnykh vozdeystviyakh poyezdnoy nagruzki [Dynamic work of bridges of high-speed railways with longitudinal forces of train load]. Saint-Petersburg: Peterburg. gos. un-t putey soobshcheniya, 2013. 65 p. (rus)

16. Ju S.H. A frictional contact finite element for wheel/rail dynamic simulations. Nonlinear Dynamics. 2016. No. 85(1). Pp. 365-374.

17. Xia H., Zhang N. Dynamic analysis of railway bridge under high-speed trains. Computers and Structures. 2005. No. 83(23-24). Pp. 1891-1901.

18. Claus H, Schiehlen W. Modeling and simulation of railway bogie structural vibrations. Vehicle Syst Dyn. 1997. No. 29(1). Pp. 538-552.

19. Kwark J.W., Choi E.S., Kim Y.J., Kim B.S., Kim S.I. Dynamic behavior of two-span continuous concrete bridges under moving high-speed train. Computers and Structures. 2004. No. 82(4-5). Pp. 463-474.

20. Iwnicki S.D. The Manchester benchmarks for rail vehicle simulation. Vehicle Syst Dyn Suppl. 1999. No. 31. Pp. 41 -48.

21. Kondratov V.V. Osobennosti dinamicheskogo vzaimodeystviya vysokoskorostnykh poyezdov s proletnymi stroyeniyami mostov [Features of dynamic interaction of high-speed trains with span structures of bridges]. Izvestiya Peterburgskogo universiteta putey soobshcheniya. 2013. No. 3(36). Pp. 28-37. (rus)

22. Olmos J.M., Astiz M.A. Analysis of the lateral dynamic response of high pier viaducts under high-speed train travel. Engineering Structures. 2013. No. 56. Pp. 13841401.

23. Chen G., Zhai W.M. A new wheel/rail spatially dynamic coupling model and its verification. Vehicle Syst Dyn. 2004. No. 41. Pp. 301-322.

24. Vatin N.I., Havula J., Martikainen L., Sinelnikov A.S., Orlova A.V., Salamakhin S.V. Thin-walled cross-sections and their joints: tests and FEM-modelling. Advanced Materials Research. 2014. Vol. 945. Pp. 1211-1215.

25. Vatin N.I., Rybakov V.A. Raschet metallokonstruktsiy: sedmaya stepen svobody [Calculation of metal structures: the seventh degree of freedom]. Stroyprofil. 2007. No. 2. Pp. 60-63. (rus)

26. Grebowski K., Zielinska M. Dynamic analysis of historic railway bridges in Poland in the context of adjusting them to Pendolino trains. International Journal of Applied Mechanics and Engineering. 2015. No. 20(2). Pp. 283-297.

27. Xia F., Cole C., Wolfs P. The dynamic wheel-rail contact stresses for wagon on various tracks. Proceeding of the 7th international conference on contact mechanics and wear of rail/wheel systems. Brisbane, Australia, 2006. Pp. 33-39.

28. Iwnicki S.D. The Manchester benchmarks for rail vehicle simulation. Vehicle Syst Dyn. 1998. No. 30. Pp. 295-313.

29. Xia H., De Roeck G., Zhang N., Maeck J. Experimental analysis of a high-speed railway bridge under Thalys trains. J Sound Vib. 2003. No. 268. Pp. 103-113.

30. Lin X., Zhang L., Guo Q., Zhang Y. Dynamic finite element model updating of prestressed concrete continuous box-girder bridge. Earthq Eng Eng Vib. 2009. No. 8. Pp. 399-407.

31. Xia H., Zhang N., Guo W. Analysis of resonance mechanism and conditions of train-bridge system. J Sound Vib. 2006. No. 297. Pp. 810-822.

32. Xia H., Zhang N., De Roeck G. Dynamic analysis of high speed railway bridge under articulated trains. Comput Struct. 2003. No. 81. Pp. 2467-2478.

dynamic simulations // Nonlinear Dynamics. 2016. № 85(1). Pp. 365-374.

17. Xia H., Zhang N. Dynamic analysis of railway bridge under high-speed trains // Computers and Structures. 2005. № 83(23-24). Pp. 1891-1901.

18. Claus H, Schiehlen W. Modeling and simulation of railway bogie structural vibrations // Vehicle Syst Dyn. 1997. № 29(1). Pp. 538-552.

19. Kwark J.W., Choi E.S., Kim Y.J., Kim B.S., Kim S.I. Dynamic behavior of two-span continuous concrete bridges under moving high-speed train // Computers and Structures. 2004. № 82(4-5). Pp. 463-474.

20. Iwnicki S.D. The Manchester benchmarks for rail vehicle simulation // Vehicle Syst Dyn Suppl. 1999. № 31. Pp. 41-48.

21. Кондратов В.В. Особенности динамического взаимодействия высокоскоростных поездов с пролетными строениями мостов // Известия Петербургского университета путей сообщения. 2013. № 3(36). С. 28-37.

22. Olmos J.M., Astiz M.A. Analysis of the lateral dynamic response of high pier viaducts under high-speed train travel // Engineering Structures. 2013. № 56. Pp. 1384-1401.

23. Chen G, Zhai W.M. A new wheel/rail spatially dynamic coupling model and its verification // Vehicle Syst Dyn. 2004. № 41. Pp. 301-322.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

24. Vatin N.I., Havula J., Martikainen L., Sinelnikov A.S., Orlova A.V., Salamakhin S.V. Thin-walled cross-sections and their joints: tests and FEM-modelling // Advanced Materials Research. 2014. Vol. 945. Pp. 1211-1215.

25. Ватин Н.И., Рыбаков В.А. Расчет металлоконструкций: седьмая степень свободы // Стройпрофиль. 2007. № 2. С. 60-63.

26. Grebowski K., Zielinska M. Dynamic analysis of historic railway bridges in Poland in the context of adjusting them to Pendolino trains // International Journal of Applied Mechanics and Engineering. 2015. № 20(2). Pp. 283-297.

27. Xia F., Cole C., Wolfs P. The dynamic wheel-rail contact stresses for wagon on various tracks // Proceeding of the 7th international conference on contact mechanics and wear of rail/wheel systems. Brisbane, Australia, 2006. Pp. 33-39.

28. Iwnicki S.D. The Manchester benchmarks for rail vehicle simulation // Vehicle Syst Dyn. 1998. № 30. Pp. 295-313.

29. Xia H., De Roeck G., Zhang N., Maeck J. Experimental analysis of a high-speed railway bridge under Thalys trains // J Sound Vib. 2003. № 268. Pp. 103-113.

30. Lin X., Zhang L., Guo Q., Zhang Y. Dynamic finite element model updating of prestressed concrete continuous box-girder bridge // Earthq Eng Eng Vib. 2009. № 8. Pp. 399-407.

31. Xia H., Zhang N., Guo W. Analysis of resonance mechanism and conditions of train-bridge system // J Sound Vib. 2006. № 297. Pp. 810-822.

32. Xia H., Zhang N., De Roeck G. Dynamic analysis of high speed railway bridge under articulated trains // Comput Struct. 2003. № 81. Pp. 2467-2478.

Vladimir Smirnov,

+7(812)572-61-19; svn193921@rambler.ru Ekaterina Shestakova,

+7(921)094-51-06; ekaterinamost6@gmail.com Sergei Chizhov,

+7(921)793-53-21; sergchizh@yandex.ru Anatoly Antonyuk,

+7(999)025-18-33; aaa.12.03.1992@mail.ru Alexander Lediaev,

+7(812)457-86-29; Tunnels@pgups.ru Andrey Indeykin,

+7(812)457-82-49; andrey.indeykin@mail.ru Sergey Evtukov,

+7(911)258-85-55; s.a.evt@mail.ru

Владимир Николаевич Смирнов,

+7(812)572-61-19;

эл. почта: svn193921@rambler.ru

Екатерина Борисовна Шестакова, +7(921)094-51-06;

эл. почта: ekaterinamost6@gmail.com

Сергей Владимирович Чижов,

+7(921)793-53-21;

эл. почта: sergchizh@yandex.ru

Анатолий Анатольевич Антонюк,

+7(999)025-18-33;

эл. почта: aaa.12.03.1992@mail.ru

Александр Петрович Ледяев, +7(812)457-86-29; эл. почта: Tunnels@pgups.ru

Андрей Викторович Индейкин,

+7(812)457-82-49;

эл. почта: andrey.indeykin@mail.ru

Сергей Аркадьевич Евтюков, +7(911)258-85-55; эл. почта: s.a.evt@mail.ru

© Smirnov V.N.,Shestakova E.B.,Chizhov S.V.,Antonyuk A.A.,Lediaev L.A.,Indeykin I.A.,

Evtukov E.S., 2017

i Надоели баннеры? Вы всегда можете отключить рекламу.