Wind Speed and Turbulence Intensity Profile for Determination of Quasi-Static Wind Loads According to Gust Concept

Technical Article on the Topic Structural Analysis Using Dlubal Software

  • Knowledge Base

Technical Article

Structures react differently to wind action depending on stiffness, mass, and damping. A basic distinction is made between buildings that are prone to vibration and those that are not.

Usually, structures are not considered to be susceptible to vibrations if the deformations under wind action by gust resonance are not increased by more than 10% [2]. In this case, the time-variable wind action can be described as a static equivalent load.

Assuming that the turbulences within the wind flow are very large in relation to the building dimensions, a statically acting distribution of pressure p on a building geometry can be calculated with RWIND Simulation according to the "quasi-stationary method" or the gust concept [3].

Basically, a stationary flow field around the analysis model is assumed for the turbulent speed fluctuation over the gust's duration [3]. The pressure fluctuation on the model surface from the inflow turbulence is thus seen as a state stationary over the certain time period t. Thus, the fluctuations follow exactly the course of the time-averaged pressure coefficients cp,mean on the model surface.

The resulting wind-induced pressure Δp(t) on the model surfaces then depends purely on the inlet velocity v(t).

Wind-Induced Pressure as Function of Time

Δp(t) = 12 · ρ · v2 (t) · cp,mean

ρ Density of air
v Inlet velocity
cp,mean Time-averaged pressure coefficient
t Time

Thus, the value of the inlet velocity vector v(t) is:

v(t)² = (vx,mean + vx,fluctuation(t))² + vy,fluctuation(t)² + vz,fluctuation(t)²

If the squared terms only make a small contribution, the effective value of the inlet velocity vector v(t) is the result:

v(t)² = vx,mean² + 2 ⋅ vx,mean ⋅ vx,fluctuation(t)

Using the effective inlet velocity in the equation of the wind-induced pressure results in:

Δp(t) = 1/2 ⋅ ρ ⋅ vx,mean² [1 + (2 ⋅ vx,fluctuation(t)) / vx,mean] ⋅ cp,mean

This transformation shows that the fluctuation of the wind pressure Δp(t) only depends on the fluctuation of the wind speed vx,fluctuation(t) in the main inflow direction x.

If you replace the time-variable speed fluctuation vx,fluctuation(t) by the maximum occurring speed fluctuation vx,fluctuation,max, you remove the temporal variability from the system.

If you then compare the term vx,fluctuation,max / vx,mean as a multiple g of the turbulence intensity Iv(z),

Turbulence Intensity as Function of Altitude

Iv(z) = δvvmean (z)

δv
Standard deviation from mean velocity vmean
vmean(z) Average velocity depending on altitude
z Height above ground

you can describe the term in square brackets as the gust factor G(z). Inserting the terms into the nominal wind load equation results in:

Nominal Wind Load

W = 12 · ρ · vmean2 (z) · G(z) · cp,mean

ρ Density of air
vmean Average inlet velocity
G(z) Gust factor depending on altitude
cp,mean Time-averaged pressure coefficient

where

Gust Factor

G(z) = 1 + 2 · g · Iv (z)

g Factor for defining gust duration
Iv(z) Turbulence intensity as a function of altitude
z Height above ground

For example, in EN 1991‑1‑4, factor g is used for describing the gust duration 3.5.

RWIND Simulation calculates the mean values of the pressures pmean on the model surface depending on an inlet velocity vx(z) by means of a stationary solution of the RANS equations using the SIMPLEC algorithm. Since the mean values of the pressure coefficients cp,mean are based on the ratio between the determined mean pressure values pmean to the undisturbed peak wind velocity pressure at the roof height q(height of roof),

cp,mean = pmean / q(height of roof)

it is possible to use the inlet velocity from the converted peak wind velocity pressure q(z) over the height to determine the nominal wind loads according to the gust concept [1].

v(z) = √(2 ⋅ q(z) / ρ)

Thus, this wind speed includes the mean wind velocity vmean and the maximum fluctuation component vfluctuation. In this case, the inflow turbulence intensity can be set constantly over the height to a very small value of about 5% [4].

When considering the effects of forces acting on the entire building or on large surface areas, this method provides a very good approximation to the natural wind loading [3]. The reason is that the small turbulence effects masked by averaging act only in partial areas and do not have a noticeable effect due to the global integration of the force values.

Furthermore, the concept reacts very well even for small partial areas with frontal inflow, since here the effective pressure fluctuations are already very well recorded in the peak wind speed profile [3].

Conversely, the system results in poorer convergence with the reality for surfaces with flow separation (side and rear walls). It is especially in these zones that the building-induced turbulence "faded away" by averaging using the gust concept has a greater effect than the inflow turbulence effect contained in the inlet velocity profile.

Author

Dipl.-Ing. (BA) Andreas Niemeier, M.Eng.

Dipl.-Ing. (BA) Andreas Niemeier, M.Eng.

Mr. Niemeier is responsible for the development of RFEM, RSTAB, and the add-on modules for tensile membrane structures. Also, he is responsible for quality assurance and customer support.

Keywords

Peak wind speed Gust wind speed Basic wind velocity Mean wind speed Turbulence intensity Wind pressure

Reference

[1]   Eurocode 1: Actions on structures - Part 1-4: General actions - Wind actions; German version EN 1991-1-4:2005 + A1:2010 + AC:2010
[2]   Albert, A. (2016). Schneider - Bautabellen für Ingenieure mit Berechnungshinweisen und Beispielen (22nd ed.). Bochum: Bundesanzeiger.
[3]   Kiefer, H. (2003). Windlasten an quaderförmigen Gebäuden in bebauten Gebieten (Dr.-Ing.). Karlsruhe Institute of Technology.
[4]   Werth, M. (2019). Vergleichende Studie zu Windlastmodellen im Hochbau: Numerische Strömungsberechnung vs. Druckmessungen im Windkanal (MSc.). Graz University of Technology.

Links

Write Comment...

Write Comment...

  • Views 2420x
  • Updated 23 March 2022

Contact Us

Contact Dlubal

Do you have any further questions or need advice? Contact us via phone, email, chat or forum or find suggested solutions and useful tips on our FAQ page, available 24/7.

+49 9673 9203-0

[email protected]

Event Invitation

International Mass Timber Conference

Conference 27 March 2023 - 29 March 2023

CFD Wind Simulation with RWIND 2

CFD Wind Simulation with RWIND 2

Webinar 29 June 2022 2:00 PM - 3:00 PM EDT

Tensile Membrane Structure Design in RFEM 6

Tensile Membrane Structure Design in RFEM 6

Webinar 17 March 2022 2:00 PM - 3:00 PM EDT

Wind Load Generation Based on CFD in RWIND 2

Wind Load Generation Based on CFD in RWIND 2

Webinar 4 March 2022 10:00 AM - 11:00 AM CET

Online Training | English

Eurocode 3 | Steel Structures According to DIN EN 1993-1-1

Online Training 9 December 2021 8:30 AM - 12:30 PM CET

Online Training | English

Eurocode 5 | Timber Structures According to DIN EN 1995-1-1

Online Training 23 September 2021 8:30 AM - 12:30 PM CEST

Online Training | English

Eurocode 3 | Steel Structures According to DIN EN 1993-1-1

Online Training 25 August 2021 8:30 AM - 12:30 PM CEST

Online Training | English

RFEM for Students | USA

Online Training 11 August 2021 1:00 PM - 4:00 PM EDT

Online Training | English

RFEM | Structural Dynamics and Seismic Design According to EC 8

Online Training 11 August 2021 8:30 AM - 12:30 PM CEST

Online Training | English

Eurocode 2 | Concrete structures according to DIN EN 1992-1-1

Online Training 29 July 2021 8:30 AM - 12:30 PM CEST

Online Training | English

RFEM | Basics

Online Training 13 July 2021 9:00 AM - 1:00 PM CEST

Online Training | English

RFEM | Basics | USA

Online Training 17 June 2021 9:00 AM - 1:00 PM EDT

Online Training | English

RFEM for Students | Part 3

Online Training 15 June 2021 2:00 PM - 4:30 PM CEST

Glass Design with Dlubal Software

Glass Design with Dlubal Software

Webinar 8 June 2021 2:00 PM - 2:45 PM CEST

Online Training | English

RFEM | Structural dynamics and seismic design according to EC 8

Online Training 2 June 2021 8:30 AM - 12:30 PM CEST

Online Training | English

Eurocode 5 | Timber structures according to DIN EN 1995-1-1

Online Training 20 May 2021 8:30 AM - 12:30 PM CEST

Online Training | English

RFEM for Students | Part 2

Online Training 17 May 2021 2:00 PM - 4:30 PM CEST

Blast Time History Analysis in RFEM

Blast Time History Analysis in RFEM

Webinar 13 May 2021 2:00 PM - 3:00 PM EDT

Online Training | English

Eurocode 2 | Concrete structures according to DIN EN 1992-1-1

Online Training 12 May 2021 8:30 AM - 12:30 PM CEST

Timber Structures | Part 2: Design

Timber Beam and Surface Structures | Part 2: Design

Webinar 11 May 2021 2:00 PM - 3:00 PM CEST

RFEM 5
RFEM

Main Program

Structural engineering software for finite element analysis (FEA) of planar and spatial structural systems consisting of plates, walls, shells, members (beams), solids, and contact elements

Price of First License
3,950.00 EUR
RSTAB 8
Structural Frame Analysis Software RSTAB

Main Program

Structural engineering software for designing frame, beam, and truss structures, as well as performing linear and nonlinear calculations of internal forces, deformations, and support reactions

Price of First License
2,850.00 EUR