John Marius Hegseth
John Marius Hegseth
Doctoral thesis
NTNU
Norwegian University of
Science and Technology
Thesis for the degree of
Philosophiae Doctor
Faculty of Engineering
Department of Marine Technology
Doctoral theses at NTNU, 2020:356
John Marius Hegseth
NO - 1598
Floating wind turbines (FWTs) are considered a promising solution for wind
energy harvesting in deep water, but are currently too expensive to compete
with other energy sources. Being a relatively new and immature technology
means that there still is a large potential for cost reductions through optim-
ization of the FWT structure. Optimized designs will bring the construction
costs of FWTs down and increase profitability, which is currently the major
challenge for the industry. Optimized designs can also result in increased
reliability, which is an important issue.
A FWT may encounter a large number of different loading situations dur-
ing its lifetime, including operational, parked, fault, start up, and shutdown
events. In each of these situations, various combinations of environmental
loads from wind and waves must be considered. Combined with computa-
tionally expensive response analyses, this results in a comprehensive design
process.
FWTs are also highly multidisciplinary systems, primarily combining the
areas of aerodynamics, hydrodynamics, structural dynamics, and control
theory. The interactions between the disciplines, and between different
components in the system, calls for integrated analysis and design. Due
to the complex dynamics, strong couplings, and potentially large number
of design variables, identifying optimal design solutions become a difficult
task.
The main purpose of this work was to improve the design process for FWTs
and thus contribute to reducing the cost of energy. This is addressed through
two overall research objectives, which consider i) increased computational
efficiency of global design analyses for FWTs, and ii) methods for numerical
iii
iv
design optimization which can help identify cost-effective and reliable design
solutions. The main focus was on the support structure and controller for
10 MW spar-type turbines, considering fatigue and ultimate loads.
A linearized aero-hydro-servo-elastic model was shown to yield good results
for the fatigue loads in the support structure, where agreement within ±
30 % was achieved for the long-term fatigue damage compared to nonlin-
ear time-domain simulations. Acceptable agreement was also observed for
short-term extreme response, especially for the support structure bending
moments, which were quite Gaussian also in harsh environmental condi-
tions. The resonant platform pitch response was overestimated by the linear
model, especially in near-rated conditions.
A gradient-based optimization approach with analytic derivatives was de-
veloped to perform integrated design optimization of the support structure,
blade-pitch controller, and mooring system for an elastic 10 MW spar FWT,
including the scantling design of the hull, where the goal was to minimize
a combination of design costs and rotor speed variation. Different con-
trol strategies were compared through integrated design of the controller
and support structure, which allowed for identification of optimal control
parameters in a lifetime perspective, and fair comparisons between different
strategies.
The impact of environmental modelling on the long-term fatigue reliability,
and associated design costs, of the support structure was assessed through
re-design of the tower and platform. Considering stochastic turbulence in-
tensity, wind-wave misalignment and the wind directional distribution re-
duced the long-term fatigue damage by approximately two-thirds along the
support structure, compared to the base model. Implications of the chosen
fatigue safety factor on the trade-offs between CAPEX and OPEX were also
assessed.
The methodologies for global response analyses and integrated design op-
timization developed in the present work have been shown to be suitable for
preliminary design of spar FWTs, where they can provide a starting design
for later and more detailed design phases. Different modelling and design
aspects for cost-effective and reliable solutions have been identified and as-
sessed. The methodologies can be further extended to account for different
FWT concepts, additional design parameters, and other load cases, and
may help identify novel design solutions.
Acknowledgements
v
vi
List of publications
The following five papers are considered part of this thesis:
vii
viii
The following paper is not regarded as part of the thesis due to scope:
Declaration of authorship
In P1, P3, P4, and P5, J. M. Hegseth had the main ideas for the papers,
developed the numerical models, performed all calculations, and wrote the
manuscripts. Profs. E. E. Bachynski, J. R. R. A. Martins and B. J. Leira
contributed with valuable suggestions and feedback during the processes.
P2 was initiated by C. E. Souza, who also developed the forced oscillation
method and wrote the manuscript. J. M. Hegseth derived the analytic
expressions for the aerodynamic damping and inertia coefficients, developed
the linearized model, and participated in discussions. Prof. E. E. Bachynski
provided guidance and feedback.
x
Contents
Abstract iii
Acknowledgements v
Publications vii
1 Introduction 1
1.1 Motivation and background . . . . . . . . . . . . . . . . . . . 1
1.2 Research objectives . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3 Main contributions . . . . . . . . . . . . . . . . . . . . . . . . 9
2 Literature Survey 11
2.1 Global response analyses . . . . . . . . . . . . . . . . . . . . . 11
2.1.1 Environmental modelling . . . . . . . . . . . . . . . . 11
2.1.2 State-of-the art simulations . . . . . . . . . . . . . . . 14
2.1.3 Linear models . . . . . . . . . . . . . . . . . . . . . . . 17
2.1.4 Structural reliability . . . . . . . . . . . . . . . . . . . 20
2.2 Design optimization . . . . . . . . . . . . . . . . . . . . . . . 21
2.2.1 Basic optimization theory . . . . . . . . . . . . . . . . 21
xi
xii CONTENTS
3 Numerical Models 35
3.1 Linearized formulation . . . . . . . . . . . . . . . . . . . . . . 35
3.2 System overview . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.3 Structural dynamics . . . . . . . . . . . . . . . . . . . . . . . 38
3.4 Hydrodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.5 Aerodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.6 Control system description . . . . . . . . . . . . . . . . . . . . 45
3.7 Response to stochastic input . . . . . . . . . . . . . . . . . . 47
3.7.1 Fatigue damage . . . . . . . . . . . . . . . . . . . . . . 47
3.7.2 Extreme response . . . . . . . . . . . . . . . . . . . . . 49
3.8 Optimization framework . . . . . . . . . . . . . . . . . . . . . 49
3.9 Model development . . . . . . . . . . . . . . . . . . . . . . . . 54
4 Research Findings 55
4.1 Linearized dynamic analyses . . . . . . . . . . . . . . . . . . . 55
4.1.1 Aerodynamic damping and inertia . . . . . . . . . . . 55
4.1.2 Long-term fatigue damage . . . . . . . . . . . . . . . . 59
4.1.3 Short-term extreme response . . . . . . . . . . . . . . 61
4.2 Integrated design optimization . . . . . . . . . . . . . . . . . 64
4.2.1 Multimodality . . . . . . . . . . . . . . . . . . . . . . 64
4.2.2 Cost and power quality . . . . . . . . . . . . . . . . . 65
4.2.3 Support structure design . . . . . . . . . . . . . . . . . 66
4.2.4 Control strategies . . . . . . . . . . . . . . . . . . . . . 72
CONTENTS xiii
A Appended Papers 99
Introduction
This chapter gives an overview of the background and motivation for this
work, and states the objectives of the research. A summary of the main
contributions to the scientific community, as well as their relations to the
research objectives and appended papers, are also presented.
1
2 Introduction
Table 1.1: Floating wind farms in Europe with expected commissioning date
within 2022 (WindEurope 2019).
(The Carbon Trust 2015). However, moving the energy production to deeper
water and larger distances from shore increases the complexity and requires
innovative technology. To date, the world’s only floating wind farm is the
Hywind Scotland pilot park, which consists of five turbines with an installed
capacity of 30 MW. However, more than 30 different FWT concepts have
been proposed (Leimeister et al. 2018), and several floating wind farms are
planned to start operation in the near future. Seven European farms, listed
in Table 1.1, are expected to be commissioned within the next three years,
with a total capacity of 274 MW.
The cost of floating wind is currently higher than for bottom-fixed turbines,
and large reductions are needed to reach a competitive level. Floating wind
turbines have larger construction costs than the bottom-fixed alternative;
however, installation costs are expected to be smaller, because the turbines
can be assembled at port using onshore cranes instead of expensive jack-up
vessels at site (Katsouris and Marina 2016). The current levelized cost of
energy (LCOE) is estimated to e180-200/MWh for pre-commercial projects,
but large reductions are expected as the technology matures and reaches
commercial-scale deployment, and may attain values of e40-60/MWh in
2030 (WindEurope 2018).
Lower LCOE can be achieved in several ways, such as reduction of capital
1.1. Motivation and background 3
Figure 1.1: Estimated breakdown of capital costs for commercial scale floating
wind farms (The Carbon Trust 2015).
Offshore wind turbines are highly dynamic systems subjected to static, peri-
odic, stochastic, and transient loads. Combined with compliant structures
and a relatively long design lifetime, this often leads to FLS being design
driving for several components (Vorpahl et al. 2013). Contrary to the design
of conventional offshore structures, which tends to be driven by ultimate
loads, this is also the case for the floating platform (Strach-Sonsalla and
Muskulus 2016). To give a realistic estimate of the fatigue life of the struc-
ture, all relevant environmental conditions throughout the lifetime of the
system must be considered, which leads to a comprehensive set of design
calculations. This is especially the case if directionality of wind and waves
is considered, which is of particular importance for fatigue design (DNV GL
2018a).
In addition to fatigue, the survival of the platform and turbine in extreme
1.1. Motivation and background 5
Sd ≤ R d , (1.1)
where Sd is the design load effect, and Rd is the design resistance. These
are usually found using the partial safety factor method, where the different
characteristic load effects and resistances are multiplied by individual load
and material factors to obtain their design values. The safety factors reflect
the uncertainty in the respective design parameters, as well as the required
reliability level.
Probability-based design is an alternative to the partial safety factor method,
where the load and resistance variables are described by their probability
distributions. The design of the system is then performed to meet a safety
level which is expressed by the failure probability:
• Power production
6 Introduction
• Start up
• Normal shutdown
• Emergency shutdown
• Parked turbine
Within each of these design situations, a subset of load cases are specified.
The subsets consider different wind and wave conditions, fault types, and
limit states, and the total number of load cases in a design process may
therefore add up to several thousands.
The goal of the design optimization process is to find the design that min-
imizes or maximizes a given performance measure, such as costs or power
production, and at the same time meets the appropriate design criteria. In
traditional design optimization, the performance of an initial design is eval-
uated through analyses or experiments. Based on the results, the designer
makes changes to the design, and the process is repeated until a satisfact-
ory design is achieved. This process is typically limited to a few iterations,
where changes to the design to a large extent depend on the experience of
the designer (Muskulus and Schafhirt 2014).
In contrast, numerical optimization techniques can efficiently explore large
parts of the design space to identify improved and possibly non-intuitive
solutions, as illustrated in Fig. 1.2. Such techniques have been widely util-
ized in the automotive and aerospace industries, but have seen limited use
in design of wind turbine structures (Muskulus and Schafhirt 2014). Simil-
arly to a conventional design process, the performance of an initial design is
evaluated based on results from analyses, and iterations are performed until
a satisfactory solution has been achieved. The key differences are that the
design changes are made automatically by an optimization algorithm, and
that the design process does not end until certain optimality criteria are
satisfied. These criteria ensure that the final design not only complies with
the relevant constraints, but that it indeed is the solution that minimizes
or maximizes the chosen performance measure, at least in that particular
neighbourhood of the design space (Martins and Ning 2020).
1.1. Motivation and background 7
Optimization prob-
lem formulation
No No
Figure 1.2: Comparison of main steps in a manual and numerical design optim-
ization process. Adapted from Muskulus and Schafhirt (2014).
integrated design process, all relevant disciplines are coupled, and the sub-
systems are designed simultaneously. An integrated optimization approach
involves a large number of design variables and load cases, and requires
sophisticated and multidisciplinary numerical tools. Although the overall
system in terms of total cost and energy production is the whole wind farm,
the studies are limited to a single turbine in the present work, and farm-level
effects are not considered.
improved and possibly novel solutions, but have seen limited use for FWT
structures. Consequently, there is currently no software or framework avail-
able for efficient FWT optimization, and there are uncertainties related to
the formulation of the optimization problem. Design methodologies with
improved accuracy and efficiency, which can help identify innovative and
optimized solutions, are expected to result in cost reductions and increased
reliability, and can help floating wind technology become commercially feas-
ible.
The second research objective addresses the potential for design improve-
ments through numerical optimization, and is formulated as:
O1 O2
Increased Numerical
computational design
efficiency optimization
P1 P2 P3 P4 P5
C1 C2 C3 C4 C5 C6
Linearized Damping Integrated Control Modelling of Inspection
model and inertia optimization strategies environment planning
Literature Survey
This chapter summarizes previous work relevant for global response ana-
lyses and numerical design optimization of FWT structures. Methods used
in state-of-the-art and simplified dynamic simulations are described, some
basic optimization theory is given, and specific aspects and challenges re-
lated to design optimization of FWTs are discussed.
11
12 Literature Survey
Figure 2.1: Relative wind-wave direction distributions for two wind speed bins
at a North Sea location 25 km west of Norway.
The wave conditions in a sea state can be separated into two categories:
wind seas, which are generated by local wind, and swell, which are waves
that have travelled out of the areas where they were generated (DNV GL
2019a). For wind seas, the Pierson–Moskowitz (PM) and JONSWAP spec-
tra are frequently applied (DNV GL 2019a). The JONSWAP spectrum,
which was developed for North Sea conditions, is used in the present work.
If the sea state is expected to contain swell components in addition to wind-
generated waves, these must also be properly represented (DNV GL 2018a).
This can be done through two-peaked spectra, such as the Torsethaugen
model (Torsethaugen 1996, Torsethaugen and Haver 2004), which take into
account both wind and swell sea. Alternatively, the wind and swell compon-
ents can be separated and described using two JONSWAP spectra, possibly
with different directions, which are superimposed to represent the total wave
elevation process (DNV GL 2018a).
2.1.2 State-of-the art simulations
Due to important interactions between aerodynamics, hydrodynamics, struc-
tural dynamics, and control dynamics, fully coupled aero-hydro-servo-elastic
models are needed to properly simulate the global response of FWTs. Dif-
ferent aspects and interactions considered in such analyses are illustrated
in Fig. 2.2. Over the last decade, several fully coupled simulation tools
for global analysis of FWTs have been developed, often based on existing
codes for either conventional offshore structures or land-based wind tur-
bines. The present thesis uses the SIMA workbench developed by SINTEF
Ocean, which couples the RIFLEX and SIMO programs (SINTEF Ocean
2016a;b). A brief summary of common structural, aerodynamic, hydro-
dynamic, and control capabilities in state-of-the-art analysis tools is given
in the following.
The structural dynamics of FWTs are usually modelled using multibody dy-
namics, modal analysis, finite elements, or a combination of these (Robertson
et al. 2014). While modal analysis is computationally efficient, finite ele-
ments offer the possibility of geometric and material nonlinearities, and
do not require pre-calculation of the system modes. Most finite element
models for FWTs are limited to beam elements, which are applicable for
long, slender structures such as the tower and blades, whereas the large
volume platform typically is considered rigid. The drivetrain is usually
modelled using a multibody dynamics approach, where the shaft flexibility
often is included through a torsional spring-damper to capture the first tor-
sional mode of the drivetrain (Smilden 2019). The mooring lines are either
2.1. Global response analyses 15
Control
system
Drive-
Wind Rotor Grid
train
Nacelle
Rotor-nacelle assembly
Tower
Mooring
system
et al. 2015).
Frequency-domain models have a long history for offshore structures, which
is reflected in the early efforts on frequency-domain modelling of FWTs
(Wayman et al. 2006, Tracy 2007, Sclavounos et al. 2008, Philippe et al.
2011), where response amplitude operators (RAOs) were established based
on diffraction and radiation loads from an external potential flow solver. As
hydrodynamic loads from numerical panel codes are computed in the fre-
quency domain, they can be readily applied in such global response models
without modification. A drawback with this approach in a design process
is that it requires an additional hydrodynamic analysis each time changes
are made to the design, which must be run prior to the response analysis,
and prevents quick examination of different design solutions. More compu-
tationally efficient options include Morison’s equation and MacCamy–Fuchs
theory; however, the validity of these methods is limited to simple geomet-
ries.
While early studies on frequency-domain methods for FWTs considered
only rigid-body modes, Bachynski and Moan (2012) compared a frequency-
domain model for a tension-leg platform (TLP) wind turbine with nonlinear
time domain simulations, and found that the linear model was insufficient
for design calculations, partly due to the rigid modelling of the tower. The
effect of tower flexibility on the global response was confirmed by full-scale
measurements from the Hywind Demo FWT, where the tower bending mo-
ment spectra clearly showed peaks around the first bending mode (Skaare
et al. 2015). Kvittem and Moan (2014) used a frequency-domain model,
which included the first tower bending mode, to evaluate the fatigue loads
for a semi-submersible FWT in stochastic wind-wave conditions, where the
shape of the bending mode was found from decay tests in calm air with a
coupled time-domain model. A stochastic linearization procedure, consid-
ering the two first tower modes, was presented by Kluger et al. (2016), who
computed the modes from the eigenvalue problem with an isolated tower.
Although more recent efforts on frequency-domain modelling have included
tower modes, platform flexibility has been little explored. As turbines grow
larger, support structure flexibility is expected to become more important,
and may need to also include the platform to properly capture the correct
global behaviour.
Due to the relatively complex aerodynamic calculations, and because the
design of the rotor usually is separated from the support structure design,
frequency-domain methods have typically focused on decoupling of the rotor
dynamics. Early work on FWTs considered only stochastic response due to
2.1. Global response analyses 19
kept unchanged. However, due to interaction with the controller and ro-
tor dynamics, the damping is frequency-dependent, and neglecting these
interaction effects leads to non-conservative damping estimates (Jonkman
2008). As damping is most important for the response around resonance,
individual values for each mode can be found from decay tests with steady
wind and an active controller (Jurado et al. 2018). Model tests have shown
that there also is a frequency-dependent inertia (or stiffness) effect from the
controller, which results in changes in the pitch natural frequency for the
system (Bachynski et al. 2016, Goupee et al. 2017). Souza and Bachynski
(2019) computed frequency-dependent inertia and damping coefficients for
the NREL 5MW turbine numerically, based on forced harmonic oscillations
of the nacelle, and were able to accurately predict pitch decay periods for
three different FWTs, using a simplified 2-DOF model with modified inertia
and damping matrices.
As linearized models assume small displacements, they are seen as unsuit-
able for extreme response. However, results from previous studies have
shown that both the tower base bending moment (Karimirad and Moan
2011) and global motions (Aggarwal et al. 2017) for spar FWTs can be quite
Gaussian in harsh environmental conditions, which suggests that a linear
model also may be used to assess the extreme response in early stages of
design, or reduce the computational burden by quickly identifying critical
load cases. The degree of nonlinearity depends on the considered response
parameter, and also varies with location, floater geometry, environmental
conditions, and design situations. Therefore, the results obtained with a
linear model cannot be generalized to any FWT structure, and such mod-
els must be continuously verified against more sophisticated tools during a
design process.
2.1.4 Structural reliability
In the design of FWTs, stochastic variables such as turbulence intensity,
wind-wave directionality, and fatigue strength are typically replaced by
design values, which are obtained from their characteristic values combined
with a given safety factor. In order to reduce conservatism and thus the
design costs caused by these safety factors, probabilistic analyses can be
performed to assess the safety level of the structure through the probabil-
ity of failure. Such analyses can also be used to evaluate the importance
of different model uncertainties on the resulting structural reliability, and
possibly calibrate appropriate safety factors, which has been the focus of
several studies for bottom-fixed wind turbines. Márquez-Domínguez and
Sørensen (2012) calibrated design fatigue factors (DFFs) for a 2.3 MW off-
2.2. Design optimization 21
shore wind turbine, and investigated the effect of inspections on the required
DFF values. Horn et al. (2019) quantified the effect of environmental load
uncertainties on the fatigue reliability of a 10 MW monopile, and derived
corresponding differences in DFFs resulting from the increased model fidel-
ity. Velarde et al. (2020) designed a 10 MW monopile with different DFFs
by varying the wall thickness, and used reliability analyses to identify the
necessary safety factor to achieve the target probability of failure during a
service life of 25 years without inspections.
Structural reliability problems consider the probability of failure for the
system, which can be expressed as
PF = P [g(x) ≤ 0] = fX (x) dx, (2.2)
g(x)≤0
eral methods have been suggested to resolve the issue, such as detuning the
controller gains (Larsen and Hanson 2007, Jonkman 2008), or to introduce
a feedback term proportional to the pitch velocity (Lackner 2009) or nacelle
velocity (van der Veen et al. 2012, Fleming et al. 2014) to manipulate the
rotor speed reference. For the latter approach, Fleming et al. (2019) sug-
gested that such controllers should not include the wave-frequency range
within their bandwidth, as doing so yielded an increase in the loads around
these frequencies. All methods lead to trade-offs between structural loads,
rotor speed tracking, and blade-pitch actuator use, and an integrated control
and structural design optimization approach is therefore needed to properly
evaluate and compare different solutions.
The shapes and natural frequencies of the tower modes, and thus the tower
response, are dependent on the platform design. In addition, as turbines
grow larger, the flexibility of the platform may become important, which
also can affect the response of the tower. In Fig. 2.5, the first bending mode
for a 10 MW spar is shown for both a rigid and flexible platform. In addition
to the modal shape being clearly affected by the platform flexibility, a 17 %
decrease in the natural frequency is observed. A similar effect was seen for
the OO-Star 10 MW semi-submersible, where the natural frequency of the
first tower mode was reduced by 25 % when the flexibility of the floater
was included (Müller et al. 2018). As the first bending mode of FWT
towers typically is located close to the 3P range, an accurate estimation of
the natural frequency is important to capture the correct modal excitation
forces. Souza and Bachynski (2020) studied the effect of flexible pontoons
for a 5 MW TLP wind turbine, and found that the flexible platform resulted
in a 12.7 % reduction in the weighted 1-h fatigue damage at the tower base,
due to the reduced natural frequency of the tower mode, which resulted in
less 3P excitation.
Commonly, wind turbine towers are designed with linearly tapered diameter
and wall thickness (Jonkman et al. 2009, Bak et al. 2013, Müller et al. 2018),
and are thus fully described by the values at the base and top. However,
the load picture for a FWT tower alters significantly from a land-based
or bottom-fixed offshore tower, with inertial and gravitational loads arising
from large platform motions, as well as possible restrictions at the tower base
due to compliance with the platform. An integrated design process which
considers the actual loads on the FWT system, and consistent reliability
levels, may therefore result in cost reductions.
Some effort has been made to perform integrated design optimization of the
blades and tower for onshore (Ning and Petch 2016) and bottom-fixed off-
28 Literature Survey
Figure 2.5: First support structure bending mode for a 10 MW spar with flexible
and rigid platform, and associated natural periods.
2.2. Design optimization 29
shore (Ashuri et al. 2014) wind turbines, which experience interaction effects
due to coupled eigenmodes, increased blade loads from tower shadow effects,
and the possibility of tower strike by the blade tips. It is expected that the
inclusion of blade design in an integrated optimization process will affect
the results also for FWTs. However, design optimization of wind turbine
blades is a separate body of literature, which requires high-fidelity aerody-
namic and structural formulations to properly capture the blade response,
and is therefore considered outside the scope of this thesis.
2.2.4 Formulation of objective and constraints
The success of a numerical optimization process depends heavily on the
formulation of objective function and constraints for the design problem.
Because the goodness of a design is measured purely through the value of
a scalar objective function, it is essential that the applied function reflects
a relevant performance measure. The constraints should ideally cover all
important limit states, in addition to any restrictions that may be necessary
to avoid non-physical or impractical solutions. However, the goal of an
optimization process is not necessarily to provide a ready-to-build design,
but may also be used to study the sensitivity of different input parameters,
or identify design trends that are non-intuitive for a human designer.
A common objective function in design optimization of offshore wind sub-
structures is the structural weight (Pasamontes et al. 2014, Chew et al.
2016, Oest et al. 2017), which is assumed to be closely related to the cost
of the design. The weight is a convenient cost parameter in an engineering
optimization problem, because it is easily calculated with small uncertainty.
However, the weight, or associated material costs, may not be a good meas-
ure of the actual cost optimality, and even in the case where only capital
costs are considered, manufacturing costs due to e.g. welding of joints should
be included (Muskulus and Schafhirt 2014). The manufacturing costs are
more uncertain, and are typically expressed as empirically fitted functions
of the structural dimensions. Cost models for manufacturing of general steel
structures are presented e.g. by Farkas and Jármai (2013); however, their
accuracy with respect to FWT structures have not been verified.
Operational costs are rarely considered in design optimization of offshore
wind turbines. For conventional offshore structures, the focus has been on
inspection planning, where a cost optimal schedule for crack inspection and
repair is sought. Madsen et al. (1989) presented an early formulation of an
optimization problem which minimized expected costs related to construc-
tion, inspection, repair, and failure, where a set of structural design para-
30 Literature Survey
Figure 2.6: Side view of the wind speed for a spar FWT without tilt (left) and
with a tilt angle of 10◦ (right). The black dots represent the vertical wake center
position (Wise and Bachynski 2020).
reduction in projected rotor area and consequently the resulting power pro-
duction. The use of a static tilt angle has computational benefits, because
it can be computed without the need of simulations, using the platform
restoring stiffness and the specified rated thrust for the turbine. Common
constraint values in the literature are 5◦ (Matha et al. 2015, Steinert et al.
2016) or 10◦ (Wayman et al. 2006, Tracy 2007, Hall et al. 2013, Karimi et al.
2017), which, given that the produced energy is proportional to the cosine of
the tilt angle to the power three (Matha et al. 2015), result in energy losses
of about 1 % and 4.5 %, respectively. Increased stability is likely to trans-
late into larger design costs, and the static tilt angle varies with the mean
wind speed. Additionally, Wise and Bachynski (2020) showed that larger
tilt reduced the wake loss for a downstream turbine due to vertical wake
deflection, as illustrated in Fig. 2.6. Concept-specific assessments of the
trade-off with energy loss, also considering farm level effects, are therefore
needed to find cost-optimal values for this constraint.
An important consideration in design is to maintain the structural integrity
of the power cable, which is designed for relatively small curvatures. The
maximum offset of the platform is often constrained in lieu of an explicit
dynamic model and design criterion for the power cable. As the relation
between offset and power cable loads is somewhat dependent on the design
of the specific system, no general offset criterion exists in design guidelines,
but limitations are typically expressed as a fraction of the water depth. The
maximum offset is also of interest in relation to the tendon connectors for
TLP wind turbines, where a maximum offset limitation of 10 % of the water
depth has been common industry practice for conventional TLPs (Bachynski
2014).
The offset is to a large degree a function of the mooring system design, which
in addition must possess sufficient breaking strength to withstand the design
32 Literature Survey
tension in the ULS and ALS (DNV GL 2018a). A minimum tension should
also be maintained at all times, as slack line events may result in high snap
loads when the line regains its tension. For drag anchors, which are not
designed to take vertical loads, a constraint may also be needed to prevent
uplift at the anchors for extreme offsets (Brommundt et al. 2012).
Mooring lines should be designed against fatigue failure; however, as the
design life of a mooring line often is shorter than that of the platform and
wind turbine, replacement of mooring lines during the design life of the
platform is an implicit design assumption (DNV GL 2018a). Consequently,
fatigue is not considered a key criterion in mooring system optimization
(Benassai et al. 2014). The effect of including a fatigue life constraint for
the mooring lines was demonstrated by Fylling and Berthelsen (2011), who
found that it resulted in a 20 % more expensive mooring system and about
5 % increase in spar buoy costs for the optimized system.
A challenge related to constraint formulations for FWTs is to derive expres-
sions based on actual design criteria for fatigue and ultimate loads. Often,
constraints are expressed in terms of response parameters which can indic-
ate the level of response, but are only applicable for relative comparison
of solutions. However, to achieve cost-effective structures with a consistent
reliability level, the constraints should reflect the actual limit state criteria.
This requires a large number of simulations, and often, no closed-form solu-
tion is available for the design loads. Steinert et al. (2016) attempted to
properly express ULS constraints, where the 50-year extreme response in
a single parked condition was considered. The constraints were based on
the assumption of Gaussian response, which simplifies the estimation of the
extreme value distribution considerably. Fatigue damage in the tower and
possibly the platform is of concern, but has rarely been assessed directly in
design optimization studies, where proxy variables such as tower top dis-
placement due to bending (Sandner et al. 2014, Lemmer et al. 2017) or
nacelle acceleration (Hall et al. 2013, Karimi et al. 2017) have been used
to assess the support structure loads. Although the large number of re-
quired load cases makes optimization for fatigue challenging, it is typically
the design-driving limit state for large parts of the support structure, and
accurate and consistent fatigue design will thus lead to more effective use
of material.
Another reason for constraining the nacelle acceleration has been the as-
sumption of strong correlation with loads on the drivetrain; however, Nejad
et al. (2019) showed that the accelerations at the tower top were not a good
indicator for either maximum or fatigue loads on the main bearings of the
2.2. Design optimization 33
Numerical Models
The results presented in Chapter 4 and the appended papers are based on
a linearized numerical model for global response analyses of floating wind
turbines, which was developed over the course of the thesis. The main
features of the model are summarized in the following. A brief description
of the evolution of the model, and the associated publications, is also given.
35
36 Numerical Models
while the system in Eq. (3.2) can be linearized by a Taylor series expansion:
∂f ∂f
ẋ ≈ f (x0 , u0 ) + Δx + Δu. (3.4)
∂x ∂u
If the operational point is chosen such that ẋ0 = f (x0 , u0 ) = 0, the linear-
ized dynamic system for the differential state and input variables can be
written in state-space form as (Chen 2013):
where A is the state matrix, and E is the input matrix. Given a set of system
outputs, y, which are described by y = g(x, u), a similar linearization can
be applied:
Δy = CΔx + DΔu. (3.6)
Here, C and D are commonly known as the output matrix and feedthrough
matrix, respectively. If the system is time-invariant, the matrices in Eq. (3.5)
and Eq. (3.6) are constant (Balchen et al. 2003).
loads are the thrust force (FT ), the tilting moment (MT ), and the aerody-
namic torque (QA ). In addition, distributed wave excitation forces (dFW )
are applied on the hull.
For each wind-wave condition, the linearization point is found from equilib-
rium when the system is subjected to the mean environmental loads. As the
model neglects current and mean wave loads, the wind forces on the rotor
and tower are the only load components with non-zero mean. The opera-
tional point for a given system design is thus solely a function of the mean
wind speed, which determines the references for the rotor speed, collective
blade pitch angle, and generator torque, as well as the equilibrium position
of the support structure and mooring system.
The model presented here considers three support structure degrees-of-
freedom (DOFs), namely surge, pitch, and the first bending mode. For
deep-draught spar platforms, the first order wave excitation in heave is
small. As a consequence, the heave response is not considered to have a
large effect on the structural loads on the system, and is thus not included in
the model. However, the natural period can be estimated, and subsequently
used in the design process to avoid heave resonance in the wave-frequency
range. Another issue relevant for spar platforms is the phenomenon known
as Mathieu instability, which can occur when the pitch restoring moment
varies harmonically due to large heave motions (Koo et al. 2004). Haslum
and Faltinsen (1999) reported that unstable solutions occur for certain ratios
between the frequency of the heave motion and the pitch natural frequency,
and can thus be prevented assuming that the heave response is dominated
by resonant motions.
In addition to the support structure DOFs described above, a rotor speed
DOF is included. The inputs to the structural system consist of both control
system outputs and disturbances due to wind and wave loads. The control
outputs, i.e. the generator torque and the collective blade pitch angle, are
described by a separate control system model, which can be combined with
the structural model to create a single closed-loop system.
The model considers coupled response of four disciplines: structural dy-
namics, hydrodynamics, aerodynamics, and control dynamics. The applied
formulations for each of the disciplines are detailed in the following subsec-
tions.
38 Numerical Models
MT QA
FT
z y
Line 2
dFW Wind
x Line 1 x
Waves
Line 3
n
w(z, t) = ψk (z)χk (t) = ψ(z) χ(t). (3.7)
k=0
The weight functions are assumed harmonic, i.e. χ(t) = χ0 eiωt , and are
found by solving the generalized equations of motion. Using the principle
of virtual work, the generalized system matrices can be established as a
combination of continuous and discrete terms. The support structure is
modelled as a slender beam undergoing pure bending deformations, and
3.3. Structural dynamics 39
z z z
l
Bχ = b(z)ψ(z)ψ (z)dz + Bp ψ(zp )ψ (zp ), (3.9)
−d p
and
l
Kχ = Kp ψ(zp )ψ (zp ) + EI(z)ψ ,zz (z)ψ
,zz (z)dz
p −d
l (3.10)
+ N (z)ψ ,z (z)ψ
,z (z)dz,
−d
where dF are distributed forces, and Fp and Qp are discrete forces and mo-
ments, respectively. The linear system of equations for the weight functions
then becomes
Mχ χ̈(t) + Bχ χ̇(t) + Kχ χ(t) = F(t). (3.12)
ID ϕ̈ = QA − Ngear QG , (3.13)
where ID is the drivetrain inertia, Ngear is the gear ratio, QG is the high-
speed shaft (HSS) generator torque, and ϕ̇ is the low-speed shaft (LSS)
rotational speed.
With the structural DOFs described above, the state vector is defined as
x = χ χ̇ ϕ̇ . (3.14)
The state matrix in Eq. (3.5) can now be derived from the total mass,
damping, and stiffness matrices for the system (Chen 2013):
0 I
A= , (3.15)
−M−1 K −M−1 B
while the input matrix E relates the system inputs to the external forces
acting on the states.
The mean configuration of the mooring system, as well as the linearized
mooring stiffness, for a given mean thrust is found using the elastic catenary
equations. For the dynamic response of the mooring lines, the simplified
dynamic model proposed by Larsen and Sandvik (1990), and later extended
by Lie and Sødahl (1993), is utilized. Here, each line is modelled as a single-
DOF spring/damper system, that consider inertia and drag forces due to
3.4. Hydrodynamics 41
TD
x
kE
cL
u
kG
Figure 3.3: Simplified dynamic mooring line model as illustrated by Lie and
Sødahl (1993).
the line motion in addition to elastic and geometric stiffness. Given the
displacement of the upper end of the line in the tangential direction, x(t),
the dynamic tension in the line can then be calculated from
where u(t) is the generalized displacement of the line, found from a quasi-
static analysis where the upper end is displaced two times the standard
deviation of the platform motion from the equilibrium position (signific-
ant motion). The geometric stiffness, kG , is found as the secant stiffness
when the line is moved from equilibrium to the position of significant mo-
tion. m∗ and c∗ are the generalized inertia and damping coefficients due to
Morison-type hydrodynamic forces, respectively, which are found using the
quasi-static line configuration as the shape function. kE is the elastic line
stiffness. The model is illustrated in Fig. 3.3.
3.4 Hydrodynamics
Assuming a circular cross section for the hull, the first order wave excit-
ation forces can be calculated using the analytic expression developed by
MacCamy and Fuchs (1954). The force per unit length, dFW , for a regular
wave with unit amplitude is then given by
1 J1 (ka)
G= , tan α = . (3.18)
(J1 (ka))2 + (Y1 (ka))2 Y1 (ka)
Here, Jn and Yn are the derivatives of the Bessel function of the first and
second kind, respectively, of order n, and a is the hull radius. The assump-
tion of a circular cross section also allows for simplification of the radiation
forces. The transverse added mass per unit length can be approximated
using the analytical 2D coefficient for a circular cylinder in infinite fluid
(Newman 1977):
a11 (z) = ρπD2 (z)/4, (3.19)
where D is the hull diameter. For a relatively slender structure, radiation
damping will be small and can be neglected, especially at the relevant nat-
ural frequencies.
The calculation of the natural period in heave requires an estimate of the
added mass in the vertical direction. Here, the vertical added mass is ap-
proximated as the value for a 3D circular disc (DNV GL 2017), with the
same diameter as the bottom of the platform (Db ), which was found to
yield good agreement with numerical linear potential flow solutions for sev-
eral considered spar designs:
The viscous drag forces on the hull are nonlinear, and must therefore be
linearized. While viscous wave excitation is small compared to the wind
excitation forces and thus not included in the model, viscous damping is
important for the low-frequency surge motions. Neglecting the wave particle
velocities, the viscous damping due to the velocity of the structure (ẋ) can
be added based on stochastic linearization of the quadratic drag term in
Morison’s equation (Borgman 1969):
8
Bvisc |ẋ|ẋ ≈ Bvisc σẋ ẋ = Bvisc,lin ẋ, (3.21)
π
where Bvisc is the dimensional quadratic damping coefficient, and Bvisc,lin is
the stochastically linearized coefficient. The standard deviation of the velo-
city, σẋ , depends on the damping forces and vice versa, and must therefore
be found using an iteration scheme.
3.5. Aerodynamics 43
3.5 Aerodynamics
Aerodynamic loads in global response analyses of FWTs are typically com-
puted using BEM theory, where induction factors, and subsequently blade
element forces, are found iteratively (Burton et al. 2011). Assuming that
the rotor design does not change during the course of the design process,
induction factors can be pre-calculated and stored in look-up tables, before
they are used in the response analyses. For each blade element, the induc-
tion factors are functions of the blade pitch angle (θ) and tip-speed ratio
(λ), where the latter is defined as
ϕ̇R
λ= . (3.22)
v
Here, R is the rotor radius, and v is the wind speed. In the present model,
the induction factors are computed from the BEM equations following the
procedure presented by Ning (2013), which includes Prandtl hub and tip
loss factors, and a modified Glauert correction as described by Buhl (2005)
for high induction factors. A quasi-static approach is used, meaning that
dynamic wake and dynamic stall effects are neglected. The calculated in-
duction factors at midspan for the DTU 10 MW blade design are shown in
Fig. 3.4.
Following quasi-steady BEM theory, the aerodynamic forces on a blade ele-
ment are nonlinear functions of (relative) wind speed, rotor speed, and
blade pitch angle, which can be found once the induction factors are known.
Linearizing the equations, the normal and tangential blade element forces
relative to the rotor plane can be expressed as:
The blade root loads of interest, i.e. the flapwise shear force Fy , the flapwise
bending moment Mz , and the edgewise bending moment My , are found by
44 Numerical Models
Figure 3.4: Calculated induction factors at midspan for the DTU 10 MW blade.
R
Fy = Fn dr, (3.24a)
0
R
Mz = r Fn dr, (3.24b)
0
R
My = r Ft dr. (3.24c)
0
Here, the subscripts point to the blade coordinate system, which is defined
as in Fig. 3.5. The resultant rotor loads can then be calculated as:
3.6. Control system description 45
3
FT = Fy,i , (3.25a)
i=1
3
MT = Mz,i cos ϕi , (3.25b)
i=1
3
QA = My,i . (3.25c)
i=1
As only the resultant loads on the rotor are of interest in the global response
model, it is desirable to simplify the wind field and express it using a single
scalar variable. This is achieved by expressing a rotor effective wind speed
for each of the resultant rotor loads in Eq. (3.25). The rotor effective wind
speed is a spatially constant wind speed which yields identical resultant
loads as the full wind field, and can be established from the incoming wind
spectrum, the spatial coherence function, and the transfer functions between
wind speed and blade element loads from Eq. (3.23). A detailed derivation
of the rotor effective wind speed, based on the thesis by Halfpenny (1998),
is provided in Appendix B.
In addition to the rotor loads, the aerodynamic drag forces on the tower are
taken into account, using a Morison drag formulation. The tower loads are
expressed as a combination of a mean force, which is added directly, and
a frequency-dependent force, which is found using stochastic linearization.
Here, only the excitation forces are considered, meaning that the tower drag
forces arising from the movement of the turbine are neglected.
torque is set proportional to the squared rotor speed, to track the optimal
tip-speed ratio and thus optimize the power output:
2
QG = Kg ϕ̇ , (3.27)
where Kg is the generator torque constant. The change in generator torque
for a rotor speed error Δϕ̇ then becomes
∂QG
ΔQG = Δϕ̇ = 2Kg ϕ̇0 Δϕ̇, (3.28)
∂ ϕ̇
where ϕ̇0 is the rotor speed reference. Alternatively, the generator torque
can be regulated using a proportional-integral (PI) controller:
t
ΔQG = kp,Q Δϕ̇ + ki,Q Δϕ̇ dτ, (3.29)
0
where kp,Q and ki,Q i are the proportional and integral gains for the rotor
speed error feedback, respectively. Above rated wind speed, the blade pitch
angle is found using a gain-scheduled PI controller:
t
Δθ = ηk kp,θ Δϕ̇ + ηk ki,θ Δϕ̇ dτ, (3.30)
0
where kp,θ and ki,θ are the proportional and integral gains, respectively,
and ηk is the gain-scheduling parameter. For the torque controller, two
common strategies exist for the above-rated regime: (1) constant power or
(2) constant torque. The variation in generator torque then becomes
⎧
⎪
⎨− N Pr ϕ̇2 Δϕ̇, constant power,
gear r
ΔQG = (3.31)
⎪
⎩ 0, constant torque,
where Pr and ϕ̇r are the rated power and rotor speed, respectively.
To avoid the potential problem of negative damping for the platform motions
above rated wind speed, a feedback term proportional to the platform pitch
velocity or nacelle velocity may also be included, to manipulate the rotor
speed reference. Here, a modified rotor speed reference, ϕ̇0 , as defined by
Lackner (2009) is used:
ϕ̇0 = ϕ̇0 (1 + kf ẋf ), (3.32)
where kf is the velocity feedback gain, and ẋf is the nacelle or pitch velocity.
An updated expression for the rotor speed error can then be established as
Δϕ̇ = ϕ̇ − ϕ̇0 = Δϕ̇ − ϕ̇0 kf ẋf . (3.33)
Here, ẋf may also be passed through a low-pass filter before it is fed back
to the blade-pitch controller.
3.7. Response to stochastic input 47
The cross spectral density matrix of the response vector y can be calculated
from
Sy (ω) = H(ω)Su (ω)H(ω)H , (3.35)
where (·)H denotes the conjugate transpose (Naess and Moan 2013). The
response spectra of y are then found along the diagonal of Sy (ω).
The input vector is expressed as
u = vFT vMT vQA F
W , (3.36)
where vFT , vMT and vQA are rotor effective wind speeds for thrust, tilting
moment and aerodynamic torque, respectively. The wave load vector, FW ,
contains the generalized wave excitation force for each support structure
DOF. Su (ω) is the cross spectral density matrix for the load process, which
has the following structure:
S (ω) 0
Su (ω) = wind , (3.37)
0 Swave (ω)
by SN curves, where the relation between the stress range (S) and number
of cycles to failure (N ) is given as (DNV GL 2019c):
N = KS −m , (3.40)
where K and m are material parameters dependent on the type of fatigue
detail. Using the Palmgren-Miner hypothesis of linear accumulation of dam-
age, the fatigue damage in time T can be expressed as (Naess and Moan
2013): ∞
νT m
D= s fS (s)ds. (3.41)
0 K
Here, the fatigue loading is described by the cycle rate, ν, and the PDF of
the stress ranges, fS (s). If the response process is narrow-banded and Gaus-
sian, the stress ranges can be described by a Rayleigh distribution; however,
because the response of FWTs in general is wide-banded, the narrow-band
formulation may result in overly conservative damage estimates. For a gen-
eral wide-banded response, cycle counting methods based on the stress time
series are typically applied, with rainflow counting being considered the
most accurate for estimation of fatigue damage (Naess and Moan 2013).
This method has no closed-form solution in the frequency domain for wide-
banded Gaussian response, and given a stress response spectrum, the stress
range distribution must thus be estimated either by transformation of the
response spectra to the time domain, or directly in the frequency domain us-
ing empirical formulae. For the latter approach, various methods have been
proposed, where the formulae presented by Dirlik (1985) and Benasciutti
and Tovo (2005) have been shown to yield accurate approximations of the
rainflow counting method for a variety of spectral shapes and bandwidths
(Gao and Moan 2008).
The Dirlik method, which is implemented in the current model, uses the
stress response spectrum and empirical factors to fit the PDF of the stress
cycles to a combination of an exponential and two Rayleigh distributions,
and the accuracy of the method is therefore dependent on how well the
rainflow-count can be represented by these distributions. The resulting
stress range distribution is given as:
2
2
G1
Q exp − Q
Z
+ G2 Z
R2
exp − 2R
Z
2 + G3 Z exp − Z2
fS (s) = √ , (3.42)
2 m0
where G1 , G2 , G3 , R, and Q are empirical weight factors found from the
spectral moments, and Z is the normalized stress range:
S
Z= √ . (3.43)
2 m0
3.8. Optimization framework 49
where νξ+ is the mean upcrossing rate of the level ξ. Assuming a stationary
Gaussian process with zero mean, the CDF can be written (Naess and Moan
2013)
ξ2
FΞ (ξ) = exp −ν0 T exp − 2
+
, (3.45)
2σy
where ν0+ is the mean zero-upcrossing rate, which can be found from the
zeroth and second order spectral moments:
1 m2
ν0+ = . (3.46)
2π m0
The p-fractile of the distribution, i.e. the value which is exceeded with prob-
ability (1 - p) during the time T , becomes (Naess and Moan 2013):
ν0+ T
ξp = σ 2 ln , (3.47)
ln(1/p)
while the most probable maximum value, ξc , also known as the characteristic
largest extreme, can be found approximately from the relation FΞ (ξc ) = e−1 :
ξc = σ 2 ln ν0+ T . (3.48)
∂R dy ∂R
=− . (3.51)
∂y dx ∂x
Given the partial derivatives of R with respect to x and y, dy/dx can then
be found by solving the linear system in Eq. (3.51). This is known as the
direct (forward) method, where the linear system must be solved once for
each model input (Martins and Ning 2020).
Alternatively, dy/dx in Eq. (3.49) can be replaced by the expression in
Eq. (3.51):
df ∂F ∂F ∂R −1 ∂R
= − . (3.52)
dx ∂x ∂y ∂y ∂x
ψ
Here, ψ is known as the adjoint vector, which can be described through the
following relation (Martins and Ning 2020):
∂R ∂F
ψ= . (3.53)
∂y ∂y
Now, ψ can be found by solving the linear system in Eq. (3.53), and used
together with Eq. (3.52) to calculate the total derivative. This is known
as the adjoint (reverse) method. The linear system in Eq. (3.53) must,
contrary to the forward method, be solved once for each model output,
and this method therefore offers increased computational efficiency when
the number of model inputs is larger than the number of model outputs.
OpenMDAO supports both direct and adjoint gradient computations, and
the method is chosen based on the specific problem structure.
To illustrate the extent of the analytic gradient implementation, the required
calculations are shown for a single model component. The component which
calculates the uncoupled heave natural period (Tn3 ), given the structural
mass of the system (M ), the added mass in heave (A33 ), and the heave
restoring stiffness (C33 ), uses the following relation:
M + A33
Tn3 = 2π . (3.54)
C33
52 Numerical Models
∂Tn3 π
= , (3.55b)
∂A33 C33 MC+A 33
33
M +A33
∂Tn3 π C33
=− . (3.55c)
∂C33 C33
As both M , A33 , and C33 are dependent on the design variables for the
optimization problem, the chain rule in Eq. (3.49) is then needed to compute
the total derivative of Tn3 with respect to an arbitrary design variable p:
dTn3 ∂Tn3 ∂Tn3 dM ∂Tn3 dA33 ∂Tn3 dC33
= + + + . (3.56)
dp ∂p ∂M dp ∂A33 dp ∂C33 dp
If there is an implicit relationship between M , A33 , or C33 , and p, obtaining
the total derivative also requires solving the linear system in Eq. (3.51) or
Eq. (3.53).
The component-level derivatives must be provided for all parts of model, and
may result in rather involved computations, especially for more advanced
components. The complexity and required effort related to the gradient
implementation are thus comparable to the development of the analysis
model itself, and makes up about 70 % of the total code volume in the
present work.
The model layout and data dependencies are illustrated using an extended
design structure matrix (XDSM) (Lambe and Martins 2012) in Fig. 3.7.
Here, the rounded blue box is the optimizer, green boxes are analysis models,
gray parallelograms are data, white parallelograms are fixed input paramet-
ers, and grey lines represent data dependencies. The data moves from top
to bottom and left to right in the upper triangular part, and from bottom
to top and right to left in the lower triangular part.
Mean wind speed,
Corner freq., SN data, Buckling resistance,
Mean wind speed Water depth Water depth RNA data sig. wave height, Drag coefficient
gain schedule DFFs pitch and offset limits
spectral peak period
Wave excitation
Wave excitation
transfer function
Control
Control matrices
model
Environmental
Load matrix
loads
Upcrossing
Extreme response
analysis
Cost
System costs
models
Figure 3.7: XDSM diagram showing the model layout and data dependencies.
3.8. Optimization framework
53
54 Numerical Models
V1, applied in P1: The structural model for the support structure using
generalized coordinates is established, but assumes a rigid platform, while
the elasticity of the tower is considered. The hydrodynamic formulations
presented in Section 3.4 are implemented, and the quasi-static mooring sys-
tem response is found from the elastic catenary equations. The rotor DOF
is not included, and aerodynamic excitation (thrust and tilting moment) is
computed from time-domain simulations with turbulent wind and an active
controller. The aerodynamic damping is found as the change in thrust due
to a small change in wind speed, assuming no reaction from the rotor or
controller. The Dirlik method for fatigue damage in the support structure
and extreme response based on upcrossing analysis are implemented.
V2, applied in P3: Building off of V1, A single-DOF model for the rotation
of the shaft is implemented, and a flexible platform is included. Linearized
aerodynamic loads are derived as functions of wind speed, rotor speed, and
blade-pitch angle, using quasi-steady BEM theory and a rotor effective wind
speed. Generator-torque and blade-pitch controllers are added to the model.
Dynamic mooring line tension is found based on a single-DOF model which
consider drag and inertia forces due to the motion of the line. Analytic
gradients are calculated for the model, which is implemented in the Open-
MDAO framework.
V3, applied in P4: More advanced control strategies for the blade-pitch
controller are implemented. The new controllers consider velocity feedback
using either the nacelle or the platform pitch velocity, with an optional low-
pass filter on the velocity signal.
Research Findings
This chapter presents the main research findings from the papers included
in the thesis, and discusses their relation to the contributions stated in
Section 1.3.
55
56 Research Findings
Here, FT is the thrust force, v is the wind speed, kP and kI are the propor-
tional and integral gains of the blade-pitch controller, and θ is the blade-
pitch angle. C̃ϕ is the complex transfer function between x and ϕ, i.e.,
The derived coefficients in Eq. (4.1) were compared to results from nu-
merical simulations using forced oscillations with AeroDyn as described by
Souza and Bachynski (2019), as well as constant damping coefficients, which
neglect the controller and rotor dynamics. In Fig. 4.1, the coefficients are
shown for the low-frequency range, for a mean wind speed of 25 m/s. The
overall agreement between the analytical and numerical solutions was good,
while the constant damping coefficient was found to significantly overes-
timate the aerodynamic damping at low frequencies. The magnitude of the
aerodynamic inertia is low compared to the total mass of typical FWTs, res-
ulting in negligible consequences for the surge response. Due to the height
of the nacelle, however, the contribution to the moment of inertia in pitch
is considerable.
To assess the importance of frequency-dependent coefficients, a frequency-
domain model of the OO-Star 10 MW semi-submersible, which was part of
the the LIFES50+ project (Müller et al. 2018), was created. Results from
the frequency-domain model, using different formulations for the aerody-
namic damping and inertia, are compared to nonlinear time-domain simu-
lations using SIMA. The response spectra for platform pitch and tower base
4.1. Linearized dynamic analyses 57
Figure 4.1: Aerodynamic inertia and damping coefficients, 25 m/s mean wind
speed.
bending moment are shown in Fig. 4.2, and relative errors in the standard
deviations are presented in Table 4.1.
While the constant damping coefficient underestimated the responses, the
opposite was observed for the frequency-dependent coefficients. There were
also notable differences between the analytical and AeroDyn models, despite
the good agreement for the damping coefficient, which suggested that the
response was sensitive to resonant motions and consequently the damping
level in the model. Although the actual numbers also depend on the con-
sidered wind-wave conditions, platform concept, and additional damping in
the model, the results demonstrated that the use of constant damping coef-
ficients is non-conservative, and that this can be avoided in a simple manner
through the frequency-dependent coefficients in Eq. (4.1).
The deviations in response between the frequency-domain models in Fig. 4.2
can be attributed to differences in aerodynamic damping coefficients. The
effect of aerodynamic inertia was much less prominent, but can be seen in
the calculated platform pitch natural period. Calculated natural periods in
pitch, with and without aerodynamic inertia, are shown in Fig. 4.3.
58 Research Findings
Figure 4.3: Pitch natural period with and without aerodynamic inertia effect.
4.1. Linearized dynamic analyses 59
Figure 4.4: Skewness and kurtosis values for tower fore-aft bending moment.
4.1. Linearized dynamic analyses 61
Figure 4.5: 1-h weighted tower base fatigue damage arranged by wind speed.
Figure 4.6: Linearized tower base bending moment response spectrum for U = 15
m/s, Hs = 3 m, and Tp = 12 s, where the natural frequency of the first bending
mode is artificially increased in steps of 0.06 rad/s (≈ 0.01 Hz) by modifying the
modal stiffness. The corresponding increase in calculated 1-h fatigue damage is
shown, together with the 3P frequency at rated speed (black line).
For the nonlinear time-domain simulations, a function for the tail of the
mean upcrossing rates can be established empirically, as described in P1.
The nonlinear and linear average upcrossing rates (AUR) for surge, pitch,
and tower base bending moment are shown for two 50-year conditions in 4.7.
The plots show the multiplication factor κ, which is related to the extreme
response, y, through the relation
y = μ + κσ, (4.4)
where μ and σ are the mean and standard deviation of the process.
Very good agreement was obtained with an extreme wind speed in the
parked condition, where the response was governed by wave loads. Lar-
ger nonlinearities are observed near the rated wind speed, due to the large
thrust forces in this condition. This was especially the case for the pitch
response. The tower base bending moment, on the other hand, was also
quite Gaussian in the operational condition, due to the large influence of
inertial forces arising from wave-frequency motions.
4.1. Linearized dynamic analyses 63
(a) Surge.
(b) Pitch.
Figure 4.7: Upcrossing rates for 11 m/s (left) and 50 m/s (right) mean wind
speed.
64 Research Findings
While these results suggest that a linear model also can be used to evaluate
extreme response in early phases of design, especially for the bending mo-
ments in the support structure, such a conclusion cannot be drawn in the
general case. The accuracy of the linear model depends on the importance
of nonlinearities in the system, which will vary with location, floater geo-
metry, and different operational conditions. The agreement with nonlinear
analyses also depends on how the system is represented in the time-domain
model, and care should be taken to make sure that relevant nonlinearities
are included.
Barrier
Figure 4.9: Optimized blade-pitch controller gains plotted against resulting rotor
speed variation.
Figure 4.10: Trade-offs between rotor speed variation and system costs in the
multiobjective optimization.
where the diameter increased from the base. The importance of integrated
design optimization was demonstrated here, as simultaneous design of the
platform and tower is needed to identify the optimal solution near the in-
terface.
Tower shape
While typical wind turbine towers are linearly tapered, the integrated op-
timization procedure employed in the present work suggested a different and
somewhat unconventional tower design. The reduction in costs achieved by
adopting such a design was quantified by optimizing the tower using three
different strategies. In all strategies, 20 years fatigue life and global buck-
ling in selected 50-year conditions were considered, while a 15◦ pitch angle
constraint was applied on the platform.
In T1, the tower and platform were optimized in an integrated fashion, as
described in previous sections. The same procedure was used in T2, but
here, the tower was constrained to a linearly tapered shape. Finally, T3
considered a procedure where only the tower was optimized with a linearly
tapered shape using a fixed platform design, taken as the optimal design
from T1. The main consequence is then that the tower base diameter is
fixed in the optimization, such that it matches the platform. The resulting
optimized tower designs are shown in Fig. 4.12.
The resulting normalized tower and platform costs are shown in Table 4.2.
The consequence of constraining the tower shape depends on the optimiz-
68 Research Findings
Figure 4.11: Optimized tower and platform design with fatigue utilization (left)
and buckling utilization (right). The wall thickness is scaled by a factor of 40
relative to the diameter for illustration purposes.
4.2. Integrated design optimization 69
Figure 4.12: Optimized tower designs using three different design strategies,
with fatigue utilization (left) and buckling utilization (right). The wall thickness
is scaled by a factor of 40 relative to the diameter for illustration purposes.
70 Research Findings
Table 4.2: Normalized platform and tower costs using different design strategies.
(a) Weighted average rotor speed variation. (b) Long-term tower base fatigue damage.
Strategy Description
CS1 PI controller
CS2 PI controller + platform pitch velocity feedback
CS3 PI controller + nacelle velocity feedback
CS4 PI controller + nacelle velocity feedback + low-pass filter
74 Research Findings
(a) Costs relative to CS1. (b) Rotor speed variation and ADC values.
Figure 4.15: Objective and constraint function values for optimized designs.
(a) System reliability index during service (b) Combined design and inspection costs
life, including the effect of inspections. comparison, r=0.05.
Figure 4.17: Reliability and total costs with different DFFs and inspection sched-
ules.
Conclusions and
Recommendations for Future
Work
5.1 Conclusions
This thesis has sought to improve the design process for FWTs through
the development of efficient methods for global dynamic response analyses,
and numerical design optimization techniques. The main focus has been
on the support structure for 10 MW spar-type turbines, considering fatigue
and ultimate loads. Some attention has also been given to the blade-pitch
controller and mooring system, particularly related to coupling effects with
the support structure.
A linearized aero-hydro-servo-elastic model with four structural DOFs was
shown to yield acceptable accuracy for the support structure response with
co-directional wind and waves. Nonlinear effects resulted in poor perform-
ance of the linearized model when the tower bending natural frequency was
close to the 3P range. Some overestimation was also observed for the res-
onant platform pitch response, which caused the model to be conservative
in general compared to nonlinear analyses. Analytical frequency-dependent
aerodynamic damping and inertia coefficients were derived, which removed
79
80 Conclusions and Recommendations for Future Work
83
84 BIBLIOGRAPHY
DNV GL. Loads and site conditions for wind turbines. Technical Report
DNVGL-ST-0437, DNV GL, 2016.
DNV GL. Support structures for wind turbines. Technical Report DNVGL-
ST-0126, DNV GL, 2018b.
M. Muskulus. Simplified rotor load models and fatigue damage estimates for
offshore wind turbines. Philosophical Transactions of the Royal Society
A, 373(20140347), 2015. doi: 10.1098/rsta.2014.0347.
The Carbon Trust. Floating offshore wind: market and technology review.
Technical report, The Carbon Trust, 2015.
J. van der Tempel. Design of support structures for offshore wind turbines.
PhD thesis, TU Delft, 2006.
Appended Papers
99
100 Appended Papers
Paper 1
101
102 Appended Papers
0DULQH6WUXFWXUHV ²
Marine Structures
journal homepage: www.elsevier.com/locate/marstruc
AR T IC L E I NFO A B S T R A C T
Keywords: A linear model for efficient design evaluation of spar floating wind turbines is presented, and
Offshore wind verified against a nonlinear time domain model with regards to long-term fatigue and short-term
Design analysis extreme response for two different spar designs. The model uses generalized displacements and a
Frequency domain semi-analytical approach to establish the equations of motion for the system, which are solved in
Fatigue
the frequency domain. The results show agreement within ± 30% for the long-term fatigue
Extreme response
considering operational conditions, however, the linear fatigue damage estimates are sensitive to
the accuracy of the estimated natural frequency of the first bending mode. The results also
suggest that a small number of environmental conditions can be simulated with a nonlinear time
domain model to verify and possibly tune the linear model, which then can be used to run the full
long-term analysis. Short-term extreme tower base bending moments and surge and pitch mo-
tions are observed to be nearly Gaussian above cut-out wind speed, as the response is dominated
by wave forces. Consequently, the linear model is able to accurately capture the upcrossing rates,
which are used to calculate the characteristic largest extreme response. For an operational case
near rated wind speed, the response is somewhat non-Gaussian, which gives larger discrepancies
between the linear and nonlinear models. However, due to large mean values in this condition,
the total error in the extreme response is reduced, and reasonable agreement is achieved.
1. Introduction
The offshore wind industry has had significant growth over the last decade, and as a majority of the global wind resources are
located in deeper waters, there has recently been an increased interest in floating wind turbines (FWTs). For floating wind farms to be
economically feasible, cost-effective and reliable designs are needed. Design optimization of FWTs is a complex task, and due to
interactions between aerodynamics, hydrodynamics, structural dynamics and control, coupled nonlinear time domain (TD) simu-
lations are usually applied. As numerous load cases need to be analysed, the design process becomes computationally very expensive.
It is therefore desirable to use simplified models, especially in preliminary design [1]. Frequency domain (FD) models provide an
efficient way of performing dynamic simulations and are frequently used for floating structures in the offshore industry, usually in
fatigue analyses [2]. Although less common in the offshore wind industry, design standards state that frequency domain analyses also
can be used to calculate the fatigue loads for FWTs, if validation against time domain analyses or experiments is performed [3].
Several studies have investigated the applicability of FD models for FWTs. Typically, these models have considered first order
hydrodynamic loads from potential theory and steady aerodynamic loads, and have been used to calculate response amplitude
operators (RAOs) for the rigid body motions of the platform [4–7].
∗
Corresponding author.
E-mail address: [email protected] (J.M. Hegseth).
https://doi.org/10.1016/j.marstruc.2018.10.015
Received 10 August 2018; Received in revised form 2 October 2018; Accepted 30 October 2018
(OVHYLHU/WG$OOULJKWVUHVHUYHG
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Although useful information about the system may be found using only rigid body motions, the elasticity of the tower has been shown
to significantly affect the global response of FWTs. Bachynski and Moan [8] compared a linear FD model with three rigid body modes
to a nonlinear TD model for different tension leg platform wind turbine designs in both wave-only and combined wind-wave con-
ditions. The study was later expanded to include turbulent wind excitation and aerodynamic damping [9], and concluded that the
linear model was insufficient for design calculations, partly due to the rigid modelling of the tower. Kvittem and Moan [10] used a
similar procedure to calculate tower base bending moments for a semi-submersible wind turbine, where the first bending mode of the
tower was found from a free decay test and included using generalized coordinates. The contribution from the first bending mode on
the tower base bending moment response was expressed as a dynamic amplification factor. The model was found to perform rea-
sonably well in the 13 considered load cases, however, the fatigue damage was underestimated by up to 60% due to its exponential
nature. A somewhat different approach was applied by Kluger et al. [11], who used statistical linearization to develop a FD model for
the OC3 spar wind turbine [12]. The two first elastic tower modes were found for a fixed foundation and included in the analysis, and
equivalent fatigue stress (EFS) at the tower base due to wave excitation and steady wind was calculated for 11 different environ-
mental conditions. The results were compared to values from TD simulations reported by Matha [13], where the FD analysis was
found to underestimate the total EFS by 12%.
Although earlier work has assessed the fatigue damage based on FD analyses for FWTs, these calculations have only been per-
formed for a limited number of environmental conditions. Thus, little information exists regarding the accuracy of such models for a
full long-term analysis. In addition, the design process must also consider the extreme response of the system, which can be a
comprehensive task. Results from previous studies have shown that both the tower base bending moment [14] and global motions
[15] for spar FWTs can be quite Gaussian in harsh environmental conditions, which suggests that a linear model also may be used to
assess the extreme response in early stages of design. For a spar FWT, in addition to the ultimate stresses in the structure, the extreme
surge and pitch response of the platform may be of interest, as these motions are important for loads on the mooring system and
nacelle components, respectively.
In the present work, a semi-analytical FD approach is presented, where generalized degrees-of-freedom (DOFs) are used to
describe the dynamic behaviour of the FWT. The system is linearized, and due to the simple geometry of the hull, closed-form
expressions may be used for the hydrodynamic loads with good accuracy. This removes the need for a separate numerical hydro-
dynamic analysis, which significantly increases the computational efficiency. The aerodynamic loads are found numerically from TD
simulations, and the forces are superimposed to find the total response. The approach is then compared to fully coupled nonlinear TD
simulations for two different 10 MW spar FWT designs, considering long-term fatigue damage at the tower base, as well as short-term
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Fig. 1. Spar designs considered in the present work. Spar 1 (left) and spar 2 (right).
extreme values for the tower base bending moment and surge and pitch motions of the platform in selected conditions along the 50-
year environmental contour. The linear FD model presented in the present work is not applicable for detailed analyses, but may
complement state-of-the-art tools in early stages of design. The model can help the designer get a quick overview of the response of
the system, and give indications of how the system responds to changes in the design parameters. It may also be used to identify
critical load cases, which then can be analysed using coupled nonlinear TD analyses.
2. System description
The two spar buoys considered in the current work are shown in Fig. 1. The designs are based on the OC3-Hywind design [12] and
consist of two cylinders with different diameters, connected by a linearly tapered section. The hull extends to a height of 10 m above
still water level (SWL), where it is fixed to the tower. Spar 1 has the same draft as the original OC3 spar of 120 m, but with increased
diameters to provide enough buoyancy to support a 10 MW wind turbine, and to match the tower base diameter. The large draft of
the spar provides good stability and hydrodynamic performance, however, it limits the use to deep water, and complicates con-
struction and towing to site. The draft of the second design (spar 2) was therefore reduced by 25%, which makes it more suitable for
intermediate water depths. In addition to the draft reduction, the lower hull diameter was increased by 25% compared to spar 1, to
achieve a sufficient amount of buoyancy.
The two designs are presented in Table 1, where values for the vertical centre of gravity (CoG) and moment of inertia are
calculated by assuming that the hulls are partially filled with concrete ballast to achieve the correct draft, using a ballast density of
2600 kg/m3, while a constant wall thickness of 6 cm is assumed for the steel hull. The water depth for both platforms is 320 m.
Table 1
Platform properties.
Parameter Spar 1 Spar 2
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Table 2
Main characteristics of the DTU 10 MW reference wind turbine [16].
Parameter Value
Rated power 10 MW
Rotor orientation and configuration Upwind, three blades
Rotor, hub diameter 178.3 m, 5.6 m
Hub height 119.0 m
Cut-in, rated, cut-out wind speed 4.0 m/s, 11.4 m/s, 25.0 m/s
Cut-in, rated rotor speed 6.0 rpm, 9.6 rpm
Overhang, shaft tilt, pre-cone 7.1 m, 5.0°, −2.5°
Rotor, nacelle, tower massa 230.7 t, 446.0 t, 628.4 t
a
Mass for original land-based tower.
The spar buoys support the DTU 10 MW reference wind turbine [16], with the tower shortened by 10 m to achieve the same hub
height as the onshore turbine at 119 m above SWL. The tower has a linearly tapered outer diameter and consists of ten sections with
constant wall thickness in each section. In the modified tower, each section is shortened by 1 m, while the top and bottom diameters
are kept unchanged. The key characteristics of the turbine are listed in Table 2.
The natural frequencies of the tower change when the turbine is placed on a floating substructure, and preliminary calculations
showed that the natural frequency of the first bending mode for spar 2 was very close to the blade passing frequency at rated rotor
speed. The wall thickness in all sections of the tower was therefore increased by 50% for spar 2, which moved the natural frequency
above the 3P range.
The mooring system, described in Table 3, consists of three catenary lines spread symmetrically about the vertical axis. As in the
OC3 project, a simplification is made in that the delta lines are removed, and lines with constant properties are instead used all the
way up to the fairleads. A rotational spring is added to the model to ensure that the yaw stiffness from the mooring system is included.
The fairleads are placed at a depth equal to the total CoG of the system including the wind turbine, in order to limit the coupling
between surge and pitch motion. This results in a slightly stiffer mooring system for spar 2.
The fully coupled nonlinear aero-hydro-servo-elastic analyses are carried out in the time domain using the simulation workbench
SIMA developed by SINTEF Ocean, which couples two computer codes: Riflex, a finite element solver developed for flexible beam
elements; and SIMO, which calculates large volume hydrodynamic loads [17,18]. The spar buoys are modelled as six-DOF rigid
bodies with first order wave forces found numerically from potential flow theory using WAMIT [19], combined with viscous forces
from the drag term in Morison's equation, which are integated up to the instantaneous wave elevation. Second order potential flow
forces are not considered in the model, as studies have shown that these loads have limited effect on the response for spar-type FWTs
[20,21].
Bar elements with only axial stiffness are used to model the mooring lines, together with hydrodynamic loads from Morison's
equation. The model thus includes the nonlinear restoring forces from the mooring system, as well as the dynamic behaviour of the
mooring lines. The tower and blades are modelled using flexible beam elements.
The aerodynamic loads are calculated using blade element/momentum (BEM) theory, including Glauert correction, Prandtl hub
and tip loss factors, dynamic stall, dynamic wake, skewed inflow and tower shadow. The code has previously been verified for FWTs
[22,23]. As wind drag forces on the tower may become important at high wind speeds, this effect is included for the extreme response
conditions using a drag force formulation with a drag coefficient of 0.7.
An external control system written in Java is used to modify the generator torque and blade pitch. In order to avoid pitch motion
Table 3
Mooring system properties.
Parameter Spar 1 Spar 2
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
instability above rated wind speed due to negative feedback [24,25], the proportional and integral gains for the blade pitch are
modified from the original controller. With the resulting gains, the controller has a natural frequency of 0.13 rad/s and a damping
ratio of 0.7.
In the frequency domain model, a linear representation of the system is created as shown in Fig. 2. For simplicity, only the
response in the xz-plane is considered in the present study. The platform is considered rigid, while the tower is modelled as a slender
flexible beam. The rotor and nacelle assembly (RNA) is replaced by a point mass and inertia at the top of the tower, with resultant
wind loads (thrust force, rotor pitching moment, and aerodynamic damping) acting on the tower top. The mooring system is re-
presented by a linear spring at the position of the fairleads. As the structure is statically determinate, internal loads in the hull can be
found by considering equilibrium between external, inertial, damping, restoring and internal reaction forces. In the present work, the
tower base bending moment is the only internal load which is considered (see also Appendix A), but the method can easily be applied
to other components.
Generalized displacements are used in combination with the principle of virtual work to establish the dynamic equations of
motion [26], which are solved in the FD. In this procedure, the total response is described by a weighted combination of an arbitrary
number of shape functions, ψ. The accuracy of the formulation thus depends on how well the actual displacement field is captured by
the shape functions, which typically are chosen to represent the most important eigenmodes of the system.
Three generalized DOFs are included in the model, namely surge, pitch, and the first bending mode, as illustrated in Fig. 3.
Orthogonal eigenmodes can be found numerically from the solution of the eigenvalue problem for a finite element (FE) model of the
linearized system. However, as the shapes of the rigid body modes are known a priori, they are instead chosen as analytical functions
for simplicity. In addition, to have global motions that are consistent with the TD model, the pitch mode is chosen as a rotation about
the SWL. Using non-orthogonal modes introduces off-diagonal terms in the system matrices, however, these coupling effects are
readily taken into account in the formulation.
The actual shape of the bending mode, on the other hand, is not known, and is therefore found from the eigenvalue solution. A
third-order spline function is then fitted to the nodal displacements to have continuous expressions for the first and second order
derivatives. The FE model used to solve the eigenvalue problem is based on the linearized system in Fig. 2 and does not include the
blades, however, added mass on the hull is included using strip theory.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
In order to predict the dynamic behaviour of the system, the inertia, damping, restoring and excitation forces on the structure
need to be determined. Due to the simple and relatively slender geometry of the hull, the hydrodynamic loads on the floater can be
simplified without significant loss of accuracy. The transverse added mass per unit length is approximated using the analytical
expression for a 2D circular cylinder with diameter D in infinite fluid [27]:
The linear wave excitation forces are taken from MacCamy-Fuchs theory [28]. The force per unit length, dFW , for a regular wave
with unit amplitude is given by
4 g cosh k (z + h)
dFW (z , ) = G ei ( t )
k cosh kh (2)
Fig. 4. Aerodynamic damping ratio for the first fore-aft mode of the land-based DTU 10 MW turbine as a function of wind speed.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
where
1 J 1 (ka)
G= , tan = .
(J1 (ka))2 + (Y1 (ka))2 Y 1 (ka) (3)
For the aerodynamic loads, the BEM equations usually have to be solved in the TD [29], and there is no straight-forward way to
establish linearized frequency-dependent wind forces. In the present study, the turbulent thrust force (FT ) and rotor moment (MT ) for
different wind speeds are extracted from TD simulations on a fixed rotor in SIMA. The resulting time series are transformed to the FD
using fast Fourier transform (FFT), and normalized with the incoming wind spectrum to create force transfer functions. In addition,
the so-called aerodynamic damping, which arises due to the change in thrust force as a result of the nacelle's velocity, needs to be
explicitly calculated, which can be done in various ways [29,30]. Here, the linearized damping is estimated following the procedure
described by Bachynski [9]. Time domain simulations with a range of constant wind speeds are run on a fixed rotor, and the damping
values are found as the change in thrust force for small variations in wind speed while the blade pitch and rotor speed are kept fixed:
dFT
Baero = .
dU (4)
The resulting aerodynamic damping values are shown as a function of mean wind speed in Fig. 4, presented as critical damping
ratio for the first fore-aft tower mode of the original land-based DTU 10 MW turbine. The linearized aerodynamic damping applies to
all generalized modes with nonzero deflection at the tower top, and the coefficient Baero is varied with the mean wind speed. Damping
caused by pure rotation of the rotor is not considered. It is assumed that any changes in aerodynamic forces that arise due to the
motions of the FWT are captured by the aerodynamic damping term, which means that variations in control system outputs caused by
platform motions are not considered in the model.
Rayleigh damping is used to model structural damping in the tower, with a stiffness-proportional coefficient β = 0.007 for both
designs. For a given β, the damping ratio is proportional to the natural frequency, and using a constant value will thus reward stiffer
designs, such as spar 2, without any physical reason. In an actual design process, the coefficient should therefore be continuously
updated to keep the damping ratio constant. However, as the scope of the current work is not to perform design calculations, but
rather to compare different modelling techniques, a constant coefficient is used for simplicity. The applied coefficient corresponds to
a damping ratio of 0.9% and 1.2% at the first bending natural frequency for spar 1 and spar 2, respectively.
Viscous damping, which is important for the low-frequency surge response, is added based on stochastic linearization of the
quadratic drag term in Morison's equation [31]:
1 1 8
Cd D (z ) x |x | Cd D (z ) (x ) ( z ) x = b visc (z ) x
2 2 (5)
with a drag coefficient Cd = 0.7. The standard deviation of the velocity, (x ) , is found using an iteration scheme. Viscous wave
excitation was found to be small compared to the wind excitation forces, and is therefore not included in the model. Wave radiation
damping is also neglected. The aerodynamic drag forces on the tower, which only are included in the extreme response conditions,
are taken into account in the linear model as a combination of a mean force, which is added directly, and a frequency-dependent
force, which is found using stochastic linearization. Here, only the excitation forces are considered, meaning that the tower drag
forces arising from the movement of the turbine are neglected.
The analytical solution for quasi-static horizontal tension in elastic catenary lines can be found in e.g. Faltinsen [32]. In order to
find the linear mooring stiffness used in the model, Kt , this equation is differentiated numerically. The stiffness is calculated for the
zero offset position of the hull in all environmental conditions, as the offset was found to have little effect on the results. However,
updated stiffness values based on the actual mean position of the platform can easily be included. Mooring line excitation, damping
and inertia are neglected in the model.
The hydrostatic restoring in pitch, which is applied as a rotational spring at the SWL, is found from
Kr = gVzB Mgz G + gIwp. (6)
Using the vector containing the shape functions and the terms derived in the previous section, the generalized mass (including
added mass), damping and stiffness matrices can be established [26]:
ztop 0
M= m (z ) (z ) (z )dz + [m (z ) + a11 (z )] (z ) (z )dz
0 L
+ MRNA (ztop) (ztop) + IRNA , z (ztop) ,z (ztop) (7)
ztop 0
B= EI (z ) , zz (z ) , zz (z )dz + b visc (z ) (z ) (z )dz + Baero (ztop) (ztop)
z TB L (8)
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
ztop ztop
K= EI (z ) , zz (z ) , zz (z )dz + N (z ) , z (z ) ,z (z )dz
z TB z TB
The transfer functions relating modal response to wave and wind input, H X ( ) and HUX ( ) respectively, are defined as
H X ( ) = HFX ( ) F ( ), HUX ( ) = HFX ( ) FU ( ) (12)
where
HFX ( ) = [ 2M + i B + K] 1
(13)
is the frequency response function matrix. Transfer functions between tower base bending moment and wind/wave input (H M ( )
and HUM ( ) ) are found by considering equilibrium between external, inertial, damping and internal forces. The complete equations
may be found in Appendix A. As both the wind speed and wave elevation are considered to be stationary Gaussian processes within
the short-term duration, the linear response will also be a stationary Gaussian process. The responses to wind and wave input are
assumed to be independent; that is, there is no interaction between the responses at different frequencies. This is an assumption
inherent in the linearization, and may not be equally applicable for all FWTs.
The response spectrum for an arbitrary response parameter ξ is then found by superimposing the wind and wave responses:
S ( ) = |H ( )|2 S ( ) + |HU ( )|2 SU ( ), (14)
while the variance is given as
2
( ) = S ( )d .
0 (15)
5. Dynamic simulations
The long-term fatigue assessment uses the joint probability distribution of metocean parameters given in Johannessen et al. [33].
Probability density functions (PDFs) for the parameters are shown in Fig. 5. IEC 61400-3 [34] prescribes a minimum bin size of 2 m/s
for the wind speed, 0.5 m for wave heights, and 0.5 s for wave periods. Kvittem and Moan [35] performed a long-term fatigue analysis
of a semi-submersible FWT, where damage sensitivities with regards to simulation length, bin size and number of samples were
addressed. Based on their findings, the bin sizes for wave heights and wave periods in the present work are increased to 1 m and 1 s,
respectively, and each condition is simulated using six 1-h realizations after removal of transients.
Long-term fatigue analyses used in design of FWTs should cover all relevant conditions over the lifetime of the structure, including
operational, fault, idling, and survival conditions. The fatigue analysis presented in this study is limited to operational cases, which
means that the mean wind speed is varied between 4 and 25 m/s, and the turbine is assumed to be operating normally. A total of 546
environmental conditions (ECs) are thus considered in the analysis, which corresponds to all ECs within 4–25 m/s with probability of
occurrence greater than 10-4. It is worth noting that, although not employed in the current fatigue comparison, the FD model also is
applicable for wind speeds below cut-in or above cut-out, as long as only steady-state conditions are considered.
Three-dimensional turbulent wind time series are generated using TurbSim [36], assuming a Kaimal spectrum for IEC61400-1
class B turbines and a normal turbulence model. A power law profile with exponent 0.14 is used to model the vertical wind shear
[34]. The same wind files are used both directly in the TD simulations, and to establish the wind thrust and rotor moment transfer
functions used in the FD calculations, which are derived by taking the average values over the six seeds. The calculated transfer
functions for each individual wind time series, as well as the average curves, are shown in Fig. 6 for 11 m/s mean wind speed, to
illustrate the degree of variability between the realizations. Wave time series are generated using JONSWAP spectra with a peak-
edness parameter of 3.3. Co-directional waves and wind travelling in the positive x-direction are applied in all simulations, and no
current is considered in the analysis.
Four ECs along the 50-year contour surface are selected to compare the short-term extreme response, as described in Table 4. In
EC 1, the mean wind speed is close to the rated speed of the turbine, which gives the maximum thrust force on the rotor. EC 2
considers an operational condition near cut-out, while ECs 3 and 4 represent wind speeds above cut-out, where the turbine is parked
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Fig. 5. Marginal distribution of 1-h wind speed, and conditional distributions of wave height and peak period. Based on Johannessen et al. [33].
Fig. 6. Calculated transfer functions for aerodynamic excitation loads at 11 m/s mean wind speed. Six different wind realizations (grey) and average
values (black).
and the blades are feathered. As for the fatigue analysis, 1-h simulations (excluding transients) are used in the extreme response
comparison, but with 20 different random seeds for each condition. No fault conditions are considered, however, previous studies on
blade pitch faults, grid loss and shutdown have shown that such events primarily affect the loads on the blades and shaft, and that the
tower base fore-aft bending moment and global motions tend to be smaller than during fault-free conditions for spar-type FWTs
[37,38].
The fore-aft axial stress at the outer radius of the tower is calculated from
N My
x = + r.
A I (16)
Stress variations due to fluctuations in the axial force are found to be negligible compared to the moment-induced stress
Table 4
Environmental conditions for extreme response calculation.
Condition 1 2 3 4
Mean wind speed at hub height, U (m/s) 11.0 25.0 36.0 50.0
Significant wave height, Hs (m) 6.9 10.3 12.9 15.1
Spectral peak period, Tp (s) 10.0 12.5 14.0 16.0
Turbulence intensity at hub height, I (−) 0.18 0.14 0.13 0.12
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
variations, and are thus not included in the fatigue calculations. In the TD, rainflow counting is used together with SN curves and the
Miner-Palmgren rule to estimate the fatigue damage, which is taken as the average value over the six realizations for each EC. In the
FD, two different methods are considered:
In the first method, 1-h stress time series are generated by performing IFT on the response spectrum with random phase angles.
The fatigue damage is then calculated using rainflow counting and the Miner-Palmgren rule. The phase angles introduce randomness
to the damage estimates, and the method is therefore performed with six different random seeds, as for the TD model. The average
value is then used in the fatigue comparison.
The Dirlik method [39] uses the stress response spectrum and empirical factors to fit the PDF of the stress cycles to a combination
of an exponential and two Rayleigh distributions, and the accuracy of the method is therefore dependent on how well the rainflow-
count can be represented by these distributions. The expression for the fatigue damage during T seconds becomes [26]:
mk
pT m t
DDK = (2 )m G1 Qm (1 + m) + ( 2 )m (G2 Rm + G3) 1+
K 2 tref (17)
where
2(xm 2
2) 1 2 G1 + G12 1.25( 2 G3 G2 R) 2 xm G12
G1 = , G2 = , G3 = 1 G1 G2, Q = ,R=
1+ 2
2 1 R G1 1 2 G1 + G12 (18)
and
m2 m m2 1 m4
2 = , xm = 1 , p = .
m 0 m4 m0 m4 2 m2 (19)
mn = nS ( )d . (20)
0
The Dirlik method has previously been shown to give accurate results over a wide range of bandwidths for a stationary Gaussian
process [40], and the closed-form expression makes it well-suited for efficient design optimization. However, the accuracy of the
method deteriorates if the fatigue loads are dominated by a few frequency components [41].
The SN curve used for the tower base is the D curve in air taken from DNV-RP-C203 [42]. The curve is bilinear, i.e. the material
parameters, and thus also the slope of the curve, change at a certain stress range. For the rainflow procedure, this is easily taken into
account, however, the Dirlik method only allows for single-slope curves. Based on initial calculations, the Dirlik damage estimation is
performed with material parameters valid for stress ranges above 52.63 MPa (fewer than 107 cycles).
Several methods have been used to estimate the short-term extreme response of FWTs, such as the Gumbel method and Weibull
tail method for global motions [15,43], and the Winterstein method for the tower base bending moment [14]. The present study uses
the average upcrossing rate (AUR) method, see e.g. Naess and Gaidai [44], which previously has been shown to perform well for a
bottom-fixed offshore wind turbine under combined wind and wave loading [45]. The AUR method is based on the assumption that
upcrossings of high levels are statistically independent, which means that the number of upcrossings during time T will be Poisson
distributed. The cumulative distribution function (CDF) of the extreme value Y may thus be written as
FY (y ) = exp{ + (y ) T }, (21)
where Y can represent an arbitrary response parameter, such as tower base bending moment, global motions of the platform, or
mooring line tension. + (y ) is the mean upcrossing rate of level y, which can be estimated empirically from j simulated time series of
length T0 :
j
1
ˆ+(y ) = ni+ (y; T0 )
jT0 i=1 (22)
where ni+ (y ;
T0) is the number of upcrossings in time series i. Assuming that the appropriate asymptotic extreme value distribution for
the considered response is the Gumbel distribution, the tail of the mean upcrossing rate may be written as
+ (y ) = q (y ) exp{ a (y b)c }, y y0 (23)
where q (y ) is slowly varying and can be approximated as a constant. The parameters a, b, c and q are then determined by fitting Eq.
(23) to the empirical data, as described in detail by Saha and Naess [46]. This is done by minimizing the mean square error function:
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
N
F (a , b , c , q ) = wi |ln ˆ+(y ) ln q + a (yi b) c |2
i=1 (24)
using the Levenberg-Marquardt least-squares optimization algorithm. wi is a weight factor which puts more emphasis on the more
reliable data points, and is found from the estimated 95% confidence interval.
The AUR method is also applicable for the linear model, where an analytical expression is available. For a stationary Gaussian
process with zero mean, the upcrossing rate can be written as [26]:
+ (y ) + (0) y2
= exp ,
2 2 (25)
where + (0) is the mean zero-upcrossing rate, which can be found from the zeroth and second order spectral moments:
+ (0) 1 m2
= .
2 m0 (26)
The characteristic largest extreme response y1h , i.e. the most probable maximum value in one hour, is used to compare the TD and
FD models. The value is found approximately from the following relation:
FY (y1h ) = e 1. (27)
6. Results
The natural periods of the linear system are estimated and compared to results from decay tests performed with the TD model in
still water, as shown in Table 5. The decay tests are performed by releasing the platform from an offset position in the considered
DOF, and the natural period is then found from the decaying motion. For the first bending mode, the hull is released together with the
tower top, such that the coupling between platform pitch and tower bending is included.
The FD model slightly overestimates the natural periods, which for the rigid modes is found to mainly be a consequence of the
simplified added mass formulation. The reason why the first bending natural frequency is better approximated for spar 1 than spar 2
can be understood by examining the eigenmodes. In Fig. 7, the bending mode from the linearized eigenvalue solution is compared to
the mode shape found from a decay test in TD, where a bandpass filter around the natural frequency is used to extract the shape
including one of the blades. Due to the increased stiffness of the tower in spar 2, the coupled blade-tower mode changes, and the
blades undergo a larger amount of bending. The simplification of replacing the RNA by a point mass and inertia is thus less accurate
than for spar 1, where the blades behave more like a rigid body in the first tower bending mode. The shape of the tower is, however,
seemingly not affected by the blade behaviour, and both designs show good agreement between the decay tests and the eigenvalue
solutions.
Table 5
Natural periods.
Mode Spar 1 Spar 2
TD FD TD FD
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Scatter plots showing the calculated fatigue damage for each environmental condition in TD and FD are presented in Fig. 9, where
the Dirlik results are used for the linear model. The solid black line indicates perfect agreement, while conservative FD predictions are
found below the line.
The accuracy of the linear model is to a large extent insensitive to the amount of fatigue damage in the individual conditions, and
the data follow a fairly straight line, with most points located within ± 10% of the line of best fit for both designs. For spar 1, this line,
when neglecting the 9 m/s conditions, is very close to the diagonal. The 9 m/s results, although being highly overestimated by the
linear model, are also seen to be reasonably consistent, as the ratio between the linear and nonlinear estimates is more or less
constant with increasing damage.
For spar 2, all points are located below the 1:1 line, and as for spar 1, the conditions with the largest discrepancies between linear
and nonlinear predictions are in general found at lower wind speeds (5–9 m/s). These wind speeds are associated with low sea states,
and little excitation of the bending mode. The tower base bending moment is thus dominated by the low-frequency wind response,
which is less accurately predicted by the linear model than the wave- and bending-frequency response. The low-frequency response in
the linear model is also somewhat sensitive to how the thrust force spectrum is estimated, as reported by Kvittem and Moan [10]. It
should be noted that these conditions are associated with relatively small bending moments, and their contribution to the total
fatigue damage over the lifetime of the structure is thus not significant.
The limited scatter in the data for both designs suggest that it is possible to simulate a small number of environmental conditions
in the TD to verify and, if necessary, calibrate a linear model which can be used to perform the full long-term analysis. The ver-
ification cases should cover the entire wind speed range, as resonance effects may have a large impact on the accuracy of the model.
In addition, a 1 m/s interval for the wind speed bins may be needed to correctly capture the response around possible resonance wind
speeds [35].
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
natural frequency for spar 2 is increased. The effect is most prominent for the rotor moment, which is about 25% lower at the TD
natural frequency than at the natural frequency found from the linearized eigenvalue solution.
To quantify the importance of this error on the fatigue damage results, the natural period of the bending mode is shifted from
1.85 s to 1.79 s in the linear model by adding an artificial stiffness term, and the fatigue analysis is rerun. The results, presented in
Fig. 12, show a large improvement at most wind speeds, and the total fatigue damage error is reduced by approximately 50%. This
highlights the importance of having accurate estimates for the natural frequency and shape of the first bending mode, which is
dependent on how well the RNA can be approximated as a rigid body. In an optimization process, the accuracy of the model
simplification may vary as changes are made to the design, however, this can be taken into account by a straightforward calibration
of the natural frequency in the linear model.
The eigenmode in tower bending depends on wind speed, due to increased blade stiffness caused by centrifugal forces and
pitching of the blades, as well as changes in the mean platform pitch angle. This is illustrated in Fig. 13, where mode shapes from
decay tests are shown for calm air and 21 m/s uniform wind. As a consequence, the natural frequency of the bending mode is slightly
increased by 2–3% at higher wind speeds. In addition, the mode shape derivative at the tower top is reduced by up to 20%.
From Fig. 11, an increase in the natural frequency of the bending mode is seen to result in higher wind loads for spar 1. However,
this effect is counteracted by the decrease in the tower top derivative, which reduces the generalized wind load caused by the rotor
moment (see Eq. (10)). For spar 2, on the other hand, both effects reduce the wind excitation loads for the bending mode, and cause
the tuned linear model to somewhat overestimate the fatigue damage at higher wind speeds, as observed in Fig. 12.
Regarding the bending-frequency response, it should also be mentioned that, even though no validation has been performed, the
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Fig. 9. Time and frequency domain fatigue damage for individual ECs.
Fig. 10. Tower base stress response spectrum for U = 21 m/s, Hs = 7.5 m, Tp = 14.5 s.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Fig. 11. Wind excitation transfer functions around first bending natural frequency, U = 21 m/s.
Fig. 12. 1-h weighted fatigue damage arranged by wind speed, spar 2 with tuned natural frequency for first bending mode.
present study assumes that the behaviour of the physical system is described correctly by the TD model. Full-scale measurements
from the Hywind Demo FWT [48] have shown that the tower base bending moment response spectra contain energy around the first
tower bending mode and that the numerical models thus capture real effects, however, the actual level of the response is sensitive to
the modelling of the system, and may therefore not be correctly represented in the simulations.
The short-term extreme response of the system in the four 50-year environmental conditions is calculated for the tower base
bending moment, as well as for the surge and pitch motions of the platform. The results are presented in the following subsections.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Fig. 13. Normalized bending mode shapes in calm air and for U = 21 m/s. Natural periods found from decay tests are 2.41 s (calm air) and 2.35 s
(21 m/s) for spar 1, and 1.79 s (calm air) and 1.75 s (21 m/s) for spar 2.
Table 6
Tower base bending moment statistics, 50-year environmental conditions.
Model Statistical quantitya Spar 1 Spar 2
Condition Condition
1 2 3 4 1 2 3 4
Time domain Mean (MNm) 247.2 123.0 50.3 105.2 305.7 152.9 63.6 132.8
Std.dev. (MNm) 89.1 101.0 106.5 118.8 103.8 106.2 115.5 129.7
Skewness −0.20 −0.03 −0.01 0.01 −0.30 −0.04 −0.01 0.01
Kurtosis 3.03 3.00 3.00 3.01 3.13 3.00 3.01 3.01
Frequency domain Mean (MNm) 260.8 127.8 52.3 101.3 329.4 161.4 66.5 128.9
Std.dev. (MNm) 89.4 100.8 106.3 115.1 104.8 110.6 118.4 128.5
a
Time domain results are averaged over 20 seeds.
4.1%. The mean moments are somewhat less accurate, with the largest errors of 5.5% (spar 1) and 7.8% (spar 2) occurring in EC 1,
near rated wind speed, where also the largest mean values are found.
Empirical upcrossing rates are plotted together with optimized curve fits and Gaussian upcrossing rates based on the FD results in
Figs. 14 and 15. The plots show the multiplication factor κ, which is related to the extreme response, y, through the relation
y=μ+ , (28)
where μ and σ are the mean and standard deviation of the process. The characteristic largest 1-h extreme for a given environmental
condition, which is used in the comparison, corresponds to an upcrossing rate of 2.78 × 10 4 .
The Gaussian upcrossing rates are seen to be conservative for EC 1, due to the slightly nonlinear behaviour which becomes more
prominent as the bending moments increase. In the parked conditions, on the other hand, the response is linear also at lower
upcrossing rates, and the empirical data are well-described by the linear model.
For spar 1, the FD solution is seen to follow the optimized curve fit closely in EC 2, however, larger discrepancies are observed for spar 2.
The fact that the shape of the fitted curve is captured quite well, suggests that the error is related to the overprediction of the bending
response discussed in the previous section. As the bending mode is located at a relatively high frequency, the response level has significant
influence on the second order spectral moment. Consequently, the mean zero-upcrossing rate (see Eq. (26)) is overestimated by 29%
compared to the TD model, which causes a right shift in the linear upcrossing rate curve as observed in Fig. 15b.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Fig. 14. Tower base bending moment upcrossing rates for spar 1.
The calculated 1-h extreme bending moments are listed in Table 7. The linear estimates give reasonably good approximations of
the extreme response, particularly for EC 2–4, as indicated by the upcrossing results.
Due to the high mean value caused by the thrust force, the largest extreme response is found in EC 1. However, this may not be the
case in general, and previous studies comparing short-term extremes in selected operational and survival conditions have shown that
the largest tower base bending moments for a spar FWT can occur in storm conditions with extreme wind speeds [49,50]. The critical
load cases may vary with floater design and turbine size, as it alters the ratio between wind and wave loads, and for long-term
extremes, the probability of the conditions must be taken into account. The most critical operational condition is also not necessarily
associated with maximum thrust force on the turbine. In a recent study by Sultania and Manuel [51], the largest 50-year tower base
bending moment for a 5 MW spar FWT, considering only operational cases, was found for a mean wind speed of 21.7 m/s. This
suggests that a wider range of environmental conditions should be considered in an actual design process.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Fig. 15. Tower base bending moment upcrossing rates for spar 2.
Table 7
Characteristic largest extreme tower base bending moment y1h , 50-year environmental conditions.
Condition Condition
1 2 3 4 1 2 3 4
Time domain (MNm) 534.2 480.6 418.7 509.8 624.1 529.6 462.6 577.4
Frequency domain (MNm) 575.7 491.7 415.4 491.3 699.3 562.0 469.9 563.3
Relative error (%) 7.8 2.3 −0.8 −3.6 12.0 6.1 1.6 −2.4
least accurate linear predictions in this condition. Large discrepancies are also seen in the most extreme wave condition (EC 4),
mostly due to underestimation of the standard deviations in FD. For the surge motions, the extremes found from FD generally agree
well with the TD simulations. Most notable are the results in EC 1, especially for spar 1, where the relative error between FD and TD is
only −1%. However, this result is somewhat coincidental, as it is a consequence of the error in κ being almost exactly balanced out by
the errors in μ and σ. Nonetheless, the results suggest that the extreme surge motions may be predicted with roughly the same
accuracy as the tower base bending moments, despite larger discrepancies in the standard deviations, as the mean values account for
a larger share of the total extreme response.
The accuracy of the linear model depends on the importance of nonlinearities in the system, which will vary with location, floater
geometry, and different operational conditions. Therefore, the results presented here cannot be generalized to any FWT, and the FD
model should be continuously verified against state-of-the-art tools during a design process. The agreement with nonlinear analyses
also depends on how the system is represented in the TD model, and care should be taken to make sure that relevant nonlinearities
are included.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Table 8
Platform surge response statistics, 50-year environmental conditions.
Model Statistical quantitya Spar 1 Spar 2
Condition Condition
1 2 3 4 1 2 3 4
Time domain Mean (m) 27.3 13.8 8.0 15.1 23.7 12.0 6.6 13.7
Std.dev. (m) 6.9 2.4 2.9 4.6 6.2 2.2 2.7 3.9
Skewness −0.14 −0.02 0.05 −0.03 −0.13 −0.01 0.01 0.06
Kurtosis 2.49 2.90 3.03 3.06 2.57 2.95 3.02 2.97
Frequency domain Mean (m) 28.6 14.1 7.9 15.3 26.1 12.8 7.0 13.6
Std.dev. (m) 5.4 2.3 2.5 3.5 5.2 2.2 2.5 3.5
a
Time domain results are averaged over 20 seeds.
Table 9
Platform pitch response statistics, 50-year environmental conditions.
Model Statistical quantitya Spar 1 Spar 2
Condition Condition
1 2 3 4 1 2 3 4
Time domain Mean (deg) 6.7 3.4 1.7 3.4 9.2 4.6 2.2 4.4
Std.dev. (deg) 1.8 1.0 1.2 1.7 2.4 1.3 1.3 2.1
Skewness −0.62 0.03 0.04 0.10 −0.72 0.03 0.03 0.10
Kurtosis 3.37 2.90 3.02 2.99 3.54 2.92 2.99 3.02
Frequency domain Mean (deg) 7.2 3.5 1.7 3.2 10.0 4.9 2.2 4.3
Std.dev. (deg) 1.7 1.0 1.0 1.4 2.3 1.2 1.2 1.7
a
Time domain results are averaged over 20 seeds.
Table 10
Characteristic largest extreme surge response y1h , 50-year environmental conditions.
Condition Condition
1 2 3 4 1 2 3 4
Time domain (m) 44.2 21.4 17.8 29.4 38.6 19.0 15.7 26.4
Frequency domain (m) 43.8 21.4 16.1 26.6 40.7 20.0 15.1 24.9
Relative error (%) −1.0 0.0 −9.5 −9.4 5.4 5.4 −3.6 −5.7
Table 11
Characteristic largest extreme pitch response y1h , 50-year environmental conditions.
Condition Condition
1 2 3 4 1 2 3 4
Time domain (deg) 10.7 6.7 5.5 9.1 14.4 8.6 6.6 11.3
Frequency domain (deg) 12.3 6.6 5.0 7.7 16.7 8.7 6.0 9.8
Relative error (%) 14.8 −0.9 −9.4 −14.9 16.2 −4.9 −8.9 −13.2
7. Conclusions
A semi-analytical FD model has been developed, which allows for efficient preliminary design optimization of spar FWTs. The
model has been verified against a fully coupled nonlinear TD model for two different 10 MW spar designs, with regards to fatigue
damage in the tower base and short-term extreme response in both operational and survival conditions. The fatigue results show that
the FD model is accurate in most environmental conditions, especially for spar 1, where the long-term fatigue damage is
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
overestimated by 7%. For spar 2, larger discrepancies are observed, and the linear model overpredicts the damage by approximately
30%. This is found to primarily be caused by errors in the response around the first bending mode, which is highly sensitive to the
estimated natural frequency. The accuracy depends on how well the RNA can be approximated as a rigid body for the first bending
mode, however, tuning the frequency by adding an artificial stiffness term can be performed in the case where the approximation
leads to unacceptable errors.
For 9 m/s mean wind speed, the 3P frequency of the rotor coincides with the natural frequency of the first bending mode for spar
1, which leads to increased damage values and causes the linear model to significantly overpredict the response. The sensitivity with
regards to wind speed also suggests that a finer discretization of the wind bins should be applied in a design process, to accurately
capture the response around possible resonances. The accuracy of the linear model is found to be fairly independent of wave height
and wave period, also for the wind speed which induces 3P resonance with the first tower bending frequency. This indicates that a
small number of environmental conditions, which cover the full wind speed range, can be simulated with a TD model to verify and
possibly tune the linear model, which then can be used to carry out the full long-term analysis.
The extreme response, both with regards to tower base bending moment and rigid body motions in surge and pitch, is seen to be
nearly Gaussian near cut-out and in parked conditions, as it is dominated by wave forces. This results in good agreement for the linear
model in these conditions. The largest extremes values are found in EC 1, near rated wind speed, which indicates that a wider range of
environmental conditions may be important for the extreme response. The large aerodynamic forces in this condition also introduce
some non-Gaussian behaviour, which causes the linear model to overestimate the multiplication factors (κ). However, as the mean
values in this condition are relatively large, the inaccuracies in the total extreme response are reduced, and the total error in y1h is
within 16% for all considered responses in this condition. Regarding the standard deviations, good agreement between the linear and
nonlinear models is obtained for the tower base bending moments, while larger discrepancies are seen for surge and pitch. The reason
for this is that the global motions are more dependent on the low-frequency (wind-induced) response, which is less accurately
captured by the linear model than the wave- and bending-frequency response. The TD model in the present study only includes first
order wave forces, and could thus be improved by also considering higher order loads, which may be important for certain FWT
concepts. The model should also be validated through comparisons with physical measurements.
The developed FD model is shown to be well-suited for preliminary design. It gives reasonable agreement with a fully coupled
nonlinear TD model, both with regards to fatigue and extreme response, and the semi-analytical approach provides an efficient way
to explore the design space. In addition, the closed-form expressions offers the possibility of using analytic gradients in a computer-
aided optimization process. However, the efficiency of the model is partly due to the simple geometry of the platform, and the
formulations are thus not valid for more complex platform designs, where more comprehensive numerical calculations may be
necessary. Also, the approach is only applicable for steady-state conditions, and transient events, which also must be considered,
cannot be analysed. Nonetheless, it can serve as a useful complement to more sophisticated models in the design process of spar-type
FWTs, and significantly reduce the needed computational effort.
A Transfer functions
The transfer function from wave elevation to tower base bending moment is given in Eq. (A.1). The equation considers bending
moments due to inertial forces in the tower and RNA, aerodynamic damping forces, and gravitational forces.
ztop
H M( ) = m (z ) H X¨ ( ) (z )(z zTB)dz MRNA H X¨ ( ) (ztop)(ztop zTB)
z TB
IRNA H X¨ ( ) , z (ztop) Baero H X ( ) (ztop)(ztop zTB)
ztop
+ m (z ) g H X ( )[ (z ) (zTB)]dz + MRNA g H X ( )[ (ztop) (zTB)]
z TB (A.1)
The transfer function from wind speed to tower base bending moment is given in Eq. (A.2). The equation considers bending
moments due to inertial forces in the tower and RNA, aerodynamic damping forces, gravitational forces, and wind excitation forces.
ztop
HUM ( ) = m (z ) HUX¨ ( ) (z )(z zTB)dz MRNA HUX¨ ( ) (ztop)(ztop zTB)
z TB
IRNA HUX¨ ( ) , z (z top ) Baero H X ( ) (ztop)(ztop zTB)
ztop
+ m (z ) g HUX ( )[ (z ) (zTB)]dz + MRNA g HUX ( )[ (ztop) (zTB)]
z TB
+ FT ( )(ztop zTB) + MT ( ) (A.2)
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
B.1 Surge
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
B.2 Pitch
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
Funding
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
References
[1] Strach-Sonsalla M, Muskulus M. Dynamics and design of floating wind turbines. Proceedings of the twenty-sixth (2016) international ocean and polar en-
gineering conference (ISOPE2016), rhodes, Greece. 2016.
[2] DNV. Environmental Conditions and Environmental Loads. 2014. Tech. rep. DNV-RP-C205.
[3] DNV. Design of Floating Wind Turbine Structures. 2013. Tech. rep. OS-J103.
[4] Wayman EN, et al. Coupled dynamic modeling of floating wind turbine systems. Offshore technology conference. 2006.
[5] Sclavounos P, Tracy C, Lee S. Floating offshore wind turbines: responses in a seastate pareto optimal designs and economic assessment”. proceedings of the ASME
27th international conference on offshore mechanics and arctic engineering (OMAE2008), estoril, Portugal. 2008.
[6] Philippe M, Babarit A, Ferrant P. Comparison of time and frequency domain simulations of an offshore floating wind turbine. Proceedings of the ASME 2011 30th
international conference on ocean, offshore and arctic engineering (OMAE2011), rotterdam, The Netherlands. 2011.
[7] Wang K, et al. Frequency domain approach for the coupled analysis of floating wind turbine system. Ships Offshore Struct 2017;12(6):767–74.
[8] Bachynski EE, Moan T. Linear and nonlinear analysis of tension leg platform wind turbines. Proceedings of the twenty-second (2012) international offshore and
polar engineering conference (ISOPE2012), rhodes, Greece. 2012.
[9] Bachynski EE. Design and dynamic analysis of tension leg platform wind turbines PhD thesis Norwegian University of Science and Technology; 2014
[10] Kvittem MI, Moan T. Frequency versus time domain fatigue analysis of a semisubmersible wind turbine tower. J Offshore Mech Arctic Eng 2014;137:1.
[11] Kluger JM, Sapsis TP, Slocum AH. A reduced-order, statistical linearization approach for estimating nonlinear floating wind turbine response statistics”.
Proceedings of the twenty- sixth (2016) international ocean and polar engineering conference (ISOPE2016), rhodes, Greece. 2016.
[12] Jonkman J. Definition of the Floating System for Phase IV of OC3. National Renewable Energy Laboratory; 2010. Tech. rep. NREL/TP-500-47535.
[13] Matha D. Model development and loads analysis of an offshore wind turbine on a tension leg platform, with a comparison to other floating turbine concepts MA
thesis University of Colorado - Boulder; 2009
[14] Karimirad M, Moan T. Extreme dynamic structural response analysis of catenary moored spar wind turbine in harsh environmental conditions. J Offshore Mech
Arctic Eng 2011;133.
[15] Aggarwal N, Manikandan R, Saha N. Nonlinear short term extreme response of spar type floating offshore wind turbines”. Ocean Eng 2017;130:199–209.
[16] Bak C, et al. Description of the DTU 10 MW Reference Wind Turbine. 2013. Tech. rep. DTU Wind Energy Report-I-0092. DTU Wind Energy.
J.M. Hegseth, E.E. Bachynski 0DULQH6WUXFWXUHV ²
128 Appended Papers
Paper 2
129
130 Appended Papers
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
1. Introduction
Frequency-domain (FD) methods can be helpful in the design of floating wind turbines (FWTs),
providing responses for a large number of loading conditions with relatively low computational
time. Previous work has indicated that the response of FWTs in moderate environmental
conditions may be estimated reasonably well using FD models [1, 2, 3]. However, their reliability
depends on an accurate prediction of loads and interactions with linearized models.
An especially interesting interaction takes place between nacelle motions and the rotor thrust.
Fluctuations in the nacelle’s horizontal velocity provoke changes in the flow through the rotor,
leading to oscillations in the rotor speed and thrust. An aerodynamic damping effect results from
this interaction, and is normally included in FD models by means of a constant aerodynamic
damping coefficient. This coefficient may be obtained from the thrust derivative w.r.t. the
relative wind speed [1, 3, 4], or by means of decay simulations of the FWT under different
incident wind speeds [2].
This method is convenient to implement due to its relative simplicity, and is normally adopted
in combination with the frequency-dependent thrust obtained for a fixed wind turbine, installed
on the top of the tower. A disadvantage of this approach is that the interactions of the nacelle
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
motions with the blade pitch control system are not considered, an effect that is known to reduce
the damping or even render it negative [5, 6, 7], besides introducing an inertia effect [8]. As a
consequence, the damping effect can be significantly overpredicted at frequencies close to the
controller bandwidth, if the constant coefficient is used.
Control system effects can be included in an FD model by means of linearization of the
thrust and torque. The rotor speed is then included as an additional degree of freedom, and the
blade-pitch controller can be written in terms of the system states [9].
The objective of the present work is to relate the linearized thrust equations to the above-
mentioned inertia and damping effects, providing a visualization of how those effects vary with
wind speed and frequency for a given control system.
The aerodynamic inertia and damping coefficients are calculated both from the linearized
thrust equations and from simulations where the nacelle is forced to oscillate with different
frequencies, under uniform wind. The obtained coefficients are then included in an FD model
of a semi-submersible 10 MW FWT. The responses predicted with the different methods for
obtaining the coefficients are compared to coupled time-domain simulations, where the thrust
is calculated with the blade element momentum (BEM) theory.
Section 2 presents the frequency domain model for the FWT; Section 3 introduces the thrust
linearization procedure and provides formulations for the frequency-dependent aerodynamic
inertia and damping coefficients; the time-domain simulations are explained in Section 4, and
the results are presented in Section 5; final discussions and conclusions are provided in Section 6.
where zhub is the hub height above still water level (SWL). Transfer functions from wave and
wind input to platform response, HζX (ω) and HU X (ω) respectively, are expressed as
where −1
HF X (ω) = −ω 2 (M + A(ω)) + iωB(ω) + K (3)
is the frequency response function matrix. After solving for the platform motions, the tower
base bending moment response is found by considering dynamic equilibrium.
Using the wind and wave spectra, the response spectrum for an arbitrary response parameter
ξ is then expressed as
2
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
∂T
T = T0 + Δv . (5)
∂v
Let the nacelle dynamics be represented by an 1-DOF, 1st -order system:
Uniform, non-turbulent wind is assumed for now. Recalling that v = Uw − ẋ, it is noted that
Δv = v − Uw = −ẋ. Replacing (5) in (6), the damping effect becomes clear:
∂T
mẍ + c + ẋ + kx = T0 , (7)
∂v
∂T
baer = . (8)
∂v
The damping effect estimated with this method considers the change in the steady-state
thrust for a small perturbation in the uniform wind speed seen by the rotor, while neglecting
the effect of the control system and rotor dynamics. Different coefficients are obtained for
each incident wind velocity, but they are constant with respect to the nacelle frequency of
oscillation. This method is relatively simple to use and normally provides satisfactory results
when the turbine oscillates in a frequency range distant from the controller bandwidth. For
lower frequencies of oscillation, interaction with the control system takes place and the damping
coefficient dependence on frequency becomes more important.
T = T0 + Tv Δv + TΩ ΔΩ + Tβ Δβ , (9)
∂A
where the index indicates partial derivative w.r.t. the indicated variable, i.e., AB = ∂B . In order
to find the damping coefficient, Equation (9) must be written in terms of ẋ only. First, the
aerodynamic torque, Q, is also given in its linearized version:
Q = Q0 + Qv Δv + QΩ ΔΩ + Qβ Δβ , (10)
3
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
Id Ω̇ = Q − Ng E , (11)
where Id is the drivetrain inertia, Ng is the gear ratio and E is the electrical (generator) torque,
here assumed as constant above rated wind speed. In equilibrium, Q0 = Ng E, so replacing (10)
in (11):
Id Ω̇ = Qv Δv + QΩ ΔΩ + Qβ Δβ . (12)
It is assumed that a PI controller commands the blade pitch angle, i.e.,
Δβ = Kp ΔΩ + Ki ΔΩdt , (13)
where Kp and Ki are the proportional and integral controller gains. Equations (9) and (12) are
now rewritten replacing Δβ as given in (13), and noting that ΔΩ = φ̇, where φ is the rotor
azimuth angle:
T = T0 − Tv ẋ + (TΩ + Kp Tβ ) φ̇ + Ki Tβ φ , (14)
Id φ̈ − (QΩ + Kp Qβ ) φ̇ − Ki Qβ φ = −Qv ẋ . (15)
Assuming harmonic oscillation, Equation (15) can be written in the frequency domain, and the
following transfer function between x(ω) and φ(ω) is obtained:
iωQv
φ(ω) = x(ω) = C(ω)x(ω) . (16)
Id ω 2 + (QΩ + Kp Qβ ) iω + Ki Qβ
The nacelle dynamics (Equation 6) are now written in terms of T as given in (14), and also
expressed in the frequency domain, with φ(ω) as given in (16):
⎧ ⎫
⎪
⎨ ⎪
⎬
−mω 2 + iωc + k + iωTv − [(TΩ + Kp Tβ ) iω + Ki Tβ ] C(ω) x(ω) = T0 . (17)
⎪
⎩ ⎪
⎭
−ω 2 aaer (ω)+iωbaer (ω)
Based on the assumption that T responds harmonically to harmonic oscillations of the nacelle,
the thrust can be written as a combination of terms proportional to the nacelle acceleration and
velocity [8]. This assumption is further discussed in Section 3.3, but an immediate consequence
is the definition of frequency-dependent, aerodynamic inertia and damping coefficients:
1
aaer (ω) = − Re {iωTv − [(TΩ + Kp Tβ ) iω + Ki Tβ ] C(ω)} (18)
ω2
1
baer (ω) = Im {iωTv − [(TΩ + Kp Tβ ) iω + Ki Tβ ] C(ω)} (19)
ω
4
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
0.8 1080
0.6 1060
1040
0.4
Thrust (kN)
1000
0
980
-0.2
960
-0.4
940
-0.6 920
-0.8 900
100 150 200 250 300 350 400 450 500
Time (s)
Figure 1: Thrust response to oscillations of the DTU 10 MW wind turbine, under constant
incident wind speed of 13.0 m/s and period of oscillation of 100.0 s.
the control system. When the turbine oscillates harmonically, the thrust also oscillates nearly
harmonically, but with a phase relative to the nacelle velocity. This can be seen in Figure 1,
where the DTU 10 MW wind turbine is forced to oscillate with period 100.0 s and under constant
incident wind speed of 13.0 m/s. The thrust can then be assumed to be composed of a mean
plus an oscillating component:
ẋ = v0 cos(ωt) , (21)
the oscillating part of the thrust can then be assumed to be given by [8]
Tosc = f0 v0 cos(ωt + α)
= f0 [v0 cos(ωt) cos(α) − v0 sin(ωt) sin(α)]
(22)
f0 sin(α)
= f0 cos(α)ẋ + ẍ .
ω
The amplifying factor f0 and the phase α between nacelle velocity and thrust depend on
the dynamic effects mentioned above, and are not straightforward to determine analytically.
However, for a given turbine and control system they vary with the incident wind speed and
frequency of oscillation, while motion amplitude does not seem to have a significant influence.
They can therefore be obtained from forced oscillations covering the ranges of interest for both
parameters. Fast Fourier transforms (FFTs) of the nacelle velocity and thrust are calculated,
and f0 and α can be directly obtained from the ratio between amplitudes or phases of both FFTs
at the period of interest. Formulations for the aerodynamic inertia and damping coefficients are
then directly determined from equation (23):
f0 sin(α)
aaer (Uw , ω) = − , baer (Uw , ω) = −f0 cos(α) , (23)
ω
5
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
obtained from the linearized BEM equations. For the latter, forced oscillations of the DTU
10 MW wind turbine were performed with a simulator which obtains aerodynamic loads from
AeroDyn [10], updates the blade pitch angle based on the same controller strategy presented
in Section 4, and updates the rotor speed based on the rotor dynamics as in Equation (11).
Oscillation periods varied from 20.0 s to 130.0 s, and uniform wind with the three velocities in
Table 3 was considered. The BEM formulation with dynamic stall is adopted for the AeroDyn
calculations.
The coefficients obtained with both approaches are shown in Figures 2 and 3. In spite of the
oscillatory character of the curves obtained with the oscillation method, the agreement between
both methods looks satisfactory, especially for higher wind velocities. The inertia effect is of
the order of 1% of the mass of typical FWTs, resulting in negligible consequences for the surge
dynamics. The contribution to the moment of inertia in pitch is however considerable, due to
the nacelle height. Important changes in the pitch natural period can then result, as shown
in [8]. The constant damping coefficients as obtained in equation (5) are also plotted in Fig.
3, showing that the discrepancy w.r.t. the frequency-dependent coefficients is larger for lower
frequencies and higher wind velocities.
3.5. The aerodynamic inertia and damping effects on the FWT dynamics
It was already shown that the fluctuations in the aerodynamic thrust resulting from the nacelle
motions may be treated as frequency-dependent inertia and damping effects. The derivation
assumed non-turbulent wind, under the hypothesis that the fluctuations in the apparent wind
flow in the rotor due to turbulence can be decoupled from those caused by the nacelle motions.
Under the same assumption, the fluctuations in the thrust due to nacelle oscillations can now
be added to the FD model presented in Section 2. This is done by noting that, from Equations
(22) and (23), the oscillating component of the thrust can be written as
The following aerodynamic and damping matrices, Aaer and Baer , are then defined:
⎡ ⎤ ⎡ ⎤
aaer (ω) 0 aaer (ω)zhub baer (ω) 0 baer (ω)zhub
Aaer (ω) = ⎣ 0 0 0 ⎦ , Baer (ω) = ⎣ 0 0 0 ⎦.
aaer (ω)zhub 2
0 aaer (ω)zhub 2
baer (ω)zhub 0 baer (ω)zhub
(25)
Aaer and Baer are now summed to A and B in equation (3), which includes the effect of nacelle
motion in the FD analysis.
4. Simulations
The platform considered is the OO-Star 10 MW [11], a semi-submersible concept designed for
the LIFES50+ project. The main properties are summarized in Table 1. The potential-theory
hydrodynamic loads are generated with WADAM, and viscosity is added in form of the Morison
drag term. The platform is installed at a water depth of 130.0 m, and the mooring system
consists on a simplified 3-catenary line arrangement, which differs from the original of taut lines
with clump weights.
The DTU 10 MW turbine (Table 2) is installed at the top of the tower, at a height of 118.4 m.
The blade-pitch controller corresponds to the PI formulation of Equation (13), i.e., the power
error contribution of the original controller [12] is not included. The controller gains are tuned
such that its bandwidth is below the pitch natural period, and a gain scheduling strategy corrects
the gains according to the current blade pitch angle.
The time-domain simulations were performed with SIMA, a workbench which allows for
coupled analyses between floating bodies and slender elements [14, 15]. The platform, nacelle
6
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
106 105
3 8
Calculated Calculated
Forced Forced
2.5 7
2 Tn,5 6 Tn,5
aaer (kg)
aaer (kg)
1.5 5
1 4
0.5 3
0 2
0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05
Frequency (Hz) Frequency (Hz)
(a) Uw = 13.0 m/s. (b) Uw = 19.0 m/s.
105
4.5
Calculated
Forced
4
3.5 Tn,5
aaer (kg)
2.5
1.5
0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05
Frequency (Hz)
(c) Uw = 25.0 m/s.
Figure 2: Aerodynamic inertia coefficients, calculated with the linearized approach or obtained
from the forced oscillations method.
and hub are modeled as rigid bodies, while the blades, tower and mooring lines are modeled as
flexible structures, with finite elements. The aerodynamics are based on the BEM theory, with
dynamic stall and dynamic wake effects based on Øye’s models [16]. Hub and tip losses are
modeled with the Prandtl factor, and Glauerts correction is applied for high induction factors.
The environmental conditions are listed in Table 3. Turbulent wind was generated with TurbSim
[17], using the IEC Kaimal spectrum with turbulence characteristic B (moderate). A JONSWAP
7
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
104 104
16 16
14
14
12
12 Tn,5
10
10
baer (N.s/m)
baer (N.s/m)
8 Tn,5
6 8
4
6
2
4
0
Calculated Calculated
Forced 2 Forced
-2
Constant Constant
-4 0
0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05
Frequency (Hz) Frequency (Hz)
(a) Uw = 13.0 m/s. (b) Uw = 19.0 m/s.
104
16
14
12 Tn,5
10
baer (N.s/m)
4
Calculated
2 Forced
Constant
0
0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05
Frequency (Hz)
(c) Uw = 25.0 m/s.
Figure 3: Aerodynamic damping coefficients, calculated with the linearized approach (both
constant and frequency-dependent) or obtained from the forced oscillations method.
Table 2: Main properties of the DTU 10MW wind turbine [13]. The controller gains refer to
β = 0◦ .
wave spectrum with a γ factor of 3.3 was used to describe the sea state. The simulations lasted
for one hour, which showed to provide a sufficient number of LF oscillations for the analysis.
Only one realization of each condition was considered.
8
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
5. Results
Figure 4 shows the frequency-domain pitch response for the three conditions in Table 3, obtained
with time-domain simulations in SIMA; and with the FD model from Section 2, using the
different approaches for including the aerodynamic damping and inertia effects, as presented in
Section 3. When the constant damping coefficient from Equation (8) is adopted, the method
is referred to as FD constant; with the frequency-dependent inertia and damping coefficients
as given by Equations (18) and (19), the approach is named FD calculated ; and when the
coefficients are obtained from the forced harmonic oscillations, as in Equation (23), the FD
forced identifier is adopted.
The curves indicate that the adoption of a constant damping coefficient overestimates the
damping effect at lower frequencies, resulting in significantly underpredicted resonant response
for frequencies below 0.05 Hz. When frequency-dependent coefficients are adopted, however,
the situation is the opposite: the damping effect is underestimated, resulting in higher response
for lower frequencies. Using the coefficients based on the forced oscillations leads to improved
results in comparison with those obtained from the linearized thrust equations, in terms of
the pitch motion standard deviation (Table 4). All three methods perform equivalently in the
wave-frequency range (0.05-0.15 Hz), where the response is underpredicted due to the absence
of viscous excitation in the FD model.
Similar conclusions can be obtained regarding the tower-base bending moment (Figure 5).
The most remarkable difference in comparison with the plots for the pitch response is for Uw =
13.0 m/s (Figure 5a), which shows a persistent underprediction for all the three FD models
for frequencies above 0.04 Hz. The errors in predicted tower base bending moment, shown in
Table 4, are in general lower than for the pitch motion, since wave-frequency and 3p loads also
influence the tower-base bending moment, and are less affected by the damping dependence of
frequency. The accuracy is satisfactory for calculated and, especially, forced methods.
Figure 6 compares the low-frequency responses in pitch, for Uw = 19.0 m/s, when the inertia
effect is included or disregarded. The FD model with damping coefficient obtained from forced
oscillations is adopted, but only the blue curve includes the inertia coefficient – resulting in a
peak period about 2.5 s longer then when the inertia effect is not considered.
6. Conclusion
A method was developed to illustrate the importance of including the interaction between
the nacelle motions and thrust in the frequency-domain representation of a FWT. The thrust
fluctuations resulting from nacelle motions can be interpreted as frequency-dependent inertia and
damping effects, which challenges the traditional approach of adopting a constant aerodynamic
damping coefficient.
Expressions were derived for the respective coefficients based both on a linearized expression
for the thrust and on forced oscillations of a wind turbine in the time domain. The obtained
coefficients were included in a 3 DOF FD model of the OOstar semi-submersible FWT.
Comparisons were made between the models using a constant aerodynamic damping coefficient,
the frequency-dependent inertia and damping coefficients obtained in two different ways, and
9
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
10-1 10-1
SIMA SIMA
FD constant FD constant
-2 -2
10 FD calculated 10 FD calculated
FD forced FD forced
10-3 10-3
Sx5 (rad2.s)
Sx5 (rad2.s)
10-4 10-4
10-5 10-5
10-6 10-6
10-7 10-7
0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2
Frequency (Hz) Frequency (Hz)
(a) Uw = 13.0 m/s. (b) Uw = 19.0 m/s.
SIMA
10-2
FD constant
FD calculated
FD forced
10-4
Sx5 (rad2.s)
10-6
10-8
0 0.05 0.1 0.15 0.2
Frequency (Hz)
(c) Uw = 25.0 m/s.
time-domain simulations. Both platform pitch motion and tower base bending moment were
analyzed. It was observed that the constant damping coefficient underestimates the responses
in the frequency range of the FWT surge and pitch natural frequencies.
The frequency-dependent coefficients, on the other hand, overestimated the response in the
same frequency range, but are closer to the time-domain predictions. The aerodynamic inertia
was shown to increase the pitch natural period by about 2.5 s, for a mean incident wind speed
of 19.0 m/s. The differences between using coefficients based on the linearized thrust or on the
forced oscillations were not very significant, but a slightly better agreement of the LF responses
with the time-domain simulations was attained with the latter.
The introduction of inertia and damping coefficients was chosen for its didactic interest, but
the adoption of a linearized thrust as in Equation (9) in a FD model is a simpler and equivalent
approach to include the interaction between nacelle motion and thrust, with due consideration
of controller and rotor dynamic effects.
In any case, a better agreement of the LF response should still be pursued. Other linearization
10
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
1017 1017
SIMA SIMA
FD constant FD constant
FD calculated FD calculated
FD forced FD forced
1016 1016
Sbm (N2m2.s)
Sbm (N2m2.s)
1015 1015
1014 1014
0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2
Frequency (Hz) Frequency (Hz)
(a) Uw = 13.0 m/s. (b) Uw = 19.0 m/s.
1017
SIMA
FD constant
FD calculated
FD forced
1016
Sbm (N2m2.s)
1015
1014
0 0.05 0.1 0.15 0.2
Frequency (Hz)
(c) Uw = 25.0 m/s.
methods could be attempted, trying to preserve effects other than those obtained with the partial
derivatives of thrust and torque w.r.t. relative wind velocity, rotor speed and blade-pitch angle.
11
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
Table 5: Error in tower base bending moment standard deviation compared to time domain
simulations.
0.02
FD forced - w. inertia
FD forced - wo. inertia
0.015
Sx5 (rad2.s)
0.01
0.005
0
20 25 30 35 40 45 50
Period (s)
Figure 6: Pitch response estimated in the frequency domain for Uw = 19.0 m/s, with
aerodynamic damping estimated using the forced oscillations method and with/without
considering the aerodynamic inertia effect.
References
[1] Kvittem MI, Moan T. Frequency Versus Time Domain Fatigue Analysis of a Semisubmersible Wind
Turbine Tower. Journal of Offshore Mechanics and Arctic Engineering. 2015 Feb;137(1):011901–011901–
11. Available from: http://dx.doi.org/10.1115/1.4028340.
[2] Pegalajar-Jurado A, Borg M, Bredmose H. An efficient frequency-domain model for quick load analysis
of floating offshore wind turbines. Wind Energy Science. 2018;3(2):693–712. Available from: https:
//www.wind-energ-sci.net/3/693/2018/.
[3] Hegseth JM, Bachynski EE. A semi-analytical frequency domain model for efficient design evaluation
of spar floating wind turbines. Marine Structures. 2019;64:186 – 210. Available from: http://www.
sciencedirect.com/science/article/pii/S0951833918303149.
[4] Bachynski EE. Design and dynamic analysis of tension leg platform wind turbines [Ph.D. thesis]. Norwegian
University of Science and Technology. Trondheim, Norway; 2014.
[5] Nielsen F, Hanson T, Skaare B. Integrated dynamic analysis of floating offshore wind turbines. In:
Proceedings of the ASME 2016 25th International Conference on Offshore Mechanics and Arctic
Engineering; 2006. .
[6] Larsen TJ, Hanson TD. A method to avoid negative damped low frequent tower vibrations for a floating,
pitch controlled wind turbine. Journal of Physics: Conference Series. 2007 jul;75:012073. Available from:
https://doi.org/10.1088%2F1742-6596%2F75%2F1%2F012073.
[7] Jonkman JM. Influence of control on the pitch damping of a floating wind turbine. In: Proceedings of the
ASME Wind Energy Symposium; 2008. .
[8] Souza CES, Bachynski EE. Changes in surge and pitch decay periods of floating wind turbines for varying
wind speed. Ocean Engineering. 2019;180:223 – 237. Available from: http://www.sciencedirect.com/
science/article/pii/S0029801818315932.
[9] Goupee AJ, Kimball RW, Dagher HJ. Experimental observations of active blade pitch and generator control
influence on floating wind turbine response. Renewable Energy. 2017;104:9 – 19.
[10] Moriarty PJ, Hansen AC. AeroDyn Theory Manual. Denver, Colorado: National Renewable Energy
12
NAWEA WindTech 2019 IOP Publishing
Journal of Physics: Conference Series 1452 (2020) 012040 doi:10.1088/1742-6596/1452/1/012040
13
144 Appended Papers
Paper 3
145
146 Appended Papers
0DULQH 6WUXFWXUHV
-
journal homepage: http://www.elsevier.com/locate/marstruc
! " #$ %
& ' # * , , . ,* *
) * *,
(* ** /0 1 ) * *
, * . . 2
+ ,
3 .
!* ) * * * 2 , ,
4 ,
* 4 $ * * .*
** . * . * *
* * $ * !
* . . * ,! $
* *, *
** * * -, ,, * ,
* * * . . * * ,
* * *
**
(* 5(1$ 6 *2 * * **
7 * * **
* ! * 8 . * * * *
* .** *.* *
* * * ** *. * * * 57)6 !
*. * 9/:
# * * 2 (1$ .
$ 9;: * 4 !, * .* * . (1$ 9<:
** ** * , * * 4.
* * 5$&=6 + ** 2 9>: 2 * , (1$
# * ,
%
( ' 3 ? 5 6
%* 9D: * ! . * * (1$ * 9B:
, * (1$ !, * $ * *
* . * * * * * ** * 9C: *,
* , 3 . ** * * ** * *
= * .* *2 * .* * * $ * 2
F * 9G: * * 4* .
* * ' 9E: , * 2 * *
* . * 4* * * * .
* * - 9/0: . .
7 , * 4 ** * 9/0:
(1$ , $ * * 4*
) 2 (** * 9//: ** * *
(1$ , * 4
$ * *
42 *
(1$ * * * 9/;: ' *
*, ** $ * * 9/<: * * 2 * (1$ # 4
& * 9/>: , * , * (1$
** ** * ! * 5&H"6 $
* ,! .* * . *
*
# * . (1$ *,. * *
* ** - *
2 * * * * * .
$ * * * ! ’
$ .* *
* * * # 2 * ** ,
* * * * * * * # **
4 * ,4* * * * , , . ,* /0 1 (1$ !
$ * * ) 7#) 9/D: * 4 ,
* $ * *
2 !*
* 3 . ( ** * . 4 * , * $
4 * * * * * $
* 9/B–/G:
- * * 4 * *
* 7 * 4 * * ! *
. * *
7 * * ! * * *
$ ** $ * (1$ * * - ;
- < * . * 3 . *
$ * - > * . 4
,4* * * ( ** 4 - D
$ * (1$ ** * * , * *2 ,
* * ( !
.* . . *. * * ** .
- * 4 * ** * * 2 *
* *
$ 7$I /0 1 9/E:
<;0 $ * ** ** 4** ** . ** ;B00 A<
$ * /0 . ** * 5-1&6 * //E . -1& #
* ** . * 2 , ( *
,** 5 * * 6 . * , * $
* * * *
)* ),*, * , * . .** . ),
* ** * # .. (1$ * * ( /
+, , - , 0DULQH 6WUXFWXUHV
$ * ! *. ! $ * 2
, * * * * ,*
, , . ,* * - ;> $ * . * x u
.* * @
x ¼ x0 þ Δx; u ¼ u0 þ Δu: 5/6
$ ! 2
Δx_ ¼ AΔx þ BΔu; 5;6
A 2 B 2 Δ .
" . 2 ! 23 !
+, , - , 0DULQH 6WUXFWXUHV
( ,. * ! * 3
. * * *
! *
#* * * * * * *
* * * * ** . *. * *
) 2 9;/: * * **,* *
!, * * * ,
. * * * * * ! 4
* *
. (1$ . * * *
* 2 * * 9;;;<: $ * *8
* * *
1 * 7)( . . 4
2 3
ξ
xs ¼ 4 ξ_ 5: 5>6
ϕ_
$ * * . * *
. u u $ * .
QG
usc ¼ ; 5D6
θ
θ ** . * * * * ! $
* * ! *
* .*
$ . . 8 . * * *
* * 2 * * . * $ . 2 .
* ( / $ . ** *
* ** 4* .
* * $ . 2
>
usd ¼ vFT vMT vQA FW;1 FW;5 FW;7 ; 5B6
4 2 . * ! .* .
- ;D; 45; * . 2 7)( %–(
$ * @
+, , - , 0DULQH 6WUXFWXUHV
ys ¼ Cs xs ; 5C6
2 * 5M * 6 5K6 5D6
0 I
As ¼ 1 1 5G6
M K M D
$ B B * * * 7)( * 2
* * C * * @
Cs ¼ ½0 0 0 0 0 0 1: 5E6
.,.,.,
$ * ** * 4 !, * &
-. 9;>: * 2 & -ø* 9;D: * ** * 7)( ( <
$ * .
0ðÞ * * * ! , * *
. * ! * 5 4 6 $ 563 6
* . ! * 4 $ * 4
5 6 4 5 6 ! , * 4 $ *
4 5/ 6 *
rffiffiffi
8
c juðtÞj
_ uðtÞ_ c _ ¼ cL uðtÞ;
σ ðuÞ_ uðtÞ _ 5//6
π
σ ð0Þ_ . 0 _ ! 5/06 . ’ ! ,+
2 * * 2
"* * * *.* 2 * * *.* 9;D:@
8
y2
>
> 1 exp 3k2 þ 1 ; 0 y y0
>
< 2
FY ðyÞ ¼ " pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi # 5/;6
>
> 3k2 þ 1 y0
>
: 1 exp y ; y y0 ;
2k 2
1
y0 ¼ pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi: 5/<6
2k 3k2 þ 1
$ 6 * .
+, , - , 0DULQH 6WUXFWXUHV
$ * * * ** . *, , ** 9;B: $
** @ , ! ** * , ** . *,
** ., $ ** ** * *
. 8 . * * * *
* ** 4 ,
x_ c ¼ Ac xc þ Bc uc ; 5/>6
yc ¼ Cc xc ; 5/>6
* ! * *@
_ QG
uc ¼ ½ϕ; yc ¼ : 5/D6
θ
( 4* 4 , * , 4* ! ω . ,! 2
*
" #
0 1 0 ϕlp
Ac ¼ ; Bc ¼ ; xc ¼ ; 5/B6
0 ωlp ωlp ϕ_ lp
ϕ_ * , 4* * ! * !
* ! * , 2
# . ! ! = ** *,
* . * *@
0 2K ϕ_ 0
Cc ¼ ðbelow ratedÞ; 5/C6
0 0
0 0
Cc ¼ ðabove ratedÞ:
ηk k i ηk kp
& ! η6 , * * 6 6 * *
= ** .*
$ * * * * * ,* ys ¼ uc u ¼ yc @
x_ ¼ Ax þ Busd ; 5/G6
y ¼ Cx; 5/G6
As Bsc Cc Bsd xs
A¼ ; B¼ ; x¼ ; 5/E6
Bc Cs Ac 0 xc
C 2 $ * ,* ! $ 2
! ω 4 @
ðÞ- 3 9;C: $ . y * * Sy ðωÞ
$ * 2 5Su ðωÞ6 * **
Swind ðωÞ 0
Su ðωÞ ¼ 5;;6
0 Swave ðωÞ
$ , * . . * 9;G:
$ * * ! * 4 *
+, , - , 0DULQH 6WUXFWXUHV
* $ ** * * 7 * 9;E:
$ 2 . 5#I"6 9<0: #
+ * . 5%7(6 2 .* 9
9;C:
y2
FY ðyÞ ¼ exp νþ ð0Þ T exp 2 ; 5;<6
2σ
νþ ð0Þ , , , * @
rffiffiffiffiffiffi
1 m2
νþ ð0Þ ¼ : 5;>6
2π m0
$ * 2 .* / / * 2 $ .*
2 @
Vt ¼ ϕ_ r ð1 þ at Þ; 5;C6
* . * * .* * ,
* ,** * * * *
$ .* 5 * φ * * 2
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
W ¼ V 2n þ V 2t ; 5;G6
Vn
φ ¼ tan1 : 5;G6
Vt
$ * * *’ *@
α ¼ φ ðθ þ βÞ: 5;E6
β * $ * * * 4
1
L ¼ ρcW 2 Cl ðαÞ; 5<06
2
1
D ¼ ρcW 2 Cd ðαÞ: 5<06
2
$ * * * . *
Fn ¼ L cosφ þ D sinφ 5</6
+, , - , 0DULQH 6WUXFWXUHV
Z R
Fy ¼ Fn dr; 5<;6
0
Z R
Mz ¼ r Fn dr; 5<;6
0
Z R
My ¼ r Ft dr: 5<;6
0
# ! 5;C6–5<;6 * * * $
* * * 4 (1$ * * * 4 $*
2 @
* . . . *
∂Ax
Ax;y ¼ : 5<>6
∂y
% ðωÞ ψ * * & ( 4
γ ( 2 *@
Z
1 π
KnðjkÞ ðωÞ ¼ γ ω; djk cosðnθÞ dθ: 5<B6
π 0
$ ! 5 6 % 9<>: 2
!
2 2 0:5
fd 0:12d
γðf ; dÞ ¼ exp 12 þ ; 5<C6
Uhub Lc
, / *
$ * * ** * *
* $ * . . * * ! 5</6
1 * * * * ** *
* * ! 5 <= B= E= 6 *
<= * $ *
$ . ** * * . @
SF ðωÞ
SvFT ðωÞ ¼ T 2 ; 5<G6
3 Fy;v
SMT ðωÞ
SvMT ðωÞ ¼ 2
; 5<G6
3
2
Mz;v
SQ ðωÞ
SvQA ðωÞ ¼ A 2 : 5<G6
3 My;v
$ * 2 ! 5;;6 S ðωÞ * .
+, , - , 0DULQH 6WUXFWXUHV
* # * 9<<: ! *
* ** ** * * $ . ,
3
MT;v ¼ Mz;v ; 5<E6
2
% * * C A 5*6 /D A 5 6
+, , - , 0DULQH 6WUXFWXUHV
$ * * . . * *8 .
** . * * ** * !** . *
* $ * * 2 @
.,@,
$ * ** *. . * 5% 6 . * 8 . *
* <0 % .* * , $
/ <0 A / A ( * .* 4 .
* $ .* ** * 3
* * 9<E: ( D $ * * . *
*
X
NEC
pi ¼ 1; 5>/6
i¼1
% % * *
- 4 . * * * *
**
+, , - , 0DULQH 6WUXFWXUHV
$ % $* / * D0, * .* 2
% / * * 2 % ;
* ., * % < . ,
*
( ** % F * % * * * * 5'$6 9<>:
$ * *@
/ - *
; 1
<
> *, *
7,, 4 6
$ (1$ * * ) 7#) 9/D: , * * *
# * * ** 5 6 *
*2 * * * * . . .
* * ** * . .
9>0: $ * * * 4 *** =
$ , * . . * * #*
! 4 * * 4
4 , 2 $ * . . 3
. * 9>/: 1 * G0 /G
* > * * ,. 2 * D * $
* * ** 2 2 5K7-6 9>;: ( / # 2
$ -')=$ * 9><: *. * -')=$ ! * ! 5-H=6
*. * * ) 7#) ) - 9>>:
, . * ** * * 3 . * ,
* .* *
7,., A!=
** (1$ * *.* 5&%)6 *
. * . &%)
* * * * 3 . #
*
* * 9/0: * * $
** * * * . (1$ #
*. !* * .
,!* * . * * .* . *
. $ 3 . @
8 * . * . . σ ðϕÞ
_
. * 2 **
% / ; <
+, , - , 0DULQH 6WUXFWXUHV
$ * . ** / ; / þ ; ¼ /; / ;; 2 ½0; /
!* 4* ! .* ,
, 3 .
$ *
Ctotal ¼ Cspar þ Ctower þ Cmoor ; 5><6
8 8 8 * .* $ *
* * ( 9>D: % * **
* $ * 5 ** 6 2
X
Cspar ¼ km Mspar þ kf Ti ; 5>>6
i
6 * * ** 6
2 $ * * * **
* * $ * 6 =6 *
. /0 * .* 1 * 9>D:
$ ** * 9>B: * * G0 ** 5 6 $
** *
* * /0 €A
5 * * 9>C:6 >D €A 5 * * 9/>:6 * 6 .
.* ;C €A 2 * .* (** * 9//: $ *
2 <–>D €A 9>G: * ">
<D €A * * D–/0J $
* . *
, 8 . * * * .
$ . . ! 5>;6 .*
<0 (&- - ;B * @
X
NEC
_ ¼
σðϕÞ _ ;
pi σðϕÞ;i 5>D6
i¼1
σðϕÞ;
_ .
-*
+, , - , 0DULQH 6WUXFWXUHV
7,7,.,
- ** * * * * * $
* 42 . * 7$I /0 1 //E
. -1& * /0D 2 * ///< $
. * ** ,
** * GD00 A< . *
9/E:
7,7,7,
$ * . * @ * * *
-1& * * * * * * * ;7
, * . * * *
$ * * $ !
2 * .* * 9>E:@
2
7 . *
7 - * ' . * I
- //
- ** //
- * / /0
" ; >
" ; >
" ! >
" >
" >
* /
( * * /
* * / /
* * ) /
+, , - , 0DULQH 6WUXFWXUHV
7,7,:, %0
% . * . * * * //0
. * $ ! * . * , *
, * ,4 * ¼ > ,
$ . * . *
* . $ * G0 . * *
$* ;
$ * * 5(&-6 * 5I&-6 *
. , * * * * * $
2 * * 3
;0 ;0
** * 7'L +& 9DB: ** * @
σ=;- ! .* . ** * $ . * * .*
! $ *
2
Astiff þ 0:06 lstiff tspar ; 5D/6
Zl2
B . #* . "
Istiff Ix þ Ixh þ Ih ; 5D;6
! * * ** 3 2 *
5") 6 5") 6 2* 5" 6 ** *
* * ** 2* * *,
* 9D>: $ . * * * *
+, , - , 0DULQH 6WUXFWXUHV
sffiffiffiffi
E
hstiff 1:35 tw;stiff : 5D<6
fy
M ’ * , * ** $ * <DD = $ . *
* ** ! * 4
2
kLC E
2:5 ; 5DD6
iC fy
6 . * 5 /06 /8 * * * 8 *
( ** * * γ 4 ¼ /:<D * . * * 9D>: * * γ .
** ! 9D>:@
8
< 1:10; λs < 0:5
γM ¼ 0:80 þ 0:60 λs ; 0:5 λs 1:0 5DB6
:
1:40; λs > 1:0;
7,,.,
$ * 7 . 7(( ¼ ;:0
9DD:
$ * < 9B>: *
* ( ** ! 4
σcr
σx ; 5DC6
γM γF
σ ) 2 * σ * * γ4 γ * *
.* $ * γ ¼ // 9BD: γ4 ¼ /<D 9D>: $
* ** !* ** * #
** F-
+, , - , 0DULQH 6WUXFWXUHV
7,,7,
$ * @ 2 * * 2
* * * $ 4 2 9D>:
Td 0:95 Smbs ; 5DG6
; ; D0, ,
.* γ ¼ /< γ ¼ /CD * * * 9D>: $
$ * ** . 2 (
* * * * $, , 3 . ! 5>;6
* - < - !* *
* * **
* ( ** . 4 * *
$ 4* * ! 4 ,
* ! ** . <= $
. ** * ** ( C 7 “ ”
* ! ** * . . , * ,
. .
* * ,
* . * * * 4 4
+, , - , 0DULQH 6WUXFWXUHV
, #* * ** (1$ , *
* 2 $ * * 4 * % * ,
** . * .
* 2* . * , * .
. * * * * ,
*
$ , !* .*
! 5>;6 . * .* . *, 7$I /0
1 * ** 9C0: * * -# $ *,
. . 00>< A
* 4 ** . 2 * *, .*
$ * ( G 4 . . * .* **
* . * . * . 2
$ . * . *, ** *
* * * , D0, 5% /6 * * $
*, ** 2 , *
., ** 5 6 * *
!* *
** * * * ( G / ¼ C
/0G ; ¼ ð/ C /0G Þ ! 5>;6
$ * * * * ( E ( ** .*
* *
$ * * .* ** * . . *
** 2 * . *
* * ** (&- ( * *
2 * /0 * .
* . ** 2 * 2 .* 2 * <;
DB * -1& .* $ * 2
* 2 * . * $
! * 4 * * .* **
3 $ * * . *2
! . * * * *
$ * . .* 4 * * $ * 3 * .
! $, . * 3 .
+, , - , 0DULQH 6WUXFWXUHV
" ) * * 5*6 * * 5 6 $ ** *
>0 * . . *
+, , - , 0DULQH 6WUXFWXUHV
5% 6 * . **
% * . 5% +6 .* *
$ 2 * ( ** 4**
** * ! * !
* ! ** * * * % + -
* . * .
.,! ) * # ,
! * 5* *6 * 2 * ! **
$ * D 1 )%<, 9C/: . *
. * ! ,** #* * **
3 . ** * . * *
** *
2 * /D 2 <; D0, .
$ * % / * $
* CC $ * * $* <
* * * * $* > 7 *
. $* /0 // # 2 #
- * 2 ;0D . -1& $
** ** ** $
. * *
* * * * * * , .
. ! . **
. * * * *
. * $ * **,, *
!* *,, * ** 5 . * 6
** 4*
$ * ** . * * *
. * I* *
% !* * 4 (&- I&-
* * * * * $* D * *.
42
( /0 )),- 1 (* - /0 1
&(-D0þ 3 9C;: $ * * ** **
)),- % .** 2 *
;>J * )),- * * * . 5GD00 . G;>< A<6 $ * !
4 , 5 2 * 6 0G 0DE )),-
.* * ** , ) *
* $ *
* * <= 2 * *
$ * ( // #** 4. . * **
* ** ! * * * .
* ** *
$ * * * !* ( /;
* ** * * .
$ ** (1$ *
* 6 6 42 * .* $* B $ .*
9;0: * * 7$I /0 1 ** 9C0: *
) * 5 6
= L* I
7 G0>
* ** /:DC /0C
7 * /:C; /0> <
% + >:BB /0E ;
L * % + * -1& BB0
L * % * -1& >;E
+, , - , 0DULQH 6WUXFWXUHV
'* (1$
L* I
- />>C
= <>>
. ;D0
/ /;D
)
= L* I
* * * ! *
$ !* * ** 5 3 . 6
$* C I * ** * * * * *
* . #* * .* * .*
* 2 * * ** * * 2 . *
*
+, , - , 0DULQH 6WUXFWXUHV
) $ * * *
) *, ** * * .
( 2 *
= L* I
)3 . .* ** *
- 9/0B €: ' * . 9:
+, , - , 0DULQH 6WUXFWXUHV
,:,
$ * $* G $ * I&- ,
% / 2 * - >< $ *
$ * * .* * * ! * $
* * * * 4* . ** * 8 .
2 . ** 9C<: * *
** * , . *
* 4 * * *
* * #* * * ,
!* *.* 2 * .
* $ *
* * * , *
* $ * * **
* * * # * * ,! *
* ** * ,! * 9C>: . *
$ * * * *
* - (1$ * 9CD:
. 4 *
* . . * **
*
$ . 4 * * * , ,
. ,* , * -# ** * ** * * ,
* * * * %–( .
’ ! . * *
* * ( –F* .
* - . * * 4 . 2 . **
*
* * 2 * * * * *
’ ! # * * *
/:>G /0G 'A $ *
- ;D> . *
$ ** * <0 (&- 2 /, * "
–=* * . .* . 2
$ 2 * 2 .* / D0, #I"
4 * 9CB: /, * ;0
D0,
$ * “- /” 9;0: B0J **
>0J * ** * 4 , $ (1$
* * $ * . 4
* * *
( * * * , * ( /< ( /> *
* .* ( * .* . * * .
* * * $ * * 7 * *,
* *, * -' . .* ,* $ *
* .* * 4 * ,* $ *
* . ** 2 * .* *
( * * * ** *
* * ** . . $ * . *
* * 2 * *
** * *
* * ( /<
$ * * * * ( />
* , * -# 1 * . !
,+ 7 * (
* + * * 7 *
(1$ 9;0:
+, , - , 0DULQH 6WUXFWXUHV
!
)
= L* I
00G/
* /0:C
/ /><<0
) />EC;
#* . ** * . * *
$ * *
* ** * **
$ 2 .* * * $* E $ * * .
* * * * * * $
* * * >J 2 .* *
. * *
$
# * , , . ,* (1$ * .* *
*, ** /0 1 (1$ $ *
. !*
$ * * * * 2 , ,
4 * * * , ** .
(1$ * 4 (&- 4 . *
. * * . . , ** * * *
*
+, , - , 0DULQH 6WUXFWXUHV
"
% 2
* * )
& <D/ C0 B:/D /0B ' <;0 /D0 D:E/ /0B '
' * ;E> DE B:BG /0B ' ;C/ //< B:BG /0B '
$ * * * . ** * . *
$ * 2 % % +
* ! !* . . * ,! * $
* * . * * * ! .,!
$ . * *, *
* ** ( * . *
* L * , 3 . * * . *
. , . 5 * 6
) 42 * * 0;J /0DJ *
. *
$ * . *
. 2 . * 4 *
* * , 3 . * . * . #
- <; 4 , 3 . ** . *
, . * $ * *
$ * * ! * *
+, , - , 0DULQH 6WUXFWXUHV
$ * * *
#* * * * * * *
# * * ! . *
* * ** 2 * *
, . *
$ * . 4 , ,, * * 4 * .
* 2 * * * . . **
* * * * $
. * * * *
$ * * * . *
* $ .* * ** *
. .** . $ * 2 * *
* .* *
N $ * . 4 * * * * .
' ()
( * " % * ' 5'("6 3 ;C>G;C , “+ -” **
* $ 4 * * * ) 7#) *
' % * & + *& )
#
7 *
7
+, , - , 0DULQH 6WUXFWXUHV
, -%./ %
, K7- * * " * 2 2 * *
*** *** 42 * $ .
* * * * * - &
9>;: * 4
&
9/: ""# & # * * @ . ### ;0/<8D/5E6@;0>E–CD @AA A/0;D/>A/0D/GED
9;: $ % = ’ $ * 8 ;00C
9<: $ 7 * * - ;0/;8;E@GE–//> @AA A/0/0/BA3
;0/;0E00/
9>: + ** 2 ,% = = * * ,* * * , @ =
, 5;0/>6 * * 5-)=;0/>6 F @ 8 ;0/>
9D: %* +( & * # * ) " /EEB8/G5>6@/DC–C/ @AA A/0/0/BA-0/>/,//GC5EB6
000;G,>
9B: F & F * ! * = ;0/;8
;>@;GE–EB @AA A/0/0/BA3 ;0/;0B///
9C: ** % % . *. @ * , @
)%#'- ;0/< $-A @ ** ' 8 ;0/< @AA A/0//0EA)%#'-,;0/<BB0G/C<
9G: F ** % % # * , 3 . )
;0/C8<5/6@BE–GC @AA A/0/00CA >0C;;,0/B,00C;,>
9E: # ' $# & ,*, * @ = , 5;0/;6
* * 5-)=;0/;6 " +8 ;0/;
9/0: * - - 7 , . ) 1 ;0/>8/5/6@/;–;;
9//: (** * =# 1'7)=$, * @ = #- ;0// <0
* # 5)#;0//6 " $ '* 8 ;0// @AA A/0///DA)#;0//,
>EEGD
9/;: -,- ** * 7 @ = , 2 5;0/B6 * *
5-)=;0/B6 " +8 ;0/B
9/<: $ * % &% ) * * * ** = % ;0/>8DDD5/;0EE6 @AA A/0/0GGA/C>;,
BDEBADDDA/A0/;0EE
9/>: & ( O** F M 1 -* 7 % =1 ) * *, ** @ =
#- ;0/C <B * # 5)#;0/C6 ' 8 ;0/C @AA A/0///DA
)#;0/C,B;0<G
9/D: + - $ ""# F$ '* # ) 7#)@ , * * *
- * ) ;0/E8DE@/0CD–/0> @AA A/0/00CA 00/DG,0/E,0;;//,
9/B: # $ P 3 ""# . * +1 . F +# * * *.*
" ;0/>8BG@GE<–E0D @AA A/0/0/BA3;0/>0;0>D
9/C: # $ ""# P 3 . F +# . * +1 # . * 4 ;0 1 *
1 ;0/B8/E@;0C/–GC @AA A/0/00;A/EC0
9/G: # P* ( -ø '' ""# * ,4* %(7, /0 1 1
- ;0/E8>@/B<–E; @AA A/0D/E>A ,>,/B<,;0/E --' ;<BBC>D/
+, , - , 0DULQH 6WUXFWXUHV
9/E: % P* ( " M # &% '3 # 7 7$I /0 1 $ * " 7$I 1
" ,,00E; 7$I 1 8 ;0/<
9;0: # ,* * ! * 4 .* - ;0/E8B>@
/GB–;/0 @AA A/0/0/BA3 ;0/G/00/D
9;/: F % * **,* =7 ' I . - $ * 8 ;0//
9;;: F. & % $ 7 * * *
= ;0/<8<D@;/0–;; @AA A/0/0/BA3 ;0/<0C/C>
9;<: P F $ 7 * * * * . 1 ;0/>8
/C@/<GD–>0E @AA A/0/00;A/B<E
9;>: & F -. =% 4 ** * @ = 4 5/EE06
' 8 /EE0
9;D: & -ø* ' - * 4 * 2 * @ ) # * ; # * * *
2 * # * 8 /EE<
9;B: 4* - * 1 - + 74 D,1 .* $ * " '"&A$=,D00,
<G0B0 ' * "* & 8 ;00E
9;C: ' # $ - % I . = 8 ;0/<
9;G: 7'L +& & $ * " 7'L+&,-$,0><C 7'L +&8 ;0/B
9;E: 7 * $ # * * =7 I . 1 8 /EGD
9<0: ' # + ) %* 2 * ;00G8/<>5G6@B;G–<B @AA A/0/0B/A
5#-%60C<<,E<EE 5;00G6/<>@G5B;G6
9</: $ ' - 7 1 1 *8 ;0//
9<;: ** ( + + " #& 1 2 * 1 *8 ;00E
9<<: * # 7 * ! =7 I . % ** & 8 /EEG
9<>: % 1 , /@ ! $ * " % B/>00,/ * * * % 8 ;00D
9<D: -'$( ) "(&K ;0/B
9<B: -'$( ) -) ;0/B
9<C: F * & $- ’ @ . /D0 $ * " '"&A$=,D00,>B/EG ' * "* & 8 ;00E
9<G: % 1 , <@ ! $ * " % B/>00,< * * * % 8 ;00E
9<E: F * $- . - . ' ' - ) = * ;00;8/;5/6@/–G
9>0: $ ""# # * * * * * . . #% $
- ;0/G8>>5>6 @AA A/0//>DA</G;<E<
9>/: ""# $ ". 4 . . * * * * ### ;0/<8D/5//6@
;DG;–EE @AA A/0;D/>A/0D;/G> --' 000//>D;
9>;: & # ""# 2 2 * * * -
* ) ;0/;8>B@;C<–G> @AA A/0/00CA 00/DG,0/;,0CB<,
9><: + ** = 1 - # -')=$@ -H= * *, * -# ) ;00;8/;5>6@ECE–/00B @AA A
/0//<CA-00<B/>>D0>>>B0EB
9>>: = " =1 ""# ) @ # = , 3, * - * ) ;0/;8
>D5/6@/0/–/G @AA A/0/00CA 00/DG,0//,0BBB,<
9>D: ( F ) * - 8 ;0/<
9>B: $ ** FO* F + * L 7 M * ( + " . - ' ( ' * L - " + " !
+ & % # ' # 7 * ,% = * + ( 1$($ . . *
* 7;<@ ** * ** * * * * 1.%
) "* 8 ;0/B $ *
9>C: # 3 % Q # ' $# &.* * * . " ;0/>8BB@
C/>–;G @AA A/0/0/BA3;0/>0/0/C
9>G: & F = * ;0/E
9>E: * ' (* @ * ) *4* =* 8 /EEG
9D0: & $ $7 # . . * ! . ** = % - ;00C8CD@
0/;0C< @AA A/0/0GGA/C>;,BDEBACDA/A0/;0C<
9D/: . L + % ") % * @ ;0/; # * * %@ ( H
* 8 ;0/; </>G–D< @AA A/0//0EA;0/;B</D/;0
9D;: & % ** * * * 1 ;00E8<<5B6@D>/–D< @AA A/0/;B0A0<0E,
D;>K<<BD>/
9D<: 7'L +& ( * $ * " 7'L+&,"=,%;0< 7'L +&8 ;0/E
9D>: 7'L +& (* $ * " 7'L+&,-$,0//E 7'L +&8 ;0/G
9DD: 7'L 7 $ * " 7'L,)-,/0< 7'L8 ;0/<
9DB: 7'L +& * ** $ * " 7'L+&,"=,%;0; 7'L +&8 ;0/E
9DC: - # * - F. 7 % , *
* @ = , 2 5;0/B6 * *
5-)=;0/B6 " +8 ;0/B >/E–;B
9DG: '3 #" $ 2 * * . , ) # ;0/E8/>/
@AA A/0///DA/>0>/EEB
9DE: 7'L +& ** * $ * " 7'L+&,"=,'/0< 7'L +&8 ;0/C
9B0: F F "** " * * ) ;00>8</@;/CD–;0G @AA A/0/0/BA3
;00>0>00D
9B/: * # (* ) #* . * * @ ) * 8 /EEE @AA A/0>0><A
/0ED<,
9B;: F * + - " - * . 2 @ * . *
, * P - *8 /ECE @AA A/0/0/BA />C>,BBC05/C6BDDG>,G
9B<: & # F + ""# # .* 2 - * ) ;0/C8DD@
;DC–CC @AA A/0/00CA 00/DG,0/B,/>ED,/
9B>: % - <@ * /,B@ * ** $ * " ' /EE<,/,
B@ ;00C % - 8 ;00C
9BD: 7'L +& - $ * " 7'L+&,-$,0/;B 7'L +&8 ;0/G
9BB: 7'L +& ) $ * " 7'L+&,)-,<0; 7'L +&8 ;0/G
9BC: 7'L +& = $ * " 7'L+&,)-,<0/ 7'L +&8 ;0/G
9BG: % %,$ & )2 I . = 8 ;0/<
+, , - , 0DULQH 6WUXFWXUHV
9BE: - ( -* 7 7 % =1 @ = #- ;0/> << *
# 5)#;0/>6 ( %* I-#8 ;0/> @AA A/0///DA)#;0/>,;>;>>
9C0: &% 7$I ** $ * " 7$I 1 " ,,00/G 7$I 1 8 ;0/<
9C/: 74 L )%< $ * " '"&A$=,D00,>CD<D ' * "* & 8 ;0/0
9C;: O** F & ( M 1 &(-D0þ 7>;@ * 4 &(-D0þ /01 $ * I . -8 ;0/G
9C<: * 1 4* - # ( * * @ ;< #- " '. I-#8 ;00>
@AA A/0;D/>AB;00>,/00C
9C>: & + P $ * * 4 * @ = )# ;00C ;B *
# 7 %* I-#8 ;00C @AA A/0///DA)#;00C,;E/DD
9CD: - ' * (+ $7 M. " .ø** ) " * # #* * 7
1 ;0/D8/G@//0D–;; @AA A/0/00;A/CD0
9CB: - ' ' # %* , 2 .* ;0/08/<B5/;6@/>E/–D0/ @AA
A/0/0B/A5#-%6/E><,CGGE0000/E>
Paper 4
175
176 Appended Papers
Design Optimization of Spar Floating Wind Turbines
Considering Different Control Strategies
John Marius Hegseth1 , Erin E. Bachynski1 , Joaquim R. R. A. Martins2
1
Department of Marine Technology, NTNU, 7491 Trondheim, Norway
2
Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109, USA
E-mail: [email protected]
Abstract. One of the challenges related to the design of floating wind turbines (FWTs) is the strong
interactions between the controller and the support structure, which may result in an unstable system.
Several control strategies have been proposed to improve the dynamic behaviour, all of which result in
trade-offs between structural loads, rotor speed variation, and blade pitch actuator use, which makes
controller design a challenging task. Due to the interactions, simultaneous design of the controller and
support structure should be performed to properly identify and compare different solutions. In the present
work, integrated design optimization of the blade-pitch controller and support structure is performed for a
10 MW spar FWT, considering four different control strategies, to evaluate the effect of the controller on the
structural design and associated costs. The introduction of velocity feedback control reduces the platform
pitch response and consequently the fatigue loads in the tower, which leads to a decrease in the tower
costs compared to a simple PI controller. Low-pass filtering of the nacelle velocity signal to remove the
wave-frequency components results in reduced rotor speed variation, but offers only small improvements
in costs, likely due to the limited wave-frequency response for the considered designs. Comparisons with
nonlinear time-domain simulations show that the linearized model is able to capture trends with acceptable
accuracy, but that significant overpredictions may occur for the platform pitch response.
1. Introduction
For floating wind turbines (FWTs), the performance of the control system generally depends on the
support structure design and vice versa. A well-known interaction between the blade pitch controller and
platform motions is the introduction of negative damping above rated wind speed [1]. While detuning
the controller gains such that the bandwidth is reduced achieves stability, this results in poorer rotor
speed tracking performance [1, 2]. Several alternative methods have also been suggested to resolve the
issue, such as introducing a feedback term proportional to the pitch velocity [3] or nacelle velocity [4, 5]
to manipulate the generator speed reference. For these types of feedback control, Fleming et al. [6] also
suggested to remove the wave-frequency components from the velocity signal, as this reduced the tower
loads. All of these control strategies result in trade-offs between structural loads, rotor speed variation,
and blade pitch actuator use, which vary with different environmental conditions. Therefore, identifying
optimal control parameters is a challenging task.
Due to the strong interactions, simultaneous design of the controller and support structure should be
performed to have a fair comparison between different solutions. A step toward simultaneous design was
made by Lemmer et al. [7], who optimized the main dimensions of a three-column semi-submersible
FWT. They minimized a combination of material costs and damage-equivalent loads in the tower for
seven operational conditions, with a constraint on the static pitch angle at rated thrust. The controller
was tuned at each design iteration using a linear quadratic regulator (LQR) approach. To properly
identify and compare optimal solutions, the integrated control and structural designs should, however,
be evaluated over the lifetime of the system, considering actual design limits. The purpose of the present
study is to perform integrated design optimization of the blade-pitch controller and support structure for
a 10 MW FWT, considering long-term fatigue damage and extreme response constraints, to evaluate the
effect of different control strategies on the structural design and associated costs.
where ξn represents generalized support structure DOF n, and ϕ̇ is the rotor speed.
Hydrodynamic excitation loads on the hull are described by MacCamy–Fuchs theory, and transverse
added mass is based on analytical 2D coefficients. Radiation damping is neglected, while viscous
damping is found from stochastic linearization of the drag term in Morison’s equation.
Wind loads on the rotor are derived from linearized BEM theory with the incoming wind field
described by the Kaimal spectrum and an exponential coherence function for the longitudinal wind
2
velocity component [11]. The blades are considered rigid in the model, and the aerodynamic forces
on the rotor are applied as resultant loads at the tower top. In addition, the static component of the
aerodynamic quadratic drag force on the tower is included.
The inputs to the structural system consist of both control system outputs and disturbances due to
environmental loads. The control input vector contains the references for the generator torque (QG ) and
the collective blade pitch angle (θ), and is defined as
The disturbance vector contains rotor-effective wind speeds for thrust, tilting moment and aerodynamic
torque, and generalized wave excitation forces for each support structure DOF, i.e.
where ϕ̇0 is the nominal rotor speed reference, kf is the velocity feedback gain, and ẋf is either the
platform pitch velocity or the nacelle velocity. An updated expression for the rotor speed error can then
be established as
Δϕ̇ = ϕ̇ − ϕ̇0 = Δϕ̇ − ϕ̇0 kf ẋf , (7)
where Δϕ̇ is the nominal rotor speed error. In CS4, the nacelle velocity signal is passed through a first
order low-pass filter to remove the wave-frequency components before it is fed back to the blade-pitch
controller.
N EC
Condition 1 2 3
Mean wind speed at hub height, U (m/s) 13.0 21.0 50.0
Significant wave height, Hs (m) 8.1 9.9 15.1
Spectral peak period, Tp (s) 14.0 15.0 16.0
Turbulence intensity at hub height, I (-) 0.17 0.14 0.12
The costs consider both material and manufacturing, using the cost models developed by Farkas and
Jármai [19]. Costs related to installation, maintenance, and decommissioning are not included. The cost
of the platform (and similarly of the tower) is expressed as
Cspar = km Mspar + kf Ti , (13)
i
where km is the steel cost factor, Mspar is the steel mass of the hull, and kf is the fabrication cost per unit
time. Ti is the time spent at fabrication stage i, expressed as a function of the geometry. The steel cost
factor, km , is assumed to have a value of 2.7 e/kg, while the ratio between the material and fabrication
cost factors, km /kf , is set to 1.0, which is a typical value for West European labour [19]. The cost of
the concrete ballast is neglected in the current work.
3.3. Constraints
The fatigue damage at each tower node is evaluated using an SN curve approach, where the D curve in
air from DNV-RP-C203 [20] is applied together with a design fatigue factor (DFF) of 2.0 [21], and the
lifetime of the FWT system is chosen to be 20 years. The fatigue design constraints are thus expressed
as
N EC
1.0
Dtot = N20 pi D i ≤ , (14)
DFF
i=1
where Dtot is the total fatigue damage in 20 years, N20 is the number of short term conditions in 20
years, and Di is the fatigue damage in condition i.
Tower buckling is assessed using Eurocode 3 [22], assuming that the tower is stiffened between each
section to reduce the buckling length. To ensure a smooth transition between the platform and tower,
the tower base diameter is set to be equal to the diameter at the platform top. Both fatigue and buckling
constraints are aggregated using Kreisselmeier–Steinhauser (KS) functions [23].
5
The maximum platform pitch angle in the considered 50-year conditions is limited to 15◦ . Although
the heave response is not included in the model, heave resonance in the wave frequency range is avoided
by placing a lower limit of 25 s on the heave natural period. The added mass in heave is approximated
as the value for a 3D circular disc with the same diameter as the platform bottom [24].
The presented model is valid strictly for hull sections with vertical walls, and a maximum taper angle
of 10◦ is therefore applied as a constraint for each section of the platform, to avoid shapes where the
physics are not captured correctly. Offset constraints are not considered, as the surge response is mostly
governed by the rotor and mooring system design.
Appropriate upper limits for the rotor speed variation and blade-pitch actuator use are difficult to
quantify. The constraints are therefore based on values taken from the land-based DTU 10 MW wind
turbine with the original controller [25], where the weighted average rotor speed standard deviation and
ADC are found from nonlinear time domain analyses using the simulation tool SIMA. Initial analyses
found that the rotor speed variation obtained with the land-based turbine was unrealistic for the floating
system with the simplified controllers considered in the present work. To enlarge the feasible region of
the design space, the constraints for both the rotor speed variation and the ADC are scaled by a factor of
1.5 compared to the land-based values, as shown in Table 2.
4. Results
4.1. Optimized designs
The optimized support structure design for CS1 is illustrated in Fig. 3. The hourglass shape
taken by the platform below the wave zone increases the
distance between the center of buoyancy and the center of
gravity, which leads to increased pitch restoring stiffness, while
the relatively large diameter at the bottom results in larger
added mass and consequently longer natural period in heave.
For the upper part of the platform and intersection with the
tower, the optimizer finds a balance between a small diameter,
which is desirable with regards to hydrodynamic loads, and a
large diameter, which (together with a small wall thickness) is
the most cost-effective way to achieve the required fatigue life.
The optimized tower diameter and wall thickness distribu-
tions for the different control strategies are plotted in Fig. 4. All
four solutions follow the same trends, and the effect of velocity
feedback control is most visible for the wall thickness, where
the values for CS2-4 are approximately 20 % lower than for
the simple PI controller (CS1) along most of the tower length.
The reduced wall thickness is enabled by a decrease in the fa-
tigue loads for these controllers, which is the design-driving
constraint for the tower. For all four control strategies, the 15◦
pitch angle constraint is also active; however, as the 50-year
storm condition with parked turbine (EC3) is found to be the
critical load case for extreme response, this constraint is not
affected by the controller.
The improved fatigue performance can be understood by
examining the response spectra for the optimized designs. In
Fig. 5a, the tower base bending moment spectra are shown for a
mean wind speed of 15 m/s. Large differences are seen around
the pitch natural frequency at 0.15 rad/s, where the velocity
feedback controllers increase the aerodynamic damping and
Figure 3: Optimized support structure thus reduce the response. The low-pass filtering of the nacelle
design for CS1. The wall thickness in velocity removes the wave-frequency range from the signal,
the tower is scaled by a factor of 40 which results in a higher optimal velocity feedback gain and
relative to the diameter for illustration nearly eliminates the tower base bending moment response
purposes.
6
Table 2: Land-based and applied constraint values for the rotor speed variation and blade-pitch actuator
use.
arising from resonant pitch motions. The wave-frequency bending moments, on the other hand, is
unaffected by the control strategy. The blade pitch spectra in Fig. 5b show how the larger PI gains in
the velocity feedback controllers result in overall increased actuator use, with the exception of the pitch
natural frequency, as well as the wave-frequency band for CS4.
The optimized costs of the tower, platform, and tower plus platform for CS2-4 are shown in Fig. 6a,
compared to the optimized costs for CS1. The majority of the cost reductions come from the tower, due
to the improved fatigue behaviour, whereas the platform is less affected due to the fixed costs related
to buckling resistance and the 15◦ pitch angle constraint. However, because a lighter tower results in a
lower overall center of gravity, the platform pitch response is somewhat improved, and a small reduction
in platform costs of about 2 % is also observed. Because the platform accounts for 70-75 % of the total
costs for the considered designs, the resulting total cost reduction is approximately 6 %.
The resulting rotor speed standard deviations and ADCs, normalized by their maximum allowable
values from Table 2, are shown in Fig. 6b. For each control strategy, there exists a limit where no further
reduction in cost can be achieved by increasing the actuator use. This limit is higher for the velocity
feedback controllers than for a controller using only the rotor speed error as input, which causes the
ADC constraint to be inactive at the optimum for CS1. A larger ADC may result in higher probability
of fatigue failure for the actuator bearings, and therefore more detailed design considering the lifetime
of the system should be performed to determine appropriate values for this constraint.
For the rotor speed variation, better performance is achieved with CS4 than with the other control
strategies. Since there is a trade-off between rotor speed variation and structural loads (and thus costs)
as previously discussed, it is expected that larger cost reductions can be achieved with CS4 if the rotor
speed constraint is tightened.
Some limitations to this work should be noted. The standard deviation in steady-state conditions is
used as the only measure of the rotor speed tracking performance, and extreme rotor speed excursions
due to gusts are not considered. The effect of the controllers on surge motions, drivetrain response,
mooring line tension, or blade response has also not been studied. In addition, the performance of the
FWT system could likely be further improved by also adding individual pitch control or modifications
to the torque controller, which have not been examined in the present work.
7
(a) Tower base bending moment. (b) Blade pitch.
(a) Costs relative to CS1. (b) Rotor speed variation and ADC values.
approximately 0.5 % reduction in the costs of the tower, or 0.2 % in total costs. Although small coupling
effects are seen for the considered optimization problem, larger differences are expected in cases where
the controller has a greater effect on the structural response.
4.3. Verification
The optimized controller designs are verified through fully coupled nonlinear time-domain simulations
using SIMA, where two different ECs above rated wind speed are simulated with each control strategy
for a specified support structure design. Comparisons with the linearized model for different response
parameters, weighted by the probability of each condition, are shown in Fig. 8.
The linearized model is seen to mostly follow the trends observed in the nonlinear simulations, but
some errors are present. The largest errors are observed for the tower base fatigue damage, which
is significantly overestimated for CS1, whereas good agreement is obtained with the other control
strategies. The reason for the poor agreement is that the aerodynamic (and thus the overall) damping
for the platform pitch mode is much lower with this control strategy. Consequently, the resonant pitch
response becomes very sensitive to the presence and amount of additional damping in the system, which
is either not considered or underpredicted in the linear model. This disagreement was also observed for
a linearized model of the 10 MW OO-Star semisubmersible in Souza et al. [28], which used a control
strategy similar to CS1. The overestimation of fatigue damage means that the optimized tower design
for CS1 is more conservative than for the other control strategies, suggesting that the cost reductions
in Fig. 6a are highly optimistic, and that CS1 may yield a fatigue design similar to that with a nacelle
velocity feedback controller. It also suggests that future optimization studies using the linearized model
should consider more advanced control strategies than CS1, to limit the pitch response error.
In addition, the rotor speed standard deviation is consistently underestimated by about 20 % for all
four control strategies. This disagreement could be taken into consideration in the optimization process
by adjusting the constraint value.
5. Concluding remarks
The design of the platform, tower, and blade-pitch control system for a 10 MW spar FWT was optimized
simultaneously using a linearized aero-hydro-servo-elastic model and gradient-based optimization with
9
analytic derivatives. The goal has been to minimize the material and manufacturing costs of the support
structure, with constraints on tower fatigue damage and buckling, extreme platform pitch motions, rotor
speed variation, and blade-pitch actuator use, considering four different strategies for the blade-pitch
controller.
The effect of the controller on the structural response was limited to the fatigue damage in the tower,
since the storm condition with parked turbine was found to govern the extreme responses of the system
considered here. The reduction in tower loads for the velocity feedback controllers compared to the
simple PI control system was a consequence of lower platform pitch response, which led to a reduction
in the wall thickness required to satisfy the long-term fatigue damage constraint. Consequently, the
tower costs were reduced, and also the platform costs due to better dynamic performance. Although
low-pass filtering of the nacelle velocity signal did not offer significant cost reductions, this control
strategy also saw a reduction in rotor speed variation, since the constraint was inactive at the optimum.
It is also expected that the effect of this filter will be more prominent for FWT concepts with larger
wave-frequency response.
Comparisons with nonlinear time-domain simulations showed that the linearized model in general is
able to capture trends with acceptable accuracy, but that the platform pitch response can be significantly
overpredicted for designs with low aerodynamic damping if contributions from other sources of damping
are small. This was the case for CS1, which indicates that the cost reductions achieved for CS2-
4 are considerably overestimated. For the velocity feedback controllers, which increases the amount
of aerodynamic damping induced by the control system, this problem was not observed, and good
agreement was achieved.
The presented approach is useful for conceptual FWT design, where it can be used to quickly explore
the design space before resorting to higher fidelity tools for detailed subsystem analysis and design.
The model captures important interactions between the controller and support structure, and enables
assessment of trade-off effects in a lifetime perspective. Further, the methodology can be extended to
account for additional design parameters and load cases, which may help identify novel design solutions.
References
[1] T. J. Larsen and T. D. Hanson. “A Method to Avoid Negative Damped Low Frequent Tower
Vibrations for a Floating, Pitch Controlled Wind Turbine”. In: Journal of Physics: Conference
Series 75 (2007), p. 012073.
[2] J. Jonkman. “Influence of Control on the Pitch Damping of a Floating Wind Turbine”. In: 46th
AIAA Aerospace Science Meeting and Exhibit, Reno, Nevada. 2008.
[3] M. Lackner. “Controlling Platform Motions and Reducing Blade Loads for Floating Wind
Turbines”. In: Wind Engineering 33.6 (2009), pp. 541–553.
[4] G. J. van der Veen, I. J. Couchman, and R. O. Bowyer. “Control of floating wind turbines”.
In: 2012 American Control Conference, Fairmont Queen Elizabeth, Montreal, Canada. 2012,
pp. 3148–3153.
[5] P. A. Fleming et al. “Evaluating Methods for Control of an Offshore Floating Turbine”. In:
Proceedings of the ASME 2014 33rd International Conference on Ocean, Offshore and Arctic
Engineering (OMAE2014), San Francisco, California, USA. 2014.
[6] P. A. Fleming, A. Peiffer, and D. Schlipf. “Wind Turbine Controller to Mitigate Structural Loads
on a Floating Wind Turbine Platform”. In: Journal of Offshore Mechanics and Arctic Engineering
141 (2019).
[7] F. Lemmer et al. “Optimization of Floating Offshore Wind Turbine Platforms With a Self-Tuning
Controller”. In: Proceedings of the ASME 2017 36th International Conference on Ocean, Offshore
and Arctic Engineering (OMAE2017), Trondheim, Norway. 2017.
[8] C. Bak et al. Description of the DTU 10 MW Reference Wind Turbine. Tech. rep. DTU Wind
Energy Report-I-0092. DTU Wind Energy, 2013.
[9] J. M. Hegseth, E. E. Bachynski, and J. R. R. A. Martins. “Integrated design optimization of spar
floating wind turbines”. In: Marine Structures 72 (2020).
10
[10] J. M. Hegseth and E. E. Bachynski. “A semi-analytical frequency domain model for efficient
design evaluation of spar floating wind turbines”. In: Marine Structures 64 (2019), pp. 186–210.
[11] IEC. Wind Turbines - Part 1: Design Requirements. Tech. rep. IEC 61400-1. 2005.
[12] T. Dirlik. “Application of Computers in Fatigue Analysis”. PhD thesis. University of Warwick,
1985.
[13] L. Kendall et al. “Application of Proportional-Integral and Disturbance Accommodating Control
to Variable Speed Variable Pitch Horizontal Axis Wind Turbines”. In: Wind Engineering 21.1
(1997), pp. 21–38.
[14] C. L. Bottasso et al. “Optimization-based study of bend-twist coupled rotor blades for passive and
integrated passive/active load alleviation”. In: Wind Energy 16 (2013), pp. 1149–1166.
[15] K. Johannessen, T. S. Meling, and S. Haver. “Joint Distribution for Wind and Waves in the
Northern North Sea”. In: International Journal of Offshore and Polar Engineering 12.1 (2002),
pp. 1–8.
[16] J. S. Gray et al. “OpenMDAO: An open-source framework for multidisciplinary design, analysis,
and optimization”. In: Structural and Multidisciplinary Optimization 59 (2019), pp. 1075–1104.
[17] P. E. Gill, W. Murray, and M. A. Saunders. “SNOPT: an SQP Algorithm for Large-Scale
Constrained Optimization”. In: SIAM Journal on Optimization 12.4 (2002), pp. 979–1006.
[18] R. E. Perez, P. W. Jansen, and J. R. R. A. Martins. “pyOpt: A Python-Based Object-Oriented
Framework for Nonlinear Constrained Optimization”. In: Structural and Multidisciplinary
Optimization 45.1 (2012), pp. 101–118.
[19] J. Farkas and K. Jármai. Optimum Design of Steel Structures. Springer, 2013.
[20] DNV. Fatigue Design of Offshore Steel Structures. Tech. rep. DNV-RP-C203. 2010.
[21] DNV. Design of Floating Wind Turbine Structures. Tech. rep. DNV-OS-J103. 2013.
[22] European Committee for Standardization. Eurocode 3: Design of Steel Structures, Part 1-6:
Strength and Stability of Shell Structures. Tech. rep. EN 1993-1-6: 2007. 2007.
[23] G. Kreisselmeier and R. Steinhauser. “Systematic Control Design by Optimizing a Vector
Performance Index”. In: International Federation of Active Controls Symposium on Computer-
Aided Design of Control Systems, Zurich, Switzerland. 1979.
[24] DNV. Modelling and Analysis of Marine Operations. Tech. rep. DNV-RP-H103. 2011.
[25] M. H. Hansen and L. C. Henriksen. Basic DTU Wind Energy controller. Tech. rep. DTU Wind
Energy Report-E-0018. 2013.
[26] J. R. R. A. Martins and A. B. Lambe. “Multidisciplinary Design Optimization: A Survey of
Architectures”. In: AIAA Journal 51.9 (2013), pp. 2049–2075.
[27] M. Muskulus and S. Schafhirt. “Design optimization of wind turbine support structures-A
review”. In: Journal of Ocean and Wind Energy 1.1 (2014), pp. 12–22.
[28] C. E. S. Souza, J. M. Hegseth, and E. E. Bachynski. “Frequency-dependent aerodynamic damping
and inertia in linearized dynamic analysis of floating wind turbines”. In: Journal of Physics:
Conference Series 1452.012040 (2020).
11
188 Appended Papers
Paper 5
189
220 Appended Papers
Appendix B
Formulation of Rotor
Effective Wind Speed
where SU (ω) is the incoming wind spectrum, ψ is the azimuth angle between
the elements, and Kn is the nth Fourier coefficient of the coherence function
γ, which is Fourier expanded in the rotor plane:
π
1
Kn(jk) (ω) = γ(ω, djk ) cos(nθ) dθ. (B.2)
π 0
The coherence function for two points with separation distance d and fre-
quency f in Hz is taken from IEC (2005) and expressed as
⎛ ! ⎞
"2 ! "2 0.5
fd 0.12d
γ(f, d) = exp ⎝−12 + ⎠, (B.3)
Uhub Lc
where Uhub is the mean hub-height wind speed and Lc is a coherence scale
parameter.
Using Eq. (B.1), the wind speed seen by each blade element can be found,
and used together with aerodynamic derivatives to calculate blade root
221
222 Formulation of Rotor Effective Wind Speed
However, if only the resultant loads are of interest, and the blades are
identical, the calculations can be simplified by defining a single ‘effective’
blade, and the cross spectral density matrix for the effective blade can be
found from the auto-spectral density for one of the blades in Eq. (B.4). For
the thrust force and aerodynamic torque, which is the sum of the blade root
shear force Fy and edgewise bending moment My respectively, the same
formulation can be used:
∞
Sv (ω) = 32 Sv,n (ω), (B.5)
n=−∞
where
⎧
⎪
⎨Kn (|ω − nϕ̇|) SU (|ω − nϕ̇|), n ∈ {..., −6, −3, 0, 3, 6, ...}
Sv,n (ω) =
⎪
⎩0, otherwise,
(B.6)
and Kn only considers elements on the same blade. For the tilting moment,
the effective blade formulation is somewhat different, as the azimuth angle
of the blade affects the contribution from the blade root moment on the
resultant load. The contribution from the flapwise blade root moment on
the rotor titling moment is proportional to cos(ϕ̇t), and the cross spectral
density matrix thus becomes:
! "2
∞
3
Sv (ω) = Sv,n (ω), (B.7)
2 n=−∞
where
⎧
⎪
⎪K (|ω − (n + 1)ϕ̇|) SU (|ω − (n + 1)ϕ̇|), n ∈ {..., −7, −4, −1, 2, 5, ...}
⎪ n
⎪
⎨
Sv,n (ω) = Kn (|ω − (n − 1)ϕ̇|) SU (|ω − (n − 1)ϕ̇|), n ∈ {..., −5, −2, 1, 4, 7, ...}
⎪
⎪
⎪
⎪
⎩
0, otherwise
(B.8)
223
As seen from the non-zero terms in Eqs. (B.6) and (B.8), the cancellation of
harmonics with a three-bladed turbine results in peaks at multiples of three
of the rotor frequency (i.e. 3P, 6P, 9P, etc.) for the resultant loads. The
tower top load spectra can be found by pre- and post-multiplying Sv (ω)
with the weight factors for each element, i.e. factors that relate wind speed
to the aerodynamic loads on the element:
Here,
⎡ ⎤ ⎡ ⎤ ⎡ ⎤
Fn,v (r1 )Δr r1 Fn,v (r1 )Δr r1 Ft,v (r1 )Δr
⎢ F (r )Δr ⎥ ⎢ r F (r )Δr ⎥ ⎢ r F (r )Δr ⎥
⎢ n,v 2 ⎥ ⎢ 2 n,v 2 ⎥ ⎢ 2 t,v 2 ⎥
WFT =⎢ ⎥ , WMT = ⎢ ⎥ , WQA = ⎢ ⎥.
⎣ · ⎦ ⎣ · ⎦ ⎣ · ⎦
Fn,v (rN )Δr rN Fn,v (rN )Δr rN Ft,v (rN )Δr
(B.10)
The spectra for the effective wind speeds can be calculated from the tower
top load spectra and the total weight factors over the rotor:
SFT (ω)
SvFT (ω) = (B.11a)
(3 Fy,v )2
SMT (ω)
SvMT (ω) = 2 (B.11b)
3
2 Mz,v
SQA (ω)
SvQA (ω) = . (B.11c)
(3 My,v )2
The cross spectral density matrix for the wind loads can be found by ob-
serving which harmonics of the wind are shared by the different load com-
ponents. As shown by Halfpenny (1998), the thrust force and aerodynamic
torque on the turbine are perfectly correlated, as they share all harmon-
ics, and completely uncorrelated from the tilting moment, as they have no
harmonics in common.
224 Formulation of Rotor Effective Wind Speed
Previous PhD theses published at the Department of Marine Technology
(earlier: Faculty of Marine Technology)
NORWEGIAN UNIVERSITY OF SCIENCE AND TECHNOLOGY
MTA-91- Løland, Geir, MH Current forces on and flow through fish farms.
78 (Dr.Ing. Thesis)
MTA-92- Mørch, Hans Jørgen Bjelke, MH Aspects of Hydrofoil Design: with Emphasis on
86 Hydrofoil Interaction in Calm Water. (Dr.Ing.
Thesis)
MTA-92- Ask, Tor Ø., MM Ignition and Flame Growth in Lean Gas-Air
91 Mixtures. An Experimental Study with a Schlieren
System. (Dr.Ing. Thesis)
MTA-96- Erikstad, Stein O., MP A Decision Support Model for Preliminary Ship
114 Design. (Dr.Ing. Thesis)
MTA-99- Andersen, Trond M., MM Short Term Maintenance Planning. (Dr.Ing. Thesis)
130
MTA-99- Tveiten, Bård Wathne, MK Fatigue Assessment of Welded Aluminium Ship
131 Details. (Dr.Ing. Thesis)
MTA- Heggelund, Svein E., MK Calculation of Global Design Loads and Load
2001-147 Effects in Large High Speed Catamarans. (Dr.Ing.
Thesis)
MTA- Mohammed, Abuu K., MK Nonlinear Shell Finite Elements for Ultimate
2001-149 Strength and Collapse Analysis of Ship Structures.
(Dr.Ing. Thesis)
MTA- Holmedal, Lars E., MH Wave-current interactions in the vicinity of the sea
2002-150 bed. (Dr.Ing. Thesis)
MTA- Rognebakke, Olav F., MH Sloshing in rectangular tanks and interaction with
2002-151 ship motions. (Dr.Ing. Thesis)
MTA- Yang, Qinzheng, MH Wash and wave resistance of ships in finite water
2002-153 depth. (Dr.Ing. Thesis)
IMT- Reite, Karl Johan Modelling and Control of Trawl Systems. (Dr.Ing.
2006-19 Thesis)
IMT- Rustad, Anne Marthine Modelling and Control of Top Tensioned Risers.
2007-23 (PhD Thesis, CeSOS)
IMT- Refsnes, Jon Erling Gorset Nonlinear Model-Based Control of Slender Body
2008-29 AUVs (PhD Thesis, IMT)
IMT- Nystad, Bent Helge Technical Condition Indexes and Remaining Useful
2008-43 Life of Aggregated Systems. PhD thesis, IMT
IMT- Pedersen, Tom Arne Bond Graph Modelling of Marine Power Systems.
2009-46 PhD Thesis, IMT
IMT- Korsvik, Jarl Eirik Heuristic Methods for Ship Routing and
2009-51 Scheduling. PhD-thesis, IMT
IMT- Bruun, Kristine Bond Graph Modelling of Fuel Cells for Marine
2009-54 Power Plants. Ph.d.-thesis, IMT
IMT Seim, Knut Sponheim Mixing Process in Dense Overflows with Emphasis
2010-65 on the Faroe Bank Channel Overflow. Ph.d.thesis,
IMT
IMT Abrahamsen, Bjørn Christian Sloshing Induced Tank Roof with Entrapped Air
2010-68 Pocket. Ph.d.thesis, CeSOS.
IMT- You, Jikun Numerical studies on wave forces and moored ship
2012-86 motions in intermediate and shallow water, CeSOS.
IMT-10- Ramìrez, Pedro Agustìn Pèrez Ageing management and life extension of technical
2013 systems-
Concepts and methods applied to oil and gas
facilities, IMT
IMT-13- Lindstad, Haakon Strategies and measures for reducing maritime CO2
2013 emissons, IMT
IMT-4- Bachynski, Erin Elizabeth Design and Dynamic Analysis of Tension Leg
2014 Platform Wind Turbines, CeSOS
IMT-8- Muliawan, Made Jaya Design and Analysis of Combined Floating Wave
2014 and Wind Power Facilities, with Emphasis on
Extreme Load Effects of the Mooring System,
CeSOS
IMT-12- Kvittem, Marit I. Modelling and response analysis for fatigue design
2014 of a semisubmersible wind turbine, CeSOS
IMT-13- Akhtar, Juned The Effects of Human Fatigue on Risk at Sea, IMT
2014
IMT-14- Syahroni, Nur Fatigue Assessment of Welded Joints Taking into
2014 Account Effects of Residual Stress, IMT
IMT-4- Jose Patricio Gallardo Canabes Numerical studies of viscous flow around bluff
2015 bodies, IMT
IMT-5- Vegard Longva Formulation and application of finite element
2015 techniques for slender marine structures subjected
to contact interactions, IMT
IMT-9- Daniel de Almeida Fernandes An output feedback motion control system for
2015 ROVs, AMOS
IMT-13- Arturo Jesùs Ortega Malca Dynamic Response of Flexibles Risers due to
2015 Unsteady Slug Flow, CeSOS
IMT-3- Mia Abrahamsen Prsic Numerical Simulations of the Flow around single
2016 and Tandem Circular Cylinders Close to a Plane
Wall, IMT
IMT-7- Øivind Kåre Kjerstad Dynamic Positioning of Marine Vessels in Ice, IMT
2016
IMT-12- Øyvind Selnes Patricksson Decision support for conceptual ship design with
2016 focus on a changing life cycle and future
uncertainty, IMT
IMT-16- Wilson Ivan Guachamin Acero Assessment of marine operations for offshore wind
2016 turbine installation with emphasis on response-
based operational limits, IMT
IMT-1- Yingguang Chu Virtual Prototyping for Marine Crane Design and
2018 Operations
IMT-10- Chenyu Luan Design and analysis for a steel braceless semi-
2018 submersible hull for supporting a 5-MW horizontal
axis wind turbine
IMT-13- Stein Melvær Nornes Guidance and Control of Marine Robotics for
2018 Ocean Mapping and Monitoring
IMT-4- Siri Mariane Holen Safety in Norwegian Fish Farming – Concepts and
2019 Methods for Improvement
IMT-12- Jørgen Bremnes Nielsen Modeling and Simulation for Design Evaluation
2019
IMT-17- Endre Sandvik Sea passage scenario simulation for ship system
2019 performance evaluation
IMT-18- Loup Suja-Thauvin Response of Monopile Wind Turbines to Higher
2019 Order Wave Loads