Combined Seismic and Scoured Numerical Model for Bucket-Supported Oﬀshore Wind Turbines

: Numerous oﬀshore wind turbines (OWTs) with bucket foundations have been installed in seismic regions. Compared to the relative development of monopiles (widely installed), seismic design guidelines for bucket - supported OWTs still need to be developed. Moreover, scour around bucket foundations induced by water –current actions also creates more challenges for the seismic design of OWTs. In this study, a simpliﬁed seismic analysis method is proposed that incorporates the soil–structure interaction (SSI) for the preliminary design of scoured bucket - supported OWTs, aiming to balance accuracy and eﬃciency. The dynamic SSI eﬀects are represented using lumped parameter models (LPMs), which are developed by ﬁtting impedance functions of the soil – bucket foundation obtained from the four - spring Winkler model. The water –structure interaction is also considered by the added mass in seismic analysis. Based on the OpenSees 3.3.0 platform, an integral model is established and validated using the three-dimensional ﬁnite element method. The results indicate that the bucket - supported OWT demonstrates greater dynamic impedance and ﬁrst -order natural frequency compared to the monopile-supported OWT, which has an increased seismic response. Seismic spectral characteristics and intensities also play an important role in the responses. Additionally, scour can change the bucket impedance functions and the frequency characteris tics of the OWT system, leading to a signiﬁcant alteration in the seismic response. Scour eﬀects may be advantageous or disadvantageous, depending on the spectral characteristics of seismic excitations. These ﬁndings provide insights into the seismic response of bucket - supported OWTs under scoured conditions.


Introduction
Recently, the offshore wind turbine industry has experienced significant growth, with the global installed capacity reaching 75.2 GW by the end of 2023.Notably, China contributed an impressive 50.5% share of this capacity [1].Additionally, bucket foundations have demonstrated a potential developmental trend in terms of cost-effectiveness and reliability compared to their counterparts [2,3].Compared to those of well-established monopile foundations [4], the static characteristics [5][6][7] and wind-wave cyclic loading performance [3,8] of bucket foundations have been extensively researched, yielding valuable insights.However, further exploration is needed regarding their seismic performance.Furthermore, many offshore wind farms in China are located in seismic regions [9,10].The seismic design of offshore wind turbine systems is critically important.
Numerous scholars have investigated the seismic responses of bucket-supported offshore wind turbines using experimental and numerical methods.Experimental methods typically include shaking table tests [11][12][13], field tests [14], and centrifuge model tests [15][16][17].However, experiments can reflect the nonlinearity of soil and the actual dynamic soil-foundation interaction.Specifically, real seismic loads and stress conditions can be replicated in centrifuge shaking table tests, which can accurately reflect structural dynamic characteristics and seismic responses [17].However, conducting experiments is usually uneconomical and time-consuming; therefore, numerical methods have been employed by many scholars to investigate the seismic analysis of bucket-supported OWTs.The three-dimensional finite element modeling method has been generally adopted [18][19][20][21][22][23][24].The method can reflect well the nonlinear dynamic behavior, incorporating soil material nonlinearity, soil-foundation gaps, sliding, and other geometric nonlinearities, through the selection of a suitable soil constitutive model and soil-foundation interface contacts.However, detailed simulations inevitably result in a decrease in computational efficiency and fail to meet the rapid and efficient design needs of preliminary offshore wind turbine system design.Therefore, simplified dynamic analysis methods need to be proposed to improve computational efficiency.Alati et al. [25] directly treated the bottom constraint of the turbine as a fixed base without considering the soil and foundation, which is commonly adopted by structural engineers.However, ignoring the SSI effect greatly overestimates the natural frequency of the turbine [26], leading to significant errors in the seismic design of OWTs.Therefore, soil-structure interactions must be considered in simplified dynamic analysis methods to balance accuracy and computational efficiency.Andersen et al. [27] obtained the impedance of bucket foundations by utilizing a three-dimensional-coupled boundary-element/finite-element scheme, including torsional, horizontal, and rocking directions.Then, a computationally efficient lumped parameter model (LPM) of the bucket foundation was developed and utilized within aeroelastic codes to rapidly evaluate the dynamic responses of OWTs [28].Suryasentana et al. [29] proposed a 1D elastoplastic Winkler model for suction caisson foundations to capture the nonlinear behavior of soil-foundation interactions.Antoniou et al. [30] also presented an enhanced Winkler-based model [31], where nonlinear hysteretic elements replaced the soil.The model was successfully used to study the performance of bucket-supported OWTs under combined seismic and wind loading.Skau et al. [32,33] adopted a macroelement approach to reflect the nonlinear behavior of soil and bucket foundations under cyclic loading.Based on substructures, Li et al. [34] proposed a finite element method to study the equivalent stiffness of bucket foundations for OWTs.
However, there are limited simplified methods available for analyzing seismic issues concerning bucket-supported offshore wind turbines.Lumped parameter models have been extensively utilized in the seismic analysis of offshore structures such as bridges, platforms, and wind turbines.Furthermore, based on linear elastic assumptions, the model offers greater computational efficiency than does the nonlinear Winkler beam model, making it suitable for preliminary seismic design requirements of OWT systems.Dezi et al. [35] employed an LPM to investigate the seismic response of a three-span bridge, which was constructed by fitting the pile impedance functions.Liang et al. [36] expanded this study and modified the impedance functions to consider the stress history effect caused by scour.Damgaard et al. [37,38] proposed a simplified analysis method for monopile-supported OWTs using a lumped parameter model under normal operating conditions while also exploring the impact of soil uncertainties on the dynamic characteristics of OWTs [39,40].Taddei et al. [41] developed a numerical model of 5 MW offshore wind turbines for seismic analysis utilizing a simplified lumped parameter model to simulate SSI.Padrón et al. [42] investigated the seismic responses of large OWTs under aerohydrodynamic loads by incorporating a lumped parameter model to account for the SSI.Similarly, Sánchez et al. [43] employed a similar model and outlined specific steps for input motion and SSI considerations.Liang et al. [44] examined the influence of marine environmental factors on the seismic responses of monopile-supported OWTs based on a lumped parameter model.
Nevertheless, the lumped parameter model is primarily employed for the seismic responses of monopile-supported OWTs.There is limited research on developing dynamic LPMs for bucket foundations.Additionally, scour phenomena are common in complex marine environments, where soil erosion around foundations is induced by water flow and wave actions [45].It has been demonstrated that soil-foundation dynamic characteristics can be impacted by scour [36,45], leading to changes in the seismic responses of OWTs [46,47].Neo et al. [48] studied the post-scour seismic responses of offshore wind turbines in sandy soils, revealing significant changes in seismic behavior due to scour.Li et al. [24] investigated the combined effects of scour and earthquakes using the ABAQUS program, revealing a notable increase in offshore wind turbine seismic responses.As a consequence, the seismic responses of bucket-supported OWTs under scoured conditions considering soil-structure interactions based on LPMs are still unclear.
This study focuses on the numerical investigation of seismic analyses of bucket-supported OWTs.To facilitate rapid and efficient preliminary seismic design, a simplified dynamic analysis method is proposed.A lumped parameter model for bucket foundations is employed to ensure a certain level of accuracy while prioritizing computational efficiency.Section 2 outlines the methodology used in this study, including the establishment of the lumped parameter model and the integral model of the wind turbine system.The earthquakes selected in this study are described in Section 3. Section 4 presents the seismic responses of the OWT under different seismic intensities, spectrum characteristics, and scouring conditions.Additionally, to investigate the seismic performance of bucket-supported OWTs, a comparison with well-established monopile foundations [4,46] is also performed.Finally, the primary conclusions are summarized in Section 5.

Materials and Methods
Wolf [49] emphasized the importance of soil-structure interactions in seismic response analysis, particularly in soft clays, which commonly exist in the southeastern seismic regions of China [50].The numerical methods for determining soil-structure interactions generally include direct and substructure methods.The former treats the soil, foundation, and superstructure as a whole for analysis, but a large number of degrees of freedom impose a significant computational burden.The latter divides them into two substructures separated by the mudline, which can significantly improve efficiency [49].For the preliminary seismic design of OWTs, this study seeks efficiency while ensuring accuracy, and the substructure method based on lumped parameter models is employed.Specifically, to facilitate time-domain analysis of the superstructure, a dynamic LPM is used to consider the SSI effects.This model consists of a series of frequency-independent springs, dampers, and mass blocks, which can be obtained from fitting the impedance functions of the soil-bucket system.Subsequently, based on the OpenSees 3.3.0finite element platform, an integral model of a bucket-supported offshore wind turbine under scoured conditions is established.It is worth noting that OpenSees 3.3.0 is object-oriented software that is commonly used to develop finite element models to simulate the responses of structural systems subjected to earthquakes.In consideration of practical applications, a simplified seismic analysis method is proposed for bucket-supported offshore wind turbines.

Numerical Model
In this paper, a 6.45 MW reference bucket-supported offshore wind turbine located in the southeastern coastal region of China is employed for this case study, as shown in Figure 1.Table 1 presents the key dimensions and physical parameters of the offshore wind turbine.   1, with a slenderness ratio (depth-to-diameter ratio) of 0.33 (<0.5), the bucket foundation for the reference wind turbine in this study is categorized as a wide-shallow bucket foundation [51].The primary deformation mode of the foundation is rigid rotation, and excessive horizontal displacement can lead to soil punching failure [52].He et al. [53] indicated that when the slenderness ratio of a bucket foundation is less than 1, the foundation can be considered to be rigid.This assumption has been widely adopted by numerous scholars to investigate the static characteristics [54] and dynamic properties [27,[55][56][57] of bucket foundations.In this study, the bucket foundation is also treated as being rigid, and the dynamic impedance functions are derived based on the four-spring Winkler model, as proposed by Gerolymos and Gazetas [58] for rigid caisson foundations.As illustrated in Figure 2, four types of complex springs are employed to simulate the interactions between the soil and foundation, which are combined with dampers: concentrated translational and rotational complex springs (  h K ,  r K ) at the bottom of the foundation and distributed translational and rotational complex springs (  x k ,  kθ ) along the sidewall of the foundation.
(1)  The definition of the dynamic impedance function is the ratio between the harmonic excitation force and the displacement at the top foundation.Therefore, the dynamic impedance matrix, referred to as the base of the bucket, can be represented by Equation (4).According to the coordinate transformation, the impedance of the bucket foundation can be expressed as Equation (5): where C is the complex stiffness matrix of the foundation. HH S and  MM S represent the horizontal and rocking impedance functions, respectively, while  HM S and  MH S are the coupled horizontal-rocking impedance functions.
According to the four-spring Winkler model depicted in Figure 1b, the complex stiffness matrix of the foundation  b K , referred to as the base of the bucket, can be derived through the integration and accumulation of the reaction forces from multiple springs: where represent the concentrated translational and rotational complex stiffness matrices, respectively, while  x k and  kθ represent the distributed translational and rotational complex stiffness matrices, respectively.l represents the embedded length of the bucket foundation.If the bucket is completely embedded, when considering the scour depth ( .The impedance functions of foundations on the half-space surface [59] are used to represent the concentrated translational and rotational complex stiffness matrices of the bucket [60,61].
where h K and h k represent the horizontal static and dynamic stiffness coefficients of foundations on the half-space surface, respectively; h C denotes the horizontal damping coefficient, which can simultaneously reflect the radiation damping and material damping of soil; D represents the diameter of the bucket foundation; G , υ , ρ , and s V repre- sent the shear modulus, Poisson's ratio, density, and shear wave velocity of soil, respectively, as shown in Table 2; and 0.03 represents the hysteretic damping (material damping), which is frequency-independent and typically ranges from 0.03 to 0.06 [62].
where r K and r k denote the rocking static and dynamic stiffness coefficients of founda- tions on the half-space surface, respectively; r C represents the rocking damping coeffi- cient, which can simultaneously reflect soil radiation damping and material damping; 0 a is the dimensionless frequency; and r c represents the radiation damping coefficient, which can be obtained from Gazetas [59].When treating the bucket foundation as an embedded rigid foundation, the dynamic impedance  emb K can also be determined as follows [59,60,62,63]: where  HH K and  MM K represent the horizontal and rocking impedances of the embedded foundations, respectively, and  HM K and  MH K are the coupled horizontal-rocking im- pedances of the embedded foundations. where

/ 1-
where MM K and MM C are the rocking dynamic stiffness and damping coefficient of the embedded foundations, respectively; w Γ is the rocking embedment factor [58,61]; and , emb r χ is the rocking dynamic stiffness coefficient of the embedded foundations.
To ensure the equality of the main diagonal elements and similarity of the off-diagonal elements of matrices  b K (Equation ( 6)) and  emb K (Equation ( 15)), the complex spring stiffness matrix along the sidewall of the bucket foundation is derived [58]: Based on Equations ( 3)-( 6), the dynamic impedance functions of the bucket foundation  S can be obtained, including the horizontal (  HH S ), rocking (  MM S ), and coupled hor- izontal-rocking impedance components (  HM S ,  MH S ).

Lumped Parameter Model Considering the SSI
To facilitate nonlinear time-domain analysis of the offshore wind turbine, a frequency-independent lumped parameter model of the bucket foundation is employed to account for the SSI effect.This model, consisting of several linear springs, dampers, and mass blocks [35,64], is obtained by fitting impedance functions of the soil-foundation system using the least squares method.Additionally, the parameters in this model are frequency-independent.Therefore, direct substitution of the soil-foundation in the time-domain analysis of the superstructure is allowed, resulting in a significant increase in computational efficiency.Various lumped parameter models have been proposed by many scholars [35,[37][38][39][40]42,64].The LPM introduced by Dezi et al. [35] features a reduced number of undetermined parameters (only nine, as depicted in Figure 3) and a simpler model architecture, which is chosen in this study for its simplicity and efficiency over alternative models, with the impedance function expressed by Equations ( 28)- (30).
Figure 3.A lumped parameter model for the bucket foundation.

Integral Model of the Wind Turbine System
Based on the OpenSees platform, a comprehensive numerical model of the bucketsupported OWT is constructed, as illustrated in Figure 1b.Due to its object-oriented nature, OpenSees requires the writing of a Tcl code file with commands to define the model, materials, boundary conditions, and seismic loading.Subsequently, the script is run in an OpenSees interpreter for seismic analysis of the OWT.The modeling environment is threedimensional with six degrees of freedom.A lumped parameter model is employed at the mudline to constrain the superstructure above the mudline, including the transition segment, tower, and rotor-nacelle assembly (RNA).To incorporate varying cross-sectional properties, the transition segment and tower are divided into 25 and 34 segments, respectively, each of which is equipped with distinct properties.These segments are established using displacement-based beam-column elements.Subsequently, the nacelle and blades are simplified as concentrated masses.Considering the geometric dimensions of the nacelle, a rigid link element is utilized to simulate the distance from the center of the mass to the top of the tower.
Additionally, water-structure interactions impact the dynamic responses of offshore wind turbines and can be represented in seismic analysis using the water-equivalent additional mass.The transition segment between the bucket foundation and the tower can be simplified as a frustum-shaped foundation, as indicated in Figure 4a.To calculate the water-equivalent additional mass, the simplified formulas proposed by Wang et al. [65], which were derived using finite element techniques, are employed.A dimensionless parameter b is introduced ( .In this study, the width-todepth ratio of the simplified transition segment is 0.643 with a slope angle of 63°, which falls within the effective range.Then, based on Equations (31)(32)(33)(34)(35)(36)(37)(38)(39)(40), the distribution curve of the additional mass density along the foundation height is calculated and illustrated in Figure 4b.d q e q e = + (36) 3. 27 6.998 0.0531

Simulation of Scoured Conditions
Scour phenomena refer to the erosion or removal of sediment from the bed induced by wave-current actions [66][67][68].They naturally occur around structures such as offshore wind turbine foundations, bridges, and piers, compromising their stability and complicating seismic analysis [46].Scour depth is widely acknowledged as the most significant factor under scoured conditions [69], with limited impacts observed when treating local scour as general scour [70].Onsite monitoring data (see Figure 5) indicate that a southeastoriented scour hole formed within a 100 m radius of the foundation center, and the scour phenomena were more pronounced within a radius of 36 m.The maximum scour depth relative to the mudline was approximately 6.05 m (0.5 L).In other studies, scholars have typically defined the range of scour depths of bucket foundations as 0.3-0.5 L [24,48,71].Therefore, this study concentrates on how scour depth affects the seismic responses of OWTs.For demonstration purposes, scour depths are set to 0, 0.1 L, 0.3 L, and 0.5 L.

Validation
The feasibility of the present study is validated by developing the same model of the 6.45 MW OWT based on the three-dimensional finite element method.This method is widely recognized as relatively accurate and has been adopted by numerous scholars for seismic analyses of bucket-supported OWTs [18][19][20][21][22][23][24].The Mohr-Coulomb model has been used to simulate soil, with the wind tower and bucket foundation modeled using linear isotropic material with C3D8R elements.Figure 6 shows the impedance functions of the bucket foundation and the corresponding LPM fittings, which are used to establish the proposed integral model in Section 2. Figure 7 displays first-order natural frequencies and seismic response envelopes of the bucket-supported wind turbine obtained from the present study and the three-dimensional finite element model (FEM) under Kobe excitation (0.15 g).Natural frequencies represent structural vibration modes, and first-order natural frequencies significantly influence system behavior [72,73].The difference in the natural frequencies between the models is 1.3%, indicating that both models have similar dynamic characteristics.Additionally, the envelope trends along the offshore wind turbine (OWT) show consistent patterns, with the maximum absolute acceleration observed at the middle tower and the maximum absolute displacement at the RNA differing by less than 1% and 2%, respectively.Therefore, the feasibility of the present study is further demonstrated.In addition, the proposed simplified dynamic analysis approach is based on linear elastic theory, while the three-dimensional FEM relies on ideal elastic-plastic theory and can consider material and geometric nonlinearity.Therefore, a higher first-order natural frequency is obtained from the proposed method, which overestimates the soil stiffness, as shown in Figure 7a.As a consequence, there are differences in the seismic responses between the two methods (Figure 7b,c).However, refined modeling (three-dimensional FEM) always demands substantial computational resources, which does not align with the rapid and efficient requirements of preliminary wind turbine design.Nevertheless, it is undeniable that the present study can ensure a certain level of accuracy while significantly improving efficiency.Furthermore, the proposed method for monopile-supported OWTs has also been validated and successfully applied in seismic analyses [44].

Seismic Excitations
Three earthquakes with different spectral characteristics were chosen as excitations in this study: Kocaeli motion, Kobe motion, and Acc100 motion.Specifically, the Kocaeli and Kobe motions are seismic records derived from significant earthquakes that occurred in Kocaeli, Turkey, and Kobe, Japan, respectively.These motions are based on ground motion recordings from actual seismic events that occurred in these regions.Conversely, Acc100 motion is a synthetic motion derived from the geological conditions of southeastern coastal regions in China [44]. Figure 8 presents the Fourier spectra and acceleration time histories of the chosen earthquakes.The frequencies of the Kocaeli motion range from approximately 0.14 Hz to 0.36 Hz, while those of the Kobe motion range from 0.586 Hz to 2.94-3.56Hz.The main frequency range of the Acc100 artificial motion is from 1 to 3 Hz.The amplitudes of each input earthquake are scaled to 0.05 g, 0.15 g, and 0.30 g to study the impact of the seismic intensity on the responses.The accelerations correspond to seismic fortification intensities of 6, 7, and 8 degrees, essentially covering the seismic fortification levels of offshore wind farm areas along the southeastern coast [74].

Effect of Seismic Intensities
Figure 9 illustrates the time history responses under Kobe excitation with different seismic intensities, including 0.05 g, 0.15 g, and 0.30 g.The seismic responses typically tend toward zero as time progresses, which is generally attributed to energy dissipation by damping [75].In comparison with the Kobe input motion shown in Figure 8, the filtering effect of the structure on the moderate-to high-frequency components of free-field motions is observed.Moreover, as the seismic intensity increases from 0.05 g to 0.15 g and 0.30 g, respectively, the maximum acceleration, displacement at the RNA, bending moment at the mudline increase by approximately 67% and 83%, respectively.Figure 10 displays response envelopes for accelerations (Figure 10a), displacements (Figure 10b), and bending moments (Figure 10c) along the OWT.Despite the seismic intensity variations, the response envelope trends remain consistent, with the maximum acceleration, displacement, and bending moment occurring at specific locations: the middle of the tower, the RNA, and the mudline, respectively.With increasing seismic amplitude from 0.05 g to 0.15 g and 0.30 g, the maximum acceleration, displacement, and bending moment of the OWT all increase by approximately 67% and 83%, respectively.The responses exhibit nearly linear increases, suggesting that the wind turbine system predominantly maintains an elastic state within the selected seismic intensity range.

Effect of Seismic Spectral Characteristics
Figure 11 presents the response envelopes under excitations with identical intensities (0.15 g) but different spectral characteristics (Kocaeli, Kobe, and Acc100).The trends of the envelopes remain consistent, with the maximum absolute acceleration (Figure 11a), displacement (Figure 11b), and bending moment (Figure 11c) occurring at specific locations: the middle of the tower, the RNA, and the mudline, respectively.Notably, the seismic responses under Kocaeli motion excitation are more pronounced due to its greater proportion of low-frequency components.Figure 5b indicates that the prominent frequencies of the Kocaeli motion range from approximately 0.14 Hz to 0.36 Hz, closer to the firstorder natural frequency of the OWT (0.305 Hz).Conversely, the prominent frequencies of the Kobe motion are relatively distant from the natural frequency of the OWT, resulting in a reduced response under the same seismic amplitude.Therefore, the seismic responses of OWTs are affected not only by the earthquake intensity but also by the seismic spectral characteristics.This finding is consistent with those of other scholars [24,44,46,48,66].Figure 12 shows the normalized maximum responses corresponding to the accelerations (Figure 12a), displacements (Figure 12b), and bending moments (Figure 12c) under the different earthquakes.Similar to the observations in Figures 10 and 11, the seismic response is amplified when the motion excitation has a higher amplitude or its prominent frequency aligns more closely with the natural frequency of OWTs.Moreover, it is worth noting that the responses exhibit nearly linear increases with increasing seismic intensity, suggesting that the wind turbine system predominantly maintains an elastic state under all selected earthquake excitations.The findings also suggest that the proposed simplified method in this paper seismic analysis of wind turbines when intensity less than g.

Effect of Foundation Types
In the initial development of offshore wind turbines, monopile foundations emerged as the predominant choice due to their simplicity and cost-effectiveness.To date, more than 70-80% of OWTs are still supported by monopiles [4,46].Seismic analyses based on lumped parameter models have been extensively conducted to study the performance of monopile-supported OWTs [37][38][39][40][41][42][43][44].However, investigations into the dynamic characteristics and seismic performance of bucket-supported OWTs still need to be conducted.In this section, numerical models for both bucket-supported and monopile-supported OWTs are constructed using proposed simplified methods based on LPMs. Figure 13 shows a schematic of both physical models.Through a comparison of foundation impedance functions, turbine system characteristics, and seismic response features, the performance of bucket-supported turbines relative to monopile-supported turbines is evaluated.To ensure consistency in parameter analysis, only the foundation type while other conditions remain unchanged, including the wind farm location (with similar geological conditions), tower height, and turbine size (with the same capacity of 6.45 MW).The dimensions and lumped parameter models of the monopile foundation are derived from Liang et al. [44].

Dynamic Characteristics
Figure 14 presents the impedance functions and LPM approximations of the bucket and monopile foundations.The positive and negative signs on the vertical axis in the figure represent directionality only.Significant differences are evident in the dynamic impedance functions of the two types of foundations, which are attributed to their distinct load characteristics and dynamic properties.Overall, the wide and shallow bucket foundations demonstrate greater dynamic impedances than monopile impedances, including both dynamic stiffness (real part of the impedance: Re) and damping (imaginary part of the impedance: Im).This indicates substantial disparities in the dynamic characteristics between the two foundation types.Figure 15 shows a comparison of the mode frequencies between the bucket-supported OWT and the monopile-supported OWT, including 16 vibration modes.The firstorder of the bucket-supported offshore wind turbine system is 0.305 Hz, while that of the monopile-supported turbine system is 0.279 Hz.The dynamic behavior of the bucket-supported OWT is relatively more rigid than that of the monopilesupported OWT under the first mode.This phenomenon may be attributed to the shallow and wide characteristics of the bucket foundation, while the monopiles are relatively slenderer in this study.

Seismic Responses
Figure 16 displays the time history responses under the same Kobe motion (0.15 g) excitation with different foundation types.The energy dissipation phenomenon of the bucket-supported OWT is more pronounced, primarily due to the distinct dynamic characteristics exhibited between the bucket and monopile foundations.As illustrated in Figure 14, the bucket foundation exhibits greater dynamic impedances, comprising dynamic stiffness (Impedance-Re) and damping (Impedance-Im).Additionally, the monopile-supported OWT experiences a percentage increase of 46.5% in the maximum acceleration and a percentage increase of 29.0% in the maximum displacement at the RNA, while the maximum bending moment at the mudline is 9.6% greater than that at the bucket-supported OWT.Due to the different foundation types, comparing foundation responses is meaningless, and only tower responses are compared, as shown in Figure 17.The trends tower response envelopes remain consistent across different foundation types, with the maximum acceleration (Figure 17a), displacement (Figure 17b), and bending moment (Figure 17c) occurring at specific locations: the middle of the tower, the RNA, and the top of the foundation, respectively.Compared to the bucket-supported OWT, the monopilesupported OWT exhibits 43.9% greater maximum tower acceleration, 30.3% greater maximum tower displacement, and 19.4% greater maximum tower bending moment.Figure 18 illustrates the impedance functions and LPM approximations of the bucket foundation before and after scour (Sd = 0.5 L).There are changes in the dynamic impedances of the soil-foundation system before and after scour.The scour effects erode the soil around the foundation, resulting in a decrease in soil stiffness, which is the primary cause of the alterations in the dynamic characteristics of the OWT.Therefore, it is essential to account for the scour effect in the dynamic analysis of the soil-foundation system.Figure 19 displays the frequency of the bucket-supported OWT before and after scour (Sd = 0.5 L).The soil stiffness decreases after scour, which impacts the stiffness of the bucket-supported offshore wind turbine system, resulting in a decrease in the firstorder natural frequency.The dynamic behavior of the bucket-supported OWT is relatively "softer" under scoured conditions.

Seismic Responses
Figure 20 presents response under same Acc100 (0.15 g) but with different scour depths, including 0.0 L, L, and 0.5 L. trends of the envelopes but the response values differ under variable scouring conditions.It is worth noting that with increasing scour depth, seismic of the OWT does not necessarily increase.In Figure 20a and 20b, the maximum acceleration increases by 8.62% when the scour depth reaches 0.5 L, while the maximum displacement decreases by 3.92%.Figure 20c shows that the maximum bending moment is still near the mudline, and it increases by 13.23%.However, under different seismic motions, the changes in seismic responses with increasing scour depth also vary. Figure 21 shows the normalized maximum responses for accelerations (Figure 21a), displacements (Figure 21b), and bending moments (Figure 21c).When the scour depth increases to 0.5 L, the maximum accelerations at the RNA increase by 8.65%, 3.79%, and 2.64% for the Kobe, Kocaeli, and Acc100 excitations, respectively.Conversely, under Kobe and Acc100 excitations, the maximum displacements at the RNA decrease by 8.10% and 9.63%, respectively, while they increase by 6.48% under Kocaeli motion.For the bending moment at the mudline, the maximum increase is 2.1%, 6.28%, and 13.23%, respectively, as shown in Figure 21c.It can be confirmed that scour alters the dynamic characteristics of the system, thereby changing the seismic responses of the OWT.This change advantageous, but in most remains disadvantageous, depending on the spectral characteristics of the seismic excitations.Therefore, it is imperative to account for the scour effect in wind turbine seismic design to mitigate potential risks.

Conclusions
Numerous offshore wind turbines (OWTs) with bucket foundations have been installed in seismic regions.Compared to the relative development of monopiles (widely installed), seismic design guidelines for bucket-supported OWTs still need to be developed.Moreover, the scour phenomenon around bucket foundations induced by watercurrent actions also creates more challenges to the seismic design of OWTs.In this study, a simplified seismic analysis method is proposed that incorporates the soil-structure interaction (SSI) for the preliminary design of scoured bucket-supported OWTs, aiming balance accuracy and efficiency.It is worth noting that this paper extends the groundwork established in earlier research [36,[44][45][46].Compared with monopiles [44], this study focuses on bucket foundations and employs a lumped parameter model (LPM) to consider the SSI effect.The model is developed by fitting frequency-dependent impedance functions obtained from the four-spring Winkler model.Additionally, the water-structure interaction is also considered by including additional mass to the transition segment in water during seismic analysis, where the transition is simplified to a frustum-shaped foundation.Then, an integral numerical model of the bucket-supported OWT is constructed based on the OpenSees platform, and the validity of the model is confirmed through comparisons of dynamic characteristics and seismic responses with those obtained from the three-dimensional finite element method.Through a series of parametric analyses, the following conclusions can be drawn:

•
As the amplitude of the earthquake decreases from 0.30 g to 0.15 g and 0.05 g, the maximum acceleration, displacement, and bending moment increase by approximately 67% and 83%, respectively.

•
The seismic responses under the Kocaeli motion (characterized by abundant low-frequency components) are more pronounced than those under the Kobe and Acc100 excitations.The seismic response is amplified when the prominent frequency of motion closely aligns with the natural frequency of the OWT.

•
Compared with the monopile-supported OWT, the bucket-supported OWT demonstrates greater dynamic impedances of foundations and first-order natural frequencies.Additionally, the maximum acceleration, displacement, and bending moment of the monopile-supported OWT increase by 43.9%, 30.3%, and 19.4%, respectively.

•
Increasing the scour depth (from 0 to 0.5 L) results in a decrease in the soil stiffness, causing alterations in the impedance functions and frequency characteristics of the OWT.Moreover, with increasing scour depth, the seismic response of the OWT changes significantly.The effect may be advantageous, but in most cases, it remains disadvantageous, depending on the spectral characteristics of the seismic excitations.Therefore, it is imperative to account for scour phenomena in wind turbine seismic design to mitigate potential risks.
In this study, a linear elastic assumption is employed to derive the dynamic impedance functions of the foundation, aiming to streamline calculations and enhance efficiency.However, nonlinear behavior actually exists not only in soil materials but also in the geometric interface between the foundation and soil, which becomes particularly pronounced during strong earthquakes.Therefore, further research and optimization of the proposed method are necessary.Subsequent adjustments, informed by experimental data, can be made to modify the stiffness and damping.This modification partially accounts for the nonlinear characteristics, thereby enhancing the accuracy of seismic analysis results for wind turbines.

Figure 1 .
Figure 1.Schematic of the bucket-supported offshore wind turbine: (a) prototype and (b) numerical model based on a lumped parameter model.

Figure 2 .
Figure 2. Four-spring Winkler model for the impedance of bucket foundations.Based on the D'Alembert principle, the vibration equation for the bucket foundation is derived through force equilibrium: b b b b b b b K denote the mass, stiffness, and damping matrices of the bucket foundation, respectively.b P is the external force vector, as expressed in Equation (2).This study focuses on horizontal seismic excitation, and responses are obtained in the horizontal direction.b u and b θ represent the displacement and rotation at the foundation base, respectively.m , L , and c J represent the mass, total length, and mass moment of iner- tia, respectively.
HH K and HH C represent the horizontal dynamic stiffness and damping coefficient of the embedded foundations, respectively; tw I represents the horizontal embedment fac- tor (Gerolymos and Gazetas, 2006); , emb h χ denotes the horizontal dynamic stiffness coeffi- cient of the embedded foundations; and should be highlighted that the simplified formulas apply to foundations with 0

m
represents the additional mass per unit volume; r denotes the foundation radius; 0 r is the top radius of the foundation; 0 θ represents the slope angle; z indi- cates the distance of the segment element from the mudline; w h represents the water depth; and w ρ is the water density.

Figure 4 .
Figure 4.Additional mass density along the foundation height: (a) schematic representation and (b) distribution curve.

Figure 5 .
Figure 5. Erosion topographic map for the bucket foundation.

Figure 7 .
Figure 7.Comparison between the results from the present study and FEM: (a) natural frequencies, (b) acceleration, and (c) displacement response envelopes along the OWT.

Figure 8 .
Figure 8.(a) Fourier spectra and (b) acceleration time histories of the earthquakes.

Figure 9 .
Figure 9.Time history responses corresponding to accelerations, displacements, and bending moments under different seismic intensities.

Figure 10 .
Figure 10.Response envelopes corresponding to accelerations, displacements, and bending moments along the OWT under different seismic intensities.

FigureFigure 12 .
Figure corresponding mo-along the OWT under different spectral characteristics.

Figure 13 .
Figure 13.Schematic of the physical models: (a) monopile-supported offshore wind turbine and (b) bucket-supported offshore wind turbine.

Figure 14 .
Figure 14.Impedance functions and LPM approximations of the bucket foundation (a-c) and the monopile (d-f).

Figure 15 .
Figure 15.Mode frequency comparison between the bucket-supported OWT and the monopile-supported OWT.

Figure 16 .Figure 17 .
Figure 16.Time history responses corresponding to accelerations, displacements at the RNA, and bending moments at the mudline for OWTs with different foundation types.

Figure 18 .
Figure 18.Impedance functions and LPM approximations of the bucket foundation (a-c) before scour and (d-f) after scour.

Figure 19 .
Figure 19.First-order natural frequency of the OWT before and after scour.

Figure 20 .Figure 21 .
Figure 20.Response envelopes corresponding to accelerations, displacements, and bending moments along towers under different scouring conditions.
Interest: Author Panpan Shen was employed by the company Shanghai Investigation, Design & Research Institute Co., Ltd.The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Table 1 .
Key parameters of the reference bucket-supported offshore wind turbine.

Position Parameter Type 6.45 MW
MSL: mean sea level.