1. Introduction
Offshore wind energy has become a globally recognized source of clean energy [
1,
2]. It boasts abundant wind resources and is not constrained by land limitations. Compared with onshore wind energy, offshore wind energy can effectively utilize larger-capacity wind turbines for power generation, resulting in higher wind energy utilization rates [
3,
4,
5]. It is projected that by 2028, the annual installed capacity of offshore wind energy may increase by 100% compared to that of 2023, thereby increasing its share of global added capacity from the current 9% to 20% by 2028. China and Europe are expected to continue leading this growth in 2024–2025 [
6].
With the expansion of the size and scale of offshore wind farms, the overall development potential of offshore wind power is limited. In sea areas with water depths exceeding 60 m, floating wind turbines have become a necessary option for building offshore wind farms. Among them, a semi-submersible foundation is one of the preferred structural forms of a floating wind turbine foundation [
7,
8,
9]. It has the advantages of a small draft, good stability, ease of transportation and installation, and relatively low cost. At present, among the wind power projects that have been put into practice in China, they have all adopted semi-submersible foundations. Therefore, the research object selected in this article was a semi-submersible floating wind turbine.
Starting in 2025, 25,000 tons of wind turbine blades are expected to be retired each year. The operation and maintenance cost of offshore floating wind power accounts for approximately 23% of the total life cycle cost of the entire project [
10,
11,
12]. Currently, if a wind turbine needs repair, the only options are to tow the wind turbine back to port for repair or to use a jack-up vessel for repair. However, as floating wind turbines increase in scale and reach deep into the sea, the jack-up installation vessels are no longer suitable for deep water environments, and towing the wind turbines back to port is time-consuming and costly. Therefore, both approaches are no longer applicable. For large floating wind turbines, crane vessels are indispensable equipment. Therefore, in this article, we selected a floating crane ship as a construction ship to study its blade replacement process for semi-submersible floating wind turbines. When facing large floating wind turbines under deep sea environmental conditions, the relative movement between a crane ship and a floating wind turbine will directly affect the safety and feasibility of the operation and maintenance process. Therefore, auxiliary measures need to be considered to make the operation and maintenance ship safer and more stable during the operation and maintenance process.
When replacing large components or performing maintenance operations in floating wind turbines, the floating crane and the floating wind turbine form a multi- body system. The surrounding flow field is different from that of a single floating body, and the hydrodynamic interaction is much more complex. Hydrodynamic analysis of multi-floating body systems mainly focuses on hydrodynamic resonance and shielding effects. Zou et al. [
13] studied the phenomenon of gap resonance between two adjacent fixed barges. A series of experiments and numerical simulations were conducted by varying the gap width, providing valuable insights for assessing the hydrodynamics and safety of offshore operations involving multiple floaters. Chen et al. [
14] studied the hydrodynamic interaction phenomena between a crane vessel and a barge arranged in both tandem and L-shaped configurations. Various hydrodynamic parameters were compared between these two arrangements, validating the observed hydrodynamic interference phenomenon. Based on the RMFC method, Zou et al. [
15] studied the coupled hydrodynamic interference problem by using a frequency–time-domain model with viscous correction. By using potential flow theory with viscous correction and introducing a state-space model to replace the convolution term in the Cummins equation, he achieved higher accuracy. In systems with more than one floating body, strong hydrodynamic interactions occur between the bodies when the system is subjected to wave loading, resulting in significant wave elevation phenomena.
When using potential flow theory for the hydrodynamic analysis of multiple floating bodies, the free surface response near resonance is often overpredicted [
16,
17]. This excessive hydrodynamic interaction can cause time-domain models to be difficult to converge. Therefore, corrections are needed. The main current method is to use artificial damping lids, that is, to introduce rigid boundary conditions or damping boundary conditions to the free surface boundary of the gap. Some scholars [
18,
19,
20] have used rigid lids to correct the slow drift forces on floating bodies at certain frequencies.
This method has been recognized by many scholars. Yao and Dong [
21] introduced an artificial damping coefficient for a free liquid surface based on three-dimensional potential flow theory to accurately simulate the hydrodynamic characteristics of multiple floating bodies. They found that, compared to other factors, the shielding effect of small-gap floating bodies was more significantly affected by transverse spacing. Li [
22] explored the hydrodynamic interactions in multi-body side-by-side models and referenced the use of artificial damping lids.
Tian et al. [
23], based on three-dimensional potential flow theory, investigated the influence of different configurations within wind wave combinations on the hydrodynamic performance of a semi-submersible platform in the frequency domain. The study examined the hydrodynamic interaction between the surrounding fluid and the semi-submersible platform, comparing results such as the Response Amplitude Operators (RAOs) to analyze the coupling effects. Degrieck et al. [
24] utilized the double-body potential flow code HYDINTER, developed in the early stages of this field of research by the University of Lisbon, Portugal, to investigate hydrodynamic ship–ship and ship–bank interactions, obtaining numerous numerical results on these interactions. Yuan [
25] developed a 3-D method based on the Rankine-type Green function to study the influence of parameters such as the oscillation frequency, forward velocity, and lateral distance between two ships on hydrodynamic interactions.
Compared to the jack-up installation vessels used in offshore wind turbine installation, floating crane vessels offer more flexible depth restrictions and significantly faster relocation speeds. They are expected to play a more prominent role in future component replacement tasks for floating wind turbines. However, wave motion poses challenges to crane operations. For crane vessels conducting component replacements, the motion of the crane tip is crucial as it directly affects the precision and safety of blade docking. Zhao et al. [
26] studied the motion response of two types of floating crane vessels, namely single-hull vessels and semi-submersible vessels, during single-blade hoisting processes and compared them with jack-up vessels. Sarker and Gudmestad [
27] Chen et al. [
28] developed a full coupled CPTDM model that can be used to consider all nonlinear problems in the float-over system and cross-validated it with AQWA. proposed a method for using floating crane vessels to install offshore wind turbine components and analyzed its feasibility. Zhu et al. [
29] compared the motion responses of jack-up crane vessels and floating crane vessels during tripod foundation lifting operations and analyzed the vertical motion of the crane tip. Ku and Roh [
30] investigated the dynamic responses of an offshore wind turbine and crane (tower–nacelle–rotor assembly) during lifting operations by a floating crane barge. They studied the motion responses of both during the lifting process of modules and analyzed the motion of the crane tip.
This study aims to investigate whether utilizing a lashing method for maintenance operations between floating crane vessels and floating wind turbines can enhance the safety and stability of the maintenance process. Initially, the frequency-domain hydrodynamic coefficients of both the crane vessel and floating wind turbine are computed using ANSYS-AQWA 2023R1 to explore hydrodynamic interference issues between them. Given the complex and transient nature of this multi-body system, comprehensive time-domain simulation calculations are necessary to obtain reliable results. Six operational conditions are selected for the time-domain simulation calculations of the berthing process between the floating crane vessel and floating wind turbine. Comparative analyses are then conducted on various aspects, including relative displacements, relative motions, and the displacements and motions of the floating crane vessel itself across these conditions. Ultimately, the feasibility of the lashing method is assessed based on the comprehensive analysis performed.
Section 2 introduces the theoretical basis used in this study.
Section 3 describes the crane vessel model, the floating wind turbine model, and the mooring line setup used in this study.
Section 4 and
Section 5 present the frequency-domain hydrodynamic calculation results and the time-domain calculation results, respectively.
Section 6 presents the conclusions of this paper.
3. Description of the Coupled System
The system studied in this paper comprises 12 degrees of freedom in total, with the floating crane and FOWT having six degrees of freedom each. The geometric model of the floating crane is depicted in
Figure 1, with the coordinate system defined in
Figure 2. The relevant parameters of the floating crane are listed in
Table 1. In
Figure 3, the FOWT is depicted. The origin of the global coordinate system is located at the waterline, above the FOWT’s center of gravity (COG). The local coordinate system of the floating crane is positioned at (−61 m, 0 m, 0 m), that is, a horizontal translation of the global coordinate system centered in the geometric center of the floating crane. In both reference frames, the Y-axes point aft and the Z-axes point upward.
Table 2 presents the COG positions for the FOWT and crane vessel.
The FOWT platform in this study is an adapted version of the OC4 DeepCWind semi-submersible for a more powerful wind turbine. Though this particular design has been developed for a 13.2 MW wind turbine [
39], this study employs the DTU-10 MW wind turbine, a three-bladed horizontal-axis wind turbine (HAWT), for it is publicly available and, therefore, can also be more easily validated and compared.
The FOWT consists of three main components: the wind turbine, the platform, and the mooring system. Since the platform design is based on scaling up of an OC4 platform, it is referred to hereafter as OC4+ [
39].
Table 3 and
Table 4 give details of the parameters of the platform and DTU-10 MW wind turbine.
The mooring system of the wind turbine consists of three anchor chains, with each pair of anchor chains forming a 120° angle. The fairleads are arranged at the base columns 21 m below the waterline, as shown in
Figure 3, and connected to the anchor points at a water depth of 200 m, as detailed in
Table 5 (OCML1 to OCML3 mooring lines).
The mooring system of the floating crane consists of eight anchor chains, as illustrated in
Figure 4. Each mooring chain is 800 m long, with the anchor points detailed in
Table 5 (QZCML1 to QZCML8 mooring lines). The connection setup between the crane vessels and floating wind turbine was referenced from the experiment conducted in [
40].
Simulations were performed in two different operational scenarios, namely, berthed and unberthed conditions. In the berthed conditions, two fenders and six lines constrained the relative motions between the OC4+ and the crane vessel for the maintenance operations. The berthed scenario is depicted in
Figure 5, with the lines defined in
Table 5 (XL1 to XL6) and the fenders defined in
Table 6. The stiffness of the XL lines used in berthing was defined as 5000 kN/m.
Last,
Figure 6 presents the fully coupled model in the Ansys AQWA interface with omitted mooring lines, so it is straightforward to check the relative position between the OC4+ FOWT and the crane vessel in the berthing condition for the purpose of maintenance operations.
4. Frequency-Domain Analysis
First, grid convergence analysis was conducted for the diffraction and radiation solver. Grid sizes of 1.6, 1.8, and 2.0 m were selected for this analysis. As shown in
Figure 7, hydrodynamic analyses of the crane vessel were conducted at a depth of 200 m and in the 180° wave direction using three different mesh sizes, resulting in comparative data across these mesh sizes.
Figure 7 presents the hydrodynamic results for the crane vessel, focusing on the added mass, radiation damping, and RAO for the heave, pitch, and roll degrees of freedom. The results for added mass and RAO are consistent across all three mesh sizes. However, while the results for radiation damping are consistent in the low-frequency range, minor discrepancies appear in the high-frequency range. Based on the mesh convergence analysis, the results for the 1.8 m and 1.6 m meshes are generally consistent. Considering computational efficiency and mesh limitations, a mesh size of 1.8 m was chosen for subsequent calculations.
In this study, the floating wind turbine was axisymmetric along its longitudinal axis. Wave directions of 0, 45, and 90 degrees were considered. The RAOs of the FOWT in single-body and multi-body scenarios with the crane vessel aside were calculated under those wave incidence angles and presented in
Figure 8. Clearly, by increasing the wave frequency, the RAOs for pitch, roll, and yaw exhibit a trend of initially increasing and then decreasing, while those for heave, sway, and surge gradually decrease. A comparison of the RAOs with and without the presence of the crane vessel (i.e., single- and multi-body conditions) indicates that the vessel has a noticeable impact only on the RAOs of pitch at wave angles of 45° and 90°, with minimal effect on other degrees of freedom. This is attributed to the shielding effect of the crane vessel on the floating wind turbine. The peak values of roll and pitch RAOs for the floating wind turbine occur approximately at the same wave frequency.
A comparison between
Figure 7 and
Figure 8 shows that the crane vessel and FOWT in single-body conditions present rather different responses to waves, which is expected from their qualitatively different geometries. Also, the crane vessel has much lower inertia, since its displacement is one order of magnitude below that of the OC4+. The results in
Figure 8 show that the barge-like geometry of the crane vessel makes it pitch considerably around the resonance frequency (ω
N ≅ 0.8 rad/s) and makes it heave with the waves for frequencies below the same value (or wave periods above 8 s). On the other hand, the OC4+ has a semi-submersible geometry with large dimensions, which makes it respond much less to waves and have relatively high natural periods, above 60 s for all vertical motions, namely, heave, roll and pitch. The effects of the presence of the crane vessel in the surroundings of the FOWT have practically no effect on the diffraction and radiation forces of the FOWT, as evident from
Figure 8. Significant differences are found only in the yaw motion, though it depends on wave incidence as well.
Figure 9 illustrates the results of added mass with six degrees of freedom for the crane vessel in single- and multi- body conditions. From this figure, it is evident that in the multi-body condition, there is no significant hydrodynamic effect of the presence of the FOWT on the surge’s added mass of the crane vessel. However, there is a considerable influence on the added mass in the heave, pitch, and yaw directions. Moreover, when the incident wave frequency exceeds 0.8 rad/s, there is significant interference in the added mass in the yaw directions.
Figure 10 illustrates the radiation damping results for the crane vessel with six degrees of freedom in single- and multi-body conditions. From this figure, it can be observed that there is no significant hydrodynamic effect from the presence of the FOWT on the radiation damping acting in the transverse sway direction. However, when the incident wave frequency exceeds 0.7 rad/s, there is significant interference in the radiation damping in the heave, pitch, and roll directions. Moreover, when the incident wave frequency is around 0.9 rad/s, the turbine has a significant impact on the radiation damping in the heave and yaw directions of the crane vessel.
Figure 11 depicts the RAOs of the crane vessel with six degrees of freedom, in three wave directions, as well as in single- and multi-body conditions. From this figure, it is observed that, with an increase in wave frequency, the RAO in the roll, pitch, and yaw directions shows a trend of initially increasing and then decreasing, while the RAO in the heave, sway, and surge directions demonstrates a gradual decrease. By comparing the RAO curves with and without the presence of the FOWT, it is evident that the yawing of the crane vessel is significantly influenced by the presence of the FOWT, whereas the impact on the other degrees of freedom is minimal.
5. Time-Domain Analysis
This section presents the multi-body time-domain simulations and analyses of crane vessel and FOWT dynamics under the combined action of wind and waves in berthed and unberthed conditions for maintenance purposes. Six conditions are selected to analyze the time-domain dynamics and compare systematically berthed and unberthed responses. Among them, CASE#3 and CASE#6 are extreme conditions. The simulation time is set to the recommended value for irregular sea states, namely, 3 h (10,800 s), though focus is placed on the time span of 6000–7000 s, when a steady state is already achieved. The calculation step is 0.1 s, and the result output step is also 0.1 s.
Considering the complex and variable nature of actual wind turbine operation, to better reflect real-world conditions and ensure applicability to Chinese offshore environments, the simulation conditions are designed with reference to the IEC 61400 standard [
41], regulations for offshore facilities operation during berthing, and relevant wind and wave statistical data for Chinese offshore regions. The selection of environmental conditions for the time-domain simulation is summarized in
Table 7, where the last column (mooring) stands for the crane vessel condition, whether it is moored to the FOWT or not.
First,
Figure 12 presents the time histories of the relative horizontal distances between the COGs for the six simulated cases, i.e., the difference between the X and Y positions of the crane vessel’s COG and the FOWT’s COG. In
Figure 13, the statistical parameters obtained from those time histories are given using bar graphs, which includes the maximum, minimum, and average values for each scenario.
A closer look into
Figure 12 and
Figure 13 makes it evident that the berthed condition results in smaller relative displacement between the center-of-gravity positions in both the X- and Y-directions after a steady state is achieved. That is due to the constrains of berthing. In unberthed conditions, the relative position in the Y-direction is very stable. On the other hand, the results in the X-direction show that the berthed condition can lead to a reduced distance between the floating bodies. Comparing the results of CASE#1 and #4 with CASE#2 and #5 reveals that berthing arrangements reduce the amplitude of the Y-direction distance between the crane vessel and the wind turbine by 60%. By comparing CASE#3 and #6, it is evident that under extreme conditions, berthed conditions lead to more significant variation in the lateral distance between the vessels.
Figure 12 depicts the motion response of the crane vessel in heave, roll, and pitch, in berthed and unberthed conditions. the Analysis reveals that the presence or absence of mooring lines has minimal impact on the vessel’s heave and pitch motions within operational environments. The roll motion is rather small in all scenarios, though the berthed condition tends to add energy to the roll motion due to coupling effects.
The time histories shown in
Figure 14 provide the basis for the statistical parameters presented in
Figure 15. These parameters are depicted using bar graphs, which include the maximum, minimum and average values for each scenario. From
Figure 15, it is evident that the berthing has minimal impact on the vessel’s motion response in terms of the heave degree of freedom, as there are minimal differences in all statistical parameters. However, the mooring configuration can reduce the surge motion response to some extent, particularly evident in unfavorable sea conditions, which is, of course, due to horizontal coupling forces arising from the mooring.
Figure 16 illustrates the mooring line tensions that ultimately act on the OC4+ under operational conditions and berthing conditions, namely, CASE#4, #5, and #6. This figure shows that under head waves, Cable2 and Cable3 experience similar forces, while Cable1 experiences the highest force. Due to the wind’s influence, the tension in Cable2 oscillates considerably around 6500 s in CASE#5, which is not an effect of the presence of the crane or berthing. The maximum tension in Cable2 for all conditions does not exceed 2.4 × 10
6 N. Observing the tension in Cable1 for the three conditions reveals that in windy conditions, the tension is more stable. In condition six, the maximum tension in Cable1 can reach 4.1 × 10
6 N, which is the maximum observed in all scenarios.
Figure 17 presents the time history of the displacement of the crane’s tip with three degrees of freedom. Overall, it can be observed that the berthing condition reduces the magnitude of crane tip displacement over time, particularly in harsh sea conditions. Notably, conditions three and six show significant differences. Around 6050 s, the X-displacement of the crane’s tip reaches 15 m without berthing, with a difference of up to 30 m between the maximum and minimum values. After berthing, this difference reduces to around 10 m. The simulations show that the crane tip movement in the X-, Y-, and Z-directions is reduced by over 30%. Large lateral displacements during blade docking are critical, as excessive lateral displacement can damage the blade root, leading to increased costs. In the Y-direction, berthing significantly reduces the variation in crane tip displacement, enhancing safety and saving time during blade docking.
6. Conclusions
In this study, a temporary berthing system between a crane vessel and a floating wind turbine is employed to better suit the maintenance operation in offshore locations. Based on ANSYS-AQWA software, coupled multi-body hydrodynamic analysis of the floating crane vessel and floating wind turbine was conducted and a systematic comparison was performed between berthed and unberthed conditions. Whereas convergence analysis and verification were performed in a frequency-domain analysis, which also gave the frequency-domain RAOs of the vessel and wind turbine in both single- and multi-body conditions, the time-domain analysis resulted comprehensive time histories of the fender forces, mooring forces, relative displacements, and velocities between the vessel and wind turbine, as well as cable tension in the wind turbine’s mooring system, and the displacement of the crane tip. However, the approach followed in this paper has some limitations. Only one type of connection method between the crane vessel and the floating wind turbine was selected. A comparison of the effects of various mooring methods should be conducted.
A summary of the conclusions drawn is given below:
1. The time-domain numerical results show that berthing the FOWT to the crane vessel reduces the relative displacements between floating bodies, with no effect on the relative surge velocity and a small impact on the relative sway velocity. On the other hand, the advantages of berthing appear in harsh sea conditions due to the reduced displacement, for it significantly reduces the motion and displacement of the crane vessel tip itself.
2. Mechanical and hydrodynamic interactions between the floating bodies in the berthed condition are observed. The hydrodynamic coupling leads to changes in the diffraction and radiation forces on the crane vessel in a multi-body configuration, though the impact on the RAOs is small, with the exception of the yaw RAO. The mechanical coupling is stronger in the surge, roll, and yaw directions, though it is definitely beneficial in surge, and roll is always limited in the simulated scenarios.
3. The adoption of the berthing method significantly reduces the magnitude of displacement in the crane tip’s position, facilitating more efficient maintenance operations on the floating wind turbine, particularly for the replacement of large components of the wind turbine. It also enhances safety during blade docking and saves a substantial amount of time, which is a major advantage given the limited weather windows these operations are subject to. This method may serve as a valuable reference for future maintenance operations of floating wind turbines.
Future research may explore the combined use of tug assistance and different mooring arrangements to further improve the safety and effectiveness of maintenance operations in berthing conditions, also simulating the exact marine operation of replacing large components of floating wind turbines. The time-domain simulations in this study were computationally intensive, using convolution integrals based on the Cummins equation. In future studies, more efficient time-domain models, such as those based on the state-space model, could be utilized to analyze the dynamic behavior of floating crane vessels during maintenance operations with a substantial improvement in time-efficiency