Commit 0d9dec59 authored by famen's avatar famen
Browse files

Dev floatstep

parent 0be69428
Pipeline #24073 passed with stage
in 3 minutes and 32 seconds
......@@ -488,10 +488,10 @@
constant 88 0.0; Frequency of BP filter [Hz] (Default 0 if filter not used)
constant 89 0.02; Damping ratio of BP filter [-] (Default 0.02)
constant 90 0.0; Time constant of BP filter [s] (Default 0)
cosntant 91 0.0; Coefficient of linear term in gain scheduling for tower-pitch loop, KK1 [1/(m/s)];
constant 91 0.0; Coefficient of linear term in gain scheduling for tower-pitch loop, KK1 [1/(m/s)];
; Kgain * (1 + KK1 *abs(WSPfilt - WSPrated) + KK2 * abs(WSPfilt - WSPrated)**2)
constant 92 0.0; Coefficient of quadratic term in gain scheduling for tower-pitch loop, KK2 [1/(m/s)^2];
cosntant 93 0.0; Coefficient of linear term in gain scheduling for tower-genTorq loop, KK1 [1/(m/s)];
constant 93 0.0; Coefficient of linear term in gain scheduling for tower-genTorq loop, KK1 [1/(m/s)];
constant 94 0.0; Coefficient of quadratic term in gain scheduling for tower-genTorq loop, KK2 [1/(m/s)^2];
constant 95 1.0; Choice of gain-scheduling variable (0: WSPfilt, 1:Pitch angle (Default))
end init;
......
......@@ -49,12 +49,14 @@ subroutine floating_controller(CtrlStatus, time, Floatingvar, TTfa_vel, wsp, Pit
TTfa_vel_filt = y(1)
!!******************
! filterd wind speed for gain-scheduling
WSPfilt = lowpass1orderfilt(deltat, stepno, wspfirstordervar, wsp)
if (Floatingvar%GSmode .gt. 0.0_mk) then
Floatingvar%GSvar = abs(WSPfilt - Floatingvar%RatedWindSpeed)
else
if (Floatingvar%GSmode .gt. 0.0_mk) then
Floatingvar%GSvar = PitchColRef
else
WSPfilt = lowpass1orderfilt(deltat, stepno, wspfirstordervar, wsp)
Floatingvar%GSvar = abs(WSPfilt - Floatingvar%RatedWindSpeed)
endif
!!******************
! Main Loop
! Tower-pitch loop
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment