Earthquake Response Analysis of Tall Reinforced Concrete Chimneys considering Eccentricity

Jiang, Yingchun; Liu, Tielin; Bai, Yikui

doi:https://doi.org/10.1155/2020/1417969

Shock and Vibration

On this page

Abstract Introduction Conclusions Appendix Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 1417969 | https://doi.org/10.1155/2020/1417969

Earthquake Response Analysis of Tall Reinforced Concrete Chimneys considering Eccentricity

Yingchun Jiang,¹Tielin Liu,²and Yikui Bai³

Academic Editor: Salvatore Caddemi

Received14 Oct 2019

Revised07 Jan 2020

Accepted29 Jan 2020

Published11 Mar 2020

Abstract

A numerical algorithm is presented to analyze earthquake response of tall reinforced concrete (RC) chimneys based on stick multidegree-of-freedom models. The algorithm considers the eccentricity phenomena between spatial discrete nodes and corresponding centroids of investigated lumps. The spatial discrete segments of the chimney are used to construct the investigated lumps. The equations of dynamic equilibrium of the investigated lumps are derived, and the numerical calculation procedure is implemented. Phenomena of eccentricity are studied for 150 m and 210 m RC chimneys. Seismic stresses and effects of vertical ground motion for the two chimneys are also studied. Numerical results show that the tensile and compressive stresses of the seismic control cross sections of the chimneys may increase under the actions of several specific earthquake waves by considering existing eccentricities. The effect of eccentricity on the earthquake responses of tall RC chimney should be considered, and stresses caused by vertical ground motion should not be neglected to obtain accurate earthquake response of chimneys.

1. Introduction

1.1. Literature Review

Many tall reinforced concrete (RC) chimneys have been built in places such as coal-fired power stations and refinery because of the requirements of air pollution control standards. Slender tall chimneys are sensitive to earthquake. Strong earthquake threatens the safety of the chimney and other buildings around it because of the nonuniformity of the diameter along the height of the chimney and the variation in dead weights at different heights. Research and calculations of earthquake response are necessary for the designs of tall RC chimneys under the action of earthquakes.

Destruction and collapse of tall chimneys bring huge economic losses. Research on seismic damage of chimneys indicates that the destructive effect of vertical earthquake on chimneys should be considered [1]. Liu and He [2] analyzed the seismic response of a chimney under the action of combined horizontal and vertical ground motions. They concluded that vertical stress cannot be neglected completely although horizontal seismic stress is the controlling part. Wang and Liao [3] developed a calculation technique for analyzing the vertical earthquake response of a chimney based on multiple reflections of traveling waves in the chimney. Liu et al. [4] conducted earthquake transient response of multidegree-of-freedom (MDOF) systems by the model superposition method. Ding et al. [5] conducted a vertical impact to the bottoms of chimney models. The test results showed that only one circle crack appeared in the chimney models, and the crack had relations with the ratio of input wave frequency to the natural frequency of the chimney model. Chen et al. [6] used a shaking table to study the fracture phenomenon of a chimney under the action of vertical earthquake. They found that vertical seismic action should be considered in the structural design of chimneys.

A 115 m high RC chimney collapsed on August 17, 1999, during the Izmit earthquake in Turkey [7–10]. Huang et al. [7, 8] investigated a collapsed chimney with methods of linear dynamic response spectrum and nonlinear static capacity spectrum. Akinci [9] used shell elements to analyze a collapsed chimney. The analysis results indicated that low ductility and inadequate capacity due to a large duct were the two significant factors for the collapsed chimney. Huang and Gould [10] further proposed a new three-dimensional (3-D) pushover analysis method to study the dynamic response of the stack due to an earthquake motion recorded at a nearby site. Lopes et al. [11] constructed a realistic numerical model to evaluate the seismic vulnerability of a chimney and simulated its damaged state. Minghini et al. [12] analyzed a brickwork chimney and a shear failure mechanism of the upper part of the chimney damaged by the 2012 Emilia earthquake. Zhou et al. selected an existing 240 m tall RC chimney as the research object and analyzed the seismic vulnerability curve of the chimney by considering structural damage [13], structural nonlinearity [14], and multidimensional earthquake [15] and evaluated and analyzed the vulnerability of the chimney structure.

The modal frequency and mode shape are the main dynamic characteristics of chimneys. Chen et al. [16] used the Adomian decomposition method to decompose dimensionless modal shape and calculated the dynamic characteristics of a high-rise chimney. Sancibrian et al. [17] used the continuum structure system to analyze the dynamic characteristics of a chimney, including natural frequencies and mode shapes. Tatara and Ratajewicz [18] evaluated the dynamic characteristics of a typical RC industrial chimney that belongs to a heating plant and determined the lowest vibration frequency of the structure.

Soil-structure interaction (SSI) has a great influence on the seismic performance of tall chimney structures [19, 20]. Halabian and Kabiri [21] studied the effect of foundation stiffness of different types of RC chimneys on the ductility reduction coefficient. A series of inelastic and elastic time history analyses up to the collapse status for a large variety of RC chimneys subjected to five scaled artificial records were conducted. Jisha et al. [22] analyzed 3-D SSI of tall RC chimneys with annular raft foundation and evaluated the effect of SSI by changing the stiffness of supporting soil and foundation [23]. Zhou et al. [24, 25] studied the effect of SSI on the seismic collapse resistance of chimney structure for a 240 m high RC chimney.

Time period is an important approach of evaluating the structural response of a structure subjected to large-scale earthquake vibrations. The seismic response of chimney structures uses time history analysis considering earthquake, which is a long-duration earthquake impulse [19]. Sarkar et al. [26] studied the effects of different base conditions on the design value of shear force and bending moment of different sections of a 275 m tall RCC multiflue chimney by using a response spectrum method. The spectrum compatibility time history of specific engineering site was considered, and the time history dynamic analysis of the chimney structure was conducted in detail. Wilson [27] studied the inelastic response of 10 RC chimneys ranging from 115 m to 301 m subject to earthquake excitation. Liu et al. [28] studied the structural dynamic performance and elastic-plastic seismic response of a 240 m high special-shaped RC chimney. The weak position of the chimney was discussed by analyzing its interlayer displacement angle.

1.2. Motivation

When the stick MDOF structural model is used to study the dynamic response of tall chimney structures, the chimney structure is usually discretized into beam elements of a certain length. Thin beam elements are good, but they are not sufficiently thin in practical analysis. The cross section area of a chimney structure is variable at different heights. Therefore, even if the discrete length of space is equal, the vertical eccentricity of mass will occur because the mass of each element will be uneven. The vertical eccentricity of mass must be considered to improve the accuracy and validity of dynamic response analysis.

A numerical algorithm is presented in this study to analyze the earthquake response of tall RC chimneys. The chimney is spatially discretized into several segments of finite length, and the investigated lumps are then constructed by using these segments. The calculating formulae of axial force, shear force, and bending moment acting on the median cross section of each segment are derived. These interior forces are associated with end displacement and angle. Not only the inertial forces of the investigated lumps are considered but also the rotary inertias of investigated lumps. Existing eccentricity phenomena between spatial discrete nodes and corresponding centroids of investigated lumps are considered. Dynamic equations and relations are applied alternately in the time domain to calculate the seismic response of the chimney. Two RC chimneys of 150 m and 210 m height are studied. The numerical results show that the effects of such eccentricities on earthquake response should not be neglected for tall chimneys, as well as the vertical ground motion.

2. Numerical Algorithm

2.1. Governing Equations of Investigated Lumps

A chimney with variable cross section can be discretized into many segments of different cross section areas along the vertical direction in practical calculation based on the stick MDOF model. These segments are used to construct the investigated lumps. Figure 1 shows a typical investigated lump (the nodal numbering is also used to describe the investigated lump) consisting of half the upper segment and half the lower segment around the discrete node . Sign is the centroid of the investigated lump. The centroid of the investigated lump is commonly not at the same point of the corresponding discrete node for the chimney with variable cross-sectional area. Here, the distance between the discrete node and the centroid of the investigated lump within the same investigated lump is called eccentricity.

The vertical equilibrium equation for the centroid of investigated lump iswhere is the mass of investigated lump , is the vertical acceleration of centroid C, and and are the axial forces acting on investigated lump .

The transverse equilibrium equation for the centroid of investigated lump iswhere is the horizontal acceleration of the centroid of investigated lump and and are the shear forces acting on investigated lump .

On the basis of mechanics theory, the rotational equilibrium equation should be established about the horizontal axis (-axis) through centroid C of investigated lump , that is,where is the moment of inertia of investigated lump about the -axis through its centroid; and are the angular acceleration and the eccentricity of investigated lump , respectively; and are the bending moments and the length of segment , respectively; and is the transverse displacement of node .

Figure 2 shows that the following relation holds:where is the horizontal acceleration of node . In Figure 2, is the vertical acceleration of node .

2.2. Forces and Bending Moment on the Central Cross Section of the Segment

Figure 3 shows the free-body diagrams of the upper and lower halves of segment by an artificial cut at the central section of segment . In Figure 3, and are the axial compressive forces at the bottom and top ends of segment , respectively, because of its vertical deformation; and are the shear forces; and are the bending moments; is the axial force coming from the dead load of the structure above node ; and is the dead load above node .

Following the beam theory considering shear deformation and axial force, we havewhere , , , and . Subscript denotes the l^th segment. is the coefficient of shear deformation in the transverse direction. , , and are the bending stiffness, shear stiffness, and tensile/compressive stiffness of segment , respectively. , where is the axial force coming from the dead weight of the chimney above node . is the cross-sectional area of segment . and are the vertical displacement and rotational angle of node , respectively. and are the vertical displacement and rotational angle of node , respectively.

Subtracting equation (8) from equation (7) yields

Given the equilibrium of the upper half of segment , we have

Given the equilibrium of the lower half of segment , we havewhere , , and are the central axial force, central shear force, and central bending moment acting on the central section of segment , respectively.

We substitute equations (10) and (12) into equation (9), and we obtain the central bending moment of segment , shown as follows:where .

2.3. Implementation of Algorithm

Horizontal and vertical earthquake waves are inputted from the foundation of the tall RC chimney. The algorithm description diagram is shown in Figure 4, and the calculating procedure is as follows. Step 1. Calculate the angular acceleration of investigated lump at time using equation (3), the vertical acceleration at time using equation (1), at time by considering being nearly equal to , the horizontal acceleration at time using equation (2), and using equation (4). Step 2. Calculate angle , vertical displacement , and horizontal displacement for investigated lump at time by integration in the time domain. Step 3. Calculate axial force at time using equation (5) and then obtain central axial force using the formula . Central shear force and central bending moment are calculated using equations (11) and (13), respectively, at time . Step 4. Calculate , , and at time as performed in Step 1.

The circulation from Step 1 to Step 4 is implemented in the time domain by a computer program. The numerical algorithm is provided to calculate the earthquake response of the tall RC chimney, considering eccentricities ().

Wilson method is used to solve the conventional dynamic equations of equilibrium in which eccentricities are not considered. Their results are called the conventional results in this study.

Here, the rotation of the foundation is not considered.

3. Stability Condition

3.1. Corresponding Relations

For simplicity in deriving the stability condition of the algorithm, we assume that , , and for all segments. Axial force is omitted in . Eccentric phenomena are neglected in equation (3); and and for all investigated lumps. On the basis of these assumptions, vertical governing equations are decoupled from the governing equations of transverse and rotational deformations. Thus, we only study the coupling equations of transverse and rotational deformations to provide the stability condition for the algorithm presented in this study.

On the basis of the above assumption, the corresponding relations between the coefficients used in Timoshenko beam theory and those used in the proposed algorithm are as follows (see Appendix A):

3.2. Stability Condition Based on Timoshenko Beam Model

The flexural wave equation of Timoshenko beam iswhere , , and .

Substituting a flexural wave into the finite difference formula of Equation (15), we obtain the following equation (see equation (A.12) in Appendix A):

Solving equation (16) with as unknown yields the following:

Given , we approximately have the following equation (see equation (A.16) in Appendix A):

3.3. Stability Condition for the Proposed Algorithm

The stability condition used in the algorithm presented in this paper can be provided by using equation (18) and the corresponding relations (14) as follows:

4. Method Validation

We study bending wave propagation of a beam to verify the effectiveness of the new algorithm. The beam with equal cross-sectional area is designed to avoid eccentricity. The result given by the finite difference method [29] is used to test the proposed numerical algorithm.

The elastic modulus of the beam is E = 209 GPa, density is ρ = 8000 kg/m³, Poisson’s ratio is = 0.3, and shear cross-sectional correction factor is k′ = 0.886. The length of the beam is X = 6.35 m, and the radius is 0.127 m. The number of discrete segments is 500 in the length direction for the uniform beam. On the basis of the calculation parameters of the beam, the time step of the beam satisfying the stability condition of the method is obtained by equation (19), and the calculation time step is .

The length from the left end of the beam is x, and the bending moment at the section is M. After normalization, the time is , the length is , and the bending moment is .

The inclination moment of the left end of the beam is

Figure 5 shows the comparison curves of the bending moments at different cross sections of the beam at .

The calculation results of the new numerical method are nearly the same as those of the finite difference method. The comparison results show that the presented algorithm is valid for calculating the earthquake responses of tall RC chimney.

5. Earthquake Responses of Tall RC Chimneys

5.1. Earthquake Records Used

Eleven earthquake acceleration records on four types of soil are used as earthquake load to calculate the structural earthquake response. Four types of soil are denoted by rock, soft rock, stiff soil, and soft soil. Soil type is characterized by shear wave velocity. The corresponding earthquake records are selected from the Pacific Earthquake Engineering Research Center website [30]. Table 1 shows the 11 earthquake records.

According to the peak value of design earthquake acceleration of 8 degrees seismic intensity in the literature [31], the peaks of the horizontal earthquake accelerations are adjusted to 0.2 g in practical calculations. The peaks of the vertical earthquake accelerations are adjusted to 0.13 g. The damping is omitted.

5.2. Two Models of RC Chimneys

Two RC chimneys [32] of 150 m and 210 m height are studied. Figure 6 shows the two models of RC chimneys. (1) The concrete strength grade of the 150 m high RC chimney is C25, and the steel strength grade is HRB335. The outer radius of the bottom is 8.46 m, and the wall thickness is 0.4 m. The outer radius of the top is 3.0 m, and the wall thickness is 0.16 m. The gradient c is 0.1 below 25 m, 0.03 ranging from 25 m to 75 m, and 0.015 ranging from 75 m to 150 m. Table 2 shows the parameters of the 150 m high RC chimney. (2) The concrete strength grade of the 210 m high RC chimney is C35, and the steel strength grade is HRB335. The outer radius of the bottom is 10.6 m, and the wall thickness is 0.6 m. The outer radius of the top is 3.85 m, and the wall thickness is 0.18 m. The gradient c is 0.1 below 30 m, 0.03 ranging from 30 m to 100 m, and 0.015 ranging from 100 m to 210 m. Table 3 shows the parameters of the 210 m high RC chimney. According to the literature [33], the standard value of the compressive strength of C25 concrete is 16.7 N/mm², the design value is 11.9 N/mm², and the elastic modulus is E_C = 2.8 × 10⁴ N/mm². The standard value of the compressive strength of C35 concrete is 23.4 N/mm², the design value is 16.7 N/mm², and the elastic modulus is E_C = 3.15 × 10⁴ N/mm². The elastic modulus of HRB335 reinforcement is E_S = 2.0 × 10⁵ N/mm², the standard value of strength (yield strength) is 335 N/mm², and the design value of strength is 300 N/mm².

(a)

(b)

The discrete segments are all 10 m long in the vertical direction for each of the two chimneys. Each chimney will become part of such calculated model consisting of several cylinders of different cross-sectional areas because an actual chimney structure has variable cross-sectional areas. Thus, the phenomena of eccentricities between the nodes and the corresponding centroids of the investigated lumps appear in the calculation of the earthquake response of the actual chimney. The mass of each investigated lump consists of half the upper segment and half the lower segment around the discrete node. The moment of inertia of each investigated lump about the horizontal axis through its centroid is calculated. On the basis of the calculation parameters of the 150 m and 210 m high chimneys, the time step satisfying the stability condition of the method is obtained by equation (19), and in the 150 m and 210 m high chimneys, respectively. The time step is 1.0 ms.

5.2.1. Effect of Eccentricity

Eleven earthquake waves are used. The effect of eccentricity on normal earthquake stress of the chimney are studied when horizontal and vertical components are inputted simultaneously.

Figure 7 shows the spatial distribution curves of the maximum earthquake normal stresses of different cross sections for the 150 m high chimney under the actions of Landers (No. 1 in Table 1), Loma Prieta (No. 2), Borrego (No. 3), and Chichi (No. 4) waves. Figure 8 shows the comparison curves of the 210 m high chimney. “New” denotes that the phenomena of eccentricities are considered. “Conventional” denotes the conventional results, in which the phenomena of eccentricities are neglected, for the earthquake responses of the real RC chimneys. Dashed lines (“tensile”) denote tensile stresses, whereas solid lines (“press”) denote compressive stresses.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

A comparison of the red solid lines and the blue solid lines in Figures 7 and 8 shows that the increased phenomena of peak values of the normal compressive stresses appear when eccentricities are considered.

A comparison of the red dashed lines and the blue dashed lines in Figures 7 and 8 shows that the increased phenomena of peak values of the normal tensile stresses appear when eccentricities are considered.

Table 4 shows the details of the increase in percentage of maximum earthquake normal tensile/compressive stresses when phenomena of eccentricities are considered.

Table 5 shows maximum earthquake compressive stress, mean value, and standard deviation of different sections of the 150 m high chimney under 11 earthquake wave actions when vertical eccentricities are considered. Table 6 shows maximum earthquake compressive stress, mean value, and standard deviation when vertical eccentricities are ignored.

In Table 5, the comparison of the data shows that the maximum earthquake compressive stress obtained under different earthquake waves is significantly different, and the degree of the difference can be expressed by the standard deviation. To quantify the degree of difference, its standard deviation S iswhere n is the number of samples, that is, the number of selected earthquake waves; is normal stress on the central section of segment , that is, the sum of bending stress and axial stress, which is ; is the bending-resistance modulus of the segment ; and is the mean normal stress of the central section.

The right side of Figure 9(a) is the comparison of the mean values of earthquake compressive stresses between Tables 5 and 6 and the corresponding standard deviation. The left side of Figure 9(a) is the comparisons of the mean values of earthquake tensile stresses and the corresponding standard deviation. The data in Tables 7 and 8 (see Appendix B) are provided regarding maximum earthquake tensile stress, mean value, and standard deviation of the 150 m high chimney when considering vertical eccentricities and neglecting vertical eccentricities. In Figure 9(a), the red and pink columns are, respectively, the mean values of earthquake compressive and tensile stresses considering vertical eccentricities. The blue and green columns are, respectively, the mean values of earthquake compressive and tensile stress neglecting vertical eccentricities. The red and pink lines are, respectively, the standard deviations of different earthquake compressive and tensile stresses considering vertical eccentricities. The blue and green lines are the standard deviations of different earthquake compressive and tensile stresses neglecting vertical eccentricities. Figure 9(b) shows the case of the 210 m high chimney (data in Tables 9–12).

(a)

(b)

Figure 9 shows that the earthquake normal stresses increase evidently when considering vertical eccentricities. Significant differences of the maximum earthquake normal stress are observed under different seismic waves actions.

Figure 10 is the boxplot of the maximum earthquake compressive stress distribution at different cross sections of 150 m and 210 m high chimneys under different seismic wave actions when considering eccentricities. They show the mean value, minimum, maximum, and 25%–75% of earthquake compressive stresses. The distribution characteristics of the earthquake compressive stress data are also shown in the figures. The stress distributions near the bottom of the chimney are relatively intensive, whereas those at the middle and upper parts of the chimney are relatively scattered.

(a)

(b)

The results studied above show that significant differences of maximum earthquake normal stress under different earthquake wave actions. The comparisons demonstrate that the conventional results (eccentricities neglected) might be less than their actual tensile/compressive stresses in the seismic control cross section of the chimneys under the actions of several specific earthquake waves. The calculation results show that eccentricities between the spatial discrete nodes and the corresponding centroids of investigated lumps should be considered in earthquake response of tall RC chimneys to obtain accurate calculation results.

5.2.2. Effects of Vertical Earthquake

The phenomena of eccentricities are considered to study the effect of vertical earthquake on earthquake normal stresses. Two types of seismic load are applied for each of the four earthquake waves. In the first case, the combined horizontal and vertical earthquake waves (see red lines in Figures 11 and 12) are inputted. In the second case, only the horizontal earthquake wave (see blue lines in Figures 11 and 12) is inputted.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

For the 150 m high chimney, Figure 11 shows that the peak normal tensile stress increases by 4% for the Landers wave (No. 1) and 10% for the Loma Prieta wave (No. 2) if vertical earthquake waves are considered. The peak normal compressive stress increases by 7% for the Loma Prieta wave (No. 2) and 10% for the Chichi wave (No. 4) if vertical earthquake waves are considered.

For the 210 m high chimney, Figure 12 shows that the peak normal tensile stress increases by 12% for the Landers wave (No. 1) and 9% for the Loma Prieta wave (No. 2) if vertical earthquake waves are considered. The peak normal compressive stress increases by 9% for Landers wave (No. 1) and 13% for Loma Prieta wave (No. 2) if vertical earthquake waves are considered.

Figure 13 shows the spatial comparison curves of the mean values of the maximum normal earthquake stresses for the 150 m and 210 m high chimneys under two types of seismic load. Figure 13(a) is obtained from the average values of Tables 5, 7, 13, and 14. Figure 13(b) is drawn from the mean values of Tables 9, 11, 15, and 16. Figure 13 shows that the normal stress of the chimney can be increased considering the vertical seismic acceleration.

(a)

(b)

The results studied above show that the vertical ground motion should not be neglected, especially for tall RC chimneys built on rock sites or soft rock sites.

6. Conclusions

A numerical algorithm that can consider the existing phenomena of eccentricities is presented to calculate the earthquake responses of tall RC chimneys. The numerical results demonstrate that the tensile and compressive stresses may increase in the seismic control cross sections of the chimneys under the actions of several specific earthquake waves when existing eccentricities are considered. The phenomena of eccentricities should be considered to obtain accurate calculation results of earthquake responses for the seismic designs of tall RC chimneys. The maximum normal stresses of earthquake response appear within the upper half of the chimneys, and their values are considerably larger than those within the lower half of the chimneys. Vertical ground motion should not be neglected, especially for chimneys on rock sites and soft rock sites.

Appendix

A. Corresponding Relations

For simplicity in deriving the stability condition, we assume that , , and for all segments and and for all investigated lumps. Axial force is omitted in , , and . Eccentric phenomena are neglected in equation (9). Let be zero. The vertical governing equations are decoupled from the governing equations of transverse and rotational deformations. Thus, we study the equations coupling transverse and rotational deformations to provide the stability condition for the method in this work.

On the basis of the aforementioned assumption, the corresponding equations between the governing equations of Timoshenko theory and the proposed method can be expressed as follows:where .

On the basis of the aforementioned corresponding relations of equations, the corresponding relations between the coefficients used in the continuum Timoshenko beam theory and those used in the proposed method are

Therefore, we havewhere is the moment of inertia of the cross-sectional area with respect to the neutral axis.

A.1. Stability Condition Based on Timoshenko Beam Model

The flexural wave equation of Timoshenko beam iswhere , , and ; and are the elastic and shear moduli of the beam, respectively; is the mass density; is the cross-sectional area; and is the effective cross-sectional area for the shear deformation of the beam.

The finite difference formula of equation (A.9) is

Assuming a flexural wave propagating in the beam and substituting it into equation (A.10), we obtain

Given and , equation (A.11) can be derived as

Equation (A.12) is a quadratic equation with unknown , and its solution can be obtained as follows:

Given , the stability condition should bewhich can be expressed approximately asThat is,

A.2. Stability Condition for the Proposed Algorithm

The stability condition used in the proposed method can be provided by using equation (A.16) and the corresponding relations (i.e., equations (A.6)–(A.8)) as follows:

B. Partial Calculation Results

Table 7 shows the maximum earthquake tensile stresses, mean values, and standard deviations of different sections of the 150 m high chimney under 11 seismic waves when horizontal and vertical seismic accelerations are simultaneously inputted and vertical eccentricities are considered. The results when vertical eccentricities are ignored are shown in Table 8.

Table 9 shows the maximum earthquake compressive stresses, mean values, and standard deviations of different sections of the 210 m high chimney under 11 seismic waves when horizontal and vertical seismic accelerations are simultaneously inputted and vertical eccentricities are considered. The results when vertical eccentricities are ignored are shown in Table 10.

Table 11 shows the maximum earthquake tensile stresses, mean values, and standard deviations of different sections of the 210 m high chimney under 11 seismic waves when horizontal and vertical seismic accelerations are simultaneously inputted and vertical eccentricities are considered. The results when vertical eccentricities are ignored are shown in Table 12.

In view vertical eccentricities, Table 13 shows the maximum earthquake compressive stresses, mean values, and standard deviations of different sections of the 150 m high chimney under 11 seismic waves when only horizontal seismic acceleration is inputted. Table 14 shows the results of maximum earthquake tensile stresses of the 150 m high chimney.

In view of vertical eccentricities, Table 15 shows the maximum earthquake compressive stresses, mean values, and standard deviations of different sections of 210 m chimney under 11 seismic waves when only horizontal seismic acceleration is inputted. Table 16 shows the results of the maximum earthquake tensile stresses of the 210 m high chimney.

Data Availability

The analysis result data used to support the findings of this study are included within the article. The calculation 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.

Acknowledgments

This study was supported in part by the National Natural Science Foundation of China (11672190 and 61673281), Liaoning Province Doctor Startup Fund (201601112), Liaoning Provincial Key Research and Development Program (2018103007), and Qinghai Provincial Key Research and Development Program (2017-NK-107).

Supplementary Materials

“Parameter data of 150 m reinforced concrete chimney.xlsx” is the RC chimney parameter data of 150 m height. “Parameter data of 210 m reinforced concrete chimney.xlsx” is the RC chimney parameter data of 210 m height. “Eleven kinds of earthquake records” are input data, which are used to study the seismic response of chimney structures. (Supplementary Materials)

References

Z. P. Liao and F. Wang, “Vertical seismic motion input to structures,” Earthquake Engineering and Engineering Vibration, vol. 3, no. 2, pp. 74–88, 1983, in Chinese.
View at: Google Scholar
J. Liu and L. M. He, “Analysis of seismic response of chimneys to combined horizontal and vertical ground motion,” Earthquake Engineering and Engineering Vibration, vol. 3, no. 2, pp. 89–98, 1983, in Chinese.
View at: Google Scholar
F. Wang and Z. P. Liao, “Wave motion response of chimneys to vertical ground motions,” China Civil Engineering Journal, vol. 17, no. 1, pp. 12–26, 1984, in Chinese.
View at: Google Scholar
Q. Liu, R. Z. Zhou, and W. J. Yuan, “Computation of earthquake transient response of skyscraper structures and study of the dynamic properties,” China Civil Engineering Journal, vol. 43, no. 5, pp. 63–69, 2010, in Chinese.
View at: Google Scholar
S. Ding, Y. J. Wang, and S. J. Zhou, “Fracture test of chimney models subjected to vertical shock,” Earthquake Engineering and Engineering Vibration, vol. 19, no. 3, pp. 34–36, 1999, in Chinese.
View at: Google Scholar
J. Y. Chen, J. Zhou, and H. C. Ma, “Study on model test of highrise chimney subjected to vertical seismic action,” Journal of Building Structures, vol. 26, no. 2, pp. 87–93, 2005, in Chinese.
View at: Google Scholar
W. Huang, P. L. Gould, R. Martinez, and G. S. Johnson, “Non-linear analysis of a collapsed reinforced concrete chimney,” Earthquake Engineering & Structural Dynamics, vol. 33, no. 4, pp. 485–498, 2004.
View at: Publisher Site | Google Scholar
P. L. Gould and W. Huang, “Higher mode effects in the nonlinear static analysis of a collapsed chimney,” in Proceedings of the 2006 Structures Congress, vol. 2006, p. 19, St. Louis, MO, USA, May 2006.
View at: Publisher Site | Google Scholar
N. O. Akinci, “An investigation on seismic resistance of reinforced concrete industrial chimneys,” in Proceedings of the 2009 Structures Congress, pp. 964–971, Austin, TX, USA, April 2009.
View at: Publisher Site | Google Scholar
W. Huang and P. L. Gould, “3-D pushover analysis of a collapsed reinforced concrete chimney,” Finite Elements in Analysis and Design, vol. 43, no. 11-12, pp. 879–887, 2007.
View at: Publisher Site | Google Scholar
V. Lopes, J. M. Guedes, E. Paupério, A. Arêde, and A. Costa, “Ambient vibration testing and seismic analysis of a masonry chimney,” Journal of Building Appraisal, vol. 5, no. 2, pp. 101–121, 2009.
View at: Publisher Site | Google Scholar
F. Minghini, G. Milani, and A. Tralli, “Seismic risk assessment of a 50 m high masonry chimney using advanced analysis techniques,” Engineering Structures, vol. 69, pp. 255–270, 2014.
View at: Publisher Site | Google Scholar
C. Zhou, X. Zeng, Q. Pan, and B. Liu, “Seismic fragility assessment of a tall reinforced concrete chimney,” The Structural Design of Tall and Special Buildings, vol. 24, no. 6, pp. 440–460, 2015.
View at: Publisher Site | Google Scholar
C. D. Zhou, X. L. Zeng, F. Zhao, and H. Zhou, “Seismic fragility analysis for high-rise reinforced concrete chimney,” Earthquake Engineering and Engineering Dynamics, vol. 36, no. 2, pp. 173–181, 2016, in Chinese.
View at: Google Scholar
C. D. Zhou, M. W. Tian, X. Zhang, P. G. Wang, and andX. Ma, “Seismic fragility analysis for high-rise RC chimney considering multi-dimensional seismic actions,” China Civil Engineering Journal, vol. 50, no. 3, pp. 54–61, 2017, in Chinese.
View at: Google Scholar
X. Chen, A. Q. Li, W. R. Cheng, and Y. Wang, “The Adomian decomposition method for dynamic characteristics of high-rise buildings,” Journal of Vibration Engineering, vol. 26, no. 4, pp. 493–499, 2013, in Chinese.
View at: Google Scholar
R. Sancibrian, I. Lombillo, E. G. Sarabia, Y. Boffill, H. Wong, and L. Villegas, “Dynamic identification and condition assessment of an old masonry chimney by using modal testing,” Procedia Engineering, vol. 199, pp. 3410–3415, 2017.
View at: Publisher Site | Google Scholar
T. Tatara and B. Ratajewicz, “The selection of a dynamic model of a RC chimney based on in Situ research,” in Proceedings of the International Conference on Experimental Vibration Analysis for Civil Engineering Structures, pp. 619–630, San Diego, CA, USA, July 2017.
View at: Google Scholar
D. Mehta and N. J. Gandhi, “Time response study of tall chimneys under the effect of soil structure interaction and long period earthquake impulse,” in Proceedings of the 14th World Conference on Earthquake Engineering, Beijing, China, October 2008.
View at: Google Scholar
R. Livaoglu, “The numerical and empirical evaluation of chimneys considering soil structure interaction and high-temperature effects,” Soil Dynamics and Earthquake Engineering, vol. 66, pp. 178–190, 2014.
View at: Publisher Site | Google Scholar
A. M. Halabian and S. Kabiri, “Effect of foundation flexibility on ductility reduction factors for R/C stack-like structures,” Earthquake Engineering and Engineering Vibration, vol. 10, no. 2, pp. 277–290, 2011.
View at: Publisher Site | Google Scholar
S. V. Jisha, B. R. Jayalekshmi, and R. Shivashankar, “3D soil-structure interaction analyses of annular raft foundation of tall RC chimneys under wind load,” Indian Geotechnical Journal, vol. 44, no. 4, pp. 409–426, 2014.
View at: Publisher Site | Google Scholar
B. R. Jayalekshmi, S. V. Jisha, and R. Shivashankar, “Analysis of foundation of tall R/C chimney incorporating flexibility of soil,” Journal of the Institution of Engineers (India): Series A, vol. 98, no. 3, pp. 211–217, 2017.
View at: Publisher Site | Google Scholar
C. D. Zhou, X. L. Zeng, J. Chen, and B. Liu, “Seismic collapse resistance analysis of tall reinforced concrete chimney,” Engineering Mechanics, vol. 33, no. 5, pp. 57–65, 2016, in Chinese.
View at: Google Scholar
C. D. Zhou, X. L. Zeng, and Q. L. Pan, “Seismic response analysis of high-rise RC chimney including soil-structure interaction,” Earthquake Engineering and Engineering Dynamics, vol. 36, no. 3, pp. 140–150, 2016, in Chinese.
View at: Google Scholar
R. Sarkar, D. Shrimal, and S. Goyal, “Seismic analysis of a 275 m tall RCC multi-flue chimney: a comparison of IS code provisions and numerical approaches,” Advances in Structural Engineering, vol. 2, pp. 1015–1025, 2015.
View at: Publisher Site | Google Scholar
J. L. Wilson, “Earthquake response of tall reinforced concrete chimneys,” Engineering Structures, vol. 25, no. 1, pp. 11–24, 2003.
View at: Publisher Site | Google Scholar
H. Q. Liu, W. B. Li, R. Y. Wu, and L. F. Yan, “Elastic-plastic seismic response analysis on the special-shaped chimney,” Applied Mechanics and Materials, vol. 711, pp. 520–524, 2014.
View at: Publisher Site | Google Scholar
N. S. Abhyankar, “Flexural waves in beams with geometrical and material longitudinal discontinuities,” International Journal of Impact Engineering, vol. 9, no. 2, pp. 205–222, 1990.
View at: Publisher Site | Google Scholar
Information on https://ngawest2.berkeley.edu/, Peer Ground Motion Database, Pacific Earthquake Engineering Research Center.
GB50010-2010, Code for Seismic Design of Buildings, China Architecture & Building Press, Beijing, China, 2010.
China Architecture & Building Press, Design Manual of Chimneys, China Architecture & Building Press, Beijing, China, 1989.
GB50010-2010, Code for Design of Concrete Structures, China Architecture & Building Press, Beijing, China, 2015.

Copyright

Copyright © 2020 Yingchun Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1331

Downloads

854

Citations