How to model a wind turbine

A wind turbine in SIMA is described by structural components, a wind turbine model component, and a controller. A recommended approach for modelling the structural components is described in this guide. For the wind turbine component and controller, see the RIFLEX wind turbine documentation.

For the structural components, RIFLEX lines are used for distributed properties, while lumped properties can be modelled using SIMO bodies and/or RIFLEX nodal components. An overview of the different components of a wind turbine and the recommended model type is given in Table 1. An example of a turbine model generated using this philosophy can be found at the SIMA GitHub repository.

Table 1. Components of a wind turbine

Component

Model type

Alternative model type

Blade

RIFLEX line

Blade eccentricity

RIFLEX line

Hub

RIFLEX nodal component

SIMO body

Shaft

RIFLEX line

Nacelle

SIMO body

RIFLEX lines

Generator

RIFLEX Flex Joint

Tower

RIFLEX line

The components are now described in more detail in the following sections, with reference to Figure 1 to Figure 3. Figure 1 shows a sketch of the wind turbine, while Figure 2 shows the supernodes and lines used in the model. This is consistent with the specification in Table 1. Figure 3 shows a close-up view of the shaft, as well as one eccentricity line and the innermost part of the corresponding blade.

Turbine outline
Figure 1. Sketch of a wind turbine
Supernodes lines
Figure 2. Lines and supernodes in model
um ii fig48
Figure 3. Lines and supernodes at nacelle and hub

1. Supernode locations and boundary conditions

Supernodes form the basis of RIFLEX models, defining the start and end points of each structural component. The supernodes given in Table 2 are required to model a wind turbine. Additional supernodes may be required to model e.g. the suport structure and foundation.

Table 2. Supernodes in a wind turbine

Node type

Number of nodes

Boundary condition

Location

Shaft (hub-end node)

1

Free

Hub centre

Shaft (tower-end node)

1

Slave to tower top/Fixed

Shaft element length behind shaft hub-end node, taking shaft tilt angle into account

Blade eccentricity

N blade

Free

Hub radius from shaft hub-end node, at blade root. Should be along blade pitch axis.

Blade root

N blade

Slave to corresponding blade eccentricity super node

Same location as blade eccentricity super nodes

Blade tip

N blade

Free

At blade tip projected onto the pitch axis. Rotor radius distance from shaft hub-end node

2. Lines and line types – blade lines, ecc lines, shaft line

Turbine blades, eccentricity lines and the shaft are modelled as lines, where a line represents each structural component, while the properties of each line are described by a line type. Typically, one line type is needed to represent the blade lines, one for the eccentricity lines, and one for the shaft line. Only in case e.g. the blades have different properties, multiple blade line types are needed. Typical properties of each line type are outlined in Table 3.

Table 3. Properties of line types

Line type

Cross section type

Properties

Blade line type

Axisymmetric cross section/general cross section

According to blade specification, see Blade modelling.

Eccentricity line type

Axisymmetric cross section

Rigid (high stiffness value), massless (negligible mass).

Shaft line type

Axisymmetric cross section

Rigid (high stiffness value). Mass as described in Shaft modelling.

Tower

Axisymmetric cross section/Thin-walled pipe cross section

According to tower specification.

The line type and supernode connections of each line are given in Table 4.

Table 4. Line properties with line types and supernodes

Line

Line type

Supernode end 1

Supernode end 2

Shaft line

Shaft line type

Shaft (hub-end node)

Shaft (tower-end node)

N blade Eccentricity lines

Eccentricity line type

Shaft (hub-end)

Blade i eccentricity supernode

N blade Blade lines

Blade line type

Blade i blade root

Blade i blade tip

Tower line

Tower line type

Tower base

Tower top

3. Blade modelling

While location of the supernodes takes the rotor cone and shaft tilt angles into account, the prebend, sweep and twist of the blades are modelled in the line type. Prebend is modelled as Offset Z, sweep is modelled as Offset Y, and twist is modelled as Twist end1 and Twist end2. It is recommended that Offset Y and Offset Z gives the distance from the pitch axis to the reference line used to model the structural properties, while Twist end1 and Twist end2 describes the angle of the chord line used to define the aerodynamic properties, relative to the element coordinate system as defined at the blade root. An example is shown in the figure below, where the dash-dot line is the pitch axis, and the solid line represents the reference line for defining the cross section properties.

sima blade
Figure 4. Blade structural model

Blade mass, stiffness and aerodynamic properties are defined for the segments in the blade line type using cross-sections. The properties defined in the cross-sections are given relative to the reference line described by Offset Y, Offset Z, and Twist end1/Twist end2. Either the general cross section or the double symmetric cross section can be used. In case of the general cross section being used, the principal axis orientation is given relative to the chord line. For the airfoil properties, the coordinates of the reference point for the aerodynamic coefficients must be given. This is defined using Foil origin Y \((Y_{FC})\) and Foil origin Z \((Z_{FC})\).

An example is shown in the figure below, using the general cross section. The origin of the local coordinate system is the reference line defined in the previous section.

BladeCS
Figure 5. Blade cross section modelled using the general cross section

4. Shaft modelling

The shaft is modelled as a RIFLEX line with two segments, with the generator torque being applied at the interface between the two segments. The shaft segment closest to the rotor is rotating and should contain the equivalent inertia of rotating parts of the shaft and generator. Note that this must contain the inertia of both the low-speed and high-speed drivetrain components, where the inertia of the latter should be scaled by the square of the gearbox ratio. The shaft segment closest to the tower transfers structural loads from the rotor and generator to the support structure.

5. Hub modelling

The hub mass and inertia are modelled using a SIMO body or nodal body attached to fore shaft node.

6. Nacelle modelling

The nacelle is modelled as a SIMO body using a slender system connection. Any property of a SIMO body can be defined, but the most important are the mass, mass moment of inertia and centre of gravity. Note that the mass moments of inertia are given relative to the local origin of the SIMO body, not the centre of gravity. Coupling terms are therefore expected.