Abstract

A safety assessment of a chimney under the action of wind and earthquake is performed on one of the reinforced concrete chimneys in the Jinchuan company in China. The assessment is based on the linear filtering method, and wind loads of several heights on one field are simulated by MATLAB. The displacement of the structure under the action of wind load and earthquake is analyzed using maximum value superposition, and the stress on the chimney is calculated in four cases using the equal curvature criterion method. The time-history responses of the tall chimney structure under linear and nonlinear elastic conditions are discussed. The results show that for a chimney under wind load or frequent earthquakes, the displacement limit can meet the code requirements; however, under the combined action of frequent earthquakes and wind load, the structural displacement is bigger than other conditions. When the chimney is acted upon by a wind load, frequent earthquake action, or the combined action of wind load and frequent earthquake action, the maximum compressive stress, maximum tensile stress, and maximum shear stress are concentrated at the bottom of the chimney. The ultimate bearing capacity of a reinforced concrete chimney under three conditions can meet the original design requirements, but the maximum compressive stress at the bottom approaches the designed tensile limit value under the action of rare earthquakes; under these conditions, the bottom of the chimney may experience partial damage.

1. Introduction

Reinforced concrete (RC) chimneys are an indispensable structure in industrial buildings, and their structural safety directly determines the safety of surrounding buildings and people. Earthquakes and wind are the main factors to consider in the design of RC chimneys. The earthquake dynamic responses of chimneys are different from those of general high-rise buildings, as they are irregular and unpredictable. Seismic action is the primary factor in chimney design. With increasing height, the dynamic response under earthquakes and wind becomes more intense and often causes fatal damage. Therefore, determining the dynamic response and assessing the safety of an RC chimney under wind and earthquakes is an urgent task.

Significant research has been carried out regarding the dynamic response and safety evaluation of RC chimneys under earthquakes and wind. Liu and He [1] analyzed the response of chimneys under the combined action of horizontal and vertical earthquake motion. Pan et al. [24] drew the conclusion that the bending moment and shear of high-rise structures obviously increased under the combined action of bidirectional seismic waves. Liu et al. [5] determined that the vertical seismic response of a high-rise structure had obvious regularity. Chen and Ren [6] carried out a simulation analysis of an RC chimney with engineering examples. Liu et al. [7] performed calculations and an analysis of the transient responses of earthquakes on tall chimney structures based on HHT transform and harmonic earthquake response theory. Shiau and Yang [8] proposed applying the vibration mode decomposition response spectrum method to the structure of concrete chimneys. Wilson [9, 10] concluded that the RC chimney possessed some ductility, and the seismic response induced by high modes was significant under complicated earthquakes with long periods. Li et al. [11] obtained actual data by measuring Taipei Building 101 and compared these data with dynamic characteristic data. Feng et al. [12] conducted an analysis of the wind-induced response on three typical high-rise buildings by wind tunnel tests and the random vibration theory method. Li and Lu [13] simulated the fluctuating wind pressure time-history of chimneys of various heights. Li and Li [14] proposed an overall optimization method combined with wind tunnel testing and the finite element method (FEM). Minghini et al. [15] proposed an accurate finite element model built by incorporating suitable tension and compression. Yang et al. [16] studied the reliability of chimney structures under the current standards requirements for long-term wind action, corrosion, hot action, lower-level construction, and lower design standards. Altunisik [17] presents a comparison of the earthquake behavior of RC minarets using fiber-reinforced polymer (FRP) composite. Karaca et al. [18] carried out the dynamic responses of chimneys after FRP strengthening. Zhou et al. [19] revealed that the RC chimney structure had considerable ductility and could withstand a strong earthquake with some structural damage using fragility analysis. Gorski [20] investigated the influence of the three components mentioned above on the total cross-wind response of the chimney. Iban et al. [21] developed a numerical technique for fluid-structure interaction. Cohen et al. [22] studied the seismic vulnerability of an elevated steel tank within a reinforced concrete chimney. Bońkowski et al. [23] analyzed the time-history response of a slender tower under translational-rocking seismic excitation.

In summary, the current research is mainly focused on the design and theoretical analysis of chimneys, and only seismic action or wind load is analyzed. There are no relevant reports on the dynamic response and safety evaluation of existing concrete chimneys under earthquakes and wind loads. Therefore, in this paper, the dynamic response of existing concrete chimneys under earthquakes and wind load is analyzed, and the safety performance of a chimney is evaluated.

2. Wind Vibration Response

For tall structures, the horizontal component of wind can make the structure vibrate in all directions, and the vertical component affects only the axial force of the structure itself. It can be seen that the horizontal component of wind significantly influences the structure. Therefore, this paper analyzes the dynamic response of the horizontal component of wind.

For any height under the action of wind, there will be forces in several directions on the structural surface. To more intuitively and effectively express the effect of wind speed on the structure, wind speed is usually transformed into wind pressure, and surface wind pressure is divided along the surface into the three forces: downward wind force, permanent wind force, and torsional moment. The three forces have the following effects on the structure:(1)The forces along the wind direction include mean wind and fluctuating wind. As a random variable vector, fluctuating wind will cause the vibration of the structure along the wind direction.(2)The horizontal wind can cause the structure itself to resonate. However, due to the wind vortex, the forces on the structure vary and are related to both the structure and the Reynolds number.(3)When the wind causes the eccentric torque of the structure, the structure will experience torsional vibration. If the structure itself cannot be restored after vibration, the structure will be unstable. The downward wind force usually causes structural vibration along the wind direction. For a relatively symmetrical structure, although the vortexes caused by the lateral effect exist in asymmetric forms while wind torque exists, the impact of vortexes is very small. Thus, wind torque is generally considered when the structure is complex and asymmetrical. Therefore, the chimney studied in this paper is used to study only the wind load response in the direction of the wind.

For tall structures [24], vibrations in all directions occur due to the horizontal component of wind, and the axial force of the structure itself is affected by the vertical component. Thus, dynamic response analysis is carried out for the horizontal component of wind.

Wind fluctuations are simulated using the linear filtering method [25], and the wind velocity time-history curve is obtained by MATLAB. The Kaimal spectrum is selected as the power spectral density function as follows:where is the power spectrum density of the Kaimal spectrum, is the frequency, is the frictional velocity, is a dimensionless coordinate, is the mean wind velocity, and is the structure height. After significant checking and testing, the average wind speed fits the formula . In this formula, , where is known as the Kaman constant, and is the surface roughness.

When multidimensional wind speed time-history is simulated by the linear filtering method, the related random wind processes are written as is calculated as

The mean value of is 0. follows a normal distribution random process with a given covariance, . is the regression coefficient matrix.

According to equation (3), the simulation of wind speed time-history is the solution of and .

2.1. Solution of

According to random vibration theory, the relationship between the power spectrum density and the correlation function follows the Wiener–Khinchin formula, that is,where is the cross-power spectrum density between points .

One space point with covariances of and can be expressed as

Because and follow the stationary random wind process with a mean value of 0 and covariance is a function of time difference only, equation (3) can be rewritten as

Both sides of equation (3) are multiplied by , and the mean value of both sides is calculated. Considering that the mean value of is 0, the relationship between covariance and the regression coefficient can be expressed aswhere

The regression coefficient matrix can be obtained by solving the linear equation group given by equation (7).

2.2. Solution of

Both sides of equation (7) are multiplied by , that is,

According to equation (9), is obtained and analyzed by the Cholesky method, that is,where . consists of independent random vectors, which has a mean value of zero and unit variance.

2.3. Generation of Random Wind Speed Time-History

The discretization of equation (9) by a time interval is carried out with and . Therefore, we can get a discrete fluctuating wind velocity time-history vector with spatial correlation time interval .

3. Earthquake Response

The modal analysis [26] of the chimney is performed by combining the vibration superposition method [27] and FEM with the lowest mass participation factor and some design requirements.

The dynamic equilibrium equation of the structural system can be expressed as a group of two-order differential equations :where represents the space vector, which is unrelated to time and represents the -th function related to time.

According to the separation variable method and assuming that the value of equation (11) can be expressed aswhere is a matrix, including the spatial vector of nontime function, can be expressed as a vector containing time function.

Taking the derivative to equation (12), there is

According to the stiffness and quality of the orthogonal condition, there iswhere is a diagonal matrix with a diagonal term of and is a unit diagonal matrix.

Substituting equation (11) into equations (12)–(14), there iswhere ( is the modal participation coefficient of load function ) and is the damping matrix. To obtain the solution, is defined as the diagonal terms of modal damping, that is,

4. Safety Assessment

4.1. Engineering Situation

The height of the RC chimney used to manufacture nickel carbonyl is 75 m. The concrete strength is 35 MPa. The bottom diameter is 5.904 m, and the top diameter is 2.082 m. The bottom of the chimney outer lining is concrete with a thickness of 250 mm, and the lining is a brick masonry insulation layer which has a thickness of 300 mm. Partial cracking occurred at the top of the chimney, but the other sections were in good condition. Thus, the effect of temperature on the structure can be excluded. The vertical and local cracks are shown in Figure 1.

There are local fractures and horizontal cracks at the top of the chimney, and significant anomalies were not found in the other sections. Additionally, vertical cracks in the surface were not found. The cause of the fractures may be the action of earthquake or wind load or the combined action of earthquake and wind load.

4.2. Analysis of the Linear Elastic Calculation Results
4.2.1. Displacement Analysis

(1) Boundary Condition. The chimney structure is only fixed at the bottom, and so the boundary conditions are as follows: the bottom is fixed and the top is free.

(2) Wind Load Action. The height of the calculation model chimney is 75 m. The concrete strength is 35 MPa. A time-history analysis with an elastic range is performed. The acceleration and displacement time-history curves at 25 m, 50 m, and 75 m are shown in Figures 24, respectively.

From Figures 2 to 4, with the increase of height, the acceleration and the displacement of the chimney obviously increase. The extreme displacement under several heights appears at 74.6 s. For both the acceleration curve and the displacement curve, the response curve of each node in the chimney is similar, but the amplitude is different. The accelerations and displacements of these extreme moments are consistent. The wind pressure acceleration is relatively stable and durative. The amplitude of the displacement response of the chimney under the wind load is obvious along the height, and it can be seen from the above figure that the amplitude of the elongated protrusions at the top of the building increases sharply under the wind load. The vertex displacement is large, which conforms to the characteristics of the flexible structure. Many displacement calculations under wind load appeared in the existing literature. In consideration of the requirement of the displacement limit, the displacement limit value of the RC chimneys of various heights under wind load did not exceed the limit value of H/100. The results show that the increasing displacement amplitudes of the upper part of the chimney are larger than those of the lower part, and the ratio of the displacement extreme value to H/100 increases with increasing altitude.

(3) Earthquake Response. Considering both frequent earthquakes and rare earthquakes, according to the building seismic design specifications, horizontal earthquake influence coefficients are as shown in Table 1.

The three main points of seismic wave selection are the characteristic and peak of the continuous time spectrum. The El-Centro wave is selected (peak value 294.0 cm/s2), and the sampling period is 0.02s. The degree of seismic fortification is intensity 8 (0.15 g). The damping ratio is 5%. The seismic wave diagram can be seen in Figure 5. The extreme displacement values of several heights were calculated and can be seen in Table 2. The displacement amplitude during frequent earthquakes and rare earthquakes can be seen in Tables 3 and 4, respectively.

From Table 2 to Table 4, the three moments with extreme values are 5.96 s, 6.00 s, and 6.02 s. With the increase of chimney height, the time of displacement extreme value gradually decreases. Moreover, the variation of the maximum displacement response of the chimney under earthquakes is obvious. The displacement extreme value of the RC chimney meets the requirements under earthquake action.

(4) The Combined Response to Frequent Earthquakes and Wind Load. The crack of the RC chimney used for the manufacture of nickel carbonyl is mainly caused by the combined action of earthquakes and wind load. Given that the acting time of the wind load is longer than that of the earthquake, the analysis of the combined action of earthquake and wind load is carried out by the method of superimposing the maximum value. The calculation results are shown in Table 5.

From Table 5, the most displacement extreme value of the RC chimney is 533 mm at 75 m under the combined action of earthquake and wind load, and the structural displacement does not exceed the specified limit. The extreme displacement value of the middle and upper parts of the chimney also satisfies the displacement limit. According to Table 5, the displacement of the RC chimney under earthquake and wind load action is larger than wind and earthquake action alone.

Given that the extreme displacement value of the middle and upper parts of the chimney under the combined action of frequent earthquakes and wind load did not approach the displacement limit and the probability of suffering rare earthquakes is small, the extreme displacement value will not be calculated under the combined action of rare earthquakes and wind load in this paper.

4.2.2. Stress Analysis

(1) Wind Load Action. Basic wind pressure equals 0.35 kN/m2. The ground roughness is Grade B. Wind load is one of the main loads that influence flexible high-rise structures, and it is necessary for the design and confirmation of the wind load. The calculation and analysis of wind load can be seen in Tables 6 and 7.

Based on Table 6, the maximum stress appears at the bottom of various sections, and it meets the design value of the shear bearing capacity. Thus, the ability of the chimney to resist shearing is adequate. From Table 7, the peak value of the bending moment is 5893.57 kN·m under wind load alone, but it is almost half the standard value of the bending moment under earthquake action. Thus, the internal force of the chimney under earthquake action is almost twice that under the wind load, which still meets the design requirements. The maximum compressive stress and the maximum tensile stress are much smaller than the compression and tensile strengths of concrete.

(2) Earthquake Response. The internal force of chimneys with heights of 25 m, 50 m, and 75 m can be seen in Figures 68 under the frequent earthquake response at 6.02s, 6.00 s, and 5.96 s, respectively.

From Figures 6 to 8, the maximum stress zone gradually moves toward the bottom of the chimney as time decreases, and the maximum value decreases when the minimum value increases. The greatest dynamic response of the chimney occurs at the bottom 10 m and the minimum occurs at the top 10 m. The maximum stress value does not exceed the structural admissible stress range; therefore, the structure meets the stiffness requirements. The stress value of the chimney is calculated by the equal curvature criterion method. The maximum stress value and the maximum shear stress value of the chimney under the action of frequent earthquakes are shown in Tables 8 and 9.

From Table 8, the coefficient of the bending section of the structure is influenced by the height of the chimney, and the maximum value appears at the bottom. As height increases, the coefficient gradually decreases. The maximum compressive stress and the tensile stress both appear at the bottom of the chimney, and the compression stress is greater than the tensile stress at various heights and at the bottom of the chimney. From Table 9, shear standard values change at various sections, and the maximum value occurs at the middle and bottom of different sections. The maximum shear stress also appears at the bottom of the chimney, but with the height change, the numerical variation amplitude is not very large.

Given that the maximum stress at the bottom of various sections is far smaller than the designed value of shear capacity during frequent earthquakes, only the maximum stress values under the action of rare earthquakes will be calculated, as shown in Table 10.

From Table 10, the most disadvantaged position is the bottom of the chimney. The maximum compressive stress at the bottom is 1865.22 kPa, and maximum tensile stress is 1607.95 kPa. The maximum compressive stress at the bottom approaches the limit tensile design value, which shows that the chimney may experience partial damage at its bottom during rare earthquakes.

(3) Combined Response to Earthquake and Wind Load. Given that the probability of suffering rare earthquakes is small, the maximum stress will not be calculated under the combined action of rare earthquakes and wind load. The combined values of frequent earthquake action and wind load are shown in Table 11.

From Table 11, the most disadvantaged position of the chimney is the bottom. The maximum compressive stress at the bottom is 3569.22 kPa, and maximum tensile stress is 695.1634 kPa. The stress under the combined response to earthquake and wind load is significantly larger than the responses to either wind load or earthquake action alone. The maximum bearing capacity occurs at the structure bottom under the combined response to earthquake and wind load, and the ultimate bearing capacity of the RC chimney meets the requirements under the combined action.

4.3. Analysis of the Nonlinear Elastic Calculation Results

When considering materials’ nonlinearity and P-∆ effect, the displacement and acceleration amplitudes under the action of wind load, earthquake, or wind and earthquake are shown in Table 12.

From Table 12, with nonlinear analysis, the displacement and acceleration amplitudes under wind load are less than earthquake. Even when wind and earthquake are acted on the chimney, the amplitude does not exceed displacement limit in the code. Compared with linear elastic analysis, the nonlinear elastic analysis results, that are displacement and acceleration amplitudes, are slightly bigger except the displacement under wind action.

4.4. Security Assessment

(1) Displacement Evaluation. A large number of relevant studies show that RC displacement does not exceed displacement limit H/100 under wind load alone. The maximum displacement of an RC chimney at 75 m under earthquake action alone meets the displacement limit.

(2) Stress Assessment. Under frequent earthquake action alone, the maximum compressive stress at the bottom of the chimney is −1311.37 kPa, and the maximum tensile stress is 1205.30 kPa. Under wind load action alone, the maximum value of the bending moment is 614.81 kPa. Under rare earthquake action alone, the maximum compressive stress at the bottom of the chimney is 1865.22 kPa, and maximum tensile stress is 1607.95 kPa. Under the combined action of earthquakes and wind load, the maximum compressive stress at the bottom of the chimney is −3569.22 kPa, and the maximum tensile stress is 695.1634 kPa.

In summary, the displacement of the RC chimney meets the specification under the action of wind load or earthquake, and the ultimate bearing capacity of the RC chimney also meets the design requirements (excluding rare earthquake). Therefore, some reinforcement measures should be taken to ensure safety. Through the calculation of structural displacement and internal force under the combined action of wind and earthquake, it is concluded that the internal force and the displacement of the chimney can meet the design requirements in this case.

5. Conclusions

Based on the linear filtering method, the wind load of several heights on one field is simulated by MATLAB, the displacement of the structure under the action of wind load and earthquake is analyzed, and the stresses of the chimney in three cases are calculated. The main conclusions are as follows:(1)With the increase of chimney height, the time at which the extreme value of displacement occurs consistently decreases. The extreme displacement value of an RC chimney meets the requirements under wind or earthquake, and the displacement also meets the specified limit under the combined action of wind and earthquake. According to the calculated results, the response of the structure under wind load is very obvious, and because wind load lasts for a long time, the vibration cycle of the structure itself lengthens. Under this condition, although the displacement of the apex of the structure reaches almost 130 millimeters, the equivalent stress at the time of the extreme displacement value satisfies the requirement of the carrying capacity because of its material properties. Because the chimney construction time is long, after a long period of vibration, the chimney may experience security risks.(2)The extreme value of stress develops progressively from the bottom of the structure to the top of the structure over time. The bearing capacity of the RC chimney meets the design requirements under wind load or earthquake action alone. However, a partial tensile failure at the bottom of the chimney may occur during rare earthquakes. Therefore, we should take reinforcement measures to ensure safety. By calculating the internal force of the structure under the combined response to wind load and earthquakes, the internal force meets the design requirements, and the effect on the displacement response is greatly influenced by the effect of the structure. By comparing the combined response with the stress response, it can be determined that the combined response has a greater impact on the displacement of the structure.(3)It is not scientific or authentic to determine the safety performance of the structure based only on the displacement from the value of maximum tensile stress under rare earthquake conditions. Additionally, the weak spots of the structure indicate that the bottom may experience partial tensile failure, and the displacement at the upper apex may not meet the specified limit. The displacement and stress should be used as the basis for the safety evaluation of the RC chimney in future structure design, detection, and reinforcement.(4)For the results of wind load and earthquake together, the extreme displacement values satisfy the displacement limits at the top of the chimney, but the extreme stress is too large, which is consistent with the location of chimney rupture and crack. Thus, the calculation is true and effective, and there are some hidden dangers in the chimney under the combined action of earthquake and wind load, and the safety performance is unreliable. Therefore, it is necessary to consider both wind load and earthquake interaction in future chimney design.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This paper is part of the project of the National Natural Science Foundation of China (Grant no. 51478212) and part of the Science and Technology Project of the China Railway 12th Bureau Group (Grant no. 14B-3).