High-Accuracy and Fast Calculation Framework for Berthing Collision Force of Docks Based on Surrogate Models

: The accurate prediction of the collision force magnitude resulting from ship berthing on docks is crucial for the rationality and safety of dock structural design. This paper presents a novel framework for the calculation of berthing collision force for ships (CBCF), which integrates field data, finite element models, and surrogate models. Based on field data and finite element analysis, the framework constructs and compares four surrogate models with low sample requirements, ultimately selecting the optimal surrogate model for predicting collision force. Furthermore, a sensitivity analysis of the parameters is conducted based on the selected model, followed by a comparison with the various methods used for collision force prediction. The results illustrate the effectiveness of the proposed framework in replacing finite element models for the rapid and accurate prediction of collision force. Comparison with existing methods also underscores the advantages of the proposed framework, including low sample requirements, high calculation accuracy, and exceptional efficiency. In summary, this study not only introduces a novel and precise surrogate model framework for the swift prediction of berthing collision force, but it also offers valuable insights into the prevention of ship collision with wharf accidents and facilitates the rational and safe design of wharf structures.


Introduction
Due to its cost-effectiveness and substantial carrying capacity, maritime shipping has emerged as the predominant method of global cargo transport, playing a pivotal role in international trade and the broader development of the global economy [1,2].Yet, with the increasing volume of shipping traffic, a range of challenges has emerged.Established port facilities may fall short of accommodating the berthing needs of increasingly larger vessels, potentially raising the likelihood of accidents.Research [3] indicates that 43.6% of waterborne transportation accidents are attributed to vessel collisions.To mitigate such occurrences, existing standards such as the AASHTO [4], European norms [5], and the Chinese railway bridge and culvert specifications [6] offer formulas to compute the collision forces in ship-bridge collisions.An accurate calculation of the collision forces can inform the design of structures or the implementation of protective measures, thus ensuring safety.However, there is a lack of a direct and suitable method for calculating the collision forces when ships berth at docks.
In 1959, Minorsky [7] conducted a comprehensive examination of 26 nuclear-powered ship collision incidents, thus pioneering the field of ship collision research.Since then, the topic of collision force has been extensively explored through methods like empirical formula fitting, impact simulation experiments, and numerical simulations.For instance, Zhou et al. [8] analyzed ship-bridge collision forces as delineated in differing standards and devised a theoretical formula supported by finite element analysis, considering cumulative damage in their calculations.Song et al. [9] performed numerical simulations that resulted in 45 time-history curves of ship-bridge collision forces to assist in parameter identification within collision models and aid in the prediction of impact loads over time.Consolazio [10] proposed the recommended bow extrusion curve and revised the AASHTO formula based on an analysis of multiple high-resolution bow-crushing simulations.Sha et al. [11] evaluated the ship impact force and wharf response through indoor impact tests and numerical simulation analysis, proposing an empirical formula for calculating the peak impact force related to barge mass and speed according to numerical calculation results.It is evident that the current research primarily focuses on ship-bridge collision issues.However, there has been limited analysis, discussion, or study on the significance of port wharves in water transport networks.Port wharves play a crucial role in undertaking berthing tasks; yet, this aspect has not received adequate attention.
In contrast to bridges, wharves are commonly equipped with rubber fenders to mitigate the impact response of ships during berthing [12,13].The current specifications do not provide a direct formula for calculating the collision force of vessels berthing at the wharf; instead, they derive the collision force manually by considering the berthing energy and combining it with the energy absorption curve of rubber fenders [14 -16].Considering that different types of fenders have varying characteristics in terms of energy absorption curves when installed at different wharves, the specification offers a more general calculation method to accommodate various specific situations.However, this method introduces a certain level of uncertainty in determining the berthing collision force, which may lead to an excessively conservative structural design for the wharf and subsequently escalate the construction costs.Therefore, accurately calculating the collision force of vessels berthing at the wharf holds significant practical value.
Among the various calculation methods available, high-resolution finite element models offer superior accuracy by considering material details, boundary conditions, and other relevant factors [17].In comparison to the limited availability of on-site data, a validated finite element model can effectively complement deficiencies in data.However, despite their ability to compensate for data shortages, performing high-precision finite element analyses not only necessitates significant computational resources [18] but also demands operators with a certain level of expertise.This scenario poses significant challenges to advancing more in-depth analyses.Thus, finding a computation method that can rapidly and accurately predict berthing collision forces is of crucial importance.When confronted with this dilemma, the emerging technology of machine learning (ML) demonstrates its superiority.
Unlike the traditionally stringent finite element simulations, surrogate models based on machine learning can perform rapid calculations in place of finite element models while ensuring accuracy, thereby significantly improving computational efficiency [19].Common types of surrogate models include artificial neural networks (ANN), polynomial chaos expansion (PCE), and response surface methodology (RSM).Surrogate models based on machine learning, with their excellent performance, have been widely applied in various fields.For example, Fan et al. [20] developed an artificial neural network-based method for calculating ship-bridge collision forces using numerical simulation data.Qin et al. [21] argued that the rise of machine learning has greatly promoted the development of the oceanography field and has used Physics-informed neural networks to address storm surge issues.Kameshwar [22] predicted the collision response of bridge piles based on the response surface methodology.In the field of engineering, using surrogate models to replace traditional finite element simulations for rapid analysis has become a new research trend.Considering the efficiency and reliability of surrogate models, applying them to calculate the collision force of ship berthing is a reasonable and feasible choice.
Addressing the existing research gaps in this particular domain, this paper proposes a novel framework for calculating the berthing collision force for ships (CBCF) by integrating finite element analysis and surrogate model technology with field data.The proposed framework can be widely applied to various wharves on-site, providing an exclusive and efficient calculation method for berthing collision force at each individual wharf.In Section 2, the CBCF framework is established in the order of "Construction-Optimization-Application".In Section 3, the effectiveness of the proposed framework is demonstrated through its application to a wharf in the Jiangsu province, China, where it accurately predicts and analyzes berthing collision forces as verified by physical model tests.Finally, Section 4 summarizes the key findings obtained from this study.

Construction of the CBCF Framework Based on Surrogate Models
The present section presents the specific composition of the proposed framework for calculating ship berthing collision force based on surrogate models.It encompasses the construction, evaluation, and optimization of the collision surrogate model, as well as the utilization of the optimized model to predict berthing collision force.

Construction of Collision Surrogate Model
According to ship-bridge collision codes [4,5], the primary factors influencing collision forces are typically the ship tonnage and the berthing speed.However, this study expands upon these factors by also considering the berthing angle and the rubber fender as additional influential parameters.It is worth mentioning that the berthing angle here is the angle between the berthing speed in the direction of ship navigation and the perpendicular line of the dock.The energy absorption capacity of rubber fenders is primarily determined by the stress-strain relationship of the material [23,24].Given that a ship may encounter multiple rubber fenders during berthing, it is assumed that each individual fender undergoes elastic deformation when compressed.Therefore, for simplicity in describing the energy absorption effect of rubber fenders in this study, we have selected the elastic modulus of rubber.In summary, this study utilizes ship tonnage t, berthing speed v, the rubber fender elastic modulus E, and the berthing angle θ as the input parameters for the collision surrogate model; the berthing collision force F serves as the output parameter; and a corresponding "Input-Output" dataset is constructed.The output parameter Y = {F} is also a random variable due to the inherent randomness of the input parameter vector t, v, E, θ, which follows a probability density function z(t, v, E, θ).The sample set is designated as (t, v, E, θ) j , Y j , j = 1, • • • , N. Considering the randomness of berthing modes, and integrating the development of computational models in recent years with the cost of sample calculations, this paper establishes four low-cost operational collision surrogate models for predicting and analyzing the collision force during ship berthing.
Col-Kriging Kriging was initially introduced by Ord [25] as an interpolation technique for stochastic spatial processes in the Statistical Encyclopedia.The collision Kriging (Col-Kriging) model developed in this study is designated as M K : where λ i , i = 0, • • • , 14 are the coefficient of the trend term, and σ 2 and Z(t, v, E, θ) denote the variance and the zero mean, the unit variance, and the stationary Gaussian process, respectively.

Col-PCE
The polynomial chaos expansion (PCE) principle aims to construct a spectral representation using appropriate polynomial functions to effectively approximate the output of a computational model, such as a finite element model [26].The collision PCE(Col-PCE) formula in this study is as follows: where λ α is the undetermined coefficient to be expanded; ψ α is the polynomial orthogonal to the joint probability density function of the input parameter t, v, E, θ; the index subscript α represents the basis function; and A 4,q denotes the truncated set of quaternion polynomials with a q-norm multi-index, encompassing any combination of t, v, E, θ.This constraint effectively limits the size of the polynomial by restricting the number of nonzero coefficients it can have. 3

Col-PCK
The PCK method combines the merits of Kriging and polynomial chaotic expansion, enabling the capture of both global behavior and local changes in computational models.The surrogate model established by Zhao et al. [27] using PCK exhibits the high flexibility and the strong nonlinear modeling ability, making it an optimal choice.Specifically, PCK consists of a general-purpose Kriging model whose trend part is a sparse set of orthogonal polynomials: where ∑ α∈A 4,q λ α ψ α (t, v, E, θ) is the weighted sum of orthogonal polynomials describing the trend of the PCK model; A 4,q corresponds to the same concept as in the PCE model; and σ 2 and Z(t, v, E, θ) share a similar interpretation as that of the Kriging model.Thus, PCK can be regarded as a specialized version of the general Kriging model, with distinct trend characteristics.
4 Col-SVR Support vector machines (SVM) originated from Cortes's research [28] and were initially employed for the binary classification and pattern recognition of test cases.This principle provides a significant generalization ability for support vector machine regression and effectively reduces the risk of overfitting.The formula of the Col-SVR model can be written as follows: where ω represents the weight coefficient vector, Φ t, v, E, θ; is the distance Gaussian kernel function, and b is the estimated offset parameter.

Evaluation and Optimization of Surrogate Models
Existing research indicates that the accuracy of surrogate models varies in different application scenarios.The aforementioned surrogate model has achieved good results in various fields [29][30][31], but whether it can maintain the same accuracy when applied to the prediction of dock ship berthing collision forces still requires further verification.The "Input-Output" dataset is split into two parts: the training set is used to build different surrogate models, and the validation set is used to evaluate the training results.Then, the most suitable surrogate model is selected to replace the on-site dock or finite element analysis to predict the dock berthing collision force.
Based on the prediction results, surrogate models are evaluated to select the optimal one for predicting the berthing collision force, which will be used for the construction of the CBCF framework.Numerous accuracy evaluation indexes are available, and, based on the studies conducted by Shao [32] and Hu [33], validation error (val) and leave-one-out crossvalidation (LOO) are selected as the evaluation metrics for their ability to effectively reflect the accuracy of surrogate models.The validation set, denoted as {X val , Y val = M(X val )}, is utilized to calculate the validation error E val according to the following formula: where μY val = 1 val ) represents the mean value of the validation set, and K denotes the total number of samples in the validation set.
The leave-one-out (LOO) error exhibits a relatively higher robustness.It works on the principle of isolating a point from the complete experimental data and constructing a new surrogate model from the remaining points.This process is repeated for every point in the design of experiments to evaluate the accuracy and applicability of the constructed model.The formula for calculating the error of E Loo is as follows: where M •\i (•) is the surrogate model obtained after removing the ith sample point from the complete experimental data, and μY is the mean value of the samples in the design of experiments.

CBCF Framework Process
By employing the surrogate model determined through evaluation in the preceding section, the computation of the berthing collision force for ships can be performed with an acceptable accuracy across various operational scenarios.The CBCF framework in this study is constructed following the "Construction-Optimization-Application" sequence, as illustrated in Figure 1.During the construction phase, a diverse "Input-Output" dataset composition is obtained primarily through two channels: field data and finite element analysis [34].The optimization phase involves training and evaluating with training and validation sets to identify the optimal collision surrogate model, taking into account both computational cost and prediction accuracy.Finally, during the application phase, the most optimal surrogate model is employed to predict berthing collision forces under different operational conditions at specific wharves.

Framework Application
In this section, the proposed CBCF framework is implemented at a wharf located in the Jiangsu province, China, and its reliability is validated through physical model testing.The "Input-Output" dataset is obtained by combining finite element analysis with field data verification.

Project Introduction
The wharf segment measured 68 m in length and 36 m in width, featuring a total of 10 beams spaced at intervals of 7 m.It was divided into upper and lower parts, with the upper part consisting of cast-in-place slabs and prefabricated longitudinal and transverse beam structures made from C40 concrete material.The lower part comprised 94 PHC tubular piles with a diameter of 1 m, utilizing C80-grade concrete.Rubber fenders were positioned every 14 m along the front edge of the wharf.The main design vessel dimensions for the high-pile wharf berth are provided in Table 1.

The Establishment of the Finite Element Model 3.2.1. Material Parameter Setting
Based on the ABAQUS (2021) finite element software, the ship-rubber fender-wharf numerical model was constructed, as shown in Figure 2. Considering that the primary research focus was on the collision force during ship berthing, simplifications were made to the ship model: its geometric dimensions remained constant, while variations in material density and the speed parameters were employed to simulate the berthing process for ships of different tonnages at various normal berthing speeds and angles.The upper structure was represented by C3D8R elements corresponding to C40 concrete material, whereas S4R elements were utilized to model the piling foundation with respect to C80 concrete material.

Material Parameter Setting
Based on the ABAQUS (2021) finite element software, the ship-rubber numerical model was constructed, as shown in Figure 2. Considering tha research focus was on the collision force during ship berthing, simplificatio to the ship model: its geometric dimensions remained constant, while material density and the speed parameters were employed to simulate process for ships of different tonnages at various normal berthing speeds an upper structure was represented by C3D8R elements corresponding to material, whereas S4R elements were utilized to model the piling foundatio to C80 concrete material.Currently, several commonly used superelastic constitutive mode fenders include the Mooney-Rivlin model, the Arriuda-Boyce model, the O the Yeoh model, etc. [35].Among these models, the strain energy function o Boyce model is excessively complex, and its parameters are challenging accurately; thus, it is unsuitable for finite element simulation.The determina parameters for the strain energy function in the Ogden model poses difficu experimental means.Additionally, the strain energy function of the Yeoh m significant errors at small deformations and fails to simulate small compression effectively.In comparison with these models, the Mooneyoffers a simpler strain energy function that allows the calculation of releva using the elastic modulus [36], making it widely employed in simula constitutive relationship.Therefore, this study utilized the Mooney-Riv simulate such a relationship, as it considered strain energy density as function of strain invariant and represented an incompressible material w term for the strain invariant-referred to as Mooney-Rivlin material [37].Th relationship of the biparametric Mooney-Rivlin model, as provided in expressed as follows: where W represents the strain energy density, 1 I and 2 I denote the fir Currently, several commonly used superelastic constitutive models for rubber fenders include the Mooney-Rivlin model, the Arriuda-Boyce model, the Ogden model, the Yeoh model, etc. [35].Among these models, the strain energy function of the Arriuda-Boyce model is excessively complex, and its parameters are challenging to determine accurately; thus, it is unsuitable for finite element simulation.The determination of related parameters for the strain energy function in the Ogden model poses difficulties through experimental means.Additionally, the strain energy function of the Yeoh model exhibits significant errors at small deformations and fails to simulate small deformation compression effectively.In comparison with these models, the Mooney-Rivlin model offers a simpler strain energy function that allows the calculation of relevant parameters using the elastic modulus [36], making it widely employed in simulating rubber's constitutive relationship.Therefore, this study utilized the Mooney-Rivlin model to simulate such a relationship, as it considered strain energy density as a polynomial function of strain invariant and represented an incompressible material with only one term for the strain invariant-referred to as Mooney-Rivlin material [37].The constitutive relationship of the biparametric Mooney-Rivlin model, as provided in ABAQUS, is expressed as follows: where W represents the strain energy density, I 1 and I 2 denote the first and second strain invariants, respectively, C 10 and C 01 are the Rivlin coefficients, E = 7.5C 10 = 30C 01 , J stands for the volume compression ratio, and D represents the incompressibility coefficient of the volume.
If the material is completely incompressible (such as rubber, Poisson's ratio 0.4997) [36], the strain energy density formula can be changed into the following: The material characteristics of each structural component in the finite element model utilized in this study are presented in Table 2.

Model Validity Verification
In order to validate the accuracy of the established finite element model in reflecting real-world conditions, this case employed three approaches to verify the efficacy of the model.

Berthing Collision Force Verification
According to the whole-life-cycle monitoring project of the wharf conducted by the China Communications Construction Company, berthing collision force sensors were installed in the rubber fender of the wharf.Representative berthing record samples from the wharf operation logs were selected to verify the material parameters of the finite element model.Taking examples with a berthing tonnage of 2100 t and a berthing speed of 0.29 m/s, as well as a berthing tonnage of 14,500 t and a berthing speed of 0.14 m/s, these cases involved ships approaching in an approximately parallel direction to the surface of the wharf, as illustrated in Figure 3.The berthing collision forces measured by the on-site monitoring system were 1792 kN and 1961 kN, respectively.The simulations carried out in the ABAQUS software indicated that the time-history data of the collision forces shown in Figure 4 correspond to the actual scenario.The finite element model effectively replicated the fluctuation of collision force during the berthing process, consistent with the data obtained from on-site monitoring.At the time of berthing of 0.43 s and 0.7 s, the collision force reached peak values of 1768 kN and 2013 kN, respectively, with errors of only 1.34% and 2.7%, indicating that the berthing collision force of the finite element model is in accordance with the actual situation.

Berthing Collision Energy Verification
Energy analysis serves as a crucial approach to ensure the accuracy of cal [38].The berthing process of a ship can be regarded as an energy conversion where the initial kinetic energy of the ship is transformed into absorbed interna by both the rubber fender and the wharf structure, along with artificial strain ene energy change diagram of the berthing collision system obtained in the finite analysis software ABAQUS is presented in Figure 5, depicting a 5000-ton ship be a normal velocity of 0.3 m/s.In the initial stage of the collision, there is a rapid de kinetic energy, followed by a gradual reduction in its rate of decline.Over time, for the internal energy exhibits an opposite behavior to that of kinetic energy, w reaching their peaks around 1.1 s.Additionally, due to the hourglass effect [39], strain energy inevitably increases slightly over time, leading to a marginal rise in system energy.To ensure that the simulation results remain unaffected by the h problem, it is generally necessary to control the proportion of artificial strai within 10% [38,40].In this study's berthing model, we successfully contro proportion of artificial strain energy within 5%, thus demonstrating the feasibili

Berthing Collision Energy Verification
Energy analysis serves as a crucial approach to ensure the accuracy of calculations [38].The berthing process of a ship can be regarded as an energy conversion process, where the initial kinetic energy of the ship is transformed into absorbed internal energy by both the rubber fender and the wharf structure, along with artificial strain energy.The energy change diagram of the berthing collision system obtained in the finite element analysis software ABAQUS is presented in Figure 5, depicting a 5000-ton ship berthing at a normal velocity of 0.3 m/s.In the initial stage of the collision, there is a rapid decrease in kinetic energy, followed by a gradual reduction in its rate of decline.Over time, the trend for the internal energy exhibits an opposite behavior to that of kinetic energy, with both reaching their peaks around 1.1 s.Additionally, due to the hourglass effect [39], artificial strain energy inevitably increases slightly over time, leading to a marginal rise in the total system energy.To ensure that the simulation results remain unaffected by the hourglass problem, it is generally necessary to control the proportion of artificial strain energy within 10% [38,40].In this study's berthing model, we successfully controlled the proportion of artificial strain energy within 5%, thus demonstrating the feasibility of our finite element simulation results.

Self-Vibration Frequency Verification
The natural frequency of a structure is an inherent characteristic th linked to the structure itself, and it serves as a crucial parameter fo accuracy of finite element models [41,42].The finite element model was study to analyze the vibration mode and obtain information on the frequ modes.Regarding the vibration data monitored on-site, the NexT-ER employed in this study for modal identification, which integrates th realization algorithm with the natural excitation technique [43].The com are compiled and presented in Table 3.The modal frequencies obtained analysis conducted in this article show a maximum error of only 1.79% to the measured results obtained on-site, indicating that they fall with margin of error.Based on the analysis results from the aforementioned three aspects, that the finite element model constructed in this study accurately depict response characteristics of the pier at the berthing site and effective mechanical behavior during berthing collisions.Consequently, this m dependable computational platform for investigating ship berthing proce simulation reliability.

Self-Vibration Frequency Verification
The natural frequency of a structure is an inherent characteristic that is intrinsically linked to the structure itself, and it serves as a crucial parameter for validating the accuracy of finite element models [41,42].The finite element model was employed in this study to analyze the vibration mode and obtain information on the frequency of vibration modes.Regarding the vibration data monitored on-site, the NexT-ERA method was employed in this study for modal identification, which integrates the Eigen-system realization algorithm with the natural excitation technique [43].The comparative results are compiled and presented in Table 3.The modal frequencies obtained from the modal analysis conducted in this article show a maximum error of only 1.79% when compared to the measured results obtained on-site, indicating that they fall within an acceptable margin of error.Based on the analysis results from the aforementioned three aspects, it can be inferred that the finite element model constructed in this study accurately depicts the mechanical response characteristics of the pier at the berthing site and effectively simulates its mechanical behavior during berthing collisions.Consequently, this model provides a dependable computational platform for investigating ship berthing processes and ensures simulation reliability.

The Acquisition of the "Input-Output" Dataset
By utilizing the finite element model established in the preceding chapter, a series of ship berthing simulations were conducted to compensate for the absence of field data, thereby acquiring the "Input-Output" dataset essential for the surrogate model.Considering the diverse working conditions at the dock, this study defined a range of ship masses from 2000 to 20,000 tons.The berthing velocities varied between 0.1 and 0.8 m per second, while the angles of berthing ranged from 0 to 15 degrees.Simultaneously, the elastic modulus for the rubber fenders was maintained within a range of 3 to 8 MPa.These crucial input-output parameters are presented in Table 4. Employing a finite element model, this research conducted simulation analyses for a total of sixty berthing scenarios, encompassing various combinations of the four key input parameters mentioned above.These simulations covered an extensive array of random permutations of the input parameters and, thus, provided substantial data support for subsequent research endeavors.

Comparative Analysis of Surrogate Models
In principle, augmenting the surrogate model with additional random samples from the training set can enhance its accuracy.However, excessively increasing the number of sampled data points is not an optimal strategy in terms of computational cost-efficiency.The surrogate model can be considered sufficiently accurate to replace the original computational model when its error reaches 10 −3 [44].When constructing surrogate models, it is crucial to achieve an optimal balance between accuracy and computational cost.In this study, four surrogate models were constructed and trained using training sets of varying sizes N tra (i.e., 10,20,30,40), all under the same validation set N val = 20.The accuracy discrepancies among the different surrogate models were discussed in order to explore the best-case scenario that minimized computational effort while maximizing precision.The detailed changes in the collision force output parameters E νal and E Loo changed in detail, as shown in Figures 6 and 7, which indicated the following: (1) The four surrogate models for E val exhibited a significant reduction in error as N tra increased.For N tra ≤ 20, both the Col-PCK and Col-SVR models demonstrated a substantial decrease in errors with increasing N tra , while, for N tra > 20, the errors of these two models remained relatively stable.Similarly, Col-PCE followed a similar trend around N tra = 30.On the other hand, the error of Col-Kriging displayed considerable fluctuations.Considering both computational cost and prediction accuracy, the Col-PCK model demonstrated optimal performance in predicting berthing collision forces within the current dataset.The Col-PCK model was preferentially selected under conditions of N tra = 20 and N val = 20 to constitute the core of the CBCF (PCK-CBCF).The predictive performance of PCK-CBCF was compared with the actual results obtained from finite element analysis (FEM).As depicted in Figure 8, the alignment between the distribution of PCK-CBCF predictions and FEM outcomes along the line of symmetry further validated the exceptional precision exhibited by the Col-PCK model.
Considering both computational cost and prediction accuracy, th demonstrated optimal performance in predicting berthing collision current dataset.The Col-PCK model was preferentially selected un to constitute the core of the CBCF (PCK-CBC performance of PCK-CBCF was compared with the actual results obt element analysis (FEM).As depicted in Figure 8, the alignment betwee of PCK-CBCF predictions and FEM outcomes along the line of symmetry the exceptional precision exhibited by the Col-PCK model.The selected optimal surrogate model now has the capability to e the finite element model for berthing collision force calculations.By parameter combinations, the rapid acquisition of collision force facilitated.In comparison, performing finite element simulation conditions on a computer equipped with an Intel Core i7-11800H CPU whereas utilizing PCK-CBCF necessitates only a total of 2.8 s, including the training set.Consequently, this approach achieves a signific computational resources and time costs.

Framework Verification through Physical Model Testing
The practical applicability of PCK-CBCF was validated through physical model test data conducted in a large harbor tank.While suitab were made to replicate dock structures, such as the length of piles, physical model to accurately reproduce the actual conditions of rubb impact force collection areas ensured an ideal validation environment f model specifically designed for this particular dock in the Jiangsu prov model settings can be found in Zhai [45], and Figure 9 depicts the simul site.A calibrated fender sensor was utilized to replicate the rubber fend modulus of 3.5 MPa, and the vessel was docked in a parallel orientation, impact force data and corresponding comparative outcomes, as present The selected optimal surrogate model now has the capability to effectively replace the finite element model for berthing collision force calculations.By inputting various parameter combinations, the rapid acquisition of collision force output results is facilitated.In comparison, performing finite element simulation under identical conditions on a computer equipped with an Intel Core i7-11800H CPU requires 35 min, whereas utilizing PCK-CBCF necessitates only a total of 2.8 s, including training time for the training set.Consequently, this approach achieves a significant reduction in computational resources and time costs.

Framework Verification through Physical Model Testing
The practical applicability of PCK-CBCF was validated through the utilization of physical model test data conducted in a large harbor tank.While suitable simplifications were made to replicate dock structures, such as the length of piles, the ability of the physical model to accurately reproduce the actual conditions of rubber fenders at key impact force collection areas ensured an ideal validation environment for the PCK-CBCF model specifically designed for this particular dock in the Jiangsu province.The specific model settings can be found in Zhai [45], and Figure 9 depicts the simulated berthing test site.A calibrated fender sensor was utilized to replicate the rubber fender with an elastic modulus of 3.5 MPa, and the vessel was docked in a parallel orientation, providing precise impact force data and corresponding comparative outcomes, as presented in Table 5.The results demonstrated that the error between the collision force calculated by the PCK-CBCF model and the experimental data was below 6.9% across the various berthing conditions, thereby effectively validating the reliability and high accuracy of the proposed calculation method in practical applications.

Sensitivity Analysis
Global sensitivity analysis (GSA) is an exceptional scientific methodology that not only elucidates the potential impact of parameter variations on research outcomes but also enables researchers to gain a more precise understanding of model robustness and predictive capabilities.The Sobol method is widely favored within the research community due to its provision of high-precision information [46].However, the calculation of the Sobol index often relies on Monte Carlo simulations with large amounts of data, which creates a bottleneck for finite element models with extensive calculations.Fortunately, leveraging the fast computation characteristics of surrogate models enables the efficient evaluation of the Sobol index without incurring significant computational costs [47].In this study, its functionality is expressed as follows: The PCK-CBCF method proposed in this study offers a practical and efficient approach for conducting sensitivity analysis on the influential factors of berthing collision force, without imposing any additional burden.By employing PCK-CBCF, a comprehensive sensitivity analysis is conducted in this paper to investigate the four key variables that affect the berthing collision force: ship tonnage, berthing speed, berthing angle, and elastic modulus of the rubber fender.The results are presented in Figure 10.The first-order Sobol index quantifies the direct impact of an individual i variable on the model output, disregarding its interactions with other input variable contrast, the full-order Sobol index considers all the potential interactions among input variables, revealing their collective influence on the output results.From analysis of the results, it is evident that the berthing speed exhibits the highest Sobol in indicating its significant impact on the berthing collision force and emphasizing importance of considering its interaction effects with other parameters.Conversely, o input variables such as the berthing angle have a minimal influence on the collision f Furthermore, a disparity between the first-order and full-order Sobol indexes for tonnage suggests that interactions involving ship tonnage and other variables may some influence on the collision force.
Due to the significant influence of the speed at which a vessel docks on the for collision, especially in the case of larger ships, it is essential to enforce strict regula on their velocity during berthing in order to reduce the risk of causing damage to structures.Although the impact of the berthing angle is relatively minor, exce berthing angles can result in significant ship rolling and compromise the stability o berthing process.Therefore, it is imperative to minimize the berthing angle du operations to ensure a stable and secure ship berth.

Comparison with Existing Methodologies
In accordance with the computational methodology suggested in this researc estimating collision force during ship berthing using surrogate modeling, adva machine learning methods like neural networks have been utilized to supp conventional finite element analysis for the prediction of ship-bridge collision fo enabling rapid and efficient calculations.For instance, Fan [20] demonstrated applicability of RBFNN in ship-bridge collision force prediction by leveragin nonlinear mapping capabilities and training/validation using finite element model Similarly, Xu [17] utilized a BPNN model along with FFT and IFFT methods to pr ship-bridge collision forces in both time and frequency domains.In this section performance of the proposed PCK-CBCF model will be compared with that of the exi RBFNN and BPNN models in terms of predicting the collision force between ships wharves.To accomplish this, 220 datasets representing random working conditions been generated by combining finite element analysis and the Monte Carlo samp method.Model training and performance evaluation will now be conducted base The first-order Sobol index quantifies the direct impact of an individual input variable on the model output, disregarding its interactions with other input variables.In contrast, the full-order Sobol index considers all the potential interactions among the input variables, revealing their collective influence on the output results.From the analysis of the results, it is evident that the berthing speed exhibits the highest Sobol index, indicating its significant impact on the berthing collision force and emphasizing the importance of considering its interaction effects with other parameters.Conversely, other input variables such as the berthing angle have a minimal influence on the collision force.Furthermore, a disparity between the first-order and full-order Sobol indexes for ship tonnage suggests that interactions involving ship tonnage and other variables may exert some influence on the collision force.
Due to the significant influence of the speed at which a vessel docks on the force of collision, especially in the case of larger ships, it is essential to enforce strict regulations on their velocity during berthing in order to reduce the risk of causing damage to dock structures.Although the impact of the berthing angle is relatively minor, excessive berthing angles can result in significant ship rolling and compromise the stability of the berthing process.Therefore, it is imperative to minimize the berthing angle during operations to ensure a stable and secure ship berth.

Comparison with Existing Methodologies
In accordance with the computational methodology suggested in this research for estimating collision force during ship berthing using surrogate modeling, advanced machine learning methods like neural networks have been utilized to supplant conventional finite element analysis for the prediction of ship-bridge collision forces, enabling rapid and efficient calculations.For instance, Fan [20] demonstrated the applicability of RBFNN in ship-bridge collision force prediction by leveraging its nonlinear mapping capabilities and training/validation using finite element model data.Similarly, Xu [17] utilized a BPNN model along with FFT and IFFT methods to predict ship-bridge collision forces in both time and frequency domains.In this section, the performance of the proposed PCK-CBCF model will be compared with that of the existing RBFNN and BPNN models in terms of predicting the collision force between ships and wharves.To accomplish this, 220 datasets representing random working conditions have been generated by combining finite element analysis and the Monte Carlo sampling method.Model training and performance evaluation will now be conducted based on these datasets.
The accuracy of three models in predicting berthing collision force under different training sets (N val = 20) is clearly demonstrated in Figure 11.The prediction accuracy of the model is quantified by the mean absolute percentage error (MAPE), which is calculated as follows: where y m represents the FEM-calculated value, y m represents the model-predicted value, and M is the number of predicted working conditions.
the model is quantified by the mean absolute percentage error (MAPE), which is calculated as follows:   6.The berthing conditions are arranged in the following order: ship tonnage, berthing speed, elastic modulus, and berthing angle ( , , , ) t v E θ .From the data in the table, it is evident that, when trained on the same dataset, the PCK-CBCF model proposed in this study not only exhibits a superior prediction accuracy but also demonstrates a higher computational efficiency for calculating impact forces under random conditions.The figure demonstrates that the prediction accuracy of the three models increases as the number of training sets N tra increases.Notably, when the training samples are limited (N tra = 10), PCK-CBCF outperforms the other two models in terms of prediction accuracy.Conversely, the BPNN-CBCF and RBFNN-CBCF models require larger training sets (N tra = 100 and N tra = 200) to achieve better prediction results.This comparison emphasizes the adaptability of the PCK-CBCF model to small-scale training sets, particularly when data constraints exist, as it maintains both efficiency and accuracy.
The three models have been selected to predict the collision force under specific berthing conditions, with N tra = 200 and N val = 20.The predicted results are presented in Table 6.The berthing conditions are arranged in the following order: ship tonnage, berthing speed, elastic modulus, and berthing angle (t, v, E, θ).From the data in the table, it is evident that, when trained on the same dataset, the PCK-CBCF model proposed in this study not only exhibits a superior prediction accuracy but also demonstrates a higher computational efficiency for calculating impact forces under random conditions.

Conclusions
Considering the complications and low accuracy plaguing current approaches to collision force calculation during ship berthing in wharves, this study presents a novel calculation framework aimed at accurately and quickly predicting the collision force during ship berthing.Furthermore, it extensively examines the suitability, effectiveness, and sample prerequisites of various surrogate models in these analyses and identifies the most appropriate model for substituting conventional finite element techniques in expedited computations.The practicality of this framework is verified through physical model tests, while a sensitivity analysis is conducted to examine the various factors influencing the collision force.Comparison with two existing models further confirms the advantages of this approach.The following conclusions can be drawn: 1.
The surrogate model based on the "Input-Output" dataset demonstrates effective substitution for the finite element model, enabling the rapid prediction of the berthing collision force.Notably, among these models, the Col-PCK model exhibits performance in predicting the berthing collision force.

2.
The sensitivity analysis conducted using the CBCF framework proposed in this study reveals that berthing speed is the most influential factor, followed by ship tonnage.This suggests that a greater emphasis should be placed on controlling berthing speed during ship berthing processes.

3.
The PCK-CBCF model demonstrates significant advantages over existing berthing collision force prediction methods in terms of reducing sample requirements, enhancing prediction accuracy, and improving computational efficiency when applied to a case study of a wharf in the Jiangsu province.This finding underscores the extensive potential of the proposed framework for future applications and its role in promoting safe wharf operations.
However, it is worth noting that current research primarily relies on validated finite element simulations rather than the direct utilization of field-measured data to enhance the frame structure.With the future acquisition and accumulation of field data on wharf impact, it is anticipated that further optimization in both the structure and performance of the CBCF frame will enable a better adaptation to collision force estimation and risk analysis for diverse wharves.

Figure 2 .
Figure 2. Finite element model for berthing analysis.

Figure 2 .
Figure 2. Finite element model for berthing analysis.

Figure 3 .
Figure 3. Picture of the berthing scene.

Figure 3 .
Figure 3. Picture of the berthing scene.

Figure 4 .
Figure 4. Comparison between FEM and measured collision force.

Figure 4 .
Figure 4. Comparison between FEM and measured collision force.

Figure 5 .
Figure 5. Energy change in the berthing collision system.

Figure 5 .
Figure 5. Energy change in the berthing collision system.

Figure 9 .
Figure 9. Pictures of the berthing test site.

y 10 traN
represents the FEM-calculated value, m y  represents the model-predicted value, and M is the number of predicted working conditions.The figure demonstrates that the prediction accuracy of the three models increases as the number of training sets tra N increases.Notably, when the training samples are limited ( = ), PCK-CBCF outperforms the other two models in terms of prediction accuracy.Conversely, the BPNN-CBCF and RBFNN-CBCF models require larger training sets ( better prediction results.This comparison emphasizes the adaptability of the PCK-CBCF model to small-scale training sets, particularly when data constraints exist, as it maintains both efficiency and accuracy.

Figure 11 .
Figure 11.Comparison of models' accuracy across diverse training datasets.

Table 1 .
Design ship type main dimension parameter table.

Table 2 .
Material characteristics of the primary structures in finite element models.

Table 3 .
Comparison of natural vibration frequencies of the wharf.

Table 3 .
Comparison of natural vibration frequencies of the wharf.

Table 4 .
Key parameters of the finite element.

Table 5 .
Comparison between physical model and PCK-CBCF.

Table 6 .
Comparison of FEM calculation results and model prediction results.

Table 6 .
Comparison of FEM calculation results and model prediction results.