NRELOffshrBsline5MW PDF
NRELOffshrBsline5MW PDF
NRELOffshrBsline5MW PDF
doc
Jason Jonkman, NREL/NWTC
June 14, 2006
Introduction:
This manuscript documents NRELs baseline wind turbine aeroelastic model for use in various
offshore analysis concept studies. This includes the aeroelastic FAST5 model, baseline control
system (torque and pitch), and rationale behind the model. The turbine parameters are fictitious
and do not represent an actual wind turbine, but do represent a good approximation of what an
actual wind turbine of the associated size would be like.
The FAST model can be used to generate an ADAMS7 model using the FAST-to-ADAMS
preprocessor built-into FAST. NREL has compared the responses between the FAST and
ADAMS model results for a number of simulationsthe results agree pretty well in general
the most notable differences are caused by the availability of mass offsets and blade torsion
DOFs in ADAMS, which are not available in FAST.
Numerous additional details and rationale may be found in the MS Excel workbook,
NRELOffshrBsline5MW.xls, which is available from:
Jason Jonkman
National Wind Technology Center (NWTC)
National Renewable Energy Laboratory (NREL)
1617 Cole Boulevard
Golden, CO 80401-3393
United States of America
Phone: +1 (303) 384 - 7026
Fax: +1 (303) 384 - 6901
E-mail: [email protected]
Next, the aeroelastic properties of the 5MW baseline had to be established. In order for the
concept studies to be defendable, the aeroelastic model had to be realistic. In order to achieve
realism, NREL gathered publicly available gross properties of several aeroelastic models and
1
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
created a composite from these systems, extracting the best available information. A summary
of the aeroelastic models utilized is given below:
FAST model of Global Energy Concepts (GECs) Wind Partnerships for Advanced
Component Development (WindPACT)14 5A0210.
FAST model of GECs WindPACT 5A0410.
BLADED model used during Garrad Hassans RECOFF project19.
BLADED and PHATAS models used during TU Delfts Dutch Offshore Wind Energy
Converter (DOWEC) project2,6,8.
and of the two largest wind turbine prototypes in the world:
Multibrid M50001,11,12
REpower 5M9,16,17
These properties are provided in the following table:
Note that the rotor diameters supplied in the table above ignore the effects of blade precone,
which reduce the actual rotor diameter and swept area.
The Multibrid M5000 machine has a significantly higher tip speed than typical onshore wind
turbines and lower tower top mass than would be expected from recently developed scaling
laws18. In contrast, the REpower 5M machine has properties that are more expected. In this
sense, the REpower 5M machine is more conventional than the Multibrid M5000, so it makes
2
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
sense to use the specifications of the REpower 5M machine as the target specifications for the
model. Thus the task was to establish an aeroelastic model that tries to resemble the REpower
5M machine as close as possible (i.e., as close as publicly known) and to fill in the data gaps
using machine properties known from the DOWEC, RECOFF, and WindPACT studies. From
the table above, the DOWEC turbine most closely matched the REpower 5M machine, which
makes sense since the DOWEC study used aeroelastic blade properties provided by LM
glasfiber9, and the REpower 5M machine also utilizes LM glasfiber blades. As such, must of the
data gaps are filled in using machine properties from the DOWEC study6.
Knowing that the rotor radius for our aeroelastic model would be about 63m (from the REpower
5M machine), and wanting the lowest reasonable hub-height possible in order to minimize the
overturning moment of the platform, NREL decided that the hub height for the baseline model
should be 90m, which would give a 15m air gap between the blade tips at their lowest point and
a 30m (15m amplitude) 50-yr extreme wave height.
Additional gross properties chosen for the NREL 5MW baseline wind turbine, most of which are
identical to the REpower 5M, are provided in the table below:
Gross Properties Chosen for the NREL 5MW Baseline Wind Turbine Model
Rating 5MW
Wind Regime IEC 61400-3 (Offshore) Class 1B / Class 6 winds
Rotor Orientation Upwind
Control Variable Speed, Collective Pitch
Rotor Diameter / Hub Diameter 126m / 3m
Hub Height 90m
Maximum Rotor / Generator Speed 12.1rpm / 1,173.7rpm
Maximum Tip Speed 80m/s
Overhang / Shaft Tilt / Precone 5m / 5 / 2.5
Rotor Mass 110,000kg Overall c.g. location:
Nacelle Mass 240,000kg (x ,y ,z ) = (-0.2m,0.0m,64.0m)
t t t
Tower Mass (Deepwater) 347,460kg
Reference Site National Data Buoy Center (NDBC) Buoy 44008
The tower mass listed in the table above is for the deepwater support structure configuration.
The overall c.g. location of this wind turbine as indicated above is in the tower base coordinate
system whose origin lies on the tower centerline at the mean sea level (MSL).
The aeroelastic FAST model, which cant support prebend, incorporates a 2.5 upwind precone.
This is used to represent the small amount of precone and larger amount of precurve built-into
the actual REpower 5M, 61.5m long, LM glasfiber blade9. The rotor and hub diameters
indicated in the table above ignore the effects of blade precone, which reduces the actual
diameters and swept area. The actual rotor diameter in the model (assuming that the blades are
undeflected) is (126m)*COS(2.5) = 125.88m and the actual swept area is (/4)*(125.88m)^2 =
12,445.3m^2.
3
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
Blade Structural Properties:
The distributed blade structural properties of the NREL 5MW baseline model are based on the
structural properties of the 62.6m long LM Glasfiber blade used in the DOWEC study, as given
in Appendix A of Ref. 8. Since the blades in the DOWEC study are 1.1m longer than the
REpower 5M blades, the structural properties of the NREL 5MW baseline blade are found by
truncating the 62.6m blades at 61.5m blade span (the structural properties of the blade tip are
found by interpolating between the 61.2m and 61.7m stations provided in Appendix A of Ref. 8).
The resulting properties are given in the table below:
The entries in the first column, Radius, are the spanwise locations along the blade pitch axis
relative to the rotor center (apex).
BlFract is the fractional distance along the blade pitch axis from the root (0.0) to the tip (1.0).
The blade root is located 1.5m along the pitch axis from the rotor center (apex).
AeroCent is a FAST variable name. By definition, the quantity (AeroCent - 0.25) is the
fractional distance to the aerodynamic center from the pitch axis along the chordline, positive
towards the trailing edge. Thus, at the root (BlFract = 0.0), AeroCent = 0.25 means that the
aerodynamic center lies on the pitch axis [since (0.25 - 0.25) = 0.0]. And at the tip (BlFract =
1.0), AeroCent = 0.125 means that the aerodynamic center lies 0.125 chordlengths toward the
leading edge from pitch axis [since (0.125 - 0.25) = -0.125].
4
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
The flapwise and edgwise section stiffness and inertia values, FlpStff, EdgStff, FlpIner and
EdgIner, are given about the principal structural axes of each cross section as oriented by the
structural twist angle, StrcTwst. The values of StrcTwst as given in the table above, are offset by
-0.09182 from the values provided in Appendix A of Ref. 8 in order to insure that the zero-twist
reference location is at the blade tip.
GJStff are the values of the blade torsional stiffness. Since the DOWEC blade data did not
contain extensional stiffness information, the blade extensional stiffness values, EAStff, given in
the table above were estimated to be 10^7 times the average mass moment of inertia at each
blade station (this rule of thumb comes from WindPACT data). This rule of thumb should get us
in the ballpark at leastthe exact values are not important due to the low rotational speed.
The flapwise c.g. offset values, FlpcgOf, and flapwise and edgewise elastic offset values,
FlpEAOf and EdgEAOf, given in the table above were assumed to be zero due to the inability to
decipher the corresponding information given in Appendix A of Ref. 8. The edgewise c.g. offset
values, EdgcgOf, are the distances in meters along the chordline from the pitch axis to the center
of gravity of the blade section, positive towards the trailing edge.
The distributed blade section mass per unit length values, BMassDen, given in the table above
are increased by 4.536% in the model in order to scale the overall (integrated) blade mass to
17,740kg, which is the mass of the blades in the REpower 5M machine. In the model, the
second mass moment of inertia, first mass moment of inertia, and the radial c.g. location of each
blade is 11,776,047 kg-m^2, 363,231 kg-m, and 20.475m with respect to the blade root,
respectively.
The NREL 5MW baseline model incorporates 0.477465% critical structural damping in all
modes of the isolated blade, which corresponds to the 3% logarithmic decrement used in the
DOWEC study from page 20 of Ref. 6.
The table below summarizes the undistributed blade structural properties discussed in this
section:
5
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
aerodynamic gradients at the blade tip. The aerodynamic properties at the blade nodes, which
are located at the center of the blade elements, are given in the table below. The blade node
locations, labeled as RNodes in the table, are directed along the pitch axis from the rotor center
(apex) to the blade cross sections.
Integrating the chord distribution along the blade span reveals that the rotor solidity is roughly
5.16%.
As indicated in the table above, the NREL 5MW baseline model incorporates eight unique airfoil
data tables. The two innermost airfoil tables represent cylinders with drag coefficients of 0.50
(Cylinder1.dat) and 0.35 (Cylinder2.dat) and no lift. The remaining six airfoil tables were
created by incorporating 3D corrections to the 2D airfoil data coefficients of the six airfoils used
in the DOWEC study as detailed in Appendix A of Ref. 6 (numerical values for these
coefficients were provided by Koert Lindenburg of ECN). The lift (Cl) and drag (Cd) values
were corrected for rotational stall delay for 0 to 90deg angle-of-attack, the Cd values were
additionally corrected using the Viterna method for 0 to 90deg angle-of-attack, and the Beddoes-
Leishman stall hysteresis parameters were found using version 2.0 of AirfoilPrep3. All
corrections were made assuming an aspect ratio of 17. The 3D corrected airfoil data coefficients
are illustrated graphically in the figures below:
6
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
3D Corrected Coefficients of the DU40 Airfoil
CL
2.5
CD
CM (not req'd.)
2.0
1.5
1.0
0.5
0.0
-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180
-0.5
-1.0
-1.5
1.5
1.0
0.5
0.0
-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180
-0.5
-1.0
-1.5
7
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
3D Corrected Coefficients of the DU30 Airfoil
CL
2.5
CD
CM (not req'd.)
2.0
1.5
1.0
0.5
0.0
-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180
-0.5
-1.0
-1.5
1.5
1.0
0.5
0.0
-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180
-0.5
-1.0
-1.5
8
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
3D Corrected Coefficients of the DU21 Airfoil
CL
2.5
CD
CM (not req'd.)
2.0
1.5
1.0
0.5
0.0
-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180
-0.5
-1.0
-1.5
1.5
1.0
0.5
0.0
-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180
-0.5
-1.0
-1.5
9
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
from Table 6 on page 26 of Ref. 6. The elevation of the yaw bearing is consequently 87.6m.
With a shaft tilt of 5, this makes the distance directed along the shaft from the hub center to the
yaw axis 5.01910m and the vertical distance along the yaw axis from the tower-top to the shaft
1.96256m. The distance directed along the shaft from the hub center to the main bearing is taken
to be 1.912m, from Table 6 on page 26 of Ref. 6.
The hub mass is taken to be 56,780 kg like the REpower 5M and is located at the hub center.
The hub inertia about the shaft is taken to be 115,926 kgm^2. This was found by assuming that
the hub spindel is a thin spherical shell with a radius 0.25m longer than the hub radius (due to the
fact that the pitch bearing is located within the spindel0.25m is chosen since the nacelle height
of the DOWEC 6MW is 3.5m based on Table 6 on page 26 of Ref. 6), which have theoretical
mass moment of inertias of: Ixx = 2/3*m*r^2.
The nacelle mass is taken to be 240,000 kg like the REpower 5M and is located 1.9m downwind
of the yaw axis like the DOWEC 6MW (from Table 7 on page 27 of Ref. 6) and 1.75m above the
yaw bearing, which is half the height of the DOWEC 6MW nacelle (from Table 6 on page 26 of
Ref. 6). The nacelle inertia about the yaw axis is taken to be 2,607.89E3 kgm^2. This was
chosen to be the same as the DOWEC 6MW's nacelle inertia about its nacelle c.g. and translated
to the yaw axis using the parallel axis theorem with the nacelle mass and downwind distance to
the nacelle c.g..
The nacelle yaw actuator is taken to have a natural frequency of 3Hz, which is equivalent to the
highest full structural system natural frequency in the FAST model, and a damping ratio of 2%
critical. This resulted in an equivalent nacelle yaw linear spring constant of 9,028.32E6 N-m/rad
and an equivalent nacelle yaw linear damping constant of 19.16E6 N-m/rad/s.
The nacelle and hub properties discussed in this section are summarized in the table below:
Drivetrain Properties:
The NREL 5MW baseline wind turbine model uses the same rated rotor speed (12.1 rpm), rated
generator speed (1173.7 rpm), and gearbox ratio (97:1) as the REpower 5M machine. The
gearbox is assumed to have no frictional lossesa requirement of the FAST-to-ADAMS
10
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
preprocessor. The electrical efficiency of the generator is taken to be 94.4%. This is chosen to
be roughly the same as the total mechanical to electrical conversion loss used by the DOWEC
6MW at ratedi.e., 3.5E5 Watts of power loss occur at about 6.25E6 Watts of aerodynamic
power from Figure 15, page 24 of Ref. 6. The generator inertia about the high-speed shaft is
taken to be 534.116 kgm^2, which is the same equivalent low-speed shaft generator inertia used
in the DOWEC studyi.e., 5025500 kgm^2 from page 36 of Ref. 6.
The drive shaft is taken to have the same natural frequency as the RECOFF model and 5%
critical damping ratio for the free-free mode with a rigid generator and rigid rotor. This resulted
in an equivalent drive shaft linear spring constant of 867.637E6 N-m/rad and a damping constant
of 6.215E6 N-m/rad/s.
The high-speed shaft brake is taken to have the same maximum brake torque to maximum
generator torque ratio and the same time lag as used in the DOWEC study (from page 29 of Ref.
6). This resulted in a fully-deployed high-speed shaft brake torque of 28,116.2 N-m and a time
lag of 0.6 seconds. This time lag is the amount of time it takes for the brake to fully engage once
deployed. FAST and ADAMS use a simple linear ramp from 0 to full over the 0.6 seconds.
The drivetrain properties discussed in this section are summarized in the table below:
Drivetrain Properties
Rated Rotor Speed 12.1 rpm
Rated Generator Speed 1173.7 rpm
Gearbox Ratio 97 :1
Electrical Generator Efficiency 94.4 %
Generator Inertia About High-Speed Shaft 534.116 kg-m^2
Equivalent Drive Shaft Torsional Spring Constant 867,637,000 N-m/rad
Equivalent Drive Shaft Torsional Damping Constant 6,215,000 N-m/rad/s
Fully-Deployed High-Speed Shaft Brake Torque 28,116.2 N-m
High-Speed Shaft Brake Time Constant 0.6 sec
Shallow Water:
The distributed tower properties of the NREL 5MW baseline model are based on the base
diameter (6m) and thickness (0.027m), top diameter (3.87m) and thickness (0.019m), and
effective mechanical steel properties of the tower used in the DOWEC study, as given in Table 9
on page 31 of Ref. 6. The radius and thickness of the tower are assumed to be linearly tapered
from base to top. The tower is connected to a monopile with a constant diameter of 6m and a
constant thickness of 0.060m. The tower base begins at an elevation of 10m above the MSL.
The monopile extends from here down to the mudline, which is at 20m below MSL. The
Youngs modulus is taken to be 210*10^9 Pa, the shear modulus is taken to be 80.8*10^9 Pa and
the effective density of the steel is taken to be 8500 kg/m^3. The 8500 kg/m^3 is meant to be an
11
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
increase above steels typical value of 7850 kg/m^3 in order to account for paint, bolts, welds,
and flanges that are not accounted for in the thickness data. The resulting distributed tower
properties are given in the table below:
The entries in the first column, Elevation, are the vertical locations along the tower centerline
relative to the MSL. HtFract is the fractional height along the tower centerline from the mudline
(0.0) to the tower top (1.0). The rest of columns are similar to those described for the distributed
blade properties.
The resulting overall (integrated) tower + monopile mass is 522,617 kg and is centered at
37.172m along the tower centerline above the mudline. This result follows directly from the
overall tower height of 87.6m.
The NREL 5MW baseline model incorporates 1% critical structural damping in all modes of the
isolated tower (without the top mass present), which corresponds to the values used in the
DOWEC study from page 21 of Ref. 6.
The table below summarizes the undistributed tower properties discussed in this section:
Deep Water:
The distributed tower properties of the NREL 5MW baseline model are based on the base
diameter (6m) and thickness (0.027m), top diameter (3.87m) and thickness (0.019m), and
effective mechanical steel properties of the tower used in the DOWEC study, as given in Table 9
on page 31 of Ref. 6. The Youngs modulus is taken to be 210*10^9 Pa, the shear modulus is
taken to be 80.8*10^9 Pa and the effective density of the steel is taken to be 8500 kg/m^3. The
8500 kg/m^3 is meant to be an increase above steels typical value of 7850 kg/m^3 in order to
account for paint, bolts, welds, and flanges that are not accounted for in the tower thickness data.
12
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
The radius and thickness of the tower are assumed to be linearly tapered from base to top. The
thickness of the tower was scaled up relative to the values indicated above in order to strengthen
the tower since the REpower 5M machine has a larger tower-top mass than the DOWEC 6MW
machine. An increase of 30% was chosen to ensure that the first tower fore-aft (FA) and side-to-
side (SS) tower frequencies lie between 1P and 3P through the operational range of the wind
turbine in a Campbell diagram. The resulting distributed tower properties are given in the table
below:
The entries in the first column, Elevation, are the vertical locations along the tower centerline
relative to the ground or MSL. HtFract is the fractional height along the tower centerline from
the tower base (0.0) to the tower top (1.0). The rest of columns are similar to those described for
the distributed blade properties.
The resulting overall (integrated) tower mass is 347,460 kg and is centered at 38.234m along the
tower centerline above the MSL. This result follows directly from the overall tower height of
87.6m.
The NREL 5MW baseline model incorporates 1% critical structural damping in all modes of the
isolated tower (without the top mass present), which corresponds to the values used in the
DOWEC study from page 21 of Ref. 6.
The table below summarizes the undistributed tower properties discussed in this section:
13
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
wind turbine is required for such analyses. As defined in equation 141 on page 149 of Ref. 15,
the body mass matrix of the rigid, undeflected 5MW baseline wind turbine in deep water to four-
significant digits is:
697,500 0 0 0 44,630,000 0
0 697,500 0 44,630,000 0 144,800
0 0 697,500 0 144,800 0
Mij =
0 44,630,000 0 3,566,000,000 0 11,670,000
44,630,000 0 144,800 0 3,551,000,000 0
0 144,800 0 11,670,000 0 25,410,000
This matrix was derived by linearizing the FAST model about a reference point located at the
tower centerline at the base of the tower (i.e., MSL). The row/column order is: 1 = surge, 2 =
sway, 3 = heave, 4 = roll, 5 = pitch, 6 = yaw and the units in the matrix are: kg, kg-m, and kg-
m^2.
Naturally, since FAST does not have a tower torsion DOF, neither does it have knowledge of the
yaw inertia of the tower about the tower centerline. So, the effects of the tower yaw inertia are
not included in the body mass matrix given above. By integration of the distributed tower inertia
data along the length of the tower, this yaw inertia is estimated to be 2,241,820 kg-m^2. This
inertia could be added to element (6,6) of the body mass matrix given above, but it is not
straightforward how this additional inertia will affect the other elements in the matrix. However,
neglecting the tower yaw inertia from the body mass matrix should not be of significant concern,
since it may be seen that the addition of the tower yaw inertia to element (6,6) will only increase
the overall turbine yaw inertia (tower, nacelle, and rotor) by a small amount (less than 10%).
Users of the body mass matrix should be aware of this issue and address it as they see fit.
14
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
Control System Flowchart
The generator speed measurement for both the torque and pitch controllers is filtered using a
recursive, single-pole, low-pass filter with exponential smoothing. The corner frequency (-3dB
point) of the low pass filter is set to 0.25 Hz, which is chosen to be roughly of the blade
edgewise natural frequency.
The generator torque is computed as a function of the filtered generator speed using a variable-
speed controller. The variable-speed controller incorporates five control regions: 1, 1, 2, 2,
and 3. Region 1 is the control region before cut-in where the generator torque is zero so that no
power is extracted from the wind; instead, the wind is used to accelerate the rotor. Region 2 is
the control region for optimizing power capture. Here, the generator torque is proportional to the
square of the filtered generator speed. In region 3, the generator power is held constant so that
the generator torque is inversely proportional to the filtered generator speed. Region 1 is the
start-up region and is a linear transition between regions 1 and 2. It is used to place a lower limit
on the generator speed in order to limit the operational speed range. Region 2 is a linear
transition between regions 2 and 3 with a torque slope corresponding to the slope of an induction
machine. Region 2 is typically needed (as is the case for this 5MW model) since a machine
does not typically reach rated torque at its rated speed using Region 2s control law [i.e., the
optimal K is typically lower than that which would make the rated torque = K*(rated speed)^2,
since the rated speed is generally limited from optimal in order to limit tip speed for noise
reasons].
15
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
The peak of the power coefficient as a function of tip-speed ratio and blade pitch surface was
found by running FAST simulations at a number of fixed rotor speeds and a number of fixed
rotor-collective blade pitch angles at a fixed wind speed of 8 m/s. From these simulations, the
peak power coefficient of 0.482 was found to occur at a tip-speed ratio of 7.55 and a rotor-
collective blade pitch of 0.0. With the 97:1 gearbox ratio, this resulted in an optimal constant of
proportionality of 0.0255764 N-m/rpm^2 in the region 2 control law. With the rated generator
speed of 1173.7 rpm, rated electric power of 5MW, and a generator efficiency of 94.4%, the
rated mechanical power is 5.296610MW and the rated generator torque is 43,093.55 N-m.
Region 1 spans the range of generator speeds between 670 rpm and 30% above this value, or
871 rpm. The minimum generator speed of 670 rpm corresponds to the minimum rotor speed of
6.9 rpm used by the REpower 5M machine16. The transitional generator speed between Regions
2 and 3 is taken to be 99% of the rated generator speed, or 1,161.963 rpm. The generator slip
percentage in region 2 is taken to be 10% in accordance with the value used in the DOWEC
studyreference page 24 of Ref. 6. The resulting generator torque versus generator speed curve
is shown in the figure below:
45000
Optimal
40000
Variable Speed Controller
35000
Generator Torque, Nm
30000
25000
20000
15000
10000
5000
0
0 200 400 600 800 1000 1200 1400
Generator Speed, rpm
Finally, a conditional statement is placed on the variable speed controller so that the torque is
computed as if it was in Region 3 regardless of the generator speed whenever the previous pitch
angle command is greater or equal to 1. This results in improved output power quality (fewer
dips below rated) at the expense of short term overloading of the generator and gearbox. To
avoid excessive overloading of both, the torque is saturated to a maximum of 10% above rated,
or 47,402.91 N-m. Also, a torque rate limit of 15,000 N-m/s is imposed.
16
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
The collective blade pitch angle commands are computed using gain-scheduled PI control on the
speed error between the filtered generator speed and the rated generator speed (1173.7 rpm).
The parameter for determining the gain scheduling correction factor is the previous pitch angle
command. The gains (see the figure below) were derived according to the technique
documented in Ref. 4. In particular, the partial derivative of the rotor power with respect to the
collective pitch angle as a function of the pitch angle was computed by performing a
linearization analysis in FAST at a number of fixed pitch angles. The frozen wake assumption
was used when computing these partial derivatives through a minor modification to AeroDyn.
1.0 0.025
0.7
Correction Factor
0.6 0.015
Gains
0.5
0.4 0.010
0.3
0.2 0.005
0.1
0.0 0.000
0 5 10 15 20
Collective Blade Pitch Angle (deg)
The pitch rate limit is set at 8/sec in absolute value. This is speculated to be the pitch rate limit
of conventional 5MW machines based on GE Winds long blade test program. The minimum
and maximum pitch settings are chosen to be 0 and 90, respectively. The integral term in the
PI controller is saturated between these limits so as to ensure fast response between the Region 2
/ Region 3 transition.
Due to limitations in the code, the FAST model does not include any pitch actuator dynamic
effects. However, pitch actuator dynamics are required in ADAMS. Hence, in the ADAMS
model, the blade pitch actuator is taken to have a very high natural frequency of 30Hz, which is
equivalent to 10x the highest full structural system natural frequency in the FAST model, and a
damping ratio of 2% critical. This resulted in an equivalent blade pitch linear spring constant of
1,069.02E6 N-m/rad and an equivalent blade pitch linear damping constant of 0.23E6 N-m/rad/s.
17
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
The variable-speed generator-torque and blade pitch control properties discussed in this section
are summarized in the table below:
Time integration occurs at a constant time step of 0.0125 seconds in FASTs fixed-step size
integration scheme and at a maximum step size of 0.0125 seconds in ADAMS variable-step size
18
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
integration scheme. Aerodynamic calculations are performed every other time step (0.025
seconds) to ensure that there are at least 200 azimuth step computations per revolution. These
time steps are the largest possible to ensure numerical stability across a range of operating
modes.
19
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
Wind Speed Relationships Legend
BldPitch1 Pitch Angle of Blade 1
GenPwr Electrical Generator Power
GenSpeed Angular Speed of the High-Speed Shaft and Generator
GenTq Electrical Generator Torque
IPDefl1 In-Plane Tip Deflection of Blade 1
OoPDefl1 Out-of-Plane Tip Deflection of Blade 1
RotPwr Mechanical Rotor Power
RotSpeed Angular Speed of the Low-Speed Shaft and Rotor
RotThrust Rotor Thrust
RotTorq Mechanical Rotor Torque
TSR Rotor Blade Tip Speed Ratio
TTDspFA Fore-Aft Displacement of the Tower Top and Yaw Bearing
TTDspSS Side-to-Side Displacement of the Tower Top and Yaw Bearing
6000
Region 1 2 2 3
GenSpeed (rpm)
RotPwr (kW)
5000
GenPwr (kW)
RotThrust (kN)
RotTorq (kNm)
4000
3000
2000
1000
0
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Wind Speed (m/s)
20
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
50
Region 1 2 2 3
45 RotSpeed (rpm)
BldPitch1 (deg)
40 GenTq (kNm)
TSR (-)
35
30
25
20
15
10
0
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Wind Speed (m/s)
6
Region 1 2 2 3
OoPDefl1 (m)
5
IPDefl1 (m)
TTDspFA (m)
4 TTDspSS (m)
-1
-2
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Wind Speed (m/s)
References:
1
De Vries, Eize, Multibrid: 'A New Offshore Wind Turbine Contender',
http://www.jxj.com/magsandj/rew/2004_05/multibird.html, Renewable Energy World, Vol. 7,
No. 5, James & James (Science Publishers) Ltd, London, United Kingdom, September-October
2004.
2
Goezinne, F., Terms of reference DOWEC, 176-FG-R0300, DOWEC 10041_000, September
2001.
21
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
3
Hansen, Craig, NWTC Design Codes: AirfoilPrep,
http://wind.nrel.gov/designcodes/preprocessors/airfoilprep/, Golden, CO: National Renewable
Energy Laboratory, October 2004.
4
Hansen, Morten H.; Hansen, Anca; Larsen, Torben J.; ye, Stig; Srensen; and Fuglsang, Peter,
Control Design for a Pitch-Regulated, Variable-Speed Wind Turbine, Ris-R-1500(EN),
Roskilde, Denmark: Ris National Laboratory, January, 2005.
5
Jonkman, J.M. and Buhl, M.L., Jr., FAST User's Guide, NREL/EL-500-29798, Golden, CO:
National Renewable Energy Laboratory, October 2004.
6
Kooijman, H.J.T., Lindenburg, C., Winkelaar, D., and van der Hooft, E.L., DOWEC 6 MW
Pre-Design: Aero-elastic modeling of the DOWEC 6 MW pre-design in PHATAS, ECN-CX--
01-135, DOWEC 10046_009, Petten, the Netherlands: Energy Research Center of the
Netherlands, September 2003.
7
Laino, David J. and Hansen, A. Craig, User's Guide to the Computer Software Routines
AeroDyn Interface for ADAMS, Salt Lake City, UT: Windward Engineering LLC, Prepared
for the National Renewable Energy Laboratory under Subcontract No. TCX-9-29209-01,
September 2001.
8
Lindenburg, C., Aeroelastic Modelling of the LMH64-5 Blade, DOWEC-02-KL-083/0,
DOWEC 10083_001, Petten, the Netherlands: Energy Research Center of the Netherlands,
December 2002.
9
LM Glasfiber Group, Wind Turbine Blades, Product Overview, Standard Products Max.
Rated Power <=5000 kW,
http://www.lmglasfiber.dk/UK/Products/Wings/ProductOverView/50000kw.htm, Lunderskov,
Denmark: LM Glasfiber Group, January 2005.
10
Malcolm, D.J. and Hansen, A.C., WindPACT Turbine Rotor Design Study, NREL/SR-500-
32495, Golden, CO: National Renewable Energy Laboratory, August 2002.
11
Multibrid Technology, Technical Data Multibrid M5000,
http://www.multibrid.com/download/Datenblatt_M5000_eng.pdf, Bremerhaven, Germany:
Multibrid Technology, March 2004.
12
Multibrid Technology, The Concept in Detail, http://www.multibrid.com/english/concept.htm,
Bremerhaven, Germany: Multibrid Technology, January 2005.
13
Musial, Walt, Butterfield, Sandy, and Boone, Andrew, Feasibility for Floating Platform
Systems for Wind Turbines, Collection of the 2004 ASME Wind Energy Symposium Technical
Papers Presented at the 42nd AIAA Aerospace Sciences Meeting and Exhibit, 58 January 2004,
Reno, NV, Washington, D.C.: American Institute of Aeronautics and Astronautics, January 2004,
pp. 476486; NREL/CP-500-36504.
22
NRELOffshrBsline5MW.doc
Jason Jonkman, NREL/NWTC
June 14, 2006
14
National Renewable Energy Laboratory, About the Program: WindPACT,
http://www.nrel.gov/wind/windpact/, Golden, CO: National Renewable Energy Laboratory,
January 2005.
15
Newman, J.N., Marine Hydrodynamics, Cambridge, MA: The Massachusetts Institute of
Technology Press, 1977.
16
REpower Systems, REpower 5M,
http://www.repower.de/typo3/fileadmin/download/produkte/5m_uk.pdf, Hamburg, Germany:
REpower Systems, January 2005.
17
REpower Systems, REpower Systems AG Renewable Energy for the Future,
http://www.repower.de/, Hamburg, Germany: REpower Systems, January 2005.
18
Smith, Kevin, WindPACT Turbine Design Scaling Studies; Technical Area 2: Turbine, Rotor,
and Blade Logistics, NREL/SR-500-29439, Golden, CO: National Renewable Energy
Laboratory, June 2001.
19
Tarp-Johansen, Neils Jacob, RECOFF Home Page, http://www.risoe.dk/vea/recoff/, Roskilde,
Denmark: Ris National Laboratory, July 2004.
23