Steady Wind Load on External Surface and Its Effect on Wind-Induced Response for a 200m High Natural-Draught Cooling Tower

Wind tunnel tests were carried out to measure the wind pressure of a 200m high natural-draught cooling tower. An analysis of the distribution characteristics of external pressure was then conducted to determine the pressure coefficients Cp(θ, z) in a given wind profile. Finally, the effect on the response of the shell and the buckling safety of the shell, applying the simplified height-constant pressure coefficient Cp(θ) and the realistic pressure Cp(θ, z), was determined. Taking the wind load specified in the code as an example, the influence of the distribution of external pressure on the wind-induced response was further analyzed. (e results indicate that the pressure distribution varies with not only the height z but also the circumferential angle θ, and the wind load of both ends of the tower is significantly greater than that of its middle. Moreover, the wind-induced static responses of the tower under the action of the realistic pressure distribution Cp(θ, z) and the simplified approach Cp(θ) are basically consistent, because the wind load distribution is more important than its magnitude for the wind-induced response of cooling tower, and the windinduced response of the cooling tower is dominated by the local shell deformation.


Introduction
Natural-draught cooling towers are high-rise and thinwalled flexible shell structures commonly used in thermal and nuclear power plants as cooling devices. e shell of a cooling tower is rather thin, and the ratio of the minimum shell thickness to the throat diameter is of the order of 1/400. In areas of negligible seismic activity, the dominating load of natural-draught cooling towers is induced by the action of turbulent wind due to the exposure of a large area to the wind (Bamu and Zingoni [1]; Nimeann and Zerna [2]; Hashish and Abu-Sitta [3]; Armitt [4]). In November 1965, three out of a group of eight reinforced-concrete cooling towers were blown down by strong wind at Ferrybridge Power Station in Yorkshire County in the UK (Central Electricity Generating Board, 1965) [5]. Since then, much attention has been paid to the wind load and wind-induced response of cooling towers (Zou et al. [6]; Ke and Ge [7]; Ke et al. [8]; Zhao et al. [9]). Following Davenport's concept, the wind forces acting on the external surface of a cooling tower may be subdivided into mean forces (averaged over a period of some 10 minutes) and turbulent wind forces. Usually, the resonant response induced by turbulent wind forces is small and an appropriate gust response factor taking into account the effects of turbulence is used to determine a static design load (Zhang et al. [10]; Zou et al. [11]). An outline of the quasistatic wind action in cooling tower shells used in the current design practice was presented by Abu-Sitta and Hashish [12]. e effects of the wind load are usually determined by applying equivalent static loads in terms of surface pressures (Kasperski and Niemann [13]; Ke et al. [14]). Once quasistatic pressure distribution is defined, the response to turbulent wind is introduced by the gust factors. Mean wind pressures are usually characterized by the pressure coefficients Cp.
where p(θ, z) is the mean wind pressure acting on the external surface of the shell depending on the circumferential angle θ and the height z above the ground level and q(z) is the reference velocity pressure at height z. As seen in equation (1), the pressure coefficients Cp depend on both the height z and the circumferential angle θ, Cp � Cp(θ, z), which is also confirmed by wind tunnel tests in shear flow. However, since the wind pressure distributions in various loading codes are based on data from on-the-spot measurements, in which sensors were usually installed at the area near the throat, and wind tunnel tests in homogeneous flow for cooling towers with heights of about 90 m∼120 m (Sun and Zhou [15]; Sun et al. [16]), it has become common practice to use a simplified height-constant pressure distribution Cp(θ) in the current design, disregarding the fact that Cp values are height-dependent, and to attribute the influence of heights on the pressures to the wind profile alone. e wind pressure p(θ, z) at any point on the outer surface is then simplified as follows: However, with the rapid development of nuclear power plants on the mainland of China, a 200 m high naturaldraught cooling tower has been proposed whose height would far exceed the 165 m height limit specified in the existing design codes for cooling towers in China (GB/T 50102 [17]; NDGJ5-88 [18]). Accordingly, the design value of the wind loading on the external surfaces of a cooling tower with a height of 200 m has attracted the extensive attention of Chinese engineers. e aim of the present study is to elaborate the differences between the simplified heightconstant pressure distribution Cp(θ) and the realistic pressures Cp(θ, z) and to find out whether the design procedure is realistic or, at least, conservative. Firstly, the wind pressure coefficients on the external surface of a 200 m high natural-draught cooling tower are obtained by synchronous pressure measurements in wind tunnel tests in shear flow, which is a closer approximation of reality, to determine the pressure coefficients Cp(θ, z) in a given wind profile. en, the resulting effect on the response of the shell and the buckling safety of the shell applying the simplified height-constant pressure coefficient Cp(θ) or the realistic pressure Cp(θ, z) was established. It should be noted that wind-induced interference effects in tower groups or from other adjacent buildings are not considered in the present study, so the results are restricted to the case of a single tower.

Details of Wind Tunnel Tests
e wind tunnel test is conducted at the HD-2 BLWT at Hunan University, whose laboratory is a closed-circuit atmospheric boundary layer.
e tests are carried out in turbulent shear flow in the high-speed test section (3.0 m in width and 2.5 m in height), and an artificial thickening of the boundary layer was achieved with the help of spires at the entrance and irregularities on the floor of the wind tunnel. e flow velocity was measured with a Cobra Probe, and the sampling time and sampling frequency were, respectively, set as 30 s and 2000 Hz. e duration of each sample was chosen to obtain an error of less than 0.5% on the mean value. e wind characteristics measured at the center of the rotary plate are shown in Figure 1. As shown in Figure 1(a), the mean wind speed profile is represented by the power law, and the profile exponent α � 0.16, which corresponds to moderately rough terrain as specified in the Chinese loading code. e turbulence intensities at z � 50 cm and z � 40 cm above the floor are 15.6% and 16.6%, corresponding to the top and throat levels of the model tower, respectively. Figure 1(b) shows that the simulated longitudinal turbulence power spectrum at 50 cm above the floor is similar to common theoretical spectra, such as those generated by Von Kármán, Kaimal, and Davenport, which illustrates that the wind field simulation is reliable. In addition, the integral length scale is about 50 cm, corresponding to a full-scale value of 200 m.

Test Model.
e total height of the prototype cooling tower is 200.00 m, the throat section above ground level is 156.70 m, and the inlet opening above ground level is 12.59 m, with a cooling area of 18,000 m 2 . e diameters of the top, throat section, and bottom of the shell are 96.60 m, 94.60 m, and 153.00 m, respectively. e shell thickness varies with the height z. e minimum shell thickness (0.25 m) is at the throat section, whereas the maximum shell thickness (1.40 m) is at the lower stiffening ring, which is supported by 52 uniformly distributed pairs of Φ1.4 m Vshaped circular columns. e testing model was made from high-quality aluminum alloy by a digitally controlled lathe to ensure that the strength and stiffness of the model are sufficient and no deformation or obvious vibration occurs under the wind speed during the pressure measuring tests to guarantee the precision of the pressure measurements. e experimental model is geometrically similar to the prototype object, and the geometric scale ratio was chosen to be 1/400 (total height of the model H � 0.5 m). e maximum blockage ratio caused by the installation of the model is 4.6%, and thus its effects are neglected in the following discussions. e shell model was supported by 52 pairs of V-shaped circular column models, permitting the reproduction of the real draught. e photo of the processed model is given in Figure 2. A total of 20 levels of measurement points were arranged on the internal surface of the model, and 36 measurement points were placed at each level along the circumferential direction (i.e., an angular distance of 10°), resulting in 720 total pressure taps. Figure 3 presents the layout of the measurement points and the definition of the circumferential angle θ.

Test Conditions.
A DTC net electronic type pressure scanning system (Pressure Systems, Inc., USA) was employed to measure the wind pressure. In this study, sixteen modules were used and a total of 1024 pressure measuring points could be monitored simultaneously. e sampling duration of each measurement was 30 s and the sampling frequency was 330 Hz. e pressure taps were linked to the transducers with 500 mm silicon tubes with an inner diameter of 1.0 mm and an outer diameter of 2.0 mm. e system obtained had a flat amplitude and a linear variation of the phase up to the sampling frequency, which could guarantee the transmission of the fluctuations without distortions.
It is important to note that the flow conditions around rounded shapes are sensitive to the Reynolds number. However, the Reynolds number in a wind tunnel is usually 2∼3 magnitudes less than that in full scale. Fortunately, the transcritical flow regime may be obtained by adding sufficient roughness to the model surface in the wind tunnel tests. It should be noted that the design curve is simulated by pasting a rough paper tape on the surface of the model to increase the surface roughness in the wind tunnel test, as shown in Figure 4. By making use of this beneficial effect, thirty-six equally spaced ribs 0.3 mm in height were applied in the present study to provide an average roughness coefficient rib height to a distance of 0.015. In a preliminary test, it was found that the pressure distribution was unaffected for wind speeds V H from 12 to 20 m/s. is was monitored using a Pitot tube at the centerline of the test section, about 3 m upstream of the tested model. e test results obtained at the free stream oncoming velocity V H � 18 m/s were taken as representative, corresponding to a Reynolds number of 6.03 × 10 5 based on V H and the mean diameter of the tower model.

Test Results of Mean Wind Load.
To enable the data obtained from the wind tunnel tests to be applied to the prototypes, the pressure coefficients presented in this study are nondimensionalized, referring the pressures at a given height z to the velocity pressure q(z) at the same height, as shown in equation (1). A positive pressure indicates the direction of the pressure force points into the surface of the structure, while a negative one indicates the direction points outwards from the surface. Since the cooling tower is axially symmetrical and the theoretical mean wind pressure distribution on the external surface should be symmetrical compared with the incoming wind direction, only half of the wind pressures (0°to 180°) are presented in this paper.  e mean wind pressure coefficients at typical levels z/H measured at the free stream oncoming velocity V H � 18 m/s are shown in Figure 5. It should be noted that the "typical curve" is the average of the pressure coefficients over the region z/H � 0.30 to 0.80, since the pressure distribution at this region is more or less independent of the height z. is is in accordance with the data reported by Sun and Zhou [15] and Harnach and Niemann [20]. It is observed that the curves of the pressure pattern can be broadly divided into three regions: area I, area II, and area III. In area I (0°≤ θ < 40°), the pressure is positive and more or less independent of the height z; in particular, at stagnation

Shock and Vibration
(θ � 0°), the coefficient is still approximately unity. In area II (40°< θ < 110°), the variation of the pressure coefficient with height is obvious. In particular, the typical and sensitive measure of the pressure distribution suction peak min Cp varies from −1.84 at z/H � 0.13 to −1.08 at z/H � 0.95. In area III (110°< θ < 180°), since the flow separates from the sides of the body (θ �110°), the base pressure coefficient Cp,b is negative and independent of the circumferential angle θ. In addition, the variation of the base pressure coefficient Cp,b with height is small, apart from an area near the tower bottom, due to the flow separating from the lower edge of the shell. Figure 6 shows typical parameters of the pressure distribution varying with height, where the pressure rise coefficient ΔCp representing the pressure rises downstream of the suction peak, min Cp, to the base pressure coefficient Cp,b at a horizontal level, is given by where Pb(z) and min P(z) are the base pressure and minimum pressure at height z, respectively. From Figure 5, it is clear that the typical parameters of the pressure distribution in the lower third of the shell are  Rough tapes on the model Pressure represents the local load acting on the surface of the shell; however, global forces such as resistance should also be considered in the design. For structures with circular cross sections, such as a cooling tower, the mean drag coefficient C D in the downwind direction and lift coefficient C L of a section obtained by the area-weighted integral of the mean pressure distributions are given by where N is the number of pressure taps arranged uniformly in the circumferential direction and Cp i and θ i represent the pressure coefficient and circumferential angle at point i, respectively. Figure 7 shows the result of the aerodynamic coefficients at each section. It can be seen that the average value of the lift coefficient C L is basically constant along the height and is mostly close to zero. In fact, for the circular cross-section structure of the cooling tower, due to the symmetrical distribution of the wind pressure on the surface, the theoretical average value of the lift coefficient is zero. e consistency of the theoretical values verifies the correctness of the test results in the present study. e three-dimensional effect of the drag coefficient C D is significant, showing the characteristics of a large end and a small middle. e drag coefficient at the top of the tower (C D � 0.72) is 41.2% larger than the middle section, and the drag coefficient at the bottom of the tower (C D � 0.90) is 76.5% larger than the middle section. It can be seen that if the design load is taken according to the "representative curve," the design load at both ends of the cooling tower will be significantly smaller than the actual load at the ends.

Details of the Numerical Calculation.
e numerical calculations of the wind-induced static response of the prototype cooling tower were carried out with the help of the finite element software ANSYS. e internal forces were calculated using the membrane theory of shells. e body of the tower was simulated using a SHELL63 element with four nodes and both membrane and flexural stiffnesses, and the columns were idealized as three-dimensional Timoshenko beams with a BEAM188 element. e finite element model and analysis result of the first-order modal are shown in Figure 8. e fundamental frequency f is 0.738 Hz, and there are 4 circular waves and 2 vertical waves in the mode. e combined effects of the dead load g, the wind load acting on the external surface of the shell Pe, and the internal suction Pi were considered. ey usually produce extreme stresses in the tower. For the wind load, a velocity pressure of q(z) � q (10) * (z/10) 2α � 0.45(z/10) 0.32 is applied, where q kN/    Shock and Vibration m 2 and z m, α � 0.16 is the dimensionless power exponent, and q (10) � 0.45 kPa. As mentioned before, two different pressure coefficients were investigated, that is, the simplified height-constant pressure coefficient Cp(θ) and the realistic pressure Cp(θ, z). Hence, two different wind loads acting on the external surface of the shell are given by In addition to the external surface, the internal surface of the cooling tower also suffers wind load due to allowing wind to pass through to the circulating cooling water. Since no regulations concerning the internal suction of the cooling towers have been included in the existing Chinese code, the value adopted in the present study is incorporated in the German code (VGB [21]). e internal suction Pi is constant with respect to the height z and the circumferential angle θ and related to the velocity pressure q(H) at the top of the tower by a pressure coefficient of −0.50. e buckling safety of the shell is one of the essential factors considered in the design of cooling towers. e Chinese codes have provided stability checking formulas for the entire cooling tower and the tower body, respectively. ey are where q cr is the critical wind pressure; E is the elasticity modulus of concrete; h is the wall thickness at the throat of the cooling tower; r 0 is the diameter at the throat; K B is the safety factor for elastic buckling; σ 1 and σ 2 refer to the circumferential and meridian pressures, respectively; c is the Poisson ratio of concrete; K 1 and K 2 are the parameters determined from the geometric parameters of the tower body; σ cr1 and σ cr2 are the critical circumferential and meridian pressures, respectively, which are given by V cr � 40 �� � q cr √ .

Shock and Vibration
According to equation (5), the unstable critical wind speed V cr is calculated by e critical circumferential and meridian pressures, σ cr1 and σ cr2 , can be calculated by equations (7) and (8), respectively. Substituting the hoop and meridional stresses calculated from ANSYS into equation (6), the safety coefficient of the local stability K B can be obtained.

Results and Analysis of Static Responses
Since the maximum positive value (tension) of the membrane forces caused by wind load and dead load is usually at the windward meridian (0°meridian) and their maximum negative value (compression) is usually at about the 70°m eridian, the membrane forces at the 0°and 70°meridian were taken as examples to analyze the influence of the pressure distribution on the static response of the cooling tower in this study. e membrane forces induced by the combined effects mentioned before are plotted in Figure 9. From Figure 9(a), it is observed that the simplified heightconstant pressure distribution Cp(θ) results in considerably higher meridional axial forces compared to the realistic pressure distribution Cp(θ, z). e increase averaged over the height of the shell is about 20%. However, as seen in Figure 7, the overall aerodynamic drag of the realistic pressure distribution (the average value over the height is 0.61) is higher than the simplified (0.51). is is because cooling towers, like spatial structures in general, are sensitive not only to load intensity but just as much to the load distributions, which will be further explained in another section below. However, the pressure distributions exert little influence on the maximum tensile and compressive hoop axial force, as can be seen in Figure 9(b). e realistic pressure distribution produces slightly higher hoop axial forces in the lower third of the shell but smaller ones in the upper fifth of the shell. is is because the hoop forces are created mainly by the local effect of wind pressure at a given level, and the peak suctions of the realistic pressure distribution at the lower third and upper fifth of the shell are larger than those in the simplified pressure distribution. is is in accordance with the results reported by Harnach and Niemann [20]. Figures 9(c) and 9(d) show the comparison of the calculation results of the bending moments. It can be seen that the meridian and hoop bending moments under the two loads are basically the same along the height change trend, and the magnitudes are similar. Figure 10 shows the comparison of the absolute value of the displacement under two wind pressure distribution curves. It can be seen that the trend of the deformation of the cooling tower is basically the same. e maximum deformation occurs near the throat, and the displacement under the action of the "sexual curve Cp(θ)" is significantly larger (up to 15%) than that of the three-dimensional curve Cp(θ, z) at the height above 1/3 H. However, at the height below 1/ 3 H, the displacement under Cp(θ) is slightly smaller than that of Cp(θ, z). In order to observe the difference between the internal forces under the two loads more clearly, the meridian stress along the hoop direction at the typical upper, middle, and lower height sections of the cooling tower is displayed in Figure 11. It can be seen that the meridian stress under different loads has the same distribution trend along the hoop direction, and the stress peak value of the "representative curve Cp(θ)" is slightly larger than that of the three-dimensional curve Cp(θ, z). Table 1 shows the comparison of the maximum response and the occurrence location under the two wind loads. It can be seen from the table that the maximum response location under the two loads is basically the same, and the maximum displacement of the "representative curve Cp(θ)" is slightly larger than that of the three-dimensional curve Cp(θ, z). e maximum principal tensile stress is quite different from the maximum principal compressive stress. e maximum principal tensile stress of the "representative curve" is 52% larger than that of the three-dimensional curve, and the absolute value of the maximum principal compressive stress is 52.7% smaller. e buckling safety is tabulated in Table 2. Comparing the simplified height-constant pressure distribution Cp(θ) and the realistic pressure distribution Cp(θ, z), it can be seen that the buckling safety is generally smaller in the first case (the difference is 6.1%). In other words, with relation to buckling safety, the simplified approach Cp(θ) is always conservative. is also means that when determining the shell thickness for a given minimum buckling safety, the required wall thickness will be somewhat reduced when the realistic pressure distribution Cp(θ, z) is taken into account. In general, the effect of the pressure distribution on the buckling mode is slight.

Wind Pressure Curves Detailed by the Specifications.
e current relevant design codes for cooling tower design in China respectively give two average wind pressure coefficient distribution curves for smooth and ribbed hyperbolic cooling towers, and the expression adopts the Fourier series as listed in equation (11) (GB/T 50102, 2003; NDGJ5-88, 2006). e German standard gives 6 curves according to the surface roughness and uses the piecewise function expression (VGB, 2010), but Gould and Kratzig [22] pointed out that the German curve can also be expressed by the Fourier series polynomials for design input. erefore, the least squares method is adopted to fit the German code curves into an eight-term Fourier series in the present study, which unifies the expression form of the Chinese and German gauge curves. e comparison between the Chinese and German standard curves and the harmonic coefficients of each order is shown in Figure 12 and Table 3.
where a i is the i th harmonic coefficient; θ is the angle with the incoming flow direction. Taking VGB-K1.0 as an example, the comparison of the wind pressure coefficients of the various harmonic components is depicted in Figure 13. It can be seen that a 0 is the 8 Shock and Vibration     direct current component, which is equivalent to the uniform pressure acting in the radial direction of the cooling tower. Since a 0 is a negative value, its effect is expressed as suction. a i is the wind pressure coefficient of the i-th harmonic. e wind pressure coefficient of the 7th harmonic fluctuates near zero, and its shape is approximately a straight line because the coefficient a 7 is small. e sum of those terms is the VGB-K1.0 wind pressure curve. Table 4 shows the comparison between the downwind drag coefficient and the total drag coefficient obtained by integrating the corresponding wind pressure coefficients of each harmonic component of the specification curve along the circumference. It can be seen from the table that the downwind drag coefficient is completely contributed by the wind pressure coefficient of the first harmonic, and the wind pressure coefficients of other harmonics contribute zero to the total drag coefficient. In fact, from the calculation expression of the harmonic resistance coefficient of each order (equation (11)), only when i � 1, Cd i is not zero, so the total resistance C D � Cd i � Cd 1 � a 1 π/2.
From the analysis of the contribution of each order of harmonics to the drag coefficient, it can be seen that the DC component a 0 , which is probably equal around the surface, induces the uniform expansion of the circular section without rigid body displacement; that is, the drag coefficient

Shock and Vibration
of a 0 is zero. e first-order component a 1 , which may be a superposition of uniform force around the surface and a whole deflection force acting on the structure, causes rigid body displacement of the cooling tower without local deformation of the circular section; that is, the drag coefficient of a 1 is equal to the drag of the entire wind pressure curve zero. e higher-order component a 2 ∼a 7 , which is not a uniform force around the surface, induces localized deformation in the circular section, but no rigid body displacement; that is, the drag coefficient of a 2 ∼a 7 is zero. A schematic diagram of the radial displacement component generated by each order of harmonics is shown in Figure 14. It can be seen that the contribution of the order of harmonics is shown in Figure 14. e cross-sectional resistance coefficient of the cooling tower to the overall displacement of the structure is only the rigid body displacement.

Influence of Different Wind Pressure Curves on the Response of Cooling Tower.
Taking the finite element model in Figure 8 as an example, the wind pressure curves specified in Figure 12 are adopted to study the influence of different wind pressure curves on the response of the cooling tower. Because the meridian stress and other response results are basically consistent with the change law of the displacement response, the present study only analyzes the wind-induced displacement response. Figure 15 shows the calculation results of the displacement response at the throat of the harmonic components of each order of the Chinese ribbed curve. It can be seen that the distribution of the displacement components generated by the harmonics of each order is consistent with the analysis conclusion of Figure 14; that is, the DC component induces the uniform expansion of the circular section, the first-order component causes the rigid body displacement of the cooling tower, and the higher-order component causes the local deformation of the circular section. e inference of Figure 15 is further confirmed. e displacement of each harmonic component at the meridian of 0°(maximum positive pressure zone) and 70°( minimum negative pressure zone) is shown in Figure 16. It can be found from Figure 16(a) that the displacement of the first-order harmonic component of each standard curve increases along the height, and the displacement of the same height increases with the increase of the drag coefficient C D of the standard curve. is is because the effect of the firstorder harmonic is resistance. e cooling tower structure undergoes a deformation similar to that of a vertical cantilever beam under the action of this component. e deformation size increases with heights; the greater the resistance, the greater the deflection induced. It can be seen from Figure 16(b) that the displacement caused by the DC and high-order harmonic components shows a distribution trend that is large near the throat and small at both ends along the height. is phenomenon can be explained as the effect of DC and high-order harmonic components. e results show that the shell is locally deformed and the deformation is related to the local stiffness of the shell. e bottom is the area with the largest cooling tower wall thickness, so the stiffness of this area is the largest in the structure. e middle section is the thinnest part of the cooling tower wall, so its rigidity is the smallest. Although the top wall is thinner and the size is the same as the middle section, the rigid ring on the top greatly enhances the rigidity of the area. erefore, it is not difficult to understand that the  local deformation of the shell occurs near the throat. e displacement of the DC and higher-order harmonic components at the same height decreases with the increase of the resistance coefficient of the gauge curve. is is because the local deformation of the shell is more closely related to the local wind pressure. In fact, as the drag coefficient C D of the wind pressure curve increases, the minimum wind pressure coefficient amplitude decreases, and the maximum local load decreases, so the local deformation of the shell also decreases. e distribution law of the overall displacement is basically consistent with the displacement of the DC and high-order harmonic components (Figure 16 that the wind-induced deformation of the cooling tower is dominated by the local deformation of the shell. Figure 17 shows the 0°and 70°meridian displacement of the first-order harmonic component as a percentage of the total displacement. It can be seen that the displacement ratio of the first-order harmonic component is mostly below 10%, which increases slightly as the resistance increases. Contrary to the distribution trend of the total displacement, the proportion of the first-order harmonic component changes along the height following a trend wherein it is small in the middle and large at both ends. Figure 18 shows the comparison of the displacement percentage of the harmonic components of each order of the standard curves. It can be seen that the displacement percentages of the harmonic components of each order of the different standard curves are basically the same. e third-order harmonic has the largest proportion, which is about 55%, followed by the second-order harmonic (about 35%), while the first-order harmonic is only about 5%.
It can be clearly seen from Figure 19 that the displacement of the first-order harmonic component increases with the increase of the resistance, but its largest proportion of the total displacement is only 7%. In spite of the fact that the displacement of the DC and higher-order harmonic components decreases with the increase of the drag coefficient C D , its minimum proportion is still above 93%. e "increase" is much smaller than the "decrease," so the total displacement decreases as the drag coefficient increases.

Conclusions
Based on rigid model pressure measurement wind tunnel tests, the three-dimensional wind pressure coefficient distribution curve C p (z, θ) of a 200 m high super large cooling tower in a nuclear power plant was obtained along the height and circumferential direction. e three-dimensional effect of the wind load on the outer surface of the cooling tower was analyzed. e static response and the buckling safety under the three-dimensional wind pressure coefficient curve C p (z, θ) and the "representative curve C p (θ)" were then calculated, and the influence mechanism of the external pressure on the wind-induced response was analyzed. e following main conclusions were drawn: (1) e wind pressure coefficient on the outer surface of the cooling tower is basically symmetrically distributed along the circumferential direction, and the average value of the lift coefficient is close to zero. e three-dimensional effect of the wind pressure distribution makes the three-dimensional effect of the drag coefficient significant. Besides, the drag coefficient is characterized by a large value at the end of the cooling tower and a small value at the middle of the cooling tower. e drag coefficients at both ends are greater than that of the middle section.
(2) e wind-induced response of the cooling tower is dominated by the local shell deformation. e rigid deformation caused by resistance does not exceed 10% of the total response, while the local shell deformation caused by other harmonics accounts for more than 90% of the total response. e overall response of the cooling tower has no absolute relationship with the drag coefficient but is closely related to the characteristics of the wind pressure distribution.
(3) e coefficient distribution curve of the three-dimensional average wind pressure C p (z, θ) along the height and the circumferential direction is simplified to the wind pressure distribution curve along the constant height C p (θ). Although the overall wind load of the structure is underestimated, the windinduced response of the cooling tower is more closely related to the distribution of the wind load. e response induced by the "representative curve C p (θ)" is basically the same as the three-dimensional curve C p (z, θ), and the meridian internal force is even larger than the three-dimensional curve. In general, replacing the three-dimensional curve with the average curve of each height in the middle section is an effective simplification method that can meet the design requirements of the project.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.  Figure 19: Percentage of first-order and high-order harmonic displacement (middle average). 16 Shock and Vibration