Fatigue Assessment of Moorings for Floating Offshore Wind Turbines by Advanced Spectral Analysis Methods

: The fatigue assessment of mooring lines for ﬂoating offshore wind turbines represents a challenging issue not only for the reliable design of the stationkeeping system but also for the economic impact on the installation and maintenance costs over the entire lifetime of the offshore wind farm. After a brief review about the state-of-art, the nonlinear time-domain hydrodynamic model of ﬂoating offshore wind turbines moored by chain cables is discussed. Subsequently, the assessment of the fatigue damage in the mooring lines is outlined, focusing on the combined-spectrum approach. The relevant fatigue parameters, due to the low-and wave-frequency components of the stress process, are estimated by two different methods. The former is based on the time-domain analysis of the ﬁltered stress process time history. The latter, instead, is based on the spectral analysis of the stress process by two advanced methods, namely the Welch and Thomson ones. Subsequently, a benchmark study is performed, assuming as reference ﬂoating offshore wind turbine the OC4-DeepCWind semisubmersible platform, equipped with the 5 MW NREL wind turbine. The cumulative fatigue damage is determined for eight load conditions, including both power production and parked wind turbine situations. A comparative analysis between time-domain and spectral analysis methods is also performed. Current results clearly show that the endorsement of advanced spectral analysis methods can be helpful to improve the reliability of the fatigue life assessment of mooring lines.


Introduction
The reliability of long-term mooring systems for floating offshore wind turbines (FOWTs), deployed on intermediate and deep-water depth, is a key factor to move the offshore wind energy sector towards the commercialization phase, provided that in the last two decades several accidents, mainly due to mooring failures, occurred throughout the world [1].Indeed, one of the expected goals of the Horizon 2020 "Secure, clean and efficient energy" work programme [2] is the development of reliable, sustainable and cost-efficient mooring systems for FOWTs, with the main aim of making the offshore wind energy sector more competitive on the international market, as regards other well-established renewable sources.
As previously said, the need of further improving the hydrodynamic modelling of FOWTs, as well as the design of long-term moorings, is due to the high number of failure events that occurred in the last decade, as stressed by Kvitrud [3], among others, who reported 15 mooring failures in the Norwegian sea area during the years 2010-2013, onethird of which were due to single or multiple fatigue failures of mooring lines.Starting from the pioneering research by de Laval [4] and van Helvoirt [5], who performed the first fatigue tests on mooring chains at the beginning of 1970s and 1980s, respectively, the fatigue analysis of mooring lines quickly became a very popular topic, following the growth of the oil and gas industry, the need of exploiting subsea oilfields and the increasing demand for floating offshore structures to be deployed on deep-water depths.Past research activities focused on a variety of topics involved in the fatigue analysis of mooring systems, among which were the following: the simultaneous presence of low-and wave-frequency components in the mooring line stress process [6], the occurrence of bimodal non-Gaussian random processes [7], the application of spectral analysis and rainflow counting methods for the fatigue analysis of mooring lines [8].
More recently, following the increasing demand for energy from renewable sources and the rapid development of the offshore wind technology, attention was also paid to the design and analysis of stationkeeping systems for FOWTs deployed on intermediate and deep-water depths.Most of past research activities in this field focused on the design and optimization of mooring systems, with the main aim of reducing the installation costs of wind farms and making the offshore wind energy sector more competitive on the international market, in terms of the levelized cost of energy.Particularly, Brommundt et al. [9] provided a frequency-domain tool for the optimization of mooring systems and applied it to a reference semisubmersible tri-floater support structure deployed on two different water depths, namely 75 and 330 m.Benassai et al. [10,11] carried out a benchmark study, devoted to investigating the incidence of the main design parameters, namely the number and type of mooring lines, on the total weight of the stationkeeping system at different water depths.Kim et al. [12] designed the stationkeeping system for an FOWT, deployed in the Jeju offshore basin.Hall and Goupee [13] introduced a lumped-mass mooring line model for FOWTs and validated it against a set of experimental tests on a tri-floater platform.Campanile et al. [14] focused on the design of mooring systems for FOWTs to be deployed on intermediate water depths and furnished some suggestions to reduce the weight of the stationkeeping system from the preliminary design phase.
Following this brief review about the state-of-art concerning the analysis of mooring systems for FOWTs, only few attempts were devoted to the fatigue analysis of mooring lines, as most of past research focused on the ultimate and accidental limit state design conditions for stationkeeping systems [10][11][12][13][14]. Really, the assessment of the fatigue damage in the mooring lines of FOWTs still presents some challenging issues, mainly related to the combined presence of both low-and wave-frequency stress components that, in turn, need to be separately evaluated to assess the cumulative fatigue damage in the mooring line.In this respect, the actual rules and guidelines for the design of stationkeeping systems [15] suggest the employment of the combined-spectrum approach when bimodal stress processes are expected to occur.Nevertheless, one of the main issues that arise when this method is applied is mainly due to the assessment of the low-and wave-frequency components, starting from the time history of the combined stress process.In this respect, starting from the pioneering work by Winterstein [16] in 1988, some years later Jiao and Moan [17] developed the so-called dual narrow-band method for the assessment of the fatigue damage for bimodal stress processes and this is still embodied in current rules and guidelines for the design and analysis of mooring systems [15].In the following years, several spectral methods were developed to improve the accuracy of the fatigue assessment [18,19], when both low-and high-frequency components occur.In the last decade, improved methods were developed for bimodal Gaussian [20][21][22] and non-Gaussian processes [23][24][25][26][27][28], with the main aim of improving the accuracy of the fatigue damage assessment.
Hence, in the current analysis, a time-domain method, based on the application of low-pass and high-pass filters to the initial stress process, is compared with two advanced spectral analysis techniques, based on the Welch and Thomson methods, with the main aim of investigating the incidence of the use of time-and frequency-domain analyses on the assessment of cumulative fatigue damage in the mooring lines.This issue is worthy of investigation not only because the stress process in the mooring lines is generally bimodal but also because the low-frequency component consists of two different contributions, due to the slowly varying drift forces and the turbulent component of the wind loads acting on the emerged part of the support structure and the wind turbine.The paper is structured as follows.Section 2 provides a brief review of the nonlinear time-domain hydrodynamic and mooring models.Section 3 focuses on the fatigue analysis of mooring lines by the combined-spectrum approach.Section 4 discusses the Welch and Thomson methods, applied in the benchmark study.Section 5 provides the main data of the reference FOWT, as well as the basic load conditions, involving both power production and parked wind turbine situations that are analysed in Section 6.Finally, Sections 7 and 8 provide the discussion of the main outcomes gathered by the benchmark study and the conclusions.All calculations were performed by a set of dedicated codes developed in MATLAB.

Model Description of Dynamic Analysis
The reference tri-floater support structure, assumed in the subsequent benchmark study, is the OC4-DeepCWind platform [29], equipped with the 5 MW NREL wind turbine [30].It is equipped with three mooring lines, spread 120 • apart, as depicted in Figure 1.investigation not only because the stress process in the mooring lines is generally bimodal but also because the low-frequency component consists of two different contributions, due to the slowly varying drift forces and the turbulent component of the wind loads acting on the emerged part of the support structure and the wind turbine.The paper is structured as follows.Section 2 provides a brief review of the nonlinear time-domain hydrodynamic and mooring models.Section 3 focuses on the fatigue analysis of mooring lines by the combined-spectrum approach.Section 4 discusses the Welch and Thomson methods, applied in the benchmark study.Section 5 provides the main data of the reference FOWT, as well as the basic load conditions, involving both power production and parked wind turbine situations that are analysed in Section 6.Finally, Sections 7 and 8 provide the discussion of the main outcomes gathered by the benchmark study and the conclusions.All calculations were performed by a set of dedicated codes developed in MATLAB.

Model Description of Dynamic Analysis
The reference tri-floater support structure, assumed in the subsequent benchmark study, is the OC4-DeepCWind platform [29], equipped with the 5 MW NREL wind turbine [30].It is equipped with three mooring lines, spread 120° apart, as depicted in Figure 1.The nonlinear time-domain motion equations, as regards the global reference system with origin on the projection of the FOWT centre of mass on the waterplane area, are provided by Equation (1), which represents the commonly embodied hydrodynamic model for FOWTs [31][32][33]: where  is the mass matrix;  is the added mass matrix at infinite frequency; ̈, ̇ and  are the 6DOF acceleration, velocity and displacement vectors;  is the retardation function matrix;   is the vector of first-and second-order wave excitation forces;   is the wind force vector, including both steady and turbulent components;   is the drag force vector and   is the mooring force vector.The retardation function matrix is determined by the inverse Fourier transform of the frequency-dependent linear radiation damping matrix , according to Equation (2), which is based on the Cummins [34] approach: The nonlinear time-domain motion equations, as regards the global reference system with origin on the projection of the FOWT centre of mass on the waterplane area, are provided by Equation (1), which represents the commonly embodied hydrodynamic model for FOWTs [31][32][33]:

q,
. q and q are the 6DOF acceleration, velocity and displacement vectors; R is the retardation function matrix; F wave is the vector of first-and second-order wave excitation forces; F wind is the wind force vector, including both steady and turbulent components; F visc is the drag force vector and F moor is the mooring force vector.The retardation function matrix is determined by the inverse Fourier transform of the frequency-dependent linear radiation damping matrix B, according to Equation (2), which is based on the Cummins [34] approach: where ω denotes the wave circular frequency.Wave forces include both first-and secondorder components.The former was randomly generated using Equation (3): where A j is the magnitude of a regular incident wave with circular frequency ω j and direction β j , while X (I) is the vector of the first-order complex excitation transfer function, corresponding to the circular frequency ω j and direction β j .As concerns the second-order wave forces, they include both difference-and sum-frequency components among pairs of incoming waves [31].Nevertheless, in the current analysis, only the difference-frequency components, namely the so-called slowly varying drift forces, were considered, according to Equation (4): where the asterisk (*) denotes the complex conjugate operator; h w is the wave direction vector and D is the difference-frequency quadratic transfer function (QTF) at the pair of frequencies ω j and ω k that, in turn, is assessed by the far-field approach [35].The magnitude of each wave component was determined using Equation (5): where S is the wave spectrum and ∆ω is the wave frequency interval, to be selected less than 2π/T to ensure the randomness of simulation in the time interval T.
The wind loads acting on the emerged part of the tri-floater support structure and wind turbine were determined by using Equation ( 6): where ρ a is the standard air density; A is the structural element projected area normal to the wind direction; C d,w is the drag coefficient; h wind is the wind direction vector and v w is the wind relative velocity, consisting of both steady and turbulent components, the latter generated according to the Norwegian Petroleum Directorate wind spectrum [36].A similar approach was followed to assess the force exerted by the wind turbine, after replacing the drag coefficient C d by the relevant thrust coefficient C t , depending on the wind speed at hub height [37].This approach is certainly slightly approximate, as regards wellestablished and more refined techniques, mainly based on the blade element momentum (BEM) theory [38], but it is accurate enough for the fatigue life assessment of mooring lines, as proved by the benchmark study recently carried out by Barrera et al. [32] on the OC4-DeepCwind platform.
The viscous forces acting on the submerged part of the tri-floater support structure were determined according to the Morison equation and the relative velocity formulation [39] using Equation (7): . q(t)•h w h w (7) where ρ w is the standard water density, u is the wave particle velocity and C d,v is the drag coefficient.The same notation, already endorsed for wind forces, was also applied for the remaining physical quantities.The quasi-static equilibrium model was applied to evaluate the nonlinear mooring forces, considering the seabed friction, the nonlinear restoring components and the line stiffness.Both inertia and damping forces, exerted by the mooring lines, were neglected in this study with no appreciable errors.Really, the validity of this assumption mainly depends on the weight of the mooring system relative to the weight of the floating platform.From a general point of view, if the inertia and damping forces exerted by the stationkeeping system are neglected, some appreciable errors in the assessment of the mooring line forces may arise [13].Anyway, if the weight of the mooring system is low in comparison with that of the floating platform, the quasi-static model does not significantly affect the hydrodynamics of the FOWT [40] and the fatigue life assessment.This outcome was recently stressed by Trubat et al. [41], who investigated the incidence of the mooring hydrodynamics on the assessment of mooring line tension and fatigue for the stationkeeping system of the OC5-DeepCWind platform, verifying that almost the same values were obtained without and with the above-mentioned components.The horizontal, H F , and vertical, V F , forces at the fairlead of each mooring line were determined according to a local reference system, having origin at the anchoring point, according to the nonlinear implicit equation systems (8a) and (8b) for fully and partly developed catenary configurations, respectively: In Equations ( 8a) and (8b), x F (z F ) is the fairlead horizontal (vertical) coordinate, as regards the local reference system; w is the mooring line unit weight in water; L is the line length; E (A c ) is the line Young modulus (cross-section); L b = L − V F /w is the line length lying on the seabed and c b is the bottom drag coefficient.The mooring force equations were preliminarily solved by the Newton-Raphson iteration scheme [42] and subsequently resembled into the global reference system.
The equation system (1) was solved by a dedicated programme developed in MAT-LAB [43], by the ode45 solver, after determining an approximate 10th-order state-space realization of each convolution integral [44] by the Hankel singular-value decomposition algorithm [45].Finally, the frequency-dependent added masses, radiation dampings and normalized wave excitation forces were determined by the open-source code NEMOH [46].

The Combined-Spectrum Approach
Mooring lines are subjected to random cyclic loadings during the entire lifetime that significantly affect the fatigue life of the stationkeeping system.Under irregular cyclic loadings, the fatigue analysis of mooring lines is generally performed according to the Miner rule [47], and the accumulated fatigue damage, d, is determined by Equation ( 9): where n j is the number of cycles, associated with a certain stress range σ j , and N j is the maximum allowable number of cycles, which, in turn, is determined by Equation (10): where K and m are the intercept and slope parameters of the mooring line S-N curve, whose values are listed in the following with reference to stud-link chain cables.Nevertheless, when the stress history in the mooring line is random, as typically occurs in a seaway, the cumulative fatigue damage needs to be replaced by the relevant expected value [19,48].
In this respect, if the stress history is a Gaussian narrow-band process, the stress range follows the Rayleigh distribution, so that the expected fatigue damage, d, is determined by Equation ( 11): where ν j and σ j are the zero up-crossing frequency and the standard deviation of the random j-th stress process, respectively; T j is the exposure time to the j-th sea state condition and Γ is the Gamma function.The previous approach is valid for unimodal stress processes, and so it needs to be properly modified if the stress history is bimodal, as generally occurs for FOWTs.In this case, the low-frequency components, due to the slowly varying drift forces and the turbulent component of wind loads, arise in addition to the wave-frequency stress process, due to first-order wave forces.In this respect, there is a wide consensus that the rainflow method provides the best estimate of the cumulative fatigue damage for bimodal stress processes, but it is very time-consuming [49]; therefore, it is generally replaced by the combined-spectrum approach.Hence, the expected cumulative fatigue damage was assessed by Equation ( 11), after applying the replacements provided by Equation ( 12): where σ j,L (σ j,W ) and ν j,L (ν j,W ) are the standard deviation and the zero up-crossing rate of the j-th low-frequency (wave-frequency) stress process.This approach, which is based on the simple summation of two narrow-band Gaussian stress processes, provides a conservative estimate of the cumulative fatigue damage if bimodal conditions occur.The validity of the assumption about the Gaussianity of the stress process is further discussed in Appendix A. As concerns the second assumption, the fatigue damage was assessed by the combined-spectrum approach, based on Equations ( 11) and (12), assuming that the lowand wave-frequency components of the stress process were narrow-banded.Really, this assumption is on the safe side, as it furnishes a slight conservative estimate of the cumulative fatigue damage in the mooring lines.In addition, it does not affect the main findings of the benchmark study, provided that the cumulative fatigue damage is always obtained by the combined-spectrum approach, after determining the main input parameters by time-history or spectral analysis, as further discussed in Sections 3.2 and 3.3.

Fatigue Assessment by Time-History Analysis
One of the main issues, arising from the application of the combined-spectrum approach, is mainly related to the assessment of the standard deviation and zero up-crossing rates of the low-and wave-frequency components of the stress process.A method, that can be easily implemented to evaluate the above-mentioned quantities, is the signal filtering of the stress process time history.In this respect, the time histories of the low-(wave-) frequency component of the stress process can be detected, after applying an infinite impulse response low-pass (high-pass) filter to the initial random signal [43], with a proper passband frequency and impulse steepness.The passband frequency, f , is generally set depending on the typical periods of the low-and wave-frequency components of the stress process.The impulse steepness, s, instead, is related to the width of the transition region interested by low-pass ∆ f L and high-pass ∆ f H filters, according to Equation (13): where f Ny is the Nyquist frequency.This method is very effective and less time-consuming, but the selection of the passband frequency and impulse steepness is a key factor for the assessment of the low-and wave-frequency components of the stress process.In the current analysis, the passband frequency, f , was set equal to 0.05 Hz and impulse steepness was assumed to be equal to 0.99.The passband frequency corresponded to a period of 20 s, which is a reasonable upper bound for typical waves in open sea conditions.The impulse steepness, instead, was selected high enough to ensure a good separation between the low-and high-frequency contents of the stress signal.The reliability of these assumptions was implicitly checked by the comparative analysis outlined in Section 7.After separately detecting the time histories of the low-and wave-frequency components, the relevant standard deviations were immediately determined.Instead, the zero up-crossing rates, ν j,L and ν j,W , were estimated using Equation ( 14): where T n is the surge natural period of the FOWT and T p is the wave peak period of the first-order wave forces.The effectiveness of these two assumptions was checked by comparative analysis with the advanced spectral analysis methods outlined in Section 4.

Fatigue Assessment by Spectral Analysis
The cumulative fatigue damage in the mooring lines can be also determined by spectral analysis, after determining the bimodal spectrum of the stress process, S σ ( f ), by the advanced spectral analysis techniques outlined in Section 4. Once the combined spectrum was determined and the separation frequency, f , was selected, the spectra of lowand wave-frequency components were determined using Equation ( 15): which implies that the two spectra are non-overlapping and the equality provided by Equation ( 16) holds: The standard deviation of the low-and wave-frequency stress processes was directly determined by the zero-order spectral moments of the low-and wave-frequency spectra, according to Equation (17): while the zero up-crossing rates were assessed using Equation (18), depending on the zeroand second-order spectral moments: This procedure allows the easy detection of all fatigue parameters required to apply the combined-spectrum method, once the spectrum of the combined stress process, S σ , is detected, as discussed in Section 4.

Welch Method
Spectrum estimation methods include two main groups, namely nonparametrical and parametrical [50].The former are based on the application of the Fourier transform (FT) to data, together with proper pre-processing and, eventually, averaging operations.The latter are based on fitting flexible models to the data, such the auto-regressive (AR) or the auto-regressing moving average (ARMA) ones, and then obtaining analytically the corresponding power spectrum.In this paper, two main nonparametrical methods were investigated, namely the Welch's overlapped segment averaging (WOSA) and Thomson multi-taper method (MTM), leaving the examination of the nonparametrical approach to a possible future development [51].
The Welch method consists in parsing the data record, corresponding to an overall observation duration, T, in smaller segments of duration T 0 , with partial overlap, typically from 20% to 50%.Each segment was pre-treated by tapering with a smooth window, to reduce the bias due to spectral leakage, then the "periodogram", i.e., the square of the FT of the series of observations, normalized in respect to the duration T 0 , was calculated for each of them.The spectrum was obtained by averaging over such periodograms.This combination of data tapering and averaging provides a reduction of bias due to spectral leakage and variance related to the randomness of the process [52].
To get some insight in the corresponding estimation procedure, let us denote the series of observations by x i = x(i∆t), where ∆t is the sampling interval, and i = 1, . . .N, with T = N∆t, and T 0 = N 0 ∆t.Let w 1 , . . ., w N 0 be a data taper, then the modified periodogram for the l-th segment is given as follows: where j is the imaginary unit.The spectral estimator is then: where n is the number of segments and m is an integer-valued shift factor, satisfying 0 < m ≤ N 0 and m(n − 1) = N − N 0 .
To optimize the performance of the analyser, two main conflicting features, namely the "effective bandwidth" and the variance (or, equivalently, the standard deviation), have to be controlled to achieve an optimum trade-off.The effective bandwidth is a measure of the minimum separation in frequency between approximately uncorrelated spectral estimates: thus, the wider such bandwidth is, the worse spectral resolution is.For Welch's method, it can be expressed as ∆ f e = α w T −1 0 , where α w is a factor that depends upon the kind of the selected taper and on the way bandwidth is defined.In the case of the Hanning window and considering a half-power bandwidth, we obtain α w = 1.44.On the other hand, standard deviation of the estimator is a measure of its uncertainty, for a given effective bandwidth.For the Welch method, with a 50% overlap, a relative standard uncertainty (standard deviation) u S ( f )/S( f ) = (11/18)N 0 N −1 can be assumed, where u S ( f ) is the absolute standard uncertainty.Once the record duration T is fixed and the kind of taper and the degree of overlap have been decided, the duration of the observation window, T 0 , remains the only design parameter to be optimized, with a trade-off between the need to have a good spectral resolution, which demands a large T 0 , and a small standard uncertainty, which requires a small T 0 .

Thomson Method
This method generalizes the tapering issue by adopting multiple orthogonal tapers to recover the information that may be lost when using a single taper.The estimator is the average of K direct spectral estimators, each acting on the whole data record (rather than on a signal segment, as happens in Welch method) and applying a different taper.Each partial estimator is defined by where h i,k is the kth data taper, usually chosen as the kth discrete prolate spheroidal sequence with parameter W, where 2W is the normalized bandwidth of the tapers, i.e., the bandwidth for ∆t = 1 s.The final estimator is thus as follows: where K is typically chosen to be equal to 2NW − 1 [53].Again, a trade-off is required between effective bandwidth and relative standard uncertainty.The effective bandwidth can be estimated as ∆ f e = 2W/∆t (Hz).Furthermore, considering that the distribution of the estimator can be approximated by S( f )χ 2 2K /2K, the relative standard uncertainty was equal to K − 1 2 .There is here much less arbitrariness than in Welch's method, since, for a fixed observation time, T, the only parameter to be chosen was the half-bandwidth W, which influences both spectral resolution and relative standard uncertainty.

The OC4-DeepCwind Platform
The OC4-DeepCwind platform [29] consists of a tri-floater support structure, equipped with the NREL 5 MW wind turbine [30].It has three offset columns, connected to each other by two sets of bracings and consisting of an upper and a base cylinder, the latter with an enlarged diameter to suppress, as far as possible, the heave motion of the FOWT.The tri-floater is also equipped with a central main column that supports the tower of the wind turbine and is connected to the offset columns by an additional set of bracings.The main properties of the FOWT are listed in Table 1, while the thrust curve of the wind turbine is provided in Figure 2, as a function of the wind speed at hub height, for both power production and parked wind turbine conditions [37,54].In this respect, it must be pointed out that the rotor control mechanism of the wind turbine was not considered here.Anyway, it is expected that this assumption does not significantly affect the analysis performed in Section 6.In fact, Barrera et al. [32] recently carried out a comparative analysis on the stress histories in the mooring lines of the OC4-DeepCWind platform, modelling the wind forces by a look-up table of thrust coefficients and by the blade element momentum theory, including, in the latter case, the rotor control mechanism of the wind turbine.They verified that the agreement between the different numerical approaches and the experimental results was always very good.The response amplitude operators (RAOs) due to first-order wave loads and the second-order difference-frequency quadratic transfer functions (QTFs) of the OC4-DeepCWind platform are listed in Appendix B.

The Stationkeeping System
The FOWT is equipped with a three-line spread mooring system, consisting of studlink chain cables.The main properties of the stationkeeping system and the soil conditions are listed in Table 2, while a sketch of the mooring layout is depicted in Figure 3.

The Stationkeeping System
The FOWT is equipped with a three-line spread mooring system, consisting of studlink chain cables.The main properties of the stationkeeping system and the soil conditions are listed in Table 2, while a sketch of the mooring layout is depicted in Figure 3.

Selection of Reference Conditions
The benchmark study, carried out in Section 6, was performed based on the reference environmental conditions recently investigated by Xu et al. [36] and listed in Table 3.They are representative of both power production and parked wind turbine situations.The sea state condition is described by the JONSWAP spectrum [39], with a significant wave height,  , ranging from 1.96 to 9.77 m and a wave peak period ranging from 9.72 to 12.95

Selection of Reference Conditions
The benchmark study, carried out in Section 6, was performed based on the reference environmental conditions recently investigated by Xu et al. [36] and listed in Table 3.They are representative of both power production and parked wind turbine situations.The sea state condition is described by the JONSWAP spectrum [39], with a significant wave height, H s , ranging from 1.96 to 9.77 m and a wave peak period ranging from 9.72 to 12.95 s.The wind climate was described by the 10 min wind speed at hub height, combined with the Norwegian Petroleum Directorate wind spectrum [36], which is representative of the turbulent component.Wave and wind loads were assumed to be codirectional, with heading angle equal to 0 • , to exert the highest loads in the mooring line #1.

Preliminary Analysis
Before carrying out the benchmark study, a preliminary analysis was performed to test the effectiveness and robustness of the hydrodynamic model developed in MATLAB.Figures 4 and 5 provide the mean values, µ, and the standard deviations, σ, of the stress process in the mooring line #1, obtained by the numerical model, currently developed in MATLAB, and by the numerical simulations performed by Xu et al. [36] by the commercial code SIMA [55].The SIMA software, currently licensed and supported by the Det Norske Veritas, is a complete tool for the simulation of marine operations from modelling to results, which include lifting of topsides and modules, offshore wind turbine installation, offshore crane operations, subsea installations, transportation and installation of TLPs and towing operations, among others.The software, widely tested and validated in the past, was recently extended to FOWTs, and it can be considered a reliable and effective mean of comparison.By the analysis of current results, it was gathered that the numerical code, developed in MATLAB allows the efficient estimation of the mean value and the standard deviation of the stress process, with errors less than 5%.
tion, offshore crane operations, subsea installations, transportation and installation o TLPs and towing operations, among others.The software, widely tested and validated in the past, was recently extended to FOWTs, and it can be considered a reliable and effectiv mean of comparison.By the analysis of current results, it was gathered that the numerica code, developed in MATLAB allows the efficient estimation of the mean value and th standard deviation of the stress process, with errors less than 5%.tion, offshore crane operations, subsea installations, transportation and installation o TLPs and towing operations, among others.The software, widely tested and validated i the past, was recently extended to FOWTs, and it can be considered a reliable and effectiv mean of comparison.By the analysis of current results, it was gathered that the numerica code, developed in MATLAB allows the efficient estimation of the mean value and th standard deviation of the stress process, with errors less than 5%.

Benchmark Study 6.1. Assessment of the Stress Process in the Mooring Lines
The stress process in mooring line#1, corresponding to a typical 3 h storm duration, is plotted in Figure 6a-h, with a sampling frequency of 5 Hz.As can be gathered from the relevant stress histories, the standard deviation of the stress process increased along with the sea state severity.In this respect, the load conditions from LC1 to LC6, which refer to the power production mode, were clearly dominated by the low-frequency component of the stress process, while the low-and wave-frequency components of the stress process were comparable in the last two conditions, which refer to the parked wind turbine situation.These outcomes will be further investigated by the fatigue analysis, performed in Sections 6.2 and 6.3, by the time-domain and spectral models, respectively.

Time-Domain Analysis of the Stress Process
Table 4 provides the fatigue assessment of mooring line #1 by time-domain analysis, considering a separation frequency equal to 0.05 Hz.The cumulative fatigue damage per unit time is reported for each load condition, together with the standard deviations and up-crossing rates of the low-and wave-frequency components of the stress process.In this respect, it must be pointed out that the up-crossing rate of the low-frequency component was constant according to the first of Equation ( 14), as it depends on the surge natural period of the FOWT, equal to 116 s.The effectiveness of this assumption is checked in Section 7. (h) load condition LC8.

Time-Domain Analysis of the Stress Process
Table 4 provides the fatigue assessment of mooring line #1 by time-domain analysis, considering a separation frequency equal to 0.05 Hz.The cumulative fatigue damage per unit time is reported for each load condition, together with the standard deviations and up-crossing rates of the low-and wave-frequency components of the stress process.In this respect, it must be pointed out that the up-crossing rate of the low-frequency component was constant according to the of Equation ( 14), as it depends on the surge natural period of the FOWT, equal to 116 s.The effectiveness of this assumption is checked in Section 7.

Spectral Analysis of the Stress Process
The spectra of the stress process in mooring line #1, determined by Welch and Thomson methods, are plotted in Figure 7a-h using black continuous and red dashed lines, respectively.In all cases, the separation frequency was set equal to 0.05 Hz.In addition, the low-and wave-frequency components of the stress process were clearly detectable.The former is due to the combined action of the slowly varying drift forces and the turbulent component of the wind forces; the latter, instead, is related to the first-order wave forces.The cumulative fatigue damage per unit time and the relevant fatigue parameters required by the combined spectrum method are reported in Tables 5 and 6.By a comparative analysis between the cumulative fatigue damage provided by the Welch and Thomson methods, the two spectral analysis techniques provide comparable results, in terms of standard deviations and up-crossing rates both low-and wave-frequency components of the stress process in the mooring line #1.Further details are provided in Section 7.

Discussion
Figure 8a-d provides a comparative study between the time-domain analysis and the advanced spectral reconstruction models, focusing on the estimated values of the standard deviations and up-crossing rates of the low-and wave-frequency components.As concerns the stress process standard deviation, a very good agreement was reached in all cases.Some differences between time-and frequency-domain analysis arose with reference to the up-crossing rate of the low-frequency component of the stress process.Anyway, these differences had a negligible impact on the cumulative fatigue damage per unit time reported in Figure 9.In fact, a substantial agreement was in all cases recognized between time-domain and spectral models, with absolute errors generally less than 5%.This outcome was partly due to the endorsement of the combined spectrum approach in both cases, but it also highlighted that both simplified time-domain and advanced spectral analysis can be applied to assess the cumulative fatigue damage in the mooring lines, with almost the same level of accuracy.
Based on current results, some main outcomes can be stressed: (i) Time-domain models are accurate enough, at least when a variety of cases need to be analysed to assess the design fatigue life of the mooring lines.In this respect, the up-crossing rate of the low-and wave-frequency components of the stress process can be assumed to be equal to the surge/sway natural period of the tri-floater platform and the wave peak period of the relevant sea state condition.These assumptions are accurate enough, as can be gathered from Figure 8c,d.(ii) Welch and Thomson methods allow accurate detection of the low-and wave-frequency components of the stress process once the partitioning frequency is provided.It was verified that the separation frequency of 0.05 Hz allowed accurate detection of the low-and wave-frequency components.(iii) It is quite challenging to identify which spectral analysis method, between the Welch and Thomson models, is the most suitable one.Presumably both methods can be applied with the same confidence, at least based on current results, provided they lead to almost the same cumulative fatigue damage per unit time in all the analysed loading conditions.Based on current results, some main outcomes can be stressed: (i) Time-domain models are accurate enough, at least when a variety of cases need to be analysed to assess the design fatigue life of the mooring lines.In this respect, the upcrossing rate of the low-and wave-frequency components of the stress process can be assumed to be equal to the surge/sway natural period of the tri-floater platform and the wave peak period of the relevant sea state condition.These assumptions are accurate enough, as can be gathered from Figure 8c,d

Conclusions
The paper focused on the assessment of the cumulative fatigue damage in the mooring lines of an FOWT, deployed on intermediate water depth, by time-domain and advanced spectral analysis methods.After a brief review about the state-of-art, the nonlinear

Conclusions
The paper focused on the assessment of the cumulative fatigue damage in the mooring lines of an FOWT, deployed on intermediate water depth, by time-domain and advanced spectral analysis methods.After a brief review about the state-of-art, the nonlinear timedomain hydrodynamic model, was provided and the various force components were separately discussed.The cumulative fatigue damage in the mooring line #1 was assessed by the combined spectrum approach, accounting for both low-and wave-frequency components of the stress process, as typically occurs for offshore structures.The fatigue analysis was carried out by time-domain analysis and advanced spectral reconstruction models, based on Welch and Thomson methods.
In order to test the effectiveness of both time-and frequency-domain models, a benchmark study was performed, assuming as the reference FOWT the OC4-DeepCwind trifloater platform, equipped with the NREL 5 MW wind turbine.Several loading conditions, representative of both power production and parked wind turbine situations, were considered to investigate the incidence of met-ocean conditions on the assessment of the fatigue cumulative damage.As stressed in Section 8, a very good agreement between time-and frequency-domain models, was in all cases recognized, confirming that the former is accurate enough when a variety of cases need to be analysed, while the latter provides a more rational approach for the assessment of the main fatigue parameters.In this respect, the Welch and Thomson methods were revealed to be very promising techniques for the reliable assessment of the mooring fatigue life, as they allow the proper estimation of all the fatigue parameters, required by the combined stress method.The current results are encouraging for future research, devoted to endorsing spectral analysis in the assessment of the mooring line fatigue life, in view of a more reliable and rational design approach of stationkeeping systems, and eventually to update the mooring-related safety factors, embodied in current rules and guidelines.

Appendix B
Figure A1 depicts the absolute values of the complex exciting forces per unit wave amplitude at 0 • , which was the heading angle considered in the benchmark study.Particularly, Figure A1a,b refers to the surge/heave and pitch motions, respectively.All the remaining forces are null, due to the symmetry condition as regards the x-axis of the global reference system.The frequency-dependent added masses and radiation dampings are provided in Figure A2.Particularly, Figure A2a provides the surge/sway and heave added mass, while Figure A2b depicts the relevant radiation dampings.Similarly, Figure A2c,d provides the same hydrodynamic properties for the roll/pitch and yaw motions.Finally, Figure A2e,f reports the surge/pitch and sway/roll added masses and radiation dampings.
In all cases, a very good agreement with the relevant values provided by Robertson et al. [29] was recognized.As concerns the difference-frequency QTF for the surge motion of the tri-floater platform, reference is made to Chuang et al. [57].
the stress process [56]: where  ( ) is the kurtosis (skewness) of the stress process and  is the slope of the S-N curve.The high-order statistics of the eight stress processes investigated in Section 6 are reported in Table A1, from which it is gathered that the non-Gaussianity of the stress process is generally quite low.

Figure 1 .
Figure 1.Reference layout of the OC4-DeepCwind platform (the symbol # means line number).

Figure 1 .
Figure 1.Reference layout of the OC4-DeepCwind platform (the symbol # means line number).

Figure 2 .
Figure 2. Thrust curve of the NREL 5 MW wind turbine.

Figure 2 .
Figure 2. Thrust curve of the NREL 5 MW wind turbine.

Figure 3 .
Figure 3. Layout of the reference stationkeeping system (the symbol # means line number).

Figure 4 .
Figure 4. Comparative analysis-mean value of stress history in line #1.

Figure 5 .
Figure 5. Comparative analysis-standard deviation of stress history in line #1.

Figure 4 .
Figure 4. Comparative analysis-mean value of stress history in line #1.

Figure 4 .
Figure 4. Comparative analysis-mean value of stress history in line #1.

Figure 5 .
Figure 5. Comparative analysis-standard deviation of stress history in line #1.

Figure 5 .
Figure 5. Comparative analysis-standard deviation of stress history in line #1.

Figure 8 .
Figure 8. Comparative analysis between time-and frequency-domain models.(a) Standard deviation of the low-frequency stress component; (b) standard deviation of the wave-frequency stress component; (c) up-crossing rate of the low-frequency stress component; (d) up-crossing rate of the wave-frequency stress component.

Figure 8 .
Figure 8. Comparative analysis between time-and frequency-domain models.(a) Standard deviation of the low-frequency stress component; (b) standard deviation of the wave-frequency stress component; (c) up-crossing rate of the low-frequency stress component; (d) up-crossing rate of the wave-frequency stress component.J. Mar.Sci.Eng.2022, 10, x FOR PEER REVIEW 19 of 24

Figure 9 .
Figure 9. Assessment of cumulative fatigue damage per unit time.

Figure 9 .
Figure 9. Assessment of cumulative fatigue damage per unit time.

Figure
FigureA1depicts the absolute values of the complex exciting forces per unit wave amplitude at 0°, which was the heading angle considered in the benchmark study.Particularly, FigureA1a,b refers to the surge/heave and pitch motions, respectively.All the remaining forces are null, due to the symmetry condition as regards the x-axis of the global reference system.The frequency-dependent added masses and radiation dampings are provided in FigureA2.Particularly, FigureA2aprovides the surge/sway and heave added mass, while FigureA2bdepicts the relevant radiation dampings.Similarly, FigureA2c,d provides the same hydrodynamic properties for the roll/pitch and yaw motions.Finally, FigureA2e,f reports the surge/pitch and sway/roll added masses and radiation dampings.

Table 2 .
Main properties of the reference stationkeeping system.

Table 3 .
Main data of reference conditions.

Table 4 .
Fatigue assessment by time-domain analysis.

Table 5 .
Fatigue assessment by spectral analysis-Welch method.

Table 6 .
Fatigue assessment by spectral analysis-Thomson method.

Table A1 .
High-order statistics of the eight stress processes.

Table A1 .
High-order statistics of the eight stress processes.