17.08.2013 Views

Publishers version - DTU Orbit

Publishers version - DTU Orbit

Publishers version - DTU Orbit

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

encompasses the storage of preview measurements. If the gains associated with the preview<br />

storage are partitioned from the gains applied to turbine measurements, then the former<br />

are correctly viewed as disturbance feedforward control and this is how they are depicted in<br />

Fig. 148. In the course of design, it is determined that the majority of benefit of preview is<br />

obtained using 0.5 s (10 samples) worth of feed-forward compensation, but this can change<br />

depending on the dynamics of the turbine involved.<br />

The FAST simulation code marches the TurbSim wind distribution past the turbine at<br />

the mean speed of 18 m/s. The pre-processing induces evolution corresponding to distances<br />

iTs×18 where Ts is the sample period and i is the number of preview samples used. Where<br />

the evolution distance is greater than 10Ts×18, the process is equivalent to taking preview<br />

measurements further in front of the turbine, and then waiting the appropriate amount of<br />

time before applying the preview gains to the measurements.<br />

10.7.1 H2 Optimal Preview Control<br />

In this section, we describe the design of the preview controller and then evaluate its performance<br />

in the presence of emulated wind evolution. Linearmodels are obtained from the FAST<br />

wind turbine modeling code (Jonkman and Buhl, 2005) developed at NREL. The model is<br />

based on a 40 m diameter, three-bladed controls advanced-research turbine (CART3) located<br />

at NREL’s National Wind Technology Center (NWTC). The nominal operating point for design<br />

(and simulation) is a uniform wind of 18 m/s, a blade pitch of 12.7 ◦ , and a rotor speed<br />

of 41.7 rpm.<br />

As in Laks et al. (2011b), the FAST linearizations are used to obtain a discrete-time statespace<br />

model representing the turbine dynamics Pt with a 20 Hz sample rate. For controller<br />

design,the linearizedmodel includes a first-ordergeneratordegree of freedom (DOF), secondorder<br />

dynamics for each blade’s out-of-plane blade flap compliance, and a second-order drivetrain<br />

compliance. This model is then augmented with simple pitch actuator dynamics that<br />

provide the pitch rates generated by individual pitch commands. During simulation, all DOFs<br />

provided by FAST are used except for yaw and teeter (the freedom of the rotor to tilt).<br />

In addition, integral control on generator speed error Ωg (= ΩRNgb where Ngb is the gearbox<br />

ratio and ΩR is the rotor speed) and rejection of 1P variation in the bending moments are<br />

obtained by augmenting the turbine model with additional dynamics Pa driven by measurements<br />

of these signals. We configure the control system to use individual point measurements<br />

of the longitudinal wind that the blades will encounter at 75% span after a fixed time delay<br />

of i samples, where the point measurements rotate in unison with the blades. As discussed in<br />

Laks et al. (2011b), this corresponds to using a blade-local model, that describes the effect<br />

on the turbine of wind speed perturbations local to each blade.<br />

Thestatefeedbackisdesignedbyformingthegeneralizedplant(theFASTlinearizationwith<br />

preview storage and augmented dynamics) with weighting on outputs that include generator<br />

speed perturbations and perturbations in out-of-plane blade bending moments. The statefeedback<br />

gains [Kx Ka Kff] are optimized to minimize the H2 norm from the preview wind<br />

measurement to weighted <strong>version</strong>s of flap and speed perturbations; this includes a penalty on<br />

excessive pitch rates so that the associated linear quadratic regulator problem is well posed. It<br />

is possible to compute the optimal state feedback independent of the amount of preview used<br />

(Hazell, 2008). This involves using the stabilizing solution of a Riccati equation to compute<br />

both the optimal state-feedback and feedforward gains. Details can be found in Laks et al.<br />

(2013).<br />

In formulating the H2 cost, emphasis is initially focused on the flap response to wind<br />

perturbations. Then, the penalty on pitch effort is increased until the (linear) closed-loop<br />

response to a step change in collective wind produces pitch rates on the order of 10 ◦ /s.<br />

Generally, the H2 performance improves with the number of samples of preview, but remains<br />

bounded below so that there is a diminishing return; this lower bound is essentially reached<br />

with the use of about 4 samples of preview at 20 Hz. However, the goal is not the precise<br />

value of the H2 norm, but the attenuation of perturbations in blade load due to perturbations<br />

212 <strong>DTU</strong> Wind Energy-E-Report-0029(EN)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!