Experimental Study on Development of Mooring Simulator for Multi Floating Cranes

: In this study, the coupled motion of a mooring system and multiﬂoating cranes were analyzed. For the motion analysis, the combined equations of motions of the mooring line and multiﬂoating cranes were introduced. The multibody equations for ﬂoating cranes were derived from the equations of motion. The ﬁnite element method (FEM) was used to derive equations to solve the stretchable catenary problem of the mooring line. To verify the function of mooring simulator, calculation results were compared with commercial mooring software. To validate the analysis results, we conducted an experimental test for offshore operation using two ﬂoating crane models scaled to 1:40. Two ﬂoating crane models and a pile model were established for operation of uprighting ﬂare towers. During the model test, the motion of the ﬂoating cranes and tensions of the mooring lines were measured. Through the model test, the accuracy of the mooring analysis program developed in this study was veriﬁed. Therefore, if this mooring analysis program is used, it will be possible to perform a mooring analysis simulation at the same time as a maritime work simulation. mean absolute deviation between simulation and model test by the maximum amplitude. This error rate, which was less than 5%, shows the reliability and validity of the simulation results.


Introduction
There has been the need to perform dynamic analysis of mechanical systems. Consequently, several simulation tools that focus on traditional mechanical systems, such as automobiles and robots, have been developed. In recent years, the need for motion analysis in the shipbuilding industry has been increasing.
However, shipbuilding and offshore industries generally differ from conventional mechanical systems in terms of purpose, size, and shape. Therefore, there are limits to the application of analysis tools developed for existing mechanical systems in the shipbuilding and offshore industries. For this reason, simulation tools were developed in several studies for the dynamic analysis of shipbuilding and offshore operation processes [1,2]. In the previously developed simulation program, the multibody dynamics theory is applied to perform a more accurate dynamic analysis of crane operations on the sea. However, mooring lines are defined as linear springs in most programs, and spring constants are determined based on some assumptions and linearization. However, for some offshore operations, such as installation and decommissioning, mooring analysis may be a vital point of analysis. Therefore, in this paper, the following research was conducted in relation to mooring analysis.

•
The mooring analysis module was implemented by referring to the existing research results [3]. • This mooring analysis module has been integrated with the multibody dynamics module that we have previously developed [1,2].
• A model test was conducted to verify the performance of the integrated program.
The other parts of this paper are organized as follows. In Section 2, previous studies related to this study are reviewed, and the scope of the study is defined. Section 3 describes the equations for motion analysis of multifloating crane and mooring line. Section 4 shows the comparison between simulation results of the developed simulator and commercial program, and Section 5 describes the model test results for the verification of the developed program. Section 6 presents a summary of this study and briefly describes a future related study.

Related Studies
Several studies on mooring line analysis have been conducted. In the research of Garrett, D.L., and Kim, B.W., the motion of the vessel and the mooring system are analyzed simultaneously, called the fully coupled method [4][5][6]. The equations of motion are formulated for both the floating body and the mooring line. Numerical models and procedures provide accurate and efficient global modeling, and the analysis results are more accurate than the simulation results obtained using the linear spring methods. However, in these studies, the floating structure is modelled as a single body, which is not capable of simulating offshore operation of floating cranes.
On the other hand, Ham, S.H., and Cha., J.H. present the simulation methods for offshore operation of floating cranes [7,8]. In these studies, the floating crane was modelled as a multibody, which has joints that can simulate the motion of offshore cranes. However, the mooring line was assumed as a linear spring line for the motion analysis. For this assumption, the determination of the spring constant is crucial. Moreover, the complicated motion of the mooring line was approximated based on the linear spring model, so the results of the analysis cannot be sufficiently precise for some aspects, in spite of the reduced simulation time [9,10].
Recently, Gutiérrez-Romero, J.E. presented the nonlinear dynamic FEM mooring model to analyze the response of moored floating wind turbines [11]. In this research, a FEM mooring module was implemented, and the floating wind turbine was modeled as a multibody, so that this research carried out an accurate simulation. However, to model the wind turbine, Gutiérrez-Romero, J.E. used FAST, which can be used only for wind turbine analysis. So, we cannot use the result of this research for analysis of floating cranes.
On the other hand, as the need for an accurate analysis of the dynamic behavior in shipbuilding and offshore industries has increased, we developed a floating crane simulation tool in previous research [1,2]. The programs developed from these studies were based on multibody dynamics, and more sophisticated motions of vessels and cranes can be simulated. Therefore, we decided that the simulation tool is suitable for simulation of offshore operations using floating cranes and planned to integrate the mooring analysis module developed in this study into the simulation tool. Figure 1 shows an example of the GUI program and its application.
Therefore, in this research, we implemented nonlinear mooring analysis module by referring to the existing research results [3]. For simultaneous analysis of mooring line and multibody floating cranes, the interface between the mooring analysis module and existing floating crane simulation tools was designed and implemented. To validate the developed simulator, a model test was conducted. The input and output data for the static analysis are listed in Table 1.

Equations Used for Simulation Program
For motion analysis of the moored floating crane, the equations of motion for the floating cranes were derived based on multibody dynamics, and the concept of the constrained force was applied. The mooring line was assumed to be stretchable catenary and modeled as a finite element. In this section, the equations of motion used for motion analysis are described.

Equations of Motion for Floating Crane
All of the floating crane systems in operation are multibody systems, in which multiple rigid bodies are joined. Thus, a dynamic analysis module for the multibody system was developed in this study. The computer methods used in the automated dynamic analysis of multibody systems that consist of rigid bodies are generally divided into two main approaches [12,13]. In the first approach, the system configuration is identified using a set of Cartesian coordinates that describe the locations and orientations of the bodies. This approach is called the "augmented formulation" [14]. In the second approach, relative joint variables are used to formulate a minimum set of differential equations of motion, The input and output data for the static analysis are listed in Table 1.

Equations Used for Simulation Program
For motion analysis of the moored floating crane, the equations of motion for the floating cranes were derived based on multibody dynamics, and the concept of the constrained force was applied. The mooring line was assumed to be stretchable catenary and modeled as a finite element. In this section, the equations of motion used for motion analysis are described.

Equations of Motion for Floating Crane
All of the floating crane systems in operation are multibody systems, in which multiple rigid bodies are joined. Thus, a dynamic analysis module for the multibody system was developed in this study. The computer methods used in the automated dynamic analysis of multibody systems that consist of rigid bodies are generally divided into two main approaches [12,13]. In the first approach, the system configuration is identified using a set of Cartesian coordinates that describe the locations and orientations of the bodies. This approach is called the "augmented formulation" [14]. In the second approach, relative joint variables are used to formulate a minimum set of differential equations of motion, and two types of formulation use relative joint variables: embedding formulation and recursive formulation. The augmented formulation has the disadvantage of increasing the equation of motion according to the number of objects, but the analysis result is intuitive and has the advantage of calculating the constrain force. Because of this advantage, in this study, the augmented formulation, expressed as Equation (1), was used [13,15]. q (4) where q is the generalized coordinate matrix, λ is the constraint force matrix, C(q) is the constraint matrix, J(q) is the velocity transformation matrix, M is the mass and inertia matrix, F is the external force matrix, and N is the three-dimensional transformation matrix. A Newton-Euler equation can be expressed simply based on the center of gravity. However, since a general mechanical system uses a generalized coordinate system in consideration of the user's convenience, coordinate transformation is necessary, and Equations (2)-(4) shows the process of performing this coordinate transformation, which is calculated by multiplying J(q) T and J(q) on M, F, and N [15].

Equations of Motion for Mooring Line
An analytical solution of the nonstretchable catenary mooring line has been provided [16]. In this study, however, the mooring line was assumed to be stretchable, and several numerical solutions of the problem are available. Among the previous studies, the equation derived by Kim was adopted in this study [3]. In this research, the mooring line was divided into several elements to model the mooring line, and Equation (5) was derived by integrating the equations of each element. The equation for static analysis using the FEM is expressed as follows: where u is the position of the node for each finite element. By solving this static equation, the static equilibrium position, tension, and mooring forces can be calculated. Each term in Equation (5) is presented as follows: where s 0 is the initial length of each finite element of the mooring line, E is Young's modulus of the mooring line, and A is the section area of the mooring line.
In Equations (8)- (11), u ge is the global position of the node for each finite element, R is the transformation matrix that transforms the global vector to a local vector, and u e is the local position of the node for each finite element. (12) [K e ] = λ[C] (13) [C] In Equations (12)- (18), λ is the Lagrange multiplier, l is the distance between the local coordinates, and w is the weight of the unit length of the mooring line. B, C, and D are symbols used to simplify the equation. B stands for Equation (10), C stands for Equation (14), and D stands for Equation (9). The Newton-Rapson method was adopted to solve Equation (5) [17].

Definition of Interface between Modules
The multibody analysis module and the mooring line analysis module are independent of each other because their governing equations are different. Therefore, the analysis modules should translate the mooring information to each other. For the translation, an interface between two modules has been developed. The interface and input/output details are depicted in Figure 2. This figure shows only limited information exchanged between the mooring analysis module and the multibody analysis module targeted in this paper.
The mooring line and multibody analysis modules are integrated in the analysis program of the moored multifloating crane, and the mooring line analysis module exchanges the data with a database, a multibody analysis module, and a postprocessing module.
The mooring line analysis module receives the following types of information from the database: physical information about the mooring line; -physical information about the seabed.
The mooring line analysis module receives the following types of information from the multibody analysis module: -position, velocity, and acceleration information about the multifloating cranes.
The mooring line analysis module provides the following types of information to the multibody analysis module and postprocessing module.
tension (mooring force) is provided to the multibody analysis module; -tension (mooring force), position, and shape of the mooring rope are provided to the postprocessing module.
J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW 6 of 2 -tension (mooring force), position, and shape of the mooring rope are provided to th postprocessing module.

Verification of Mooring Line Analysis Module
As mentioned in the previous chapter, the mooring line and multibody analysis mod ules are integrated in the analysis program. Among them, the multibody analysis module were developed and verified in previous research [15]. To verify the mooring line analysi module, we conducted several simulations of mooring examples. The example in Sectio 4.1 is referenced in [3], and the example in Section 4.2 is an example of pulling test o floating mass made up in this study. The calculation results were compared with [3] an a commercial program for mooring analysis called Orcaflex, respectively, to verify th calculation accuracy of the module developed in this paper.

Verification of Mooring Line Analysis Module
As mentioned in the previous chapter, the mooring line and multibody analysis modules are integrated in the analysis program. Among them, the multibody analysis modules were developed and verified in previous research [15]. To verify the mooring line analysis module, we conducted several simulations of mooring examples. The example in Section 4.1 is referenced in [3], and the example in Section 4.2 is an example of pulling test of floating mass made up in this study. The calculation results were compared with [3] and a commercial program for mooring analysis called Orcaflex, respectively, to verify the calculation accuracy of the module developed in this paper.

Verification Examples of Static Analysis
In the previous section, equations for analysis were presented. In this section, a static analysis solver was developed first, and for verification, the simulation results and the nonstretchable analytic solution, which is described in research [3], were compared. An example of a mooring line is shown in Figure 3. To make the mooring line nonstretchable in the developed solver, Young's modulus was set as infinite number. Figure 4 shows the comparison results. In the graph, the line represents the simulation results obtained using the developed solver, and the dots represent the analytical solution-they are the same.

Verification Examples of Static Analysis
In the previous section, equations for analysis were presented. In this section, a static analysis solver was developed first, and for verification, the simulation results and the nonstretchable analytic solution, which is described in research [3], were compared. An example of a mooring line is shown in Figure 3. To make the mooring line nonstretchable in the developed solver, Young's modulus was set as infinite number. Figure 4 shows the comparison results. In the graph, the line represents the simulation results obtained using the developed solver, and the dots represent the analytical solution-they are the same.     The slope angle q varies from 0° to 60°, where q is the angle between the x-axis and dotted linear line L. The analysis results depict a parabolic catenary curve shape. As the axial stiffness, which is proportional to EA, decreases, the sag increases quadratically, and the tension of the mooring line decreases, which indicates that the mooring force also decreases. These results were compared with other research [3], and they were the same. Therefore, we could confirm that the stretchable mooring line analysis function is accurate enough to use.  The slope angle q varies from 0 • to 60 • , where q is the angle between the x-axis and dotted linear line L. The analysis results depict a parabolic catenary curve shape. As the axial stiffness, which is proportional to EA, decreases, the sag increases quadratically, and the tension of the mooring line decreases, which indicates that the mooring force also decreases. These results were compared with other research [3], and they were the same. Therefore, we could confirm that the stretchable mooring line analysis function is accurate enough to use.     The calculation process described above is performed without considering the bottom reaction of the sea. However, the mooring line is generally in contact with the seabed. In addition, mooring lines in contact with the seabed have a significant effect on mooring performance; hence, the bottom condition must be considered. The method of considering the ground conditions in the analysis module was applied as follows: where is the height of the mooring node, is the seabed depth, is the seabed stiffness, is the weight of the mooring node, is the penetration depth of the mooring line, is the damping coefficient of the seabed, is the damping ratio, and is the The calculation process described above is performed without considering the bottom reaction of the sea. However, the mooring line is generally in contact with the seabed. In addition, mooring lines in contact with the seabed have a significant effect on mooring performance; hence, the bottom condition must be considered. The method of considering the ground conditions in the analysis module was applied as follows: where z is the height of the mooring node, z c is the seabed depth, k c is the seabed stiffness, w c is the weight of the mooring node, δ c is the penetration depth of the mooring line, c c is the damping coefficient of the seabed, ε c is the damping ratio, and m c is the weight of the mooring node. The contact force of the seabed was computed using the above equation. An example of the mooring line is shown in Figures 9-12, which depict the static equilibrium position and tension of the numerical example. The slope angle q was varied from 20° to 60°. The analysis results show a parabolic catenary curve shape. The sag increases quadratically with a decrease in the axial stiffness (the axial stiffness is proportional to EA). Moreover, the tension of the mooring line decreases with the axial stiffness, which suggests that the mooring force also decreases.  The contact force of the seabed was computed using the above equation. An example of the mooring line is shown in Figures 9-11, which depict the static equilibrium position and tension of the numerical example. The slope angle q was varied from 20 • to 60 • . The analysis results show a parabolic catenary curve shape. The sag increases quadratically with a decrease in the axial stiffness (the axial stiffness is proportional to EA). Moreover, the tension of the mooring line decreases with the axial stiffness, which suggests that the mooring force also decreases. The contact force of the seabed was computed using the above equation. An example of the mooring line is shown in Figures 9-12, which depict the static equilibrium position and tension of the numerical example. The slope angle q was varied from 20° to 60°. The analysis results show a parabolic catenary curve shape. The sag increases quadratically with a decrease in the axial stiffness (the axial stiffness is proportional to EA). Moreover, the tension of the mooring line decreases with the axial stiffness, which suggests that the mooring force also decreases.

Verification Examples of Pulling Test
In this section, a performed simulation of pulling test for evaluating the function of the mooring line analysis module is described. The example of the simulation is shown in Figure 13. This example was made up for comparison with commercial programs. It is possible to calculate the restoration force generated by the mooring line by moving the

Verification Examples of Pulling Test
In this section, a performed simulation of pulling test for evaluating the function of the mooring line analysis module is described. The example of the simulation is shown in Figure 13. This example was made up for comparison with commercial programs. It is possible to calculate the restoration force generated by the mooring line by moving the

Verification Examples of Pulling Test
In this section, a performed simulation of pulling test for evaluating the function of the mooring line analysis module is described. The example of the simulation is shown in Figure 12. This example was made up for comparison with commercial programs. It is possible to calculate the restoration force generated by the mooring line by moving the point mass, and through this, we can assess whether it is the same as the calculation performance of the verified commercial program.
For each of the two cases, the floating bodies were forced to move from −55 to +55 m. The tension of wire rope was calculated, a restoring force curve was plotted, and the simulation results were compared with those of OrcaFlex [18]. The simulation conditions for Case 1 and Case 2 are outlined as follows.
A. Case 1 • Vertical distance between anchor and floating body: 300 m.

Verification Examples of Pulling Test
In this section, a performed simulation of pulling test for evaluating the function o the mooring line analysis module is described. The example of the simulation is shown in Figure 13. This example was made up for comparison with commercial programs. It i possible to calculate the restoration force generated by the mooring line by moving th   For each of the two cases, the floating bodies were forced to move from −55 to +55 m. The tension of wire rope was calculated, a restoring force curve was plotted (Figure 13), and the simulation results were compared with those of OrcaFlex [18]. The simulation conditions for Case 1 and Case 2 are outlined as follows.

Experimental Study
To validate the developed program, we performed an experimental test for offshore operation in the ocean engineering basin of the Korea Research Institute of Ships and Ocean Engineering (KRISO). The ocean engineering basin is 56 × 30 × 4.5 m and is equipped with a wave generator capable of regular, irregular and multidirectional waves, a current generator, and a wind generator as shown in Figure 14. For this experimental test, two floating crane models scaled to 1:40 were prepared. Detailed specification of the crane model is shown in Table 2. Two floating crane models and one pile model were used to conduct operation of uprighting flare towers. Figure 15 shows a scene of the model test. A 1200-ton floating crane was moored with a spring type mooring line; on the other hand, 2000-ton floating crane was moored with a catenary type mooring line. For more accurate performance verification of pure mooring line and crane dynamics analysis, this experiment was conducted under the sea state 0 condition-the so called calm water condition.   Figure 16 shows mooring points for the 1200-ton and 2000-ton crane models. The mooring coordinates of each crane model are listed in Table 3. Positions of fairlead and anchor are expressed based on local coordinate system and global coordinate system, respectively.    Figure 16 shows mooring points for the 1200-ton and 2000-ton crane models. The mooring coordinates of each crane model are listed in Table 3. Positions of fairlead and anchor are expressed based on local coordinate system and global coordinate system, respectively.   The connection points of the mooring lines are listed in Table 3. The properties and pretension details of each mooring line for 1200-ton crane are as follows.  The connection points of the mooring lines are listed in Table 3. The properties and pretension details of each mooring line for 1200-ton crane are as follows. The connection points of the mooring lines are listed in Table 3. The properties and pretension details of each mooring line for 1200-ton crane are as follows.

•
Spring constant: 15.625 N/m. The procedure of the model test is shown in Figure 17. Figure 18 shows the control input values for the model test and simulation. The 1200-ton crane's backstay and boom started to control at 280 s, the 2000-ton crane's backstay did not rotate and the boom started to rotate at 100 and 280 s. The 1200-ton crane wire started to control at 270 s and the 2000-ton crane wire started at 100 s. The procedure of the model test is shown in Figure 17. Figure 18 shows the control input values for the model test and simulation. The 1200-ton crane's backstay and boom started to control at 280 s, the 2000-ton crane's backstay did not rotate and the boom started to rotate at 100 and 280 s. The 1200-ton crane wire started to control at 270 s and the 2000-ton crane wire started at 100 s. The control input values of the simulation and the model test were matched, and the motion of the flare tower, hoisting wire rope tension, and mooring line tension were compared. Figure 19 shows a comparison of the motion of the flare tower. The yellow and red lines represent the results of the model test and simulation, respectively. The trends of model tests and simulations were similar. In the case of surge motion, as it was erected, the flare tower moved in the direction of the 2000-ton crane. In the case of sway motion and yaw motion, the flare tower moved out of the center and then came back while the backstay of the 1200-ton crane was rotating. Roll motion rarely occurred. Since it is an erect work, the pitch angle and heave elevation increased. In the case of the model test, the roll and yaw motions can be seen as measurement errors as the flare tower was erected The control input values of the simulation and the model test were matched, and the motion of the flare tower, hoisting wire rope tension, and mooring line tension were compared. Figure 19 shows a comparison of the motion of the flare tower. The yellow and red lines represent the results of the model test and simulation, respectively. The trends of model tests and simulations were similar. In the case of surge motion, as it was erected, the flare tower moved in the direction of the 2000-ton crane. In the case of sway motion and yaw motion, the flare tower moved out of the center and then came back while the backstay of the 1200-ton crane was rotating. Roll motion rarely occurred. Since it is an erect work, the pitch angle and heave elevation increased. In the case of the model test, the roll and yaw motions can be seen as measurement errors as the flare tower was erected after 350 s. The maximum error rate between simulation and model test was 0.56% for surge, 0.19% for sway, 1.14% for heave, 5.72% for roll, 0.23% for pitch, and 0.34% for yaw.  In addition, the maximum error rate for the mooring line of the 2000-ton crane and was 3.9%, and the maximum error rate for the mooring line of the 1200-ton crane was 4.7%. The error rate was calculated by dividing the error of the mean absolute deviation between simulation and model test by the maximum amplitude. This error rate, which was less than 5%, shows the reliability and validity of the simulation results.  In addition, the maximum error rate for the mooring line of the 2000-ton crane and was 3.9%, and the maximum error rate for the mooring line of the 1200-ton crane was 4.7%. The error rate was calculated by dividing the error of the mean absolute deviation between simulation and model test by the maximum amplitude. This error rate, which was less than 5%, shows the reliability and validity of the simulation results. The free decay test in the model experiment was performed to obtain the natural period of the floating crane. In the case of the 2000-ton crane, the natural period was 18.11 s for surge, 19.81 s for sway, 1.36 s for roll, 1.60 s for pitch and 12.56 s for yaw. In the case of a 1200-ton crane, the natural period is 13.41 s for surge, 14.73 s for sway, 1.38 s for roll, 1.25 s for pitch and 7.69 s for yaw. In the case of heave motion, it was difficult to measure because the damping was too large. As shown in Figure 20, the tension of the 2000-ton crane vibrates for a period of approximately 20 s, which is the natural frequency of the 2000-ton crane model. Akin to the general model test, the operation was started in static equilibrium in this experiment. However, it was not easy to maintain a perfect static equilibrium in actual experiments. Due to this, very small natural vibration was measured, but the value was not large and did not significantly affect the experimental results. This phenomenon did not occur in the simulation because the initial state of the simulation was static equilibrium.
In addition, the amplitudes of the mooring tension in the model test and simulation are different, which is a limitation of the resolution performance of the measurement equipment. The free decay test in the model experiment was performed to obtain the natural period of the floating crane. In the case of the 2000-ton crane, the natural period was 18.11 s for surge, 19.81 s for sway, 1.36 s for roll, 1.60 s for pitch and 12.56 s for yaw. In the case of a 1200-ton crane, the natural period is 13.41 s for surge, 14.73 s for sway, 1.38 s for roll, 1.25 s for pitch and 7.69 s for yaw. In the case of heave motion, it was difficult to measure because the damping was too large. As shown in Figure 20, the tension of the 2000-ton crane vibrates for a period of approximately 20 s, which is the natural frequency of the 2000-ton crane model. Akin to the general model test, the operation was started in static equilibrium in this experiment. However, it was not easy to maintain a perfect static equilibrium in actual experiments. Due to this, very small natural vibration was measured, but the value was not large and did not significantly affect the experimental results. This phenomenon did not occur in the simulation because the initial state of the simulation was static equilibrium.
In addition, the amplitudes of the mooring tension in the model test and simulation are different, which is a limitation of the resolution performance of the measurement equipment.

Conclusions
A few programs are available for the multibody analysis of offshore operations, such as installation using a multicrane, considering catenary mooring analysis. Therefore, the aim of this study was to develop a nonlinear mooring analysis code that can be applied to the floating multibody analysis code.
In this study, the equations of motion for calculating the mooring force connected to the multibody system were introduced. The multibody equations for floating cranes were

Conclusions
A few programs are available for the multibody analysis of offshore operations, such as installation using a multicrane, considering catenary mooring analysis. Therefore, the aim of this study was to develop a nonlinear mooring analysis code that can be applied to the floating multibody analysis code.
In this study, the equations of motion for calculating the mooring force connected to the multibody system were introduced. The multibody equations for floating cranes were

Conclusions
A few programs are available for the multibody analysis of offshore operations, such as installation using a multicrane, considering catenary mooring analysis. Therefore, the aim of this study was to develop a nonlinear mooring analysis code that can be applied to the floating multibody analysis code.
In this study, the equations of motion for calculating the mooring force connected to the multibody system were introduced. The multibody equations for floating cranes were derived from the equations of motion. The finite element method (FEM) was used to derive equations to solve the stretchable catenary problem of the mooring line.
To validate the developed program, we conducted an experimental test for offshore operation. For this experimental test, two floating crane models scaled to 1:40 were prepared. Two floating crane models and one pile model were used to conduct the operation of uprighting flare towers. During the model test, the motion of the floating cranes and tensions of the mooring lines were measured. By comparing the simulation results and the model test results, we confirmed an error rate of 5% or less in the case of the mooring line. Therefore, we confirmed that the developed mooring simulator could be effectively used for multifloating crane simulations.