POLItecnico
di MIlano
ESTIMATION OF DAMPING FOR
WIND TURBINES OPERATING IN
CLOSED LOOP
C.L. Bottasso, S. Cacciola, A. Croce
Politecnico di Milano, Italy
S. Gupta
Clipper Windpower Inc., USA
EWEC 2010
Warsaw, Poland, April 20-23, 2010
Damping Estimation of Wind Turbines
Outline
• Introduction and motivation
• Approach: modified Prony’s method for linear time
periodic systems
• Applications and results:
- Simulation models
- Library of procedures for modes of interest
- Examples: tower, rotor and blade modes
• Conclusions and outlook
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Introduction and Motivation
Focus of present work: estimation of damping in a wind turbine
Applications in wind turbine design and verification:
•
•
•
•
Explaining the causes of observed vibration phenomena
Assessing the proximity of the flutter boundaries
Evaluating the efficacy of control laws for low-damped modes
…
Highlights of proposed approach:
•
•
•
Closed loop: damping of coupled wind turbine/controller system
Applicable to arbitrary mathematical models (e.g., finite element
multibody models, modal-based models, etc.)
In principle applicable to a real wind turbine in the field
POLITECNICO di MILANO
Poli-Wind Research Lab
Introduction and Motivation
Damping Estimation of Wind Turbines
Previous work:
•
Linear Time Invariant (LTI) systems:
Hauer et al., IEEE TPS, 1990; Trudnowski et al., IEEE TPS 1999
•
However: wind turbines are characterized by periodic coefficients
(vertical/horizontal shear layer, up-tilt, yawed flow, blade-tower
interaction, etc.)
Linear Time Periodic (LTP) systems:
Bittanti & Colaneri, Automatica 2000; Allen IDETC/CIE 2007
However: methods well suited only when
characteristic time τ (time to half/double)
much larger than period T (1rev): τ ≫T
τ2 0,96 sec, 2nd fore-aft tower mode
τ1 3.45 sec, 1st fore-aft tower mode
T
5.5 sec
Typically not the case for WT problems
E.g.: damping of tower fore-aft modes ▶
Proposed approach: transform LTP in equivalent/approximate LTI, then use
Prony’s method (standard for LTI analysis)
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Outline
• Introduction and motivation
• Approach: modified Prony’s method for linear time
periodic systems
• Applications and results:
- Simulation models
- Library of procedures for modes of interest
- Examples: tower, rotor and blade modes
• Conclusions and outlook
POLITECNICO di MILANO
Poli-Wind Research Lab
Approach
LTP system:
Damping Estimation of Wind Turbines
.
x = A(ψ)x + B(ψ)u
A(ψ) = closed-loop matrix (accounts for pitch-torque controller)
u = exogenous input (wind), constant in steady conditions
Fourier reformulation (Bittanti & Colaneri 2000):
A(ψ) = A0+Σi(Aissin(i ψ)+Aiccos(i ψ))
B(ψ) = B0+Σi(Bissin(i ψ)+Biccos(i ψ))
1. Approximate state matrix: A(ψ) ≈ A0
2. Transfer periodicity to input term (remark: arbitrary amplitude)
Obtain linear time invariant (LTI) system:
.
x = A0x + Ub(ψ)
where b(ψ) = exogenous periodic input
Remark: no need for model generality, just good fit with measures
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Approach
Given reformulated LTI system
.
x = A0x + Ub(ψ)
use standard Prony’s method (Hauer 1990; Trudnowski 1999):
1. Trim and perturb with doublet (or similar, e.g. 3-2-1-1) input
2. Identify discrete time ARX model (using Least Squares or Output
Error method) with harmonic input
3. Compute discrete poles, and transform to continuous time (Tustin
transformation)
4. Obtain frequencies and damping factors
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Outline
• Introduction and motivation
• Approach: modified Prony’s method for linear time
periodic systems
• Applications and results:
- Simulation models
- Library of procedures for modes of interest
- Examples: tower, rotor and blade modes
• Conclusions and outlook
POLITECNICO di MILANO
Poli-Wind Research Lab
Simulation Models
Cp-Lambda highlights:
Damping Estimation of Wind Turbines
• Geometrically exact composite-ready
beam models
• Generic topology (Cartesian
coordinates+Lagrange multipliers)
• Dynamic wake model (Peters-He,
yawed flow conditions)
• Efficient large-scale DAE solver
• Non-linearly stable time integrator
• Fully IEC 61400 compliant (DLCs,
wind models)
Compute
sectional
stiffness
• Rigid body
ANBA (Anisotropic Beam Analysis)
cross sectional model
• Geometrically exact beam
Recover cross
sectional
stresses/strains
POLITECNICO di MILANO
• Revolute joint
• Flexible joint
• ActuatorPoli-Wind Research Lab
Damping Estimation of Wind Turbines
Simulation Environment
Measurement
noise
POLITECNICO di MILANO
Wind
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Applications and Results
Definition of best practices for the identification
of modes of interest:
For each mode:
Excitations (inputs)
• Consider possible excitations (applied loads,
pitch and/or torque inputs) and outputs (blade,
shaft, tower internal reactions)
• Verify presence of modes in response (FFT)
• Verify linearity of response
• Perform model identification
• Verify quality of identification (compare
measured response with predicted one)
Compiled library of mode id procedures:
In this presentation:
• Tower fore-aft mode
• Rotor in-plane, blade first edge modes
POLITECNICO di MILANO
Response (outputs)
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Example: Damping Estimation of
Fore-Aft Tower Modes
Doublets of varying intensity to verify linearity
Excitation: doublet of hub
force in fore-aft direction
Verification of linearity of response
Output: tower root
fore-aft bending
moment
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Example: Damping Estimation of
Fore-Aft Tower Modes
Verification of linearity of response and presence of modes
First tower mode
1P
POLITECNICO di MILANO
Second tower mode
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Example: Damping Estimation of
Fore-Aft Tower Modes
◀ Time domain
▼ Frequency domain
• Excellent quality of identified models
(supports hypothesis A(ψ) ≈ A0)
• Necessary for reliable estimation
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Example: Damping Estimation of
Fore-Aft Tower Modes
Estimated damping ratios for varying wind speed
POLITECNICO di MILANO
Poli-Wind Research Lab
Example: Damping Estimation of
Blade Edge and Rotor In-Plane Modes
Damping Estimation of Wind Turbines
Quality of identified model, using blade root bending
First blade
edgewise mode
Excitation: doublet of
• In-plane blade tip force
• Generator torque
Rotor in-plane
mode
Quality of identified model, using shaft torque
Outputs:
• Blade root bending moment
• Shaft torque
Rotor in-plane
mode
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Example: Damping Estimation of
Blade Edge and Rotor In-Plane Modes
Rotor in-plane mode
◀ Little sensitivity to used output
(blade bending or shaft torque)
Blade edge mode
POLITECNICO di MILANO
Poli-Wind Research Lab
Damping Estimation of Wind Turbines
Outline
• Introduction and motivation
• Approach: modified Prony’s method for linear time
periodic systems
• Applications and results:
- Simulation models
- Library of procedures for modes of interest
- Examples: tower, rotor and blade modes
• Conclusions and outlook
POLITECNICO di MILANO
Poli-Wind Research Lab
Conclusions
Damping Estimation of Wind Turbines
Proposed a method for the estimation of damping in wind turbines:
• Modified Prony’s method (accounts for periodic nature of wind turbine
models)
• Good quality model identification is key for reliable damping
estimation
• Compiled library of mode id procedures (need specific inputs/outputs
for each mode)
• Fast and robust
Outlook:
• Riformulation leading to Periodic ARX, and comparison
• Effect of turbulence (simulation study):
- Turbulence as an excitation
- Turbulence as process noise (filter error method)
• Verify applicability in the field (theoretically possible)
POLITECNICO di MILANO
Poli-Wind Research Lab
Scarica

estimation of damping for wind turbines operating in closed loop