Research onBridge Structure Reliability Evaluation due toVessels Collison Based on a Statistical Moment Method

Damage to bridge structures caused by vessel collision is a risk for bridges crossing water traffic routes. ,erefore, safety around vessel collision of existing and planned bridges is one of the key technical problems that must be solved by engineering technicians and bridge managers. In the evaluation of the reliability of the bridge structure, the two aspects of vessel-bridge collision force and structural resistance need to be considered. As there are many influencing parameters, the performance function is difficult to express by explicit function. ,is paper combines the moment method theory of structural reliability with finite element analysis and proposes a statistical moment method based on finite element analysis for the calculation of vessel-bridge collision reliability, which solves the structural reliability problemwith a nonlinear implicit performance function. According to the probability model based on current velocity, vessel velocity, and vessel collision tonnage, the estimate points in the standard normal space are converted into estimate points in the original state space through the Rosenblatt reverse transform. According to the estimate points in the original state space and the simplified dynamic loadmodel of vessel-bridge collision, the sample time-history curve of random vessel-bridge collision force is generated, the dynamic response of the bridge structure and the structural resistance of the bridge are calculated by establishing a finite element model, and the failure probability and reliability index of the bridge structure is calculated according to the fourth-moment method. ,e statistical moment based on the finite element analysis is based on the finite element analysis and the moment method theory of structural reliability. ,e statistical moment of the limited performance function is calculated through a quite small amount of confirmatory finite element analysis, and the structural reliability index and failure probability are obtained. ,e method can be widely used in existing finite element analysis programs, greatly reducing the number of finite element analyses needed and improving the efficiency of structural reliability analysis.


Introduction
According to data [1], vessel-bridge collision accidents occur every year around the world and bridge structures are seriously damaged or even collapse, such as the Sunshine Bridge and Queen Isabella Bridge in the United States and the Tasman Bridge in Australia [2]. erefore, vessel-bridge collision is a factor that cannot be ignored in bridge design, and it is necessary to calculate the reliability of the bridge structure when bridge-vessel collision occurs. For the calculation of the reliability of the bridge structure, it is usually necessary to consider two aspects: vessel-bridge collision force and structural resistance [3]. Because of its many influencing parameters, it is difficult to display the performance function. In response to this problem, this paper combines the moment method theory of structural reliability and finite element analysis, according to the probability model of each influencing factor that affects the failure probability, and the finite element-statistical moment method is used to solve the structural reliability problem with a nonlinear implicit performance function.
Since the 1960s, bridges around the world have collapsed or been seriously damaged by vessel collision. Vessel collision accidents not only cause huge economic losses and casualties but can also have negative political impacts and cause serious environmental losses. e safety assessment of vessel collision of existing and planned bridges is one of the key technical problems that must be solved carefully by engineering technicians and bridge managers.
Reliability analysis of bridge structures under vessel collision is an important aspect of vessel-bridge collision safety assessment. e calculation of vessel collision reliability of a bridge structure involves two aspects: vessel collision force and structural resistance.
ere are many influencing factors, and the performance function is difficult to express by explicit function. At present, the methods of structural reliability analysis mainly include the first-order second-moment central point method [4], JC method [5], Monte Carlo method [6], and so on. e traditional structural reliability methods need to expand the performance function. It usually requires complex partial derivative calculation to solve the checking point by iterative calculation, which is difficult to use for structural reliability analysis with implicit performance function.
e Monte Carlo method has the disadvantages of requiring large amounts of calculation and being time-consuming. e AASHTO vessel collision design guidelines of the United States [7] provide a method to determine the bridge collapse probability PC by using the ratio of the limit lateral resistance H of the pier to the vessel collision force P, as shown in Figure 1. In addition, Proske and Curbach [8], Manuel et al. [9], Geng [10], Consolazio et al. [11], and other scholars used the first-order second-moment method, Latin hypercube sampling Monte Carlo method, and Markov chain Monte Carlo method to study the vessel collision reliability of bridge structures, mainly using Monte Carlo simulation based on deterministic finite element analysis, which requires large amounts of calculation and is time-consuming. In this paper, the vessel collision reliability of a bridge structure is analyzed by combining the structural reliability moment method theory and finite element analysis.  [12]; therefore, we consider the following: When the pier column undergoes shear failure, the limit state performance function of the bridge pier column can be expressed as

Limit
When the pier column undergoes bending failure, the limit state performance function of the bridge pier column can be expressed as In the above two formulas, V represents the vessel collision velocity, DWT represents the vessel collision deadweight tonnage, and θ represents the vessel collision angle.
(2) Failure Mode of Pile Group Foundation. For pile group foundations, according to the deformation failure criterion, in order to ensure that the pile group foundation is not destroyed, the rotation θ S(Pile) and limited values of component rotation θ R(Pile) of the most dangerous pile of the pile group foundation under the vessel collision need to meet the following requirements [13]: (3) us, the limit state performance function of the bridge pier column can be expressed as For a single pile-column pier system composed of pile group, pile caps, and piers (or tower columns), the failure of any one part will lead to the failure of the entire pier system. erefore, the entire pier system forms a series system; then the limit state performance function of a single pile-column pier under vessel collision can be expressed as

Open-End Caisson Foundation (Sunk Shaft Foundation).
e open-end caisson foundation sustains stability damage under vessel collision, mainly in the form of overturning and sliding. erefore, the damage index can be expressed by the antioverturning stability coefficient K 0 and the antisliding stability coefficient K c .
When the open-end caisson foundation undergoes overturning failure, the corresponding limit state performance function can be expressed as When the open-end caisson foundation undergoes sliding failure, the corresponding limit state performance function can be expressed as 2 Advances in Materials Science and Engineering In the above two formulas, [K 0 ] represents the permissible antioverturning stability coefficient, and [K c ] represents the permissible antisliding stability coefficient. us, the limit state performance function of the openend caisson foundation under the vessel collision can be expressed as e limit state performance function of pile-column pier and open-end caisson foundation is unclosed to basic random variables; that is, the performance function is implicit, and the performance function itself has nonlinear characteristics. Because the functional relation of the bridge structure's response to basic variables under vessel collision is relatively complicated, the function in the result is often difficult to express explicitly [14].

Limit State Performance Function of Bridge Structure
System. When any pier in the bridge structure system reaches the corresponding performance level state, the entire bridge structure system is considered to have reached the corresponding performance level state [8]. erefore, the entire bridge structure can be regarded as a series system, so the calculation of the failure probability of vessel-bridge collision can be simplified to the calculation of the failure probability of a series structure system. If a bridge system has k piers, the failure probability of the entire bridge system under vessel collision force can be expressed as P F � Prob min g 1 , g 2 , . . . , g k > 0 . (9) erefore, the limit state performance function of the entire bridge structure system under the action of vessel collision force can be expressed as the minimum value of all the bridge pier performance functions; that is, rough the performance function of the bridge structure system established above, the reliability of the vessel-bridge collision system of the bridge structure under the action of vessel collision can be solved.
Aiming at the problem of vessel-bridge collision reliability analysis and solving failure probability of implicit performance function, this paper proposes a statistical moment method for calculating failure probability of vesselbridge collision based on finite element analysis, which is especially suitable for addressing the problem that the performance function cannot be expressed explicitly.

Probability Model of Random Factors in Vessel-Bridge Collision Performance Function
is paper selects the main factors that have a greater impact on the failure probability of vessel-bridge collision, including vessel collision tonnage, collision velocity, and collision angle.

Probability Model of Vessel Tonnage.
e determination of vessel collision tonnage firstly depends on the determination of the deadweight tonnage of each type of vessel. For each type of vessel, its deadweight tonnage is not a fixed value but is distributed within a certain range.
ere is no unified view on the probability distribution of each type of vessel tonnage interval. is paper uses an unknown distribution model to simulate the distribution of vessel collision tonnage. Since the statistical probability characteristics of vessels passing through the bridge area can be described by statistical moments (mean value, standard deviation, skewness, and kurtosis), the probability distribution can be described by the first four order moments of statistical variables. e probability density function can be expressed as e distribution function is Among them μ is the mean value of the statistical variable; σ is the standard deviation of the statistical variable; and a 1 , a 2 , a 3 , a 4 are all coefficients, which are obtained from the first four order moments of the statistical variable.
Using this method to describe the tonnage of navigable vessels and directly reflecting the statistical characteristics of random variables according to the first four order statistical moments of random variables is more accurate.

Probability Model of Collision Velocity.
e collision velocity of a vessel is one of the important factors that affect the collision force of a vessel. Typical upward and downward navigation velocities of various representative vessels are obtained by investigating vessel flow [15]. e AASHTO [16] specification adopts the broken line form shown in Figure 2 for the distribution of vessel collision velocity.
ere are two random variables affecting vessel collision velocity: e navigation velocity of the vessel and the mean value of current velocity are all subject to normal distribution, so the vessel collision velocity also follows a certain state distribution. is paper uses normal distribution to describe its probability characteristics.
Under the condition of x ≤ x c and x > x L , the mean value and standard deviation of the collision velocity are consistent with the mean value and standard deviation of V T and V MIN .
Under the condition of x c < x < x L , the average collision velocity can be calculated as follows: Advances in Materials Science and Engineering e standard deviation can be calculated as follows:

Probability Model of Collision
Angle. e angle of the vessel impacting a pier is not only related to wind direction, water flow, and other factors, but also related to the geometric shape of the bridge substructure and the relative channel position. Due to the lack of existing observation data on vessel yaw angle, the probability model of collision angle is described differently.
is paper uses extreme I-type distribution to describe the probability model of collision angle, the mean value of collision angle on the straight route can be taken as 10-15°, and the standard deviation is 4-6°on the nonstraight route, the mean value should consider the angle between the channel and the normal direction of the bridge axis [17]. e probability density function of extreme I-type distribution is In the formula,

Calculation Model of Vessel-Bridge Collision Dynamic Response
To simplify the analysis and calculation method, when analyzing the dynamic response of vessel-bridge collision, the vessel collision can be applied to the finite element model of bridge structure as a forced dynamic load [18], and the time-history response of the structure can be obtained by solving the basic dynamic equation. e analysis model is shown in Figure 3: e vessel's forced dynamic load model adopts the modified half-wave sinusoidal load model. e time-history curve of vessel collision force is calculated according to the following formula: In the formula, F(t) refers to the time course of impact force (MN); I refers to the momentum of vessels (MNS); T refers to the duration of vessel collision (s); t refers to the time (s); and k, m, and n are dimensionless model parameters.

Statistical Moment Method for Calculating
Failure Probability of Vessel-Bridge Collision e probability distribution of the performance function is a function of random variables. e shape and position of the distribution curve of the probability density function of the performance function are closely related to the moment of the performance function [19]. e failure probability of the structure can be calculated by finding the relationship between the failure probability and the central moment of the performance function.
erefore, the first several order central moments of the performance function need to be solved.

Point Estimate of Statistical Moments of Performance
Function. For the point estimate of performance function with multiple random variables, when calculating the moment of performance function g(X), the performance function can be replaced by the following function formula: Since u i is an independent random variable and g i is only a function of the random variable, u i is also independent, and the first four order statistical moments of g i can be calculated by the following formula:

Fourth-Moment Reliability Index Based on Pearson
System. For standardized random variable Z s , the probability density function satisfies the following partial differential equation: In the formula, respectively, a, b, c, and d are erefore, the reliability index based on the fourth moment can be expressed as

Statistical Moment Method for Vessel-Bridge Collision Failure Probability Calculation Based on Finite Element Analysis
For the nonlinear implicit performance function, the traditional reliability analysis method cannot be solved. In this paper, the structural reliability moment method is combined with the finite element method, and the finite element analysis method is used to calculate the dynamic response of the vessel collision and the performance index boundary value of the structure. rough the finite element calculation, the failure probability of the structure is obtained, which greatly improves the calculation efficiency of the failure probability analysis. e calculation process is shown in Figure 4 below. e concrete ideas are as follows: (1) In view of the factors affecting the failure probability of vessel-bridge collision, vessel tonnage, collision velocity, and collision angle, probability models are selected, respectively, and the limit state performance function of vessel-bridge collision is established according to the analysis model. (2) According to the probability model of each influencing factor, the estimate points in the standard normal space are converted into estimate points in the original state space through the Rosenblatt reverse transform. According to the estimate points in the original state space and the simplified dynamic load model of vessel-bridge collision, the sample time-history curve of random vessel-bridge collision force can be generated. (3) Establish the finite element model of bridge structure. According to the sample time-history curve of vessel collision force, the dynamic response of bridge under vessel collision is solved by finite element analysis. (4) is paper considers that when the most unfavorable angle of pile top is greater than the limit angle, the bridge structure is destroyed. According to the actual reinforcement at the top of the pile foundation, the limit value of the performance index of the bridge structure, namely, the limit rotation angle, is calculated by finite element analysis.  Advances in Materials Science and Engineering the box girder are rigid, and the foundation is rock-socketed pile; No. 3 transition pier is a single-row double-column pier, the foundation is rock-socketed pile, and each pier is connected by a tie beam.

7.2.
e Finite Element Model of Bridge Structure. Midas/Civil was used for finite element modeling analysis. e main beam, pier, cover beam, tie beam, cap, and pile foundation were simulated by spatial beam element. According to the specific form, the bridge bearing is simulated by coupling the degree of freedom between the node on the pier or tower and the corresponding node on the main beam, and the bottom of the cap and pile foundation adopt a master-slave constraint; the main beam and the middle-pier are rigidly connected, the main beam and the side pier are elastically connected, and the pile foundation is embedded at 4 times the pile diameter below the scour line. e model takes the longitudinal direction as the X axis, the transverse direction as the Y axis, and the vertical direction as the Z axis (Figure 7). e corresponding relationship between each component and the model is as follows:    6 Advances in Materials Science and Engineering e material parameters of the bridge structure are shown in Table 1 (the plastic hinge is set at the top of the pile foundation, and the limit rotation angle is determined by the M-φ curve according to the actual reinforcement of the section).

Sample Value of Vessel Collision Velocity.
In this paper, two random variables of collision tonnage and collision velocity are considered as influencing factors of vessel collision force. As shown in Figure 4, vessel collision velocity is determined by two random variables of typical navigation velocity and average current velocity. According to the actual hydrological situation of the bridge, the probability models and parameters of the two random variables of typical navigation velocity and average current velocity are shown in Table 2. According to the mean value and coefficient of variation of these two random factors, combined with the structural layout and navigation requirements of the bridge, the mean value of collision velocity is 3.47 m/s, and the coefficient of variation is 0.043. According to the above theoretical analysis, this paper uses five-point estimation, and the estimation points in the standard normal space are shown in Table 3. e collision velocity in this paper adopts the normal distribution. According to the estimation points in the standard normal space and the mean value and coefficient of variation of the collision velocity, five grades of collision velocity are obtained through the Rosenblatt inverse transform, which are, respectively, 3

e Sample Time-History Curve of the Vessel Collision
Force. In this paper, the vessel tonnage adopts unknown distribution, and the impact velocity adopts normal distribution. e corresponding relationships of specific samples are shown in Table 4. is paper adopts the modified half-wave sinusoidal load model. According to the collision tonnage and collision velocity under various working conditions, 10 sample timehistory curves of the vessel collision force are obtained. e sample time-history curves of vessel collision force under partial working conditions are shown in Figure 8.

Calculation of Bridge Failure Probability.
Bridges are prone to bending failure at the top of pile foundation under vessel collision. When bending failure occurs at the top of the most unfavorable pile foundation, it is considered that the entire bridge system is damaged.
As shown in Figure 4, the failure probability of bridge pier is calculated by ultimate rotation angle of pile top and the most unfavorable rotation angle of pile top. In this instance, the ultimate rotation angle is calculated by the actual reinforcement at the pile top according to the M-Φ curve, and the most unfavorable rotation angle at the pile top is obtained by analyzing the finite element model. 10 different sample time-history curves of vessel collision force are applied to the bridge cap, elastic-plastic hinges are set on the top of the pile, and through nonlinear dynamic time-history analysis, the maximum rotation response of the plastic hinge zone at the top of the most unfavorable pile foundation of the main pier group pile foundation is shown in Table 5.
According to the actual reinforcement of the pile top, the ultimate rotation angle of the pile top of No. 1 pier is 2.432E-02, and the ultimate rotation angle of the pile top of No. 2 pier is 2.475E-02. According to the above theory, the failure probability of middle-pier calculated by the program is shown in Table 6.

Conclusions
In this paper, a method and process of reliability analysis of bridge structure under vessel collision are established. In view of the implicit and nonlinear characteristics of the limit state performance function of the vessel-bridge collision and the limitations of the traditional reliability analysis method to solve the problem, a statistical moment method based on finite element analysis is proposed, which greatly improves the computational efficiency of structural failure probability with implicit performance function. e method is applied to the reliability analysis of vessel-bridge collision. e reliability analysis model and performance function of vesselbridge collision failure probability calculation are established, and the probability model of each influencing factor affecting the failure probability is established. Due to the characteristic of large discreteness of vessel tonnage in bridge areas, the unknown distribution model is used to simulate the distribution of vessel collision tonnage, and the first four order statistical moments of statistical variables are used to describe its probability distribution. e method established in this paper is used to analyze the reliability of vessel-bridge collision in practical bridge engineering, which provides a reference for the specific application of this method in practical engineering.

Data Availability
e underlying data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.   Advances in Materials Science and Engineering 9