3D-Multi Body Simulation of Wind Turbines With Flexible Components
3D-Multi Body Simulation of Wind Turbines With Flexible Components
3D-Multi Body Simulation of Wind Turbines With Flexible Components
A new method and an example of a 3-dimensional simulation of the dynamic behaviour of wind tur-
bine generator systems (WTGS) is presented. The essential feature is a multi body simulation of a
complete 3-D model of the wind turbine including flexible elements (for example rotor and tower).
By reducing the analysis to standard methods in multi body simulations (MBS) and finite element
methods (FEM) the simulation model avoids programming errors. Both, motion and forces at any
point of interest within the turbine can be represented by 3D animation, graphs or result files for
advanced analysis. The 3-D animation guarantees an easy check of the reliability. The calculated
example takes into account the FEM structures of the entire rotor blades and the tower and their
flexible properties. The inertia effects of all structural elements (rotor, drive train, gearbox and gen-
erator) are implemented. A pitch and yawing control as well as a characteristic torque of the genera-
tor are modelled. The induced wind velocities and the local interaction of airflow and vibrations of
the blades and the tower are calculated. The developed 3-D simulation permits a safe realistic fore-
cast of normal operation- and extreme loads. The paper is part of a project supported by the gov-
ernment department for education and research (BMBF) of the Federal Republic of Germany.
1. Object of 3D-Simulation
In present simulation programs different types of wind
turbines, especially different foundations of off-shore
facilities, are modelled either by rough simplification
or by additional programming susceptible to faults.
Using flexible FEM-structures (for example founda-
tion, tower, blades) within the simulation gives a bet-
ter approach to reality. Local strain concentrations at
discontinuities (end of struts, ribs) are covered realis-
tically. In dynamic simulations the flexible blades are
substituted up to now only by beams. The detailed
structure of the blades and the local interaction of air-
flow and vibrations of the blades are often neglected.
Most of the present simulation programs for WTGS
permit only a limited optical check of the used calcu-
lation model, expressive 3D-animations are missing.
A 3D-multi body simulation of WTGS can solve these
problems (Fig. 1).
It is not possible to model wind turbine generator systems using geometric based tools only.
Aerodynamic loads must be calculated by additional programming. For that purpose undisturbed wind
fields (IEC-61400-1 [11]) are processed from input data.
2. Flexible Components and Rigid Bodies
A flexible body in MBS is defined by a linear dynamic finite element model whose normal modes and
static costraint modes have to be taken into account. Flexible properties are calculated by the finite ele-
ment program ANSYS [10]. In order to build the multi body simulation model DADS [9] is used.
2.1 Multi body simulation, rigid bodies
A mechanism within a multi body simulation is described geometrically, the associated equations of
motion are set up automatically and are solved in time domain. The results are represented graphical-
ly. A mechanism consists of (rigid) bodies, whose relative motions are restricted by constraints. Every
body has a body fixed co-ordinate system. Its origin and orientation is described by independent gen-
eralised co-ordinates {k}. The joints connecting the bodies are modelled by MBS standard elements.
They define kinematic constraints in the following form:
63
3D-Multi Body Simulation of Wind Turbines with Flexible
Components
F. Baumjohann, M. Hermanski, R. Diekmann, University of Applied Sciences Bielefeld;
J. Krning, TV Nord
DEWI Magazin Nr. 21, August 2002
Fig. 1: Partial view of the wind turbine
Externer Artikel
64
The program calculates the mass matrix M, Jacobian Z,k, a term (t,k,l) resulting from the constraints
for k and uses the applied forces Q
A
(t,k,l) and Lagrange multipliers for constraint forces in order to
formulate the equations of dynamics
2.2 Flexible bodies, normal and constraint modes
In general all important flexible components of a WTGS are available as FEM models. The presented
example of a blade [2] has been created parametrically. Local flexible displacements u are combina-
tions of fixed interface normal modes n and constraint modes c , which are determined within the
FEM code. Normal modes (Fig.2) with fixed displacements at the interfaces are calculated from the
eigenvalue problem
Constraint modes c are static deformations with normalized displacements at nodes coupling to the
rigid bodies. The modes have to be transformed, so that they are orthogonal and decouple the equa-
tions of motion [6]. Moreover rigid body modes have to be removed. Both is done in the multi body
code. The modal coordinates qn(t), qc(t) define the local displacements of the flexible components
where unchanged terms for the transformed modes were used.
2.3 Flexible bodies in a moving mechanism
The descrition of flexible bodies in MBS will be illus-
trated with the example of a beam (Fig. 3). The
motion xN(t) of a node N of the flexible beam is
described by superposing its rigid body motion and
flexible relative displacement uN represented by the
modal coordinates q. The rigid body motion of the
node is expressed by generalised coodinates of posi-
tion R and orientation p of the undeformed structure
and the relative position rN in the body-fixed reference
frame.
I denotes the identity, E(rN) p describes a rotation using Euler parameters p and A is a transformation
to global co-ordinates.
The rigid-body and flexible-element motions influence each other according to the following equation
of dynamics:
DEWI Magazin Nr. 21, August 2002
Z(k,t)} =
) t , (
h
z
) t , (
1
z
k
k
M = 0.
]
]
]
,
,
0 Z,
Z, M
k
T
k
k
& &
=
A
Q
.
[-
2
M + K]
n
= 0.
u(t) = [
n
,
c
]
(t)
(t)
c
n
q
q
,
Fig. 2: Third normal mode of the blade
x
N
(t) = R(t) + r
N
(t) + u
N
(t) = [I; E(r
N
); A]
) t (
) t (
) t (
q
p
R
[ ] +
]
]
]
]
,
,
T
k
A
, Z Q
q
p
R
qq
K 0 0
0 0 0
0 0 0
q
p
R
M
& &
& &
& &
Externer Artikel
3 Modelling Steps
The necessary steps to work out a simulation of a wind
turbine under operating conditions are described in the
following chapters.
3.1 Mechanical system
The blades and the tower have been built as parame-
terised finite element models. The fixed interface nor-
mal modes and the constraint modes are calculated by
ANSYS. The finite element model and the computed
modes are transferred via the interface implemented in
the MBS code. The rigid components (hub, shaft,
nacelle), their joints and connections are modelled with
standard DADS elements. The aerodynamic coeffi-
cients are fed into the program and linked with the
blade data.
3.2 Control and electrical system
The pitch control of the example is modelled by standard elements of DADS. It is also possible to cre-
ate individual control algorithms with user defined subroutines. To simulate the electrical load of the
wind turbine a file containing the generator characteristic curve is used as input to a standard control
element.
3.3 Loads
The undisturbed turbulent wind field is computed in a
separated program and has been made available for
the MBS model. In an internal user defined subrou-
tine the induced wind is calculated for each segment
of the blade. The relative local wind speed takes into
account the vibration of the tower and the blades. The
forces calculated by the aerodynamic coefficients are
transferred automatically to the corresponding finite
element nodes and are used at the next time step.
Thereby a direct aeroelastic coupling is coming
about. The time step of the solver is automatically
chosen by the program, the print interval is deter-
mined by the user.
An input file controls different load cases, e.g. normal
operation, start up and shut down, gusts or opera-
tional conditions resulting from a failure.
4 Evaluation of Results
4.1 3D-animation, graphical evaluation
The dynamic behaviour can be displayed as anima-
tion. The view point can be any position in a local or
global co-ordinate system. The viewer can move with
any body of the wind turbine model. The reaction
forces and moments at any position of the model are
saved in a result file and can be used in further analy-
sis.
4.2 Stress evaluation
Strain or stress in the flexible elements of the wind turbine are calculated by the orthogonalised deflec-
tion shapes used in the MBS formulation. The deflections are transferred back to ANSYS as static load
cases. Additionally the time dependent modal co-ordinates are transferred to the FE-code and are used
65
DEWI Magazin Nr. 21, August 2002
X
G
Y
G
Y
C
Global co-ordinate system
X
C
u
N
R
x
N
r
N
flexible beam with
6 finite elements
local FEM
-ordinate system
Fig. 3: Superposing rigid body motion and flexible rela-
tive displacements
Fig. 4: Snap-shot of the deflection; left side scaled by fac-
tor 25; right side scaled by factor 3
Fig. 5: Force components at the top of the tower
Externer Artikel
66
to superimpose the shapes describing the deflection
of the finite element nodes. An ANSYS specific
APDL program has been developed to calculate
stresses and strains at arbitrary time steps. The
resulting time histories can be used to analyse
extreme loads and operational stability.
5 Summary
A new technique for analysing the dynamics and
stresses of wind turbines is presented. The simula-
tion is carried out with a complete 3-dimensional
model of the wind turbine within a multi body simu-
lation program using flexible elements. Models of
different types of blades or towers need not to be
simplificated. Modelling is supported by optical
checks. The analysis is reduced to standard meth-
ods in multi body simulations and finite element cal-
culation. Strain concentrations at discontinuities
within the flexible parts are covered realistically. The
systems behaviour can be represented by 3D ani-
mation of the moving bodies. Forces or stresses (in
the flexible parts) may be displayed at any point of
interest within the turbine. The new method is to be
extended to a software package for the simulation of
wind turbines.
References:
[1] Krning J., Baumjohann F.: Zur Absicherung der Betriebsfestigkeit von Rotorblttern,
DEWEK98, 4. Deutsche Windenergie-Konferenz, Wilhelmshaven, 1998, 215-219
[2] Baumjohann F., Korjukin E., Krning J.: Lastsimulation und parametrisierte FE-Berechnung
eines Windenergierotors. Proc. 2nd. Int. Scient. Colloquium CAE Techniques, Bielefeld, 1995.
[3] Edward J. Haug: Computer-Aided Kinematics and Dynamics of Mechanical Systems, Volume I,
Allyn And Bacon, ISBN 0-205-11669-8.
[4] Wan S. Yoo: Dynamics of flexible mechanical systems using finite element lumped mass approx-
imation and static correction modes, Abstract of doctor Thesis, University of Iowa 1985
[5] Olof Friberg; Vesa Karhu: Use of mode orthogonalization and modal damping in flexible multi-
body dynamics, Department of Structural Mechanics; Chalmers University of Technology; S-412
96 Gteborg, Sweden
[6] Roy R. Craig: Structural Dynamics, John Wiley & Sons, Inc. 1981
[7] R.E Roberson, R. Schwertassek: Dynamics of multibody systems, Springer Verlag 1988
[8] Richard Schwertassek, Oskar Wallrapp: Dynamik flexibler Mehrkrpersysteme, Vieweg & Sohn
Braunschweig/Wiesbaden 1998, ISBN 3-528-06629-6
[9] DADS Manual Revision 9.0: Dads/ Flex Manual, CADSI, 1999
[10] ANSYS Release 5.7.1, ANSYS, Inc.
[11] IEC Standard 61400-1, Wind turbine generator systems, 1999-02
DEWI Magazin Nr. 21, August 2002
Fig. 6: Stress and deformation of a blade
Externer Artikel