7904x
001696
3.3.2021

Profil rychlosti větru a intenzity turbulence pro stanovení kvazistatického zatížení větrem v souladu s koncepcí poryvu větru

Odezva konstrukcí na působení větru se liší v závislosti na jejich tuhosti, hmotnosti a na tlumení. V zásadě rozlišujeme mezi konstrukcemi náchylnými a nenáchylnými ke kmitání.

Obvykle se konstrukce nepovažují za náchylné ke kmitání, pokud se deformace při působení větru při poryvové rezonanci nezvýší o více než 10 % [2]. V tomto případě lze časově proměnné působení větru popsat jako statické ekvivalentní zatížení.

Za předpokladu, že turbulence v proudění větru jsou vzhledem k rozměrům budovy velmi velké, lze v programu RWIND Simulation vypočítat staticky působící rozdělení tlaku p na geometrii budovy podle „kvazistacionární metody“ nebo takzvaného poryvu větru. koncept [3].

V zásadě se předpokládá stacionární pole proudění okolo analytického modelu pro turbulentní kolísání rychlosti po dobu trvání poryvu [3]. Kolísání tlaku na povrchu modelu vlivem turbulence větru se tak považuje za stav, který setrvává po určitou dobu t. Fluktuace tak přesně sledují průběh časově středovaných součinitelů tlaku cp,mean na povrchu modelu.

Výsledný tlak na povrch modelu vyvolaný větrem Δp(t) pak závisí čistě na rychlosti proudění v(t).

Hodnota vektoru rychlosti proudění v(t) je tedy:
v(t)² = (vx,mean + vx,fluctuation(t))² + vy,fluctuation(t)² + vz,fluctuation(t)²

Pokud je příspěvek kvadratických členů nepatrný, je efektivní hodnota vektoru rychlosti proudění v(t):
v(t)² = vx,mean² + 2 ⋅ vx,mean ⋅ vx,fluctuation(t)

Po dosazení efektivní rychlosti proudění do rovnice pro tlak vyvolaný větrem dostaneme:
Δp(t) = 1/2 ⋅ ρ ⋅ vx,mean² [1 + (2 ⋅ vx,fluctuation(t)) / vx,mean] ⋅ cp,mean

Tato transformace ukazuje, že kolísání tlaku větru Δp(t) závisí pouze na kolísání rychlosti větru vx,fluctuation(t) v hlavním směru proudění x.

Pokud nahradíme časově proměnné kolísání rychlosti vx,fluctuation(t) maximálním kolísáním rychlosti vx,fluctuation,max, odstraníme ze systému časovou proměnnost.

Jestliže poté porovnáme výraz vx,fluctuation,max / vx,mean jako násobek g intenzity turbulence Iv(z),


lze tento výraz v hranatých závorkách popsat jako součinitel nárazu větru G(z). Při dosazení výrazů do rovnice pro jmenovité zatížení větrem dostaneme:

kde:

Například v EN 1991-1-4 je součinitel g udávající trvání poryvu větru 3,5.

Program RWIND Simulation stanoví střední hodnoty tlaků pmean na povrchu modelu v závislosti na rychlosti proudění vx(z) stacionárním řešením RANS rovnic pomocí algoritmu SIMPLEC. Protože střední hodnoty součinitelů tlaku cp,mean vycházejí z poměru mezi stanovenými středními hodnotami tlaku pmean a maximálním dynamickým tlakem nerušeného větru ve výšce střechy q (výška střechy),
cp,mean = pmean / q(výška střechy)
lze pro stanovení jmenovitých zatížení větrem podle poryvové koncepce použít rychlost proudění z přepočteného maximálního dynamického tlaku větru q(z) podle výšky [1].
v(z) = √(2 ⋅ q(z) / ρ)

Tato rychlost větru tak zahrnuje střední rychlost větru vmean a maximální fluktuační složku vfluctuation. Intenzitu turbulence proudění lze přitom nastavit konstantně po výšce na velmi malou hodnotu asi 5 % [4].

Při zohlednění účinků sil působících na celou budovu nebo na velké plochy poskytuje tato metoda velmi dobrou aproximaci k přirozenému zatížení větrem [3]. Důvodem je, že průměrováním potlačené malé turbulentní účinky působí pouze v dílčích oblastech a při globální integraci sil nemají znatelný efekt.

Kromě toho tato koncepce velmi dobře reaguje i na malé dílčí oblasti s čelním prouděním, protože zde jsou již velmi dobře zaznamenány efektivní kolísání tlaku v profilu maximální rychlosti větru [3].

Naopak u ploch s oddělenými toky větru (boční a zadní stěna) vykazuje systém horší přiblížení realitě. Zejména v těchto zónách má turbulence indukovaná stavbou a „odstraněná“ průměrováním podle poryvové koncepce větší účinek, než je turbulence proudění zahrnutá v profilu rychlosti proudění.


Autor

Ing. Niemeier je zodpovědný za vývoj hlavních programů RFEM, RSTAB, RWIND a oblast membránových konstrukcí. Zároveň má na starosti řízení jakosti a podporu zákazníkům.

Odkazy
Reference
  1. Eurokód 1: Zatížení konstrukcí - Část 1-4: Obecné účinky, Zatížení větrem; BS EN 1991-1-4:2010-12
  2. Albert, A. Schneider - Bautabellen für Ingenieure mit Berechnungshinweisen und Beispielen, 22. vydání. Bochum: Bundesanzeiger, 2016
  3. Kiefer, H. Windlasten an quaderförmigen Gebäuden in bebauten Gebieten, 2003
  4. Werth, M. Vergleichende Studie zu Windlastmodellen im Hochbau: Numerische Strömungsberechnung vs. Druckmessungen im Windkanal, 2019


;