Artiﬁcial Neural Network-Based Prediction of the Extreme Response of Floating Offshore Wind Turbines under Operating Conditions

: The development of ﬂoating offshore wind turbines (FOWTs) is gradually moving into deeper offshore areas with more harsh environmental loads, and the corresponding structure response should be paid attention to. Safety assessments need to be conducted based on the evaluation of the long-term extreme response under operating conditions. However, the full long-term analysis method (FLTA) recommended by the design code for evaluating extreme response statistics requires signiﬁcant computational costs. In the present study, a power response prediction method for FOWT based on an artiﬁcial neural network algorithm is proposed. FOWT size, structure, and training algorithms from various artiﬁcial neural network models to determine optimal network parameters are investigated. A publicly available, high-quality operational dataset is used and processed by the Inverse First Order Reliability Method (IFORM), which signiﬁcantly reduces simulation time by selecting operating conditions and directly yielding extreme response statistics. Then sensitivity analysis is done regarding the number of neurons and validation check values. Finally, the alternative dataset is used to validate the model. Results show that the proposed neural network model is able to accurately predict the extreme response statistics of FOWT under realistic in situ operating conditions. A proper balance was achieved between prediction accuracy, computational costs, and the robustness of the model.


Introduction
FOWT has consistently advanced in recent years, bringing with it not only more environmentally friendly engineering applications but also financial rewards and scientific advancements [1].China, Japan, Portugal, Norway, and the UK will make up the top five markets for floating wind projects as of 2022 [2].Since most wind energy is found in deep water, FOWT technology is advancing quickly [3].Tension leg platform (TLP), Spar, semi-submersible (Semi), and barges are the three primary types of floating wind foundations now in use [4].However, as offshore development activities move further out into the ocean, FOWT will produce stronger structural vibrations when subjected to more demanding environmental loads, which could cause instability issues and significant fatigue damage to the FOWT structure and anchor cables [5].The motion response of a FOWT happens in six degrees of freedom (6DOF), which presents difficult design and evaluation problems [6].Effectively lowering the likelihood of accidents in the system is real-time or advance prediction of the dynamic response of the FOWT system, which provides early warning of risky circumstances.
It is necessary to calculate the structure's long-term severe reaction or loading for a specific recurrence period (once every 20 or 50 years) for any form of offshore structure, including the FOWTs [7].The approach known as Full Long-term Analysis (FLTA) is thought to be the most accurate one for determining how maritime constructions will react to extremely high loads [8].However, FLTA mandates simulation computations for at least 3000 different environmental scenarios, even if only a small number of them significantly affect how the structure responds.The Inverse First-Order Reliability Method (IFORM) is a well-liked technique to evaluate the ultimate loads of maritime structures with multidimensional variables and complicated nonlinear systems in order to minimize the number of simulations.IFORM was used to evaluate an oil tanker's dependability in various sea situations [9].It does not take a lot of calculations to estimate the joint long-term probability distributions of various FOWT severe loads using the IFORM [10].IFORM is also featured among the dependability techniques outlined and assessed for the design of ships and offshore buildings [11].
The FOWT is a strongly nonlinear system, and its dynamic response is usually forecasted using a time-domain coupled analysis method.This method is computationally intensive, inefficient, requires expensive storage space and computer time, and is inconvenient for quickly obtaining the system response or forecasting ahead of time.In order to obtain the extreme response of FOWT quickly, some researchers have improved the time-domain coupling analysis method [12].Research is being done on simplified and equivalent wave condition approaches that can decrease the number of computational conditions and increase the computational effectiveness of time domain evaluation methods [13,14].Additionally, real-time simulation methods are employed to track how the FOWT's structural performance responds to varied wind and wave loads [15,16].Although the computational effectiveness or reaction times of these systems have increased, they are still unable to accommodate the need for brief delays in actual use.
Some researchers have used artificial neural networks in the discipline of marine engineering as a result of the development of big data and artificial intelligence [17].The ANN algorithm, which has great computational efficiency and produces precise results, simply has to concentrate on the significant nonlinear relationship between the inputs and outputs.As a result, it is frequently utilized.An enhanced radial basis function (RBF) neural network strategy is provided for torque management of the FOWT system at speeds lower than the rated wind speed.This approach targets the external disturbances and nonlinear uncertain dynamic systems of the FOWT due to the proximity to load centers and strong wave coupling [18].Artificial neural network-based fuzzy logic control enhances the dynamic behavior of the FOWT platform and boosts its stability under a variety of wind and wave situations [19].An intelligent damage detection framework on FOWT was created using a multi-scale deep convolutional neural network model.The framework's performance is at least 10% better than the diagnostic precision attained with the majority of popular industrial-grade diagnostic tools [20].The results demonstrate the possibility of including artificial neural networks in the mooring design process.Artificial neural networks were used to simulate the time course of FOWT mooring tension in a range of irregular wave sea conditions [21].
As a result of their high inference accuracy, flexibility, resilience, fault tolerance, adaptability, scalability, generalization, and extrapolation capabilities in several successful applications, artificial neural network-based models have generally been shown to be more accurate in prediction [22].In research fields including FOWT behavioral control, fault analysis, generating power prediction, and dynamics modeling, machine learning is presently being employed more.There are, however, few techniques for accurately forecasting the FOWT systems' severe dynamic response values.As a result, a model for predicting extreme responses for FOWT enterprises is built in this research using artificial neural networks.To improve the accuracy of predicting FOWT extreme responses, the performance of shallow and deep ANN models with various dimensions is investigated, as is the impact of various network structures and the use of various training algorithms on ANN model performance and computational complexity.Artificial neural networks require a lot of input data to be trained, and the quality and quantity of that data directly impact the training outcomes.So, by processing marine environmental data from the TRL+ National Research Program (Spain) using IFORM, a sizable amount of training data was acquired.Additionally, a sensitivity study was carried out to look at the impact of changing the validation check values or the number of neurons on the model's ability to predict outcomes accurately.In order to explore the generated ANN model's generalizability and accuracy, the forecasting model's dependability was also confirmed.

Aerodynamic Load
The Blade Element Momentum Theory (BEM), which consists of the momentum theory and the blade element theory, is one of the primary techniques for estimating the aerodynamic stresses on wind turbine blades.According to the blade element hypothesis, the blade is divided into foliation units that are all infinitely thin and independent of one another [23].Figure 1 shows an illustration of an airfoil with the speeds and angles that define the forces acting on the element as well as the induced speeds brought on by the effect of the wake.Mathematical modeling of aerodynamic loads will be undertaken by AeroDyn [24].
neural networks.To improve the accuracy of predicting FOWT extreme res performance of shallow and deep ANN models with various dimensions is in as is the impact of various network structures and the use of various training on ANN model performance and computational complexity.Artificial neura require a lot of input data to be trained, and the quality and quantity of that d impact the training outcomes.So, by processing marine environmental da TRL+ National Research Program (Spain) using IFORM, a sizable amount of tr was acquired.Additionally, a sensitivity study was carried out to look at th changing the validation check values or the number of neurons on the mode predict outcomes accurately.In order to explore the generated ANN model's bility and accuracy, the forecasting model's dependability was also confirmed

Aerodynamic Load
The Blade Element Momentum Theory (BEM), which consists of the mom ory and the blade element theory, is one of the primary techniques for estimat odynamic stresses on wind turbine blades.According to the blade element hyp blade is divided into foliation units that are all infinitely thin and independen other [23].Figure 1 shows an illustration of an airfoil with the speeds and ang fine the forces acting on the element as well as the induced speeds brought on b of the wake.Mathematical modeling of aerodynamic loads will be undertak Dyn [24].
where ϕ denotes the angle of incidence; ∞ U the speed of the entering wind a the axial induction factor, and a′ the tangential induction factor; Ω sta wind turbine's angular velocity; r for its foliation to hub center distance; and regional tip speed ratio.
By incorporating the vane's aerodynamic properties throughout the radiu turbine's aerodynamic properties are determined.Figure 2 depicts the result namic forces on the element and its parts, both perpendicular to the rotor pla allel to it.The axial thrust T d and torque Q d on a blade element micro-se distance r from the leaf root and thickness r d are as follows [26,27]: The two components of the local velocity vector yield the angle of inflow [25].
where ϕ denotes the angle of incidence; U ∞ the speed of the entering wind at infinity; a the axial induction factor, and a the tangential induction factor; Ω stands for the wind turbine's angular velocity; r for its foliation to hub center distance; and λ r for its regional tip speed ratio.By incorporating the vane's aerodynamic properties throughout the radius, the wind turbine's aerodynamic properties are determined.Figure 2 depicts the resulting aerodynamic forces on the element and its parts, both perpendicular to the rotor plane and parallel to it.The axial thrust dT and torque dQ on a blade element micro-segment with distance r from the leaf root and thickness dr are as follows [26,27]:  The relationship between the forces and moments acting on the wind turb momentum and angular momentum is established by the momentum theory sumes that the loss of pressure or momentum in the plane of the wind turbin by the work done by the airflow in the plane of the wind turbine on the blad The following describes the axial thrust T d and torque M d on a micro-seg vane located at a point r away from the vane's root and having a thickness In this scenario, 1 v denotes the entering velocity, a the axial induction b the tangential induction factor.The iterative technique may be used to ca force, torque, and power operating on the wind turbine's wind wheel in conju the Blade Element Theory and momentum theory.

Hydrodynamic Load
In particular, the 3D potential flow theory is utilized to solve for the wave the floating foundation for large-scale constructions ( , D represe ameter of the circular member, and L signifies the wavelength) [28], as ill Equation (6).For small-scale structures ( ), Morrison's equations a solve the wave loads for slender structures, as shown in Equation (7) [29].Ma modeling of hydrodynamic loads will be undertaken by HydroDyn [30].
is hydrostatic force.In this equation, B stands for the number of blades, ρ for air density, c for air density, v 0 for relative incoming velocity, C L for lift coefficient, ϕ for incidence angle, and C D for drag coefficient.
The relationship between the forces and moments acting on the wind turbine and the momentum and angular momentum is established by the momentum theory, which assumes that the loss of pressure or momentum in the plane of the wind turbine is caused by the work done by the airflow in the plane of the wind turbine on the blade element.The following describes the axial thrust dT and torque dM on a micro-segment of the vane located at a point r away from the vane's root and having a thickness dr [26,27].
In this scenario, v 1 denotes the entering velocity, a the axial induction factor, and b the tangential induction factor.The iterative technique may be used to calculate the force, torque, and power operating on the wind turbine's wind wheel in conjunction with the Blade Element Theory and momentum theory.

Hydrodynamic Load
In particular, the 3D potential flow theory is utilized to solve for the wave stresses on the floating foundation for large-scale constructions (D/L > 0.2, D represents the diameter of the circular member, and L signifies the wavelength) [28], as illustrated in Equation ( 6).For small-scale structures (D/L < 0.2), Morrison's equations are used to solve the wave loads for slender structures, as shown in Equation (7) [29].Mathematical modeling of hydrodynamic loads will be undertaken by HydroDyn [30]. where ξ j is radiation force and ρgV 0 δ i3 − C h ij ξ j is hydrostatic force.
where A represents the cross-sectional area of the elongated structure, D represents the diameter, ρ represents the seawater density, u represents the velocity of the wave mass, C m represents the added mass, C d represents the drag coefficient, and ξ is the amplitude of the wave height.

Mooring Load
As shown in Figure 3, the lumped mass method discretizes the mooring anchor chain into n segments of uniformly equal length with massed nodes and massless spring damping, which are connected through n + 1 nodes to form the mooring anchor chain [31].The mathematical modeling of the lumped mass method will be undertaken by MoorDyn [32].
mass, m C represents the added mass, d C represents the dra the amplitude of the wave height.

Mooring Load
As shown in Figure 3, the lumped mass method discretizes th into n segments of uniformly equal length with massed nod damping, which are connected through 1 + n nodes to form th [31].The mathematical modeling of the lumped mass method MoorDyn [32].This technique has the benefit of accounting for the cable's se as well as the dynamic impacts of current drag and inertia forces.T a thorough explanation of how the equations were derived [33] equations of each node are given here as follows: where i m and i a stand for the mass and added mass of e for the internal stiffness and dam gravity and the force of contact with the seafloor; and pi D F F + each node.

Description of the FOWT Model
The NREL 5MW horizontal axis wind turbine model is dep developed using data from numerous wind power research and d is now widely used in the field of research pertaining to wind po wind turbine parameters.This technique has the benefit of accounting for the cable's self-weight and buoyancy as well as the dynamic impacts of current drag and inertia forces.The publication provides a thorough explanation of how the equations were derived [33], and only the dynamic equations of each node are given here as follows: where m i and a i stand for the mass and added mass of each node, respectively; for the internal stiffness and damping force; W i + F B for gravity and the force of contact with the seafloor; and F D pi + F D qi for the drag force at each node.

Description of the FOWT Model
The NREL 5MW horizontal axis wind turbine model is depicted in Figure 4.It was developed using data from numerous wind power research and development projects and is now widely used in the field of research pertaining to wind power.Table 2 lists the key wind turbine parameters.Semi-submersible FOWTs will be the focus of future development due to their short draft and outstanding stability, according to the U.S. Offshore Wind Technology Market Report [3].The support platform, as seen in Figure 4, primarily consists of a central column, three offset outside columns, and several thin supports.The primary characteristics of the DeepCwind semi-submersible platform are listed in Table 2.The three mooring anchor chains are evenly spaced along the axial direction of the middle short columns, and one end of each mooring chain is connected to the edge short columns through the guide holes and anchored to the seabed at the other end.The mooring anchor chain uses a suspension chain line type to provide the restoring force.In this study, more information on the mooring chain is provided [35].Semi-submersible FOWTs will be the focus of future development due to their short draft and outstanding stability, according to the U.S. Offshore Wind Technology Market Report [3].The support platform, as seen in Figure 4, primarily consists of a central column, three offset outside columns, and several thin supports.The primary characteristics of the DeepCwind semi-submersible platform are listed in Table 1.The three mooring anchor chains are evenly spaced along the axial direction of the middle short columns, and one end of each mooring chain is connected to the edge short columns through the guide holes and anchored to the seabed at the other end.The mooring anchor chain uses a suspension chain line type to provide the restoring force.In this study, more information on the mooring chain is provided [35].

Turbine Properties Value
Rotor

Marine Environmental Data
The open-ocean test facility BiMEP, which is 2 km offshore of the Spanish coastal town of Armintza, provided the data utilized in this simulation.The BiMEP test center, which has a total surface area of 5.2 km 2 and depths ranging from 50 to 90 m, is depicted in Figure 5a.The WAVESCAN directional buoys set up at BiMEP will gather marine meteorological data (currents, winds, and waves) and send it in real time to coastal stations through satellite communications, as illustrated in Figure 5b.

Marine Environmental Data
The open-ocean test facility BiMEP, which is 2 km offshore of the Spanish coastal town of Armintza, provided the data utilized in this simulation.The BiMEP test center, which has a total surface area of 5.2 km 2 and depths ranging from 50 to 90 m, is depicted in Figure 5a.The WAVESCAN directional buoys set up at BiMEP will gather marine meteorological data (currents, winds, and waves) and send it in real time to coastal stations through satellite communications, as illustrated in Figure 5b.In this work, the Inverse First Order Reliability Method (IFORM) is used to identify the severe sea state, and greater emphasis is placed on the accurate prediction of wind turbine loads in severe sea conditions [39], which can obtain maritime meteorological conditions over a long period.IFORM is used to create environmental contours on which design circumstances are chosen to compute the reaction, greatly cutting down on the simulation time [40].IFORM was used to analyze data from the TRL+ National Research Project to determine the effects of long-term severe sea conditions on FOWTs [37,38], and Figure 6 displays the wave characteristics and the distributions of wind speed parameters (90 m above sea level).In this work, the Inverse First Order Reliability Method (IFORM) is used to identify the severe sea state, and greater emphasis is placed on the accurate prediction of wind turbine loads in severe sea conditions [39], which can obtain maritime meteorological conditions over a long period.IFORM is used to create environmental contours on which design circumstances are chosen to compute the reaction, greatly cutting down on the simulation time [40].IFORM was used to analyze data from the TRL+ National Research Project to determine the effects of long-term severe sea conditions on FOWTs [37,38], and Figure 6 displays the wave characteristics and the distributions of wind speed parameters (90 m above sea level).
Considering the offshore standard [41], the four loading conditions were chosen from the 50-year environmental contour (yellow) in Figure .6 because the maritime environment chosen for the FOWT must have a return duration of 50 years.Separating the datasets correctly can increase the artificial neural network's dependability [42].The neural network model must be trained on a huge volume of data, and by discretizing the aforementioned four working conditions, additional environmental factors may be gained.
ditions over a long period.IFORM is used to create environmental contours on which design circumstances are chosen to compute the reaction, greatly cutting down on the simulation time [40].IFORM was used to analyze data from the TRL+ National Research Project to determine the effects of long-term severe sea conditions on FOWTs [37,38], and Figure 6 displays the wave characteristics and the distributions of wind speed parameters (90 m above sea level).For example, in the meaningful wave height interval H S ∈ (8.10, 8.81), the step size t = 0.05 is taken to generate 14 groups of meaningful wave height working conditions.In order to obtain 14 sets of wind speed conditions matching the meaningful wave height, the appropriate step size is taken through the wind speed interval W ∈ (7.93, 11.93).By analogy, another 29 groups of wind and wave loading data are obtained at H S ∈ (8.10, 9.50), W ∈ (11.93, 16.08) and H S ∈ (9.50, 10.30), W ∈ (16.08, 25.69).The total number of wind and wave conditions is 43.
To cover the spectral peak periods of short, medium, and long bands in the marine environment, three spectral peak periods were designed, namely 6 s, 12 s, and 24 s.In order to investigate the effect of different turbulence intensities on the load response of the floating wind turbine, three turbulence intensities were designed with reference to the requirements of the International Electrotechnical Commission (IEC) [43], which are 0.16, 0.14, and 0.12.With 43 wind and wave loading data available, 43 × 3 × 3 environmental parameters were obtained by considering three different spectral peak periods and turbulence intensities and were divided into three data sets, as shown in Table 3.

Load Simulation Calculation
A tool called OpenFAST [44] can simulate the coupled dynamic response of wind turbines using several physics models and levels of accuracy.To enable coupled nonlinear aero-hydro-servo-elastic modeling in the time domain, OpenFAST joins computational modules for aerodynamics, hydrodynamics for offshore structures, control and electrical system (servo) dynamics, and structural dynamics.As shown in Figure 7, AeroDyn [24], HydroDyn [30], and MoorDyn [32] in OpenFAST are responsible for calculating blade aerodynamic loads, hydrodynamic loads, and mooring loads, respectively.The FOWT coupled system dynamics model is established based on the above functional modules.In fact, other functional modules of OpenFAST are also involved in the computation, and the details can be found in this paper [44].A stochastic, full-field turbulent wind model is TurbSim [45].It simulates a threecomponent windspeed time series numerically at locations in a fixed, two-dimensional vertical rectangular grid using a statistical model.The turbulent wind field was generated using TurbSim.Referring to the turbulence intensity mentioned in the preceding section, choose the Kaimal turbulent wind spectrum and the NTM wind load model.
The full turbine dynamics simulation took ten minutes.To avoid being impacted by the FOWT's transitory behavior during startup, the simulation time was extended to 630 s, omitting the first 30 s, and the time step was changed to 0.0125 s.As shown in Figure 8, the blade root out-of-plane bending moment, tower base fore-aft bending moment, and mooring line tension of the FOWT in Condition 1 were obtained according to the above setup.All sea states are calculated, and the maximum value of the dynamic response is counted.Later, a neural network prediction model is developed using the environmental parameters as inputs and the maximum value of the dynamic response as outputs.A stochastic, full-field turbulent wind model is TurbSim [45].It simulates a threecomponent windspeed time series numerically at locations in a fixed, two-dimensional vertical rectangular grid using a statistical model.The turbulent wind field was generated using TurbSim.Referring to the turbulence intensity mentioned in the preceding section, choose the Kaimal turbulent wind spectrum and the NTM wind load model.
The full turbine dynamics simulation took ten minutes.To avoid being impacted by the FOWT's transitory behavior during startup, the simulation time was extended to 630 s, omitting the first 30 s, and the time step was changed to 0.0125 s.As shown in Figure 8, the blade root out-of-plane bending moment, tower base fore-aft bending moment, and mooring line tension of the FOWT in Condition 1 were obtained according to the above setup.All sea states are calculated, and the maximum value of the dynamic response is counted.Later, a neural network prediction model is developed using the environmental parameters as inputs and the maximum value of the dynamic response as outputs.A stochastic, full-field turbulent wind model is TurbSim [45].It simulates a threecomponent windspeed time series numerically at locations in a fixed, two-dimensional vertical rectangular grid using a statistical model.The turbulent wind field was generated using TurbSim.Referring to the turbulence intensity mentioned in the preceding section, choose the Kaimal turbulent wind spectrum and the NTM wind load model.
The full turbine dynamics simulation took ten minutes.To avoid being impacted by the FOWT's transitory behavior during startup, the simulation time was extended to 630 s, omitting the first 30 s, and the time step was changed to 0.0125 s.As shown in Figure 8, the blade root out-of-plane bending moment, tower base fore-aft bending moment, and mooring line tension of the FOWT in Condition 1 were obtained according to the above setup.All sea states are calculated, and the maximum value of the dynamic response is counted.Later, a neural network prediction model is developed using the environmental parameters as inputs and the maximum value of the dynamic response as outputs.

Artificial Neural Network Framework
In all varieties of artificial neural networks, the input, hidden, and output layers are all combined into a multi-layer feed-forward (MLFF) structure.It is the neural network structure of this paper, as seen in Figure 9.The maximum value of the blade root out-of-plane bending moment, the maximum value of the tower base fore-aft bending moment, and the maximum value of mooring line tension are selected as the model's output parameters.Wind speed, turbulence intensity, meaningful wave height, and spectral peak period are chosen as the model's input parameters.

Artificial Neural Network Framework
In all varieties of artificial neural networks, the input, hidden, and output layers are all combined into a multi-layer feed-forward (MLFF) structure.It is the neural network structure of this paper, as seen in Figure 9.The maximum value of the blade root out-ofplane bending moment, the maximum value of the tower base fore-aft bending moment, and the maximum value of mooring line tension are selected as the model's output parameters.Wind speed, turbulence intensity, meaningful wave height, and spectral peak period are chosen as the model's input parameters.

Let the input
.The activation function from the input layer to the hidden layer is f , and the activation func- tion from the hidden layer to the output layer is g .The process of signal smooth propa- gation is essentially the process of obtaining the output result by calculating the activation function, and the output k z of the k neuron in the output layer is as follows: One of the mechanical learning systems that uses strong learning from samples is an Artificial Neural Network (ANN), which can comprehend complicated nonlinear interactions [46].In this study, a neural network known as a back propagation (BP) neural network that exhibits both forward and backward propagation of error correction is used.Currently, complicated nonlinear systems of marine structures frequently employ this feed-forward back-propagation neural network [47].BP neural network uses error correction backward propagation, which first corrects the error from the output layer to the implicit layer before utilizing the gradient descent approach to change the weights and thresholds until the prediction of an artificial neural network closely resembles the desired outcome.Minimizing the mean square error of the desired value k t and the network output value k z is taken as the objective, and the mean square error (MSE) of the output layer is calculated as follows: Let the input x = (x 1 , x 2 , • • • , x n ) T , totaling n layer, and the implicit layer output be the output layer input y = (y 1 , y 2 , • • • , y h ) T , totaling h layer.Output layer z = (z 1 , z 2 , • • • , z m ) T , totaling m layers, w ij and θ ij denote the weights and thresholds from the i neuron in the input layer to the j neuron in the hidden layer, and w jk and θ jk denote the weights and thresholds from the j neuron in the hidden layer to the k neuron in the output layer, respectively The activation function from the input layer to the hidden layer is f , and the activation function from the hidden layer to the output layer is g.The process of signal smooth propagation is essentially the process of obtaining the output result by calculating the activation function, and the output z k of the k neuron in the output layer is as follows: One of the mechanical learning systems that uses strong learning from samples is an Artificial Neural Network (ANN), which can comprehend complicated nonlinear interactions [46].In this study, a neural network known as a back propagation (BP) neural network that exhibits both forward and backward propagation of error correction is used.Currently, complicated nonlinear systems of marine structures frequently employ this feedforward back-propagation neural network [47].BP neural network uses error correction backward propagation, which first corrects the error from the output layer to the implicit layer before utilizing the gradient descent approach to change the weights and thresholds until the prediction of an artificial neural network closely resembles the desired outcome.Minimizing the mean square error of the desired value t k and the network output value z k is taken as the objective, and the mean square error (MSE) of the output layer is calculated as follows: The accuracy of the output results, training duration, and network convergence speed are all impacted by wind, waves, and other maritime environment factors, as well as by the vast numerical range of the eigenvalues of the dynamic response of the floating wind turbine system.As a result, before training the neural network, the data must be normalized.The bipolar activation function's normalization formula is shown below.
with x being the normalized data, and the data set's minimum and maximum values are represented, respectively, by the letters x min and x max .
The entire dataset was randomly separated into three subsets: a training set (70%), a validation set (15%), and a test set (15%) in order to prevent the overfitting issue [48,49].Only the training dataset is used to train the network; the validation dataset is used to test the network's generalizability and to stop the training process before it becomes overfit.The artificial neural network was independently tested using the test dataset.
A floating wind turbine is a complex dynamical system, and it is time-consuming to obtain the power response from modeling to calculation.As shown in Figure 10, artificial neural networks may be useful to accelerate FOWT power response predictions.The inputs to the extreme response forecasting model are the environmental parameters acting on the FOWT, and the output is the extreme response of the FOWT.
The accuracy of the output results, training duration, and network convergence speed are all impacted by wind, waves, and other maritime environment factors, as well as by the vast numerical range of the eigenvalues of the dynamic response of the floating wind turbine system.As a result, before training the neural network, the data must be normalized.The bipolar activation function's normalization formula is shown below.x and max x .
The entire dataset was randomly separated into three subsets: a training set (70%), a validation set (15%), and a test set (15%) in order to prevent the overfitting issue [48,49].Only the training dataset is used to train the network; the validation dataset is used to test the network's generalizability and to stop the training process before it becomes overfit.The artificial neural network was independently tested using the test dataset.
A floating wind turbine is a complex dynamical system, and it is time-consuming to obtain the power response from modeling to calculation.As shown in Figure 10, artificial neural networks may be useful to accelerate FOWT power response predictions.The inputs to the extreme response forecasting model are the environmental parameters acting on the FOWT, and the output is the extreme response of the FOWT.

Training of Neural Networks
To make an accurate forecast, choosing the appropriate network topology and parameters is crucial [50,51].In order to better understand how different artificial neural network characteristics, particularly the hidden layer network topology and training technique, affect prediction performance.Based on MATLAB, a BP neural network was created [52], and a series of activation functions and training functions were provided in the MATLAB toolbox.BP neural networks require the activation function to be derivable, so a bipolar activation function (tanning-sigmoid) is used as an implicit layer [53].The output layer is a pure linear activation function (purelin) [54].The training procedures for the other datasets are identical; however, in the following, training for dataset 1 will be used as an example.Twenty computations will be performed using the same network configuration, and the accuracy of the forecasted results will be assessed using the mean square error average.After the neural network's fundamental structure has been established, it must be trained.
Various training techniques will be used, and the results will be examined.

Hidden Layer Network Structure
The number of hidden layers and the number of neurons in each hidden layer together form the network structure.The network structure will have a significant impact on prediction accuracy.For the BP neural network, the number of hidden layers can be set arbitrarily, but in fact, the more hidden layers, the more complex the network structure will be, and the more complex the calculation will be.For most engineering problems, the neural network with one or two hidden layers can complete the fitting task of almost all mapping relations [55].By changing the number of hidden layers and the number of neurons in each hidden layer, different network structures are obtained, and the accuracy of the prediction model under the network structure is calculated.In order to determine the number of hidden layers, neural networks with 1, 2, 3, and 4 hidden layers were tried, respectively.As shown in Figure 11, a neural network with four hidden layers has a maximum of nine neurons per layer, and there are a total of 6561 different network structures.In order to avoid contingency, each network structure needs to run 20 times to take the average, and finally, the network structure with the smallest error needs to be taken.The number of calculations is too large, so a batch file is written.The data set 1 was used for training.After calculation, the number of neurons with the smallest error under different hidden layers was obtained, as shown in Table 4.
MATLAB toolbox.BP neural networks require the activation function to be deriv a bipolar activation function (tanning-sigmoid) is used as an implicit layer [53].Th layer is a pure linear activation function (purelin) [54].The training procedure other datasets are identical; however, in the following, training for dataset 1 will as an example.Twenty computations will be performed using the same network ration, and the accuracy of the forecasted results will be assessed using the mea error average.After the neural network's fundamental structure has been establ must be trained.Various training techniques will be used, and the results will b ined.

Hidden Layer Network Structure
The number of hidden layers and the number of neurons in each hidden gether form the network structure.The network structure will have a significan on prediction accuracy.For the BP neural network, the number of hidden layer set arbitrarily, but in fact, the more hidden layers, the more complex the network s will be, and the more complex the calculation will be.For most engineering probl neural network with one or two hidden layers can complete the fitting task of a mapping relations [55].By changing the number of hidden layers and the numbe rons in each hidden layer, different network structures are obtained, and the acc the prediction model under the network structure is calculated.In order to determ number of hidden layers, neural networks with 1, 2, 3, and 4 hidden layers we respectively.As shown in Figure 11, a neural network with four hidden layers ha imum of nine neurons per layer, and there are a total of 6561 different network str In order to avoid contingency, each network structure needs to run 20 times to average, and finally, the network structure with the smallest error needs to be ta number of calculations is too large, so a batch file is written.The data set 1 was training.After calculation, the number of neurons with the smallest error under hidden layers was obtained, as shown in Table 4.By increasing the number of implied layers, the mean square error of the computed results of the power response prediction model is significantly reduced.Too many hidden layers can lead to the negative effect of overfitting.The maximum value of the mooring chain tension is much more accurately represented in the output response than are the maximum values of the external bending moments at the leaf root surface and the bending moments in front of and behind the tower base.The mean square error of the computed results of the power response prediction model is greatly decreased by increasing the number of assumed layers.Overfitting might be a problem if there are too many covert layers.The accuracy of the maximum value of the mooring line tension in the output response is much greater than that of the blade root out-of-plane bending moment and the tower base fore-aft bending moment.
As shown in Figure 12, the output of the forecasting model under the 2(4-9) network structure is compared with the expectation.The training, validation, and test sets' respective regression coefficients are 0.99424, 0.98969, and 0.98614.A better fit is denoted by a bigger R. The created neural network model has a significant capacity for generalization, as evidenced by the regression coefficient R of all the dataset training reaching 0.99275 and the good agreement between the neural network projected value and the anticipated value.
By increasing the number of implied layers, the mean square error of the computed results of the power response prediction model is significantly reduced.Too many hidden layers can lead to the negative effect of overfitting.The maximum value of the mooring chain tension is much more accurately represented in the output response than are the maximum values of the external bending moments at the leaf root surface and the bending moments in front of and behind the tower base.The mean square error of the computed results of the power response prediction model is greatly decreased by increasing the number of assumed layers.Overfitting might be a problem if there are too many covert layers.The accuracy of the maximum value of the mooring line tension in the output response is much greater than that of the blade root out-of-plane bending moment and the tower base fore-aft bending moment.As shown in Figure 12, the output of the forecasting model under the 2(4-9) network structure is compared with the expectation.The training, validation, and test sets' respective regression coefficients are 0.99424, 0.98969, and 0.98614.A better fit is denoted by a bigger R. The created neural network model has a significant capacity for generalization, as evidenced by the regression coefficient R of all the dataset training reaching 0.99275 and the good agreement between the neural network projected value and the anticipated value.

Comparison of Training Algorithms
By using the network structure of three hidden layers as an example, different training algorithms' effects on the accuracy of prediction models may be further examined.The first, second, and third layers of this network each contain 3, 6, and 9 neurons, respectively.The selection of different BP neural network training algorithms for prediction.The adaptive learning rate gradient descent method [56] and the scaled conjugate gradient backpropagation [57] are significantly more efficient compared to the most basic gradient descent method, but the computational accuracy does not meet the requirements; see Table 5. Bayesian regularization [58] and Levenberg-Marquardt [59] with significantly higher computational accuracy.In reality, as rapid backpropagation techniques are highly advised as first-choice supervised algorithms, the Levenberg-Marquardt and Bayesian regularization backpropagation algorithms are considered the finest options for such nonlinear dynamics.For this simulation, Levenberg-Marquardt has no additional training time and is recommended as the neural network training algorithm for dynamic response prediction models.

Results and Discussion
First, the sensitivity of the model's performance metrics to the number of neurons and the number of validations will be discussed.Then the network structure with the highest computational accuracy will be selected as the structure of the power response prediction model, and the reliability of the model will be verified with a new data set.

Sensitivity Analysis
Determination of the number of neurons is a more complex problem.In order to study the sensitivity of the neural network model performance to the number of neurons and the number of validations, three additional evaluation metrics, RMSE, MAE, and Coefficient of Determination R 2 , were added.The three responses were evaluated as a whole without calculating the evaluation metrics for each response variable separately.
The root mean square error (RMSE) can be obtained by taking the square root of the MSE, which is more sensitive to larger error values.
Mean Absolute Error (MAE) does not take into account the square of the variance, so it does not amplify the squared difference in the variance values and is more intuitive relative to MSE, which reflects the absolute magnitude of the prediction error rather than the squared magnitude of the error.
Coefficient of Determination R 2 , which indicates the proportion of variability in the dependent variable that can be explained by the model, and how well the model fits the data.
where t k is the mean of the expected value t k .Smaller MSE, RMSE, and MAE represent better model performance.On the contrary, larger t k represents better model performance.Using neural network with two hidden layers, the number of neurons in each layer is continuously increased to a maximum of 50.By changing the number of neurons, the neural network model can make better predictions as shown in Figure 13.For example, as the number of neurons increases to 15, lower prediction errors and higher coefficients of determination occur.However, as the number of neurons exceeds 15, the performance of the artificial neural network model begins to gradually degrade, with higher error and risk of overfitting.Therefore, it is recommended to have no more than 15 neurons per layer.
Coefficient of Determination R 2 , which indicates the proportion of variability in the dependent variable that can be explained by the model, and how well the model fits the data.

 
where k t is the mean of the expected value k t .Smaller MSE, RMSE, and MAE represent better model performance.On the contrary, larger k t represents better model perfor- mance.Using neural network with two hidden layers, the number of neurons in each layer is continuously increased to a maximum of 50.By changing the number of neurons, the neural network model can make better predictions as shown in Figure 13.For example, as the number of neurons increases to 15, lower prediction errors and higher coefficients of determination occur.However, as the number of neurons exceeds 15, the performance of the artificial neural network model begins to gradually degrade, with higher error and risk of overfitting.Therefore, it is recommended to have no more than 15 neurons per layer.Under the optimal network structure, the neural network with two hidden layers can reduce the MSE, RMSE, and MAE by about 29.55%, 16.57%, and 19.86%, respectively, Under the optimal network structure, the neural network with two hidden layers can reduce the MSE, RMSE, and MAE by about 29.55%, 16.57%, and 19.86%, respectively, compared to the neural network with a single layer.In addition, the R 2 can be improved by 1.44%; see Table 6.The validation count is used to validate the dataset.Specifically, if the current validation count exceeds the set value, the neural network will terminate the training to avoid overfitting, thus improving the generalization ability of the network.The proposed neural network with two hidden layers is trained by varying the set value of the validation count to obtain the performance metrics of the neural network with different validation counts.As shown in Figure 14d, the number of validations increases from 6 to 50, and the training time increases by 95.42%.However, there is almost no improvement in performance metrics; see Figure 14a-c.In addition, using a validation count lower than 6 leads to unacceptable errors.Therefore, a validation check number of 6 is appropriate.compared to the neural network with a single layer.In addition, the R 2 can be improved by 1.44%; see Table 6.The validation count is used to validate the dataset.Specifically, if the current validation count exceeds the set value, the neural network will terminate the training to avoid overfitting, thus improving the generalization ability of the network.The proposed neural network with two hidden layers is trained by varying the set value of the validation count to obtain the performance metrics of the neural network with different validation counts.As shown in Figure 14d, the number of validations increases from 6 to 50, and the training time increases by 95.42%.However, there is almost no improvement in performance metrics; see Figure 14a-c.In addition, using a validation count lower than 6 leads to unacceptable errors.Therefore, a validation check number of 6 is appropriate.

Reliability of Predictive Models
Extreme sea states will bring more dangers to offshore floating wind turbines in service, and accurate and fast prediction of the dynamic response of FOWT in extreme sea states is crucial to ensuring the safety of floating systems during service.In order to verify the reliability of the trained artificial neural network model, the environmental loading conditions used to verify the reliability of the model were prepared, which are different from the environmental loading conditions used for neural network training; see Table 7.A total of 12 environmental loading conditions used to validate the reliability of the model were categorized into four datasets.Data sets 1, 2, and 3 have different spectral peak periods, which represent the short, medium, and long wavelengths of the ocean environment, and the three conditions in each data set have the same meaningful wave height, spectral peak period, wind speed, and different turbulence intensities.For dataset 4, which represents the extreme ocean environment, the wind speed is close to the cut-out speed of the NREL 5MW wind turbine, and the dataset takes into account both different spectral peak periods and turbulence intensities.The environmental load conditions used to verify the reliability of the model are substituted into the trained artificial neural network model, and the predicted values of the maximum blade root out-of-plane bending moment, the maximum tower base fore-aft bending moment, and the maximum mooring line tension are obtained.The comparison between the predicted values of the artificial neural network and the OpenFAST simulation values is shown in Figure 15a-c.Data sets 1, 2, and 3 have the same meaningful wave height, spectral peak period, and wind speed.But the turbulence intensity is different.From Figure 15a, it can be seen that the effect of turbulence intensity on the maximum blade root out-of-plane bending moment is obvious, and the larger the turbulence intensity, the larger the maximum blade root out-of-plane bending moment.For data set 4, the average wind speed of the data set is close to the cut-out wind speed of 25 m/s.For the safety of FOWT, the control system will automatically adjust the pitch angle to reduce the blade root out-of-plane bending moment.As shown in Figure 14a, although the wind speed of conditions 10 to 12 is larger than that of conditions 1 to 9, the blade root out-of-plane bending moment is decreasing.It can be seen that the trained neural network model successfully predicts the variation of load with wind speed.
Using OpenFAST simulation to obtain the extreme response for 12 operating conditions takes 12 × 630 s.If the proposed artificial neural network prediction model is used, it takes 45 s (Laptop system: Windows 10 home, processor: Intel (R) Core (TM) i7-7500U CPU @ 2.70 GHz 2.90 GHz).The 99.4% time reduction and 167-fold efficiency improvement compared to the OpenFAST simulation demonstrate the computational efficiency advantages of this extreme response prediction model.The MSE, RMSE, MAE, and R 2 of the artificial neural network prediction value and the OpenFAST simulation value are shown in Figure 15d.The MSE, RMSE, and MAE of the maximum mooring line tension are lower than the MSE, RMSE, and MAE of the maximum blade root out-of-plane bending moment of the blade root and the maximum tower base fore-aft bending moment, and the R 2 of the maximum mooring line tension is higher than the R 2 of the maximum blade root out-of-plane bending moment of the blade root and the maximum tower base fore-aft bending moment.Smaller MSE, RMSE, and MAE represent higher accuracy.Conversely, a larger R 2 represents higher accuracy.Therefore, among the three output responses, the prediction accuracy of the maximum mooring line tension is significantly higher than that of the maximum blade root out-of-plane bending moment and the maximum tower base fore-aft bending moment.
Using OpenFAST simulation to obtain the extreme response for 12 operating conditions takes 12 × 630 s.If the proposed artificial neural network prediction model is used, it takes 45 s (Laptop system: Windows 10 home, processor: Intel (R) Core (TM) i7-7500U CPU @ 2.70 GHz 2.90 GHz).The 99.4% time reduction and 167-fold efficiency improvement com-pared to the OpenFAST simulation demonstrate the computational efficiency advantages of this extreme response prediction model.

Summary and Conclusions
The experimental physical quantities (EPQs) that can be obtained are very limited.Therefore, there is a practical concern for designers and industry to obtain accurate non-EPQs.Traditional theoretical and experimental approaches to studying FOWT extreme response statistics are often costly (i.e., model testing) or require time-consuming computational methods.
In order to accelerate FOWT reliability assessment, this study suggests the application of artificial neural networks.Marine environmental data from the TRL+ National Research Project (Spain) was used.The training process of the artificial neural network is highly dependent on the input data, and the number and accuracy of the input data directly affect the training results.Therefore, a large amount of training data was obtained by processing the marine environmental data of the TRL + National Research Project (Spain) through IFORM.TurbSim numerical tool was used to generate turbulent winds as an input to the OpenFAST (Releases v3.5.0) software tool.
The above data were normalized and used for the training of the neural network.The training is completed, and the artificial neural network model is validated, tested, and evaluated by four statistical features.In addition, a sensitivity analysis was performed to show the effect of different numbers of neurons and different numbers of validation checks on the performance of different artificial neural network models.Finally, the new dataset was used to validate the reliability of the proposed neural network extreme response prediction model.
Conclusions drawn through numerical simulation and analysis are as follows: (1) Application of IFORM to assess response for data-selected design cases of the TRL + National Research Project can significantly reduce simulation time.(2) Neuron sensitivity analysis has shown that deep neural networks with double hidden layers yield better predictions, compared to shallow neural networks with single hidden layers, however further increases in neuron numbers lead to a gradual decrease in model performance.(3) The sensitivity analysis of the number of validations shows that setting too few validations when training the neural network will lead to larger MSE, RMSE, MAE, and smaller R 2 , which indicates that the accuracy of the model is not high.However, too many verifications will take more time, and the improvement in accuracy is almost negligible.(4) Reliability verification has shown that a reasonable neural network, as well as appropriate parameter settings of the neural network given in situ wind conditions, may become a useful design strategy.(5) Blade element momentum theory, applied to key FOWT parts, wave kinematics, and hydrodynamics calculations, forms a solid mechanical model.Finally, given all the above-mentioned steps, an accurate extreme response statistic for FOWT key parts was obtained.
The developed neural network model is based on marine environment data from bimep (Biscay Marine Energy Platform), so it is not applicable to other wind farms.Each wind farm has a unique marine environment, and therefore, the neural network model needs to be designed for different wind farms.However, there is no need to develop separate grids for individual wind or sea conditions for a single wind farm, as all conditions can be included in the training of the neural network in the early stages, and in fact, the more types of data that are available, the more accurate the predictions will be.For developing a neural network model for a new wind farm, the biggest cost comes from the acquisition of environmental data for the wind farm.
Combining the TRL+ project, the inverse first-order reliability method, and machine learning is a new technical attempt.Future work could consider combining AI with the FOWT digital twin, which requires accurate and instantaneous display and monitoring of engineering structures in the real environment.The technique of combining the digital model of the wind turbine with the actual operating data of the wind turbine to achieve real-time output of the wind turbine loads can provide effective information to wind turbine researchers and developers to reduce the cost of wind turbine design, increase the means of wind turbine control, and optimize the wind turbine operation process.The neural network prediction model proposed in this paper has high accuracy, comparable to the OpenFAST simulation results, and is also fast, which can realize the instant calculation of the FOWT system response and the over-warning of the dangerous state, so as to take measures in advance to guarantee the operational safety of the FOWT system.

2
In this equation, B stands for the number of blades, ρ for air density density, 0 v for relative incoming velocity, L C for lift coefficient, ϕ for in gle, and D C for drag coefficient.

Figure 6 .
Figure 6.Marine Environment of Different Return Periods.(a) IFORM results of the significant wave height HS and wind speed W; (b) IFORM results of the HS and peak wave period TP (right) for the BiMEP area [34].
n layer, and the implicit layer output be the output layer input m layers, ij w and ij θ denote the weights and thresholds from the i neuron in the input layer to the j neuron in the hidden layer, and jk w and jk θ denote the weights and thresholds from the j neuron in the hidden layer to the k neuron in the output layer, respectively,
x ˆ being the normalized data, and the data set's minimum and maximum values are represented, respectively, by the letters min

Figure 11 .
Figure 11.Flowchart for determining the optimal network structure with four implicit lay

Figure 11 .
Figure 11.Flowchart for determining the optimal network structure with four implicit layers.

Figure 12 .
Figure 12.Comparison of Forecast Model Output and Expectations.

Figure 12 .
Figure 12.Comparison of Forecast Model Output and Expectations.

Figure 13 .
Figure 13.Performance indicators of neural networks.(a) Mean squared error; (b) Root mean square error; (c) Mean absolute error; (d) Coefficient of determination.

Figure 13 .
Figure 13.Performance indicators of neural networks.(a) Mean squared error; (b) Root mean square error; (c) Mean absolute error; (d) Coefficient of determination.

Figure 14 .
Figure 14.Impact of the number of validations during neural network training on performance metrics.(a) Mean squared error and Root mean square error; (b) Mean absolute error; (c) Coefficient of determination; (d) Training time.

Figure 15 . 15 .
Figure 15.Verify the reliability of the model, RootMyc1 (blade root out-of-plane bending moment), TwrBsMyt (tower base fore-aft bending moment).(a) Maximum blade root out-of-plane bending Figure 15.Verify the reliability of the model, RootMyc1 (blade root out-of-plane bending moment), TwrBsMyt (tower base fore-aft bending moment).(a) Maximum blade root out-of-plane bending moment; (b) Maximum tower base fore-aft bending moment; (c) Maximum mooring line tension; (d) Performance indicators.

Table 3 .
Environmental loading conditions for neural network training.

Table 4 .
The impact of hidden layer network structure on computational accuracy.

Table 4 .
The impact of hidden layer network structure on computational accuracy.

Table 5 .
Effect of different neural network algorithms on computational accuracy.

Table 6 .
Performance indicators of dual hidden layer neural networks and single hidden layer neural networks.

Table 6 .
Performance indicators of dual hidden layer neural networks and single hidden layer neural networks.

Table 7 .
New environmental parameters for verifying reliability.