Windturbineblad: analyse van vloeistof–structuurinteractie
Het doel was te zien hoe een composiet windturbineblad van 41 m doorbuigt onder wind en rotatie: de tip buigt 400,9 mm door onder de gecombineerde belasting, en een handcontrole van de kracht in de bladwortel klopt tot op 0,1%. Opgebouwd uit 20 profieldoorsneden in SolidWorks, doorgerekend in ANSYS Fluent en daarna vertaald naar ANSYS Mechanical (eenzijdig gekoppelde vloeistof–structuurinteractie).
- Project
- Casestudy, Cornell University (onlinecursus)
- Mijn rol
- CAD, CFD, FEA en verificatie
- Software
- SolidWorks, ANSYS Fluent, ANSYS Mechanical
- Methode
- Eenzijdig gekoppelde vloeistof–structuurinteractie (CFD naar FEA)
- 41,25maerodynamische spanwijdte van het blad, uit 20 profieldoorsneden
- 7,7Melementen in de fijnste van 5 CFD-meshes
- 400,9mmmaximale tipdoorbuiging onder gecombineerde belasting
- ≤0,1%grootste afwijking van tipsnelheid en reactiekracht in de wortel ten opzichte van de handberekening
Aangetoonde vaardigheden
Aanpak
De studie volgt dezelfde keten die een bladontwerp doorloopt vóór de constructieve vrijgave: geometrie en materiaalsysteem vastleggen, de aerodynamische belasting bepalen, die belasting overbrengen op een sterktemodel, en elke stap toetsen aan een onafhankelijk getal in plaats van alleen op de solver te vertrouwen.
-
Het blad modelleren
SolidWorks
Het blad is gedimensioneerd naar de algemene klasse van een rotor van 1,5 MW van het type GE 1.5xLE: een aerodynamische spanwijdte van 41,25 m, opgebouwd uit 20 stations over de spanwijdte met drie NREL-profielfamilies. Elke profieldoorsnede is op een eigen vlak geschetst, om de voorrand gedraaid tot de plaatselijke verdraaiingshoek en verbonden met een afzonderlijke loft van de boven- en onderzijde. Zo blijven de lijnen van voor- en achterrand eenduidig voor de CFD- en FEA-meshers die erop volgen.
Huid en ligger zijn gemodelleerd als dunne referentieschalen in plaats van volumelichamen, verdeeld in vijf secties over de spanwijdte zodat de aerodynamische belasting goed langs het blad kan variëren, en samengevoegd tot één verlijmd onderdeel uit meerdere lichamen, passend bij een gezamenlijk uitgeharde composietconstructie.
- Aerodynamische spanwijdte
- 41,25 m
- Verdraaiing, van wortel tot tip
- 42° → 0°
- Huid
- Laminaat van E-glas/epoxy
- Ligger
- Laminaat van koolstofvezel/epoxy
Huid en ligger in de doorsnede bij de wortel. Eén ligger op een kwart van de koorde is een gangbare vereenvoudiging in het ontwerp van composietbladen. -
De luchtstroming simuleren
ANSYS Fluent
Eén periodieke sector van 120° met één blad is doorgerekend in een roterend referentiestelsel bij een instroming onder nominale omstandigheden (12 m/s, ω = 2,22 rad/s, snellopendheid ≈ 8), met het turbulentiemodel SST k-ω: de standaardkeuze voor het voorspellen van grenslaagloslating op een liftgevend oppervlak.
De drukcontour vertoont een duidelijke asymmetrie tussen beide zijden van het blad. Het donkerblauwe gebied aan de zuigzijde bereikt een aanzienlijk grotere waarde dan de rode en groene gebieden, zodat die zuigpiek het netto drukverschil bepaalt dat lift opwekt.
- Instroming
- 12 m/s
- Hoeksnelheid
- 2,22 rad/s
- Tipsnelheid
- ≈98 m/s
- Basismesh
- 356.628 elementen
Snelheid ten opzichte van het blad, doorgetrokken naar alle drie bladen: de snelheid neemt toe met de straal en bereikt ≈98 m/s nabij de tip. 
Statische druk rond een representatieve doorsnede: de zuigpiek die lift opwekt. 
Inflatielagen in de grenslaag en meshverfijning in het nabije veld bij de bladtip. -
De windbelasting overbrengen
Van ANSYS Fluent naar ANSYS Mechanical
Het geconvergeerde CFD-drukveld is knoop voor knoop overgebracht op een afzonderlijk gemesht sterktemodel van schaal en ligger, met 100% getrouwheid van de overdracht (eenzijdige FSI).
- Overgebrachte knopen
- 100%
- Drukbereik
- +2.897 tot −5.466 Pa
Aerodynamisch drukveld uit ANSYS Fluent geïmporteerd op de FEA-schaalmesh; 100% van de doelknopen kreeg een overgebrachte waarde. -
Doorbuiging en spanning berekenen
ANSYS Mechanical
Het model is doorgerekend onder gecombineerde aerodynamische druk en centrifugale belasting, met spanningsverstijving door grote vervormingen ingeschakeld, om vast te leggen hoe de centrifugale trekkracht in een draaiend blad de effectieve buigstijfheid vergroot. Een lokaal coördinatenstelsel legt in elk element de orthotrope vezelrichting langs de spanwijdte-as van het blad.
- Tipdoorbuiging
- 400,9 mm
- Maximale spanning
- 27,8 MPa
- Reactiekracht in de wortel
- 1.578,1 kN
- Schaalelementen
- 4.831
De constructieve belasting komt uit de basismesh met 356.628 elementen. De vervorming en spanning laten dus zien dat de FSI-werkwijze correct is uitgevoerd, maar vormen geen volledig mesh-geconvergeerde ontwerpbelasting. Zie de controles.
Animatie van de vervorming: het blad beweegt van onvervormd naar de maximaal doorgebogen vorm onder belasting. 
Equivalente (Von Mises-)spanning: maximaal 27,8 MPa, bij de overgang van de wortel naar het buitendeel. Selecteer om te vergroten.
Controle van de resultaten
Elk hoofdresultaat in deze studie is getoetst aan een onafhankelijke handberekening, een trend uit de meshgevoeligheid of externe referentiegegevens, en wordt niet als kale solveruitvoer gerapporteerd.
≤0,1%
grootste afwijking van de tipsnelheid (CFD) en de reactiekracht in de wortel (FEA) ten opzichte van onafhankelijke handberekeningen
De tipsnelheid van het blad en de reactiekracht in de wortel komen tot op 0,1% overeen met onafhankelijke handberekeningen, en de massabalans van de CFD sluit tot op ongeveer een miljoenste van de doorstroming. Dat bevestigt sterk dat randvoorwaarden, massa-eigenschappen en belastingsoverdracht consistent zijn gedefinieerd en toegepast. Het belangrijkste open punt is de meshonafhankelijkheid van de CFD: de vermogenscoëfficiënt van de rotor is bij de fijnste geteste mesh nog niet gestabiliseerd (zie de grafiek hieronder), en dat wordt expliciet meegenomen in plaats van verbloemd.
| Controle | Resultaat | Criterium | Status |
|---|---|---|---|
| Massabalans van de CFD | −0,755 kg/s (≈0,0001% van de doorstroming) | ≈ 0 | Voldoet |
| Convergentie van de residuen in de CFD | Alle residuen ≤ 1×10−6 | ≤ 1×10−6 | Voldoet |
| Tipsnelheid van het blad (CFD tegenover handberekening) | 98,05 tegenover 98,12 m/s | Goede overeenstemming | Voldoet (0,07%) |
| Reactiekracht in de wortel (FEA tegenover handberekening) | 1.578,1 tegenover 1.576,3 kN | Goede overeenstemming | Voldoet (0,1%) |
| Overdracht van de belasting CFD → FEA | 100% van de knopen overgebracht | 100% | Voldoet |
| Vermogenscoëfficiënt tegenover de Betz-limiet | 13–42% (meshafhankelijk) | < 59,3% | Voldoet, alle meshes |
| Vermogenscoëfficiënt tegenover de referentieturbine | 41,5% (fijnste mesh) tegenover streefwaarde 30% | Overeenstemming in orde van grootte | Open, niet mesh-geconvergeerd |
| Meshonafhankelijkheid van de CFD | Niet bereikt bij 7,7M elementen | Stabilisatie van Cp | Open |
Meshverfijning: vermogenscoëfficiënt tegenover meshdichtheid
Vermogenscoëfficiënt van de rotor (Cp), opnieuw berekend bij vijf meshdichtheden. Beweeg over een punt voor details.
De vermogenscoëfficiënt verdrievoudigt bijna tussen de grofste en de fijnste geteste mesh en is bij 7,7 miljoen elementen nog niet gestabiliseerd. Hij blijft wel steeds ruim onder de Betz-limiet en sluit de streefwaarde van 30% in, die onafhankelijk is berekend uit het nominale vermogen van de echte turbine.
De constructieve belasting in stap 4 komt uit de basismesh met 356.628 elementen. De vervorming en spanning die volgen laten dus zien dat de FSI-werkwijze correct is uitgevoerd, maar vormen geen volledig mesh-geconvergeerde ontwerpbelasting; dat gat dichten is het belangrijkste vervolgpunt in de aanbevelingen hieronder.
Technisch oordeel en aanbevelingen
Naast de controles voldoet/voldoet niet laat de studie zien waar het model solide is en waar het nog werk nodig heeft voordat het een productiebeslissing kan dragen.
-
Aandachtspunt
De meshafhankelijkheid van de CFD wegwerken
De vermogenscoëfficiënt van de rotor verdrievoudigt bijna tussen de grofste en de fijnste geteste mesh (13,4% → 41,5%) zonder duidelijk te stabiliseren. Omdat de constructieve belasting uit de grovere basismesh komt, laten de gerapporteerde vervorming en spanning een correct uitgevoerde FSI-werkwijze zien, maar geen definitieve, mesh-geconvergeerde ontwerpbelasting. Herhaal de sterkteberekening zodra de CFD meshonafhankelijk is.
-
Spanningsverstijving doet ertoe bij een draaiend blad
Door Large Deflection in te schakelen wordt meegenomen hoe de centrifugale trekkracht de effectieve buigstijfheid van het blad verhoogt: hetzelfde mechanisme waardoor de eigenfrequentie van een draaiend turbine- of helikopterblad boven die in stilstand komt te liggen. Zonder deze optie wordt de stijfheid onderschat en de doorbuiging overschat.
-
De constructieve marge is ruim, maar slechts indicatief
Een indicatieve veiligheidsfactor van ongeveer 20 ten opzichte van een representatieve toelaatbare waarde voor composiet weerspiegelt de bewust conservatieve schaaldiktes in dit model, die ook de centrifugale belasting van het blad vergroten. Een productieontwerp vraagt om een laminaatcontrole per laag en een vermoeiingsbeoordeling voordat die marge iets zegt voor certificering.
-
Kans: een op massa geoptimaliseerd dikteverloop
Het huidige lineaire dikteverloop levert een blad op dat meerdere keren zwaarder is dan een productieblad met deze spanwijdte. Omdat de centrifugale belasting rechtstreeks met de massa schaalt, zou een taps, laminaatgeoptimaliseerd dikteverloop massa en belasting samen verlagen: een logische volgende iteratie voor dit model.
Technische details
Volledige instellingen, materiaalgegevens en verificatie, voor wie zich verder wil verdiepen.
Opzet van CFD en FEA Vergelijkingen, domein, mesh, turbulentiemodel, randvoorwaarden, elementen
Aerodynamische analyse (ANSYS Fluent)
Eén periodieke sector van 120° met één blad is doorgerekend in een roterend referentiestelsel bij een instroming onder nominale omstandigheden (12 m/s, ω = 2,22 rad/s, snellopendheid ≈ 8), met het turbulentiemodel SST k-ω: de standaardkeuze voor het voorspellen van grenslaagloslating op een liftgevend oppervlak.
∇·(ρvrvr) + ρ(2ω×vr + ω×ω×r) = −∇p + ∇·τr
- Domein: periodieke sector van 120° met één blad, radiaal ongeveer 3–5× de bladstraal en axiaal 6× de spanwijdte.
- Mesh: ongestructureerde tetraëders, 5 inflatielagen aan de wand, elementgrootte 0,30 m op het bladoppervlak, 356.628 elementen in de basismesh.
- Turbulentie: SST k-ω, productiebegrenzer ingeschakeld, automatische wandbehandeling.
- Solver: gekoppeld druk-snelheidsschema, tweede-orde-upwind, residuen tot 1×10−6.
- Randvoorwaarden: snelheidsinlaat van 12 m/s met 5% turbulentie-intensiteit, drukuitlaat, no-slipwand op het blad, rotatieperiodieke zijvlakken.
De drukcontour vertoont ook een duidelijke asymmetrie tussen beide zijden van het blad. Op de kleurschaal zijn de rode (hoge positieve druk) en groene (matig negatieve druk) gebieden ongeveer even groot, maar het donkerblauwe gebied aan de zuigzijde bereikt een aanzienlijk grotere waarde dan beide. Omdat die zuigpiek zwaarder weegt dan de positieve druk aan de windzijde, bepaalt ze het netto drukverschil dat lift opwekt.
Sterkteanalyse (ANSYS Mechanical)
Het geconvergeerde CFD-drukveld is knoop voor knoop overgebracht op een afzonderlijk gemesht sterktemodel van schaal en ligger (100% getrouwheid van de overdracht) en daarna doorgerekend onder gecombineerde aerodynamische druk en centrifugale belasting, met spanningsverstijving door grote vervormingen ingeschakeld om vast te leggen hoe de centrifugale trekkracht in een draaiend blad de effectieve buigstijfheid vergroot.
[K]{d} = {F}
- Elementen: schaalelementen met 6 vrijheidsgraden (type SHELL181/281), 4.831 elementen en 4.671 knopen.
- Randvoorwaarde in de wortel: alle wortelknopen star gekoppeld aan één extern punt in de naaf, dat alle zes reactiecomponenten levert en niet alleen krachten.
- Belasting: geïmporteerd CFD-drukveld plus de centrifugale volumekracht bij een hoeksnelheid van 2,22 rad/s, met Large Deflection (spanningsverstijving) ingeschakeld.
- Materiaaloriëntatie: een lokaal coördinatenstelsel legt in elk element de orthotrope vezelrichting langs de spanwijdte-as van het blad, een noodzakelijke stap voor een fysisch zinvol anisotroop resultaat.
Materiaalmodel van het composiet Gehomogeniseerde orthotrope eigenschappen in het FE-model
De buitenhuid is gemodelleerd als een laminaat van E-glas/epoxy en de inwendige ligger als een laminaat van koolstofvezel/epoxy: de gangbare verdeling in de industrie, met een kosteneffectieve glashuid die de aerodynamische vorm en de afschuiving draagt en koolstof geconcentreerd in de liggerflens, waar het voordeel in stijfheid per gewicht het meest telt voor buiging in klaprichting (flapwise). Beide zijn in het FE-model weergegeven door één gehomogeniseerd orthotroop materiaal.
| Eigenschap | Waarde | Eenheid |
|---|---|---|
| Dichtheid | 1.550 | kg/m³ |
| Elasticiteitsmodulus, X (spanwijdte / vezel) | 113,75 | GPa |
| Elasticiteitsmodulus, Y (koorderichting) | 7,583 | GPa |
| Elasticiteitsmodulus, Z (dikterichting) | 7,583 | GPa |
| Dwarscontractiecoëfficiënt, XY | 0,32 | — |
| Dwarscontractiecoëfficiënt, YZ | 0,37 | — |
| Dwarscontractiecoëfficiënt, XZ | 0,35 | — |
| Glijdingsmodulus, XY | 5,446 | GPa |
| Glijdingsmodulus, YZ / XZ | 2,964 | GPa |
Vergeleken met handboekgegevens voor laminaten van E-glas en koolstof/epoxy vallen deze gehomogeniseerde eigenschappen binnen de gangbare bereiken, en de resulterende stijfheidsverhouding tussen spanwijdte- en koorderichting (≈15:1) past bij een door de ligger gedomineerde doorsnede met veel koolstof.
Details van de analyse Elke verificatiecontrole uitgewerkt: wat ze vaststelt en wat ze niet bewijst
De sectie “Controle van de resultaten” vat elke verificatiecontrole in één regel samen. Hier wordt elke controle afzonderlijk uitgewerkt: wat ze vaststelt, de maatgevende vergelijking, en wat ze wel en niet bewijst.
V1. Convergentie en behoudswetten op solverniveau
Bevestigt dat de gediscretiseerde CFD-vergelijkingen correct zijn opgelost, ongeacht of het onderliggende model met de werkelijkheid overeenkomt. Drie fysisch verschillende indicatoren worden samen gecontroleerd in plaats van één afzonderlijk: de globale massabalans, de geschaalde residuen en de over het oppervlak geïntegreerde druk op het blad, de grootheid die de overgebrachte constructieve belasting werkelijk bepaalt.
De netto massaonbalans over de domeingrenzen is −0,755 kg/s bij een doorstroming van ruim 660.000 kg/s; de residuen voor continuïteit, impuls en turbulentietransport bereiken alle de streefwaarde; en de monitor van de oppervlaktedruk blijft over de volledige berekening van 1.500 iteraties vlak, in plaats van alleen een residugrens te halen terwijl de geïntegreerde belasting nog verloopt.
V2. Kwaliteit van de CFD-mesh
Een geconvergeerde oplossing zegt alleen iets als de onderliggende mesh goed gevormd is. De scheefheid wordt getoetst aan de algemene kwaliteitsklassen uit de ANSYS Meshing User’s Guide; dat kan een residuplot alleen niet bevestigen.
Het overgrote deel van de basismesh met 356.628 elementen valt in de scheefheidsklassen “uitstekend” (0–0,25) tot “zeer goed” (0,25–0,50), zonder elementen in de problematische bereiken: een mesh die goed geschikt is voor de SST k-ω-oplossing, zonder plaatselijke numerieke kwaliteitsproblemen die anders de druk- en snelheidsgradiënten voor de constructieve belastingsoverdracht zouden kunnen verstoren.
V3. Controle van de kinematica in het roterende stelsel
Het roterende referentiestelsel en de randvoorwaarde voor de hoeksnelheid worden getoetst aan elementaire kinematica van een star lichaam: voor een punt op straal r dat met ω roteert, is de tangentiële snelheid v = rω, onafhankelijk van wat de CFD-solver zelf berekent.
vtip = rω = (44,2 m)(2,22 rad/s) = 98,12 m/s
De door de CFD berekende maximale snelheid ten opzichte van het blad komt tot op 0,07% overeen met de handberekening: een sterke bevestiging dat de MRF-rotatierandvoorwaarde (ω = −2,22 rad/s om +z) correct is geïmplementeerd. De snelheidsdriehoek toont dezelfde geometrie van de relatieve wind die deze controle valideert: de axiale instroming (U) en de tangentiële bladsnelheid (ΩR) vormen samen de relatieve wind die de plaatselijke invalshoek bepaalt.
V4. Toereikendheid van de domeingrootte in de CFD
Een formele studie naar domeinonafhankelijkheid (blokkering), waarbij op een groter domein opnieuw wordt gerekend en wordt bevestigd dat het resultaat gelijk blijft, is niet uitgevoerd. In plaats daarvan is de gemodelleerde vrije ruimte rond het blad kwalitatief getoetst aan algemene richtlijnen voor CFD van externe aerodynamica.
Het domein reikt radiaal tot ongeveer 120–240 m (circa 1,4–2,7 rotordiameters, bij een tipstraal van 44,25 m) en axiaal tot ongeveer 270 m (circa 3 rotordiameters). Dat is ruim minder dan de 5–10 rotordiameters die doorgaans worden aanbevolen voor de externe aerodynamica van rotoren, zodat blokkering de berekende belastingen en vermogenscoëfficiënt enigszins kan opdrijven. Een specifieke studie naar domeinonafhankelijkheid blijft open als vervolgwerk.
V5. Getrouwheid van de belastingsoverdracht bij eenzijdige FSI
De geconvergeerde oppervlaktedruk uit de CFD moet volledig en consistent worden overgebracht op de afzonderlijk gegenereerde FEA-schaalmesh, zonder niet-overgebrachte (geëxtrapoleerde of standaard ingevulde) doelknopen, wil men de twee meshes geometrisch consistent kunnen noemen.
De profielbehoudende overdracht van ANSYS meldt 100% getrouwheid: elke FEA-doelknoop kreeg een geldige, geïnterpoleerde drukwaarde uit de CFD-bronmesh. Het overgebrachte bereik komt overeen met de lage druk aan de zuigzijde en de hogere druk aan de drukzijde die in de aerodynamische resultaten kwalitatief zijn beschreven, en weerspiegelt rechtstreeks het netto liftgevende drukverschil tussen beide zijden van het blad. Het geïmporteerde veld is weergegeven in stap 3.
V6. Kwaliteit van de FEA-mesh en controle van de reactiekracht
Net als aan de CFD-kant wordt de kwaliteit van de sterktemesh op scheefheid gecontroleerd. Daarnaast wordt de reactiekracht in de wortel getoetst aan een handberekening voor een puntmassa op vaste straal onder centrifugale belasting: een controle die alleen van de kinematica van een star lichaam afhangt en niet van de FEA-solver zelf.
Fr = mrω² = (22.473 kg)(14,232 m)(2,22 rad/s)² = 1.576,3 kN
De FEA-schaalmesh met 4.831 elementen heeft geen elementen in het problematische scheefheidsbereik. Daarnaast komt de door ANSYS berekende reactiekracht in de wortel tot op 0,1% overeen met de handberekening: een sterke bevestiging dat de bladmassa, de ligging van het zwaartepunt, de hoeksnelheid en de randvoorwaarde met extern punt in de wortel correct zijn gedefinieerd en consistent zijn toegepast.
