1. Introduction
At present, offshore structures tend to increase in dimension and complexity. The so-called large multi-body floating offshore structure shown in
Figure 1 and
Figure 2 is a new type of floating offshore structure which has appeared in recent years, which is mainly composed of a substructure, connecting structure (truss or column), upper box structure, and superstructure. The substructure is composed of a row of distributed floating bodies. The length and width of the floating offshore structure are very large, while the height is very small. It is an extremely flat and flexible structure. Its stiffness and vibration natural frequency are low, and the elastic deformation cannot be neglected. The calculation of the floating body as a rigid body may greatly underestimate the bending moment [
1]. In addition, under the wave action of multiple floating bodies, the reflected wave of a floating body will affect the motions and loads on other floating bodies, and waves will be amplified or sheltered in some areas. The fluid flow around the floating bodies is very complex, especially when the spacing between dispersed floating bodies is small and the fluid resonates inside the spacing. In addition, the viscous effect is also obvious, which makes it particularly difficult to predict the motions and loads of floating bodies. Since the size of a single floating body in the substructure is relatively small, strong nonlinear phenomena, such as wave overtopping the floating body and the water entry and exit of the floating body, will occur under severe sea conditions. Therefore, it is necessary to analyze the deformation of this large-scale offshore structure under fluid loads and the influence of floating body deformation using nonlinear hydroelastic methods accounting for the fluid viscosity, so as to obtain the wave loads and structural responses. According to the results of the tank model test, the deformation and displacement of each floating body are consistent at each time in head seas. Therefore, in order to reasonably reduce the difficulty and complexity and also to focus on the main problems, the present study chose to analyze the hydroelastic behavior of a single body instead of an entire floating structure and verify the feasibility of the CFD-FEA coupling method. The large multi-body floating offshore structure has many usages. The floating offshore structure of a length of 300 m is an independent and complete module, i.e., a module of a 300 m long structure alone can also be used as a floating port; in this case, a single module can be analyzed. If a larger length is required, such as that used as a floating airport, several modules can be connected by connectors. In the case of a large floating structure comprising of several modules, it is still reasonable and appropriate to start with the hydroelastic analysis of a single module.
Since the late 1970s, the potential-flow-based theory and methods for hydroelasticity analysis have made remarkable progress. Betts et al. [
2] proposed a two-dimensional hydroelasticity theory. Wu [
3] devised a fully three-dimensional approach using a generalized boundary condition of the fluid–solid interface, and the method has been very successful and widely applied in the community. However, the limitations of the methods based on potential flow theory are also obvious. Compared with potential-flow-based theory, CFD (computational fluid dynamics) methods are more accurate in describing the changes in velocity and pressure in the fluid domain, capturing the nonlinear phenomenon of waves on the free surface, and simulating the fine flow field with hydrodynamic interactions between multiple bodies. Xie et al. [
4], Wang and Guedes Soares [
5], Jiao and Huang [
6], Huang et al. [
7], and Jiao et al. [
8] applied CFD-based methods to analyze the seakeeping and wave loads of ships and offshore structures. It has been found that, compared with potential-flow-based methods, CFD methods can more realistically simulate the flow field and obtain more accurate numerical results.
In recent years, with the development of CFD technology and the improvement of computer performance, some researchers have begun to study how to couple CFD (computational fluid dynamics) and FEA (finite element analysis) for hydroelasticity analysis. El Moctar et al. [
9] and Oberhagemann [
10] combined the CFD method and the Timoshenko beam model for a two-way coupling simulation of ship slamming loads, and validation against the experimental results was also conducted. Because the Timoshenko beam model is relatively simple, only vertical bending and horizontal bending were studied, thus it is not applicable to the torsion problem. Wilson et al. [
11] used a CFD method to study the flow of a ship in waves, and then applied the motions and pressure distribution on the ship’s hull to the finite element model of the ship’s structure for a structural analysis. In order to reduce the amount of calculation, only a one-way coupling method was adopted, i.e., only the influence of fluid on the structure was considered, while the influence of structural deformation on the fluid was neglected. Liu et al. [
12] developed a CFD-FEM (computation fluid dynamics–finite element method) coupling simulation method for estimating the hydroelastic responses of the floating elastic plate. For offshore structures with complex forms, the finite element model for the entire hull structure would also be very complex, thus the amount of time required for the fluid–structure coupling calculation would be very large and result in great difficulties in practical applications. Lakshmynarayanana and Hirdaris [
13] and Lakshmynarayanana and Temarel [
14] studied the symmetrical motions and loads of the S175 container ship in severe waves and under slamming impacts using one-way and two-way CFD-FEA coupling methods. It was shown that the two-way coupling method is able to more accurately estimate the hydroelastic effects of ships in waves. Jiao et al. [
15] used the two-way CFD-FEA coupling method to predict the motions, wave loads, and hydroelastic responses of the S175 container ship in regular waves. All the methods adopted in [
13,
14,
15] are based on the linear Timoshenko beam element to simplify the simulation of the hull structure, neglecting the influence of rotational inertia on the motions and wave loads of offshore structures. This simplification is acceptable for the case of head seas, but not for quartering and beam seas as the effect of rotational inertia is significant and cannot be completely neglected.
To sum up, the existing hydroelastic analysis methods have certain limitations, and there are still some difficulties in calculating the wave loads on the large multi-body floating offshore structure shown in
Figure 1. Previous studies focused on traditional ships, such as barges or container ships. Although traditional ship types can also reach 300 m in length, due to the differences in structural forms, the hydrodynamic and structural response characteristics are different. The hydroelastic analysis method for traditional ship types may not be well suited for large multi-body floating offshore structures. For instance, wave overtopping on the structure is a phenomenon that may behave differently between the two types of structures. To capture such behaviors with complex phenomena, the CFD-FEA coupling method, for its inherent natures of methodology, is a possible solution.
In this study, a CFD-FEA coupling simulation approach is devised, where the fluid domain is solved by the RANS (Reynolds-averaged Navier–Stokes)-based commercial code STAR-CCM+ [
16], and the structure domain is solved by the commercial FE code ABAQUS [
17]. The advantage of this method is that it can fully consider the interaction between elastic structure and fluid through two-way fluid–structure coupling and it can effectively solve strong nonlinear problems such as wave overtopping the floating body and water entry and exit of the floating body. At the same time, it can significantly reduce the computational time of the fluid–structure coupling calculation, and it is applicable to the case of head seas, quartering seas, and beam seas. It is suitable for solving the nonlinear wave loads of large multi-body floating offshore structures accounting for the influence of hydroelasticity. In the present study, a single floating body of the large multi-body floating offshore structure as shown in
Figure 1 is studied, and the hydrodynamic and structural response characteristics of a single floating body of the large multi-body floating offshore structure are studied with a CFD-FEA coupling method. By comparing the potential-flow-based numerical results for wave loads with and without accounting for the elasticity, the CFD-FEA-based method for hydroelastic analysis of flexible floating structures is validated. It is found that the CFD-FEA method has the potential to capture the complex phenomena of the wave–structure interaction.
The paper is organized as follows: the problem and method of solution used in this study are described in
Section 2. Then, in
Section 3, the numerical modeling and mesh sensitivity analysis are carried out with the verification of the numerical method. In
Section 4, the numerical method is validated by comparing the potential-flow-based numerical results for wave loads with and without accounting for the elasticity. Finally, the conclusion is drawn in
Section 5.
2. Problem Statement and Method of Solution
2.1. Problem Statement
In the wave loads on a single floating body of the large multi-body floating offshore structure, the nonlinear factors of the hydrodynamic problem are significant, such as the instantaneous position variation in the body surface, the nonlinearity of the free surface, the elastic deformation of the floating body. These nonlinear factors have a great influence on the loads and motions of the floating body. In order to account for the nonlinear factors, a two-way fluid–structure simulation approach is adopted by coupling CFD and FEA. The fluid domain is solved using the finite volume method (FVM) using STAR-CCM+, and the structure domain is solved by the finite element method (FEM) using ABAQUS.
For the fluid domain, according to Park et al. [
18], the contribution of the compressibility and surface tension of the fluid is small and can be neglected. The mass conservation equation and Navier–Stokes equation are used to describe the flow of viscous fluid:
where
is the fluid pressure,
is the velocity vector,
is the fluid density and
is the dynamic viscosity, and
is the gravity.
For the structure domain, assuming that the structure is a linear elastic material, the motion equation of the structure is obtained using the finite element method:
where
is the mass matrix of the structure,
is the damping matrix of structure,
is the stiffness matrix of the structure,
is the external fluid loads from the CFD solution in the fluid domain, including the fluid pressure and shear force, and
is the motions and deformations of the structure.
2.2. Coupling Approach
The coupling approach is shown in
Figure 3, where
is the initial time and
is the time step. At the initial time, the finite volume method (FVM)-based code STAR-CCM+ is used to calculate the fluid pressure and shear force on the floating body surface, and then the quantities obtained by the FVM are applied to the floating body surface in the finite element model in ABAQUS by means of data mapping. Under the action of fluid pressure and shear force in the fluid domain, the velocity and acceleration of the finite element nodes of the floating body will vary, which leads to the motions and deformations of the fluid–structure coupling interface. Then, the node displacements are transferred to STAR-CCM+ to update the coupling interface. The pressure field and velocity field obtained by STAR-CCM+ are then transferred to the next time step, as well as the velocity and acceleration of the nodes obtained by ABAQUS.
2.3. Coordinate Systems
Three right-handed orthogonal coordinate systems are used to describe the floating body motions: earth-fixed coordinate system , body-fixed coordinate system , and steadily translating coordinate system .
The origin of the earth-fixed coordinate system
lies on the still water surface. The positive
axis is in the direction of the wave propagation, the positive
axis is directed upwards, and
is the wave direction, as shown in
Figure 4.
The body-fixed coordinate system is connected to the ship, with its origin placed at the center of gravity . The positive axis is in the longitudinal forward direction. The positive axis is pointing in the port side direction. The positive axis is perpendicular to the waterline of the hull and is pointing upward.
The steadily translating coordinate system is moving forward at the nominal speed. If the floating body is stationary, the origin of the coordinate system is located at the intersection of the still water surface, the midship cross-section and the center plane, and the directions of the axes are the same as those of the .
2.4. Backbone Beam
A backbone beam is used to simulate the stiffness of the whole floating body. The stress and deformation of the micro-segment of the backbone beam are analyzed. The micro-segment of the backbone beam is shown in
Figure 5 [
19]. The beam is located in the
coordinate system. As shown in
Figure 5c, the reference point
is the intersection of the
axis and the cross-section of the micro-segment when floating still, i.e., the intersection of
axis and section of the micro-segment. In
Figure 5c,
is the mass center of the micro-segment;
is the shear center of the micro-segment;
is the distance between the mass centroid
and
;
is the distance between the shear center
and
; and
is the distance between the center of mass
and the shear center
,
.
For the vertical bending, as shown in
Figure 5a, the differential equations of the motion of the backbone beam are:
For the horizontal bending and torsion, as shown in
Figure 5b,c, the differential equations of the motion of the backbone beam are:
where
is the mass distribution per unit length of the beam;
is the moment of inertia per unit length of the beam;
is the torsional moment of inertia per unit length passing through the mass center
and parallel to the
axis, and its relationship with the moment of inertia
about the
axis is
;
is the vertical bending deformation of the cross-section center;
is the rotation angle of the section due to bending;
is the rotation angle of the cross-section due to torsion;
is the force acting on the unit length of the beam;
is the cross-section bending moment;
is the cross-section shear force;
is the cross-section disturbing force; and
is the roll moment acting on the unit length of the beam. The subscript
denotes the
direction or rotating about the
axis, and the subscript
stands for the
direction or rotating about the
axis. The dot over the variable denotes the time derivative, and the prime indicates the partial derivative of the length
.
According to the Timoshenko beam theory, the relationship between the force and deformation on the section can be obtained as below. For the vertical bending,
For the horizontal bending and torsion,
where
is the shear stiffness of the section;
is the bending stiffness of the section;
is the torsional stiffness of the section;
is the shear damping coefficient;
is the bending damping coefficient;
is torsional damping coefficient;
is the shear angle; and
is the warping stiffness. Subscripts
and
indicate that the quantities are related to a vertical or horizontal motion, respectively.
The geometric relationship of the deformation of the backbone beam is:
In the above Equations (4), (6) and (8) constitute the governing equations of vertical motion, and (5), (7) and (9) constitute the horizontal motion and torsional motion of the backbone beam.
2.5. Section Loads Calculation
Through the above calculation, the pressure distribution, velocity, and acceleration on the surface of the floating body can be obtained, and then the forces and moments at the cross-section of the floating body can be calculated according to the following formula [
20].
where
is the six-component of the section loads,
is the wetted surface area of a floating body over a partial length,
is the total fluid pressure and shear force,
is the additional item,
is the mass matrix of the partial length floating body, and
is the 6-DOF acceleration of the floating body.
3. Numerical Modeling and Simulation
In this section, the numerical modeling is introduced. Then, the extent of the computational domain, boundary condition, time step, and mesh sensitivity analysis are discussed and analyzed, and the numerical method is verified.
3.1. Numerical Modeling
The length of a single floating body of large multi-body floating offshore structure is 300 m, the height is 7.5 m, the width is 5 m, the draft is 5 m, and the mass is evenly distributed along the length of the single floating body. The transverse cross-section of the single floating body is shown in
Figure 6, and the thickness of the single floating body plate is 20 mm. The hollow rectangular beam is used as the backbone beam to simulate the stiffness of the whole floating body. The floating body has 21 stations (0–20 stations) in total, and is divided into 20 sections, as shown in
Figure 7. The red line in
Figure 7 denotes the backbone beam.
The vertical bending stiffness and longitudinal torsional stiffness of the actual floating body cross-section and backbone beam are shown in
Table 1.
The single floating body of a large multi-body floating offshore structure is shown in
Figure 8. The structural model is composed of two parts, the hull surface and the backbone beam, as shown in
Figure 9. The hull surface of the floating body is modeled using shell elements, denoted by the blue surface in
Figure 9. The hull surface is the fluid–structure coupling interface, where the external fluid loads are mapped from the fluid domain, which makes no contribution to the overall stiffness and mass of the floating body. The backbone beam is simulated by three-dimensional Timoshenko beam elements, denoted by the red line in
Figure 9, which contributes to the stiffness, mass, and moment of the inertia of the floating body. The beam elements are connected to the shell elements by kinematic coupling constraints. The CFD solutions for the external fluid loads in the fluid domain are transferred to the FEA model in the structure domain through the hull surface for structural dynamic analysis. Once the structural dynamic analysis is done, the solutions for the motions and deformations of the backbone beam are obtained. Because of the kinematic coupling constraints between the beam elements and the shell elements, the motions and deformations of the backbone beam will also be fed back to the hull surface for the subsequent hydrodynamic analysis in the fluid domain.
The fluid domain modeling and mesh generation are carried out in STAR-CCM+. The wavelength is set to be
L, the length and width of the fluid domain are 4
L, the distance between the static water surface and the bottom of the domain is 2
L, and the distance between the static water surface and the top is 2
L. The boundary conditions are set [
8,
21] as shown in
Figure 10. A pressure outlet is applied at the top boundary. The nonslip wall boundary condition is applied at the floating body surface. Since the floating body is a symmetrical structure, and so is the flow about it in head seas, only half of the flow field needs to be simulated. The symmetry plane boundary condition is applied at the center plane of the structure. Regarding the other four boundary planes, the velocity inlet is applied, where the velocity of the wave is prescribed to avoid the gradient generated from the wall and flow [
21]. The volume of fluid (VOF) method is used to track free surfaces at the interface between air and water [
22]. The fluid domain is meshed as shown in
Figure 11. The extent of the whole computational domain is determined by referring to the recommendations of ITTC [
23]. The distance between the bow of the floating body and the inlet boundary is 1.5
L, and the distance between the stern of the floating body and the outlet boundary is also 1.5
L. The FVM is used to solve the fluid domain, and the fluid pressure and shear force on the surface of the floating body are obtained. The direction of the force due to the pressure is perpendicular to the surface of the floating body, and the direction of the shear force is tangent to the surface of the floating body.
The Courant number is used to evaluate the time step requirements of a transient simulation for a given mesh size and flow velocity.
where
is the Courant number,
is the flow velocity,
is the representative time step of the simulation, and
is the characteristic size of the mesh cell.
In order to ensure the accuracy and stability of the numerical results, the Courant number should be less than one [
24]. A constant time step of 0.01 s is chosen for the simulations.
The long-crest regular waves are created using linear airy wave theory. Since the offshore structure is in a moored condition during the operation, the forward speed is set to zero in the simulation. In head seas, the heave and pitch motions of the floating body are released, while the surge motion is constrained. Additionally, because of the symmetry, only half of the flow field is simulated, and a symmetrical boundary condition is applied at the center plane, that is, there are no sway, roll, and yaw motions.
3.2. Mesh Sensitivity Analysis
For the mesh convergence analysis, three sets of meshes, as shown in
Table 2, are generated for the structure domain and fluid domain mentioned above. A regular wave of 300 m in wave length and 16 m in wave height is used for the mesh sensitivity analysis.
The accuracy of the wave simulation is very important for the accuracy of the simulation of floating body’s motions and loads. According to the recommendation by STAR-CCM+, which has also been validated by many researchers [
25], a minimum of 40 cells per wavelength and 20 cells per wave height on the free surface is necessary to produce a stable wave with an acceptable dissipation. The mesh created according to the principle is Mesh B, as shown in
Table 2, and the mesh which is relatively coarser is Mesh A, while the mesh which is relatively coarser is Mesh C.
The numerical results for the amplitudes of the heave and pitch motion of the single floating body are presented in
Figure 12. The values of the heave and pitch motion are only rigid body displacements and do not include the elastic deformations of the floating body. For the heave motion, the results obtained with the three sets of meshes are very close to each other. For the pitch motion, the results obtained with Mesh B are very close to that obtained with Mesh C, with an average difference of 1.7%, but there is a relatively large difference, 21%, between the results obtained with Mesh A and Mesh B.
The numerical results for the amplitudes of the vertical bending moment at midship are presented in
Figure 13. Mesh B and Mesh C produce similar results, and the average difference between them is 1.6%. In comparison, there is a relatively large average difference, 35%, between the results obtained with Mesh A and Mesh B.
From the comparison of the results of three sets of meshes, it is found that the results from Mesh B and Mesh C are very close to each other, which shows the good mesh convergence property of the numerical model. Mesh B, with the meshing principle introduced above, is selected in this paper considering the accuracy and efficiency of the calculation.
3.3. The Influence of Shear Force along the Wetted Surface
Because the fluid is assumed to be inviscid in potential flow theory, the shear force acting on the wetted surface of the floating body is neglected in SESAM/WADAM [
26] and COMPASS-WALCS-NE [
27]. As described in
Section 2.5, the CFD-FEA method can account for the contribution of the shear force to the global solutions. In order to study the influence of shear force on the vertical bending moment of the cross-section, the numerical results obtained with and without shear force by the CFD-FEA method are compared, as shown in
Figure 14 and
Figure 15. From the comparison, it can be seen that the influence of shear force on the vertical bending moment of the cross-section is very small.
5. Conclusions
In this study, the hydrodynamic and structural responses of an unconventional large floating structure are simulated and analyzed using a CFD-FEA method. The results generated by the CFD-FEA method are compared with those generated by a rigid body linear frequency domain potential flow theory (WADAM) and by a three-dimensional potential flow nonlinear hydroelastic theory (COMPASS-WALCS-NE). Based on the analysis of the numerical results, the following concluding remarks can be drawn:
(1) The proposed CFD-FEA method can effectively exchange the information of the fluid domain and structure domain through the coupling interface, and a two-way coupling computation of the fluid domain and structure domain is realized. The numerical results of the elastic deformation of the floating body are reasonable.
(2) In the case of large wave heights where wave overtopping occurs, the rate of an increase in the amplitude of the vertical bending moment obtained by the CFD-FEA method slowly decreases with the increase in the wave height. The amplitude of the vertical bending moment obtained by the CFD-FEA method is smaller than that by WADAM. Hydroelastic responses of elastic floating structures and nonlinear behaviors are observed in large wave heights, for which case COMPASS-WALCS-NE and the CFD-FEA method are recommended to solve such problems.
(3) For large wave heights, the peak values of the vertical bending moment obtained by the CFD-FEA method are in a good agreement with those by COMPASS-WALCS-NE, which shows that the CFD-FEA method is valid. Moreover, the CFD-FEA method can also capture complex phenomena, such as the details of the change in the vertical bending moment , that the potential-flow-based methods cannot.
In this study, the CFD-FEA method for the hydroelastic response of large floating offshore structures is validated. Not only can the behavior of the hydroelastic responses of elastic floating structures be well estimated by the CFD-FEA method, but also the significant advantages of the CFD-FEA method over the traditional methods in capturing complex nonlinear behaviors are demonstrated.