Dynamic Characteristics of an Offshore Wind Turbine with Tripod Suction Buckets via Full-Scale Testing

.e dynamic characteristics of an offshore wind turbine with tripod suction buckets are investigated through finite element analysis and full-scale experiments. In finite element analysis, an integrated framework is suggested to create a simple yet accurate high fidelity model. .e integrated framework accounts for not only the strain dependency of the soil but also for all dynamics in the seabed, including those of the soil, suction bucket skirt, and cap. Hence, themodel accurately describes the coupling effect of translational and rotational motions of the seabed..e prediction results are compared to the experimental results obtained via full-scale testing in four stages during construction and in several operational conditions. .e comparison shows that the stiffness of the suction bucket cap and strain dependency of the soil play a significant role in predicting natural frequency, suggesting that these two factors should be considered in finite element analysis for the accurate prediction of dynamic responses of an offshore wind conversion system. Moreover, dynamic analysis of the strain and acceleration measured during operational conditions shows that strain is more robust than acceleration with regard to the characterization of the overall dynamics of an offshore wind conversion system because the natural frequency of an offshore wind turbine is very low. It can be inferred that the measurement of strain is a more effective way to monitor the long-term evolution of dynamic characteristics. .e suggested integrated framework and measurement campaign are useful not only to avoid conservatism that may incur additional costs during load calculation and design phases but also to establish an intelligent operation and maintenance strategy with a novel sensing technique.


Introduction
Renewable and sustainable energies are attractive alternatives for mitigating issues related to the consumption of fossil fuels, including the volatility of oil prices, emission of carbon dioxide, and aggravation of air quality due to fine dust emission. Among renewable and sustainable energies, offshore wind power shows a variety of advantages including high energy density, low turbulence, and low wind shear [1]. Technological advances are also making wind energy competitive from the economic perspective [2][3][4][5][6][7][8][9]. In 2018, 409 new offshore wind turbines were commissioned in Europe, which provided an additional capacity of 2,649 MW, and the cumulative capacity of wind farms was 18,499 MW [10]. ere are two trends for offshore wind power. One trend is the increase of turbine rated capacity. Most wind farms under construction are using turbines with capacities higher than 6 MW. e other trend is increasing the distance from the shore to mitigate concerns raised by civil complaints as well as to achieve better wind potential. For offshore wind farms under construction in Europe in 2018, the mean distance from the shore is 33 km [10]. However, the farther that a wind farm is from the shore, the deeper the water depth of its candidate site. e average water depth of European offshore wind farms under construction in 2018 is 27.1 m [10]. e long distance from the shore and deep water depth hinder the economic feasibility of wind farms.
To increase the economic feasibility of offshore wind farms, many studies have focused on substructures and foundations because the cost for these components accounts for 20-40% of the total cost, depending on water depth [11]. With regard to substructures, space frame substructures, including tripods and jacket structures, are the most promising. ey provide sufficient strength and stiffness at a competitive cost in transitional water depths at which most new wind farms are permitted for installation nowadays. Regarding foundations, bearing capacity is the most important factor in supporting an offshore wind turbine. To estimate bearing capacity, soil-pile interaction should be exactly elucidated depending on types of foundations: gravitational type, pile type, and rigid bucket.
For the gravitational type, the deformation of bulk soil on the surface is of interest in the sense that the most interactions between the foundation and soil occur at the contact surface between the foundation and soil. However, this type is only feasible for shallow water because the construction cost is proportional to water depth and thus expensive for transitional water depths. Consequently, this type is no longer used for MW-class wind turbines installed in transitional water depths. For the pile type, the ratio of flexibility between soil and a pile is of interest because the ratio of flexibility governs soil-structure interaction [12]. Hence, the ratio of flexibility between soil and the pile is important to determine the specifications of the pile, including size and height [13]. Many current wind farms have been installed in transitional water depths with the pile type because of its proven reliability. Regarding the deep and rigid caisson type, interactions occur only below and beside the foundation. e shape of the foundation (aspect ratio) and soil properties are important to evaluate the influence on the interaction. If the foundation is very rigid compared to the soil, coupling of translational and rotational stiffness should be considered.
Interestingly, a suction bucket has the characteristics of both flexible pile and rigid caisson regarding the interaction with soil [14], making the design and construction more difficult and complex. Several factors affect the interaction between the soil and the suction bucket, including dimension, aspect ratio of a suction bucket, soil properties, and the ratio of flexibility between the soil and the suction bucket. Although designing suction buckets is difficult and complex, suction buckets become more promising than pile foundations in the transitional water depth because suction buckets are economic and offer simple and fast installation [15]. Moreover, suction buckets can help overcome several issues associated with the installation of the pile type including low level of vibration, noise, and suspended sediments.
Several studies have been conducted to deploy suction buckets on the foundations of offshore wind turbines [16,17]. Design procedures for the installation of suction buckets in sand and clay have been suggested [18,19]. Experimental studies on the installation of suction buckets have been conducted for sand and clay to validate suggested models and provide appropriate parameters in models [20,21]. Bearing capacities and resistances of suction buckets under a variety of operational conditions have also been studied [22,23]. Moreover, the effects of plasticity and cyclic loads on the behavior of suction buckets have been analyzed to observe operational performance and degradation of stiffness [24,25]. ese intensive studies have elaborated the model accuracy of suction buckets and provided useful information for understanding soil-pile interactions. However, many of the suggested high fidelity models require heavy computational efforts to calculate the static and dynamic responses of offshore wind conversion systems (OWCSs) comprising offshore wind turbines, substructures, and foundations. ese heavy calculation efforts hinder the use of these models in the design phases of commercial wind farms because hundreds of load calculations should be conducted in the actual design of substructures and foundations to ensure the reliability and safety of the OWCSs at a variety of operational conditions. Previous experimental studies have been conducted with laboratory model tests or reduced-scale field trials. Hence, uncertainty still exists with regard to the interapplicability of the results obtained from laboratory model tests, field trials with a reduced scale, and full-scale installations. Moreover, significant site-to-site variation of soil properties and the strain dependence of modulus of elasticity cause numerous difficulties and complications in the actual design and construction of a commercial wind farm.
us, it is essential to bridge the gap between theory and experiments as well as that between academic research and field applications and thereby mitigate concerns regarding uncertainty and increase prediction accuracy in the development phase of commercial offshore wind farms with multipod suction buckets. To this end, this study suggests a simple yet effective design procedure and its validation method. In particular, an efficient integrated framework is proposed to predict the dynamic responses of OWCSs. is framework also suggests a method of collaboration between civil engineering and mechanical engineering to design a commercial offshore wind farm.
From the civil engineering perspective, a suction bucket was designed in detail to provide sufficient bearing capacity for a selected offshore wind turbine based on an actual soil investigation, including the cone penetration test (CPT) and standard penetration test (SPT). en, design parameters of the foundation were fed into a three-dimensional (3D) finite element model (FEM) to estimate the equivalent stiffness matrix of the seabed. is matrix accounts for the coupling effect between the translational and rotational motions of the seabed to enhance prediction accuracy. From the mechanical engineering perspective, dynamic analysis was conducted via a high fidelity 3D FEM of the OWCSs with the equivalent stiffness matrix of the seabed that was obtained via civil engineering as the boundary condition.
Two full-scale measurement campaigns were also conducted during the construction and operation. e dynamic characteristics were measured to validate the prediction of the integrated simulation framework when significant progress in construction was achieved. e forces of impact caused by a ship or aero/hydrodynamic excitation forces from winds and waves were used for the excitation of structures in four stages during construction. Moreover, a structural health monitoring (SHM) system was installed to measure the dynamic characteristics during operation. e data measured using the SHM system were also used for the validation of the integrated simulation framework and evolution of dynamic responses of the OWCS. A comparison between predicted and experimental results suggests that several important features should be accounted for to accurately predict the dynamic responses of the OWCSs, including the stiffness of the cap of the suction bucket and strain dependency of the modulus of elasticity of the soil. Figure 1 shows a schematic of an integrated simulation framework for predicting the dynamic characteristics of an OWCS. Dynamic characteristics of the OWCS should be predicted to determine the specifications of the substructure in the design and load calculation phases because fundamental force frequencies (1P and 3P) should be avoided in the first natural frequency of the OWCS during operation [26,27]. Hence, prediction accuracy is very important to secure the reliability and safety of OWCSs, given that the design lifespan of an OWCS is more than 20 years. e integrated framework proposed herein is separated into two parts. e first part is estimating a 6 × 6 stiffness matrix K Seabed that represents the stiffness of the seabed including the soil and suction buckets. e method for estimating the stiffness matrix K Seabed is described in detail through experiments and simulation in Section 2.1. e second part is predicting the dynamic characteristics of the OWCS and its substructure with a simple yet effective modeling method constrained by boundary conditions from the stiffness matrix K Seabed . e modeling of the wind turbine system and substructure is described in detail in Section 2.2. Figure 2 shows the procedure for estimating the stiffness matrix K Seabed , accounting for the stiffness of the soil and suction buckets. Hence, the stiffness matrix K Seabed represents the dynamic characteristics of the seabed connected to a substructure as a unified model, as shown in Figure 1.

Soil and Suction Buckets.
In step 1, soil characteristics of the installation location of the OWCS were investigated. e candidate site is located near a breakwater in a southwestern sea area of the Korean Peninsula. e exact location is 35°58′19.43″ North (latitude) and 126°30′53.39″ East (longitude). e mean sea level of the site is 13.6 m. is site was selected as a testbed of the OWCS to establish the design procedure and criteria for the commercial development of an offshore wind farm with a multipod substructure with suction buckets. e soil investigation, including the CPT and SPT, was conducted through a site survey. Note that site surveys for conducting the CPT and SPT are indispensable in the design phase for the characterization of the actual soil properties of the site because soil properties are highly nonlinear and vary by site. is nonlinearity includes strain dependency, different dynamic behavior depending on type of soil (i.e., sand or clay), grain size, and pore water pressure. is complexity is the main reason for developing a standard procedure and code of analysis in the design phase of an offshore wind farm. Bearing resistance, friction resistance, and pore water pressure were measured through CPT, whereas soils were sampled at depth through SPT. e tests were conducted up to a depth of 10.7 m.
In step 2, soil properties of each layer were estimated through laboratory experiments with soil samples collected at depth. Table 1 shows soil properties estimated through laboratory experiments.
e results show that the site consists of three layers, with each layer featuring different soil properties. Sand is dominant in the first and third layers, whereas clay is dominant in the second layer (see Table 1).
Using an empirical equation, Young's modulus E, which is a key parameter that describes soil dynamics, can be estimated with the measured cone tip bearing resistance q c over depth. Different empirical equations should be applied to estimate Young's modulus E for sand and clay because of inherent differences between various types of soils.
To estimate the initial tangent shear modulus G max for the sand layers, 14 empirical models were used. ree models [28][29][30] directly calculated the initial tangent shear modulus G max .
ree models [29,31,32] first calculated Young's modulus E; then, this value was converted to the initial tangent shear modulus G max as . (1) Eight models [33][34][35][36][37][38][39][40] first calculated the velocity of a shear wave traveling V s through soil with the cone tip bearing resistance q c ; then, this value was converted to the initial tangent shear modulus G max as ese values were calibrated with the ratio of the secant shear modulus and initial tangent shear modulus, i.e., (G/G max ), with respect to the shear strain considering that Young's (shear) modulus of the soil significantly depends on (shear) strain. e ratio of (G/G max ) with respect to shear strain was obtained from laboratory experiments on samples from the Saemangeum Embankment, which is located near the installation site of OWCS [41] (Figure 3(a)).
To estimate Young's modulus E for the clay layer, six empirical models were used. One model [43] calculated the initial tangent shear modulus G max , whereas the remaining five models [44][45][46][47][48] calculated the velocity of a shear wave traveling through soil V s with the cone tip bearing resistance q c ; subsequently, the estimated velocity was converted into the initial tangent shear modulus G max using equation (2). All values were calibrated with the ratio of the secant shear modulus and the initial tangent shear modulus (G/G max ) with respect to the shear strain to account for strain dependency. Shear strain-shear modulus curves suggested by Vucetic and Dobry [42] were used in this study to calculate Young's modulus E of the clay layer ( Figure 3(b)). A plasticity index (PI) of 30 was used for calibrating the modulus because the PI of the clay layer was 29.7 (see Table 1).

Complexity 3
Finally, all estimated values from different empirical equations, excluding the maximum and minimum values for each layer, were averaged and converted into Young's modulus E with respect to strain ( Table 2). Table 2 shows that Young's modulus E significantly decreases with an increasing strain ε, suggesting that the strain dependency of Young's modulus E should be considered for soil that supports OWCSs because the large blades in OWCSs intensify the aerodynamic force acting on supporting structures, through which significant deformation might occur; other structures installed in the ocean are designed such that the aerodynamic force is minimized.
In step 3, 3D finite element analyses (FEAs) were conducted to elucidate soil-structure interaction between the soil and the skirt of suction buckets. is 3D FEA accounts for multilayer characteristics with different types of soil over depth. Different values of Young's modulus E were also used in case studies to account for the strain dependency of the soil. en, a flexibility matrix A soil+suc. s was calculated for an integrated dynamic analysis, where the subscript suc. s denotes the skirt of the suction bucket. Hence, A soil+suc. s is the flexibility matrix of the skirt of the suction bucket embedded in soil. Note that this matrix does not account for the flexibility of the cap of suction buckets. e flexibility of the cap is estimated in the next step.
Regarding 3D FEA, an elaborate model describing the interaction between soil and the skirt was built with hexagonal elements ( Figure 4). e diameter and length of the suction buckets are 6 m and 12 m, respectively, and the thickness of the skirt is 19 mm. e following properties of Step 1. CPT and SPT for a candidate site Step 3. estimate a flexibility matrix U of soil and a skirt of a suction bucket through 3D FEA Step 4. estimate a flexibility matrix U of a cap of a suction bucket through 3D FEA Step 5. calculate a stiffness matrix K representing soil and suction buckets Step 2. characterize soil properties over the depth: ρ, S u , v, E(ε) Figure 2: Procedure for estimating a stiffness matrix K Seabed representing soil and suction buckets.  [49], which is defined as where a ij is a flexibility influence coefficient. e flexibility of a suction bucket embedded in soil is similar to that of passive joints [50], suggesting that the flexibility matrix A includes a number of zero elements as A basic assumption in this approach is that force in the x (y) direction and moment in the y (x) direction mainly contribute to deformation in the x (y) direction. In a similar manner, force in the x (y) direction and moment in the y (x) direction mainly contribute to angular displacement in the y (x) direction. Note that a coupling effect between translational and rotational degrees of freedom (DOF) is an actual physical phenomenon in the response of the suction bucket.
erefore, the coupling effect should be considered to accurately evaluate equivalent soil stiffness. Each element of the flexibility matrix A is calculated by where U and F denote a displacement vector expressed as u x u y u z θ x θ y θ z T and a force vector expressed as f x f y f z M x M y M z T . Each element of the flexibility matrix A is calculated with the displacement vector U and force vector F when a unit force (or unit moment) is applied. Specifically, a 11 and a 51 are calculated with responses u x and θ y , respectively, by 3D FEA when a unit force (f x � 1) is applied with other zero forces and moments. In a similar manner, a 24 and a 44 are calculated with responses u y and θ x by 3D FEA when a unit moment (M x � 1) is applied with other zero forces and moments. In step 4, 3D FEAs were conducted to estimate the flexibility matrix A suc. c of the cap, where the subscript suc. c denotes the cap of the suction bucket. e diameter of the cap experiences significant changes over depth, whereas the diameter of a mono pile remains unchanged. is significant change of shape occurs in the suction bucket to provide sufficient bearing capacity with relatively small length of the skirt, whereas the mono pile generates a large bearing force with a long length of the pile. is significant change in diameter may induce stress concentration and local flexibility even if the cap is reinforced by many stiffeners. 3D FEA is an efficient solution for calculating the flexibility matrix A suc. c because of the complex design of the cap. e cap was modeled as combination of 2D plates and an edge of the cap was assumed to be clamped in the sense that the edge connecting the cap with the skirt is much stiffer than others ( Figure 5(a)). en, a unit force or a unit moment in three directions was applied to calculate the 6 × 6 flexibility matrix A suc. c with deformation of the cap (Figures 5(b) and 5(c)).
In step 5, an equivalent flexibility matrix A Seabed representing the flexibility of soil and the suction bucket is calculated as is calculation can be explained similar to the case of two spring elements connected in series. One spring represents the stiffness of the soil and the skirt, and the other represents that of the cap. Finally, the equivalent stiffness matrix K Seabed of the seabed can be calculated by matrix inversion as is stiffness matrix was used as boundary conditions for a substructure as mentioned earlier.

Wind Turbine and Substructure.
is study modeled the WinDS3000 TC-2 (Doosan Heavy Industries & Construction, South Korea) wind turbine with a tripod substructure. WinDS3000 TC-2 is a 3 MW offshore wind turbine with blades having a diameter of 100 m. Based on the specifications of WinDS3000, which are provided by Doosan Heavy Industry & Construction, the substructure was designed to prevent natural frequencies of the entire system from operational frequencies. Note that specifications of the wind turbine cannot be changed depending the site because turbine manufacturers require reissue of certificates when any specification of a wind turbine is changed [27]. Table 3 describes the dimensions of the tower and the tripod substructure in detail.
GH Bladed (DNV GL) was used to predict natural frequencies of the OWCS. A tube-shaped tower and a substructure were modeled using beam elements. An assembly of a nacelle, a hub, and blades can be modeled as a point mass. e assembly weighs 186 tons; the nacelle, hub, and blades weigh 128, 28, and 10 tons, respectively. Local flexibility of the blades, tower, and nacelle is negligible because a shape change in the tower or nacelle is small and blades of MW-class OWCS have sufficient stiffness with a sole natural frequency of above 0.6 Hz [51]. e substructures were designed as tripod suction buckets because the cost analysis for the offshore meteorological tower, Herald of the Meteorological and Oceanographic Special Research Unit (HeMOSU)-2, confirmed that the construction cost of tripod suction buckets was half of that of jacket piles, for the same seabed geology, construction, transportation, and supply chain in South Korea [26]. Note that the mean sea level for the candidate site of a 100 MW wind farm is around 20 m, meaning that gravity or monopole types are not economically feasible [52]. e pile type also gives rise to several environmental issues during construction, including vibration, noise, and suspended sediments. Moreover, the construction of HeMOSU-2 took only 6 h, whereas that of HeMOSU-1, in which jacket piles were used, took two months for the installation [53]. e properties of steel described in Section 2.1 were used for the tower and substructure. e effect of added mass was considered for the submerged substructure, which is defined as the substructure from the depth of the seabed to the depth of the mean sea level because the water added mass has a significant influence on the natural frequencies and mode shapes of a multipod substructure [54]. e inner space of the tubeshaped element below the sea level is filled by seawater and the outer surface is also in contact with seawater, suggesting that the amount of added mass is exactly related to the deformation shape of an element. However, the lowest mode deformation of an element in the finite element analysis can be assumed to be only translational. erefore, added mass of the inner and the outer water can be set by the area of the element [55].

Measurement under Construction.
e evolution of dynamic responses was measured during the construction of the 3 MW-class OWCS. e measurement framework consisted of four stages, as shown in Figure 6. e natural frequency was measured at each stage when significant progress was achieved during construction. Specifically, the natural frequency was measured after completing the installation of tripod suction buckets at stage I (Figure 6(a)) and after completing the installation of bottom and middle towers at stage II ( Figure 6(b)). In stage III, the natural frequency was measured after completing the installation of 6 Complexity the nacelle on top of the tower (Figure 6(c)). Finally, the natural frequency was measured at stage IV after completing the installation of the hub and three blades (Figure 6(d)). Two accelerometers with a sensitivity of 10 V/g were installed on the top of the transition piece to measure accelerations along two directions perpendicular to each other. Figure 7(d) shows the installed accelerometers on the transition piece.
In the experiments, the structures were excited in two ways. In one, a ship was used to provide impact force on only the substructure at stage I, whereas in the other, natural excitation from winds and waves was used to structures under construction from stage II to stage IV. Note that natural exciting forces from winds and waves were not sufficient to excite only the substructure in stage I. Hence, a small ship of 6.6 ton was used to provide impact force to the substructure (Figure 7(c)) because impact force in the time domain is white noise in the frequency domain [56]. e ship slowly approached on landing ladder to excite a substructure (Figures 7(a) and 7(b)); then, the dynamic response was measured. In contrast, structures at stages II, III, and IV were excited sufficiently by natural wind and waves. Hence, the natural frequencies at stages II, III, and IV were determined by measuring acceleration and through subsequent spectral analysis.

Measurement during Operation.
A SHM system was installed to observe long-term variations of the dynamic responses of the OWCS during operation, such as acceleration, strain, temperature, and sea level. e natural frequency of the OWCS with the substructure was calculated using measured accelerations and strains. To measure the accelerations, 10 accelerometers were installed at 5 locations on the tower of the OWCS along the vertical direction (Figure 8(a)). Two accelerometers were installed perpendicular to each other at different heights in each location. For each set, one accelerometer was aligned north-northwest (NNW), the main direction of wind at the site, whereas the  other accelerometer was aligned east-northeast (ENE), perpendicular to the main direction of wind (Figure 8(b)) [57]. In addition, 4 sets of strain gauges were installed near the transition piece between the offshore wind turbine system and the substructure to measure bending deformation. Two sets of strain gauges were installed +0.5 m from the transition piece; one set was aligned NNW and the other was aligned ENE. Two other sets of strain gauges were installed at −0.5 m from the transition piece. Each set of strain gauges is a full bridge type.

Stage I.
e dynamic response of the tripod suction bucket measured at stage I was compared to that predicted from the integrated framework. To calculate the natural frequency based on an experiment performed at stage I, the measured acceleration (Figure 9(a)) with a record length of 11.8 s and sampling frequency of 256 Hz was converted to the frequency domain via fast Fourier transform (FFT). en, the maximum peak in the spectrum was determined (Figure 9(b)). Hence, the resolution in the frequency domain for this data is 0.085 Hz. Figure 9(a) shows the measured acceleration in the time domain, which follows the conventional transient response of a second-order system. is observation can be explained by the fact that the first mode is dominant for the tripod suction buckets embedded in the seabed; thus, it governs the overall dynamic behavior of the OWCS. e natural frequency of the tripod substructure determined through spectral analysis is 3.30 Hz, as shown in Figure 9(b).
Four different cases were studied in the integrated framework to compare not only the effect of the cap of the suction bucket but also the strain dependency of the modulus (Table 4). In each case, the flexibility matrix A soil+suc. s accounts for only the stiffness of the soil and the suction skirt, as in many previous studies, whereas the flexibility matrix A Seabed accounts for all stiffness components, including the soil, suction skirt, and suction cap. In case 1, strain is assumed to be less than 10 −3 for all layers, whereas in case 2, it is assumed to be approximately 0.01 for all layers. In case 3, the two uppermost layers are considered to experience a large strain of approximately 0.1, whereas the deepest layer experiences a deformation of approximately 0.01, assuming that a smaller deformation occurs in the deeper layer because of relatively high stiffness. In case 4, a large deformation of approximately 0.1 is considered to occur for all layers.
e last two columns of Table 4 show the predicted natural frequencies (and the corresponding errors compared to the measured values that are given in parentheses) calculated with the corresponding flexibility matrix. Remarkably, relatively large errors were encountered when only the flexibility matrix A soil+suc. s was considered. e errors ranged from 8.18% to 30.3% when only the flexibility matrix A soil+suc. s was considered, whereas they ranged from 1.82% to 14.2% when the flexibility matrix A Seabed was considered. e errors attributed to this difference were at least double, i.e., 30.3% vs. 14.2% (case 1) or 8.18% vs. 1.82% (case 4), when the same modulus for each layer was applied. e analysis suggests that the stiffness of the cap plays a critical role in the prediction of dynamic responses of the substructure. erefore, the contribution of the cap should be considered.
From the perspective of strain dependency, the error could be reduced by considering larger deformations in soil. Specifically, the error between the measured and predicted values was 1.82% when the strain for all layers is assumed to be approximately 0.1 (case 4 with the flexibility matrix A Seabed ), whereas the error between measurement and prediction was 14.2% when strain is assumed to be less than 10 −3 (case 1 with the flexibility matrix A Seabed ). e analysis of the dynamic characteristics of the tripod substructure in stage I clearly shows that the stiffness of the cap and Young's modulus, which accounts for strain dependency [58], should be considered to accurately predict the dynamic response of the OWCS and its substructure.

Stages II to IV.
e natural frequencies under construction during stages II to IV were also obtained from the measured acceleration under natural excitation from wind and waves. In stages II to IV, acceleration could be measured over a long period because these experiments were conducted under natural excitation, whereas stage I was conducted using an impact test with a ship. erefore, analysis could be conducted with high resolution in the frequency domain. e accelerations excited by wind and wave were measured for a period of 125 s with a sampling frequency of 1707 Hz, meaning that the frequency resolution at the analysis is 0.008 Hz. Figure 10 shows results of spectral analyses for stages II to IV. ese frequency spectra show one clear peak even though the structure is under natural excitation, suggesting that natural excitation from wind and waves was sufficient for determining the natural frequencies of the structure under construction. Table 5 shows a comparison between the measured and predicted values of natural frequency. e measured values were compared with the predicted values for cases that showed the largest and smallest errors in stage I, namely, case 1 and case 4, respectively. Similar to the results obtained in stage I, case 4 showed smaller errors than case 1. Moreover, the percentage errors of case 1 with the flexibility matrix A soil+suc. s in stage II were much smaller than those in stage I. is observation suggests that the contribution of the equivalent stiffness of the seabed to the natural frequency of the OWCS during stage II was lesser than that at stage I. In a similar manner, the errors of case 4 with the flexibility matrix A soil+suc. s in stages II, III, and IV decreased compared to those in stage I. e predicted results for case 4 also show  that the flexibility matrix A Seabed shows more accurate predictions compared to the flexibility matrix A soil+suc. s . In conclusion, the model that accounts for the stiffness of the cap and strain dependency demonstrates increased prediction accuracy. e suggested elaborate model and method of predicting the dynamics of a substructure and foundation would be useful in predicting the long-term degradation of soil supporting the OWCS because soil degrades faster than the suction bucket and the substructure [59]. Figure 11 shows the results of a spectral analysis of data obtained from a SHM system under idling (stop) and operational conditions. In the idling or stop conditions, the first natural frequency calculated was 0.318 Hz for the spectra of both acceleration and strain measured at various locations, further validating that the first natural frequency measured at stage IV is accurate. Moreover, the spectra did not exhibit location dependency between NNW and ENE. e spectra for ENE are omitted  e natural frequency does not appear for the measured acceleration (Figure 11 11(c)). e spectra of acceleration simply show various peaks related to 1P (rotational speed), 3P (blade passing component), and their harmonics. In contrast, the natural frequency clearly appears along with operational forcing frequencies for the measured strain ( Figure 11(d)). e same phenomena were observed in the spectra measured using the accelerometers and strain gauges at ENE. is observation can be explained by the fact that various auxiliary machines, including oil pumps, cooling fans, and electric equipment, are switched on during operation and these machines increase the level of background vibration noise. is could be the reason for previous studies measuring the natural frequency with accelerometers after applying emergency stops [60,61]. In other words, a significant impact force would be induced in measuring the dynamic responses with the accelerometers [62]. In contrast, the spectra of strain show the natural frequency of the OWCS as well as frequencies associated with operation. Hence, it can be inferred that the measurement of strain using SHM is more effective than that of acceleration in characterizing the dynamic response or monitoring the evolution of natural frequencies of the OWCS because the frequencies of interest in this study are low frequencies. Note   that acceleration is proportional to the second derivative of the displacement with respect to time. is means that the acceleration is proportional to the square of the frequency in the frequency domain. Hence, the strain would be more apparent than the acceleration at low frequencies, but the acceleration would be larger than the strain at high frequencies.

Dynamic Responses during Operation.
In addition, the natural frequency under the rated operating condition (0.33 Hz) was slightly higher than that under the idling condition (0.32 Hz). is result is similar to that in [63], in which the first natural frequency predicted under the idling or stop conditions was slightly lesser than the measured frequency.
Several analytical approaches have been applied to elucidate the load dependency of the dynamic response, suggesting that the load dependency is primarily relevant to soil-pile interactions. One reason is that the static thrust force induced on the OWCS is significantly larger than the dynamic and cyclic force, as shown in Figure 12. is static load changes the friction angle of soil because of soil-pile interaction [64]. Another reason is that the equivalent stiffness of soil depends on the ratio of the diameter and length of embedding [60]. A relatively large ratio increases the bending stiffness of soil and piles, thereby increasing the first natural frequency with increasing load induced on the OWCS. is phenomenon can explain the slight increase in the natural frequency measured via SHM because the key characteristic of the suction bucket is the relatively large ratio of the diameter with respect to the length of embedding. Sediment transportation at the seabed also results in changes in the dynamic characteristics of the OWCS [61]. e scour depth depends on time because of aero-and hydrodynamic loads and thereby changes the eigenfrequency and modal damping. Compaction of the backfilled material occurs in the presence of waves, thereby inducing depth compaction. is phenomenon increases the friction angle and thereby shifts the eigenfrequency of the OWCS [64].
is interesting phenomenon might also be related to the pull in and out mechanism of the suction bucket. Under operational conditions, a significant thrust force, which is larger than the cyclic force, as shown in Figure 12, acts against the tip bearing force, which is a major supporting force of the suction bucket. erefore, the equivalent stiffness of soil, including that of the suction bucket, increases. Note that the suction bucket provides two supporting forces for the structure. One is skin force between the soil and the skirt. e skin force is related to contact friction, which is influenced by pressure and type of soil. e contact friction is independent of the friction angle of the soil, which is defined by the inherent characteristics of the soil. e other is tip bearing force, which is the supporting force from the soil and is transmitted to the perimeter of the suction bucket. e contribution of the tip bearing force would be significant in the sense that the stiffness of the cap plays a critical role, as already shown in case studies. Moreover, the tip bearing force is governed by the friction angle, which is the same as shear strength.

Conclusion
is study evaluates the dynamic characteristics of an OWCS through simulation and experiments. An integrated framework was suggested to accurately model the seabed, including the soil, the skirt of a suction bucket, and the cap of a suction bucket. Full-scale experiments were also conducted for various stages under construction as well as under operation. e following conclusions can be drawn: (i) e stiffness of the cap of a suction bucket plays a critical role in predicting the natural frequencies of the OWCS. In particular, the contribution of the stiffness of the cap of the suction bucket is more significant when predicting dynamic responses of the substructure embedded in the seabed with suction buckets [65]. (ii) e strain dependency of the modulus of elasticity of the soil significantly affects the dynamic responses of the substructure. However, the effect of the strain dependency is limited with regard to predicting the dynamic responses of the entire structure. Hence, it can be inferred that the suction bucket embedded in soil should be designed with a focus on providing sufficient bearing capacity, whereas the substructure submerged in seawater should account for the dynamic responses of the OWCS and its substructure. Future work will include a systematic study using a 3D finite element model of the soil and the suction bucket to elucidate the strain with respect to depth and soil-pile interactions. (iii) e natural frequency during the operational stage slightly differs from that during stop or idling conditions. is may be attributed to soil-pile interactions in the seabed. Specifically, the tip bearing force governed by thrust forces contributes to the evolution of the natural frequency depending on operational conditions such as stop and rated operation. (iv) During operation, strain is more sensitive than acceleration for measuring the dynamic characteristics or long-term evolution of dynamic responses of the OWCS and its substructure. Hence, strain measurement can be used to predict the remaining useful life of the OWCS with the SHM system. Shear wave velocity f * : Force in * direction q c : Cone tip bearing resistance u * : Displacement in * direction θ * : Angular displacement in * direction ϕ: Internal friction angle ε: Strain ]: Poisson's ratio c d : Unit weights for zero water content.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest.