Published: 25 April 2019

Two span bridge under moving load

Daniela Kuchárová1
Jozef Melcer2
1, 2University of Zilina, Zilina, Slovakia
Corresponding Author:
Jozef Melcer
Views 75
Reads 29
Downloads 1142

Abstract

The dynamic effect of the moving load on the bridge construction is the subject of the solution in this article. The bridge is modeled as two span continuous beam with two degrees of freedom. The assumption describing the dynamic deflection curve and the assumption describing the load distribution on individual lumped masses are adopted at the creation of bridge computational model. The plane computational model of heavy vehicle with five degrees of freedom is adopted. The problem is described by ordinary differential equations which are solved numerically by using MATLAB. The results are presented by graphical and numerical form.

1. Introduction

The problems of vehicle bridge interaction belongs to the oldest solved problems of structural dynamics. The works of the civil engineer R. Willis [1] and mathematician G. G. Stoks [2] in which they tried to clarify the breakdown of Chester Rail Bridge in England in 1847 are not only considered the first attempts to solve the problems of vehicle runway interaction but also the first works in the field of structural dynamics. Slovak and Czech Republic are world-known by the high level of bridge engineering and theoretical approach to the solution of dynamical problems of bridges. Basic knowledge concerning of the dynamic investigation of highway and railway bridges was published in monographs [3, 4]. Numerical methods offer an effective tool for the solution of this problem. Current state of computers enables to solve all the problems in real time. The results obtained from numerical analyses are used in the process of the design of optimal parameters of bridges with respect to its lifetime and reliability.

2. Computational vehicle and bridge model

This contribution works with half planar computational model of the lorry Tatra 815, Fig. 1. Equations of motion for 5 unknown functions r1t÷r5t describing the vehicle vibration are derived as ordinary differential Eq. (1), [5]:

1
r¨1t=-k1d1t+b1d˙1t+k2d2t+b2d˙2t+f2d˙2td˙cm1,
r¨2t=--ak1d1t-ab1d˙1t+bk2d2t+bb2d˙2t+f2d˙2td˙cIy1,
r¨3t=--k1d1t-b1d˙1t+k3d3t+b3d˙3tm2,
r¨4t=--k2d2t-b2d˙2t-f2d˙2td˙c+k4d4t+b4d˙4t+k5d5t+b5d˙5tm3,
r¨5t=--ck4d4t-cb4d˙4t+ck5d5t+cb5d˙5tIy3.

The contact forces Fint,i(i= 6, 7, 8) belong to individual contact points [5] are expressed as:

2
Fint,6t=-G6+k3d3t+b3d˙3t,
Fint,7t=-G7+k4d4t+b4d˙4t,
Fint,8t=-G8+k5d5t+b5d˙5t.

In the above equations d represents the deformation of joining elements, k, b, m are stiffness, damping, mass characteristics and G is gravity force. The dot above the symbol marks the derivative with respect to time t.

Fig. 1Planar vehicle computational model

Planar vehicle computational model

For the bridge the beam computational model with two degrees of freedom was adopted. The shape of deflection curve vx,t in time moment t and load distribution function ϕ(x) were adopted in the computation as sine function, Fig. 2.

Fig. 2Beam computational model of a bridge with two degrees of freedom and load distribution function

Beam computational model of a bridge with two degrees of freedom and load distribution function

Than by the [6]:

3
vx,t=ϕxvt+hx, vxt=ϕxtvt+ht,

where:

4
ϕx=sinπxl, ϕxt=sinπetl=sinωt, ω=πel,
5
Ft=Finttϕxt.

Equations of motion can be written as:

6
[m]D{v¨(t)}+2ωb[m]Dv˙t+kvt=Ft.

In the above expressions t is time coordinate, x length coordinate, dt deformation of joining elements, hx, ht road roughness, e speed of vehicle motion in [m/s], ω angular frequency, ωb damping angular frequency. Derivations of functions with respect to time are denoted by dot above the symbol.

3. Numerical solution

For the numerical simulation of vehicle motion along the bridge structure the computer program in the program language MATLAB was created. The following input data for vehicle TATRA 815 and for the bridge were used during computation.

Vehicle input data and initial conditions: k1= 287433 N/m, k2= 1522512 N/m, k3= 2550600 N/m, k4=k5= 5022720 N/m, b1= 19228 kg/s, b2= 260197 kg/s, b3= 2746 kg/s, b4=b5= 5494 kg/s, m1= 22950 kg, m2= 910 kg, m3= 2140 kg, Iy1= 62298 kg∙m2, Iy3= 932 kg∙m2, a= 3.135 m, b= 1.075 m, s= 4.210 m, c= 0.660 m, r1(0)= –0.02 m, r˙1(0)= 0.0 m/s, r2(0)= 0.00 rad, r˙2(0)= 0.0 rad/s, r3(0)= –0.002 m, r˙3(0)= 0.0 m/s, r4(0)= –0.003 m, r˙4(0)=0.0 m/s, r5(0)= 0.00 rad, r˙5(0)= 0.0 rad/s.

Bridge input data and initial conditions: μ= 19680.0 kg/m, I= 1.60622 m4, E= 3.85e10 N/m2, ωb= 0.1 rad/s, mm1= 285360.0 kg, l1= 29.0 m, mm2= 285360.0 kg, l2= 29.0 m, v1(0)= 0.0 m, v˙1(0)= 0.0 m/s, v2(0)= 0.0 m, v˙2(0)= 0.0 m/s.

The results of the numerical simulation are the time course of the vehicle and bridge vibration in graphical or numerical form. For the vehicle the time courses of the displacement components of characteristic points, corresponding the individual degrees of freedom, can be displayed. For the bridge the time courses of vertical mid span deflections of the individual bridge spans can be displayed. It is possible to change vehicle and bridge parameters, initial conditions and vehicle speed. The road profile can be smooth or a random variable.

As an example of possible outputs the time courses of bridge and vehicle vibration while the lorry Tatra 815 passes the bridge at the speed 70 km/h along smooth road surface are presented. Fig. 3 and 4 show the time courses of vertical deflections in the middle of the 1st and 2nd spans.

Fig. 3Vertical bridge deflection in the middle of the 1st span. Vehicle speed V= 70 km/h

Vertical bridge deflection in the middle of the 1st span. Vehicle speed V= 70 km/h

In Fig. 5 the time course of vertical deflection of the vehicle sprung mass gravitational center is displayed. In Fig. 6 the time course of rotation of the vehicle sprung mass is displayed.

Fig. 4Vertical bridge deflection in the middle of the 2nd span. Vehicle speed V= 70 km/h

Vertical bridge deflection in the middle of the 2nd span. Vehicle speed V= 70 km/h

Fig. 5Vertical deflection of the vehicle sprung mass gravitational center. Vehicle speed V= 70 km/h

Vertical deflection of the vehicle sprung mass gravitational center. Vehicle speed V= 70 km/h

Fig. 6Rotation of the vehicle sprung mass. Vehicle speed V= 70 km/h

Rotation of the vehicle sprung mass. Vehicle speed V= 70 km/h

4. Conclusions

Numerical modeling of the problems of vehicle and bridge interaction is an effective tool for the solution of real tasks of engineering practice. It allows numerically to obtain results that were previously obtained only by experimental test on the finished structure. Quality of obtained results is dependent on the quality of input data. Current state of computers enables the numerical processing of solved problems in real time. From the practical point of view the influence of various vehicle and bridge parameters is interested. The results obtained from numerical analyses are used in the process of design of optimal bridge parameters with respect to lifetime and reliability of the bridge structure.

References

  • Willis R. Report of the Commissioners Appointed to Inquire into the Application of Iron to Railway Structures. Stationary Office, London, 1849.
  • Stokes G. G. Discussion of a differential equation relating to the breaking of railway bridges. Transactions Cambridge Philosophic Society, Vol. 8, 1849, p. 707-220.
  • Melcer J. Dynamic Computation of Highway Bridges. Edis, University of Žilina, Žilina, 1997, (in Slovak).
  • Frýba L. Dynamics of Railway Bridges. Academia, Prague, 1992, (in Czech).
  • Kuchárová D., Lajčáková G. Modelling of moving load effect on concrete pavements. MATEC Web of Conferences, 2017.
  • Melcer J. Experimental verification of a computing model. Applied Mechanics and Materials, Vol. 732, 2015, p. 345-348.

About this article

Received
26 February 2019
Accepted
20 March 2019
Published
25 April 2019
SUBJECTS
Vibration in transportation engineering
Keywords
moving load
bridges
dynamics
numerical solution
computational models