Determination of limit vehicle bump height in dynamic detection of simply supported bridge

Tan, Guojin; Liu, Ziyu; Jiao, Yubo; Wang, Wensheng

Vibroengineering Procedia

Browse Procedia

Published: 08 December 2016

Check for updates

Determination of limit vehicle bump height in dynamic detection of simply supported bridge

Guojin Tan¹

Ziyu Liu²

Yubo Jiao³

Wensheng Wang⁴

^{1, 2, 3, 4}College of Transportation, Jilin University, Changchun, China

Corresponding Author:

Yubo Jiao

Cite the article Download PDF

Downloads 1177

Abstract

Based on the D’Alembert principle and the hypothesis of Euler-Bernoulli beam, the vehicle-bridge coupling vibration equations of vehicle bump test were established using mode decomposition method. As bump test applies an impact pulse excitation to bridge structure, the initial conditions to calculate the system vibration equations were derived using momentum conservation law. Then the system vibration equations were solved by using Runge-Kutta method and the dynamic response of bridge was obtained through mode decomposition method. By comparing the maximum bending moment of each section at h height with the design bending moment of bridge, the limit height in bump test was determined. Finally, a simply supported hollow slab bridge was taken as an example, and the dynamic response of this bridge was calculated to determine the limit bump height in dynamic test. The correctness and effectiveness of this method was validated.

1. Introduction

The excitation methods of bridge dynamic test include self-vibration method, forced vibration method, pulsation method, etc. Vehicle bump test is the most commonly used in self-vibration method for short to medium span bridge [1-3]. The researchers have conducted lots of useful exploratory works in the field of vehicle test with obstacles. Huang et al. let the rear wheels of the truck fall down from the concrete block of 20 cm height at the deck to generate the impulsive force in dynamic test [4]. Kwasniewski et al. placed a wooden plank across the deck to simulate major deterioration of the deck surface [5]. Racanel artificially produced an obstacle at the deck in order to induce an impulse at the superstructure level in the dynamic loading of bridge [6]. Cantieni performed dynamic load test on the undisturbed pavement which placed a thick plank on the roadway [7]. The excitation imposed to bridges by vehicle test with obstacles is forced excitation, whose purpose is to obtain the impact factor of the bridge. However, the so called vehicle bump test of Chinese code suggests that the front wheels stay at bridge in static and the rear wheels fall from the obstacle of $h$ height. Then vehicle-bridge coupling system starts to free-decay vibration under impact loads which generated by the process of vehicle bump to bridge structure. Through analyzing free-decay response to obtain the free vibration characteristics (natural frequency, damping ratio, modal shape, etc.) of bridge structure. Thus, the limit bump height for the vehicle bump test is the key to obtain more obvious signal-to-noise ratio of dynamic response and not to cause damage for bridge. Literature [8] suggests that the height of vehicle bump is 5-15 cm. Literature [9] regulates that the height is 15 cm in the case of a single vehicle which is close to the standard load heavy vehicle is adopted as the bump excitation. However, the theoretical basis of the height specified is not given. Therefore, it is necessary to establish a theoretical calculation model of vehicle bump height and determine it accurately.

2. Vibration equations of vehicle-bridge coupling after the vehicle falling down

The process of vehicle falling down is shown in Fig. 1(a). Let the rear wheels of vehicle suddenly fall from the vehicle bump device without horizontal velocity and the vehicle remains stationary level after it falling down. The impulse effect is applied to bridge by utilizing the momentum generated during the falling process of vehicle’s rear wheels.

Fig. 1The vehicle-bridge coupling model of vehicle bumping process

a) The half-car planar and bridge model in the process of falling

b) The half-car planar and bridge model after the contact with bridge

In this study, a half-car planar model with four degrees of freedom is adopted.

When the rear wheels of vehicle fallen to the deck, the whole mechanical analysis of vehicle can be performed according to Fig. 1(b). The force of vehicle acts on bridge fall into two parts of static and dynamic force. Let $P_{1} (t)$ and $P_{3} (t)$ be the dynamic load and static load of front wheel respectively, $P_{2} (t)$ and $P_{4} (t)$ be those of rear wheel:

1

P_{1} (t) = \emptyset [m_{t 1} {\ddot{y}}_{t 1} + \frac{m_{b}}{2} {\ddot{y}}_{b} - \frac{I_{b}}{a} \ddot{θ}], P_{2} (t) = \emptyset [m_{t 2} {\ddot{y}}_{t 2} + \frac{m_{b}}{2} {\ddot{y}}_{b} + \frac{I_{b}}{a} \ddot{θ}],

P_{3} (t) = \emptyset [m_{t 1} g + \frac{m_{b}}{2} g], P_{4} (t) = \emptyset [m_{t 2} g + \frac{m_{b}}{2} g] .

As shown in Fig. 1(b), the equation of motion for bridge structure can be established based on D’Alembert principle and the hypothesis of Euler-Bernoulli beam when the excitation of vehicle is exerted on bridge:

2

m \ddot{y} (x, t) + c \dot{y} (x, t) + E I \frac{\partial^{4} (x, t)}{\partial x^{4}} = P_{1} (t) δ_{1} (x - (c + a)) + P_{2} (t) δ_{2} (x - c),

where $δ_{1} (x - (c + a))$ and $δ_{2} (x - c)$ are Dirac’s functions.

The dynamic deflection can be set according to mode decomposition method as follows:

3

y (x, t) \approx \sum_{i = 1}^{\infty} φ_{i} (x) q_{i} (t),

where $φ_{i} (x) = s i n (i π x / l)$ is the $i$ th order vibration mode function of simply supported bridge and $q_{i} (t)$ is the corresponding general coordinate.

Combining Eqs. (2-3) with the orthogonal property of mode shape can obtain:

4

\int_{0}^{l} m φ_{n}^{2} (x) {\ddot{q}}_{n} (t) d x + \int_{0}^{l} 2 β m φ_{n}^{2} (x) {\ddot{q}}_{n} (t) d x + ω_{n}^{2} \int_{0}^{l} m φ_{n}^{2} (x) q_{n} (t) d x

= \int_{0}^{l} φ_{n} (x) [P_{1} (t) δ_{1} (x - (c + a)) + P_{2} (t) δ_{2} (x - c)] d x,

where $β = c / 2 m$ ; $ω_{n}$ is the $n$ th natural frequency of bridge. Let $s_{n} = \int_{0}^{l} m φ_{n}^{2} (x) d x$ , then Eq. (4) can be expressed by:

5

s_{n} {\ddot{q}}_{n} (t) + 2 β {\dot{q}}_{n} (t) + ω_{n}^{2} s_{n} q_{n} (t) = φ_{n} (c + a) P_{1} (t) + φ_{n} (c) P_{2} (t) .

Equations of motion for wheel masses $m_{t 1}$ and $m_{t 2}$ :

6

m_{t 1} {\ddot{y}}_{t 1} + c_{a 1} ({\dot{y}}_{t 1} + {\dot{y}}_{1}) + k_{a 1} (y_{t 1} + y_{1}) + c_{t 1} ({\dot{y}}_{t 1} - {\dot{y}}_{b} + \frac{a}{2} \dot{θ})

+ k_{t 1} (y_{t 1} - y_{b} + \frac{a}{2} θ) = 0,

m_{t 2} {\ddot{y}}_{t 2} + c_{a 2} ({\dot{y}}_{t 2} + {\dot{y}}_{2}) + k_{a 2} (y_{t 2} + y_{2}) + c_{t 2} ({\dot{y}}_{t 2} - {\dot{y}}_{b} - \frac{a}{2} \dot{θ})

+ k_{t 2} (y_{t 2} - y_{b} - \frac{a}{2} θ) = 0 .

Equations of motion for sprung mass $m_{b}$ :

7

m_{b} {\ddot{y}}_{b} + c_{t 1} ({\dot{y}}_{b} - {\dot{y}}_{t 1} - \frac{a}{2} \dot{θ}) + k_{t 1} (y_{b} - y_{t 1} - \frac{a}{2} θ)

+ c_{t 2} ({\dot{y}}_{b} - {\dot{y}}_{t 2} + \frac{a}{2} \dot{θ}) + k_{t 2} (y_{b} - y_{t 2} + \frac{a}{2} θ) = 0,

I_{b} \ddot{θ} - \frac{a}{2} [c_{t 1} ({\dot{y}}_{b} - {\dot{y}}_{t 1} - \frac{a}{2} \dot{θ}) + k_{t 1} (y_{b} - y_{t 1} - \frac{a}{2} θ)]

+ \frac{a}{2} [c_{t 2} ({\dot{y}}_{b} - {\dot{y}}_{t 2} + \frac{a}{2} \dot{θ}) + k_{t 2} (y_{b} - y_{t 2} + \frac{a}{2} θ)] = 0,

where $y_{1}$ and $y_{2}$ are dynamic deflections at positions of front and rear wheel for bridge.

Eq. (2) and Eqs. (5-7) can be expressed as the generalized motion equation in matrix form:

8

M \ddot{u} + C \dot{u} + K u = 0,

where:

M = [\begin{matrix} m_{t 1} & 0 & 0 & 0 & 0 & \dots & 0 \\ 0 & m_{t 2} & 0 & 0 & 0 & \dots & 0 \\ 0 & 0 & m_{b} & 0 & 0 & \dots & 0 \\ 0 & 0 & 0 & I_{b} & 0 & \dots & 0 \\ A_{1} & B_{1} & C_{1} & D_{1} & s_{1} & \dots & 0 \\ ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋱ & 0 \\ A_{n} & B_{n} & C_{n} & D_{n} & 0 & 0 & s_{n} \end{matrix}],

A_{n} = - \emptyset φ_{n} (c + a) m_{t 1}, B_{n} = - \emptyset φ_{n} (c) m_{t 2},

C_{n} = - \emptyset (φ_{n} (c + a) \frac{m_{b}}{2} + φ_{n} (c) \frac{m_{b}}{2}),

D_{n} = \emptyset (φ_{n} (c + a) \frac{I_{b}}{a} - φ_{n} (c) \frac{I_{b}}{a}),

C = [\begin{matrix} {c_{a 1} + c}_{t 1} & 0 & - c_{t 1} & c_{t 1} \frac{a}{2} & φ_{1} (c + a) c_{a 1} & \dots & φ_{n} (c + a) c_{a 1} \\ 0 & c_{a 2} {+ c}_{t 2} & - c_{t 2} & - c_{t 2} \frac{a}{2} & φ_{1} (c) c_{a 2} & \dots & φ_{n} (c) c_{a 2} \\ - c_{t 1} & - c_{t 2} & c_{t 1} + c_{t 2} & - c_{t 1} \frac{a}{2} + c_{t 2} \frac{a}{2} & 0 & \dots & 0 \\ c_{t 1} \frac{a}{2} & - c_{t 2} \frac{a}{2} & - c_{t 1} \frac{a}{2} + c_{t 2} \frac{a}{2} & c_{t 1} \frac{a^{2}}{4} + c_{t 2} \frac{a^{2}}{4} & 0 & \dots & 0 \\ 0 & 0 & 0 & 0 & 2 β s_{1} & \dots & 0 \\ ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋱ & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 2 β s_{n} \end{matrix}],

K = [\begin{matrix} {k_{a 1} + k}_{t 1} & 0 & - k_{t 1} & k_{t 1} \frac{a}{2} & φ_{1} (c + a) k_{a 1} & \dots & φ_{n} (c + a) k_{a 1} \\ 0 & k_{a 2} {+ k}_{t 2} & - k_{t 2} & - k_{t 2} \frac{a}{2} & φ_{1} (c) k_{a 2} & \dots & φ_{n} (c) k_{a 2} \\ - k_{t 1} & - k_{t 2} & k_{t 1} + k_{t 2} & - k_{t 1} \frac{a}{2} + k_{t 2} \frac{a}{2} & 0 & \dots & 0 \\ c_{t 1} \frac{a}{2} & - c_{t 2} \frac{a}{2} & - c_{t 1} \frac{a}{2} + c_{t 2} \frac{a}{2} & c_{t 1} \frac{a^{2}}{4} + c_{t 2} \frac{a^{2}}{4} & 0 & \dots & 0 \\ 0 & 0 & 0 & 0 & ω_{1}^{2} s_{1} & \dots & 0 \\ ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋱ & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & ω_{n}^{2} s_{n} \end{matrix}],

u = {\{\begin{matrix} y_{t 1} & y_{t 2} & \begin{matrix} y_{b} & θ & \begin{matrix} q_{1} (t) & \dots & q_{n} (t) \end{matrix} \end{matrix} \end{matrix}\}}^{T} .

3. Solution of vibration equations and determination of bump height

According to conservation of energy, the potential energy of vehicle at the highest point of vehicle bump device is equal to kinetic energy at the moment of vehicle fallen to the deck. Let $ω$ be the rotation angular velocity of vehicle body at the moment of rear wheel falling to the deck:

9

\frac{1}{2} (I_{b} + m_{t 2} a^{2}) ω^{2} = m_{b} g \frac{h}{2} + m_{t 2} g h .

Let $γ$ be the distance between front wheel and mass center of particle system composed of $m_{b}$ and $m_{t}$ , it has:

10

m_{b} \frac{a}{2} + m_{t 2} a = (m_{b} + m_{t 2}) γ .

The total momentum $p_{0}$ of vehicle before the contact of vehicle and bridge can be written:

11

p_{0} = (m_{b} + m_{t 2}) γ ω = (m_{b} \frac{a}{2} + m_{t 2} a) ∙ \sqrt{\frac{2 (m_{b} g h / 2 + m_{t 2} g h)}{(I_{b} + m_{t 2} a^{2})}} .

Let $p_{1}$ be the momentum of bridge at the moment of contact of vehicle and bridge. Assuming that the momentum $p_{1}$ along the beam length is evenly distributed, it has:

12

p_{1} = \int_{0}^{l} \frac{p_{1}}{l} d x = \int_{0}^{l} m ∙ \dot{y} (x, t) d x .

Differentiating both sides of Eq. (12) with respect to $x$ yields:

13

\frac{p_{1}}{l} = m ∙ \dot{y} (x, t) .

Substituting Eq. (3) into Eq. (13) at the time $t = 0$ with the orthogonal property of mode shape can obtain:

14

\int_{0}^{l} \frac{p_{1}}{l} φ_{n} (x) d x = \int_{0}^{l} m φ_{n}^{2} (x) {\dot{q}}_{n} (0) d x .

Let $s_{n} = \int_{0}^{l} m φ_{n}^{2} (x) d x$ , then Eq. (14) can be expressed by:

15

{\dot{q}}_{n} (0) = \frac{\int_{0}^{l} φ_{n} (x) d x ∙ p_{1}}{s_{n} ∙ l} .

According to momentum conservation law, it has:

16

p_{1} = \emptyset [p_{0} + m_{t 1} {\dot{y}}_{t 1} + m_{t 2} {\dot{y}}_{t 2} + m_{b} {\dot{y}}_{b} + I_{b} \dot{θ}] .

Combining Eqs. (6-7) and Eqs. (11, 15, 16), the initial velocities of all degrees of freedom in $u$ can be obtained and all the initial conditions required to solve the theoretical calculation model are determined. Then the values of $u$ of each moment can be calculated by using fourth-order Runge-Kutta method. Consequently, the time-history curve of dynamic displacement of bridge $y (x, t)$ can be obtained according to mode decomposition method.

The moment value of each section of certain moment $M (t_{i}, x)$ can be expressed as:

17

M (t_{i}, x) = E I \frac{d^{2} [y (t_{i}, x) + y_{j} (x)]}{d x^{2}} .

Assuming that the design bending moment of each section of bridge is $M_{s} (x)$ , if the following equation is satisfied, $h$ can be considered as limit bump height:

18

M_{m a x} (t_{i}, x) = M_{s} (x) .

4. Numerical simulation

A simply supported hollow slab bridge with uniform section and 20 m span is used for the following numerical simulations and the bridge is composed of 5 slabs. The general layout and the lateral position of vehicle bump are shown in Fig. 2. The cross section of simply supported bridge is shown in Fig. 3. The geometrical and material parameters of vehicle and bridge are listed in Table 1.

Fig. 2The general layout and the lateral position of vehicle bump (unit: cm)

The maximum bending moment of mid-span in 1# slab at the bump height of 0, 5, 10, 15, 20 cm are calculated respectively and the calculation results are shown in Fig. 4. When the bump height of vehicle is 20 cm, the maximum bending moment of mid-span in 1# slab has reached the designed value of bending moment 1580.8 kN·m through analysis. Then it can be considered that the limit bump height of simply supported bridge is 20 cm in the case of vehicle and bump position given in this paper. Fig. 5 provides the displacement response of mid-span. Fig. 6 gives the displacement of each point in 1# slab at the moment of maximum displacement in mid-span. Fig. 7 shows that the bending moment of each point in 1# slab is at the moment of maximum bending moment in mid-span.

Fig. 3The cross section of simply supported bridge (unit: cm)

Table 1The calculation parameters of vehicle and bridge

Vehicle parameters		Bridge parameters
$m_{t 1}$ , $m_{t 2}$ / kg	4330	$l$ / m	20
$m_{b}$ / kg	3.85×10⁴	$a$ / m	4.2
$I_{b}$ / (kg∙m²)	2.446×10⁵	$c$ / m	10
$k_{a 1}$ , $k_{a 2}$ / (N∙m^-1)	4.28×10⁵	Concrete $ρ$ / (kg∙m^-3)	2500
$k_{t 1}$ , $k_{t 2}$ / (N∙m^-1)	2.535×10⁵	$E$ / Pa	3.25×10¹⁰
$c_{a 1}$ , $c_{a 2}$ / (kg∙s^-1)	0.98×10⁴	$\emptyset$	0.265
$c_{t 1}$ , $c_{t 2}$ / (kg∙s^-1)	1.96×10⁴	$β$	10 %

Fig. 4The relationship between bump height and bending moment of mid-span in 1# slab

Fig. 5The displacement response of mid-span in 1# slab

Fig. 6The displacement of each point in 1# slab at the moment of maximum displacement in mid-span

Fig. 7The bending moment of each point in 1# slab at the moment of maximum bending moment in mid-span

5. Conclusions

1) The theoretical research on the height of vehicle bump is carried out, and the feasible and effective theoretical calculation model and determination method are given, which provides a theoretical upper limit value for the selection of vehicle bump height.

2) When the theoretical limit height of vehicle bump is determined, a higher bump height can be used to obtain a higher signal-to-noise ratio of dynamic response and improve the measurement accuracy without causing damage to the bridge.

References

Farrar C. R., Duffey T. A., Cornwell P. J., et al. Excitation Methods for Bridge Structures. Office of Scientific and Technical Information Technical Reports, 1999.

Search CrossRef
Peng D. W., Hong J. X., Guo A. M., et al. Dynamic analysis and field-test of jointless bridge. Earthquake Engineering and Engineering Vibration, Vol. 25, Issue 2, 2005, p. 72-76.

Search CrossRef
Zhang J. P. Bridge Detection and Maintenance and Reinforcement. Second Edition, China Communications Press, Beijing, 2011, (in Chinese).

Search CrossRef
Huang C. S., Yang Y. B., Lu L. Y., et al. Dynamic testing and system identification of a multi-span highway bridge. Earthquake Engineering and Structural Dynamics, Vol. 28, Issue 8, 1999, p. 857-878.

Search CrossRef
Kwasniewski L., Li H., Wekezer J., et al. Finite element analysis of vehicle-bridge interaction. Finite Elements in Analysis and Design, Vol. 42, Issue 2, 2006, p. 950-959.

Search CrossRef
Racanel I. R. Static and dynamic testing of a concrete bridge in Bucharest. International Multidisciplinary Scientific Geoconference, 2011.

Search CrossRef
Cantieni R. Dynamic Load Testing of Highway Bridge. Transportation Research Record, 1984.

Search CrossRef
Test Method for Long Span Concrete Bridges (Final Recommendations). Industry Standards of Ministry of Transport, 2013, (in Chinese).

Search CrossRef
Specification for Inspection and Evaluation of Load-Bearing Capacity of Highway Bridges. Highway Science Institution of the Ministry of Transportation, 2011, (in Chinese).

Search CrossRef

About this article

Received

07 November 2016

Accepted

08 November 2016

Published

08 December 2016

SUBJECTS

Vibration in transportation engineering

Keywords

bridge engineering

dynamic detection

limit bump height

Runge-Kutta method

vehicle-bridge coupling vibration

Acknowledgements

This work was funded by the National Natural Science Foundation of China (Grant Number 51478203), Training Program for Outstanding Young Teachers of Jilin University.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.