Vibroengineering PROCEDIA

Published: 29 June 2020

# Influence of bridge geometric nonlinearity on dynamic response of vehicle-bridge coupling system

Zhijun Ni1
Liangliang Zhang2
Tao Wu3
Binbin Zhang4
Teng Wang5
1, 2School of Civil Engineering, Chongqing University, Chongqing, 400045, China
3Central and Southern China Municipal Engineering Design and Research Institute, Wuhan, 430010, China
4, 1China Academy of Building Research, Beijing, 100013, China
5University Facilities Horticulture Laboratory in Shandong, Shouguang, 262700, China
2City College of Science and Technology, Chongqing University, Chongqing, 400045, China
Corresponding Author:
Liangliang Zhang
Views 76

#### Abstract

Based on the theory of vehicle-bridge coupling vibration, the differential equation of vehicle-bridge coupling system is set up according to different conditions. The differential equation of the system is converted into matrix form using mode decomposition method and is solved using MATLAB. The system equation has a non-linear matrix term when the geometric nonlinearity of the bridge is considered. The influence of wheel acceleration on the dynamic response of the bridge is analyzed without simplification under four speeds. The results show that it is acceptable to neglect the influence of wheel acceleration at low speeds, but it has a significant influence which must be considered at high speeds.

#### Highlights

• The system equation has a non-linear matrix term when the geometric nonlinearity of the bridge is considered. The non-linear matrix is related only to the nature of the bridge structures.
• The system equation of vehicle-bridge coupling system without simplifying wheel acceleration is established.
• The influence of wheel acceleration on the bridge dynamic response is analyzed without simplification under four speeds.

## 1. Introduction

With the remarkable and increasing span and vehicle load of bridges, as well as the gradually decreasing mass and stiffness of bridge structures, much attention has been focused on the complexity and diversity of vehicle-bridge coupling vibration.

There are numerous methods of analyzing vehicle-bridge coupling vibration. Xia  presented a method of simplifying the bridge into a modal model to analyze the vehicle-bridge-pier system. Based on the finite element method, Chen  analyzed the resonance of a rigid-frame bridge under vehicle-bridge coupling vibration. Li  conducted in-depth analyses of the vehicle-bridge interaction under stochastic vibration. Shen  used an ODE function based on the Runge-Kutta method to compile a secondary development function which can solve the differential equation of the vehicle-bridge coupling system.

The influence of the geometric nonlinearity of the structures on the vehicle-bridge coupling vibration is important, especially for long-span and low-mass bridges. Considering the geometric nonlinearity of the bridge, it can be more realistic to analyze the dynamic response of the structure, which is of great theoretical and practical significance in the design and construction of the bridge structure.

## 2. Establishment of system equation considering the geometric nonlinearity of the bridge

Fig. 1 shows the simply supported beam subjected to sprung mass, the vehicle system passes through the beam at a uniform speed $v$, $w\left(x,t\right)$ and $Z\left(t\right)$ denotes the dynamic deflection of the beam and the mass ${M}_{2}$, respectively. Meanwhile, the displacement of the mass ${M}_{1}$ is consistent with the deflection of the beam at its position expressed as $w\left(x,t\right)$.

The main derivation process of the differential equation is as follows. First, the dynamic equilibrium equation of a simply-supported beam can be written as Eq. (1):

1
$EI\frac{{\partial }^{4}w\left(x,t\right)}{\partial {x}^{4}}+\rho A\frac{{\partial }^{2}w\left(x,t\right)}{\partial {t}^{2}}-\frac{3EA}{2}\frac{{\partial }^{2}w\left(x,t\right)}{\partial {x}^{2}}{\left(\frac{\partial w\left(x,t\right)}{\partial x}\right)}^{2}+c\frac{\partial w\left(x,t\right)}{\partial t}$

Fig. 1Simply supported beam subjected to sprung mass model The dynamic equilibrium equation of ${M}_{2}$ can be written as:

2
${M}_{2}Z\left(t\right)+{k}_{1}\left[\stackrel{¨}{Z}\left(t\right)-{w\left(x,t\right)|}_{x=vt}\right]+{c}_{1}\left[Z\left(t\right)-{\frac{\partial \stackrel{˙}{w}\left(x,t\right)}{\partial t}|}_{x=vt}\right]=0.$

The external load acting on the beam can be expressed as Eq. (3):

3
$P\left(x,t\right)=\delta \left(x-vt\right)\left[\begin{array}{c}\left({M}_{1}+{M}_{2}\right)g-{M}_{1}\frac{{\partial }^{2}w\left(x,t\right)}{\partial {t}^{2}}\\ +{k}_{1}\left[Z\left(t\right)-w\left(x,t\right)\right]+{c}_{1}\left[\stackrel{˙}{Z}\left(t\right)-\frac{\partial w\left(x,t\right)}{\partial t}\right]\end{array}\right].$

According to the mode decomposition method, the right end of Eq. (3) can be transformed into:

4
${P}_{n}\left(t\right)={P}_{n1}\left(t\right)+{P}_{n2}\left(t\right),$
${P}_{n1}\left(t\right)=\frac{2}{ml}\left[\left({M}_{1}+{M}_{2}\right)g\mathrm{sin}\frac{n\pi vt}{l}-{M}_{1}\sum _{i=1}^{\infty }{\stackrel{¨}{q}}_{i}\left(t\right)\mathrm{sin}\frac{i\pi vt}{l}\mathrm{sin}\frac{n\pi vt}{l}\right],$
${P}_{n2}\left(t\right)=\frac{2}{ml}\left[{k}_{1}Z\left(t\right)+{c}_{1}\stackrel{˙}{Z}\left(t\right)\right]\mathrm{sin}\frac{n\pi vt}{l}-\frac{2}{ml}\sum _{i=1}^{\infty }\left[{k}_{1}{q}_{i}\left(t\right)+{c}_{1}{\stackrel{˙}{q}}_{i}\left(t\right)\right]\mathrm{sin}\frac{i\pi vt}{l}\mathrm{sin}\frac{n\pi vt}{l}.$

The left end of the Eq. (1) can be transformed into:

5
${\stackrel{¨}{q}}_{n}\left(t\right)+2{\xi }_{n}{\omega }_{n}{\stackrel{˙}{q}}_{n}\left(t\right)+{{\omega }_{n}}^{2}{q}_{n}\left(t\right)+\frac{3A}{8I}{{\omega }_{n}}^{2}{{q}_{n}}^{3}\left(t\right).$

Therefore, Eq. (3) can be expressed as Eq. (6):

6
$\left[{\stackrel{¨}{q}}_{n}\left(t\right)+\frac{2{M}_{1}}{ml}\sum _{i=1}^{\infty }{\stackrel{¨}{q}}_{i}\left(t\right)\mathrm{sin}\frac{i\pi vt}{l}\mathrm{sin}\frac{n\pi vt}{l}\right]+\left[2{\xi }_{n}{\omega }_{n}{\stackrel{˙}{q}}_{n}\left(t\right)+\frac{2{c}_{1}}{ml}\sum _{i=1}^{\infty }\stackrel{˙}{{q}_{i}}\left(t\right)\mathrm{sin}\frac{i\pi vt}{l}\mathrm{sin}\frac{n\pi vt}{l}\right]$

Similarly, using the mode decomposition method, the Eq. (2) can be written as:

7
${M}_{2}\stackrel{¨}{Z}\left(t\right)+{c}_{1}\stackrel{˙}{Z}\left(t\right)+{k}_{1}Z\left(t\right)-{c}_{1}\sum _{i=1}^{\infty }{\stackrel{˙}{q}}_{i}\left(t\right)\mathrm{sin}\frac{i\pi vt}{l}-{k}_{1}\sum _{i=1}^{\infty }{q}_{i}\mathrm{sin}\frac{i\pi vt}{l}=0.$

By combining Eq. (6) and Eq. (7), the system dynamic equilibrium equations of the vehicle coupling with the simply-supported beam are obtained. The $N$th order matrix expression of the system motion equation is:

8
$M\left\{\stackrel{¨}{X}\right\}+C\left\{\stackrel{˙}{X}\right\}+K\left\{X\right\}=\left\{F\right\},$

where $\left\{X\right\}={\left[{q}_{1},{q}_{2},\cdots {q}_{N,}Z\right]}^{T}$ represents the generalized displacement vector, $\left\{F\right\}={\left[{\rho }_{F}{\phi }_{1},{\rho }_{F}{\phi }_{2},\cdots {\rho }_{F}{\phi }_{N},0\right]}^{T}$ represents the generalized force vector, $M$ represents the generalized mass matrix, $C$ represents the generalized damping matrix and $K$ represents the generalized stiffness matrix, $H$ is used to describe the nonlinear matrix. This can be described by the following matrices:

$C=\left[\begin{array}{ccc}\begin{array}{cc}2{\xi }_{1}{\omega }_{1}+{\rho }_{c}{\phi }_{11}& {\rho }_{c}{\phi }_{12}\\ {\rho }_{c}{\phi }_{21}& 2{\xi }_{2}{\omega }_{2}+{\rho }_{c}{\phi }_{22}\end{array}& \cdots & \begin{array}{cc}{\rho }_{c}{\phi }_{1N}& {-\rho }_{c}{\phi }_{1}\\ {\rho }_{c}{\phi }_{2N}& {-\rho }_{c}{\phi }_{2}\end{array}\\ ⋮& \ddots & ⋮\\ \begin{array}{cc}{\rho }_{c}{\phi }_{N1}& {\rho }_{c}{\phi }_{N2}\\ {-c}_{1}{\phi }_{1}& {-c}_{1}{\phi }_{2}\end{array}& \cdots & \begin{array}{cc}2{\xi }_{N}{\omega }_{N}+{\rho }_{c}{\phi }_{NN}& {-\rho }_{c}{\phi }_{N}\\ -{c}_{1}{\phi }_{N}& {c}_{1}\end{array}\end{array}\right],$
$K=\left[\begin{array}{ccc}\begin{array}{cc}{{\omega }_{1}}^{2}+{\rho }_{k}{\phi }_{11}& {\rho }_{k}{\phi }_{12}\\ {\rho }_{k}{\phi }_{21}& {{\omega }_{2}}^{2}+{\rho }_{k}{\phi }_{22}\end{array}& \cdots & \begin{array}{cc}{\rho }_{k}{\phi }_{1N}& {-\rho }_{k}{\phi }_{1}\\ {\rho }_{k}{\phi }_{2N}& -{\rho }_{k}{\phi }_{2}\end{array}\\ ⋮& \ddots & ⋮\\ \begin{array}{cc}{\rho }_{k}{\phi }_{N1}& {\rho }_{k}{\phi }_{N2}\\ -{k}_{1}{\phi }_{1}& -{k}_{1}{\phi }_{2}\end{array}& \cdots & \begin{array}{cc}{{\omega }_{N}}^{2}+{\rho }_{N}{\phi }_{NN}& {-\rho }_{k}{\phi }_{N}\\ -{k}_{1}{\phi }_{N}& {k}_{1}\end{array}\end{array}\right],$

where ${\rho }_{M}=2{M}_{1}/ml$, ${\rho }_{C}=2{C}_{1}/ml$, ${\rho }_{k}=2{k}_{1}/ml$ and ${\phi }_{n}=\mathrm{s}\mathrm{i}\mathrm{n}n\pi vt/l$.

It can be seen from the above derivations that in cases where the geometric nonlinearity of the bridge is taken into consideration, the system equation has a non-linear matrix term. It can be found that the non-linear matrix is related only to the nature of the bridge itself, and not to the parameters of the vehicle.

## 3. Establishment of system equation without simplifying wheel acceleration

The external load acting on the bridge is deduced by Eq. (6). In fact, the overall analysis of the vehicle shows that the external loads on the bridge can be expressed as Eq. (9):

9
$P\left(x,t\right)=\delta \left(x-vt\right)\left[\left({M}_{1}+{M}_{2}\right)g-{M}_{1}\frac{{d}^{2}w\left(x,t\right)}{d{t}^{2}}-{M}_{2}\stackrel{¨}{Z}\left(t\right)\right],$
10
$\frac{{d}^{2}w\left(x,t\right)}{d{t}^{2}}=\left(i\right)\frac{{\partial }^{2}w\left(x,t\right)}{\partial {t}^{2}}+\left(ii\right)2\frac{{\partial }^{2}w\left(x,t\right)}{\partial x\partial t}v+\left(iii\right)\frac{{\partial }^{2}w\left(x,t\right)}{\partial {x}^{2}}{v}^{2},$

where ($i$) represents the vertical acceleration of the vibration of the bridge at the position of the wheel, regardless of the change in the position of the load; ($ii$) represents the vertical acceleration of the bridge caused by the movement of the load; ($iii$) represents the vertical acceleration of the bridge caused by the curvature during structural vibration. The external loads on the bridge can be described by Eq. (11), which is calculated by substituting Eq. (10) into Eq. (9):

11
$P\left(x,t\right)=\delta \left(x-vt\right)\left[\left({M}_{1}+{M}_{2}\right)g\right$

According to the mode decomposition method, the right end of Eq. (6) can be obtained and converted into the following equations:

12
${P}_{n}\left(t\right)={P}_{n1}\left(t\right)+{P}_{n2}\left(t\right)+{P}_{n3}\left(t\right),$
13
${P}_{n1}\left(t\right)=\frac{2}{ml}\left[\left({M}_{1}+{M}_{2}\right)g\mathrm{sin}\frac{n\pi vt}{l}-{M}_{1}\sum _{i=1}^{\infty }{\stackrel{¨}{q}}_{i}\left(t\right)\mathrm{sin}\frac{i\pi vt}{l}\mathrm{sin}\frac{n\pi vt}{l}\right],$
14
${P}_{n2}\left(t\right)=-\frac{2{M}_{1}}{ml}2v\frac{n\pi v}{l}\mathrm{cos}\frac{n\pi vt}{l}{c}_{1}{\stackrel{˙}{q}}_{i}\left(t\right)-\left[-\frac{2{M}_{1}}{ml}{v}^{2}{\left(\frac{n\pi v}{l}\right)}^{2}\mathrm{sin}\frac{n\pi vt}{l}{k}_{1}{q}_{i}\left(t\right)\right],$
15
${P}_{n3}\left(t\right)=\frac{2{M}_{2}}{ml}\stackrel{¨}{Z}\left(t\right)\mathrm{sin}\frac{n\pi vt}{l}.$

The system dynamic equilibrium equations of a simply-supported beam, wheel and spring (damper) mass system can be obtained by combining Eq. (5), Eq. (13), Eq. (14) and Eq. (15).

The $N$th-order matrix expression can be expressed as follows:

16
$M=\left[\begin{array}{ccccc}1+{\rho }_{m}{\phi }_{1}{\phi }_{1}& {\rho }_{m}{\phi }_{1}{\phi }_{2}& \cdots & {\rho }_{m}{\phi }_{1}{\phi }_{N}& {{\rho }_{M}\phi }_{1}\\ {\rho }_{m}{\phi }_{2}{\phi }_{1}& 1+{\rho }_{m}{\phi }_{2}{\phi }_{2}& \cdots & {\rho }_{M}{\phi }_{2N}{\phi }_{1}& {{\rho }_{M}\phi }_{2}\\ ⋮& ⋮& \ddots & ⋮& ⋮\\ {\rho }_{m}{\phi }_{N}{\phi }_{1}& {\rho }_{m}{\phi }_{N}{\phi }_{2}& \cdots & 1+{\rho }_{m}{\phi }_{N}{\phi }_{N}& {{\rho }_{M}\phi }_{N}\\ 0& 0& \cdots & 0& {M}_{2}\end{array}\right],$
$H=\frac{3A}{8I}\left[\begin{array}{ccccc}{{\omega }_{1}}^{2}& 0& \cdots & 0& 0\\ 0& {{\omega }_{2}}^{2}& \cdots & 0& 0\\ ⋮& ⋮& \ddots & ⋮& ⋮\\ 0& 0& \cdots & {{\omega }_{N}}^{2}& 0\\ 0& 0& \cdots & 0& 0\end{array}\right],$
$C=\left[\begin{array}{ccccc}2{\xi }_{1}{\omega }_{1}+{\rho }_{m}{\phi }_{1}{\stackrel{˙}{\phi }}_{1}& {\rho }_{m}{\phi }_{1}{\stackrel{˙}{\phi }}_{1}& \cdots & {\rho }_{m}{\phi }_{1}{\stackrel{˙}{\phi }}_{1}& 0\\ {\rho }_{m}{\phi }_{2}{\stackrel{˙}{\phi }}_{1}& 2{\xi }_{2}{\omega }_{21}+{\rho }_{m}{\phi }_{2}{\stackrel{˙}{\phi }}_{2}& \cdots & {\rho }_{m}{\phi }_{1}{\stackrel{˙}{\phi }}_{1}& 0\\ ⋮& ⋮& \ddots & ⋮& ⋮\\ {\rho }_{m}{\phi }_{N}{\stackrel{˙}{\phi }}_{1}& {\rho }_{m}{\phi }_{1}{\stackrel{˙}{\phi }}_{1}& \cdots & 2{\xi }_{N}{\omega }_{N}+{\rho }_{m}{\phi }_{N}{\stackrel{˙}{\phi }}_{N}& 0\\ {-c}_{1}{\phi }_{1}& {-c}_{1}{\phi }_{2}& \cdots & -{c}_{1}{\phi }_{N}& {c}_{1}\end{array}\right],$
$K=\left[\begin{array}{ccccc}{{\omega }_{1}}^{2}+{\rho }_{m}{v}^{2}{\phi }_{1}{\stackrel{¨}{\phi }}_{1}& {\rho }_{m}{v}^{2}{\phi }_{1}{\stackrel{¨}{\phi }}_{2}& \cdots & {\rho }_{m}{v}^{2}{\phi }_{1}{\stackrel{¨}{\phi }}_{N}& 0\\ {\rho }_{m}{v}^{2}{\phi }_{2}{\stackrel{¨}{\phi }}_{1}& {{\omega }_{1}}^{2}+{\rho }_{m}{v}^{2}{\phi }_{2}{\stackrel{¨}{\phi }}_{2}& \cdots & {\rho }_{m}{v}^{2}{\phi }_{1}{\stackrel{¨}{\phi }}_{N}& 0\\ ⋮& ⋮& \ddots & ⋮& ⋮\\ {\rho }_{m}{v}^{2}{\phi }_{1}{\stackrel{¨}{\phi }}_{1}& {\rho }_{m}{v}^{2}{\phi }_{1}{\stackrel{¨}{\phi }}_{1}& \cdots & {{\omega }_{1}}^{2}+{\rho }_{m}{v}^{2}{\phi }_{N}{\stackrel{¨}{\phi }}_{N}& 0\\ -{k}_{1}{\phi }_{1}& -{k}_{1}{\phi }_{2}& \cdots & -{k}_{1}{\phi }_{N}& {k}_{1}\end{array}\right],$
$\left\{F\right\}={\left[{\rho }_{F}{\phi }_{1},{\rho }_{F}{\phi }_{2},\cdots ,{\rho }_{F}{\phi }_{N},0\right]}^{T},$

where:

## 4. Effect of velocity on dynamic response of bridge without simplifying wheel acceleration

For the dynamic response of the wheel mass coupled acceleration of a high-speed vehicle on the bridge, it is necessary to analyze the influence of the last two terms of Eq. (12) ($ii$) and $\left(iii$) on the dynamic response of the vehicle-bridge coupling system [1, 5].

In this section, the influence of wheel acceleration on the bridge’s dynamic response is analyzed without simplification under four speeds: 20 m/s, 25 m/s, 30 m/s and 35 m/s. The calculating model of Eq. (12) is defined as Model A, and the model neglecting ($ii$) and ($iii$) is defined as Model B. Fig. 2 shows the vertical displacement response of the two models at different speeds. The maximum displacement of the structure and the maximum displacement ratio of the two models are shown in Table 1.

Table 1The dynamic responses of the bridge structure

 $v$ velocity(m/s) 20 79.1 74.2 1.0660 25 86.3 74.4 1.1599 30 99.1 73.1 1.3557 35 123.6 70.5 1.7532 ${D}_{Amax}$and${D}_{Bmax}$represent the maximum displacement in the span of Model A and Model B

Fig. 2Comparisons of dynamic response of the mid-span at different speeds a)20 m/s b)25 m/s c)$v=$ 30 m/s d)$v=$ 35 m/s

It can be observed that ${D}_{Amax}$ increases continuously with the increasing velocity. However, ${D}_{Bmax}$ changes slightly with the increase in velocity, varying between 70 mm and 75 mm. With the increase in velocity, the position of the maximum mid-span displacement gradually recedes for both Model A and Model B. When the speed is 35 m/s, the maximum mid-span displacement of Model A and Model B occurs when the vehicle is about to leave the bridge.

When the speed is 20 m/s, Model A and Model B exhibit little difference in their mid-span displacement response. However, when the speed is 30 m/s, the ratio is 1.3557, and model A should be adopted. When the speed is 20 m/s 30 m/s, the appropriate model can be selected according to the calculation accuracy requirements.

## 5. Conclusions

In this paper, the effect of the geometric nonlinearity of bridges on the dynamic response of the vehicle-bridge coupling system is taken into account, and the following conclusions are obtained:

(1) The system equation has a non-linear matrix term when the geometric nonlinearity of the bridge is considered. The non-linear matrix is related only to the nature of the bridge structures, and the parameters of the vehicles have no influence on the matrix.

(2) The value of the non-linear matrix will increase as the inertia moment decreases when the other conditions remain unchanged.

(3) The influence of wheel acceleration on the bridge dynamic response is analyzed without simplification under four speeds. It is acceptable to neglect the influence of wheel acceleration at 20 m/s, while it is necessary to consider the significant influence of the wheel acceleration at 30 m/s.

#### Acknowledgements

This study was funded by the Natural Science Foundation of Chongqing under grant number cstc2018jscx-msyb1299, the Science and Technology Research Program of Chongqing Municipal Education Commission under the grant number KJZD-K201802501, the China Academy of Building Research under the grant number 20190112470730020.

#### Cited by

Stochastic non-stationary characteristics of vehicle-induced bridge vibrations