﻿ Theoretical Study and Parameter Optimization of CFST Pillar in Cooling Tower

ISSN: 2472-0437

## Journal of Steel Structures & Construction

Reach Us +44-1522-440391

# Theoretical Study and Parameter Optimization of CFST Pillar in Cooling Tower

Xiaoxiong Zha*, Peichang Xu and Dejin Chen
School of Civil and Environment Engineering, Shenzhen Graduate School, Harbin Institute of Technology, Shenzhen, China
*Corresponding Author: Xiaoxiong Zha, School of Civil and Environment Engineering, Shenzhen Graduate School, Harbin Institute of Technology, Shenzhen, China, Tel: +86-158-1685-3765, Email: [email protected]

Received Date: Nov 19, 2017 / Accepted Date: Nov 26, 2018 / Published Date: Dec 03, 2018

### Abstract

The application of concrete pillar in cooling towers has become more common. However, the CFST (Concrete Filled Steel Tube) pillar has no relevant applications and researches in cooling towers. In order to research the CFST pillar, we will analyze the optimization of the parameters of the CFST pillar first. In CFST pillar, the parameters that can be optimized are: The position of the intersection (above, below, or in the middle), the cross angle, the tilt angle of the CFST pillar, the form of the intersection, the form of the pillar section, etc. Changes in these parameters can also cause a large change in the CFST pillar internal force. The purpose of this paper is to obtain the optimal value range of “intersection position, cross angle, and CFST pillar tilt angle” in the CFST pillar parameters through theoretical calculations, and to provide parameters selection suggestions for cooling tower engineering design. The main theories adopted are: classical structural mechanics (e.g., force method, displacement method, moment distribution method, graph multiplication, etc.), moment distribution method of space rigid frame, and transformation rule of space coordinates. We analyze the CFST pillar by three-dimensional structural mechanics and use MatLab to obtain the optimized results of CFST pillar parameters.

Keywords: CFST pillar; Structural mechanics; MatLab; Intersection; Configuration; Horizontal; Load

#### Background

The parameters related to the construction of the CFST X-pillar include: the cross angle of the CFST X-pillar is α, the angle between the pillar and the ground is β, the length of the pillar is L, and the distance between the pillar intersection and the upper end is λL.

This article will use the method of classical structural mechanics, combined with the transformation rules of spatial coordinates, to derive the CFST X-pillar internal force calculation formula and then use calculation formula to optimize the parameters of CFST X-pillars.

This paper only considers the static load of the CFST X-pillar, mainly including the gravity load G of the superstructure, and the equivalent static wind load P (Figure 1).

Figure 1: CFST X-pillar configuration parameters.

#### Structural Analysis

Gravity load distribution: CFST X-pillars are generally spaced at regular intervals, then the gravity load G will be evenly distributed to each CFST X-pillar. Suppose that the number of CFST X-pillars in the cooling tower is n, then G1=G2=…=Gm=…=Gn (Figure 2).

Figure 2: Cooling tower overall model.

Horizontal load distribution: The lateral stiffness of different CFST X-pillar is different along the direction of the horizontal load P, so the horizontal load P cannot be simply distributed evenly to each CFST X-pillar (Figure 3).

Figure 3: A rigid frame model consisting of CFST X-pillars and a ring beam.

Considering that the upper structure is larger and the deformation is allowed to be smaller, the ring beam can be regarded as a rigid body. Then under the action of the horizontal load P, the ring beam will move in its entirety, and it will also bring about the same horizontal displacement Δp at the top of all CFST X-pillars. Suppose the CFST X-pillar closest to the horizontal load P is the first pillar. Counterclockwise, the lateral stiffness of the i-th CFST X-pillar is KSi and the horizontal load it is assigned to is Pi.

There is the equation,

Therefore, the load distribution value of the horizontal load P at the top of each CFST X-pillar is determined according to the ratio between KSi. In addition, the lateral stiffness Ks of the X pillar is mainly composed of two parts: the in-plane stiffness Kin and the out-of-plane stiffness Kout (Figure 4) the in-plane stiffness is Kin along the x-direction, and the out-of-plane stiffness is Kout along the y-direction.

Figure 4: In-plane stiffness and Out-of-plane stiffness of CFST X-pillar.

Similarly, the horizontal load P can be decomposed into an inplane horizontal load Pin and an out-plane horizontal load Pout. The final load condition of a single CFST X-pillar can be represented as shown in Figure 5.

Figure 5: CFST X-pillar load conditions.

Conclusion: In order to obtain the lateral stiffness Ks in any direction of the CFST X-pillar, the in-plane rigidity Kin and the out-ofplane stiffness Kout of the CFST X-pillar are required to be calculated. And then according to Ks, the load distribution of the horizontal load P at the top of each CFST X-pillar can be obtained in the overall structure.

Boundary condition analysis

When the stiffness of the ring beam is assumed to be infinite, the boundary condition at the top of the CFST X-pillar can be considered as the slip bearing (Figure 6).

Figure 6:CFST X-pillar boundary condition.

Basic assumption

Basic assumptions for theoretical analysis in this paper:

• The ring beam above the CFST X-pillars does not allow large deformations. So it can be assumed to be a rigid body.

• Suppose EA=∞ in the analysis of internal forces.

#### In-plane Stiffness Kin (Two-dimensional calculation)

According to the simplified calculation model, the X pillar is fixed at the bottom and the top of it is the sliding bearing. Because the force and deformation of the CFST X-pillars are both in-plane, the calculation of the in-plane stiffness Kin can be reduced to a two-dimensional problem (Figure 7).

Figure 7: CFST X-pillar in-plane stiffness calculation model.

AECD is a CFST X-pillar with intersection B. Points C and E are sliding bearings. Points A and D are fixed bearings. A rigid body CE is connected between points C and E at the top of the pillar. The length of the CD is L. CB/CD=λ. The load of the pillar can be equivalent as shown in Figure 8.

Figure 8: CFST X-pillar in-plane stiffness equivalent model.

The expression of internal force of bending moment can be obtained by moment distribution method and displacement method [1]. Here only the internal force of the BC rod end is given:

Then calculate the internal force of the bending moment of the basic structure under the force of unit “1”. Here only the internal force of the BC rod end is given:

Finally, CFST X-pillar top displacement Δin can be obtained by the method of graph multiplication under in-plane force Pin.

So the in-plane stiffness is,

#### Moment Distribution Method for Space Rigid Frames

When calculating the out-of-plane bending stiffness Kout of the CFST X-pillar, we should use the moment distribution method for space rigid frames. Regarding to the moment distribution method for space rigid frames, Xiling Ma et al. has done relevant research [2]. Their research results are only suitable when structural members are perpendicular to each other. Here we will use the law of spatial coordinate transformation to derive the case for structural members that are not perpendicular to each other.

As shown in Figure 9, O′x′y′z′ is the local coordinate system and Oxyz is the global coordinate system. The angular displacement components of the member ij in the local coordinate system are: θx′i, θy′i, θz′i, θx′j, θy′j, θz′j. Among them, θx′i and θx′j are the twist angles of the member, and θy′i, θy′j and θz′i, θz′i are the rotation angles of the end surfaces around the y-axis and the z-axis, respectively, when the member ij is bent. Figure 10 shows the six moments acting on the ends of the member in the local coordinate system: Mx′ij, My′ij, Mz′ij, Mx′ji, My′ji, Mz′ji. Among them, My′ij, Mx′ji are torques; My ij, Mx′ji and Mz′ij, Mz′ji are bending moments about the y-axis and the z-axis, respectively. In addition, the above-mentioned angular displacement and couple moment vector are positive in the positive direction of the coordinate axis according to the right-handed helix rule [3].

Figure 9: Member end angular displacement in local coordinate system.

Figure 10: Member end moment in local coordinate system.

The stiffness matrix of the member ij in the local coordinate system is:

Marked as: M′=K′θ′;

Let the relationship between the local coordinate system and the global coordinate system be:

The conversion matrix from the local coordinate system to the global coordinate system is:

In the global coordinate system, there are: M=K; M=TM; θ=Tθ′

So, M=TK′TTθ; K=TK′TT

The global stiffness matrix is:

When considering the distribution coefficient, it is not necessary to consider the coupling between the bending moments at the member end, but only the diagonal elements (rotational stiffness) of the stiffness matrix K [4].

Therefore, the rotational stiffness of the member in the global coordinate system is:

The distribution coefficient is:

The distribution coefficients include bending moment distribution coefficients and torsional distribution coefficients. There are three groups of distribution coefficients at each rigid node, namely one group around the x-axis, one group around the y-axis, and one group around the z-axis. The sum of the distribution coefficients in each group is 1.

#### Out-plane Stiffness Kout (Three-dimensional calculation)

AECD is a CFST X-pillar with intersection B. Points C and E are sliding bearings. Points A and D are fixed bearings. A rigid body CE is connected between points C and E at the top of the pillar. The length of the CD is L. CB/CD=λ. The load of the pillar can be equivalent as shown in Figures 11 and 12.

Figure 11: CFST X-pillar out-of-plane force model.

Figure 12: CFST X-pillar out-of-plane force equivalent model.

Due to the symmetrical structure of the CFST X-pillar, only half of the CFST X-pillars have been studied. As shown in the Figure 13, AB-BC is a CFST X-pillar calculation structure, in which point A is a fixed branch, and point C is a slip support that is movable only in the y-direction [5].

Figure 13: CFST X-pillar out-of-plane force equivalent calculation model.

Define the local coordinate system of the AB member: Take the AB direction as the x′ direction, and the direction perpendicular to x′ to the left as the y′ direction in the plane of the CFST X-pillar. And define the z′ direction according to the right-handed spiral rule [6].

Define the local coordinate system of the BC member: Take the BC direction as the x′′ direction, and the direction perpendicular to x′′ to the left as the y′′ direction in the plane of the CFST X-pillar. And define the z′′ direction according to the right-handed spiral rule (Figure 14).

Figure 14: Global coordinate system and local coordinate system.

The relationship between the local coordinate system and the global coordinate system can be obtained through geometric relations [7].

Using the previous torque distribution method, the internal force expression of the CFST X-pillar can be calculated. Finally, CFST X-pillar top displacement Δout can be obtained by the method of graph multiplication under in-plane force Pout.

among them,

#### Horizontal Load Distribution Calculation

The in-plane stiffness Kin is along the x-direction and the outof- plane stiffness Kout is along the y-direction (Figure 15) [8]. The s-direction is an arbitrary direction, the angle between which and the y-direction is φ. When unit “1” displacement occurs in the s-direction, that is, Δ=1, it can be decomposed into displacements in x and y directions.

Figure 15: CFST X-pillar top view (stiffness).

Δ = sinϕ , y Δ = cosϕ

The force generated along the s-direction under displacement is:

Suppose that there are n CFST X-pillars, and the stiffness of each CFST X-pillar is denoted by Ksi (i=1,2,…..,n), which is numbered counterclockwise.

The overall lateral stiffness of all CFST X-pillars is

When a horizontal force P acts vertically on the ring beam above the first CFST X-pillar, the horizontal force assigned by the i-th CFST X-pillar is given as (Figure 16):

Figure 16: CFST X-pillar top view (horizontal load).

Decomposes Pi into an in-plane load Pin and an out-of-plane load Pout.

#### CFST X-pillar Internal Force Expression under Horizontal Load

Here only the internal forces of the BC rod end are given:

• Under in-plane load Pin

• Under out-of-plane load Pout

Among them,

#### CFST X-pillar Internal Force Expression under Gravity Load

The load of the pillar can be equivalent as shown in Figure 17.

Figure 17: Calculation model of CFST X-pillar under gravity load.

Here only the internal forces of the BC rod end are given:

Among them,

#### X-pillar Total Internal Force Expression

Previously, the internal forces of the CFST X-pillar under horizontal loads (including in-plane and out-of-plane) and gravity loads have been calculated. They are MPin, MPout and MG respectively. And the total internal force expression is M= MPin+MPout+MG [9].

Here only the internal forces of the BC rod end are given:

#### Specific Project Optimization Examples

The cooling tower bottom radius is 70 meters and the height is 173 meters. The pillars used are concrete-filled steel tube X pillars. The steel pipe is Q420 steel. The concrete strength grade is C60 and the shell concrete grade is C60 (Figure 18).

Figure 18: CFST X-pillar configuration parameters.

The height of the CFST X-pillar is H=30 m, then the length of the corresponding CFST X-pillar is . And the number of CFST X-pillars is n=44. The equivalent static wind load P=6155 kN, and the total gravity load of the ring beam and superstructure is G=419832 kN. The outer diameter of the CFST X-pillar section is D=1000 mm, the thickness of the steel tube is t=20 mm, and the combined elastic modulus Esc=5.245×104 N/mm2 [10].

• Optimization result of α

It can be seen that the stress reaches the minimum when α=76°. Considering the calculation error, the optimal value range of cross angle α is 17°~23° (Figures 19 and 20).

Figure 19: Cross angle α-stress relationship diagram.

Figure 20: Cross angle α-stress relationship envelope diagram.

• Optimization result of β

It can be seen that the stress reaches the minimum when β=76°. Considering the calculation error, the optimal value range of slope angle β is 73°~79° (Figures 21 and 22).

Figure 21: Slope angle β-stress relationship diagram.

Figure 22: Slope angle β-stress relationship envelope diagram.

• Optimization result of λ

It can be seen that the stress reaches the minimum when λ=0.32. Considering the calculation error, the optimal value range of determining the node ratio λ is 0.29~0.35 (Figures 23 and 24).

Figure 23: Node ratio λ-stress relationship diagram.

Figure 24: Node ratio λ-stress relationship envelope diagram.

#### Conclusion

The position of the intersection point is better above the middle of CFST pillar, the range of ratio λ is 0.2~0.35; the optimal value of cross angle α is 18°~22°; the optimal range of CFST pillar tilt angle β is 74°~78°.

#### References

Citation: Zha X, Xu P, Chen D (2018) Theoretical Study and Parameter Optimization of CFST Pillar in Cooling Tower. J Steel Struct Constr 4: 148. DOI: 10.4172/2472-0437.1000148

Copyright: © 2018 Zha X, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Select your language of interest to view the total content in your interested language

##### Recommended Journals
Viewmore
###### Article Usage
• Total views: 169
• [From(publication date): 0-0 - Dec 19, 2018]
• Breakdown by view type
• HTML page views: 149