Drijvende Offshore Windturbine: Actieve Golfdemping

Drijvende offshore windturbines (FOWT) vormen het speerpunt van duurzame energie op diepe kustwateren. In tegenstelling tot turbines met een vaste fundering die in ondiep water (< 50 m) op de zeebodem worden geheid, drijven spar-buoy platforms zoals de NREL OC3-Hywind op waterdieptes van 100 tot meer dan 300 meter.

Op open zee zijn de windsnelheden hoger en constanter, maar de dynamische uitdagingen zijn aanzienlijk. Een drijvende spar-buoy gedraagt zich in rotatie (pitchhoek $\theta$ ) als een enorme onderwaterslinger. Doordat de golfstralingsdemping van zeewater uiterst gering is ( $\zeta \approx 0{,}03$ ), veroorzaakt golfexcitatie bij of rond de eigenfrequentie gevaarlijke opslingering die de toren en constructie overbelast.

Deze casestudy onderzoekt actieve golfdemping via een gereduceerd 1-DOF spar-buoy pitchmodel, toont aan waarom een naïeve PD-regelaar faalt, en demonstreert hoe een tweede-orde notchfilter de fundamentele frequentiekloof overbrugt.

De fysische opstelling

Het model is gebaseerd op het benchmark NREL OC3-Hywind concept. In rotatierichting (bewegingsgraad 5 in de scheepsbouw) wordt het evenwicht bepaald door de Newton-Euler rotatiebalans:

I_{\text{tot}} \ddot{\theta} + B_{55} \dot{\theta} + C_{55} \theta = M_{\text{wave}} - M_{\text{ctrl}}

Parameter	Betekenis	Waarde
$I_{55}$	Traagheidsmoment stalen constructie	$5{,}00 \times 10^{10}\ \text{kg}\cdot\text{m}^2$
$A_{55}$	Hydrodynamische toegevoegde massa (meebewegend water)	$1{,}70 \times 10^{10}\ \text{kg}\cdot\text{m}^2$
$I_{\text{tot}}$	Totale virtuele traagheid ( $I_{55} + A_{55}$ )	$6{,}70 \times 10^{10}\ \text{kg}\cdot\text{m}^2$
$B_{55}$	Hydrodynamische stralingsdemping	$8{,}65 \times 10^8\ \text{N}\cdot\text{m}\cdot\text{s/rad}$
$C_{55}$	Hydrostatische herstelstijfheid	$3{,}111 \times 10^9\ \text{N}\cdot\text{m/rad}$
$f_n$	Natuurlijke pitch-eigenfrequentie	$0{,}0343\ \text{Hz}\ (\omega_n = 0{,}2155\ \text{rad/s},\ T_n \approx 29{,}2\ \text{s})$
$\zeta$	Onbelaste dempingsratio	$\approx 0{,}03$ ( $Q \approx 16{,}7$ , $+24{,}5\ \text{dB}$ resonantiepiek)

De vrijheidsgraden ("55" conventie)

In de maritieme techniek en hydrodynamica heeft elk drijvend voorwerp 6 standaard vrijheidsgraden:

Surge (langsscheepse verplaatsing)
Sway (dwarsscheepse verplaatsing)
Heave (verticale verplaatsing / dompen)
Roll (slingeren om de lengteas)
Pitch (stampen / kantelen om de dwarsas)
Yaw (gieren om de verticale as)

De index $55$ staat voor de reactiekracht in richting 5 (pitch) veroorzaakt door beweging in richting 5 (ongekoppelde pitch-dynamica).

Overdrachtsfunctie & De Spectrale Scheiding

De open-lus overdrachtsfunctie van moment naar pitchhoek $\theta$ luidt:

G(s) = \frac{\Theta(s)}{M(s)} = \frac{1}{I_{\text{tot}} s^2 + B_{55} s + C_{55}}

In het frequentiedomein zien we:

De constructie-eigenfrequentie van de spar ligt op $\omega_n = 0{,}2155\ \text{rad/s}$ ( $f_n = 0{,}0343\ \text{Hz}$ ).
Dominante zeegang concentreert zich rond $f_p \approx 0{,}10\ \text{Hz}$ ( $\omega_p = 0{,}628\ \text{rad/s}$ , periode $T_p \approx 10\ \text{s}$ ).

Omdat de golffrequentie bijna een factor 3 hoger ligt dan de structuurresonantie ( $\omega_p / \omega_n \approx 2{,}91$ ), is er sprake van een duidelijke frequentiekloof:

Constructieresonantie is traag ( $0{,}034\ \text{Hz}$ ).
Golfexcitatie is sneller ( $0{,}10\ \text{Hz}$ ).

Waarom een naïeve PD-regelaar faalt

Een conventionele PD-regelaar gebruikt differentiërende actie ( $K_d s$ ). De versterking van de D-actie groeit met $+20\ \text{dB/decade}$ :

|K_d \cdot j\omega| = K_d \cdot \omega

Bij de golffrequentie ( $0{,}10\ \text{Hz}$ ) beantwoordt de D-actie elke radiaal golfgeïnduceerde kanteling met $\sqrt{K_p^2 + (K_d\omega)^2} \approx 561\ \text{MN}\cdot\text{m}/\text{rad}$ tegenkoppel, terwijl de golven niet datgene zijn wat hij moet dempen. Met de golven van deze les ( $20\ \text{MN}\cdot\text{m}$ ) is het piekkoppel slechts ongeveer $1{,}25\ \text{MN}\cdot\text{m}$ — $2{,}5\%$ van de grens van $\pm 50\ \text{MN}\cdot\text{m}$ , die pas bij een kanteling op golffrequentie van ongeveer $5{,}1^\circ$ bereikt wordt. Bijna al die inspanning gaat op aan de golffrequentie, waar zij niets doet voor de resonantie; de PD verhoogt de resonantiedemping wel van $\zeta = 0{,}030$ naar ongeveer $0{,}054$ .

Actieve golfdemping met een notchfilter

De oplossing bestaat uit het plaatsen van een 2e-orde notchfilter (bandsperfilter) $N(s)$ rond de golffrequentie ( $f_0 = 0{,}10\ \text{Hz}$ ) in de terugkoppellus:

N(s) = \frac{s^2 + \omega_0^2}{s^2 + 2\zeta_n \omega_0 s + \omega_0^2} \quad \text{met } \omega_0 = 2\pi f_0

Bij de golffrequentie ( $s = j\omega_0$ ) wordt de teller exact nul:

|N(j\omega_0)| = 0 \quad (-\infty\ \text{dB onderdrukking})

Het filter maakt de regelaar blind voor de eerste-orde golfbeweging, zodat het platform natuurlijk meebeweegt met de kammen: in deze les dalen het piekkoppel en de gedissipeerde energie met ruwweg 70%. Gratis is dat niet: bij de constructieresonantie ( $0{,}034\ \text{Hz}$ ) is $|N| \approx 0{,}93$ , maar het filter vertraagt de terugkoppeling daar $21^\circ$ , waardoor de demping die de PD toevoegt met ongeveer 70% krimpt (gesloten lus $\zeta \approx 0{,}037$ met notch, $0{,}054$ zonder, $0{,}030$ open lus).

Reduced-Order afbakening & Visuele context

Deze les modelleert de ongekoppelde pitchbeweging (1-DOF gereduceerd spar-buoy model). De windturbine-rotor, meertrossen en deining dienen als visuele context om fysisch inzicht te vergroten, zonder de didactische essentie te vertroebelen met zware hydro-elastische simulaties.

Zie ook

StateSpace: compacte 2-toestandsweergave van de spar-buoy hydrodynamica.
TransferFunction: implementatie van het 2e-orde notchfilter.
Integrator: integratieketen van hoekversnelling naar pitchhoek.
Saturation: realistische koppelbegrenzing van de actuator.
BodePlot: spectrale analyse en frequentiekloof in het frequentiedomein.
uPlotDisplay: real-time telemetrie voor hoek, hoeksnelheid en stuurkoppel.