Dynamic Responses for WindFloat Floating Offshore Wind Turbine at Intermediate Water Depth Based on Local Conditions in China

: Offshore wind energy, a clean energy resource, is considered to be a possible alternative to fossil energy. Floating offshore wind technology is considered to be a proper concept to develop abundant wind energy in deep water. Considering the reality of offshore wind energy development in China, the ﬂoating offshore wind turbine concept is expected to be developed at moderate water depths. In this paper, a mooring system of the WindFloat semisubmersible ﬂoating offshore wind turbine (SFOWT) at a water depth of 60 m is designed. The dynamic responses of the WindFloat SFOWT under different wind–wave combination conditions are investigated using the coupled method and the simpliﬁed method, which do not include the effect of the tower top motion in the aerodynamic calculation. The results show that the dynamic responses of the WindFloat SFOWT, including the platform motions, tower loads, and mooring line tensions, perform fairly well at a moderate water depth. A comparison between the coupled method and simpliﬁed method shows that the calculated results are slightly different between the different conditions for the time domain results, response spectra results, and fatigue results. In addition, mooring line 1 (ML 1) suffers higher fatigue damage than ML2, which should be paid more attention. tension responses. The results show that the tower base force and mooring line tension have similar STD values under the extreme sea conditions (LC 2 and LC 3), while the STD value of the simpliﬁed method is slightly lower than those of the coupled method for the over-rated wind speed condition (LC 1). The means and maximums of the tower base shear force and moment are larger in the calculation of the


Introduction
Fossil energy shortages and environmental pollution caused by the consumption of fossil fuels have become one of the most critical problems around the world [1,2]. Great efforts have been made to pursue clean and sustainable energy. Offshore wind energy, as an excellent substitute energy of fossil energy, is trying to move from shallow water to deep water to exploit the abundant wind energy far from the shore [3]. China is now the world's third largest offshore market in total installation as of the end of 2019 (after the UK and Germany) [4]. As offshore wind technology moves to deeper water (over 50 m water depth), floating foundations are expected to be an active research area. Floating offshore wind turbine foundations can be divided into three types: semisubmersible [5,6], tension leg platforms [7,8], and spar [9,10]. Among all floating types, semisubmersible floating foundations are considered to be very efficient design configurations due to their low construction cost and widely applicable range of water depths [11]. Considering the energy distribution and technical reality, the floating offshore wind technology is expected to be developed mainly at moderate water depths with the application of SFOWTs.
A variety of numerical simulation and tank tests have been conducted for investigating the dynamic responses of floating an offshore wind turbine with combined wind-wave interaction. Based on the designed Windfloat SFOWT, Roddier et al. [12] investigated the dynamic responses of WindFloat SFOWT, which was equipped with a National Renewable Energy Laboratory (NREL) 5 MW offshore wind turbine at 200 m water depth through an experimental approach. Comparison of WindFloat SFOWT is made between the experiment and the numerical simulation via TimeFloat-Fast software. The response amplitude operators (RAOs) in the surge, heave, and pitch show good agreement between the numerical results and experimental results. To investigate the damping effect of the heave plate on the SFOWT, Moreno et al. [13] carried out a model test based on WindFloat SFOWT. The hexagonal and circular structure of the heave plates is explored. The results show that different shapes of the heave plates show similar trends in the hydrodynamic performance, and the maximum difference of the damping coefficient is 8%. Based on the DeepCwind project in the US [14], the University of Maine (UoM) and NREL proposed a floating offshore wind turbine concept named OC4-DeepCwind SFOWT. Additionally, the OC4-DeepCwind SFOWT was selected as the reference model for IEA Wind Task 30 (OC4 to OC6), which was used for the validation of the response quantities of the system by a comparison with the experimental data. The numerical simulations of Braceless SFOWT by Luan et al. [15] were carried out to analyze its intact stability, natural period, and dynamic global responses. Moreover, a scale model of the Braceless SFOWT has been investigated in a hybrid model test by SINTEF Ocean (Marintek).
Several papers have discussed the utilization of floating offshore wind turbines at moderate water depths. Xu et al. [16] designed the mooring systems of Braceless SFOWT at 50 m water depth. A second order wave force is considered, and a fully coupled method was used in the study to investigate the response performance of Braceless SFOWT under parked conditions. The results show that it is feasible to mount the semisubmersible floating offshore wind turbine at shallow water depths and that the second order wave force has great impacts on the dynamic responses. Zhang et al. [17] compared different floating offshore wind turbine configurations at moderate and deep water depths. The results revealed that second-order difference wave loads can excite the resonance of motion, especially for the pitch motion. Bayati et al. [18] analyzed the hydrodynamic performance of OC4-semisubmersible offshore wind turbine from deep water depths to shallow water depths, illustrating that the water depth has a great impact on the low frequency motion of the platform. Zhang et al. [19] designed a semisubmersible floating offshore wind foundation at a water depth of 60 m. Combined wind and wave environments were used to investigate the dynamic responses of the SFOWT. The motions and mooring line tensions show good performance in the ocean environment.
In this paper, the WindFloat SFOWT is used to investigate the feasibility of floating offshore wind turbines at moderate water depths. Most of papers discuss the dynamic responses of WindFloat SFOWT at deep water depth. Considering the reality of floating offshore wind turbine development in China, it is very valuable to investigate the dynamic response of WindFloat SFOWT at shallow water depths. Therefore, mooring system design for WindFloat SFOWT at 60 m water depth is carried out in this paper. A comparison between the coupled analysis and decoupled analysis is made to investigate the dynamic responses of WindFloat SFOWT at a moderate water depth. This paper may help resolve the fundamental design trade-offs between the floating-system concepts in areas with moderate water depths. In the present simulation, blade element momentum theory (BEMT) [20] is used to derive the aerodynamic load on the blade. The blade is divided into a set of elements, and the blade elements are assumed to be independent from each other. As shown in Figure 1, the aerodynamic loads are calculated by evaluating the relative inflow wind speed. The relative wind speed can be calculated by the Equation (1):

Theory Background
where ρ is the air density, c is the chord length, and C L and C D are the lift coefficient and drag coefficient of each air foil, respectively. The NREL 5 MW wind turbine of the lift coefficient and drag coefficient can be found in [21]. It should be noted that the unsteady effects are not considered in Equation (2), which are illustrated in [22,23].
where ρ is the air density, c is the chord length, and L C and D C are the lift coefficient and drag coefficient of each air foil, respectively. The NREL 5 MW wind turbine of the lift coefficient and drag coefficient can be found in [21]. It should be noted that the unsteady effects are not considered in Equation (2), which are illustrated in [22,23].

Aerodynamic Calculation by Coupled Method and Simplified Method
For the calculation of the rotor thrust force [24], the coupled method ( Figure 2a) considers the tower top motion. When the wind turbine moves along the direction of the wind, the aerodynamic loads decrease. For the one-dimensional model, the rotor thrust thrust F with the assumption of an ideal rotor can be calculated by:

Aerodynamic Calculation by Coupled Method and Simplified Method
For the calculation of the rotor thrust force [24], the coupled method (Figure 2a) considers the tower top motion. When the wind turbine moves along the direction of the wind, the aerodynamic loads decrease.
where ρ is the air density, c is the chord length, and L C and D C are the lift coefficient and drag coefficient of each air foil, respectively. The NREL 5 MW wind turbine of the lift coefficient and drag coefficient can be found in [21]. It should be noted that the unsteady effects are not considered in Equation (2), which are illustrated in [22,23].

Aerodynamic Calculation by Coupled Method and Simplified Method
For the calculation of the rotor thrust force [24], the coupled method ( Figure 2a) considers the tower top motion. When the wind turbine moves along the direction of the wind, the aerodynamic loads decrease. For the one-dimensional model, the rotor thrust thrust F with the assumption of an ideal rotor can be calculated by: For the one-dimensional model, the rotor thrust F thrust with the assumption of an ideal rotor can be calculated by: 4 of 17 where C T is the thrust coefficient, ρ is the air density, A is the rotor disc area, and V is the wind speed. The wind speed is separated into two parts, as shown below: where V mean is the mean wind speed and V turb is the turbulent wind speed. Considering the tower top motion of the wind turbine, the relative wind speed V rel can be obtained by subtracting the velocity of the tower top .
x t : Then, the rotor thrust force F total can be divided into the mean thrust force F mean , the turbulent thrust force F turb and the aerodynamic damping force F damping : For the simplified method ( Figure 2b) used in the present paper, the precomputed aerodynamic load from a rigid wind turbine model includes mean and turbulent part. The aerodynamic damping caused by the interaction between the wind speed and tower top motion is not considered. Firstly, the precomputed rotor thrust force is calculated from a rigid wind turbine. Secondly, the rotor thrust force is added at the hub height for the simplified method.

Morison Equation
For the hydrodynamic calculation of braces, the Morison equation [25] is used. The form of the Morison equation accounts for the wave-induced excitation load, added mass and, drag force. The wave loads are calculated by using the Morison equation with the following formula: ..
x 1 represent the floating body motion, velocity and acceleration, respectively; C D is the drag coefficient; C M is the inertia coefficient; u and . u are the velocity and acceleration of water particles, respectively; D is the diameter of cylinder and dz is the length of a strip.

Potential Flow Theory
For a large volume floating structure, the Morison equation is not capable of calculating the hydrodynamic loads. The potential flow theory [26], assuming that the fluid is inviscid, incompressible, and irrotational, is used to calculate the hydrodynamic loads of the large floating structures. Three components of hydrodynamic loads are accounted for: hydrostatic force, radiation wave force, and wave-induced excitation force. ..
x k k = 1, 2 . . . 6 j = 1, 2 . . . 6 (10) → F r is the radiation wave load; ω is the wave frequency; → n is the normal direction vector of the wet surface; s is the wet surface area;φ i is the incident potential of the wave; φ d is the diffraction potential of the wave; ρ is the density of the water; φ Re i and φ Im i are the real and imaginary parts of the incident potential of the wave, respectively; A jk and B jk are the added mass and radiation damping coefficients, respectively; and the indices k and j refer to the degrees of freedom (DOFs) of the platform.

Mooring System Theory
We use the lumped-mass model [27] (Figure 3) to calculate the mooring system connected to the floating platform. For the lumped-mass method, the mooring line is divided into a set of segments, which are represented by the connected nodes and spring-damper systems. Each segment is divided into two components, and the properties are assigned and lumped to the two nodes at each end of that segment. The connections between adjacent nodes are represented by damper-spring systems. Only the axial properties of the mooring lines are accounted for, whereas the torsional and bending properties are not considered. are the added mass and radiation damping coefficients, respectively; and the indices k and j refer to the degrees of freedom (DOFs) of the platform.

Mooring System Theory
We use the lumped-mass model [27] (Figure 3) to calculate the mooring system connected to the floating platform. For the lumped-mass method, the mooring line is divided into a set of segments, which are represented by the connected nodes and spring-damper systems. Each segment is divided into two components, and the properties are assigned and lumped to the two nodes at each end of that segment. The connections between adjacent nodes are represented by damper-spring systems. Only the axial properties of the mooring lines are accounted for, whereas the torsional and bending properties are not considered.

Fatigue Damage Calculation
In this paper, the tower base section is simplified as a thin-walled cylinder structure that does not consider the welding effect. Fatigue damage includes shear fatigue and axial fatigue, while shear fatigue is less than axial fatigue. Therefore, shear fatigue is not considered in this paper. The procedure of the short-term fatigue damage analysis for the tower base and mooring lines is displayed in Figure 4. The Mlife code [28] developed by the NREL is used to calculate the fatigue of the tower base and mooring lines.

Fatigue Damage Calculation
In this paper, the tower base section is simplified as a thin-walled cylinder structure that does not consider the welding effect. Fatigue damage includes shear fatigue and axial fatigue, while shear fatigue is less than axial fatigue. Therefore, shear fatigue is not considered in this paper. The procedure of the short-term fatigue damage analysis for the tower base and mooring lines is displayed in Figure 4. The Mlife code [28] developed by the NREL is used to calculate the fatigue of the tower base and mooring lines.
 is the radiation wave load; ω is the wave frequency; n  is the normal direction vector of the wet surface; s is the wet surface area; i φ is the incident potential of the wave; d φ is the diffraction potential of the wave; ρ is the density of the water; Re i φ and Im i φ are the real and imaginary parts of the incident potential of the wave, respectively; jk A and jk B are the added mass and radiation damping coefficients, respectively; and the indices k and j refer to the degrees of freedom (DOFs) of the platform.

Mooring System Theory
We use the lumped-mass model [27] (Figure 3) to calculate the mooring system connected to the floating platform. For the lumped-mass method, the mooring line is divided into a set of segments, which are represented by the connected nodes and spring-damper systems. Each segment is divided into two components, and the properties are assigned and lumped to the two nodes at each end of that segment. The connections between adjacent nodes are represented by damper-spring systems. Only the axial properties of the mooring lines are accounted for, whereas the torsional and bending properties are not considered.

Fatigue Damage Calculation
In this paper, the tower base section is simplified as a thin-walled cylinder structure that does not consider the welding effect. Fatigue damage includes shear fatigue and axial fatigue, while shear fatigue is less than axial fatigue. Therefore, shear fatigue is not considered in this paper. The procedure of the short-term fatigue damage analysis for the tower base and mooring lines is displayed in Figure 4. The Mlife code [28] developed by the NREL is used to calculate the fatigue of the tower base and mooring lines. According to [29], the tower base axial stress can be calculated by the following equation: where σ is the nominal stress axial stress; A is the area of the cross section; F z , M x and M y are the axial force, moment about the x axis, and moment about the y axis, respectively; I x and I y are the sections of the area to x axis, y axis, respectively; and θ is the angle from the point to the calculated point ( Figure 5), Point 1 (P1) to point 12 (P12) are used to defined the location peripheral direction where the stress was calculated.
According to [29], the tower base axial stress can be calculated by the following equation: where σ is the nominal stress axial stress; A is the area of the cross section; z F ,  To evaluate the mooring line fatigue damage [30], the nominal stress of the mooring line can be calculated by the following equation: where mooring F is the mooring line force.
Then, the rainflow method is used to process the time series results of the axial stress analysis. After that, the fatigue damage can be obtained using the S-N curve. Based on the S-N curve, the predicted number of cycles to failure for the stress range i N can be given as: where ult L is the ultimate stress; mf L is the fixed mean stress value; m is the negative inverse slope of S-N curve; m is defined as 3 in the present paper; and , i RF L is the stress value corresponding to the number of cycles.

Model Description
In this paper, the NREL 5 MW wind turbine mounted on the WindFloat semisubmersible floating foundation is modeled and simulated, as shown in Figure 6. To evaluate the mooring line fatigue damage [30], the nominal stress of the mooring line can be calculated by the following equation: where F mooring is the mooring line force. Then, the rainflow method is used to process the time series results of the axial stress analysis. After that, the fatigue damage can be obtained using the S-N curve. Based on the S-N curve, the predicted number of cycles to failure for the stress range N i can be given as: where L ult is the ultimate stress; L m f is the fixed mean stress value; m is the negative inverse slope of S-N curve; m is defined as 3 in the present paper; and L i,RF is the stress value corresponding to the number of cycles.

Model Description
In this paper, the NREL 5 MW wind turbine mounted on the WindFloat semisubmersible floating foundation is modeled and simulated, as shown in Figure 6.

Wind Turbine Model
The wind turbine model [30] was developed by the National Renewable Energy Laboratory (NREL), USA, which is a conventional three-bladed, upwind, variable-speed, collective-pitch controlled horizontal axis wind turbine. Its main properties are listed in

Wind Turbine Model
The wind turbine model [30] was developed by the National Renewable Energy Laboratory (NREL), USA, which is a conventional three-bladed, upwind, variable-speed, collective-pitch controlled horizontal axis wind turbine. Its main properties are listed in Table 1.

WindFloat Floating Platform Properties
In this paper, the WindFloat semisubmersible platform was adopted to support the NREL 5-MW reference wind turbine. The platform consists of one center column and two side columns connected by braces. It should be noted that the active ballast method was used for the WindFloat semisubmersible platform. However, this was not considered in the present simulation. It should be noted that all the parameters are given about the reference system at the mean sea level (MSL). The properties of the WindFloat semisubmersible floating platform are summarized in Table 2.

Mooring System Design
Considering the geographical and technical strategy in China, the mooring system of WindFloat SFOWT at 60 m water depth is designed (Figure 7). The objective of the design of the mooring system at moderate water depth is to achieve the reasonable natural frequency of the structure, pretension, and the safety of the mooring system. For each mooring line, a clump mass is used. Compared to the original mooring system of the WindFloat SFOWT at 200 m water depth, a heavier mooring line and clump are used to obtain a reasonable mooring tension.

Mooring System Design
Considering the geographical and technical strategy in China, the mooring system of WindFloat SFOWT at 60 m water depth is designed (Figure 7). The objective of the design of the mooring system at moderate water depth is to achieve the reasonable natural frequency of the structure, pretension, and the safety of the mooring system. For each mooring line, a clump mass is used. Compared to the original mooring system of the WindFloat SFOWT at 200 m water depth, a heavier mooring line and clump are used to obtain a reasonable mooring tension.  Table 3. The static configuration and axial tension for each mooring line element are presented in Figure 8. The results show that the mooring line tension at the fairlead position is the largest, which is 1400 kN in the static equilibrium position and a large drop in the mooring line tension at the clump position. Figure 9 presents the mooring line force and restoring force with offset in surge. As shown in Figures 7 and 8, the characteristics of the nonlinearity appear in the mooring line force and restoring force with offset in surge.  The fairlead, anchoring positions, and mooring line properties are listed in Table 3. The static configuration and axial tension for each mooring line element are presented in Figure 8. The results show that the mooring line tension at the fairlead position is the largest, which is 1400 kN in the static equilibrium position and a large drop in the mooring line tension at the clump position. Figure 9 presents the mooring line force and restoring force with offset in surge. As shown in Figures 7 and 8, the characteristics of the nonlinearity appear in the mooring line force and restoring force with offset in surge.

Mooring System Design
Considering the geographical and technical strategy in China, the mooring system of WindFloat SFOWT at 60 m water depth is designed (Figure 7). The objective of the design of the mooring system at moderate water depth is to achieve the reasonable natural frequency of the structure, pretension, and the safety of the mooring system. For each mooring line, a clump mass is used. Compared to the original mooring system of the WindFloat SFOWT at 200 m water depth, a heavier mooring line and clump are used to obtain a reasonable mooring tension.  Table 3. The static configuration and axial tension for each mooring line element are presented in Figure 8. The results show that the mooring line tension at the fairlead position is the largest, which is 1400 kN in the static equilibrium position and a large drop in the mooring line tension at the clump position. Figure 9 presents the mooring line force and restoring force with offset in surge. As shown in Figures 7 and 8, the characteristics of the nonlinearity appear in the mooring line force and restoring force with offset in surge.

Load Case Definition
Considering the metocean conditions in the South China Sea, typical load cases (LCs) at 60 m water depth, including the operating conditions and parked conditions, are included in the present simulation. The wind and wave are collinear along the positive X axis. The misaligned wind-wave conditions are not considered in the simulation. It should be noted that turbulent wind is used in this paper. The turbulent intensity (TI) is 0.12. The load cases are listed in Table 4.

Natural frequencies of WindFloat SFOWT
The natural frequencies of the WindFloat SFOWT for different modes are obtained by free decay test simulation in SIMA by giving an initial displacement to the platform. Then, the time series of the free decay for each DOF can be obtained for the output results. The natural frequency for each DOF can be calculated based on the fast Fourier transform (FFT). The natural frequencies of the WindFloat SFOWT for all DOFs are shown in Table 5.

Comparison of Dynamic Responses Using Coupled and Simplified Method
In this section, the dynamic response of the coupled method and simplified method are compared and investigated. For the simplified method, the wind force is considered as a time series signal. In the present simulation, the dynamic damping force is not yet considered. The other two components (mean part and the dynamic turbulent part), can be calculated from the monopole offshore wind turbine with A rigid foundation. The total simulation time is 4600 s, and, for the calculation of the statistical results, the first 1000 s has been excluded to avoid the transient effect. It should be noted that the time series results of the motion and force responses between 2000 and 2500 s are displayed to better see the difference between the coupled method and simplified method.  Table 6 shows a comparison of statistics, including the maximum, minimum, average, and standard deviation (STD) values in surge, heave, and pitch motion. The motion statistics of the WindFloat SFOWT at moderate water depth are quite acceptable for both normal operating condition (LC 1) and extreme sea conditions (LC 2 and LC 3). For surge motion, the amplitude is 9.07 m under normal operation, while it is 20.86 m under extreme condition (LC 3). For the heave and pitch motion, their amplitudes are 3.5 m and 4 • , respectively. With respect to the difference of the motion response between the coupled method and the simplified method, the calculated response of the simplified method is larger than that of the coupled method. For example, under normal operational condition, the surge and pitch motion maxima are increased by 11% and 16%, respectively. It seems that the STDs between the different motions are quite close. For example, the STDs of the surge motion of the two methods are 1.09 and 1.10 under normal operating condition (LC 1).  Table 7 shows the statistical comparison of tower base loads and mooring line tensions for different load cases. It should be noted that only the tension forces of ML 1 and 2 are displayed in the present paper, since the ML 2 and 3 are in the upwind condition and show the similar mooring tension responses. The results show that the tower base force and mooring line tension have similar STD values under the extreme sea conditions (LC 2 and LC 3), while the STD value of the simplified method is slightly lower than those of the coupled method for the over-rated wind speed condition (LC 1). The means and maximums of the tower base shear force and moment are larger in the calculation of the simplified method. For the mooring line tension, the forces acting on the hub under the parked conditions are lower than those under the normal operating conditions. The mean surge motion shown in Table 6 is lower under the parked condition. Therefore, the ML 1 tension (downwind direction) is expected to be larger in extreme sea conditions, while the ML 2 tension is not.  Table 7 shows the statistical comparison of tower base loads and mooring line tensions for different load cases. It should be noted that only the tension forces of ML 1 and 2 are displayed in the present paper, since the ML 2 and 3 are in the upwind condition and show the similar mooring tension responses. The results show that the tower base force and mooring line tension have similar STD values under the extreme sea conditions (LC 2 and LC 3), while the STD value of the simplified method is slightly lower than those of the coupled method for the over-rated wind speed condition (LC 1). The means and maximums of the tower base shear force and moment are larger in the calculation of the simplified method. For the mooring line tension, the forces acting on the hub under the parked conditions are lower than those under the normal operating conditions. The mean surge motion shown in Table 6 is lower under the parked condition. Therefore, the ML 1 tension (downwind direction) is expected to be larger in extreme sea conditions, while the ML 2 tension is not.  Table 7 shows the statistical comparison of tower base loads and mooring line tensions for different load cases. It should be noted that only the tension forces of ML 1 and 2 are displayed in the present paper, since the ML 2 and 3 are in the upwind condition and show the similar mooring tension responses. The results show that the tower base force and mooring line tension have similar STD values under the extreme sea conditions (LC 2 and LC 3), while the STD value of the simplified method is slightly lower than those of the coupled method for the over-rated wind speed condition (LC 1). The means and maximums of the tower base shear force and moment are larger in the calculation of the simplified method. For the mooring line tension, the forces acting on the hub under the parked conditions are lower than those under the normal operating conditions. The mean surge motion shown in Table 6 is lower under the parked condition. Therefore, the ML 1 tension (downwind direction) is expected to be larger in extreme sea conditions, while the ML 2 tension is not.

Comparison of Frequency-Domain Response
In this section, a frequency-domain comparison of motions and mooring line tensions is made between the coupled method and simplified method. To investigate the contribution from different frequency regions, the integrals of the power spectral density (PSD) plots at low frequency and wave frequency regions are displayed. It should be noted that the first 1000 s of the total simulation time (4600 s) was removed for the frequency domain analysis to avoid the transient effect. Figure 13 shows the PSD plots of the surge, heave, and pitch motion of the WindFloat SFOWT under normal operation (LC 1) and extreme sea conditions (LC 2 and LC 3). For the surge motion, the wave frequency and natural frequency of the surge motion are indicated in the surge motion PSD plot. The low frequency response contribution is larger than the wave frequency response under both the normal operating condition and extreme sea conditions, as shown in Figure 13b. Compared to the normal operation condition (LC 1), the contribution to surge motion from the wave excitation force increases. Additionally, a comparison between the coupled and simplified methods indicates that the simplified method has a relatively similar frequency response in both the low frequency and wave frequency areas. Being different with the surge motion response, the heave motion response (Figure 13c, d) is mainly dominated by the wave excitation force.

Comparison of Frequency-Domain Response
In this section, a frequency-domain comparison of motions and mooring line tensions is made between the coupled method and simplified method. To investigate the contribution from different frequency regions, the integrals of the power spectral density (PSD) plots at low frequency and wave frequency regions are displayed. It should be noted that the first 1000 s of the total simulation time (4600 s) was removed for the frequency domain analysis to avoid the transient effect. Figure 13 shows the PSD plots of the surge, heave, and pitch motion of the WindFloat SFOWT under normal operation (LC 1) and extreme sea conditions (LC 2 and LC 3). For the surge motion, the wave frequency and natural frequency of the surge motion are indicated in the surge motion PSD plot. The low frequency response contribution is larger than the wave frequency response under both the normal operating condition and extreme sea conditions, as shown in Figure 13b. Compared to the normal operation condition (LC 1), the contribution to surge motion from the wave excitation force increases. Additionally, a comparison between the coupled and simplified methods indicates that the simplified method has a relatively similar frequency response in both the low frequency and wave frequency areas. Being different with the surge motion response, the heave motion response (Figure 13c,d) is mainly dominated by the wave excitation force. Compared to the coupled method, the low frequency response is slightly larger for the simplified method under normal operating conditions, while the wave frequency response is larger under extreme sea conditions. The response spectra of the pitch motion shows that the pitch natural frequency response has a larger contribution than the wave excitation frequency in response to the pitch motion under the normal operating condition, while the motion response is mainly dominated by the wave excitation force under the parked conditions. For both normal operating conditions and extreme sea conditions, small peaks at the surge natural frequency are observed in the pitch motion response spectra, showing a strong coupling effect between the surge and pitch motion. Under the normal operation, the low frequency response of the pitch motion using the simplified method is larger than that using the coupled method.
J. Mar. Sci. Eng. 2021, 9, 1093 13 of 17 Compared to the coupled method, the low frequency response is slightly larger for the simplified method under normal operating conditions, while the wave frequency response is larger under extreme sea conditions. The response spectra of the pitch motion shows that the pitch natural frequency response has a larger contribution than the wave excitation frequency in response to the pitch motion under the normal operating condition, while the motion response is mainly dominated by the wave excitation force under the parked conditions. For both normal operating conditions and extreme sea conditions, small peaks at the surge natural frequency are observed in the pitch motion response spectra, showing a strong coupling effect between the surge and pitch motion. Under the normal operation, the low frequency response of the pitch motion using the simplified method is larger than that using the coupled method. As mentioned in the time domain analysis, ML 2 and ML 3 share similar responses due to the symmetry of the mooring configuration. Therefore, only the ML 1 and ML 2 tension response spectra are displayed in the present section. Figure 14 presents a comparison of tension spectra for ML 1 and ML 2 between the two methods, including the As mentioned in the time domain analysis, ML 2 and ML 3 share similar responses due to the symmetry of the mooring configuration. Therefore, only the ML 1 and ML 2 tension response spectra are displayed in the present section. Figure 14 presents a comparison of tension spectra for ML 1 and ML 2 between the two methods, including the coupled and simplified method under different wind-wave combination cases. The consequences reveal that the mooring forces are mainly dominated by responses at the surge natural frequency for the normal operation. Wave actions also contribute to the mooring tension responses. When the WindFloat SFOWT operates in the extreme sea conditions, the contribution of the ML 1 tension caused by the wave excitation force increases. Compared with the ML 1 tension response, it seems that the low frequency has a much larger impact on the mooring line tension response. A comparison of the mooring line tension between the two methods indicates that, the calculated results of the mooring line tension response using the simplified are lower than those using the coupled method. For ML 2 tension response spectra, the calculated response using the coupled method is slightly larger under the normal operating conditions, while the low frequency response of the mooring line is lower than those using the simplified method. A comparison of the mooring line tension responses shows that the different methods have an impact on the mooring tension performance. Therefore, fatigue damage estimation is needed to further check the difference between the two methods.
J. Mar. Sci. Eng. 2021, 9, 1093 14 of 17 coupled and simplified method under different wind-wave combination cases. The consequences reveal that the mooring forces are mainly dominated by responses at the surge natural frequency for the normal operation. Wave actions also contribute to the mooring tension responses. When the WindFloat SFOWT operates in the extreme sea conditions, the contribution of the ML 1 tension caused by the wave excitation force increases. Compared with the ML 1 tension response, it seems that the low frequency has a much larger impact on the mooring line tension response. A comparison of the mooring line tension between the two methods indicates that, the calculated results of the mooring line tension response using the simplified are lower than those using the coupled method. For ML 2 tension response spectra, the calculated response using the coupled method is slightly larger under the normal operating conditions, while the low frequency response of the mooring line is lower than those using the simplified method. A comparison of the mooring line tension responses shows that the different methods have an impact on the mooring tension performance. Therefore, fatigue damage estimation is needed to further check the difference between the two methods.

Comparison of Fatigue Damage
As shown in Figure 15, judging by the axial stress and damage rate of the tower base position at the different points (P1 to P12) shown in Figure 5, the points along the windwave direction have the maximum fatigue damage, while the points perpendicular to the wind-wave direction have the lowest fatigue damage. Meanwhile, P7 experiences greater fatigue damage compared to P1. When the floating wind turbine operates in the parked condition, the fatigue damage rate difference between P1 and P7 fatigue damage rate decreases. Compared to the normal operating condition (LC 2), the fatigue load at the tower base is larger under the extreme conditions (LC 2 and LC 3). According to the comparison between the coupled method and simplified method under the normal operating condition (LC 1), there is a large difference between the P1 and P7, while they are smaller under extreme sea conditions (LC 2 and LC 3). For example, the damage rates calculated based

Comparison of Fatigue Damage
As shown in Figure 15, judging by the axial stress and damage rate of the tower base position at the different points (P1 to P12) shown in Figure 5, the points along the wind-wave direction have the maximum fatigue damage, while the points perpendicular to the wind-wave direction have the lowest fatigue damage. Meanwhile, P7 experiences greater fatigue damage compared to P1. When the floating wind turbine operates in the parked condition, the fatigue damage rate difference between P1 and P7 fatigue damage rate decreases. Compared to the normal operating condition (LC 2), the fatigue load at the tower base is larger under the extreme conditions (LC 2 and LC 3). According to the comparison between the coupled method and simplified method under the normal operating condition (LC 1), there is a large difference between the P1 and P7, while they are smaller under extreme sea conditions (LC 2 and LC 3). For example, the damage rates calculated based on the coupled method and decoupled method of P1 at tower base are 1.64 × 10 −6 and 1.1710 −6 under the normal operating condition (LC 1), while the damage rates are 2.18 × 10 −6 and 2.10 × 10 −6 under extreme sea condition (LC 2), respectively. This is because the wave excitation force dominates the tower base responses under the extreme sea condition.  Figure 16 presents the axial stress and damage rate of the mooring system. The results show that the mooring system experiences a smaller fatigue load under normal operating condition than that under extreme sea conditions (LC 2 and LC 3). According to the time domain response of ML 2, although the mean value of ML 2 is larger than that of ML 1, the STD value of the mooring line tension is larger for ML 1, which has a significant impact on fatigue damage. A comparison between ML 1 and ML 2 indicates that ML 1 suffers more fatigue damage than ML 2, which should be paid more attention. Moreover, under extreme sea conditions (LC 2 and LC 3), the mooring lines experience more fatigue damage than under normal conditions. The time series statistics in Table 7 show that the mean value and STD value of ML 1 are larger using the coupled method, while those of ML 2 are larger using the simplified method. The comparison of fatigue damage between the coupled method and simplified method shows that the fatigue damage using the coupled method is larger than that using the simplified method for ML 1, while the fatigue of ML 2 is the opposite.

Conclusions
In this paper, the dynamic response analysis of the WindFloat SFOWT at a moderate water depth was carried out. A comparison of the dynamic responses between the coupled method and simplified method is made. The coupled method considers the motion interaction between the wind turbine and tower top and the aerodynamic damping, while the simplified method considered the wind load as an external load. The following conclusion can be drawn from the results:  Figure 16 presents the axial stress and damage rate of the mooring system. The results show that the mooring system experiences a smaller fatigue load under normal operating condition than that under extreme sea conditions (LC 2 and LC 3). According to the time domain response of ML 2, although the mean value of ML 2 is larger than that of ML 1, the STD value of the mooring line tension is larger for ML 1, which has a significant impact on fatigue damage. A comparison between ML 1 and ML 2 indicates that ML 1 suffers more fatigue damage than ML 2, which should be paid more attention. Moreover, under extreme sea conditions (LC 2 and LC 3), the mooring lines experience more fatigue damage than under normal conditions. The time series statistics in Table 7 show that the mean value and STD value of ML 1 are larger using the coupled method, while those of ML 2 are larger using the simplified method. The comparison of fatigue damage between the coupled method and simplified method shows that the fatigue damage using the coupled method is larger than that using the simplified method for ML 1, while the fatigue of ML 2 is the opposite.  Figure 16 presents the axial stress and damage rate of the mooring system. The results show that the mooring system experiences a smaller fatigue load under normal operating condition than that under extreme sea conditions (LC 2 and LC 3). According to the time domain response of ML 2, although the mean value of ML 2 is larger than that of ML 1, the STD value of the mooring line tension is larger for ML 1, which has a significant impact on fatigue damage. A comparison between ML 1 and ML 2 indicates that ML 1 suffers more fatigue damage than ML 2, which should be paid more attention. Moreover, under extreme sea conditions (LC 2 and LC 3), the mooring lines experience more fatigue damage than under normal conditions. The time series statistics in Table 7 show that the mean value and STD value of ML 1 are larger using the coupled method, while those of ML 2 are larger using the simplified method. The comparison of fatigue damage between the coupled method and simplified method shows that the fatigue damage using the coupled method is larger than that using the simplified method for ML 1, while the fatigue of ML 2 is the opposite.

Conclusions
In this paper, the dynamic response analysis of the WindFloat SFOWT at a moderate water depth was carried out. A comparison of the dynamic responses between the coupled method and simplified method is made. The coupled method considers the motion interaction between the wind turbine and tower top and the aerodynamic damping, while the simplified method considered the wind load as an external load. The following conclusion can be drawn from the results:

Conclusions
In this paper, the dynamic response analysis of the WindFloat SFOWT at a moderate water depth was carried out. A comparison of the dynamic responses between the coupled method and simplified method is made. The coupled method considers the motion interaction between the wind turbine and tower top and the aerodynamic damping, while the simplified method considered the wind load as an external load. The following conclusion can be drawn from the results: (1) The WindFloat SFOWT motions and force responses are quite reasonable at a moderate water depth, showing the possibility of semisubmersible floating wind turbines installing at moderate depths. The amplitude of the surge motion response is 9.07 m under a normal operating condition. For the extreme conditions, the amplitude of the surge motion response is approximately 20 m. The maximum pitch motion for all load cases is below 5 • . For the mooring line tension responses, the maximum value for ML 1 and ML 2 are 1337.7 kN and 1791.6 kN, respectively, under normal operating conditions. The PSD analysis of the WindFloat SFOWT under different load cases shows that the low frequency responses are the main part of the surge motion response for all load cases, while the heave and pitch motion are mainly dominated by the wave excitation force. The fatigue damage results indicate that the points on the leeward side experience greater fatigue damage than those points on the heading wave side, to which should be paid more attention. For the mooring line tension fatigue, the fatigue damage on ML 1 is larger than that of ML 2 for all load cases. (2) A comparison between the coupled and simplified methods shows that the simplified method overestimates the maximum motion responses of the WindFloat SFOWT, while the STD value does not change a lot. The simplified method predicts a close PSD response in surge motion. The response spectra in the low frequency region of the heave and pitch motion are larger when using the simplified method. For the mooring line tension responses, the calculated response spectra of the simplified method are underestimated under different load cases for both the low frequency and wave frequency parts, while the response spectra of ML 2 are the opposite. Under normal operating conditions, the difference between the coupled method and simplified method is larger, while it is close under extreme sea conditions because the wave excitation force dominates the tower base responses.