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

### 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 c_{p,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).

$\mathrm{\Delta p}\left(\mathrm{t}\right)=\frac{1}{2}\xb7\mathrm{\rho}\xb7{\mathrm{v}}^{2}\left(\mathrm{t}\right)\xb7{\mathrm{c}}_{\mathrm{p},\mathrm{mean}}$

ρ |
density of air |

v |
inlet velocity |

c_{p,mean} |
time-averaged pressure coefficient |

t |
time |

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

v(t)² = (v_{x,mean} + v_{x,fluctuation}(t))² + v_{y,fluctuation}(t)² + v_{z,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)² = v_{x,mean}² + 2 ⋅ v_{x,mean} ⋅ v_{x,fluctuation}(t)

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

Δp(t) = 1/2 ⋅ ρ ⋅ v_{x,mean}² [1 + (2 ⋅ v_{x,fluctuation}(t)) / v_{x,mean}] ⋅ c_{p,mean}

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

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

And if you then compare the term v_{x,fluctuation,max} / v_{x,mean} as a multiple g of the turbulence intensity I_{v}(z),

${\mathrm{I}}_{\mathrm{v}}\left(\mathrm{z}\right)=\frac{{\mathrm{\delta}}_{\mathrm{v}}}{{\mathrm{v}}_{\mathrm{mean}}\left(\mathrm{z}\right)}$

δ_{v} |
standard deviation from mean velocity v _{mean} |

v_{mean}(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:

$\mathrm{W}=\frac{1}{2}\xb7\mathrm{\rho}\xb7{\mathrm{v}}_{\mathrm{mean}}^{2}\left(\mathrm{z}\right)\xb7\mathrm{G}\left(\mathrm{z}\right)\xb7{\mathrm{c}}_{\mathrm{p},\mathrm{mean}}$

ρ |
density of air |

v_{mean} |
average inlet velocity |

G(z) |
gust factor depending on altitude |

c_{p,mean} |
time-averaged pressure coefficient |

where

$\mathrm{G}\left(\mathrm{z}\right)=1+2\xb7\mathrm{g}\xb7{\mathrm{I}}_{\mathrm{v}}\left(\mathrm{z}\right)$

g |
factor for defining gust duration |

I_{v}(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 p_{mean} on the model surface depending on an inlet velocity v_{x}(z) by means of a stationary solution of the RANS equations using the SIMPLEC algorithm. Since the mean values of the pressure coefficients c_{p,mean} are based on the ratio between the determined mean pressure values p_{mean} to the undisturbed peak wind velocity pressure at the roof height q(height of roof),

c_{p,mean} = p_{mean} / 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 v_{mean} and the maximum fluctuation component v_{fluctuation}. 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 a 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.

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

#### Links

Write Comment...

Write Comment...

Contact Us

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

Wind Simulation Using RWIND Simulation and Transferring Wind Forces to RFEM or RSTAB

Computer technology has a firm grip on digital structural analysis and design. With each new development, the planners involved are able to increase the limits of what is feasible.

New

RWIND 2 | Features of RWIND Pro

- Calculation of transient incompressible turbulent wind flow with the BlueDyMSolver solver
- LES SpalartAllmarasDDES turbulence model
- Consideration of stationary solution as initial state for transient calculation
- Automatic determination of analysis period and time steps
- Use of intermediate results during the calculation
- Organized display of time-varying results via time step units
- Diagram of drag force and point probe results over analysis time
- Display of line probe results for any time steps in a diagram

- Is it possible to keep the results from RWIND Simulation when deleting the FE mesh if I only adjust materials or support conditions?
- I have defined a wind profile for RWIND Simulation up to a height of 100 m. Is the wind profile dependent on the size of the wind tunnel?
- How can I import an IFC model into the RWIND Simulation environment?
- I would like to use the export of visual objects to RWIND Simulation. How do I scale my visual object for the import?
- Where are the RWIND Simulation projects created by RFEM/RSTAB saved?

- Where can I define a wind velocity diagram according to a specific wind standard for the numerical wind tunnel analysis in RWIND Simulation?
- Where can I adjust the settings of my wind simulation in RFEM 6 / RSTAB 9?
- Is it possible to use load transfer surfaces in RFEM 6 for the wind load generation using RWIND?
- Can I also use RWIND 2 with RFEM 5/RSTAB 8?
- How can I open RWIND using the load wizard in RFEM and RSTAB?

Associated Products