A Stereolithographic Model-Based Dense Body Plan Generation Method to Construct a Ship Hydrodynamic Coefficients Database

A ship’s body plan is a vital data resource of ship hydrodynamics analysis, especially for time-domain simulations. Motivated by 3D printing technology, a novel dense body plan generation method is developed in this study. The slicing algorithm is adopted to generate dense 2D body plans from ship stereolithographic models. The dense body plan can be produced automatically under arbitrary rotational angles and slices. Moreover, a section redistribution algorithm is integrated to eliminate the non-uniform distribution features in sliced data inherited from the stereolithographic models. The benchmark ship models are selected to validate the accuracy of the method. The hull volumes of three ship models are calculated based on the produced data. The calculated results show satisfactory agreement with the published values. Furthermore, the estimation formulas of wetted surface area (WSA) are reviewed and utilized for validation. The calculated WSAs by slice integration turn out to be adaptive and accurate. The time costs of different slices are provided to illustrate the computational efficiency. A ship hydrodynamic coefficients database is constructed based on a 2D strip method and the produced data. The proposed method aims to improve the generation process of the body plan, which could meet the accuracy requirements of the strip method. As a result, hydrodynamic coefficients utilized in time-domain simulations could be obtained smoothly from the database.


Introduction
The maneuvering behaviors of a ship in an actual ocean environment are complicated. Large roll and pitch angles may occur in adverse weather conditions [1]. Ships may even sink in extremely rough seas [2]. However, simulation-based approaches provide a flexible and resource-saving way to predict and analyze ship performance in the ocean environment [3]. Therefore, numerical simulations are applied to enhance the safety of navigation in the maritime sector [4]. In addition, the simulation of ship motion in waves is a key component of marine simulators [5]. Precise motion prediction would ensure the physical reality of the simulators. Moreover, fast simulations of trajectory and motion status allow deck officers to judge whether the maneuvering actions have the correct tendency or indicate the need for corrections [6]. The simulation of ship motions could also contribute to the development of the model's predictive control algorithm [7]. Since the wave disturbance is the most dangerous effect in navigation, many investigations have focused on simulating ship maneuvering behaviors in waves [8][9][10][11]. The pioneering work in this field was undertaken by Cummins [12] and Ogilvie [13]. The impulse response function (IRF) is introduced and developed in their work, and the memory effect The most-commonly used 3D printing process (46% as of 2018) is a material extrusion technique called fused deposition modeling (FDM) [35]. The FDM process builds a three-dimensional object from a CAD model by successively adding material layer by layer. The flowchart of a general FDM process is shown in Figure 1. The initial part of a 3D printing pipeline is to create digital prototypes using various CAD software. After establishing the digital 3D model, the integrated CAM software is utilized to generate the G-code file. Technically, G-code is a machine language used to drive the motion of motors. The main task of the CAM software is to slice the digital 3D model into layers and generate a G-code subsequently. The G-code will be transferred to the control module of a 3D printer after slicing. The 3D physical model of the KVLCC2 printed by the FDM process is shown in Figure 2.

Stereolithographic Model
The prerequisite for producing any 3D printed part is a digital 3D model or a CAD file. The most widely used format of a digital 3D model in 3D printing is stereolithographic format due to its compact structure. The triangular facet is adopted to describe closed complex geometries in a stereolithographic model. The two format types of the model, the ASCII, and the binary are illustrated in Table 1. An ASCII stereolithography (STL) file begins with a solid name and concludes with the end solid name. It contains any number of triangles, where each n or v is a floating-point number. Since the ASCII STL files can become very large, a binary version of STL exists. A binary STL file has an 80-character header, where a four-byte little-endian unsigned integer indicating the number of triangular facets is present. Each triangle is described by twelve 32-bit floating-point numbers: three for the normal and then three for the X/Y/Z coordinate of each vertex. The file simply ends after the last triangle. The initial part of a 3D printing pipeline is to create digital prototypes using various CAD software. After establishing the digital 3D model, the integrated CAM software is utilized to generate the G-code file. Technically, G-code is a machine language used to drive the motion of motors. The main task of the CAM software is to slice the digital 3D model into layers and generate a G-code subsequently. The G-code will be transferred to the control module of a 3D printer after slicing. The 3D physical model of the KVLCC2 printed by the FDM process is shown in Figure 2.
J. Mar The most-commonly used 3D printing process (46% as of 2018) is a material extrusion technique called fused deposition modeling (FDM) [35]. The FDM process builds a three-dimensional object from a CAD model by successively adding material layer by layer. The flowchart of a general FDM process is shown in Figure 1. The initial part of a 3D printing pipeline is to create digital prototypes using various CAD software. After establishing the digital 3D model, the integrated CAM software is utilized to generate the G-code file. Technically, G-code is a machine language used to drive the motion of motors. The main task of the CAM software is to slice the digital 3D model into layers and generate a G-code subsequently. The G-code will be transferred to the control module of a 3D printer after slicing. The 3D physical model of the KVLCC2 printed by the FDM process is shown in Figure 2.

Stereolithographic Model
The prerequisite for producing any 3D printed part is a digital 3D model or a CAD file. The most widely used format of a digital 3D model in 3D printing is stereolithographic format due to its compact structure. The triangular facet is adopted to describe closed complex geometries in a stereolithographic model. The two format types of the model, the ASCII, and the binary are illustrated in Table 1. An ASCII stereolithography (STL) file begins with a solid name and concludes with the end solid name. It contains any number of triangles, where each n or v is a floating-point number. Since the ASCII STL files can become very large, a binary version of STL exists. A binary STL file has an 80-character header, where a four-byte little-endian unsigned integer indicating the number of triangular facets is present. Each triangle is described by twelve 32-bit floating-point numbers: three for the normal and then three for the X/Y/Z coordinate of each vertex. The file simply ends after the last triangle.

Stereolithographic Model
The prerequisite for producing any 3D printed part is a digital 3D model or a CAD file. The most widely used format of a digital 3D model in 3D printing is stereolithographic format due to its compact structure. The triangular facet is adopted to describe closed complex geometries in a stereolithographic model. The two format types of the model, the ASCII, and the binary are illustrated in Table 1. An ASCII stereolithography (STL) file begins with a solid name and concludes with the end solid name. It contains any number of triangles, where each n or v is a floating-point number. Since the ASCII STL files can become very large, a binary version of STL exists. A binary STL file has an 80-character header, where a four-byte little-endian unsigned integer indicating the number of triangular facets is present. Each triangle is described by twelve 32-bit floating-point numbers: three for the normal and then three for the X/Y/Z coordinate of each vertex. The file simply ends after the last triangle. In this study, the benchmark models from the SIMMAN workshop are selected for investigation and validation [36]. The stereolithographic models and corresponding facets amount of the KVLCC2, KCS, and DTMB 5415 are shown in Figure 3.  In this study, the benchmark models from the SIMMAN workshop are selected for investigation and validation [36]. The stereolithographic models and corresponding facets amount of the KVLCC2, KCS, and DTMB 5415 are shown in Figure 3. The vast number of facets describes the geometry of the hull precisely, as can be seen from the figures above. The facets distributed on the stereolithographic models are non-uniform. The smaller these facets are, the smoother the model will be, and the bigger the size of the model. The middle portion of a body plan generated from the original stereolithographic model is provided in Figure 4. The portion with small curvature is marked by A, and the portion with large curvature is marked by B in the left figure of Figure 4. An undeniable fact can be observed: the number of triangular facets in part A is less than part B. The large, flat plane in the stereolithographic models is easy to describe by a few large triangular facets, which could save storage space. In contrast, a curved plane needs a huge amount of small triangular facets due to the significant curvature. The middle sections of the generated body plan are provided in the right figure. Note that the body plan is produced by planetriangular intersection processing; the points are scattered in the body plan due to the non-uniform distribution of triangular facets in the stereolithographic model. That is, the points of each section will accumulate in the portions with large curvature, as well as influencing the seakeeping calculation results negatively. Besides, the non-uniform distribution is a common feature in stereolithographic models, which can be observed in other models as shown in Figure 3. The vast number of facets describes the geometry of the hull precisely, as can be seen from the figures above. The facets distributed on the stereolithographic models are non-uniform. The smaller these facets are, the smoother the model will be, and the bigger the size of the model. The middle portion of a body plan generated from the original stereolithographic model is provided in Figure 4.  In this study, the benchmark models from the SIMMAN workshop are selected for investigation and validation [36]. The stereolithographic models and corresponding facets amount of the KVLCC2, KCS, and DTMB 5415 are shown in Figure 3. The vast number of facets describes the geometry of the hull precisely, as can be seen from the figures above. The facets distributed on the stereolithographic models are non-uniform. The smaller these facets are, the smoother the model will be, and the bigger the size of the model. The middle portion of a body plan generated from the original stereolithographic model is provided in Figure 4. The portion with small curvature is marked by A, and the portion with large curvature is marked by B in the left figure of Figure 4. An undeniable fact can be observed: the number of triangular facets in part A is less than part B. The large, flat plane in the stereolithographic models is easy to describe by a few large triangular facets, which could save storage space. In contrast, a curved plane needs a huge amount of small triangular facets due to the significant curvature. The middle sections of the generated body plan are provided in the right figure. Note that the body plan is produced by planetriangular intersection processing; the points are scattered in the body plan due to the non-uniform distribution of triangular facets in the stereolithographic model. That is, the points of each section will accumulate in the portions with large curvature, as well as influencing the seakeeping calculation results negatively. Besides, the non-uniform distribution is a common feature in stereolithographic models, which can be observed in other models as shown in Figure 3. The portion with small curvature is marked by A, and the portion with large curvature is marked by B in the left figure of Figure 4. An undeniable fact can be observed: the number of triangular facets in part A is less than part B. The large, flat plane in the stereolithographic models is easy to describe by a few large triangular facets, which could save storage space. In contrast, a curved plane needs a huge amount of small triangular facets due to the significant curvature. The middle sections of the generated body plan are provided in the right figure. Note that the body plan is produced by plane-triangular intersection processing; the points are scattered in the body plan due to the non-uniform distribution of triangular facets in the stereolithographic model. That is, the points of each section will accumulate in the portions with large curvature, as well as influencing the seakeeping calculation results negatively.
Besides, the non-uniform distribution is a common feature in stereolithographic models, which can be observed in other models as shown in Figure 3.
A body plan (KCS) of 20 slices is selected to further illustrate the non-uniform feature inherited from the stereolithographic model. For each slice, the distance interval in space between two adjacent points is calculated. The statistical characteristics of these intervals are presented in Figure 5. The top figure in Figure 5 is the boxplot of the intervals. In each box, the central red line mark indicates the median, and the bottom and top edges of the blue box indicate the 25th and 75th percentiles, respectively. The whiskers extend to the most extreme data points not considered outliers, and the outliers are plotted individually using the red "+" symbol. In addition, the standard normal distribution is utilized to evaluate the non-uniform feature of the intervals. The probability density function f(x) is calculated based on the mean value µ and the standard deviation σ as in Equation (1), where the sample data X refer to the intervals. The f(x) of the intervals is provided in the bottom figure of Figure 5, where the dimension of the horizontal axis could be in centimeters or meters, depending on the stereolithographic model. Mar A body plan (KCS) of 20 slices is selected to further illustrate the non-uniform feature inherited from the stereolithographic model. For each slice, the distance interval in space between two adjacent points is calculated. The statistical characteristics of these intervals are presented in Figure 5. The top figure in Figure 5 is the boxplot of the intervals. In each box, the central red line mark indicates the median, and the bottom and top edges of the blue box indicate the 25th and 75th percentiles, respectively. The whiskers extend to the most extreme data points not considered outliers, and the outliers are plotted individually using the red "+" symbol. In addition, the standard normal distribution is utilized to evaluate the non-uniform feature of the intervals. The probability density function f(x) is calculated based on the mean value μ and the standard deviation σ as in Equation (1), where the sample data X refer to the intervals. The f(x) of the intervals is provided in the bottom figure of Figure 5, where the dimension of the horizontal axis could be in centimeters or meters, depending on the stereolithographic model. The probability density indicates that the intervals below five make up the majority of all sections. Several outliers could be observed in each section from the boxplot. These outliers attribute to the non-uniform features in a triangular mesh. Therefore, these features cause similar nonuniformity in the produced body plan. The re-triangulation and optimization of the stereolithographic model will refine the mesh by rearranging the triangular facets. However, the nonuniformity features will remain [37][38][39]. Since the non-uniform features cannot be fully eliminated by optimization methods. A section redistribution algorithm is developed and integrated into the model slicing pipeline, which is explained in the model slicing algorithm.

Model Slicing Pipeline
The proposed method is developed based on the model slicing algorithm from the FDM process. Since the essential element of stereolithographic models is a triangular facet, the major task of the slicing operation is therefore converted into a space intersection problem. The intersection lines between a ship model and a moving cutting plane are computed from bow to stern in order to obtain the longitudinal offset data. Furthermore, similar operations can be performed from starboard to port side or each waterline to generate a transverse or vertical lines plan under a set of uniform or non- The probability density indicates that the intervals below five make up the majority of all sections. Several outliers could be observed in each section from the boxplot. These outliers attribute to the non-uniform features in a triangular mesh. Therefore, these features cause similar non-uniformity in the produced body plan. The re-triangulation and optimization of the stereolithographic model will refine the mesh by rearranging the triangular facets. However, the non-uniformity features will remain [37][38][39]. Since the non-uniform features cannot be fully eliminated by optimization methods. A section redistribution algorithm is developed and integrated into the model slicing pipeline, which is explained in the model slicing algorithm.

Model Slicing Pipeline
The proposed method is developed based on the model slicing algorithm from the FDM process. Since the essential element of stereolithographic models is a triangular facet, the major task of the slicing operation is therefore converted into a space intersection problem. The intersection lines between a ship model and a moving cutting plane are computed from bow to stern in order to obtain the longitudinal offset data. Furthermore, similar operations can be performed from starboard to port side or each waterline to generate a transverse or vertical lines plan under a set of uniform or non-uniform cutting planes. The schematic diagram of the slicing process between a KCS stereolithographic model and a slice plane is shown in Figure 6. uniform cutting planes. The schematic diagram of the slicing process between a KCS stereolithographic model and a slice plane is shown in Figure 6. The set of cutting planes is prepared at the beginning of the slicing pipeline. The constant layer thickness along the slicing direction is adopted in uniform slicing and non-uniform slice planes and can be easily produced by a density distribution function or adaptive method based on surface curvature. Examples based on both uniform and non-uniform slicing are provided in Figure 7. In the time-domain strip method, the instantaneous wetted surface is calculated at each time step [19]. Particularly, the wetted surface varies during the six degrees of freedom (6-DOF) motions of a ship. The dense body plan under arbitrary roll and pitch angles can be conveniently produced by rotating the stereolithographic model. Based on the incident wave elevation around the ship, an instantaneous wetted surface could be integrated by limiting the body plan with the elevation in each longitudinal location. The slicing process of the rotated model is presented in Figure 8. Although the instantaneous free surface in the time-domain calculation is fully non-linear, an average steady water plane is adopted here to illustrate the effectiveness. The present method shows the applied potential in time-domain strip theory. In detail, the intersection lines between each triangular facet and the moving cutting plane need to be calculated in advance. The triangle-plane intersection problem can be considered as a specific uniform cutting planes. The schematic diagram of the slicing process between a KCS stereolithographic model and a slice plane is shown in Figure 6. The set of cutting planes is prepared at the beginning of the slicing pipeline. The constant layer thickness along the slicing direction is adopted in uniform slicing and non-uniform slice planes and can be easily produced by a density distribution function or adaptive method based on surface curvature. Examples based on both uniform and non-uniform slicing are provided in Figure 7. In the time-domain strip method, the instantaneous wetted surface is calculated at each time step [19]. Particularly, the wetted surface varies during the six degrees of freedom (6-DOF) motions of a ship. The dense body plan under arbitrary roll and pitch angles can be conveniently produced by rotating the stereolithographic model. Based on the incident wave elevation around the ship, an instantaneous wetted surface could be integrated by limiting the body plan with the elevation in each longitudinal location. The slicing process of the rotated model is presented in Figure 8. Although the instantaneous free surface in the time-domain calculation is fully non-linear, an average steady water plane is adopted here to illustrate the effectiveness. The present method shows the applied potential in time-domain strip theory. In detail, the intersection lines between each triangular facet and the moving cutting plane need to be calculated in advance. The triangle-plane intersection problem can be considered as a specific In the time-domain strip method, the instantaneous wetted surface is calculated at each time step [19]. Particularly, the wetted surface varies during the six degrees of freedom (6-DOF) motions of a ship. The dense body plan under arbitrary roll and pitch angles can be conveniently produced by rotating the stereolithographic model. Based on the incident wave elevation around the ship, an instantaneous wetted surface could be integrated by limiting the body plan with the elevation in each longitudinal location. The slicing process of the rotated model is presented in Figure 8. Although the instantaneous free surface in the time-domain calculation is fully non-linear, an average steady water plane is adopted here to illustrate the effectiveness. The present method shows the applied potential in time-domain strip theory. uniform cutting planes. The schematic diagram of the slicing process between a KCS stereolithographic model and a slice plane is shown in Figure 6. The set of cutting planes is prepared at the beginning of the slicing pipeline. The constant layer thickness along the slicing direction is adopted in uniform slicing and non-uniform slice planes and can be easily produced by a density distribution function or adaptive method based on surface curvature. Examples based on both uniform and non-uniform slicing are provided in Figure 7. In the time-domain strip method, the instantaneous wetted surface is calculated at each time step [19]. Particularly, the wetted surface varies during the six degrees of freedom (6-DOF) motions of a ship. The dense body plan under arbitrary roll and pitch angles can be conveniently produced by rotating the stereolithographic model. Based on the incident wave elevation around the ship, an instantaneous wetted surface could be integrated by limiting the body plan with the elevation in each longitudinal location. The slicing process of the rotated model is presented in Figure 8. Although the instantaneous free surface in the time-domain calculation is fully non-linear, an average steady water plane is adopted here to illustrate the effectiveness. The present method shows the applied potential in time-domain strip theory. In detail, the intersection lines between each triangular facet and the moving cutting plane need to be calculated in advance. The triangle-plane intersection problem can be considered as a specific In detail, the intersection lines between each triangular facet and the moving cutting plane need to be calculated in advance. The triangle-plane intersection problem can be considered as a specific case in which the two planes intersect. Figure 9 presents the intersection between a triangle facet and a cutting plane. case in which the two planes intersect. Figure 9 presents the intersection between a triangle facet and a cutting plane.
The facet intersects the cutting plane if the value of di appears not only positive but negative. For the sake of generality, V1 is above the cutting plane and V2 and V3 are below the cutting plane, as shown in Figure 6. The orientation of the intersection vector is determined by the cross product of two planes' normal, as in Equation (3).
The intersection segment satisfies both the plane equation and the linear equation.
Substituting Equation (5) into Equation (4) gives The ratio t in the linear equation is solved by Equation (7).
Based on the dot products obtained by Equation (2), one of the endpoints of the intersection segment S1 can be conveniently calculated using Equation (8), which saves the computational cost. In addition, the other endpoint S2 can be derived from a similar approach. Figure 9. The intersection of a triangular facet with a cutting plane.
The plane α is the cutting plane with a normal n a . A triangular facet consists of three vertices (V 1 , V 2 , V 3 ) coincident with the plane β. The corresponding normal of plane β is n b . The vector → S 1 S 2 is the intersection segment between the facet and the cutting plane. Besides, the vector → I 1 I 2 is the intersection segment of the two planes. The situations of intersection can be determined by the geometric positional relationship between the three vertices and the plane. By selecting an arbitrary point P 0 on the plane α, vectors are formed by the three vertices and the selected point. The dot products between the three vectors and the normal n a will be utilized to determine the intersection as shown in Equation (2).
The facet intersects the cutting plane if the value of d i appears not only positive but negative. For the sake of generality, V 1 is above the cutting plane and V 2 and V 3 are below the cutting plane, as shown in Figure 6. The orientation of the intersection vector is determined by the cross product of two planes' normal, as in Equation (3).
The intersection segment satisfies both the plane equation and the linear equation.
Substituting Equation (5) into Equation (4) gives The ratio t in the linear equation is solved by Equation (7).
Based on the dot products obtained by Equation (2), one of the endpoints of the intersection segment S 1 can be conveniently calculated using Equation (8), which saves the computational cost. In addition, the other endpoint S 2 can be derived from a similar approach.
Consequently, all intersection segments are obtained by sweeping all the facets of the hull stereolithographic model with all cutting planes. Note that the produced points on each slice are non-uniformly distributed, as described in Figure 5. Redundancy appears in the produced body plans, especially the portions where the curvature changes rapidly. Therefore, a redistribution algorithm (Algorithm 1) is implemented to obtain a uniform, low-redundancy section.  The developed algorithm improves the uniformity by redistributing the original sliced data. The algorithm consists of three basic steps. At first, the minimum value of the spacing of a slice is searched by sweeping the spacing of each two adjacent points. Then, the points are densified according to the minimum value by linear interpolation. Finally, uniform data under a specified amount of points are obtained by extracting the points from the densified data at a constant interval. The final dense body plan is generated based on the processes above. The pipeline of the stereolithographic model slicing is illustrated in Figure 10.

Dense Body Plan Results
The produced body plans are shown in the following figures based on the present method. Initially, the data of the models under a rolling and pitching angle are presented in Figure 11. The sliced data along the transverse and vertical directions are provided in Figure 12. Besides, the body plans under different slices are given in Figure 13. Finally, the original sliced data are compared with the redistributed data in Figure 14.

Dense Body Plan Results
The produced body plans are shown in the following figures based on the present method. Initially, the data of the models under a rolling and pitching angle are presented in Figure 11. The sliced data along the transverse and vertical directions are provided in Figure 12. Besides, the body plans under different slices are given in Figure 13. Finally, the original sliced data are compared with the redistributed data in Figure 14.

Dense Body Plan Results
The produced body plans are shown in the following figures based on the present method. Initially, the data of the models under a rolling and pitching angle are presented in Figure 11. The sliced data along the transverse and vertical directions are provided in Figure 12. Besides, the body plans under different slices are given in Figure 13. Finally, the original sliced data are compared with the redistributed data in Figure 14.

Validation of the Produced Body Plan
Accuracy and efficiency are major factors to evaluate the present method. Two cases are provided to validate the accuracy of the produced data. The hull volume is widely used in the hydrostatic analysis. In addition, the WSA is often utilized for the estimation of the ship resistance in maneuvering simulations. Thus, the hull volumes and WSAs are calculated and compared with the published value. The published value of hull volumes and WSAs are given in Table 2. The hull volumes and WSAs can be integrated by produced body plan using Equation (9). Ci and Ai are the wetted perimeter length and area of each section, respectively.

Validation of the Produced Body Plan
Accuracy and efficiency are major factors to evaluate the present method. Two cases are provided to validate the accuracy of the produced data. The hull volume is widely used in the hydrostatic analysis. In addition, the WSA is often utilized for the estimation of the ship resistance in maneuvering simulations. Thus, the hull volumes and WSAs are calculated and compared with the published value. The published value of hull volumes and WSAs are given in Table 2. The hull volumes and WSAs can be integrated by produced body plan using Equation (9). C i and A i are the wetted perimeter length and area of each section, respectively.
The calculated volumes are compared with the published value. The relative errors of the calculated results under different slice numbers are given in Table 3. The maximum relative error of the calculated hull volumes is 0.915%, which occurred in the DTMB 5415 model of 20 slices. The underwater geometry cannot be accurately described by only 20 slices due to the rapid change around the bulbous bow. In the case of 500 slices, the relative error decreases to 0.057%, 0.001%, and 0.003%. The relative errors are converged when the number of slices exceeds 200. Twelve regression-based approaches of WSA are reviewed and listed in Table 4 in order to validate the adaptability of the present method. The calculated WSAs the 20 and 200 slices of the body plans are selected to perform the comparison.  [43] No.14 200 slices Figure 15 illustrates the relative error of each WSA compared with the published value. The results show that the traditional methods have difficulties covering all the ship types. Note that the estimation methods (1-12) of WSA are mostly based on the regression analysis from a group of ships, or the approximate calculation using simplified ship geometries. Thus, the estimation approaches have limited accuracy, especially across the diversity of possible vessel types. However, the relative error of DTMB at 20 slices is greater than the other two ships, which demonstrates the limitation of a body plan at a small number of slices. The accuracy and adaptability based on the dense body plan of 200 slices are higher than all the estimation methods. Moreover, the present method could handle the ships equipped with special underwater appendages such as a bulbous bow. or the approximate calculation using simplified ship geometries. Thus, the estimation approaches have limited accuracy, especially across the diversity of possible vessel types. However, the relative error of DTMB at 20 slices is greater than the other two ships, which demonstrates the limitation of a body plan at a small number of slices. The accuracy and adaptability based on the dense body plan of 200 slices are higher than all the estimation methods. Moreover, the present method could handle the ships equipped with special underwater appendages such as a bulbous bow. The computational cost is another vital factor for evaluating this method. The number of facets varies according to the size of stereolithographic models, which leads to different time costs. Time costs of the three models are provided in Figure 16. The time cost is linearly increased from 20 slices

Hydrodynamic Coefficients Database
The speed and heading of a ship change during maneuvering motion, which is not only caused by the operating engine and rudder but also by the environmental disturbances. However, ship motion in waves is commonly predicted at a constant forward speed, wave direction, and wave frequency in the seakeeping simulations. Both the forward speed and wave parameters vary in the actual sea. The different combination of the parameters leads to a different motion performance. Therefore, a hydrodynamic coefficients database is adopted to simulate the ship maneuvering motion in short-crested irregular waves. It is a powerful tool to handle the varying wave-induced forces/moments due to the changing ship speed and wave encounter angle. Such forces/moments could be practically superposed using coefficients of discrete frequency and phase based on the linear hypothesis. For instance, a unified model expanded from Fossen [9] is shown in Equation (10).
where M is the mass matrix, and H, P, R refer to the hull, propeller, and rudder according to the modular concept. The wind, current, excitation, and drift subscripts represent the wind, current, wave excitation forces/moments, and wave steady drift forces/moments. NL-Hull is the non-linear viscous forces/moments in maneuvering. Note that the IRF (K), the added mass (A), and the damping (B) are in a 6 × 6 matrix at each speed and wave direction. The convolution integral of IRFs captures the so-called fluid memory effects. Therefore, it is necessary to consider the varying IRFs, added

Hydrodynamic Coefficients Database
The speed and heading of a ship change during maneuvering motion, which is not only caused by the operating engine and rudder but also by the environmental disturbances. However, ship motion in waves is commonly predicted at a constant forward speed, wave direction, and wave frequency in the seakeeping simulations. Both the forward speed and wave parameters vary in the actual sea. The different combination of the parameters leads to a different motion performance. Therefore, a hydrodynamic coefficients database is adopted to simulate the ship maneuvering motion in short-crested irregular waves. It is a powerful tool to handle the varying wave-induced forces/moments due to the changing ship speed and wave encounter angle. Such forces/moments could be practically superposed using coefficients of discrete frequency and phase based on the linear hypothesis. For instance, a unified model expanded from Fossen [9] is shown in Equation (10).
where M is the mass matrix, and H, P, R refer to the hull, propeller, and rudder according to the modular concept. The wind, current, excitation, and drift subscripts represent the wind, current, wave excitation forces/moments, and wave steady drift forces/moments. NL-Hull is the non-linear viscous forces/moments in maneuvering. Note that the IRF (K), the added mass (A), and the damping (B) are in a 6 × 6 matrix at each speed and wave direction. The convolution integral of IRFs captures the so-called fluid memory effects. Therefore, it is necessary to consider the varying IRFs, added mass, and damping during the time-domain simulation. These IRFs are pre-calculated offline based on the frequency-domain hydrodynamic coefficients using Equation (11).
By utilizing a database, the IRFs are interpolated by instant speed and wave encounter angle in the time-domain simulations. The first-order wave forces/moments in random seas can be easily computed based on a Response Amplitude Operator (RAO) database and 2D wave spectra. In addition, second-order wave load coefficients could be interpolated from a similar database. The databases are made as a function of ship speed, wave frequency and wave direction. For example, the steady wave-induced drift forces/moments could be obtained by Equation (12). Here, C is the wave drift force coefficient at a specific combination of the ship speed U, the wave frequency ω k , and the wave encounter angle χ e . Besides, the symbol T refers to the average wave period, S ζζ is the wave spectra, θ n is the wave spreading angle, H 1/3 represents the wave significant height, and F refers to the wave drift forces/moments.
A ship hydrodynamic coefficients database is constructed in this study based on the produced dense body plan and strip method. The strip theory-based Kochin-function method (SKFM) is used to calculate the hydrodynamic coefficients. This method has been proved effective in maneuvering simulation in waves [20]. The generated dense body plan of a real bulk carrier is shown in Figure 17.
A ship hydrodynamic coefficients database is constructed in this study based on the produced dense body plan and strip method. The strip theory-based Kochin-function method (SKFM) is used to calculate the hydrodynamic coefficients. This method has been proved effective in maneuvering simulation in waves [20]. The generated dense body plan of a real bulk carrier is shown in Figure 17. The hydrodynamic coefficients are organized by ship speed, wave frequency and wave encounter angle in the database. Part of the coefficients in the database are provided in Figure 18   The hydrodynamic coefficients are organized by ship speed, wave frequency and wave encounter angle in the database. Part of the coefficients in the database are provided in Figure 18 to 21. The added mass, damping coefficients, motion RAOs and wave load RAOs are matrix values in the database. The ship speed varies from 0 to 12 knots, and the wave encounter angle varies from 0° to 180°. The frequency interval is 0-0.3 Hz.      The hydrodynamic coefficients vary due to different combinations of ship speed, wave frequency and wave encounter angle. However, the second-order drift force coefficients appear to have poor stability under the condition of a following sea, which needs further research. The IRFs can be computed in advance of simulations to save the time cost. Parts of the IRF (K22) are calculated based on Equation (11) and shown in Figure 22. The contribution of first-order wave load to ship motion can be computed based on the IRFs and linear convolution in time-domain simulation. The hydrodynamic coefficients vary due to different combinations of ship speed, wave frequency and wave encounter angle. However, the second-order drift force coefficients appear to have poor stability under the condition of a following sea, which needs further research. The IRFs can be computed in advance of simulations to save the time cost. Parts of the IRF (K 22 ) are calculated based on Equation (11) and shown in Figure 22. The contribution of first-order wave load to ship motion can be computed based on the IRFs and linear convolution in time-domain simulation.
The hydrodynamic coefficients vary due to different combinations of ship speed, wave frequency and wave encounter angle. However, the second-order drift force coefficients appear to have poor stability under the condition of a following sea, which needs further research. The IRFs can be computed in advance of simulations to save the time cost. Parts of the IRF (K22) are calculated based on Equation (11) and shown in Figure 22. The contribution of first-order wave load to ship motion can be computed based on the IRFs and linear convolution in time-domain simulation.

Conclusions
In this paper, motivated by 3D printing technology, a fast and accurate dense body-plan generation method based on the stereolithographic model is developed. Not only the longitudinal

Conclusions
In this paper, motivated by 3D printing technology, a fast and accurate dense body-plan generation method based on the stereolithographic model is developed. Not only the longitudinal but the transverse and vertical sections can be produced. The number of sections is fully customized to meet the requirement of different calculations. Moreover, the model under arbitrary rotation is also covered in the present method, which could contribute to computing the instantaneous wetted surface area in the strip method. A section redistribution algorithm is purposed for post-processing the non-uniform distribution features inherited from the triangular facets. The redistributed results indicate that the non-uniformity is eliminated in each section based on the purposed algorithm. Two cases are provided to validate the produced data. Firstly, the hull volumes calculated by using the created body plan are compared with the published value. Note that the KVLCC2 model has 24,586 facets, and the relative error of the hull volume is 0.003% under 500 slices. Secondly, the calculated WSAs are compared with the existing methods. The calculated WSAs based on the produced body plan are highly consistent with the published value at 200 slices. Consequently, the dense body plan is accurate and adaptive for such calculations. Furthermore, time costs under different slices are presented. The slicing operation of the three models is less than 1 s under 200 slices, which implies the applied potential of the method in time-domain strip theory due to its computational efficiency. A ship hydrodynamic coefficients database is constructed based on the produced data and 2D strip method. The first-order wave memory effect can be computed by Cummins equations based on the interpolation of IRFs. The second-order wave drift effect in the time-domain can be calculated based on the interpolation of drift force coefficients in the database. It is necessary to emphasize that the present method turned the time-consuming task of generating a body plan manually into building the 3D model of a ship hull, and the latter is only a one-time task. Accordingly, the present hydrodynamic database is convenient for simulating the ship maneuvering motion in a seaway. This paper only used the body-plan generation method. Sufficient validation of time-domain simulation should be done in the next step. However, it is reasonable to believe that the method could provide the high-quality data fundamental for strip theory to perform maneuvering simulations in waves.