002059

Mahyar Kazemian, M.Sc.

29.6.2026

Výpočet lokálních a celkových součinitelů tlaku (cp,1 a cp,10) pomocí extrémní analýzy v časové oblasti

Koeficienty tlaku větru jsou základními parametry při navrhování obvodových plášťů budov, střešních plášťů a konstrukčních prvků. Tradiční normy pro zatížení větrem poskytují tabelované hodnoty pro místní součinitele tlaku (\(c_{p,1}\)) a celkové součinitele tlaku (\(c_{p,10}\)), které jsou obvykle odvozeny z měření v aerodynamických tunelech a statistického zpracování. Moderní výpočetní techniky však umožňují odhadnout tyto součinitele přímo z časových průběhů tlaku získaných z měření, CFD simulací nebo syntetické generace turbulence.

1. Úvod

Tlak větru působící na povrchy budov je vysoce proměnlivý v důsledku atmosférické turbulence a jevů odtržení proudění. Zatímco okamžité tlaky mohou vykazovat velké špičky, konstrukční návrh vyžaduje reprezentativní extrémní hodnoty spojené s konkrétními zatěžovacími plochami. Dva běžně používané součinitele jsou [1]:

* cp,1: Představuje lokální špičkové tlaky působící na malé prvky opláštění s účinnou plochou přibližně 1 m²

* cp,10: Představuje tlaky působící na větší konstrukční oblasti s účinnou plochou přibližně 10 m²

Protože větší plochy průměrují struktury malé turbulence, hodnoty cp,10 jsou obecně méně extrémní než hodnoty cp,1. Úkolem je transformovat vysoce proměnlivé tlakové signály na statisticky významné návrhové součinitele. Tento článek představuje kompletní metodiku v časové oblasti pro odvození cp,1 a cp,10 pomocí:

Kvazistacionárního modelování tlaku
Plošného průměrování
Filtrování po dobu trvání zatížení
Extrakce extrémních hodnot
Přizpůsobení Gumbelovu rozdělení
Korekce součinitelem zatížení nárazem větru

Metodika vychází z rámce teorie pseudostacionárního zatížení větrem a poskytuje fyzikálně konzistentní postup pro převod přechodných tlakových fluktuací na návrhové součinitele tlaku.

1 Algoritmus Cp,10 a Cp,1 z přechodové CFD

Vstupní parametry, doporučené rozsahy hodnot a odvozené veličiny použité v analýze extrémního součinitele tlaku

2 Vstupní parametry, doporučené rozsahy hodnot a odvozené veličiny použité v analýze extrémního tlakového součinitele

2. Model kvazistacionárního tlaku větru

Výchozím bodem je kvazistacionární předpoklad, který dává do přímé souvislosti tlakové fluktuace a fluktuace rychlosti. Okamžitá rychlost větru je:

V (t) = U + u^{'} (t)

U	Střední rychlost větru
u'(t)	Turbulentní fluktuace

Kvazistacionární model tlaku větru

Odpovídající tlak se vypočítá jako:

w (t) = 0.5 ρ V^{2} (t) \cdot {\bar{c}}_{p} = 0.5 ρ U^{2} \cdot {[1 + \frac{u^{'} (t)}{U}]}^{2} \cdot {\bar{c}}_{p}

cp,bar	Střední součinitel tlaku
ρ	Hustota vzduchu

Kvazi-stacionární tlak

Tato formulace vystihuje nelineární vztah mezi fluktuacemi rychlosti a aerodynamickým tlakem. Pro oblasti sání na střeše je cp,bar typicky záporný.

3 Okamžitá rychlost větru

3. Generování časových průběhů tlaku

Signál turbulentní rychlosti může být získán z:

Měření v aerodynamickém tunelu
CFD simulací
Generátorů syntetické turbulence
Terénních měření

Výsledným časovým průběhem součinitele tlaku je cp(t). V praxi vznikají další tlakové fluktuace z:

Odtržení proudění
Odtrhávání vírů
Lokálních rohových vírů
Interakcí v turbulentním úplavu

Tyto efekty jsou superponovány na kvazistacionární tlakovou složku, aby se vytvořil realistický tlakový signál.

4 Quasi-Steady Pressure Time History

4. Plošné průměrování a účinky doby trvání zatížení

Normy pro zatížení větrem uznávají, že tlakové špičky jsou méně extrémní, jak se zatěžovaná plocha zvětšuje. Aby se tento jev zohlednil, součinitele tlaku jsou filtrovány po ekvivalentní dobu trvání zatížení:

c_{p, A} (t) = \frac{1}{τ} \int_{0}^{τ} c_{p} (t - \bar{t}) d \bar{t}, (τ = \frac{4.5 \sqrt{A}}{V_{m}})

A	Přidružená plocha
Vm	Průměrná rychlost větru

Filtrované součinitele tlaku

Typické hodnoty jsou běžně definovány na základě zatěžovací plochy zatěžovaného povrchu, kde cp,1 představuje součinitel tlaku spojený s účinnou plochou 1 m², zatímco cp,10 odpovídá účinné ploše 10 m². Jak se zatěžovací plocha zvětšuje, související doba trvání zatížení a doba průměrování se také prodlužují, což vede k silnějšímu vyhlazovacímu efektu na tlakový signál. V důsledku toho delší doba průměrování spojená s cp,10 odfiltruje větší část vysokofrekvenčních tlakových fluktuací způsobených lokální turbulencí a odtržením proudění. Toto snížení krátkodobých tlakových špiček vede k nižším extrémním hodnotám sání ve srovnání s cp,1, což činí cp,10 obecně méně extrémním a reprezentativnějším pro tlaky působící na hlavní konstrukční prvky spíše než na lokální části opláštění.

Vyhlazování tlakových součinitelů pomocí klouzavého průměru pro hodnocení cp,1 a cp,10

5 Filtrování klouzavým průměrem součinitelů tlaku pro vyhodnocení cp,1 a cp,10

Vliv Doby Trvání Zatížení na Časové Historie Součinitele Tlaků

6 Vliv Délky Trvání Zatížení na Časové Průběhy Součinitele Tlaků

5. Extrakce extrémních tlaků

Po plošném průměrování jsou filtrované záznamy tlaku rozděleny do bloků o stejné délce trvání. Pro každý blok je extrahován minimální součinitel tlaku cp,min. Tato bloková minima představují statisticky nezávislé extrémní události. Tento postup transformuje tisíce vzorků tlaku do zvládnutelného souboru extrémních hodnot vhodných pro pravděpodobnostní analýzu.

7 Extrahované extrémní hodnoty sání

6. Teorie extrémních hodnot a Gumbelovo rozdělení

Extrémy tlaku větru jsou běžně modelovány pomocí Gumbelova rozdělení. Pro výběrový soubor extrémních hodnot: X=−cp. Parametry Gumbelova rozdělení jsou odhadnuty z výběrového průměru a směrodatné odchylky. Distribuční funkce je:

F (x) = \exp [- \exp (- \frac{x - μ}{β})]

Kumulativní distribuční funkce

kde:

μ = parametr polohy
β = parametr měřítka

Gumbelův model poskytuje spojitou reprezentaci pravděpodobnosti extrémních tlakových událostí.

7. Vyhodnocení 78% fractilu

Normy pro zatížení větrem často definují charakteristické součinitele tlaku pomocí 78% fractilu rozdělení extrémních hodnot. Charakteristická extrémní hodnota je:

c_{p, (78)} = μ + 1.4 β

Charakteristická extrémní hodnota

78% fractil poskytuje robustní odhad návrhových tlakových špiček, přičemž se vyhýbá nadměrné citlivosti na izolované odlehlé hodnoty.

Gumbelova rozdělení extrémních součinitelů tlaku a odpovídající 78% kvantilové hodnoty

8 Gumbelovo rozdělení extrémních součinitelů tlaku a odpovídající hodnoty 78% kvantilu

8. Součinitel zatížení nárazem větru

Součinitel zatížení nárazem větru zohledňuje zesílení v důsledku atmosférické turbulence. Součinitel špičkové rychlosti a součinitel zatížení nárazem větru se získají z:

g_{u} = \sqrt{2 \ln (ν_{u} T)} + \frac{0.5772}{\sqrt{2 \ln (ν_{u} T)}}

G_{w, l i n} (A) = 1 + 2 g_{u} I_{u} \sqrt{B^{2} (A)}

νu	Průměrná četnost vzestupných překročení
T	Průměrné období
Iu	intenzita turbulence
B²(A)	Faktor pozadí

Součinitel zatížení větrem

Tato formulace zahrnuje charakteristiky turbulence a efekty prostorového průměrování.

9. Výpočet konečných součinitelů tlaku

Konečné součinitele tlaku cp,1 a cp,10 se získají vydělením charakteristických špičkových součinitelů tlaku odpovídajícím součinitelem zatížení nárazem větru. Tím se odstraní zesílení způsobené turbulencí větru a převedou se špičkové tlaky na ekvivalentní pseudostacionární součinitele tlaku. Výsledné hodnoty představují statické aerodynamické zatěžovací účinky, které lze přímo použít pro konstrukční návrh a výpočty zatížení větrem podle norem.

c_{p, 1} = \frac{c_{p, (78), 1}}{G_{w, lin} (A_{1})}, c_{p, 10} = \frac{c_{p, (78), 10}}{G_{w, lin} (A_{10})}

Součinitele Konečného Tlakového Zatížení

Výsledné součinitele představují ekvivalentní statické součinitele tlaku vhodné pro konstrukční návrh.

Konečné pseudostacionární tlakové součinitele ∣cp,1∣ a ∣cp,10∣ získané analýzou extrémních hodnot

9 Konečné pseudostacionární součinitele tlaku ∣cp,1∣ a ∣cp,10∣ získané z analýzy extrémních hodnot

10. Ekvivalentní zatížení větrem

Konečné ekvivalentní zatížení větrem působící na zatěžovací plochu je

W_{l i n} = \frac{1}{2} ρ U^{2} c_{p} G_{w, l i n} A

Ekvivalentní zatížení větrem

Tato rovnice transformuje součinitele tlaku na návrhová zatížení přímo použitelná v modelech pro statickou analýzu konstrukce.

11. Interpretace cp,1 a cp,10

Metodika přirozeně reprodukuje dobře známý efekt redukce plochou:

cp,1

Citlivý na lokální víry
Vyšší špičky sání
Používá se pro návrh opláštění a spojovacích prostředků
Krátká doba průměrování

cp,10

Reprezentuje větší zatížené oblasti
Redukovaná špička sání
Používá se pro návrh hlavních konstrukčních prvků
Delší doba průměrování

Jak se zatěžovací plocha zvětšuje, vysokofrekvenční tlakové fluktuace jsou postupně odfiltrovány, čímž se snižují extrémní součinitele tlaku.

12. Výhody metodiky v časové oblasti

Ve srovnání s tradičními normovými přístupy nabízí prezentovaný postup několik výhod:

Použitelný pro tlaková pole generovaná CFD
Použitelný pro měření v aerodynamickém tunelu
Explicitní zpracování intenzity turbulence
Fyzikální reprezentace plošných efektů
Konzistentní statistika extrémních hodnot
Vhodný pro libovolné zatěžovací plochy
Přímé generování pseudostacionárních návrhových zatížení

Metodologie překlenuje propast mezi daty přechodného aerodynamického tlaku a praktickými parametry pro konstrukční návrh.

Závěr

Byl představen komplexní rámec v časové oblasti pro odvození lokálních a konstrukčních součinitelů tlaku, cp,1 a cp,10, ze záznamů přechodného tlaku větru. Metodologie kombinuje teorii kvazistacionárního zatížení větrem, plošné průměrování, statistiku extrémních hodnot a korekce součinitelem zatížení nárazem, aby generovala fyzikálně smysluplné návrhové součinitele. Tento přístup je zvláště cenný pro CFD inženýrství větru, kde jsou k dispozici kompletní časové průběhy tlaku a tradiční normové součinitele nemusí adekvátně reprezentovat složité geometrie. Integrací účinků doby trvání zatížení a Gumbelovy analýzy extrémních hodnot poskytuje tento postup rigorózní cestu od přechodných tlakových fluktuací k ekvivalentním statickým návrhovým zatížením.

Autor

Mahyar je zodpovědný za vývoj produktů a marketing společnosti Dlubal Software, zejména za RWIND 2. Spojuje technický vývoj s dalším rozvojem komunikace produktu.

Reference

Picozzi, V., et al. (2022). The relationship between wind pressure and pressure coefficients for the definition of wind loads on buildings. Buildings 12(2). p. 225.

;