Free vibration analysis of a cracked simply supported bridge considering bridge-vehicle interaction

This paper presents an analytical approach to investigate the free vibration of simply supported bridge with cracks under arbitrary number of vehicles. Calculation methods for natural frequencies and mode shapes are proposed based on Euler-Bernoulli beam theory, transfer matrix method and numerical assembly method. The vehicle is modeled as a half-car planar model. Equations of motion and displacement functions for bridge and vehicle are derived, respectively. The undermined coefficient matrices for wheels, vehicles and boundary conditions are obtained based on equilibrium and continuity conditions. Numerical assembly technique is adopted to construct the overall matrix of coefficients for bridge-vehicle vibration system. And natural frequencies and corresponding mode shapes are determined based on iterative method and overall matrix solution. Numerical simulation is presented to verify the effectiveness of the proposed method. The results reveal that solutions of the proposed method have favorable reliability. Natural frequencies and associate modal shapes of simply supported multi-girder bridge under the effects of crack and vehicle are investigated. The influences of crack and vehicle parameters on dynamic characteristics are also demonstrated. Meanwhile, a practical simply supported box-girder bridge model is analyzed by the proposed method and an effective crack identification algorithm is proposed.


Introduction
Vibration characteristics of bridge structures under the effects of vehicles have been an interesting topic.Kim [1] investigated the effect of vehicle on variation of dynamic characteristics through traffic-induced vibration tests on three bridges and found that the vehicle could cause the changes of dynamic characteristics.Besides, cracks often appear in bridges under the effects of external loads and environmental factors, which will lead to the variations of dynamic characteristics [2].Therefore, investigation on dynamic characteristics of cracked bridge considering bridge-vehicle interaction is important for health monitoring and condition assessment of bridge structures.
Dimarogonas [3] proposed a state of review on the dynamics of cracked structures.Many works in this fields deal with the cracked beam subjected to various boundary conditions.In some papers, the beam was subdivided into several parts and separated by cracks, which were simulated by massless rotational springs [4,5].The cracks lead to energy release and additional deformation of beam.Based on the simulation of crack using massless rotational spring, many methods including energy approach, transfer matrix method, etc. were proposed.The methods with exponentially decaying crack disturbance functions were proposed to develop vibration equations for continuous models [6][7][8].The exponential function was firstly presented by Christides and Barr to model the stress-strain variation around the crack zone for one or more pairs of symmetric cracks [6].
Saavedra and Cuitino [9] proposed a modeling method for multi-beams systems containing a transverse crack.The flexibility matrix of cracked element given by strain energy density function is directly added to that of corresponding intact element to obtain the total flexibility matrix.However, results of this method are not quite accurate.Transfer matrix method is an efficient tool for free vibration analysis of beams with cracks.This method was firstly introduced by Pestel and Leckie [10].Modified transfer matrix methods are also developed to investigate the dynamic behaviors of beams with various attachments.Ostachowicz et al. [11] presented the relationships of parameters (deformations, moments, slope, etc.) at the position of cracks.Four integration constants of eigenfunctions between adjacent sub-beams can be determined.Lin et al. [12] studied the vibrations of a beam with arbitrary number of cracks using transfer matrix methods.An analytical form of eigenvalue problem was introduced.However, the effect of vehicle has not been taken into account in most of the studies for cracked structure.
Vibration analysis of a beam-like structure subjected to moving loads has been a subject of interest in many fields [13,14].Bridges are generally modeled as elastic beams, while the models for vehicles can be divided into three categories [15]: the so-called moving load [16,17], moving mass [18,19] and moving sprung-mass models [20,21].Esmailzadeh et al. [22] studied the dynamic interaction between high-speed train and simply supported girders based on rigid-body dynamics theory, finite element analysis and wheel-rail displacement corresponding assumption.Song et al. [23] proposed a novel finite element model for simulating the interactions between high-speed train and bridge.Vehicle model devised for high-speed train is employed, which has an articulated bogie system.Nasrellah et al. [24] presented a strategy for structural system identification considering bridge-vehicle interaction by particle filtering algorithms, which can take into account measurement noise, guideway unevenness and incomplete information.Zhang et al. [25] proposed an inter-system iteration method for dynamic analysis of coupled bridge-vehicle system, which can avoid the interaction within time-step and improve computation efficiency.Yang et al. [26] investigated the vertical responses of bridge and moving vehicle using modal superposition, and corresponding closed-form solutions were determined.
From above literature review, one can obtain that most of the researches on vibration analysis of bridge considering bridge-vehicle interaction are focused on dynamic response.Little attention is paid on free vibration characteristics, which are important for health monitoring of bridge.And most of the objects in these researches are undamaged beams, damage of beam is needed to be considered in practical studies.Bilello et al. [27] modeled the damage of beam through rotational springs and carried out theoretical and experimental studies on the response of damaged Euler-Bernoulli beam traversed by a moving mass.Lin and Chang [28] obtained an analytical solution of forced response for a cantilever beam with a crack subjected to a concentrated moving load by using equivalent rotational spring model, transfer matrix method and modal series expansion technique.Yoon and Son [29] investigated the effects of open crack and the moving mass on dynamic behavior of simply supported pipe conveying fluid.However, previous studies still have certain insufficiencies.Vehicle models are relatively simple and could not reflect their vibration characteristics.Therefore, free vibration analysis of simply supported cracked bridges considering bridge-vehicle interaction need to be investigated.
In this paper, free vibration characteristics of cracked bridge considering bridge-vehicle interaction were investigated.Corresponding calculation method is proposed based on Euler-Bernoulli beam theory, transfer matrix method and numerical assembly method.Numerical simulations for simply supported bridges under the effects of crack and bridge are used to verify the reliability and accuracy of the proposed method.Changes of natural frequencies and modal shapes are used to demonstrate the influences of crack and vehicle parameters on simply supported bridge.Moreover, a practical simply supported box-girder bridge model is also analyzed using the proposed method and an algorithm is proposed for crack identification.

Simplified models for bridge and vehicle
Short-to medium-span bridges are the most widely used types for simply supported bridge in structural design.They can be divided into two categories as integral bridge and multi-girder bridge according to the composition forms of girders.For integral bridge, there is only one piece of girder.Its characteristics of cross section are determined by this girder.For multi-girder one, some horizontal connections between girders are adopted to bear loads as shown in Fig. 1.Its characteristics of cross section are determined comprehensively by the girders.Assuming each girder in multi-girder bridge possesses the same geometry dimensions, and also the same crack position and depth.Multi-girder bridge can be simplified as integral one for vertical vibration analysis.Equivalent principles are listed as follows: (a) Height of simplified integral bridge is equal to that of multi-girder one.
(b) Moment of inertia of simplified integral bridge is the sum of values for all girders in multi-girder bridge.It can be expressed as where is the moment of inertia of simplified integral bridge, is the moment of inertia for th girder of multi-girder bridge.(c) Mass per unit length of simplified integral bridge is the sum of values for all girders in multi-girder bridge.It is calculated by: where is mass per unit length of simplified integral bridge, is mass per unit length for th girder of multi-girder bridge.
(d) The position and depth of cracks for simplified integral bridge are the same with multi-girder bridge.
Based on above calculation, dynamic analysis of a multi-girder bridge can be simplified as corresponding integral one.
There are several vehicle models, which possesses two degrees of freedom, four degrees of freedom and six degrees of freedom, respectively.In this study, a half-car planar model with four degrees of freedom is adopted.
, and ( , ) are sprung mass, rotatory mass and wheel mass for the th vehicle, respectively; and are suspension spring constants; and are tyre stiffness coefficients; is distance between the left and right wheels of the th vehicle.

Equation of motion and displacement function for bridge
Taking section for example, corresponding equation of motion and displacement function are established.Section is divided into ( + 1) parts by cracks, which can be denoted by 1, 2,…, , + 1, and 1 , 2 ,…, +1 are the lengths of small parts (shown in Fig. 3).
where is flexural rigidity.is Young's modulus of bridge, is moment of inertia, is mass per unit length.( , ) is vertical deflection of section at position and time .Assuming the whole vibrating system shown in Fig. 2 performs harmonic free vibration at equilibrium position, it has: where is angular frequency for bridge-vehicle system, ( ) is corresponding vibration mode where , , , are undetermined coefficients for section , and = ⁄ .

Equation of motion and displacement function for vehicle
Equations of motion for wheel masses and of the th vehicle are given by: The force balance equations for the th vehicle body (sprung mass, rotatory mass) are given by: The equilibrium equations of motion for the th vehicle is composed of Eqs.(6)(7)., are deformations for wheel masses and , respectively; is deformation for sprung mass; is rotation angle for vehicle body; , are deformations at positions and of bridge, respectively.
Because the whole vibrating system performs harmonic free vibration, it can obtain: where , , and are the amplitudes of , , and respectively.They are undetermined coefficients for amplitudes.

Equations of bridge sections
Taking section for example (shown in Fig. 3), continuity of deformations and equilibrium of moments and forces are satisfied at the right end of th part and left end of ( + 1)th part at the th crack of section .It has: The slope has a jump at the crack location: is called the flexibility function for open crack.For a rectangular-sectional beam, ( ℎ ⁄ ) could be expressed as [30]: in which ℎ is height of cross-section for bridge at crack position, is depth of the th crack in the section .The Eq. ( 10) can also be applicable for other sections, such as box-girder, circular-sectional beam, etc. Corresponding flexibility function in Eq. ( 11) varies according to the type of cross-section.In this paper, a numerical example of box-girder is shown in Section 5.3.
Eq. ( 13) can be changed into: where is transfer matrix for undetermined coefficients at the right end of th part and left end of ( + 1)th part for section : where: The transfer relationship from the 1st part to ( + 1)th part in section can be given by: Assuming: Eq. ( 19) can be expressed by: where is transfer matrix for undetermined coefficients of section .If section has no crack, division is not necessary.One group of underdetermined coefficients is able to represent the mode function.In order to obtain the uniform representation of coefficients in each section, Eq. ( 21) is also used to represent its coefficient equation even if the section of has no crack.Here, has the same meaning with , and is identity matrix.

Equations at wheels
As shown in Fig. 3, sections , and are divided into ( + 1), ( + 1) and ( + 1) parts, respectively.The last parts for , and can be denoted by The substitution of Eq. ( 5) into Eq.( 22), it has: Eq. ( 23) can be represented by matrix form, it is: here: When the left wheel is just at the crack position, , , , and are not equal to zero; on the contrary, , , , and are equal to zero.
The substitution of Eq. ( 21) into Eq.( 24), it has: where = , which is the undetermined coefficient vector for mode shape of the first part of section ; is transfer matrix for undetermined coefficients of section .
At the right wheel of th vehicle ( ), it can also obtain that: Eq. ( 29) can be represented by matrix form, it is: here: When the right wheel is just at the crack position, , , , and are not equal to zero; on the contrary, , , , and are equal to zero.
The substitution of Eq. ( 21) into Eq.( 30), it has: where is transfer matrix for undetermined coefficients of section .

Equations from boundary conditions
At the left end, displacement and moment are zero, and it has: where is mode function for the first part of section .The substitution of Eq. ( 5) into Eq.(40), one obtains: where: , , and are undetermined coefficients for mode shape of the first part of section .
Section is divided into ( + 1) parts, the length for ( + 1)th part is .According to the boundary condition of right end of bridge, it has: The substitution of Eq. ( 5) into Eq.( 45), one obtains: where , , and are undetermined coefficients for mode shape of the ( + 1)th part of section .Eq. ( 46) is represented by matrix form, it obtains: where: The substitution of Eq. ( 21) into Eq.(47), one obtains: GUOJIN TAN, WENSHENG WANG, YUBO JIAO of section ; is transfer matrix for coefficients of section .

Determination of natural frequencies and modal shapes
For cracked simply supported bridge under s vehicles (shown in Fig. 2), the bridge is divided into 2 + 1 sections by vehicles, which can be represented by , , , … , , , , ,…, , .At one section, undetermined coefficients of mode shape for each part can be represented by four coefficients of the first part (shown in Eq. ( 21)).Therefore, there are four undetermined coefficients for each section.For 2 + 1 sections, there are totally 4 × (2 + 1) coefficients.Moreover, there are four undetermined coefficients for each vehicle.The whole system possesses 4 × (2 + 1) + 4 × = 12 + 4 undetermined coefficients.
There are four equations for undetermined coefficients at the left wheel of vehicle (shown in Eq. ( 28)); and four equations for undetermined coefficients at the right wheel of vehicle (shown in Eq. ( 34)).Four equations for coefficients can be obtained by equation of motion of vehicle (shown in Eq. ( 36)).Therefore, there are totally twelve equations for each vehicle, and 12 equations for vehicles.At the left end and right end of bridge, there are two equations for undetermined coefficients, respectively (shown in Eqs.(42, 50)).In general, there are totally 12 + 4 equations for undetermined coefficients.Numerical assembly method is adopted to obtain the matrix equation of all undetermined coefficients, one has: where: (53) Non-trivial solution of Eq. ( 52) requires that: Half-interval method is used to determine the natural frequencies ( = 1, 2, ...) of cracked The accurate values of are obtained respectively using the half-interval method.The substitution of the obtained frequency into Eq.( 5) will determine the corresponding mode shape of the bridge.

Cracked simply supported beam without vehicle
In order to illustrate the proposed method, numerical example will show how to determine the natural frequencies of a simply supported beam with uniform cross-section and with cracks.
The first two order natural frequencies are obtained through the proposed method in this paper.Corresponding reduction coefficient is calculated by Eq. (57), which can be used to demonstrate the effect of crack on natural frequencies of simply supported beam without vehicle: As can be seen from Fig. 4, and decrease with the increasing of crack width for simply supported beam with the same crack location.It indicates that the bigger the crack depth is, the greater the frequency reduction is.As for the effect of crack location, and present different variations.For , the minimum value is at the mid-span ( = 400 mm) for simply supported beam with the same crack depth.For , the minimum value is at the 1/4 span ( = 200 mm).The phenomena coincide with the change rules for 1st and 2nd mode shapes.The results are consistent with those obtained by Liang et al. [31], which can validate the effectiveness of the proposed method for dynamic analysis of simply supported cracked beam without vehicle.

Intact simply supported beam with vehicle
Two numerical examples are used to verify the accuracy of proposed method for simply supported beam with vehicles through the comparison of results calculated by proposed method and FEM (shown in Fig. 5).For finite element model, the beam is divided into 40 beam elements with length 0.5 m and each element has two nodes, in which each node has rotational and vertical displacements.The mass matrix and stiffness matrix of beam elements are formed using Lagrange interpolation function, and equation of vibration is derived as follows: where = , = , is the 2nd order derivation of modal coordinates for beam, is the first th order mode shape matrix of free vibration for beam, is the beam-vehicle interaction force: here, is the location matrix of external force point; and are the forces of front and back wheel on the beam, respectively.Eq. ( 59) is substituted into Eq.( 58).Equation of beam-vehicle system can be obtained combining with Eqs.(6-7), which is shown as follows [32]: where and are the mass and stiffness matrices of vehicle, respectively.is coupling stiffness matrix of the beam-vehicle interaction system.Based on the finite element program constructed using MATLAB, modal properties of beam-vehicle system can be obtained by solving Eq. (60).
For the parameters of vehicle model, wheel masses = = 1500 kg, sprung mass = 1.77×10 4  The first three natural frequencies for cases 1 and 2 are calculated by the proposed method and FEM, respectively.The results are listed in Tables 1-2.As can be seen from the calculation results, relative errors for the first three natural frequencies for cases 1 and 2 are lower than 0.21 %.It reveals that the proposed method possesses favorable accuracy for vibration analysis of simply supported beams with vehicle.

Simply supported bridge model
A simply supported bridge with uniform cross-section and 30 m span is adopted for numerical simulation and the bridge is composed of 4 girders (shown in Fig. 6(a)).According to the equivalent principle listed in section 2.1, cross section for this model can be simplified as a rectangular with width 0.8 m, height 1.5 m (shown in Fig. 6(b)).Elastic modulus is 3×10 10 Pa, material density is 2500 kg/m 3 .The first three natural frequencies are calculated by the proposed method and listed in Table 3.

Cracked simply supported bridge without vehicle
Cracked simply supported bridge without vehicle is shown in Fig. 7. Its cross section is consistent with Fig. 6(b) and crack locations are 2, 4, …, 28 m from the left end of bridge, respectively.In order to demonstrate the effects of crack depth, crack position and crack number on natural frequencies of cracked simply supported bridge without vehicle, change ratio of natural frequencies considering cracks (∆ ) is used, which can be calculated by: where and are natural frequencies calculated by the proposed method for cracked and intact simply supported bridge without vehicle, respectively.In order to demonstrate the effect of crack depth, the case of crack at position (1) (mid-span) with ratio of crack depth ( ) to cross section height (ℎ) ( /ℎ) varying from 10 % to 70 % (shown in Fig. 7) is adopted.The first three natural frequencies of bridge are calculated by the proposed method.Corresponding results are shown in Fig. 8.
In order to demonstrate the effect of crack position, the cases of crack with /ℎ = 15 % located at 2 m, 4 m, …, 26 m, 28 m from the left end of the bridge are used, respectively.Corresponding results are shown in Fig. 9.
As for the effect of crack number, /ℎ = 15 %, number of crack is 1 to 15, respectively.The results are shown in Fig. 10.
As can be seen from Fig. 8, ∆ of the first and third frequencies decrease nonlinearly with the increasing of crack depth.The greater crack depth is, the larger ∆ presents.However, crack depth has little influence on the second order frequency because crack is located at zero point for this mode.
In Fig. 9, ∆ is the largest if crack position is closer to the antinodal points of corresponding mode shape.Meanwhile, when the crack is located at the zero point of a certain mode, ∆ of corresponding frequencies are close to zero, which indicates that the crack has little influence.As can be seen from Fig. 12, mode shapes of cracked bridge have significant changes at the crack positions, which is because the slopes of mode shapes at crack positions are discontinuous.Mode shapes calculated by the proposed method are coincident with the factual data, and it can be used for damage detection of bridge.

Simply supported bridge with one vehicle at different positions
Different vehicle parameters are determined and listed in Table 5.For these vehicles, each has four freedoms.Therefore, four natural frequencies ( , , and ) can be calculated for each vehicle model.Calculation results of natural frequencies for vehicles are listed in Table 6.As for bridge model, it is consistent with Fig. 6(b) (shown in Fig. 13).The first natural frequencies are calculated by the proposed method for cases that the centerline of vehicle moves from support A to B (shown in Fig. 13).Change ratio of the first natural frequency considering vehicle (∆ ) is used to demonstrate the effect of vehicle, which can be calculated by: where is the first natural frequency of simply supported bridge with one vehicle calculated by the proposed method, is the first natural frequency of intact bridge without vehicle calculated by the proposed method.
Results of ∆ for vehicle moving from A to B are shown in Fig. 14.As can be seen from Fig. 14, ∆ > 0, when vehicle 1 or 2 is used.It indicates that the first natural frequency of simply supported bridge under vehicle 1 or 2 is higher than that without vehicle.However, ∆ 0 , when vehicle 3 or 4 is adopted.For all vehicle parameters, = 16.4493rad/s, which is between and .When parameters of vehicle 1 or 2 is adopted, ( − ) ( − ); while ( − ) > ( − ) for vehicle 3 or 4. Based on above analysis, one conclusion can be obtained for assessing relative size of and .The closest vehicle frequency ( ) with should be determined firstly.Then, and are compared.If > , ; if , > .Other cases are also simulated to verify above results and that is effective.Moreover, the influence of vehicle on natural frequency decreases with the increasing of distance between the vehicle centerline and mid-span.

Cracked simply supported bridge with one vehicle at different positions
Model for cracked simply supported bridge with one vehicle at different positions is shown in Fig. 15.The first natural frequencies are calculated for cracked simply supported bridge with one vehicle in which the crack is located at mid-span and centerline of vehicle moves from A to B. In order to demonstrate the effect of bridge-vehicle interaction, change ratio of natural frequency considering crack and vehicle (∆ ) is calculated by: where is the first natural frequency of cracked simply supported bridge with one vehicle calculated by the proposed method.
In Fig. 16, coupled effects of vehicle and crack are presented for ∆ of simply supported bridge.Increasing of crack depth leads to the reduction of ∆ regardless of vehicle 1 or 3.If vehicle 3 is used, ∆ 0, ∆ 0 and ∆ 0. If vehicle 1 is adopted, ∆ 0, ∆ > 0 and ∆ (> or ) 0. As can be seen from Fig. 16, there are specific points for ∆ = 0 when crack depth /ℎ = 15 %, 20 % and 25 %, respectively.Therefore, vehicle effect must be considered for crack identification of bridge.
In order to demonstrate the effect of vehicle on crack identification, change ratio of the first natural frequency (∆ ) for simply supported bridge under vehicle 1 at mid-span, under vehicle 1 and with crack depth /ℎ = 15 % at mid-span, only with crack depth ℎ ⁄ = 15 % at mid-span are calculated by the proposed method, respectively.Corresponding results are shown in Fig. 17.
As shown in Fig. 17, difference of ∆ caused by ( 1) and ( 2) (( 1)-( 2)) is equal to that of 15 % crack depth without vehicle.Therefore, bridge-vehicle vibration system can be considered as a linear one, and crack-induced ∆ is comparable for bridge identification under the same parameters of vehicle.The local flexibility due to a sectional crack can be calculated according to reference [33].
For a shallow open crack, the penetration depth ' ' is contained within the bottom solid-sectional region [0 ≤ ≤ ]: For a deeper open crack, depth ' ' goes into the middle hollow-sectional region [ ≤ ≤ (ℎ − )]: where = /ℎ , = ℎ ⁄ , = ℎ ⁄ , and is a function of the relative depth : Substituting the local flexibility coefficient ' ' into Eq.(10), transfer matrix can be calculated by the proposed method in this paper.Therefore, natural frequency and modal shape can be further obtained.

Cracked simply supported box-girder bridge without vehicle
Cases of various crack depths, positions and numbers are studied.Changes of crack depths, positions and numbers are the same as Section 5.2.2 and the results of ∆ are shown in Figs.19-21.As can be seen from Figs. [19][20][21], relationships between ∆ and crack depths, positions, numbers are similar with the conclusions obtained in Section 5.2.2.It reveals that the proposed method is general for different bridge type.For ∆ in Fig. 19, mutation point can be found due to crack going into the middle hollow-sectional region.As shown in Fig. 24, the variation laws of normalized mode shapes are similar to that in Section 5.2.3 and there are also significant changes at crack positions.

Diagnosis of crack for a cracked simply supported box-girder with one vehicle
As discussed above, presence of cracks changes the characteristic equation and has influences on mode shapes and natural frequencies of the beam.In an inverse problem, crack positions and depths are unknown parameters, however, natural frequencies and mode shapes of damaged beam could be obtained through various testing techniques, and the 1st natural frequency and corresponding normalized mode shape can easily be determined experimentally.Because no experiment is carried out in this paper, a numerical example shown in Fig. 25 is used to identify the crack positions and depths.The 1st natural frequency ' ' and corresponding normalized mode shape ' ' in Section 5.3.3 are adopted as input data for crack identification.
is the random generator function in MATLAB with a zero mean and a standard deviation of , is error lever.
According to the matrix equation of all undetermined coefficients in Eq. ( 53), = = 1, one can obtain the assembly matrix [ ] with dimension of 16×16.The undetermined coefficients correspond to , , and vehicle, respectively.Finally, there are 16 equations and 16 undetermined coefficients.Assuming crack positions and depths are known and is a constant ' ', other undetermined coefficients could be solved by using any other fifteen equations.The undetermined coefficients of and are obtained by Eq. ( 21), the substitution of the undetermined coefficients and ' ' into Eq.( 5) will determine the normalized mode shape of the bridge ' '.The iterative algorithm is used to obtain the crack positions and depths by comparing normalized ' ' and ' ', and calculation process for this algorithm is as follows.(b) According to ' ', ' ' and ' ', the normalized mode shape ' ' can be obtained, ∆ = − , note that all normalized mode shapes ' ' and ' ' are consist of fixed position points; (c) Sensitivity matrix [ ] can be calculated by Perturbation method [34]: where: (d) Solve the residuals Δ and Δ : where [ ] is the generalized inverse matrix of [ ] using singular value decomposition.
(e) Update the crack parameters: where is an underrelaxation parameter and can be used to avoid overshoots and increase reliability of the algorithm in early steps of iterations.
(f) Iterate the procedures (b)-(e) until the residuals of Δ and Δ become sufficiently small and satisfactorily.
The actual and predicted crack parameters calculated from the proposed identification algorithm are listed in Table 9.From Table 9, the maximum relative error for crack position identification is 0.45 %, and it is 0.41 % for crack depth identification when = 0 %.Relative error increases with increasing of , and it is within 5.5 % for both position and depth identifications when = 10 %.The results reveal that the proposed algorithm for crack identification possesses favorable accuracy and robustness to noise.

Conclusions
In this paper, free vibration analysis of cracked simply supported bridge under arbitrary number of vehicles is presented.The exact solutions for natural frequencies and mode shapes of cracked bridge considering bridge-vehicle interaction are derived based on Euler-Bernoulli beam theory and numerical assembly method.Numerical simulation on simply supported bridge is used to verify its feasibility.Following conclusions can be obtained: 1) The first two natural frequencies of simply supported bridge with different crack depths are calculated by the proposed method.Reduction coefficients for natural frequencies are consistent with the present study.The first three natural frequencies of simply supported bridge under one and two vehicles are analyzed, respectively.The results are in good agreement with FEM.It reveals that the proposed method possesses favorable reliability.
2) Simply supported multi-girder bridge is adopted and simplified.Change ratio of natural frequencies considering cracks (∆ ) is presented to investigate the influences of crack depth, position and number, respectively.For crack depth, the greater depth is, the larger ∆ presents for the first and third frequencies.However, it has little influence on the second order frequency.For crack position, ∆ is the largest if crack position is closer to the antinodal points of corresponding mode shape.For crack numbers, ∆ decreases with the increasing of crack number.Moreover, significant changes of modal shapes occur at the crack positions.
3) The influences of vehicle parameters on modal properties are also investigated.For intact bridge under one vehicle with four different parameters, the first natural frequencies present different change rules, which are closely associated with vehicles parameters.Corresponding judgment method is also proposed based on frequencies of vehicle and bridge without vehicle.Moreover, the influence of vehicle on natural frequency decreases with the increasing of distance between the vehicle centerline and mid-span.For cracked simply supported bridge with one vehicle at different positions, coupled effects of vehicle and crack are presented for frequency changes, and vehicle effect must be considered for crack identification of bridge.Meanwhile, bridge-vehicle vibration system can be considered as a linear one, and crack-induced frequency changes are comparable for bridge identification under the same parameters of vehicle.
4) The proposed method is applicable for simply supported box-girder bridge.An effective approach for crack identification is proposed, which presents favorable accuracy.

Fig. 2 .
Fig. 2. Simply supported bridge with cracks and under multiple vehicles

Fig. 3 .
Fig. 3. Sketch for division of sections , and Assuming = − ( + ∑ ); for the section , 0 ≤ ≤ .Equation of motion for section can be established based on Euler-Bernoulli beam theory, which is given by: represent the mode functions of the first part for sections , and ., and are used to represent the mode functions of the last part for sections , and .Continuity of deformations and equilibrium of moments and forces are satisfied at the position of left wheel of the th vehicle ( ), it has: are vehicle numbers.In the process of assembling matrix [H], position for each element can be determined by identification number shown at the right side and top of matrix shown in Eq. (53).Identification numbers for undetermined coefficients are shown in Eq. (54).

4 .
57) where and are the th natural frequency of the cracked and un-cracked beams, respectively, is the reduction coefficient of the th natural frequency.In this paper, four different crack depths ( = 10, 15, 20 and 25 mm) at four different locations ( = 100, 200, 300 and 400 mm, is the distance between crack location and the left end of the beam) are investigated.Corresponding and are calculated and shown in Fig. 4. The results obtained by Liang et al. [31] are also presented in Fig. 4, which are denoted by "•".a) Reduction coefficient for 1st natural frequency ( ) b) Reduction coefficient for 2nd natural frequency ( ) Fig. Crack-induced eigen-frequency changes for various crack locations and depths

6 .
a) Sketch of cross section for simply supported bridge b) Simply supported bridge model and cross section Fig. Simply supported bridge for numerical simulation (unit: m)

Fig. 8 .Fig. 9 .
Fig. 8. Relationships between ∆ and crack depths for the first three natural frequencies

Fig. 10 .
Fig. 10.Relationships between ∆ and crack numbers for the first natural frequency 5.2.3.Cracked simply supported bridge with one vehicle Cracked simply supported bridge with one vehicle is shown in Fig. 11, in which cross section is consistent with Fig. 6(b).For the parameters of vehicle model, wheel masses = = 1500 kg, sprung mass = 1.7735×10 4 kg, rotatory mass = 2.4×10 5 kg•m 2 , suspension spring constants = = 2.0×10 6 N/m, tire stiffness coefficients = = 1.4×10 6 N/m, wheel distance = 4.2 m and = = 0.5.Crack positions are shown in Fig. 11, and crack depth is the same, /ℎ = 30 %.The first three natural frequencies for intact bridge without vehicle and cracked bridge with vehicle are calculated by the proposed method and listed in Table 4. Corresponding normalized modal shapes are shown in Fig. 12.

Fig. 18 .
Simply supported box-girder bridge: a) longitudinal section; b) cross section

Fig. 21 .
Fig. 21.Relationship between ∆ and crack numbers for box-girder bridge 5.3.3.Cracked simply supported box-girder bridge with one vehicle The parameters of vehicle are the same with Section 5.1.2.The depths of cracks are 30 %. Layouts of vehicle and cracks are shown in Figs.22-23.The first three natural frequencies calculated by the proposed method are listed in Table 8, and normalized modal shapes are shown in Fig. 24.
2148.FREE VIBRATION ANALYSIS OF A CRACKED SIMPLY SUPPORTED BRIDGE CONSIDERING BRIDGE-VEHICLE INTERACTION.
GUOJIN TAN, WENSHENG WANG, YUBO JIAO 2148.FREE VIBRATION ANALYSIS OF A CRACKED SIMPLY SUPPORTED BRIDGE CONSIDERING BRIDGE-VEHICLE INTERACTION.GUOJIN TAN, WENSHENG WANG, YUBO JIAO simply supported bridge under multiple vehicles.For each order of natural frequency, it satisfies Eq. (56).Mode shapes can be obtained by substituting natural frequencies ( = 1, 2, ...) into Eq.(5).The detailed process for determining lowest natural frequencies is as follows: Firstly, an initial value of frequency coefficient is assumed, the values of coefficient matrix | | corresponding to are obtained and the values of | | are calculated, = .Then a new value = + ∆ with ∆ (e.g.∆ = 0.5) representing the increment of is assumed, and the same calculations are repeated to determine the new value of determinant corresponding to , ] is obtained.The same calculations are repeated to determine the new interval [ , ] until the rank of matrix | | is equal to − 1.
) 2148.FREE VIBRATION ANALYSIS OF A CRACKED SIMPLY SUPPORTED BRIDGE CONSIDERING BRIDGE-VEHICLE INTERACTION.GUOJIN TAN, WENSHENG WANG, YUBO JIAO

Table 1 .
Calculation results of natural frequencies for case 1

Table 2 .
Calculation results of natural frequencies for case 2

Table 3 .
Natural frequencies for simply supported bridge without vehicle

Table 4 .
The first three natural frequencies for intact bridge without vehicle and cracked bridge with vehicle

Table 5 .
Four different vehicle parameters

Table 6 .
First four natural frequencies of vehicle models

Table 9 .
Crack identification of a simply supported bridge with vehicle and double cracks Crack parameters Actual Initial value Predicted Relative error (%)