# Research on Wind Load Calculation Based on Identical Guarantee Rate Method

Authors: Tao Ye; Ledong Zhu; Zhongxu Tan; Lanlan Li
Journal: Shock and Vibration
 (2021)

---

## Abstract

Wind load on building surface is one of the main loads for structural design; scholars in this field have put forward some methods to calculate wind load, such as Simiu method and Kasperski method. Based on the basic theory of probability and the systematic analysis of the surrounding environment and turbulence, a random variable model for calculating wind load is established. According to the model, through the analysis of the relationship between guarantee rate and wind load, a numerical calculation method to calculate wind load is proposed based on extreme value analysis and polynomial fitting theory. To verify the performance of the algorithm, wind tunnel experiments were carried out to obtain a large number of first-hand measured data of high-rise building (Shanghai World Financial Center). Based on the measured data, the algorithm is simulated, and calculated results are analyzed, including wind pressure distribution on building and probability distribution of fluctuating wind pressure of some measuring points. The validity and accuracy of the proposed model and algorithm are verified by the comparative analysis and theoretical analysis of the calculation results.


---

## Body

## 1. Introduction

The wind load on the building surface is one of the main loads during the stage of the structure design. For high-rise buildings structure, the response caused by wind load accounts for a considerable proportion of the total load and even plays a decisive role. Accurate estimation of the wind load on the building surface envelope is critical to reducing the wind disaster loss [1, 2]. Reasonable wind-resistant design is of great significance to ensure the function of building structures. In the wind-resistant design of structures, to ensure the safety and reliability of structures, it is necessary to evaluate the extreme wind load reasonably. The extreme wind load is related to the wind characteristics and the aerodynamic characteristics of the structure. The former is in the form of wind speed, and the latter is in the form of the wind pressure coefficient. Targeting the wind-resistant design of high-rise buildings, in recent research studies, many wind load calculation and wind pressure extremum estimation methods are proposed. Zhang [3] provides a novel wind load calculation for the structural design concerning the directionality and uncertainty effect. Xue et al. [4] proposed a novel method of identifying structural state, parameters, and unknown wind load estimation from incomplete measurements. Chowdhury et al. [5] give a better understanding of wind load effects on roof-to-wall connections of a typical low-rise gable roof residential structure subjected to combined impacts of wind and a potential breach of the building envelope. Yao et al. [6] built a joint probability density function of along-wind direction and across-wind direction on wind load effect based on the high-frequency force balance test and Copula Frank function, and taking the rule of stress, a design basis solved the combination coefficient with a certain assurance of wind load component. Zhuang et al. [7] researched the non-Gaussian features of wind pressure fluctuations on a rectangular high-rise building when the wind was normal to the short and long edges through a wind tunnel test of the rigid model. Meng et al. [8] analyzed impacts of various parameters like turbulence model, approaching-flow speed, and grid type on wind pressure coefficients over CAARC buildings. Luo [9] regarded wind load as a function of average wind and load coefficient and analyzed the relationship among wind load, wind speed, and load coefficient. Bhattacharyya and Dalui [10] conducted a comprehensive study on the average pressure coefficient of “E” plan-shaped high-rise buildings using wind tunnel test and computational fluid dynamics (CFD) numerical simulation. Luo et al. [11] proposed a simplified formula for calculating wind pressure extremum based on Hermite polynomial and wind tunnel test data during the long sampling time intervals. Ma [12] proposed the wind pressure extremum algorithm combining the experimental data with numerical simulation. Malick et al. [13] put forward the wind pressure coefficient model of “C” plan-shaped tall building model and process of obtaining its solution. Quan et al. [14] proposed a wind pressure extreme value estimation method based on generalized extreme value theory. Quan et al. [15] proposed extreme value estimating method of non-Gaussian wind pressure.In terms of wind load calculation, this topic has been studied a long time ago. The widely used methods include the extreme value estimation method based on the assumption that the wind pressure or the extreme wind pressure on the building surface follows a certain probability distribution [16]. Based on the Davenport method for the Gaussian process, several scholars have proposed different improvements to the calculation model and determination method of the non-Gaussian peak factor of fluctuating wind pressure. Kareem and Zhao [17] expressed the non-Gaussian process as Hermite polynomials of the Gaussian process considering higher-order statistics, thus extending the peak factor method to the non-Gaussian process. The process conversion method of Simiu and Filliben [18] uses the principle of equivalent probability to transform non-Gaussian process into Gaussian process and proposes a method for calculating the extreme value of non-Gaussian process. Simiu method [18, 19] transforms the most unfavorable annual wind load sequences into the annual equivalent wind speed sequences from the point of view of the wind pressure, where the equivalent wind speed reflects implicitly the influence of wind direction. By assuming that the annual equivalent wind speed series obeys a certain probability distribution, the wind load value with a certain return period is obtained. Kasperski and Hoxey [20–22] improved the abovementioned method and proposed a more accurate extreme value calculation method (Kasperski method). Diniz et al. method [23] is based on Monte Carlo simulation and introduces a series of factors of cognitive uncertainty to study the wind load estimation method considering wind direction more comprehensively and accurately. In their methods, the peak factor method is used to obtain the most unfavorable wind pressure coefficient extreme value under the whole wind direction, and the wind pressure coefficient time history of the actual measuring point is obtained through wind tunnel test, and the wind load estimation value is obtained according to Bernoulli expression. The existing calculation method requires a large amount of data and calculation, a complex calculation process, but inaccurate results are obtained. Therefore, it is not suitable for engineering applications. It is valuable to design a method with “simplifying calculation and ensuring accuracy” to estimate wind pressure.Based on the basic probability theory and combined with the systematic analysis of the surrounding environment and turbulence of buildings, a random variable model is established to calculate the wind load considering wind direction. Through the analysis of the relationship between the guarantee rate and the wind load, the “asymmetric quantile truncation method” is proposed to determine the quantile of the wind pressure extreme value. On this basis, a numerical solution for the design wind load is proposed. From the point of view of structural wind resistance design, the numerical model of wind load calculation and the design structure wind load are obtained. Taking the wind pressures on claddings of the Shanghai World Financial Center as an example, the accuracy and effectiveness of the proposed algorithm are systematically verified.


## 2. Wind Load Calculation Model

At present, to estimate wind pressure, researchers first use the peak factor method to calculate the most unfavorable extreme value of the wind pressure coefficient under the whole wind direction and obtain the wind speed extreme value of the annual wind speed extreme series of the whole wind direction from the meteorological data and then take the calculated quantile value corresponding to a certain guarantee rate (return period) as the design wind speed extreme value. The influence of wind direction is not considered in the above calculation of the two parameters. The calculated wind pressure coefficient and wind speed extremum are substituted into Bernoulli formula to obtain the wind pressure design value (hereinafter referred to as the traditional method). The traditional calculation method of design wind load is based on model:(1)Wdes=12ρV^des2C^des,where Wdes is the design wind load, V^des2 is the extreme design wind speed, and C^des is the extreme design wind pressure coefficient. The design extreme wind pressure coefficient can be calculated by peak factor method.In model (1), the influence of wind direction is not fully considered, so the calculated wind load based on this model lacks the results caused by wind direction effect. In the process of calculating wind load, the probability distribution that annual wind speed discrete time history under all wind direction follows is obtained by assuming based on the original wind speed data, and then the quantile value corresponding to the specified guarantee rate is obtained, which is the extreme wind speed. For the wind pressure coefficient, their most unfavorable value under the whole wind direction is regarded as the extreme wind pressure coefficient, and the synchronous influence caused by wind speed and wind pressure coefficient on wind load is ignored. In this way, the actual guarantee rate corresponding to the design wind load obtained is quite different from the guarantee rate that the engineering application needs to meet.For a certain measuring point on the building surface, there are great differences in the extreme values of wind pressure under different wind directions. There are two reasons for the difference: one is that the extreme wind speeds at different wind directions are different, which is mainly related to the macro meteorology in the region, such as climate or cyclone type; the other is that the extreme wind pressure coefficients at different wind directions are different, which is mainly related to the environment around buildings and wind flow [24, 25]. The impact of wind on buildings and the surrounding environment is as follows: (1) in the area with dense high-rise buildings, the original wind field is changed by the buildings; under the same conditions, the local wind speed around the buildings increases; (2) wind load is a kind of random load, which is greatly affected by building height, wind direction, wind intensity, and duration; (3) the outer contour of a building is generally nonmanifold, so the flow field is inevitably accompanied by separated flow, vortex shedding and oscillation, and even more seriously, the coupling oscillation of structure and fluid; (4) the roughness of the building surface also affects the magnitude of the wind force, and the roughness of the building surface will also increase the effect of the wind force. The interaction between buildings and surrounding environment will be caused by wind. The shape, size, and number of buildings, the relative position of buildings, wind direction angle, wind field, and wind speed are all related to the mutual interference. The model described by formula (1) and the existing wind load calculation methods [18, 22, 23] do not fully involve these factors. In this paper, a series of factors of cognitive uncertainty are considered in the design of the algorithm model, so as to more comprehensively and accurately study the wind load estimation method considering wind direction, in which a series of factors of cognitive uncertainty include the wind-induced effect in short duration τ (bending moment, shear force, and displacement), the roughness degree z0 of the representative landform of the upwind location of the building, the wind speed conversion factor u, the wind speed v, the site type γ, shape coefficient Tx and wind pressure point correlation coefficient r, local shape coefficient of wind pressure μp, and variation of wind pressure coefficient. With the progress of technique for the measurement of unsteady aerodynamic force [26], wind pressure distributions and unsteady aerodynamic characteristics are investigated and analyzed [27], and the local aerodynamic force coefficient [28] plays an important role in studying the characteristics of unsteady aerodynamic forces acting on a structure. Based on the above analysis, the obtained wind load calculation model is more comprehensive than (1). The wind load model proposed in this paper is given by(2)Wdes=12ρV^des2v,θi,hC^desz0,θi,h,γ,r,uθi,h,C,μp/Cp,i=1,…,N,Nis number of wind angles,where Wdes is the design wind load, V^des2 is the extreme design wind speed related to wind angle θi, and C^des is the extreme design wind pressure coefficient related to wind angle θi. C is the aerodynamic force coefficient. μp/Cp is the local shape coefficient of wind pressure related to coherence of wind pressure.The algorithm proposed in this paper is based on the above model to calculate the wind load with the premise of uniform guarantee rate, so as to ensure a more safe and economical design.


## 3. Calculation of Wind Load Based on an Identical Guarantee Rate

### 3.1. Algorithm Overview

The random variables in the model (2) follow the nonlinear distribution, and the calculation of wind load based on this model involves the calculation of the square and product of random variables. To calculate the wind load according to model (2), it is necessary to solve the calculation problem of nonlinear random variables, which is difficult in probability theory. In addition, due to the complexity of wind field disturbance, it is difficult to find a certain probability distribution that wind pressure time history follows, which is not conducive to the application of practical engineering. Therefore, the numerical calculation method not only obtains the required wind load results, but also meets the accuracy requirements. This paper proposes a discrete numerical calculation method to get wind load. Based on the analysis of the relationship between guarantee rate and wind load, the “asymmetric quantile truncation method” is used to calculate the wind load related to a specific guarantee rate, and the maximum and minimum wind load sequences are obtained. Then, the required design wind load is calculated by using the algorithm proposed in this paper. Figure 1 is the schematic diagram of the algorithm model proposed in this paper.Figure 1
Modeling diagram for calculating wind load.


### 3.2. Determination of Asymmetric Quantile

In order to solve the aforementioned nonlinear probability model and calculate the wind load value in the model, a bilateral guarantee rate model is proposed in this paper. The wind pressure time history of high-rise building envelope surface is mostly non-Gaussian distribution, and there is a one-to-one correspondence between a specified guarantee rate and its corresponding peak factor. Therefore, the bilateral guarantee rate model can be determined according to the maximum sequenceC^P and minimum sequence C∨P of wind pressure time history. The guarantee rate corresponding to a quantile between the positive and negative extreme wind pressure coefficients symmetrical to the mean value is Pfix. It is defined as follows:The establishment process of the model can be seen in Figure2, and the model can be expressed by(3a)PC^P−PCP∨=Pfix,(3b)fk=PCPmean+k⋅CPrms−PCPmean−k⋅CPrms−Pfix=0,where k is the peak factor, while CPmean and CPrms are the mean pressure coefficient and root mean square of the fluctuating pressure coefficient, respectively. This algorithm is based on analyzing the relationship between the guarantee rate and wind load. In general, the corresponding guarantee rates for 50 years and 100 years are 98% and 99%, respectively P=1−1/T, and the corresponding guarantee rate for 1000 years is 99.9%. The principle of the algorithm is illustrated by taking the guarantee rate of 98%. When the guarantee rate is 98%, the right quantile of the extreme wind pressure probability density curve in Figure 3 is the wind pressure extreme value corresponding to probability 99%, and the left quantile is the extreme value of wind pressure corresponding to a probability of 1%.In the calculation process described below, the calculation symbols are defined as follows:W99%: extreme design wind load at quantile (value is 99%) on the right sideW1%: extreme design wind load at the quantile (value is 1%) on the left sideFigure 2
Bilateral guarantee rate model.Figure 3
Diagram of relationship between quantile and guarantee rate. (a) Illustration of design extreme wind pressure controlled byW99%. (b) Illustration of design extreme wind pressure controlled by W99%. (c) Illustration of design extreme wind pressure controlled by W1%. (d) Illustration of design extreme wind pressure controlled byW1%.
(a)(b)(c)(d)For a single random variable, there are four cases of the probability distribution and location quantile value in coordinate system, as shown in Table1. The four cases correspond to Figures 3(a)–3(d), respectively:Table 1
Quantile value in different cases.
Case numberCorresponding figure numberWind load is positive or negativeWind load controlling quantile valueQuantile 99%Quantile 1%1Figure3(a)PositivePositiveW99%2Figure3(b)PositiveNegativeW99%3Figure3(c)PositiveNegativeW1%4Figure3(d)NegativeNegativeW1%According to the above four extreme wind pressure probability density distributions, to meet the guarantee rate requirements, bothW99% and W1% need to be calculated, and then select the one with the larger absolute value. That is,(4)Wdesign=maxW0.99,design,W0.01,design.The meanings of the above symbols are as follows:Wdesign: design extreme wind load with at least 98% guarantee rateW0.99,design: design extreme wind load meeting 99% guarantee rateW0.01,design: design extreme wind load meeting 1% guarantee rateFor multiple extreme wind pressure probability distributions, at a certain measuring point underN wind angles, there are N extreme wind pressure coefficient time histories, from which N S-shaped probability distribution can be obtained. For these N wind pressure coefficient time histories, there is design wind load ≥load effect. The design wind load Wdesign can be obtained by several probability distribution curves. In the following figures, for simplicity, two S-shaped probability distribution curves are used to illustrate the idea of solving the guarantee rate (in fact, there are N S-shaped probability distribution curves; the solving method is the same as that for two S-shaped curves), as shown in Figures 4(a) and 4(b). In these two figures, W1 and W2 are the wind pressure extremum corresponding to a certain probability on the probability distribution curve of wind pressure extremum.Figure 4
Diagram of probability distribution curve corresponding to guarantee rate. (a) Probability distribution curve with 99% assurance rate. (b) Probability distribution curve with 1% assurance rate.
(a)(b)DefineW0.99,design as design wind load. Here, W0.99,design=maxW1,W2; it can be seen from Figure 4(a) that when guarantee rate P<0.99, wind load W (load effect) corresponding to red and blue curves is both less than or equal to W2; here,W2 is design wind load, and W2=W0.99,design=maxWi,i=1,2.DefineW0.01,design as design wind load. Here, W0.01,design=maxW1,W2, and it can be seen from Figure 4(b) that when guarantee rate P>0.01, wind load W (load effect) corresponding to red and blue curves is both less than or equal to W1; here, W1 is design wind load, and W1=W0.01,design=maxWi,i=1,2.The above analysis shows that if wind loadW≤Wdesign, thenwhenP=0.99, W0.99,design satisfies W0.99,design=maxWi,i=1,2,…N andwhenP=0.01, W0.01,design satisfies W0.01,design=maxWi,i=1,2,…N.DetermineWdesign and enable it satisfying Wdesign=maxW0.99,design,W0.01,design.The symbols in the above formulas are defined the same as those in (4). The word “at least” in the above definition of Wdesignmeans that the area of non-shaded area in Figures 5(a) and 5(b) is greater than 98%. The design extreme wind loads correspond to the quantiles W0.99,design and W0.01,design. Refer to Figures 5(a) and 5(b) below.The range of guarantee rate values in Figures5(a) and 5(b) is described in Tables 2 and 3:Figure 5
Diagram of the relationship between two design wind load:W0.99,design and W0.01,design. (a)
Gaussian distribution case; (b) non-Gaussian distribution case.
(a)(b)Table 2
Guarantee rate range (1).
Range of guarantee rateCondition>98%WhenW0.99,design and W0.01,design are not symmetric with respect to μW=98%WhenW0.99,design and W0.01,design are symmetric with respect to μWTable 3
Guarantee rate range (2).
Range of guarantee rateCondition>98%WhenW0.01,design+2μW−W0.01,design>W0.99,designWhenW0.01,design+2μW−W0.01,design<W0.99,design=98%WhenW0.01,design+2μW−W0.01,design=W0.99,designCombining with (4), it can be inferred that the guarantee rate corresponding to Wdesign is more than 98%. As shown in Figures 5(a) and 5(b), because the distance μW−W0.01,design between the left quantile and the mean μW is not necessarily equal to the distance W0.99,design−μW between the right quantile and the mean μW, this method is called “asymmetric quantile truncation method” with 98% guarantee rate in this case. It is similar to other guarantee rate values.


### 3.3. Numerical Algorithm

#### 3.3.1. Monte Carlo Sampling Simulation of Original Wind Speed Data

In the process of extreme value theory modeling and parameter fitting, problems such as missing wind speed sample data and data singularity are often encountered, which will decrease the accuracy of parameter fitting and the estimation of extreme wind speed, and Monte Carlo sampling of original data can solve these problems [29]. Based on the basic principle of Monte Carlo, the specific calculation steps of the Algorithm 1 are as follows:Algorithm 1: Process of original wind speed data.
(1)
The original wind speed sampling data is classified and counted according to the wind direction. There are 6000 wind speed samples at a certain wind direction in this numerical calculation case. The kernel density estimation method is used to fit these 6000 wind speed samples to obtain 6000 probability valuesPi=i=1,…,6000 corresponding to the probability distribution curve.(2)
The coordinates of 6000 points on the probability distribution curve are obtained：U1,P1,U2,P2,…,U6000,P6000.(3)
According to the fitting results calculated in Step 2, a mathematical model can be obtained:P=Fx. Then, U is substituted as an independent variable.(4)
ForP=fx in Step 3, find its inverse function: F−1x.(5)
SubstituteN random numbers P1,P2…,PN generated by uniform distribution U0,1 into F−1Pi=xi.(6)
In this way, the sample of Monte Carlo simulation with sample size ofN is obtained: x1,x2,…,xN. These are the new extreme wind speed samples with sample size of N: U1,U2,…,UN.


#### 3.3.2. Numerical Algorithm

Based on model (8) and the above theoretical analysis of the guarantee rate, this paper proposes an algorithm to calculate design wind load according to the wind speed data of wind tunnel test. The algorithm in this paper can be implemented by numerical calculation method, and the specific calculation process can be divided into the following steps:Step 1: Determine nonlinear wind load probability modelW^=1/2U^2C^p (or W^=1/2U^2Cˇp), C^p is maximum value, and Cˇp is minimum value.Step 2: Determine the time interval of two subsamples in the process of determining the probability model. This step is divided into the following two substeps:(1)
The time interval of extreme wind speed subsample is defined as:TU^′=TU¯/m. The length of duration of wind speed time history U¯ is TU¯. The time history of the extreme wind speed U^ (i.e., sub-sample) and its time interval TU^′ are determined.(2)
The time interval of extreme wind pressure coefficient is defined asTCp′=TCp/n. In this paper, the value can be calculated: 10 = 6000/600. The time history of extreme wind pressure coefficient C^pCˇp determined by the wind pressure coefficient time history Cp is 600.Step 3: SinceU^ and C^pCˇp are independent of each other, find the probability distribution that they obey.Step 4: For a certain measuring point, find probability distributionfC^pθi,h,i=1,2,…,N of C^pθi,h. Follow the same step for Cˇpθi,h.Step 5: For a certain measuring point, calculate the converted wind speedU^2 according to height from the original extreme wind speed U^1, and then do the Monte Carlo sampling simulation to find the simulated extreme wind speed U^ out of large sample, i.e., U^θi,h,i=1,2,…,N.Step 6: For a certain measuring point, calculate the probability distributionfU^θi,h,i=1,2,3…N of U^θi,h, and the probability distribution fU^2θi,h, i=1,2,3…N of U^2θi,h.Step 7: For a certain measuring point, find the probability distributionfZ1θi,h,i=1,2,3…N of Z1θi,h=U^2⋅C^p.Step 8: For a certain measuring point, find the probability distributionfZ1θi,h,i=1,2…N of Z1θi,h=1/2ρU^2⋅C^p.Step 9: For a certain measuring point, calculate the design wind loadWj,design, j=1,2,3,…,m with identical guarantee rate, and m is the total number of measuring points.For different measurement points, the same guarantee rate corresponds to different quantile values, so the algorithm is referred to as “based on identical guarantee rate.” The maximum and minimum wind load sequences can be calculated as follows:(1) Calculation of wind load maximum sequence: the wind load corresponding to 99% (or 99.9%) guarantee percentile of each distribution is maximum design wind loadW^0.99θI,h. The wind load value corresponding 1% (or 0.1%) guarantee rate quantile of each distribution is minimum design wind load W^0.01θI,h (negative maximum). If most of the maximum series samples are negative or all of them are negative, the design wind load can be calculated by W^0.01θI,h. For example, the wind load time history is shown in Figure 6.Figure 6
Most of time history are negative, or all of them are negative.(2) Calculation of wind load minimum sequence: the wind load corresponding to 99% (or 99.9%) guarantee percentile of each distribution is maximum design wind loadWˇ0.99θI,h. If most of the samples of the minimum series are positive, or all are positive, the design wind load can be calculated by Wˇ0.99θI,h. For example, the wind load time history is shown in Figure 7.Figure 7
Most of time history are positive, or all of them are positive.The minimum design wind loadW is taken as the value corresponding to the 1% (or 0.1%) guarantee rate quantile of each distribution Wˇ0.99θI,h (negative maximum). Then, the maximum or minimum wind load mentioned above has(5)W^0.99θI,hsatisfying:W^0.99θI,h=maxW^0.99θI,h,I=1,…N,I∈1,N,W^0.01θI,hsatisfying:W^0.01θI,h=maxW^0.01θi,h,i=1,…N,I∈1,N,Wˇ0.99θI,hsatisfying:Wˇ0.99θI,h=maxWˇ0.99θI,h,i=1,…N,I∈1,N,Wˇ0.01θI,hsatisfying:Wˇ0.01θI,h=maxWˇ0.01θI,h,i=1,…N,I∈1,N.According to the wind load maximum sequence and minimum sequence, design wind load of a measuring point is(6)Wj,design=maxW^0.99θI,h,W^0.01θI,h,Wˇ0.99θI,h,Wˇ0.01θI,h,where j=1,2…m,I∈1,N, m is the total number of measuring points, N is the number of wind directions, and it is a positive integer.For each measuring point, 4N wind load values were calculated by the above method based on U^ and C^p (Cˇp) at N wind directions. The maximum absolute value of these 4 N values is the design extreme wind load of this measuring point. The design extreme wind load of 360° in any direction can be obtained by harmonic function fitting. If the layers in height direction and N wind directions in horizontal direction are considered at the same time, this wind load calculation method can be applied to each small block (each measuring point).With the above method, the overall distribution of the maximum wind pressure of the equal probability design of the exterior envelope of the whole building can be obtained.


## 3.1. Algorithm Overview

The random variables in the model (2) follow the nonlinear distribution, and the calculation of wind load based on this model involves the calculation of the square and product of random variables. To calculate the wind load according to model (2), it is necessary to solve the calculation problem of nonlinear random variables, which is difficult in probability theory. In addition, due to the complexity of wind field disturbance, it is difficult to find a certain probability distribution that wind pressure time history follows, which is not conducive to the application of practical engineering. Therefore, the numerical calculation method not only obtains the required wind load results, but also meets the accuracy requirements. This paper proposes a discrete numerical calculation method to get wind load. Based on the analysis of the relationship between guarantee rate and wind load, the “asymmetric quantile truncation method” is used to calculate the wind load related to a specific guarantee rate, and the maximum and minimum wind load sequences are obtained. Then, the required design wind load is calculated by using the algorithm proposed in this paper. Figure 1 is the schematic diagram of the algorithm model proposed in this paper.Figure 1
Modeling diagram for calculating wind load.


## 3.2. Determination of Asymmetric Quantile

In order to solve the aforementioned nonlinear probability model and calculate the wind load value in the model, a bilateral guarantee rate model is proposed in this paper. The wind pressure time history of high-rise building envelope surface is mostly non-Gaussian distribution, and there is a one-to-one correspondence between a specified guarantee rate and its corresponding peak factor. Therefore, the bilateral guarantee rate model can be determined according to the maximum sequenceC^P and minimum sequence C∨P of wind pressure time history. The guarantee rate corresponding to a quantile between the positive and negative extreme wind pressure coefficients symmetrical to the mean value is Pfix. It is defined as follows:The establishment process of the model can be seen in Figure2, and the model can be expressed by(3a)PC^P−PCP∨=Pfix,(3b)fk=PCPmean+k⋅CPrms−PCPmean−k⋅CPrms−Pfix=0,where k is the peak factor, while CPmean and CPrms are the mean pressure coefficient and root mean square of the fluctuating pressure coefficient, respectively. This algorithm is based on analyzing the relationship between the guarantee rate and wind load. In general, the corresponding guarantee rates for 50 years and 100 years are 98% and 99%, respectively P=1−1/T, and the corresponding guarantee rate for 1000 years is 99.9%. The principle of the algorithm is illustrated by taking the guarantee rate of 98%. When the guarantee rate is 98%, the right quantile of the extreme wind pressure probability density curve in Figure 3 is the wind pressure extreme value corresponding to probability 99%, and the left quantile is the extreme value of wind pressure corresponding to a probability of 1%.In the calculation process described below, the calculation symbols are defined as follows:W99%: extreme design wind load at quantile (value is 99%) on the right sideW1%: extreme design wind load at the quantile (value is 1%) on the left sideFigure 2
Bilateral guarantee rate model.Figure 3
Diagram of relationship between quantile and guarantee rate. (a) Illustration of design extreme wind pressure controlled byW99%. (b) Illustration of design extreme wind pressure controlled by W99%. (c) Illustration of design extreme wind pressure controlled by W1%. (d) Illustration of design extreme wind pressure controlled byW1%.
(a)(b)(c)(d)For a single random variable, there are four cases of the probability distribution and location quantile value in coordinate system, as shown in Table1. The four cases correspond to Figures 3(a)–3(d), respectively:Table 1
Quantile value in different cases.
Case numberCorresponding figure numberWind load is positive or negativeWind load controlling quantile valueQuantile 99%Quantile 1%1Figure3(a)PositivePositiveW99%2Figure3(b)PositiveNegativeW99%3Figure3(c)PositiveNegativeW1%4Figure3(d)NegativeNegativeW1%According to the above four extreme wind pressure probability density distributions, to meet the guarantee rate requirements, bothW99% and W1% need to be calculated, and then select the one with the larger absolute value. That is,(4)Wdesign=maxW0.99,design,W0.01,design.The meanings of the above symbols are as follows:Wdesign: design extreme wind load with at least 98% guarantee rateW0.99,design: design extreme wind load meeting 99% guarantee rateW0.01,design: design extreme wind load meeting 1% guarantee rateFor multiple extreme wind pressure probability distributions, at a certain measuring point underN wind angles, there are N extreme wind pressure coefficient time histories, from which N S-shaped probability distribution can be obtained. For these N wind pressure coefficient time histories, there is design wind load ≥load effect. The design wind load Wdesign can be obtained by several probability distribution curves. In the following figures, for simplicity, two S-shaped probability distribution curves are used to illustrate the idea of solving the guarantee rate (in fact, there are N S-shaped probability distribution curves; the solving method is the same as that for two S-shaped curves), as shown in Figures 4(a) and 4(b). In these two figures, W1 and W2 are the wind pressure extremum corresponding to a certain probability on the probability distribution curve of wind pressure extremum.Figure 4
Diagram of probability distribution curve corresponding to guarantee rate. (a) Probability distribution curve with 99% assurance rate. (b) Probability distribution curve with 1% assurance rate.
(a)(b)DefineW0.99,design as design wind load. Here, W0.99,design=maxW1,W2; it can be seen from Figure 4(a) that when guarantee rate P<0.99, wind load W (load effect) corresponding to red and blue curves is both less than or equal to W2; here,W2 is design wind load, and W2=W0.99,design=maxWi,i=1,2.DefineW0.01,design as design wind load. Here, W0.01,design=maxW1,W2, and it can be seen from Figure 4(b) that when guarantee rate P>0.01, wind load W (load effect) corresponding to red and blue curves is both less than or equal to W1; here, W1 is design wind load, and W1=W0.01,design=maxWi,i=1,2.The above analysis shows that if wind loadW≤Wdesign, thenwhenP=0.99, W0.99,design satisfies W0.99,design=maxWi,i=1,2,…N andwhenP=0.01, W0.01,design satisfies W0.01,design=maxWi,i=1,2,…N.DetermineWdesign and enable it satisfying Wdesign=maxW0.99,design,W0.01,design.The symbols in the above formulas are defined the same as those in (4). The word “at least” in the above definition of Wdesignmeans that the area of non-shaded area in Figures 5(a) and 5(b) is greater than 98%. The design extreme wind loads correspond to the quantiles W0.99,design and W0.01,design. Refer to Figures 5(a) and 5(b) below.The range of guarantee rate values in Figures5(a) and 5(b) is described in Tables 2 and 3:Figure 5
Diagram of the relationship between two design wind load:W0.99,design and W0.01,design. (a)
Gaussian distribution case; (b) non-Gaussian distribution case.
(a)(b)Table 2
Guarantee rate range (1).
Range of guarantee rateCondition>98%WhenW0.99,design and W0.01,design are not symmetric with respect to μW=98%WhenW0.99,design and W0.01,design are symmetric with respect to μWTable 3
Guarantee rate range (2).
Range of guarantee rateCondition>98%WhenW0.01,design+2μW−W0.01,design>W0.99,designWhenW0.01,design+2μW−W0.01,design<W0.99,design=98%WhenW0.01,design+2μW−W0.01,design=W0.99,designCombining with (4), it can be inferred that the guarantee rate corresponding to Wdesign is more than 98%. As shown in Figures 5(a) and 5(b), because the distance μW−W0.01,design between the left quantile and the mean μW is not necessarily equal to the distance W0.99,design−μW between the right quantile and the mean μW, this method is called “asymmetric quantile truncation method” with 98% guarantee rate in this case. It is similar to other guarantee rate values.


## 3.3. Numerical Algorithm

### 3.3.1. Monte Carlo Sampling Simulation of Original Wind Speed Data

In the process of extreme value theory modeling and parameter fitting, problems such as missing wind speed sample data and data singularity are often encountered, which will decrease the accuracy of parameter fitting and the estimation of extreme wind speed, and Monte Carlo sampling of original data can solve these problems [29]. Based on the basic principle of Monte Carlo, the specific calculation steps of the Algorithm 1 are as follows:Algorithm 1: Process of original wind speed data.
(1)
The original wind speed sampling data is classified and counted according to the wind direction. There are 6000 wind speed samples at a certain wind direction in this numerical calculation case. The kernel density estimation method is used to fit these 6000 wind speed samples to obtain 6000 probability valuesPi=i=1,…,6000 corresponding to the probability distribution curve.(2)
The coordinates of 6000 points on the probability distribution curve are obtained：U1,P1,U2,P2,…,U6000,P6000.(3)
According to the fitting results calculated in Step 2, a mathematical model can be obtained:P=Fx. Then, U is substituted as an independent variable.(4)
ForP=fx in Step 3, find its inverse function: F−1x.(5)
SubstituteN random numbers P1,P2…,PN generated by uniform distribution U0,1 into F−1Pi=xi.(6)
In this way, the sample of Monte Carlo simulation with sample size ofN is obtained: x1,x2,…,xN. These are the new extreme wind speed samples with sample size of N: U1,U2,…,UN.


### 3.3.2. Numerical Algorithm

Based on model (8) and the above theoretical analysis of the guarantee rate, this paper proposes an algorithm to calculate design wind load according to the wind speed data of wind tunnel test. The algorithm in this paper can be implemented by numerical calculation method, and the specific calculation process can be divided into the following steps:Step 1: Determine nonlinear wind load probability modelW^=1/2U^2C^p (or W^=1/2U^2Cˇp), C^p is maximum value, and Cˇp is minimum value.Step 2: Determine the time interval of two subsamples in the process of determining the probability model. This step is divided into the following two substeps:(1)
The time interval of extreme wind speed subsample is defined as:TU^′=TU¯/m. The length of duration of wind speed time history U¯ is TU¯. The time history of the extreme wind speed U^ (i.e., sub-sample) and its time interval TU^′ are determined.(2)
The time interval of extreme wind pressure coefficient is defined asTCp′=TCp/n. In this paper, the value can be calculated: 10 = 6000/600. The time history of extreme wind pressure coefficient C^pCˇp determined by the wind pressure coefficient time history Cp is 600.Step 3: SinceU^ and C^pCˇp are independent of each other, find the probability distribution that they obey.Step 4: For a certain measuring point, find probability distributionfC^pθi,h,i=1,2,…,N of C^pθi,h. Follow the same step for Cˇpθi,h.Step 5: For a certain measuring point, calculate the converted wind speedU^2 according to height from the original extreme wind speed U^1, and then do the Monte Carlo sampling simulation to find the simulated extreme wind speed U^ out of large sample, i.e., U^θi,h,i=1,2,…,N.Step 6: For a certain measuring point, calculate the probability distributionfU^θi,h,i=1,2,3…N of U^θi,h, and the probability distribution fU^2θi,h, i=1,2,3…N of U^2θi,h.Step 7: For a certain measuring point, find the probability distributionfZ1θi,h,i=1,2,3…N of Z1θi,h=U^2⋅C^p.Step 8: For a certain measuring point, find the probability distributionfZ1θi,h,i=1,2…N of Z1θi,h=1/2ρU^2⋅C^p.Step 9: For a certain measuring point, calculate the design wind loadWj,design, j=1,2,3,…,m with identical guarantee rate, and m is the total number of measuring points.For different measurement points, the same guarantee rate corresponds to different quantile values, so the algorithm is referred to as “based on identical guarantee rate.” The maximum and minimum wind load sequences can be calculated as follows:(1) Calculation of wind load maximum sequence: the wind load corresponding to 99% (or 99.9%) guarantee percentile of each distribution is maximum design wind loadW^0.99θI,h. The wind load value corresponding 1% (or 0.1%) guarantee rate quantile of each distribution is minimum design wind load W^0.01θI,h (negative maximum). If most of the maximum series samples are negative or all of them are negative, the design wind load can be calculated by W^0.01θI,h. For example, the wind load time history is shown in Figure 6.Figure 6
Most of time history are negative, or all of them are negative.(2) Calculation of wind load minimum sequence: the wind load corresponding to 99% (or 99.9%) guarantee percentile of each distribution is maximum design wind loadWˇ0.99θI,h. If most of the samples of the minimum series are positive, or all are positive, the design wind load can be calculated by Wˇ0.99θI,h. For example, the wind load time history is shown in Figure 7.Figure 7
Most of time history are positive, or all of them are positive.The minimum design wind loadW is taken as the value corresponding to the 1% (or 0.1%) guarantee rate quantile of each distribution Wˇ0.99θI,h (negative maximum). Then, the maximum or minimum wind load mentioned above has(5)W^0.99θI,hsatisfying:W^0.99θI,h=maxW^0.99θI,h,I=1,…N,I∈1,N,W^0.01θI,hsatisfying:W^0.01θI,h=maxW^0.01θi,h,i=1,…N,I∈1,N,Wˇ0.99θI,hsatisfying:Wˇ0.99θI,h=maxWˇ0.99θI,h,i=1,…N,I∈1,N,Wˇ0.01θI,hsatisfying:Wˇ0.01θI,h=maxWˇ0.01θI,h,i=1,…N,I∈1,N.According to the wind load maximum sequence and minimum sequence, design wind load of a measuring point is(6)Wj,design=maxW^0.99θI,h,W^0.01θI,h,Wˇ0.99θI,h,Wˇ0.01θI,h,where j=1,2…m,I∈1,N, m is the total number of measuring points, N is the number of wind directions, and it is a positive integer.For each measuring point, 4N wind load values were calculated by the above method based on U^ and C^p (Cˇp) at N wind directions. The maximum absolute value of these 4 N values is the design extreme wind load of this measuring point. The design extreme wind load of 360° in any direction can be obtained by harmonic function fitting. If the layers in height direction and N wind directions in horizontal direction are considered at the same time, this wind load calculation method can be applied to each small block (each measuring point).With the above method, the overall distribution of the maximum wind pressure of the equal probability design of the exterior envelope of the whole building can be obtained.


## 3.3.1. Monte Carlo Sampling Simulation of Original Wind Speed Data

In the process of extreme value theory modeling and parameter fitting, problems such as missing wind speed sample data and data singularity are often encountered, which will decrease the accuracy of parameter fitting and the estimation of extreme wind speed, and Monte Carlo sampling of original data can solve these problems [29]. Based on the basic principle of Monte Carlo, the specific calculation steps of the Algorithm 1 are as follows:Algorithm 1: Process of original wind speed data.
(1)
The original wind speed sampling data is classified and counted according to the wind direction. There are 6000 wind speed samples at a certain wind direction in this numerical calculation case. The kernel density estimation method is used to fit these 6000 wind speed samples to obtain 6000 probability valuesPi=i=1,…,6000 corresponding to the probability distribution curve.(2)
The coordinates of 6000 points on the probability distribution curve are obtained：U1,P1,U2,P2,…,U6000,P6000.(3)
According to the fitting results calculated in Step 2, a mathematical model can be obtained:P=Fx. Then, U is substituted as an independent variable.(4)
ForP=fx in Step 3, find its inverse function: F−1x.(5)
SubstituteN random numbers P1,P2…,PN generated by uniform distribution U0,1 into F−1Pi=xi.(6)
In this way, the sample of Monte Carlo simulation with sample size ofN is obtained: x1,x2,…,xN. These are the new extreme wind speed samples with sample size of N: U1,U2,…,UN.


## 3.3.2. Numerical Algorithm

Based on model (8) and the above theoretical analysis of the guarantee rate, this paper proposes an algorithm to calculate design wind load according to the wind speed data of wind tunnel test. The algorithm in this paper can be implemented by numerical calculation method, and the specific calculation process can be divided into the following steps:Step 1: Determine nonlinear wind load probability modelW^=1/2U^2C^p (or W^=1/2U^2Cˇp), C^p is maximum value, and Cˇp is minimum value.Step 2: Determine the time interval of two subsamples in the process of determining the probability model. This step is divided into the following two substeps:(1)
The time interval of extreme wind speed subsample is defined as:TU^′=TU¯/m. The length of duration of wind speed time history U¯ is TU¯. The time history of the extreme wind speed U^ (i.e., sub-sample) and its time interval TU^′ are determined.(2)
The time interval of extreme wind pressure coefficient is defined asTCp′=TCp/n. In this paper, the value can be calculated: 10 = 6000/600. The time history of extreme wind pressure coefficient C^pCˇp determined by the wind pressure coefficient time history Cp is 600.Step 3: SinceU^ and C^pCˇp are independent of each other, find the probability distribution that they obey.Step 4: For a certain measuring point, find probability distributionfC^pθi,h,i=1,2,…,N of C^pθi,h. Follow the same step for Cˇpθi,h.Step 5: For a certain measuring point, calculate the converted wind speedU^2 according to height from the original extreme wind speed U^1, and then do the Monte Carlo sampling simulation to find the simulated extreme wind speed U^ out of large sample, i.e., U^θi,h,i=1,2,…,N.Step 6: For a certain measuring point, calculate the probability distributionfU^θi,h,i=1,2,3…N of U^θi,h, and the probability distribution fU^2θi,h, i=1,2,3…N of U^2θi,h.Step 7: For a certain measuring point, find the probability distributionfZ1θi,h,i=1,2,3…N of Z1θi,h=U^2⋅C^p.Step 8: For a certain measuring point, find the probability distributionfZ1θi,h,i=1,2…N of Z1θi,h=1/2ρU^2⋅C^p.Step 9: For a certain measuring point, calculate the design wind loadWj,design, j=1,2,3,…,m with identical guarantee rate, and m is the total number of measuring points.For different measurement points, the same guarantee rate corresponds to different quantile values, so the algorithm is referred to as “based on identical guarantee rate.” The maximum and minimum wind load sequences can be calculated as follows:(1) Calculation of wind load maximum sequence: the wind load corresponding to 99% (or 99.9%) guarantee percentile of each distribution is maximum design wind loadW^0.99θI,h. The wind load value corresponding 1% (or 0.1%) guarantee rate quantile of each distribution is minimum design wind load W^0.01θI,h (negative maximum). If most of the maximum series samples are negative or all of them are negative, the design wind load can be calculated by W^0.01θI,h. For example, the wind load time history is shown in Figure 6.Figure 6
Most of time history are negative, or all of them are negative.(2) Calculation of wind load minimum sequence: the wind load corresponding to 99% (or 99.9%) guarantee percentile of each distribution is maximum design wind loadWˇ0.99θI,h. If most of the samples of the minimum series are positive, or all are positive, the design wind load can be calculated by Wˇ0.99θI,h. For example, the wind load time history is shown in Figure 7.Figure 7
Most of time history are positive, or all of them are positive.The minimum design wind loadW is taken as the value corresponding to the 1% (or 0.1%) guarantee rate quantile of each distribution Wˇ0.99θI,h (negative maximum). Then, the maximum or minimum wind load mentioned above has(5)W^0.99θI,hsatisfying:W^0.99θI,h=maxW^0.99θI,h,I=1,…N,I∈1,N,W^0.01θI,hsatisfying:W^0.01θI,h=maxW^0.01θi,h,i=1,…N,I∈1,N,Wˇ0.99θI,hsatisfying:Wˇ0.99θI,h=maxWˇ0.99θI,h,i=1,…N,I∈1,N,Wˇ0.01θI,hsatisfying:Wˇ0.01θI,h=maxWˇ0.01θI,h,i=1,…N,I∈1,N.According to the wind load maximum sequence and minimum sequence, design wind load of a measuring point is(6)Wj,design=maxW^0.99θI,h,W^0.01θI,h,Wˇ0.99θI,h,Wˇ0.01θI,h,where j=1,2…m,I∈1,N, m is the total number of measuring points, N is the number of wind directions, and it is a positive integer.For each measuring point, 4N wind load values were calculated by the above method based on U^ and C^p (Cˇp) at N wind directions. The maximum absolute value of these 4 N values is the design extreme wind load of this measuring point. The design extreme wind load of 360° in any direction can be obtained by harmonic function fitting. If the layers in height direction and N wind directions in horizontal direction are considered at the same time, this wind load calculation method can be applied to each small block (each measuring point).With the above method, the overall distribution of the maximum wind pressure of the equal probability design of the exterior envelope of the whole building can be obtained.


## 4. Research Methods Based on Wind Tunnel Test

### 4.1. Wind Field Simulation

This study takes the Shanghai World Financial Center as an engineering background, and the pressure tests were carried out on a rigid model in the TJ-2 Wind Tunnel at Tongji University. The TJ-2 Wind Tunnel is a boundary layer tunnel of closed-circuit-type, and the working section of the tunnel is 3 m wide, 2.5 m high, and 15 m long. The achievable mean wind speed ranges from 0.5 m/s to 68.0 m/s, adjusted continuously. The mean wind speed profile together with the turbulence intensity for terrain category C was simulated mainly by using hybrid passive devices including roughness blocks, spirelets, and vertical-bar fences. The simulation result and the theoretical formula according to the Load code for the design of building structures by Ministry of Construction of the People’s Republic of China (GB50009－2001) are shown in Figure8.Figure 8
Experimental wind field characteristics. (a) Mean wind speed profile. (b) Turbulence profile of along-wind. (c) Fluctuating wind spectrum of along-wind.
(a)(b)(c)


### 4.2. Test Model

The rigid model for wind tunnel pressure test is shown in Figure9, together with the definitions of model orientation, wind direction angle, and coordinate axis of this rigid model. Each pressure signal was sampled with 6000 data points with a frequency of 312.5 Hz, which means that the sampling time is 19.2 s. Because of the limitations of the channels of the electronic scan valve, the signals of measuring taps could not be collected simultaneously at a time and should be divided into 4 groups of A, B, C, and D. The way of grouping situation, the characteristic size of model section, and the layout of surface pressure measurement taps are shown in Figure 7. It can be seen that pressure measurement taps were distributed on 22 levels, which were numbered from bottom to top of the model. There are 28 to 40 measurement taps in each level, changing with height. In this study, “L-i” is used to indicate the measurement taps, where L represents the number of levels, and i represents the measurement point number at this level. For example, “5–17” represents the 17 measurement points of the 5th level.Figure 9
Setup, facade, wind direction, and coordinate axis of the rigid model. (a) Model test of rigid body pressure measurement. (b) Definition of model orientation, elevation, wind direction, and coordinate axis.
(a)(b)


### 4.3. Test Conditions

Corresponding to the 1:350 length scales, the total height of this model is 1.4 m. In this test, the reference wind speed was 14 m/s, which was monitored by a Pitot tube at 1.2 m height in wind tunnel (above the gradient height, equivalent to the actual 420 m height), about 1.0 m away from the side wall of the wind tunnel. Considering a recurrence period of 50 years, the actual wind speed corresponding to the reference point is 52.93 m/s according to the conversion of basic wind pressure of the C-type landform in Shanghai, and then the wind velocity ratio can be deduced to beλV=0.2645. The time ratio is λt=λL/λV=0.0108, which can be deduced by the dimensional analysis method.


## 4.1. Wind Field Simulation

This study takes the Shanghai World Financial Center as an engineering background, and the pressure tests were carried out on a rigid model in the TJ-2 Wind Tunnel at Tongji University. The TJ-2 Wind Tunnel is a boundary layer tunnel of closed-circuit-type, and the working section of the tunnel is 3 m wide, 2.5 m high, and 15 m long. The achievable mean wind speed ranges from 0.5 m/s to 68.0 m/s, adjusted continuously. The mean wind speed profile together with the turbulence intensity for terrain category C was simulated mainly by using hybrid passive devices including roughness blocks, spirelets, and vertical-bar fences. The simulation result and the theoretical formula according to the Load code for the design of building structures by Ministry of Construction of the People’s Republic of China (GB50009－2001) are shown in Figure8.Figure 8
Experimental wind field characteristics. (a) Mean wind speed profile. (b) Turbulence profile of along-wind. (c) Fluctuating wind spectrum of along-wind.
(a)(b)(c)


## 4.2. Test Model

The rigid model for wind tunnel pressure test is shown in Figure9, together with the definitions of model orientation, wind direction angle, and coordinate axis of this rigid model. Each pressure signal was sampled with 6000 data points with a frequency of 312.5 Hz, which means that the sampling time is 19.2 s. Because of the limitations of the channels of the electronic scan valve, the signals of measuring taps could not be collected simultaneously at a time and should be divided into 4 groups of A, B, C, and D. The way of grouping situation, the characteristic size of model section, and the layout of surface pressure measurement taps are shown in Figure 7. It can be seen that pressure measurement taps were distributed on 22 levels, which were numbered from bottom to top of the model. There are 28 to 40 measurement taps in each level, changing with height. In this study, “L-i” is used to indicate the measurement taps, where L represents the number of levels, and i represents the measurement point number at this level. For example, “5–17” represents the 17 measurement points of the 5th level.Figure 9
Setup, facade, wind direction, and coordinate axis of the rigid model. (a) Model test of rigid body pressure measurement. (b) Definition of model orientation, elevation, wind direction, and coordinate axis.
(a)(b)


## 4.3. Test Conditions

Corresponding to the 1:350 length scales, the total height of this model is 1.4 m. In this test, the reference wind speed was 14 m/s, which was monitored by a Pitot tube at 1.2 m height in wind tunnel (above the gradient height, equivalent to the actual 420 m height), about 1.0 m away from the side wall of the wind tunnel. Considering a recurrence period of 50 years, the actual wind speed corresponding to the reference point is 52.93 m/s according to the conversion of basic wind pressure of the C-type landform in Shanghai, and then the wind velocity ratio can be deduced to beλV=0.2645. The time ratio is λt=λL/λV=0.0108, which can be deduced by the dimensional analysis method.


## 5. Numerical Results and Analysis

### 5.1. Analysis of Calculated Wind Load

In an appearance in Figure10 the Shanghai world financial center is a 57 m × 57 m × 492 m square cylinder, and there is a diagonal line from northeast to southwest at the top. Take one point at each of the two vertical edges from northwest to southeast 52 meters above the ground and cut an oblique section through the diagonal line and the two points, forming the shape of Shanghai World Financial Center. Here, the front elevation in Figure 11 is divided into left elevation I, middle-inclined section 2, and right elevation III. Figure 9(b) is the schematic diagram of front elevation zoning.Figure 10
Layout of pressure measuring points on the front elevation of experimental measuring points (only A and C groups).
(a)(b)(c)Figure 11
Layout of measuring points for pressure measurement model of Shanghai World Financial Center. (a) Side elevation. (b) Front elevation. (c) Cross section on typical height.Through the wind tunnel test described in Section4, the test data are obtained. According to the test data, the numerical method is applied to calculate the wind load with an identical guarantee rate. Based on the wind tunnel experiment described in Section 4.1, the wind speed sample data of each layer are obtained by sampling the wind speed in the wind tunnel. According to these data, the numerical calculation is carried out according to the algorithm implementation steps in Section 3.3, and the wind load and wind pressure are analyzed according to the intermediate and final results of the calculation process. Due to the complex building shape and too many measuring points in this example, in order to grasp the core idea and highlight the trend law, the measuring points on the seventh floor are selected to analyze the calculation results of wind load. The arrangement of measuring points and the relationship between wind direction and angle on the seventh floor are shown in Figure 12, and the calculation results are listed in Table 4.Figure 12
Cross section of the 7th measurement level and definition of wind direction angle.Table 4
Statistical value of wind load on different floors.
Measuring point no.7-17-27-37-47-57-67-77-87–97–107–117–127–137–147–157–167–177–187–19Wind load−0.612−0.769−0.357−0.382−0.445−0.468−0.47−0.563−0.482−0.508−0.444−0.417−0.517−0.593−0.686−0.594−0.53−0.394−0.355Measuring point no.7–207–217–227–237–247–257–267–277–287–297–307–317–327–337–347–357–367–377–38Wind load−0.381−0.487−0.479−0.4880.881.0290.9570.9380.9460.950.960.9650.8952.0970.9941.1621.0081.2681.13According to Figure10(b), the wind load values of the measuring points deployed at 16th, 18th, and 20th floors are negative, and they are shown in Table 5. The results show that the direction of wind load is away from the building surface, which is suction.Table 5
Wind loads at 16th, 18th and 20th floors.
Measuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind load20-1−1.25320-2−1.20620-3−1.12120-4−0.89720-13−0.82720-14−0.96720-15−1.11420-16−1.12218-1−1.2518-2−1.17518-3−1.04418-4−0.87818-13−0.77318-14−0.91618-15−1.01618-16−1.10216-1−0.9616-2−1.07716-3−0.93616-4−0.79216-13−0.72716-14−0.84416-15−0.95216-16−0.497According to Figures9(b) and 10(b), the wind load values of the measuring points deployed at the junction of the left elevation I and the middle-inclined section 2 of the front elevation are shown in Table 6. The results show that the wind loads change greatly at the junction of the two outer inclined planes I and II (see Table 6 and Figure 13(c)), which is mainly due to the influence of incoming flow separation, reattachment, and vortex shedding. In the junction area of two inclined planes, the wind pressure exhibits strong non-Gaussian property, and the wind pressure of the measuring point in the structure surface area with strong non-Gaussian property is more sensitive to the change of wind direction angle and the surface structure of the building object.Table 6
Wind loads at junction of plane I and II.
Measuring point no.28-528-625-525-623-623-722-622-720-620-718-718-816-716-8Wind load1.044−1.4120.959−1.3180.921−1.0910.877−1.0130.819−0.9490.787−0.7870.706−0.701Figure 13
Wind load calculation results analysis.
(a)(b)(c)(d)(e)(f)The wind load of the measuring point at the edge of the building opening is relatively small (see Table7 and Figure 1(d)), and the wind load value at the edge of the high-rise opening is close to or even smaller than that of the lower floor in the same column. It is because the building opening can effectively reduce the wind pressure coefficient.Table 7
Wind load of measuring points near opening.
Measuring point no.39-1939-2039-2139-2239-2339-2439-2539-2639-2732-432-532-632-732-1032-1132-1232-1332-14Wind load−1.345−1.216−1.204−1.283−1.192−1.147−1.123−1.123−1.385−1.239−1.149−1.139−1.102−1.106−1.414−1.478−1.398−1.691The above analysis shows that the wind load calculated by this algorithm is consistent with the actual situation.


### 5.2. Fitting of Probability Distribution Function

#### 5.2.1. Minimal Polynomial Fitting

In this paper, the polynomial fitting method and kernel density estimation method are used to find the probability distribution function of discrete wind speed sample data. In this section, the two algorithms and parameter selection are briefly introduced. Polynomial fitting is to use a polynomial expansion to fit all the observation points in a small analysis area, which contains several points in analysis grid to obtain the objective analysis field of the observation data. The expansion coefficient is determined by least-square fitting. In the numerical experiment of polynomial fitting, it is found that the fitting results do not change with increasing of the order of polynomial, and the fitting result is better when the order of the polynomial is 20. When the order is more or less than 20, the fitting curve begins to oscillate up and down.


#### 5.2.2. Kernel Density Estimation Fitting

Kernel density estimation is a density estimation method proposed in the 1950s and 1960s. It is an effective nonparametric estimation method. Supposing thatx1,x2,…,xN are samples of random variable x, and the probability density function of variable x is fx, then the kernel density estimation of f can be expressed as(7)f^x=1Nh∑i=1NKx−xih,where h is the bandwidth, N is the sample function, K⋅ is the kernel function, the selection principle of K⋅ is Borel measurable in its measure space, and it has the following asymptotically unbiased properties:(8)∫Kudu=1,∫uKudu=0,∫u2Kudu=μ2K>0.It can be proved by (8) that when the number of samplesn⟶∞, the window width h⟶0 and nh⟶0, f^v converges to fv in probability. Therefore, when the sample data is known, the accuracy of f^v depends entirely on the choice of kernel function K⋅ and window width h. When the window width h is fixed, different kernel functions have different effects on f^v. The influence of the two is equivalent. The commonly used kernel functions are Uniform, Triangle, Epanechnikov, and Gaussian. Different kernel functions have different contributions to each sample point to be estimated according to the distance. The choice of kernel function is usually not the key of density estimation, and the bandwidth has a great influence on the smoothness of the model. If the bandwidth is very large, more points will have an impact on the density at x, and the farther points will weaken the weight of the closer points, the deviation of kernel estimation is large, and the variance is small, so the peak values of the main part is covered, and the curve is too smooth. On the contrary, if the bandwidth is too small, the weight of each point has a large drop due to the influence of distance, and the deviation of kernel estimation is small, the variance is large, and the curve presents the characteristics of multipeak instability. In this paper, Gaussian function is chosen as the kernel function.


### 5.3. Analysis of Wind Pressure Distribution

The wind pressure coefficient can be calculated according to the experimental wind pressure sample data. Figure14 shows the average wind pressure coefficient on the surface of the experimental model when the incoming wind direction angle is several typical angles. The mean (cpmean), root mean square (cprms), skewness (cpsk), and kurtosis (cpku) are calculated. In this paper, the wind direction angle 180° is taken as an example. At 180° wind direction angle, the opposite left elevation I is equivalent to that at the wind direction angle of 135o, the inclined section 2 is equivalent to that at the wind direction angle of 180°, and the right elevation III is equivalent to that at the wind direction angle of 225° (135°).Figure 14
Mean wind pressure coefficient.It can be seen from Figure15 that the wind pressure coefficients range from −1 to −0.2 on the left side elevation I. The wind pressure coefficients on the 135° leeward side are negative. The wind pressure coefficients at inclined section 2 range from −1 to −0.7, which is equivalent to a 180° wind angle on the leeward side of the wake area and bears the suction. On the leeward side, the middle and lower parts bear the least suction, and the suction gradually increases spreading around. The deviation at left elevation I ranges from −5 to −1−1. Skewness at inclined section 2 ranges from −2 to −1. The deviations at the right elevation range from −5 to −1. The kurtosis at the left elevation ranges from 0 to 40. The kurtosis at inclined section 2 ranges from 0 to 7. The kurtosis of the right elevation III ranges from 0 to 50.Figure 15
Distribution of statistics of building facade wind pressure coefficient at 180° wind direction.From the experimental data and results, it is known that wind loads at the oblique section2 follow non-Gaussian distribution. The left elevation I and the right elevation III belong to a strong non-Gaussian distribution area. The vortex movement such as the separation and reattachment of the incoming flow on the building surface changes the disorganization and uniformity of the turbulence movement in the original wind field, resulting in a non-Gaussian wind field. The wind pressure on the envelope surface is very sensitive to the change of the incoming flow turbulence intensity, and the non-Gaussian characteristics increase with the increase of the turbulence intensity. The results show that the variation of wind pressure with wind direction angle is more sensitive in the region with strong non-Gaussian characteristics.


### 5.4. Probability Distribution of Fluctuating Wind Pressure

The tail of wind pressure probability distribution directly determines the calculation result of wind pressure extreme value, so the fitting of tail region in probability fitting diagram is the focus of attention. In order to check the fitting of the tail of wind pressure probability distribution more accurately, the longitudinal probability density is set as logarithmic coordinate system, and the time history of wind pressure coefficient on building surface is fitted by normal distribution, Gamma distribution, and Weibull distribution function respectively. Figure16 shows the detailed fitting of wind pressure probability distribution. Figures 16(a) and 16(f) show that, in the region of positive wind pressure, fluctuating wind pressure is in good agreement with Gaussian distribution, and its probability distribution can be well described by normal distribution and Gamma distribution, while Weibull distribution is relatively poor in describing these probability distributions (see Figures 16(a), 16(b), 16(f), and 16(j)). Figures 16(c)–16(e) show that the probability distribution of wind pressure time history at each measuring point in the separation zone of incoming flow on the roof is obviously non-Gaussian and has obvious negative skewness. As for the characteristics of spatial flow field, due to the influence of airflow separation, vortex shedding, and the formation of characteristic turbulence, the correlation between the measuring points in the crosswind area parallel to the incoming flow is strong. On the windward side, the correlation between the measuring points decays rapidly, and the location between the measuring points is weak.Figure 16
Probability density fitting diagram of wind pressure coefficient at typical measurement points. (a) 0° measurement point 28–39. (b) 45° measurement point 28–39. (c) 90° measurement point 28–39. (d) 135° measurement point 28–39. (e) 180° measurement point 28–39. (f) 45° measurement point 5-3. (g) 90° measurement point 5-3. (h) 135° measurement point 5-3.
(a)(b)(c)(d)(e)(f)(g)(h)


### 5.5. Probability Distribution Fitting in the Process of Calculation

In the process of numerical calculation, polynomial fitting and kernel density estimation are used to fit the discrete samples to obtain the probability distribution function of wind pressure extreme value, and the curve consistent with the expectation is obtained. The calculation example used in this paper includes 808 measuring points, for which experiments are carried out at 16 wind directions, and includes two cases of maximum sequence and minimum sequence in the process of calculation. Therefore, a total of808×16×12×2=310272 figures of calculation results can be obtained in the calculation process. This paper takes measuring point 41-1 as an example. The method of calculating the wind load with identical guarantee rate at 0° (due north wind) extreme wind pressure coefficient at measuring point 41-1 obtains the relevant probability density and probability distribution graph. This paper lists 12 graphs of calculating the maximum wind load at 0° at measuring point 41-1 to illustrate the calculation results. As shown in Figures 17(b), 17(d), and 17(f), the fitting probability distribution curve of discrete points is smooth, and the value of probability gradually approaches 1, which is consistent with the distribution of wind load. The polynomial curves of order 20 in Figures 17(i) and 17(j) are completely fitted with the test points, which ensures the accuracy of the algorithm.Figure 17
The probability density and probability distribution diagram of random variables in the intermediate process of the algorithm. (a) Probability density diagram of wind pressure coefficientCp time history of measuring point 41-1 at 0° wind direction. (b) Probability distribution diagram of wind pressure coefficient Cp time history of measuring point 41-1 at 0° wind direction. (c) Probability density diagram of wind pressure coefficient U time history of measuring point 41-1. (d) Probability distribution diagram of wind pressure coefficient U time history of measuring point 41-1 at 0° wind direction. (e) Probability density diagram of U2 time history of measuring point 41-1 at 0° wind direction. (f) Probability distribution diagram of U2 time history of measuring point 41-1 at 0° wind direction. (g) Probability density diagram of U2Cp time history of measuring point 41-1 at 0° wind direction. (h) Probability distribution diagram of U2Cp time history of measuring point 41-1 at 0° wind direction. (i) 20-order polynomial fitting diagram of probability distribution of time history of measuring point 41-1 at 0o wind direction. (j) 20-order polynomial fitting diagram of inverse function of probability distribution of time history of measuring point 41-1 at 0° wind direction. (k) Diagram of superimposed probability density curves of points in top 4 of maximum wind pressure coefficient series (64 curves). (l) Diagram of superimposed probability distribution curves of points in top 4 of maximum wind pressure coefficient series (64 curves).
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)


## 5.1. Analysis of Calculated Wind Load

In an appearance in Figure10 the Shanghai world financial center is a 57 m × 57 m × 492 m square cylinder, and there is a diagonal line from northeast to southwest at the top. Take one point at each of the two vertical edges from northwest to southeast 52 meters above the ground and cut an oblique section through the diagonal line and the two points, forming the shape of Shanghai World Financial Center. Here, the front elevation in Figure 11 is divided into left elevation I, middle-inclined section 2, and right elevation III. Figure 9(b) is the schematic diagram of front elevation zoning.Figure 10
Layout of pressure measuring points on the front elevation of experimental measuring points (only A and C groups).
(a)(b)(c)Figure 11
Layout of measuring points for pressure measurement model of Shanghai World Financial Center. (a) Side elevation. (b) Front elevation. (c) Cross section on typical height.Through the wind tunnel test described in Section4, the test data are obtained. According to the test data, the numerical method is applied to calculate the wind load with an identical guarantee rate. Based on the wind tunnel experiment described in Section 4.1, the wind speed sample data of each layer are obtained by sampling the wind speed in the wind tunnel. According to these data, the numerical calculation is carried out according to the algorithm implementation steps in Section 3.3, and the wind load and wind pressure are analyzed according to the intermediate and final results of the calculation process. Due to the complex building shape and too many measuring points in this example, in order to grasp the core idea and highlight the trend law, the measuring points on the seventh floor are selected to analyze the calculation results of wind load. The arrangement of measuring points and the relationship between wind direction and angle on the seventh floor are shown in Figure 12, and the calculation results are listed in Table 4.Figure 12
Cross section of the 7th measurement level and definition of wind direction angle.Table 4
Statistical value of wind load on different floors.
Measuring point no.7-17-27-37-47-57-67-77-87–97–107–117–127–137–147–157–167–177–187–19Wind load−0.612−0.769−0.357−0.382−0.445−0.468−0.47−0.563−0.482−0.508−0.444−0.417−0.517−0.593−0.686−0.594−0.53−0.394−0.355Measuring point no.7–207–217–227–237–247–257–267–277–287–297–307–317–327–337–347–357–367–377–38Wind load−0.381−0.487−0.479−0.4880.881.0290.9570.9380.9460.950.960.9650.8952.0970.9941.1621.0081.2681.13According to Figure10(b), the wind load values of the measuring points deployed at 16th, 18th, and 20th floors are negative, and they are shown in Table 5. The results show that the direction of wind load is away from the building surface, which is suction.Table 5
Wind loads at 16th, 18th and 20th floors.
Measuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind loadMeasuring point no.Wind load20-1−1.25320-2−1.20620-3−1.12120-4−0.89720-13−0.82720-14−0.96720-15−1.11420-16−1.12218-1−1.2518-2−1.17518-3−1.04418-4−0.87818-13−0.77318-14−0.91618-15−1.01618-16−1.10216-1−0.9616-2−1.07716-3−0.93616-4−0.79216-13−0.72716-14−0.84416-15−0.95216-16−0.497According to Figures9(b) and 10(b), the wind load values of the measuring points deployed at the junction of the left elevation I and the middle-inclined section 2 of the front elevation are shown in Table 6. The results show that the wind loads change greatly at the junction of the two outer inclined planes I and II (see Table 6 and Figure 13(c)), which is mainly due to the influence of incoming flow separation, reattachment, and vortex shedding. In the junction area of two inclined planes, the wind pressure exhibits strong non-Gaussian property, and the wind pressure of the measuring point in the structure surface area with strong non-Gaussian property is more sensitive to the change of wind direction angle and the surface structure of the building object.Table 6
Wind loads at junction of plane I and II.
Measuring point no.28-528-625-525-623-623-722-622-720-620-718-718-816-716-8Wind load1.044−1.4120.959−1.3180.921−1.0910.877−1.0130.819−0.9490.787−0.7870.706−0.701Figure 13
Wind load calculation results analysis.
(a)(b)(c)(d)(e)(f)The wind load of the measuring point at the edge of the building opening is relatively small (see Table7 and Figure 1(d)), and the wind load value at the edge of the high-rise opening is close to or even smaller than that of the lower floor in the same column. It is because the building opening can effectively reduce the wind pressure coefficient.Table 7
Wind load of measuring points near opening.
Measuring point no.39-1939-2039-2139-2239-2339-2439-2539-2639-2732-432-532-632-732-1032-1132-1232-1332-14Wind load−1.345−1.216−1.204−1.283−1.192−1.147−1.123−1.123−1.385−1.239−1.149−1.139−1.102−1.106−1.414−1.478−1.398−1.691The above analysis shows that the wind load calculated by this algorithm is consistent with the actual situation.


## 5.2. Fitting of Probability Distribution Function

### 5.2.1. Minimal Polynomial Fitting

In this paper, the polynomial fitting method and kernel density estimation method are used to find the probability distribution function of discrete wind speed sample data. In this section, the two algorithms and parameter selection are briefly introduced. Polynomial fitting is to use a polynomial expansion to fit all the observation points in a small analysis area, which contains several points in analysis grid to obtain the objective analysis field of the observation data. The expansion coefficient is determined by least-square fitting. In the numerical experiment of polynomial fitting, it is found that the fitting results do not change with increasing of the order of polynomial, and the fitting result is better when the order of the polynomial is 20. When the order is more or less than 20, the fitting curve begins to oscillate up and down.


### 5.2.2. Kernel Density Estimation Fitting

Kernel density estimation is a density estimation method proposed in the 1950s and 1960s. It is an effective nonparametric estimation method. Supposing thatx1,x2,…,xN are samples of random variable x, and the probability density function of variable x is fx, then the kernel density estimation of f can be expressed as(7)f^x=1Nh∑i=1NKx−xih,where h is the bandwidth, N is the sample function, K⋅ is the kernel function, the selection principle of K⋅ is Borel measurable in its measure space, and it has the following asymptotically unbiased properties:(8)∫Kudu=1,∫uKudu=0,∫u2Kudu=μ2K>0.It can be proved by (8) that when the number of samplesn⟶∞, the window width h⟶0 and nh⟶0, f^v converges to fv in probability. Therefore, when the sample data is known, the accuracy of f^v depends entirely on the choice of kernel function K⋅ and window width h. When the window width h is fixed, different kernel functions have different effects on f^v. The influence of the two is equivalent. The commonly used kernel functions are Uniform, Triangle, Epanechnikov, and Gaussian. Different kernel functions have different contributions to each sample point to be estimated according to the distance. The choice of kernel function is usually not the key of density estimation, and the bandwidth has a great influence on the smoothness of the model. If the bandwidth is very large, more points will have an impact on the density at x, and the farther points will weaken the weight of the closer points, the deviation of kernel estimation is large, and the variance is small, so the peak values of the main part is covered, and the curve is too smooth. On the contrary, if the bandwidth is too small, the weight of each point has a large drop due to the influence of distance, and the deviation of kernel estimation is small, the variance is large, and the curve presents the characteristics of multipeak instability. In this paper, Gaussian function is chosen as the kernel function.


## 5.2.1. Minimal Polynomial Fitting

In this paper, the polynomial fitting method and kernel density estimation method are used to find the probability distribution function of discrete wind speed sample data. In this section, the two algorithms and parameter selection are briefly introduced. Polynomial fitting is to use a polynomial expansion to fit all the observation points in a small analysis area, which contains several points in analysis grid to obtain the objective analysis field of the observation data. The expansion coefficient is determined by least-square fitting. In the numerical experiment of polynomial fitting, it is found that the fitting results do not change with increasing of the order of polynomial, and the fitting result is better when the order of the polynomial is 20. When the order is more or less than 20, the fitting curve begins to oscillate up and down.


## 5.2.2. Kernel Density Estimation Fitting

Kernel density estimation is a density estimation method proposed in the 1950s and 1960s. It is an effective nonparametric estimation method. Supposing thatx1,x2,…,xN are samples of random variable x, and the probability density function of variable x is fx, then the kernel density estimation of f can be expressed as(7)f^x=1Nh∑i=1NKx−xih,where h is the bandwidth, N is the sample function, K⋅ is the kernel function, the selection principle of K⋅ is Borel measurable in its measure space, and it has the following asymptotically unbiased properties:(8)∫Kudu=1,∫uKudu=0,∫u2Kudu=μ2K>0.It can be proved by (8) that when the number of samplesn⟶∞, the window width h⟶0 and nh⟶0, f^v converges to fv in probability. Therefore, when the sample data is known, the accuracy of f^v depends entirely on the choice of kernel function K⋅ and window width h. When the window width h is fixed, different kernel functions have different effects on f^v. The influence of the two is equivalent. The commonly used kernel functions are Uniform, Triangle, Epanechnikov, and Gaussian. Different kernel functions have different contributions to each sample point to be estimated according to the distance. The choice of kernel function is usually not the key of density estimation, and the bandwidth has a great influence on the smoothness of the model. If the bandwidth is very large, more points will have an impact on the density at x, and the farther points will weaken the weight of the closer points, the deviation of kernel estimation is large, and the variance is small, so the peak values of the main part is covered, and the curve is too smooth. On the contrary, if the bandwidth is too small, the weight of each point has a large drop due to the influence of distance, and the deviation of kernel estimation is small, the variance is large, and the curve presents the characteristics of multipeak instability. In this paper, Gaussian function is chosen as the kernel function.


## 5.3. Analysis of Wind Pressure Distribution

The wind pressure coefficient can be calculated according to the experimental wind pressure sample data. Figure14 shows the average wind pressure coefficient on the surface of the experimental model when the incoming wind direction angle is several typical angles. The mean (cpmean), root mean square (cprms), skewness (cpsk), and kurtosis (cpku) are calculated. In this paper, the wind direction angle 180° is taken as an example. At 180° wind direction angle, the opposite left elevation I is equivalent to that at the wind direction angle of 135o, the inclined section 2 is equivalent to that at the wind direction angle of 180°, and the right elevation III is equivalent to that at the wind direction angle of 225° (135°).Figure 14
Mean wind pressure coefficient.It can be seen from Figure15 that the wind pressure coefficients range from −1 to −0.2 on the left side elevation I. The wind pressure coefficients on the 135° leeward side are negative. The wind pressure coefficients at inclined section 2 range from −1 to −0.7, which is equivalent to a 180° wind angle on the leeward side of the wake area and bears the suction. On the leeward side, the middle and lower parts bear the least suction, and the suction gradually increases spreading around. The deviation at left elevation I ranges from −5 to −1−1. Skewness at inclined section 2 ranges from −2 to −1. The deviations at the right elevation range from −5 to −1. The kurtosis at the left elevation ranges from 0 to 40. The kurtosis at inclined section 2 ranges from 0 to 7. The kurtosis of the right elevation III ranges from 0 to 50.Figure 15
Distribution of statistics of building facade wind pressure coefficient at 180° wind direction.From the experimental data and results, it is known that wind loads at the oblique section2 follow non-Gaussian distribution. The left elevation I and the right elevation III belong to a strong non-Gaussian distribution area. The vortex movement such as the separation and reattachment of the incoming flow on the building surface changes the disorganization and uniformity of the turbulence movement in the original wind field, resulting in a non-Gaussian wind field. The wind pressure on the envelope surface is very sensitive to the change of the incoming flow turbulence intensity, and the non-Gaussian characteristics increase with the increase of the turbulence intensity. The results show that the variation of wind pressure with wind direction angle is more sensitive in the region with strong non-Gaussian characteristics.


## 5.4. Probability Distribution of Fluctuating Wind Pressure

The tail of wind pressure probability distribution directly determines the calculation result of wind pressure extreme value, so the fitting of tail region in probability fitting diagram is the focus of attention. In order to check the fitting of the tail of wind pressure probability distribution more accurately, the longitudinal probability density is set as logarithmic coordinate system, and the time history of wind pressure coefficient on building surface is fitted by normal distribution, Gamma distribution, and Weibull distribution function respectively. Figure16 shows the detailed fitting of wind pressure probability distribution. Figures 16(a) and 16(f) show that, in the region of positive wind pressure, fluctuating wind pressure is in good agreement with Gaussian distribution, and its probability distribution can be well described by normal distribution and Gamma distribution, while Weibull distribution is relatively poor in describing these probability distributions (see Figures 16(a), 16(b), 16(f), and 16(j)). Figures 16(c)–16(e) show that the probability distribution of wind pressure time history at each measuring point in the separation zone of incoming flow on the roof is obviously non-Gaussian and has obvious negative skewness. As for the characteristics of spatial flow field, due to the influence of airflow separation, vortex shedding, and the formation of characteristic turbulence, the correlation between the measuring points in the crosswind area parallel to the incoming flow is strong. On the windward side, the correlation between the measuring points decays rapidly, and the location between the measuring points is weak.Figure 16
Probability density fitting diagram of wind pressure coefficient at typical measurement points. (a) 0° measurement point 28–39. (b) 45° measurement point 28–39. (c) 90° measurement point 28–39. (d) 135° measurement point 28–39. (e) 180° measurement point 28–39. (f) 45° measurement point 5-3. (g) 90° measurement point 5-3. (h) 135° measurement point 5-3.
(a)(b)(c)(d)(e)(f)(g)(h)


## 5.5. Probability Distribution Fitting in the Process of Calculation

In the process of numerical calculation, polynomial fitting and kernel density estimation are used to fit the discrete samples to obtain the probability distribution function of wind pressure extreme value, and the curve consistent with the expectation is obtained. The calculation example used in this paper includes 808 measuring points, for which experiments are carried out at 16 wind directions, and includes two cases of maximum sequence and minimum sequence in the process of calculation. Therefore, a total of808×16×12×2=310272 figures of calculation results can be obtained in the calculation process. This paper takes measuring point 41-1 as an example. The method of calculating the wind load with identical guarantee rate at 0° (due north wind) extreme wind pressure coefficient at measuring point 41-1 obtains the relevant probability density and probability distribution graph. This paper lists 12 graphs of calculating the maximum wind load at 0° at measuring point 41-1 to illustrate the calculation results. As shown in Figures 17(b), 17(d), and 17(f), the fitting probability distribution curve of discrete points is smooth, and the value of probability gradually approaches 1, which is consistent with the distribution of wind load. The polynomial curves of order 20 in Figures 17(i) and 17(j) are completely fitted with the test points, which ensures the accuracy of the algorithm.Figure 17
The probability density and probability distribution diagram of random variables in the intermediate process of the algorithm. (a) Probability density diagram of wind pressure coefficientCp time history of measuring point 41-1 at 0° wind direction. (b) Probability distribution diagram of wind pressure coefficient Cp time history of measuring point 41-1 at 0° wind direction. (c) Probability density diagram of wind pressure coefficient U time history of measuring point 41-1. (d) Probability distribution diagram of wind pressure coefficient U time history of measuring point 41-1 at 0° wind direction. (e) Probability density diagram of U2 time history of measuring point 41-1 at 0° wind direction. (f) Probability distribution diagram of U2 time history of measuring point 41-1 at 0° wind direction. (g) Probability density diagram of U2Cp time history of measuring point 41-1 at 0° wind direction. (h) Probability distribution diagram of U2Cp time history of measuring point 41-1 at 0° wind direction. (i) 20-order polynomial fitting diagram of probability distribution of time history of measuring point 41-1 at 0o wind direction. (j) 20-order polynomial fitting diagram of inverse function of probability distribution of time history of measuring point 41-1 at 0° wind direction. (k) Diagram of superimposed probability density curves of points in top 4 of maximum wind pressure coefficient series (64 curves). (l) Diagram of superimposed probability distribution curves of points in top 4 of maximum wind pressure coefficient series (64 curves).
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)


## 6. Conclusion

Based on the basic theory of probability theory and the systematic analysis of the surrounding environment and turbulence of buildings, a random variable model is established to calculate the wind load. Based on this model, a numerical method is proposed. The algorithm is verified by the wind tunnel test data, and the following conclusions are obtained:(1)
It is not safe to use peak factor method to calculate extreme wind load of enclosure components in these areas. The investigation of wind disaster shows that it is the local damage caused by the separation of air flow and vortex shedding that leads to the larger damage of the envelope. Based on this, this paper studies the non-Gaussian characteristics of local wind pressure and its extreme value in high-rise buildings, so as to determine the extreme wind load reasonably and accurately.(2)
In this paper, the wind direction effect is considered in calculation method of the extreme wind load with the equivalent guarantee rate. For the windward side, because the wind pressure is close to Gaussian distribution, the extreme wind load based on identical guarantee rate is close to that based on traditional peak factor method. On the crosswind side, because the wind pressure deviates from the Gaussian distribution, the extreme wind load with the identical guarantee rate is generally significantly larger than that based on the traditional peak factor method. On the leeward side, the results are between the windward side and the crosswind side.(3)
In the traditional method of calculating wind load, the “50-year return period” extreme value of the annual extreme value of wind speed at all wind direction and the wind pressure coefficient extreme value calculated by the peak factor method are used. Through the analysis of the test results, the proposed algorithm can reflect the situation of wind field, and its main advantages are accuracy and practicability. This method makes full use of the overall information of wind speed and wind pressure coefficient, so it is a reasonable numerical calculation method.



---

*Source: 1831330-2021-11-28.xml*