Winddruk op gebouwen: de windstudie volgens Eurocode 1 met CFD-simulatie

Windstudie van gebouwen: waarom de wind de constructie bepaalt

Op een gebouw is wind geen oppervlakkige hinder: het is een mechanische belasting die van de gebouwomhulling tot in de fundering wordt doorgegeven en die vaak maatgevend is voor slanke, lichte of ruim overspannen bouwwerken.

De windbelasting op gebouwen is in wezen driedimensionaal, instationair en sterk afhankelijk van de omgeving. Zware wind wekt aerodynamische belastingen op die niet alleen de dragende constructie kunnen aantasten, maar ook de gebouwomhulling (vliesgevels, beglazing, gevelbekleding), de technische installaties, de bevestigingssystemen en de buiteninrichting.

Naast de globale krachten zijn het de plaatselijke effecten die de belangrijkste risicofactor vormen. De combinatie van hoogte, geometrie (scherpe hoeken, terugliggende delen, complexe daken), de eventuele doorlaatbaarheid van de gevels en de directe omgeving (obstakels, andere gebouwen, reliëf) brengt geavanceerde aerodynamische verschijnselen voort: overdruk en onderdruk van grote intensiteit, hoge drukgradiënten, plaatselijke versnellingen van de wind, venturi-effecten tussen bouwvolumes, hoekwervels en zones met recirculatie die instabiel zijn.

Deze verschijnselen kunnen leiden tot extreme belastingen op zeer plaatselijke zones, die de globale normatieve methoden vaak slecht weergeven. In dicht bebouwde stedelijke omgevingen versterken de interacties tussen gebouwen deze effecten en maken ze de windbeoordeling bijzonder gevoelig voor de rekenaannames. Een verkeerde inschatting van deze mechanismen kan leiden tot:

Het beheersen van de windeffecten is dus een kwestie van veiligheid, duurzaamheid en kostenbeheersing van het project. Het bereik van onze studies is breed:

Windkrachten op een gebouw: globale stabiliteit en plaatselijk afwaaien

Elke serieuze windstudie is opgebouwd rond één fundamenteel onderscheid, want beide vragen andere methoden, andere coëfficiënten en andere toetsingen. Wij scheiden de globale krachten, die op de hoofddraagconstructie werken en tot in de fundering worden doorgegeven, van de plaatselijke krachten, die maatgevend zijn voor de gebouwomhulling, de aansluitingen en de bevestigingen.

Dit leeskader stuurt de hele verdere studie. De globale krachten worden uitgedrukt als resultante van krachten en momenten die aan de constructeur wordt doorgegeven voor de stabiliteitstoetsingen; de plaatselijke krachten worden uitgedrukt in drukkaarten en extreme waarden per zone, die worden gebruikt voor de berekening van de gebouwomhulling.

In tegenstelling tot benaderingen met gemiddelde coëfficiënten brengt de simulatie fijnmazig in kaart welke over- en onderdruk op alle wanden optreedt, spoort de drukpieken op die met hoek- en randeffecten samenhangen, analyseert de plaatselijke versnellingen van de wind en de afschuifzones, en karakteriseert de aerodynamische interacties tussen naburige gebouwen.

Berekening van de winddruk volgens Eurocode 1 (NEN-EN 1991-1-4)

De studies van EOLIOS blijven strikt binnen het normkader van Eurocode 1, windbelasting (NEN-EN 1991-1-4) en de bijbehorende Franse nationale bijlage. CFD wordt daarbij niet gebruikt als alternatief voor de normen, maar als aanvullend, consistent en verantwoord hulpmiddel. De logica van de norm is een keten van omzettingen: men start bij een meteorologische referentiewind, zet die om in wind op de locatie via ruwheid en orografie, vervolgens in extreme stuwdruk waarin de windvlagen zijn verwerkt, en tot slot in aangebrachte belastingen via de drukcoëfficiënten. Elke schakel heeft zijn parameters:

1. Referentiewindsnelheid: Eurocode-zonering en coëfficiënten

De basissnelheid v_b,0 hoort bij een extreme maar zeldzame windgebeurtenis (gemiddelde herhalingstijd van ongeveer 50 jaar). Zij is gedefinieerd over een periode van 10 minuten, op 10 m boven het maaiveld, in open vlak terrein. Zij hangt af van het geografische gebied en wordt vastgelegd door de zonering van de Eurocode: het Franse vasteland is door de nationale bijlage in 4 windzones verdeeld. Met twee coëfficiënten wordt vervolgens de referentiesnelheid v_b.

Het richtingscoëfficiënt c_dir houdt er rekening mee dat de sterkste wind niet altijd uit de meest ongunstige richting komt: hij maakt het mogelijk de referentiesnelheid te verlagen wanneer de kritische richtingen weinig waarschijnlijk zijn. De seizoenscoëfficiënt c_season geeft de spreiding van extreme wind over het jaar weer: hij is 1 (maximum) voor permanente bouwwerken en kan lager worden gekozen voor tijdelijke bouwwerken zoals steigers.

Vgl. 1

v_b = c_dir · c_season · v_b,0

v_b: referentiesnelheid van het project · c_dir: richtingscoëfficiënt (waarschijnlijkheid van de kritische richtingen) · c_season: seizoenscoëfficiënt (1 voor permanente bouwwerken) · v_b,0: regionale basissnelheid uit de Eurocode-zonering.

2. Terreincategorieën en ruwheid: het windsnelheidsprofiel

De referentiewind is een wind in open vlak terrein; de wind die het bouwwerk werkelijk treft, hangt af van het doorstroomde terrein. De nationale bijlage geeft in tabellen de schijnbare ruwheidshoogte z₀ naar gelang de aard van de locatie, in terreincategorieën 0 tot IV. Bij stedelijke studies hoort de aangehouden ruwheid bij een stedelijke ruwheid van type IV: zij bepaalt het profiel van de windsnelheden, een fundamenteel gegeven voor de inlaatvoorwaarden van het CFD-model.

Het ruwheidscoëfficiënt c_r(z) geeft de gemiddelde verandering van de snelheid met de hoogte weer. Daarin komen z₀, de onderste geldigheidsgrens z_min, de maximale hoogte van het onderzoeksgebied z_max en de terreinfactor k_r:

Vgl. 2

c_r(z) = k_r · ln(z / z₀)

c_r(z): ruwheidscoëfficiënt op hoogte z · k_r: terreinfactor · z₀: ruwheidslengte van de terreincategorie · z: hoogte boven het maaiveld. Geldig voor z_min ≤ z ≤ z_max; onder z_min, c_r(z) = c_r(z_min).

3. Orografie en gemiddelde windsnelheid v_m(z)

Het orografiecoëfficiënt c₀(z) houdt rekening met de invloed van het terreinreliëf (heuvels, ruggen, steilranden) op de windsnelheid: het reliëf versnelt de stroming plaatselijk en verhoogt de wind op de locatie. De effecten mogen worden verwaarloosd wanneer de gemiddelde helling aan loefzijde kleiner is dan 3°; dan geldt c₀(z) = 1. De gemiddelde snelheid v_m(z) combineert de referentiesnelheid met de twee coëfficiënten; door de vergelijking over de hoogte te itereren verkrijgt men het logaritmisch profiel van de gemiddelde wind:

4. Extreme stuwdruk q_p(z): het effect van windvlagen

Natuurlijke wind is niet gelijkmatig: hij fluctueert in vlagen waarvan de relatieve intensiteit toeneemt nabij het maaiveld en op ruwe locaties. De Eurocode vat dit effect met de turbulentie-intensiteit I_v(z) en definieert de extreme stuwdruk q_p(z), waarin het gemiddelde en het fluctuerende deel van de druk zijn samengevoegd. Dit is de spilgrootheid van de hele norm: alle aangebrachte belastingen volgen eruit.

5. Drukcoëfficiënten Cpe en Cpi: van tabellen naar windbelastingen

De druk op een wand volgt uit de vermenigvuldiging van q_p met een drukcoëfficiënt die de vorm van het gebouw en de plaats van de beschouwde zone weergeeft. De Eurocode geeft in zijn tabellen de coëfficiënt uitwendig Cpe per genormeerde zone (zones A tot E op de gevel, F tot I op het dak, waarbij de hoekzones F het zwaarst zijn) en de coëfficiënt inwendig Cpi, die afhangt van de luchtdoorlatendheid van de gebouwomhulling en bij de uitwendige zuiging kan worden opgeteld. De netto rekenbelasting combineert beide:

De grenzen van de Eurocode-berekening: een vrijstaand gebouw, zonder omgeving

De Eurocode leidt zijn coëfficiënten af uit geïdealiseerde geometrieën: vrijstaande balkvormen, eenvoudige zadeldaken, homogene omgeving. Dat kader is robuust, bewezen en voldoet in de meeste gangbare gevallen. Maar het rust op twee sterke aannames: dat de vorm van het project op de tabellen lijkt, en dat de omgeving er niet mee in wisselwerking staat. Zodra het project daarvan afwijkt, worden de forfaitaire waarden ofwel te ruim aan de veilige kant (die voorzichtigheid wordt betaald in staal en bevestigingen), ofwel plaatselijk ontoereikend (het risico verschuift naar het bouwwerk). CFD-simulatie maakt het mogelijk deze grenzen te overstijgen door het volgende mee te nemen:

De norm voorziet die situatie zelf: voor bouwwerken buiten haar geldigheidsgebied verwijst zij naar de windtunnelproef of naar bewezen numerieke methoden. CFD speelt dan de rol van digitale windtunnel: dezelfde windprofielen aan de inlaat, dezelfde grootheden aan de uitgang, met daarbij toegang tot het volledige veld in elk punt van het domein en de vrijheid om varianten te testen zonder een nieuw model te bouwen.

CFD-simulatie van de wind: de methode van EOLIOS

Een simulatie van windbelastingen is niets waard zonder een strikte opbouw. De methodologie van EOLIOS volgt een beheerste keten in vijf stappen: geometrie, windcondities, rekenrooster, berekening voor meerdere richtingen, controle van de convergentie. Elke stap is in het studierapport gedocumenteerd, zodat de keuringsinstantie ze kan controleren.

Atmosferische grenslaag: het juiste windprofiel opleggen

De werkelijke wind is niet uniform: hij is opgebouwd als een atmosferische grenslaag die uit drie sublagen bestaat:

1. 1

   De buitenlaag (inertiële sublaag)

   Een dikte in de orde van een kilometer (0,5 tot 3 km), waar de stroming geleidelijk overgaat in de onverstoorde wind.

2. 2

   De oppervlaktegrenslaag

   Van 10 tot 100 m (ongeveer 10 % van de totale dikte), met sterke snelheids- en temperatuurgradiënten; de windrichting blijft er met de hoogte betrekkelijk constant.

3. 3

   De ruwe sublaag

   Enkele meters boven het maaiveld: driedimensionale en wanordelijke stromingen, sterk beïnvloed door obstakels. Dit is de laag van de voetgangers en de lage geveldelen.

Bij het maaiveld wordt de wind vertraagd door obstakels en ruwheid. Op hoogte, in de onverstoorde lagen van de geostrofische wind (ongeveer 5 km), wordt hij niet meer beïnvloed door de toestand van het aardoppervlak. Daartussen verloopt de snelheid volgens een logaritmisch profiel: dat is de verticale windafschuiving. Het profiel dat aan de inlaat van het domein wordt opgelegd, is het profiel van de maximale vlaagsnelheden, afgeleid uit het gemiddelde profiel en de turbulentie-intensiteit, strikt in lijn met de Eurocode:

Vgl. 6

v_r(z) = √[ 1 + 7 · k_l / ln(z / z₀) ] · v_m(z)

z: hoogte · z₀: ruwheidslengte van de terreincategorie · k_l: turbulentiecoëfficiënt · v_m(z): gemiddelde snelheid. Dit tweede profiel dient als inlaat van het domein voor de simulaties: de stuwdruk die het gebouw ondervindt, is exact de extreme stuwdruk q_p(z) van het normkader.

CFD-rekenrooster: de drukpieken van de wind vatten

Het rekenrooster wordt verfijnd op de plaatsen waar de extremen ontstaan: randen, hoeken, dakranden en dakopstanden, zones van loslating en herhechting. Nabij de wanden vatten prismatische roosterlagen de snelheidsgradiënt; ver van het gebouw wordt het rooster geleidelijk grover om de rekenkosten te beperken. Een gevoeligheidsstudie van het rekenrooster controleert of de drukcoëfficiënten niet meer veranderen bij verdere verfijning: dat is de voorwaarde om de extreme waarden als fysische resultaten te beschouwen en niet als numerieke artefacten.

Acht windrichtingen om de maatgevende belastinggevallen te bepalen

De simulaties worden uitgevoerd voor de acht hoofdrichtingen van de windroos, elk overeenkomend met een extreme vlaagwind volgens Eurocode 1. Deze aanpak met meerdere gevallen is essentieel: de richting die de globale weerstand maximaliseert, is bijna nooit de richting die de hoekzuiging op het dak maximaliseert; scheve wind onder 45° is doorgaans het zwaarst voor afwaaien, zoals bij panelen en zonneparken die aan extreme wind blootstaan, omdat die de conische dakwervels opwekken. De waarden voor de tussenliggende richtingen (NNO, ONO, OZO…) worden geschat met lineaire interpolatie tussen de dichtstbijzijnde gesimuleerde richtingen, en elke zone van de gebouwomhulling wordt uiteindelijk gekarakteriseerd door haar richtingsafhankelijk maximum.

Drukcoëfficiënten, afwaaien: het gebruik van de CFD-resultaten

Stromingsanalyse: begrijpen waar de drukpieken vandaan komen

Naast de resultaatvelden zit de expertise van EOLIOS in de fysische interpretatie van de stromingen. De simulaties maken het mogelijk nauwkeurig de oorsprong van de drukpieken te bepalen, de mechanismen te begrijpen (loslating, wisselwerking wervel/constructie, kanalisering) en optimalisatiemogelijkheden voor te stellen. Deze fijnmazige analyse geldt voor eenvoudige én complexe geometrieën, met een driedimensionale bepaling van de belastingen en de locatie-effecten.

Van drukkaart naar rekenbelastingen

De simulatie geeft in elk punt de snelheids-, druk- en turbulentievelden. De expertise bestaat erin daaruit direct bruikbare grootheden af te leiden voor de constructeur en de gevelbouwer: kaarten van de coëfficiënten Cpe en Cpi per richting en per zone van de gebouwomhulling, vergelijkende richtingsanalyses, plaatselijke extreme waarden voor de bevestigingen, de resultante van krachten en momenten op globaal niveau (drie krachten, drie momenten) voor de stabiliteitstoetsingen. De resultaten worden geleverd in de zone-indeling van het project, niet in die van de software.

Toetsing op afwaaien van daken en gevels

In de zuigzones wordt de maximale plaatselijke onderdruk omgezet in een opwaartse kracht per onderdeel: dakpaneel, baan waterdichting, afzonderlijke bevestiging. De drukcoëfficiënt varieert met de grootte van het oppervlak dat wordt beschouwd: zijn maximale waarde is Cpe,1, de coëfficiënt voor een oppervlak kleiner dan 1 m², waarmee in de praktijk de kleine bevestigingsonderdelen worden berekend. De maximale over 1 m² gemiddelde druk wordt P = q_p · c_pe. Deze belasting wordt vergeleken met de karakteristieke weerstand van de bevestiging en de ondergrond (bepaald door proeven of technische goedkeuring), met de bijbehorende partiële veiligheidsfactoren. Daaruit volgt de vereiste bevestigingsdichtheid, zone per zone:

Kruisvalidatie: CFD-simulatie en analytische Eurocode 1-berekening

Om de uitgevoerde berekeningen te valideren worden de CFD-resultaten systematisch vergeleken met analytische Eurocode-berekeningen op sterk vereenvoudigde geometrieën, die een orde van grootte van de drukken geven. De schatting koppelt de extreme stuwdruk aan de standaard drukcoëfficiënten uit de tabellen. Voor een toren met hoogte h en basis b waarvoor geldt h > 2b houdt de norm q_p(z) = q_p(h) aan voor h − b < z < h; daaruit volgt, voor z > h − b:

Deze dubbele lezing, numeriek en normatief, waarborgt dat geen enkel resultaat wordt gebruikt zonder door het normkader te zijn begrensd. De deliverables van de studie zijn verdeeld over de constructie en de gebouwomhulling:

Een CFD-windstudie ten dienste van beheerste projecten

Voor wie is de windstudie bedoeld?

De windstudie van EOLIOS voegt zich in de ontwerpketen van het project en spreekt elke betrokkene in zijn eigen taal aan:

Waarom Eurocode 1 aanvullen met een CFD-windstudie?

Eurocode 1 blijft de regelgevende basis: hij legt de referentiesnelheid, het windprofiel en de veiligheidsfactoren vast en waarborgt de conformiteit. CFD vervangt hem niet, maar verlengt hem daar waar zijn geïdealiseerde geometrieën hun grenzen bereiken: waar de forfaitaire berekening een vrijstaand gebouw beschouwt, plaatst de simulatie het project terug in de werkelijke omgeving van de locatie en geeft zij de werkelijke verdeling van de krachten weer, inclusief omgeving, windschaduw en interacties.

Door grondige beheersing van de normeisen, geavanceerde CFD-expertise en vermogen om de resultaten fysisch te interpreteren te combineren begeleidt EOLIOS haar klanten bij het veiligstellen en optimaliseren van hun projecten tegenover windrisico's: zuigzones veiligstellen voordat ze schade worden, constructie en bevestigingen zo scherp mogelijk dimensioneren in plaats van zo ruim mogelijk, en elke aanname verantwoorden. CFD wordt zo een strategisch instrument voor risicobeperking, technische verantwoording en kostenbeheersing, ten dienste van veilige, duurzame en goed presterende gebouwen. Een windstudie uitgevoerd binnen het kader van Eurocode 1, van de berekening van de extreme stuwdruk tot de toetsing op afwaaien.