Application of a Statistical Regression Technique for Dynamic Analysis of Submarine Pipelines

: This study employs a statistical regression technique to investigate the maximum displacement, stress, and natural vibration frequencies of a submarine pipeline subjected to hydrodynamic wave forces. Eighteen pipeline models are designed, varying in wall thickness from 10 mm to 30 mm and diameter from 500 mm to 1000 mm. The hydrodynamic drag and inertia forces are performed by using the Morison equation. Computer-aided Finite Element Analysis is employed to simulate the complex interactions between the fluid and structure in 18 pipelines. Multiple Regression technique is used to evaluate the reliability metrics, considering uncertainties in geometrical properties affecting pipeline performance. Full Quadratic models are developed for expressing more effective and concise mathematical equations. Analysis of Variance (ANOVA) is performed to determine the adequacy of the model in representing the observed data. The Coefficient of Determination (R 2 ), Mean Square Error (MSE), and Mean Absolute Error (MAE) are calculated to assess the equation’s predictive accuracy and reliability. The results confirm the suitability of the suggested regression technique for analyzing the relationships between predictor variables and the response variable.


Introduction
The design of submarine pipelines involves a comprehensive approach that integrates knowledge from multiple disciplines, aiming to ensure the long-term integrity and functionality of critical components in the oil and gas industry's transportation infrastructure.Environmental loads, including wave forces, are carefully analyzed, with particular attention to the physical properties of the water and the characteristics of the pipeline itself, such as diameter (D), thickness (t), and material composition.Additionally, factors like wave height (H), wave period (T), and water depth (d) play a significant role in influencing hydrodynamic processes and the dynamic behaviors of pipelines.Researchers have proposed different methodologies, such as numerical modeling, experimental testing, field measurements, and fault detection techniques in subsea environments, to analyze and mitigate stability concerns [1][2][3][4][5][6][7].Repeated and complex analyses are performed under variable conditions to determine the structural design.
In recent years, where the importance of achieving results in a shorter time and with fewer resources has increased, the use of statistical prediction methods in the analysis of submarine pipelines has also increased.Youssef et al. employ a statistical method to create a response surface model that predicts the maximum horizontal displacement of a pipeline during storm conditions.The Monte Carlo simulation technique is applied alongside the established response surface model to determine extreme response statistics [8].Xu and Sinha introduce a thorough analytical structure for statistically analyzing field performance data for water pipelines.The paper provides a detailed explanation of the methodology's implementation steps, along with preliminary analyses conducted on datasets from two different water utility systems [9].The study by Zhang and Weng presents a Bayesian network model for analyzing buried gas pipeline failures caused by corrosion and external interference.The model offers a probabilistic framework for assessing the risk of pipeline failure, taking into account various contributing factors and their uncertainties [10].Zhang et al. develop an optimal statistical regression model capable of predicting the wave-induced equilibrium scour depth beneath pipelines in sandy and silty seabeds [11].Most of the studies conducted on offshore structures are based on ground movements occurring on the seabed [12,13].In this study, a suspended pipeline model is examined to avoid the impact of seabed movement.
This study aims to provide a more comprehensive understanding of the factors influencing pipeline dynamic behavior and to enhance the accuracy of stability assessments by using statistical methods.This study focuses on the structural model rather than seabed movement and material properties.This approach allows for a comprehensive assessment of pipeline stability under hydrodynamic wave forces, enabling informed decisionmaking in design processes.
The multiple regression model, which is most commonly used for complex decision problems in many fields [14], is employed to demonstrate the general tendencies of the relationship between the pipeline design parameters and stability criteria.The hydrodynamic wave forces are obtained by using Airy Wave Theory.Le Méhauté's diagram is used to determine the applicability theories of water waves according to the values of normalized wave height H and the water depth d [15].The hydrodynamic drag force (FD) and inertia force (FM) are performed by using the Morison equation.Since the ratio of the gap between the pipe and the seabed to the pipe diameter is considered to be greater than 1.0, the hydrodynamic lift force (FL) has been neglected [16].The vibration and stability of a simply supported steel pipeline are analyzed when lateral hydrodynamic forces (FH) act on it.For this, the dynamic behavior of the pipeline is modeled by computer-aided Time History Analyses based on Newmark β [17].The pipe wall thickness ranges from 10 to 30 mm, and the pipe outer diameter ranges from 500 to 1000 mm.The initial three natural frequencies, maximum values of stress, and displacement are used to create a dataset to explore relationships between variables for 18 pipe models.The multiple regression method is employed to assess the reliability metrics while accounting for uncertainties in the geometric properties that influence pipeline performance.Full quadratic models are constructed to present more efficient and concise mathematical formulations.An Analysis of Variance is conducted to ascertain the model's adequacy in accurately representing the collected data.The performance of the regression methods is evaluated based on different metrics: Mean Squared Error (MSE), Mean Absolute Error (MAE), and Coefficient of Determination (R 2 ).The results highlighted that the proposed regression methods can be an alternative practical way that provides the correct solution in a shorter time in the preliminary design of submarine pipelines.The relationship between the variables can be implemented in any computational calculation software.

Analyzing Procedure
In this paper, a 2D structural model is performed to predict the dynamic behavior of submarine pipelines under wave forces according to displacements, maximum stresses, and natural vibration frequencies.Airy Wave Theory is used to obtain wave forces, which are adopted as environmental loads.Numerical simulations are also applied using the Time History Tool of SAP 2000.The natural vibration frequencies, maximum values of stress, and deflection are derived to create a dataset for the regression process.The multiple regression model is used to learn relationships between four categories of input and five categories of output from example data and to predict novel inputs.Moreover, 144 input values are used for estimating the relationships between dependent and independent variables for each case.The performance metrics of full quadratic equations have demonstrated that regression analyses can be used for forecasting and predicting submarine pipeline reactions based on geometrical properties and effective hydrodynamic forces.

Structural Model Application
In application, the predominant configuration of suspended submarine pipelines consists of multiple spans (refer to Figure 1), with support provided at intermediate locations through various support mechanisms.In this study, the pipe is selected as similar to the one used in Gücüyen [18], having Lp = 10 m span length and consisting of uniform and homogeneous material.The structural behavior of the 18 pipe models is achieved under regular waves.The outside diameters of the pipe range from ϕ500 to ϕ1000, and the wall thicknesses are chosen as 10 mm, 20 mm, and 30 mm to verify the effect of cross-section on frequency, stress, and displacement values.The properties of the pipe models are given in Table 1.The material properties of the pinned-pinned pipe are chosen to represent the steel, with Young's modulus of 21 × 10 7 kN/m 2 , Poisson's ratio of 0.3, and density of 78.50 kN/m 3 .Various combinations of pipe geometrical properties are performed using Time History analyses.The outer diameter of pipe models influences the value of hydrodynamic forces [19,20].However, the wall thickness changes the behavior of the pipe in response to the force.

Environmental Conditions and Load Assessment
The submarine pipeline is designed by considering a combination of dead loads and wave loads.The environmental conditions are modeled in the numerical analysis based on Airy Wave Theory, referring to Le Méhauté's diagram, which indicates the validity domain of each theory-type of water waves [21].
Water depth, wave period, and height are three major parameters in the design of marine facilities [22].The employed parameters d = 10.00 m, T = 8.00 s, and H = 0.31 m are considered for the applicability criteria of Airy Wave Theory, which is widely used to model a marine environment [23].The wave loads are calculated with water particle velocities (u) and accelerations (̇) by two different approaches presented as follows: where the z-axis is directed vertically upward from the still water level, which is positive.The x-axis is along the direction of the propagation of waves.When the value of the relative water depth (water depth/wavelength, d/L) is between 0.005 and 0.05, it is classified as an intermediate depth wave, as mentioned in this study [24].The free surface elevation η during the wave generation process is obtained based on the potential flow approach and is given in Figure 2. The seawater characteristics are assigned as fluid properties, with a salty water density () of 1025 kg/m 3 and a dynamic viscosity (μ) of 0.0015 Ns/m 2 .The wavelength (L) is determined as ≈70.90 m by considering wave parameters.

Description of Hydrodynamic Loads Acting on Pipeline
Drag, lift, and inertia forces are hydrodynamic loads caused by the motion of the surrounding water in a pipe.Drag and inertia are defined as the components that are parallel to the flow direction, and lift is the component of the hydrodynamic load that is perpendicular to the flow direction.In this study, the magnitude of the gap ratio between the pipe and the seabed is assumed to be e/D ≈ 1.2.It is observed that when e/D > 1, the importance of the lift force vanishes compared to the in-line forces.This means that the importance of the cross-flow vibration becomes negligible, and in-line vibrations dominate the behavior of the pipe [25].Hydrodynamic stability is established by using the Morison equation, which relates hydraulic inertial (FI), drag (FD) forces to local values of particle velocity and acceleration.The total force (FH) can be derived by integrating Equation (3) along the z-axis as follows: 1 24 where D is represented as the diameter of the pipe.As seen in Equation ( 3), the force components are required to determine the drag force coefficient CD and the inertia force coef-ficient CM.In this paper, the unknown coefficients are assumed to be 2.0 and 0.7, respectively [26].The dynamic behavior related to changes in hydrodynamic forces over a wave period is analyzed by collecting data at 1 s intervals.The examinations of total hydrodynamic forces are made on the basis of the definition of Airy Wave Theory, and the results are given in Figure 3.The hydrodynamic forces are assigned as time-varying external loads on the pipeline in the finite element analysis.

Time History Analyses
Time history analysis is a sophisticated technique used in structural engineering to assess the dynamic response of structures under various time-dependent loads [27].This method involves inputting a detailed time series of loading conditions into a computer model of the structure, which then calculates the resulting stresses and displacements at each time step.
The horizontal hydrodynamic forces (FH) resulting from FD and FM are performed in computer-aided Time History analyses based on Newmark β [17].Time History analyses provide for both linear and nonlinear structural dynamic responses and sinusoidal and non-sinusoidal loading.The time step method, which is supported by computer-aided Time History analyses, determines the solution of an equation at a succession of values of t, t + ∆t, t + 2∆t, etc.For finite element formulation, it was decided to use the well-known commercial software SAP2000, which is based on the implicit time step method.Various implicit methods are available.The widely used Newmark β method is used in the method given below.Starting from a Taylor series [28] as follows: Here, m is mass, k is stiffness, c is damping, and F is the loading term.The accuracy of the solution depends on the length of the step interval (∆t).It must be short enough for the load time history, the response time history, and, in many cases, the shortest natural periods to be well defined [29].The analysis is executed over a duration of 30 s, employing a step interval of 0.01 s and a β value of 0.5.Numerical analysis uses the time-varying external loads as time history functions presented in Figure 3. Various combinations of pipe geometrical properties and hydrodynamic loading conditions are performed using Time History analyses, which provide a time-dependent function for obtaining dynamic structural response under wave forces.The natural frequencies for the initial three modes, the maximum lateral displacement of the pipeline, and the maximum values of stresses are used to generate a set of 144 input-output data pairs.

Regression Techniques Implemented
Regression, as denoted in the statistical literature, contains the prediction or learning of numeric features [30].Regression holds significance across numerous engineering applications, as a multitude of real-world issues can be effectively addressed through the modeling of regression problems.Extensive research has been conducted in the field of maritime transportation [31][32][33].
Multiple regression (MR) is commonly used to investigate the causal relationship between multiple independent variables and a dependent variable simultaneously based on mathematical equations [34].Each independent variable has its own coefficient, representing its unique contribution to the prediction of the dependent variable, while holding other variables constant.The algorithm calculates coefficients for the equation in a way that minimizes the sum of squared errors between the predicted values and the actual values [35].The general mathematical form of multiple regression with "n" predictors is expressed as follows [36]: Here, y is the dependent variable, xj is an independent variable, j represents the regression coefficient associated with each regressor, ₀ is the constant term, and  denotes the random error term.
The prediction accuracy of the regression model is investigated by performance metrics containing Coefficient of Determination (R 2 ), Mean Square Error (MSE), and Mean Absolute Error (MAE).
The performance of all regression models is also investigated based on the Coefficient of Determination (R 2 ) as given by the following expression: R 2 measures the effectiveness of a model in fitting its data and presents variances in the measured variable.Its scale spans from 0 (indicating no explanatory power) to 1 (representing a perfect fit).
MSE determines the variations between predicted and actual values by averaging their squared disparities.Smaller values of MSE signify higher algorithm performance; however, it is susceptible to the influence of outliers [33].
( ) MAE evaluates the mean variation between predicted and actual values, treating all errors uniformly.It presents lower sensitivity to outlier data, rendering it advantageous in scenarios where substantial errors are undesirable [37].
where the number of observations is represented by n.Within this manuscript, six main stages are used in the system architecture, as presented in Figure 4.These steps are assigned distinct and crucial responsibilities within the system.

Natural Vibration Frequencies, Displacements, and Stresses of Pipeline Models
The analysis is conducted for a pipeline with a single opening and simple supported ends using the Euler-Bernoulli beam theory.The basis of the Euler-Bernoulli beam theory relies on three assumptions of continuum mechanics: first, that straight lines initially normal to the beam's axis retain their straightness after deformation; second, that these lines are inextensible; and third, that they rotate as rigid lines perpendicular to the bent axis of the beam [38].The axis x of the pipe is defined as a longer dimension, and the transverse vibration of the pipe that occurred on this axis is calculated.
The pipeline system's dynamic performance can be described by identifying one or more natural vibration frequencies [39].The high-frequency modes exhibit a more rapid decay in amplitude compared to the slower decay observed in the low-frequency modes.The dominant behavior is caused by the lower frequency modes [40].Hence, the natural vibration frequencies of the pipe have been calculated for the first three modes in this study.Figure 5 shows the specific displacements required for the first three modes, initiating harmonic vibrations in the pipe.The comparative results in Table 2 present the obtained natural vibration frequency values from the analysis conducted for a single opening with simply supported ends of the pipeline.In assessments where the internal flow in the pipe is neglected, it is observed that the natural vibration frequency values increase with an increase in the pipe diameter.Additionally, for each diameter, three different values of wall thickness also affect the frequency.It is evident that as the wall thickness increases, the results obtained for the same diameter decrease.The highest natural vibration frequency occurs when the wall diameter is 1000 mm and the wall thickness is 10 mm.The maximum values of stress and displacement are used to create a dataset to explore relationships between variables for 18 pipe models.The maximum stress values contain the most significant strains occurring at the section where the maximum displacement is observed.The hydrodynamic forces vary depending on the velocity of fluid particles changing direction throughout the wave period.This leads to stresses calculated at the cross-section of the pipeline having different effects, with positive (+) and negative (−) pressures, as given in Figure 6."BS EN 1993-1-1:2005+A1:2014 Eurocode 3. Design of steel structures General rules and rules for buildings" [41] and "BS EN 10025-2:2019 Hot rolled products of structural steels Technical delivery conditions for non-alloy structural steels" [42] are selected as steel standards in this study.The stress levels for all models remain within permissible limits.The blue, red, and black colors represent the results for wall thicknesses of 10 mm, 20 mm, and 30 mm, respectively.As observed in the figure, results obtained at t = T/2, T/4, 3T/4, and T are compared.The stress values in the pipeline are highest at 2 s and 6 s, where wave crests and troughs occur, representing the peak hydrodynamic forces.An increase in wall thickness enhances strength, leading to a decrease in stress.However, the decrease in stress is not significant due to the increase in hydrodynamic force according to the Morison equation resulting from the increase in diameter.The maximum displacement that can occur in the cross-section of the pipeline is as crucial as stress.Particularly in systems where the gap between the pipeline and the seabed is small, erroneous calculations can lead to the deformation of the pipeline and significant financial losses [43].Comparative results of maximum values are presented in Figure 7, considering the time intervals as t = T/2, T/4, 3T/4, and T. Maximum displacements are observed again at the wave crest and trough positions.However, unlike stress values, it is evident that both the diameter and wall thickness values contribute to the decrease in displacement.Evaluating the scenario where displacement takes small values aims to ensure safety within the limits under the internal flow conditions, which will be considered in subsequent studies.The largest displacement, calculated as 0.113 mm, occurs at t = 6 s with D = 500 mm and tw = 10 mm.The smallest value, 0.003 mm, is determined at t = 8 s when the surface profile returns to a calm sea surface, with D = 1000 mm and t = 30 mm.The relationship between datasets consisting of maximum stress and displacement values, independent of time, is presented in Figure 8.

Comparative Results of Regression Models
In this section of this study, ANOVA (Analysis of Variance) results are generated.The effect magnitude of input parameters (D, tw, FH, and t) on the outcomes (wmax, σmax, ω1, ω2, and ω3) is determined.Thus, the impact ratios of pipe diameter D, pipe wall thickness tw, hydrodynamic force FH, and time t are evaluated.In analyses where time is considered as input, the t/T approach is used to reach a general correlation.By using the dimensionless time term, the mathematical expression obtained ensures validity for wave models with different wave periods.Full quadratic models that include not only linear terms but also quadratic terms are developed for expressing more effective and concise mathematical equations [44].The squares of the independent variables appear as part of the quadratic terms.These quadratic terms are included to capture nonlinear relationships between the predictor variables.It is expressed mathematically as follows: where Y is the dependent variable, and X1, X2, X3, …, Xk are the predictor variables.β0, β1, β2, β3, …, βk are the coefficients for the linear terms; β11, β22, β33, …, βkk are the coefficients for the quadratic terms; β12, β13, …, β(k−1)k are the coefficients for the interaction terms; and ε represents the error term.
Additionally, the effects of input parameters on the outcome parameters are demonstrated with 3D graphics.

The Effect of Input Variables on Displacement Values
Displacement is an important parameter in the strength of the pipe.Knowing the extent to which each input affects displacement is of great importance.New mathematical models have been created using statistical methods to solve this problem in pipeline transportation systems.The results of the ANOVA analysis are provided in Table 3.As seen in Table 3, within a 95% confidence interval, FH, D × FH, and tw × FH are the most influential parameters on the maximum displacement occurring throughout the independent period.At the same time frame, both the diameter and wall thickness of the pipe have a significant impact on displacement values.However, when a general regression model independent of time is constructed, hydrodynamic force and the conditions under which this parameter is multiplied become important.
The full quadratic model equation obtained using these parameters is given in Equation (10).wmax = −0.0745+ 0.000061D + 0.001261tw + 0.003184FH + 0.1330t − 0.000002D × FH is the most influential parameter on maximum displacement, accounting for 29.476% of the effect.The effect ratio for D × FH is 7.727%, while the effect ratio for tw × FH is 2.511%.The Pareto chart related to effective parameters is given in Figure 9.The regression model achieves a high R 2 of 8.952 × 10 −1 , which is close to the 1.0 perfect fit.Additionally, the model's MSE value is 1.182 × 10 −2 and its MAE value is 8.081 × 10 −2 .The results show that the model can accurately predict the maximum displacements of the pipeline.The relationship between maximum displacement values and input data is presented in Figure 10.It is observed that hydrodynamic forces reach their maximum values under the crest and trough of the wave surface profile.When the hydrodynamic force is at its maximum, displacement values also increase.The increase in displacement values is more pronounced when the wall thickness is taken as 10 mm.Displacement values decrease as the pipe diameter increases.For the smallest selected diameter of 500 mm, the pipe experiences more oscillation due to the hydrodynamic force, resulting in a displacement change figure resembling the wave surface profile.
For the largest selected diameter of 1000 mm, the displacement values are maximum at the (−) and (+) maximum points of the hydrodynamic force, but the change curve reaches a more linear form.This indicates that increasing the pipe diameter and thickness reduces the destabilizing effect on stability in the transition regions of the hydrodynamic force causing tension and compression.The results obtained also highlight the importance of independent variable effect coefficients determined by ANOVA analysis in the design process.

The Effect of Input Variables on Stress Values
Stress values serve as key indicators of the structural performance of submarine pipelines, and their proper evaluation and management are essential for ensuring the safety, reliability, and sustainability of these critical infrastructure systems [45].The statistical analysis results, aimed at succinctly and practically illustrating the variation of this important parameter with respect to the geometric properties of the submarine pipeline, are presented in Table 4.As seen in Table 4, FH is the most influential parameter on maximum stress values, accounting for 29.490% of the effect.The effect ratio for D × FH is 4.540%, while the effect ratio for tw × FH is 4.224%.In the mathematical model, the properties of the pipeline are represented by 2-Way Interaction terms: D × FH and tw × FH.The full quadratic model equation is presented in Equation (11).σmax = 0.2530 − 0.000188D − 0.00503tw − 0.016386FH − 0.450t + 0.000010D × FH + 0.000334D × t + 0.000164tw × FH + 0.00894tw × t (11) The Pareto chart related to effective parameters is given in Figure 11.The sudden changes in maximum stress values occurring under the wave crest and trough decrease with increasing diameter values.As the cross-sectional area of the pipeline increases, stability is enhanced, resulting in smaller oscillations in frequency.While the increase in pipe diameter leads to a higher hydrodynamic inertia force at t = 0, the resulting reductions in maximum stress values are relatively lower.However, over time, a significant decrease in stress becomes apparent.
The time-dependent maximum stress variation figure is obtained to be consistent with the water surface profile.The maximum stress values decrease with increasing wall thickness.The greatest difference between positive and negative stress values is obtained under the conditions of D = 500 and tw = 10 mm.Significant variations in maximum tensile and compressive stresses can lead to fatigue of the pipeline over time and the occurrence of additional stresses.

The Effect of Input Variables on Natural Vibration Frequencies
The natural vibration frequencies of submarine pipelines are essential considerations in ensuring their structural integrity and resilience against dynamic forces in underwater environments [46].Hence, in this study, a variation analysis of the natural vibration frequencies for the initial three modes is conducted.The results of the Analysis of Variance (ANOVA) are presented in Table 5.Typically, the lowest frequency mode, termed the first mode, is characterized by a uniform bending motion along the length of the pipeline.As seen in Table 5, the effect of D in the regression model is calculated as 99.768%, while the effect of tw is computed as 0.214%.In contrast, the second mode presents a single-loop oscillation, resembling a halfwave pattern.Therefore, in addition to the effects of D at 99.718% and tw at 0.215%, the term D × D also exhibits an effect of 0.066%.The third mode introduces an additional loop, leading to a more intricate vibration pattern.This also increases the number of terms required to define the mathematical model.In addition to D being the most effective parameter at 90.280% and tw at 1.868%, the square (D × D, tw × tw) and 2-Way Interaction (D × tw) parameters also gain importance.
The equations for the full quadratic model, obtained to describe the relationships between variables, are presented for the initial three natural vibration frequencies in Equations ( 12)- (14).
The model given by Equation ( 12) exhibits a strong R 2 value of 9.998 × 10 −1 , indicating that it explains a significant portion of the variance in the data.Moreover, MSE is low at 1.207 × 10 −3 and MAE is 2.849 × 10 −2 .
The second vibration frequency model accounts for 100% of the variance in the dependent variable.It means that the model explains all the variability of the dependent variable around its mean using the independent variables.Additionally, both the MSE and MAE are low, recorded at 9.283 × 10 −3 and 8.667 × 10 −2 , respectively.ω3 = 11.39 + 0.0463D + 0.153t + 0.000041D × D + 0.01720t × t − 0.001441D × t (14) The model obtained using the third natural vibration frequency value exhibits a strong R 2 value of 9.437 × 10 −1 , indicating that it explains a significant portion of the variance in the data.Even though defining the equation with more variables may decrease the actual prediction percentage, the results are within acceptable limits.The MSE value is determined to be 5.461, accompanied by a MAE value of 1.758.
The pareto charts, which are used to identify the most significant parameters for natural vibration frequencies, and 3D graphics are presented in Figure 13.As seen in Figure 13, the geometry of the pipeline, including diameter and wall thickness, affects its natural vibration frequencies.Thicker pipelines tend to have lower frequency values.The effect of pipe thickness on natural vibration frequency increases with larger pipe diameters.Developing explicit expressions that define the relationship between the inputs and the natural vibration frequencies of the pipe would be particularly advantageous in terms of time and computational load, especially in the presence of internal flow.
Table 6 summarizes the previously obtained results, where the models are sorted in descending order, in terms of performance criteria, depending on the train data set used.Table 6.Performance comparison of the prediction equations.

Performance Criteria
Prediction Performance As seen in Table 6, the suggested equations for the full quadratic model exhibit the best performance for the initial two natural vibration frequencies.However, the relationship for displacement estimation, while remaining within the validity bounds, demonstrates weaker performance in terms of the R 2 value.The natural vibration frequency values for the third mode yield numerically larger results compared to other outputs.Therefore, MSE and MAE values are correspondingly larger.

Conclusions
The primary objective of this investigation is to establish the correlation among the factors influencing submarine pipelines without an extensive and complex analysis process.In this study, the widely used Airy Wave Theory is employed.Hydrodynamic wave forces acting on submarine pipeline models, each characterized by 18 different geometric properties, are calculated using the Morison Equation.
Dynamic behaviors are investigated through Finite Element-based computer-aided Time History analyses, with consideration given to constant clearance and environmental factors for each pipeline.A dataset including maximum displacement, stress, and the initial three natural vibration frequencies is generated.Multiple Regression analysis is performed using 144 data pairs with 4 inputs and 5 outputs.The relationship between inputs and outputs is defined by mathematical equations using a Full Quadratic model.Generally, it is observed that variations in diameter had a greater effect on displacement values than on stress values.The largest variation between predicted and actual values for maximum displacement and stress values occurs under conditions where the diameter is 1000 mm at T/4 and 3T/4.This finding supports the notion that hydrodynamic forces are the predominant variables in stress and displacement prediction equations.Examination of natural vibration frequency values reveals that increasing the number of variables required to describe the relationship between inputs and outputs reduced prediction performance by 5.630%.Therefore, this study aims to obtain the simplest mathematical equation possible, yielding the best results.Considering the influence of submarine pipeline geometry on hydrodynamic force variations, this study is significant in elucidating and interpreting the relationship between maximum stress, displacement, and natural vibration frequency.Furthermore, performance criteria are determined for each equation.These results show that the proposed models provide a highly accurate depiction of the data and give a great deal of confidence in the accuracy and reliability of the model.The results are observed in the case of neglecting the internal flow of the pipeline to simplify the analysis.When the internal flow of the pipeline is taken into account, correlations that include the effects on the outputs, especially on the natural vibration frequency, will become more significant.Future research will incorporate the effects of internal and external fluid flow.From this perspective, this study provides valuable insights into the underlying dynamics of submarine pipeline systems and analyzes inflection points in the relationships between variables.

Figure 2 .
Figure 2. Surface elevations based on Airy Wave Theory.

Figure 4 .
Figure 4. Main steps of the system architecture.

Figure 6 .
Figure 6.Comparison of maximum values of stress.

Figure 7 .
Figure 7.Comparison of maximum values of displacement.

Figure 8 .
Figure 8.The maximum stress value versus the maximum displacement.It is observed that the relationship between maximum stress and displacement values is significant based on the Coefficient of Determination.Deviation in the values occurs when the hydrodynamic forces change from positive to negative.In this wave position, while stress increases, displacement values numerically decrease due to the change in direction.

Figure 9 .
Figure 9. Pareto chart of the standardized effects for wmax.

Figure 11 .
Figure 11.Pareto chart of the standardized effects for σmax.The regression model exhibits a strong R 2 value of 9.518 × 10 −1 , indicating that it explains a significant portion of the variance in the data.Moreover, MSE is low at 2.691 × 10 −3 and MAE is 3.903 × 10 −2 .The results of quality metrics showed that the model's predictions are generally close to the actual values.The relationship between the independent variables and maximum stress is given in Figure12.

Table 1 .
Properties of the pipe models.

Table 2 .
Comparison of natural vibration frequency values.

Table 3 .
Analysis of Variance for wmax values.

Table 4 .
Analysis of Variance for σmax values.