Lifetime assessment of semi-submersible wind turbines by Gaidai risk evaluation method

As the global agenda turns more towards the so-called challenge of climate change and lowering carbon emissions, research into green, renewable energy sources becoming nowadays more and more popular. Offshore wind power, produced by FOWTs (i


Introduction
One of the most important ecologically friendly renewable green energy sources, that might meet the world's growing energy demands is wind energy.To capture wind energy and produce electricity, offshore WT farms are being built on a regular basis.These wind farms are usually constructed on the continental shelf.Since windspeeds offshore being often higher, than those on land, the offshore wind power sector plays a major industrial role in the electricity supply.Greater in situ windspeeds are often the result of a relatively smooth sea surface.The effective production of renewable energy is dependent on the latest enhancements to the FOWT's modern design.Due to FOWTs' innate vulnerability to ambient wind turbulence and nonlinear aerodynamic and hydrodynamic loads, their extreme dynamics characteristics should be accounted for in both initial design and operations.
There are typically 2 distinct approaches for calculating the design loads of WTs (i.e., Wind Turbine): (a) simulating rare, high-load-level occurrence, e.g., 10,000-year wave; (b) simulating WT reaction under usual operating conditions, then extrapolating joint PDF (i.e., Probability Density Function) of WT failure probability (International Energy Agency, 2020; Poujol et al. 2020;Ferraz de Paula and Carmo 2022;Maienza et al. 2022).For life cycle assessment under environmental impacts for FOWTs in deep-sea areas, see (Brussa et al. 2023).IEC (i.e., International Electrotechnical Commission) 61,400-1 standard, (GWEC 2021) recommends using both techniques to assess characteristic design loads.FOWT in offshore operations has been illustrated in Fig. 1.This study contributes to the 2nd approach (b), (Chen et al. 2022) , by advocating state-of-art methodology, which has already been shown successful for a wide range of naval and marine structures, including FOWTs, offshore platforms, and vessels.The internal bending moments and anchor tensions of the FOWT were numerically calculated to evaluate the reaction of the structural dynamic system.NREL (i.e., National Renewable Energy Laboratory, USA) created the FAST code, a potent nonlinear aero-hydroservo-elastic modelling toolbox.Modal dynamics is incorporated into the FAST software package to simulate multibody structural dynamics, (Jonkman and Buhl 2005;NOAA 2023;IEC 2009;Haid et al. 2013;Jain et al. 2012).To examine critical FOWT's bending moments and anchor stresses, load data from 2550× 10-min FAST simulations, which spanned about 18 days in total, were assessed.These simulations represented a range of potential in situ meteorological circumstances.Since the 10-min single simulation duration mandated by current regulations for both land-based and offshore-based turbines may not be adequate for combined wind and wave loads on FOWT, extensive numerical Monte Carlo (MC) simulations have been carried out.
Developing more accurate tools to assess critical structural stresses that FOWTs could encounter when operating in choppy waters was the aim of numerous recent studies.(Ronold et al. 1999;Ronold and Larsen 2000;Manuel et al. 2001;Fitzwater and Cornell 2002;Moriarty et al. 2002;Agarwal and Manuel 2008;Barreto et al. 2022;McCluskey et al. 2021).In (Fogle et al. 2008), authors have utilized block maxima extrapolation, to assess critical structural loads.In (Ernst and Seume 2012) authors have accounted for ambient turbulence intensity, while performing a WT risks evaluation study.MC method has been widely utilized by various authors in (Freudenreich and Argyriadis 2007;Gaidai 2022;Peeringa 2009;Abdallah 2015) to assess long-term FOWT fatigue loadings.For predictions of design WT loads, based on limited datasets, see e.g., (Stewart et al. 2015).The techniques described above fit an Extreme Value Distribution, or EVD, to the underlying dataset; nevertheless, Fig. 1 Example of the Deep-C-wind semi-submersible-type FOWT, (Xu et al. 2022) this may not meet the asymptotic condition.The latter becomes a crucial matter since it raises the possibility of less conservative outcomes when evaluating design/ characteristic values.To allow improved FOWT design and its control system, it is imperative to examine nonlinear extreme structural dynamics using more accurate and efficient risks evaluation models and numerical approaches.With improved WT design, extreme climatic loadings will be less likely to cause mechanical damage to FOWTs; precise design load predictions would also enable more efficient sizing of important FOWT components, (Gaidai et al. 2022).
Contemporary risks evaluation approaches that are currently in use in offshore and wind engineering, suggest that modern structural risks evaluation methods are not well-suited to handle high-dimensional dynamic systems, particularly when there are intricate non-linear cross-correlations between the FOWT system's essential components.Primary purpose and objective of the current study was to fill the later research gap, in particular what regards systems high-dimensionality.The operational lifespan evaluation of complex dynamic systems, which is crucial for the robust and long-lasting design of a variety of offshore structures exposed to in situ environmental pressures, is carried out in this study by employing a novel Gaidai risks evaluation methodology.For recent studies on WT safety and reliability, see (Cuesta 2024;Sharar and Hoover 2024;Kenworthy, et al. 2024).In the following Sect. 2 to introduce WT system, Sect. 3 to outline novel reliability methodology, Sect. 4 to report results and numerical findings.

System's description
Obtaining excellent, high-resolution met-ocean data may be challenging.The NOAA (i.e., National Oceanic and Atmospheric Administration, USA) provided the dataset that was utilized in this investigation.Extensive network of floating data collecting buoys being kept up to date by NOAA in both domestic and foreign waters.The NOAA Data Buoy Center offers online access to data from these buoys, (NOAA 2023).Wind directions and observed windspeeds for the specific dataset under analysis in this study have been averaged across 8-min intervals.Significant wave heights may be identified based on 20-min measurements periods.Throughout each 20-min wave measurement interval, the peak-spectral period, along with the wave period with the highest wave energy had been assessed.The measuring buoy station at Cape Elizabeth, located 45 nautical miles northwest of Aberdeen, state of Washington, close to the US continental shelf line, has been selected for this investigation.Combined wind-wave statistics for the indicated site have been produced using hourly historical met-ocean in situ measurements from2010 to 2017.For an MC-based statistical long-term analysis flowchart see Fig. 2. Authors will further refer to all scatter diagram's environmental in situ sea/ocean conditions, including local windspeeds, as "sea states".
Windspeeds should then be extrapolated to a standard FOWT's hub height of 90 m.Anemometers have been installed at Cape Elizabeth buoy, five meters above MSL (i.e., Mean Sea Level).PL (i.e., Power Law) wind shear equation to be utilized with α being PL constant, U (z) , U (z r ) being windspeed at height z and reference windspeed at elevation z r , respec- tively.In this study, PL constant, specified by Eq. (2) has been fixed to α = 0.14, (IEC 2009).The following nota- tions have been employed.
Windspeed U Significant wave height H s Peak spectral wave period T p In this work, an empirical multi-dimensional joint PDF was formed from the observed buoy dataset through post-processing using only in-situ scatter diagrams, without any assumptions or simplifications.Joint 3D (i.e., 3-Dimensional) PDF p U , H s T p in-situ scatter diagrams have been straightforwardly derived from underlying (1) in-situ met-ocean dataset.In the FAST numerical simulation, 3D PDF U , H s T p has been used as an input; wind/wave misalignments have not been accounted for.
Gaidai risks evaluation approach being particularly efficient for high-dimensional PDFs.In the next, predictions, obtained from a direct MC-based long-term numerical simulation approach, will be compared with those by the Gaidai risks evaluation method.The direct MC-based approach is schematically displayed in Fig. 2. The OC4 semi-submersible Deep-C-wind floating dynamic system, supporting the FOWT platform, has been selected in this study.Three external offset columns and one central main column make up the semi-submersible type FOWT platform.The FOWT bottom has heave plates affixed to it to lessen excessive heave motions.Table 1 displays the principal dimensions of the semi-submersible type FOWT platform under examination.NREL 5-MW baseline FOWT had been mounted on top of the Deep-C-wind semi-submersible-type FOWT platform.Hub's height of the FOWT cylindrical tower was 90 m, 126 m for the FOWT 3-bladed rotor.The 5-MW FOWT baseline's summary properties have been presented in Table 2.
FOWT structural dynamics, as well as applied aerodynamic and gravitational stresses, have been modelled using FAST and Aero-Dyn.The FOWT dynamics encompass tower elasticity and dynamic interaction between the FOWT platform and tower movements.FAST utilizes a numerical model that combines modal dynamics along with structural dynamics for multibody dynamics, (Jonkman and Buhl 2005).
All relevant/key FOWT system's DOFs (Degrees Of Freedom) have been well accounted for in FAST numerical simulation, including 2 flap-wise and 1 edge-wise moment, 1 torsion DOF for FOWT powertrain, 1 variable generator-speed DOF for FOWT nacelle.The FOWT platform had 3 DOFs for translation (surge, sway, heave), 3 for rotation (roll, pitch, yaw), 2 fore-aft and 2 side-toside FOWT tower modes DOFs.To numerically simulate the FOWT dynamic system's cross-correlated key/critical components, the aero-hydro-servo-elastic simulation software package FAST, (Jonkman and Buhl 2005) has been utilized.Turb-Sim software has been utilised to generate stochastic wind speed fields on a 145 m wide, 31× 31 square grid.When modelling aerodynamics of baseline WTs, the Aero-Dyn module of the FAST code, which utilizes the BEM (i.e., Blade Element Momentum) technique, in order to describe FOWT aerodynamics, takes into consideration FOWT rotor wake effects, in addition to dynamic stall.To model FOWT hydrodynamic loads, the FAST's Hydro-Dyn module had been utilized, combining potential flow equations for large diameter structure, given Morison's hydrodynamic force equation F D = 1 2 C D ρAv 2 for offshore slender structure, with C D being drag coefficient, ρ being water density, A tubular cross-sectional area, vbeing water particle veloc- ity.Morison's equation effectively takes into account the viscous drag external loads operating on FOWT via the drag force component.FAST numerical models have done a good job of accounting for second-order hydrodynamic forces.Standard IEC-61400-1 states that to assess ultimate loads over 50 years, based on standard IEC-61400-3, at least 15 short-term MC simulations, 10 min duration each, are required, (IEC 2009).Note the latter 10-min MC sample's duration being typical for FOWT aerodynamic modelling, as opposed to 3-h direct MC runs for the hydrodynamic type of problems, involving 3-h stationary sea/ocean states.Figure 3 presents a timeseries example for the FOWT tower's fore-aft bending moment, affecting the FOWT tower's base.
Extensive experimental effort of FOWT modelling methodologies being required to give experimental data that is essential in verifying the completion of OC4, OC5, and similar projects.The OC5 FAST numerical data, the experimental/lab findings, and the numerical results all showed an average underprediction of the FOWT tower top ultimate shear loads, tower-base loads, and upwind mooring/hawser tensions of around 10%, 14%, 20%,  respectively, (Robertson et al. 2017).The authors relied on FAST as an accurate numerical simulation tool, based on the latter experimental validation findings.

Gaidai risks evaluation method
LTD (i.e., Life Time Distribution) estimation approach for complex, nonstationary, nonlinear FOWT dynamic systems, subjected to various risk/hazard modes over specified service timelapse, is briefly discussed in this section.FOWT dynamic system successive lifetimes L i , i = 1,2, . . ., measured in-between consecutive failure/ hazard/risk/damage threshold crossings, generating sequential L values.The lifetime of the FOWT dynamic system is thus represented by L , which is a random vari- able.LTD of complex nonstationary FOWT systems cannot be easily calculated, using current technical risks evaluation methodologies, particularly when there are many important dimensions and components (failure/hazard modes) in the system.It is now possible to evaluate the target CDF (i.e., Cumulative Distribution Function) using extensive direct type MC simulations, or sufficient measurement data, if available for the FOWT dynamic system, (Gaidai et al. 2022;Gaidai et al. 2022;Gaidai et al. 2022;Gaidai 2022;Gaidai et al. 2022;Gaidai and Xing 2022;Gaidai et al. 2022).The expense of computational resources and lab testing may frequently be prohibitive (2) LDT(L) ≡ Prob(Lifetime ≤ L) for the majority of complex dynamic energy/environmental systems.The authors of this study advocate a novel LTD assessment methodology for FOWT systems that lowers the cost of computation and measurement during the FOWT system's design phase, Fig. 4.
Consider the vector (X(t), Y (t), Z(t), . . . ) , which rep- resents the MDOF (Multi-Degree Of Freedom) system's overall dynamic that has been either simulated or measured within a representative timelapse (0, T ) .This vector is made up of the key/critical components/dimensions of the FOWT system, which are inter-correlated nonlinear FOWT system's critical components (loads or responses The 1D (1-Dimensional) critical/key global maxima of the FOWT system are represented by the notationsX max T = X(t), Y max T = Y (t) , Z max T = Z(t), . . .The authors define "representative" timelapse (0, T ) as a big enough value of T in relation to the autocorrelation of the dynamic FOWT system and the relaxation temporal scales.
The crucial process X = X(t) a dynamic system has local maxima X 1 , . . ., X N X that happen at discrete, tempo- rally increasing instants t X 1 < • • • < t X N X inside (0, T ) .For additional MDOF crucial components, such as Y (t), Z(t) , and others, specifically Y 1 , . . ., Y N Y ;Z 1 , . . ., Z N Z similar definitions must be used.For the purpose of simplicity, all local maxima for the critical elements of the FOWT system have been assumed to be positive.
Next, Fig. 3 FOWT tower base's fore-aft bending moment timeseries, U hub = 11m/s, (Xu et al. 2022) being the target probability of the dynamic system survival, andη X ,η Y ,η Z ,… being the critical/hazard/ limit values of the FOWT system's component parts; p X max T ,Y max T ,Z max T ,... being the system's joint PDF of each global maxima of the system's component parts.It could not be practically possible to predict directly the joint PDFp X max T ,Y max T ,Z max T ,... , and therefore the FOWT system's goal survival chances/probability P if the dynamic FOWT system's NDOF (number of degrees of freedom) is enormous.Based on the underlying system's lifespan CDF distributionLTD(L) , provided by Eq. ( 1), it is simple to evaluate the hazard/failure risk/probability P failure of the FOWT system as well as the associated system's expected lifetime L expected having the likelihood, risk, and failure of the sys-temP failure = 1 − P , which is supplementary to the survival probability of the FOWT systemP .If any of the 1D (unidimensional) critical components of the dynamic FWT system X(t) exceedingη X , Y (t) sur- passingη Y , Z(t) exceedingη Z , etc.-are exceeded, the system is said to have failed (or been damaged).Individually defined fixed failure, danger, and damage levels ( X, Y , Z, . . .etc.) for each system's 1D critical and cru- cial componentsX max N X = {X j ; j = 1, . . ., 1D key components of a dynamic system, X, Y, Z,… being rescaled, and then non-dimensionalized resulting in the crucial components of the dynamic FOWT system being non-dimensional, with target failure probability P = P(1) and target hazard/failure/ risk/damage limits = 1 .P( ) is a non-dimensional level smooth C 1 function.The 1D critical compo- nents local maxima of the dynamic FOWT system are merged into a single temporally increasing vector R(t) ≡ − → R = (R 1 , R 2 , . . ., R N ) in accordance with the cor- responding temporal vector of the merged FOWT system t 1 < • • • < t N , N ≤ N X + N Y + N Z + . . . .Every local maxima of R j came into contact with the local maxima of the dynamic FOWT system critical component, which corresponded X(t) or Y (t) , or Z(t) , or other dynamic sys- tem components, (Gaidai et al. 2022;Gaidai et al. 2022;Gaidai et al. 2022).For constructed synthetic − → R -vector thus having no data loss, see Fig. 5.In simple words, the − → R vector has been obtained by the temporal coalescing (5) In the case of clustering of neighbouring local maxima, the de-clustering procedure has to be done.The non-exceedance (system's survival) probability P( ) must be reevaluated for k < j ≤ N (level of conditioning k ) in order to eliminate clustering between the local maxima of the temporally neighbouring system's components.The approximation that was first proposed by Eq. ( 6) may now be further stated as The objective is to maintain a record of every individual damage, danger, and failure that happened locally and originally in time.This will lessen the likelihood of intercorrelated local exceedances and cascade (i.e., de-clustering).Given that the MDOF (6 5 Example of how 2 components, X, Y being merged into 1 new synthetic vector − → R .Red ellipse highlights the case of simultaneous maxima for different system's components.Nondimensional parameter on the vertical axis FOWT system is thought to be piecewise ergodic, the chance of the system failing or sustaining damage p k ( ) := Prob{R j > η j |R j−1 ≤ η j−1 , R j−k+1 ≤ η j−k+1 } for j ≥ k is independent of j and only dependent on the degree of conditioning k .Now, the likelihood of non-exceedance may be evaluated as Note that Eq. ( 7) follows from Eq. ( 6) if neglecting Prob(R 1 ≤ η 1 , . . ., R k ≤ η k ) ≈ 1 , as design failure/dam- age/hazard probability being of small order of magnitude, with N ≫ k , (Gaidai et al. 2022;Gaidai et al. 2022;Gaidai et al. 2022).Regarding the conditioning parameterk , convergence is required The synthetic − → R -vector of the built system exhibits a net 0 data loss, (Gaidai et al. 2023;Gaidai et al. 2023;Gaidai et al. 2023;Zhang et al. 2023;Liu et al. 2023;Gaidai et al. 2023).Now that the synthetic vector − → R has been presented, the stochastic longevity variable L may be defined along with the corresponding temporally increasing occurrence temporal instants t 1 < • • • < t N where i = 2, . . ., N .The integrated synthetic process R(t) the FOWT system thus contains crucial informa- tion about the target LTD PDF.Currently, target FOWT survival chances/probability can be evaluated using the MUR (i.e., Mean Up-crossing Rate) function P = P(1) with ν + ( ) denoting the MUR function of the hazard/ failure/risk/damage threshold level for the non-dimensional synthetic vector R(t) of the FOWT system as pre- viously discussed, (Liu et al. 2023;Gaidai 2023;Gaidai et al. 2023;Yakimov et al. 2023;Gaidai et al. 2023;Sun et al. 2023;Gaidai et al. 2023;Yakimov et al. 2023).For the MUR function, Rice's formula is represented by Eq. ( 5), where the joint PDF corresponding to R, Ṙ is represented by p R Ṙ and the time-derivative is given by Ṙ = R′(t) .The dynamic FOWT system approaches the pre-established threshold of → 1 for failure, risk, haz- ard or damage (7) The MDOF dynamic FOWT system had been assumed to be jointly quasi-stationary.A novel short-term environmental, sea, or ocean state is given an in-situ environmental scatter diagram with m = 1, .., M environmental/ sea/ocean states, each with its own occurrence probability/chances probability/chances q m , with M m=1 q m = 1 , (Gaidai et al. 2023;Gaidai et al. 2023;Gaidai et al. 2023;Gaidai et al. 2023;Sun et al. 2023;Gaidai et al. 2024;Gaidai et al. 2024;Gaidai et al. 2024).Currently, a longterm probability equation indicates that Applying analogous functions to Eq. ( 6) for p k ( , m) , which represent each unique short-term ambient sea/ocean condition/state with number m.The following Section will show how to evaluate the relevant LTD q quantiles.using underlying time series that are either MC simulated or measured, where ° denotes functional composition, q ∈ (0, 1) and LTD −1 are the inverse of LTD, and LTD°LTD −1 = 1 is the identity-operator.The target lifes- pan distribution will resemble Poisson's distribution with the parameter ν + ( ) Tsince the failure, hazard, or dam- age events of the FOWT dynamic system become almost independent at high failure, hazard, risk, or damage levels.Current study primarily focused on the Poisson-type process, but it may have broader implications if any failure, hazard, risk, or damage of any of the critical components or dimensions of the FOWT system does not indicate that the system is about to collapse, but instead merely consists of a sequence of interdependent local maxima of the FOWT system's critical components, (Gaidai 2024;Gaidai 2024;Gaidai et al. 2024;Gaidai et al. 2024) (Yuan 2023;Gaidai 2022;Gaidai. 2023;Gaidai 2023)

Results
In this section, the efficacy of the previously mentioned Gaidai risks evaluation methodology is illustrated through an analysis of three dynamic key/critical components of the FOWT system using the Gaidai risks evaluation method: 1) the FOWT tower's fore-aft bending moment; 2) the FOWT blade's root's out-of-plane bending moment; and 3) the FOWT anchor's tension.( 12) FOWT dynamics is a difficult dynamic system to assess due to its cross-correlated, multidimensional, highly nonlinear characteristics.Moreover, system risks evaluation research is crucial for FOWTs working in extreme weather conditions.This section will demonstrate the Gaidai risks evaluation technique in practice.3 FOWT system's critical components X, Y , Z(internal forces in our study's case), defined above, had been selected as an engineering example of a 3D dynamic FOWT system.For each FOWT system's critical component, its failure/hazard/risk limit had been set equal to its recorded global maximum, but increased twice; note that this being a purely subjective choice, done only for an illustration, (Sun et al. 2023;Gaidai et al. 2024;Gaidai et al. 2024;Gaidai et al. 2024;Gaidai 2024;Gaidai 2024;Gaidai et al. 2024.The three recorded time series' local maxima were then integrated into a single synthetic time series, with the rising periods of occurrences of the local maxima of the three important components of these systems being utilized to sort each set X j , Y j , Z j .As a result, the criti- cal dimensions and components of all three systems become non-dimensional, with the same = 1 hazard, risk, failure, and risk limitations.The non-dimensional constructed synthetic vector − → R in (Fig. 6a) is derived from assembled local maxima of FOWT critical components (in this case, internal forces), with lower values ≥ 0 being irrelevant for target hazard/failure/dam- age PDF tail's extrapolation towards the target threshold/level = 1 .The > 0.3 cut-on limit has been used purely illustratively.Since it has been made up of various FOWT critical components, having distinct measurement units-in our example, Nm for bending moments and kN for forces, the FOWT system's vector − → R does not possess any particular physical meaning.Figure 6a) presents non-dimensional system's synthetic vector − → R .Fig- ure 6b) presents extrapolation towards a 1-year return period of interest, marked with the red star, along with PDF plot illustrating LTD.CI (i.e., Confidence Intervals) have been estimated empirically based on inverse normal distribution.
In Fig. 6b), the extended 95% CI is highlighted with 2 dotted lines.Due to convergence, about conditioning parameter k , namely lim k→∞ p k ( ) = P( ) , with condition- ing parameter k = 6found to be convergent, for further evidence, see (Gaidai et al. 2023).Figure6b) displays a relatively low 95% CI, which is beneficial.As demonstrated in Fig. 6b), (Gaidai et al. 2023), MUR from the equationP( ) ≈ exp (−ν + ( ) T ) had been extrapolated in order to attain desired failure/hazard/risk threshold = 1 .With the target dynamic FOWT system's failure/ risk/hazard/damage rate, expressed in years −1 , and an expected FOWT system lifetime of L expected = T p( )N , according to Eq. ( 9), it can be concluded that the target dynamic FOWT system's LTD will follow Poisson-type PDF, with temporal parameter ν + (1).Hence, critical factor to assess expected FOWT system's lifetime being the MUR ν + , presented by Fig. 6b).Note that this is being methodological study, thus not a case study, hence authors did not aim to present exact numbers, diagrams and tables, but rather advocated novel Gaidai risks evaluation methodology, (Gaidai et al. 2023;Yakimov et al. 2023;Gaidai et al. 2023;Sun et al. 2023;Gaidai et al. 2023;Yakimov et al. 2023).Fig- ure6b) presents the PDF tail's extrapolation of 3 decimal orders of magnitude, meaning that the MC-type method's computational costs could be saved by about 10 3 times.With its robust statistical features and relatively simple approach, the Gaidai multivariate risks evaluation methodology appears to be a compelling approach.It often agrees with data-intensive POT (i.e., Peaks Over the Threshold), when examining extended return levels.Although the data's extreme value distribution may be fully represented nonparametrically by the Gaidai multivariate risks evaluation methodology, a specific class of parametric distributions was established in order to facilitate extrapolation to extended return period levels.The Gumbel-type PDF is the suitable asymptotic extreme value distribution in certain instances, which is the focus of this class of PDFs.This assumption is supported in the context of this chapter by the underlying statistics on water level, which always fall inside the zone of attraction of an asymptotic Gumbel-type PDF for the extremes.

Conclusions
Traditional timeseries risks evaluation design methods may not always be able to include complex cross-correlations across various system's critical components, due to system's high-dimensionality.Gaidai multivariate risks evaluation methodology's primary benefit lies within its capacity to assess lifetime distribution of multi-dimensional nonlinear nonstationary dynamic system.Current study investigated FOWT systems' performance, given realistic environmental settings, resulting in accurate assessment of FOWT lifetime distribution.The proposed approach's theoretical underpinnings have been presented in detail.In order to optimally utilize underlying datasets-which can be measured or MC numerically simulated-novel, accurate, dependable approaches have to be further developed.Evaluating structural failure or damage risk, based on direct measurements or extensive MC simulations is not always affordable.Current study's primary goal was to offer a multifunctional, trustworthy, easy-to-use solution for multidimensional offshore energy systems design.The proposed strategy has produced reasonably narrow confidence intervals; thus, the suggested approach can be helpful for nonlinear dynamic systems design.The examined FOWT example does not in any way limit potential use of the suggested methodology.Primary benefit of the advocated risk assessment methodology being is ability to treat reliability of multidimensional systems with practically unlimited number of dimensions or components.For non-stationary systems, the Gaidai risks evaluation approach can still be applied, utilizing an in-situ scatter diagram, composed of multiple short-term stationary sea states.In the latter case the underlying trend should be identified first -an example of an underlying trend could be a dynamic system with deterioration, where initial manufacturing faults cause local fractures to gradually grow, e.g., in welded hot spots.Material properties have also to be taken into account.

Fig. 6
Fig. 6 a) Non-dimensional system's vector − → R , index on the x-axis.b Extrapolation of p k ( ) towards failure/risk/damage level (star).An extrapolated 95% CI being indicated with 2 dotted lines.Vertical PDF plot illustrates LTD

Table 1
Main dimensions of semi-submersible-type FOWT platform