Stochastic Dynamic Response Analysis of a 10 MW Tension Leg Platform Floating Horizontal Axis Wind Turbine

: The dynamic response of ﬂoating horizontal axis wind turbines (FHWATs) are affected by the viscous and inertia effects. In free decay motion, viscous drag reduces the amplitude of pitch and roll ﬂuctuation, the quasi-static mooring system overestimates the resonant amplitude values of pitch and roll. In this paper, the quasi-static mooring system is modiﬁed by introducing linear damping and quadratic damping. The dynamic response characteristics of the FHAWT modiﬁed model of the DTU 10 MW tension leg platform (TLP) were studied. Dynamic response of the blade was mainly caused by wind load, while the wave increased the blade short-term damage equivalent load. The tower base bending moment was affected by inclination of the tower and the misaligned angle β wave between wind and wave. Except the yaw motion, other degrees of freedom motions of the TLP were substantially affected by β wave . Ultimate tension of the mooring system was related to the displacement caused by pitch and roll motions, and standard deviation of the tension was signiﬁcantly affected by the wave frequency response. Under the action of wave load, the viscous drag would stimulate the mooring system and increase the resonance of the platform motion.


Introduction
In order to cope with the increasingly serious problems of fossil energy shortages and environmental pollution, renewable energy technologies have developed rapidly in recent years. Among the various types of renewable energy, wind energy has the most promising prospects and has received wide public attention [1,2]. So far, most wind farms have been deployed onshore due to the relatively lower cost of construction, operation, and maintenance [3]. However, offshore wind farms have higher wind speeds and smaller turbulence, and they are also less sensitive to space utilization, noise constraints, visual pollution and regulation issues. What is more, they can provide more and better green energy than onshore wind farms [4,5]. Although the cost of construction, operation and maintenance of offshore wind farms is relatively higher compared to onshore wind farms, as more and more research is being invested in offshore projects, as the advantages of offshore wind farms are gradually emerging. The substantial increment of offshore wind power capacity in the future may meet the need of the demand for green energy [4].
In shallow water areas (<50 m), offshore wind turbines, such as monopiles, gravity-based structures, jackets and buckets, are usually supported by bottom-fixed foundations [6,7]; when the offshore wind turbines are constructed in deep water areas (>50 m), the floating support platform has more advantages from an economic point of view [8,9]. According to the static stability method, there are mainly three types of FHAWTs: spar type, semi-submersibles, and TLP type [9][10][11][12]. The spar-type platform is usually a stable ballast in a small horizontal plane area; the semi-submersible floater neceds a large water line restoring moment to achieve sufficient stability; the TLP type is stabilized by tendons or tightened mooring lines to resists the wave-induced motions, thus reducing generator power variations and reducing the impact to the grid [6]. In the case of wind-wave alignment, both the spar platform and the semi-submersible platform showed large fatigue damage. In contrast, the limited platform motion of the TLP FHAWT could reduce the structural loads on the tower and the blades, resulting in less fatigue damage [13].
FHAWT operates in a complex environment which withstands loads from a variety of sources, including wind, wave, and ocean currents. Considering the mobility of floating platforms, FHAWT is a much more complex system than an onshore wind turbine or a stationary offshore wind turbine. Bachynski et al. conducted a series of studies in order to understand the dynamic response of the FHAWT more accurately [11,13]. According to the design considerations of a TLP FHAWT in deep sea water, they studied the design parameters of the TLP wind turbine which used full coupling simulation to evaluate the motion of the platform and the structural loads of the wind turbine components and mooring line. Significant works on TLP wind turbines had been done by the National Renewable Energy Laboratory and Massachusetts Institute of Technology [14]. Li et al. studied effects of the vertical heave motion, the aerodynamic forces and the quadratic wave forces on a TLP-type wind turbine [15]. Glosten Associates Company has designed PelaStar TLP for water depths ranging from 50 m to 200 m. This platform has the advantages of small displacement, economical to transport and install, low energy cost and capable of supporting 5-10 MW wind turbines [16]. Sant et al. conducted an experiment on TLP wind turbines, compared the power coefficient measured by the experiment with the power coefficient predicted by three independent models, and studied influences of the platform surge motion on the time average power coefficient [17]. Nihei et al. conducted a number of model tests to determine stochastic dynamic performances of three floating wind turbines on the Spar platform, TLP and semi-submersible platforms, respectively, and concluded that the TLP model showed the best performance [18]. However, due to the scale effect, the viscosity phenomenon cannot be accurately simulated in the model test.
The dynamic response of the TLP FHAWT depends on several factors such as the aerodynamic loads on the wind rotor, the hydrodynamic loads on the floating platform, the recovery effects of the mooring system, and the structural characteristics of the wind turbine and floating platform [19]. Viscous drag and inertia effects are two important factors affecting the dynamic response of FHAWT, approximate mooring systems simulated by quasi-static mooring lines (i.e., ignoring inertial effects and viscous drag) overestimates the pitch resonant amplitude values to some extent [20]. Farrugia's research showed a linear relationship between the amplitude of the surge velocity and the amplitude of the aerodynamic load response [21], which indicated that the aerodynamic load was affected by the surge motion [22]. The experimental results of Stewart et al. showed that the quadratic damping caused by viscous drag played an important role in the low-frequency zones where the radiation damping disappears [23]. Most studies assumed that the mooring system is a quasi-static system, in order to better predict the dynamic response of FHAWT, this paper considered the inertial effect and viscous drag in the mooring system and corrected the quasi-static mooring system. The stochastic dynamic response analysis of the modified 10 MW TLP FHAWT is carried out to understand the operating law and load characteristics of the TLP FHAWT more accurately.
Most of the studies on FHAWT dynamic analysis are limited to the environmental conditions in which wind waves are assumed to have the same direction of propagation. However, according to observational data, although there is less than 5% of the incident that the angle between wind and wave is greater than 60 degrees, the angle can usually reach 30 degrees [24]. Wind and waves may be significantly misaligned, especially under stable atmospheric conditions [25]. Barj et al. [26] observed that the spar-type FHAWT wind-wave misalignment could increase the extreme side-tower loads in numerical simulations. In the monopile installed at Bockstigen by Trumars, it was observed that the lateral load was increased under a 90 • misalignment conditions, but it was difficult to interpret the results due to the lack of information on environmental conditions (ECs) [27]. Aerodynamic and viscous damping are keys to determine the pitch resonance amplitude [20]. When the wave comes from a different direction with the wind direction, the aerodynamic damping is less effective in relieving the wave load [13]. Especially in the extreme case where the wind and wave direction are different from 90 • , because there is no aerodynamic damping, the inertial force generated by the sway motion is likely to cause roll resonance response. It is estimated that the viscous drag and inertia effect of the mooring system would be more obvious at this time. In addition, due to the high stiffness of the mooring system, TLP is susceptible to high frequency excitation, which generates heave resonance and pitch motion, and causes fatigue damage to the mooring system. When the wind and waves misalignment, the hydrodynamic load exceeds the aerodynamic load, which dominates the system response. Therefore, the modified model including viscous drag and inertia effects is used to study the wind and wave misalignment.
The structure of this paper is as follows: Section 2 introduces the numerical method of fully coupled dynamic analysis and the environmental conditions of wind turbine operation used in TLP FHAWT modeling simulation. Aiming at the problem that the quasi-static mooring system overestimates the surge and pitch motion, a method for modifying the model is proposed. Section 3 introduces the structure of the FHAWT and the data parameters used in the modeling and simulation, also corrects the simulation model. Section 4 discusses the dynamic response characteristics of TLP FHAWT based on statistical methods and frequency domain analysis methods: (1) Compared the differences in dynamic characteristics of TLP FHAWT and equivalent onshore HAWT; (2) Analyzed the response characteristics of wind turbines with and without waves. Based on the above analysis, the impact of wind-waves misalignment on the turbine is analyzed. The effects of inertial and viscous effects on the TLP wind turbine response are also demonstrated by typical wind and wave misalignment. The conclusions are summarized in Section 5.

Aerodynamic Modeling
Aerodynamic load on the blade is calculated based on the blade element momentum (BEM) theory and the wake influence was calculated by the induction factor [28]. Since BEM theory is originally developed for wind-stabilized wind turbines, the transient wind field should be corrected. Assuming that the pressure/momentum loss in the plane of the wind rotor is caused by the air flowing through the rotor blade planar element, the induced velocity resulting from loss of axial and tangential momentum in the flow field can be computed, and the aerodynamic and induced velocity near the rotor can be iteratively determined. The correction factor F is: Tip-Loss Model: Hub-Loss Model: The dynamic stall model [29] is commonly used to solve the transient aerodynamic problems caused by turbulent wind, yawing, rotational speed regulation and pitch regulation. The dynamic stall phenomena of airfoils have been extensively experimentally analyzed and studied [30,31]. In the literature, the B-L mode is widely used in helicopter and wind turbine analysis [3]. When using the semi-empirical Gonzalez and Minnema/Pierce models to calculate aerodynamic loads, a reasonable hysteresis can be generated in the normal force, tangential force, and pitching moment coefficient if the model parameters are set appropriately [32]. The potential-flow model uses the analytical potential-flow solution for the surrounding flow of a cylinder to model the tower dam effect on upwind rotors. Aerodynamic loads on the tower are calculated based on tower diameter, drag coefficient, and local relative velocity between free incoming flow and analytical node structure of each tower [32].

Hydrodynamic Modeling
In the simulation calculation, it is assumed that the TLP is a rigid structure and the mooring line is represented by a finite beam element. The loads acting on the TLP include first-order wave forces, quadratic sum-frequency forces and viscous forces. The first-order wave force is processed by the linear potential flow model. The quadratic sum-frequency forces is solved according to the complete quadratic potential flow. The viscous force is provided by the semi-empirical Morison equation [33]. The potential flow model in Sesam/Wadam deduces a first-order wave transfer function, additional mass, and radiation damping [34]. The general form of the complete time-domain equation of motion for coupled wind turbines and supporting platform systems is [35]: where M ij is the (i, j) component of the volumetric mass (inertia) matrix, which depends nonlinearly on the set of system degree of freedom q, control inputs u and time t; .. q j is the second time derivative of degree of freedom q j ; f i is the force function component associated with the degree of freedom q i ; The force function fi depends nonlinearly on the system degree of freedom and the first time derivatives (q and . q respectively), as well as the control inputs u and time t. The total external load on TLP are: where A ij is the (i, j) component of the hydrodynamic-added-mass matrix; F Hydro i is the ith component of the applied hydrodynamic load on the TLP associated with everything but A ij ; F Lines i is the ith component of the applied load on the TLP from the contribution of all mooring lines; For each TLP degree of freedom q i , where q 1 is surge, q 2 is sway, q 3 is heave, q 4 is roll, q 5 is pitch and q 6 is yaw.
In hydrodynamics problems, the form of F Hydro i is as follows: The expression of the first term F Waves i on the right-hand side of Equation (6) is: where W(ω) represents a Fourier transform of a realization of the white Gaussian noise time-series process with zero mean and unit variance; S 2−Sided ζ represents the expected two-sided power spectral density of the wave elevation per unit time; X i (ω, β) represents the wave-excitation force normalized per unit wave amplitude on the TLP.
The combination of the second and third terms on the right-hand side of Equation (6) ρgV 0 δ i3 − C Hydrostatic ij q j represents the load contribution from hydrostatics. ρ is the density of seawater; g is the gravitational acceleration; V 0 is the displaced volume of fluid when the TLP is in the undisplaced position; δ i3 is the (i,3) component (water depth) of the Kronecker-Delta function (unit matrix), and is the (i, j) component of the linear hydrostatic-restoring matrix. The last term on the right-hand side of Equation (6), − t 0 K ij (t − τ)q j (τ)dτ is a convolution integral representing the load contribution from wave-radiation damping.
For mooring systems approximated by quasi-static mooring lines, the total loads on the TLP from the contribution of all mooring lines F Lines i would be: where C Lines ij is the (i, j) component of the linearized restoring matrix from all mooring lines; F Lines,0 i is the ith component of the mooring system load acting on the TLP in the undisplaced position.
It can be seen from the above equations that the mooring system fitted with a quasi-static mooring line ignores the inertia effect and the viscous drag, which overestimates the pitch and the surge resonance amplitude to some extent. Because the coupling between axial and transversal motions contributes a lot to the pitch damping [20]. In order to correct the inertia effect neglected in the mooring system, two-thirds of the total dry weight of the mooring mass is added to the quality of the supporting platform. The quadratic damping is calculated using the method proposed by Hoff [36] to correct the viscous drag that is neglected in the mooring system. The pitch motion can be considered as a single degree of freedom system, and is assumed governed by the following form of an equation: where φ is the pitch angle; F( . φ) is the nonlinear pitch damping characteristics; G(φ) is the restoring moment characteristic; β is the assumed small scaling parameter.
For linear damping and quadratic damping terms: For small to moderate pitch angles, a positive cubic term restoring moment characteristic is usually appropriate: where ω 0 is the undamped natural pitch frequency. The average energy loss rate per cycle, namely the energy loss function L(V) is: where T(V) is the periodic time of an undamped oscillation.
In the calculation of L(V), V is held constant, and the variation of φ(t) thus corresponds to the case of undamped free oscillations. In this case, we can write: where U(φ) = φ 0 G(ξ)dξ; θ 0 (t) is a time-dependent phase angle. The final energy loss function L(V) is: where Perform dimensionless processing on the loss function L(V) to get the function Q(V) In view of the linear and quadratic damping term format, we may convert the loss function to a non-dimensional Q(V) function as: where: To find the appropriate value for b * 1 and b * 2 , assuming that there are N estimates Q(V 1 ), Q(V 2 ), · · · , Q(V N ), according to the least squares method: In order to get the best b * 1 and b * 2 , the e is derived and to obtain the simultaneous equation of b * 1 and b * 2 : The final solutions b * 1 and b * 2 : where:

Environmental Conditions
Normal wind profile model are used in the operational conditions. In the normal wind profile conditions, the wind profile U(z) is the average wind speed as a function of the height z above the free water level: (32) where U ref is the reference wind speed at hub height; z ref is hub height; α is the power law index, typically between 0.07 and 0.15, and the α value of the TLP FHAWT is set to 0.14 according to IEC 61400-3 [37].
Turbulent winds and irregular waves are considered related. For normal turbulence model, based on the Kaimal turbulence model for IEC Class C, turbulent wind fields are generated using the TurbSim program of NREL. For the irregular wave, the significant wave height (Hs) and Peak-spectral period (Tp) were set based on their correlation with the wind speed at the Statfjord site in the North Sea [38], and the irregular wave time history was generated by using the JONSWAP wave model.

Modeling
The 10 MW TLP FHAWT model included a 10 MW DTU wind turbine, a TLP, and a mooring system. As shown in Table 1, the wind turbine adopts the DTU 10 MW reference wind turbine developed by DTU Wind Energy and Vestas in the Light Rotor project [39]; Table 2 lists the TLP and mooring system characteristics. In order to study the impact of wind-wave misalignment, the rotor was aligned with the wind and the direction of the wave was defined relative to the wind direction under the conditions of EC1-EC7. During the process of research, one wind direction (0 • ) and four-wave directions (0 • , 30 • , 60 • , and 90 • ) were considered. The wind-wave misalignment angle was denoted by β wave , as shown in Figure 1. In order to study the impact of wind-wave misalignment, the rotor was aligned with the wind and the direction of the wave was defined relative to the wind direction under the conditions of EC1-EC7. During the process of research, one wind direction (0°) and four-wave directions (0°, 30°, 60°, and 90°) were considered. The wind-wave misalignment angle was denoted by βwave, as shown in

Validation
The fully coupled aero-hydro-servo-elastic time domain simulations performed on the DTU 10 MW reference wind turbine can be carried out by using the code FAST [40] which is a wind turbine design code developed by the National Renewable Energy Laboratory (NREL). The FAST simulation model was compared with the parameters of the DTU reference wind turbine to assess the accuracy. Under the same operating conditions, the rotor speed w and the blade pitch angle θ calculated by the calculation model must be consistent with the corresponding parameters in the DTU model for a given wind speed V m/s. Thrust is the most important aerodynamic factor in the overall dynamics of a wind turbine, so the calculation model needs to match the thrust of the DTU 10 MW reference wind turbine at each wind speed section. By matching the above key system features, the test and calibration of the calculation model were completed. As can be seen from Figure 2, except for the low wind speed region of 4-7 m/s, the remaining positions are well matched. The larger error at low wind speeds is due to the fact that the simulation model uses simple PI control and does not pitch to 5° as in the DTU reference wind turbine. Therefore, a wind speed of 7 m/s and above were selected for simulation.
Surge and pitch motions are the most important in the six motions of rigid-body floating platforms because they dominate the dynamic response of the system. The experimental data of the DTU 10MW reference wind turbine scaling model was used to test and verify the simulation results of the surge free decay motion [41]. All the experimental results were extrapolated to the full size.

Validation
The fully coupled aero-hydro-servo-elastic time domain simulations performed on the DTU 10 MW reference wind turbine can be carried out by using the code FAST [40] which is a wind turbine design code developed by the National Renewable Energy Laboratory (NREL). The FAST simulation model was compared with the parameters of the DTU reference wind turbine to assess the accuracy. Under the same operating conditions, the rotor speed w and the blade pitch angle θ calculated by the calculation model must be consistent with the corresponding parameters in the DTU model for a given wind speed V m/s. Thrust is the most important aerodynamic factor in the overall dynamics of a wind turbine, so the calculation model needs to match the thrust of the DTU 10 MW reference wind turbine at each wind speed section. By matching the above key system features, the test and calibration of the calculation model were completed. As can be seen from Figure 2, except for the low wind speed region of 4-7 m/s, the remaining positions are well matched. The larger error at low wind speeds is due to the fact that the simulation model uses simple PI control and does not pitch to 5 • as in the DTU reference wind turbine. Therefore, a wind speed of 7 m/s and above were selected for simulation.
Surge and pitch motions are the most important in the six motions of rigid-body floating platforms because they dominate the dynamic response of the system. The experimental data of the DTU 10 MW reference wind turbine scaling model was used to test and verify the simulation results of the surge free decay motion [41]. All the experimental results were extrapolated to the full size.
The relationships between the length, the time, and the mass of the full-size model and the experimental model were shown in Equation (33). In the surge free decay motion test, the full-scale displacement of the platform in the surge direction was 6.72 m. The comparison among experimental data, FAST simulation results, and modified model simulation results including inertial effects and viscous drag were shown in Figure 3. It can be seen from the figure that the free decay motion obtained by the modified model is closer to the experimental data than the FAST simulation. And the quasi-static Energies 2018, 11, 3341 9 of 24 mooring system of FAST simulation overestimates the amplitude of the turbulence to some extent due to the neglected of viscous drag: where λ = 60 is the geometric scaling factor; ρ sw and ρ f w are the densities of seawater and fresh water, respectively; L, T and M are the length, time and mass for the full-size model, respectively; l, t and m are the length, time and mass of the scaled experimental model, respectively. The relationships between the length, the time, and the mass of the full-size model and the experimental model were shown in Equation (33). In the surge free decay motion test, the full-scale displacement of the platform in the surge direction was 6.72 m. The comparison among experimental data, FAST simulation results, and modified model simulation results including inertial effects and viscous drag were shown in Figure 3. It can be seen from the figure that the free decay motion obtained by the modified model is closer to the experimental data than the FAST simulation. And the quasi-static mooring system of FAST simulation overestimates the amplitude of the turbulence to some extent due to the neglected of viscous drag: where λ = 60 is the geometric scaling factor; sw ρ and fw ρ are the densities of seawater and fresh water, respectively; L, T and M are the length, time and mass for the full-size model, respectively; l, t and m are the length, time and mass of the scaled experimental model, respectively. As shown in Figure 4, the pitch resonance frequency response predicted by the modified model is smaller than it predicted by FAST. The peak value of the modified model pitch response amplitude operator (RAO) is approximately 2/3 of the peak value predicted by FAST, which is caused by the damping effects of the viscous drag on the mooring system; The natural frequency of the modified model changes from 0.1918 Hz to 0.1896 Hz because of the neglecting of viscous damping and inertia effects of the mooring system. The difference could be explained by both the viscous damping effect and the inertial effect reducing the natural frequency. As shown in Figure 4, the pitch resonance frequency response predicted by the modified model is smaller than it predicted by FAST. The peak value of the modified model pitch response amplitude operator (RAO) is approximately 2/3 of the peak value predicted by FAST, which is caused by the damping effects of the viscous drag on the mooring system; The natural frequency of the modified model changes from 0.1918 Hz to 0.1896 Hz because of the neglecting of viscous damping and inertia effects of the mooring system. The difference could be explained by both the viscous damping effect and the inertial effect reducing the natural frequency.  Through the surge free decay motion verification and pitch response amplitude operator test, the results showed that the modified model could provide more accurate results for the stochastic dynamic response of the 10 MW FHAWT.

Simulation Research
For the modified model that considered the inertial effects and viscous drag in the mooring system, a series of representative working conditions EC1-EC7 as shown in Table 3 were selected to perform a 600 s fully coupled time domain simulation. The stochastic dynamic response of FHAWT under various wind and wave conditions was studied by the statistical analysis (mean, standard deviation, maximum/minimum value) and frequency domain analysis. Firstly, the differences between TLP FHAWT and equivalent onshore HAWT responses were compared and analyzed; then the effect of waves on TLP FHAWT was studied.  As shown in Figure 4, the pitch resonance frequency response predicted by the modified model is smaller than it predicted by FAST. The peak value of the modified model pitch response amplitude operator (RAO) is approximately 2/3 of the peak value predicted by FAST, which is caused by the damping effects of the viscous drag on the mooring system; The natural frequency of the modified model changes from 0.1918 Hz to 0.1896 Hz because of the neglecting of viscous damping and inertia effects of the mooring system. The difference could be explained by both the viscous damping effect and the inertial effect reducing the natural frequency.  Through the surge free decay motion verification and pitch response amplitude operator test, the results showed that the modified model could provide more accurate results for the stochastic dynamic response of the 10 MW FHAWT.

Simulation Research
For the modified model that considered the inertial effects and viscous drag in the mooring system, a series of representative working conditions EC1-EC7 as shown in Table 3 were selected to perform a 600 s fully coupled time domain simulation. The stochastic dynamic response of FHAWT under various wind and wave conditions was studied by the statistical analysis (mean, standard deviation, maximum/minimum value) and frequency domain analysis. Firstly, the differences between TLP FHAWT and equivalent onshore HAWT responses were compared and analyzed; then the effect of waves on TLP FHAWT was studied. Through the surge free decay motion verification and pitch response amplitude operator test, the results showed that the modified model could provide more accurate results for the stochastic dynamic response of the 10 MW FHAWT.

Simulation Research
For the modified model that considered the inertial effects and viscous drag in the mooring system, a series of representative working conditions EC1-EC7 as shown in Table 3 were selected to perform a 600 s fully coupled time domain simulation. The stochastic dynamic response of FHAWT under various wind and wave conditions was studied by the statistical analysis (mean, standard deviation, maximum/minimum value) and frequency domain analysis. Firstly, the differences between TLP FHAWT and equivalent onshore HAWT responses were compared and analyzed; then the effect of waves on TLP FHAWT was studied.

Comparison of TLP FHAWT and Onshore HAWT
In this section, the equivalent onshore HAWT used the same wind turbine blade as the TLP FHAWT with the difference that the former was mounted on a fixed base on the ground. Therefore, comparing and analyzing the stochastic dynamic responses of the two could provide a solution of designing the TLP FHAWT. In the comparative analysis process, the same turbulent wind loads were applied to the TLP FHAWT and onshore HAWT while the wave load was only applied to the TLP FHAWT.
The blade root is the connecting part between the blade and the turbine hub, and the load on it can reflect the load characteristics of the blade under the coupled action of wind and wave. Figure 5 shows the statistical analysis of the blade root bending moment and the corresponding short-term damage equivalent load of the TLP FHAWT and those of the equivalent onshore HAWT. It can be seen from Figure 5a,b that the two mean values of the blade root bending moments are almost the same, which indicates that the dynamic response of the blade is mainly caused by the wind load; The difference between the two standard deviations of blade root bending moment is quite large, which is related to the effects of wave load. The maximum difference is up to 3400 kNm when the wind speed is 14.7 m/s. It can be seen from Figure 5c that short-term damage equivalent load of the out-of-plane moment My at the blade root of the TLP FHWAT is larger than that of the onshore HAWT in the mid-high wind speed. The difference reaches a maximum value of 10,480 kNm when the wind speed is 14.7 m/s, which indicates that the wave load will increase the blade fatigue damage; It can be seen from Figure 5d that short-term damage equivalent load of the in-plane moment Mx at the blade root of the TLP FHAWT is basically consistent with that of the onshore HAWT, which indicates that the pitch control can reduce the influence of the wave load on the torque fluctuation.

Comparison of TLP FHAWT and Onshore HAWT
In this section, the equivalent onshore HAWT used the same wind turbine blade as the TLP FHAWT with the difference that the former was mounted on a fixed base on the ground. Therefore, comparing and analyzing the stochastic dynamic responses of the two could provide a solution of designing the TLP FHAWT. In the comparative analysis process, the same turbulent wind loads were applied to the TLP FHAWT and onshore HAWT while the wave load was only applied to the TLP FHAWT.
The blade root is the connecting part between the blade and the turbine hub, and the load on it can reflect the load characteristics of the blade under the coupled action of wind and wave. Figure 5 shows the statistical analysis of the blade root bending moment and the corresponding short-term damage equivalent load of the TLP FHAWT and those of the equivalent onshore HAWT. It can be seen from Figure 5a,b that the two mean values of the blade root bending moments are almost the same, which indicates that the dynamic response of the blade is mainly caused by the wind load; The difference between the two standard deviations of blade root bending moment is quite large, which is related to the effects of wave load. The maximum difference is up to 3400 kNm when the wind speed is 14.7 m/s. It can be seen from Figure 5c that short-term damage equivalent load of the out-ofplane moment My at the blade root of the TLP FHWAT is larger than that of the onshore HAWT in the mid-high wind speed. The difference reaches a maximum value of 10,480 kNm when the wind speed is 14.7 m/s, which indicates that the wave load will increase the blade fatigue damage; It can be seen from Figure 5d that short-term damage equivalent load of the in-plane moment Mx at the blade root of the TLP FHAWT is basically consistent with that of the onshore HAWT, which indicates that the pitch control can reduce the influence of the wave load on the torque fluctuation. The stability of the nacelle is very important for the safe operation of the wind turbine. Therefore, the fore-aft displacement and corresponding accelerations at the tower-top of the TLP FHAWT and the onshore HAWT were compared and analyzed. As shown in Figure 6a, the average fore-aft displacement at the tower-top of the TLP FHAWT is slightly larger than that of the equivalent onshore HAWT, and the difference is 0.02 m at the rated wind speeds. This was caused by the surge motion of the TLP. The standard deviation of the fore-aft displacement at the tower-top of the TLP FHAWT is greater than that of the equivalent onshore HAWT, and the difference between the two increases significantly as the wave increases. It can be explained by the wave load effects from Figure  6b, which indicates that the wave load will affect the stability of the nacelle. There was a similar conclusion about the maximum value of the fore-aft displacement at the tower-top. The spectrum analysis (under the rated conditions) of Figure 6d shows that the linear acceleration at the tower-top of the onshore HAWT is more sensitive to the 3P frequency, which indicates that the tower of the TLP FHAWT is less sensitive to the passing frequency of the blades. The stability of the nacelle is very important for the safe operation of the wind turbine. Therefore, the fore-aft displacement and corresponding accelerations at the tower-top of the TLP FHAWT and the onshore HAWT were compared and analyzed. As shown in Figure 6a, the average fore-aft displacement at the tower-top of the TLP FHAWT is slightly larger than that of the equivalent onshore HAWT, and the difference is 0.02 m at the rated wind speeds. This was caused by the surge motion of the TLP. The standard deviation of the fore-aft displacement at the tower-top of the TLP FHAWT is greater than that of the equivalent onshore HAWT, and the difference between the two increases significantly as the wave increases. It can be explained by the wave load effects from Figure 6b, which indicates that the wave load will affect the stability of the nacelle. There was a similar conclusion about the maximum value of the fore-aft displacement at the tower-top. The spectrum analysis (under the rated conditions) of Figure 6d shows that the linear acceleration at the tower-top of the onshore HAWT is more sensitive to the 3P frequency, which indicates that the tower of the TLP FHAWT is less sensitive to the passing frequency of the blades. The stability of the nacelle is very important for the safe operation of the wind turbine. Therefore, the fore-aft displacement and corresponding accelerations at the tower-top of the TLP FHAWT and the onshore HAWT were compared and analyzed. As shown in Figure 6a, the average fore-aft displacement at the tower-top of the TLP FHAWT is slightly larger than that of the equivalent onshore HAWT, and the difference is 0.02 m at the rated wind speeds. This was caused by the surge motion of the TLP. The standard deviation of the fore-aft displacement at the tower-top of the TLP FHAWT is greater than that of the equivalent onshore HAWT, and the difference between the two increases significantly as the wave increases. It can be explained by the wave load effects from Figure  6b, which indicates that the wave load will affect the stability of the nacelle. There was a similar conclusion about the maximum value of the fore-aft displacement at the tower-top. The spectrum analysis (under the rated conditions) of Figure 6d shows that the linear acceleration at the tower-top of the onshore HAWT is more sensitive to the 3P frequency, which indicates that the tower of the TLP FHAWT is less sensitive to the passing frequency of the blades. Moreover, the bending moment at tower base is an important factor to consider when designing the tower, so the bending moments at tower base of TLP FHAWT and equivalent onshore HAWT were investigated. It is observed in Figure 7a that the mean tower base fore-aft moment My of the TLP FHAWT is slightly greater than that of the equivalent onshore HAWT at low wind speeds, the difference reaches a maximum value of 4702.37 kNm at the rated wind speed. This could be explained by the non-zero tilt angle of the TLP FHAWT tower causes gravity to generate bending moments. The standard deviation of fore-aft bending moment at tower base of the TLP FHAWT is greater than that of the equivalent onshore HAWT, and the difference between the two increases significantly as the wave increases. The difference reaches a maximum value of 46,525.48 kNm when the wind speed is 24.8 m/s. There was also a similar conclusion about the maximum value of fore-aft bending moment at tower base. The standard deviation could reflect the fatigue load of the structure to a certain extent and the maximum value indicated the limit load that the structure faced. Therefore, compared with the onshore HAWT, the tower of the TLP FHAWT needed an enhanced design. It is observed in Figure 7b that the tower base fore-aft moment My of Moreover, the bending moment at tower base is an important factor to consider when designing the tower, so the bending moments at tower base of TLP FHAWT and equivalent onshore HAWT were investigated. It is observed in Figure 7a that the mean tower base fore-aft moment My of the TLP FHAWT is slightly greater than that of the equivalent onshore HAWT at low wind speeds, the difference reaches a maximum value of 4702.37 kNm at the rated wind speed. Moreover, the bending moment at tower base is an important factor to consider when designing the tower, so the bending moments at tower base of TLP FHAWT and equivalent onshore HAWT were investigated. It is observed in Figure 7a that the mean tower base fore-aft moment My of the TLP FHAWT is slightly greater than that of the equivalent onshore HAWT at low wind speeds, the difference reaches a maximum value of 4702.37 kNm at the rated wind speed. This could be explained by the non-zero tilt angle of the TLP FHAWT tower causes gravity to generate bending moments. The standard deviation of fore-aft bending moment at tower base of the TLP FHAWT is greater than that of the equivalent onshore HAWT, and the difference between the two increases significantly as the wave increases. The difference reaches a maximum value of 46,525.48 kNm when the wind speed is 24.8 m/s. There was also a similar conclusion about the maximum value of fore-aft bending moment at tower base. The standard deviation could reflect the fatigue load of the structure to a certain extent and the maximum value indicated the limit load that the structure faced. Therefore, compared with the onshore HAWT, the tower of the TLP FHAWT needed an enhanced design. It is observed in Figure 7b that the tower base fore-aft moment My of This could be explained by the non-zero tilt angle of the TLP FHAWT tower causes gravity to generate bending moments. The standard deviation of fore-aft bending moment at tower base of the TLP FHAWT is greater than that of the equivalent onshore HAWT, and the difference between the two increases significantly as the wave increases. The difference reaches a maximum value of 46,525.48 kNm when the wind speed is 24.8 m/s. There was also a similar conclusion about the maximum value of fore-aft bending moment at tower base. The standard deviation could reflect the fatigue load of the structure to a certain extent and the maximum value indicated the limit load that the structure faced. Therefore, compared with the onshore HAWT, the tower of the TLP FHAWT needed an enhanced design. It is observed in Figure 7b that the tower base fore-aft moment My of the TLP FHAWT has a amplitude near the wave frequency of 0.12 Hz, again indicating that the wave load has a significant effect on the tower. There was a similar result observed for the tower base side-to-side moment Mx of TLP FHAWT.

Influence of Wave on TLP FHAWT
Many studies on the TLP FHAWT dynamic analysis were limited to the environmental conditions of wind load and wave load in the same propagation direction. However, the wind and wave can be significantly misaligned, and it is expected that such misalignments would affect the motions and structural responses of the wind turbine.

The Influence of the Wave
The mean value of out-of-plane moment My at the blade root under the waving conditions was equal to that under the non-waving conditions. It means that the wave loads had very limited effects on out-of-plane moment My at the blade root. There was a similar result observed from the in-plane moment Mx at the blade root. Those were consistent with the conclusion from the comparative analysis of TLP FHAWT and equivalent onshore HAWT in Section 4.1 that the loads on the wind turbine blade were mainly caused by turbulent wind. The blade structure design mainly considered the wind load, but it should be noted that the wave would increase the blade fatigue load.
Tower base bending moments were significantly affected by wave loads and it was up to the nacelle movement relative to the TLP, the thrust, the inertia transmitted to the platform and the mooring line tension. According to the Figure 8a, the mean value of tower base fore-aft moment My under the waves conditions is consistent with that under the non-waves conditions, indicating that the mean value is mainly affected by the wind; The large difference in standard deviations can be explained based on the power spectrum (under the rated conditions) of Figure 8b that the wave load causes the tower to generate wave frequency response. When the wind speed was 24.8 m/s, the standard deviation of tower base fore-aft moment My reached the maximum value of 46,130.38 kNm. There was a similar conclusion about tower base side-to-side moment Mx. Therefore, the tower fatigue load caused by waves must be taken into consideration during the design phase. the TLP FHAWT has a amplitude near the wave frequency of 0.12 Hz, again indicating that the wave load has a significant effect on the tower. There was a similar result observed for the tower base sideto-side moment Mx of TLP FHAWT.

Influence of Wave on TLP FHAWT
Many studies on the TLP FHAWT dynamic analysis were limited to the environmental conditions of wind load and wave load in the same propagation direction. However, the wind and wave can be significantly misaligned, and it is expected that such misalignments would affect the motions and structural responses of the wind turbine.

The Influence of the Wave
The mean value of out-of-plane moment My at the blade root under the waving conditions was equal to that under the non-waving conditions. It means that the wave loads had very limited effects on out-of-plane moment My at the blade root. There was a similar result observed from the in-plane moment Mx at the blade root. Those were consistent with the conclusion from the comparative analysis of TLP FHAWT and equivalent onshore HAWT in Section 4.1 that the loads on the wind turbine blade were mainly caused by turbulent wind. The blade structure design mainly considered the wind load, but it should be noted that the wave would increase the blade fatigue load.
Tower base bending moments were significantly affected by wave loads and it was up to the nacelle movement relative to the TLP, the thrust, the inertia transmitted to the platform and the mooring line tension. According to the Figure 8a, the mean value of tower base fore-aft moment My under the waves conditions is consistent with that under the non-waves conditions, indicating that the mean value is mainly affected by the wind; The large difference in standard deviations can be explained based on the power spectrum (under the rated conditions) of Figure 8b that the wave load causes the tower to generate wave frequency response. When the wind speed was 24.8 m/s, the standard deviation of tower base fore-aft moment My reached the maximum value of 46,130.38 kNm. There was a similar conclusion about tower base side-to-side moment Mx. Therefore, the tower fatigue load caused by waves must be taken into consideration during the design phase.  Figure 9 shows the statistical analysis and spectral analysis (under the rated conditions) of the surge, heave and pitch motions of the TLP FHAWT with and without waves. Motions of sway, roll and yaw were not discussed because they were very small due to the symmetry of structures and external loads [42]. It can be seen from Figure 9a that the two mean values (one with wave and another without wave) of surge motion of the TLP are substantially identical. And similar results was  Figure 9 shows the statistical analysis and spectral analysis (under the rated conditions) of the surge, heave and pitch motions of the TLP FHAWT with and without waves. Motions of sway, roll and yaw were not discussed because they were very small due to the symmetry of structures and external loads [42]. It can be seen from Figure 9a that the two mean values (one with wave and another without  Figure 9c,e for the heave and pitch motions. As can be seen from Figure 9b,d, the standard deviations of the surge motion and heave motion are mainly composed of the wave frequency response component and the surge resonance response component, and the former is much smaller than latter induced by low-frequency wind. Since the heave resonance frequency was high, there was no necessary to consider the heave resonance response. The pitch natural frequency is reduced to 1.85 Hz as shown in Figure 9f due to the shift of the pitch natural frequency caused by the TLP-tower coupling effect [43]. The standard deviation of pitch motion was the combined actions of wind and wave loads, i.e., the inertia of the surge motion and the hydrodynamic coupling. The pitch motion consisted of surge resonance response component, wave response component, and pitch resonance response component. Based on the above analysis, it could be inferred that the surge resonance and the pitch resonance affected the motions of the TLP platform deeply, and they should be considered in the TLP design phase as key issues. observed from the Figure 9c,e for the heave and pitch motions. As can be seen from Figure 9b and Figure 9d, the standard deviations of the surge motion and heave motion are mainly composed of the wave frequency response component and the surge resonance response component, and the former is much smaller than latter induced by low-frequency wind. Since the heave resonance frequency was high, there was no necessary to consider the heave resonance response. The pitch natural frequency is reduced to 1.85 Hz as shown in Figure 9f due to the shift of the pitch natural frequency caused by the TLP-tower coupling effect [43]. The standard deviation of pitch motion was the combined actions of wind and wave loads, i.e., the inertia of the surge motion and the hydrodynamic coupling. The pitch motion consisted of surge resonance response component, wave response component, and pitch resonance response component. Based on the above analysis, it could be inferred that the surge resonance and the pitch resonance affected the motions of the TLP platform deeply, and they should be considered in the TLP design phase as key issues. Since the mooring system was related to the safety of the wind turbine, the mooring line tension was also investigated. While the mooring line 3 was in the upwind position, with greater tension than the downwind mooring lines, the mooring line 3 was selected to be analyzed. It can be seen from the Figure 10a that the wave load has almost no effect on the mean value of mooring line tension, but it has a great influence on the standard deviation. According to the Figure 10b, the standard deviation response of the tension is affected by the wave-frequency excitation and the pitch resonance response, which is similar to the pitch motion in Figure 9. Therefore, reducing the surge resonance response and pitch resonance response were effective ways to ensure the safety of the TLP FHAWT. Since the mooring system was related to the safety of the wind turbine, the mooring line tension was also investigated. While the mooring line 3 was in the upwind position, with greater tension than the downwind mooring lines, the mooring line 3 was selected to be analyzed. It can be seen from the Figure 10a that the wave load has almost no effect on the mean value of mooring line tension, but it has a great influence on the standard deviation. According to the Figure 10b, the standard deviation response of the tension is affected by the wave-frequency excitation and the pitch resonance response, which is similar to the pitch motion in Figure 9. Therefore, reducing the surge resonance response and pitch resonance response were effective ways to ensure the safety of the TLP FHAWT. (e) (f) Since the mooring system was related to the safety of the wind turbine, the mooring line tension was also investigated. While the mooring line 3 was in the upwind position, with greater tension than the downwind mooring lines, the mooring line 3 was selected to be analyzed. It can be seen from the Figure 10a that the wave load has almost no effect on the mean value of mooring line tension, but it has a great influence on the standard deviation. According to the Figure 10b, the standard deviation response of the tension is affected by the wave-frequency excitation and the pitch resonance response, which is similar to the pitch motion in Figure 9. Therefore, reducing the surge resonance response and pitch resonance response were effective ways to ensure the safety of the TLP FHAWT.

The Influence of the Wind-Wave Misalignment
From the results of the study with and without wave loading in Section 4, it was expected that wind-wave misalignment would affect the dynamic response of the TLP FHAWT, including the structural loads of the blade and tower, platform motion and mooring line tension. The mean of the dynamic response of TLP FHAWT was mainly induced by wind, so the dynamic standard deviation (STD) of TLP FHAWT was mainly analyzed [44].

The Influence of the Wind-Wave Misalignment
From the results of the study with and without wave loading in Section 4, it was expected that wind-wave misalignment would affect the dynamic response of the TLP FHAWT, including the structural loads of the blade and tower, platform motion and mooring line tension. The mean of the dynamic response of TLP FHAWT was mainly induced by wind, so the dynamic standard deviation (STD) of TLP FHAWT was mainly analyzed [44].
The standard deviations of the blade root bending moments, the tower base bending moments and the tension of the mooring lines were shown as functions of the wind-wave misalignment angle β wave in Figure 11. As shown in Figure 11a,b, the standard deviations of the blade in-plane moment Mx at the blade root and out-of-plane moment My at the blade root are basically not affected by β wave at the low wind speed. While under the high wind speed conditions, the standard deviation of the in-plane moment Mx at the blade root increases slightly as β wave increases and the standard deviation of the out-of-plane moment My at the blade root decreases as β wave increases. According to Figure 11c,d, it can be seen that the tower base moment changes significantly as β wave increases. This was because the tilting of the platform shifts the center of gravity of the wind turbine nacelle thus changed the tower base bending moment. When β wave changed, the weight of the rotor and nacelle components changed the contribution to the tower base side-to-side moment Mx and tower base fore-aft moment My. The standard deviation of the tower base side-to-side moment Mx reaches a maximum value of 58,219.14 kNm when the FHAWT operates in EC7, β wave is 90 • , and the standard deviation of tower base fore-aft moment My reaches a maximum value of 73,157.88 kNm when the FHAWT operates in EC7, β wave is 0 • . Under the low wind speed conditions EC1-EC2, the fore-aft tilt angles of the tower are very small, the wind and wave misalignment has little influence on the fore-aft bending moments. While at mid-high wind speeds, the tilt angle of the tower increases, resulting in a significant influence. The mooring line 2 and the mooring line 3 were respectively located on the rotating tension side and the downwind side of the rotor, which meant they would bear larger loads than other mooring lines. It can be seen from Figure 11e,f that the standard deviation of tension of the mooring line 2 increases with the increment of β wave , while the standard deviation of the tension of the mooring line 3 is opposite. This performance is more pronounced at high wind speed conditions. The standard deviation of the tension of mooring line 2 reaches a maximum value of 329.55 kN when the wind speed is 24.8 m/s, β wave is 90 • , and the standard deviation of the tension of mooring line 3 reaches a maximum value of 552.12kN when the wind speed is 24.8 m/s, β wave is 0 • . The mooring line 3 had the same tendency as the pitch motion, which meant that the pitch motion had a great influence on the tension of the mooring line 3.
of the out-of-plane moment My at the blade root decreases as βwave increases. According to Figure  11c,d, it can be seen that the tower base moment changes significantly as βwave increases. This was because the tilting of the platform shifts the center of gravity of the wind turbine nacelle thus changed the tower base bending moment. When βwave changed, the weight of the rotor and nacelle components changed the contribution to the tower base side-to-side moment Mx and tower base fore-aft moment My. The standard deviation of the tower base side-to-side moment Mx reaches a maximum value of 58,219.14 kNm when the FHAWT operates in EC7, βwave is 90°, and the standard deviation of tower base fore-aft moment My reaches a maximum value of 73,157.88 kNm when the FHAWT operates in EC7, βwave is 0°. Under the low wind speed conditions EC1-EC2, the fore-aft tilt angles of the tower are very small, the wind and wave misalignment has little influence on the fore-aft bending moments. While at mid-high wind speeds, the tilt angle of the tower increases, resulting in a significant influence. The mooring line 2 and the mooring line 3 were respectively located on the rotating tension side and the downwind side of the rotor, which meant they would bear larger loads than other mooring lines. It can be seen from Figure 11e,f that the standard deviation of tension of the mooring line 2 increases with the increment of βwave, while the standard deviation of the tension of the mooring line 3 is opposite. This performance is more pronounced at high wind speed conditions. The standard deviation of the tension of mooring line 2 reaches a maximum value of 329.55 kN when the wind speed is 24.8 m/s, βwave is 90°, and the standard deviation of the tension of mooring line 3 reaches a maximum value of 552.12kN when the wind speed is 24.8 m/s, βwave is 0°. The mooring line 3 had the same tendency as the pitch motion, which meant that the pitch motion had a great influence on the tension of the mooring line 3.  The standard deviation of the platform motion was shown as a function of the wind-wave misalignment angle βwave in Figure 12. Figure 12a shows that the standard deviation of the surge motion increases as βwave increases at low to medium wind conditions EC1-EC5. This was because the aerodynamic damping decreases as βwave increases, resulting in an increment in low frequency surge motion. It can also be observed that the standard deviation of surge motion under high wind speed condition EC6-7 decreases as βwave increases. This was because the viscous hydrodynamic load in the TLP surge direction was greatly reduced as the wind speed increases, so that the wave-frequency surge motion decreased significantly as βwave increased. As shown in Figure 12b, the standard deviation of the sway motion at a low wind speed is the minimum when βwave is 30°, and in the case of mid-high wind speed, the sway motion increases as βwave increases due to the larger effective wave The standard deviation of the platform motion was shown as a function of the wind-wave misalignment angle β wave in Figure 12. Figure 12a shows that the standard deviation of the surge motion increases as β wave increases at low to medium wind conditions EC1-EC5. This was because the aerodynamic damping decreases as β wave increases, resulting in an increment in low frequency surge motion. It can also be observed that the standard deviation of surge motion under high wind speed condition EC6-7 decreases as β wave increases. This was because the viscous hydrodynamic load in the TLP surge direction was greatly reduced as the wind speed increases, so that the wave-frequency surge motion decreased significantly as β wave increased. As shown in Figure 12b, the standard deviation of the sway motion at a low wind speed is the minimum when β wave is 30 • , and in the case of mid-high wind speed, the sway motion increases as β wave increases due to the larger effective wave height and the longer wave period. As shown in Figure 12c, the standard deviation of the heave motion at low and medium wind speeds increases slightly with the increment of β wave , because the viscous fluid dynamic load decreases slightly as the increment of β wave . While at high wind speeds, the dynamic load of viscous fluid decreases significantly as β wave increases, and the standard deviation of the heave motion decreases with the increase of β wave . According to Figure 12d, it can be seen that the standard deviation of the roll motion increases as β wave increases.
From Figure 12e, it can be seen that the standard deviation of the pitch motion decreases as β wave increases. This was because as β wave increases, the contribution of the wave load to the Y direction increased and the contribution to the X direction decreased.
As shown in Figure 12f, the standard deviation of yaw motion is not sensitive to the wave directions of all ECs, indicating that it is mainly affected by changes in wind speed. From the above analysis, we could conclude that, with the exception of yaw motion, the other degrees of freedom motions of TLP were significantly affected by the misalignment of wind and waves. The influence of wind and wave misalignment on the platform could be reduced by increasing the aerodynamic damping of the wind turbine.
From the above analysis, it could be concluded that the wind-wave misalignment had obvious influence on the surge, pitch, sway and roll motion of the TLP. In order to demonstrate the influence of the viscous drag and inertial effect of quasi-static mooring system on the stochastic dynamic response of wind turbines, the power spectrum comparison analysis was carried out for the standard deviation of TLP motions of the original model and the correction model under wind-wave misalignment. Figure 13 shows the spectral analysis of the standard deviation of surge and pitch motion at the wind speed of 24.8 m/s (the maximum standard deviation of the pitch motion is occurred at this time). Among them, Figure 13a,b show the situation when β wave is 0 • , and Figure 13c,d show the situation when β wave is 90 • . It can be seen from Figure 13a,c that, for the surge motion, the surge component in the correction model is larger, indicating that the viscous drag has an effect on the fluctuation of the TLP motion; according to the power spectrum analysis of pitch motion in Figure 13b,d, for the modified model, the inertial force generated by the surge motion causes larger pitch resonance response, which indicates that the viscous resistance and inertia effect of the mooring system should be considered. It can be seen by comparing Figure 13a,c that the fluctuations of the surge motion are different for different β wave angles. It has bigger effect on the surge motion when β wave is 0 • , while the effect is smaller when β wave is 90 • . There are similar conclusions about the pitch motion. (e) (f) From Figure 12e, it can be seen that the standard deviation of the pitch motion decreases as βwave increases. This was because as βwave increases, the contribution of the wave load to the Y direction increased and the contribution to the X direction decreased.
As shown in Figure 12f, the standard deviation of yaw motion is not sensitive to the wave directions of all ECs, indicating that it is mainly affected by changes in wind speed. From the above analysis, we could conclude that, with the exception of yaw motion, the other degrees of freedom From the above analysis, it could be concluded that the wind-wave misalignment had obvious influence on the surge, pitch, sway and roll motion of the TLP. In order to demonstrate the influence of the viscous drag and inertial effect of quasi-static mooring system on the stochastic dynamic response of wind turbines, the power spectrum comparison analysis was carried out for the standard deviation of TLP motions of the original model and the correction model under wind-wave misalignment.  Figure 13 shows the spectral analysis of the standard deviation of surge and pitch motion at the wind speed of 24.8 m/s (the maximum standard deviation of the pitch motion is occurred at this time).
Among them, Figure 13a,b show the situation when βwave is 0°, and Figure 13c,d show the situation when βwave is 90°. It can be seen from Figure 13a,c that, for the surge motion, the surge component in the correction model is larger, indicating that the viscous drag has an effect on the fluctuation of the TLP motion; according to the power spectrum analysis of pitch motion in Figure 13b,d, for the modified model, the inertial force generated by the surge motion causes larger pitch resonance response, which indicates that the viscous resistance and inertia effect of the mooring system should be considered. It can be seen by comparing Figure 13a,c that the fluctuations of the surge motion are  Figure 14 shows the sway and roll spectrum analysis of the standard deviation of the sway and roll motion at the wind speed of 24.8 m/s (the maximum standard deviation of the roll motion occurs at this time). Among them, Figure 14a,b show the situation when β wave is 0 • , and Figure 14c,d show the situation when β wave is 90 • . It can be seen from Figure 14a,c that, when wind and wave in the same direction, the fluctuation of the sway motion is not much different. While for β wave is 90 • , due to the inclusion of viscous drag, a larger deviation of the sway motion is produced under the wave loads in the corrected model. According to the power spectrum analysis of the roll motion in Figure 14b,d, when β wave is 90 • , for the modified model, the inertial force generated by the sway motion causes a larger roll resonance response, which indicates that the viscous drag and inertia effects of the mooring system should be considered. Comparing Figure 14a,c, the fluctuations of the sway motion are different for different β wave angles. Which is exactly the opposite of the surge motion, it has smaller effects on the sway motion when β wave is 0 • , while the effect is bigger when β wave is 90 • . There are similar conclusions about the roll motion.
Energies 2018, 11, x FOR PEER REVIEW 21 of 24 different for different βwave angles. It has bigger effect on the surge motion when βwave is 0°, while the effect is smaller when βwave is 90°. There are similar conclusions about the pitch motion.  Figure 14 shows the sway and roll spectrum analysis of the standard deviation of the sway and roll motion at the wind speed of 24.8 m/s (the maximum standard deviation of the roll motion occurs at this time). Among them, Figure 14a,b show the situation when βwave is 0°, and Figure 14c,d show the situation when βwave is 90°. It can be seen from Figure 14a,c that, when wind and wave in the same direction, the fluctuation of the sway motion is not much different. While for βwave is 90°, due to the inclusion of viscous drag, a larger deviation of the sway motion is produced under the wave loads in the corrected model. According to the power spectrum analysis of the roll motion in Figure 14b,d, when βwave is 90°, for the modified model, the inertial force generated by the sway motion causes a larger roll resonance response, which indicates that the viscous drag and inertia effects of the mooring system should be considered. Comparing Figure 14a,c, the fluctuations of the sway motion are different for different βwave angles. Which is exactly the opposite of the surge motion, it has smaller effects on the sway motion when βwave is 0°, while the effect is bigger when βwave is 90°. There are similar conclusions about the roll motion.

Conclusions
In this paper, the quasi-mooring system was modified by calculating the quadratic damping, and an improved model simulating the 10 MW TLP FHAWT was constructed based on the modified

Conclusions
In this paper, the quasi-mooring system was modified by calculating the quadratic damping, and an improved model simulating the 10 MW TLP FHAWT was constructed based on the modified mooring system. Then the response characteristics of the TLP FHAWT were analyzed and compared with the equivalent onshore HAWT to determine effects of the wave on dynamic responses of the TLP FHAWT. Moreover, the influence of turbulent wind loads and wave loads misalignment on FHAWT under typical operating conditions was studied. Finally, the differences between the correction model and the original model TLP response are analyzed. On the basis of those works, it can be concluded that: (1) Both mean values of blade root bending moments of the TLP FHAWT and of the equivalent onshore HAWT were almost the same. Loads on the blade root were mainly caused by turbulent wind, but short-term damage equivalent loads on the blade were dramatically affected by the wave. Due to inclination of the tower, mean value of tower bending moment of the TLP FHAWT was greater than it of the equivalent onshore HAWT. At the rated wind speed, My value of the TLP FHAWT exceeded it of the equivalent Onshore HAWT My value by 4702.37 kNm. (2) Vibration of the TLP FHAWT tower was substantially affected by the wind-wave misalignment.
The aerodynamic damping generated by rotation of the rotor could significantly reduce vibration of the tower caused by wave loads. As β wave increased, inclination of the platform caused the center of gravity of the nacelle to shift which further increased bending moments of the tower base. Therefore, reinforcement design for the TLP FHAWT tower should supposed to be enhanced. (3) Ultimate tension of the mooring system was mainly related to the displacement caused by pitch motion. Standard deviation of the tension was significantly affected by the wave frequency response. Therefore, ultimate tension of the mooring line should be decreased by reducing the pitch motion. (4) Amplitudes of the pitch and roll motion were overestimated in the free decay motion due to the neglected of viscous drag of the mooring system. Viscous drag would increase the resonance of the platform under the wave loads, while the quasi-static mooring system correction model could capture those resonances better. Wave loads mainly affected the pitch and surge motion when β wave is 0 • . The inertial force caused by the surge motion would generate high-frequency moments which would cause pitch resonance. When the wind-wave misalignment, the inertial force caused by the sway motion would generate high-frequency moments without the inhibition of aerodynamic damping, which would cause roll resonance.