Kriging Model for Reliability Analysis of the Offshore Steel Trestle Subjected to Wave and Current Loads

: Offshore steel trestles (OSTs) are exposed to severe marine environments with stochastic wave and current loads, making structural safety assessment challenging and difﬁcult. Reliability analysis is a suitable way to consider both wave and current loading intensity uncertainties, but the implicit and complex limit state functions of the reliability analysis usually imply huge computational costs. This paper proposes an efﬁcient reliability analysis framework for OST using the kriging model of optimal linear unbiased estimation. The surrogate model is built with stochastic waves, current parameters, and the corresponding load factors. The framework is then used to evaluate the reliability of an example OST subjected to wave and current loads at three limit states of OST, including ﬁrst yield (FY), full plastic (FP), and collapse initiation (CI). Three different distributions are used for comparison of the results of failure probability and reliability index. The results and the computational cost by the proposed framework are compared with that from the Monte Carlo sampling (MCS) and Latin hypercube sampling (LHS) method. The inﬂuences of sample number on the prediction accuracy and reliability index are investigated. The inﬂuence of marine growth on the reliability analysis of the OST is discussed using MCS and the kriging model. The results show that the reliability analysis based on the kriging model can obtain the reliability index for the OST efﬁciently with less calculation time but similar results compared with MCS and LHS. With the increase of the number of samples, the prediction accuracy of the kriging model increases, and the corresponding failure probability ﬂuctuates greatly at ﬁrst and then tends to be stable. The reliability of the example OST is reduced with the increase of marine growth, regardless of the limit state.


Introduction
In recent decades, many long-span sea-crossing bridges have been constructed or are in the planning phase to accelerate worldwide ocean development and promote coastal economics [1][2][3]. As a type of auxiliary facility for construction that can be quickly assembled on the sea, the offshore steel trestles (OSTs) have been widely used to construct sea-crossing bridges. These OSTs play an essential role in transporting materials, equipment, and labor in the construction of a sea-crossing bridge.
An OST is usually composed of bridge deck, bailey beam, and jacket support structures and located in a marine environment consisting of wave and current [4]. Many studies have been conducted to investigate the wave, current, or both of them, applied on different components of an offshore bridge, including the deck, pile-cap, pier and the piles [5][6][7][8][9][10]. They state that calculating both the applied loads and the reaction forces for coastal decks and bridge substructures are complex due to the random nature and challenging characteristics of the hydrodynamic loads [11][12][13]. For example, there are significant aleatory uncertainties in the wave slamming and breaking process [14][15][16], influences of trapped air [17,18], and variable temporal and spatial evolution of hydrodynamic pressures on different structures [19]. Moreover, when the relative angle between the wave propagation and structural orientation is not normal, the applied load and structural reaction are sensitive to the three-dimensional (3D) influence, either due to the obliqueness of the wave or the skewness of the structure [20].
Pushover analysis is a common nonlinear approach in evaluating the ultimate behavior and capacity of offshore jackets under environmental loads [21]. Raheem and Shehata [22] made a nonlinear response analysis of a fixed jacket type platform using fifth-order Stokes wave theory and accounting for direction spreading and irregularity of the wave profile. Naess et al. [23] predicted both short-term and long-term extreme response statistics of a jacket structure based on the time histories of deck response under waves. Wei et al. [24] developed a nonlinear analysis framework based on the pushover approach and applied it to assess the nonlinear response of OSTs subjected to wave and current loads. All these studies emphasize that the combination of high wave and strong current aggravate the nonlinear response of offshore jacket structures significantly, and hence threaten the safety of OSTs [25]. Therefore, the accurate calculation of wave, current loads, and the precise determination of structural performance are necessary to utilize OSTs successfully during their overall life cycle.
As the nonlinear analysis of OST is usually time-consuming and requires professional knowledge in structure engineering, many design practices recommend deterministic approaches in the design of OSTs, for example, by multiplying the safety factor, to ensure the safety of the OSTs under wave and current [26]. However, the reliability levels achieved by the partial safety factors are not optimal and clear. As an approach considers the randomness of resistance and the load effect, reliability analysis can control the occurrence of failure at an acceptable level according to the structure's reliability requirements to seek the balance between safety and economic benefit of the structure. To avoid the computational difficulty and the lack of accuracy caused by the strong nonlinear equation, scholars use the Monte Carlo sampling (MCS) method to solve offshore structures' reliability [27]. However, MCS usually generates duplicate samples in the random sampling. To obtain more accurate reliability results, one needs to generate a large amount of calculation (the total number of samples N ≥ 100/P f ) by MCS for a better assessment of structural reliability [28]. Moreover, every sample needs to be calculated by pushover analyses, which consumes too much time and limits its practical application scope. To improve the sampling efficiency and reduce the amount of calculation, researchers developed some advanced sampling methods; for example, the Latin hypercube sampling (LHS) method, which aims to use as few samples as possible to construct the structure's function in order to predict the rest of the unknown information [29]. The results show that LHS can reduce the sampling number significantly under the same calculation precision as MCS. Besides improving sampling methods, one other solution is to increase the calculation efficiency of the reliability analysis. Kriging models have been widely used in offshore structural analyses. The advantage of the kriging model is that it is an efficient method that approximates the high-order nonlinear function by the surrogate model, and avoids directly solving the high-order nonlinear function. The establishment of a surrogate model only needs a small number of samples, and can greatly reduce the reliability calculation time. Yang et al. [30] constructed a finite element (FE) model to calculate the dynamic response of a 5 MW wind turbine, then used the kriging model to replace the original time-consuming dynamic response analysis of the FE model and gain the final best design point, considering uncertainties by a global optimization algorithm. Morató et al. [31] proposed a computationally efficient reliability framework for an offshore wind turbine jacket by the kriging model to approximate the system's response, capturing both the dynamic behavior of the structural and inherent uncertainties. Benefiting from the advantages mentioned above, the kriging method has been widely used in the reliability analysis of offshore wind turbines. However, the existing research has not researched the reliability of OSTs using the kriging model. This paper proposed a reliability analysis method for OST under wave and current loads. Kriging model of optimal linear unbiased estimation and LHS with good spatial sampling are employed in the method. First, the reliability analysis method framework of OST based on the kriging model is elaborated in detail. Second, the engineering background, finite element modeling, and wave and current load calculation methods of the calculation case are introduced. Finally, based on the reliability analysis method proposed in this paper, the reliability index of the calculation case is solved. Three different distributions are used for comparison of the results of failure probability and reliability index. The results of the reliability calculation and analysis time of this method are compared with those of MCS and LHS, proving the method's efficiency and accuracy. The influence of the sample number on the kriging model of OST is discussed in detail. In addition, the influence of marine growth on the reliability analysis of the OST is discussed using MCS and the kriging model.

The Framework of the Reliability Analysis Method for Offshore Steel Trestles
A reliability analysis method for OST under wave and current loads is proposed in this section. Both kriging model and LHS sampling are employed to reduce the amount of calculation of the proposed procedure. The framework includes the wave and current load model, nonlinear assessment of OST subjected to wave and current, kriging model and reliability analysis process.

Calculation of Wave and Current Loads
Considering the shallow water depth of the OST and the apparent nonlinear characteristics of waves, the stream function wave theory [32] is used to simulate the waves in this paper. Stream function wave theory has the advantages of fewer undetermined coefficients and is easy to extend to any order. It is very suitable for nonlinear wave simulation in shallow water where the OST is located. At the same time, according to the relevant evaluation standards, the wave theory applicable to the working conditions (wave height H 5% and water depth d) is identified in this paper. The 10th order stream function theory is selected to generate nonlinear regular wave kinematics in this paper. The study shows that the conditions considered in this paper meet the stream function wave theory's applicable scope, so it is reasonable to use the stream function wave theory to simulate the OST waves. The reason for judging the wave theory can be observed in Wei et al. [24]. Since none of the used wave heights exceeds the breaking limit of 0.78d [33], the slamming force of breaking wave is not included.
For the beams and substructure of the OST, such as a transverse brace and steel pipe pile, the ratio of component diameter D to wavelength L meets D/L < 0.2, which belongs to a small-scale structure. The wave force is calculated by the Morison equation. The basic idea of the Morison equation is to divide the wave force f H into two parts: one is the inertial force term f D, which is proportional to the acceleration, and the other is the resistance term f I, which is proportional to the square of velocity, that is, acting on any height of the vertical column. The horizontal wave force on the unit column height at z is f H : where ρ is the density of seawater; D is the diameter of the cylinder; u x is the horizontal velocity of the water particle; and the drag coefficient and inertia force coefficient of C D and C M are 1.2 and 2.0, respectively. Ocean currents usually refer to the large-scale horizontal or vertical flow of seawater from one area to another. The causes can be divided into three types: tidal currents, windinduced currents, and wave currents. For the OST studied in this paper, the current of its marine environment mainly considers the tidal current, while the influence of wave currents and wind-induced currents can be ignored. The following equation can calculate the current profile along with the water depth: where u c,sub is the subsurface current velocity at the still water level; d is the water depth; z u is the distance from the still water level. The above Morison equation represents only the wave action result, without considering the combined action of waves and current. When there is current flow in the wave field, the combined action of waves and current should be considered. The current is often regarded as a constant current, and the force on the cylinder is only the drag force in the Morison equation. The combined action of the current velocity and the water particle velocity inevitably causes a change in the cylinder's drag force. As the drag force is proportional to the square of the velocity, the effect is significant. The combined action of waves and current is very complicated. In the current engineering design, the horizontal wave current force acting on the unit height of the vertical cylinder at any height z is often calculated as f : where v is the current velocity. The highest wave height set in this article does not act on the deck in the hypothesis of this article. Therefore, the load acting on the deck is not listed here. However, if the waves reach the deck of the OST, the uplift and overturning moment cannot be neglected [34].

Nonlinear Assessment of Offshore Steel Trestles
This paper evaluates the capacity of the OST under the combined action of waves and currents. Wei et al. have defined three limit states, named first yield (FY), full plastic (FP), and collapse initiation (CI), according to the performance characteristics of the OST under stress [24]. The corresponding damage conditions are given in Figure 1. Wei et al. have investigated the dynamic effect of an OST subjected to wave and current loads in the term of dynamic amplification factor (DAF) and found the DAF of the large wave is always smaller than 1.05 [24]. Neglection dynamic effect results in tiny influence on the nonlinear response of the OST. Therefore, the capacity of the three limit states of the OST, first yield (FY), full plastic (FP), and collapse initiation (CI), are calculated by the pushover analyses. The process of the pushover analyses is described as follows:

1.
Generating a regular wave train with the given wave height H 5% and the corresponding wave period T 5% , which is calculated as a function of wave height according to Valamanesh et al. [35], where g defines the acceleration of gravity: 2.
Using Equation (3) to calculate wave and current loads in the process of a complete phase of a regular wave (through-crest-through) stepping through the OST with a notably small-time interval. Define the wave and current loads corresponding to the maximum base shear, abbreviated as F demand, in the process of wave stepping through the OST as a lateral load pattern.

3.
Pushover analyses of the OST are conducted based on the lateral load pattern. The pushover analyses will not stop until the OST collapses, taking the lateral displacement at the top of the OST as the classification index for three limit states. The base shear of the OST in three limit states corresponds to the capacity of the OST in three limit states, abbreviated as F capacity. See Wei et el. [24,36] for more details of the process of pushover analyses.
This paper introduces the dimensionless load factor into the performance evaluation method to measure the capacity of the OST, which is defined as follows: The load factors corresponding to the first yield (FY), full plastic (FP), and collapse initiation (CI) are γ C1 , γ C2 , and γ C3 , respectively. Figure 2 shows the relationship curve between the load factors and the lateral displacement at the top of the OST. When γ C = 1, the structure's capacity to resist damage F capacity (FY, FP, and CI) is just equal to the maximum base shear F demand generated by a wave and current loads, and the structure is in the critical state of corresponding damage. When γ C > 1, the capacity of the structure F capacity is greater than the maximum base shear F demand caused by wave and current loads, and the structure has no corresponding damage. When γ C < 1, the corresponding damage occurs before the structure fully bears the given wave and current loads. Therefore, this paper uses the load factor to replace the response value corresponding to the OST limit state, which is more conducive to establishing the OST limit state function. The specific analysis process of this method is shown in Figure 3. Because the evaluation method of the capacity of the OST under wave and current loads can realize the determination of the initial state of first yield (FY), full plastic (FP), and collapse initiation (CI), three limit state functions are still established based on these three limit states. Additionally, since the load factors corresponding to the three limit states can be obtained by this method, the three limit state functions Z can be rewritten as follows: First yield: Full plastic: Collapse initiation: where R is the structural resistance; S is the load effect corresponding to wave height H 5% and current velocity v.
According to the basic form of a limit state function, we only need to build the kriging model, which characterizes the relationship between wave height, current velocity, and load factor. By judging whether the predicted value of the load factor corresponding to the sample points randomly selected by LHS is less than 1, the failure point in the sample points can be obtained, and the failure probability and reliability index of the structure can be obtained.

Kriging Model for Offshore Steel Trestles Subjected to Wave and Current Loads
The kriging model was initially derived and used to calculate deposit reserves and error estimates. The kriging method is a statistical prediction method based on the stochastic process. The corresponding relationship between the response value and the kriging model's independent variable consists of the regression part F(β, x) and the random process z(x). The specific model is as follows: where F(β, x) is used to provide the global approximation of the kriging model, the mathematical expectation of y(x); the mean value of z(x) is 0, and the variance is σ z 2 . A random process with nonzero covariance provides an approximation to simulate local deviations, that is, local changes of y(x). In engineering practice, the regression part generally adopts a polynomial form.
where β is the regression coefficient of each function; p is the number of f i (x); and f i (x) is a polynomial function of variable x, which is similar to the polynomial form of the response surface method. For structures with a high nonlinear degree, such as OST, using a quadratic function as a regression model can obtain calculation results with high accuracy. The covariance matrix z(x) is: where x i and x j are the i-th and j-th components of the training sample, respectively, and R θ, x i , x j is the correlation function with a super parameter θ, representing the spatial correlation between training sample points, which plays a decisive role in the simulation's accuracy. There are many correlation functions, such as the exponential correlation function, Gaussian correlation function, Matérn-3/2, and Matérn-5/2. The most used function is the Gaussian correlation function, which is as follows: where n dv is the number of known design variables and θ k is the k-th element of vector θ. Given the training sample S = [x 1 , x 2 , · · ·, x n ] and its real response value Y = [y 1 , y 2 , · · ·, y n ] (n is the number of samples), the correlation matrix R of n × n and the regression partial matrix F of n × p are as follows: The estimated valueŷ(x) of any point x can be expressed as: where r(x) is the correlation vector between test point x and the sample point, which can be expressed as: The maximum likelihood estimation factor is as follows: Under the Gaussian process assumption, the optimal kriging model is constructed by solving the unknown quantity θ. According to the maximum likelihood estimation, it can be concluded that: Among them, θ can be obtained by solving the following optimization problems: A certain amount of available information is required to develop the kriging model. Therefore, it is necessary to select a certain number of sample points in a specific design space range and conduct numerical experiments. LHS is a multi-dimensional stratified sampling method. Compared with MCS, LHS is designed to accurately reconstruct the input distribution by sampling with fewer iterations, providing more comprehensive sample information to establish the kriging model. The basic principles of LHS are as follows: 1.
Define the sampling number N participating in the computer operation; 2.
Divide each input equally into N columns, where Select only one sample for each column, and each column's bin position is random.
LHS avoids the collapse of sample data and can use fewer samples to represent the overall distribution of parameters, which has a good filling for the sampling space, as shown in Figure 4. Suppose a parameter A obeys a Gaussian normal distribution, the average value of the ideal distribution is 1.5, and the standard deviation is 2.0. LHS and MCS were used to obtain samples with a capacity of 1800, and the mean and variance of samples were calculated. A total of 20 simulations were performed. Figure 5a,b show the fluctuation of the sample mean and variance with simulation times, respectively. The results show that LHS's statistical parameters are more accurate and stable than those obtained by MCS.

Reliability Analysis Process
The analysis can proceed as follows: 1.
According to the practical problems of OSTs, the influencing factors of the capacity of OSTs are determined; namely, random variable X 1 , X 2 , · · · , X j and the limit state of OST are defined; 2.
According to the distribution function of the random variable X 1 , X 2 , · · · , X j , LHS is used to sample random variables and m×j dimension random variable samples are generated; 3.
The random variable samples are substituted into the finite element model to generate the corresponding analysis conditions and the corresponding response values of the OST's limit state are obtained; 4.
The kriging model of the response value in accordance with a certain accuracy is constructed by using the random variable samples and the response values; 5.
Again, LHS is used to extract n×j dimensional random variable samples, and the predicted response values of each random variable sample are obtained by substituting them into the kriging model; 6.
The predicted response value is substituted into the limit state equation Z. If Z is less than 0, it is the failure point. Assuming that the number of failure points is k, the failure probability is P f = k/n. The failure probability can be expressed as: When Z < 0, I(Z) = 1; when Z < 0, I(Z) = 0. Then, the reliability index of OST is: where Φ is the standard normal distribution function.
As an approximate fitting model, the kriging model can have some errors with the real model. To test the fit effect of the kriging model, it is mainly judged by the Theil inequality coefficient, which can be calculated as follows: where n is the number of samples; y i is the calculated actual value of the i-th sample point; y i is the predicted value of the kriging model of the i-th sample point; and τ is the Theil inequality coefficient, which is between 0 and 1, and the closer to 0, the better the fitting effect is. This paper's method only pays attention to the randomness of load parameters, that is, the randomness of wave and current parameters (X wave , X current ) in the reliability analysis of OST. The kriging models proposed above are compiled using ooDACE toolbox of the software MATLAB (Mathworks Inc., Natick, MA, USA) [37].

Example Structure
The OST of the Pingtan Strait Rail-cum-road Bridge is taken as the example structure, as shown in Figure 6. The OST is a two-span continuous beam structure with a length of 54 + 54 m, a width of 6 m, and a water depth of 25 m. The height from the bridge deck to the seabed is 42 m, and the clearance of the OST is 15 m. The bridge deck is fixed by bailey beams, which are supported by beams. The steel pipe pile and the supporting frame are welded together to form a jacket-type support system. The steel pipe pile is embedded in the seabed rock. The bailey beam's height is 1.5 m, and the distance between adjacent bailey plates is 1.5 m. The height of the steel pipe pile is 40 m, the diameter is 1.5 m, and the spacing of the steel pipe pile is 12 m. The supporting frame's height is 6.5 m, the diameter of the upper and lower transverse bracing and diagonal bracing is 0.63 m, and the diameter of the vertical brace is 0.35 m. The thicknesses of the steel pipe pile and supporting frame are 16 mm and 10 mm, respectively. The bailey beam is made of Q345 steel, the steel pipe pile is made of Q235 steel, and the bridge deck is made of 16 Mn steel. The specific dimensions are shown in Figure 6. The measured wave and current parameters in the bridge site's sea area are shown in Table 1.

Establishment of the Finite Element Model
The three-dimensional numerical model of the OST is established using the finite element software USFOS (SINTEF marintek, Oslo, Norway and the Norwegian University of Science and Technology, Trondheim, Norway) as shown in Figure 7. Nonlinear beam elements are used to simulate the bailey beam, jacket, and steel pipe pile. The plastic hinge model is used to simulate the nonlinearity of the bar material. The elastic-plastic model considering strain strengthening is adopted for the mechanical behavior of the plastic hinge. The plastic hinge is automatically inserted into the two ends or the beam element's middle point according to the element's plastic development. The four-point shell element model is used to establish the bridge deck. A rigid connection is adopted between the bailey frame and steel pipe column, and a rigid connection is still adopted between the superstructure and substructure. The mass of guardrail and vehicle load on the bridge deck is considered by increasing the bridge deck's thickness. The bottom of the steel pipe pile is consolidated and constrained. The nonlinear geometric effects of large displacement and strain are considered. In this paper, the waves and the currents are set in the same direction and along the transverse direction of the OST. The local elastic-plastic deformation of tubular joints occurs when subjected to large external forces. However, the traditional rigid joint connection method assumes that the local elastic deformation and plastic deformation after the yielding of tubular joints are not considered in the stress analyses. This paper uses MSL joint model, developed by MSL Engineering Ltd. in County Cork, Ireland [38] to simulate the plastic behavior in the OST joints. The simulation method models the pipe joint's nonlinear deformation after being stressed by a nonlinear spring between the chord skin and the chord centerline.

Stochastic Model of Wave and Current Conditions
Since there are only three sets of wave and current conditions, three distributions are used for the comparison of the probability distribution of wave and current parameters of the OST. They are the three-parameter Weibull extreme value distribution, the Gaussian distribution, and the Rayleigh distribution.
The three-parameter Weibull extreme value distribution model was proposed by the South China Sea Institute of Oceanography, Chinese Academy of Sciences. The threeparameter Weibull extreme value distribution function is as follows: where α is the scale parameter, µ is the position parameter, and ξ is the shape parameter.
The Gaussian distribution function is as follows: where µ is the mu parameter, and σ is the sigma parameter.
The Rayleigh distribution function is as follows: where σ is the sigma parameter. The extreme value x p of environmental load parameters with the n-year return period can be determined according to the probability distribution function of load parameters. The probability of occurrence of environmental load parameters with the n-year return period is 1/n.
Then, the extremum quantile x p can be expressed as follows: Referring to the measured wave and current parameters in the bridge site's sea area (Table 1), the distribution of wave height and current velocity is shown in Figure 8. The three-parameter Weibull extreme value distribution function, Gaussian distribution function and Rayleigh distribution function of wave height and the current velocity are as follows:

Verification of the Proposed Method
When constructing the kriging model, it is assumed that the wave height and current velocity parameters are uniformly distributed. The sampling range of wave height and the current velocity is based on the actual distribution range of wave height and current velocity; that is, the sampling range of wave height is (2 m, 12 m), and the sampling range of current velocity is (1 m/s, 4 m/s). To verify the proposed method, this paper selects 500 samples of wave height and current velocity according to the uniform distribution type and distribution range. The predicted values of 500 samples are obtained by using the kriging model. The Theil inequality coefficient of the kriging model is calculated from the actual values and the predicted values obtained by the kriging model. When the Theil inequality coefficient is less than 0.01, the kriging model's prediction accuracy can be considered to meet the requirements.
The kriging models corresponding to three limit states with specific sample numbers of 100 are given in Figures 9-11, as well as the comparison between actual and predicted load factor values. The figure shows that the kriging model is consistent with the original surface and the predicted value is equal to the actual value. The Theil inequality coefficient of the kriging model is less than 0.01, which verifies the accuracy of the method.    Table 2 gives the failure probability and reliability index of the three limit states of the OST under three different distributions of wave and current. Table 2 shows that the failure probability of the OST gradually decreases with the increase of the limit state and its corresponding reliability index increases with the increase of the limit state. Different distributions of wave and current result in different failure probability and reliability index of the three limit states of the OST. The result of using the Rayleigh distribution will have the highest failure probability, and the use of Gaussian distribution will have the lowest failure probability. This is due to the use of Rayleigh distribution for fitting of wave and current distribution, resulting in more large waves and rapid current velocities. The failure probability obtained for limit state FY is relatively high; because the partial yielding in the definition of FY is more likely to occur, it does not affect the normal use of the OST. As the distributions of both wave and current conditions affect the reliability significantly, it is necessary to have more measured data to fit the distributions and make the results of reliability more realistic. Table 2. Failure probability and reliability index corresponding to three limit states of OST using the three-parameter Weibull extreme value, Gaussian, and Rayleigh distribution for fitting of wave and current parameters.

Comparison of MCS, LHS, and Kriging Model
The failure probability and reliability index of OST are calculated using the other two widely used reliability analysis methods, MCS and LHS, to reflect the accuracy and efficiency of the reliability analysis method for OST based on the kriging model. Threeparameter Weibull extreme value distribution is used for sampling. Three methods are defined as follows: 1.
MCS is used to randomly select a large number of samples of wave height and current velocity that obey the three-parameter Weibull extreme value distribution, all samples are analyzed by the finite element method one by one, and the load factors corresponding to the three limit states of the OST are obtained. The number of failure samples whose load factor is less than one is counted for each limit state. Finally, the ratio of failure events to the whole sample is calculated. The failure probability and reliability index corresponding to the three limit states of the OST are obtained.

2.
The second method is LHS. The only difference between LHS and MCS is that the sampling method is different. LHS is used in the sampling of wave height and current velocity samples.

3.
For the kriging model, it is worth noting that the kriging model constructed in this paper uses LHS twice with different purposes. The first time LHS is used to extract a certain sample number to establish a kriging model, the second time LHS is used to extract the load factors to calculate reliability index and failure probability. Table 3 summarizes the reliability results obtained by the three methods and compares the total analysis time spent by the three methods in solving the reliability results. The whole analysis time is mainly the calculation time of all finite element analyses. The time spent in other processes is generally less than ten minutes, which can be ignored, such as sample sampling and failure sample quantity statistics.  Table 3 shows that when MCS is used to calculate the reliability of OST directly, 1.15×10 5 finite element analyses are needed to obtain an accurate approximate solution of the structural reliability. The analysis time is approximately 1437.5 h. When the LHS is used instead of the MCS, the calculation amount of finite element analyses can be effectively reduced, the time saved is approximately 26%, and the relative error of the reliability index is within 1%. The LHS is better than the MCS in the reliability calculation. The reliability analysis method of OST based on the kriging model proposed in this paper only needs 100 finite element analyses to obtain more accurate reliability results. The analysis time is approximately 1.25 h. Compared with MCS and LHS, the calculation time is reduced by three orders of magnitude, and the relative error of the reliability index is within 0.01%.

Influence of Sample Number on the Prediction Accuracy of the Kriging Model
In theory, the more samples that are extracted when constructing the kriging model, the higher the accuracy of the kriging model. However, if the number of samples is too large, the number of finite element analyses inevitably increases, which affects the calculation efficiency and fails to reflect the efficiency of the dynamic kriging model. To improve the calculation efficiency, it is necessary to determine the appropriate sample size using LHS. This section analyzes the impact of sample size on the accuracy of the kriging model.
Previous studies have shown that the number of samples required for constructing the kriging model should not be less than 30. Therefore, the number of samples is selected to be 40, 60, 80, 100, and 120. Table 4 shows the kriging model's Theil inequality coefficient corresponding to three limit states under different sample numbers. Table 4 clearly shows that with the increase in sample number, the Theil inequality coefficient of the kriging model corresponding to the three limit states gradually decreases; that is, the prediction accuracy increases continuously. When the number of samples is 100, the kriging model's Theil inequality coefficients corresponding to three limit states are less than 0.01, which meets the prediction accuracy. Therefore, a 100 sample number is used to construct the kriging model. After obtaining the kriging model of the relationship between wave height, current velocity, and load factor, to obtain the failure probability and reliability index of OST, we need to extract a certain number of load factors using LHS again. Although LHS can represent the random distribution characteristics of parameters with only a few samples, it may not obtain an accurate reliability index and failure probability with a few samples. Therefore, it is necessary to explore the influence of the sample number on the failure probability of the OST to determine the appropriate sample number in advance. The sample size ranges from 5000 to 200,000, and the step size is 5000. The sample of wave height and current velocity obey the three-parameter Weibull extreme value distribution. Figure 12 shows the variation in failure probability with the number of samples corresponding to the three limit states of the OST. Figure 12 shows that when the number of samples is small, the failure probability corresponding to each limit state fluctuates greatly. It is mainly because of the large wave height and current velocity required for the structure's FY, FP, and CI. Suppose the number of samples is too small. In that case, the number of Latin hypercube delamination intervals is less. The wave height or current velocity range corresponding to the equal cumulative probability interval is large, resulting in the randomness of large wave height or high current velocity being sampled. In addition, when the number of samples reaches 35,000, 55,000, and 85,000, the failure probability corresponding to the FY, FP, and CI of OST fluctuates slightly, and the fluctuation ranges are ±0.5%, ±1%, ±1.5%, respectively. As demonstrated by the blue dotted lines in the figure, the failure probability tends to be stable. It should be noted that the size of the failure probability is inversely proportional to the number of samples corresponding to the stable failure probability. That is, the greater the failure probability is, the fewer samples need to be calculated. To ensure that the failure probability corresponding to the three limit states can reach high accuracy and reduce the sampling process, the sample size is 85,000.

Influence of Marine Growth on the Reliability Analysis of the OST
A measurement of marine growth on a platform in the China South Sea demonstrates that the annual marine growth thickness can be as large as 10 mm, when the water depth is smaller than 50 m [39]. Studies have shown that marine growth increase wave and current loads on the OST by increasing the hydrodynamic radius and the roughness of the steel pipe pile and changing the hydrodynamic coefficient, and thereby affecting the capacity of the OST [40]. To understand the influence of marine growth on the reliability of the OST, the proposed kriging model is used to analyze the reliability of the OST under different thicknesses of marine growth.
The load factors γ C1,C2,C3 (t m , H 5% , v) of three limit states corresponding different thickness of marine growths have an approximately linear relationship with the load factors γ corresponding to three limit states of the OST. This relationship does not include the case where the thickness of the marine growth t m = 0 m, because the hydrodynamic coefficient of t m = 0 m is set up as 1.2 for C D and 2.0 for C M, while the hydrodynamic coefficient is set up as 1.4 for C D and 1.8 for C M when t m is not 0 m [41]. Therefore, based on the condition of t m = 0.02 m, the load factors γ C1 , γ C2 , and γ C3 with any thickness of marine growth of the OST under wave and current loads can be calculated according to the following equation: where λ γ C1 ,γ C2 ,γ C3 (t m ) represents the relationship between γ C1,C2,C3 (t m , H 5% , v) and γ C1,C2,C3 (0.02, H 5% , v), the correlation equation is as follows [42]: λ γ C1 ,γ C2 ,γ C3 = −1.45t m + 1.027 (35) As the above relationship does not include the case of t m = 0 m, when calculating the failure probability and reliability index corresponding to the thickness of different marine growths, it is necessary to construct the kriging model corresponding to the three limit states of the OST when the thickness of the marine growth is at first 0.02m. Based on the kriging model of t m = 0.02 m combined with Equation (35), the kriging model corresponding to the three limit states of the OST under different thicknesses (0.04 m, 0.06 m, 0.08 m, and 0.1 m) of marine growths is constructed, and then the corresponding reliability index is obtained using the LHS. Figure 13 shows the reliability index and failure probability of the OST under different marine growth thicknesses obtained by the kriging method and MCS. It can be observed from Figure 13 that the reliability indexes of the OST under different marine growth thicknesses, calculated by the kriging method and MCS, are in good agreement with the results, and the relative error is within 1.5%. In the process of the OST from having no marine growth to having marine growths, its failure probability and reliability indexes are significantly increased and decreased, respectively. As the thickness of marine growths increases, the failure probability of the OST increases linearly, and the corresponding reliability indicators gradually decrease linearly. Table 5 shows the change of the reliability index of the OST with the increase of the thickness of the marine growth from 0 to 0.1m. It can be observed from Table 5 that with the influence of marine growths, the failure probability corresponding to three limit states of the OST will be significantly increased, and the corresponding reliability index will be significantly reduced. Among them, the degree of reliability index reduction is as follows: First yield (31.23%) > Full plastic (17.79%) > Collapse initial (10.15%). Marine growths will affect the fatigue behavior of the OST, and thus affect the capacity of the OST. However, this aspect is not considered in this paper and is worth considering in future work.

Conclusions
Based on the kriging model and LHS, a reliability analysis method for OST is proposed in this paper. The reliability index of the OST used in the Pingtan Strait Rail-cum-road Bridge is solved based on this method. Three different distributions are used for comparison of the results of failure probability and reliability index. The reliability calculation results and analysis time of this method are compared with MCS and LHS. The sample number's influence on the prediction accuracy and reliability index of the kriging model for the reliability analysis of OST is discussed. In addition, the influence of marine growth on the reliability analysis of the OST is discussed using the kriging model and MCS. The main conclusions are described as follows: 1.
The kriging model is constructed assuming that the wave height and current velocity parameters are uniformly distributed. The sampling range of wave height is (2 m, 12 m), and the sampling range of current velocity is (1 m/s, 4 m/s). A total of 500 samples of wave height and current velocity are selected according to the uniform distribution type and distribution range to verify the proposed method. The Theil inequality coefficient of the kriging model is less than 0.01, which verifies the accuracy of the method. The failure probability of the OST gradually decreases with the increase of the limit state, all less than 0.600%, and its corresponding reliability index gradually increases with the increase of the limit state, all greater than 2.5. Different distributions of wave and current result in different failure probability and reliability index of the three limit states of the OST. The result of using the Rayleigh distribution will have the highest failure probability, and the use of Gaussian distribution will have the lowest failure probability.

2.
Compared with MCS and LHS, the reliability analysis method based on the kriging model can obtain the reliability index of OST efficiently and accurately. The analysis time is approximately 1437.5 h and 1062.5 h when using MCS and LHS, while the calculation time is approximately 1.25 h when using the kriging model. Compared with MCS and LHS, the calculation time is reduced by three orders of magnitude. Compared with MCS, the relative error of the reliability index using the kriging model is within 0.01%, which shows the accuracy of the kriging model.

3.
When constructing the kriging model, the sample number is selected to be 40, 60, 80, 100, and 120. With the increase in sample number, the Theil inequality coefficient of the kriging model corresponding to the three limit states gradually decreases; that is, the prediction accuracy increases continuously. When the number of samples is 100, the kriging model's Theil inequality coefficients corresponding to three limit states are less than 0.01, which meets the prediction accuracy. When extracting a certain number of load factors using LHS to obtain an accurate reliability index and failure probability, the sample size ranges from 5000 to 200,000 with the step size of 5000 are discussed. When the number of samples reaches 35,000, 55,000, and 85,000, the failure probability corresponding to the FY, FP, and CI of OST fluctuates slightly, and the fluctuation ranges are ±0.5%, ±1%, and ±1.5%, respectively. To ensure that the failure probability corresponding to the three limit states can reach high accuracy and reduce the sampling process, the sample size is 85,000.

4.
The influence of marine growths on the reliability of the OST is discussed using MCS and the kriging model. The reliability indexes of the OST under different marine growth thicknesses calculated by the kriging method and MCS are in good agreement with the results, and the relative error is within 1.5%. With the influence of marine growth, the failure probability corresponding to three limit states of the OST will be significantly increased, and the corresponding reliability index will be significantly reduced. Among them, the degree of reliability index reduction is as follows: First yield (31.23) %) > Full plastic (17.79%) > Collapse initial (10.15%).
It should be noted that the above conclusions are drawn based on the reliability analyses of OST, which only considers the randomness of waves and currents. The dynamic effect of wave load on the structural response is neglected as well [24]. The influences of more factors, such as geometric and material uncertainties [43], wave direction [44] and ship/debris impact [45] on the reliability deserve to be investigated in further study. The correlations between wave, current, and marine growth are not considered due to the lack of measured data in this research. The correlations among these parameters are worth being included in the further study [46,47]. In terms of the reliability analysis approach, some methods in the framework of active learning [48] or Bayesian experimental design [49] may achieve better performance than kriging and LHS method. The active learning method selects a small amount of experimental design points to build the initial surrogate model and identify the best next point by active learning equation until the required failure probability is obtained. Some possible active learning functions include but are not limited to, H learning function [50], least improvement function [51], reliability based expected improvement function [52], and folded normal based expected improvement function [53]. Bayesian experimental design selecting the next best point considers not only the potentially "dangerous" points to be close to the limit state, but also the updating process of the next best point [54][55][56]. These advanced approaches are worth considering for future reliability analysis frameworks of offshore structures.