SlideShare a Scribd company logo
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 1
Vibration Suppression of Vehicle-Bridge-Interaction System using
Multiple Tuned Mass Dampers
Mehmet Akif KOÇ1
, İsmail ESEN2
, Yusuf ÇAY3
, Ömer ÇERLEK3
, Muhammed Asım
KESERCİOĞLU3
, RecepKILIÇ3
, Hüseyin DAL3
, Mustafa EROĞLU3
1
(Mechatronics Engineering, Sakarya University, Turkey)
2
(Mechanical Engineering, Karabuk University, Turkey)
3
(Mechanical Engineering, Sakarya University, Turkey)
ABSTRACT :In this study, vibration suppression of vehicle-bridge interaction system has been studied in
terms of bridge dynamics. In vibration suppression, some configurations of tuned passive mass dampers have
been considered. Motion equations of vehicle and absorbers are coupled with the motion equation of the bridge
beam using first four mode of beam and Lagrange equations. Then the coupled equation of vehicle-bridge-
absorbers (VBA) has been solved using the fourth order Raunge-Kutta algorithm. It is proved that the
configuration of the absorbers that are placed on the anti-nodes of first and second mode shapes of the bridge
beam is the most effective one of all others.
KEYWORDS -Anti-nodes,Tuned mass damper, Vehicle bridge interaction,Vibration suppression
I. INTRODUCTION
Vibrations of bridges under moving loads
are vital in engineering and have been studied over a
century. Early studies of the moving load problem
have been carried out by structural engineers who
accepted the moving load as a moving point force.
Later, with the implementation of the increased
transport speed, the subject has also been examined
by the automotive and railway engineers. Now the
subject, as vehicle-bridge-interaction (VBI), has
been studied in terms of both vehicle and bridge
dynamics with the new direction such as vibration
suppression of the both systems. In this field, some
of analytical studies[1,2] are valuable for the
dynamics of the bridge beams under different
loading cases and end support conditions. When
higher transportation speeds are implemented, the
dynamic amplification factor (DAF) is considered as
important due to the inertial effects of the moving
loads in rail and other transportation systems.
Vibration suppression and control of the bridge
beams have been studied by [3–6]. It has been
reported that using Passive Tuned Mass Dampers
(PTMD), the vibration of a high-speed train bridge
at 300 km/h had been reduced by 21 percent.
Another application of PTMD has been reported by
[4] for a thin bridge under successive moving loads.
Under a periodic excitation [10] has been studied
vibration reduction of a thin beam by PTMD. It has
been reported in[6] that it is impossible to suppress
all the vibration due to the VBI, since the excitation
frequency of the moving load is not constant. For a
constant velocity of the load [7], for vibration
absorbers [8] , for Timoshenko beams [9], for
optimized damped-absorbers [10] are the other
valuable studies in this field. Plate structures under
moving loads are also research interests and some
FEM of such system can be found in [11]. Resonant
response of a beam plate of a high-speed a railway
bridge has been optimized by [12] using passive
viscous dampers. Vibration suppression of a plate
structure under random excitation using PTMD and
its H∞ and H2 optimizations have been given in [13]
and it is the early study of PTMD. It has reported
that the most effective use of PTMD is to place it at
a point where the mode shape is maximum at an
anti-node. One of the other implementation field of
the PTDM is defence science applications such as
projectile barrel interaction in various type of heavy
and light gun systems. For example, studies [14,15]
have reported that the muzzle vibrations of barrels
can be decreased at 50 percent using proper
optimized PTMD.
One of the biggest challenges in PTDM
use, the influence frequency of the moving vehicle is
not constant and when the vehicle velocity is altered,
it is changed. Secondly, there can be more than one
vehicle on the bridge system with different
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 2
characteristics. Thirdly, the moving mass case,
taking into account the effects of mass inertia, the
natural frequencies of the bridge beam are changed
depending on the position of the moving load.
Further, a VBI system is coupled and its
linearization is difficult.
In order to suppress vibrations of the
vehicle and bridge beam, in this study, the idea of
the usage of the anti-nodes for the location of the
PTDM in [13] has been extended, and developed
such that at least the anti-nodes of the first and
second mode shapes may help more effectively.
Since a general vibration, response of a vibrating
system can be approximated by superimposing of
the effect of each vibration mode. Using this idea, in
this study the vibration response of the beam has
been approximated including the effect of the first
four modes. Then tuned mass dampers are
considered for the anti-node of the first mode at the
mid-point of the bridge beam, and the anti-nodes of
the second mode that are the 0.25 and 0.75 of the
bridge length L. The tuning of the absorber at the
mid-point is carried out using the fundamental
frequency while the tunings of the absorbers at the
0.25 and 0.75 locations, the second mode frequency
is used. Except this, the effects of one undamped
absorber and damped absorber are also studied as a
comparison. Moreover, the effect of the other
position of the PTDM on beam has been analysed
widely. From this perspective, this study can help
both structural and automotive engineers in order to
extend the service life of the bridges and predict the
dynamic forces applied to the vehicle from VBI, and
to achieve desired ride-safety and passenger
comfort.
II. MATHEMATICAL MODELLING
In order to reduce vibrations of the vehicle
and bridge due to vehicle-bridge interaction, a
special system shown in Fig. 1 has been studied.
The system consists of a single degree of freedom
vehicle, a simply supported Euler- Bernoulli beam
on which the vehicle moves with constant velocity v,
and pasive vibration absorbers that are placed at
0.25, 0.5 and 0.75 of the length of the beam. Where
m is the mass of the vehicle, k and c spring constant
and the damping coefficient of the suspension
system, while y is the vertical displacement of the
vehicle body; and parameters ki and mi(i = 1, 2, 3),
respectively, are the stiffness and mass of the
absorbers at the given locations i=1, 2, and 3, as
shown in Fig. 1. Symbol ycis the displacement of
the wheel at the contact point of the beam and wb(x,
t) is the displacement of the beam at point x and time
t. The displacements of the absorbers are
represented by y1, y2, y3from left to rigth.
In the formulation for the vehicle-bridge-
absorbers (VBA) analysis following assumptions
will be adopted:
 The bridge is modelled as a simple
supported beam based on Euler-Bernoulli
theory.
 The vehicle is modelled as a single DOF
lumped parameter system.
 Only one vehicle is accepted moving on the
bridge with constant velocity v.
 The wheel is always in contact with the
bridge.
 The effect of road roughness upon vehicle
and bridge dynamic is not considered
during analysis in this study.
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 3
Figure 1. Model of Euler-Bernoulli Bridge beam with attached PTMDs and subjected to moving vehicle.
With all these assumptions, for the system
of vehicle-bridge- absorbers (VBA) shown in Fig.
1 the kinetic and potential energy are expressed as
follows, respectively:
2 2
0
2 2 2
1 1 2 2 3 3
( , ) ( )1
,
2
( ) ( ) ( )

 
       
 
  
  
  
L
b
k
w x t dx my t
E
m y t m y t m y t
(1a) (1a)
 
   
 
22
0
2 2
1 1 2 2
2
3 3
( , ) ( ) ( ( ), ) +
1
( ) ( / 4, ) ( ) ( / 2, ) ( , ( )),
2
( ) (3 / 4, )


 
     
 
 
    
 
  
 
 

L
b b
p b b
b
EI w x t dx k y t w t t
E k y t w L t k y t w L t H x t
k y t w L t
(1b)
whereμ and EI, refer to the mass of the
unit length and flexural rigitidy of the bridge
girder, respectively.
For the system shown in Fig.1, in order to
obtain the equations of motion one can use the
virtual work principle, Hamilton's principle and
D'Alembert's principle. In this study, motion
equation of the VBA integrated system is optained
using Langrange’s equations and mode expansion
method. For any point x on the beam at time tthe
deflection function wb (x, t) can be aproximated
using the Galerkin’s method:
1
1
1
( , ) ( ) ( ),
( , ) ( ) ( ),
( , ) ( ) ( ),
2
( ) sin , 1,2,..., .
 
 
 







 
 
  
 




n
b i bi
i
n
b i bi
i
n
b i bi
i
i
w x t x t
w x t x t
w x t x t
i x
x i n
L L
(2)
whereηbi is the i-th time dependent
generalized nodal coordinate,φi is i-th mode shape
function. The applied static axle load of the
vehicle at the contact point of the wheel can be
determined as:
( , ) ( ( ( ))),  cf x t mgH x t (3)
whereH(x-ξ(t))is the Heaviside function.
For the vehicle- bridge system the Rayleigh’s
dissipation function can be expressed as below:
 2 21
( , ) [ ( ) ( ( ), )] ( ( )) ,
2
      eq b bD c w x t c y t w t t H x t (4)
In Eq. (4) ceq is the equivalent damping
coefficient of bridge girder. In addition, for the
given system, the Lagrangian (L = Ek- Ep) is equal
to the difference of the kinetic and potential
energies. In such a case, for the five independent
coordinates, the Lagrange equations can be as
follows:
x
y
EI,c ,
vt
L
0.5L
ck
k2k1 k3
0.75L
t
0.25L
m1 m2 m3
m
v=const.
F
w(tt)b
F
y(t)
y (t)1
y (t)2
y (t)3
eq
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 4
d
0 , 1,
d ( ) ( ) ( )
   
    
    k k k
L L R
k
t p t p t p t
(5a)
d
, i 1,2,3,4,
d ( ) ( ) ( )  
   
    
     i
i i i
L L D
Q
t t t t
(5b)
With the state variables that are:
   1 2 3 4( ) ( ) , (t) (t) (t) (t) (t) ,     
T T
p t y t (6)
And the corresponding generalized forces
are:
0
( ) ( , ) , 1,2,3,4, 
L
i i cQ x f x t dx i (7)
Using the orthogonality conditions and
Galerkin approach for displacements of the beam,
(Eq. (2)) the equation of motion is obtained for
VBA interaction model given in Fig. 1. In this
case, the equation for the acceleration of the
vehicle body and absorbers are written as follows,
respectively:
 
1
- [ ( ) - ( , )]- [ ( ) - ( , )] ,  b by c y t w x t k y t w x t
m
(8a) (
 1 1
1
1
- [ ( ) - ( 4, )] , by k y t w L t
m
(8b) (
 2 2
2
1
- [ ( ) - ( 2, )] , by k y t w L t
m
(8c) (
 3 3
3
1
- [ ( ) - (3 4, )] , by k y t w L t
m
(8d)
Dynamic equations of the bridge as the n
second order ordinary differential equations can be
expressed as follows:
 
 
 
   
 
2
1
1
1 1 2 2
3 3
( ) ( ) ( ) ( )
( ), - ( )
( )
( ), - ( )
( 4, ) - ( ) ( 2, ) - ( )
(3 4, ) - ( ) 0 1,2,3,4.
   

 

  
      
  
      
 
  
 
 
i i eq i i i i
c b
i
b b
b
N t c x t t
f c w t t y t
t
k w t t y t
k w L t y t k w L t y t
k w L t y t i
(
(9)
whereΛ is:
11, 0
0 ,
for t t
elsewhere
ì £ <ïïL = í
ïïî
(10)
Where t1is the time when the vehicle
leaves the bridge. The coefficient Λis time
dependent and is used with Eigen function
φi(ξ(t))in order to determine motion equation.
In Eqs. (8a), (8b-c) and (9), shows the
total eight second order differential equations for
the vehicle, bridge and absorbers. The first four
equations are for vehicle and the absorbers and
they are transferred into first order equations by
[16]. In addition, the bridge dynamic is expressed
by using the second order differential equations in
Eq. (9). In this study, the bridge dynamic is
calculated by considering the first four-vibration
mode and it is represented with four differential
equations of second order. These four differential
equations are transferred to eight equations by [16].
In order to solve this equation system, comprising
of sixteen equations, a fourth order Runge-Kutta
method has been used[16].
III. NUMERICAL ANALYSIS
In this section, using four different VBA
models (Model A, B, C, D) the effects of the
absorbers on the vehicle and bridge dynamics have
been studied extensively. System equations of
motion Eq. (8a-d) for the VBA is solved by using
Runge-Kutta algorithm of the fourth order, with a
special m.file prepared MATLAB © environment.
The integration time step size for numerical
analysis during the study Δt = 0.01 s, the solution
final time is taken to be t = L / v. In addition, the
parameters used in this study for the systems VBA
are presented in Table 1.
When the vehicle moves over the bridge,
it forces the bridge to vibrate. At certain vehicle
speeds, the bridge enters into resonance and the
amount of oscillation increases significantly. The
velocity of the vehicle that causes resonance is
called the critical velocity vcr and it is calculated
using the circular natural frequencies of the beam.
For a simply supported beam the natural
frequencies of jth
vibration mode is expressed
as[17]:
4 4
2
4
( / )j
j EI
rad s
L
p
w
m
= (11)
When Eq. (11) and the excitation
frequency of a moving vehicle on the bridge are
rearranged, the following is obtained[16]:
1/2
2
,
2j j cr
vL v
f j EI v
w w m
a
w p p
æ ö÷ç= = = =÷ç ÷÷çè ø
(12)
Using Eq. (12), Fig. 2 shows the first three
modes of vibration of the bridge beam. Frequency
values of the first three modes of vibration of the
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 5
bridge beam are, respectively, f1=0.2108 Hz,
f2=0.8432 Hz, and f3=1.8972 Hz. Using these
values and Eq. (22), the critical speeds for vehicles
moving on the bridge are calculated, respectively,
as vcr1=42.15 m/s, vcr2=168.62 m/s, and vcr3=252.82
m/s.
Figure2.Mode shapes of the bridge.
In order to analyze the effects of the
absorbers on the dynamics of vehicles and bridges
four different models of the VBA shown in Fig. 3.
In the model A as shown in Fig. 3a, an undamped
mass spring system is suspended from the bridge at
the mid-point. The main vibration displacement of
the bridge system is affected mostly from the
fundemantal mode of the beam, for such reason the
absorber is placed at the mid point where the
maximum deflection occurs at the midpoint when
the beam vibrates with only the first mode. The
first natural frequency of the bridge is f1 = 0.2108
Hz. Therefore, the absorber is tuned to resonate at
the first natural frequency
1
1 1 1(2 ) 0.2108f k mp -
= = . Thus, it is expected
that this absorber can reduce the vibration energy
of the beam caused from the first vibration mode of
the bridge. In the second model, model B, a
viscous damper added to the spring mass system.
The damping coefficient of the absorber was
chosen as c2=10000 Nsm-1
, not on purpose but to
show the effect of any damping in the absorber.
The critical damping value
2 2 2 22 2 26476.46cr nc m k mw= = = Nsm-1
, and
the damping ratio of the damper is
( )2 2 2 2 2 2 22 (2 ) 0.377cr nc c c m c k mz w= = = =
obtained using the parameters listed in Table 1.
When this damper is added, the natural frequency
of the absorber has reduced from 0.2108 Hz to
0.195 Hz that is the damped natural frequency,
where 2
1 1.226 / 0.195 Hz.d n rad sw w z= - = =
it is expected that adding damping may not be
useful in order to reduce the vibration energy of the
first mode. However, it should be known that in
some case the damping could be considered in
terms of design criteria. These issues are examined
in detail in the following sections.
Figure 3. Model of a vehicle-bridge-absorber interaction system (a) an undamped absorber at middle
of the bridge; (b) A damped absorber at middle of the bridge; (c) three undamped absorbers at 0.25L, 0.5L, and
0.75L of bridge length; (d) without absorber.
The Models A and B can reduce oscillations originating from the first mode of the
0 20 40 60 80 100
-1
-0.5
0
0.5
1
x 10
-3
Length of the bridge (m)
Mode 1 ( 0.2108 Hz. )
Mode 2 ( 0.8432 Hz. )
Mode 3 (1.8972 Hz. )
x
y
x=vt
L
L/2
ck
k c2 2
x
y
x=vt
L
L/2
ck
x
y
x=vt
L
ck
x
y
x=vt
L
L/2
ck
k2
k2k1 k3
Model A Model B
Model C Model D
3L/4
L/4
m1 m2 m3
v=const.
m2
m2
m m
m m
mode 1 mode 1
mode 1
mode 1
mode 2
v=const.
v=const.
v=const.
EI,c ,eq
EI,c ,eq
EI,c ,eq
EI,c ,eq
a-) b-)
c-) d-)
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 7
bridge during vibration. Because these absorbers
prepared at a middle point of maximum amplitude
of the first mode of the bridge and their resonance
frequency is set to the first oscillation frequency of
the bridge. However, structures are not only
affected from the first mode of vibration the
contribution of the higher modes are also involved.
Therefore, in addition to the first vibration mode
for the second vibration mode of the bridge a
second and a third absorber are placed at the
locations 0.25 L and 0.75 L of the bridge; and they
are tuned to the second natural frequency of the
bridge beam. As seen from Figure 3 the peaks of
the second mode shape of the beam are at these
locations. In order to compare with the effects of
the three previous models, the Model D has no
absorber. The parameters of the bridge and the
vehicle used in the study are presented in Table 1.
Table 1.Prpoperties of the vehicle and bridge.
Bridge Vehicle parameters
L (m)
E
(Gpa)
I (m4
)
μ
(kg/m)
ceq
(Ns/m)
100
207
0.174
20000
1750
m (kg)
k (kg)
c (kg)
k1 ( N/m
)
k2 ( N/m
)
2172
85439.6
2219.6
280401.
59
17525.0
8
c2 ( Ns/m)
m1( Ns/m)
m2( Ns/m)
m3( Ns/m)
k3 ( N/m )
10000
10000
10000
10000
280401.
59
Based on the four different VBA models
in Figure 3, Figure 4, displays the midpoint
displacements wb (x=L/2, t) of the bridge beam for
a vehicle speed of v=90 km/h = 25 m/s that this
speed is chosen so that it is applicable in
transportation.
For Model D model that there is no
absorber, the maximum midpoint displacement is
20.34 mm when the vehicle is at about 0.75 L,
while for Model A, the maximum displacement is
observed as 15.8 mm, which means a decrease of
22.32 percent when compared with no absorber
case. After a damped absorber in Model B is used,
the maximum midpoint displacement is 17.6 mm
when the vehicle has reached to 81.2 percent of the
bridge length. When a damping added to the
absorber, the occurrence time of the maximum
bridge midpoint displacement shifted to forward,
because of the added phase difference from
damping. Damping in the passive absorbers could
not reduce the displacements much when compared
to the undamped absorbers. For Model C, a special
case when three absorbers are used, the maximum
displacement of the midpoint of the bridge is 15
mm and; that means the displacement has been
reduced by 26.25%. All the results are presented in
Table 2.
Table 2.Comparation of the four different models according to
bridge midpoint maximum displacement (m).
Model Max. Mid-point
disp.(mm)
Location of
the vehicle
(%)
According to
the Model D,
the reduction
rate (%)
A 15.8 73 22.32
B 17.6 81.2 13.47
C 15 87.5 26.25
D 20.34 75 ----
Figure 4. The effect of the absorber upon bridge midpoint dynamic response for vehicle constant velocity v=25
m/s and time step size ∆t=0.01 s (a) Bridge midpoint displacement; (b) Absorber displacement.
In traditional highway or railroad bridge
0 20 40 60 80 100
-20
-15
-10
-5
0
x 10
-3
Position of the vehicle on the bridge ( m )
Midpointdisplacementofthebridge(m)
model A
model B
model C
model D
0 20 40 60 80 100
-0.03
-0.025
-0.02
-0.015
-0.01
-0.005
0
0.005
Position of the vehicle on the bridge ( m )
Displacementoftheabsorbers(m)
Model A
Model B
Model C 0.25L
Model C 0.5L
Model C 0.75L
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 7
design, the Dynamic Amplification Factor (DAF) is
considered. DAF means the ratio of the maximum
dynamic displacement to the static displacement
when the load is located at the midpoint of the
bridge; and can be expressed as given below:
2 2
   
    
   
d s
L L
DAF R R (13)
The DAF of any bridge system can be
affected by many parameters such as vehicle
velocity, mass of the vehicle, natural frequency of
the bridge, and road roughness, etc. considering the
Models A, B, C and D, the DAFs of the different
velocities of the vhecile are presented in Figure 5.
As shown in Figure 5, for the Model D, that is no
absorber on the bridge the maximum value of the
DAF is 1.65 for the vehicle speed v = 26 m/s, while
for the Model B, the maximum value of the DAF is
1.47 for the vehicle speed v = 18.5 m / s; and for
the Model C the maximum value of the DAF is
1.45 for the vehicle speed v = 13.5 m / s; and for
the model A the DAF=1.37 at v=12.35 m/s. Up to
a speed of 31 m/s, the Model A showed the best
performance, but after this speed, the Model C was
the best in reducing the vibrations. Considering
this situation it can genarally be accepted that
velocity of the vehicle determine the behaviour of
the absorbers. This is because; the forcing
frequency of the VBI system is determined by the
velocity of the vehicle.
Figure 5. DAF response of the bridge mid-point for
four different vehicle-bridge-absorber interaction
system.
IV. CONCLUSION
In this study, in order to reduce the
vibration of the vehicle that is caused by vehicle –
bridge- interaction (VBI), the effects of various
dynamic vibration absorbers on the bridge
dynamics have been investigated. Due to the
complexity of the loading in moving load problem,
that the interaction is a non-linear coupled system
of the two multiple degree of freedom subsystems.
In the vibration of any beam system one can
consider the general forced vibration response of
the beam can be calculated by super positioning the
mode shapes of consecutive modes at least
including the first-four modes. However, in this
assumption the effect of the first fundamental mode
and the second mode is compromises the 80
percent of the response even in very flexible
systems. In such a case, the anti-node points of the
two first nodes can be used as the location of the
absorbers. The anti-node is the midpoint in the
first mode, and they are at 0.25L and 0.75 L in the
second mode shape. From the results of the
analysis, the Model C in this study, represents this
situation where three absorbers have been placed
on the beam at the location of the anti-nodes of the
first and second modes, is the best design for the
reduction of the vibrations of the bridge.
REFERENCES
[1] L. Fryba, Vibration solids and structures under
moving loads, Thomas Telford House, 1999.
[2] C.I. Bajer, B. Dyniewicz, Numerical Analysis of
Vibrations of Structures under Moving Inertial
Load, Springer, New York, 2012. doi:10.1007/
978-3-642-29548-5.
[3] H.C. Kwon, M.C. Kim, I.W. Lee, Vibration
control of bridges under moving loads, Comput.
Struct. 66 (1998) 473–480. doi:10.1016/S0045-
7949(97)00087-4.
[4] J.F. Wang, C.C. Lin, B.L. Chen, Vibration
suppression for high-speed railway bridges
using tuned mass dampers, Int. J. Solids Struct.
40 (2003) 465–491. doi:10.1016/S0020-
7683(02)00589-9.
[5] R.G. Jacquot, Optimal dynamic vibration
absorbers for general beam systems, J. Sound
Vib. 60 (1978) 535–542. doi:10.1016/S0022-
460X(78)80090-X.
[6] Y. Kajikawa, M. Okino, Y. Matsuura, Control
of traffic vibration on urban viaduct with tuned
mass dampers, J. Struct. Eng. (1989).
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90
0.4
0.6
0.8
1
1.2
1.4
1.6
Vehicle velocity (m/s)
DAF
1=2:
p
K a2=M a2 = 0:2108H z:
Model A
Model B
Model C
Model D
max. DAF
value 1.65
Vehicle
velocity
18.5
m/s,
Vehicle velocity 26
m/s,
max DAF
value 1.47
Vehicle velocity 12.35
m/s, max DAF 1.37
Vehicle velocity 13.5
m/s, max DAF 1.45
Vehicle velocity
31
m/s,DAF value
1.25
International Journal of Modern Research in Engineering and Technology (IJMRET)
www.ijmret.org Volume 2 Issue 7 ǁ November 2017
w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 8
[7] J.S. Issa, Vibration absorbers for simply
supported beams subjected to constant moving
loads, 226 (2012) 398–404.
doi:10.1177/1464419312450652.
[8] L. Öztürk, Vibration Absorbers as Controllers,
Bogazici University, 1997.
[9] M. Moghaddas, E. Esmailzadeh, R. Sedaghati,
P. Khosravi, Vibration control of Timoshenko
beam traversed by moving vehicle using
optimized tuned mass damper, J. Vib. Control.
18 (2012) 757–773.
doi:10.1177/1077546311404267.
[10] M. Tursun, E. Eskinat, H 2 Optimization of
Damped-Vibration Absorbers for Suppressing
Vibrations in Beams With Constrained
Minimization, J. Vib. Acoust. 136 (2014) 1–14.
doi:10.1115/1.4026246.
[11] İ. Esen, A new FEM procedure for transverse
and longitudinal vibration analysis of thin
rectangular plates subjected to a variable
velocity moving load along an arbitrary
trajectory, Lat. Am. J. Solids Struct. 12 (2015)
808–830.
[12] M.D. Martínez-Rodrigo, P. Museros, Optimal
design of passive viscous dampers for
controlling the resonant response of orthotropic
plates under high-speed moving loads, J. Sound
Vib. 330 (2011) 1328–1351.
doi:10.1016/j.jsv.2010.10.017.
[13] Y.L. Cheung, W.O. Wong, H∞ and H2
optimizations of a dynamic vibration absorber
for suppressing vibrations in plates, J. Sound
Vib. 320 (2009) 29–42.
doi:10.1016/j.jsv.2008.07.024.
[14] İ. Esen, M.A. Koç, Optimization of a passive
vibration absorber for a barrel using the genetic
algorithm, Expert Syst. Appl. 42 (2015) 894–
905. doi:10.1016/j.eswa.2014.08.038.
[15] A. Littlefield, E. Kathe, R. Messier, K. Olsen,
Gun barrel vibration absorber to increase
accuracy, New York, 2002.
[16] M.A. Koç, İ. Esen, Modelling and analysis of
vehicle-structure-road coupled interaction
considering structural flexibility , vehicle
parameters and road roughness †, J. Mech. Sci.
Technol. 31 (2017) 1–18. doi:10.1007/s12206-
017-0913-y.
[17] İ. Esen, Hareketli Yükler Altındaki Köprülü
Kren Kirişlerinin Dinamik (Mukavemet)
Analizi, Istanbul Technical University, 2009.

More Related Content

PDF
D011132935
PDF
SLAM of Multi-Robot System Considering Its Network Topology
PDF
Design and analysis of x y- positioning stage based on redundant parallel li...
PDF
Design and analysis of x y-q positioning stage based on
PDF
ADAPTIVE TREADMILL CONTROL BY HUMAN WILL
PDF
The Research of Induced Current in Coils When the Process of Motion of Magnet...
PDF
EXPERIMENTAL STUDY ON CONCRETE BOX GIRDER BRIDGE UNDER TRAFFIC INDUCED VIBRATION
PDF
The Effect of Arm Stiffness on the Elasto-Kinematic Properties of Single-Axle...
D011132935
SLAM of Multi-Robot System Considering Its Network Topology
Design and analysis of x y- positioning stage based on redundant parallel li...
Design and analysis of x y-q positioning stage based on
ADAPTIVE TREADMILL CONTROL BY HUMAN WILL
The Research of Induced Current in Coils When the Process of Motion of Magnet...
EXPERIMENTAL STUDY ON CONCRETE BOX GIRDER BRIDGE UNDER TRAFFIC INDUCED VIBRATION
The Effect of Arm Stiffness on the Elasto-Kinematic Properties of Single-Axle...

What's hot (19)

PDF
Measuring Axle Weight of Moving Vehicle Based on Particle Swarm Optimization
PDF
E013152730
PDF
D04452233
PDF
Coletti Dual Stage Ppt
PDF
Integrated inerter design and application
PDF
Quantum Current in Graphene Nano Scrolls Based Transistor
PDF
I012326065
PDF
An Implementation Mechanisms of SVM Control Strategies Applied to Five Levels...
PDF
To compare different turbulence models for the simulation of the flow over NA...
PDF
Extended abstract isps
PDF
3d fem
PDF
jkps.66.323
PDF
Robust Control of Rotor/AMB Systems
DOC
Analysis of Automobile Suspension
PDF
2006 Green Incorporating Pulse to Pulse Motion Effects anto Side Looking Arra...
PDF
SUSPENSION
PDF
International Journal of Computational Engineering Research(IJCER)
PDF
Fatigue Life Assessment for Power Cables in Floating Offshore Wind Turbines
Measuring Axle Weight of Moving Vehicle Based on Particle Swarm Optimization
E013152730
D04452233
Coletti Dual Stage Ppt
Integrated inerter design and application
Quantum Current in Graphene Nano Scrolls Based Transistor
I012326065
An Implementation Mechanisms of SVM Control Strategies Applied to Five Levels...
To compare different turbulence models for the simulation of the flow over NA...
Extended abstract isps
3d fem
jkps.66.323
Robust Control of Rotor/AMB Systems
Analysis of Automobile Suspension
2006 Green Incorporating Pulse to Pulse Motion Effects anto Side Looking Arra...
SUSPENSION
International Journal of Computational Engineering Research(IJCER)
Fatigue Life Assessment for Power Cables in Floating Offshore Wind Turbines
Ad

Similar to Vibration Suppression of Vehicle-Bridge-Interaction System using Multiple Tuned Mass Dampers (20)

PDF
VKI_RVAD_2005_Application_of_a_Lattice_Boltzmann_Code
PDF
Design and analysis of x y-q positioning stage based on
PDF
T130403141145
PPTX
DYNAMIC RESPONSE OF SIMPLE SUPPORTED BEAM VIBRATED UNDER MOVING LOAD
PDF
The International Journal of Engineering and Science (The IJES)
PDF
COMPARISON OF RESPONSE TO UNBALANCE OF OVERHUNG ROTOR SYSTEM FOR DIFFERENT SU...
PDF
Damage detection in cfrp plates by means of numerical modeling of lamb waves ...
PDF
E1102032834
PDF
International Journal of Engineering and Science Invention (IJESI)
PDF
Free vibration analysis of laminated composite beams using fem
PDF
Iaetsd design of a robust fuzzy logic controller for a single-link flexible m...
PDF
Effectiveness Evaluation of Seismic Protection Devices for Bridges in the PB...
PDF
PERFORMANCE OF HYBRID ELECTROMAGNETIC DAMPER FOR VEHICLE SUSPENSION
PDF
1 s2.0-s0997753817308902-main
PDF
Study on Vibration Reduction Effect of the Rail Vehicle with Axle Dynamic Vib...
PDF
Design of Compensator for Roll Control of Towing Air-Crafts
PDF
Comparative study of a flexible and rigid two degree of
PDF
A Modal Pushover Analysis on Multi-Span Bridge to Estimate Inelastic Seismic ...
PDF
C012131116
PDF
Analysis of rc bridge decks for selected national a nd internationalstandard ...
VKI_RVAD_2005_Application_of_a_Lattice_Boltzmann_Code
Design and analysis of x y-q positioning stage based on
T130403141145
DYNAMIC RESPONSE OF SIMPLE SUPPORTED BEAM VIBRATED UNDER MOVING LOAD
The International Journal of Engineering and Science (The IJES)
COMPARISON OF RESPONSE TO UNBALANCE OF OVERHUNG ROTOR SYSTEM FOR DIFFERENT SU...
Damage detection in cfrp plates by means of numerical modeling of lamb waves ...
E1102032834
International Journal of Engineering and Science Invention (IJESI)
Free vibration analysis of laminated composite beams using fem
Iaetsd design of a robust fuzzy logic controller for a single-link flexible m...
Effectiveness Evaluation of Seismic Protection Devices for Bridges in the PB...
PERFORMANCE OF HYBRID ELECTROMAGNETIC DAMPER FOR VEHICLE SUSPENSION
1 s2.0-s0997753817308902-main
Study on Vibration Reduction Effect of the Rail Vehicle with Axle Dynamic Vib...
Design of Compensator for Roll Control of Towing Air-Crafts
Comparative study of a flexible and rigid two degree of
A Modal Pushover Analysis on Multi-Span Bridge to Estimate Inelastic Seismic ...
C012131116
Analysis of rc bridge decks for selected national a nd internationalstandard ...
Ad

More from International Journal of Modern Research in Engineering and Technology (20)

PDF
Numerical Simulations of the Bond Stress-Slip Effect of Reinforced Concrete o...
PDF
Building an integrated vertical chain - a factor for sustainable construction
PDF
Applicability Study on the Optical Remote Sensing Techniques in a River
PDF
There is Always A Better Way: The Argument for Industrial Engineering
PDF
Study on the LandCover Classification using UAV Imagery
PDF
Comparative Analysis between Five Level Conventional and Modified Cascaded H-...
PDF
Cytotoxicity Studies of TiO2/ZnO Nanocomposites on Cervical Cancer Cells
PDF
Investigation of Performance Properties of Graphene Coated Fabrics
PDF
Effects of bagasse ash additive on the physiochemical and biological paramete...
PDF
Production and Analysis of Bioresin From Mango (Mangifera Indica) Kernel Oil
PDF
Particle Swarm Optimization Algorithm Based Window Function Design
PDF
Computed Tomography Image Reconstruction in 3D VoxelSpace
PDF
Antimicrobial Activity of Capsicum Essential Oil of Peppers
PDF
Design of Window Function in LABVIEW Environment
PDF
A study of the temporal flow of passenger and cargo transport in a Brazilian ...
PDF
Determination of Linear Absorption Coefficient for Different Materials
PDF
Evaluation of Naturally Occurring Radionuclide in Soil Samples from Ajiwei Mi...
PDF
Kinematics Modeling and Simulation of SCARA Robot Arm
PDF
Strength and durability assessment of concrete substructure in organic and hy...
Numerical Simulations of the Bond Stress-Slip Effect of Reinforced Concrete o...
Building an integrated vertical chain - a factor for sustainable construction
Applicability Study on the Optical Remote Sensing Techniques in a River
There is Always A Better Way: The Argument for Industrial Engineering
Study on the LandCover Classification using UAV Imagery
Comparative Analysis between Five Level Conventional and Modified Cascaded H-...
Cytotoxicity Studies of TiO2/ZnO Nanocomposites on Cervical Cancer Cells
Investigation of Performance Properties of Graphene Coated Fabrics
Effects of bagasse ash additive on the physiochemical and biological paramete...
Production and Analysis of Bioresin From Mango (Mangifera Indica) Kernel Oil
Particle Swarm Optimization Algorithm Based Window Function Design
Computed Tomography Image Reconstruction in 3D VoxelSpace
Antimicrobial Activity of Capsicum Essential Oil of Peppers
Design of Window Function in LABVIEW Environment
A study of the temporal flow of passenger and cargo transport in a Brazilian ...
Determination of Linear Absorption Coefficient for Different Materials
Evaluation of Naturally Occurring Radionuclide in Soil Samples from Ajiwei Mi...
Kinematics Modeling and Simulation of SCARA Robot Arm
Strength and durability assessment of concrete substructure in organic and hy...

Recently uploaded (20)

PPTX
CYBER-CRIMES AND SECURITY A guide to understanding
PPTX
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
PPTX
Welding lecture in detail for understanding
PDF
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
PPT
Mechanical Engineering MATERIALS Selection
PDF
Enhancing Cyber Defense Against Zero-Day Attacks using Ensemble Neural Networks
PDF
Model Code of Practice - Construction Work - 21102022 .pdf
PPTX
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
PPTX
Geodesy 1.pptx...............................................
PPTX
Lecture Notes Electrical Wiring System Components
PPT
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
PDF
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
PDF
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
PPTX
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
PPTX
Sustainable Sites - Green Building Construction
PPT
Project quality management in manufacturing
PDF
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
PPTX
Foundation to blockchain - A guide to Blockchain Tech
PDF
R24 SURVEYING LAB MANUAL for civil enggi
PPTX
CH1 Production IntroductoryConcepts.pptx
CYBER-CRIMES AND SECURITY A guide to understanding
MCN 401 KTU-2019-PPE KITS-MODULE 2.pptx
Welding lecture in detail for understanding
SM_6th-Sem__Cse_Internet-of-Things.pdf IOT
Mechanical Engineering MATERIALS Selection
Enhancing Cyber Defense Against Zero-Day Attacks using Ensemble Neural Networks
Model Code of Practice - Construction Work - 21102022 .pdf
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
Geodesy 1.pptx...............................................
Lecture Notes Electrical Wiring System Components
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
Sustainable Sites - Green Building Construction
Project quality management in manufacturing
TFEC-4-2020-Design-Guide-for-Timber-Roof-Trusses.pdf
Foundation to blockchain - A guide to Blockchain Tech
R24 SURVEYING LAB MANUAL for civil enggi
CH1 Production IntroductoryConcepts.pptx

Vibration Suppression of Vehicle-Bridge-Interaction System using Multiple Tuned Mass Dampers

  • 1. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 1 Vibration Suppression of Vehicle-Bridge-Interaction System using Multiple Tuned Mass Dampers Mehmet Akif KOÇ1 , İsmail ESEN2 , Yusuf ÇAY3 , Ömer ÇERLEK3 , Muhammed Asım KESERCİOĞLU3 , RecepKILIÇ3 , Hüseyin DAL3 , Mustafa EROĞLU3 1 (Mechatronics Engineering, Sakarya University, Turkey) 2 (Mechanical Engineering, Karabuk University, Turkey) 3 (Mechanical Engineering, Sakarya University, Turkey) ABSTRACT :In this study, vibration suppression of vehicle-bridge interaction system has been studied in terms of bridge dynamics. In vibration suppression, some configurations of tuned passive mass dampers have been considered. Motion equations of vehicle and absorbers are coupled with the motion equation of the bridge beam using first four mode of beam and Lagrange equations. Then the coupled equation of vehicle-bridge- absorbers (VBA) has been solved using the fourth order Raunge-Kutta algorithm. It is proved that the configuration of the absorbers that are placed on the anti-nodes of first and second mode shapes of the bridge beam is the most effective one of all others. KEYWORDS -Anti-nodes,Tuned mass damper, Vehicle bridge interaction,Vibration suppression I. INTRODUCTION Vibrations of bridges under moving loads are vital in engineering and have been studied over a century. Early studies of the moving load problem have been carried out by structural engineers who accepted the moving load as a moving point force. Later, with the implementation of the increased transport speed, the subject has also been examined by the automotive and railway engineers. Now the subject, as vehicle-bridge-interaction (VBI), has been studied in terms of both vehicle and bridge dynamics with the new direction such as vibration suppression of the both systems. In this field, some of analytical studies[1,2] are valuable for the dynamics of the bridge beams under different loading cases and end support conditions. When higher transportation speeds are implemented, the dynamic amplification factor (DAF) is considered as important due to the inertial effects of the moving loads in rail and other transportation systems. Vibration suppression and control of the bridge beams have been studied by [3–6]. It has been reported that using Passive Tuned Mass Dampers (PTMD), the vibration of a high-speed train bridge at 300 km/h had been reduced by 21 percent. Another application of PTMD has been reported by [4] for a thin bridge under successive moving loads. Under a periodic excitation [10] has been studied vibration reduction of a thin beam by PTMD. It has been reported in[6] that it is impossible to suppress all the vibration due to the VBI, since the excitation frequency of the moving load is not constant. For a constant velocity of the load [7], for vibration absorbers [8] , for Timoshenko beams [9], for optimized damped-absorbers [10] are the other valuable studies in this field. Plate structures under moving loads are also research interests and some FEM of such system can be found in [11]. Resonant response of a beam plate of a high-speed a railway bridge has been optimized by [12] using passive viscous dampers. Vibration suppression of a plate structure under random excitation using PTMD and its H∞ and H2 optimizations have been given in [13] and it is the early study of PTMD. It has reported that the most effective use of PTMD is to place it at a point where the mode shape is maximum at an anti-node. One of the other implementation field of the PTDM is defence science applications such as projectile barrel interaction in various type of heavy and light gun systems. For example, studies [14,15] have reported that the muzzle vibrations of barrels can be decreased at 50 percent using proper optimized PTMD. One of the biggest challenges in PTDM use, the influence frequency of the moving vehicle is not constant and when the vehicle velocity is altered, it is changed. Secondly, there can be more than one vehicle on the bridge system with different
  • 2. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 2 characteristics. Thirdly, the moving mass case, taking into account the effects of mass inertia, the natural frequencies of the bridge beam are changed depending on the position of the moving load. Further, a VBI system is coupled and its linearization is difficult. In order to suppress vibrations of the vehicle and bridge beam, in this study, the idea of the usage of the anti-nodes for the location of the PTDM in [13] has been extended, and developed such that at least the anti-nodes of the first and second mode shapes may help more effectively. Since a general vibration, response of a vibrating system can be approximated by superimposing of the effect of each vibration mode. Using this idea, in this study the vibration response of the beam has been approximated including the effect of the first four modes. Then tuned mass dampers are considered for the anti-node of the first mode at the mid-point of the bridge beam, and the anti-nodes of the second mode that are the 0.25 and 0.75 of the bridge length L. The tuning of the absorber at the mid-point is carried out using the fundamental frequency while the tunings of the absorbers at the 0.25 and 0.75 locations, the second mode frequency is used. Except this, the effects of one undamped absorber and damped absorber are also studied as a comparison. Moreover, the effect of the other position of the PTDM on beam has been analysed widely. From this perspective, this study can help both structural and automotive engineers in order to extend the service life of the bridges and predict the dynamic forces applied to the vehicle from VBI, and to achieve desired ride-safety and passenger comfort. II. MATHEMATICAL MODELLING In order to reduce vibrations of the vehicle and bridge due to vehicle-bridge interaction, a special system shown in Fig. 1 has been studied. The system consists of a single degree of freedom vehicle, a simply supported Euler- Bernoulli beam on which the vehicle moves with constant velocity v, and pasive vibration absorbers that are placed at 0.25, 0.5 and 0.75 of the length of the beam. Where m is the mass of the vehicle, k and c spring constant and the damping coefficient of the suspension system, while y is the vertical displacement of the vehicle body; and parameters ki and mi(i = 1, 2, 3), respectively, are the stiffness and mass of the absorbers at the given locations i=1, 2, and 3, as shown in Fig. 1. Symbol ycis the displacement of the wheel at the contact point of the beam and wb(x, t) is the displacement of the beam at point x and time t. The displacements of the absorbers are represented by y1, y2, y3from left to rigth. In the formulation for the vehicle-bridge- absorbers (VBA) analysis following assumptions will be adopted:  The bridge is modelled as a simple supported beam based on Euler-Bernoulli theory.  The vehicle is modelled as a single DOF lumped parameter system.  Only one vehicle is accepted moving on the bridge with constant velocity v.  The wheel is always in contact with the bridge.  The effect of road roughness upon vehicle and bridge dynamic is not considered during analysis in this study.
  • 3. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 3 Figure 1. Model of Euler-Bernoulli Bridge beam with attached PTMDs and subjected to moving vehicle. With all these assumptions, for the system of vehicle-bridge- absorbers (VBA) shown in Fig. 1 the kinetic and potential energy are expressed as follows, respectively: 2 2 0 2 2 2 1 1 2 2 3 3 ( , ) ( )1 , 2 ( ) ( ) ( )                       L b k w x t dx my t E m y t m y t m y t (1a) (1a)         22 0 2 2 1 1 2 2 2 3 3 ( , ) ( ) ( ( ), ) + 1 ( ) ( / 4, ) ( ) ( / 2, ) ( , ( )), 2 ( ) (3 / 4, )                              L b b p b b b EI w x t dx k y t w t t E k y t w L t k y t w L t H x t k y t w L t (1b) whereμ and EI, refer to the mass of the unit length and flexural rigitidy of the bridge girder, respectively. For the system shown in Fig.1, in order to obtain the equations of motion one can use the virtual work principle, Hamilton's principle and D'Alembert's principle. In this study, motion equation of the VBA integrated system is optained using Langrange’s equations and mode expansion method. For any point x on the beam at time tthe deflection function wb (x, t) can be aproximated using the Galerkin’s method: 1 1 1 ( , ) ( ) ( ), ( , ) ( ) ( ), ( , ) ( ) ( ), 2 ( ) sin , 1,2,..., .                           n b i bi i n b i bi i n b i bi i i w x t x t w x t x t w x t x t i x x i n L L (2) whereηbi is the i-th time dependent generalized nodal coordinate,φi is i-th mode shape function. The applied static axle load of the vehicle at the contact point of the wheel can be determined as: ( , ) ( ( ( ))),  cf x t mgH x t (3) whereH(x-ξ(t))is the Heaviside function. For the vehicle- bridge system the Rayleigh’s dissipation function can be expressed as below:  2 21 ( , ) [ ( ) ( ( ), )] ( ( )) , 2       eq b bD c w x t c y t w t t H x t (4) In Eq. (4) ceq is the equivalent damping coefficient of bridge girder. In addition, for the given system, the Lagrangian (L = Ek- Ep) is equal to the difference of the kinetic and potential energies. In such a case, for the five independent coordinates, the Lagrange equations can be as follows: x y EI,c , vt L 0.5L ck k2k1 k3 0.75L t 0.25L m1 m2 m3 m v=const. F w(tt)b F y(t) y (t)1 y (t)2 y (t)3 eq
  • 4. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 4 d 0 , 1, d ( ) ( ) ( )              k k k L L R k t p t p t p t (5a) d , i 1,2,3,4, d ( ) ( ) ( )                 i i i i L L D Q t t t t (5b) With the state variables that are:    1 2 3 4( ) ( ) , (t) (t) (t) (t) (t) ,      T T p t y t (6) And the corresponding generalized forces are: 0 ( ) ( , ) , 1,2,3,4,  L i i cQ x f x t dx i (7) Using the orthogonality conditions and Galerkin approach for displacements of the beam, (Eq. (2)) the equation of motion is obtained for VBA interaction model given in Fig. 1. In this case, the equation for the acceleration of the vehicle body and absorbers are written as follows, respectively:   1 - [ ( ) - ( , )]- [ ( ) - ( , )] ,  b by c y t w x t k y t w x t m (8a) (  1 1 1 1 - [ ( ) - ( 4, )] , by k y t w L t m (8b) (  2 2 2 1 - [ ( ) - ( 2, )] , by k y t w L t m (8c) (  3 3 3 1 - [ ( ) - (3 4, )] , by k y t w L t m (8d) Dynamic equations of the bridge as the n second order ordinary differential equations can be expressed as follows:             2 1 1 1 1 2 2 3 3 ( ) ( ) ( ) ( ) ( ), - ( ) ( ) ( ), - ( ) ( 4, ) - ( ) ( 2, ) - ( ) (3 4, ) - ( ) 0 1,2,3,4.                                      i i eq i i i i c b i b b b N t c x t t f c w t t y t t k w t t y t k w L t y t k w L t y t k w L t y t i ( (9) whereΛ is: 11, 0 0 , for t t elsewhere ì £ <ïïL = í ïïî (10) Where t1is the time when the vehicle leaves the bridge. The coefficient Λis time dependent and is used with Eigen function φi(ξ(t))in order to determine motion equation. In Eqs. (8a), (8b-c) and (9), shows the total eight second order differential equations for the vehicle, bridge and absorbers. The first four equations are for vehicle and the absorbers and they are transferred into first order equations by [16]. In addition, the bridge dynamic is expressed by using the second order differential equations in Eq. (9). In this study, the bridge dynamic is calculated by considering the first four-vibration mode and it is represented with four differential equations of second order. These four differential equations are transferred to eight equations by [16]. In order to solve this equation system, comprising of sixteen equations, a fourth order Runge-Kutta method has been used[16]. III. NUMERICAL ANALYSIS In this section, using four different VBA models (Model A, B, C, D) the effects of the absorbers on the vehicle and bridge dynamics have been studied extensively. System equations of motion Eq. (8a-d) for the VBA is solved by using Runge-Kutta algorithm of the fourth order, with a special m.file prepared MATLAB © environment. The integration time step size for numerical analysis during the study Δt = 0.01 s, the solution final time is taken to be t = L / v. In addition, the parameters used in this study for the systems VBA are presented in Table 1. When the vehicle moves over the bridge, it forces the bridge to vibrate. At certain vehicle speeds, the bridge enters into resonance and the amount of oscillation increases significantly. The velocity of the vehicle that causes resonance is called the critical velocity vcr and it is calculated using the circular natural frequencies of the beam. For a simply supported beam the natural frequencies of jth vibration mode is expressed as[17]: 4 4 2 4 ( / )j j EI rad s L p w m = (11) When Eq. (11) and the excitation frequency of a moving vehicle on the bridge are rearranged, the following is obtained[16]: 1/2 2 , 2j j cr vL v f j EI v w w m a w p p æ ö÷ç= = = =÷ç ÷÷çè ø (12) Using Eq. (12), Fig. 2 shows the first three modes of vibration of the bridge beam. Frequency values of the first three modes of vibration of the
  • 5. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 5 bridge beam are, respectively, f1=0.2108 Hz, f2=0.8432 Hz, and f3=1.8972 Hz. Using these values and Eq. (22), the critical speeds for vehicles moving on the bridge are calculated, respectively, as vcr1=42.15 m/s, vcr2=168.62 m/s, and vcr3=252.82 m/s. Figure2.Mode shapes of the bridge. In order to analyze the effects of the absorbers on the dynamics of vehicles and bridges four different models of the VBA shown in Fig. 3. In the model A as shown in Fig. 3a, an undamped mass spring system is suspended from the bridge at the mid-point. The main vibration displacement of the bridge system is affected mostly from the fundemantal mode of the beam, for such reason the absorber is placed at the mid point where the maximum deflection occurs at the midpoint when the beam vibrates with only the first mode. The first natural frequency of the bridge is f1 = 0.2108 Hz. Therefore, the absorber is tuned to resonate at the first natural frequency 1 1 1 1(2 ) 0.2108f k mp - = = . Thus, it is expected that this absorber can reduce the vibration energy of the beam caused from the first vibration mode of the bridge. In the second model, model B, a viscous damper added to the spring mass system. The damping coefficient of the absorber was chosen as c2=10000 Nsm-1 , not on purpose but to show the effect of any damping in the absorber. The critical damping value 2 2 2 22 2 26476.46cr nc m k mw= = = Nsm-1 , and the damping ratio of the damper is ( )2 2 2 2 2 2 22 (2 ) 0.377cr nc c c m c k mz w= = = = obtained using the parameters listed in Table 1. When this damper is added, the natural frequency of the absorber has reduced from 0.2108 Hz to 0.195 Hz that is the damped natural frequency, where 2 1 1.226 / 0.195 Hz.d n rad sw w z= - = = it is expected that adding damping may not be useful in order to reduce the vibration energy of the first mode. However, it should be known that in some case the damping could be considered in terms of design criteria. These issues are examined in detail in the following sections. Figure 3. Model of a vehicle-bridge-absorber interaction system (a) an undamped absorber at middle of the bridge; (b) A damped absorber at middle of the bridge; (c) three undamped absorbers at 0.25L, 0.5L, and 0.75L of bridge length; (d) without absorber. The Models A and B can reduce oscillations originating from the first mode of the 0 20 40 60 80 100 -1 -0.5 0 0.5 1 x 10 -3 Length of the bridge (m) Mode 1 ( 0.2108 Hz. ) Mode 2 ( 0.8432 Hz. ) Mode 3 (1.8972 Hz. ) x y x=vt L L/2 ck k c2 2 x y x=vt L L/2 ck x y x=vt L ck x y x=vt L L/2 ck k2 k2k1 k3 Model A Model B Model C Model D 3L/4 L/4 m1 m2 m3 v=const. m2 m2 m m m m mode 1 mode 1 mode 1 mode 1 mode 2 v=const. v=const. v=const. EI,c ,eq EI,c ,eq EI,c ,eq EI,c ,eq a-) b-) c-) d-)
  • 6. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 7 bridge during vibration. Because these absorbers prepared at a middle point of maximum amplitude of the first mode of the bridge and their resonance frequency is set to the first oscillation frequency of the bridge. However, structures are not only affected from the first mode of vibration the contribution of the higher modes are also involved. Therefore, in addition to the first vibration mode for the second vibration mode of the bridge a second and a third absorber are placed at the locations 0.25 L and 0.75 L of the bridge; and they are tuned to the second natural frequency of the bridge beam. As seen from Figure 3 the peaks of the second mode shape of the beam are at these locations. In order to compare with the effects of the three previous models, the Model D has no absorber. The parameters of the bridge and the vehicle used in the study are presented in Table 1. Table 1.Prpoperties of the vehicle and bridge. Bridge Vehicle parameters L (m) E (Gpa) I (m4 ) μ (kg/m) ceq (Ns/m) 100 207 0.174 20000 1750 m (kg) k (kg) c (kg) k1 ( N/m ) k2 ( N/m ) 2172 85439.6 2219.6 280401. 59 17525.0 8 c2 ( Ns/m) m1( Ns/m) m2( Ns/m) m3( Ns/m) k3 ( N/m ) 10000 10000 10000 10000 280401. 59 Based on the four different VBA models in Figure 3, Figure 4, displays the midpoint displacements wb (x=L/2, t) of the bridge beam for a vehicle speed of v=90 km/h = 25 m/s that this speed is chosen so that it is applicable in transportation. For Model D model that there is no absorber, the maximum midpoint displacement is 20.34 mm when the vehicle is at about 0.75 L, while for Model A, the maximum displacement is observed as 15.8 mm, which means a decrease of 22.32 percent when compared with no absorber case. After a damped absorber in Model B is used, the maximum midpoint displacement is 17.6 mm when the vehicle has reached to 81.2 percent of the bridge length. When a damping added to the absorber, the occurrence time of the maximum bridge midpoint displacement shifted to forward, because of the added phase difference from damping. Damping in the passive absorbers could not reduce the displacements much when compared to the undamped absorbers. For Model C, a special case when three absorbers are used, the maximum displacement of the midpoint of the bridge is 15 mm and; that means the displacement has been reduced by 26.25%. All the results are presented in Table 2. Table 2.Comparation of the four different models according to bridge midpoint maximum displacement (m). Model Max. Mid-point disp.(mm) Location of the vehicle (%) According to the Model D, the reduction rate (%) A 15.8 73 22.32 B 17.6 81.2 13.47 C 15 87.5 26.25 D 20.34 75 ---- Figure 4. The effect of the absorber upon bridge midpoint dynamic response for vehicle constant velocity v=25 m/s and time step size ∆t=0.01 s (a) Bridge midpoint displacement; (b) Absorber displacement. In traditional highway or railroad bridge 0 20 40 60 80 100 -20 -15 -10 -5 0 x 10 -3 Position of the vehicle on the bridge ( m ) Midpointdisplacementofthebridge(m) model A model B model C model D 0 20 40 60 80 100 -0.03 -0.025 -0.02 -0.015 -0.01 -0.005 0 0.005 Position of the vehicle on the bridge ( m ) Displacementoftheabsorbers(m) Model A Model B Model C 0.25L Model C 0.5L Model C 0.75L
  • 7. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 7 design, the Dynamic Amplification Factor (DAF) is considered. DAF means the ratio of the maximum dynamic displacement to the static displacement when the load is located at the midpoint of the bridge; and can be expressed as given below: 2 2              d s L L DAF R R (13) The DAF of any bridge system can be affected by many parameters such as vehicle velocity, mass of the vehicle, natural frequency of the bridge, and road roughness, etc. considering the Models A, B, C and D, the DAFs of the different velocities of the vhecile are presented in Figure 5. As shown in Figure 5, for the Model D, that is no absorber on the bridge the maximum value of the DAF is 1.65 for the vehicle speed v = 26 m/s, while for the Model B, the maximum value of the DAF is 1.47 for the vehicle speed v = 18.5 m / s; and for the Model C the maximum value of the DAF is 1.45 for the vehicle speed v = 13.5 m / s; and for the model A the DAF=1.37 at v=12.35 m/s. Up to a speed of 31 m/s, the Model A showed the best performance, but after this speed, the Model C was the best in reducing the vibrations. Considering this situation it can genarally be accepted that velocity of the vehicle determine the behaviour of the absorbers. This is because; the forcing frequency of the VBI system is determined by the velocity of the vehicle. Figure 5. DAF response of the bridge mid-point for four different vehicle-bridge-absorber interaction system. IV. CONCLUSION In this study, in order to reduce the vibration of the vehicle that is caused by vehicle – bridge- interaction (VBI), the effects of various dynamic vibration absorbers on the bridge dynamics have been investigated. Due to the complexity of the loading in moving load problem, that the interaction is a non-linear coupled system of the two multiple degree of freedom subsystems. In the vibration of any beam system one can consider the general forced vibration response of the beam can be calculated by super positioning the mode shapes of consecutive modes at least including the first-four modes. However, in this assumption the effect of the first fundamental mode and the second mode is compromises the 80 percent of the response even in very flexible systems. In such a case, the anti-node points of the two first nodes can be used as the location of the absorbers. The anti-node is the midpoint in the first mode, and they are at 0.25L and 0.75 L in the second mode shape. From the results of the analysis, the Model C in this study, represents this situation where three absorbers have been placed on the beam at the location of the anti-nodes of the first and second modes, is the best design for the reduction of the vibrations of the bridge. REFERENCES [1] L. Fryba, Vibration solids and structures under moving loads, Thomas Telford House, 1999. [2] C.I. Bajer, B. Dyniewicz, Numerical Analysis of Vibrations of Structures under Moving Inertial Load, Springer, New York, 2012. doi:10.1007/ 978-3-642-29548-5. [3] H.C. Kwon, M.C. Kim, I.W. Lee, Vibration control of bridges under moving loads, Comput. Struct. 66 (1998) 473–480. doi:10.1016/S0045- 7949(97)00087-4. [4] J.F. Wang, C.C. Lin, B.L. Chen, Vibration suppression for high-speed railway bridges using tuned mass dampers, Int. J. Solids Struct. 40 (2003) 465–491. doi:10.1016/S0020- 7683(02)00589-9. [5] R.G. Jacquot, Optimal dynamic vibration absorbers for general beam systems, J. Sound Vib. 60 (1978) 535–542. doi:10.1016/S0022- 460X(78)80090-X. [6] Y. Kajikawa, M. Okino, Y. Matsuura, Control of traffic vibration on urban viaduct with tuned mass dampers, J. Struct. Eng. (1989). 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 0.4 0.6 0.8 1 1.2 1.4 1.6 Vehicle velocity (m/s) DAF 1=2: p K a2=M a2 = 0:2108H z: Model A Model B Model C Model D max. DAF value 1.65 Vehicle velocity 18.5 m/s, Vehicle velocity 26 m/s, max DAF value 1.47 Vehicle velocity 12.35 m/s, max DAF 1.37 Vehicle velocity 13.5 m/s, max DAF 1.45 Vehicle velocity 31 m/s,DAF value 1.25
  • 8. International Journal of Modern Research in Engineering and Technology (IJMRET) www.ijmret.org Volume 2 Issue 7 ǁ November 2017 w w w . i j m r e t . o r g I S S N : 2 4 5 6 - 5 6 2 8 Page 8 [7] J.S. Issa, Vibration absorbers for simply supported beams subjected to constant moving loads, 226 (2012) 398–404. doi:10.1177/1464419312450652. [8] L. Öztürk, Vibration Absorbers as Controllers, Bogazici University, 1997. [9] M. Moghaddas, E. Esmailzadeh, R. Sedaghati, P. Khosravi, Vibration control of Timoshenko beam traversed by moving vehicle using optimized tuned mass damper, J. Vib. Control. 18 (2012) 757–773. doi:10.1177/1077546311404267. [10] M. Tursun, E. Eskinat, H 2 Optimization of Damped-Vibration Absorbers for Suppressing Vibrations in Beams With Constrained Minimization, J. Vib. Acoust. 136 (2014) 1–14. doi:10.1115/1.4026246. [11] İ. Esen, A new FEM procedure for transverse and longitudinal vibration analysis of thin rectangular plates subjected to a variable velocity moving load along an arbitrary trajectory, Lat. Am. J. Solids Struct. 12 (2015) 808–830. [12] M.D. Martínez-Rodrigo, P. Museros, Optimal design of passive viscous dampers for controlling the resonant response of orthotropic plates under high-speed moving loads, J. Sound Vib. 330 (2011) 1328–1351. doi:10.1016/j.jsv.2010.10.017. [13] Y.L. Cheung, W.O. Wong, H∞ and H2 optimizations of a dynamic vibration absorber for suppressing vibrations in plates, J. Sound Vib. 320 (2009) 29–42. doi:10.1016/j.jsv.2008.07.024. [14] İ. Esen, M.A. Koç, Optimization of a passive vibration absorber for a barrel using the genetic algorithm, Expert Syst. Appl. 42 (2015) 894– 905. doi:10.1016/j.eswa.2014.08.038. [15] A. Littlefield, E. Kathe, R. Messier, K. Olsen, Gun barrel vibration absorber to increase accuracy, New York, 2002. [16] M.A. Koç, İ. Esen, Modelling and analysis of vehicle-structure-road coupled interaction considering structural flexibility , vehicle parameters and road roughness †, J. Mech. Sci. Technol. 31 (2017) 1–18. doi:10.1007/s12206- 017-0913-y. [17] İ. Esen, Hareketli Yükler Altındaki Köprülü Kren Kirişlerinin Dinamik (Mukavemet) Analizi, Istanbul Technical University, 2009.