168x
006036
22.7.2025

H.1.5.1. Skupina 1: Součinitel síly a momentu

Příběh uživatele

V tomto příkladu spočítáme průměrný součinitel síly a součinitel momentu pro Tokijská polytechnická univerzita (TPU) jako experimentální příklad (M0/S0), například ty, které se vztahují na proces statických výpočtů na základě WTG-Merkblatt-M3.

Tato část výsledků patří do „skupiny 1“ podle „obrázku 2.2“ ve WTG-Merkblatt-M3:

  • G2: Absolutní hodnoty se středními požadavky na přesnost. Oblast použití může zahrnovat parametry nebo předběžné studie, pokud jsou plánovány pozdější studie s vyšší přesností (např. zkouška ve větrném tunelu třídy G3).
  • R2: Samostatné, všechny relevantní směry větru s dostatečně jemným rozlišením směru.
  • Z1: Statistické střední hodnoty, pokud se týkají stacionárních procesů proudění, kde lze kolísání (např. v důsledku přibližující se turbulence proudění) dostatečně zachytit jinými opatřeními.
  • S1: Statické účinky. Stačí znázornit statický model s nezbytnými mechanickými detaily, ale bez hmoty a tlumicích vlastností.

Popis

Toto srovnání, které zkoumá aerodynamické chování v různých rovinách, má za cíl poskytnout komplexnější pochopení rozložení sil. Tyto poznatky jsou zásadní pro hodnocení celkové stability modelu, aerodynamického výkonu a odezvy konstrukce za podmínek proudění. Průměrné profily rychlosti a intenzity turbulence proudění reprodukované ve větrném tunelu odpovídají profilům terénní kategorie IV podle normy Japonského architektonického institutu (AIJ). Exponent 0,25 uvedený v některých číslech byl po dotazu klasifikován jako typografická chyba. Předpokládaný model je znázorněn na obrázku 1:

Předpoklady, z nichž vychází analýza a simulace, jsou shrnuty v tabulce 1, která nabízí jasný přehled parametrů a podmínek zohledněných v celé studii.

Tabulka 1: Vstupní data 3D obdélníkového modelu

Parametr Symbol Hodnota Jednotka
Referenční rychlost větru UH 11 m/s
Výška střechy Href 0,4 m
Exponent profilu α 0,25 -
Kategorie terénu - IV -
Hustota vzduchu – RWIND ρ 1,25 kg/m³
Turbulentní model – RWIND RANS K-Omega - -
Kinematická viskozita – RWIND ν 1,5×10⁻⁵ m²/s
Pořadí schématu – RWIND Druhé - -
Reziduální cílová hodnota – RWIND 10⁻⁴ - -
Reziduální typ – RWIND Tlak - -
Minimální počet iterací – RWIND 800 - -
Mezní vrstva – RWIND NL 10 -
Typ stěnové funkce – RWIND Vylepšená / Smíšená - -

Studie výpočetní sítě

Obrázek 2 shrnuje studii citlivosti sítě a ukazuje, že s nárůstem hustoty sítě z 20 % na 55 % se součinitel síly (Cf) zvyšuje z 0,95 na 1,05 a poté se stabilizuje. To naznačuje, že simulace dosahuje nezávislosti na síti při 55 %, což znamená, že další zahuštění sítě již významně nemění výsledky. Obrázek zdůrazňuje význam kvality sítě pro zajištění spolehlivých výsledků simulace větru.

Studii vytvoření sítě je nutné provést podle následujícího odkazu:

Obrázek 3 představuje diagram, který porovnává výsledky RWIND s experimentálními daty, pokud jde o vykreslení normalizovaných rychlostních profilů podle výšky. Výsledky vykazují velmi dobrou shodu.

Požadavky na přesnost WTG-Merkblatt

WTG-Merkblatt M3 poskytuje dvě klíčové metody pro ověření výsledků simulace. Metoda „Hit Rate“ vyhodnocuje, kolik simulovaných hodnot Pi správně odpovídá referenčním hodnotám Oi v rámci definované tolerance, a to pomocí binární klasifikace (hit nebo miss). Tento přístup hodnotí spolehlivost simulace výpočtem míry úspěšnosti q, podobně jako funkce spolehlivosti používané v teorii spolehlivosti. Naopak metoda „normalizované střední kvadratické chyby“ (e2) nabízí podrobnější hodnocení přesnosti kvantifikací průměrné kvadratické odchylky mezi simulovanými a referenčními hodnotami, normalizované tak, aby zohledňovala rozdíly ve stupnici. Společně tyto metody poskytují jak kvalitativní, tak kvantitativní měřítka pro ověření simulace.

Výsledky a diskuse

Tabulka 3 ukazuje velmi dobrou shodu mezi výsledky RWIND a referenčními údaji WTG pro normalizované hodnoty rychlosti. Všechny odchylky jsou v přijatelném rozmezí (pod 10 %), takže je dosaženo „míry úspěšnosti“ „q = 100 %“ a „normalizovaná střední kvadratická chyba“ „e2 = 0,00001“ je extrémně nízká. Tyto výsledky potvrzují, že program RWIND přesně reprodukuje referenční profil větru a splňuje přísná validační kritéria.

Tabulka 3: Metrika validace pro rychlostní profil

WTG - u / uref RWIND - u / uref Odchylka (%) n
0,394 0,390 1,045 1,00
0,478 0,473 1,093 1,00
0,566 0,564 0,291 1,00
0,629 0,629 0,081 1,00
0,675 0,674 0,183 1,00
0,713 0,712 0,072 1,00
0,750 0,753 0,422 1,00
0,784 0,783 0,098 1,00
0,822 0,819 0,362 1,00
0,869 0,865 0,463 1,00
0,897 0,901 0,496 1,00
0,939 0,940 0,113 1,00
1,010 1,005 0,493 1,00
1,065 1,070 0,540 1,00
1,123 1,124 0,044 1,00
1,161 1,162 0,121 1,00
1,195 1,195 0,051 1,00
1,237 1,240 0,263 1,00
1,266 1,270 0,291 1,00
1,299 1,302 0,223 1,00

Obrázek 4 ilustruje, jak jsou v programu RWIND definovány zóny plochy pro přesný výpočet součinitelů síly na každé ploše modelu budovy vystavené působení větru. 3D model s jednoduchým obdélníkovým tvarem má každé ploše (přední, zadní, levé, pravé a horní) přiřazenu samostatnou zónu označenou barevným kódem. Tyto zóny umožňují programu RWIND samostatně analyzovat a vypočítat aerodynamické síly a rozložení tlaku nezávisle pro každý povrch. Místo umístění četných bodových sond jsme se rozhodli definovat povrchové zóny, protože program RWIND může přímo určit součinitel síly pro každou definovanou zónu.

Ve středu obrázku 4 je otevřeno okno „Upravit zónu“ pro zónu č. 2 – přední jako příklad. Tento podrobný výstup obsahuje klíčové aerodynamické parametry. Promítnutá plocha této zóny ve směru větru je 0,04 m². Vypočtená odporová síla na této ploše je 1,7 N a odpovídající součinitel síly (Cx) je 0,56. V následujícím vzorci je jako příklad vypočítán součinitel síly větru (aerodynamický součinitel síly) pro konkrétní tlakovou zónu. Tuto hodnotu lze poté porovnat s hodnotou zobrazenou na záložce Info v dialogovém okně Upravit data modelu.

Obrázek 5 zobrazuje výsledky následného zpracování simulace větru provedené pro vysokopodlažní budovu pomocí modelu turbulence k-omega. Na pravé straně obrázku je srovnávací tabulka, která obsahuje koeficienty síly získané z CFD simulace spolu s referenčními hodnotami z pokynů WTG. Přední plocha má součinitel síly 0,562, zatímco referenční hodnota WTG je 0,540, což odpovídá odchylce 4,07 %. Zadní plocha, která je vystavená podtlaku v důsledku jevů v úplavu, vykazuje mírně větší odchylku 6,09 %. Boční plochy také vykazují malé rozdíly s odchylkami v rozmezí 3,68 % až 5,38 %. Je pozoruhodné, že celkový součinitel síly ve směru větru, který představuje globální odporový účinek, vykazuje pouze odchylku 1,95 %, což podtrhuje přesnost a spolehlivost CFD simulace při zachycení celkového aerodynamického chování konstrukce.

Skutečné hodnoty síly v newtonech pro každou stranu budovy jsou uvedeny pod touto tabulkou. Přední plocha je vystavena kladné odporové síle 1,7 N, zatímco zadní, pravá a levá plocha jsou vystaveny záporným silám v důsledku sání a bočního odporu. Celková čistá síla ve směru větru činí 3,2 N, což odpovídá celkové síle uvedené v souhrnu sítě. V simulaci byla použita referenční rychlost větru 11 m/s a plocha čela 0,04 m².

Tabulka 4 ukazuje, že program RWIND přesně replikuje referenční průměrné součinitele síly WTG na všech plochách, s odchylkami v rozmezí od 1,95 % do 6,09 % a ‚‘'q=100 % jako míra shody‘‘. Nízká normalizovaná střední kvadratická chyba e2=0,0015 potvrzuje výraznou shodu mezi simulací a měřeními, čímž účinně splňuje validační standardy.

Metrika validace pro součinitele síly mezi WTG a RWIND

Součinitel síly Cf,průměr – WTG Cf – RWIND – k-omega Odchylka (%) n
Přední plocha 0,540 0,562 4,07 1,00
Zadní plocha –0,528 –0,496 6,09 1,00
Pravá boční plocha –0,829 –0,860 3,68 1,00
Levá boční plocha –0,847 –0,893 5,38 1,00
Celková síla ve směru větru 1,070 1,058 1,95 1,00
Celková síla ve směru bočního větru 0,018 0,020 10,00 1,00
Moment ve směru větru 0,582 0,520 10,65 1,00
Moment ve směru bočního větru 0,010 0,0099 1,00 1,00