A Coupled Hydrodynamic–Structural Model for Flexible Interconnected Multiple Floating Bodies

: Evaluating the structural safety and seakeeping performance of very large ﬂoating structures (VLFS) using the rigid module ﬂexible connector (RMFC) method remains challenging due to the complexity of the coupled hydrodynamic–structural responses in this system. In this study, a coupled hydrodynamic–structural frequency–time domain model is developed based on the RMFC method employing the planar Euler–Bernoulli beam elements to investigate the dynamic responses of multi-module ﬂoating systems. To reveal the dynamic characteristics of the systems, the coupled hydrodynamic–structural responses are investigated using a frequency–time-domain numerical model with viscous correction, in which the mass and stiffness attributes of connectors are incorporated into the system. Given the effects of hydrodynamic interaction, consideration is given to the case of three modular boxes connected by ﬂexible beams aligned in series in shallow water to validate the present model. Higher efﬁciency and accuracy can be found in the system using viscous correction in potential ﬂow theory and introducing state–space model to replace the convolution terms in the Cummins equation for the time domain. Moreover, this model can be extended to a considerable number of ﬂoating modules, which provides possibilities to analyze N-module ﬂoating systems.


Introduction
With the acceleration of the exploitation of ocean resources, multiple floaters near offshore operations are becoming more common, such as multi-module floating systems [1], offloading operations from FLNG to LNG carrier [2,3], offshore platform float-over deck installations [4][5][6] and floating wind-wave power-generation platforms [7,8].These systems may experience strong hydrodynamic interactions, resulting in enormous wave loads generated in the gaps between two floaters.On the other hand, the multi-module floating system has become increasingly prosperous in offshore activities due to the need for large size and the functional integration of floating structures, e.g., the MOB (mobile offshore base) [9][10][11] and floating airports [12].Such structures require a high degree of integrity, due to the complex hydrodynamic interactions and dynamic interactions between different modules via the connectors [13,14].Additional attention should therefore be paid to mitigating the relative movements of adjacent modules in the system subjected to environmental load.Therefore, the assessment of the coupling dynamic effects of the interconnected multi-module floating system is of paramount importance to improving the safety and the seakeeping performance of offshore structures operating in the form of adjacent multiple floating bodies.Since much literature has been dedicated to investigating the effect of hydrodynamic interactions, such as gap resonance [6,[15][16][17][18], this study will concentrate on developing a coupled hydrodynamic-structural frequency-time-domain model based on the previous work by Chen et al. [1].
Connectors can be classified into rigid ones and flexible ones according to their characteristics.Rigid connectors between modules are usually welded and hinged, which are characterized by simple structure and easy manufacture.However, rigid connectors are unable to release the bending stress and shear stress induced by the large deflection of the system, resulting in prominent fatigue problems and short service life.To avoid excessive stress caused by rigid connectors, research on flexible connectors has gained prosperity in recent years.The structural dynamic responses of the interconnected multi-module system are generally investigated using the rigid module flexible connector (RMFC) method.In this method, it is assumed that the stiffness of connectors is much smaller than that of the floating modules; meanwhile, structural deformation occurs on the connector while the modules themself do not deform.
Ever since the concept of very large floating structures (VLFS) was put forward, various connector designs have been proposed by scholars worldwide.Some previous connector designs for VLFS were summarized by Jiang et al. [19].In some sea conditions, the natural frequency of the overall system is close to that of the connectors, which may result in large relative responses of the floating bodies and further dangerous overloading of the connectors [20].To reduce loads on the connectors, Haney [21] proposed several compliant connector configurations.On top of this, the designs using rubber cushions and cables [22,23] and using springs [24][25][26] were generally applied to provide flexible connections.On this basis, Qi et al. [27] and Zhang et al. [28] investigated, respectively, the use of ball-spring devices and nylon-rubber materials in the design of a flexible connector.
Apart from proposing new designs for a connector, scholars are also dedicated to investigating the characteristics of dynamic responses of the multi-module system using the RMFC method both experimentally and numerically.Three-dimensional experiments were conducted by Loukogeorgaki et al. [29] to investigate flexible connectors' internal forces of a pontoon-type multi-module floating system.Ding et al. [30] carried out both experimental and numerical analyses to investigate the RMFC model.Michailides et al. [31] developed a numerical analysis framework to evaluate connectors' internal loads and further optimize the configuration of a modular pontoon-type floating structure.Yang et al. [32] established a three-step time-domain method to assess the motions of each module and connector load based on the Boussinesq equation and Cummins equation.Bispo et al. [33,34] investigated the wave interactions of multi-floating systems with articulated and hinged connections using the potential-flow-based model and analytical models.Chen et al. [5] developed a hydrodynamic-structural frequency-domain model by combining the dynamic substructuring method and static condensation as well as considering the flexibility of connectors and hydrodynamic interactions to analyze the catamaran tow operation in the Spar float-over deck installation scenario.The frequency-time-domain method was first proposed by Cummins [35], and this method has been adopted in standard software, such as AQWA [36][37][38].
The connector load and module motion are two key indicators to evaluate the safety and seakeeping performance of the system.Wang et al. [39] found that connector stiffness is the most important factor affecting relative motion and load conditions.Fu et al. [13] further proved that the stiffness of flexible connectors has a great influence on the hydroelastic responses of the multi-floating system.By comparing connectors with different stiffness, Gao et al. [14] suggested that proper stiffness can effectively reduce hydroelastic responses in the system.Wang et al. [40] carried out a parametric study to investigate the connector characteristics, and appropriate stiffness and damping coefficients have been demonstrated to release the load while retaining reasonable relative motion between adjacent modules.On the other hand, the high-efficiency time-domain analysis model is one of the research interests.The State-Space Model (SSM) is widely used in analyzing marine operations, which has proven to be more efficient than directly solving the Cummins equation with a convolution term.Chen et al. [41][42][43] and Zou et al. [44] proposed a constant-parameter time-domain model (CPTDM) by replacing the convolution terms in the Cummins equation with SSMs to investigate the complex dynamics in the float-over deck installation for a jacket platform with a single barge.This method was then further developed to analyze other marine operation systems, such as multi-floater float-over systems [6,45] and marine renewable energy devices [46].However, due to the limitation of current analysis methods, the RMFC model usually requires simplification, such as using springs and dampers to replace the connection component.From the aforementioned study, a brief conclusion can be drawn that connector loads and module motions are significantly influenced by connector stiffness and damping coefficients.Therefore, it is of great necessity to develop a robust and efficient coupled frequency-time-domain model for a multi-module floating system, which can consider connector flexibility.
In this study, a hydrodynamic-structural model for an interconnected multi-module floating system is developed, which is based on the hydrodynamic results obtained from the common hydrodynamic analysis code, e.g., the commercial code AQWA.At the same time, to eliminate the adverse effects of distortion of hydrodynamic resonance obtained in the potential flow theory system on the results, the damping lid method [1,6,47,48] is used for viscous correction.The frequency-domain simulations are first verified with the results calculated by AQWA in the free-floating status.When further considering the connector system, the effects of the bending stiffness of the connectors are investigated, since the corresponding frequency-domain model has been validated in previous studies [49][50][51].Based on the frequency-domain results, the constant-parameter hydrodynamic-structural time-domain model is established for a multi-module floating system connected end to end and the accuracy of the developed time-domain model is comprehensively discussed.
The paper is arranged as follows.Section 2 introduces the methodology for developing the hydrodynamic-structural model in the frequency domain and time domain.Section 3 discusses the application of the developed numerical model with a three-module system connected end to end, in which the influences of the connector parameters on the overall system behavior are also investigated.Finally, some concluding remarks are given in Section 4.

Frequency-Domain Model
In this research, a hydrodynamic-structural analysis model has been used to analyze the dynamics characteristics in the series of multiple floating systems based on the twostage method proposed by Sun et al. [49][50][51].In the first stage, a commercial boundary element method (BEM) code, ANSYS-AQWA, is employed to calculate the first-order wave excitation forces and hydrodynamic coefficients.In the second stage, the dynamic model for the multi-module floating system with flexible connections is established, which is schematically shown in Figure 1.As illustrated in Figure 1, the motions of the barges are indicated as X b , the motions of the action point between the flexible substructures and barges are defined as master DOF (degrees of freedom), and the nodal motions on the substructures are deemed as slave DOF, which are designated as X m and X s , respectively.The establishment of this dynamic model can help to condense the flexible substructures to the barges through master nodes in the form of a condensed mass and stiffness matrix using static condensation; therefore, the wave-induced dynamics of the interconnected system can be conveniently obtained.By assuming that the only external forces on the substructure are the constraint forces Fm(t) at the connections, the motion equation for the substructure can be written in the following form [52,53]: where the subscripts m and s represent master and slave DOF, respectively, and the superscript T denotes transpose.
In this state, the displacement and acceleration of the flexible substructure can be expressed as: Since the resonance responses of the flexible structure itself are not taken into account, the damping terms can be further ignored.Therefore, the second term in the equation can be written as: Further ignoring the influence of inertia force, the slave DOF s X and the accelera- tion vector s X  can be expressed by the master DOF m X and m X  in the following form: ( ) ( ) ( ) where is the coordinate transformation matrix.
On this basis, the equation of the displacement and acceleration of the substructure can be rewritten as: where I denotes the identity matrix and TG is the transpose matrix.Substituting Equation (6) into Equation (1) leads to:  By assuming that the only external forces on the substructure are the constraint forces F m (t) at the connections, the motion equation for the substructure can be written in the following form [52,53]: where the subscripts m and s represent master and slave DOF, respectively, and the superscript T denotes transpose.In this state, the displacement and acceleration of the flexible substructure can be expressed as: Since the resonance responses of the flexible structure itself are not taken into account, the damping terms can be further ignored.Therefore, the second term in the equation can be written as: Further ignoring the influence of inertia force, the slave DOF X s and the acceleration vector .. X s can be expressed by the master DOF X m and .. X m in the following form: ..
where T T ms is the coordinate transformation matrix.On this basis, the equation of the displacement and acceleration of the substructure can be rewritten as: ..
where I denotes the identity matrix and T G is the transpose matrix.Substituting Equation (6) into Equation (1) leads to: Therefore, an analytical equation involving the master DOF and the constraining force on the master nodes only is constructed.The mass and stiffness matrix on the master nodes of the substructures can be expressed as follows: The harmonic responses of multiple hydrodynamic interacting structures can be obtained by solving linear equations in linear potential flow theory, which are commonly referred to as Response Amplitude Operators (RAOs).The motions of N floating structures under free-floating conditions can be expressed by frequency-dependent hydrodynamic coefficients as follows: where X(iω) represents the 6N × 1 complex vector of multiple bodies' RAOs; N is the number of bodies; M, A H (ω) and B H (ω) are the (6N) × (6N) structural mass, hydrodynamic added mass, and added damping matrices, respectively; C bb is the hydrostatic stiffness matrix assembled by each 6 × 6 sub-hydrostatic stiffness of individual structure along the diagonal; FE (iω) is the 6N × 1 wave excitation force at frequency ω.On this basis, by taking the Fourier transform of both sides of Equation ( 7), the frequency-domain response equation of rigid floating modules with a flexible connection can be written as: The above equation provides the ability to solve the motion response and constraint force of the master DOF on the flexible structure while solving the RAOs of the rigid barges with connections.The expressions are as follows: Fm (iω To develop the structural dynamic model in the frequency domain, the connectors are considered to be Euler-Bernoulli beams.When considering the motions in a twodimensional plane projected on the side of the system, the motions of the two nodes at both ends of the space beam element can be simplified as the translation of the x and z axes and the rotation about the y axis, a total of 6 degrees of freedom, which is shown in Figure 2. According to structural mechanics, the 6-DOF stiffness matrix of a beam element moving in the two-dimensional plane is defined as: Since the Euler-Bernoulli beams used in this model are based on the plane-section assumption, the effect of transverse shear force is therefore ignored.ν Φ and z Φ in the matrix above can be viewed as zero, and it can be therefore further simplified as: Considering that the mass distribution of the beam element is uniform, its mass matrix form can be obtained according to structural dynamics: Since the Euler-Bernoulli beams used in this model are based on the plane-section assumption, the effect of transverse shear force is therefore ignored.Φ ν and Φ z in the matrix above can be viewed as zero, and it can be therefore further simplified as: Considering that the mass distribution of the beam element is uniform, its mass matrix form can be obtained according to structural dynamics: Furthermore, the Euler-Bernoulli beam element can be used to obtain the condensed mass and stiffness matrix of the horizontal beam at each node using static condensation, in which the mass matrix contained in Equation ( 8) can be expressed as: M ms = mL 420 In the same way, the stiffness matrices contained in Equation ( 9) have the following form: When the stiffness of the joint is non-zero, it corresponds to the case of a semi-rigid connection, while when the stiffness C is zero, it can be considered to be the case of a hinged connection.

Coupled Time-Domain Model
Since the frequency-domain method is based on linear potential theory, the nonlinearities of the multi-module floating system with connectors are not able to be evaluated.The time-domain model described by the Cummins equation [35] has been widely used in analyzing the nonlinear dynamics of various types of marine multi-floater systems, such as the floating wind-wave power generation platform, catamaran float-over deck installation and multi-module floating system.For the case of the freely floating condition with zero forward speed, the Cummins equation can be written as: where F E (t) denotes the time-domain wave excitation forces, K(t) and A H (∞) is the matrix of the impulse response functions and added mass at infinite frequency, respectively, which can be obtained from the frequency-domain results [54]: Further consideration is given to the aforementioned cohesion stiffness and mass, the coupling equation considering the flexible connectors in the frequency domain can be transformed to the time domain, which can be written as: When solving the Cummins equation, the convolution term with kernel function K(t) within the time range from 0 to t is calculated at each time step, therefore the integrating process is time-consuming and easy to accumulate and amplify errors.Since the convolution term is a linear time-invariant system, one can replace it with either a state-space model or a transfer function to improve its accuracy and efficiency [55].A single-input-single-output (SISO) state-space model can be easily established to be equivalent to each convolution term [42].Once the SISO state-space models for all the impulse response functions are obtained, the multi-input-multi-output (MIMO) state-space model replaces the matrix of the convolution terms [43]: where z(t) represents the state vector; P, Q, and S are the parametric matrices of the MIMO state-space model, which have the following forms [4]: where q denotes the number of degrees of freedom, which is 18 for the three-barge model; k is the order of the individual SISO state-space model.By replacing the convolution term with the MIMO state-space model, a coupled constant parameter hydrodynamic-structural time-domain model (CPHSTDM) can be established based on the Cummins equation: It should be noted that the hydrodynamic-structural system described by Equation ( 30) is significantly more efficient compared to AQWA [6].In addition, the multi-module floating system to be analyzed in this study involves multi-body hydrodynamic interactions with strong resonant phenomena, which will have a significant effect on the hydrodynamic coefficients, resulting in error generated in the frequency-time-domain transformation.To ensure the error lies within an acceptable threshold, the damping lid method is employed in this study.Detailed information about the damping lid method can be found in works by Zou et al. [6] and Chen et al. [1].To quantify the fitting quality of the SSMs, the coefficient R 2 is applied [56]: where K i,j (t l ) denotes the output of the SISO state-space model at time t l for the input being the Dirac delta function, and K i,j is the mean value of Ki,j (t).

Discussion of the Development of CPHSTDM and Parametric Study on Connector Parameters
A three-module floating system aligned in a longitudinal direction is applied to this study in this section, which is based on the numerical study carried out by Chen et al. [1].The results of the hydrodynamic analysis required in the coupled constant parameter hydrodynamic-structural time-domain model (CPHSTDM) are obtained by the standard panel code ANSYS-AQWA; on this basis, the influence of connector parameter characteristics between modules is further considered, and the characteristics of motion responses and loads of the multi-floating system are analyzed by CPHSTDM, which is subjected to wave forces and connection constraints between modules simultaneously.

Particulars of the Analyzed Model
The selected floating body module is a box structure with an aspect ratio of 2:1, and its mass distribution is uniform.The particulars of the single floating module in the prototype are summarized in Table 1.The model to be analyzed in this study remains consistent with Chen et al. [1], in which a series of analyses on a three-module system to investigate the gap resonance phenomenon in shallow water is conducted.Although the previous work mainly investigates the gap resonance phenomenon, this study further investigates the dynamic characteristics of the interconnected three-module floating system with a narrow gap, whose connectors are considered to be Euler-Bernoulli beams.To verify the coupled CPHSTDM, a three-module model, in which adjacent modules are connected end to end, has been analyzed under head sea conditions as illustrated in Figure 3.The water depth is 50 m.As shown in Figure 4, one connector consists of two identical Euler-Bernoulli beams combined through a single joint.The length of each beam is 2.5 m and it is considered that the mass of the beam is uniformly distributed, and the mass of each unit length is taken as 10 kg/m.As mentioned in Section 2, the motion responses at the connecting points on the module are known as X m , while the motion at the joint is defined as X s .

Frequency-Domain Simulations of the Interconnected Three-Module System
First, a full-scale three-module floating system is analyzed using the panel code ANSYS-AQWA in the frequency domain.To meet the accuracy of the subsequent timedomain model calculation, the wave frequencies selected in the process of frequencydomain hydrodynamic analysis cover a wide-enough range from low frequency to high frequency, i.e., from 0.01 rad/s to 2.41 rad/s with a step of 0.05 rad/s to reduce calculation error.According to the analyses conducted by Chen et al. [1] and Zou et al. [6], the overestimated gap resonances can be observed in the multi-module system; therefore, applying the external damping lid on the free surface of the gap fluid is essential, which has been proved to be capable of effectively mitigating the unrealistic resonances.To highlight the multi-body hydrodynamic interactions, the single-module case and three-module case are both calculated in this study.
Figures 5 and 6 show part of the hydrodynamic coefficients.It can be seen that, compared with the single-module condition, due to the existence of adjacent floating bodies, the hydrodynamic coefficient of the wave-facing module in the multi-module system has changed significantly, revealing the characteristics of pumping resonance [17].When the damping factor is zero, it means no viscous correction is introduced into the hydrodynamic analysis.Under the influence of the hydrodynamic interactions, the added mass experiences first a maximum and then a minimum value (sometimes becoming negative); meanwhile, the damping term appears at a certain frequency near the phenomenon of the obvious resonance peak.Considering the recommended damping lid factor, which ranges from 0 to 0.2 [57], and the computability problem of the time-domain model, a damping factor with an intermediate value of 0.1 was selected for analysis in this study to tune the numerical results more realistic while reserving the hydrodynamic interaction characteristics.With the introduction of the damping lid factor, the exorbitant oscillations hydrodynamic coefficient curves are effectively controlled, and the intensity has been improved.
system has changed significantly, revealing the characteristics of pumping resonance [17].When the damping factor is zero, it means no viscous correction is introduced into the hydrodynamic analysis.Under the influence of the hydrodynamic interactions, the added mass experiences first a maximum and then a minimum value (sometimes becoming negative); meanwhile, the damping term appears at a certain frequency near the phenomenon of the obvious resonance peak.Considering the recommended damping lid factor, which ranges from 0 to 0.2 [57], and the computability problem of the time-domain model, a damping factor with an intermediate value of 0.1 was selected for analysis in this study to tune the numerical results more realistic while reserving the hydrodynamic interaction characteristics.With the introduction of the damping lid factor, the exorbitant oscillations hydrodynamic coefficient curves are effectively controlled, and the intensity has been improved.After introducing the damping lid factor, the hydrodynamic coefficients calculated by ANSYS-AQWA can be further input into the coupled CPHSTDM for post-processing and analysis.The frequency-domain RAOs of the three-module system in free-floating conditions by solving Equation ( 10) are first compared with the ones obtained by AN-SYS-AQWA to verify the accuracy of the developed program.Figure 7 shows the comparisons of the surge, heave, and pitch RAOs of three different modules of the three-module system in a free-floating state with a gap width of 5 m under the head sea.It can be observed that the motion responses obtained by two different computational tools are highly consistent, which verifies that the hydrodynamic coefficients to be used in the CPHSTDM, i.e., wave force, added mass, and additional damping, are the same as those used in ANSYS-AQWA.[17].When the damping factor is zero, it means no viscous correction is introduced into the hydrodynamic analysis.Under the influence of the hydrodynamic interactions, the added mass experiences first a maximum and then a minimum value (sometimes becoming negative); meanwhile, the damping term appears at a certain frequency near the phenomenon of the obvious resonance peak.Considering the recommended damping lid factor, which ranges from 0 to 0.2 [57], and the computability problem of the time-domain model, a damping factor with an intermediate value of 0.1 was selected for analysis in this study to tune the numerical results more realistic while reserving the hydrodynamic interaction characteristics.With the introduction of the damping lid factor, the exorbitant oscillations hydrodynamic coefficient curves are effectively controlled, and the intensity has been improved.After introducing the damping lid factor, the hydrodynamic coefficients calculated by ANSYS-AQWA can be further input into the coupled CPHSTDM for post-processing and analysis.The frequency-domain RAOs of the three-module system in free-floating conditions by solving Equation ( 10) are first compared with the ones obtained by AN-SYS-AQWA to verify the accuracy of the developed program.Figure 7 shows the comparisons of the surge, heave, and pitch RAOs of three different modules of the three-module system in a free-floating state with a gap width of 5 m under the head sea.It can be observed that the motion responses obtained by two different computational tools are highly consistent, which verifies that the hydrodynamic coefficients to be used in the CPHSTDM, i.e., wave force, added mass, and additional damping, are the same as those used in ANSYS-AQWA.After introducing the damping lid factor, the hydrodynamic coefficients calculated by ANSYS-AQWA can be further input into the coupled CPHSTDM for post-processing and analysis.The frequency-domain RAOs of the three-module system in free-floating conditions by solving Equation ( 10) are first compared with the ones obtained by ANSYS-AQWA to verify the accuracy of the developed program.Figure 7 shows the comparisons of the surge, heave, and pitch RAOs of three different modules of the three-module system in a free-floating state with a gap width of 5 m under the head sea.It can be observed that the motion responses obtained by two different computational tools are highly consistent, which verifies that the hydrodynamic coefficients to be used in the CPHSTDM, i.e., wave force, added mass, and additional damping, are the same as those used in ANSYS-AQWA.
Based on the free-floating model in the frequency domain, the stiffness and mass attributes of the connecting components are introduced into the equation of motion in the frequency domain using an additional matrix, and the three-module floating system under the constraint of the connector is now considered.A type of connector consisting of two beams joined end to end was used, which is illustrated schematically in Figure 4.The tensile rigidity, bending rigidity, and mass characteristics of the connector can be considered simultaneously due to the introduction of the concept of the Euler-Bernoulli beam.By solving Equation ( 11), the frequency-domain responses of the modules at any wave frequency can be obtained, and the effects of the stiffness of the connections can be analyzed.Heave and pitch responses of three modules under flexible connection constraints are plotted in Figure 8. Axial stiffness EA and mass per unit length of the connector are set to 1.0 × 10 10 N and 10 kg/m, respectively, and the kinematic characteristics of the system with different bending stiffness EI are to be discussed.There is no restoring force in the surge direction considered in CPHSTDM by now, whether from mooring force or hydrostatic force.To remain consistent with the following time-domain analysis, axial stiffness EA is set to be a constant.As for the mass of the connector, it should be noted that the mass of the connected component, consisting of four Euler-Bernoulli beams, accounts for a very small proportion of the total mass of the system, therefore the response to the change can also be ignored.Moreover, the stiffness of the connection joint C is selected as the same value as bending stiffness.Based on the free-floating model in the frequency domain, the stiffness and attributes of the connecting components are introduced into the equation of motion frequency domain using an additional matrix, and the three-module floating s under the constraint of the connector is now considered.A type of connector con of two beams joined end to end was used, which is illustrated schematically in Fi The tensile rigidity, bending rigidity, and mass characteristics of the connector considered simultaneously due to the introduction of the concept of the Euler-Be beam.By solving Equation ( 11), the frequency-domain responses of the modules wave frequency can be obtained, and the effects of the stiffness of the connections analyzed.Heave and pitch responses of three modules under flexible connectio straints are plotted in Figure 8. Axial stiffness EA and mass per unit length of th nector are set to 1.0 × 10 10 N and 10 kg/m, respectively, and the kinematic characte of the system with different bending stiffness EI are to be discussed.There is no res force in the surge direction considered in CPHSTDM by now, whether from m force or hydrostatic force.To remain consistent with the following time-domain an axial stiffness EA is set to be a constant.As for the mass of the connector, it sho noted that the mass of the connected component, consisting of four Euler-Be beams, accounts for a very small proportion of the total mass of the system, therefo response to the change can also be ignored.Moreover, the stiffness of the connectio C is selected as the same value as bending stiffness.
As illustrated in Figure 8, it can be seen that the bending stiffness of the con As illustrated in Figure 8, it can be seen that the bending stiffness of the connector would significantly affect the motion of the module, especially at low frequency, within 1 rad/s.Moreover, when EI is relatively small, the responses of each module exhibit behavior similar to those under free-floating conditions.It can be observed that the influence on the maximum pitch of the module is neglectable when EI is less than or equal to the order of 1.0 × 10 9 Nm 2 .With the increase in the bending stiffness, the differences between different modules decreased, because the overall system is more rigid due to the strengthening of the constraint effect.In addition, it should be noted that both the heave and pitch motions of Module 2 showed an obvious trend of weakening around 0.5 rad/s, while this phenomenon did not occur in Module 1 and Module 3.This is because it was in the center of the system and received the connection constraints from both the front and back sides.However, when the bending stiffness EI changes from 1.0 ×   In addition to the frequency-domain responses, as one of the most important characteristics of the RMFC model, the constraint forces of the interconnected structure and loads of the connector itself caused by the relative motions of the adjacent modules are of vital importance for connector design in both academic research and practice.Some typical results are given in Figure 9, showing the absolute values of the constraint forces on each module, respectively, plotted against wave frequency in the mode of surge, heave, and pitch.The peak value of the constraint force acting on each module are summarized in Table 2.As the bending stiffness EI of the connector increases, the constraints imposed on each floating module become progressively more severe.The heave and pitch modes experience increasing constraint forces as the bending stiffness increases, although the rate of increase slows down.On the other hand, changes in stiffness barely affect the constraint forces in the surge direction.Figure 10 presents the effects of stiffness on the connector load and the peak value of the connector loads on common joints are summarized in Table 3.As expected, the load at the common joint increases with stiffness.Notably, the force in the surge direction on the connector is significantly lower than that on the floating body.Nevertheless, the connector bears a larger load which is roughly similar to the load on the floating body in the heave motion.In addition to the frequency-domain responses, as one of the most important characteristics of the RMFC model, the constraint forces of the interconnected structure and loads of the connector itself caused by the relative motions of the adjacent modules are of vital importance for connector design in both academic research and engineering practice.Some typical results are given in Figure 9, showing the absolute values of the constraint forces on each module, respectively, plotted against wave frequency in the mode of surge, heave, and pitch.The peak value of the constraint force acting on each module are summarized in Table 2.As the bending stiffness EI of the connector increases, the constraints imposed on each floating module become progressively more severe.The heave and pitch modes experience increasing constraint forces as the bending stiffness increases, although the rate of increase slows down.On the other hand, changes in stiffness barely affect the constraint forces in the surge direction.Figure 10 presents the effects of stiffness on the connector load and the peak value of the connector loads on common joints are summarized in Table 3.As expected, the load at the common joint increases with stiffness.Notably, the force in the surge direction on the connector is significantly lower than that on the floating body.Nevertheless, the connector bears a larger load which is roughly similar to the load on the floating body in the heave motion.

Discussion of the Development of CPHSTDM and Parametric Study on Connector Parameters
After obtaining the hydrodynamic coefficients in the frequency domain, one can establish a time-domain model based on the Cummins equation.By replacing the convolution terms in the Cummins equation with state-space models (SSMs) and further considering the constraints between modules, the constant-parameter hydrodynamic-structural time-domain model (CPHSTDM) is then developed.
Since the multi-module floating system would experience strong hydrodynamic interactions, the damping lid method is employed to obtain reasonable hydrodynamic coefficients.The damping lid introduced in the multi-floating system not only affects the attenuation efficiency of the impulse response function [1], but also affects the accuracy of the constructed SSM, which has been verified by Zou et al. [6].In this study, the frequencydomain results with damping lid factor α = 0.1 are used to develop the time-domain model.As the key step in CPHSTDM, appropriate parametric matrices of SSMs obtained via system identification based on the calculated impulse response function can be vital for the subsequent time-domain analysis.The established SSMs are quantified by the value of R 2 defined in Equation (31), where the value of 0.97 for R 2 is proved by Duarte et al. [52] to be sufficient to satisfy the time-domain accuracy in the single floating system.
Figure 11 shows that excellent fitting quality of the MDOF SSMs for the calculated K(t) of the windward module is achieved when the order of the SSMs is selected to be 40.
of the constructed SSM, which has been verified by Zou et al. [6].In this study, the fr quency-domain results with damping lid factor α = 0.1 are used to develop t time-domain model.As the key step in CPHSTDM, appropriate parametric matrices SSMs obtained via system identification based on the calculated impulse response fun tion can be vital for the subsequent time-domain analysis.The established SSMs a quantified by the value of R 2 defined in Equation (31), where the value of 0.97 for R 2 proved by Duarte et al. [52] to be sufficient to satisfy the time-domain accuracy in t single floating system.
Figure 11 shows that excellent fitting quality of the MDOF SSMs for the calculate K(t) of the windward module is achieved when the order of the SSMs is selected to be 4 In addition, the corresponding R 2 of the SSMs for the coupled terms between d ferent modules in the selected motions are summarized in Table 4.The R 2 values of t system shown in Figure 11 and Table 4 are greater than 0.99, surpassing the suggeste value of 0.97 by Duarte et al. [56].In addition, the corresponding R 2 of the SSMs for the coupled terms between different modules in the selected motions are summarized in Table 4.The R 2 values of the system shown in Figure 11 and Table 4 are greater than 0.99, surpassing the suggested value of 0.97 by Duarte et al. [56].In this study, surge, heave, and pitch are the only modes of interest when the system is subjected to head sea.Based on the identified SSMs, the coupled nine-DOF time-domain model is established, which can take the constraint effect of the connectors into account.To analyze the three-module interconnected system under the head sea, the stiffness and mass matrix of the connecting system obtained in the frequency domain also play fundamental roles in the CPHSTDM.By solving Equation ( 30), the time-domain dynamic response of each module subjected to specified incident waves can be obtained, in which the fourthorder Runge-Kutta (RK4) integration method is used to solve the Cummins equation with a time step of 0.1 s.The numerical simulation costs about 100 s for a simulation duration of 500 s (AMD Ryzen™ 7 2700 Eight-Core Processor 3.20 GHz).
The application of the substructure method based on the static condensation in a floating-body system makes the CPHSTDM by Equation ( 30) equivalent to the frequencydomain model by Equation (11) in a free-floating state.Therefore, the widely used method of frequency-time-domain verification in free floating can be used for reference, i.e., the accuracy of time-domain model calculation can be estimated by comparing the time-domain results with the corresponding frequency-domain RAOs.
Figure 12 shows a comparison between the RAO-based response and the results of the CPHSTDM for the heave and pitch motions of the windward module.Several connector stiffness values are considered in this comparison.The wave amplitude is 2 m, and the wave frequency is 1.26 rad/s, which is close to the frequency at which hydrodynamic resonance occurs.The results of the CPHSTDM show high accuracy compared to the RAO-based responses, which indicates that the time-domain results are reliable and can be used for further analysis.However, it should be noted that the connector stiffness EI is 1.0 × 10 9 Nm 2 , and a resonance phenomenon occurs, making it impossible to obtain stable time-domain results.This phenomenon may have something to do with the hydrostatic restoring force in pitch-pitch mode, with the pitch-pitch stiffness being 4.1732 × 10 10 Nm/rad.Therefore, for the time-domain comparison, the value of EI equal to 1.0 × 10 9 Nm 2 , is not considered.
Due to the constraint effect of connectors between adjacent modules, the three floating bodies form a whole in the mode of the surge and the motion response of each unit module is consistent.Therefore, only the heave and pitch motions of each module are plotted here to reveal the influence of different connector parameters on the dynamic characteristics of the whole system.The connector between the two beam units is semi-rigid, and therefore, the overall pitch response results demonstrate uniformity for different bending stiffness EI values, in the order of 10 −3 .However, it is crucial to note that the motion amplitude of each module is not necessarily the same.On the other hand, it can also be seen that the heave response amplitude of the floating system under various connection stiffness is different, but there is no significant difference between the results when the stiffness EI is 1.0 × 10 6 Nm 2 and 1.0 × 10 7 Nm 2 , which remain consistent with the conclusion drawn in the frequency-domain model.Due to the constraint effect of connectors between adjacent modules, th floating bodies form a whole in the mode of the surge and the motion response unit module is consistent.Therefore, only the heave and pitch motions of each m are plotted here to reveal the influence of different connector parameters on the d characteristics of the whole system.The connector between the two beam u semi-rigid, and therefore, the overall pitch response results demonstrate uniform The periodic response characteristics can be illustrated from regular wave analysis in CPHSTDM. Figure 13 demonstrates the comparison of the heave motions of the three modules under a regular wave of ω = 1.26 rad/s and wave amplitude of 2 m in the head sea for different connector stiffness.The typical time is chosen when the windward module is at the peak of its own heave motion, which is shown in the right image of the figure with the colors of the modules corresponding to their own motion curves in the left.The areas in slash shadow are the wet surfaces of the modules under calm water while those in dotted shadow are the dry surfaces.The results clearly show the position relationship between the three modules in the system.In addition to the individual motions of each floating body, the relative motion between adjacent modules is also an important concern in multi-floating system application scenarios, such as a floating airport.According to the time-domain results obtained by the CPHSTDM, a box plot is demonstrated to reveal the statistics of the relative heave motion of adjacent modules under different connector parameters in Figure 14.Due to the shielding effect, the relative heave motion between the two leeward bodies, Module 1 and Module 2 is weaker than that of the two windward bodies, Module 2 and Module 3.Moreover, when the bending stiffness is set to 1.0 × 10 8 Nm 2 , the relative motion between Module 2 and Module 3 is significantly greater than those when the value of EI is lower.

Wet surface in hydrostatic environment
Wet surface in hydrostatic environment In addition to the individual motions of each floating body, the relative motion between adjacent modules is also an important concern in multi-floating system application scenarios, such as a floating airport.According to the time-domain results obtained by the CPHSTDM, a box plot is demonstrated to reveal the statistics of the relative heave motion of adjacent modules under different connector parameters in Figure 14.Due to the shielding effect, the relative heave motion between the two leeward bodies, Module 1 and Module 2 is weaker than that of the two windward bodies, Module 2 and Module 3.Moreover, when the bending stiffness is set to 1.0 × 10 8 Nm 2 , the relative motion between Module 2 and Module 3 is significantly greater than those when the value of EI is lower.

Concluding Remarks
A coupled constant-parameter hydrodynamic-structural model (CPHTDM) ha been established in this study, which can be used to investigate the dynamics of a mul ti-module floating system with flexible connectors between adjacent floating bodies.A three-module floating system, which has been studied in the previous study, is furthe analyzed after increasing the gap width to 5 m using the standard panel code AQWA With the employment of the dynamic substructure method, structural parameters such as stiffness and mass of flexible connectors can be fully considered and equipped in the whole system, which provides a basis for the analysis of the interconnected multi-module floating system in both the frequency domain and time domain.In this study, the fre

Concluding Remarks
A coupled constant-parameter hydrodynamic-structural model (CPHTDM) has been established in this study, which can be used to investigate the dynamics of a multi-module floating system with flexible connectors between adjacent floating bodies.A three-module floating system, which has been studied in the previous study, is further analyzed after increasing the gap width to 5 m using the standard panel code AQWA.With the employment of the dynamic substructure method, structural parameters such as stiffness and mass of flexible connectors can be fully considered and equipped in the whole system, which provides a basis for the analysis of the interconnected multi-module floating system in both the frequency domain and time domain.In this study, the frequency-domain model, which has been verified by published papers [5,49,50,58], is further extended to the time domain, and the behaviors of the system are investigated by the developed CPHSTDM.Considering that strong hydrodynamic interactions are a severe problem in the multi-module floating system in which the floating bodies operate with a small gap between each other, the external damping lid method is applied to ensure the precision of the analysis in AQWA.On this basis, the modified frequency-domain results with viscosity correction are transformed into the time domain using the Ogilvie relationship, and the Cummins equation is further used to construct the CPHSTDM, in which the SSM is used to replace the time-consuming convolution term.After finishing the construction of the CPHSTDM, the time-domain response is further demonstrated.Some findings can be confirmed from this investigation.

1.
The effects of hydrodynamic interactions in the multi-module floating system play a crucial role in the response of the system and directly affect the feasibility of the numerical model calculation, especially when transforming the frequency-domain results into the time domain.When the distorted hydrodynamic interaction results in potential theory are corrected, it would lead to stable time-domain results.

2.
Connector parameters such as bending stiffness have a significant impact on overall system performance.Both the response of the floating body and the internal load such as forces and moments obtained through the connections between adjacent modules are affected to varying degrees.For the semi-rigid connection system analyzed here, since the tension-compression characteristics are mainly determined by the axial stiffness EA, the force along the connector layout shows signs of insensitivity to different bending stiffness EI.With the exponential increase of the bending stiffness, loads in the surge direction seem to be insensitive for both connectors and the joint, while the loads in the heave and pitch direction of the connectors exhibit first exponential growth and then a slow-down at EI = 1.0 × 10 9 Nm 2 .

3.
The development of CPHSTDM makes it possible not only to analyze the system behaviors such as the specific motion state at a certain time and the relative motion between adjacent modules under the connection constraints, but also to judge the influence of stiffness selection on the whole system so that the phenomenon of global resonance caused by inappropriate stiffness selection can be avoided in the analysis of practical problems.
In addition to analyzing the three-module floating system aligned in series, the developed numerical model can be extended to N floating bodies provided that the hardware performance is not limited.Meanwhile, the application of the CPHSTDM provides the probability to analyze complex connector configurations.For example, connectors with bending stiffness or hyperstatic marine systems are not available in AQWA, but it is possible by employing the Euler-Bernoulli beam in the CPHSTDM.Therefore, the CPH-STDM is suitable for analyzing complex and flexible connections, which makes the model more applicable to various marine operation scenarios, such as float-over installation and wave energy converters.The use of the RMFC model can be expanded to other scenarios with CPHSTDM.
It should be noted that the coupled numerical model developed in this study only considers the linear wave excitation forces and only regular waves are considered in the timedomain simulation.However, according to the previous works on the SSM [4,6,7,43,46,59], it is recommended to extend the analysis into irregular wave scenarios, and second-order wave forces on the mooring model, such as MoorDyn [60], can be further introduced to complete the whole system's nonlinear dynamic characteristics of the mooring state in the future study.Moreover, it is a limitation that the Euler-Bernoulli beam model used in this study is based on the vertical planar motion assumptions.The 6-DOF response of each module in an environment such as oblique waves cannot be accurately captured by the present model considering only the coupled surge-heave-pitch motions.Hence, future studies are recommended to use the space flexible beam model to achieve full DOF simulation, and therefore many marine operation problems involving complex wave conditions, such as the floating breakwater and floating bridge, can be analyzed in the time domain based on the developed CPHSTDM.In addition, the developed CPHSTDM can also be applied to other configurations, e.g., the floating photovoltaic system as analyzed by Li and Choung [61].

Figure 3 .
Figure 3. Notation of the three-module hydrodynamic-structural model.

Figure 4 .
Figure 4. Plan view of the system and connector modeling.

Figure 5 .
Figure 5.Comparison of the added mass coefficients.

Figure 6 .
Figure 6.Comparison of the added damping coefficients.

Figure 7 .
Figure 7.Comparison between the developed numerical model and the well-proven pan AQWA.

Figure 7 .
Figure 7.Comparison between the developed numerical model and the well-proven panel code AQWA.
10 6 Nm 2 to 1.0 × 10 7 Nm 2 , it has almost no effect on the heave motion of the floating body, and only has a slight change in the pitch motion at some frequencies.
(a) Heave motion of Module 1 (b) Heave motion of Module 2 (c) Heave motion of Module 3 (d) Pitch motion of Module 1 (e) Pitch motion of Module 2 (f) Pitch motion of Module 3

Figure 8 .
Figure 8. Heave and pitch motions of modules under different connection parameters

Figure 8 .
Figure 8. Heave and pitch motions of modules under different connection parameters.

Figure 10 .
Figure 10.Connector loads on common joints between beam elements.

Figure 10 .
Figure 10.Connector loads on common joints between beam elements.

Figure 11 .
Figure 11.Fitting quality of the MDOF SSM for K(t) based on the frequency-domain results wi damping lid factor α = 0.1.

Figure 11 .
Figure 11.Fitting quality of the MDOF SSM for K(t) based on the frequency-domain results with damping lid factor α = 0.1.

Figure 12 .
Figure 12.Comparisons of heave and pitch motions of Module 3 between the RAO-based r and CPHSTDM-based response under the regular wave of ω = 1.26 rad/s and wave amplit m in the head sea for different connector stiffness.

Figure 12 .
Figure 12.Comparisons of heave and pitch motions of Module 3 between the RAO-based response and CPHSTDM-based response under the regular wave of ω = 1.26 rad/s and wave amplitude of 2 m in the head sea for different connector stiffness.

J
. Mar. Sci.Eng.2023, 11, x FOR PEER REVIEW 20 of 25 increase compared to the former two cases.Notably, the response amplitude of Module 1 experiences a slight decrease.No phase shift occurs on the responses of both Module 2 and Module 3, while Module 1 seems to have a positive phase shift when the value of EI is set to 1.0 × 10 8 Nm 2 , which indicates that the heave response of Module 1 lags more compared to the previous cases.(a) EI = 1.0 × 10 6 Nm 2 (b) EI = 1.0 × 10 7 Nm 2 (c) EI = 1.0 × 10 8 Nm 2

Figure 13 .
Figure 13.The system's instantaneous motion state under different connector stiffness when the windward module's heave motion is at the peak at ω = 1.26 rad/s and wave amplitude of 2 m in a head sea.

Figure 13 .
Figure 13.The system's instantaneous motion state under different connector stiffness when the windward module's heave motion is at the peak at ω = 1.26 rad/s and wave amplitude of 2 m in a head sea.With the value of EI increasing from 1.0 × 10 6 Nm 2 to 1.0 × 10 7 Nm 2 , the phases of the responses of all three modules are hardly affected.However, when the EI is set to 1.0 × 10 8 Nm 2 , the response amplitudes of Module 2 and Module 3 exhibit differing extents of increase compared to the former two cases.Notably, the response amplitude of Module 1 experiences a slight decrease.No phase shift occurs on the responses of both Module 2 and Module 3, while Module 1 seems to have a positive phase shift when the

Figure 14 .
Figure 14.Relative motion between adjacent modules in a heave direction with different connecto stiffness.

Figure 14 .
Figure 14.Relative motion between adjacent modules in a heave direction with different connector stiffness.

Table 1 .
The particulars of the single floating module in the prototype.

Table 2 .
The peak value of the constraint force acting on each module.

Table 3 .
The peak value of the connector loads on common joints.

Table 3 .
The peak value of the connector loads on common joints.

Table 4 .
R 2 values for SSMs identified based on the frequency-domain results for the coupling term between different modules.