3D Numerical Modeling and Quantification of Oblique Wave Forces on Coastal Bridge Superstructures

<jats:p>Simply supported bridges comprise the majority of bridge systems in coastal communities and are susceptible to severe damage from extreme waves induced by storms or tsunamis. However, the effects of oblique wave impacts have been less investigated due to the lack of appropriate numerical models. To address this issue, this study investigates the effects of wave incident angles on coastal bridge superstructures by developing an advanced computational fluid dynamics (CFD) model. Different wave scenarios, including wave height, relative clearance, incident angle, and wavelength are tested. It is found that the maximum wave forces in the horizontal and longitudinal directions could reach 1901 and 862 kN under extreme conditions, respectively, destroying bearing connections. Three surrogate models, i.e., the Gaussian Kriging surrogate model, the Artificial Neural Network (ANN), and the Polynomial Chaos Expansion (PCE), are established by correlating the wave parameters with the maximum wave forces. Through comparisons among the three surrogate models, it is found that the 3-order PCE model has better performance in predicting loads in vertical and horizontal directions, while the ANN model is more suitable for results in the longitudinal direction. This study contributes to the optimized design of coastal bridges and also offers an opportunity for future studies to investigate hazard damage-mitigation measures.</jats:p>

section_content: The prosperous economy in the bay area triggers the construction of more coastal bridges, forming comprehensive transportation networks. The harsh environment, particularly extreme waves and surges induced by tropical cyclones or tsunamis, poses considerable challenges to the structural safety of these infrastructures. According to a post-hazard report, Hurricanes Ivan (2004) and Katrina (2005) led to a USD 2. 5 billion loss for the damaged bridge networks [1]. More recently, the 2011 Great East Japan Tsunami destroyed more than 250 coastal bridges [2]. The bearing connections could be damaged by extreme waves, leading to the bridge superstructure being washed away [3] [4] [5] [6]. However, although a design methodology for the tsunami design of skewed bridges is available [7, 8], the estimation methods of storm-induced oblique wave load on coastal bridges are still not clear, and more investigations are necessary to improve existing design standards [9]. 

The wave incident angle is one of the critical issues affecting structural safety, which induces uneven load distributions and large overturning moments on the bridge superstructure. Previous studies on tsunami effects on coastal bridges have proven that large yaw and roll moments will be generated under extreme solitary wave impacts, but relevant studies on storm-induced wave effects are still limited. The hindcast analyses for Hurricane Katrina showed that most of the bridge damages were caused by oblique waves [10]. Among the few relevant studies [11] conducted experiments to measure the oblique wave loads on two typical bridge decks with slab-and-beam and box-culvert cross-sections, but only the wave with a period of 3. 0 s was tested. In recent years [12] performed numerical simulations of tsunami-induced forces on skewed bridges using the OpenFOAM software and found that the angle between wave direction and bridge deck would significantly increase the horizontal force [13] conducted laboratory experiments to investigate the influences of wave incident angles. It was also highlighted that the limited observations in that study could not provide a comprehensive understanding of such a complex problem, and further studies on this issue are required [14, 15]. 

Moreover, existing studies on skewed bridges [7, 8, 16] have proven that the angle of a bridge will lead to a much more complicated phenomenon, since it reduces the lateral force in the direction of the flow propagation (beneficial) but increases it in the other direction (detrimental) [14]. Another significant detrimental effect is the generation of yaw and roll moments, which result in an unequal distribution of both the horizontal and uplift forces in the two bent-caps and in the individual connections and bearings [15, 16]. Although the terms "skew angle of a bridge" and "incident angle of a wave" are not exactly the same things, they both generate a complex three-dimensional (3D) wave-bridge interaction, sharing similar difficulties when calculating the maximum loads on the bridge. In the case of a skewed bridge, the bent-caps and abutments are generally designed to be aligned with the flow (so the abutments are skewed relative to the direction of the span length) [14], while in the case of an 'oblique wave impact on a straight bridge' the bent caps are not aligned with the flow. Due to the misalignment between the direction of the substructure (bent-caps or abutments) with the flow, the wave-bridge interaction is more complex causing larger impulsive loads and different trends in the steady-state loads than a respective skew bridge. 

The oblique wave effects have not been thoroughly investigated, as appropriate numerical models have not been developed until recent years. There are some promising methods that have been used in the past to investigate the wave-bridge interaction such as the Arbitrary Lagrangian-Eulerian (ALE) method [17], the Smoothed Particle Hydrodynamics [18], the Particle Finite Element Method (PFEM) method [19], and the coupled SPH-FEM method [20]. However, the investigations on wave incident angles are still limited, which is mainly because of the fact that (i) it requires the development of sophisticated 3D models allowing waves from different angles; (ii) validation of these numerical models requires prior experimental data of such phenomena, which are very limited; and (iii) the computational fluid dynamics (CFD) process is computationally expensive requiring high performance computer systems. Nonetheless, in recent years, ref. [14] has demonstrated that an incompressible hydrodynamic solver based on the FEM can simulate accurately the solitary wave impacts on skewed bridges and capture all the 3D effects [21]. utilized a Finite Volume Method (FVM) to investigate the tsunami impact on skewed bridges. Moreover, during the wave-bridge interaction, there is an aleatory uncertainty in the wave breaking process [22]. A minor change in the initial water level or height of an unbroken wave can result in a major variability in the impact pressures and forces on a bridge [23], which adds complexity to these investigations. Inspired by these works, this study develops a 3D computational fluid dynamics (CFD) model to investigate oblique wave forces acting on coastal bridges. Different wave heights, wavelengths, relative clearances, and wave incident angles are tested. The maximum wave forces are quantified by using different surrogate models. These results can facilitate the following structural reliability analyses [24] [25] [26] [27] [28]. 

Overall, this study develops a 3D CFD model capable of computing oblique wave forces on coastal bridges. Time histories of wave loads in the x, y, and z directions under different wave heights, wavelengths, relative clearances, and incident angle conditions are computed. The maximum wave forces are quantified by different surrogate models. The remainder of this study is organized as follows: The 3D numerical modeling method is introduced in Section 2. A brief review of different surrogate models that are tested is given in Section 3. Oblique wave impacts and examinations of quantification metamodels are presented in Section 4. Finally, conclusions are drawn in Section 5.
section_name: Introduction

section_name: 3D Numerical Modeling for Wave-Bridge Interactions

section_content: In this study, the 3D CFD numerical domain was established in the ANSYS Fluent software, as shown in Figure 1a. After building the bridge superstructure model, subtract Boolean was generated with the fluid domain set as target bodies and structure domain set as tool bodies. Waves were generated in the velocity inlet plane ABCD and flowed toward the pressure outlet plane EFGH. Plane BCGF was also set as a pressure outlet with an atmosphere. Other planes were set as stationary walls. In the pressure outlet planes EFGH and BCGF, the free surface level model was utilized as the pressure specification and density interpolation method. The volume-of-fluid method was used as the multiphase model, and the open-channel flow model and open-channel wave boundary conditions were adopted. The SST k-omega model was used as the viscous model, and the turbulence damping factor was set as 10. The turbulence damping factor was set as 2%, and the turbulence viscosity ratio was set as 10%. 

given in Section 3. Oblique wave impacts and examinations of quantification metamodels are presented in Section 4. Finally, conclusions are drawn in Section 5.
section_name: Model Setups and Boundary Conditions

section_content: In this study, the 3D CFD numerical domain was established in the ANSYS Fluent software, as shown in Figure 1a. After building the bridge superstructure model, subtract Boolean was generated with the fluid domain set as target bodies and structure domain set as tool bodies. Waves were generated in the velocity inlet plane ABCD and flowed toward the pressure outlet plane EFGH. Plane BCGF was also set as a pressure outlet with an atmosphere. Other planes were set as stationary walls. In the pressure outlet planes EFGH and BCGF, the free surface level model was utilized as the pressure specification and density interpolation method. The volume-of-fluid method was used as the multiphase model, and the open-channel flow model and open-channel wave boundary conditions were adopted. The SST k-omega model was used as the viscous model, and the turbulence damping factor was set as 10. The turbulence damping factor was set as 2%, and the turbulence viscosity ratio was set as 10%. After performing mesh sensitivity analyses to satisfy the requirements of the Courant number [29], the tetrahedron meshes with a maximum size of 0. 6 m and a minimum size of 0. 07 m were utilized to form the whole numerical domain. Under this settlement, there are a total of 2,905,470 meshes in the numerical domain. The corresponding time step was set as 0. 0025 s. The CFL number during the CFD process is around 0. 25. The results of mesh sensitivity analysis are shown in Figure 1b. The numerical modeling method was validated by the authors' previous studies [30]. Compared with previous studies [31], in this study, we aimed to improve the structural details of the CFD model by modeling the I-shaped girders, diaphragms, and overhangs of the bridge. To reach this purpose, the meshes in the CFD domain were refined by adding wall layers along the fluid-solid interfaces, producing a smooth distribution along all the structure surfaces. It should be noted that the air and water are assumed as incompressible flows in the established numerical model. The trapped air has a critical role and increases the uplift forces [32] [33] [34] [35] [36], as well as the overturning moment [34] on the bridges. Moreover, refs. [35] and [36] indicated the generation of three-dimensional effects even in the case of a normal wave impact on a straight bridge (without any obliqueness) due to the tendency of the air to move in the 3rd direction. The effects of trapped air may be mitigated by setting air venting holes. This study mainly focuses on a 3D numerical modeling method for investigating storm-induced oblique waves. Additional experimental and numerical studies are encouraged to explore the 3D effects of trapped air in future work. 

The investigated bridge model in this paper is a typical simply supported bridge designed to carry two-lane traffic on the deck (as Figure 1a ), which is commonly used to connect island communities and suffered from severe damage during previous hurricanes After performing mesh sensitivity analyses to satisfy the requirements of the Courant number [29], the tetrahedron meshes with a maximum size of 0. 6 m and a minimum size of 0. 07 m were utilized to form the whole numerical domain. Under this settlement, there are a total of 2,905,470 meshes in the numerical domain. The corresponding time step was set as 0. 0025 s. The CFL number during the CFD process is around 0. 25. The results of mesh sensitivity analysis are shown in Figure 1b. The numerical modeling method was validated by the authors' previous studies [30]. Compared with previous studies [31], in this study, we aimed to improve the structural details of the CFD model by modeling the I-shaped girders, diaphragms, and overhangs of the bridge. To reach this purpose, the meshes in the CFD domain were refined by adding wall layers along the fluid-solid interfaces, producing a smooth distribution along all the structure surfaces. It should be noted that the air and water are assumed as incompressible flows in the established numerical model. The trapped air has a critical role and increases the uplift forces [32] [33] [34] [35] [36], as well as the overturning moment [34] on the bridges. Moreover, refs. [35] and [36] indicated the generation of three-dimensional effects even in the case of a normal wave impact on a straight bridge (without any obliqueness) due to the tendency of the air to move in the 3rd direction. The effects of trapped air may be mitigated by setting air venting holes. This study mainly focuses on a 3D numerical modeling method for investigating storm-induced oblique waves. Additional experimental and numerical studies are encouraged to explore the 3D effects of trapped air in future work. 

The investigated bridge model in this paper is a typical simply supported bridge designed to carry two-lane traffic on the deck (as Figure 1a ), which is commonly used to connect island communities and suffered from severe damage during previous hurricanes [37]. The bridge superstructure is 10. 45 m in width, 0. 3 m in depth, and 15. 85 m in length. There are six I-shaped girders evenly distributed along the width of the deck, and the height of each girder component is 1. 05 m. A total of three sets of diaphragms are settled at the two ends and the middle part of the deck, each with a thickness of 0. 5 m. 

The overall density of the concrete and steel is taken as 2. 6 × 10 3 kg/m 3. Hence, the total weight of the superstructure is around 3063 kN/span. The total wave loads on the bridge superstructure are calculated by accumulating the pressures on the structural surface of the bridge model. Detailed information on the bridge model could be found in ref. [14]. This study mainly focuses on the oblique wave impacts on the bridge superstructure and assumes that the substructure would not affect the forces on the superstructure. Effects of the substructures will be explored in future studies.
section_name: Model Setups and Boundary Conditions

section_content: Different water and wave scenarios were tested in this study to explore the oblique wave impacts. Several definitions are given first for the ease of following discussions. As shown in Figure 2a, the distance from the bottom of the bridge girder to the seabed was set as 16. 5 m. Different initial water depths before the wave arrived were tested, covering submerged (initial water level exceeds the bottom of girders, clearance Z c < 0 m) and unsubmerged (initial water level does not exceed the bottom of girders, clearance Z c ≥ 0 m) conditions. In the CFD simulations of wave-bridge interactions, linear waves were generated at the velocity inlet domain by setting multiphase open-channel wave boundary conditions. By rotating the bridge model in the computational domain, the wave incident angle θ can be adjusted, as shown in Figure 2b. The coordinate x 0 -y 0 -z 0 is the original coordinate system of the cuboid fluid domain, and the coordinate x-y-z is set for the measurements of oblique wave-induced forces on the bridge model in the horizontal (x), vertical (y), and longitudinal (z) directions. [37]. The bridge superstructure is 10. 45 m in width, 0. 3 m in depth, and 15. 85 m in length. There are six I-shaped girders evenly distributed along the width of the deck, and the height of each girder component is 1. 05 m. A total of three sets of diaphragms are settled at the two ends and the middle part of the deck, each with a thickness of 0. 5 m. The overall density of the concrete and steel is taken as 2. 6 × 10 3 kg/m 3. Hence, the total weight of the superstructure is around 3063 kN/span. The total wave loads on the bridge superstructure are calculated by accumulating the pressures on the structural surface of the bridge model. Detailed information on the bridge model could be found in ref. [14]. This study mainly focuses on the oblique wave impacts on the bridge superstructure and assumes that the substructure would not affect the forces on the superstructure. Effects of the substructures will be explored in future studies.
section_name: Tested Wave Cases

section_content: Different water and wave scenarios were tested in this study to explore the oblique wave impacts. Several definitions are given first for the ease of following discussions. As shown in Figure 2a, the distance from the bottom of the bridge girder to the seabed was set as 16. 5 m. Different initial water depths before the wave arrived were tested, covering submerged (initial water level exceeds the bottom of girders, clearance Zc < 0 m) and unsubmerged (initial water level does not exceed the bottom of girders, clearance Zc ≥ 0 m) conditions. In the CFD simulations of wave-bridge interactions, linear waves were generated at the velocity inlet domain by setting multiphase open-channel wave boundary conditions. By rotating the bridge model in the computational domain, the wave incident angle θ can be adjusted, as shown in Figure 2b. The coordinate x0-y0-z0 is the original coordinate system of the cuboid fluid domain, and the coordinate x-y-z is set for the measurements of oblique wave-induced forces on the bridge model in the horizontal (x), vertical (y), and longitudinal (z) directions. In addition, determining investigated wave cases is also an important step in quantification since it is rarely practical to study whole wave conditions. The tested water and wave cases should cover a wide range of scenarios that could occur in natural disasters, In addition, determining investigated wave cases is also an important step in quantification since it is rarely practical to study whole wave conditions. The tested water and wave cases should cover a wide range of scenarios that could occur in natural disasters, which is conducive to analyzing how the oblique wave impacts can be attributed to the influences of various input factors (water depth, wave height, etc. ), as well as the effects of interactions between these inputs [38]. To serve this purpose well, a random sampling method that provides good coverage of the parameter space and avoids an unconscious bias reflected in the collected data is required [39]. A superior input group is also helpful to facilitate the metamodels to capture data characteristics, optimize the training speed, and improve the applicability of the trained model [40, 41]. Recognizing these, the Sobol sequence method [42] was adopted to sample different wave parameters to form the input group, which is a high-efficient quasi-Monte Carlo sequence and has been proven effective in reducing model learning costs [43, 44]. 

While the established numerical model can be used to calculate oblique wave impacts on coastal bridges, these structural demands are not deterministic quantities in real hazards because of aleatory and epistemic uncertainties in the modeling parameters [4]. Probability distributions are, therefore, applied for different parameters considered to be random variables. In this study, we used the simulated sea surface directional wave spectra under hurricane wind forcing [45], which was validated based on buoy observations and NASA Scanning Radar Altimeter (SRA) data obtained during the 1998 Hurricane Bonnie. Accordingly, uniform distributions were used to sample clearance Z c (-2. 5~3. 5 m), wave height H (4~10 m), and wavelength λ (100~250 m). In addition, because of the symmetry of the bridge structure, the incident angle θ was (0~90 • ). The sampled water and wave cases are listed in Table 1, and the distributions of each parameter are illustrated in Figure 2c.
section_name: Tested Wave Cases

section_content: This section provides an introduction to different surrogate models examined in this study for the prediction of oblique wave forces. The CFD computational results were used as data sources for the model training. These advanced mathematical models can help to approximate and extrapolate the hydrodynamic features during wave-bridge interactions.
section_name: Surrogate Models for Quantification of Oblique Wave Forces

section_content: The Gaussian Kriging surrogate model is a stochastic process model furnished with features of uncertainty quantification. It has been widely applied in many engineering disciplines [46]. Supposing a function to be approximated has n sampled points and k variables, (as x i = (x i (1),. . ., x i (k) ) for i = 1,. . ., n), and the associated observation value y i = y(x i ), the predicted value Y i can be expressed as

where each f l (x i ) is a linear or nonlinear function of x i ; β l is the coefficient corresponding to the lth variable; ε(x i ) is randomly distributed with a mean zero and a variance equals σ.
section_name: Gaussian Kriging Surrogate Model

section_content: The artificial neural network (ANN) model is a powerful computing system based on a collection of connected nodes, called artificial neurons, connecting the inputs X and outputs Y [47]. The general structure of an ANN comprises an input layer, an output layer, and several hidden layers formed by neurons. The nonlinear relations between X and Y are modeled through the connections between the neurons, forming a comprehensive network. The output of the jth neuron in the (I + 1)th layer a j (i+1) is a function of the linear combinations of the outputs from the neurons in the previous layer

where h() is the activation function (the commonly used sigmoid function); a m (i) is the output from the previous layer; w m (i) is the weight; b (i) is the bias term.
section_name: Artificial Neural Network (ANN)

section_content: The polynomial chaos expansion (PCE) is a widely used metamodel constructed based on the spectral representation of the performance function [48, 49]. The key concept of PCE is that the outputs can be expanded to a basis composed of several multivariate polynomials [50]. A typical PCE model can be expressed as

where α is a set of indices mapping to the components of H α (X); C α is coefficients; φ αi

is the univariate orthogonal polynomial with respect to the ith variable in degree α i. For practical purposes, the original PCE is truncated to a finite sum by setting the total degree of all the polynomials associated with the input variables smaller than or equal to p as

where A p,k is the truncation set. In this study, a three-order PCE model (i. e., p = 3) was adopted after several calculations and comparisons, and the total number of coefficients in a PCE (P) could thus be calculated as

The unknown coefficients C α (or C) are commonly determined by the projection approach or by the regression approach. In this paper, the least-square minimization method (regression approach) was adopted to calculate the unknown coefficients C by minimizing the mean-square truncation error [50].
section_name: Polynomial Chaos Expansion (PCE)

section_name: Oblique Wave-Induced Forces and Metamodeling

section_content: Figure 3 shows one typical example of wave-bridge interaction result for the case of clearance Z c = -2. 39 m, wave height H = 9. 27 m, wavelength λ = 226 m, and incident angle θ = 77. 31 •. In this case, the bridge is already submerged beneath the initial water level by storm-induced surges, and massive waves continue to hit the bridge superstructure, generating periodically changing loads. Additionally, there is a large incident angle between the bridge model and the wave forwarding direction, so that it would generate changing wave forces in both horizontal and longitudinal directions. The wave-bridge interactions at representative moments with t = 8, 10, 12, and 14 s are illustrated in Figure 3a-d From Figure 3, it is observed that as the water surface rises, the vertical wave force F y gradually increases to 3192 kN, as shown in Figure 3e. F y reaches its positive peak at around t = 10 s, while the maximum values of F x and F z occur a little later. Although the bridge is fully submerged at the initial stage and subjected to a large buoyance, F y reduces to a negative value (around -1565 kN) as the water level drops to the trough (not the static weight of the structure). Such wave-induced downward load is relatively small, compared with normal traffic load, which means the bridge is less likely to be crushed by wave impacts. However, it is worth noting that the maximum longitudinal force F z-max is larger than the maximum horizontal force F x-max (288 kN and 506 kN, respectively). In the general design of such kinds of simply supported bridges, the constraints for the superstructure are mainly contributed by the bearing connections beneath the girders [51]. The connections are assumed as elastomeric bearings for the investigated bridge. The limit state is defined as the state where the failure of the bearing occurs. Once the bearing is damaged, it needs to be repaired immediately, otherwise, it will affect the normal traffic use. These bearings usually have a large capacity in the vertically downward direction (y) and horizontal direction (x), but the resistance in the longitudinal direction (z) is relatively small. This means the bearing connections could be damaged by the wave forces in the longitudinal direction (z), and once damaged, they will lose their constraints in the other directions (x and y), leading to the deck overturning or unseating failure [51]. Therefore, it is of vital importance to quantify the varying loads on coastal bridges subjected to oblique wave impacts and consider these effects in the design stage.
section_name: Time Histories of Oblique Wave Forces

section_content: For an in-depth understanding of the oblique wave impacts, the maximum forces (Fxmax, Fy-max, and Fz-max) were obtained from the time histories of the wave forces and are plotted in Figures 5 and 6. The computed results of the tested cases are listed in Table 2. The colors of the circular marks in the figures correspond to the wave incident angles, and the sizes are proportional to the values of wave forces. Several typical results are marked with arrows. Figure 5 shows the results under unsubmerged conditions. Generally, Fz-max
section_name: Maximum Wave Forces under Different Scenarios

section_content: For an in-depth understanding of the oblique wave impacts, the maximum forces (F x-max, F y-max, and F z-max ) were obtained from the time histories of the wave forces and are plotted in Figures 5 and 6. The computed results of the tested cases are listed in Table 2. The colors of the circular marks in the figures correspond to the wave incident angles, and the sizes are proportional to the values of wave forces. Several typical results are marked with arrows. Figure 5 shows the results under unsubmerged conditions. Generally, F z-max has a much larger magnitude than F x-max and F y-max in most cases and could exceed the static weight of the bridge superstructure under extreme conditions, which is around 3063 kN per span [51]. The values of F x-max and F z-max highly depend on the wave incident angle. F z-max is larger than F x-max in most cases when θ > 60 •, as marked by the yellow and orange circles in Figure 5c. Results under submerged conditions are presented in Figure 6. Similar characteristics are observed except that the peak values of F y-max are much smaller than those under unsubmerged conditions, while the results of F x-max and F z-max are very close in these two scenarios. Although the bridge is affected by buoyancy when it is submerged, the smaller F y-max indicates that the vertical uplift force is mainly contributed by the impulse impact under unsubmerged conditions. Conversely, the effects of F z-max could be more critical to structural safety under oblique wave impacts considering the insufficient constraints in the longitudinal direction (z). Besides, the bridge is assumed as a rigid body in the established numerical model, which means the structural dynamics are not considered. The dynamic characteristics have been proven significant in previous studies on tsunami waves and breaking waves [52] [53] [54]. These dynamic structural characteristics should also be investigated for oblique wave impacts in future studies. The established numerical model provides reliable results under different wave conditions but may lead to expensive or even prohibitive computational costs when it requires large-scale calculations, such as fragility and risk assessment. Thus, how to quantify the complex relationships between the maximum wave forces and these wave parameters to the greatest extent is also a critical step. Conventional methods estimated wave loads using empirical formulas summarized based on experimental or numerical data; however, this method cannot accurately measure the effects of each parameter and their interactions. Additionally, the introduction of the incident angle makes it more complicated. Recognizing this, this study used advanced stochastic techniques to predict the maximum wave forces under different wave scenarios. In this section, the performance of different surrogate models (the ANN model, Kriging model, and three-order PCE model) are examined and compared. The coefficient of determination (R 2 ) and root-mean-square error (RMSE) were adopted as the goodness of prediction. 

In this study, the ANN metamodel was trained using the stochastic gradient descent method, which uses randomness to find the optimal weights and bias from the inputs to the outputs. The network weights were first initialized randomly. To this end, 70% of the datasets were used to train the neural model, 15% were used as a validation set, and the remaining 15% were retained to test the trained data. After several calculations and comparisons, 3 hidden layers were adopted, and 10 neurons are used in each layer. The Bayesian Regularization algorithm was employed as the transfer function. A typical example of the gradient descent during the training process for F y-max is shown in Figure 7. The gradient gradually converged after 100 epochs. After several calculations and comparisons, the fitting results and prediction performance for Fx-max, Fy-max, and Fz-max by using different metamodels were determined, which are listed in Table 3. For most of the cases, the Kriging surrogate model has a smaller R 2 than the other two models. The Kriging model fails to capture the characteristics of the sampling data well, which also means a larger volume of data is required for a better prediction performance, increasing the computational cost. On the other hand, the threeorder PCE has the largest R 2 in most cases, and the results of the ANN model are in between. The RMSE values from the three-order PCE and ANN models are close when predicting Fx-max and Fy-max. For Fz-max, the ANN model results in a much smaller RMSE, compared with the PCE model, indicating a smaller prediction bias and better performance. Thus, the ANN model is adopted for the prediction of maximum wave loads in the following fragility analyses. It should be noted that the trained surrogate models are only applicable to the investigated bridge model, which mainly considers the effects of different wave parameters. Previous studies have proven that the bridge dimensions also have significant effects on the wave-induced loads [4, 32]. Thus, more investigations on different structural dimensions are recommended in future studies, to ultimately reach a flexible prediction model which could be applied for different bridge and hazard scenarios. After several calculations and comparisons, the fitting results and prediction performance for F x-max, F y-max, and F z-max by using different metamodels were determined, are listed in Table 3. For most of the cases, the Kriging surrogate model has a smaller R 2 than the other two models. The Kriging model fails to capture the characteristics of the sampling data well, which also means a larger volume of data is required for a better prediction performance, increasing the computational cost. On the other hand, the three-order PCE has the largest R 2 in most cases, and the results of the ANN model are in between. The RMSE values from the three-order PCE and ANN models are close when predicting F x-max and F y-max. For F z-max, the ANN model results in a much smaller RMSE, compared with the PCE model, indicating a smaller prediction bias and better performance. Thus, the ANN model is adopted for the prediction of maximum wave loads in the following fragility analyses. It should be noted that the trained surrogate models are only applicable to the investigated bridge model, which mainly considers the effects of different wave parameters. Previous studies have proven that the bridge dimensions also have significant effects on the wave-induced loads [4, 32]. Thus, more investigations on different structural dimensions are recommended in future studies, to ultimately reach a flexible prediction model which could be applied for different bridge and hazard scenarios.
section_name: Maximum Wave Forces under Different Scenarios

section_content: In this study, an elaborate 3D CFD numerical model was developed to investigate oblique wave impacts on coastal bridges. Time histories of wave-induced forces in the horizontal (x), vertical (y), and longitudinal (z) directions were computed, and the maximum values were quantified by using different surrogate models. The developed 3D CFD model has superior performance in calculating wave load distributions on coastal bridges for its spatial advantages and structural details. More relevant studies are required to study the responses of bearing constraints under changing loading distributions. The maximum oblique wave forces were quantified by using ANN surrogate models with different wave parameters, including wave height, relative clearance, wave incident angle, and wavelength. 

Overall, this study provides a preliminary numerical investigation of oblique wave impacts on coastal bridges, providing data reference for future studies and design standards. More experimental and parameterizing studies on oblique wave impacts are encouraged to improve existing design specifications.
section_name: Conclusions

section_name: Acknowledgments:

section_content: The opinions and conclusions presented in this paper are those of the authors and do not necessarily reflect the views of the sponsoring organizations.
section_name: Acknowledgments:

section_name: Institutional Review Board Statement: Not applicable.

section_content: Funding: This research was funded by the National Key R&D Program of China (No. 2019YFB2102703 ).

section_content: Informed Consent Statement: Not applicable.
section_name: Institutional Review Board Statement: Not applicable.

section_content: Data Availability Statement: Not applicable.

section_name: Conflicts of Interest:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: