Hydrodynamic Response of a Combined Wind–Wave Marine Energy Structure

: Due to the energy crisis and greenhouse e ﬀ ect, o ﬀ shore renewable energy is attracting increasing attention worldwide. Various o ﬀ shore renewable energy systems, such as ﬂoating o ﬀ shore wind turbines (FOWTs), and wave energy converters (WECs), have been proposed and developed so far. To increase power output and reduce related costs, a combined marine energy structure using FOWT and WEC technologies has been designed, analyzed and presented in the present paper. The energy structure combines a 5-MW braceless semisubmersible FOWT and a heave-type WEC which is installed on the central column of the semisubmersible. Wave power is absorbed by a power take-o ﬀ (PTO) system through the relative heave motion between the central column of the FOWT and the WEC. A numerical model has been developed and is used to determine rational size and draft of the combined structure. The e ﬀ ects of di ﬀ erent PTO system parameters on the hydrodynamic performance and wave energy production of the WEC under typical wave conditions are investigated and a preliminary best value for the PTO’s damping coe ﬃ cient is obtained. Additionally, the e ﬀ ects of viscous modeling used during the analysis and the hydrodynamic coupling on the response of the combined structure are studied.


Introduction
Due to the exhaustion of fossil fuels and their pollution on the environment, people are trying to develop various clean and renewable energy technologies to replace conventional energy resources. Offshore wind energy has experienced rapid development [1] during the past ten years however it faces significant challenges to ensure commercial viability due to the policy of grid parity [2]. The ocean contains a variety of clean and renewable resources such as wind energy, wave energy and tidal energy, etc. Wave energy is still in the very early stages of development. The integration of offshore wind energy and wave energy power generation by sharing space, transmission infrastructure, operation and maintenance (O&M) activities and supporting structure will improve the utilization efficiency of marine renewable energy, but also reduce their levelized cost of energy.
Many scholars tried to develop different wind-wave combined systems. The EU funded Marine Renewable Integrated Application Platform (MARINA Platform) program [3] that proposed possible combined systems for use. For example, two combined systems have been proposed, which respectively J. Mar. Sci. Eng. 2020, 8,253 2 of 23 consist of the same WindFloat semisubmersible turbine platform with a point-absorber-type wave energy converter (WEC) and an oscillating-water-column type WEC. Peiffer et al. [4] and Aubault et al. [5] investigated the coupled dynamic analysis by numerical and laboratory models. Muliawan et al. [6,7] proposed an integrated concept combining a spar-type floating wind turbine with a torus WEC (named as STC). The dynamic response characteristics of the combined structure in typical operational and extreme sea states were studied. Further, Wan et al. [8][9][10] studied the dynamic response characteristics of the STC under typical operational and extreme sea states based on numerical and experimental methods. Another combined structure combining a floating semisubmersible wind turbine and flap-type WECs (named as SFC) was proposed. Michailides et al. [11][12][13] systematically studied the integrated operation of the SFC by numerical simulation and physical model tests and it was found that the total produced power can be increased. Moreover, Ren et al. [14] studied a new combined structure consisting of a monopile type wind turbine and a heave-type wave energy converter (denoted as MWWC). However, the conceptual designs of combined structure based a floating semi-submersible wind turbine and WECs in deep water have been insufficient so far.
In the present paper, we introduce a combined marine energy structure consisting of a 5-MW braceless wind turbine developed in [15] with a heave-type WEC inspired by the STC system [16,17], which was proposed first by our team. At the first stage, it aims to investigate the dynamic response and its geometry and various parameters. Based on potential flow theory by means of ANSYS/AQWA [18], we developed a numerical model, and the best diameter and draft of the WEC are selected. Moreover, the effect of different power take-off (PTO) parameters on the performance of the WEC's produced power is investigated. A rational PTO damping and stiffness coefficient related with the power production is proposed. The effect of viscous load modeling and hydrodynamic coupling on the response of the combined structure is studied. Motions and dynamic responses of the combined structure under typical operational and extreme sea states are presented. The relative heave, semisubmersible pitch, horizontal contact force and PTO damping force of the combined structure, have an almost linear relationship with the wave height, while WEC surge and PTO produced power have an almost quadratic relationship under regular waves. Finally, a survival mode is examined in order to ensure combined concepts survivability in extreme wave conditions.

Potential Flow Theory
For large-diameter elements, the existence of the structure affects the characteristics in the surrounding flow field. For the calculation of wave load, the three-dimensional potential flow theory is adopted [19,20], which generally assumes an ideal fluid and incompressible flow. The Laplace equation can be obtained with a velocity potential Φ based on the conservation of fluid mass and momentum, and is expressed as follows: The free wave surface, structure surface and subsea conditions of Φ can be expressed as below: where U j and n j represent the generalized velocity and generalized normal vector of the corresponding motion mode of a point on the structure surface, respectively.
In the potential flow theory, the control equations and boundary conditions can be linearized, and the linear superposition principle can be used to decompose the velocity potential as below: φ(x, y, z, t) = φ(x, y, z, t) i + φ(x, y, z, t) d + φ(x, y, z, t) r (5) where φ(x, y, z,t) i represents the incident wave velocity potential, φ(x, y, z,t) d is the diffraction potential and φ(x, y, z,t) r represents the radiation potential. If the floating platform exhibits simple harmonic motion at the balanced position in static water, the velocity potential of the floating platform can be expressed as follows: φ(x, y, z, t) = Re ϕ(x, y, z) · exp(−iωt) (6) ϕ(x, y, z, t) = ϕ i (x, y, z, t) + ϕ d (x, y, z, t) + ϕ r (x, y, z, t) The incident wave velocity potential can be expressed as follows: cosh kh · exp(ik(x · cos β + y · sin β)) (8) where A, g, h, k and β represent the wave amplitude, gravitational acceleration, water depth, wavenumber and wave propagation direction, respectively. The radiation potential that can be solved by the following governing equations and boundary conditions can be expressed as: The control equation, free wave surface condition, structure surface condition, subsea conditions and r → ∞ can be expressed as follows: where x j is the motion amplitude in the motion mode of j, and ϕ r j represents the radiation potential in the motion mode of j.
Finally, the Green's function method is used to calculate the diffraction potential and is no elaborated further here.

Morison Equation (Cylinders of the Semisubmersible Platform)
For small-diameter offshore structures, the Morrison equation is mainly used to calculate the hydrodynamic loads [19,21]. When the diameter of the structure is less than or equal to 0.15 times the wavelength, the existence of the structure has small effect on the wave field. The wave force applied on the structure is the sum of the drag force and the inertia force, which can be expressed as follows: where C d represents the drag force coefficient. In the paper, C d = 1.2 is selected due to H/d ≤ 0.2 and d/L ≥ 0.2; H, d and L represent wave height, water depth and wave length, respectively; µ, D, and ρ are the incoming flow velocity, structure diameter and fluid density, respectively; and f d represents drag force on a unit height of the structure.

Heave Mode Viscous Load Effects on the WEC
The viscosity can significantly affect the dynamic responses of the combined structure (especially the estimation of the PTO produced power), but the effect of viscous damping on the WEC cannot be calculated based on potential flow theory, which has a great difference with physical truth. Therefore, the viscous effect of the heave mode of the WEC has been additionally included, which is calculated as follows [22,23]: where M 33 , m 33 and K 33 represent the mass, added mass at heave natural frequency and hydrostatic stiffness coefficient of the WEC in the heave direction and their values are 540.982 t, 422.3 t and 1515.4 kN/m, respectively; F viscous , D critical and η are the viscous force, the critical damping and damping coefficient based on the viscous effects of the heave mode of the WEC. Based on the reference [24], a damping coefficient of 6-10% in heave direction is recommended using experimental study for different cylindrical semisubmersible platforms. Therefore, 8% of damping coefficient for WEC in heave direction is used in our simulation.

Equation of Motion of the Combined Structure
The simulation of the combined structure is carried out in the time domain. The hydrodynamic properties of the numerical modeling, such as the linear excitation forces, added mass and potential damping as well as the coupling terms of the two bodies, are obtained. In addition, the mechanical coupling (PTO of WEC) between the two bodies is modeled with both linear damper and spring. The mooring lines are also regarded as linear springs. Therefore, the equation of motion for the analysis of the combined structure can be expressed as follows [9]: .. ..
where the subscripts 1 and 2 are denoting the semisubmersible platform and the WEC, respectively; M and m represent the structural mass matrix and the added mass, respectively; x, .
x and ..
x represent the displacement, velocity and acceleration of the bodies; h(t-τ) represents the retardation function; R ij is the restoring coefficient matrix; f wind (t) and f wave (t) represent wind loads on the rotor and wave forces, respectively; f viscous 1 (t) is the viscous force based on Morison equation applied on the semisubmersible and f viscous 2 (t) is the viscous load applied on the WEC (Equation (21)); f inter f ace (t) is the interface forces between the semisubmersible platform and the WEC including horizontal contact forces and vertical friction forces; and F PTO represents the PTO force.

Design Parameters of the 5-MW Semisubmersible Wind Turbine and of the WEC
The combined model proposed in the present paper consists of a braceless 5-MW semisubmersible wind turbine and a heave-type WEC. The conceptual model and the main design parameters of the 5-MW semisubmersible wind turbine and the WEC are shown in Figure 1 and Table 1, respectively. In the model, the WEC can be moved up and down relative to the center column of the semisubmersible, and the wave energy is produced by the PTO at the same time. The relative movement is constrained by the mechanical connection between the center column and the WEC in a direction perpendicular to the center column [25]. As a result, the utilization efficiency of marine renewable energy is improved through shared power generation platforms and shared transmission systems, effectively reducing the unit power generation cost [6,13]. For the numerical simulation, viscous forces of the columns of the semisubmersible are considered as viscous drag loads using Morison equation in the horizontal direction. Center of mass (0, 0, −1 m) Figure 1. Conceptual sketch of the combined structure.

Power-Take-Off (PTO) System
The combined structure has been modeled as two rigid bodies consisting of a semisubmersible platform and a WEC, which are connected by a PTO. The PTO has been simplified as a system of heave-direction linear dampers (BPTO) and linear springs (KPTO) and is modeled with the use of corresponding fender features in ANSYS/AQWA ( Figure 2).  Center of mass (0, 0, −1 m)

Power-Take-Off (PTO) System
The combined structure has been modeled as two rigid bodies consisting of a semisubmersible platform and a WEC, which are connected by a PTO. The PTO has been simplified as a system of heave-direction linear dampers (B PTO ) and linear springs (K PTO ) and is modeled with the use of corresponding fender features in ANSYS/AQWA ( Figure 2). For the PTO system [9], the PTO force includes the linear damping force and the linear spring force, which can be expressed as: where BPTO and KPTO are the linear damping stiffness coefficient and the linear spring stiffness coefficient, respectively. Moreover, the produced power by the relative heave motion between the column of the semisubmersible and the WEC is calculated as follows:

Preliminary determination of structural design parameters of the WEC
The modification of the structural parameters of the WEC will directly affect its hydrodynamic performance, resulting in changes in the produced wave energy efficiency. By evaluating the motion responses of the combined structure, a rational outer diameter and draft can be obtained.
Five different outer diameters of the WEC ranging from 12.0 m to 20.0 m are examined. To obtain the best WEC in terms of PTO produced power, WEC's cost and structure's safety, some factors such as gravity of the WEC and the horizontal contact forces between the column of the semisubmersible and the WEC have been considered. The values of those parameters are calculated under a regular wave with wave height 2 m and wave period 9 s. It is noted that wind loads are not considered here. To reduce the cost and increase safety, the best selection of outer diameter of the WEC can be obtained by Equation (24) and Figure 3. It is noted that we did not solve the mathematical optimization problem in this paper. Based on the curves in Figure 3, a rational value of 16 m is selected for the rest of the paper since all the curves are close to their peak for 16 m.
where θ reflects a performance indicator of the WEC; P, Fh and G represent the PTO produced power, horizontal contact force and WEC weight, respectively, and α and β are related to the cost and safety of the WEC. With the comparison of different parameters with various groups (α = 0.4, β = 0.6; α = 0.3, β = 0.7; α = 0.2, β = 0.8), it is found that the best selection will be achieved when the outer diameter is 16 m. In this paper, we assume that cost reduction (related to β) is For the PTO system [9], the PTO force includes the linear damping force and the linear spring force, which can be expressed as: where B PTO and K PTO are the linear damping stiffness coefficient and the linear spring stiffness coefficient, respectively. Moreover, the produced power by the relative heave motion between the column of the semisubmersible and the WEC is calculated as follows:

Preliminary Determination of Structural Design Parameters of the WEC
The modification of the structural parameters of the WEC will directly affect its hydrodynamic performance, resulting in changes in the produced wave energy efficiency. By evaluating the motion responses of the combined structure, a rational outer diameter and draft can be obtained.
Five different outer diameters of the WEC ranging from 12.0 m to 20.0 m are examined. To obtain the best WEC in terms of PTO produced power, WEC's cost and structure's safety, some factors such as gravity of the WEC and the horizontal contact forces between the column of the semisubmersible and the WEC have been considered. The values of those parameters are calculated under a regular wave with wave height 2 m and wave period 9 s. It is noted that wind loads are not considered here. To reduce the cost and increase safety, the best selection of outer diameter of the WEC can be obtained by Equation (24) and Figure 3. It is noted that we did not solve the mathematical optimization problem in this paper. Based on the curves in Figure 3, a rational value of 16 m is selected for the rest of the paper since all the curves are close to their peak for 16 m.
where θ reflects a performance indicator of the WEC; P, F h and G represent the PTO produced power, horizontal contact force and WEC weight, respectively, and α and β are related to the cost and safety of the WEC. With the comparison of different parameters with various groups (α = 0.4, β = 0.6; α = 0.3, β = 0.7; α = 0.2, β = 0.8), it is found that the best selection will be achieved when the outer diameter is 16 m. In this paper, we assume that cost reduction (related to β) is more important and the safety is less important (related to α) for the combined structure. Therefore, α = 0.3, β = 0.7 is assumed in this study. more important and the safety is less important (related to α) for the combined structure. Therefore, α= 0.3, β = 0.7 is assumed in this study. Moreover, the effect of different drafts of the WEC on the responses of the semisubmersible are presented in Figure 4 and are slightly different when the draft is 3.0 m and 3.5 m. However, considering the possible slamming effects due to water entry and exit, which can destroy the WEC, choosing the value of the draft (3.5 m) is more reasonable. Therefore, the selection of diameter 16 m and draft 3.5 m is followed in the rest of the paper. Moreover, the effect of different drafts of the WEC on the responses of the semisubmersible are presented in Figure 4 and are slightly different when the draft is 3.0 m and 3.5 m. However, considering the possible slamming effects due to water entry and exit, which can destroy the WEC, choosing the value of the draft (3.5 m) is more reasonable. Therefore, the selection of diameter 16 m and draft 3.5 m is followed in the rest of the paper.

Determination of rational BPTO and KPTO parameters
Before determination of linear damping and stiffness coefficients, the effect of different wave periods (BPTO = 2000 kNs/m, KPTO = 100 N/m was arbitrary selected) on the PTO produced power has been investigated under regular waves with wave height 2 m. The values of produced power are shown in Figure 5. It indicates that there is a maximum produced power when the wave period is around 9 s because the WEC was designed to have a natural period in heave around 7-2 s.

Determination of Rational B PTO and K PTO Parameters
Before determination of linear damping and stiffness coefficients, the effect of different wave periods (B PTO = 2000 kNs/m, K PTO = 100 N/m was arbitrary selected) on the PTO produced power has been investigated under regular waves with wave height 2 m. The values of produced power are shown in Figure 5. It indicates that there is a maximum produced power when the wave period is around 9 s because the WEC was designed to have a natural period in heave around 7-2 s.
The B PTO and K PTO parameters can significantly affect the dynamic responses of the combined structure. To illustrate mutual effects of linear damping and stiffness coefficients on the PTO produced power of the WEC, different B PTO and K PTO values are simulated with wave height 2 m and wave period 9 s. The results of produced power are shown in Figure 6; it can be seen that the value of produced power will be the largest when the B PTO coefficient is 1500 kNs/m with an invariable K PTO . The value of produced power increases when the B PTO is smaller than 1500 kNs/m, while the values will decrease when the B PTO is larger than 1500 kNs/m. Moreover, for same B PTO values, produced power decreases as K PTO increases, which is more significant when its value is over 1000 N/m. Therefore, a best B PTO coefficient of 1500 kNs/m and a K PTO coefficient of 1 N/s are selected for the rest of the paper, respectively.

Determination of rational BPTO and KPTO parameters
Before determination of linear damping and stiffness coefficients, the effect of different wave periods (BPTO = 2000 kNs/m, KPTO = 100 N/m was arbitrary selected) on the PTO produced power has been investigated under regular waves with wave height 2 m. The values of produced power are shown in Figure 5. It indicates that there is a maximum produced power when the wave period is around 9 s because the WEC was designed to have a natural period in heave around 7-2 s. The BPTO and KPTO parameters can significantly affect the dynamic responses of the combined structure. To illustrate mutual effects of linear damping and stiffness coefficients on the PTO produced power of the WEC, different BPTO and KPTO values are simulated with wave height 2 m and wave period 9 s. The results of produced power are shown in Figure 6; it can be seen that the value of produced power will be the largest when the BPTO coefficient is 1500 kNs/m with an invariable KPTO. The value of produced power increases when the BPTO is smaller than 1500 kNs/m, while the values will decrease when the BPTO is larger than 1500 kNs/m. Moreover, for same BPTO values, produced power decreases as KPTO increases, which is more significant when its value is over 1000 N/m. Therefore, a best BPTO coefficient of 1500 kNs/m and a KPTO coefficient of 1 N/s are selected for the rest of the paper, respectively.

Viscous effects on the heave motion of the WEC
By ignoring viscous effects, the produced power may be overestimated when there is a large relative motion between the semisubmersible and the WEC. The dynamic responses of the combined structure with or without considering viscous effects have been presented for regular waves with a wave height 2 m and a wave period 9 s. In Figure 7, the dynamic responses considering viscous effects are smaller than those without considering them. The relative heave motion with viscous effect has a peak value when wave period is 11 s, which is different from that without viscous effect. However, the radiation force, PTO damping and produced power, a peak value exists when the period is 9 s. Moreover, the ratio which can quantitatively describe viscous effects on the combined structure between those responses (Equation 25) always increases with increasing period.
where ε is a ratio, and λ and ν represent the dynamic responses with and without viscous effects.

Viscous Effects on the Heave Motion of the WEC
By ignoring viscous effects, the produced power may be overestimated when there is a large relative motion between the semisubmersible and the WEC. The dynamic responses of the combined structure with or without considering viscous effects have been presented for regular waves with a wave height 2 m and a wave period 9 s. In Figure 7, the dynamic responses considering viscous effects are smaller than those without considering them. The relative heave motion with viscous effect has a peak value when wave period is 11 s, which is different from that without viscous effect. However, the radiation force, PTO damping and produced power, a peak value exists when the period is 9 s. Moreover, the ratio which can quantitatively describe viscous effects on the combined structure between those responses (Equation (25)) always increases with increasing period.
where ε is a ratio, and λ and ν represent the dynamic responses with and without viscous effects.

The hydrodynamic coupling effect on the combined structure
In the flow field, in addition to the direct interaction between the semisubmersible and WEC, there is also an interaction between the two bodies due to the radiation and diffraction effects of the combined structure on the flow field. The radiation and diffraction potentials acting on the semisubmersible and WEC can interact when considering hydrodynamic coupling [26,27], which is proved by comparing diffraction forces and radiation forces of the semisubmersible and the WEC with considering hydrodynamic coupling with those without considering hydrodynamic coupling ( Figure 8). It can be seen that the diffraction forces ( Figure 8a) and radiation forces (Figure 8b) acting on the semisubmersible have same trends for different wave periods when considering hydrodynamic coupling. There is a little increase of the diffraction force and radiation force acted on the semisubmersible considering hydrodynamic coupling. The largest difference for diffraction force and radiation force is approximately 3.8% and 3.6% larger than those without considering hydrodynamic coupling at wave period of 11 s. However, for WEC, the effects of the hydrodynamic coupling are significant on the diffraction forces ( Figure 8c) and radiation forces (Figure 8d). And the diffraction force acted on the WEC considering hydrodynamic coupling is approximately 52.2% larger than that without considering hydrodynamic coupling at the wave period of 8 s and the radiation force is 55.4% larger than that without considering hydrodynamic coupling at the period of 11 s. Therefore, motion and dynamic responses of the WEC should be put more emphasis compared with the semisubmersible for the combined structure.
As hydrodynamic loads with considering hydrodynamic coupling are different from those without considering hydrodynamic coupling the relevant responses of the combined structure with wave height 2 m are presented and studied. The WEC surge motion (Figure 9a) and the relative heave motion (Figure 9b) considering hydrodynamic coupling are not always larger compared to

The Hydrodynamic Coupling Effect on the Combined Structure
In the flow field, in addition to the direct interaction between the semisubmersible and WEC, there is also an interaction between the two bodies due to the radiation and diffraction effects of the combined structure on the flow field. The radiation and diffraction potentials acting on the semisubmersible and WEC can interact when considering hydrodynamic coupling [26,27], which is proved by comparing diffraction forces and radiation forces of the semisubmersible and the WEC with considering hydrodynamic coupling with those without considering hydrodynamic coupling ( Figure 8). It can be seen that the diffraction forces ( Figure 8a) and radiation forces (Figure 8b) acting on the semisubmersible have same trends for different wave periods when considering hydrodynamic coupling. There is a little increase of the diffraction force and radiation force acted on the semisubmersible considering hydrodynamic coupling. The largest difference for diffraction force and radiation force is approximately 3.8% and 3.6% larger than those without considering hydrodynamic coupling at wave period of 11 s. However, for WEC, the effects of the hydrodynamic coupling are significant on the diffraction forces ( Figure 8c) and radiation forces (Figure 8d). And the diffraction force acted on the WEC considering hydrodynamic coupling is approximately 52.2% larger than that without considering hydrodynamic coupling at the wave period of 8 s and the radiation force is 55.4% larger than that without considering hydrodynamic coupling at the period of 11 s. Therefore, motion and dynamic responses of the WEC should be put more emphasis compared with the semisubmersible for the combined structure.
As hydrodynamic loads with considering hydrodynamic coupling are different from those without considering hydrodynamic coupling the relevant responses of the combined structure with wave height 2 m are presented and studied. The WEC surge motion (Figure 9a) and the relative heave motion (Figure 9b) considering hydrodynamic coupling are not always larger compared to those without considering hydrodynamic coupling. Additionally, the amplitudes of the relative heave motion (Figure 9b) and WEC pitch motion (Figure 9c) have similar trends, but WEC pitch motion considering hydrodynamic coupling is always larger than the case without considering hydrodynamic coupling. Figure 9d shows the wave period considering hydrodynamic coupling (approximately 9 s) is larger than that without considering hydrodynamic coupling (approximately 8 s). The largest produced power considering hydrodynamic coupling is smaller than that without considering hydrodynamic coupling. Same conclusions with PTO damping force (Figure 9e)  those without considering hydrodynamic coupling. Additionally, the amplitudes of the relative heave motion (Figure 9b) and WEC pitch motion (Figure 9c) have similar trends, but WEC pitch motion considering hydrodynamic coupling is always larger than the case without considering hydrodynamic coupling. Figure 9d shows the wave period considering hydrodynamic coupling (approximately 9 s) is larger than that without considering hydrodynamic coupling (approximately 8 s). The largest produced power considering hydrodynamic coupling is smaller than that without considering hydrodynamic coupling. Same conclusions with PTO damping force (Figure 9e) are presented.   Figure  10) without wind excitation loads [28]. In Figure 10a, the maximum relative heave velocity of the combined structure is smaller than 0.375 m/s while the minimum absolute value is larger than 0.375 m/s, which is an explanation why the peaks of produced power (Figure 10c) are not obtaining same values (Equation 23). In Figure 10b, the PTO damping force induced by fenders in AQWA mode is presented. The phase of the forces is similar with relative heave, which means the damping forces  Figure 10) without wind excitation loads [28]. In Figure 10a, the maximum relative heave velocity of the combined structure is smaller than 0.375 m/s while the minimum absolute value is larger than 0.375 m/s, which is an explanation why the peaks of produced power (Figure 10c) are not obtaining same values (Equation (23)). In Figure 10b, the PTO damping force induced by fenders in AQWA mode is presented. The phase of the forces is similar with relative heave, which means the damping forces are negative or positive when the WEC is up or down. Moreover, the capture width ratio, which is used to evaluate the WEC capability of absorbed wave power, is calculated by Equations (26) and (27) [29,30] as follows: where P wave and P PTO are the average wave power per unit width and average PTO produced power, respectively, and σ denotes the capture width ratio.
structure proposed in the present paper. Figure 12 shows that the amplitudes of the responses, such as the relative heave, semisubmersible pitch, horizontal contact forces and PTO damping forces of the combined structure, almost linearly increase with the wave height, while the amplitudes of semisubmersible surge and PTO produced power have an approximately quadratic relationship with the wave height. To obtain the characteristics of the combined structure, the responses of the proposed combined structure are compared with similar work done by Wan et al. [9], Gao el al. [31] and Ren et al. [32] ( Figure 11). Attention should be paid that in our work a semisubmersible platform is used while a spar platform was used in Wan and Gao's work and a tension leg platform (TLP) was used Ren's work. In Figure 11a,b, platform surge and WEC heave have similar tendency compared with Wan et al.'s work. The wave energy capture width ratio is compared with STC [9,31] and [32] in Figure 11d. The difference of the capture width ratio is caused by two reasons. The first one is the difference of PTO damping and stiffness coefficient, which result in the difference for the best wave period and produced power. The second reason is that the relative heave motions between the semisubmersible platform and WEC are much less than those between the spar and the TLP platforms and their WECs. This is caused by the heave motion difference of the spar, TLP and semisubmersible platforms. According to Equation (23), produced power of the combined structure will be smaller compared with those power produced by STC [9,31] and [32]. Special attention should be made that Wan's results are based on the irregular wave case. It is stated that the comparisons are presented for having an idea about the well comparable response of the combined structure proposed in the present paper. Figure 12 shows that the amplitudes of the responses, such as the relative heave, semisubmersible pitch, horizontal contact forces and PTO damping forces of the combined structure, almost linearly increase with the wave height, while the amplitudes of semisubmersible surge and PTO produced power have an approximately quadratic relationship with the wave height.

Irregular wave and wind conditions
The sea conditions of certain sea areas in Guangdong Province, China whose location of the offshore wind farm is presented in Figure 13 have been selected, as shown in Table 2. In the paper, the aerodynamic loads from wind flow are simplified as a mean constant thrust load based on NREL 5 MW wind turbine thrust force-wind speed curve ( Figure 14). The JONSWAP spectrum is used to simulate the irregular wave in this section. Comparisons of the time series and power spectral density (PSD) of motion and dynamic responses for LC1 (Hs = 2 m, Tp = 9 s) and LC2 (Hs = 2 m, Tp = 9 s and Uwind = 17 m/s) are presented in Figures 15 and 16, respectively. Figure 15 shows that the semisubmersible surge and pitch motions have big differences due to wind loads. However, the values of the relative heave, PTO damping force and PTO produced power have slight differences.

Irregular Wave and Wind Conditions
The sea conditions of certain sea areas in Guangdong Province, China whose location of the offshore wind farm is presented in Figure 13 have been selected, as shown in Table 2. In the paper, the aerodynamic loads from wind flow are simplified as a mean constant thrust load based on NREL 5 MW wind turbine thrust force-wind speed curve ( Figure 14). The JONSWAP spectrum is used to simulate the irregular wave in this section. Comparisons of the time series and power spectral density (PSD) of motion and dynamic responses for LC1 (H s = 2 m, T p = 9 s) and LC2 (H s = 2 m, T p = 9 s and U wind = 17 m/s) are presented in Figures 15 and 16, respectively. Figure 15 shows that the semisubmersible surge and pitch motions have big differences due to wind loads. However, the values of the relative heave, PTO damping force and PTO produced power have slight differences.

Irregular wave and wind conditions
The sea conditions of certain sea areas in Guangdong Province, China whose location of the offshore wind farm is presented in Figure 13 have been selected, as shown in Table 2. In the paper, the aerodynamic loads from wind flow are simplified as a mean constant thrust load based on NREL 5 MW wind turbine thrust force-wind speed curve (Figure 14). The JONSWAP spectrum is used to simulate the irregular wave in this section. Comparisons of the time series and power spectral density (PSD) of motion and dynamic responses for LC1 (Hs = 2 m, Tp = 9 s) and LC2 (Hs = 2 m, Tp = 9 s and Uwind = 17 m/s) are presented in Figures 15 and 16, respectively. Figure 15 shows that the semisubmersible surge and pitch motions have big differences due to wind loads. However, the values of the relative heave, PTO damping force and PTO produced power have slight differences. Figure 13. Location of the reference offshore site. Figure 13. Location of the reference offshore site.   The natural frequencies (surge, heave and pitch) of the semisubmersible could be identified through free decay tests (Table 3). In Figure 16, PSD of motion and dynamic responses of the combined structure under LC1 have same trends compared with those under LC2. In Figure 16a, a resonance is around 0.073 rad/s, which is contributed by surge natural frequency. At this frequency, the surge motion in LC2 is enlarged due to additional wind load compared with LC1, which can also be confirmed from Figure 15a. The second peak (around 0.7 rad/s) in Figure 16a presents the contribution from the wave frequency to the surge response. Figure 16c shows that the pitch motion response is dominant by pitch resonance frequency (0.218 rad/s). Resonances from surge natural frequency (0.073 rad/s) and wave frequency (0.7 rad/s) could also be identified from pitch motion response. For heave motion in Figure 16b, the wave frequency response is dominated in the range between 0.500 rad/s to 0.800 rad/s, in which WEC is designed to capture the energy from incident wave. For PTO damping force (Figure 16d), the surge natural frequency and wave natural frequency could be clearly identified from load responses. Table 3. Natural frequencies (rad/s) of the semisubmersible.

Surge
Heave Pitch 0.07545 0.2543 0.2108 Figure 14. Rotor thrust curves at different wind speeds. The natural frequencies (surge, heave and pitch) of the semisubmersible could be i through free decay tests (Table 3). In Figure 16, PSD of motion and dynamic response combined structure under LC1 have same trends compared with those under LC2. In Figu resonance is around 0.073 rad/s, which is contributed by surge natural frequency. At this fr the surge motion in LC2 is enlarged due to additional wind load compared with LC1, w also be confirmed from Figure 15a. The second peak (around 0.7 rad/s) in Figure 16a pre contribution from the wave frequency to the surge response. Figure 16c shows that the pitc wave. For PTO damping force (Figure 16d), the surge natural frequency and wave natural frequency could be clearly identified from load responses.  The statistical values of the motion and dynamic responses of the combined structure under LC3 are shown in Figure 17. In Figure 17a, we can see that the maximum surge is less than 3.5 m and maximum pitch is less than 3.5 degrees. The maximum relative heave is larger than 1.25 m, while the absolute value of the minimum relative heave is smaller than 1.25 m, which may be The statistical values of the motion and dynamic responses of the combined structure under LC3 are shown in Figure 17. In Figure 17a, we can see that the maximum surge is less than 3.5 m and maximum pitch is less than 3.5 degrees. The maximum relative heave is larger than 1.25 m, while the absolute value of the minimum relative heave is smaller than 1.25 m, which may be caused by the lift effect induced by strong wave action [22]. Compared with LC1, it finds that large mean value of surge and pitch motion are caused by wind loads and positive surge and pitch moment presents in LC3. Figure 17b shows that the horizontal contact force between the semisubmersible and the WEC is larger than the wave force acting on the WEC due to the relative acceleration of the two bodies. The mooring line 1 tension force is larger than that of mooring line 2 because the drift of the combined structure in downwind direction due to the wind load. The mean PTO produced power is around 50 kW with maximum power around 1245 kW under LC3. Moreover, in Figure 18, PSD of motion responses (surge, relative heave and pitch) and dynamic responses (horizontal contact force, PTO damping force, mooring line 1 tension and mooring line 2 tension) of the combined structure with considering hydrodynamic coupling have similar trend compared with those without considering hydrodynamic coupling. For surge motion (Figure 18a), the surge resonant motion is dominant at 0.073 rad/s, which is the surge natural frequency, followed by the wave frequency response in the range 0.400 rad/s to 0.800 rad/s. A large decrease of the response at the surge frequency exhibits if the hydrodynamic coupling is considered. Similar as LC1 and LC2 (Figure 16b), the most significant contribution to the relative heave motion comes from the wave frequency (Figure 18b). Considering hydrodynamic coupling will give a small increase in heave motion response near wave frequency, which is similar as in surge motion response. This enhancement of the surge and pitch responses near wave frequency can also be identified with wave period of 10 s in Figure 9a,b. For pitch motion (Figure 18c), the wave frequency response and pitch resonant response are very comparable, meanwhile, surge resonance response also visible in this motion. The pitch response considering hydrodynamic coupling at the surge natural frequency and wave frequency is larger than that without considering hydrodynamic coupling, opposite influence can be identified at pitch natural frequency. Figure 18d,e shows that the wave frequency responses are also dominant in the same range for horizontal contact force and damping force and the hydrodynamic coupling will contribute the increase of the response in the wave frequency range. For mooring line tension force (Figure 18g,h)  Moreover, in Figure 18, PSD of motion responses (surge, relative heave and pitch) and dynamic responses (horizontal contact force, PTO damping force, mooring line 1 tension and mooring line 2 tension) of the combined structure with considering hydrodynamic coupling have similar trend compared with those without considering hydrodynamic coupling. For surge motion (Figure 18a), the surge resonant motion is dominant at 0.073 rad/s, which is the surge natural frequency, followed by the wave frequency response in the range 0.400 rad/s to 0.800 rad/s. A large decrease of the response at the surge frequency exhibits if the hydrodynamic coupling is considered. Similar as LC1 and LC2 (Figure 16b), the most significant contribution to the relative heave motion comes from the wave frequency (Figure 18b). Considering hydrodynamic coupling will give a small increase in heave motion response near wave frequency, which is similar as in surge motion response. This enhancement of the surge and pitch responses near wave frequency can also be identified with wave period of 10 s in Figure 9a,b. For pitch motion (Figure 18c), the wave frequency response and pitch resonant response are very comparable, meanwhile, surge resonance response also visible in this motion. The pitch response considering hydrodynamic coupling at the surge natural frequency and wave frequency is larger than that without considering hydrodynamic coupling, opposite influence can be identified at pitch natural frequency. Figure 18d,e shows that the wave frequency responses are also dominant in the same range for horizontal contact force and damping force and the hydrodynamic coupling will contribute the increase of the response in the wave frequency range. For mooring line tension force (Figure 18g,h)

Dynamic Responses of the Combined Structure in Extreme Sea Conditions
For the LC4 with extreme sea conditions (wave height 8.6 m, wave period 15 s and wind speed of 31.2 m/s), the PTO system is released, the wind turbine is parked, and the WEC is locked mechanically to the semisubmersible at the mean sea level as a survival mode. In this case, the two bodies are locked and can move together, which can be confirmed from the relative heave motion between WEC and semisubmersible in Figure 19. The statistical value from LC4 (mean, maxi mum, minimum and STD) are presented in Figure 20. In this load case, maximum heave motion of the combined structure is approximate 6 m (Figure 20a), in which WEC may bear huge loads due to slamming forces. The mean surge and pitch responses of the combined structure in LC4 is much less than those responses in LC3 (Figure 17a). Particularly, the maximum and minimum pitch response is much less than those values in LC3 for operating conditions because the thrust load in LC4 is much smaller than the thrust load in LC3 ( Figure 14). Figure 20b shows that the maximum wave force of the combined structure is smaller than the vertical contact force between the semisubmersible column and the WEC, while the horizontal contact force is far smaller than the wave force and vertical contact force, which proves the locked mode of the combined structure results to a reasonable response.

Dynamic responses of the combined structure in extreme sea conditions
For the LC4 with extreme sea conditions (wave height 8.6 m, wave period 15 s and wind speed of 31.2 m/s), the PTO system is released, the wind turbine is parked, and the WEC is locked mechanically to the semisubmersible at the mean sea level as a survival mode. In this case, the two bodies are locked and can move together, which can be confirmed from the relative heave motion between WEC and semisubmersible in Figure 19. The statistical value from LC4 (mean, maxi mum, minimum and STD) are presented in Figure 20. In this load case, maximum heave motion of the combined structure is approximate 6 m (Figure 20a), in which WEC may bear huge loads due to slamming forces. The mean surge and pitch responses of the combined structure in LC4 is much less than those responses in LC3 (Figure 17a). Particularly, the maximum and minimum pitch response is much less than those values in LC3 for operating conditions because the thrust load in LC4 is much smaller than the thrust load in LC3 ( Figure 14). Figure 20b shows that the maximum wave force of the combined structure is smaller than the vertical contact force between the semisubmersible column and the WEC, while the horizontal contact force is far smaller than the wave force and vertical contact force, which proves the locked mode of the combined structure results to a reasonable response.

Dynamic responses of the combined structure in extreme sea conditions
For the LC4 with extreme sea conditions (wave height 8.6 m, wave period 15 s and wind speed of 31.2 m/s), the PTO system is released, the wind turbine is parked, and the WEC is locked mechanically to the semisubmersible at the mean sea level as a survival mode. In this case, the two bodies are locked and can move together, which can be confirmed from the relative heave motion between WEC and semisubmersible in Figure 19. The statistical value from LC4 (mean, maxi mum, minimum and STD) are presented in Figure 20. In this load case, maximum heave motion of the combined structure is approximate 6 m (Figure 20a), in which WEC may bear huge loads due to slamming forces. The mean surge and pitch responses of the combined structure in LC4 is much less than those responses in LC3 (Figure 17a). Particularly, the maximum and minimum pitch response is much less than those values in LC3 for operating conditions because the thrust load in LC4 is much smaller than the thrust load in LC3 ( Figure 14). Figure 20b shows that the maximum wave force of the combined structure is smaller than the vertical contact force between the semisubmersible column and the WEC, while the horizontal contact force is far smaller than the wave force and vertical contact force, which proves the locked mode of the combined structure results to a reasonable response.    Figure 21 presents the PSD of the motion response and mooring line forces of the combined structure and pure semisubmersible platform under the extreme sea condition LC4. PSD of surge of the combined structure have same trends compared with pure semisubmersible platform and the response is dominant by the wave frequency (Figure 21a). For heave motion (Figure 21b), the response of pure semisubmersible is dominated by the wave frequency. However, for the combined structure, the response around the frequency of 0.35 rad/s is significantly larger because this frequency is the heave natural frequency of the combined structure and is near to the range of wave frequency. The shifting of the heave natural frequency of the combined structure from the pure semisubmersible (0.254 rad/s) is due to the attachment of the WEC to the semisubmersible. The surge resonant motion of the combined structure (Figure 21c) is dominant at 0.21 rad/s, which is the pitch natural frequency, followed by the wave frequency response. But the surge motion of the pure semisubmersible is dominated by the wave frequency response. Figure 21d,e shows that PSD of mooring line 1 tension and mooring line 2 tension of the combined structure are not always identical with only pure semisubmersible platform in high peak frequency region, which is caused by the coupling effect of heave resonance and incident wave for the combined structure. The most significant contribution to the mooring line 1 tension for both the combined structure and the semisubmersible platform comes from the wave frequency range. A clear peak response from the heave natural frequency for the combined structure could be identified for mooring line 1 tension. For mooring line 2 (Figure 21e), the most significant contribution to the mooring line tension for the pure semisubmersible comes from the wave frequency range from 0.300 rad/s to 0.600 rad/s, while the low-frequency response does not contribute significantly to the total response. However, for the combined structure, the heave resonance response at approximately 0.35 rad/s dominates the total response for mooring line 2.  Figure 21 presents the PSD of the motion response and mooring line forces of the combined structure and pure semisubmersible platform under the extreme sea condition LC4. PSD of surge of the combined structure have same trends compared with pure semisubmersible platform and the response is dominant by the wave frequency (Figure 21a). For heave motion (Figure 21b), the response of pure semisubmersible is dominated by the wave frequency. However, for the combined structure, the response around the frequency of 0.35 rad/s is significantly larger because this frequency is the heave natural frequency of the combined structure and is near to the range of wave frequency. The shifting of the heave natural frequency of the combined structure from the pure semisubmersible (0.254 rad/s) is due to the attachment of the WEC to the semisubmersible. The surge resonant motion of the combined structure (Figure 21c) is dominant at 0.21 rad/s, which is the pitch natural frequency, followed by the wave frequency response. But the surge motion of the pure semisubmersible is dominated by the wave frequency response. Figure 21d,e shows that PSD of mooring line 1 tension and mooring line 2 tension of the combined structure are not always identical with only pure semisubmersible platform in high peak frequency region, which is caused by the coupling effect of heave resonance and incident wave for the combined structure. The most significant contribution to the mooring line 1 tension for both the combined structure and the semisubmersible platform comes from the wave frequency range. A clear peak response from the heave natural frequency for the combined structure could be identified for mooring line 1 tension. For mooring line 2 (Figure 21e), the most significant contribution to the mooring line tension for the pure semisubmersible comes from the wave frequency range from 0.300 rad/s to 0.600 rad/s, while the low-frequency response does not contribute significantly to the total response. However, for the combined structure, the heave resonance response at approximately 0.35 rad/s dominates the total response for mooring line 2.

Conclusions
In this paper, a combined marine renewable energy structure consisting of a 5-MW br floating wind turbine and a heave-type wave energy converter (WEC) is presented. A num model of the combined structure that is capable of simulating its motion and dynamic res under different sea states was developed and used. Different types of analysis have been per for examining the hydrodynamic response of the combined structure. The main resu summarized as follows: The outer diameter of the WEC is selected to be 16 m, in term of its economy, safety an produced power. In order to avoid the slamming phenomena of the WEC as far as possible response, the draft of the WEC with 3.5 m is selected.
The PTO K parameter has a significant impact on the performance of the WEC when its v over 1000 N/m. A preliminary best value of the PTO B parameter for the PTO of the com structure is selected to be 1500 kNs/m with the perspective of PTO produced power. The trends of hydrodynamic loads of the WEC with considering hydrodynamic coupli change compared to those without considering hydrodynamic coupling, but responses combined structure have the same trends. The viscosity of the heave of the WEC will h influence on the PTO produced power and dynamic responses of the combined structure.
Under regular waves, the amplitudes of the responses, such as the relative semisubmersible pitch, horizontal contact force and PTO damping force of the combined str have an almost linear relationship with the wave height, while the amplitudes of WEC sur PTO produced power have an almost quadratic relationship.
Under irregular waves, the surge and pitch of the semisubmersible with considering c thrust wind have great effects compared to those without considering constant thrust win there are slight effects for relative heave, damping force and produced power of the com structure. Moreover, the maximum relative heave of the combined structure is larger th absolute value of the minimum relative heave, and the wave forces acting on the WEC transferred to the semisubmersible by means of their contacted fenders. The PSD of responses and dynamic responses of the combined structure with considering hydrod coupling have similar trends compared with those without considering hydrodynamic coupl Under the extreme sea conditions, the WEC was mechanically locked to the semisubm platform to withstand the large impact force on WEC as a survival mode. Heave natural fre shifting could be identified and may coincide with wave frequency. The resonance in heave and mooring line 2 tension could be identified near the heave natural frequency. Other s models should be proposed and investigated in future work. Experiments and validation integrated structure will also be carried out in the future.

Conclusions
In this paper, a combined marine renewable energy structure consisting of a 5-MW braceless floating wind turbine and a heave-type wave energy converter (WEC) is presented. A numerical model of the combined structure that is capable of simulating its motion and dynamic responses under different sea states was developed and used. Different types of analysis have been performed for examining the hydrodynamic response of the combined structure. The main results are summarized as follows: The outer diameter of the WEC is selected to be 16 m, in term of its economy, safety and PTO produced power. In order to avoid the slamming phenomena of the WEC as far as possible motion response, the draft of the WEC with 3.5 m is selected.
The K PTO parameter has a significant impact on the performance of the WEC when its value is over 1000 N/m. A preliminary best value of the B PTO parameter for the PTO of the combined structure is selected to be 1500 kNs/m with the perspective of PTO produced power.
The trends of hydrodynamic loads of the WEC with considering hydrodynamic coupling will change compared to those without considering hydrodynamic coupling, but responses of the combined structure have the same trends. The viscosity of the heave of the WEC will have an influence on the PTO produced power and dynamic responses of the combined structure.
Under regular waves, the amplitudes of the responses, such as the relative heave, semisubmersible pitch, horizontal contact force and PTO damping force of the combined structure, have an almost linear relationship with the wave height, while the amplitudes of WEC surge and PTO produced power have an almost quadratic relationship.
Under irregular waves, the surge and pitch of the semisubmersible with considering constant thrust wind have great effects compared to those without considering constant thrust wind, but there are slight effects for relative heave, damping force and produced power of the combined structure. Moreover, the maximum relative heave of the combined structure is larger than the absolute value of the minimum relative heave, and the wave forces acting on the WEC can be transferred to the semisubmersible by means of their contacted fenders. The PSD of motion responses and dynamic responses of the combined structure with considering hydrodynamic coupling have similar trends compared with those without considering hydrodynamic coupling.
Under the extreme sea conditions, the WEC was mechanically locked to the semisubmersible platform to withstand the large impact force on WEC as a survival mode. Heave natural frequency shifting could be identified and may coincide with wave frequency. The resonance in heave motion and mooring line 2 tension could be identified near the heave natural frequency. Other survival models should be proposed and investigated in future work. Experiments and validation for this integrated structure will also be carried out in the future. This work is also supported by the international collaboration and exchange program from the NSFC-RCUK/EPSRC with grant No. 51761135011. This work is also partially supported by LiaoNing Revitalization Talents Program (XLYC1807208), State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University and the Fundamental Research Funds for the Central University (DUT19GJ209).