Fatigue assessment of flange connections in offshore wind turbines under the initial flatness divergence

Bolted ring flange connections are widely utilized in offshore wind turbines to connect steel tubular segments. After the massive production and installation of offshore wind turbines in the past decade, flatness divergence is regarded as one of the most important initial imperfections for the fatigue design of flange connections. Offshore wind turbines are subjected to wind, wave, and current loads. This initial imperfection may alter the structural response and accelerate the fatigue crack growth. This paper aims to analyse the impact of the initial flatness divergence on the structural response of flange connections and evaluate its consequences on fatigue damage. Two different offshore wind turbines with fixed foundations and floating foundations are modelled to simulate their global responses to environmental loads. Based on a superposition method, local finite-element models of flange connections are established with three types of flatness divergence. Using the same bolt pretension and external loads from global modelling, the impact of these geometric imperfections is further examined by comparing the structural responses of flanges under different radial and peripheral opening lengths. Then, the fatigue assessments on flange connections in both fixed wind turbines and floating wind turbines are conducted, and the impacts of initial flatness divergence on these two different wind turbines are analysed.


Introduction
During mass production and installations of offshore wind turbines (OWTs), in the wind industry, the bolted ring flange connections for assembling the OWT substructure and tubular steel towers have been widely used in the last few years (Pavlović et al., 2015a;Pavlović et al., 2015b). Throughout the service life of OWTs, the variable and cyclic loads from wind and waves sustained by wind turbine towers make them vulnerable to fatigue damage. Therefore, it has often been reported on fatigue cracks in wind turbine towers (Chou and Tu, 2011).
According to the location of openings, flatness divergence can be classified in flange-sided gapping, tower-sided gapping, and parallel gapping, as shown in Figure 1B. When preloading, these gaps close initially, and then, clamp solids emerge. Ji et al. (2021) and Long et al., 2021) established three types of tower flange gap to discuss the effects of the flange gap on the cumulative fatigue damage. They revealed that the special shape of the gap possesses a significant impact on fatigue damage. Jakubowski and Schmidt (2003) performed systematic experimental and numerical analyses on the influence of flange imperfections on the load transfer function (LTF), and the result showed that it was most difficult to close openings with parallel gaps, while the tower-sided gaps appeared to be less difficult. In contrast, the flange-sided gaps had a positive effect on the LTF, since it caused the resulting clamp force to move closely towards the tower wall and, thus, the resulting axis of external forces. This conclusion was also confirmed by Feldmann et al. (2011).
In design practice, the LTF according to Schmidt-Neuper does not consider the various geometric flange imperfections shown in Figure 1B explicitly, but it covers a certain (but unknown) range of gap sizes and shapes implicitly. In order to account for these gap sizes and shapes, according to DNVGL- ST-0126 (DNVGL, 2018) and the DIN 18088-3 (DNVGL, 2019), the flatness divergence of one flange must be less than 2 mm in length over the entire flange circumference and less than 1 mm over a 30°segment after completing the tubular sections. However, Weijtjens et al. (2021) proposed that this requirement is semi-empirical and it did not consider those geometric imperfections explicitly. With bolted flange connections being widely used within OWT structures, which are quite large in diameter and thickness, it becomes additionally challenging to fulfil these requirements without a major effort during the OWT installations.
This research presents a numerical investigation of the effect of the initial flatness divergence on maximum stress and fatigue damage in the ring flange connections within fixed-bottom wind turbines (FBWTs) and floating offshore wind turbines (FOWTs). These two different OWTs are modelled to simulate their global responses to environmental loads in operating conditions, and external loads on the tower shell sections are calculated. Then, local finite-element models with the initial flatness divergence are applied to the same bolt pretension and external loads from global modelling. The local flanges' responses are calculated based on LTFs. Finally, the fatigue damage of flange connections in both FBWTs and FOWTs is compared, and the impacts of the initial flatness divergence on these two different wind turbines are discussed.

Global responses of offshore wind turbines
In the past decades, significant research focused on the global response simulations of OWTs, as it was the foundation for any other further analysis. In this section, two different OWTs for both FBWTs and FOWTs are modelled to simulate their global responses to environmental loads using SESAM-SIMA software. Then, the cyclic loadings on the FOWT and FBWT tower shell sections under the operating condition are obtained. The workflow of the global responses of OWTs is shown in Figure 1A.

Description of offshore wind turbines
OWTs include FBWTs and FOWTs. FBWTs are usually operating in shallow water with a depth less than 60 m, and the foundations are inserted into the seabed. FBWTs are mainly subjected to environmental loads, such as wind and waves. Liu et al. (2010) considered that wind loads played a leading role in load combination for FBWTs' tower. Different from FBWTs, FOWTs consist of a floating platform and restrained with various mooring systems that typically undergo displacement in the six degrees of freedom (Zou et al., 2021). Furthermore, they may undergo remarkable dynamic responses when operating in the hostile marine environment. FOWTs studied in this paper are the 5 MW OC4-DeepCWind semi-submersible wind turbine, as shown in Figure 2. The floater includes a centre column connecting the tower and three offset columns which are connected with the main column through a number of smaller pontoons and braces (details of the platform are given in Figure 2B). Three mooring lines are arranged symmetrically above the platform vertical axis with 120°between them. The radius to fairleads and anchors from the platform's centreline is 40.868 m and 837.6 m, respectively. Each of the three lines has an unstretched length of 835.5 m with 0.0776 m diameter and an equivalent mass per unit length of 113.35 kg/m. The detailed layout and structural property are available in Robertson et al. (2014).

Calculating dynamic responses
OWTs have to resist winds, waves, and currents during the operation stage. Therefore, the environmental loads should be estimated before calculating the dynamic responses of OWTs (Hu et al., 2020). All these environmental loads induce an aerodynamic and hydrodynamic impact on the structures. Ye and Ji (2019) investigated the effects of both hydrodynamic and aerodynamic excitations which include current, wave, and wind excitations, as well as buoyant forces, along with the dynamic interaction between the drive-train system and tower structure on the dynamic response of the spar-type floating platform under different sea conditions. Chen and Basu (2018) considered the effects of current and wave-current interactions in the fatigue analysis of FOWTs. Their results showed that the current and the wave-current interaction could have significant influences on FOWT cable responses. Philippe et al. (2013) performed a coupled dynamic analysis of a floating wind turbine system to investigate the effect of the wave direction relative to wind on the wind turbine system. They considered that the natural modes of the tower system are excited differently regarding the wave's direction.
In this research, wind and wave loads are considered as primary environmental external forces. Despite that the current may also affect the FOWTs' motions, especially on supporting the structures of an extensive draft, it mainly acts on the mooring systems and contributes less to fatigue damage on the towers (Zou et al., 2022). Therefore, the current force is not discussed in this paper.
The external loads on the ring flange connections are calculated, as shown in Figure 2C. The wind and wave spectra are set as stationary uniform and JONSWAP for both FOWTs and FBWTs, respectively. The significant wave height is set as 3 m, spectral peak period is set as 10 s, and horizontal wind speed is set as 11.4 m/s.
The coupled global dynamic analysis carried out the time domain using the module SIMA in software package SESAM. The SESAM-SIMA calculates the internal forces and moments in the cross sections of the FOWTs' and FBWTs' tower under the same operation condition. Figure 3 shows the time series of the axial force and the bending moment in a cross section of the tower. The results show that the ranges of fluctuations in axial force and those of bending moment of the FOWTs' towers increase compared to FBWTs' towers, and FBWTs' towers bear a greater axial force and bending moment compared to FOWTs' towers under operating conditions. It is clearly observed that FOWTs' substructures undergo oscillating motions in six degree of freedoms (DOFs) with larger magnitudes compared to fixed-bottom wind turbines. Their compliant foundations can reduce the

Frontiers in Energy Research
frontiersin.org 03 counterforce around the tower and lower the stress in the structural response. Moreover, the substructure's motions may influence the aerodynamic condition of the rotor, which makes the aerodynamic loadings on FOWTs much more complicated than those on FBWTs. de VaalHansen and Moan (2012) studied the effect of a periodic surge motion on the integrated loads and induced velocity on a wind turbine rotor through numerical simulations. Among all motions of the six DOFs, surge and pitch are regarded as the most important ones affecting the FOWTs' aerodynamics (Bayati et al., 2016;Fang et al., 2020;Fang et al., 2021). Once the aerodynamic conditions are affected by the substructures' motions, the external forces on the wind tower are also changed. Hence, it is necessary to consider the impact of FOWTs' motions in the fatigue assessment.

Finite-element modelling of the connection
The finite-element method (FEM) model is established to analyse the stress distribution around flange connections. To simplify the modelling and minimize the computational cost, a segment of the connection that included one bolt and experienced the maximum load was isolated and modelled through ANSYS parametric design language (APDL) in this research. The following sections describe the structure of the FE models.

Material properties and geometry
The FE analysis was performed based on linear material models. The materials were modelled considering the steel's elastic

Finite elements and the contact model
A finite-element model should not only have a low computational requirement but also reflect the key characteristics, such as irregular geometries and large deformations. In this study, SOLID186-type finite elements are used for the flange. SOLID186 is a 3D element with 20 nodes, and three degrees of freedom can accurately simulate irregular geometries and structural phenomena, namely, plasticity, creep, and large deformations, among others (System, 2005). The mesh density has an effect on finite-element results. A good grid layout is one of the key factors to improve the reliability of simulation when using the finite-element method for structure simulation. The coarse mesh density can decrease the computation time, but it may lead to a larger finite-element stress result that is not consistent with the reality. On the contrary, fine finite elements are used to improve the accuracy of the finiteelement results with more computation time. Li et al. (2009) put forward the concept of key areas, and it was concluded that the precision of the calculation result of the whole finite-element model was determined by the meshing level of all the key areas in the finiteelement model. In this section, the mesh density to ESIZE is set as 6 in the mesh tool initially, and Figure 5 shows that the stress value of the flange is too high. With the further use of small mesh density, the tendency of the flange stress is clearly shown approaching a stable value indicating that the mesh density of 0-1 is reasonable. In this paper, the mesh size was chosen to be a minimum of 0.9 mm and a maximum of 16 mm, with a total of 107,242 nodes and 14,997 elements. Figure 6 (left) shows a localised fine mesh around the areas with high stress concentration, while a relatively coarser mesh size is adopted for other areas. In addition, a pretension element (PRETS179 element in ANSYS) is used to apply the preload force on the bolts. The pretension elements are located in the middle of the bolt shaft. Then, ANSYS internally cuts the bolt into two parts and translates the cut planes axially so that the sum of the nodal reaction forces is equal to the preload force.
The local model was created with the configuration of the contacts between elements and the application of boundary conditions, which consisted of a frictional contact between the surface of the upper and lower flanges, washers and flanges, inner surface of flanges and the bolt shank, and a bonded contact between the surface of nuts and washers and inner surface of nuts and the bolt head. The friction coefficient applied was 0.15, adopting an augmented Lagrangian equation (Simo and Laursen, 1992).

Loads and boundary conditions
To assess the fatigue damage on OWTs' tower flange, the calculation of the cyclic load amplitudes ΔF b requires the approximation of the load transfer between the tower shell loads and the resulting bolt force inside the flange. According to the current design standards, this is performed using a load transfer function (LTF). Figure 7 (left) schematically shows the LTF F b (Z), with the external load Z acting in the tower shell on the abscissa and the bolt load F b on the ordinate. Due to the eccentric load introduction and the presence of the preload F p of the bolts, the LTF is non-linear.
For the reduction of fatigue loads, bolt pretension F p is necessary for structural integrity. The preload F p is induced via either torquecontrolled tightening or direct tensioning methods. With the broadly used torque-controlled tightening procedure, the bolt shall be preloaded up to the reduced nominal preload, F p,c* = 0.7R p0.2 ·A sp, where R p0.2 is the material plastic strain limit and A sp is the stress area of the bolt thread.
Under perfect geometric conditions, as the surface pressure between the upper and lower flange increases due to the tightening, a clamp solid emerges centrically around the longitudinal bolt axis as a counterpart of the preload, as shown in Figure 7 (middle). With the clamp solid distributed symmetrically around the bolt axis, external forces acting in the tower shell first diminish the resulting clamp force, which reduces cyclic bolt forces ΔF b . However, under imperfect flatness divergence conditions, the flange contact surface is not ideally plane due to a flange inclination α s , a flange opening length L i , and a flange opening height h, which is shown in Figure 7 (right). The two flange parts may open, and the clamp solid may develop eccentrically from the bolt's longitudinal axis.
The flange connections are mainly subjected to external loads, such as moments, pressure, and shear. The axial stress applied to the upper tower flange shell sections is calculated as follows (Eq. 1). The shear is neglected due to its minor contribution to fatigue damage.

Frontiers in Energy
Here, Z is the external load, β is the angle between the location of the bolt and the x-axis, F axial is the maximum axial force (plotted in Figure 3), M y is the y-axis bending moment (plotted in Figure 3) above the horizontal axis perpendicular to the wind direction, A t is the area of the tower shell, R is the radius of the flange, and N is the number of high strength bolts installed in the flange (plotted in Figure 1).
A preload force of 510 KN is applied on the bolt in order to simulate the tightening of the bolt (Ajaei and Soyoz, 2020) Figure 9 shows that bolts' and flanges' stress concentrations occur at the bolt head under the perfect condition. Figure 10 shows that imperfect models' stress concentration mainly occurs at the bolt head and the centre of the bolt shank. With the increase in opening length, the locations of stress concentrations are transferred from the bolt head to the centre of the bolt shank. For flanges with flangesided gapping and tower-sided gapping (Figure 11 (1) (2)), stress concentrations initially occur at the bolt head and transfer to the opening gaps with the increase in opening length. For parallel gapping, a three-sector local finite-element model is established, as the contact surface between the two flange parts opens completely in the middle and the middle bolt loses its preload completely. Hence, the cyclic loadings on the FOWTs' tower shells are mainly sustained by other two bolts, and the locations of stress concentrations occurring in the bolt head and gaps are shown in Figure 11 (3). These locations may play a critical role in fatigue assessments. Figure 12 schematically shows the variation tendency of the bolt and flange stress increment under flatness divergence. It is evident that the magnitude of the bolt and flange stress increased considerably when the opening length increases. The main findings are drawn as follows:

Effect of gap length on bolt and flange stress
1. When the tower shell is subjected to tensile stress, bolt and flange stresses significantly increase in flange-sided gapping compared to the other two types of flatness divergence. The maximum stress of bolts is 2763 MPa (Figure 10), where (1 with a maximum flange opening length, and the stress increases to 1787 MPa ( Figure 12A) for flange-sided gapping compared to the perfect

FIGURE 8
Details of the imperfect model.

FIGURE 9
Von Mises stress state of the OWTs' tower flange under the perfect condition.
Frontiers in Energy Research frontiersin.org condition (Figure 9, right). The bolt stress increment is about 5% for tensile when the gap length is 49 mm (point A in Figure 8), and the bolt stress increment is about 183% for tensile when the gap length is 115 mm (point E in Figure 8) for flange-sided gapping. For flanges, the maximum stress is 9858 MPa [ Figure 11 (1)] when the location of openings occurs in the centre of the bolt shank, i.e., point C (gap length 82 mm) in Figure 8, and the stress increases to 8610 MPa for flange-sided gapping (flange stress increment is about 690%) compared to the perfect condition (Figure 9, right).

FIGURE 10
Von Mises stress state of the FOWTs' tower bolts under the initial flatness divergence when subjected to tensile (left) and compressive stresses (right) (details of the location of openings are given in Figure 8).

FIGURE 11
Von Mises stress state of the FOWTs' tower flanges under initial flatness divergence when subjected to tensile (left) and compressive stresses (right) (details of the location of openings are given in Figure 8).

Frontiers in Energy Research
frontiersin.org 2. The bolts and flanges are more sensitive to tower-sided gapping when the tower shell sections are subjected to compressive stress. Compared to the perfect condition shown in Figure 9 (right), the increase of bolts' stress is 2380 MPa (Figure 12B), and the flanges' increase is 9816 MPa ( Figure 12B). The bolt stress increment is about 11% for compression when the gap length is 35 mm (point E in Figure 8), and the bolt stress increment is about 246% for compression when the gap length is 101 mm (point A in Figure 8) for tower-sided gapping. For flanges, the maximum stress increment is about 789% when the location of openings occurs in the centre of the bolt shank, i.e., point C (gap length 68 mm) in Figure 8 compared to the perfect condition (Figure 9, right). 3. When the openings become larger, the contact surfaces of the upper and lower flanges are finally completely separated. The flange connection appears less sensitive to parallel gapping; this is due to the fact that the bolt in the middle segment has failed, and the cyclic loadings are mainly sustained by flanges. The increase of flanges is 541 MPa ( Figure 12C) compared to the perfect condition. These tendencies indicate that initial flatness divergences can result in a significant increase of stresses, which should be taken into account in fatigue assessments. Figures 9-11 show that the stress of both the bolt and the flange exceeds the materials' yield strength. In the case of the model with flatness divergence, the structure is subjected to deformation or even destruction. The plastic deformation may occur at this bolt, and it will lose its resistance against external forces. The external force on this bolt will be transferred to the surrounding bolts. In this paper, the structure response after plastic deformation is not considered. Instead, the influence of the imperfections on the structural stress is mainly discussed. When using the finite-element method for meshing, improper meshing may lead to stress concentration. The nonconvergence of stress can result in an extremely higher stress compared to the real structural stress (as shown in Figure 6). Therefore, the stress extrapolation method is generally used in finiteelement calculation. Alonso-Martinez et al. (2019) applied a 520 KN preload force on the M36 bolt and calculated the maximum von Mises stress of local finite-element models of flange connections with the perfect model by ANSYS Workbench. The result showed that the maximum stress of the flange was about 4700 MPa. It also exceeds the yield strength of the materials. Ji et al., 2021 andLong et al. (2021) established three types of tower flange gap to discuss the effects of the flange gap on the bolts stress. The result revealed that the bolt stress was more sensitive to flange-sided gapping when the tower shell was subjected to tensile stress, and the bolt stress was more sensitive to tower-sided gapping when the tower shell was subjected to compressive stress. This result is consistent with the findings in Figure 12.

Effect of gap sizes on different OWTs' bolts
The stress difference of bolts is not significant between FOWTs and FBWTs at the beginning. When the location of openings occurs at the centre of the bolt shank, i.e., point C in Figure 8, the bolt stress of FBWTs is observed to be significantly affected by the tower-sided gapping and flange-sided gapping compared to FOWTs (as shown in Figure 13). This trend is due to external forces Z acting on the FBWTs' tower shell larger than those on the FOWTs' tower shell. This large external force Z leads to a higher bolt load F b to some extent. With the openings becoming large until the flange contact surface completely opens, the bolts cannot exert an adequate clamping force to keep the joint together. Hence, parallel gapping  Figure 8).

Frontiers in Energy Research
frontiersin.org has a weak influence in the middle bolt but has a greater impact on the neighbouring bolts for both FOWTs and FBWTs. For the tower-sided gapping condition, the maximum stress of FBWTs' bolts is about 1200 MPa and the maximum stress of FOWTs' bolts is about 1800 MPa when the tower shell is applied with tensile stress. In addition, when the tower shell is subjected to compressive stress, the maximum stress of FBWTs' bolts is about 4300 MPa and the maximum stress of FOWTs' bolts is about 3300 MPa. For the flange-sided gapping condition, the maximum stress of FBWTs' bolts is about 3000 MPa and the maximum stress of FOWTs' bolts is about 2700 MPa when the tower shell is applied with tensile stress. When the tower shell is subjected to compressive stress, the maximum stress of FBWTs' bolts is about 2200 MPa and the maximum stress of FOWTs' bolts is about 2400 MPa.

Discussion
As can be observed, these geometric imperfections with the initial flatness divergence have significant effects on both flange and bolt stresses. Typical load transfer and potential reasons for the impacts between flanges and bolts are as follows: a) For ring flanges, stress concentrations occur in the bolt head and the location of geometric imperfections. When geometric imperfections occur in the bolt centre, the flange stress may have a sharp rise due to the bolt failure. b) For bolts, stress concentrations often occur centrically around the longitudinal bolt axis and bolt head. With the increase in opening length, the locations of stress concentrations are transferred from the bolt head to the centre of the bolt shank. c) For offshore wind turbines, the case with tower-sided gapping has a more significant impact on FBWTs' bolts, while FOWTs' bolts are more sensitive to flange-sided gapping.

Fatigue analysis and results
Fatigue damage is a cumulative process caused by cyclic loads (Zou and Kaminski, 2020). Initial imperfections are one of the most critical factors for fatigue strength in bolted connections, except for a vulnerable material, cyclic tensile loads, and higher stress levels (Mehmanparast et al., 2020). This study focuses on the impact of flatness divergences on bolted connections' fatigue damage, including flange-sided gapping, tower-sided gapping, and parallel gapping. Therefore, it is assumed to have a 20-year design life of bolts, i.e., 1e7 number of cycles, and conducts the fatigue assessments for flange connections in both fixed wind turbines and floating wind turbines based on the accumulation of fatigue damage computed with nominal stresses and S-N curves.

Global-local methodology
In this section, a global-local methodology of the fatigue assessment procedure is utilized. This method emphasises the theoretical concepts associated with the modal superposition and fatigue analysis. A workflow for the global-local fatigue assessment methodology is shown in Figure 14. It introduces two different phases of analysis from the global to the local scale with increasing details of geometrical, material, and contact properties. The two phases of fatigue analysis are systematized as follows. In the first phase, a global approach based on the accumulation of damage computed with nominal stresses and S-N curves is introduced. The adopted global numerical model should be accurate enough to obtain nominal stresses to be used as input loading for the linear damage accumulation method. Since the available S-N curves with nominal stress may not properly reflect the local geometrical and material characteristics of critical details, a conservative S-N curve for the fatigue resistance is, thus, considered. Based on the fatigue analysis result, the critical details should be identified, which require a more refined assessment by implementing local scale models. Then, in the second phase, submodelling techniques leveraged by modal superposition concepts are used to accurately evaluate local fatigue damage.

S-N curves and Miner's rule
For an offshore structure design, fatigue damage is usually calculated based on the S-N curve approach. The fatigue resistance is represented by S-N curves (Gao et al., 2021). S-N curves define the relationships between the stress ranges and the numbers of cycles to failure in those stress ranges. GL 2010 (Wind guideline, 2010) provides S-N curves for different structural materials, including a bolt with a diameter larger than 30 mm (M36 bolts) used in the bolted connection, which is considered in this paper. The basic design of the S-N curve is given as where N = predicted number of cycles to failure for stress range where n i is the number of the existing stress cycles in each range and N i is the number of stress cycles in each range, which can cause fatigue failure. It is usually assumed that failure occurs when the damage index D reaches 1.0. Figure 15A shows the fatigue damage on FOWTs' and FBWTs' bolts with a maximum flange opening length within 20 years. It is observed that the existence of flatness divergence in bolted flange connection increases the fatigue damage severely. Moreover, it makes a significant impact to consider the different types of flatness divergence in fatigue assessments, especially under towersided gapping and flange-sided gapping conditions. On the contrary, the fatigue damage of the middle bolt appears not so sensitive to parallel gapping which may significantly affect the neighbouring bolts. Compared to FOWTs, FBWTs' bolts are more vulnerable to fatigue damage due to higher external loads and stress range

FIGURE 14
Workflow of the global-local fatigue assessment methodology for offshore wind turbines.

Frontiers in Energy Research
frontiersin.org (the stress range of bolts is calculated by stress extrapolation method), as shown in Figure 15B. The flexible foundations of FOWTs can reduce the counterforce around the tower against the environmental loads and result in lower structural responses. The fatigue damage of FBWTs' bolts is 0.86, which is about 1% larger than that of FOWTs, which is 0.85, when the tower is subjected to tensile stress and the fatigue damage of FOWTs' bolts is 0.65, which is about 2% larger than that of FBWTs, which is 0.63, when the tower shell is subjected to compressive stress for flange-sided gapping. Furthermore, the fatigue damage of FOWTs' bolts is 0.5, which is about 2% larger than that of FBWTs, which is 0.48, when the tower shell section is subjected to tensile stress and the fatigue damage of FBWTs' bolts is 1.16, which is about 18% larger than that of FOWTs, which is 0.98, when the tower shell is subjected to compressive stress for tower-sided gapping. Hence, more attention to FBWTs' bolts with tower-sided gapping under the compressive stress is necessary during the design process for safe operation.

Conclusion
This paper analysed the impact of the initial flatness divergence on the structural response of flange connections in OWTs and evaluated its consequences on fatigue damage.
The effects of different opening lengths on flanges' and bolts' stresses were compared by establishing local finite element models with three types of flatness divergence. The fatigue assessments for flange connections in both fixed wind turbines and floating wind turbines were conducted. The impacts of the initial flatness divergence on these two different wind turbines were compared. According to the results of the finite-element analysis, the main conclusions were drawn as follows: 1. The ranges of force fluctuations in the FOWTs' towers are larger than in FBWTs. FOWTs' substructures were not fixed to the seabed; they had oscillating motions in the six DOFs and would experience larger motions than the bottom-fixed wind turbines. The flexible foundations of FOWTs could reduce the counterforce around the tower against the environmental loads, resulting in lower structural responses. The maximum tensile stress of the FBWTs' tower shell can be 50% larger than that of FOWTs, and the maximum compressive stress of the FBWTs' tower shell is shown to be 51% larger than that of FOWTs. 2. The stress of bolted flange connection can increase significantly with the increase in the opening length, and stress concentrations may occur in the bolt head, the centre of the bolt shank, and gaps. Thus, particular attention should be paid to these locations for safe operation. For the flange-sided gapping, the stress of bolts and flanges increase significantly when the tensile stress is applied to the tower shell. Compared to the perfect condition, the flangesided gapping is markedly increased and demonstrated for the maximum bolt stress and the maximum flange stress. For the tower-sided gapping, the compressive stress has a significant influence on the bolt and flange with significant increases in the maximum bolt stress and the maximum flange stress, respectively, compared to the perfect condition. The bolt stress increases by 40 MPa and the flange stress increases by 541 MPa, which is not sensitive to parallel gapping. 3. The very significant negative effect of the initial flatness divergence could be confirmed based on fatigue analysis. The local finite-element model with tower-sided gapping and flangesided gapping demonstrated a significant increase in the fatigue damage of bolts compared to parallel gapping. For the flangesided gapping, the fatigue damage of the FOWTs and FBWTs is approximately the same. For the tower-sided gapping, however, the fatigue damage of the FOWTs is markedly lower than that of FBWTs. Compared with parallel gapping, the fatigue damage of FBWTs is almost the same as that of FOWTs.

Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions
XJ: conceptualisation, software, investigation, formal analysis, and writing-original draft. TZ: first corresponding author, methodology, writing-review and editing, and validation. XB: second corresponding author and writing-review and editing. XN: picture making. LT: writing-review and editing. All authors contributed to the article and approved the submitted version.