# A Conditional Simulation Method for Predicting Wind Pressure Fields of Large-Span Spatial Structures

Authors: Fangjin Sun; Tiantian Liu; Daming Zhang; Zhonghao Xu
Journal: Shock and Vibration
 (2021)

---

## Abstract

Wind load is among the control loads for large-span spatial structures. Wind tunnel test is one of the commonly used methods for measuring wind pressure fields of different kinds of structures. However, due to the limited wind pressure data obtained from wind tunnel testing, it is quite meaningful to employ the limited measured data to predict the unknown wind pressure at target points. Considering the complexity of wind pressure fields of large-span spatial structures, a simplified nonparametric method based on conditional simulation is proposed to predict the unknown pressures using the existing data. The Karhunen–Loève (KL for short) expansion is employed to represent wind pressure random variants as eigenfunctions of the covariance operator. To reduce the variant dimensionality, the nearest neighboring estimator is given for the transition distribution of the KL expansion. The targeted wind pressure fields are obtained by expanding the Fourier basis of the eigenfunction and estimating its expansion coefficients. The proposed method is applied to estimate wind pressures on a gable roof building. The relevant parameters of the wind pressure field are obtained, and the results compare well with those from wind tunnel testing, with higher efficiency. The proposed method effectively reduces the dimensionality of the predicted wind pressures, with reduced errors, higher accuracy, and increased efficiency.


---

## Body

## 1. Introduction

Wind load is one of the control loads for large-span spatial structures. The current structural design specifications [1] are not specified for the calculation of the wind pressure of this type of structure. Usually, the most reliable method for measuring the wind pressure of large-span roofs is to conduct wind tunnel testing. But the surface of large-span roof structures generally needs to be arranged 800 to 1000 measurement points, sometimes even more. However most research institutions rarely have electronic scanner equipment that can measure such a large number of measurement points [2], the measured wind pressure coefficient is often limited. Then, under the conditions of limited and known wind pressure field data, the wind pressure field without test points (locations) is accurately and efficiently simulated; that is, conditional simulation, for accurately predicting the wind-induced response of large-span roofs, has important value and significance for its wind-resistant design.At present, there are mainly two kinds of conditional simulation methods [2]: the Kriging method and the conditional probability density function method. The Kriging method is mainly through the random variable to be estimated expressing as a linear combination of several known reference variables and the corresponding correlation function. The conditional probability density function method was first proposed by Kamenda and Morikawa, applied to earthquake engineering [3]. This method mainly uses the joint probability distribution function between variables to estimate unknown variables. Because the wind pressure field of a large-span roof is more complicated, if the above two methods are directly adopted, it may cause a huge amount of calculation and the solution is easily divergent. In recent years, a technique called “function principal component analysis” [4] can be used to limit the value of random functions in separable Hilbert spaces and perform orthogonal basis functions (which is the eigenfunction containing the covariance operators); the extended method reduces the function dimension, to reduce the amount of calculation. The so-called principal components are those corresponding to the selected first few basis functions. Much of the current research uses the “function principal component analysis” technique to estimate the principal components [5–7], but a few studies use the estimated principal components for prediction. If it is used for prediction, due to the dimension of the random function decreases, it will effectively reduce the amount of calculation and improve the calculation accuracy and efficiency. Domestic scholars mainly use the artificial neural network (ANN) [8, 9] and proper orthogonal decomposition (POD) method [10, 11] to predict the wind pressure field of the structure in terms of wind pressure field simulation. Samaniego and Anitescu explored Deep Neural Networks (DNNs) [12] to solve boundary value problems. The authors considered two approaches: a collocation method and the Deep Energy Method. An important implementation aspect that emerges is that the loss function for the DNN can be naturally devised from the energy of the body. Several very relevant examples from computational mechanics have been solved using DNNs to build the approximation space, which shows that it is possible to tackle the solution of very relevant BVPs using concepts and tools coming from deep learning.Considering the complexity of the random variables of the wind pressure field of large-span spatial structures, this paper proposes a conditional simulation method of wind pressure using the nearest neighbor estimation of the function transition distribution, which is also a nonparametric conditional simulation method. This method first uses Karhunen–Loève (abbreviated as KL) to expand the wind pressure random variable as the eigenfunction of the covariance operator, and then the nearest neighbor estimate of the KL extended transition distribution is given. Then, the expansion coefficient is estimated by extending the Fourier basis of the eigenfunction to realize the estimation of the eigenfunction and obtain the target wind pressure value. Finally, this method is applied to the estimation of the wind pressure field of the large-span structures.


## 2. Conditional Simulation Method

The assumption of Gaussian wind loading is not however always correct as pointed out by previous researchers. In those cases, the conditional simulation of a non-Gaussian wind velocity field is necessary [13]. Combine the conditional simulation technique of Gaussian random fields by Behrouz [14] and the iterative procedure for the unconditional simulation of non-Gaussian random fields [15] to conditionally simulate the time-independent non-Gaussian random fields. Gurley and Kareem [16] developed a procedure for the conditional simulation of multivariate non-Gaussian velocity/pressure fields. For mapping a Gaussian process to a non-Gaussian process and vice versa, they used a modified Hermite transformation using the Hermite polynomial function. The simulated process generated by this method compares well with the target one up to the fourth-order of statistics.Based on the known (i − 1) wind pressure test records Pi−1, the conditional simulation to predict the wind pressure Pi at the ith point can be written as [17](1)Piv=∑n=1i−1∫vΛv,uPnv+Kiv,where u,v∈Φ, Φ is a compact subset of ℜkk≥1; Ki is a random variable whose value is on a separable Hilbert space, and Λi is a sequence of the linear random operator. It is assumed that the wind pressure at each measurement point of i is independent and uncorrelated, and then(2)Λiv,u=∑s∈NWi,sδsvδsu,where Wi,s;s∈N is an independent random variable, δs is an orthogonal real-valued eigenfunction of a set of integral operators, its kernel function is written as(3)C0u,v=∑s∈NγSδSuδsv,γs≥γs+1…,where γs is the real-valued order eigenvalue.Then, in theL2 separable domain of Hilbert space, the wind pressure Pi can be expanded by Karhunen–Loève (KL) as(4)Piv=∑s∈NQi,sδsv,where Qi,s is a series of unrelated random variables, and there are [17](5)Qi,s=∫vPivδsvdv=Pi,δs.


## 3. Conditional Simulation for Wind Pressure Fields of Large-Span Roofs Based on Consistent Estimation

### 3.1. Transition Distribution and Its Nearest Neighboring Estimation

To improve the efficiency and accuracy of wind pressure field simulation, the method of reducing the dimension of wind pressure field variables is employed, which is to map high-dimensional data to low-dimensional space and maintain the properties of the original data. The optimization of these high-dimensional data usually falls into the category of eigenvalue estimation. The KL expansion of random wind pressure field variable is represented by adopting the eigenfunction of covariance operator; then, a simple nonparametric model is proposed. By considering one-dimensional nearest neighboring estimation [18] of arbitrary coefficient transition distribution in KL expansion, wind pressure field variable dimension is reduced. Therefore, the transition distribution of the wind pressure field Piv and its nearest neighboring distribution have to be found firstly, and then eigenfunction will be estimated. Reducing the dimension of the wind pressure field will be discussed from the two aspects, that is, the nearest neighboring estimation and eigenfunction estimation.Estimation of thePivi=1n transition distribution can be written as(6)Unpv|x≔∫ℜ∑s=1Sδ^svqs≤pv∏s=1SdUsnqs|Ks<x,δ^s>,rsa,where δ^s⋅,s=1,…,S is the orthogonal estimator of δs⋅,s=1,…,S, Piv is the wind pressure field, qs is the element of the random variable, and rs is a family radius established around xm=Yn,δ^s for each as s=1,…,S. The equation shows that when estimating the wind pressure field, the number of s one-dimensional estimated value is needed instead of a multidimensional variable, to reduce the calculation time and improve the efficiency. Next is how to construct the one-dimensional estimation. The one-dimensional nearest neighboring estimator of the transition distribution is used to reduce dimensions since the nearest neighboring estimator can maintain the characteristics of the original data well and easy to get.Since general combinations of linear functions are adopted in equation (6), for any Lipschitz function g:ℜ⟶ℜ can be written as(7)∫gpvUnpv|x=∑i=2n∏s=1sQi−1,sn∈Ksxs,n,rs∑i=2nQi−1,sn∈Ksxs,n,rsg∑s=1sQi,snδ^sv,where xs,n≔x,δ^s, rs is a family radius established around xm=Yn,δ^s for each as s=1,…,S.The nearest neighboring estimator is an anticlustering method to reduce the variable dimensions. The basic idea of seeking the nearest neighboring estimator is to calculate the distance between the target center and the neighboring point and control it within the expected value; then, the nearest neighboring estimator ofUsqs|xs can be defined as(8)Usnqs|Ksx,δ^s,rs≔∑i=2nδ^s,Pi−1∈Kx,δ^x,rs,δ^x,Pi≤qs∑i=2nδ^s,Yi−1∈Kx,δ^s,rs,where qs is the element of the random variable and Qi,s is the wind pressure field at the known measuring point.Next, estimating the orthogonal eigenfunction of the covariance operatorδs⋅,s=1,…,S will be discussed in Section 3.2.


### 3.2. Eigenfunction Estimation

It can be seen from equation (4) that the key to simulate the wind pressure field conditions is to seek the integral operator eigenvalue δsv of core function C0u,v in the entire Φ domain. Arbitrary variables are needed.Assumingδ^s⋅,s=1,…,S is the orthogonal estimators of δs⋅,s=1,…,S, then(9)δ^s−δs=Op∏s−1n−τ,where τ>0∏s:s=1,…,S depends on the spacing between eigenvalues.Based on the above conditions, the wind pressure simulation method proposed mainly involves Fourier basis expansion of the eigenfunctions and the expansion coefficient estimation, that is, the estimation of the eigenfunction. Generally, based on the intervalΦ=0,1, the covariance function based on the binary Fourier basis can be written as(10)C0u,v=∑r,l∈Zcr,lexpiru+lv,where cr,l is the r,l-th Fourier coefficient of the kernel function C0u,v. The estimation of the Fourier basis expansion coefficients of the eigenfunctions can be obtained by truncating between −S and S through the above Fourier expansion. Based on the Mercer theory, the Fourier basis expansion coefficient estimation can be written as(11)δs,nSv=∑r=−SSas,nrexpirv,where as,nr is the r-th input of the s-th eigenvector Cns of the Fourier coefficient cnsr,p, and cnsr,p is obtained by the truncated two-dimensional Fourier expansion of the estimated covariance function, an S×S order matrix, expressed as(12)cnSu,v=∑R=−SS∑P=−SScnSr,pexpiru+pv.


### 3.3. Conditional Simulation Procedures for Wind Pressure Prediction of Large-Span Roofs

According to the above, conditional simulation procedures for wind pressure prediction of large-span roofs can be summarized as follows:(1)
According to equations (10) and (11), the first s orthogonal eigenfunctions δs⋅,s=1,…,S of the covariance operator are estimated and to make sure equation (9) is satisfied.(2)
According to wind pressure data measured in the experiment, the estimated eigenfunctionPi,δ^s,i=1,…,n;s=1,…,S is obtained.(3)
Equations (6) and (7) are adopted to calculate the transition distribution of each known wind pressure measuring point, and then equation (8) is used for the nearest neighboring estimator. That is, δs⋅,s=1,…,S is estimated firstly, and the known wind pressure of measuring points is transformed onto two-dimensional space Qi,sni=1n: through Pi,δΛs,i=1,…,n;s=1,…,S.(4)
The estimated value of the random variableQi,sni∈1,…,n;s=1,…,S is obtained by equation (12), using the eigenfunction δ^s:s=1,…,S, where Qi,sn≔Qi,s,n=Pi,δ^s, and finally the wind pressure Pn+1 is obtained by equations (4) and (5).


## 3.1. Transition Distribution and Its Nearest Neighboring Estimation

To improve the efficiency and accuracy of wind pressure field simulation, the method of reducing the dimension of wind pressure field variables is employed, which is to map high-dimensional data to low-dimensional space and maintain the properties of the original data. The optimization of these high-dimensional data usually falls into the category of eigenvalue estimation. The KL expansion of random wind pressure field variable is represented by adopting the eigenfunction of covariance operator; then, a simple nonparametric model is proposed. By considering one-dimensional nearest neighboring estimation [18] of arbitrary coefficient transition distribution in KL expansion, wind pressure field variable dimension is reduced. Therefore, the transition distribution of the wind pressure field Piv and its nearest neighboring distribution have to be found firstly, and then eigenfunction will be estimated. Reducing the dimension of the wind pressure field will be discussed from the two aspects, that is, the nearest neighboring estimation and eigenfunction estimation.Estimation of thePivi=1n transition distribution can be written as(6)Unpv|x≔∫ℜ∑s=1Sδ^svqs≤pv∏s=1SdUsnqs|Ks<x,δ^s>,rsa,where δ^s⋅,s=1,…,S is the orthogonal estimator of δs⋅,s=1,…,S, Piv is the wind pressure field, qs is the element of the random variable, and rs is a family radius established around xm=Yn,δ^s for each as s=1,…,S. The equation shows that when estimating the wind pressure field, the number of s one-dimensional estimated value is needed instead of a multidimensional variable, to reduce the calculation time and improve the efficiency. Next is how to construct the one-dimensional estimation. The one-dimensional nearest neighboring estimator of the transition distribution is used to reduce dimensions since the nearest neighboring estimator can maintain the characteristics of the original data well and easy to get.Since general combinations of linear functions are adopted in equation (6), for any Lipschitz function g:ℜ⟶ℜ can be written as(7)∫gpvUnpv|x=∑i=2n∏s=1sQi−1,sn∈Ksxs,n,rs∑i=2nQi−1,sn∈Ksxs,n,rsg∑s=1sQi,snδ^sv,where xs,n≔x,δ^s, rs is a family radius established around xm=Yn,δ^s for each as s=1,…,S.The nearest neighboring estimator is an anticlustering method to reduce the variable dimensions. The basic idea of seeking the nearest neighboring estimator is to calculate the distance between the target center and the neighboring point and control it within the expected value; then, the nearest neighboring estimator ofUsqs|xs can be defined as(8)Usnqs|Ksx,δ^s,rs≔∑i=2nδ^s,Pi−1∈Kx,δ^x,rs,δ^x,Pi≤qs∑i=2nδ^s,Yi−1∈Kx,δ^s,rs,where qs is the element of the random variable and Qi,s is the wind pressure field at the known measuring point.Next, estimating the orthogonal eigenfunction of the covariance operatorδs⋅,s=1,…,S will be discussed in Section 3.2.


## 3.2. Eigenfunction Estimation

It can be seen from equation (4) that the key to simulate the wind pressure field conditions is to seek the integral operator eigenvalue δsv of core function C0u,v in the entire Φ domain. Arbitrary variables are needed.Assumingδ^s⋅,s=1,…,S is the orthogonal estimators of δs⋅,s=1,…,S, then(9)δ^s−δs=Op∏s−1n−τ,where τ>0∏s:s=1,…,S depends on the spacing between eigenvalues.Based on the above conditions, the wind pressure simulation method proposed mainly involves Fourier basis expansion of the eigenfunctions and the expansion coefficient estimation, that is, the estimation of the eigenfunction. Generally, based on the intervalΦ=0,1, the covariance function based on the binary Fourier basis can be written as(10)C0u,v=∑r,l∈Zcr,lexpiru+lv,where cr,l is the r,l-th Fourier coefficient of the kernel function C0u,v. The estimation of the Fourier basis expansion coefficients of the eigenfunctions can be obtained by truncating between −S and S through the above Fourier expansion. Based on the Mercer theory, the Fourier basis expansion coefficient estimation can be written as(11)δs,nSv=∑r=−SSas,nrexpirv,where as,nr is the r-th input of the s-th eigenvector Cns of the Fourier coefficient cnsr,p, and cnsr,p is obtained by the truncated two-dimensional Fourier expansion of the estimated covariance function, an S×S order matrix, expressed as(12)cnSu,v=∑R=−SS∑P=−SScnSr,pexpiru+pv.


## 3.3. Conditional Simulation Procedures for Wind Pressure Prediction of Large-Span Roofs

According to the above, conditional simulation procedures for wind pressure prediction of large-span roofs can be summarized as follows:(1)
According to equations (10) and (11), the first s orthogonal eigenfunctions δs⋅,s=1,…,S of the covariance operator are estimated and to make sure equation (9) is satisfied.(2)
According to wind pressure data measured in the experiment, the estimated eigenfunctionPi,δ^s,i=1,…,n;s=1,…,S is obtained.(3)
Equations (6) and (7) are adopted to calculate the transition distribution of each known wind pressure measuring point, and then equation (8) is used for the nearest neighboring estimator. That is, δs⋅,s=1,…,S is estimated firstly, and the known wind pressure of measuring points is transformed onto two-dimensional space Qi,sni=1n: through Pi,δΛs,i=1,…,n;s=1,…,S.(4)
The estimated value of the random variableQi,sni∈1,…,n;s=1,…,S is obtained by equation (12), using the eigenfunction δ^s:s=1,…,S, where Qi,sn≔Qi,s,n=Pi,δ^s, and finally the wind pressure Pn+1 is obtained by equations (4) and (5).


## 4. Example Analysis

In this paper, the wind pressure field of a double-sloped large-span roof structure is simulated by the above method. The roof sketch and the arrangement of some measuring points are shown in Figure1 [19]. A total of 84 pressure measuring points are arranged on the roof in the wind tunnel test. In the experiment, the roof is divided into four areas noted 1 to 4 (as shown in Figure 2). According to the division of the roof areas in the experiment [19], the wind pressure coefficients of the measuring points 74–80 (except point 79) in areas 1 and 2 are used as the known data of the conditional simulation in the experiment. The wind pressure-related parameters of the measuring points 79, 81, to 84 in areas 3 and 4 are predicted through the proposed method. The results are compared with the wind tunnel experimental data to validate the correctness of the proposed method.Figure 1
Wind tunnel experiment model and pressure tap layout for a gable roof [19].
(a)(b)Figure 2
Roof divided zones.When using the method introduced in this article to calculate, mainly follow the steps in Section3.3. In the specific implementation process, for the convenience of calculation, first, establish a set of finite grid points for each known wind pressure and sum all the grid points, that is, the integral of equation (7) is replaced by equation (13). That is, for each s, s=1,…,S establishes a circle with a radius rs=csαsn−1/4 around xsn=Pn,δ^s, and cs is an adjustment constant, where cs=1.8[14], S=7. Then, equation (8) is used to establish the nearest neighbor estimation of the transfer distribution for each known wind pressure value, and the wind pressure value of the predicted point can be obtained by replacing Qn+1,ss=1S in the KL expansion. This method is also applicable to other examples of wind pressure prediction. The number of samples can be increased according to the required accuracy. However, the subsequent calculation results of this paper show that even under the condition of limited samples, the calculation accuracy can be guaranteed:(13)∫pvdUnPv|x=∑s=1sQ¯snδ^sv.Among them,(14)Q¯sn=∑i=2nQi,snQi−1,sn∈Ksxs,n,rs∑i=2nQi−i,sn∈Ksxs,n,rs.Table1 shows the average wind pressure coefficient of the predicted point at 0° wind direction angle calculated by this method and compared with the experimental data [19].Table 1
Comparison of mean wind coefficients between conditional simulation and experiment at 0° wind azimuth.
Prediction pointsExperiment valueAnalog valueError (%)79−0.483 3−0.501 83.8381−0.122 4−0.117 34.1782−0.528 4−0.498 75.6283−0.019 0−0.017 95.7984−0.029 4−0.027 85.44Table2 shows the root mean square value of the pulsating wind pressure coefficient at the predicted point at the 0° wind direction angle calculated by the method in this paper and compares it with the experimental data [19].Table 2
Comparison of fluctuating wind coefficients RMS between conditional simulation and experiment at 0° wind azimuth.
Prediction pointsExperiment valueAnalog valueError (%)79−0.154 5−0.164 94.1581−0.097 8−0.095 22.6582−0.172 6−0.180 44.5283−0.038 0−0.040 15.5284−0.033 0−0.031 83.63As shown in Tables1 and 2, the mean wind pressure coefficient and root mean square value of fluctuating wind pressure coefficient of the predicted taps simulated by this method are close to the experimental value, and the error is small, indicating the accuracy of this simulation method is good. It is should be noticed that despite the limited number of samples, the simulation accuracy is still satisfactory. It can be concluded that, if the samples increase, the simulation accuracy will be improved further. To compare the effects of the number of samples on the accuracy of the calculation, all the taps from 73 to 83 are used as the known data for the conditional simulation. The average wind pressure coefficients and the root mean square values of the fluctuating wind pressure coefficients from the conditional simulation of tap 84 are compared with the results when taps 74 to 80 (except for tap 79) are used as conditional simulation samples. It is found that the accuracy of the predicted values is increased by 3% and 5.3%, respectively, while the calculation time increased by 11.9% and 9.6%, respectively. Therefore, while increasing conditional simulation samples to improve the computation accuracy, it is also necessary to consider the computation time cost and strive to make a balance between the two.To illustrate the superiority of the method in the work here, the mean absolute error (MAE) and the mean square error (SAE) of the predicted points at 90° wind direction angle are calculated using the proposed method and compared with the calculation errors of two other methods (LA and PF), commonly used in conditional simulation, the number of simulations in the three methods isn=100, and the results are shown in Table 3. Among them,(15)MAE=1N∑i=1Nεi′SAE=Eεi2Eyi−Ey^i2,where εi is the output error, εi=yi−y^i, i=1,…,N, yi is the experiment wind pressure coefficient, and y^i is the conditional simulation results.Table 3
Comparison of errors by different conditional simulation methods at 90° wind azimuth.
Prediction pointsAdapting methodsAverage wind pressure coefficientPulsating wind pressure coefficientMAESAEMAESAE79LA2.6520.2862.7390.336PF2.2960.2272.5270.277Article method1.0310.1521.2240.14281LA3.5470.3223.3140.353PF2.2340.2692.5310.291Article method1.2270.1261.4710.12682LA2.4590.3323.7530.416PF2.6310.3413.3540.316Article method1.1620.1371.3760.22583LA2.9370.3263.1440.397PF2.3580.4022.5280.326Article method1.2650.1131.1760.20684LA2.9410.4282.7560.376PF2.4380.3312.1730.313Article method1.2960.2011.1430.176The LA method [20] refers to the linear autoregressive estimation method, and the PF method [21] refers to the prediction factor estimation method. The disadvantages of these two methods are that the simulation results are accurate only when the model that generates the simulated data meets the linear assumption. In other words, the accuracy of the simulation results needs to be improved for the presence of nonlinear factors. Considering the nonlinear factors, the kernel function is introduced in the method in this paper. The kernel function is a key factor for implicitly processing nonlinear data in the feature space. The calculation in the feature space will be simple and efficient by using the kernel function.It can be seen from Table3 that the calculation error of this method is much lower than that of the two commonly used conditional simulation methods. Among them, the calculation accuracy of the average absolute error of this method is about 51% higher than that of the LA method, and about 36% higher than that of the PF method. The calculation accuracy of the mean square error is improved by about 58% on average compared with the LA method, and about 47% on average compared with the PF method. It shows that the method in this paper is correct and reliable under different wind direction angles.At the same time, time histories of the predicted target taps are shown with the conditional simulation through wind pressure coefficients of the known pressure taps, as shown in Figure3.Figure 3
Time histories of wind pressure at predicted points. (a) 81 points. (b) 82 points. (c) 83 points. (d) 84 points.
(a)(b)(c)(d)To compare the predicted wind pressure with the experimental data and to validate the accuracy of the conditional simulation, the average error and root mean square error of the simulated average wind pressure coefficients are given, as shown in Figure4.Figure 4
Variation of mean wind pressure coefficient errors. (a) Mean error change. (b) Root mean square error variation.
(a)(b)It can be seen from Figure4 that the mean error and root mean square error of the wind pressure coefficients simulated varied little with time, and the errors are basically below 8%, which validates the accuracy of the simulation.To illustrate the efficiency of the conditional simulation, the residuals of the two commonly used conditional simulation methods are compared, as shown in Figure5.Figure 5
Comparison of residuals by different conditional simulation methods. (a) Calculate the residual variation with time. (b) Calculated residual varies with the number of simulations.
(a)(b)It can be seen from Figure5 that, with the same convergence accuracy achieved, the time and the number of simulations required by the proposed method in this work are less than those of the other two conditional simulation methods, which indicates the efficiency of the proposed method. The reason why the residual decreases with the number of simulations is that as the number of simulations increases, the computation convergence accelerates. This is because the nearest neighbor estimation reduces the originally high-dimensional variables to one-dimensional variables, which improves the calculation efficiency; on the one hand, the use of the kernel function speeds up the calculation of the eigenfunction δs:s=1,…,S, since the use of the kernel function improves the two-dimensional space Qi,sni=1n:s=1,…,S “similarity” identification of the variables, which speeds up the calculation and convergence.


## 5. Conclusions

This paper proposes a nonparametric conditional simulation method of wind pressure using the nearest neighbor estimation of the function transfer distribution. The proposed method is applied to the prediction of wind pressure field on large-span roofs, and the main conclusions obtained are as follows:(1)
The average wind pressure coefficient and pulsating wind pressure coefficient of the large-span roof predicted by the method in this paper are in good agreement with the experimental values. The mean error and root mean square error of the wind pressure coefficient have little change with time and are basically below 8%, which proves that the method proposed in this paper has high accuracy in simulating the wind pressure field.(2)
In this method, the error of the prediction result of the wind pressure field of a large-span roof is much smaller than that of the commonly used conditional simulation method applicable to linear models. The calculation accuracy of the average absolute error of this method is about 51% higher than that of the LA method and approximately 36% higher than the PF method; the calculation accuracy of the mean square error is about 58% higher than the LA method, and about 47% higher than the PF method; and the accuracy of the simulation method is not affected by the number of samples.(3)
The calculation accuracy and convergence speed of the method in this paper are higher than the general conditional simulation method. This is due to the fact that the use of nearest neighbor estimation reduces the variable dimension and the kernel function speeds up the eigenfunction computation.



---

*Source: 8829509-2021-07-07.xml*