TLD Tuned Liquid Dampers

FEA for TLD Tuned Liquid Dampers

One of the greatest challenges when it comes to design steel stack structures such as chimneys, is to deal with oscillating loads as dynamic wind forces and seismic spectrum, as commented in this article. These loads threaten the design with coupling the excitation with its natural frequency, leading the structure to failure through resonance due to the high stresses produced by the oscillation amplitude.

Two main strategies can be followed to address this issue: either increasing the natural mode of the structure or damping its response to the excitation.

The objective of the first way (natural mode improvement) is to run away from the dangerous resonance frequencies. 

Figure 1 - Vortex induced oscillating loads on a cylindrical profile due to the wind action

This can be achieved by geometrical changes on the design such as plate thickness enhancement or a gain of chimney’s diameter, or by the addition of new structural elements as stiffeners or supports. Although this method is conceptually simple and does not require an extra designing effort other than adding some more material, is has other limitations. Usually, geometrical changes as wider diameter or the addition of new elements are restricted by other features of the whole project. 
 

Figure 2 - TLD pools coupled oscillation

Neighbour structures limiting the available volume or chimney’s exhaust properties are examples of these limitations. Other inconvenience for this solution, and the one who uses to have the last word, is the cost-efficiency. An increment on plate’s thickness involve a not negligible increment on the project’s expense moving it away further than the structure’s natural mode.

The other path to follow in order to guard the structure against oscillations is dramatically increasing the structure’s damping by the addition of a damping system. 

Even though there are different methods to achieve this, this article will be focused on the tuned liquid dampers (TLD).

The mission of the TLD is to dissipate the energy absorbed by the structure from oscillating loads, through a mass-spring-damper equivalent system that moves desynchronized with the chimney at its natural frequency. 

It is, in essence, a series of pools connected to the structure that moves against its pendulum movement. The phenomenon responsible of the energy mitigation effect inside the pools is called sloshing, and is defined as the movement of liquid inside another object. An example of sloshing familiar to almost everyone can be seen every time you walk with a soup plate. 

The human gait frequency (acting as the excitation) is close to the plate mode, so each step amplifies the soup’s wave motion, finally causing it to ruin your menu by splitting out of the rim. 

Figure 3 - Amplitude vs Frequency graph comparing a chimney with and without damper

On the other hand, if you are carrying a cup of coffee at a standard walking speed, despite it creates waves on the surface, the excitation frequency remains lower that the coffee’s natural mode, and as long as you don’t speed up, you breakfast will be safe.

Figure 4 - Plate & cup first oscillating mode shapes

The parameters that define geometrically a TLD, controlling its behaviour, are its interior [Rint] and exterior radius [Rext], the number of pools in which the TLD is divided, and the liquid height [H] that controls the liquid volume and thus the liquid mass. This liquid mass is divided in two distinguishable masses: the impulsive mass, the one that moves creating the sloshing phenomenon, and the convective mass, that has no relative movement with the TLD. Note that despite it has its influence on the sloshing frequency, the amount of convective mass also acts as dead load on the chimney, moving its natural frequency, so it has to be taken into consideration when tuning the assembly.

Figure 5 - TLD main dimensions and water mass distribution

The interior radius is used as the start point to calculate the previously mentioned variables [as it is more or less set by the chimneys geometry itself] along with the chimneys natural frequency and modal mass, acting as the focus result for the tune and as a reference for the total liquid mass respectively. Then, using the linear wave theory and the Graham & Rodríguez equations, the first group of parameters is obtained. Because of the influence of the TLD’s convective mass on the base chimney’s natural mode, an interpolation process has to be performed in order to obtain a fine tuning for the liquid damper.

Figure 6 - TLD sloshing mode shape at tuned frequency

To optimize this method’s accuracy and asses its results, SDEA’s advanced engineering team assembled a coupled FEA modal analysis plus harmonic response computational model. A fluid-structure interaction model along with a fluid free surface boundary condition is needed to capture the sloshing effect on the TLD’s pools. The next image shows the pools from a TLD where it can be appreciated that the liquid from all pools move in the same direction when they get coupled with the chimney’s response to external excitation.

Although this solution for the problem introduced in this article has a bigger technical load than enhancing the structure plate’s thickness, the cost-efficiency and adaptability of the damper installation compensate these cons by relieving the projects budget and providing certain design flexibility to overcome hypothetical installation issues.

FSI STEEL STACK STRUCTURAL FEA ANALYSIS

Steel stack wind contours

Chimney miniature

When designing steel stack structures, several rules from standards as “ASME STS-1” or “ASCE 7” have to be taken into consideration to assure security during its useful life. This is particularly important for those located in areas subjected to heavy environmental loads, such as high activity seismic zone 4, due to the slender nature of this kind of structures thus low resonant primary modes prone to be easily excited by seismic or wind actions. In this article a brief outline about the designing process and considerations is provided based on the experience in FEA use gathered by the SDEA engineering team along different projects.

The main goal to achieve for the structure is to keep its design stresses under certain load combinations below a maximum allowable value. That is a common design strategy through many sectors; the difference here is the forces acting over the structure are somehow easy to misunderstand.

The first step is to determine the load combinations that the structure is going to be subjected to. Standards as “ASCE 7” or national regulations define the basic group of combinations to be considered during the assessment of the design, for example:

Load Combinations

The whole set of load combinations should be taken from the standards, from project’s characteristics or from customer’s requirements.

Once the load combinations are defined and the most relevant loads are identified, it’s time to calculate each of these values separately by using a finite element analysis (FEA) method over a detailed 3D model. Next, a summary of three of the most common loads shared for the most part of the projects is presented:

The first one is the dead load, and can be defined as the weight of all material incorporated into the structure, including its fixed permanent equipment such as ladders, platforms, etc. A constant when calculating stresses for different loads is that the worst situation is the one to be considered. Thus, in this case unlike the following ones, effects that reduce the total mass such as corrosion have not to be taken into consideration to not underestimate the stresses induced by the dead load expected in a real situation. A static structural FEA model can be used to precisely calculate the stresses induced by this load on every section of the stack.

The second common load is the wind load, or better said, are the wind loads. Wind load can be divided in two different loads: static and dynamic. The first one represents the action of the constant part of the wind pressure and is calculated using a static FEA model where a height dependant pressure load is applied to upwind and leeward surfaces. This pressure is obtained using standard’s equations, and usually depends on factors as the local wind speed, the geometry of the structure and the height of the point where the pressure is being applied.

The other wind load is the dynamic one, and it has to be considered due to the lightweight and flexible nature of steel stacks. It is defined as the load produced by the oscillations generated by the wind-structure interaction. “ASME STS-1” standard contains a set of rules to calculate the dynamic wind response. A more complex CFD fluid structure interaction model should be used in this case to assess the structure-vortex interaction in order to then evaluate the stresses produced. More details about the dynamic wind can be seen in this article.

Static wind profile

Last but not least, we have the earthquake load, defined as the horizontal and vertical loads related to the response of the structure to seismic motions, and are calculated using the response spectrum method. In terms of the FEA model, a coupled modal analysis plus response spectrum model is required. The pseudo-spectral acceleration (Sa) used as input for the seismic response spectrum is calculated according to the following equation:

Spectral acceleration

Where g is the gravity acceleration. R and U are factors that depend on the structure itself, being R the seismic force reduction factor that depends on the geometric characteristics of the structure and U the use factor, which defines the category of the building depending on its function and importance. The ground environment is defined by the factors Z and S. The seismic zone factor Z depends on the magnitude of the seismic forces presented on the region where the stack structure is going to be standing, and the ground factor S is a value used to define the geological conditions of the terrain.

U.S. Seismic Zones Map

The last factor C is the seismic amplification factor. It is the non-constant value that provides Sa with the frequency dependant behaviour. An example of a real pseudo-spectral acceleration is showed in the next graph:

Seismic pseudo spectral acceleration

Both the dynamic and earthquake loads have a paramount importance when designing steel stack structures due to its slender and flexible nature. Frequency excitations close to the structure’s natural mode of vibration can lead it to failure because of the appearance of resonance phenomenon. To mitigate this effect, cost-effective solutions as tuned liquid dampers (TLD, as explained in this article) are widely used.

Once every relevant load is calculated and combined according to appropriate normative, the maximum allowed values have to be obtained to assess the design, with the idea of improving it in case of not meeting the requirements, or optimizing it reduce the project expenses! These value’s calculation methodology is described in standards such as “ASME STS-1”, and are directly related to the structure’s parameters such as diameter, plate thickness, material properties, etc. Various values have to be obtained depending on the project conditions, being the most common maximum allowable values for longitudinal compression and the longitudinal compression and bending combination.